Mats Larsbo

Applicability of models to predict phosphorus losses in drained fields: a review.

Most phosphorus (P) modeling studies of water quality have focused on surface runoff loses. However, a growing number of experimental studies have shown that P losses can occur in drainage water from artificially drained fields. In this review, we assess the applicability of nine models to predict this type of P loss. A model of P movement in artificially drained systems will likely need to account for the partitioning of water and P into runoff, macropore flow, and matrix flow. Within the soil profile, sorption and desorption of dissolved P and filtering of particulate P will be important. Eight models are reviewed (ADAPT, APEX, DRAINMOD, HSPF, HYDRUS, ICECREAMDB, PLEASE, and SWAT) along with P Indexes. Few of the models are designed to address P loss in drainage waters. Although the SWAT model has been used extensively for modeling P loss in runoff and includes tile drain flow, P losses are not simulated in tile drain flow. ADAPT, HSPF, and most P Indexes do not simulate flow to tiles or drains. DRAINMOD simulates drains but does not simulate P. The ICECREAMDB model from Sweden is an exception in that it is designed specifically for P losses in drainage water. This model seems to be a promising, parsimonious approach in simulating critical processes, but it needs to be tested. Field experiments using a nested, paired research design are needed to improve P models for artificially drained fields. Regardless of the model used, it is imperative that uncertainty in model predictions be assessed.

MACRO (v5.2): Model Use, Calibration, and Validation

MACRO is a one-dimensional dual-permeability model of variably saturated water flow and reactive solute transport in soil that has been used since the early 1990s as a research tool to investigate the effects of macropore flow on soil hydrology and contaminant transport under transient field conditions. It is also widely used in the form of bespoke versions in pesticide exposure and risk assessments for groundwater and surface waters, e.g., in registration procedures in the European Union (EU). Macropore flow is a highly episodic, fast, non-equilibrium process that can dominate the leaching of reactive solutes in structured soils. This has important consequences for model calibration and validation procedures. Firstly, it means that in addition to water contents and resident solute concentrations, water flows and flux concentrations measured at high time resolution are required. Secondly, it implies that more weight must be placed on flux data obtained during these important but short-lived episodic flow events if parameters controlling macropore flow are to be reliably estimated. Although the choice of approach will vary with the purpose, automatic or hybrid automatic/manual calibration procedures are generally recommended for MACRO, especially global methods that account for uncertainty within a multi-objective framework. Despite the complexity of the processes it attempts to represent, MACRO is a parsimonious model, requiring only five additional parameters to simulate water flow and reactive solute transport, compared to the use of the Richards equation and the advection-dispersion equation. Nevertheless, for practical reasons, the size of the parameter space that can be explored by calibration is often quite limited. This emphasizes the importance of measuring those parameters that can be measured and the role of sensitivity analyses in supporting the choice of parameters to calibrate. This will vary between applications, but at least for structured soils, all five parameters regulating the generation and strength of non-equilibrium flow and transport are usually rather influential.

Pesticide leaching from two Swedish topsoils of contrasting texture amended with biochar

The use of biochar as a soil amendment has recently increased because of its potential for long-term soil carbon sequestration and its potential for improving soil fertility. The objective of this study was to quantify the effects of biochar soil incorporation on pesticide adsorption and leaching for two Swedish topsoils, one clay soil and one loam soil. We used the non-reactive tracer bromide and the pesticides sulfosulfuron, isoproturon, imidacloprid, propyzamid and pyraclostrobin, substances with different mobility in soil. Adsorption was studied in batch experiments and leaching was studied in experiments using soil columns (20 cm high, 20 cm diameter) where 0.01 kg kg(-1) dw biochar powder originating from wheat residues had been mixed into the top 10 cm. After solute application the columns were exposed to simulated rain three times with a weekly interval and concentrations were measured in the effluent water. The biochar treatment resulted in significantly larger adsorption distribution coefficients (Kd) for the moderately mobile pesticides isoproturon and imidacloprid for the clay soil and for imidacloprid only for the loam soil. Relative leaching of the pesticides ranged from 0.0035% of the applied mass for pyraclostrobin (average Kd=360 cm3 g(-1)) to 5.9% for sulfosulfuron (average Kd=5.6 cm3 g(-1)). There were no significant effects of the biochar amendment on pesticide concentrations in column effluents for the loam soil. For the clay soil concentrations were significantly reduced for isoproturon, imidacloprid and propyzamid while they were significantly increased for the non-mobile fungicide pyraclostrobin suggesting that the transport was facilitated by material originating from the biochar amendment.

Simulation of pharmaceutical and personal care product transport to tile drains after biosolids application.

Pharmaceuticals and personal care products (PPCPs) carried in biosolids may reach surface waters or ground water when these materials are applied as fertilizer to agricultural land. During preferential flow conditions created by land application of liquid municipal biosolids (LMB), the residence time of solutes in the macropores may be too short for sorption equilibration. The physically based dual-permeability model MACRO is used in environmental risk assessments for pesticides and may have potential as an environmental risk assessment tool for PPCPs. The objective of this study was to evaluate MACRO and an updated version of MACRO that included non-equilibrium sorption in macropores using data from experiments conducted in eastern Ontario, Canada on the transport of three PPCPs (atenolol, carbamazepine, and triclosan), the nicotine metabolite cotinine, and the strongly sorbing dye rhodamine WT applied in LMB. Results showed that the MACRO model could not reproduce the measured rhodamine WT concentrations (Nash-Sutcliffe coefficient [NS] for the best simulation = -0.057) in drain discharge. The updated version resulted in better fits to measured data for PPCP (average NS = 0.97) and rhodamine WT (NS = 0.84) concentrations. However, it was not possible to simulate all compounds using the same set of hydraulic parameters, which indicates that the model does not fully account for all relevant processes. The results presented herein show that non-equilibrium sorption in macropores has a large impact on simulated solute transport for reactive compounds contained in LMB. This process should be considered in solute transport models that are used for environmental risk assessments for such compounds.

Spatial variation in herbicide leaching from a marine clay soil via subsurface drains.

Background Subsurface transport via tile drains can significantly contribute to pesticide contamination of surface waters. The spatial variation in subsurface leaching of normally applied herbicides was examined together with phosphorus losses in 24 experimental plots with water sampled flow-proportionally. The study site was a flat, tile-drained area with 60% marine clay in the topsoil in southeast Sweden. The objectives were to quantify the leaching of frequently used herbicides from a tile drained cracking clay soil and to evaluate the variation in leaching within the experimental area and relate this to topsoil management practices (tillage method and structure liming). Results In summer 2009, 0.14, 0.22 and 1.62%, respectively, of simultaneously applied amounts of MCPA, fluroxypyr and clopyralid were leached by heavy rain five days after spraying. In summer 2011, on average 0.70% of applied bentazone was leached by short bursts of intensive rain 12 days after application. Peak flow concentrations for 50% of the treated area for MCPA and 33% for bentazone exceeded the Swedish no-effect guideline values for aquatic ecosystems. Approximately 0.08% of the glyphosate applied was leached in dissolved form in the winters of 2008/2009 and 2010/2011. Based on measurements of glyphosate in particulate form, total glyphosate losses were twice as high (0.16%) in the second winter. The spatial inter-plot variation was large (72–115%) for all five herbicides studied, despite small variations (25%) in water discharge. Conclusions The study shows the importance of local scale soil transport properties for herbicide leaching in cracking clay soils.

Simulating Sulfadimidine Transport in Surface Runoff and Soil at the Microplot and Field Scale

To prevent residues of veterinary medicinal products (VMPs) from contaminating surface waters and ground water, an environmental impact assessment is required before a new product is allowed on the market. Physically based simulation models are advocated for the calculation of predicted environmental concentrations at higher tiers of the assessment process. However, the validation status of potentially useful models is poor for VMP transport. The objective of this study was to evaluate the dual-permeability model MACRO for simulation of transport of sulfonamide antibiotics in surface runoff and soil. Special focus was on effects of solute application in liquid manure, which may alter the hydraulic properties at the soil surface. To this end we used data from a microplot runoff experiment and a field experiment, both conducted on the same clay loam soil prone to preferential flow. Results showed that the model could accurately simulate concentrations of sulfadimidine and the nonreactive tracer bromide in runoff and in soil from the microplot experiments. The use of posterior parameter distributions from calibrations using the microplot data resulted in poor simulations for the field data of total sulfadimidine losses. The poor results may be due to surface runoff being instantly transferred off the field in the model, whereas in reality re-infiltration may occur. The effects of the manure application were reflected in smaller total and micropore hydraulic conductivities compared with the application in aqueous solution. These effects could easily be accounted for in regulatory modeling.

Effects of surfactant use and peat amendment on leaching of fungicides and nitrate from golf greens

Soil water repellency in golf putting greens may induce preferential “finger flow”, leading to enhanced leaching of surface applied agrochemicals such as fungicides and nitrate. We examined the effects of root zone composition and the use of the non-ionic surfactant Revolution on soil water repellency, soil water content distributions, infiltration rates, turf quality, and fungicide and nitrate leaching from April 2007 to April 2008. The study was made on 4-year-old experimental green seeded with creeping bentgrass (Agrostis stolonifera L.) ‘Penn A-4’ at Landvik in southeast Norway. Eight lysimeters with two different root zone materials: (i) straight sand (1% gravel, 96% sand, 3% silt and clay, and 4 g kg−1 organic matter) (SS) and (ii) straight sand mixed with Sphagnum peat to an organic matter content of 25 g kg−1 (SP) were used in this study. Surfactant treatment reduced the spatial variability of water contents, increased infiltration rates and reduced water drop penetration times (WDPTs) by on average 99% in and just below the thatch layer. These effects were most evident for SS lysimeters. Surfactant treatment resulted on average in an 80% reduction of total fungicide leaching, presumably due to reduced preferential finger flow facilitated by decreased soil water repellency. Peat amendment reduced fungicide leaching by 90%, probably due to increased sorption of the fungicides to organic matter. Nitrate leaching was also smaller from surfactant-treated straight-sand root zones, but this effect was not significant.

Spatial and temporal patterns of pesticide concentrations in streamflow, drainage and runoff in a small Swedish agricultural catchment

A better understanding of the dominant source areas and transport pathways of pesticide losses to surface water is needed for targeting mitigation efforts in a more cost-effective way. To this end, we monitored pesticides in surface water in an agricultural catchment typical of one of the main crop production regions in Sweden. Three small sub-catchments (88-242ha) were selected for water sampling based on a high-resolution digital soil map developed from proximal sensing methods and soil sampling; one sub-catchment had a high proportion of clay soils, another was dominated by coarse sandy soils while the third comprised a mix of soil types. Samples were collected from the stream, from field drains discharging into the stream and from within-field surface runoff during spring and early summer in three consecutive years. These samples were analyzed by LC-MS/MS for 99 compounds, including most of the polar and semi-polar pesticides frequently used in Swedish agriculture. Information on pesticide applications (products, doses and timing) was obtained from annual interviews with the farmers. There were clear and consistent differences in pesticide occurrence in the stream between the three sub-catchments, with both the numbers of detected compounds and concentrations being the largest in the area with a high proportion of clay soils and with very few detections in the sandy sub-catchment. Macropore flow to drains was most likely the dominant loss pathway in the studied area. Many of the compounds that were detected in drainage and stream water samples had not been applied for several years. This suggests that despite the predominant role of fast flow pathways in determining losses to the stream, long-term storage along the transport pathways also occurs, presumably in subsoil horizons where degradation is slow.

Sequential use of the STICS crop model and of the MACRO pesticide fate model to simulate pesticides leaching in cropping systems

The current challenge in sustainable agriculture is to introduce new cropping systems to reduce pesticides use in order to reduce ground and surface water contamination. However, it is difficult to carry out in situ experiments to assess the environmental impacts of pesticide use for all possible combinations of climate, crop, and soils; therefore, in silico tools are necessary. The objective of this work was to assess pesticides leaching in cropping systems coupling the performances of a crop model (STICS) and of a pesticide fate model (MACRO). STICS-MACRO has the advantage of being able to simulate pesticides fate in complex cropping systems and to consider some agricultural practices such as fertilization, mulch, or crop residues management, which cannot be accounted for with MACRO. The performance of STICS-MACRO was tested, without calibration, from measurements done in two French experimental sites with contrasted soil and climate properties. The prediction of water percolation and pesticides concentrations with STICS-MACRO was satisfactory, but it varied with the pedoclimatic context. The performance of STICS-MACRO was shown to be similar or better than that of MACRO. The improvement of the simulation of crop growth allowed better estimate of crop transpiration therefore of water balance. It also allowed better estimate of pesticide interception by the crop which was found to be crucial for the prediction of pesticides concentrations in water. STICS-MACRO is a new promising tool to improve the assessment of the environmental risks of pesticides used in cropping systems.

Surface Runoff of Pesticides from a Clay Loam Field in Sweden

Pesticides stored at or close to the soil surface after field application can be mobilized and transported off the field when surface runoff occurs. The objective of our study was to quantify the potential pesticide losses in surface runoff from a conventionally managed agricultural field in a Swedish climate. This was achieved by measuring surface runoff volumes and concentrations in runoff of six spring-applied pesticides and autumn-applied glyphosate and its metabolite aminomethylphosphonic acid (AMPA). Measurements were performed for 3 yr both during the growing seasons and during intervening winter snowmelt periods on a clay loam field close to Uppsala. During growing seasons, surface runoff was generated on only five occasions during one 25-d period in 2012 when the infiltration capacity of the soil may have been reduced by structural degradation due to large cumulative rainfall amounts after harrowing. Concentrations in surface runoff exceeded Swedish water quality standards in all samples during this growing season for diflufenican and pirimicarb. Surface runoff was generated during three snowmelt periods during the winter of 2012-2013. All of the applied pesticides were found in snowmelt samples despite incorporation of residues by autumn plowing, degradation, and leaching into the soil profile during the period between spraying and sampling. Concentrations of glyphosate ranged from 0.12 to 7.4 μg L, and concentrations of AMPA ranged from 0 to 2.7 μg L. Our results indicate that temporal changes in hydraulic properties during the growing season and when the soil freezes during winter affect pesticide losses through surface runoff.