Application of a geographically-weighted regression analysis to estimate net primary production of Chinese forest ecosystems

Aim The objective of this paper is to obtain a net primary production (NPP) regression model based on the geographically weighted regression (GWR) method, which includes spatial non‐stationarity in the parameters estimated for forest ecosystems in China. Location We used data across China. Methods We examine the relationships between NPP of Chinese forest ecosystems and environmental variables, specifically altitude, temperature, precipitation and time‐integrated normalized difference vegetation index (TINDVI) based on the ordinary least squares (OLS) regression, the spatial lag model and GWR methods. Results The GWR method made significantly better predictions of NPP in simulations than did OLS, as indicated both by corrected Akaike Information Criterion (AIC_c) and R². GWR provided a value of 4891 for AIC_c and 0.66 for R², compared with 5036 and 0.58, respectively, by OLS. GWR has the potential to reveal local patterns in the spatial distribution of a parameter, which would be ignored by the OLS approach. Furthermore, OLS may provide a false general relationship between spatially non‐stationary variables. Spatial autocorrelation violates a basic assumption of the OLS method. The spatial lag model with the consideration of spatial autocorrelation had improved performance in the NPP simulation as compared with OLS (5001 for AIC_c and 0.60 for R²), but it was still not as good as that via the GWR method. Moreover, statistically significant positive spatial autocorrelation remained in the NPP residuals with the spatial lag model at small spatial scales, while no positive spatial autocorrelation across spatial scales can be found in the GWR residuals. Conclusions We conclude that the regression analysis for Chinese forest NPP with respect to environmental factors and based alternatively on OLS, the spatial lag model, and GWR methods indicated that there was a significant improvement in model performance of GWR over OLS and the spatial lag model.     

GENETIC BASIS OF ACTIVITY METABOLISM. I. INHERITANCE OF SPEED,STAMINA, AND ANTIPREDATOR DISPLAYS IN THE GARTER SNAKE THAMNOPHIS SIRTALIS

Recent conceptual advances in physiological ecology emphasize the potential selective importance of whole-animal performance. Empirical studies of locomotor performance in reptiles have revealed surprising amounts of individual variation in speed and stamina. The present study is the first in a series examining the genetic basis of variation in locomotor performance, activity metabolism, and associated behaviors in garter snakes. Maximal sprint crawling speed, treadmill endurance, and antipredator displays (Arnold and Bennett, 1984; exhibited as snakes reached exhaustion on the treadmill) were measured for approximately six offspring (presumed to be full siblings) from each of 46 wild-caught gravid garter snakes (Thamnophis sirtalis). Each character was measured on two days; all were individually repeatable. Correlations of these characters with body mass, snout–vent length, age at testing, litter size, dam mass, and dam snout–vent length were removed by computing residuals from multiple-regression equations. These residuals were used in subsequent genetic analyses. Approximate coefficients of variation of residuals were 17% for speed, 48% for endurance, and 31% for antipredator displays. Broad-sense heritabilities were significant for all characters: speed h² = 0.58; stamina h² = 0.70; antipredator display h² = 0.42. All three residual characters showed positive and statistically significant phenotypic correlations (r = 0.19–0.36). Genetic correlations (estimated and tested by restricted maximum likelihood) among residuals were positive and highly significant between speed and endurance (0.58), but nonsignificant between speed and antipredator display (0.43), and between endurance and antipredator display (0.26). All environmental correlations were nonsignificant. These data suggest that, contrary to expectations based on previous physiological studies, there may be no necessary evolutionary trade-off between speed and stamina in these animals. This tentative conclusion will have important implications for future theoretical studies of the evolution of locomotor performance and associated antipredator behaviors.     

Effect of natural atmospheric CO2 fertilization suggested by open-grown white spruce in a dry environment

Evidence of an atmospheric CO₂‐fertilization effect on radial growth rates was uncovered for open‐grown white spruce in a mixed‐grass prairie of southwestern Manitoba, Canada. Consistent upward trends of the residuals from dendroclimatic models indicated a decreased ability for climatic parameters to predict radial growth. Despite that a similar amount (61%) of the total variation in radial growth index was explained by climate for both young and old trees, residuals from young trees for the period of 1955–1999 demonstrated a stronger upward trend (R²=0.551, P<0.0001) than old trees for the period of 1900–1996 (R²=0.020, P=0.167). Similar to young trees, the residuals from the early growth period (1900–1929) of old trees also demonstrated a stronger upward trend (R²=0.480, P<0.0001) than the period of 1900–1996. Likewise, a comparable period (1970–1999) of young trees also demonstrated a stronger upward trend (R²=0.619, P<0.0001) than the early growth period (1900–1929) of old trees. In addition, postdrought growth response was much stronger for young trees (1970–1999) compared with old trees at the same development stage (1900–1929) (P=0.011) or within the same time period (1970–1999) (P=0.014). There was no difference (P=0.221) in drought recovery between the early (1900–1929) period and the late (1970–1999) period within old trees. Together, our results suggest that (1) open‐grown white spruce trees improved their growth with time at the early developmental stage, and (2) at the same developmental stage, a greater growth response occurred in the late period when atmospheric CO₂ concentration, and the rate of atmospheric CO₂ increase were both relatively high. While it is impossible to rule out other factors, these results are consistent with expectations for CO₂‐fertilization effects.     

Validation of an Elliptical Anthropometric Model to Estimate Visceral Compartment Area

Objective: The visceral compartment is a surrogate for visceral adipose tissue. Cross‐sectional visceral compartment area (VCA) has been approximated from waist circumference using a circular model. However, the two‐dimensional shape of the abdomen is rarely circular. This study validated an elliptical model of cross‐sectional total abdominal area (TAA), subcutaneous adipose tissue (SAT) area, and VCA at the L₄–L₅ level. Research Methods and Procedures: We analyzed magnetic resonance images (MRIs) at the level of the L₄–L₅ intervertebral space from 35 subjects with a wide range of abdominal adiposity. Waist circumference, abdominal thickness (midline sagittal diameter), abdominal width (coronal diameter at one‐half of abdominal thickness), and abdominal SAT thickness at four sites (front, back, right, and left) were measured from MRI images using an image analysis software. The same anatomical regions were also estimated from anthropometrics purely by geometric formulae of circular and elliptical models. A simple linear regression model was used to interpret the association strength between anthropometric estimates and MRI measures. Results: Estimated TAA by either model was strongly related to MRI TAA (r² = 0.98, p < 0.0001). The SAT and VCA by MRI analysis showed a stronger association with calculation from an elliptical model (r² = 0.95 and 0.88, respectively; p < 0.001) than a circular model (r² = 0.69 and 0.25, respectively; p < 0.001). The absolute prediction residuals and variances were significantly smaller with an elliptical model than a circular model (p < 0.0001). Discussion: An elliptical anthropometric model might be superior to a circular model to estimate abdominal SAT and VCA.     

On the relationship of anthranilic derivatives structure and the FXR (Farnesoid X receptor) agonist activity

Farnesoid X receptor (FXR) is a nuclear receptor related to lipid and glucose homeostasis and is considered an important molecular target to treatment of metabolic diseases as diabetes, dyslipidemia, and liver cancer. Nowadays, there are several FXR agonists reported in the literature and some of it in clinical trials for liver disorders. Herein, a compound series was employed to generate QSAR models to better understand the structural basis for FXR activation by anthranilic acid derivatives (AADs). Furthermore, here we evaluate the inclusion of the standard deviation (SD) of EC₅₀ values in QSAR models quality. Comparison between the use of experimental variance plus average values in model construction with the standard method of model generation that considers only the average values was performed. 2D and 3D QSAR models based on the AAD data set including SD values showed similar molecular interpretation maps and quality (Q²_LOO, Q²_(F2), and Q²_(F3)), when compared to models based only on average values. SD-based models revealed more accurate predictions for the set of test compounds, with lower mean absolute error indices as well as more residuals near zero. Additionally, the visual interpretation of different QSAR approaches agrees with experimental data, highlighting key elements for understanding the biological activity of AADs. The approach using standard deviation values may offer new possibilities for generating more accurate QSAR models based on available experimental data.     

Factors of between- and within-plant distribution of Metopolophium dirhodum (Hom., Aphididae) on small grain cereals

Abstract: The relationship between abundance of rose‐grain aphid Metopolophium dirhodum (Walker) and leaf chlorophyll content of spring wheat, spring oats and winter barley was investigated. Within production stands of each crop 18–25 plots were established, located at places with different plant quality where aphids were counted on particular leaves of 50–100 tillers, and the leaf chlorophyll content and area were determined. In all stands aphid numbers × tiller^?1 increased exponentially with chlorophyll content (r² = 0.783–0.933). This parallel variation may be explained by increased nitrogen content and assimilate production of vigorous chlorophyll‐rich plants. The aphid numbers on leaves of particular order (within‐plant distribution) were also correlated with leaf chlorophyll content (r² = 0.373–0.827). However, in oats and barley the analysis of variance of residuals of log aphids × leaf^?1 versus leaf chlorophyll regression revealed a significant effect of leaf position (the order of the leaf from the top of the plant). The magnitude of residuals was positively related to leaf size and may be related to the intensity of phloem transport from the source leaves to sink organs.     

Stability of Adiposity Phenotypes from Childhood and Adolescence into Young Adulthood with Contribution of Parental Measures

Objective: The stability of several indicators of body composition and adipose tissue distribution over 12 years was quantified. Research Methods and Procedures: The participants were 77 boys and 76 girls who were evaluated along with their parents at baseline as children and adolescents (8 to 18 years of age) and remeasured as young adults 12 years later. Indicators of body composition included the body mass index, fat mass, fat free mass, percentage of body fat, sum of six skinfolds (SF6), and the first principal component of six age‐adjusted skinfold residuals. Relative adipose tissue distribution was represented by the second principal component of skinfold residuals and a trunk‐to‐extremity skinfold ratio, adjusted for SF6. Results: Partial interage correlations, controlling for initial age and length of follow‐up, were 0.65 and 0.59 for the body mass index, 0.59 and 0.64 for fat mass, 0.65 and 0.57 for fat free mass, 0.50 and 0.57 for percentage of body fat, 0.66 and 0.44 for SF6, 0.64 and 0.42 for the first principal component of six age‐adjusted skinfold residuals, 0.19 and 0.31 for the second principal component of skinfold residuals, and 0.41 and 0.47 for trunk‐to‐extremity skinfold ratio, adjusted for SF6, in men and women, respectively. Multiple regression analyses indicated that the significant partial R² values of parental measurements on the prediction of their offspring in young adulthood ranged from 2% to 9%. Discussion: The results indicate moderately high stability of indicators of body composition and somewhat lower stability of measures of adipose tissue distribution. Overall, parental measures offer less predictive value than do measures of childhood and adolescent body composition and adipose tissue distribution.     

Detecting biodiversity hotspots using species–area and endemics–area relationships: the case of butterflies

Using data of the Red Data Book of European Butterflies we establish the species–area relationship (S = 8.5 A^0.23) and the endemics–area relationship (S = 0.5 A^0.18) of European butterflies. Applying confidence limits as tools for the identification of hotspot countries we show that in the case of butterflies hotspots of endemism and hotspots of overall species richness do not coincide. We introduce plots of residuals from species–area relationships shifted upwards by the 95% confidence limits of the intercept (α-values) as a new tool for identifying and ranking of hotspots.     

The cost of living slowly: metabolism,Q10 and repeatability in a South American harvestman

Abstract Resources are usually considered to be limited in caves and underground habitats. Therefore, lower metabolic rates of ectotherms in these environments should be advantageous. The standard metabolic rate (SMR) of the common harvestman Pachylus paessleri Roewer, 1913 (Arachnida, Opiliones) is determined at two ambient temperatures. Repeatability of SMR is estimated by two methods: (i) product–moment correlation on residuals of body mass and (ii) variance components. Estimations of the thermal sensitivity of metabolic rate (individual Q₁₀). SMR of P. paessleri at 20 °C show a mean of 19.01 ± 6.04 and 14.99 ± 3.64 μL CO₂ h^?1 for males and females, respectively. Thermal sensitivity of SMR is 2.11 ± 0.23 with a substantial coefficient of variation of 26.7%. In males, residuals of CO₂ release are significantly repeatable (r_p = 0.61) between measurements realized 5 months apart, which reflects the consistency of this trait over time. As typical soil inhabitants, the harvestman P. paessleri present a lower metabolic rate in comparison with arthropods of similar body mass (e.g. arachnids). Although the coefficient of thermal sensitivity is within of the range reported for arthropods inhabiting Mediterranean zones, it is lower than expected for an organism with such a low SMR. This appears to be the first report of repeatability of metabolic rate in a harvestman species.     

Validation of a DNA methylation reference panel for the estimation of nucleated cells types in cord blood

Cord blood is widely used as surrogate tissue in epigenome-wide association studies of prenatal conditions. Cell type composition variation across samples can be an important confounder of epigenome-wide association studies in blood that constitute a mixture of cells. We evaluated a newly developed cord blood reference panel to impute cell type composition from DNA methylation levels, including nucleated red blood cells (nRBCs). We estimated cell type composition from 154 unique cord blood samples with available DNA methylation data as well as direct measurements of nucleated cell types. We observed high correlations between the estimated and measured composition for nRBCs (r = 0.92, R² = 0.85), lymphocytes (r = 0.77, R² = 0.58), and granulocytes (r = 0.72, R² = 0.52), and a moderate correlation for monocytes (r = 0.51, R² = 0.25) as well as relatively low root mean square errors from the residuals ranging from 1.4 to 5.4%. These results validate the use of the cord blood reference panel and highlight its utility and limitations for epidemiological studies.     

Reproductive allocation in Daphnia exposed to toxic cyanobacteria

We investigated experimentally how resources were allocated to reproduction in Daphnia pulex and Daphnia longispina when varying levels of toxic Microcystis were added to higher quality food. We used multiple regression models to estimate mean offspring size and clutch size in relation to maternal size and clutch number, and analysed effects of treatments on residuals from the models. We also measured variation in per offspring investment. At a high cyanobacterial level, D.pulex was virtually unable to reproduce. At a lower level, D.pulex produced small clutches. However, the regression model residuals indicated that the presence of cyanobacteria increased the portion of available resources allocated to reproduction. The observed allocation may be a means to maximize reproduction under diminished longevity. Effects on mean offspring size were marginal in D.pulex but variation in per offspring investment sometimes decreased in cyanobacterial exposures. Daphnia longispina was affected by a higher cyanobacterial level only, where offspring sized was reduced. Deviations from the regression model indicated that effects on maternal size alone do not explain this effect. Clutch size residuals and per offspring investment were unaffected by treatments in D.longispina. The observed responses differ from theoretical models on reproductive allocation under food imitation.      

The use of model-fitting in the interpretation of 'dual' uptake isotherms

Abstract. Published data of the concentration dependence of the uptake rate (uptake isotherms) of K⁺ Na⁺, Cl^?, SO2^?₄, and L-lysine in barley roots, and glucose and 3-O-methylglucose in potato tuber tissue, were re-examined. In as much as these isotherms yield non-linear, concave upward Eadie-Hofstee plots, they might have been termed ‘dual’ isotherms. In addition, all these isotherms have been considered to display discontinuous transitions in gradient. The following models that yield continuous isotherms were fitted to the isotherms: (1) the sum of a Michaelis-Menten term and a linear term; (2) the sum of two Michaelis-Menten terms; (3) the sum of two Michaelis-Menten terms and a linear one. Goodness of fit was judged from: (i) the weighted mean square of deviates; (ii) the standard errors of the kinetic parameters; (iii) the algebraic significance of the terms; (iv) a Rankits plot of the residuals; (v) a Runs test on the residuals. For the precise and detailed isotherms of SO^2? uptake, only model (3) gave a fit that was satisfactory in all respects. There appeared to be no reason to consider these isotherms as multiphasic. The same conclusion was reached for the L-lysine uptake isotherms. For the other isotherms the results were less conclusive. Thai for K⁺ and Na+ could, at any rate, be described satisfactorily by a continuous model, the best fit being obtained with model (2). The uptake isotherms of Cl^? and 3-0-methylglucose could be best described by model (2), and that of glucose by model (3), only the result of the Runs test being unsatisfactory. It is concluded that there is hardly any evidence that the presumed ‘jumps’ or discontinuities or inflections in the gradient of uptake isotherms are not due to experimental error in the data. It is suggested that many uptake isotherms may be described by model (3), although the reason for this is still incompletely understood.     

Modeling and Predicting Hemorrhagic Fever with Renal Syndrome Trends Based on Meteorological Factors in Hu County,China

Background

Hu County is a serious hemorrhagic fever with renal syndrome (HFRS) epidemic area, with notable fluctuation of the HFRS epidemic in recent years. This study aimed to explore the optimal model for HFRS epidemic prediction in Hu.

Methods

Three models were constructed and compared, including a generalized linear model (GLM), a generalized additive model (GAM), and a principal components regression model (PCRM). The fitting and predictive adjusted R² of each model were calculated. Ljung-Box Q tests for fitted and predicted residuals of each model were conducted. The study period was stratified into before (1971–1993) and after (1994–2012) vaccine implementation epochs to avoid the confounding factor of vaccination.

Results

The autocorrelation of fitted and predicted residuals of the GAM in the two epochs were not significant (Ljung-Box Q test, P>.05). The adjusted R² for the predictive abilities of the GLM, GAM, and PCRM were 0.752, 0.799, and 0.665 in the early epoch, and 0.669, 0.756, and 0.574 in the recent epoch. The adjusted R² values of the three models were lower in the early epoch than in the recent epoch.

Conclusions

GAM is superior to GLM and PCRM for monthly HFRS case number prediction in Hu County. A shift in model reliability coincident with vaccination implementation demonstrates the importance of vaccination in HFRS control and prevention.     

Changes in stomatal conductance and net photosynthesis during phenological development in spring wheat: implications for gas exchange modelling

Gas exchange was measured from 1 month before the onset of anthesis until the end of grain filling in field-grown spring wheat, Triticum aestivum L., cv. Vinjett, in southern Sweden. Two g _s models were parameterised using these data: one Jarvis-type multiplicative g _s model (J-model), and one combined stomatal-photosynthesis model (L-model). In addition, the multiplicative g _s model parameterisation for wheat used within the European Monitoring and Evaluation Programme (EMEP-model) was tested and evaluated. The J-model performed well (R ²=0.77), with no systematic pattern of the residuals plotted against the driving variables. The L-model explained a larger proportion of the variation in g _s data when observations of A _n were used as input data (R ²=0.71) compared to when A _n was modelled (R ²=0.53). In both cases there was a systematic model failure, with g _s being over- and underestimated before and after anthesis, respectively. This pattern was caused by the non-parallel changes in g _s and A _n during plant phenological development, with A _n both peaking and starting to decline earlier as compared to g _s . The EMEP-model accounted for 41% of the variation in g _s data, with g _s being underestimated after anthesis. We conclude that, under the climatic conditions prevailing in southern Scandinavia, the performance of the combined stomatal-photosynthesis approach is hampered by the non-parallel changes in g _s and A _n, and that the phenology function of the EMEP-model, having a sharp local maximum at anthesis, should be replaced by a function with a broad non-limiting period after anthesis.     

Application of artificial neural network model for the development of optimized complex medium for phenol degradation using <Emphasis Type="Italic">Pseudomonas pictorum</Emphasis> (NICM 2074)

Biodegradation of phenol using Pseudomonas pictorum (NICM 2074) a potential biodegradant of phenol was investigated for its degrading potential under different operating conditions. The neural network input parameter set consisted of the same set of four levels of maltose (0.025, 0.05, 0.075 g/l), phosphate (3, 12.5, 22 g/l), pH (7, 8, 9) and temperature (30°C, 32°C, 34°C) on phenol degradation was investigated and a Artificial Neural Network (ANN) model was developed to predict the extent of degradation. The learning, recall and generalization characteristic of neural networks was studied using phenol degradation system data. The efficiency of the model generated by the ANN, was tested and compared with the results obtained from an established second order polynomial multiple regression analysis (MRA). Further, the two models (ANN and MRA) were used to predict the percentage of degradation of phenol for blind test data. Performance of both the models were validated in the cases of training and test data, ANN was recommended based on the following higher coefficient of determination R ²; lower standard error of residuals and lower mean absolute percentage deviation.     

基于不同空间模型的马尾松林生态系统碳密度研究

以江西省马尾松林生态系统为研究对象,基于样地调查及样品碳含量测定结果计算其碳密度,并选取立地、植被及气象等方面的15个因子,采用多元线性逐步回归方法筛选出对生态系统碳密度影响显著的因子,然后分别利用最小二乘模型(OLS)、空间误差模型(SEM)、空间滞后模型(SLM)和地理加权回归模型(GWR)构建生态系统碳密度与其影响因子之间的关系模型,筛选出最优的拟合模型。结果表明:对马尾松林生态系统碳密度影响显著的因子分别为海拔、坡度、土层厚度、胸径、年均温度和年均降水量。4种模型拟合结果均显示碳密度与坡度呈负相关,与海拔、土层厚度、胸径呈正相关。模型的决定系数(R~2)由大到小分别为GWR(0.8043)SEM(0.6371)SLM(0.6364)OLS(0.6321),模型均方误差(MSE)与赤池信息准则(AIC)最大的均为OLS模型,最小的均为GWR模型;残差检验表明GWR模型能有效降低模型残差的空间自相关性。综合分析得出GWR模型的拟合效果最优,更适用于江西省马尾松林生态系统碳密度的估测。     

单木生物量模型估计区域尺度生物量的不确定性

采用系统抽样体系江西省固定样地杉木连续观测数据和生物量数据,通过Monte Carlo法反复模拟由单木生物量模型推算区域尺度地上生物量的过程,估计了江西省杉木地上总生物量。基于不同水平建模样本量n及不同决定系数R~2的设计,分别研究了单木生物量模型参数变异性及模型残差变异性对区域尺度生物量估计不确定性的影响。研究结果表明:2009年江西省杉木地上生物量估计值为(19.84±1.27)t/hm~2,不确定性占生物量估计值约6.41%。生物量估计值和不确定性值达到平稳状态所需的运算时间随建模样本量及决定系数R~2的增大而减小;相对于模型参数变异性,残差变异性对不确定性的影响更小。     

Applying Canonical Analysis to a Kind of Generalized Regression

This paper is concerned with the generalized model E(φ(Y)! X) =φ(X) involving the transformations on both the predictor vector X and the response variable Y. For this purpose, Taylor expansions and canonical analysis are applied. For optimizing the expansions, it is shown by a simulation study that not only prediction error, the combination of model error and noise error, is an important index, but the distribution of the residuals and the t-values of the coefficients also must be considered. Furthermore, the results of penicillin titrition show that the practical situations often need to be considered in selecting an appropriate model for a real-life problem.     

Comparative histochemical analysis of intestinal glycoconjugates in the blunthead pufferfish Sphoeroides pachygaster and grey triggerfish Balistes capriscus (Teleostei: Tetraodontiformes)

The localization of intestinal glycoconjugates of the blunthead pufferfish Sphoeroides pachygaster and the grey triggerfish Balistes capriscus from the north-western Ionian Sea was analysed by histochemical methods (PAS, AB pH 2.5, HID) and lectin binding experiments (WGA, LFA, SBA, sialidase-SBA, PNA, sialidase-PNA, ConA, AAA, UEA-I, LTA) to assess how evolutionary loss of a functional stomach in S. pachygaster affects intestinal secretions relative to the B. capriscus, which retains the plesiomorphic gastric condition. Sphoeroides pachygaster had a lower content of acid mucins but more complex sialylation patterns than B. capriscus. GalNAc and GlcNAc residuals were present in both, but GalNAc residuals in S. pachygaster were subterminal to sialic acid. Balistes capriscus lacked galactosylated residuals and its enterocytes had a glycocalyx that differed in composition between the small intestine and the rectum and was missing from S. pachygaster. Functional and ecological implications of these findings are discussed.     

Estimation of Differently Bound Water Molecules for the Gel Phase of Dimyristoylphosphatidylethanolamine-Water System as Studied by DSC and 2H-NMR Spectroscopy

A measurement of ²H spin-lattice relaxation time, T ₁, forD₂O was performed with a high resolution liquid NMR apparatus fortwo samples of dimyristoylphosphatidylethanolamine (DMPE)-D₂Osystem in a full hydration at varying temperatures of –20, –10, and 5 ^°C, and both components and compositions of differently boundfreezable water molecules were estimated from a best-fitted curve toexperimental inversion recovery data. A choice of the best-fitted curve wasbased on a distribution of weighted residuals for the experimental data. Asingle component was found for a temperature of –20 ^°C. At 5 ^°C, where all the freezable water exists in the liquid state, threecomponents were observed to be characterized by T ₁ values ofapproximately 20, 100, and 200 ms, respectively. By comparingcompositions of these individual components with those obtained in ourprevious DSC study, it was revealed that the first and secondarycomponents are members of freezable interlamellar water and the last oneis comparable to bulk water.