Methane and nitrous oxide fluxes in annual and perennial land-use systems of the irrigated areas in the Aral Sea Basin

Land use and agricultural practices can result in important contributions to the global source strength of atmospheric nitrous oxide (N2O) and methane (CH4). However, knowledge of gas flux from irrigated agriculture is very limited. From April 2005 to October 2006, a study was conducted in the Aral Sea Basin, Uzbekistan, to quantify and compare emissions of N2O and CH4 in various annual and perennial land-use systems: irrigated cotton, winter wheat and rice crops, a poplar plantation and a natural Tugai (floodplain) forest. In the annual systems, average N2O emissions ranged from 10 to 150 μg N2O-N m−2 h−1 with highest N2O emissions in the cotton fields, covering a similar range of previous studies from irrigated cropping systems. Emission factors (uncorrected for background emission), used to determine the fertilizer-induced N2O emission as a percentage of N fertilizer applied, ranged from 0.2% to 2.6%. Seasonal variations in N2O emissions were principally controlled by fertilization and irrigation management. Pulses of N2O emissions occurred after concomitant N-fertilizer application and irrigation. The unfertilized poplar plantation showed high N2O emissions over the entire study period (30 μg N2O-N m−2 h−1), whereas only negligible fluxes of N2O (<2 μg N2O-N m−2 h−1) occurred in the Tugai. Significant CH4 fluxes only were determined from the flooded rice field: Fluxes were low with mean flux rates of 32 mg CH4 m−2 day−1 and a low seasonal total of 35.2 kg CH4 ha−1. The global warming potential (GWP) of the N2O and CH4 fluxes was highest under rice and cotton, with seasonal changes between 500 and 3000 kg CO2 eq. ha−1. The biennial cotton–wheat–rice crop rotation commonly practiced in the region would average a GWP of 2500 kg CO2 eq. ha−1 yr−1. The analyses point out opportunities for reducing the GWP of these irrigated agricultural systems by (i) optimization of fertilization and irrigation practices and (ii) conversion of annual cropping systems into perennial forest plantations, especially on less profitable, marginal lands.

Greenhouse gas emissions from sub-tropical agricultural soils after addition of organic by-products

As the cost of mineral fertilisers increases globally, organic soil amendments (OAs) from agricultural sources are increasingly being used as substitutes for nitrogen. However, the impact of OAs on the production of greenhouse gases (CO2 and N2O) is not well understood. A 60-day laboratory incubation experiment was conducted to investigate the impacts of applying OAs (equivalent to 296 kg N ha−1 on average) on N2O and CO2 emissions and soil properties of clay and sandy loam soils from sugar cane production. The experiment included 6 treatments, one being an un-amended (UN) control with addition of five OAs being raw mill mud (MM), composted mill mud (CM), high N compost (HC), rice husk biochar (RB), and raw mill mud plus rice husk biochar (MB). These OAs were incubated at 60, 75 and 90% water-filled pore space (WFPS) at 25°C with urea (equivalent to 200 kg N ha−1) added to the soils thirty days after the incubation commenced. Results showed WFPS did not influence CO2 emissions over the 60 days but the magnitude of emissions as a proportion of C applied was RB < CM < MB < HC

Toward a model for predicting magnetic susceptibility of bedrock regolith and soils

One of the Department of Defense's most pressing environmental problems is the efficient detection and identification of unexploded ordnance (UXO). In regions of highly magnetic soils, magnetic and electromagnetic sensors often detect anomalies that are of geologic origin, adding significantly to remediation costs. In order to develop predictive models for magnetic susceptibility, it is crucial to understand modes of formation and the spatial distribution of different iron oxides. Most rock types contain iron and their magnetic susceptibility is determined by the amount and form of iron oxides present. When rocks weather, the amount and form of the oxides change, producing concomitant changes in magnetic susceptibility. The type of iron oxide found in the weathered rock or regolith is a function of the duration and intensity of weathering, as well as the original content of iron in the parent material. The rate of weathering is controlled by rainfall and temperature; thus knowing the climate zone, the amount of iron in the lithology and the age of the surface will help predict the amount and forms of iron oxide. We have compiled analyses of the types, amounts, and magnetic properties of iron oxides from soils over a wide climate range, from semi arid grasslands, to temperate regions, and tropical forests. We find there is a predictable range of iron oxide type and magnetic susceptibility according to the climate zone, the age of the soil and the amount of iron in the unweathered regolith.

Mineralogy and Magnetic Properties of Basaltic Substrate Soils: Kaho'olawe and Big Island, Hawaii

Magnetic behavior of soils can seriously hamper the performance of geophysical sensors. Currently, we have little understanding of the types of minerals responsible for the magnetic behavior, as well as their distribution in space and evolution through time. This study investigated the magnetic characteristics and mineralogy of Fe-rich soils developed on basaltic substrate in Hawaii. We measured the spatial distribution of magnetic susceptibility (χlf) and frequency dependence (χfd%) across three test areas in a well-developed eroded soil on Kaho'olawe and in two young soils on the Big Island of Hawaii. X-ray diffraction spectroscopy, x-ray fluorescence spectroscopy (XFCF), chemical dissolution, thermal analysis, and temperature-dependent magnetic studies were used to characterize soil development and mineralogy for samples from soil pits on Kaho'olawe, surface samples from all three test areas, and unweathered basalt from the Big Island of Hawaii. The measurements show a general increase in magnetic properties with increasing soil development. The XRF Fe data ranged from 13% for fresh basalt and young soils on the Big Island to 58% for material from the B horizon of Kaho'olawe soils. Dithionite-extractable and oxalate-extractable Fe percentages increase with soil development and correlate with χlf-and χfd%, respectively. Results from the temperature-dependent susceptibility measurements show that the high soil magnetic properties observed in geophysical surveys in Kaho'olawe are entirely due to neoformed minerals. The results of our studies have implications for the existing soil survey of Kaho'olawe and help identify methods to characterize magnetic minerals in tropical soils.

Are we getting accurate measurements of Ksat for sodic clay soils?

In this paper we discuss the use of a series of column experiments to improve understanding of the effect irrigation water chemistry (saline solutions) has on measurements of saturated hydraulic conductivity (Ksat) of a sodic clay soil. We highlight in particular the use of extended leaching periods to determine whether the duration of leaching affects the results. In the experiments, mixed cation solutions of two different salinity levels, 50 meq/L and 100 meq/L, were applied under constant head to columns of a repacked sodic clay soil using three replicates for each treatment. The maximum Ksat measured during leaching with the 100 meq/L solution was approximately double the maximum Ksat measured during leaching with the 50 meq/L solution. Measured flow rates were found to increase rapidly after flow commenced then decrease gradually until flow rates became stable. The final, stable flow rate was roughly 80% less than the maximum flow rate measured. Reasons for these changes in saturated hydraulic conductivity are discussed. The key finding from these experiments is that long term leaching, involving significantly more pore volumes than is commonly reported in the literature, is required to obtain a ‘stable’ Ksat. We recommend that further studies be carried out to (1) determine whether similar behaviour in Ksat occurs in a wide range of sodic clay soils and (2) to help build a better understanding of the causes and implications of the observed behaviour in Ksat.

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.

Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

Consolidation Behavior of Soils

Based on Terzaghi's consolidation theory, percent of consolidation, U, versus the time factor, T, relationship for constant/linear excess pore water pressure distribution, it is possible to generate theoretical log10(H2/t) versus U curves where H is the length of the drainage path of a consolidating layer, and t is the time for different known values of the coefficient of consolidation, cν. A method has been developed wherein both the theoretical and experimental behavior of soils during consolidation can be simultaneously compared and studied on the same plot. The experimental log10(H2/t) versus U curves have been compared with the theoretical curves. The deviations of the experimental behavior from the theory are explained in terms of initial compression and secondary compression. Analysis of results indicates that the secondary compression essentially starts from about 60% consolidation. A simple procedure is presented for calculating the value of cv from the δ-t data using log10(H2/t) versus U plot.

Critical Reappraisal Of Plasticity Index Of Soils

A fundamental approach, based on Gouy-Chapman theory of double layer, has been provided to micromechanistically interpret the plasticity index of soils and their relationship with liquid limit. The relationships between plasticity index and liquid limit, developed earlier, through statistical approaches and critical state concepts, have been reexamined. The statistical analysis of extensive published data has resulted in the relationship, IP = 0.74 (wL - 8). On comparison with other relationships in vogue the proposed equation has been found to give better agreement. From the reappraisal of critical state approaches consistent with the micromechanistic interpretation, the possible range of parameters have been computed and compared with those obtained by statistical means to enhance the credibility of the proposed relationship.

Risk management strategies using seasonal climate forecasting in irrigated cotton production: a tale of stochastic dominance.

Decision-making in agriculture is carried out in an uncertain environment with farmers often seeking information to reduce risk. As a result of the extreme variability of rainfall and stream-flows in north-eastern Australia, water supplies for irrigated agriculture are a limiting factor and a source of risk. The present study examined the use of seasonal climate forecasting (SCF) when calculating planting areas for irrigated cotton in the northern Murray Darling Basin. Results show that minimising risk by adjusting plant areas in response to SCF can lead to significant gains in gross margin returns. However, how farmers respond to SCF is dependent on several other factors including irrigators’ attitude towards risk.

Pig effluent-P application can increase risk of P transport: two case studies

Land application of piggery effluent (containing urine, faeces, water, and wasted feed) is under close scrutiny as a potential source of water resource contamination with phosphorus (P). This paper investigates two case studies of the impact of long-term piggery effluent-P application to soil. A Natrustalf (Sodosol) at P1 has received a net load of 3700 kg effluent P/ha over 19 years. The Haplustalf (Dermosol) selected (P2) has received a net load of 310 000 kg P/ha over 30 years. Total, bicarbonate extractable, and soluble P forms were determined throughout the soil profiles for paired (irrigated and unirrigated) sites at P1 and P2, as well as P sorption and desorption characteristics. Surface bicarbonate (PB, 0 - 0.05 m depth) and dilute CaCl2 extractable molybdate-reactive P (PC) have been significantly elevated by effluent irrigation (P1: PB unirrigated 23±1, irrigated 290±6; PC unirrigated 0.03±0.00, irrigated 23.9±0.2. P2: PB unirrigated 72±48, irrigated 3950±1960; PC unirrigated 0.7±0.0, irrigated 443±287 mg P/kg; mean±s.d.). Phosphorus enrichment to 1.5 m, detected as PB, was observed at P2. Elevated concentrations of CaCl2 extractable organic P forms (POC; estimated by non-molybdate reactive P in centrifuged supernatants) were observed from the soil surface of P1 to a depth of 0.4 m. Despite the extent of effluent application at both of these sites, only P1 displayed evidence of significant accumulation of POC. The increase in surface soil total P (0 - 0.05 m depth) due to effluent irrigation was much greater than laboratory P sorption (>25 times for P1; >57 times for P2) for a comparable range of final solution concentrations (desorption extracts ranged from 1-5 mg P/L for P1 and 50-80 mg P/L for P2). Precipitation of sparingly soluble P phases was evidenced in the soils of the P2 effluent application area.

Direct shear testing on unsaturated silty soils to investigate the effects of drying and wetting on shear strength parameters at low suction

A modified conventional direct shear device was used to measure unsaturated shear strength of two silty soils at low suction values (0 ~ 50 kPa) that were achieved by following drying and wetting paths of soil water characteristic curves (SWCCs). The results revealed that the internal friction angle of the soils was not significantly affected by either the suction or the drying wetting SWCCs. The apparent cohesion of soil increased with a decreasing rate as suction increased. Shear stress-shear displacement curves obtained from soil specimens subjected to the same net normal stress and different suction values showed a higher initial stiffness and a greater peak stress as suction increased. A soil in wetting exhibited slightly higher peak shear stress and more contractive volume change behavior than that of soil in drying at the same net normal stress and suction.

Spectral and thermal sensing for nitrogen and water status in rainfed and irrigated wheat environments

Variable-rate technologies and site-specific crop nutrient management require real-time spatial information about the potential for response to in-season crop management interventions. Thermal and spectral properties of canopies can provide relevant information for non-destructive measurement of crop water and nitrogen stresses. In previous studies, foliage temperature was successfully estimated from canopy-scale (mixed foliage and soil) temperatures and the multispectral Canopy Chlorophyll Content Index (CCCI) was effective in measuring canopy-scale N status in rainfed wheat (Triticum aestivum L.) systems in Horsham, Victoria, Australia. In the present study, results showed that under irrigated wheat systems in Maricopa, Arizona, USA, the theoretical derivation of foliage temperature unmixing produced relationships similar to those in Horsham. Derivation of the CCCI led to an r2 relationship with chlorophyll a of 0.53 after Zadoks stage 43. This was later than the relationship (r2 = 0.68) developed for Horsham after Zadoks stage 33 but early enough to be used for potential mid-season N fertilizer recommendations. Additionally, ground-based hyperspectral data estimated plant N (g kg)1) in Horsham with an r2 = 0.86 but was confounded by water supply and N interactions. By combining canopy thermal and spectral properties, varying water and N status can potentially be identified eventually permitting targeted N applications to those parts of a field where N can be used most efficiently by the crop.