Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth Systems”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Evaluation of Stratocumulus Evolution Under Contrasting Temperature Advections in CESM2 Through a Lagrangian Framework

Abstract This study leveraged a Lagrangian framework to examine the evolution of stratocumulus clouds under cold and warm advections (CADV and WADV) in the Community Earth System Model 2 (CESM2) against observations. We found that CESM2 simulates a too rapid decline in low‐cloud fraction (LCF) and cloud liquid water path (CLWP) under CADV conditions, while it better aligns closely with observed LCF under WADV conditions but overestimates the increase in CLWP. Employing an explainable machine learning approach, we found that too rapid decreases in LCF and CLWP under CADV conditions are related to overestimated drying effects induced by sea surface temperature, whereas the substantial increase in CLWP under WADV conditions is associated with the overestimated moistening effects due to free‐tropospheric moisture and surface winds. Our findings suggest that overestimated drying effects of sea surface temperature on cloud properties might be one of crucial causes for the high equilibrium climate sensitivity in CESM2.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Potential of Stratospheric Aerosol Injection to Reduce the Climatic Risks of Explosive Volcanic Eruptions

Sulfur-rich volcanic eruptions happen sporadically. If Stratospheric Aerosol Injection (SAI) were to be deployed, it is likely that explosive volcanic eruptions would happen during such a deployment. Here we use an ensemble of Earth System Model simulations to show how changing the injection strategy post-eruption could be used to reduce the climate risks of a large volcanic eruption; the risks are also modified even without any change to the strategy. For a medium-size eruption (10 Tg-SO 2 ) comparable to the SAI injection rate, the volcanic-induced cooling would be reduced if it occurs under SAI, especially if artificial sulfur dioxide injections were immediately suspended. Alternatively, suspending injection only in the eruption hemisphere and continuing injection in the opposite would reduce shifts in precipitation in the tropical belt and thus mitigate eruption-induced drought. Finally, we show that for eruptions much larger than the SAI deployment, changes in SAI strategy would have minimal effect.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Influence of Climate Variability and Future Climate Change on Atlantic Hurricane Season Length

Abstract Atlantic hurricane season length is important for emergency management preparation, motivating the need to understand its variability and change. We investigated the influence of ocean variability on Atlantic hurricane season length in observations and a future climate simulated by the Energy Exascale Earth System Model (E3SM). We found that multiple factors influence hurricane season length, through their influence on season start and end. Warm western subtropical Atlantic sea‐surface temperature anomalies (SSTAs) during boreal spring (before the official hurricane season start) drive early starts to the hurricane season, and vice versa for cool SSTAs. Meanwhile, La Niña in autumn (before the official hurricane season end) drives late ends to the hurricane season, and vice versa for El Niño. E3SM projects a 27‐day increase in future Atlantic hurricane season length given La Niña and warm northern tropical Atlantic SSTAs. This research documents sources of predictability for Atlantic hurricane season length.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Enhancing Extreme Precipitation Predictions With Dynamical Downscaling: A Convection‐Permitting Modeling Study in Texas and Oklahoma

Precipitation in the Southern Plains of the United States is relatively well depicted by the Community Earth System Model (CESM). However, despite its ability to capture seasonal mean precipitation anomalies, CESM consistently underestimates extreme pluvial and drought events, rendering it an insufficient tool for extending simulation lead times for exceptional events, such as the abnormally dry May 2011, which helped drive Texas into its worst period of drought in more than a century, and the abnormally wet May 2015, which led to widespread flooding in that state. Ensemble‐based regional climate experiments are completed for the two extreme years using Weather Research and Forecasting model (WRF) and downscaled from CESM. WRF simulations are at convection‐permitting grid resolution for improved physical representation of simulated precipitation over the Southern Great Plains. By integrating convection‐permitting models (CPMs) into each individual member of a CESM ocean data assimilation ensemble, this study demonstrates that high‐resolution dynamical downscaling can improve model skillfulness at capturing these two events and is thus a potentially useful tool for forecasting extremely high and extremely low precipitation events at subseasonal or even seasonal lead times.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Characterizing Ice Nucleating Particles Over the Southern Ocean Using Simultaneous Aircraft and Ship Observations

Abstract Supercooled liquid clouds are ubiquitous over the Southern Ocean (SO), even to temperatures below −20°C, and comprise a large fraction of the marine boundary layer (MBL) clouds. Earth system models and reanalysis products have struggled to reproduce the observed cloud phase distribution and occurrence of cloud ice in the region. Recent simulations found the microphysical representation of ice nucleation and growth has a large impact on these properties, however, measurements of SO ice nucleating particles (INPs) to validate simulations are sparse. This study presents measurements of INPs from simultaneous aircraft and ship campaigns conducted over the SO in austral summer 2018, which include the first in situ observations in and above cloud in the region. Our results confirm recent observations that INP concentrations are uniformly lower than measurements made in the late 1960s. While INP concentrations below and above cloud are similar, higher ice nucleation efficiency above cloud supports model simulations that the dominant INP composition varies with height. Model parameterizations based solely on aerosol properties capture the mean relationship between INP concentration and temperature but not the observed variability, which is likely related to the only modest correlations observed between INPs and environmental or aerosol metrics. Including wind speed in addition to activation temperature in a marine INP parameterization reduces bias but does not explain the large range of observed INP concentrations. Direct and indirect inference of marine INP size suggests MBL INPs, at least during Austral summer, are dominated by particles with diameters smaller than 500 nm.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Surface and Atmospheric Heating Responses to Spectrally Resolved Albedo of Frozen and Liquid Water Surfaces

Abstract Multiple Earth system models (ESMs) discretize surface albedo into two semi‐broadbands comprising the UV/visible and near‐infrared (NIR) wavelengths. Here, we use an offline single‐column radiative transfer model to investigate the radiative effects of spectrally resolving the surface albedo. We use the Snow, Ice, and Aerosol Radiative model, extended to simulate liquid water, to calculate snow, ice, and liquid water albedo. We flux‐weight the hyperspectral albedo into the coarser spectral bands used by the atmospheric shortwave radiative transfer model. We establish representative atmospheric profiles for the three surface types and compare their shortwave fluxes and atmospheric warming rates with the spectrally resolved albedo to those calculated with the semi‐broadband approximation. Spectrally resolved surface albedo over snow and ice reduces atmospheric warming by darkening the albedo of NIR bands, correcting the too‐strong surface absorption in visible bands, and too‐weak surface absorption in shortwave infrared bands caused by the semi‐broadband approximation. We explore the effects on surface and atmospheric warming rates of varying solar zenith angle, cloud cover, relative humidity, and snow grain/air bubble radii. The semi‐broadband albedo biases can exceed 10% and 2% for the surface and atmospheric net flux respectively, being particularly strong under conditions which alter the distributions of surface insolation (i.e., cloud cover or increased atmospheric water vapor). These results show that transmitting a higher resolution spectral radiation field between the atmosphere and surface reduces biases in surface absorption and atmospheric heating present in ESMs that currently use the semi‐broadband approximation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Role of Cloud Radiative Effects in the Propagating Southern Annular Mode

Abstract The Southern Annular Mode (SAM) is the most dominant natural mode of variability in the mid‐latitudes of the Southern Hemisphere (SH). However, both the sign and magnitude of the feedbacks from the diabatic processes, especially those associated with clouds, onto the SAM remain elusive. By applying the cloud locking technique to the Energy Exascale Earth System Model (E3SM) atmosphere model, this study isolates the positive feedback from the cloud radiative effect (CRE) to the SAM. Feedback analysis based on a wave activity‐zonal momentum interaction framework corroborates this weak but positive feedback. While the magnitude of the CRE feedback appears to be secondary compared to the feedbacks from the dry and other diabatic processes, the indirect CRE effects through the interaction with other dynamical and thermodynamical processes appear to play as important a role as the direct CRE in the life cycle of the SAM. The cross‐EOF analysis further reveals the obstructive effect of the interactive CRE on the propagation mode of the SH zonal wind directly through the CRE wave source and/or indirectly through modulating other diabatic processes. As a result, the propagation mode becomes more persistent and the SAM it represents becomes more predictable when the interactive CRE is disabled by cloud locking. Future efforts on inter‐model comparisons of CRE‐denial experiments are important to build consensus on the dynamical feedback of CRE.

58 GEOSCIENCES↗
DOE OSTI2023

Microbial Models for Simulating Soil Carbon Dynamics: A Review

Abstract Soils store the largest amount of carbon (C) in the biosphere, and the C pool in soil is critical to the global C balance. Numerous microbial models have been developed over the last few decades to represent microbial processes that regulate the responses of soil organic carbon (SOC) to climate change. However, the representation of microbial processes varies, and how microbial processes are incorporated into SOC models has not been well explored. Here, we reviewed 71 microbial models to characterize the microbial processes incorporated into SOC models and analyzed variations in mechanistic complexity. We revealed that (a) four processes (microbial‐mediated decomposition, mineral interaction, microbial necromass recycling, and active and dormant microbial dynamics) are commonly incorporated in microbial models, (b) ∼48% of models simulate only one microbial process (i.e., microbial‐mediated decomposition) and 35% of models simulate two microbial processes: for example, microbial‐mediated decomposition and mineral interaction, (c) more than 80% microbial models use nonlinear equations, such as forward Michaelis‐Menten kinetics, to represent SOC decomposition, (d) the concept of persistence of SOC due to its intrinsic properties has been replaced by organo‐mineral interaction (∼39% of microbial models) that protects SOC from decomposition, and (e) various temperature and moisture modifiers and pH effects have been used to explain the environmental effect on microbial processes. In the future, to realistically incorporate microbial processes into Earth System Models, it is imperative to identify experimental evidence on rate limitation processes and firmly ground model structure on the field and laboratory data.

Environmental Sciences & Ecology↗
DOE OSTI2024

Upscaling Soil Organic Carbon Measurements at the Continental Scale Using Multivariate Clustering Analysis and Machine Learning

Abstract Estimates of soil organic carbon (SOC) stocks are essential for many environmental applications. However, significant inconsistencies exist in SOC stock estimates for the U.S. across current SOC maps. We propose a framework that combines unsupervised multivariate geographic clustering (MGC) and supervised Random Forests regression, improving SOC maps by capturing heterogeneous relationships with SOC drivers. We first used MGC to divide the U.S. into 20 SOC regions based on the similarity of covariates (soil biogeochemical, bioclimatic, biological, and physiographic variables). Subsequently, separate Random Forests models were trained for each SOC region, utilizing environmental covariates and SOC observations. Our estimated SOC stocks for the U.S. (52.6 ± 3.2 Pg for 0–30 cm and 108.3 ± 8.2 Pg for 0–100 cm depth) were within the range estimated by existing products like Harmonized World Soil Database, HWSD (46.7 Pg for 0–30 cm and 90.7 Pg for 0–100 cm depth) and SoilGrids 2.0 (45.7 Pg for 0–30 cm and 133.0 Pg for 0–100 cm depth). However, independent validation with soil profile data from the National Ecological Observatory Network showed that our approach ( R 2 = 0.51) outperformed the estimates obtained from Harmonized World Soil Database ( R 2 = 0.23) and SoilGrids 2.0 ( R 2 = 0.39) for the topsoil (0–30 cm). Uncertainty analysis (e.g., low representativeness and high coefficients of variation) identified regions requiring more measurements, such as Alaska and the deserts of the U.S. Southwest. Our approach effectively captures the heterogeneous relationships between widely available predictors and the current SOC baseline across regions, offering reliable SOC estimates at 1 km resolution for benchmarking Earth system models.

58 GEOSCIENCES↗
DOE OSTI2024

Functionally Assembled Terrestrial Ecosystem Simulator (FATES) for Hurricane Disturbance and Recovery

Tropical cyclones are an important cause of forest disturbance, and major storms caused severe structural damage and elevated tree mortality in coastal tropical forests. Model capabilities that can be used to understand post-hurricane forest recovery are still limited. We use a vegetation demography model, the Functionally Assembled Terrestrial Ecosystem Simulator, coupled with the Energy Exascale Earth System Model Land Model (ELM-FATES) to study the processes and the key factors regulating post-hurricane forest recovery. We implemented hurricane-induced forest damage, including defoliation, structural biomass reduction, and tree mortality, performed ensemble model simulations, and used random forest feature importance. For the simulation in the Luquillo Experimental Forest, Puerto Rico, we identified factors controlling the post-hurricane forest recovery, and quantified the sensitivity of key model parameters to the post-hurricane forest recovery. The results indicate a tendency for the Bisley forests to shift toward the light demanding plant functional type (PFT) when the pre-hurricane biomass between the light demanding and shade tolerant PFTs is nearly equal and forests experience hurricane disturbance with mortality >60% for both the two PFTs. Under more realistic conditions where the shade tolerant PFT is initially dominant, mortality >80% is required for a shift toward dominance of the light demanding PFT at Bisley. Hurricane mortality and background mortality are the two major factors regulating post-hurricane forest recovery in simulations. This research improves understanding of the ELM-FATES model behavior associated with hurricane disturbance and provides guidance for dynamic vegetation model development in representing hurricane induced forest damage with varied intensities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Sensitivity of Arctic Surface Temperature to Including a Comprehensive Ocean Interior Reflectance to the Ocean Surface Albedo Within the Fully Coupled CESM2

Abstract Almost all current climate models simplify the ocean surface albedo (OSA) by assuming the reflected solar energy without the ocean interior contribution. In this study, an improved ocean surface albedo scheme is incorporated into the Community Earth System Model version 2 (CESM2) to assess the sensitivity of Arctic surface temperature to including ocean interior reflectance to the OSA. Fully coupled CESM2 simulations with and without ocean interior reflectance are subsequently performed, we focus on the analysis of Arctic surface temperature responses. Incorporating ocean interior reflectance increases absorbed solar radiation and warms the ocean, enhancing seasonal heat storage and release across the Arctic Ocean, and increasing sea ice reduction and positive climate feedbacks that elevates Arctic surface temperature. Seasonal variations in air‐surface temperature differences induce changes in turbulent heat flux patterns, concurrently modifying dynamic advection and moisture processes that affect boundary layer humidity and low clouds, especially in winter. Based on partitioning physical processes in the thermodynamic energy equation, surface air warming is induced primarily through positive heating anomalies of vertical advection, latent heat release, and longwave radiative forcing. Through an examination of the surface energy budget, skin temperature warming is driven predominantly by increased downward longwave radiation, positive surface albedo feedback in summer, and increased conductive heat transport from the ocean particularly in winter. Significant effects of ocean interior reflectance on the Arctic Ocean, including sea surface warming and sea ice reduction, justify the importance of ocean interior reflectance in climate models for better understanding of ongoing Arctic climate changes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

A Refined Understanding of the Ice Cloud Longwave Scattering Effects in Climate Model

Abstract Because longwave (LW) absorption by greenhouse gases and clouds is more significant than the LW scattering effect by clouds, most climate models neglect cloud LW scattering to save computational costs. Ignoring cloud LW scattering directly overestimates outgoing longwave radiation (OLR). This study included ice‐cloud LW scattering treatment in the Exascale Energy Earth System Model (E3SM) version 2 and ran fully‐coupled simulations, prescribed sea surface temperature simulations, and offline radiative transfer calculations to comprehensively assess the impact of ice‐cloud LW scattering on global climate simulation. The instantaneous effect due to ice‐cloud LW scattering reduces the OLR by ∼1 W/m 2 on the global average and 2 W/m 2 on the tropical average. Tropospheric warming and high cloud amount reduction act to partially compensate for such instantaneous OLR reduction caused by the inclusion of LW scattering. When the simulation reaches the equilibrium, the surface warms by 0.66 K on average with respect to the simulation without LW scattering, with the Arctic surface temperature differences more than twice as large as that of the global mean. The impact of including LW scattering on the simulated climate change in response to 4 × CO 2 is also assessed. While including the cloud LW scattering does not significantly modify radiative forcing and total radiative feedback under such a scenario, it results in a 10% more positive cloud feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

On the Sensitivity of Coastal Hypoxia to Its External Physical Forcings

Abstract The development of low‐oxygen zones threatening marine life (hypoxia) occurs annually in multiple coastal regions of the world. The largest estuary of the continental United States, the Chesapeake Bay, typically has ≈10 km 3 of water with dioxygen concentrations <3 mg L −1 in July. As numerical methods for refining model resolutions in targeted areas are becoming common, there is interest in assessing the feasibility of simulating coastal hazards such as hypoxia in Earth System Models (ESMs). These coupled models are typically not constrained by observations and thus likely to feature systematic biases in their land, atmosphere, or ocean components. This study relies on four numerical experiments to evaluate the sensitivity of Chesapeake Bay hypoxia to changes or biases in its external physical forcings. Hypoxia exhibits only a minor decrease (−1.6%) after reducing the Bay's terrestrial freshwater discharge by 9.5% (but keeping terrestrial nutrient loadings the same). Changes in freshwater discharges have their largest impact on hypoxia during one extreme event (−37% during 2011 tropical storm Lee). Similarly, changing oceanic conditions on the shelf or their temporal frequency impact hypoxia by only 5%–6%, indicating that the latter is predominantly dictated by local conditions. Although these results are promising from the perspective of ESMs, additional components of ESMs will need to be evaluated before general conclusions can be reached. We notably speculate that the Bay's hypoxia would exhibit higher sensitivity to other forcings not examined here, notably air temperatures and nutrient loadings.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Causes of Reduced Climate Sensitivity in E3SM From Version 1 to Version 2

The effective climate sensitivity in the Department of Energy's Energy Exascale Earth System Model (E3SM) has decreased from 5.3 K in version 1 to 4.0 K in version 2. This reduction is mainly due to a weaker positive cloud feedback that leads to a stronger negative radiative feedback. Present-day atmosphere-only experiments with uniform 4 K sea surface temperature warming are used to separate the contributions of individual model modifications to the reduced cloud feedback. We find that the reduced cloud feedback is mostly driven by changes over the tropical marine low cloud regime, mainly related to a new trigger function for the deep convection scheme and modifications in the cloud microphysics scheme. The new trigger function helps weaken the low cloud reduction by increasing the cloud water detrainment at low levels from deep convection under warming. Changes to the formula of autoconversion rate from liquid to rain and an introduced minimum cloud droplet number concentration threshold in cloud microphysical calculations help sustain clouds against dissipation by suppressing precipitation generation with warming. In the midlatitudes, the increased Wegener-Bergeron-Findeisen (WBF) efficiency strongly reduces present-day liquid water and leads to a stronger negative cloud optical depth feedback. The reduced trade cumulus cloud feedback in v2 is closer to estimates from recent observational and large-eddy modeling studies but might not be due to the right physical reasons. The reduced mid-latitude cloud feedback may be more plausible because more realistic present-day mixed-phase clouds are produced through the change in the WBF efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Representing Anthropogenic Dust From Agricultural Sources in E3SMv1: Implementation, Evaluation, and Assessment of the Radiative Forcing

Dust emissions related to anthropogenic activities (i.e., anthropogenic dust) is not represented in most global climate models and its radiative impact remains unassessed. In this study, we develop a new and physically based method to parameterize anthropogenic dust emission from agricultural sources (i.e., agricultural dust, or AD) based on the DOE's Energy Exascale Earth System Model version 1 (E3SMv1). This method relates AD emission to the crop land use fraction in the E3SMv1 land component. Major AD sources simulated by our parameterization include those over Central America, the Sahel, North India, and North China. The annual averaged AD emission with diameter less than 10 μm is 567 Tg yr −1 in present-day (year 2000), which contributes to 13.3% of total dust emission. Model evaluation against satellite and ground-based observations shows that the new parameterization can represent AD emissions and global dust cycle reasonably well. We find that the total dust emission increases by 13.1% (495 Tg yr −1 ) from 1850 to 2000 due to the cropland land use fraction changes without considering the impact from climate changes. This induces a net dust direct effective radiative forcing of −0.041 W m −2 at top of the atmosphere. This dust-induced cooling exceeds 10% of the total anthropogenic aerosol direct effective radiative forcing from 1750 to 2014 estimated by the Intergovernmental Panel on Climate Change Sixth Assessment Report. Our findings indicate an important role of anthropogenic dust in global climate change, which should be included in future climate change assessments.

agricultural dust↗
DOE OSTI2024

Autocalibration of the E3SM Version 2 Atmosphere Model Using a PCA-Based Surrogate for Spatial Fields

Global Climate Model tuning (calibration) is a tedious and time-consuming process, with high-dimensional input and output fields. Experts typically tune by iteratively running climate simulations with hand-picked values of tuning parameters. Many, in both the statistical and climate literature, have proposed alternative calibration methods, but most are impractical or difficult to implement. We present a practical, robust, and rigorous calibration approach on the atmosphere-only model of the Department of Energy's Energy Exascale Earth System Model (E3SM) version 2. Our approach can be summarized into two main parts: (a) the training of a surrogate that predicts E3SM output in a fraction of the time compared to running E3SM, and (b) gradient-based parameter optimization. To train the surrogate, we generate a set of designed ensemble runs that span our input parameter space and use polynomial chaos expansions on a reduced output space to fit the E3SM output. We use this surrogate in an optimization scheme to identify values of the input parameters for which our model best matches gridded spatial fields of climate observations. To validate our choice of parameters, we run E3SMv2 with the optimal parameter values and compare prediction results to expertly-tuned simulations across 45 different output fields. This flexible, robust, and automated approach is straightforward to implement, and we demonstrate that the resulting model output matches present day climate observations as well or better than the corresponding output from expert tuned parameter values, while considering high-dimensional output and operating in a fraction of the time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Machine Learning Driven Sensitivity Analysis of E3SM Land Model Parameters for Wetland Methane Emissions

Methane (CH 4 ) is globally the second most critical greenhouse gas after carbon dioxide, contributing to 16%–25% of the observed atmospheric warming. Wetlands are the primary natural source of methane emissions globally. However, wetland methane emission estimates from biogeochemistry models contain considerable uncertainty. One of the main sources of this uncertainty arises from the numerous uncertain model parameters within various physical, biological, and chemical processes that influence methane production, oxidation, and transport. Sensitivity Analysis (SA) can help identify critical parameters for methane emission and achieve reduced biases and uncertainties in future projections. This study performs SA for 19 selected parameters responsible for critical biogeochemical processes in the methane module of the Energy Exascale Earth System Model (E3SM) land model (ELM). The impact of these parameters on various CH 4 fluxes is examined at 14 FLUXNET- CH 4 sites with diverse vegetation types. Given the extensive number of model simulations needed for global variance-based SA, we employ a machine learning (ML) algorithm to emulate the complex behavior of ELM methane biogeochemistry. We found that parameters linked to CH 4 production and diffusion generally present the highest sensitivities despite apparent seasonal variation. Comparing simulated emissions from perturbed parameter sets against FLUXNET-CH 4 observations revealed that better performances can be achieved at each site compared to the default parameter values. This presents a scope for further improving simulated emissions using parameter calibration with advanced optimization techniques.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impacts of Topography‐Driven Water Redistribution on Terrestrial Water Storage Change in California Through Ecosystem Responses

Abstract Lateral subsurface flow plays an essential role in sustaining the terrestrial ecosystem, but it is not explicitly represented in most Earth System Models. In this study, we implemented an explicit lateral saturated flow model into the E3SM land model (ELM). The model explicitly describes lateral flow in the saturated zone by representing, for each model grid, an idealized hillslope consisting of five hydrologically connected soil columns. We conducted three model experiments driven by 0.125° atmospheric forcing data during 1980–2015 over California using models of the default ELM, a modified version of ELM to enhance infiltration, and the model with the lateral saturated flow model. The simulated runoff, evapotranspiration, and terrestrial water storage anomaly (TWSA) from the three simulations were evaluated against available observations, and the model explicitly representing lateral flow performs best. The new model produces greater gridcell‐averaged evapotranspiration especially over the mountainous regions with moderate relief and seasonally dry climates. Most importantly, it improves the modeled seasonal variations, interannual variabilities, and the recent decadal decline of TWSA. Many of these improvements can be attributed to the enhanced ecosystem resilience to droughts as demonstrated by transpiration increases caused by lateral flow. Model sensitivity experiments suggest that subsurface runoff is most sensitive to the ratio between horizontal and vertical saturated hydraulic conductivity, followed by hillslope planforms (convergent, divergent, and uniform), number of columns, and lower boundary conditions. Future work should effectively characterize hillslopes in global models and explore the long‐term influences of lateral water movement on modeled biogeochemical cycle.

54 ENVIRONMENTAL SCIENCES↗