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 release of silicon from amorphous silica and rice straw in Sri Lankan soils

Silicon release from rice straw and amorphous silica when shaken in solution with five Sri Lankan soils was studied indirectly using sorption isotherms and changes in concentration and directly using straw in dialysis bags examined using electron microscopy. The aim was to further our understanding of the processes and factors affecting the release of straw-Si in soils and its availability to rice. The soils (alfisols and ultisols) shaken with 0.1 M NaCl (5 g per 125 mL for 250 days) produced concentrations of 1 - 4 mg L-1 of monosilicic acid-Si. Amorphous silica added to these suspensions (36.5 mg, containing 17 mg Si) raised the concentrations to 20 - 40 mg L-1, and added rice straw (0.5 g, containing 17 mg Si) gave 10 - 25 mg L-1. Sorption isotherms (7 days equilibrations) were used to calculate from the concentrations the amounts of Si released ( 24 - 38% and 8 - 21%, respectively). Both materials gave about 40 mg L-1 of monosilicic acid-Si plus 30 mg L-1 of disilicic acid-Si when shaken in solution alone (5 g per 125 mL). Straw in dialysis bags ( 0.5 g per 25 mL in 0.1 M NaCl) was shaken in soil suspension ( 5 g per 100 mL) for 60 days. Similar concentrations and releases were measured to those obtained above. About one fifth of the mass of straw was lost by decomposition in the first 15 days. A chloroform treatment prevented decomposition, but Si release was unaffected. Disintegration continued throughout the experiments, with phytoliths being exposed and dissolved. Compared to the rate of release from straw into solution without soil, the release of Si into soil suspensions was increased during the first 20 days by adsorption on the soil, but was then reduced probably through the effect of Fe and Al on the phytolith surfaces. The extent of this blocking effect varied between soils and was not simply related to soil pH.

Quantifying time scales of pedogenic calcrete formation using U-series disequilibria

Time scales of pedogenic calcrete development are quantified by subsampling carbonate from within a mature (stage V) pedogenic calcrete profile from southeast Spain and dating the material by U-series disequilibria. The location of the earliest and latest cements can be estimated by comparing previous studies of calcrete morphological development with micromorphological analysis of the study profile. Carbonate was sampled and dated from three locations within the profile: (1) below the lower surface of clasts within the hardpan (representing the earliest cement present-207 +/- 11 ka), (2) from the centre of cement filled pores within the hardpan (reflecting the final plugging of the calcrete hardpan-155 +/- 9 ka) and (3) from the laminar calcrete overlying the hardpan (representing the latest cement-112 +/- 15 ka). These results show that the hardpan took between 73 and 31 ka to form, whilst the mature stage V profile took between 121 and 69 ka to form. This is the first time that rates of mature calcrete development have been established by direct radiometric dating of the authigenic carbonate. The technique is appropriate for dating mature calcretes in dryland regions worldwide and offers the opportunity of increasing our understanding of the spatial and temporal variability in rates of pedogenic calcrete development. (C) 2004 Elsevier B.V. All rights reserved.

Scales and voices in World Historical Archaeology

This paper briefly sets the scene for the articles that follow, introducing some key debates that have characterized the recent practice of historical archaeology. The definition of historical archaeology is explored according to parameters of chronology and methodology, drawing a distinction between New World traditions that define the subject as 'post-Columbian' and Old World approaches that establish broader connections with the 'documentary archaeology' of all literate societies. Current issues in European and American historical archaeology are highlighted, including the gradual breakdown of the medieval/post-medieval divide and the call for a global 'modern-world archaeology' to address the 'grand historical narratives' of the period, such as capitalism, economic improvement, and consumerism. The resistance to this global research agenda is explored with reference to archaeologies of diaspora and postcolonialism, which demand local perspectives to explore diversity and meaning. Finally, the innovative use of community archaeology and multi-vocality is introduced, with particular reference to the experimental narratives pursued by American historical archaeologists, in their new role as 'storytellers'.

Introduction: scales and voices in world historical archaeology

This paper briefly sets the scene for the articles that follow, introducing some key debates that have characterized the recent practice of historical archaeology. The definition of historical archaeology is explored according to parameters of chronology and methodology, drawing a distinction between New World traditions that define the subject as 'post-Columbian' and Old World approaches that establish broader connections with the 'documentary archaeology' of all literate societies. Current issues in European and American historical archaeology are highlighted, including the gradual breakdown of the medieval/post-medieval divide and the call for a global 'modern-world archaeology' to address the 'grand historical narratives' of the period, such as capitalism, economic improvement, and consumerism. The resistance to this global research agenda is explored with reference to archaeologies of diaspora and postcolonialism, which demand local perspectives to explore diversity and meaning. Finally, the innovative use of community archaeology and multi-vocality is introduced, with particular reference to the experimental narratives pursued by American historical archaeologists, in their new role as 'storytellers'.

Characteristics of silica-rich pedofeatures in a buried paleosol

Whitish and whitish-light brown milky-like textural pedofeatures and impregnations were found in the voids and the matrix of buried paleosols older than 2.7 million years in a site in Sardinia, Italy. The pedofeatures were described and analysed using micromorphology, X-ray diffraction and microprobe techniques, and their spatial distribution correlated with field evidence. The suite of analyses showed that the main components of the pedofeatures were more or less ordered silica phases. As well as forming a matrix cement, these pedofeatures also occurred as coatings and infillings in pores. Significant amounts of alumina and, less, Mg, Ca and Fe were also present in the pedofeatures, possibly in the form of silicate coatings and inclusions/impurities, or alumino-silicates of the adjacent soil matrix. A number of hypotheses are drawn on the possible mechanisms of formation of these silica-rich pedofeatures, including the possibility of prolonged weathering of volcanic materials and the resulting formation of colloids and more or less ordered silica phases, with successive dehydration and progressive ordering of phases during the at least 2.5 million years. (C) 2003 Elsevier B.V. All rights reserved.

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.

Evaluating the risk of non-point source pollution from biosolids: integrated modelling of nutrient losses at field and catchment scales

A semi-distributed model, INCA, has been developed to determine the fate and distribution of nutrients in terrestrial and aquatic systems. The model simulates nitrogen and phosphorus processes in soils, groundwaters and river systems and can be applied in a semi-distributed manner at a range of scales. In this study, the model has been applied at field to sub-catchment to whole catchment scale to evaluate the behaviour of biosolid-derived losses of P in agricultural systems. It is shown that process-based models such as INCA, applied at a wide range of scales, reproduce field and catchment behaviour satisfactorily. The INCA model can also be used to generate generic information for risk assessment. By adjusting three key variables: biosolid application rates, the hydrological connectivity of the catchment and the initial P-status of the soils within the model, a matrix of P loss rates can be generated to evaluate the behaviour of the model and, hence, of the catchment system. The results, which indicate the sensitivity of the catchment to flow paths, to application rates and to initial soil conditions, have been incorporated into a Nutrient Export Risk Matrix (NERM).

The release of silicon from amorphous silica and rice straw in Sri Lankan soils

Silicon release from rice straw and amorphous silica when shaken in solution with five Sri Lankan soils was studied indirectly using sorption isotherms and changes in concentration and directly using straw in dialysis bags examined using electron microscopy. The aim was to further our understanding of the processes and factors affecting the release of straw-Si in soils and its availability to rice. The soils (alfisols and ultisols) shaken with 0.1 M NaCl (5 g per 125 mL for 250 days) produced concentrations of 1 - 4 mg L-1 of monosilicic acid-Si. Amorphous silica added to these suspensions (36.5 mg, containing 17 mg Si) raised the concentrations to 20 - 40 mg L-1, and added rice straw (0.5 g, containing 17 mg Si) gave 10 - 25 mg L-1. Sorption isotherms (7 days equilibrations) were used to calculate from the concentrations the amounts of Si released ( 24 - 38% and 8 - 21%, respectively). Both materials gave about 40 mg L-1 of monosilicic acid-Si plus 30 mg L-1 of disilicic acid-Si when shaken in solution alone (5 g per 125 mL). Straw in dialysis bags ( 0.5 g per 25 mL in 0.1 M NaCl) was shaken in soil suspension ( 5 g per 100 mL) for 60 days. Similar concentrations and releases were measured to those obtained above. About one fifth of the mass of straw was lost by decomposition in the first 15 days. A chloroform treatment prevented decomposition, but Si release was unaffected. Disintegration continued throughout the experiments, with phytoliths being exposed and dissolved. Compared to the rate of release from straw into solution without soil, the release of Si into soil suspensions was increased during the first 20 days by adsorption on the soil, but was then reduced probably through the effect of Fe and Al on the phytolith surfaces. The extent of this blocking effect varied between soils and was not simply related to soil pH.

Linear Contributions of Different Time Scales to Teleconnectivity

The contributions of different time scales to extratropical teleconnections are examined. By applying empirical orthogonal functions and correlation analyses to reanalysis data, it is shown that eddies with periods shorter than 10 days have no linear contribution to teleconnectivity. Instead, synoptic variability follows wavelike patterns along the storm tracks, interpreted as propagating baroclinic disturbances. In agreement with preceding studies, it is found that teleconnections such as the North Atlantic Oscillation (NAO) and the Pacific–North America (PNA) pattern occur only at low frequencies, typically for periods more than 20 days. Low-frequency potential vorticity variability is shown to follow patterns analogous to known teleconnections but with shapes that differ considerably from them. It is concluded that the role, if any, of synoptic eddies in determining and forcing teleconnections needs to be sought in nonlinear interactions with the slower transients. The present results demonstrate that daily variability of teleconnection indices cannot be interpreted in terms of the teleconnection patterns, only the slow part of the variability.

Impacts of a pesticide on pollinator species richness at different spatial scales

Pesticides are an important potential cause of biodiversity and pollinator decline. Little is known about the impacts of pesticides on wild pollinators in the field. Insect pollinators were sampled in an agricultural system in Italy with the aim of detecting the impacts of pesticide use. The insecticide fenitrothion was over 150 times greater in toxicity than other pesticides used in the area, so sampling was set up around its application. Species richness of wild bees, bumblebees and butterflies were sampled at three spatial scales to assess responses to pesticide application: (i) the ‘field’ scale along pesticide drift gradients; (ii) the ‘landscape’ scale sampling in different crops within the area and (iii) the ‘regional’ scale comparing two river basins with contrasting agricultural intensity. At the field scale, the interaction between the application regime of the insecticide and the point in the season was important for species richness. Wild bee species richness appeared to be unaffected by one insecticide application, but declined after two and three applications. At the landscape scale, the species richness of wild bees declined in vine fields where the insecticide was applied, but did not decline in maize or uncultivated fields. At the regional scale, lower bumblebee and butterfly species richness was found in the more intensively farmed basin with higher pesticide loads. Our results suggest that wild bees are an insect pollinator group at particular risk from pesticide use. Further investigation is needed on how the type, quantity and timing of pesticide application impacts pollinators.

Insights into processes and deposits of hazardous vulcanian explosions at Soufrière Hills Volcano during 2008 and 2009 (Montserrat, West Indies)

During the Soufrière Hills eruption, vulcanian explosions have generally occurred 1) in episodic cycles; 2) isolated during pauses in extrusion, and 3) after major collapses of the dome. In a different eruptive context, significant vulcanian explosions occurred on 29 July 2008, 3 December 2008, and 3 January 2009. Deposits are pumiceous except for the 3 December event. We reconstructed the dispersal pattern of the deposits and their textural characteristics to evaluate erupted volume and vesicularity of the magma at fragmentation. We discuss the implications of these explosions in terms of eruptive processes and chronology, and the hazards posed by their sudden and often unheralded occurrence. We suggest that overpressurization of the conduit can develop over time-scales of months to weeks by a process of self-sealing of conduit walls and/or the cooling dome by silica polymorphs. This work provides new insights for understanding the generation of hazardous vulcanian explosions at andesitic volcanoes.

Maize yield and rainfall on different spatial and temporal scales in Southern Brazil

This study aimed to establish relationships between maize yield and rainfall on different temporal and spatial scales, in order to provide a basis for crop monitoring and modelling. A 16-year series of maize yield and daily rainfall from 11 municipalities and micro-regions of Rio Grande do Sul State was used. Correlation and regression analyses were used to determine associations between crop yield and rainfall for the entire crop cycle, from tasseling to 30 days after, and from 5 days before tasseling to 40 days after. Close relationships between maize yield and rainfall were found, particularly during the reproductive period (45-day period comprising the flowering and grain filling). Relationships were closer on a regional scale than at smaller scales. Implications of the crop-rainfall relationships for crop modelling are discussed.