Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy systems 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 703 records · Page 39

ELM model simulations of Plum Island Ecosystems LTER low marsh site 2018-2020

Model simulations using the Department of Energy's Energy Exascale Earth System Model (E3SM) land model (ELM) with improved capabilities to represent vegetation response to salinity and inundation. The simulations were conducted for a tidal salt marsh at Plum Island Ecosystems Long Term Ecological Research (LTER) site near Rowley, Massachusetts, USA; the site is a low marsh dominated by Spartina alterniflora. The model was forced with site-specific meteorology, salinity and tidal cycles from 2018-2020. Four sets of model simulations are included and described below:1. Parameterization of the salinity response function. These simulations tested different combinations of values for optimal salinity and salinity tolerance.2. Model evaluation. This comparison conducted simulations using the default model, the salinity function only, the submergence function only, and both the salinity and submergence functions. 3. Salinity scenarios. These simulations used the 2018 salinity input data varied by -5 to +10 ppt salinity.4. Water level scenarios. These simulations used the tide height varied by -10 to +50 cm. These simulations were used to demonstrate how the salinity and submergence functions better represent carbon uptake by tidal salt marshes.The data package includes netCDF files used as forcing files for tide height and salinity, one for each year 2018-2020 at observed salinity concentrations, and an additional three forcing files in which salinity concentrations were varied 5 ppt lower, 5 ppt higher, and 10 ppt higher than the measured 2018 time series. Also included are python scripts for creating forcing files, plain text parameter and command files for running simulations, model outputs in netCDF format, and python scripts for visualizing outputs. Code for the modified E3SM model is archived in Sulman et al 2023 at doi:10.15485/1991625. More detail about files is provided in the README.md file.

54 ENVIRONMENTAL SCIENCES↗

E3SM: Improved Climate Prediction with Exascale Capability

The Energy Exascale Earth System Model (E3SM) project is an ongoing, state-of-the-science earth system modeling, simulation, and prediction effort that optimizes Department of Energy (DOE) computing resources to meet the science needs of the nation and the agency’s mission objectives. Climate simulation has become a proven tool for identifying and quantifying the impacts of climate change, but even greater accuracy is required at all levels to improve forecast precision. Understanding the impact of climate change on global and regional water cycles is one of the highest priorities and most difficult challenges in climate change prediction. As part of a subproject of DOE’s Exascale Computing Project, a multidisciplinary team including geophysical and computational scientists developed a multiscale modeling framework (MMF) to refine cloud representation in E3SM climate simulation on GPU accelerated supercomputers, making higher resolution, more computationally efficient predictions possible.

54 ENVIRONMENTAL SCIENCES↗

CERES ERBE-like Instantaneous TOA Estimates (ES-8) in HDF (CER_ES8_TRMM-PFM_Edition1)

The ES-8 archival data product contains a 24-hour, single-satellite, instantaneous view of scanner fluxes at the top-of-atmosphere (TOA) reduced from spacecraft altitude unfiltered radiances using Earth Radiation Budget Experiment (ERBE) scanner Inversion algorithms and the ERBE shortwave (SW) and longwave (LW) Angular Distribution Models (ADMs). The ES-8 also includes the total (TOT), SW, LW, and window (WN) channel radiometric data; SW, LW, and WN unfiltered radiance values; and the ERBE scene identification for each measurement. These data are organized according to the CERES 3.3-second scan into 6.6-second records. As long as there is one valid scanner measurement within a record, the ES-8 record will be generated. The following CERES ES8 data sets are currently available: CER_ES8_TRMM-PFM_Edition1 CER_ES8_TRMM-PFM_Edition2 CER_ES8_TRMM-PFM_Transient-Ops2 CER_ES8_Terra-FM1_Edition1 CER_ES8_Terra-FM2_Edition1 CER_ES8_Terra-FM1_Edition2 CER_ES8_Terra-FM2_Edition2 CER_ES8_Aqua-FM3_Edition1 CER_ES8_Aqua-FM4_Edition1 CER_ES8_Aqua-FM3_Edition2 CER_ES8_Aqua-FM4_Edition2 CER_ES8_Aqua-FM3_Edition1-CV CER_ES8_Aqua-FM4_Edition1-CV CER_ES8_Terra-FM1_Edition1-CV CER_ES8_Terra-FM1_Edition1-CV. [Location=GLOBAL] [Temporal_Coverage: Start_Date=1997-12-27; Stop_Date=1998-08-31] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Temporal_Resolution=1 day; Temporal_Resolution_Range=Daily - < Weekly].

LONGWAVE FLUX↗

CERES ERBE-like Instantaneous TOA Estimates (ES-8) in HDF (CER_ES8_Aqua-FM4_Edition1)

The ES-8 archival data product contains a 24-hour, single-satellite, instantaneous view of scanner fluxes at the top-of-atmosphere (TOA) reduced from spacecraft altitude unfiltered radiances using Earth Radiation Budget Experiment (ERBE) scanner Inversion algorithms and the ERBE shortwave (SW) and longwave (LW) Angular Distribution Models (ADMs). The ES-8 also includes the total (TOT), SW, LW, and window (WN) channel radiometric data; SW, LW, and WN unfiltered radiance values; and the ERBE scene identification for each measurement. These data are organized according to the CERES 3.3-second scan into 6.6-second records. As long as there is one valid scanner measurement within a record, the ES-8 record will be generated. The following CERES ES8 data sets are currently available: CER_ES8_TRMM-PFM_Edition1 CER_ES8_TRMM-PFM_Edition2 CER_ES8_TRMM-PFM_Transient-Ops2 CER_ES8_Terra-FM1_Edition1 CER_ES8_Terra-FM2_Edition1 CER_ES8_Terra-FM1_Edition2 CER_ES8_Terra-FM2_Edition2 CER_ES8_Aqua-FM3_Edition1 CER_ES8_Aqua-FM4_Edition1 CER_ES8_Aqua-FM3_Edition2 CER_ES8_Aqua-FM4_Edition2 CER_ES8_Aqua-FM3_Edition1-CV CER_ES8_Aqua-FM4_Edition1-CV CER_ES8_Terra-FM1_Edition1-CV CER_ES8_Terra-FM1_Edition1-CV. [Location=GLOBAL] [Temporal_Coverage: Start_Date=1997-12-27; Stop_Date=2005-03-29] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Temporal_Resolution=1 day; Temporal_Resolution_Range=Daily - < Weekly].

SHORTWAVE FILTERED RADIANCE↗

Improving Aerosol Radiative Forcing and Climate in E3SM: Impacts of New Cloud Microphysics and Improved Wet Removal Treatments

Abstract Numerous Earth system models exhibit excessive aerosol effective forcing at the top of the atmosphere (TOA), including the Department of Energy's Energy Exascale Earth System Model (E3SM). Here, in the context of the E3SM version 3 effort, the predicted particle property (P3) stratiform cloud microphysics scheme and an enhanced deep convection parameterization suite (ZM_plus) are implemented into E3SM. The ZM_plus includes a convective cloud microphysics scheme, a multi‐scale coherent structure parameterization for mesoscale convective systems, and a revised cloud base mass flux formulation considering impacts of the large‐scale environment. The P3 scheme improved cloud and radiation particularly over the Northern Hemisphere and the frequency of heavy precipitation over the tropics, and the ZM_plus improved clouds in the tropics. P3 decreases aerosol effective forcing by 0.15 W m −2 , while the ZM_plus increases it by 0.27 W m −2 , resulting from excessive direct (0.31 W m −2 ) and indirect forcing (−1.79 W m −2 ). The excessive aerosol forcings are due to aerosol overestimation associated with insufficient aerosol wet removal. By improving the physical treatments in the aerosol wet removal, we effectively mitigate anthropogenic aerosol overestimation and thus attenuate direct (0.09 W m −2 ) and indirect aerosol forcing (−1.52 W m −2 ). Adjustment to primary organic matter hygroscopicity reduces direct and indirect forcing to more reasonable values: −0.13 W m −2 and −1.31 W m −2 , respectively. On climatology, improved aerosol treatments mitigate overestimation of aerosol optical depth.

54 ENVIRONMENTAL SCIENCES↗

EMPIRICAL VALIDATION OF MULTI-ZONE BUILDING AND HVAC SYSTEM MODELS UNDER UNCERTAINTY

This study implemented a framework of empirical validation of building energy models under uncertainty to a set of controlled experiments that aim to validate multi-zone building and HVAC system models. Energy models were created through iterative acquisitions of information and data and uses measurement data from various types of sensors as both inputs and as observations to validate predictions. Experimental and modeling uncertainties were quantified and propagated accordingly, and probabilistic accuracy metrics were used to evaluate the agreement between model predictions and observations under uncertainty. Sensitivity analysis was performed to identify the most influential uncertainties that will be prioritized to be addressed in the next steps. Current results of two cooling tests show an overall good agreement between predictions and observations on a set of HVAC system outputs despite considerable and influential uncertainty in DX cooling coil COP. Agreements on zone-level responses vary notably among individual rooms, likely because of significant uncertainties in room radiation heat gain and system supply air.

Li, Qi↗

CERES ERBE-like Instantaneous TOA Estimates (ES-8) in HDF (CER_ES8_TRMM-PFM_Transient-Ops2)

The ES-8 archival data product contains a 24-hour, single-satellite, instantaneous view of scanner fluxes at the top-of-atmosphere (TOA) reduced from spacecraft altitude unfiltered radiances using Earth Radiation Budget Experiment (ERBE) scanner Inversion algorithms and the ERBE shortwave (SW) and longwave (LW) Angular Distribution Models (ADMs). The ES-8 also includes the total (TOT), SW, LW, and window (WN) channel radiometric data; SW, LW, and WN unfiltered radiance values; and the ERBE scene identification for each measurement. These data are organized according to the CERES 3.3-second scan into 6.6-second records. As long as there is one valid scanner measurement within a record, the ES-8 record will be generated. The following CERES ES8 data sets are currently available: CER_ES8_TRMM-PFM_Edition1 CER_ES8_TRMM-PFM_Edition2 CER_ES8_TRMM-PFM_Transient-Ops2 CER_ES8_Terra-FM1_Edition1 CER_ES8_Terra-FM2_Edition1 CER_ES8_Terra-FM1_Edition2 CER_ES8_Terra-FM2_Edition2 CER_ES8_Aqua-FM3_Edition1 CER_ES8_Aqua-FM4_Edition1 CER_ES8_Aqua-FM3_Edition2 CER_ES8_Aqua-FM4_Edition2 CER_ES8_Aqua-FM3_Edition1-CV CER_ES8_Aqua-FM4_Edition1-CV CER_ES8_Terra-FM1_Edition1-CV CER_ES8_Terra-FM1_Edition1-CV. [Location=GLOBAL] [Temporal_Coverage: Start_Date=1997-12-27; Stop_Date=1999-07-17] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Temporal_Resolution=1 day; Temporal_Resolution_Range=Daily - < Weekly].

TOTAL FILTERED RADIANCE↗

Machine learning feature analysis illuminates disparity between E3SM climate models and observed climate change

In September of 2020, Arctic sea ice extent was the second-lowest on record. State of the art climate prediction uses Earth system models (ESMs), driven by systems of differential equations representing the laws of physics. Previously, these models have tended to underestimate Arctic sea ice loss. The issue is grave because accurate modeling is critical for economic, ecological, and geopolitical planning. We use machine learning techniques, including random forest regression and Gini importance, to show that the Energy Exascale Earth System Model (E3SM) relies too heavily on just one of the ten chosen climatological quantities to predict September sea ice averages. Furthermore, E3SM gives too much importance to six of those quantities when compared to observed data. Finally, identifying the features that climate models incorrectly rely on should allow climatologists to improve prediction accuracy.

54 ENVIRONMENTAL SCIENCES↗

CERES Single Scanner Satellite Footprint, TOA, Surface Fluxes and Clouds (SSF) data in HDF (CER_SSF_Terra-FM1-MODIS_Edition2A)

The Single Scanner Footprint TOA/Surface Fluxes and Clouds (SSF) product contains one hour of instantaneous Clouds and the Earth's Radiant Energy System (CERES) data for a single scanner instrument. The SSF combines instantaneous CERES data with scene information from a higher-resolution imager such as Visible/Infrared Scanner (VIRS) on TRMM or Moderate-Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua. Scene identification and cloud properties are defined at the higher imager resolution and these data are averaged over the larger CERES footprint. For each CERES footprint, the SSF contains the number of cloud layers and for each layer the cloud amount, height, temperature, pressure, optical depth, emissivity, ice and liquid water path, and water particle size. The SSF also contains the CERES filtered radiances for the total, shortwave (SW), and window (WN) channels and the unfiltered SW, longwave (LW), and WN radiances. The SW, LW, and WN radiances at spacecraft altitude are converted to Top-of-the-Atmosphere (TOA) fluxes based on the imager defined scene. These TOA fluxes are used to estimate surface fluxes. Only footprints with adequate imager coverage are included on CER_SSF_TRMM-PFM-VIRS_Subset_Edition1the SSF which is much less than the full set of footprints on the CERES ES-8 product. The following CERES SSF data sets are currently available: CER_SSF_TRMM-PFM-VIRS_Edition1 CER_SSF_TRMM-PFM-VIRS_Subset_Edition1 CER_SSF_TRMM-PFM-VIRS_Edition2A CER_SSF_TRMM-SIM-VIRS_Edition2_VIRSonly CER_SSF_TRMM-PFM-VIRS_Edition2A-TransOps CER_SSF_TRMM-PFM-VIRS_Edition2B-TransOps CER_SSF_TRMM-PFM-VIRS_Edition2B CER_SSF_Terra-FM1-MODIS_Edition1A CER_SSF_Terra-FM1-MODIS_Edition1A CER_SSF_Terra-FM1-MODIS_Edition2A CER_SSF_Terra-FM2-MODIS_Edition2A CER_SSF_Terra-FM1-MODIS_Edition2B CER_SSF_Terra-FM2-MODIS_Edition2B CER_SSF_Aqua-FM4-MODIS_Beta1 CER_SSF_Aqua-FM3-MODIS_Beta2 CER_SSF_Aqua-FM4-MODIS_Beta2. [Location=GLOBAL] [Temporal_Coverage: Start_Date=1998-01-01; Stop_Date=2003-12-31] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Temporal_Resolution=1 hour; Temporal_Resolution_Range=Hourly - < Daily].

ICE WATER PATH↗

Ocean physical‐biogeochemical interactions in the CMIP6 and E3SM Earth System Models

The primary objectives of this project are to quantify the spatio‐temporal linkages between physical climate and ocean biogeochemical variables in the U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SM), in CMIP6 Earth System Models (ESMs), and in historical observations, to evaluate the E3SM’s ability to reproduce the observations, and to develop a solid diagnostic package to do so. It will revolve around four hypotheses and will focus on key physical and biogeochemical processes that regulate the upper ocean carbon and oxygen cycling: water mass distribution, ventilation, and biological production.

54 ENVIRONMENTAL SCIENCES↗

CERES Single Scanner Satellite Footprint, TOA, Surface Fluxes and Clouds (SSF) data in HDF (CER_SSF_Terra-FM2-MODIS_Edition2A)

The Single Scanner Footprint TOA/Surface Fluxes and Clouds (SSF) product contains one hour of instantaneous Clouds and the Earth's Radiant Energy System (CERES) data for a single scanner instrument. The SSF combines instantaneous CERES data with scene information from a higher-resolution imager such as Visible/Infrared Scanner (VIRS) on TRMM or Moderate-Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua. Scene identification and cloud properties are defined at the higher imager resolution and these data are averaged over the larger CERES footprint. For each CERES footprint, the SSF contains the number of cloud layers and for each layer the cloud amount, height, temperature, pressure, optical depth, emissivity, ice and liquid water path, and water particle size. The SSF also contains the CERES filtered radiances for the total, shortwave (SW), and window (WN) channels and the unfiltered SW, longwave (LW), and WN radiances. The SW, LW, and WN radiances at spacecraft altitude are converted to Top-of-the-Atmosphere (TOA) fluxes based on the imager defined scene. These TOA fluxes are used to estimate surface fluxes. Only footprints with adequate imager coverage are included on CER_SSF_TRMM-PFM-VIRS_Subset_Edition1the SSF which is much less than the full set of footprints on the CERES ES-8 product. The following CERES SSF data sets are currently available: CER_SSF_TRMM-PFM-VIRS_Edition1 CER_SSF_TRMM-PFM-VIRS_Subset_Edition1 CER_SSF_TRMM-PFM-VIRS_Edition2A CER_SSF_TRMM-SIM-VIRS_Edition2_VIRSonly CER_SSF_TRMM-PFM-VIRS_Edition2A-TransOps CER_SSF_TRMM-PFM-VIRS_Edition2B-TransOps CER_SSF_TRMM-PFM-VIRS_Edition2B CER_SSF_Terra-FM1-MODIS_Edition1A CER_SSF_Terra-FM1-MODIS_Edition1A CER_SSF_Terra-FM1-MODIS_Edition2A CER_SSF_Terra-FM2-MODIS_Edition2A CER_SSF_Terra-FM1-MODIS_Edition2B CER_SSF_Terra-FM2-MODIS_Edition2B CER_SSF_Aqua-FM4-MODIS_Beta1 CER_SSF_Aqua-FM3-MODIS_Beta2 CER_SSF_Aqua-FM4-MODIS_Beta2. [Location=GLOBAL] [Temporal_Coverage: Start_Date=1998-01-01; Stop_Date=2003-12-31] [Spatial_Coverage: Southernmost_Latitude=-90; Northernmost_Latitude=90; Westernmost_Longitude=-180; Easternmost_Longitude=180] [Data_Resolution: Temporal_Resolution=1 hour; Temporal_Resolution_Range=Hourly - < Daily].

SURFACE (RADIATIVE) FLUX↗

Evaluation of the interactive stratospheric ozone (O3v2) module in the E3SM version 1 Earth system model

Stratospheric ozone affects climate directly as the predominant heat source in the stratosphere and indirectly through chemical reactions controlling other greenhouse gases. The U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) implemented a new ozone chemistry module that improves the simulation of the sharp tropopause gradients, replacing a version based partly on long-term average climatologies that poorly represented heating rates in the lowermost stratosphere. The new O3v2 module extends seamlessly into the troposphere and preserves the naturally sharp cross-tropopause gradient, with 20 %–40 % less ozone in this region. Additionally, O3v2 enables the diagnosis of stratosphere–troposphere exchange flux of ozone, a key budget term lacking in E3SMv1. Here, we evaluate key features in ozone abundance and other closely related quantities in atmosphere-only E3SMv1 simulations driven by observed sea surface temperatures (SSTs, years 1990–2014), comparing them with satellite observations of ozone and also with the University of California, Irvine chemistry transport model (UCI CTM) using the same stratospheric chemistry scheme but driven by European Centre forecast fields for the same period. In terms of stratospheric column ozone, O3v2 shows reduced mean bias and improved northern midlatitude variability, but it is not quite as good as the UCI CTM. As expected, SST-forced E3SMv1 simulations cannot synchronize with observed quasi-biennial oscillations (QBOs), but they do show the typical QBO pattern seen in column ozone. This new O3v2 E3SMv1 model mostly retains the same climate state and climate sensitivity as the previous version, and we recommend its use for other climate models that still use ozone climatologies.

54 ENVIRONMENTAL SCIENCES↗

Evaluation and modification of ELM seasonal deciduous phenology against observations in a southern boreal peatland forest

Phenological transitions determine the timing of changes in land surface properties and the seasonality of exchanges of biosphere-atmosphere energy, water, and carbon. Accurate mechanistic modeling of phenological processes is therefore critical to understand and correctly predict terrestrial ecosystem feedbacks with changing atmospheric and climate conditions. However, the phenological components in the land model of the US Department of Energy's (DOE) Energy Exascale Earth System Model (ELM of E3SM) were previously unable to accurately capture the observed phenological responses to environmental conditions in a well-studied boreal peatland forest. In this research, we introduced new seasonal-deciduous phenology schemes into version 1.0 of ELM and evaluated their performance against the PhenoCam observations at the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experiment in northern Minnesota from 2015 to 2018. We found that phenology simulated by the revised ELM (i.e., earlier spring onsets and stronger warming responses of spring onsets and autumn senescence) was closer to observations than simulations from the original algorithms for both the deciduous conifer (Larix laricina) and mixed shrub layers. Moreover, the revised ELM generally produced higher carbon and water fluxes (e.g., photosynthesis and evapotranspiration) during the growing season and stronger flux responses to warming than the default ELM. A parameter sensitivity analysis further indicated the significant contribution of phenology parameters to uncertainty in key carbon and water cycle variables, underscoring the importance of precise phenology parameterization. Furthermore, this phenological modeling effort demonstrates the potential to enhance the E3SM representation of land-climate interactions at broader spatiotemporal scales, especially under anticipated elevated CO 2 and warming conditions.

54 ENVIRONMENTAL SCIENCES↗

Pyrocumulonimbus Events Over British Columbia in 2017: An Ensemble Model Study of Parameter Sensitivities and Climate Impacts of Wildfire Smoke in the Stratosphere

Abstract The Pyrocumulonimbus (pyroCb) events over British Columbia in 2017 were observed in the lower stratosphere for about 8–10 months after the smoke injections. Several previous studies used global climate models to investigate the physical parameters for the 2017 pyroCb events, but the conclusions show strong model dependency. In this study, we use the Energy Exascale Earth System Model (E3SM) and complete an ensemble of runs exploring three injection parameters: smoke aerosol mass, the percentage of black carbon within the smoke aerosols, and plume injection height. Additionally, we consider the heterogeneous reaction of ozone and primary organic matter. According to the satellite daily observed aerosol optical depth (AOD), we find that the best ensemble member is the simulation with 0.4 Tg of smoke, 3% of which is black carbon, a 13.5 km smoke injection height, and a 10 −5 probability factor of the heterogeneous reaction. Besides AOD, we examine the ensemble score based on the metrics used in previous studies: the metrics of the extinction coefficient at 18 km altitude and the maximum plume height. The conclusion of the best estimate of the injection parameters for 2017 pyroCb events shows strong not only model but also evaluation metric dependency. We use the Random Forest machine learning technique to quantify the relative importance of each parameter in accurately simulating the 2017 pyroCb events and find that the injection height is the most critical feature, no matter which metric is used to score the ensemble members.

54 ENVIRONMENTAL SCIENCES↗

Using Satellite and ARM Observations to Evaluate Cold Air Outbreak Cloud Transitions in E3SM Global Storm‐Resolving Simulations

Abstract This study examines marine boundary layer cloud regime transition during a cold air outbreak (CAO) over the Norwegian Sea, simulated by a global storm‐resolving model (GSRM) known as the Simple Cloud‐Resolving Energy Exascale Earth System Model Atmosphere Model (SCREAM). By selecting observational references based on a combination of large‐scale conditions rather than strict time‐matched comparisons, this study finds that SCREAM qualitatively captures the CAO cloud transition, including boundary layer growth, cloud mesoscale structure, and phase partitioning. SCREAM also accurately locates the greatest ice and liquid in the mesoscale updrafts, however, underestimates supercooled liquid water in cumulus clouds. The model evaluation approach adopted by this study takes advantages of the existing computational‐expensive global simulations of GSRM and the available observations to understand model performance and can be applied to assessments of other cloud regimes in different regions. Such practice provides valuable guidance on the future effort to correct and improve biased model behaviors.

54 ENVIRONMENTAL SCIENCES↗

Vertical canopy gradients of respiration drive plant carbon budgets and leaf area index

Despite its importance for determining global carbon fluxes, leaf respiration remains poorly constrained in land surface models (LSMs). We tested the sensitivity of the Energy Exascale Earth System Model Land Model – Functionally Assembled Terrestrial Ecosystem Simulator (ELM-FATES) to variation in the canopy gradients of leaf maintenance respiration (R dark ). We ran global and point simulations varying the canopy gradient of R dark to explore the impacts on forest structure, composition, and carbon cycling. In global simulations, steeper canopy gradients of R dark lead to increased understory survival and leaf biomass. Leaf area index (LAI) increased up to 77% in tropical regions compared with the default parameterization, improving alignment with remotely sensed benchmarks. Global vegetation carbon varied from 308 Pg C to 449 Pg C across the ensemble. In tropical forest simulations, steeper gradients of R dark had a large impact on successional dynamics. Results show the importance of canopy gradients in leaf traits and fluxes for determining plant carbon budgets and emergent ecosystem properties such as competitive dynamics, LAI, and vegetation carbon. The high-model sensitivity to canopy gradients in R dark highlights the need for more observations of how leaf traits and fluxes vary along light micro-environments to inform critical dynamics in LSMs.

59 BASIC BIOLOGICAL SCIENCES↗

Representing lateral groundwater flow from land to river in Earth system models

Lateral groundwater flow (LGF) is an important hydrologic process in controlling water table dynamics. Due to the relatively coarse spatial resolutions of land surface models, the representation of this process is often overlooked or overly simplified. In this study, we developed a hillslope-based lateral groundwater flow model. Specifically, we first developed a hillslope definition model based on an existing watershed delineation model to represent the subgrid spatial variability in topography. Building upon this hillslope definition, we then developed a physical-based lateral groundwater flow using Darcy’s equation. This model explicitly considers the relationships between the groundwater table along the hillslope and the river water table levels. We coupled this intra-grid model to the land component (E3SM Land Model: ELM) and river component (MOdel for Scale Adaptive River Transport: MOSART) of the Energy Exascale Earth System Model (E3SM). We tested both the hillslope definition model and the lateral groundwater flow model and performed sensitivity experiments using different configurations. Simulations for a single grid cell at 0.5°×0.5° within the Amazon basin show that the definition of hillslope is the key to modeling lateral flow processes and the runoff partition between surface and subsurface can be dramatically changed using the hillslope approach. Although our method provides a pathway to improve the lateral flow process, future improvements are needed to better capture the subgrid structure to account for the spatial variability in hillslopes within the simulated grid of land surface models.

54 ENVIRONMENTAL SCIENCES↗