Results for “Model Coupling” — Page 50

NASA NTRS2011

Methane and Environmental Change during the Paleocene-Eocene Thermal Maximum (PETM): Modeling the PETM Onset as a Two-stage Event

An atmospheric CH4 box model coupled to a global carbon cycle box model is used to constrain the carbon emission associated with the PETM and assess the role of CH4 during this event. A range of atmospheric and oceanic emission scenarios representing different amounts, rates, and isotopic signatures of emitted carbon are used to model the PETM onset. The first 3 kyr of the onset, a pre-isotope excursion stage, is simulated by the atmospheric release of 900 to 1100 Pg C CH4 with a delta C-13 of -22 to - 30 %. For a global average warming of 3 deg C, a release of CO2 to the ocean and CH4 to the atmosphere totalling 900 to 1400 Pg C, with a delta C-13 of -50 to -60%, simulates the subsequent 1 -kyr isotope excursion stage. To explain the observations, the carbon must have been released over at most 500 years. The first stage results cannot be associated with any known PETM hypothesis. However, the second stage results are consistent with a methane hydrate source. More than a single source of carbon is required to explain the PETM onset.

Carozza, David A.↗

NASA NTRS2015

Downscaling GISS ModelE Boreal Summer Climate over Africa

The study examines the perceived added value of downscaling atmosphere-ocean global climate model simulations over Africa and adjacent oceans by a nested regional climate model. NASA/Goddard Institute for Space Studies (GISS) coupled ModelE simulations for June- September 1998-2002 are used to form lateral boundary conditions for synchronous simulations by the GISS RM3 regional climate model. The ModelE computational grid spacing is 2deg latitude by 2.5deg longitude and the RM3 grid spacing is 0.44deg. ModelE precipitation climatology for June-September 1998-2002 is shown to be a good proxy for 30-year means so results based on the 5-year sample are presumed to be generally representative. Comparison with observational evidence shows several discrepancies in ModelE configuration of the boreal summer inter-tropical convergence zone (ITCZ). One glaring shortcoming is that ModelE simulations do not advance the West African rain band northward during the summer to represent monsoon precipitation onset over the Sahel. Results for 1998-2002 show that onset simulation is an important added value produced by downscaling with RM3. ModelE Eastern South Atlantic Ocean computed sea-surface temperatures (SST) are some 4 K warmer than reanalysis, contributing to large positive biases in overlying surface air temperatures (Tsfc). ModelE Tsfc are also too warm over most of Africa. RM3 downscaling somewhat mitigates the magnitude of Tsfc biases over the African continent, it eliminates the ModelE double ITCZ over the Atlantic and it produces more realistic orographic precipitation maxima. Parallel ModelE and RM3 simulations with observed SST forcing (in place of the predicted ocean) lower Tsfc errors but have mixed impacts on circulation and precipitation biases. Downscaling improvements of the meridional movement of the rain band over West Africa and the configuration of orographic precipitation maxima are realized irrespective of the SST biases.

Regional climate models↗

NASA NTRS2014

A Practical Philosophy of Complex Climate Modelling

We give an overview of the practice of developing and using complex climate models, as seen from experiences in a major climate modelling center and through participation in the Coupled Model Intercomparison Project (CMIP).We discuss the construction and calibration of models; their evaluation, especially through use of out-of-sample tests; and their exploitation in multi-model ensembles to identify biases and make predictions. We stress that adequacy or utility of climate models is best assessed via their skill against more naive predictions. The framework we use for making inferences about reality using simulations is naturally Bayesian (in an informal sense), and has many points of contact with more familiar examples of scientific epistemology. While the use of complex simulations in science is a development that changes much in how science is done in practice, we argue that the concepts being applied fit very much into traditional practices of the scientific method, albeit those more often associated with laboratory work.

complex simulation↗

NASA NTRS2016

GEM-CEDAR Challenge: Poynting Flux at DMSP and Modeled Joule Heat

Poynting flux into the ionosphere measures the electromagnetic energy coming from the magnetosphere. This energy flux can vary greatly between quiet times and geomagnetic active times. As part of the Geospace Environment Modeling-coupling energetics and dynamics of atmospheric regions modeling challenge, physics-based models of the 3-D ionosphere and ionospheric electrodynamics solvers of magnetosphere models that specify Joule heat and empirical models specifying Poynting flux were run for six geomagnetic storm events of varying intensity. We compared model results with Poynting flux values along the DMSP-15 satellite track computed from ion drift meter and magnetic field observations. Although being a different quantity, Joule heat can in practice be correlated to incoming Poynting flux because the energy is dissipated primarily in high latitudes where Poynting flux is being deposited. Within the physics-based model group, we find mixed results with some models overestimating Joule heat and some models agreeing better with observed Poynting flux rates as integrated over auroral passes. In contrast, empirical models tend to underestimate integrated Poynting flux values. Modeled Joule heat or Poynting flux patterns often resemble the observed Poynting flux patterns on a large scale, but amplitudes can differ by a factor of 2 or larger due to the highly localized nature of observed Poynting flux deposition that is not captured by the models. In addition, the positioning of modeled patterns appear to be randomly shifted against the observed Poynting flux energy input. This study is the first to compare Poynting flux and Joule heat in a large variety of models of the ionosphere.

Rastaetter, Lutz↗

NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

NASA NTRS2022

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

NASA NTRS2022

Attribution of Stratospheric and Tropospheric Ozone Changes Between 1850 and 2014 in CMIP6 Models

We quantify the impacts of halogenated ozone-depleting substances (ODSs), greenhouse gases (GHGs), and short-lived ozone precursors on ozone changes between 1850 and 2014 using single-forcing perturbation simulations from several Earth system models with interactive chemistry participating in the Coupled Model Intercomparison Project Aerosol and Chemistry Model Intercomparison Project. We present the responses of ozone to individual forcings and an attribution of changes in ozone columns and vertically resolved stratospheric and tropospheric ozone to these forcings. We find that whilst substantial ODS-induced ozone loss dominates the stratospheric ozone changes since the 1970s, in agreement with previous studies, increases in tropospheric ozone due to increases in short-lived ozone precursors and methane since the 1950s make increasingly important contributions to total column ozone (TCO) changes. Increases in methane also lead to substantial extra-tropical stratospheric ozone increases. Impacts of nitrous oxide and carbon dioxide on stratospheric ozone are significant but their impacts on TCO are small overall due to several opposing factors and are also associated with large dynamical variability. The multi-model mean (MMM) results show a clear change in the stratospheric ozone trends after 2000 due to now declining ODSs, but the trends are generally not significantly positive, except in the extra-tropical upper stratosphere, due to relatively small changes in forcing over this period combined with large model uncertainty. Although the MMM ozone compares well with the observations, the inter-model differences are large primarily due to the large differences in the models' representation of ODS-induced ozone depletion.

Halogenated ozone-depleting substances↗

NASA NTRSDate not supplied

Beginning and Ending Droughts in a Warmer World

Droughts are typically caused by below-normal precipitation and typically terminate in response to high precipitation. In many regions, soil moisture is projected to decrease in a warmer world. At the same time, changes in both mean and extreme precipitation are expected as the global mean temperature increases. How will these precipitation changes interact with soil moisture changes to affect drought onset and termination? In this study, we use climate model projections from Phase Six of the Coupled Model Intercomparison Project (CMIP6) to investigate how changes in mean precipitation affect soil moisture drought in the context of future warming in the Mediterranean and Western North America. We show that the same precipitation deficits at a global warming level of 2 degrees Celsius increase the likelihood of both moderate and severe drought compared to the pre-industrial baseline. Further, we evaluate the effectiveness of a given precipitation anomaly to terminate a soil moisture drought under different warming levels. These analyses offer a better understanding of the cascading impacts of, and relationships between, precipitation, precipitation extremes, and drought termination under warmer global temperatures.

Drought↗

NASA NTRSDate not supplied

Beginning and Ending Droughts in a Warmer World

Droughts are typically caused by below-normal precipitation and typically terminate in response to high precipitation. In many regions, soil moisture is projected to decrease in a warmer world. At the same time, changes in both mean and extreme precipitation are expected as the global mean temperature increases. How will these precipitation changes interact with soil moisture changes to affect drought onset and termination? In this study, we use climate model projections from Phase Six of the Coupled Model Intercomparison Project (CMIP6) to investigate how changes in mean precipitation affect soil moisture drought in the context of future warming in the Mediterranean and Western North America. We show that the same precipitation deficits at a global warming level of 2 degrees Celsius increase the likelihood of both moderate and severe drought compared to the pre-industrial baseline. Further, we evaluate the effectiveness of a given precipitation anomaly to terminate a soil moisture drought under different warming levels. These analyses offer a better understanding of the cascading impacts of, and relationships between, precipitation, precipitation extremes, and drought termination under warmer global temperatures.

Drought↗

NASA NTRS2024

The Changing Influence of Precipitation on Soil Moisture Drought With Warming in the Mediterranean and Western North America

Anthropogenic climate change has already affected drought severity and risk across many regions, and climate models project additional increases in drought risk with future warming. Historically, droughts are typically caused by periods of below-normal precipitation and terminated by average or above-normal precipitation. In many regions, however, soil moisture is projected to decrease primarily through warming-driven increases in evaporative demand, potentially affecting the ability of negative precipitation anomalies to cause drought and positive precipitation anomalies to terminate drought. Here, we use climate model simulations from Phase Six of the Coupled Model Intercomparison Project (CMIP6) to investigate how different levels of warming (1, 2, and 3°C) affect the influence of precipitation on soil moisture drought in the Mediterranean and Western North America regions. We demonstrate that the same monthly precipitation deficits (25th percentile relative to a preindustrial baseline) at a global warming level of 2°C increase the probability of both surface and rootzone soil moisture drought by 29% in the Mediterranean and 32% and 6% in Western North America compared to the preindustrial baseline. Furthermore, the probability of a dry (25th percentile relative to a preindustrial baseline) surface soil moisture month given a high (75th percentile relative to a preindustrial baseline) precipitation month is 6 (Mediterranean) and 3 (Western North America) times more likely in a 2°C world compared to the preindustrial baseline. For these regions, warming will likely increase the risk of soil moisture drought during low precipitation periods while simultaneously reducing the efficacy of high precipitation periods to terminate droughts.

drought↗

NASA NTRS2018

Tropospheric Jet Response to Antarctic Ozone Depletion: An Update with Chemistry-Climate Model Initiative (CCMI) Models

The Southern Hemisphere (SH) zonal-mean circulation change in response to Antarctic ozone depletion is re-visited by examining a set of the latest model simulations archived for the Chemistry-Climate Model Initiative (CCMI) project. All models reasonably well reproduce Antarctic ozone depletion in the late 20th century. The related SH-summer circulation changes, such as a poleward intensification of westerly jet and a poleward expansion of the Hadley cell, are also well captured. All experiments exhibit quantitatively the same multi-model mean trend, irrespective of whether the ocean is coupled or prescribed. Results are also quantitatively similar to those derived from the Coupled Model Intercomparison Project phase 5 (CMIP5) high-top model simulations in which the stratospheric ozone is mostly prescribed with monthly- and zonally-averaged values. These results suggest that the ozone-hole-induced SH-summer circulation changes are robust across the models irrespective of the specific chemistry-atmosphere-ocean coupling.

chemistry-climate models↗

DOE OSTI2021

Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6

Abstract. The Scenario Model Intercomparison Project (ScenarioMIP) defines and coordinates the main set of future climate projections, based on concentration-driven simulations, within the Coupled Model Intercomparison Project phase 6 (CMIP6). This paper presents a range of its outcomes by synthesizing results from the participating global coupled Earth system models. We limit our scope to the analysis of strictly geophysical outcomes: mainly global averages and spatial patterns of change for surface air temperature and precipitation. We also compare CMIP6 projections to CMIP5 results, especially for those scenarios that were designed to provide continuity across the CMIP phases, at the same time highlighting important differences in forcing composition, as well as in results. The range of future temperature and precipitation changes by the end of the century (2081–2100) encompassing the Tier 1 experiments based on the Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and SSP1-1.9 spans a larger range of outcomes compared to CMIP5, due to higher warming (by close to 1.5 ∘C) reached at the upper end of the 5 %–95 % envelope of the highest scenario (SSP5-8.5). This is due to both the wider range of radiative forcing that the new scenarios cover and the higher climate sensitivities in some of the new models compared to their CMIP5 predecessors. Spatial patterns of change for temperature and precipitation averaged over models and scenarios have familiar features, and an analysis of their variations confirms model structural differences to be the dominant source of uncertainty. Models also differ with respect to the size and evolution of internal variability as measured by individual models' initial condition ensemble spreads, according to a set of initial condition ensemble simulations available under SSP3-7.0. These experiments suggest a tendency for internal variability to decrease along the course of the century in this scenario, a result that will benefit from further analysis over a larger set of models. Benefits of mitigation, all else being equal in terms of societal drivers, appear clearly when comparing scenarios developed under the same SSP but to which different degrees of mitigation have been applied. It is also found that a mild overshoot in temperature of a few decades around mid-century, as represented in SSP5-3.4OS, does not affect the end outcome of temperature and precipitation changes by 2100, which return to the same levels as those reached by the gradually increasing SSP4-3.4 (not erasing the possibility, however, that other aspects of the system may not be as easily reversible). Central estimates of the time at which the ensemble means of the different scenarios reach a given warming level might be biased by the inclusion of models that have shown faster warming in the historical period than the observed. Those estimates show all scenarios reaching 1.5 ∘C of warming compared to the 1850–1900 baseline in the second half of the current decade, with the time span between slow and fast warming covering between 20 and 27 years from present. The warming level of 2 ∘C of warming is reached as early as 2039 by the ensemble mean under SSP5-8.5 but as late as the mid-2060s under SSP1-2.6. The highest warming level considered (5 ∘C) is reached by the ensemble mean only under SSP5-8.5 and not until the mid-2090s.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Manipulating electromagnetic responses of metal wires at the deep subwavelength scale via both near- and far-field couplings

A hybrid coupling model containing both near- and far-field couplings is developed for radiating two-resonator structures. We demonstrate that the near- and far-field couplings make distinguished contributions to electromagnetic responses. Compared to the classical electromagnetically induced transparency configurations, the presence of far-field coupling provides more flexibility in tuning lineshapes. Planar metamaterials composed of metal wires are designed based on this model, and various electromagnetic responses are experimentally observed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS2017

Approach to Modeling Boundary Layer Ingestion Using a Fully Coupled Propulsion-RANS Model

Airframe-propulsion integration concepts that use boundary layer ingestion have the potential to reduce aircraft fuel burn. One concept that has been recently explored is NASA's Starc-ABL aircraft configuration, which offers the potential for 12% mission fuel burn reduction by using a turbo-electric propulsion system with an aft-mounted electrically driven boundary layer ingestion propulsor. This large potential for improved performance motivates a more detailed study of the boundary layer ingestion propulsor design, but to date, analyses of boundary layer ingestion have used uncoupled methods. These methods account for only aerodynamic effects on the propulsion system or propulsion system effects on the aerodynamics, but not both simultaneously. This work presents a new approach for building fully coupled propulsive-aerodynamic models of boundary layer ingestion propulsion systems. A 1D thermodynamic cycle analysis is coupled to a RANS simulation to model the Starc-ABL aft propulsor at a cruise condition and the effects variation in propulsor design on performance are examined. The results indicates that both propulsion and aerodynamic effects contribute equally toward the overall performance and that the fully coupled model yields substantially different results compared to uncoupled. The most significant finding is that boundary layer ingestion, while offering substantial fuel burn savings, introduces throttle dependent aerodynamics effects that need to be accounted for. This work represents a first step toward the multidisciplinary design optimization of boundary layer ingestion propulsion systems.

Boundary Layer Ingestion↗

NASA NTRS1998

Development of a Coupled Hydrological/Sediment Yield Model for a Watershed at Regional Level

Development of a hydrologic model for the study of environmental conservation requires a comprehensive understanding of individual-storm affecting hydrologic and sedimentologic processes. The hydrologic models that we are currently coupling are the Simulator for Hydrology and Energy Exchange at the Land Surface (SHEELS) and the Distributed Runoff Model (DRUM). SHEELS runs continuously to estimate surface energy fluxes and sub-surface soil water fluxes, while DRUM operates during and following precipitation events to predict surface runoff and peak flow through channel routing. The lateral re-distribution of surface water determined by DRUM is passed to SHEELS, which then adjusts soil water contents throughout the profile. The model SHEELS is well documented in Smith et al. (1993) and Laymen and Crosson (1995). The model DRUM is well documented in Vieux et al. (1990) and Vieux and Gauer (1994). The coupled hydrologic model, SHEELS/DRUM, does not simulate sedimentologic processes. The simulation of the sedimentologic process is important for environmental conservation planning and management. Therefore, we attempted to develop a conceptual frame work for coupling a sediment yield model with SHEELS/DRUM to estimate individual-storm sediment yield from a watershed at a regional level. The sediment yield model that will be used for this study is the Universal Soil Loss Equation (USLE) with some modifications to enable the model to predict individual-storm sediment yield. The predicted sediment yield does not include wind erosion and erosion caused by irrigation and snow melt. Units used for this study are those given by Foster et al. (1981) for SI units.

Rajbhandaril, Narayan↗

DOE OSTI2021

Coupling particle scale model and SuperDEM–CFD for multiscale simulation of biomass pyrolysis in a packed bed pyrolyzer

An efficient biomass pyrolysis process requires a comprehensive understanding of the chemical and physical phenomena that occur at multi-length and time scales. In this study, a multiscale computational approach was developed and validated for biomass pyrolysis in a packed-bed reactor by integrating pyrolysis kinetics, a particle scale model, and SuperDEM-CFD in open-source code MFiX. A one-dimensional particle-scale model that discretizes the characteristic length of biomass particle into layers was developed to predict the intraparticle phenomena inside a single particle. The 1D model was validated by comparing it with a single biomass particle pyrolysis experiment. Here, a recently developed SuperDEM-CFD model was employed to simulate the non-spherical particle-particle contact and fluid-particle interaction. The coupled model was applied to simulate the pyrolysis of cubic biomass particles in a packed bed and validated by comparing with experimental data. Simulation with and without particle-scale model was compared, and the effect of the gas-solid heat transfer models was also investigated.

09 BIOMASS FUELS↗

DOE OSTI2023

A coupled equilibrium boundary layer model with stable water isotopes and its application to local water recycling

The contribution of evapotranspiration (ET) to regional precipitation, known as “local water recycling”, is a key process in the water cycle that can affect water management. However, the amount of planetary boundary layer (PBL) moisture that arises from ET is highly uncertain due to complex atmosphere and land surface conditions. In this study, an idealized two-layer equilibrium boundary layer model was coupled with a stable water isotope module including HDO and H 2 18 O to constrain PBL growing season water transport processes. The model was validated using turbulent heat fluxes and isotope ratios of water vapor (δ m ) and precipitation (δ P ) measured at a cropland site and a nearby tall tower in the Upper Midwest, United States. The results show that the PBL equilibrium features of δ m and δ P are well-constrained by thermal and moisture equilibrium in the PBL. For this study region, the summer values of rain evaporation fraction (f) and local water recycling ratio (LRR) are estimated to be 0.09 and 0.29 ± 0.12, respectively. The summer LRR values for the years 2006–2010 were 0.35, 0.36, 0.17, 0.29, and 0.29, respectively. The small value of LRR in 2008 corresponded to a drought condition with the lowest precipitation and second lowest ET among the five years. The summer magnitude of the amount effect is –2.8‰ (mm day –1 ) –1 and –0.8‰ (mm day –1 ) –1 for HDO and H 2 18 O, respectively. The local water recycling is identified as a significant factor influencing the continental isotope effect. Cropland has likely changed the regional LRR by –7.6 to 19.5% under different pre-agriculture land use scenarios. In conclusion, the feedback processes revealed here indicate that local water recycling is expected to be weakened under drought conditions, but it will be enhanced if irrigation is applied more intensely with more frequent drought events as the climate continues to warm.

54 ENVIRONMENTAL SCIENCES↗