Nitrous oxide (N2O) is a potent and persistent greenhouse gas, with rising atmospheric concentrations driven in part by inefficient use of synthetic nitrogen (N) fertilizers in agriculture. Predicting soil N2O emissions is challenging due to high spatial and temporal variability arising from complex soil biogeochemical processes. Process-based ecosystem models and standalone machine learning (ML) approaches without extensive site-specific calibration often miss high-emission episodes. Here, we show how an Ensemble Modeling System (EMS) based on outputs from an ensemble of ecosystem models coupled to an ensemble of ML models can improve predictions and understanding of N2O fluxes from US cropland. Trained and validated on ~12,000 N2O chamber measurements at 17 US Midwest sites (six crops, 35 management practices), the EMS accurately predicted daily fluxes of N2O at both training (R2= 0.84, RMSE = 16.4 g N ha−1d−1) and held-out testing sites (R2= 0.84, RMSE = 6.2 g N ha−1d−1). Analyses identified six dominant N2O drivers: soil organic carbon (SOC), NH4+, NO3-, water-filled pore space, temperature, and aboveground biomass production. Wet, warm soils produced large N2O peaks only with sufficient SOC and mineral N; in low-SOC soils, fluxes remained low. Incorporating these drivers into process-based models might significantly improve their predictive capacity. The EMS demonstrates a strong potential to predict N2O fluxes at unseen sites, enabling more reliable regional inventories, improved gap-filling where measurements are sparse, and enhanced understanding of mechanisms to advance targeted mitigation strategies in food, feed, and bioenergy crops.
more »
« less
Data from: Coupled machine learning-ecosystem ensemble models substantially improve predictions of nitrous oxide (N2O) fluxes from US croplands
Nitrous oxide (N₂O) is a potent and persistent greenhouse gas, with rising atmospheric concentrations driven in part by inefficient use of synthetic nitrogen (N) fertilizers in agriculture. Predicting soil N₂O emissions is challenging due to high spatial and temporal variability arising from complex soil biogeochemical processes. Process-based ecosystem models and standalone machine learning (ML) approaches without extensive site-specific calibration often miss high emission episodes. Here, we show how an Ensemble Modeling System (EMS) based on outputs from an ensemble of ecosystem models coupled to an ensemble of ML models can improve predictions and understanding of N2O fluxes from US cropland. Trained and validated on approximately 12,000 N2O chamber measurements at 17 U.S. Midwest sites (six crops, 35 management practices), the EMS accurately predicted daily fluxes of N2O at both training (R² = 0.84, RMSE = 16.4 g N ha⁻¹ d⁻¹) and held-out testing sites (R² = 0.84, RMSE = 6.2 g N ha⁻¹ d⁻¹). Analyses identified six dominant N₂O drivers: soil organic carbon (SOC), NH₄⁺, NO₃⁻, water-filled pore space (WFPS), soil temperature, and biomass production. Wet, warm soils produced large N₂O peaks only with sufficient SOC and mineral N; in low-SOC soils, fluxes remained low. Incorporating these drivers into process-based models might significantly improve their predictive capacity. The EMS demonstrates a strong potential to predict N₂O fluxes at unseen sites, enabling more reliable regional inventories, improved gap-filling where measurements are sparse, and enhanced understanding of mechanisms to advance targeted mitigation strategies in food, feed, and bioenergy crops. # Data from: Coupled machine learning-ecosystem ensemble models substantially improve predictions of nitrous oxide (N~2~O) fluxes from US croplands Dataset DOI: [10.5061/dryad.pvmcvdnzx](10.5061/dryad.pvmcvdnzx) ## Description of the data and file structure We present here the data that were used for the analysis presented in: Coupled machine learning-ecosystem ensemble models substantially improve predictions of nitrous oxide (N~2~O) fluxes from US croplands. ### Files and variables Files: Dataset_S1.xlsx, Dataset_S2.csv, Dataset_S3.csv, #### **Description:** **Description of data sheets** **Dataset S1A columns** * **Site_ID:** Numeric identifier for the experimental site. * **Treatment_ID:** Numeric code for the management treatment applied at that site * **DataUse:** To assign data to be used for model training (development) and testing (held-out evaluation) * **State/Province:** State acronym * **Latitude decimal deg:** Site location latitude * **Longitude decimal deg:** Site location longitude * **Start Data Year:** Starting year of data used * **End Data Year:** Ending year of data used * **Cover crop:** Type of cover crop used within the treatment * **Rotation Descriptor:** Describe the rotation of crops within the treatment * **Tillage Descriptor:** Describe tillage type within the treatment * **Residual Removal:** Describe residual management within the treatment * **Irrigation:** Describe if irrigation was applied or not within the treatment * **N Treatment Descriptor:** Describe nitrogen amendments within the treatment * **Reference:** Reference for the data **Dataset S1B**: This sheet contains the reference list for the data used **Dataset S2 columns** * **Date:** Gas sampling days * **Site_ID:** Numeric identifier for the experimental site. * **Treatment_ID:** Numeric code for the management treatment applied at that site * **DataUse:** To assign data to be used for model training (development) and testing (held-out evaluation) * **Observed N2O:** Daily average N2O flux measured (g N2O-N ha-1d-1) * **Predicted N2O:** Daily average N2O flux predicted by multimodel hybrid framework (g N2O-N ha-1d-1) * **NH4:** Process-based models simulated daily NH4-N content in the top 30-cm soil layer (kg ha-1) * **SOC:** Process-based models simulated daily soil organic carbon in the top 30-cm soil layer (kg ha-1) * **NO3:** Process-based models simulated daily NO3-N content in the top 30-cm soil layer (kg ha-1) * **ST:** Process-based models simulated daily average soil temperature in the top 30-cm soil layer (°C) * **WFPS:** Process-based models simulated daily water-filled pore space in the top 30-cm soil layer (fraction) * **ABG:** Process-based models simulated daily above-ground biomass (kg ha-1) * **BG:** Process-based models simulated daily below-ground biomass (kg ha-1) * **SRAD:** Average solar radiation for the last five days before gas sampling (Watt m-2) * **Tmax:** Average maximum air temperature for the last three days before gas sampling (°C) * **APrecip:** Average precipitation in the last fifteen days before gas sampling (mm) * **Wspd:** Average wind in the last fifteen days before gas sampling (m s-1) * **LAI:** Process-based models simulated daily leaf area index (m2 m2) * **Nstress:** Process-based models simulated the daily nitrogen stress factor (fraction) * **Wstress:** Process-based models simulated the daily water stress factor (fraction) * **PET:** Process-based models simulated daily potential evapotranspiration (mm) * **SE:** Process-based models simulated daily soil evaporation (mm) * **SPrecip:** Cumulative precipitation in the last two days before gas sampling (mm) * **SH:** Average specific humidity in the last three days before gas sampling (g kg-1) * **RH:** Average relative humidity in the last fifteen days before gas sampling (%) **Dataset S3 columns** * **Date:** Gas sampling days * **Site_ID:** Numeric identifier for the experimental site. * **Treatment_ID:** Numeric code for the management treatment applied at that site * **DataUse:** To assign data to be used for model training (development) and testing (held-out evaluation) * **SD:** Monte Carlo standard deviation of the simulated daily N₂O flux distribution (g N2O-N ha-1d-1) * **CV**: Monte Carlo coefficient of variation of the simulated daily N₂O flux distribution (%) * **CI05:** 5th percentile (lower 90 % confidence bound) of the Monte Carlo flux distribution(g N2O-N ha-1d-1) * **CI95:** 95th percentile (upper 90 % confidence bound) of the Monte Carlo flux distribution(g N2O-N ha-1d-1)
more »
« less
- Award ID(s):
- 2224712
- PAR ID:
- 10670465
- Publisher / Repository:
- Dryad
- Date Published:
- Edition / Version:
- 6
- Subject(s) / Keyword(s):
- FOS: Earth and related environmental sciences FOS: Earth and related environmental sciences FOS: Agricultural sciences FOS: Agricultural sciences nitrous oxide emission ensemble modeling system multi-model ensemble
- Format(s):
- Medium: X Size: 5237226 bytes
- Size(s):
- 5237226 bytes
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Soil nitrous oxide (N2O) emissions exhibit high variability in intensively managed cropping systems, which challenges our ability to understand their complex interactions with controlling factors. We leveraged 17-years (2003-2019) of measurements at the Kellogg Biological Station LTER/LTAR site to better understand controls of N2O emissions in four corn–soybean–winter wheat rotations employing Conventional, No-till, Reduced input, and Biologically-based/organic inputs. We used a Random Forest machine learning model to predict daily N2O fluxes, trained separately for each system with 70% of observations, using variables such as crop species, daily air temperature, cumulative 2-day precipitation, water-filled pore space, and soil nitrate and ammonium concentrations. The model explained 29 to 42% of daily N2O flux variability in test data, with greater predictability for the corn phase in each system. The long-term rotations showed different controlling factors and threshold conditions influencing N2O emissions. In the Conventional system, the model identified ammonium (>15 kg N ha-1) and daily temperature (>23 °C) as the most influential variables; in the No-till system, climate variables, precipitation, and temperature were important variables. In low input and organic systems, where red clover (Trifolium repens L.; before corn) and cereal rye (Secale cereale L.; before soybean) cover crops were integrated, nitrate was the predominant variable, followed by precipitation and temperature. In low input and biologically-based systems, red clover residues increased soil nitrogen availability to influence N2O emissions. Long-term data facilitated machine learning for predicting N2O emissions in response to differential controls and threshold responses to management, environmental, and biogeochemical drivers.more » « less
-
Abstract. Nitrogen (N) fertilizer inputs to agricultural soils area leading cause of nitrous oxide (N2O) emissions. Legume cover cropsare an alternative N source that can reduce agricultural N2O emissionscompared to fertilizer N. However, our understanding of episodic N2Oflux following cover crop incorporation by tillage is limited and hasfocused on single-species cover crops. Our study explores whether increasingcover crop functional diversity with a legume–grass mixture can reduce pulseemissions of N2O following tillage. In a field experiment, we plantedcrimson clover (Trifolium incarnatum L.), cereal rye (Secale cereal L.), a clover–rye mixture, and a no-covercontrol at two field sites with contrasting soil fertility properties inMichigan. We hypothesized that N2O flux following tillage of the covercrops would be lower in the mixture and rye compared to the clovertreatment because rye litter can decrease N mineralization rates. Wemeasured N2O for approximately 2 weeks following tillage to capturethe first peak in N2O emissions in each site. Across cover croptreatments, the higher-fertility site, CF, had greater cover crop biomass,2-fold-higher aboveground biomass N, and higher cumulative N2Oemissions than the lower-fertility site, KBS (413.4±67.5 vs. 230.8±42.5 g N2O-N ha−1; P=0.004). Therewas a significant treatment effect on daily emissions at both sites. AtCF, N2O fluxes were higher following clover than the control 6 d aftertillage. At KBS, fluxes from the mixture were higher than rye 8 and 11 dafter tillage. When controlling for soil fertility differences betweensites, clover and mixture led to approximately 2-fold-higher N2Oemissions compared to rye and fallow treatments. We found partial supportfor our hypothesis that N2O would be lower following incorporation ofthe mixture than clover. However, treatment patterns differed by site,suggesting that interactions between cover crop functional types andbackground soil fertility influence N2O emissions during cover cropdecomposition.more » « less
-
Abstract Soil nitrous oxide (N2O) emissions exhibit high variability in intensively managed cropping systems, which challenges our ability to understand their complex interactions with controlling factors. We leveraged 17 years (2003–2019) of measurements at the Kellogg Biological Station Long‐Term Ecological Research (LTER)/Long‐Term Agroecosystem Research (LTAR) site to better understand the controls of N2O emissions in four corn–soybean–winter wheat rotations employing conventional, no‐till, reduced input, and biologically based/organic inputs. We used a random forest machine learning model to predict daily N2O fluxes, trained separately for each system with 70% of observations, using variables such as crop species, daily air temperature, cumulative 2‐day precipitation, water‐filled pore space, and soil nitrate and ammonium concentrations. The model explained 29%–42% of daily N2O flux variability in the test data, with greater predictability for the corn phase in each system. The long‐term rotations showed different controlling factors and threshold conditions influencing N2O emissions. In the conventional system, the model identified ammonium (>15 kg N ha−1) and daily air temperature (>23°C) as the most influential variables; in the no‐till system, climate variables such as precipitation and air temperature were important variables. In low‐input and organic systems, where red clover (Trifolium repensL.; before corn) and cereal rye (Secale cerealeL.; before soybean) cover crops were integrated, nitrate was the predominant predictor of N2O emissions, followed by precipitation and air temperature. In low‐input and biologically based systems, red clover residues increased soil nitrogen availability to influence N2O emissions. Long‐term data facilitated machine learning for predicting N2O emissions in response to differential controls and threshold responses to management, environmental, and biogeochemical drivers.more » « less
-
Climate change and nitrogen (N) application significantly influence agricultural productivity and soil greenhouse gas emissions. However, the interactive effects of interannual climate variability and N fertilization legacy on corn yield and soil nitrous oxide (N2O) emissions remain inadequately understood. In this study, we employed the DeNitrification-DeComposition (DNDC) model to simulate corn yield and soil N2O emissions over a 40-year period (1981~2020). We designed a series of experiments by adjusting climate year data to quantify interannual variability in corn yield and soil N2O emissions, while also disentangling the contributions of climate variability and N legacy effects. The results revealed substantial interannual variability in both corn yield and soil N2O emissions. Corn yield was primarily driven by changes in growing season precipitation, while soil N2O emissions were influenced by precipitation, exchangeable ammonium N (NH4+), and nitrification-denitrification processes. Severe drought strongly reduced corn yield, while soil N2O emissions exhibited a gradual yet pronounced legacy effect of N application, increasing from 1.69 to 7.85 kg N ha-1 over the 40-year period. This study highlights the relatively weak influence of interannual climate variability compared to the stronger legacy effects of N application on crop yield and soil N2O emissions, providing valuable insights for sustainable agricultural and environmental management.more » « less
An official website of the United States government
