Modeling & SimulationCompleted2026

Mechanistic Modeling of AMPA Persistence Across Soil Types

Meta-analysis using TESFO mechanistic modeling to evaluate glyphosate metabolite (AMPA) behavior across four soil types. Identified critical mobility-persistence trade-off requiring site-specific agricultural risk management.

Mechanistic Modeling of AMPA Persistence Across Soil Types

Gallery

Problem

Glyphosate is the world's most widely used herbicide, but its primary metabolite AMPA (aminomethylphosphonic acid) is highly persistent and potentially harmful to soil ecosystems. Regulatory agencies apply universal application limits without accounting for soil-specific environmental behavior. This one-size-fits-all approach ignores fundamental soil physics: sandy soils allow rapid leaching to groundwater, while clay and humus-rich soils retain AMPA for extended periods. Without mechanistic understanding of how soil composition affects AMPA mobility and degradation, risk assessment remains incomplete.

Approach

Conducted systematic meta-analysis using the TESFO (Two-stage Equilibrium Sorption with First-Order degradation) mechanistic model to simulate AMPA environmental fate across four representative soil types: Sand, Loam/Reference, Clay, and Humus-rich. Parameterized each soil scenario with literature-derived physical properties: bulk density (1.2-1.6 g/cm³), volumetric water content (0.15-0.45 cm³/cm³), sorption capacity (Freundlich and Langmuir isotherms with Kf = 0.5-8.2 L/kg), and first-order degradation kinetics (k = 0.0069-0.0231 d⁻¹). Simulated 200-day temporal dynamics of AMPA concentration in soil solution following standard glyphosate application. Analyzed peak concentration timing, dissipation curves, and calculated DT90 (time to 90% dissipation) for leaching risk assessment.

Data

Mechanistic simulation outputs for each soil type: temporal concentration profiles (0-200 days), peak AMPA levels (mg/L), time-to-peak (days), and leaching potential indices. Reference scenario based on University of Bonn experimental loam soil. Soil parameters sourced from published pedotransfer functions and glyphosate degradation literature (Vereecken et al. 2005, Wimmer et al. 2022). All simulations assume identical initial glyphosate application rate to isolate soil-specific effects.

Validation

Model outputs validated against published field studies reporting AMPA detection in agricultural soils across Europe. DT90 values (>100 days across all soil types) matched long-term monitoring data showing AMPA persistence in soil profiles 6-12 months post-application. Sorption coefficient ranges verified against laboratory batch equilibrium experiments from multiple independent studies. Sensitivity analysis confirmed that sorption capacity (Kf parameter) is the dominant control on leaching vs. retention behavior—higher than degradation rate or water content variability.

Results

Identified strict mobility-persistence trade-off governed by soil texture and organic matter: (1) Sand—rapid peak concentration (10-15 mg/L at 5-10 days) followed by fast dissipation, HIGH leaching risk to groundwater due to low sorption capacity (Kf ≈ 0.5-1.2), short retention time but mobile plume; (2) Humus-rich—strong sorption (Kf ≈ 6-8.2) dramatically reduces mobility, LOW leaching risk, but HIGH long-term persistence with extended tail in concentration profile; (3) Clay and Loam—intermediate behavior balancing moderate retention with gradual release. Critical finding: AMPA DT90 exceeded 100 days in ALL modeled scenarios, confirming metabolite persistence dominates system dynamics regardless of soil type. This persistence asymmetry (AMPA >> glyphosate parent compound) drives long-term environmental accumulation risk.

My Role

Environmental Modeler. Executed comprehensive meta-analysis across four soil scenarios, parameterized TESFO mechanistic model with literature-derived sorption isotherms and degradation kinetics, analyzed temporal dynamics to quantify leaching risk vs. persistence trade-offs, and synthesized findings into soil-specific risk profiles for regulatory context.

Next Steps

Extend analysis to include rainfall scenarios (infiltration-driven leaching pulses) and groundwater transport modeling. Validate predictions with field lysimeter data from long-term monitoring sites. Incorporate spatial heterogeneity using 3D reactive transport models (e.g., HYDRUS, MIN3P) for landscape-scale risk mapping. Propose soil-texture-based application rate adjustments for precision agriculture implementation.

Key Outcomes

  • Quantified mobility-persistence trade-off across 4 soil types
  • Sand: HIGH leaching risk (low Kf, rapid peak)
  • Humus-rich: HIGH persistence (strong sorption, slow release)
  • AMPA DT90 > 100 days in ALL soils (metabolite dominates)
  • Demonstrated inadequacy of universal regulatory limits

Tech Stack

PythonTESFO ModelFreundlich IsothermLangmuir IsothermNumerical Simulation

Tags

mechanistic-modelingenvironmental-sciencesoil-physicsrisk-assessmentpython