Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Nitrogen fertilizer (15N) leaching in a central pivot fertigated coffee crop

Nitrogen has a complex dynamics in the soil-plant-atmosphere system. N fertilizers are subject to chemical and microbial transformations in soils that can result in significant losses. Considering the cost of fertilizers, the adoption of good management practices like fertigation could improve the N use efficiency by crops. Water balances (WB) were applied to evaluate fertilizer N leaching using 15N labeled urea in west Bahia, Brazil. Three scenarios (2008/2009) were established: i) rainfall + irrigation the full year, ii) rainfall only; and iii) rainfall + irrigation only in the dry season. The water excess was considered equal to the deep drainage for the very flat area (runoff = 0) with a water table located several meters below soil surface (capillary rise = 0). The control volume for water balance calculations was the 0 - 1 m soil layer, considering that it involves the active root system. The water drained below 1 m was used to estimate fertilizer N leaching losses. WB calculations used the mathematic model of Penman-Monteith for evapotranspiration, considering the crop coefficient equal to unity. The high N application rate associated to the high rainfall plus irrigation was found to be the main cause for leaching, which values were 14.7 and 104.5 kg ha-1 for the rates 400 and 800 kg ha-1 of N, corresponding to 3.7 and 13.1 % of the applied fertilizer, respectively.

Effects of Soil Compaction and Organic Carbon Content on Preferential Flow in Loamy Field Soils

Abstract: Preferential flow and transport through macropores affect plant water use efficiency and enhance leaching of agrochemicals and the transport of colloids, thereby increasing the risk for contamination of groundwater resources. The effects of soil compaction, expressed in terms of bulk density (BD), and organic carbon (OC) content on preferential flow and transport were investigated using 150 undisturbed soil cores sampled from 15 × 15–m grids on two field sites. Both fields had loamy textures, but one site had significantly higher OC content. Leaching experiments were conducted in each core by applying a constant irrigation rate of 10 mm h−1 with a pulse application of tritium tracer. Five percent tritium mass arrival times and apparent dispersivities were derived from each of the tracer breakthrough curves and correlated with texture, OC content, and BD to assess the spatial distribution of preferential flow and transport across the investigated fields. Soils from both fields showed strong positive correlations between BD and preferential flow. Interestingly, the relationships between BD and tracer transport characteristics were markedly different for the two fields, although the relationship between BD and macroporosity was nearly identical. The difference was likely caused by the higher contents of fines and OC at one of the fields leading to stronger aggregation, smaller matrix permeability, and a more pronounced pipe-like pore system with well-aligned macropores.

Soil phosphorus fractions in sandy soils amended with cattle manure for long periods

Phosphorus fractions were determined in soil samples from areas fertilized or not with farmyard cattle manure (FYM) and in samples of FYM used in the semi-arid region of Paraiba state, Brazil. Soil samples were taken from the 0-20; 20-40 and 40-60 cm layers of 18 cultivated areas, which, according to interviews with farmers, had been treated with 12 to 20 t ha-1 FYM annually, for the past 2 to 40 years. Soil samples were also collected from four unfertilized pasture areas as controls. Phosphorus in the soil samples was sequentially extracted with water (Pw), resin (Pres), NaHCO3 (Pi bic and Po bic), NaOH (Pi hid and Po hid), H2SO4 (Pacid) and, finally, by digestion with H2SO4/H2O2 (Presd). Nine FYM samples were extracted with water, resin, Mehlich-1, H2SO4, NaOH or digestion with H2SO4/H2O2, not sequentially, and the extracts analyzed for P. The sampled areas had homogeneous, sandy and P-deficient soils; increases in total soil P (Pt) above the mean value of the control areas (up to 274 mg kg-1 in the 0-20 cm layer of the most P-enriched samples) were therefore attributed to FYM applications, which was the only external P input in the region. Regression analysis was used to study the relationship between soil P fractions and Pt. The Pacid fraction, related to Ca-P forms, showed the greatest increases (p < 0.01) as a result of FYM applications, rising from 8.4 mg kg-1 in a non-fertilized sample to 43.8 mg kg-1 in the sample with the highest Pt content. The sum of Pw, Pres and Pi bic, considered as labile P, showed comparable increases with Pacid, while Pi hid showed the smallest increase due to FYM applications. Organic P forms also increased, more so the fraction Po hid, considered less labile, than the more labile one, Po bic. The residual P fraction was practically half of Pt, independently of the Pt value. Increases in labile P, Pacid and organic P were justified by the high average concentration of Pw (36 %), Pacid (34 %), and Po hid (30 %) in the FYM. Significant changes in the proportion of P forms among soil layers indicated the downward movement of P in organic forms.

Mobility of inorganic and organic phosphorus forms under different levels of phosphate and poultry litter fertilization in soils

The eutrophication of aquifers is strongly linked to the mobility of P in soils. Although P mobility was considered irrelevant in a more distant past, more recent studies have shown that P, both in organic (Po) and inorganic forms (Pi), can be lost by leaching and eluviation through the soil profile, particularly in less weathered and/or sandier soils with low P adsorption capacity. The purpose of this study was to determine losses of P forms by leaching and eluviation from soil columns. Each column consisted of five PVC rings (diameter 5 cm, height 10 cm), filled with two soil types: a clayey Red-Yellow Latosol and a sandy loam Red-Yellow Latosol, which were exposed to water percolation. The soils were previously treated with four P rates (as KH2PO4 ) to reach 0, 12.5, 25.0 and 50 % of the maximum P adsorption capacity (MPAC). The P source was homogenized with the whole soil volume and incubated for 60 days. After this period the soils were placed in the columns; the soil of the top ring was mixed with five poultry litter rates of 0, 20, 40, 80, and 160 t ha-1 (dry weight basis). Treatments consisted of a 4 x 5 x 2 factorial scheme corresponding to four MPAC levels, five poultry litter rates, two soils, with three replications, arranged in a completely randomized block design. Deionized water was percolated through the columns 10 times in 35 days to simulate about 1,200 mm rainfall. In the leachate of each column the inorganic P (reactive P, Pi) and organic P forms (unreactive P, Po) were determined. At the end of the experiment, the columns were disassembled and P was extracted with the extractants Mehlich-1 (HCl 0.05 mol L-1 and H2SO4 0.0125 mol L-1) and Olsen (NaHCO3 0.5 mol L-1; pH 8.5) from the soil of each ring. The Pi and Po fractions were measured by the Olsen extractant. It was found that under higher poultry litter rates the losses of unreactive P (Po) were 6.4 times higher than of reactive P (Pi). Both the previous P fertilization and increasing poultry litter rates caused a vertical movement of P down the soil columns, as verified by P concentrations extracted by Mehlich-1 and NaHCO3 (Olsen). The environmental critical level (ECL), i.e., the P soil concentration above which P leaching increases exponentially, was 100 and 150 mg dm-3 by Mehlich-1 and 40 and 60 mg dm-3 by Olsen, for the sandy loam and clay soils, respectively. In highly weathered soils, where residual P is accumulated by successive crops, P leaching through the profile can be significant, particularly when poultry litter is applied as fertilizer.

Cation dynamics in soils with different salinity levels growing irrigated rice

Salinity levels in soils of the Outer Coastal Plain of Rio Grande do Sul, Brazil, can be high, due to excess of Na in the irrigation water, evapotranspiration and soil development from marine sediments. The cultivation of irrigated rice could be an alternative, since ion uptake as well as leaching by the establishment of a water layer could mitigate the effects of soil salinity. This study aimed to evaluate the dynamics of basic cations in the solution of Albaqualf soils with different salinity levels growing irrigated rice. The plow layer contained exchangeable Na percentages (ESP) of 5.6, 9.0, 21.2 and 32.7 %. The plant stand, dry matter, Na, K and Ca + Mg uptake at full flowering and grain yield were evaluated. The levels of Na, K, Ca + Mg and electrical conductivity (EC) in the soil solution were also measured weekly during the rice cycle at four soil depths, in the water layer and irrigation water. The Na, K and Ca + Mg uptake by rice at full flowering was used to estimate ion depletion from the layer under root influence. Soil salinity induced a reduction in the rice stand, especially in the soil with ESP of 32.7 %, resulting in lower cation uptake and very low yield at that site. As observed in the water layer and irrigation water, the Na, K, Ca + Mg and EC levels in the soil solution decreased with time at depths of 5, 10 and 20 cm, regardless of the original soil salinity, showing that cation dynamics in the plow layer was determined by leaching and root uptake, rather than by the effect of evapoconcentration of basic cations in the soil surface layer.

Physical and chemical attributes of archaeological soils developed from shell middens in the Região dos Lagos, Rio de Janeiro, Brazil

In prehistoric times, innumerous shell middens, called "sambaquis", consisting mainly of remains of marine organisms, were built along the Brazilian coast. Although the scientific community took interest in these anthropic formations, especially since the nineteenth century, their pedological context is still poorly understood. The purpose of this study was to characterize and identify the physical and chemical changes induced by soil-forming processes, as well as to compare the morphology of shell midden soils with other, already described, anthropogenic soils of Brazil. Four soil profiles developed from shell middens in the Região dos Lagos - RJ were morphologically described and the physical and chemical properties determined. The chemical analysis showed that Ca, Mn, Mg, and particularly P and Zn are indicators of anthropic horizons of midden soils, as in the Amazon Dark Earths (Terras Pretas de Índio). After the deposition of P-rich material, P reaction and leaching can mask or disturb the evidence of in situ man-made strata, but mineralogical and chemical studies of phosphate forms can elucidate the apparent complexity. Lower phosphate-rich strata without direct anthropic inputs indicate P leaching and precipitation in secondary forms. The total and bioavailable contents of Ca, Mg, Zn, Mn, Cu, P, and organic C of midden soils were much higher than of regional soils without influence of ancient human settlements, demonstrating that the high fertility persisted for long periods, at some sites for more than 4000 years. The physical analysis showed that wind-blown sand contributed significantly to increase the sand fraction in the analyzed soils (texture classes sand, sandy loam and sandy clay loam) and that the aeolian sand accumulation occurred simultaneously with the midden formation.

Nitrate role in basic cation leaching under no-till

Especially under no-tillage, subsuface soil acidity has been a problem, because it depends on base leaching, which has been associated with the presence of low molecular weigth organic acids and companion anions. The objective of this study was to evaluate exchangeable base cation leaching as affected by surface liming along with annual urea side-dressing of maize and upland rice. Treatments consisted of four lime rates (0, 1500, 3000, and 6000 kg ha-1) combined with four nitrogen rates (0, 50, 100, and 150 kg ha-1) applied to maize (Zea mays) and upland rice (Oryza sativa), in two consecutive years. Maize was planted in December, three months after liming. In September of the following year, pearl millet (Pennisetum glaucum) was planted without fertilization and desiccated 86 days after plant emergence. Afterwards, upland rice was grown. Immediately after upland rice harvest, 18 months after surface liming, pH and N-NO3-, N-NH4+, K, Ca, and Mg levels were evaluated in soil samples taken from the layers 0-5, 5-10, 10-20 and 20-40 cm. Higher maize yields were obtained at higher N rates and 3000 kg ha-1 lime. Better results for upland rice and pearl millet yields were also obtained with this lime rate, irrespective of N levels. The vertical mobility of K, Ca and Mg was higher in the soil profiles with N fertilization. Surface liming increased pH in the upper soil layers causing intense nitrate production, which was leached along with the base cations.

Chemical properties of soils treated with biological sludge from gelatin industry

The impact of agro-industrial organic wastes in the environment can be reduced when used in agriculture. From the standpoint of soil fertility, residue applications can increase the organic matter content and provide nutrients for plants. This study evaluated the effect of biological sludge from gelatin industry on the chemical properties of two Ultisols (loamy sand and sandy clay) and an Oxisol (clay). The experiment lasted 120 days and was carried out in laboratory in a completely randomized design with factorial arrangement, combining the three soils and six biological sludge rates (0, 100, 200, 300, 400, and 500 m³ ha-1), with three replications. Biological sludge rates of up to 500 m³ ha-1 decreased soil acidity and increased the effective cation exchange capacity (CEC) and N, Ca, Mg, and P availability, without exceeding the tolerance limit for Na. The increase in exchangeable base content, greater than the effective CEC, indicates that the major part of cations added by the sludge remains in solution and can be lost by leaching.

Estimating subsoil resistance to nitrate leaching from easily measurable pedological properties

Leaching of nitrate (NO3-) can increase the groundwater concentration of this anion and reduce the agronomical effectiveness of nitrogen fertilizers. The main soil property inversely related to NO3- leaching is the anion exchange capacity (AEC), whose determination is however too time-consuming for being carried out in soil testing laboratories. For this reason, this study evaluated if more easily measurable soil properties could be used to estimate the resistance of subsoils to NO3- leaching. Samples from the subsurface layer (20-40 cm) of 24 representative soils of São Paulo State were characterized for particle-size distribution and for chemical and electrochemical properties. The subsoil content of adsorbed NO3- was calculated from the difference between the NO3- contents extracted with 1 mol L-1 KCl and with water; furthermore, NO3- leaching was studied in miscible displacement experiments. The results of both adsorption and leaching experiments were consistent with the well-known role exerted by AEC on the nitrate behavior in weathered soils. Multiple regression analysis indicated that in subsoils with (i) low values of remaining phosphorus (Prem), (ii) low soil pH values measured in water (pH H2O), and (iii) high pH values measured in 1 moL L-1 KCl (pH KCl), the amounts of surface positive charges tend to be greater. For this reason, NO3- leaching tends to be slower in these subsoils, even under saturated flow condition.

Boron Leaching Decreases withIncreases on Soil pH

ABSTRACT Management of boron fertilization depends on the magnitude of B leaching in the soil profile, which varies proportionally with the concentration of B in the soil solution, which, in turn, decreases as the soil pH increases due to the higher sorption of B on soil solid surfaces. The objective of this study was to quantify the effect of liming and rates of B applied to the soil on B leaching. The experiment was carried out in the laboratory in 2012, and treatments consisted of a factorial combination of two rates of liming (without and with lime to raise the soil pH to 6.0) and five rates of B (0, 10, 20, 50 and 100 mg kg-1, as boric acid). A Typic Rhodudalf was used, containing 790 g kg-1 clay and 23 g kg-1 organic matter; the pH(H2O) was 4.6. Experimental units were composed of PVC leaching columns (0.10 m in diameter) containing 1.42 kg of soil (dry base). Boron was manually mixed with the top 0.15 m of the soil. After that, every seven days for 15 weeks, 300 mL of distilled water were added to the top of each column. In the percolated solution, both the volume and concentration of B were measured. Leaching of B decreased with increased soil pH and, averaged across the B rates applied, was 58 % higher from unlimed (pH 4.6) than from limed (pH 6.6) samples as a result of the increase in B sorption with higher soil pH. In spite of its high vertical mobility, the residual effect of B was high in this oxisol, mainly in the limed samples where 80 % of B applied at the two highest rates remained in the soil, even after 15 water percolations. Total recovery of applied B, including leached B plus B extracted from the soil after all percolations, was less than 50 %, showing that not all sorbed B was quantified by the hot water extraction method.

Metal leaching from a limestone quarry minesoil made with thermally-dried sewage sludge

Peer Reviewed

Environmental behavior of soils and mixtures of soil-whitewash mud

The present study sought to observe the behavior of soils in natural state and in mixtures, in different ratios, with the industrial solid residue called whitewash mud. The work was conducted with samples of typical soils from the region of Alagoinhas, Bahia-Brazil. Wet chemical analysis and atomic absorption spectrophotometry were used in order to obtain the classification of the industrial solid residue. Solubilization and leaching tests were performed and X-ray diffraction and electron microscopy techniques were carried out. The results showed that the whitewash mud was classified as non-inert, but with great capacity of heavy metal retention largely owed to the kaolinite and goethite presence in the clay fraction of the soils, making it difficult to have heavy metals readily available for exchange.

Leaching of Picloram in ultisol under different rainfall volumes

The present study aimed to evaluate the leaching potential of Picloram in Ultisol columns under different rainfall amounts. For such, 30 treatments were evaluated (one soil associated with three levels of rainfall and ten depths).The experiments were arranged in a split-plot design, in a completely randomized design, with four replications. PVC columns of 10 cm in diameter and 50 cm in length were filled with these soils, moistened, and placed upright for 48 hours to drain the excess water. The herbicide was applied and rainfall simulations were carried out at specified intensities, according to the treatments, to force Picloram leaching. After 72 hours, all the columns were arranged in a horizontal position and opened lengthwise. Then, soil sampling was carried out every 5 cm of depth for subsequent herbicide extraction and quantification and analysis by high performance liquid chromatography. The remaining soil samples were placed in plastic pots, and, at the respective depths, the indicator species Cucumis sativus was sown. Twenty-one days after the emergence (DAE) of the indicator plants, evaluations were conducted to verify the symptoms of toxicity caused by Picloram in the plants. It was concluded that Picloram leaching is directly dependent on the volume of rain applied. The herbicide reached the deepest regions in the soil with the highest intensity of rain. The results obtained by bioassay were in agreement with those found by liquid chromatography.

Application of biogas slurries from energy crops to arable soils and their impact on carbon and nitrogen dynamics

The use of renewable primary products as co-substrate or single substrate for biogas production has increased consistently over the last few years. Maize silage is the preferential energy crop used for fermentation due to its high methane (CH4) yield per hectare. Equally, the by-product, namely biogas slurry (BS), is used with increasing frequency as organic fertilizer to return nutrients to the soil and to maintain or increase the organic matter stocks and soil fertility. Studies concerning the application of energy crop-derived BS on the carbon (C) and nitrogen (N) mineralization dynamics are scarce. Thus, this thesis focused on the following objectives: I) The determination of the effects caused by rainfall patterns on the C and N dynamics from two contrasting organic fertilizers, namely BS from maize silage and composted cattle manure (CM), by monitoring emissions of nitrous oxide (N2O), carbon dioxide (CO2) and CH4 as well as leaching losses of C and N. II) The investigation of the impact of differences in soil moisture content after the application of BS and temperature on gaseous emissions (CO2, N2O and CH4) and leaching of C and N compounds. III) A comparison of BS properties obtained from biogas plants with different substrate inputs and operating parameters and their effect on C and N dynamics after application to differently textured soils with varying application rates and water contents. For the objectives I) and II) two experiments (experiment I and II) using undisturbed soil cores of a Haplic Luvisol were carried out. Objective III) was studied on a third experiment (experiment III) with disturbed soil samples. During experiment I three rainfall patterns were implemented including constant irrigation, continuous irrigation with periodic heavy rainfall events, and partial drying with rewetting periods. Biogas slurry and CM were applied at a rate of 100 kg N ha-1. During experiment II constant irrigation and an irrigation pattern with partial drying with rewetting periods were carried out at 13.5°C and 23.5°C. The application of BS took place either directly before a rewetting period or one week after the rewetting period stopped. Experiment III included two soils of different texture which were mixed with ten BS’s originating from ten different biogas plants. Treatments included low, medium and high BS-N application rates and water contents ranging from 50% to 100% of water holding capacity (WHC). Experiment I and II showed that after the application of BS cumulative N2O emissions were 4 times (162 mg N2O-N m-2) higher compared to the application of CM caused by a higher content of mineral N (Nmin) in the form of ammonium (NH4+) in the BS. The cumulative emissions of CO2, however, were on the same level for both fertilizers indicating similar amounts of readily available C after composting and fermentation of organic material. Leaching losses occurred predominantly in the mineral form of nitrate (NO3-) and were higher in BS amended soils (9 mg NO3--N m-2) compared to CM amended soils (5 mg NO3--N m-2). The rainfall pattern in experiment I and II merely affected the temporal production of C and N emissions resulting in reduced CO2 and enhanced N2O emissions during stronger irrigation events, but showed no effect on the cumulative emissions. Overall, a significant increase of CH4 consumption under inconstant irrigation was found. The time of fertilization had no effect on the overall C and N dynamics. Increasing temperature from 13.5°C to 23.5°C enhanced the CO2 and N2O emissions by a factor of 1.7 and 3.7, respectively. Due to the increased microbial activity with increasing temperature soil respiration was enhanced. This led to decreasing oxygen (O2) contents which in turn promoted denitrification in soil due to the extension of anaerobic microsites. Leaching losses of NO3- were also significantly affected by increasing temperature whereas the consumption of CH4 was not affected. The third experiment showed that the input materials of biogas plants affected the properties of the resulting BS. In particular the contents of DM and NH4+ were determined by the amount of added plant biomass and excrement-based biomass, respectively. Correlations between BS properties and CO2 or N2O emissions were not detected. Solely the ammonia (NH3) emissions showed a positive correlation with NH4+ content in BS as well as a negative correlation with the total C (Ct) content. The BS-N application rates affected the relative CO2 emissions (% of C supplied with BS) when applied to silty soil as well as the relative N2O emissions (% of N supplied with BS) when applied to sandy soil. The impacts on the C and N dynamics induced by BS application were exceeded by the differences induced by soil texture. Presumably, due to the higher clay content in silty soils, organic matter was stabilized by organo-mineral interactions and NH4+ was adsorbed at the cation exchange sites. Different water contents induced highest CO2 emissions and therefore optimal conditions for microbial activity at 75% of WHC in both soils. Cumulative nitrification was also highest at 75% and 50% of WHC whereas the relative N2O emissions increased with water content and showed higher N2O losses in sandy soils. In summary it can be stated that the findings of the present thesis confirmed the high fertilizer value of BS’s, caused by high concentrations of NH4+ and labile organic compounds such as readily available carbon. These attributes of BS’s are to a great extent independent of the input materials of biogas plants. However, considerably gaseous and leaching losses of N may occur especially at high moisture contents. The emissions of N2O after field application corresponded with those of animal slurries.