Determining the effect of asymmetric data on the variogram. I. Underlying asymmetry

Matheron's usual variogram estimator can result in unreliable variograms when data are strongly asymmetric or skewed. Asymmetry in a distribution can arise from a long tail of values in the underlying process or from outliers that belong to another population that contaminate the primary process. This paper examines the effects of underlying asymmetry on the variogram and on the accuracy of prediction, and the second one examines the effects arising from outliers. Standard geostatistical texts suggest ways of dealing with underlying asymmetry; however, this is based on informed intuition rather than detailed investigation. To determine whether the methods generally used to deal with underlying asymmetry are appropriate, the effects of different coefficients of skewness on the shape of the experimental variogram and on the model parameters were investigated. Simulated annealing was used to create normally distributed random fields of different size from variograms with different nugget:sill ratios. These data were then modified to give different degrees of asymmetry and the experimental variogram was computed in each case. The effects of standard data transformations on the form of the variogram were also investigated. Cross-validation was used to assess quantitatively the performance of the different variogram models for kriging. The results showed that the shape of the variogram was affected by the degree of asymmetry, and that the effect increased as the size of data set decreased. Transformations of the data were more effective in reducing the skewness coefficient in the larger sets of data. Cross-validation confirmed that variogram models from transformed data were more suitable for kriging than were those from the raw asymmetric data. The results of this study have implications for the 'standard best practice' in dealing with asymmetry in data for geostatistical analyses. (C) 2007 Elsevier Ltd. All rights reserved.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - II. Diurnal shape of soil heat flux

For vegetated surfaces, calculation of soil heat flux, G, with the Exact or Analytical method requires a harmonic analysis of below-canopy soil surface temperature, to obtain the shape of the diurnal course of G. When determining G with remote sensing methods, only composite (vegetation plus soil) radiometric brightness temperature is available. This paper presents a simple equation that relates the sum of the harmonic terms derived for the composite radiometric surface temperature to that of belowcanopy soil surface temperature. The thermal inertia, Gamma(,) for which a simple equation has been presented in a companion paper, paper I, is used to set the magnitude of G. To assess the success of the method proposed in this paper for the estimation of the diurnal shape of G, a comparison was made between 'remote' and in situ calculated values from described field sites. This indicated that the proposed method was suitable for the estimation of the shape of G for a variety of vegetation types and densities. The approach outlined in paper I, to obtain Gamma, was then combined with the estimated harmonic terms to predict estimates of G, which were compared to values predicted by empirical remote methods found in the literature. This indicated that the method proposed in the combination of papers I and II gave reliable estimates of G, which, in comparison to the other methods, resulted in more realistic predictions for vegetated surfaces. This set of equations can also be used for bare and sparsely vegetated soils, making it a universally applicable method. (C) 2007 Elsevier B.V. All rights reserved.

Droplet shape of an anisotropic liquid

We investigate how a droplet of a complex liquid is modified by its internal nanoscale structure. As the liquid passes from an isotropic disordered state to an anisotropic layered morphology, the droplet shape switches from a smooth spherical cap to a terraced hyperbolic profile, which can be modeled as a stack of thin concentric circular disks with a repulsion between adjacent disk edges. Our ability to resolve the detailed shape of these defect-free droplets offers a unique opportunity to explore the underlying physics.

Reactive uptake coefficients for heterogeneous reaction of N2O5 with submicron aerosols of NaCl and natural sea salt

The kinetics of uptake of gaseous N2O5 on submicron aerosols containing NaCl and natural sea salt have been investigated in a flow reactor as a function of relative humidity (RH) in the range 30-80% at 295±2K and a total pressure of 1bar. The measured uptake coefficients, γ, were larger on the aerosols containing sea salt compared to those of pure NaCl, and in both cases increased with increasing RH. These observations are explained in terms of the variation in the size of the salt droplets, which leads to a limitation in the uptake rate into small particles. After correction for this effect the uptake coefficients are independent of relative humidity, and agree with those measured previously on larger droplets. A value of γ=0.025 is recommended for the reactive uptake coefficient for N2O5 on deliquesced sea salt droplets at 298K and RH>40%.

Parameter estimation using a particle method: inferring mixing coefficients from sea level observations

This paper presents a first attempt to estimate mixing parameters from sea level observations using a particle method based on importance sampling. The method is applied to an ensemble of 128 members of model simulations with a global ocean general circulation model of high complexity. Idealized twin experiments demonstrate that the method is able to accurately reconstruct mixing parameters from an observed mean sea level field when mixing is assumed to be spatially homogeneous. An experiment with inhomogeneous eddy coefficients fails because of the limited ensemble size. This is overcome by the introduction of local weighting, which is able to capture spatial variations in mixing qualitatively. As the sensitivity of sea level for variations in mixing is higher for low values of mixing coefficients, the method works relatively well in regions of low eddy activity.

Modelling the lava dome extruded at Soufrière Hills Volcano, Montserrat, August 2005–May 2006. Part I: Dome shape and internal structure

Lava domes comprise core, carapace, and clastic talus components. They can grow endogenously by inflation of a core and/or exogenously with the extrusion of shear bounded lobes and whaleback lobes at the surface. Internal structure is paramount in determining the extent to which lava dome growth evolves stably, or conversely the propensity for collapse. The more core lava that exists within a dome, in both relative and absolute terms, the more explosive energy is available, both for large pyroclastic flows following collapse and in particular for lateral blast events following very rapid removal of lateral support to the dome. Knowledge of the location of the core lava within the dome is also relevant for hazard assessment purposes. A spreading toe, or lobe of core lava, over a talus substrate may be both relatively unstable and likely to accelerate to more violent activity during the early phases of a retrogressive collapse. Soufrière Hills Volcano, Montserrat has been erupting since 1995 and has produced numerous lava domes that have undergone repeated collapse events. We consider one continuous dome growth period, from August 2005 to May 2006 that resulted in a dome collapse event on 20th May 2006. The collapse event lasted 3 h, removing the whole dome plus dome remnants from a previous growth period in an unusually violent and rapid collapse event. We use an axisymmetrical computational Finite Element Method model for the growth and evolution of a lava dome. Our model comprises evolving core, carapace and talus components based on axisymmetrical endogenous dome growth, which permits us to model the interface between talus and core. Despite explicitly only modelling axisymmetrical endogenous dome growth our core–talus model simulates many of the observed growth characteristics of the 2005–2006 SHV lava dome well. Further, it is possible for our simulations to replicate large-scale exogenous characteristics when a considerable volume of talus has accumulated around the lower flanks of the dome. Model results suggest that dome core can override talus within a growing dome, potentially generating a region of significant weakness and a potential locus for collapse initiation.

Time-varying coefficients in the almost ideal demand system: An empirical appraisal

This paper provides a generalisation of the structural time series version of the Almost Ideal Demand System (AIDS) that allows for time-varying coefficients (TVC/AIDS) in the presence of cross-equation constraints. An empirical appraisal of the TVC/AIDS is made using a dynamic AIDS with trending intercept as the baseline model with a data set from the Italian Household Budget Survey (1986-2001). The assessment is based on four criteria: adherence to theoretical constraints, statistical diagnostics on residuals, forecasting performance and economic meaningfulness. No clear evidence is found for superior performance of the TVC/AIDS, apart from improved short-term forecasts.

Octanol-water partition coefficients for predicting the effects of tannins in ruminant nutrition

Tannins can cause beneficial or harmful nutritional effects, but their great diversity has until now prevented a rational distinction between tannin structures and their nutritional responses. An attempt has been made to study this problem by examining the octanol-water solubilities of tannins. A relatively simple HPLC method has been developed for screening mixtures of plant tannins for their octanol-water partition coefficients (K-ow coefficients). Tannins were isolated from the fruits and leaves of different Acacia, Calliandra, Dichrostachys, and Piliostigma species, which are known to produce beneficial or harmful effects. The K-ow coefficients of these tannins ranged from 0.061 to 13.9, average coefficients of variation were 9.2% and recoveries were 107%. Acacia nilotica fruits and leaves had the highest K-ow coefficients, that is, 2.0 and 13.9, respectively. These A. nilotica products also have high concentrations of tannins. The combined effects of high octanol solubilities and high tannin concentrations may explain their negative effects on animal nutrition and health. It is known that compounds with high octanol solubilities are more easily absorbed into tissues, and it is, therefore, proposed that such compounds are more likely to cause toxicity problems especially if consumed in large quantities. According to the literature, tannins in human foods tend to have low K-ow coefficients, and this was confirmed for the tannins in Piliostigma thonningii fruits. Therefore, unconventional feeds or browse products should be screened not only for their tannin concentrations but also for low octanol-water partition coefficients in order to identify nutritionally safe feeds and to avoid potentially toxic feeds.

The effects of environmental perturbation and measurement error on estimates of the shape parameter in the theta-logistic model of population regulation

The theta-logistic is a widely used generalisation of the logistic model of regulated biological processes which is used in particular to model population regulation. Then the parameter theta gives the shape of the relationship between per-capita population growth rate and population size. Estimation of theta from population counts is however subject to bias, particularly when there are measurement errors. Here we identify factors disposing towards accurate estimation of theta by simulation of populations regulated according to the theta-logistic model. Factors investigated were measurement error, environmental perturbation and length of time series. Large measurement errors bias estimates of theta towards zero. Where estimated theta is close to zero, the estimated annual return rate may help resolve whether this is due to bias. Environmental perturbations help yield unbiased estimates of theta. Where environmental perturbations are large, estimates of theta are likely to be reliable even when measurement errors are also large. By contrast where the environment is relatively constant, unbiased estimates of theta can only be obtained if populations are counted precisely Our results have practical conclusions for the design of long-term population surveys. Estimation of the precision of population counts would be valuable, and could be achieved in practice by repeating counts in at least some years. Increasing the length of time series beyond ten or 20 years yields only small benefits. if populations are measured with appropriate accuracy, given the level of environmental perturbation, unbiased estimates can be obtained from relatively short censuses. These conclusions are optimistic for estimation of theta. (C) 2008 Elsevier B.V All rights reserved.

Osmotic pressure and aggregate shape in BSA/poly (ethylene glycol)-lipid/Dextran solutions

In this work we study the colloidal osmotic pressure (COP) and aggregate shape in phosphate saline buffer solutions (PH 7.4) containing bovine serum albumin (BSA), poly(ethylene glycol) lipid (PEG(2000)-PE) and Dextran (Dx). Dx was added to the BSA/PEG(2000)-PE system in order to increase the COP of the solution to levels comparable to the COP of healthy adults, with the aim of using the solution as a blood COP regulator. Dynamic light scattering and small angle X-ray scattering results shown the formation of BSA/PEG(2000)-PE/Dx aggregates in the solution. Osmometry results shown that the addition of Dx to the BSA/PE2000-PE system could successfully increase the COP, through the formation of BSA/PEG(2000)-PE/Dx aggregates. The BSA/PEG(2000)-PE/Dx solutions attained COP= 15 mm Hg, representing 60% of COP measured for healthy adults. (c) 2008 Elsevier B.V. All rights reserved.

New magnetic nano-absorbent for the determination of n-octanol/water partition coefficients

A novel and generic miniaturization methodology for the determination of partition coefficient values of organic compounds in noctanol/water by using magnetic nanoparticles is, for the first time, described. We have successfully designed, synthesised and characterised new colloidal stable porous silica-encapsulated magnetic nanoparticles of controlled dimensions. These nanoparticles absorbing a tiny amount of n-octanol in their porous silica over-layer are homogeneously dispersed into a bulk aqueous phase (pH 7.40) containing an organic compound prior to magnetic separation. The small size of the particles and the efficient mixing allow a rapid establishment of the partition equilibrium of the organic compound between the solid supported n-octanol nano-droplets and the bulk aqueous phase. UV-vis spectrophotometry is then applied as a quantitative method to determine the concentration of the organic compound in the aqueous phase both before and after partitioning (after magnetic separation). log D values of organic compounds of pharmaceutical interest (0.65-3.50), determined by this novel methodology, were found to be in excellent agreement with the values measured by the shake-flask method in two independent laboratories, which are also consistent with the literature data. It was also found that this new technique gives a number of advantages such as providing an accurate measurement of log D value, a much shorter experimental time and a smaller sample size required. With this approach, the formation of a problematic emulsion, commonly encountered in shake-flask experiments, is eliminated. It is envisaged that this method could be applicable to the high throughput log D screening of drug candidates. (c) 2005 Elsevier B.V. All rights reserved.

Gas-phase rate coefficients for reactions of NO3, OH, O-3 and O(P-3) with unsaturated alcohols and ethers: correlations and structure-activity relations (SARs)

Experimental difficulties sometimes force modellers to use predicted rate coefficients for reactions of oxygenated volatile organic compounds (oVOCs). We examine here methods for making the predictions for reactions of atmospheric initiators of oxidation, NO3, OH, O-3 and O(P-3), with unsaturated alcohols and ethers. Logarithmic correlations are found between measured rate coefficients and calculated orbital energies, and these correlations may be used directly to estimate rate coefficients for compounds where measurements have not been performed. To provide a shortcut that obviates the need to calculate orbital energies, structure-activity relations (SARs) are developed. Our SARs are tested for predictive power against compounds for which experimental rate coefficients exist, and their accuracy is discussed. Estimated atmospheric lifetimes for oVOCs are presented. The SARs for alkenols successfully predict key rate coefficients, and thus can be used to enhance the scope of atmospheric models incorporating detailed chemistry. SARs for the ethers have more limited applicability, but can still be useful in improving tropospheric models. (C) 2008 Elsevier Ltd. All rights reserved.

Gas-phase rate coefficients for the reactions of NO3, OH and O-3 with alpha,beta-unsaturated esters and ketones: structure-activity relations (SARs)

Gas-phase rate coefficients for the atmospherically important reactions of NO3, OH and O-3 are predicted for 55 alpha,beta-unsaturated esters and ketones. The rate coefficients were calculated using a correlation described previously [Pfrang, C., King, M.D., C. E. Canosa-Mas, C.E., Wayne, R.P., 2006. Atmospheric Environment 40, 1170-1179]. These rate coefficients were used to extend structure-activity relations for predicting the rate coefficients for the reactions of NO3, OH or O-3 with alkenes to include alpha,beta-unsaturated esters and ketones. Conjugation of an alkene with an alpha,beta-keto or alpha,beta-ester group will reduce the value of a rate coefficient by a factor of similar to 110, similar to 2.5 and similar to 12 for reaction with NO3, OH or O-3, respectively. The actual identity of the alkyl group, R, in -C(O)R or -C(O)OR has only a small influence. An assessment of the reliability of the SAR is given that demonstrates that it is useful for reactions involving NO3 and OH, but less valuable for those of O-3 or peroxy nitrate esters. (c) 2006 Elsevier Ltd. All rights reserved.

Designs for an adaptive tuned vibration absorber with variable shape stiffness element

An adaptive tuned vibration absorber (ATVA) with a smart variable stiffness element is capable of retuning itself in response to a time-varying excitation frequency., enabling effective vibration control over a range of frequencies. This paper discusses novel methods of achieving variable stiffness in an ATVA by changing shape, as inspired by biological paradigms. It is shown that considerable variation in the tuned frequency can be achieved by actuating a shape change, provided that this is within the limits of the actuator. A feasible design for such an ATVA is one in which the device offers low resistance to the required shape change actuation while not being restricted to low values of the effective stiffness of the vibration absorber. Three such original designs are identified: (i) A pinned-pinned arch beam with fixed profile of slight curvature and variable preload through an adjustable natural curvature; (ii) a vibration absorber with a stiffness element formed from parallel curved beams of adjustable curvature vibrating longitudinally; (iii) a vibration absorber with a variable geometry linkage as stiffness element. The experimental results from demonstrators based on two of these designs show good correlation with the theory.