Comparing sampling needs for variograms of soil properties computed by the method of moments and residual maximum likelihood

It has been generally accepted that the method of moments (MoM) variogram, which has been widely applied in soil science, requires about 100 sites at an appropriate interval apart to describe the variation adequately. This sample size is often larger than can be afforded for soil surveys of agricultural fields or contaminated sites. Furthermore, it might be a much larger sample size than is needed where the scale of variation is large. A possible alternative in such situations is the residual maximum likelihood (REML) variogram because fewer data appear to be required. The REML method is parametric and is considered reliable where there is trend in the data because it is based on generalized increments that filter trend out and only the covariance parameters are estimated. Previous research has suggested that fewer data are needed to compute a reliable variogram using a maximum likelihood approach such as REML, however, the results can vary according to the nature of the spatial variation. There remain issues to examine: how many fewer data can be used, how should the sampling sites be distributed over the site of interest, and how do different degrees of spatial variation affect the data requirements? The soil of four field sites of different size, physiography, parent material and soil type was sampled intensively, and MoM and REML variograms were calculated for clay content. The data were then sub-sampled to give different sample sizes and distributions of sites and the variograms were computed again. The model parameters for the sets of variograms for each site were used for cross-validation. Predictions based on REML variograms were generally more accurate than those from MoM variograms with fewer than 100 sampling sites. A sample size of around 50 sites at an appropriate distance apart, possibly determined from variograms of ancillary data, appears adequate to compute REML variograms for kriging soil properties for precision agriculture and contaminated sites. (C) 2007 Elsevier B.V. All rights reserved.

Variation in dynamic elastic shear modulus of sandstone upon fluid saturation and substitution

Experimental acoustic measurements on sandstone rocks at both sonic and ultrasonic frequencies show that fluid saturation can cause a noticeable change in both the dynamic bulk and shear elastic moduli of sandstones. We observed that the change in dynamic shear modulus upon fluid saturation is highly dependent on the type of saturant, its viscosity, rock microstructure, and applied pressures. Frequency dispersion has some influence on dynamic elastic moduli too, but its effect is limited to the ultrasonic frequency ranges and above. We propose that viscous coupling, reduction in free surface energy, and, to a limited extent, frequency dispersion due to both local and global flow are the main mechanisms responsible for the change in dynamic shear elastic modulus upon fluid saturation and substitution, and we quantify influences.

Estimating the spatial scale of herbicide and soil interactions by nested sampling, hierarchical analysis of variance and residual maximum likelihood

An unbalanced nested sampling design was used to investigate the spatial scale of soil and herbicide interactions at the field scale. A hierarchical analysis of variance based on residual maximum likelihood (REML) was used to analyse the data and provide a first estimate of the variogram. Soil samples were taken at 108 locations at a range of separating distances in a 9 ha field to explore small and medium scale spatial variation. Soil organic matter content, pH, particle size distribution, microbial biomass and the degradation and sorption of the herbicide, isoproturon, were determined for each soil sample. A large proportion of the spatial variation in isoproturon degradation and sorption occurred at sampling intervals less than 60 m, however, the sampling design did not resolve the variation present at scales greater than this. A sampling interval of 20-25 m should ensure that the main spatial structures are identified for isoproturon degradation rate and sorption without too great a loss of information in this field.

Estimating the spatial scales of regionalized variables by nested sampling, hierarchical analysis of variance and residual maximum likelihood

The variogram is essential for local estimation and mapping of any variable by kriging. The variogram itself must usually be estimated from sample data. The sampling density is a compromise between precision and cost, but it must be sufficiently dense to encompass the principal spatial sources of variance. A nested, multi-stage, sampling with separating distances increasing in geometric progression from stage to stage will do that. The data may then be analyzed by a hierarchical analysis of variance to estimate the components of variance for every stage, and hence lag. By accumulating the components starting from the shortest lag one obtains a rough variogram for modest effort. For balanced designs the analysis of variance is optimal; for unbalanced ones, however, these estimators are not necessarily the best, and the analysis by residual maximum likelihood (REML) will usually be preferable. The paper summarizes the underlying theory and illustrates its application with data from three surveys, one in which the design had four stages and was balanced and two implemented with unbalanced designs to economize when there were more stages. A Fortran program is available for the analysis of variance, and code for the REML analysis is listed in the paper. (c) 2005 Elsevier Ltd. All rights reserved.

The rise and fall of open solar flux during the current grand solar maximum

We use geomagnetic activity data to study the rise and fall over the past century of the solar wind flow speed VSW, the interplanetary magnetic field strength B, and the open solar flux FS. Our estimates include allowance for the kinematic effect of longitudinal structure in the solar wind flow speed. As well as solar cycle variations, all three parameters show a long-term rise during the first half of the 20th century followed by peaks around 1955 and 1986 and then a recent decline. Cosmogenic isotope data reveal that this constitutes a grand maximum of solar activity which began in 1920, using the definition that such grand maxima are when 25-year averages of the heliospheric modulation potential exceeds 600 MV. Extrapolating the linear declines seen in all three parameters since 1985, yields predictions that the grand maximum will end in the years 2013, 2014, or 2027 using VSW, FS, or B, respectively. These estimates are consistent with predictions based on the probability distribution of the durations of past grand solar maxima seen in cosmogenic isotope data. The data contradict any suggestions of a floor to the open solar flux: we show that the solar minimum open solar flux, kinematically corrected to allow for the excess flux effect, has halved over the past two solar cycles.

Engineering Pt in ceria for a maximum metal-support interaction in catalysis

Conventional supported metal catalysts are metal nanoparticles deposited on high surface area oxide supports with a poorly defined metal−support interface. Typically, the traditionally prepared Pt/ceria catalyzes both methanation (H2/CO to CH4) and water−gas shift (CO/H2O to CO2/H2) reactions. By using simple nanochemistry techniques, we show for the first time that Pt or PtAu metal can be created inside each CeO2 particle with tailored dimensions. The encapsulated metal is shown to interact with the thin CeO2 overlayer in each single particle in an optimum geometry to create a unique interface, giving high activity and excellent selectivity for the water−gas shift reaction, but is totally inert for methanation. Thus, this work clearly demonstrates the significance of nanoengineering of a single catalyst particle by a bottom-up construction approach in modern catalyst design which could enable exploitation of catalyst site differentiation, leading to new catalytic properties.

Young's modulus varies with differential orientation of keratin in feathers

Feathers are composed of a structure that, whilst being very light, is able to withstand the large aerodynamic forces exerted upon them during flight. To explore the contribution of molecular orientation to feather keratin mechanical properties, we have examined the nanoscopic organisation of the keratin molecules by X-ray diffraction techniques and have confirmed a link between this and the Young's modulus of the feather rachis. Our results indicate that along the rachis length, from calamus to tip, the keratin molecules become more aligned than at the calamus before returning to a state of higher mis-orientation towards the tip of the rachis. We have also confirmed the general trend of increasing Young's modulus with distance along the rachis. Furthermore, we report a distinct difference in the patterns of orientation of beta-keratin in the feathers of flying and flightless birds. The trend for increased modulus along the feathers of volant birds is absent in the flightless ostrich.

Effects of high pressure on functional properties of soy protein

Simultaneous measurement of the effects of low soy protein concentration, pH and high pressure treatment at room temperature on solubility, emulsifying properties and rheological properties (loss modulus, G '') of soy protein isolate (SPI) were evaluated. Central composite rotatable designs (2(3)) were employed over two pH ranges (2.66-4.34 and 5.16-6.84) with SPI concentration (0.32-3.68%) and pressure (198-702 MPa) as the other independent variables. The surface responses were obtained for protein solubility, emulsifying activity index (EAI) and G ''. The samples with the highest effect on protein solubility, EAI and G '' values were evaluated, as well, by electrophoresis and free sulphydryl determination. The pH was the main factor that affected protein solubility, with solubility at a maximum at pH < 3 or pH > 6. Increasing SPI concentration and decreasing/increasing the pH away from the isoelectric point both caused a reduction in EAI. Loss modulus (G '') was found to increase with SPI concentration in both pH ranges. (c) 2006 Elsevier Ltd. All rights reserved.

Bode's maximum available feedback and phase margin

'Maximum Available Feedback' is Bode's term for the highest possible loop gain over a given bandwidth, with specified stability margins, in a single loop feedback system. His work using asymptotic analysis allowed Bode to develop a methodology for achieving this. However, the actual system performance differs from that specified, due to the use of asymptotic approximations, and the author[2] has described how, for instance, the actual phase margin is often much lower than required when the bandwidth is high, and proposed novel modifications to the asymptotes to address the issue. This paper gives some new analysis of such systems, showing that the method also contravenes Bode's definition of phase margin, and shows how the author's modifications can be used for different amounts of bandwidth.