The effect of preferential sampling on sampling variance

We examine some variations of standard probability designs that preferentially sample sites based on how easy they are to access. Preferential sampling designs deliver unbiased estimates of mean and sampling variance and will ease the burden of data collection but at what cost to our design efficiency? Preferential sampling has the potential to either increase or decrease sampling variance depending on the application. We carry out a simulation study to gauge what effect it will have when sampling Soil Organic Carbon (SOC) values in a large agricultural region in south-eastern Australia. Preferential sampling in this region can reduce the distance to travel by up to 16%. Our study is based on a dataset of predicted SOC values produced from a datamining exercise. We consider three designs and two ways to determine ease of access. The overall conclusion is that sampling performance deteriorates as the strength of preferential sampling increases, due to the fact the regions of high SOC are harder to access. So our designs are inadvertently targeting regions of low SOC value. The good news, however, is that Generalised Random Tessellation Stratification (GRTS) sampling designs are not as badly affected as others and GRTS remains an efficient design compared to competitors.

Resilience of a high-conservation-value, semi-arid grassland on fertile clay soils to burning, mowing and ploughing.

In grassland reserves, managed disturbance is often necessary to maintain plant species diversity. We carried out experiments to determine the impact of fire, kangaroo grazing, mowing and disc ploughing on grassland species richness and composition in a nature reserve in semi-arid eastern Australia. Vegetation response was influenced by winter-spring drought after establishment of the experiments, but moderate rainfall followed in late summer-autumn. Species composition varied greatly between sampling times, and the variability due to rainfall differences between seasons and years was greater than the effects of fire, kangaroo grazing, mowing or disc ploughing. In the fire experiment, species richness and composition recovered more rapidly after spring than autumn burning. Species richness and composition were similar to control sites within 12 months of burning and mowing, suggesting that removal of the dominant grass canopy is unnecessary to enhance plant diversity. Two fires (separated by 3 years) and post-fire kangaroo grazing had only minor influence on species richness and composition. Even disc ploughing caused only a small reduction in native richness. The minor impact of ploughing was explained by the small areas that were ploughed, the once-off nature of the treatment, and the high degree of natural movement and cracking in these shrink-swell soils. Recovery of the composition and richness of these grasslands was rapid because of the high proportion of perennial species that resprout vegetatively after fire and mowing. There appears to be little conservation benefit from fire, mowing or ploughing ungrazed areas, as we could identify no native plant species dependent on frequent disturbance for persistence in this grassland community. However, the ability of the Astrebla- and Dichanthium-dominated grasslands to recover quickly after disturbance, given favourable seasonal conditions, suggests that they are well adapted to natural disturbances (e.g. droughts, fire, flooding and native grazing).

Defining critical soil nutrient concentrations in soils supporting grains and cotton in Northern NSW and Queensland

This project focussed on the phosphorus (P) and potassium (K) status of northern cropping soils. Stores of P and K have been depleted by crop removal and limited fertiliser application, with depletion most significant in the subsoil. Soil testing strategies are confounded by slowly available mineral reserves with uncertain availability. The utility of new soil tests was assessed to measure these reserves, their availability to plants quantified and a regional sampling strategy undertaken to identify areas of greatest P and K deficit. Fertiliser application strategies for P and K have been tested and the interactions between these and other nutrients have been determined in a large field program.

Phosphorus in agricultural soils of Finland - characterization of reserves and retention in mineral soil profiles

An overwhelming majority of all the research on soil phosphorus (P) has been carried out with soil samples taken from the surface soils only, and our understanding of the forms and the reactions of P at a soil profile scale is based on few observations. In Finland, the interest in studying the P in complete soil profiles has been particularly small because of the lack of tradition in studying soil genesis, morphology, or classification. In this thesis, the P reserves and the retention of orthophosphate phosphorus (PO4-P) were examined in four cultivated mineral soil profiles in Finland (three Inceptisols and one Spodosol). The soils were classified according to the U.S. Soil Taxonomy and soil samples were taken from the genetic horizons in the profiles. The samples were analyzed for total P concentration, Chang and Jackson P fractions, P sorption properties, concentrations of water-extractable P, and for concentrations of oxalate-extractable Al and Fe. Theoretical P sorption capacities and degrees of P saturation were calculated with the data from the oxalate-extractions and the P fractionations. The studied profiles can be divided into sections with clearly differing P characteristics by their master horizons Ap, B and C. The C (or transitional BC) horizons below an approximate depth of 70 cm were dominated by, assumingly apatitic, H2SO4-soluble P. The concentration of total P in the C horizons ranged from 729 to 810 mg kg-1. In the B horizons between the depths of 30 and 70 cm, a significant part of the primary acid-soluble P has been weathered and transformed to secondary P forms. A mean weathering rate of the primary P in the soils was estimated to vary between 230 and 290 g ha-1 year-1. The degrees of P saturation in the B and C horizons were smaller than 7%, and the solubility of PO4-P was negligible. The P conditions in the Ap horizons differed drastically from those in the subsurface horizons. The high concentrations of total P (689-1870 mg kg-1) in the Ap horizons are most likely attributable to long-term cultivation with positive P balances. A significant proportion of the P in the Ap horizons occurred in the NH4F- and NaOH-extractable forms and as organic P. These three P pools, together with the concentrations of oxalate-extractable Al and Fe, seem to control the dynamics of PO4-P in the soils. The degrees of P saturation in the Ap horizons were greater (8-36%) than in the subsurface horizons. This was also reflected in the sorption experiments: Only the Ap horizons were able to maintain elevated PO4-P concentrations in the solution phase − all the subsoil horizons acted as sinks for PO4-P. Most of the available sorption capacity in the soils is located in the B horizons. The results suggest that this capacity could be utilized in reducing the losses of soluble P from excessively fertilized soils by mixing highly sorptive material from the B horizons with the P-enriched surface soil. The drastic differences in the P characteristics observed between adjoining horizons have to be taken into consideration when conducting soil sampling. Sampling of subsoils has to be made according to the genetic horizons or at small depth increments. Otherwise, contrasting materials are likely to be mixed in the same sample; and the results of such samples are not representative of any material present in the studied profile. Air-drying of soil samples was found to alter the results of the sorption experiments and the water extractions. This indicates that the studies on the most labile P forms in soil should be carried out with moist samples.

Soybean rotation and crop residue management to reduce nitrous oxide emissions from sugarcane soils

NITROUS OXIDE (N2O) IS a potent greenhouse gas and the predominant ozone-depleting substance in the atmosphere. Agricultural nitrogenous fertiliser use is the major source of human-induced N2O emissions. A field experiment was conducted at Bundaberg from October 2012 to September 2014 to examine the impacts of legume crop (soybean) rotation as an alternative nitrogen (N) source on N2O emissions during the fallow period and to investigate low-emission soybean residue management practices. An automatic monitoring system and manual gas sampling chambers were used to measure greenhouse gas emissions from soil. Soybean cropping during the fallow period reduced N2O emissions compared to the bare fallow. Based on the N content in the soybean crop residues, the fertiliser N application rate was reduced by about 120 kg N/ha for the subsequent sugarcane crop. Consequently, emissions of N2O during the sugarcane cropping season were significantly lower from the soybean cropped soil than those from the conventionally fertilised (145 kg N/ha) soil following bare fallow. However, tillage that incorporated the soybean crop residues into soil promoted N2O emissions in the first two months. Spraying a nitrification inhibitor (DMPP) onto the soybean crop residues before tillage effectively prevented the N2O emission spikes. Compared to conventional tillage, practising no-till with or without growing a nitrogen catch crop during the time after soybean harvest and before cane planting also reduced N2O emissions substantially. These results demonstrated that soybean rotation during the fallow period followed with N conservation management practices could offer a promising N2O mitigation strategy in sugarcane farming. Further investigation is required to provide guidance on N and water management following soybean fallow to maintain sugar productivity.

Sampling frequency affects estimates of annual nitrous oxide fluxes

Quantifying nitrous oxide (N(2)O) fluxes, a potent greenhouse gas, from soils is necessary to improve our knowledge of terrestrial N(2)O losses. Developing universal sampling frequencies for calculating annual N(2)O fluxes is difficult, as fluxes are renowned for their high temporal variability. We demonstrate daily sampling was largely required to achieve annual N(2)O fluxes within 10% of the best estimate for 28 annual datasets collected from three continents, Australia, Europe and Asia. Decreasing the regularity of measurements either under- or overestimated annual N(2)O fluxes, with a maximum overestimation of 935%. Measurement frequency was lowered using a sampling strategy based on environmental factors known to affect temporal variability, but still required sampling more than once a week. Consequently, uncertainty in current global terrestrial N(2)O budgets associated with the upscaling of field-based datasets can be decreased significantly using adequate sampling frequencies.

Country-scale carbon accounting of the vegetation and mineral soils of Finland

The increase in global temperature has been attributed to increased atmospheric concentrations of greenhouse gases (GHG), mainly that of CO2. The threat of severe and complex socio-economic and ecological implications of climate change have initiated an international process that aims to reduce emissions, to increase C sinks, and to protect existing C reservoirs. The famous Kyoto protocol is an offspring of this process. The Kyoto protocol and its accords state that signatory countries need to monitor their forest C pools, and to follow the guidelines set by the IPCC in the preparation, reporting and quality assessment of the C pool change estimates. The aims of this thesis were i) to estimate the changes in carbon stocks vegetation and soil in the forests in Finnish forests from 1922 to 2004, ii) to evaluate the applied methodology by using empirical data, iii) to assess the reliability of the estimates by means of uncertainty analysis, iv) to assess the effect of forest C sinks on the reliability of the entire national GHG inventory, and finally, v) to present an application of model-based stratification to a large-scale sampling design of soil C stock changes. The applied methodology builds on the forest inventory measured data (or modelled stand data), and uses statistical modelling to predict biomasses and litter productions, as well as a dynamic soil C model to predict the decomposition of litter. The mean vegetation C sink of Finnish forests from 1922 to 2004 was 3.3 Tg C a-1, and in soil was 0.7 Tg C a-1. Soil is slowly accumulating C as a consequence of increased growing stock and unsaturated soil C stocks in relation to current detritus input to soil that is higher than in the beginning of the period. Annual estimates of vegetation and soil C stock changes fluctuated considerably during the period, were frequently opposite (e.g. vegetation was a sink but soil was a source). The inclusion of vegetation sinks into the national GHG inventory of 2003 increased its uncertainty from between -4% and 9% to ± 19% (95% CI), and further inclusion of upland mineral soils increased it to ± 24%. The uncertainties of annual sinks can be reduced most efficiently by concentrating on the quality of the model input data. Despite the decreased precision of the national GHG inventory, the inclusion of uncertain sinks improves its accuracy due to the larger sectoral coverage of the inventory. If the national soil sink estimates were prepared by repeated soil sampling of model-stratified sample plots, the uncertainties would be accounted for in the stratum formation and sample allocation. Otherwise, the increases of sampling efficiency by stratification remain smaller. The highly variable and frequently opposite annual changes in ecosystem C pools imply the importance of full ecosystem C accounting. If forest C sink estimates will be used in practice average sink estimates seem a more reasonable basis than the annual estimates. This is due to the fact that annual forest sinks vary considerably and annual estimates are uncertain, and they have severe consequences for the reliability of the total national GHG balance. The estimation of average sinks should still be based on annual or even more frequent data due to the non-linear decomposition process that is influenced by the annual climate. The methodology used in this study to predict forest C sinks can be transferred to other countries with some modifications. The ultimate verification of sink estimates should be based on comparison to empirical data, in which case the model-based stratification presented in this study can serve to improve the efficiency of the sampling design.

Microbiological characterisation of soils : Evaluation of some critical steps in data collection and experimental design

The study of soil microbiota and their activities is central to the understanding of many ecosystem processes such as decomposition and nutrient cycling. The collection of microbiological data from soils generally involves several sequential steps of sampling, pretreatment and laboratory measurements. The reliability of results is dependent on reliable methods in every step. The aim of this thesis was to critically evaluate some central methods and procedures used in soil microbiological studies in order to increase our understanding of the factors that affect the measurement results and to provide guidance and new approaches for the design of experiments. The thesis focuses on four major themes: 1) soil microbiological heterogeneity and sampling, 2) storage of soil samples, 3) DNA extraction from soil, and 4) quantification of specific microbial groups by the most-probable-number (MPN) procedure. Soil heterogeneity and sampling are discussed as a single theme because knowledge on spatial (horizontal and vertical) and temporal variation is crucial when designing sampling procedures. Comparison of adjacent forest, meadow and cropped field plots showed that land use has a strong impact on the degree of horizontal variation of soil enzyme activities and bacterial community structure. However, regardless of the land use, the variation of microbiological characteristics appeared not to have predictable spatial structure at 0.5-10 m. Temporal and soil depth-related patterns were studied in relation to plant growth in cropped soil. The results showed that most enzyme activities and microbial biomass have a clear decreasing trend in the top 40 cm soil profile and a temporal pattern during the growing season. A new procedure for sampling of soil microbiological characteristics based on stratified sampling and pre-characterisation of samples was developed. A practical example demonstrated the potential of the new procedure to reduce the analysis efforts involved in laborious microbiological measurements without loss of precision. The investigation of storage of soil samples revealed that freezing (-20 °C) of small sample aliquots retains the activity of hydrolytic enzymes and the structure of the bacterial community in different soil matrices relatively well whereas air-drying cannot be recommended as a storage method for soil microbiological properties due to large reductions in activity. Freezing below -70 °C was the preferred method of storage for samples with high organic matter content. Comparison of different direct DNA extraction methods showed that the cell lysis treatment has a strong impact on the molecular size of DNA obtained and on the bacterial community structure detected. An improved MPN method for the enumeration of soil naphthalene degraders was introduced as an alternative to more complex MPN protocols or the DNA-based quantification approach. The main advantage of the new method is the simple protocol and the possibility to analyse a large number of samples and replicates simultaneously.

Analysis of Compressibility of Sensitive Soils

Sensitive soils, in general, are prone to mechanical disturbances while sampling, handling, and testing. This necessitates the prediction of true field behavior. The compressibility response of such soils is typical of having three zones, mechanistically explained as nonparticulate, transitional, and particulate. Such zoning has enabled the development of a simple method to predict the field compressibility response of the sample. The field compression curve with sigmact act as the most probable yield stress is considered to reflect 0% disturbance. By a comparison of experimentally determined sigmac and sigmact, it is possible to estimate the degree of sample disturbance. When the value of sigmac is closer to sigmact, the sampling disturbance approaches zero. As the value of sigmac reduces, the degree of sampling disturbance increases. The possibility of using this degree of sample disturbance from compressibility data to obtain other true properties from laboratory results of the sampled specimens has been examined.

Radioactivity Concentrations in Soils in the Qingdao Area, China

The specific activity concentrations of radionuclides U-238, Th-212, and K-40 of 2300 sampling points in the Qingdao area were measured by an FD-3022 gamma-ray spectrometer. The radioactivity concentrations of U-238, Th-232, and K-40 ranged from 3.3 to 185.3, from 6.9 to 157.2, and from 15.8 to 7834.4 Bq kg(-1), respectively. The air-absorbed dose at I meter above ground, effective annual dose, external hazard index, and radium equivalent activity were also calculated to systematically evaluate the radiological hazards of the natural radioactivity in Qingdao. The air-absorbed dose, effective annual dose, external hazard index, and radium equivalent activity in the study area were 98.6 nGy h(-1), 0.12 mSv, 0.56, 197 Bq kg(-1), respectively. Compared with the worldwide value, the air-absorbed dose is slightly high, but the other factors are all lower than the recommended value. The natural external exposure will not pose significant radiological threat to the population. In conclusion, the Qingdao area is safe with regard to the radiological level and suitable for living.

Organic Matter-Solid Phase Interactions Are Critical for Predicting Arsenic Release and Plant Uptake in Bangladesh Paddy Soils

Agroecological zones within Bangladesh with low levels of arsenic in groundwater and soils produce rice that is high in arsenic with respect to other producing regions of the globe. Little is known about arsenic cycling in these soils and the labile fractions relevant for plant uptake when flooded. Soil porewater dynamics of field soils (n = 39) were recreated under standardized laboratory conditions to investigate the mobility and interplay of arsenic, Fe, Si, C, and other elements, in relation to rice grain element composition, using the dynamic sampling technique diffusive gradients in thin films (DGT). Based on a simple model using only labile DGT measured arsenic and dissolved organic carbon (DOC), concentrations of arsenic in Aman (Monsoon season) rice grain were predicted reliably. DOC was the strongest determinant of arsenic solid-solution phase partitioning, while arsenic release to the soil porewater was shown to be decoupled from that of Fe. This study demonstrates the dual importance of organic matter (OM), in terms of enhancing arsenic release from soils, while reducing bioavailability by sequestering arsenic in solution.

Arsenic levels in the soils and macrophytes of the 'Entremuros' after the Aznalcóllar mine spill

In April 1998, a holding lagoon containing pyrite ore processing waste rich in arsenic, failed and released 5-6 million m(3) of highly polluting sludge and acidic water. Over 2700 ha of the internationally important Doñana National and Natural Parks were contaminated. The area of Natural Park to sustain the greatest impact was known as the Entremuros. This paper presents 0-5 cm soil monitoring data from the Entremuros, from sampling campaigns 6 and 18 months after the disaster; as well as macrophyte root, rhizome and stem data from samples taken 18 months after the spill. Results show a clear, decreasing, north-south arsenic soil pollution trend, both 6 and 18 months after the spill, and suggest a small reduction in total soil arsenic levels occurred over time; although a significant increase in extractable arsenic is also noted. The two macrophytes (Typha dominguensis and Scirpus maritimus) studied herein are not accumulating arsenic in stem parts, however, accumulation of arsenic on iron plaque on the roots of these plants may be occurring. Further work is recommended in order to determine the ecotoxicological significance of this process in relation to the avian food-chains of Doñana, and elsewhere.

Modelling local-scale determinants and the probability of microarthropod species occurrence in Antarctic soils

Soil fauna in the extreme conditions of Antarctica consists of a few microinvertebrate species patchily distributed at different spatial scales. Populations of the prostigmatic mite Stereotydeus belli and the collembolan Gressittacantha terranova from northern Victoria Land (Antarctica) were used as models to study the effect of soil properties on microarthropod distributions. In agreement with the general assumption that the development and distribution of life in these ecosystems is mainly controlled by abiotic factors, we found that the probability of occurrence of S. belli depends on soil moisture and texture and on the sampling period (which affects the general availability of water); surprisingly, none of the analysed variables were significantly related to the G. terranova distribution. Based on our results and literature data, we propose a theoretical model that introduces biotic interactions among the major factors driving the local distribution of collembolans in Antarctic terrestrial ecosystems. (c) 2007 Elsevier Ltd. All rights reserved.

Compositonal Analysis of Potentially harmful Elements in Soils: Implications for Epidemiology and Kidney Disease

A substantial proportion of aetiological risks for many cancers and chronic diseases remain unexplained. Using geochemical soil and stream water samples collected as part of the Tellus Project studies, current research is investigating naturally occurring background levels of potentially toxic elements (PTEs) in soils and stream sediments and their possible relationship with progressive chronic kidney disease (CKD). The Tellus geological mapping project, Geological Survey Northern Ireland, collected soil sediment and stream water samples on a grid of one sample site every 2 km2 across the rural areas of Northern Ireland resulting in an excess of 6800 soil sampling locations and more than 5800 locations for stream water sampling. Accumulation of several PTEs including arsenic, cadmium, chromium, lead and mercury have been linked with human health and implicated in renal function decline. The hypothesis is that long-term exposure will result in cumulative exposure to PTEs and act as risk factor(s) for cancer and diabetes related CKD and its progression. The ‘bioavailable’ fraction of total PTE soil concentration depends on the ‘bioaccessible’ proportion through an exposure pathway. Recent work has explored this bioaccessible fraction for a range of PTEs across Northern Ireland. In this study the compositional nature of the multivariate geochemical PTE variables and bioaccessible data is explored to augment the investigation into the potential relationship between PTEs, bioaccessibility and disease data.

Magnitude, form and bioavailability of fluvial carbon exports from Irish organic soils under pasture

Organic soils are widespread in Ireland and vulnerable to degradation via drainage for agriculture. The soil-landuse combination of pasture on organic soils may play a disproportionate role in regional C dynamics but is yet to receive study. Fluvial C fluxes and labile organic fractions were determined for two such sites at nested field (c.4 ha) and subcatchment scales (>40 ha); one relatively dry and nutrient rich, the other wetter and nutrient poor. Field scale flux from the nutrient poor site over 2 years was 38.9 ± 6.6 g C m⁻² yr⁻¹ with DIC > DOC > POC at 57, 32 and 11 % respectively, and 72 % DIC was comprised of above equilibrium CO₂. At the nutrient rich site, which overlies limestone geology, field scale export over an individual year was 90.4 g C m⁻² with DIC > DOC > POC at 49, 42 and 9 %, but with 90 % DIC as bicarbonate. By comparison with the nutrient poor site, the magnitude and composition of inorganic C exports from the nutrient rich site implied considerable export of soil-respiratory C as bicarbonate, and lower evasion losses due to carbonate system buffering. Labile DOC determined using dark incubations indicated small fractions (5–10 %) available for remineralisation over typical downstream transit times of days to weeks. These fractions are probably conservative as photolysis in the environment can increase the proportion of labile compounds via photocleavage and directly remineralise organic matter. This study demonstrates that monitoring at soil–water interfaces can aid capture of total landscape fluvial fluxes by precluding the need to incorporate prior C evasion, although rapid runoff responses at field scales can necessitate high resolution flow proportional, and hydrograph sampling to constrain uncertainty of flux estimates.