Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Biogeochemical modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

Responses of Boreal Forest Ecosystems and Permafrost to Climate Change and Disturbances: A Modeling Perspective

The north circumpolar region contains a large amount of carbon. This carbon storage is vulnerable due to permafrost degradation and wildfire disturbances under ongoing and projected climate change. Climate warming and wildfires change soil organic horizons gradually or abruptly, and modify permafrost thermal-hydrology and biogeochemistry, ecosystem structures, functions, and capability of sequestrating rising atmospheric CO2. Land models do not fully take accounts of these interactions and its complexity in the high latitude. This chapter describes a terrestrial ecosystem model with dynamic organic soil module (DOS-TEM) and its unique freezing-thawing algorithm, and presents key results of its applications mainly in boreal forests of Alaska. The DOS-TEM explicitly considers interactions of soil thermal and hydrological processes, permafrost degradation and the direct and indirect effects of wildfire disturbances, in addition to soil–plant C and N cycles. We first introduce four modules of DOS-TEM, focusing on its disturbance module and coupling with a dynamic organic soil module. Then we describe and validate DOS-TEM’s freezing-thawing algorithm and development based on two-directional Stefan algorithm (TDSA). Finally, we apply the DOS-TEM at site and region scales, with a focus on model ability to dynamically simulate soil organic thickness under warming and wildfires, and consequent impacts on permafrost in the Yukon River Basin. We conclude that land surface model development is urgently needed to include other critical landscape processes, such as thermalkarst and other disturbances, to synchronize thermal-hydrological-biogeochemical processes, and to incorporate an advanced understanding of biospheric feedbacks to atmosphere and ecosystems. Such a complexity of modeling scope is plausible with advancement of high performance computing.

Yi, Shuhua↗

Can Biogeochemical Tracer Observations Constrain Southern Ocean Diapycnal Mixing Rates?

Abstract Direct observations of background diapycnal mixing rates in the Southern Ocean (SO) are limited spatially and temporally, making the choice of an appropriate value to parameterize this mixing in Earth system models a challenge. However, the deployment of Argo floats throughout the SO has provided an extensive range of observations of both physical and biogeochemical parameters. We use an ocean state estimate run with various background diapycnal mixing coefficients to assess if biogeochemical tracer observations can be used to better constrain SO diapycnal mixing rates. We find that vertical tracer distributions in the SO are highly sensitive to the rate of background diapycnal mixing and can provide an upper limit on background mixing rates. This demonstrates the importance of biogeochemical tracer observations throughout the full depth of the water column to validate ocean models.

Ellison, Elizabeth [The Australian National Univer↗

Seasonal variability of light availability and utilization in the Sargasso Sea

A 2 year time series of optical, biogeochemical, and physical parameters, taken near the island of Bermuda, is used to evaluate the sources of temporal variability in light avaliability and utilization in the Sargasso Sea. Integrated assessments of light availability are made by examining the depth of constant percent incident photosynthetically available radiation (% PAR) isolumes. To first order, changes in the depth %PAR isolumes were caused by physical processes: deep convection mixing in the winter which led to the spring bloom and concurrent shallowing of %PAR depths and the occurrence of anomalous thermohaline water masses during the summer and fall seasons. Spectral light availability variations are assessed using determinations of diffuse attenuation coefficient spectra which illustrates a significant seasonal cycle in colored detrital particulate and/or dissolved materials that is unrelated to changes in chlorophyll pigment concentrations. Temporal variations in the photosynthetic light utilization index Psi are used to assess vertically intergrated light utilization variations. Values of Psi are highly variable and show no apparent seasonal pattern which indicates that Psi is not simply a 'biogeochemical constant.' Determinations of in situ primary production rates and daily mean PAR fluxes are used to diagnose the relative role of light limitation in determining vertically integrated rates of primary production integral PP. The mean depth of the light-saturated zone (the vertical region where the daily mean PAR flux was greater than or equal to the saturation irradiance) is only approximately 40 m, although more than one half of interal PP occurred within this zone. Production model results illustrate that accurate predictions of integral PP are dependent upon rates of light-saturated production rather than upon indices of light limitation. It seems unlikely that significant improvements in simple primary production models will come from the partitioning of the Earth's seas into biogeochemical provinces.

Siegel, David A.↗

Simulating Marine Ecosystem Dynamics and Biogeochemical Cycling With Multiple Plankton Functional Types

Current representations of marine ecosystems in Earth System Models are greatly simplified, neglecting key interactions between dynamic food webs, biogeochemistry, and climate change. We use the Marine Biogeochemistry Library code base within the Community Earth System Model 2.2.2 to create an expanded ecosystem model with eight phytoplankton groups and four zooplankton size classes (MARBL-8P4Z). Incorporating more specific plankton types and size classes has the potential to capture a wider range of possible behaviors of the ecosystem, its complex interactions with biogeochemistry, and its feedback to climate change. It also permits stronger observational constraints, including in situ group-specific biomass and various observational estimates of plankton community composition. MARBL-8P4Z broadly captures observed global-scale patterns in biomass and community composition for both phytoplankton and zooplankton, with a good performance in simulating broad biogeochemistry fields. The model shows comparable spatial patterns and magnitudes to the observed picophytoplankton biomass (Prochlorococcus, Synechococcus, picoeukaryotes), and captures the seasonal cycle of mesozooplankton biomass. Picophytoplankton groups and microzooplankton dominate biomass and production in oligotrophic, subtropical regions, while nano-phytoplankton, diatoms and the larger zooplankton groups prevail at higher latitudes and within upwelling zones. The model simulates reasonable energy transfer efficiency through the food web, with tight linkages between the phytoplankton community composition, zooplankton grazing, and carbon export, with the potential to link to fisheries models. Thus, MARBL-8P4Z has the potential to account for key climate-driven ecological shifts in the plankton that will modify ocean biogeochemistry in the future.

54 ENVIRONMENTAL SCIENCES↗

Recent Progresses in Incorporating Human Land-Water Management into Global Land Surface Models Toward Their Integration into Earth System Models

The global water cycle has been profoundly affected by human land-water management. As the changes in the water cycle on land can affect the functioning of a wide range of biophysical and biogeochemical processes of the Earth system, it is essential to represent human land-water management in Earth system models (ESMs). During the recent past, noteworthy progress has been made in large-scale modeling of human impacts on the water cycle but sufficient advancements have not yet been made in integrating the newly developed schemes into ESMs. This study reviews the progresses made in incorporating human factors in large-scale hydrological models and their integration into ESMs. The study focuses primarily on the recent advancements and existing challenges in incorporating human impacts in global land surface models (LSMs) as a way forward to the development of ESMs with humans as integral components, but a brief review of global hydrological models (GHMs) is also provided. The study begins with the general overview of human impacts on the water cycle. Then, the algorithms currently employed to represent irrigation, reservoir operation, and groundwater pumping are discussed. Next, methodological deficiencies in current modeling approaches and existing challenges are identified. Furthermore, light is shed on the sources of uncertainties associated with model parameterizations, grid resolution, and datasets used for forcing and validation. Finally, representing human land-water management in LSMs is highlighted as an important research direction toward developing integrated models using ESM frameworks for the holistic study of human-water interactions within the Earths system.

biogeochemistry↗

Long-term leaf C:N ratio change under elevated CO 2 and nitrogen deposition in China: Evidence from observations and process-based modeling

Climate change, elevating atmosphere CO 2 (eCO 2 ) and increased nitrogen deposition (iNDEP) are altering the biogeochemical interactions between plants, microbes and soils, which further modify plant leaf carbon-nitrogen (C:N) stoichiometry and their carbon assimilation capability. Many field experiments have observed large sensitivity of leaf C:N ratio to eCO 2 and iNDEP. However, the large-scale pattern of this sensitivity is still unclear, because eCO 2 and iNDEP drive leaf C:N ratio toward opposite directions, which are further compounded by the complex processes of nitrogen acquisition and plant-and-microbial nitrogen competition. Here, we attempt to map the leaf C:N ratio spatial variation in the past 5 decades in China with a combination of data-driven model and process-based modeling. These two approaches showed consistent results. Over different regions, we found that leaf C:N ratio had significant but uneven changes between 2 time periods (1960-1989 and 1990-2015): a 5% ± 8% increase for temperate grasslands in northern China, a 3% ± 6% increase for boreal grasslands in western China, and by contrast, a 7% ± 6% decrease for temperate forests in southern China, and a 3% ± 5% decrease for boreal forests in northeastern China. Additionally, the structural equation models indicated that the leaf C:N change was sensitive to ΔNDEP, ΔCO 2 and ΔMAT rather than ΔMAP and ecosystem types. In this work, process-based modeling suggested that iNDEP was the main source of soil mineral nitrogen change, dominating leaf C:N ratio change in most areas in China, while eCO 2 led to leaf C:N ratio increase in low iNDEP area. This study also indicates that the long-term leaf C:N ratio acclimation was dominated by climate constraint, especially temperature, but was constrained by soil N availability over decade scale.

54 ENVIRONMENTAL SCIENCES↗

The JGOFS North Atlantic Bloom Experiment: An overview

The North Atlantic Bloom Experiment (NABE) of JGOFS presents a unique opportunity and challenge to the data management community because of the diversity and large size of biogeochemical data sets collected. NABE was a pilot study for JGOFS and has also served as a pilot study within the U.S. NODC for management and archiving of the data sets. Here I present an overview to some of the scientific results of NABE, which will be published as an Introduction to a special volume of NABE results in Deep-Sea Research later this year. An overview of NABE data management is given elsewhere in the present report. This is the first collection of papers from the Joint Global Ocean Flux Study (JGOFS). Formed as an international program in 1987, JGOFS has four principal elements: modelling and data management, multidisciplinary regional process studies, a global survey of biogeochemical properties and long-term time series observatories. In 1989-1990 JGOFS conducted a pilot process study of the spring phytoplankton bloom, the North Atlantic Bloom Experiment (NABE). JGOFS decided to conduct a large scale, internationally-coordinated pilot study in the North Atlantic because of its proximity to the founding nations of the project, the size and predictability of the bloom and its fundamental impact on ocean bio-geochemistry (Billett et al., 1983; Watson and Whitfield, 1985; Pfannkuche, 1992). In 1989, six research vessels from Canada, Germany, The Netherlands, the United Kingdom and the USA and over 200 scientists and students from more than a dozen nations participated in NABE. Some of their initial results are reported in this volume.

Ducklow, Hugh W.↗

Multiplexed Knockouts in the Model Diatom Phaeodactylum by Episomal Delivery of a Selectable Cas9

Marine diatoms are eukaryotic microalgae that play significant ecological and biogeochemical roles in oceans. They also have significant potential as organismal platforms for exploitation to address biotechnological and industrial goals. In order to address both modes of research, sophisticated molecular and genetic tools are required. We presented here new and improved methodologies for introducing CRISPRCas9 to the model diatom Phaeodactylum tricornutum cells and a streamlined protocol for genotyping mutant cell lines with previously unknown phenotypes. First, bacterialconjugation was optimized for the delivery of Cas9 by transcriptionally fusing Cas9 to a selectable marker by the 2A peptide. An episome cloning strategy using both negative and positive selection was developed to streamline CRISPR-episome assembly. Next, cell line picking and genotyping strategies, that utilize manual sequencing curation, TIDE sequencing analysis, and a T7 endonuclease assay, were developed to shorten the time required to generate mutants. Following this new experimental pipeline, both singlegene and two-gene knockout cell lines were generated at mutagenesis efficiencies of 48% and 25%, respectively. Lastly, a protocol for precise gene insertions via CRISPRCas9 targeting was developed using particle-bombardment transformation methods. Overall, the novel Cas9 episome design and improved genotyping methods presented here allow for quick and easy genotyping and isolation of Phaeodactylum mutant cell lines (less than 3 weeks) without relying on a known phenotype to screen for mutants.

59 BASIC BIOLOGICAL SCIENCES↗

Corrinoids as model nutrients to probe microbial interactions in a soil ecosystem

Earth’s soils are habitats for microbial communities that drive biogeochemical cycling, plant growth, and carbon storage and persistence. The thousands of microbial species living in soil form an intricate web of interactions involving the exchange of molecules produced by different microbes. Understanding in detail how these molecular exchanges occur and how they shape microbial communities may lead to new methods to improve soil health, bioremediation efforts, and better understanding of biogeochemical processes. The overall goal of this research is to gain a deeper knowledge of the microbial interactions that drive soil community structure. However, the high functional and genomic diversity in soil microbiomes has posed a challenge for current microbiology methods to achieve this goal. This research leverages a model group of key metabolites related to cobalamin (vitamin B 12 ), known as corrinoids, to investigate microbial interactions.

59 BASIC BIOLOGICAL SCIENCES↗

Interactive Gas Chemistry for Enhanced Science Capabilities of the Energy Exascale Earth System Model Version 3

Atmospheric chemistry plays a crucial role in Earth system models (ESMs), controlling atmospheric composition and radiative balance; it is highly interactive with the physical climate, biogeochemical cycles, and human systems. However, it often imposes computational challenges in an ESM. Here we develop a full troposphere‐stratosphere interactive chemistry module for the US Department of Energy's Energy Exascale Earth System Model (E3SM). We intentionally build a streamlined module based on E3SM version 2 that interacts with other components and maintains all of major chemical and chemistry‐climate feedbacks. The module incorporates a new, highly efficient tracer advection scheme; linearization of stratospheric chemistry; and abridged tropospheric chemical mechanism with 28 reactive tracers. This new model, E3SM‐chem, can readily perform century‐long climate simulations of ozone, methane, and nitrous oxide based on emission scenarios as well as provide hourly budgets for the gas‐phase radicals that drive aerosol chemistry. We evaluate E3SM‐chem with an atmosphere‐only simulation as in the recent climate model intercomparison project (CMIP6) finding results similar to the other CMIP6 models. For the present‐day, E3SM‐chem matches the standard measurement metrics for stratospheric and tropospheric ozone, surface air quality, other key reactive gases like carbon monoxide, and the methane lifetime. Overall, E3SM‐chem maintains the climate fidelity of the baseline model while adding at most 20% to the computational cost of the atmosphere model. Hence, interactive chemistry can be a default configuration for long climate simulations at resolutions of 1° or finer, which is crucial for producing self‐consistent chemistry‐climate feedbacks that alter the climate system.

54 ENVIRONMENTAL SCIENCES↗

Optimal interpolation of global dissolved oxygen: 1965–2015

Oxygen inventory of the global ocean has declined in recent decades potentially due to the warming-induced reduction in solubility as well as the circulation and biogeochemical changes associated with ocean warming and increasing stratification. Earth System Models predict continued oxygen decline for this century with profound impacts on marine ecosystem and fisheries. Observational constraint on the rate of oxygen loss is crucial for assessing the ability of models to accurately simulate these changes. There are only a few observational assessments of the global oceanic oxygen inventory reporting a range of oxygen loss. This study develops a gridded data set of dissolved oxygen for the global oceans using optimal interpolation method. The resulting gridded product includes full-depth map of dissolved oxygen as 5-year moving average from 1965 to 2015 with uncertainty estimates. The uncertainty can come from unresolved small-scale and high-frequency variability and mapping errors. The multi-decadal trend of global dissolved oxygen is in the range of -281 to -373 Tmol/decade. This estimate is more conservative than previous works. In this study, the grid points far from the observations are essentially set equal to zero anomaly from the climatology. Calculating global inventory with this approach produces a relatively conservative estimate; thus, the results from this study likely provide a useful lower bound estimate of the global oxygen loss.

54 ENVIRONMENTAL SCIENCES↗

Development of a Land Use Mapping and Monitoring Protocol for the High Plains Region: A Multitemporal Remote Sensing Application

The purpose of this project was to develop a practical protocol that employs multitemporal remotely sensed imagery, integrated with environmental parameters to model and monitor agricultural and natural resources in the High Plains Region of the United States. The value of this project would be extended throughout the region via workshops targeted at carefully selected audiences and designed to transfer remote sensing technology and the methods and applications developed. Implementation of such a protocol using remotely sensed satellite imagery is critical for addressing many issues of regional importance, including: (1) Prediction of rural land use/land cover (LULC) categories within a region; (2) Use of rural LULC maps for successive years to monitor change; (3) Crop types derived from LULC maps as important inputs to water consumption models; (4) Early prediction of crop yields; (5) Multi-date maps of crop types to monitor patterns related to crop change; (6) Knowledge of crop types to monitor condition and improve prediction of crop yield; (7) More precise models of crop types and conditions to improve agricultural economic forecasts; (8;) Prediction of biomass for estimating vegetation production, soil protection from erosion forces, nonpoint source pollution, wildlife habitat quality and other related factors; (9) Crop type and condition information to more accurately predict production of biogeochemicals such as CO2, CH4, and other greenhouse gases that are inputs to global climate models; (10) Provide information regarding limiting factors (i.e., economic constraints of pumping, fertilizing, etc.) used in conjunction with other factors, such as changes in climate for predicting changes in rural LULC; (11) Accurate prediction of rural LULC used to assess the effectiveness of government programs such as the U.S. Soil Conservation Service (SCS) Conservation Reserve Program; and (12) Prediction of water demand based on rural LULC that can be related to rates of draw-down of underground water supplies.

Price, Kevin P.↗

Recent advances in integrated hydrologic models: Integration of new domains

Over the past several decades, hydrologic models have advanced from independent models of the surface and subsurface to integrated models that can capture the terrestrial hydrologic cycle within one framework. In recent years, these coupled frameworks have seen the inclusion of biogeochemical processes, ecohydrology, sedimentation and erosion, cold region hydrology, anthropogenic activities, and atmospheric processes. This expansion is the result of increased computational, data, and modeling capabilities and capacities, as well as improved understanding of the processes that drive these integrated systems. Here, in this study, we review these recent advances to integrate new processes and systems into existing terrestrial hydrologic models and highlight the significant challenges and opportunities that remain. We identify that with so many models currently available and in development, selecting the most appropriate model is difficult, and we suggest a path for new or novice modelers to find the most appropriate code based on their needs. In addition, data required to parameterize and calibrate these models can often constrain their applicability and usefulness. However, advances in environmental sensors and measurement technology, in addition to data assimilation of non-traditional data (e.g. remote sensing, qualitative data) are providing new ways of addressing this issue. As we expand hydrologic models to integrate more processes and systems, our computational demands also increase. Recent and emerging advances in computational platforms, including cloud and quantum computing, in addition to the use of machine learning to capture some processes, will continue to support the use of increasingly larger and more complex, process-based models. Finally, we highlight that it is critical to develop state-of-the-science models that are accessible to all model users, not just those applied for research and development. We encourage continued development of diverse modeling platforms, considering the user needs, data availability, and computational resources.

54 ENVIRONMENTAL SCIENCES↗

Representing Global Soil Erosion and Sediment Flux in Earth System Models

Abstract Soil erosion produces enormous amounts of sediment, carbon, and nutrient fluxes from land to rivers, thus playing crucial roles in global biogeochemical cycles and food security. To predict soil erosion in the context of climate and land use changes, we explicitly parameterize cropland management actions (i.e., irrigation, conserved agriculture, and crop residue management) and geological factors (i.e., lithology and glacier) in the Energy Exascale Earth System Model (E3SM) soil erosion module. The erosion model is calibrated using a global‐scale regionalized parameter calibration method. The spatial variabilities of the modeled and the Revised Universal Soil Loss Equation (RUSLE) based benchmark soil erosion are consistent across vegetation, climate and soil properties. Compared with independent data, our model shows a bias reduction in 59% of the observations relative to the RUSLE‐based soil erosion, with 53% of the bias reduction exceeding 50%. This improvement is mainly due to a better representation of the topographic effect on soil erosion. Our results indicate that conserved agriculture practices have effectively reduced soil erosion in cropland by over 25% in the United States and Argentina. In contrast, irrigation has increased soil erosion in many Asian countries. For upland sediment flux, our model is consistent with the WBMsed benchmark data in inter‐basin variability but could be more skillful in simulating intra‐basin variability because it couples soil erosion and sediment flux explicitly. The developed model provides useful skills for more realistic predictions of soil erosion and river sediment dynamics under environmental changes.

58 GEOSCIENCES↗

Fire, dust, air and water: Improving aerosol biogeochemistry interactions in ACME. Final report

Aerosol-biogeochemistry interactions will have important consequences for human health, ecosystem services, and climate feedbacks over the next several decades. In this proposal we seek to link how human activities on land, including direct industrial emissions and landscape fires, impact the atmosphere and the biogeochemistry of the ocean and land. Our overall goal was to complement existing work at DOE labs to improve the capability and performance of the E3SM model in short-term, high-resolution projections of climate. Projections of ocean uptake of anthropogenic carbon are highly dependent on the details of the ocean physical and biogeochemistry parameterizations. Potentially equally important are the estimates of increases in iron and soluble nitrogen deposition in aerosols to the open ocean, but there are substantial uncertainties in these feedbacks. The E3SM land model simulates natural and anthropogenic fires, but improvement is needed to more accurately simulate fire responses to climate variability as well as fire contributions to terrestrial climate feedbacks. This proposal addressed four of the six priorities within the ACME/E3SM call by improving the representation of fire, industrial combustion sources, and desert dust and adding iron and soluble iron to the atmospheric chemistry (1). Our work focused on improving ocean biogeochemistry and the response of the oceans to anthropogenic atmospheric iron and nitrogen (3). We also improved the representation of land disturbance on fires and fire mediated carbon cycle feedbacks in the land model (4). Finally we had planned to couple the land, atmosphere and ocean biogeochemical systems (6) and explore how this coupling changes the feedbacks in the climate system. Our work under this proposal was separated into the following 5 tasks: i) improve representation of fire dynamics in both boreal and tropical ecosystems, ii) improve estimates of industrial emissions of iron and phosphorus, iii) improve representation of atmospheric iron, phosphorus and nitrogen transport, chemistry and deposition, iv) improve representation of iron, phosphorus and nitrogen deposition impacts on ocean biogeochemistry, and v) couple emissions into atmosphere, to deposition onto land and ocean biogeochemistry. For all 5 tasks, we conducted offline experiments of the historical period to enable comparison with available aerosol and land surface observations. We had planned to conduct fully coupled simulations, however, these will be conducted as part of the coupling in version 3 and 4 of the E3SM model, and thus have not been completed.

54 ENVIRONMENTAL SCIENCES↗

Enhanced prediction of terrestrial feedbacks to the coastal carbon cycle:using machine learning to improve sub-grid biogeochemical processes

Focal Areas: The paper aims to improve predictability of the coastal carbon cycle through improved model representation and quantification of the terrestrial feedbacks to aquatic ecosystems. The paper covers two focus areas, 1. Improved predictive modeling of terrestrial and aquatic biogeochemistry through ML-enabled surrogate models, and 2. Employing machine learning to integrate multimodal data sets collected across terrestrial and aquatic ecosystems at high spatial and temporal resolution.

54 ENVIRONMENTAL SCIENCES↗

Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

Modeling↗

Soil oxygen dynamics: a key mediator of tile drainage impacts on coupled hydrological, biogeochemical, and crop systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

60 APPLIED LIFE SCIENCES↗