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.

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.

Recent glacier retreat in the Caucasus Mountains, Russia, and associated increase in supraglacial debris cover and supra-/proglacial lake development

This paper reports changes in supraglacial debris cover and supra-/proglacial lake development associated with recent glacier retreat (1985-2000) in the central Caucasus Mountains, Russia. Satellite imagery (Landsat TM and ETM+) was used to map the surface area and supraglacial debris cover on six neighbouring glaciers in the Adylsu valley through a process of manual digitizing on a false-colour composite of bands 5, 4, 3 (red, green, blue). The distribution and surface area of supraglacial and proglacial lakes was digitized for a larger area, which extended to the whole Landsat scene. We also compare our satellite interpretations to field observations in the Adylsu valley. Supraglacial debris cover ranges from < 5% to > 25% on individual glaciers, but glacier retreat between 1985 and 2000 resulted in a 3-6% increase in the proportion of each glacier covered by debris. The only exception to this trend was a very small glacier where debris cover did not change significantly and remote mapping proved more difficult. The increase in debris cover is characterized by a progressive upglacier migration, which we suggest is being driven by focused ablation (and therefore glacier thinning) at the up-glacier limit of the debris cover, resulting in the progressive exposure of englacial debris. Glacier retreat has also been accompanied by an increase in the number of proglacial and supraglacial lakes in our study area, from 16 in 1985 to 24 in 2000, representing a 57% increase in their cumulative surface area. These lakes appear to be impounded by relatively recently lateral and terminal moraines and by debris deposits on the surface of the glacier. The changes in glacier surface characteristics reported here are likely to exert a profound influence on glacier mass balance and their future response to climate change. They may also increase the likelihood of glacier-related hazards (lake outbursts, debris slides), and future monitoring is recommended.

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.

Why the western Pacific subtropical high has extended westward since the late 1970s

The western Pacific subtropical high (WPSH) is closely related to Asian climate. Previous examination of changes in the WPSH found a westward extension since the late 1970s, which has contributed to the inter-decadal transition of East Asian climate. The reason for the westward extension is unknown, however. The present study suggests that this significant change of WPSH is partly due to the atmosphere's response to the observed Indian Ocean-western Pacific (IWP) warming. Coordinated by a European Union's Sixth Framework Programme, Understanding the Dynamics of the Coupled Climate System (DYNAMITE), five AGCMs were forced by identical idealized sea surface temperature patterns representative of the IWP warming and cooling. The results of these numerical experiments suggest that the negative heating in the central and eastern tropical Pacific and increased convective heating in the equatorial Indian Ocean/ Maritime Continent associated with IWP warming are in favor of the westward extension of WPSH. The SST changes in IWP influences the Walker circulation, with a subsequent reduction of convections in the tropical central and eastern Pacific, which then forces an ENSO/Gill-type response that modulates the WPSH. The monsoon diabatic heating mechanism proposed by Rodwell and Hoskins plays a secondary reinforcing role in the westward extension of WPSH. The low-level equatorial flank of WPSH is interpreted as a Kelvin response to monsoon condensational heating, while the intensified poleward flow along the western flank of WPSH is in accord with Sverdrup vorticity balance. The IWP warming has led to an expansion of the South Asian high in the upper troposphere, as seen in the reanalysis.

Dynamic scaling of bred vectors in spatially extended chaotic systems

We unfold a profound relationship between the dynamics of finite-size perturbations in spatially extended chaotic systems and the universality class of Kardar-Parisi-Zhang (KPZ). We show how this relationship can be exploited to obtain a complete theoretical description of the bred vectors dynamics. The existence of characteristic length/time scales, the spatial extent of spatial correlations and how to time it, and the role of the breeding amplitude are all analyzed in the light of our theory. Implications to weather forecasting based on ensembles of initial conditions are also discussed.

Galactic cosmic ray hazard in the unusual extended solar minimum between solar cycle 23 and 24

Galactic cosmic rays (GCRs) are extremely difficult to shield against and pose one of the most severe long-term hazards for human exploration of space. The recent solar minimum between solar cycles 23 and 24 shows a prolonged period of reduced solar activity and low interplanetary magnetic field strengths. As a result, the modulation of GCRs is very weak, and the fluxes of GCRs are near their highest levels in the last 25 years in the fall of 2009. Here we explore the dose rates of GCRs in the current prolonged solar minimum and make predictions for the Lunar Reconnaissance Orbiter (LRO) Cosmic Ray Telescope for the Effects of Radiation (CRaTER), which is now measuring GCRs in the lunar environment. Our results confirm the weak modulation of GCRs leading to the largest dose rates seen in the last 25 years over a prolonged period of little solar activity.

Discriminative pattern mining in software fault detection

We present a method to enhance fault localization for software systems based on a frequent pattern mining algorithm. Our method is based on a large set of test cases for a given set of programs in which faults can be detected. The test executions are recorded as function call trees. Based on test oracles the tests can be classified into successful and failing tests. A frequent pattern mining algorithm is used to identify frequent subtrees in successful and failing test executions. This information is used to rank functions according to their likelihood of containing a fault. The ranking suggests an order in which to examine the functions during fault analysis. We validate our approach experimentally using a subset of Siemens benchmark programs.

Sonogashira coupling on an extended gold surface in vacuo: reaction of phenylacetylene with iodobenzene on Au(111)

Temperature-programmed reaction measurements supported by scanning tunneling microscopy have shown that phenylacetylene and iodobenzene react on smooth Au(111) under vacuum conditions to yield biphenyl and diphenyldiacetylene, the result of homocoupling of the reactant molecules. They also produce diphenylacetylene, the result of Sonogashira cross-coupling, prototypical of a class of reactions that are of paramount importance in synthetic organic chemistry and whose mechanism remains controversial. Roughened Au(111) is completely inert toward all three reactions, indicating that the availability of crystallographically well-defined adsorption sites is crucially important. High-resolution X-ray photoelectron spectroscopy and near-edge X-ray absorption fine structure spectroscopy show that the reactants are initially present as intact, essentially flat-lying molecules and that the temperature threshold for Sonogashira coupling coincides with that for C−I bond scission in the iodobenzene reactant. The fractional-order kinetics and low temperature associated with desorption of the Sonogashira product suggest that the reaction occurs at the boundaries of islands of adsorbed reactants and that its appearance in the gas phase is rate-limited by the surface reaction. These findings demonstrate unambiguously and for the first time that this heterogeneous cross-coupling chemistry is an intrinsic property of extended, metallic pure gold surfaces: no other species, including solvent molecules, basic or charged (ionic) species are necessary to mediate the process.

Mining, grievance and conflict in rural Tanzania

This paper provides an extended analysis of the tensions that have surfaced between large-scale mine operators and artisanal miners in gold-rich areas of rural Tanzania. The literature on grievance is used to contextualise, these disputes, the underlying cause of which is artisanal miners' mounting frustration over not being able to secure viable concessions to work. Newly implemented legislation has, for the most part, empowered foreign large-scale mine operators, while simultaneously disempowering indigenous small-scale miners. In many cases, the former have addressed mounting security and community problems on their own. Until the country's major mine operators extend assistance to marginalised small-scale mining groups, the likelihood of violent conflict unfolding between these parties will increase.

The use of early lactation milk protein content to predict subsequent fertility performance and likelihood of culling, in commercial dairy cows

A dataset of 1,846,990 completed lactation record,; was created Using milk recording data from 8,967 commercial dairy farms in the United Kingdom over a five year period. Herd-specific lactation curves describing levels of milk, Cat and protein by lactation number and month of calving were generated for each farm. The actual yield of milk and protein proportion at the first milk recording of individual cow lactations were compared with the levels taken from the lactation curves. Logistic regression analysis showed that cows production milk with a lower percentage of protein than average had a significantly lower probability of being in-calf at 100 days post calving and it significantly higher probability of being culled at the end of lactation. The culling rates derived from the studied database demonstrate the current high wastage rate of commercial dairy cows. Well of this wastage is due to involuntary culling as a result of reproductive failure.

Meta-analysis: a unifying meta-likelihood approach framing unobserved heterogeneity, study covariates, publication bias, and study quality

OBJECTIVES: This contribution provides a unifying concept for meta-analysis integrating the handling of unobserved heterogeneity, study covariates, publication bias and study quality. It is important to consider these issues simultaneously to avoid the occurrence of artifacts, and a method for doing so is suggested here. METHODS: The approach is based upon the meta-likelihood in combination with a general linear nonparametric mixed model, which lays the ground for all inferential conclusions suggested here. RESULTS: The concept is illustrated at hand of a meta-analysis investigating the relationship of hormone replacement therapy and breast cancer. The phenomenon of interest has been investigated in many studies for a considerable time and different results were reported. In 1992 a meta-analysis by Sillero-Arenas et al. concluded a small, but significant overall effect of 1.06 on the relative risk scale. Using the meta-likelihood approach it is demonstrated here that this meta-analysis is due to considerable unobserved heterogeneity. Furthermore, it is shown that new methods are available to model this heterogeneity successfully. It is argued further to include available study covariates to explain this heterogeneity in the meta-analysis at hand. CONCLUSIONS: The topic of HRT and breast cancer has again very recently become an issue of public debate, when results of a large trial investigating the health effects of hormone replacement therapy were published indicating an increased risk for breast cancer (risk ratio of 1.26). Using an adequate regression model in the previously published meta-analysis an adjusted estimate of effect of 1.14 can be given which is considerably higher than the one published in the meta-analysis of Sillero-Arenas et al. In summary, it is hoped that the method suggested here contributes further to a good meta-analytic practice in public health and clinical disciplines.