Applying HYDRUS to flow in a sodic clay soil with solution composition–dependent hydraulic conductivity

Management of sodic soils under irrigation often requires application of chemical ameliorants to improve permeability combined with leaching of excess salts. Modeling irrigation, soil treatments, and leaching in these sodic soils requires a model that can adequately represent the physical and chemical changes in the soil associated with the amelioration process. While there are a number of models that simulate reactive solute transport, UNSATCHEM and HYDRUS-1D are currently the only models that also include an ability to simulate the impacts of soil chemistry on hydraulic conductivity. Previous researchers have successfully applied these models to simulate amelioration experiments on a sodic loam soil. To further gauge their applicability, we extended the previous work by comparing HYDRUS simulations of sodic soil amelioration with the results from recently published laboratory experiments on a more reactive, repacked sodic clay soil. The general trends observed in the laboratory experiments were able to be simulated using HYDRUS. Differences between measured and simulated results were attributed to the limited flexibility of the function that represents chemistry-dependent hydraulic conductivity in HYDRUS. While improvements in the function could be made, the present work indicates that HYDRUS-UNSATCHEM captures the key changes in soil hydraulic properties that occur during sodic clay soil amelioration and thus extends the findings of previous researchers studying sodic loams.

Hydraulic conductivity increases in a sodic clay soil in response to gypsum applications impacts of bulk density and cation exchange

Amelioration of sodic soils is commonly achieved by applying gypsum, which increases soil hydraulic conductivity by altering soil chemistry. The magnitude of hydraulic conductivity increases expected in response to gypsum applications depends on soil properties including clay content, clay mineralogy, and bulk density. The soil analyzed in this study was a kaolinite rich sodic clay soil from an irrigated area of the Lower Burdekin coastal floodplain in tropical North Queensland, Australia. The impact of gypsum amelioration was investigated by continuously leaching soil columns with a saturated gypsum solution, until the hydraulic conductivity and leachate chemistry stabilized. Extended leaching enabled the full impacts of electrolyte effects and cation exchange to be determined. For the columns packed to 1.4 g/cm3, exchangeable sodium concentrations were reduced from 5.0 ± 0.5 mEq/100 g to 0.41 ± 0.06 mEq/100 g, exchangeable magnesium concentrations were reduced from 13.9 ± 0.3 mEq/100 g to 4.3 ± 2.12 mEq/100 g, and hydraulic conductivity increased to 0.15 ± 0.04 cm/d. For the columns packed to 1.3 g/cm3, exchangeable sodium concentrations were reduced from 5.0 ± 0.5 mEq/100 g to 0.51 ± 0.03 mEq/100 g, exchangeable magnesium concentrations were reduced from 13.9 ± 0.3 mEq/100 g to 0.55 ± 0.36 mEq/100 g, and hydraulic conductivity increased to 0.96 ± 0.53 cm/d. The results of this study highlight that both sodium and magnesium need to be taken into account when determining the suitability of water quality for irrigation of sodic soils and that soil bulk density plays a major role in controlling the extent of reclamation that can be achieved using gypsum applications.

Effect of nitrification inhibitors (DMPP and 3MP+ TZ) on soil nitrous oxide emissions from a sub-tropical vegetable production system in Queensland, Australia

The use of nitrification inhibitors, in combination with ammonium based fertilisers, has been promoted recently as an effective method to reduce nitrous oxide (N2O) emissions from fertilised agricultural fields, whilst increasing yield and nitrogen use efficiency. Vegetable cropping systems are often characterised by high inputs of nitrogen fertiliser and consequently elevated emissions of nitrous oxide (N2O) can be expected. However, to date only limited data is available on the use of nitrification inhibitors in sub-tropical vegetable systems. A field experiment investigated the effect of the nitrification inhibitors (DMPP & 3MP+TZ) on N2O emissions and yield from a typical vegetable production system in sub-tropical Australia. Soil N2O fluxes were monitored continuously over an entire year with a fully automated system. Measurements were taken from three subplots for each treatment within a randomized complete blocks design. There was a significant inhibition effect of DMPP and 3MP+TZ on N2O emissions and soil mineral N content directly following the application of the fertiliser over the vegetable cropping phase. However this mitigation was offset by elevated N2O emissions from the inhibitor treatments over the post-harvest fallow period. Cumulative annual N2O emissions amounted to 1.22 kg-N/ha, 1.16 kg-N/ha, 1.50 kg-N/ha and 0.86 kg-N/ha in the conventional fertiliser (CONV), the DMPP treatment, the 3MP+TZ treatment and the zero fertiliser (0N) respectively. Corresponding fertiliser induced emission factors (EFs) were low with only 0.09 - 0.20% of the total applied fertiliser lost as N2O. There was no significant effect of the nitrification inhibitors on yield compared to the CONV treatment for the three vegetable crops (green beans, broccoli, lettuce) grown over the experimental period. This study highlights that N2O emissions from such vegetable cropping system are primarily controlled by post-harvest emissions following the incorporation of vegetable crop residues into the soil. It also shows that the use of nitrification inhibitors can lead to elevated N2O emissions by storing N in the soil profile that is available to soil microbes during the decomposition of the vegetable residues over the post-harvest phase. Hence the use of nitrification inhibitors in vegetable systems has to be treated carefully and fertiliser rates need to be adjusted to avoid excess soil nitrogen during the postharvest phase.

A novel 15N tracer approach for the quantification of N2 and N2O emissions from soil incubations in a completely automated laboratory set up

The microbial mediated production of nitrous oxide (N2O) and its reduction to dinitrogen (N2) via denitrification represents a loss of nitrogen (N) from fertilised agro-ecosystems to the atmosphere. Although denitrification has received great interest by biogeochemists in the last decades, the magnitude of N2lossesand related N2:N2O ratios from soils still are largely unknown due to methodical constraints. We present a novel 15N tracer approach, based on a previous developed tracer method to study denitrification in pure bacterial cultures which was modified for the use on soil incubations in a completely automated laboratory set up. The method uses a background air in the incubation vessels that is replaced with a helium-oxygen gas mixture with a 50-fold reduced N2 background (2 % v/v). This method allows for a direct and sensitive quantification of the N2 and N2O emissions from the soil with isotope-ratio mass spectrometry after 15N labelling of denitrification N substrates and minimises the sensitivity to the intrusion of atmospheric N2 at the same time. The incubation set up was used to determine the influence of different soil moisture levels on N2 and N2O emissions from a sub-tropical pasture soil in Queensland/Australia. The soil was labelled with an equivalent of 50 μg-N per gram dry soil by broadcast application of KNO3solution (4 at.% 15N) and incubated for 3 days at 80% and 100% water filled pore space (WFPS), respectively. The headspace of the incubation vessel was sampled automatically over 12hrs each day and 3 samples (0, 6, and 12 hrs after incubation start) of headspace gas analysed for N2 and N2O with an isotope-ratio mass spectrometer (DELTA V Plus, Thermo Fisher Scientific, Bremen, Germany(. In addition, the soil was analysed for 15N NO3- and NH4+ using the 15N diffusion method, which enabled us to obtain a complete N balance. The method proved to be highly sensitive for N2 and N2O emissions detecting N2O emissions ranging from 20 to 627 μN kg-1soil-1hr-1and N2 emissions ranging from 4.2 to 43 μN kg-1soil-1hr-1for the different treatments. The main end-product of denitrification was N2O for both water contents with N2 accounting for 9% and 13% of the total denitrification losses at 80% and 100%WFPS, respectively. Between 95-100% of the added 15N fertiliser could be recovered. Gross nitrification over the 3 days amounted to 8.6 μN g-1 soil-1 and 4.7 μN g-1 soil-1, denitrification to 4.1 μN g-1 soil-1 and 11.8 μN g-1 soil-1at 80% and 100%WFPS, respectively. The results confirm that the tested method allows for a direct and highly sensitive detection of N2 and N2O fluxes from soils and hence offers a sensitive tool to study denitrification and N turnover in terrestrial agro-ecosystems.

Soil N2O emissions under N2-fixing legumes and N-fertilised canola: A reappraisal of emissions factor calculations

Introducing nitrogen (N)-fixing legumes into cereal-based crop rotations reduces synthetic fertiliser-N use and may mitigate soil emissions of nitrous oxide (N2O). Current IPCC calculations assume 100% of legume biomass N as the anthropogenic N input and use 1% of this as an emission factor (EF)—the percentage of input N emitted as N2O. However, legumes also utilise soil inorganic N, so legume-fixed N is typically less than 100% of legume biomass N. In two field experiments, we measured soil N2O emissions from a black Vertosol in sub-tropical Australia for 12 months after sowing of chickpea (Cicer arietinum L.), canola (Brassica napus L.), faba bean (Vicia faba L.), and field pea (Pisum sativum L.). Cumulative N2O emissions from N-fertilised canola (624 g N2O-N ha−1) greatly exceeded those from chickpea (127 g N2O-N ha−1) in Experiment 1. Similarly, N2O emitted from canola (385 g N2O-N ha−1) in Experiment 2 was significantly greater than chickpea (166 g N2O-N ha−1), faba bean (166 g N2O-N ha−1) or field pea (135 g N2O-N ha−1). Highest losses from canola were recorded during the growing season, whereas 75% of the annual N2O losses from the legumes occurred post-harvest. Legume N2-fixation provided 37–43% (chickpea), 54% (field pea) and 64% (faba bean) of total plant biomass N. Using only fixed-N inputs, we calculated EFs for chickpea (0.13–0.31%), field pea (0.18%) and faba bean (0.04%) that were significantly less than N-fertilised canola (0.48–0.78%) (P < 0.05), suggesting legume-fixed N is a less emissive form of N input to the soil than fertiliser N. Inputs of legume-fixed N should be more accurately quantified to properly gauge the potential for legumes to mitigate soil N2O emissions. EF’s from legume crops need to be revised and should include a factor for the proportion of the legume’s N derived from the atmosphere.

Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe

Soil microorganisms are critical to ecosystem functioning and the maintenance of soil fertility. However, despite global increases in the inputs of nitrogen (N) and phosphorus (P) to ecosystems due to human activities, we lack a predictive understanding of how microbial communities respond to elevated nutrient inputs across environmental gradients. Here we used high-throughput sequencing of marker genes to elucidate the responses of soil fungal, archaeal, and bacterial communities using an N and P addition experiment replicated at 25 globally distributed grassland sites. We also sequenced metagenomes from a subset of the sites to determine how the functional attributes of bacterial communities change in response to elevated nutrients. Despite strong compositional differences across sites, microbial communities shifted in a consistent manner with N or P additions, and the magnitude of these shifts was related to the magnitude of plant community responses to nutrient inputs. Mycorrhizal fungi and methanogenic archaea decreased in relative abundance with nutrient additions, as did the relative abundances of oligotrophic bacterial taxa. The metagenomic data provided additional evidence for this shift in bacterial life history strategies because nutrient additions decreased the average genome sizes of the bacterial community members and elicited changes in the relative abundances of representative functional genes. Our results suggest that elevated N and P inputs lead to predictable shifts in the taxonomic and functional traits of soil microbial communities, including increases in the relative abundances of faster-growing, copiotrophic bacterial taxa, with these shifts likely to impact belowground ecosystems worldwide.

Impacts of sodic soil amelioration on hydraulic conductivity and deep drainage in the Lower Burdekin

An understanding of the influence of soil chemistry on soil hydraulic properties is of critical importance for the management of sodic soils under irrigation. The hydraulic conductivity of sodic soils has been shown to be affected by properties of the applied solution including pH (Suarez et al. 1984), sodicity and salt concentration (McNeal and Coleman 1966). The changes in soil hydraulic conductivity are the result of changes in the spacing between clay layers in response to changes in soil solution chemistry. While the importance o f soil chemistry in controlling hydraulic conductivity is known, the exact impacts of sodic soil amelioration on hydraulic conductivity and deep drainage at a given location are difficult to predict. This is because the relationships between soil chemical factors and hydraulic conductivity are soil specific and because local site specific factors also need to be considered to determine the actual impacts on deep drainage rates.

Performance of soft-hard-soft (SHS) cement based composite subjected to blast loading with consideration of interface properties

This paper presents a combined experimental and numerical study on the damage and performance of a soft-hard-soft (SHS) multi-layer cement based composite subjected to blast loading which can be used for protective structures and infrastructures to resist extreme loadings, and the composite consists of three layers of construction materials including asphalt concrete (AC) on the top, high strength concrete (HSC) in the middle, and engineered cementitious composites (ECC) at the bottom. To better characterize the material properties under dynamic loading, interface properties of the composite were investigated through direct shear test and also used to validate the interface model. Strain rate effects of the asphalt concrete were also studied and both compressive and tensile dynamic increase factor (DIF) curves were improved based on split Hopkinson pressure bar (SHPB) test. A full-scale field blast test investigated the blast behavior of the composite materials. The numerical model was established by taking into account the strain rate effect of all concrete materials. Furthermore, the interface properties were also considered into the model. The numerical simulation using nonlinear finite element software LS-DYNA agrees closely with the experimental data. Both the numerical and field blast test indicated that the SHS composite exhibited high resistance against blast loading.

Determination of imidacloprid in paddy water and soil by liquid chromatography electrospray ionization-tandem mass spectrometry

A method for the determination of imidacloprid in paddy water and soil was developed using liquid chromatography electrospray ionization-tandem mass spectrometry (LC/ESI-MS/MS). Separation of imidacloprid was carried out on a Shimadzu C18 column (150 mm × 4.6 mm, 4.6 μm) with an acetonitrile?water (50 : 50, v/v) mobile phase containing 0.1% of acetic acid. The flow rate was 0.3 mL/min in isocratic mode. The product ion at 209 m/z was selected for quantification in multiple-reaction monitoring scan mode. Imidacloprid residues in soil were extracted by a solid-liquid extraction method with acetonitrile. Water samples were filtered and directly injected for analysis without extraction. Detection limits of 0.5 μg/kg and 0.3 μg/L were achieved for soil and water samples, respectively. The method had recoveries of 90 ± 2% (n = 4) for soil samples and 100 ± 2% (n = 4) for water samples. A linear relationship was observed throughout the investigated range of concentrations (1-200 μg/L), with the correlation coefficients ranging from 0.999 to 1.000. © Pleiades Publishing, Ltd., 2010.

Impacts of sodic soil amelioration on deep drainage

Groundwater tables are rising beneath irrigated fields in some areas of the Lower Burdekin in North Queensland, Australia. The soils where this occurs are predominantly sodic clay soils with low hydraulic conductivities. Many of these soils have been treated by applying gypsum or by increasing the salinity of irrigation water by mixing saline groundwater with fresh river water. While the purpose of these treatments is to increase infiltration into the surface soils and improve productivity of the root zone, it is thought that the treatments may have altered the soil hydraulic properties well below the root zone leading to increased groundwater recharge and rising water tables. In this paper we discuss the use of column experiments and HYDRUS modelling, with major ion reaction and transport and soil water chemistry-dependent hydraulic conductivity, to assess the likely depth, magnitude and timing of the impacts of surface soil amelioration on soil hydraulic properties below the root zone and hence groundwater recharge. In the experiments, columns of sodic clays from the Lower Burdekin were leached for extended periods of time with either gypsum solutions or mixed cation salt solutions and change s in hydraulic conductivity were measured. Leaching with a gypsum solution for an extended time period, until the flow rate stabilised, resulted in an approximately twenty fold increase in the hydraulic conductivity when compared with a low salinity, mixed cation solution. HYDRUS modelling was used to high light the role of those factors which might influence the impacts of soil treatment, particularly at depth, including the large amounts of rain during the relatively short wet season and the presence of thick low permeability clay layers.

Establishing turf grass increases soil greenhouse gas emissions in peri-urban environments

Urbanization is becoming increasingly important in terms of climate change and ecosystem functionality worldwide. We are only beginning to understand how the processes of urbanization influence ecosystem dynamics and how peri-urban environments contribute to climate change. Brisbane in South East Queensland (SEQ) currently has the most extensive urban sprawl of all Australian cities. This leads to substantial land use changes in urban and peri-urban environments and the subsequent gaseous emissions from soils are to date neglected for IPCC climate change estimations. This research examines how land use change effects methane (CH4) and nitrous oxide (N2O) fluxes from peri-urban soils and consequently influences the Global Warming Potential (GWP) of rural ecosystems in agricultural use undergoing urbanization. Therefore, manual and fully automated static chamber measurements determined soil gas fluxes over a full year and an intensive sampling campaign of 80 days after land use change. Turf grass, as the major peri-urban land cover, increased the GWP by 415 kg CO2-e ha 1 over the first 80 days after conversion from a well-established pasture. This results principally from increased daily average N2O emissions of 0.5 g N2O ha-1 d-1 from the pasture to 18.3 g N2O ha-1 d-1 from the turf grass due to fertilizer application during conversion. Compared to the native dry sclerophyll eucalypt forest, turf grass establishment increases the GWP by another 30 kg CO2-e ha 1. The results presented in this study clearly indicate the substantial impact of urbanization on soil-atmosphere gas exchange in form of non-CO2 greenhouse gas emissions particularly after turf grass establishment.

Optimal adaptive sequential spatial sampling of soil using pair-copulas

A spatial sampling design that uses pair-copulas is presented that aims to reduce prediction uncertainty by selecting additional sampling locations based on both the spatial configuration of existing locations and the values of the observations at those locations. The novelty of the approach arises in the use of pair-copulas to estimate uncertainty at unsampled locations. Spatial pair-copulas are able to more accurately capture spatial dependence compared to other types of spatial copula models. Additionally, unlike traditional kriging variance, uncertainty estimates from the pair-copula account for influence from measurement values and not just the configuration of observations. This feature is beneficial, for example, for more accurate identification of soil contamination zones where high contamination measurements are located near measurements of varying contamination. The proposed design methodology is applied to a soil contamination example from the Swiss Jura region. A partial redesign of the original sampling configuration demonstrates the potential of the proposed methodology.

Denitrification losses from an intensively managed sub-tropical pasture: Impact of soil moisture on the partitioning of N2 and N2O emissions

Intensively managed pastures in subtropical Australia under dairy production are nitrogen (N) loaded agro-ecosystems, with an increased pool of N available for denitrification. The magnitude of denitrification losses and N2:N2O partitioning in these agro-ecosystems is largely unknown, representing a major uncertainty when estimating total N loss and replacement. This study investigated the influence of different soil moisture contents on N2 and N2O emissions from a subtropical dairy pasture in Queensland, Australia. Intact soil cores were incubated over 15 days at 80% and 100% water-filled pore space (WFPS), after the application of 15N labelled nitrate, equivalent to 50 kg N ha−1. This setup enabled the direct quantification of N2 and N2O emissions following fertilisation using the 15N gas flux method. The main product of denitrification in both treatments was N2. N2 emissions exceeded N2O emissions by a factor of 8 ± 1 at 80% WFPS and a factor of 17 ± 2 at 100% WFPS. The total amount of N-N2 lost over the incubation period was 21.27 kg ± 2.10 N2-N ha−1 at 80% WFPS and 25.26 kg ± 2.79 kg ha−1 at 100% WFPS respectively. N2 emissions remained high at 100% WFPS, while related N2O emissions decreased. At 80% WFPS, N2 emissions increased constantly over time while N2O fluxes declined. Consequently, N2/(N2 + N2O) product ratios increased over the incubation period in both treatments. N2/(N2 + N2O) product ratios responded significantly to soil moisture, confirming WFPS as a key driver of denitrification. The substantial amount of fertiliser lost as N2 reveals the agronomic significance of denitrification as a major pathway of N loss for sub-tropical pastures at high WFPS and may explain the low fertiliser N use efficiency observed for these agro-ecosystems.

Non-linear response of soil N2O emissions to nitrogen fertiliser in a cotton-fallow rotation in sub-tropical Australia

Nitrogen fertiliser is a major source of atmospheric N2O and over recent years there is growing evidence for a non-linear, exponential relationship between N fertiliser application rate and N2O emissions. However, there is still high uncertainty around the relationship of N fertiliser rate and N2O emissions for many cropping systems. We conducted year-round measurements of N2O emission and lint yield in four N rate treatments (0, 90, 180 and 270 kg N ha-1) in a cotton-fallow rotation on a black vertosol in Australia. We observed a nonlinear exponential response of N2O emissions to increasing N fertiliser rates with cumulative annual N2O emissions of 0.55 kg N ha-1, 0.67kg N ha-1, 1.07 kg N ha-1 and 1.89 kg N ha-1 for the four respective N fertiliser rates while no N response to yield occurred above 180N. The N fertiliser induced annual N2O EF factors increased from 0.13% to 0.29% and 0.50% for the 90N, 180N and 270N treatments respectively, significantly lower than the IPCC Tier 1 default value (1.0 %). This non-linear response suggests that an exponential N2O emissions model may be more appropriate for use in estimating emission of N2O from soils cultivated to cotton in Australia. It also demonstrates that improved agricultural N management practices can be adopted in cotton to substantially reduce N2O emissions without affecting yield potential.

Mapping the risk of soil-transmitted helminthic infections in the Philippines

Background In order to increase the efficient allocation of soil-transmitted helminth (STH) disease control resources in the Philippines, we aimed to describe for the first time the spatial variation in the prevalence of A. lumbricoides, T. trichiura and hookworm across the country, quantify the association between the physical environment and spatial variation of STH infection and develop predictive risk maps for each infection. Methodology/Principal Findings Data on STH infection from 35,573 individuals across the country were geolocated at the barangay level and included in the analysis. The analysis was stratified geographically in two major regions: 1) Luzon and the Visayas and 2) Mindanao. Bayesian geostatistical models of STH prevalence were developed, including age and sex of individuals and environmental variables (rainfall, land surface temperature and distance to inland water bodies) as predictors, and diagnostic uncertainty was incorporated. The role of environmental variables was different between regions of the Philippines. This analysis revealed that while A. lumbricoides and T. trichiura infections were widespread and highly endemic, hookworm infections were more circumscribed to smaller foci in the Visayas and Mindanao. Conclusions/Significance This analysis revealed significant spatial variation in STH infection prevalence within provinces of the Philippines. This suggests that a spatially targeted approach to STH interventions, including mass drug administration, is warranted. When financially possible, additional STH surveys should be prioritized to high-risk areas identified by our study in Luzon.