Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth system 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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Accelerated Climate Model for Energy (ACME): Ocean and Sea-Ice Processes

The primary goal of this project was to understand the role of the ocean and sea-ice components in the evolution of the climate in both a high-resolution fully-coupled 1850 preindustrial control and present day transient climate simulations run using the Energy Exascale Earth System Model (E3SM), formerly known as the Accelerated Climate Model for Energy (ACME). Two distinct approaches to the initialization of the present day transients were tested: in the first, an ensemble of three transients used atmospheric reanalysis forced 0.1-degree ocean/sea-ice states, while a fourth transient was branched off the preindustrial control. Biases from the control were inherited by its successor transient, in particular, sea surface temperatures (SSTs) at high latitudes were too warm and too much northward meridional heat transport occurred in the North Atlantic Ocean. SSTs in the Northern Hemisphere and Atlantic meridional heat transport increased in the transients initialized from the ocean/sea-ice simulation over a 40-50 year period to largely resemble the transient branched off the control simulation. Excessive absorbed solar radiation (ASR) is considered to be responsible for the SST biases as well as the formation of mode waters in lighter density classes than in reality both in the Southern Ocean and the North Atlantic Ocean subtropical gyre. These mode waters ventilate the thermocline waters to the north of their formation sites explaining the excessive Atlantic meridional heat transport and some of the warming in the high-latitude North Atlantic Ocean after multiple decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of Tropical Cyclone Wind Forcing on the Global Climate in a Fully Coupled Climate Model

Tropical cyclones (TCs) alter upper-ocean temperature and influence ocean heat content via enhanced turbulent mixing. A better understanding of the role of TCs within the climate system requires a fully coupled modeling framework, where TC-induced ocean responses feed back to the atmosphere and subsequently to the climate mean state and variability. Here, we investigate the impacts of TC wind forcing on the global ocean and the associated feedbacks within the climate system using the fully coupled Community Earth System Model version 1.3 (CESM1.3). Using the low-resolution version of CESM1.3 (1° atmosphere and ocean grid spacing) with no intrinsic TCs, we conduct a suite of sensitivity experiments by inserting TC winds extracted from a high-resolution (0.25° atmosphere grid spacing) TC-permitting simulation into the low-resolution model. Results from the low-resolution TC experiment are compared to a low-resolution control simulation to diagnose TCs’ impact. We found that the added TC winds can increase ocean heat content by affecting ocean vertical mixing, air–sea enthalpy fluxes, and cloud amount. Here, the added TCs can influence mean SST, precipitation, ocean subsurface temperature, and ocean mixed layer depth. We found a strengthening of the wind-driven subtropical cells and a weakening of the Atlantic meridional overturning circulation due to the changes of surface buoyancy fluxes. TCs in the model cause anomalous equatorward ocean heat convergence in the deep tropics and an increase of poleward ocean heat transport out of the subtropics. Our modeling results provide new insights into the multiscale interactions between TCs and the coupled climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Global Temperature Responses to Large Tropical Volcanic Eruptions in Paleo Data Assimilation Products and Climate Model Simulations Over the Last Millennium

Abstract Large volcanic eruptions are one of the dominant perturbations to global and regional atmospheric temperatures on timescales of years to decades. Discrepancies remain, however, in the estimated magnitude and persistence of the surface temperature cooling caused by volcanic eruptions, as characterized by paleoclimatic proxies and climate models. We investigate these discrepancies in the context of large tropical eruptions over the Last Millennium using two state‐of‐the‐art data assimilation products, the Paleo Hydrodynamics Data Assimilation product (PHYDA) and the Last Millennium Reanalysis (LMR), and simulations from the National Center for Atmospheric Research Community Earth System Model‐Last Millennium Ensemble (NCAR CESM‐LME). We find that PHYDA and LMR estimate mean global and hemispheric cooling that is similar in magnitude and persistence once effects from eruptions occurring in short succession are removed. The estimates also compare well to Northern‐Hemisphere reconstructions based solely or partially on tree‐ring density, which have been proposed as the most accurate proxy estimates of surface cooling due to volcanism. All proxy‐based estimates also agree well with the magnitude of the mean cooling simulated by the CESM‐LME. Differences remain, however, in the spatial patterns of the temperature responses in the PHYDA, LMR, and the CESM‐LME. The duration of cooling anomalies also persists for several years longer in the PHYDA and LMR relative to the CESM‐LME. Our results demonstrate progress in resolving discrepancies between proxy‐ and model‐based estimates of temperature responses to volcanism, but also indicate these estimates must be further reconciled to better characterize the risks of future volcanic eruptions.

Tejedor, E.↗
DOE OSTI2024

Impacts of spatial heterogeneity of anthropogenic aerosol emissions in a regionally refined global aerosol–climate model

Abstract. Emissions of anthropogenic aerosol and their precursors are often prescribed in global aerosol models. Most of these emissions are spatially heterogeneous at model grid scales. When remapped from low-resolution data, the spatial heterogeneity in emissions can be lost, leading to large errors in the simulation. It can also cause the conservation problem if non-conservative remapping is used. The default anthropogenic emission treatment in the Energy Exascale Earth System Model (E3SM) is subject to both problems. In this study, we introduce a revised emission treatment for the E3SM Atmosphere Model (EAM) that ensures conservation of mass fluxes and preserves the original emission heterogeneity at the model-resolved grid scale. We assess the error estimates associated with the default emission treatment and the impact of improved heterogeneity and mass conservation in both globally uniform standard-resolution (∼ 165 km) and regionally refined high-resolution (∼ 42 km) simulations. The default treatment incurs significant errors near the surface, particularly over sharp emission gradient zones. Much larger errors are observed in high-resolution simulations. It substantially underestimates the aerosol burden, surface concentration, and aerosol sources over highly polluted regions, while it overestimates these quantities over less-polluted adjacent areas. Large errors can persist at higher elevation for daily mean estimates, which can affect aerosol extinction profiles and aerosol optical depth (AOD). We find that the revised treatment significantly improves the accuracy of the aerosol emissions from surface and elevated sources near sharp spatial gradient regions, with significant improvement in the spatial heterogeneity and variability of simulated surface concentration in high-resolution simulations. In the next-generation E3SM running at convection-permitting scales where the resolved spatial heterogeneity is significantly increased, the revised emission treatment is expected to better represent the aerosol emissions as well as their lifecycle and impacts on climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluating Northern Hemisphere Growing Season Net Carbon Flux in Climate Models Using Aircraft Observations

Understanding terrestrial ecosystems and their response to anthropogenic climate change requires quantification of land-atmosphere carbon exchange. However, top-down and bottom-up estimates of large-scale land-atmosphere fluxes, including the northern extratropical growing season net flux (GSNF), show significant discrepancies. We developed a data-driven metric for the GSNF using atmospheric carbon dioxide concentration observations collected during the High-Performance Instrumented Airborne Platform for Environmental Research Pole-to-Pole Observations and Atmospheric Tomography Mission flight campaigns. This aircraft-derived metric is bias-corrected using three independent atmospheric inversion systems. We estimate the northern extratropical GSNF to be 5.7 ± 0.3 Pg C and use it to evaluate net biosphere productivity from the Coupled Model Intercomparison Project phase 5 and 6 (CMIP5 and CMIP6) models. While the model-to-model spread in the GSNF has decreased in the CMIP6 models relative to that of the CMIP5 models, there is still disagreement on the magnitude and timing of seasonal carbon uptake with most models underestimating the GSNF and overestimating the length of the growing season relative to the observations. We also use an emergent constraint approach to estimate annual northern extratropical gross primary productivity to be 56 ± 17 Pg C, heterotrophic respiration to be 25 ± 13 Pg C, and net primary productivity to be 28 ± 12 Pg C. The flux inferred from these aircraft observations provides an additional constraint on large-scale gross fluxes in prognostic Earth system models that may ultimately improve our ability to accurately predict carbon-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluating E3SM Global Storm‐Resolving Model Simulations of Deep Convection: Insights From DP‐SCREAM During TRACER

Global Storm-Resolving Models (GSRMs) are becoming increasingly vital for advancing climate modeling and improving the prediction of extreme weather events. Houston, a coastal region frequently affected by deep convective storms, offers an ideal setting to evaluate the ability of GSRMs to simulate deep convection. This study assesses the performance of the Doubly Periodic Simple Cloud-Resolving E3SM (Energy Exascale Earth System Model) Atmosphere Model (DP-SCREAM) using observations from the TRacking Aerosol Convection interactions ExpeRiment (TRACER) campaign. DP-SCREAM effectively reproduces the diurnal cycles of clouds and precipitation, demonstrating much greater skill than the E3SM single column model. The DP-SCREAM is demonstrated to be applicable to coastal regions, partially due to the forcing data sets already capturing the influence of breezes. DP-SCREAM also replicates biases persistent in the global version of SCREAM: the underrepresentation of boundary layer shallow clouds, a lack of mid-level congestus clouds, and the popcorn convection, characterized by small and disorganized convective cells generating the strongest precipitation. To investigate these issues, two sensitivity experiments were conducted: increasing the mixing length and scaling up the buoyancy flux within the Simplified Higher Order Closure scheme. Increasing the mixing length improved mid-level congestus representation and reduced unrealistic early morning fog occurrence. Enhancing buoyancy flux only marginally improved the bias of underproduced big convective cells. In conclusion, an additional resolution sensitivity test at 0.5 km grid spacing demonstrated that a refined horizontal resolution alone is insufficient to resolve these biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Microbial dormancy and its impacts on northern temperate and boreal terrestrial ecosystem carbon budget

A large amount of soil carbon in northern temperate and boreal regions could be emitted as greenhouse gases in a warming future. However, lacking detailed microbial processes such as microbial dormancy in current biogeochemistry models might have biased the quantification of the regional carbon dynamics. Here the effect of microbial dormancy was incorporated into a biogeochemistry model to improve the quantification for the last century and this century. Compared with the previous model without considering the microbial dormancy, the new model estimated the regional soils stored 75.9 Pg more C in the terrestrial ecosystems during the last century and will store 50.4 and 125.2 Pg more C under the RCP8.5 and RCP2.6 scenarios, respectively, in this century. This study highlights the importance of the representation of microbial dormancy in earth system models to adequately quantify the carbon dynamics in the northern temperate and boreal natural terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

PAVC Gridded 20m Alaska NGEE Tier3 PFTs v1.0

These 20-meter spatial resolution gridded products provide per-pixel fractional cover (%) of Next Generation Ecosystem Experiments (NGEE) Arctic Plant Functional Types (PFTs) Tier 3 across Alaska, north of the boreal treeline. The products were developed for the NGEE Arctic project, which is improving Arctic vegetation representation and parameterization of the E3SM Land Model. This dataset includes 8 files containing fractional cover for NGEE Tier 3 PFTs (https://data.ess-dive.lbl.gov/view/doi:10.15485/2529470): (1) bryophytes; (2) lichens; (3) non-vascular plants, i.e., the sum of lichens and bryophytes; (4) deciduous shrubs, (5) evergreen shrubs, (6) forbs, (7) graminoids, and a non-PFT class, (8) litter. Each pixel contains the percent cover (expressed as a fraction of total ground cover) that was predicted by random-forest regression models. The random-forest models were trained on cover data collected at 978 plots from 2010 to 2021, of which are archived in the Pan-Arctic Vegetation Cover (PAVC) database (https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2483557). The plot cover was linked to 20-meter spatial resolution, satellite-derived predictor variables: Sentinel-2 spectra and Sentinel-1 polarizations averaged over the 2019 growing season, as well as topographical features derived from ArcticDEM. Then, spatio-temporally anomalous plot data that introduced large variability to the regression outcomes were dropped using the Cook’s distance outlier detection method, and the models were re-created using high-quality plots and their associated satellite derived explanatory variables per each PFT. The correlations between plot-observed and satellite-derived fractional cover for all PFTs were well correlated (R2 = 0.69–0.95 and 0.5 for litter) and had low RMSE bias (0.02–0.11). This research was performed as a part of the NGEE Arctic project. The NGEE Arctic project was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Understanding the Drivers of Atlantic Multidecadal Variability using a Stochastic Model Hierarchy

The relative importance of ocean and atmospheric dynamics in generating Atlantic Multidecadal Variability (AMV) remains an open question. Comparisons between climate models with SLAB and fully-dynamic (FULL) ocean components are often used to explore this question, but cannot reveal how individual ocean processes generate these differences. We build a hierarchy of physically interpretable stochastic models to investigate the contribution of two upper-ocean processes to AMV: the role of seasonal variation and mixed-layer entrainment. This interpretability arises from the stochastic model’s simplified representation of sea surface temperature (SST), considering only the local upper ocean response to white-noise atmospheric forcing and its impact on surface heat exchange. We focus on understanding differences between SLAB and FULL non-eddy resolving pre-industrial control simulations of the Community Earth System Model 1 (CESM), and estimate the stochastic model parameters from each respective simulation. Despite its simplicity, the stochastic model reproduces temporal characteristics of SST variability in the SPG, including reemergence, seasonal-to-interannual persistence and power spectra. Furthermore, unrealistically persistent SST of the CESM-SLAB ocean simulation is reproduced in the equivalent stochastic model configuration where the mixed-layer depth (MLD) is constant. The stochastic model also reveals that vertical entrainment primarily damps SST variability, thus explaining why SLAB exhibits larger SST variance than FULL. Here, the stochastic model driven by temporally stochastic, spatially coherent forcing patterns reproduces the canonical AMV pattern. However, the amplitude of low-frequency variability remains underestimated, suggesting a role for ocean dynamics beyond entrainment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Limited increases in Arctic offshore oil and gas production with climate change and the implications for energy markets

Climate change impacts on sea ice thickness is opening access to offshore Arctic resources. The degree to which these resources are exploited will depend on sea-ice conditions, technology costs, international energy markets, and the regulatory environment. We use an integrated human-Earth system model, GCAM, to explore the effects of spatial–temporal patterns of sea-ice loss under climate change on future Arctic offshore oil and gas extraction, considering interactions with global energy markets and emission reduction scenarios. We find that under SSP5, a “fossil-fueled development” scenario, the effects of sea-ice loss are larger for Arctic offshore oil production than gas. Under SSP5, future extraction of Arctic offshore oil and gas through 2100 adds roughly 0.8–2.6 EJ/year to oil and gas markets but does not have large impacts on global oil and gas markets. Surprisingly, a low-carbon scenario results in greater Arctic offshore oil production to offset the more emissions-intensive unconventional oil production.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Evaluating the impact of peat soils and snow schemes on simulated active layer thickness at pan-Arctic permafrost sites

Abstract Permafrost stability is significantly influenced by the thermal buffering effects of snow and active-layer peat soils. In the warm season, peat soils act as a barrier to downward heat transfer mainly due to their low thermal conductivity. In the cold season, the snowpack serves as a thermal insulator, retarding the release of heat from the soil to the atmosphere. Currently, many global land models overestimate permafrost soil temperature and active layer thickness (ALT), partially due to inaccurate representations of soil organic matter (SOM) density profiles and snow thermal insulation. In this study, we evaluated the impacts of SOM and snow schemes on ALT simulations at pan-Arctic permafrost sites using the Energy Exascale Earth System Model (E3SM) land model (ELM). We conducted simulations at the Circumpolar Active Layer Monitoring (CALM) sites across the pan-Arctic domain. We improved ELM-simulated site-level ALT using a knowledge-based hierarchical optimization procedure and examined the effects of precipitation-phase partitioning methods (PPMs), snow compaction schemes, and snow thermal conductivity schemes on simulated snow depth, soil temperature, ALT, and CO 2 fluxes. Results showed that the optimized ELM significantly improved agreement with observed ALT (e.g. RMSE decreased from 0.83 m to 0.15 m). Our sensitivity analysis revealed that snow-related schemes significantly impact simulated snow thermal insulation levels, soil temperature, and ALT. For example, one of the commonly used snow thermal conductivity schemes (quadratic Sturm or SturmQua) generally produced warmer soil temperatures and larger ALT compared to the other two tested schemes. The SturmQua scheme also amplified the model’s sensitivity to PPMs and predicted deeper ALTs than the other two snow schemes under both current and future climates. The study highlights the importance of accurately representing snow-related processes and peat soils in land models to enhance permafrost dynamics simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Rapid SACR Observations of Convection at Bankhead National Forest (RAPID) Field Campaign Report

Improving our representation of convective cell processes requires better quantification of convective clouds throughout their entire life cycle. This includes gaining a clearer understanding of the controls on key convective cloud properties, such as updraft intensity, particle size distributions, rainfall rates, and hydrometeor species. Our inability to improve convective cloud process modeling stems, in part, from a limited understanding of convective cell properties, particularly given how rapidly these storms evolve. This lack of detailed observations in the most intense and organized convective storms is especially significant, as large errors remain in representing these clouds, which are critical for severe weather prediction and Earth system model performance. Cloud and precipitation radars are essential tools for studying cloud microphysics and dynamics, particularly in deeper convective clouds. The recent U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s third Mobile Facility (AMF3) Bankhead National Forest (BNF) deployment provides a unique opportunity to investigate important land–atmosphere interactions, as well as the environmental controls on deep convective cloud processes, in a location favorable for frequent convection. We operate the X/Ka-band Scanning ARM Cloud Radar (X/Ka SACR) using a scan strategy optimized to capture these rapidly evolving clouds and their properties, thereby improving studies of deep convective cloud processes. This effort is strengthened by a complementary and coordinated partnership with ongoing university and multi-agency radar activities collocated in north Alabama—a unique opportunity to examine clouds and precipitation from a lifetime-centric perspective.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Wetter California Projected by CMIP6 Models With Observational Constraints Under a High GHG Emission Scenario

As the world's fifth-largest economy entity, California (CA) is vulnerable to climate changes especially the magnitude of winter (wet season) precipitation that is closely linked to regional drought severity, vegetation growth, and wildfire activities. However, the fate of CA precipitation in the future remains highly uncertain, for example, the state-of-the-art Earth system models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) projected –3% to +42% and –27% to +63% precipitation changes in northern and central-southern CA, respectively, under the high greenhouse gas (GHG) emission scenario (SSP585). In this work, we applied the Pareto optimality concept and used observed teleconnection patterns (causal relationships between ocean regions and CA precipitation) to mechanistically constrain CMIP6 projected CA precipitations. We estimated that precipitation will robustly increase by 0.4–1.3 mm d –1 (10–34%) over northern CA and increase by 0.1–0.5 mm d –1 (7–32%) over central-southern CA by the end of the 21st century compared with present-day. Up to 71% of ESM projection uncertainties were reduced mainly due to the strong and consistent causal relationship between North American west coast sea level pressure and CA precipitation in both observations and CMIP6 models. Our results suggest that teleconnection patterns are powerful mechanistic constraints that can help explain and reduce uncertainties in ESM projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Assessment of Storm-Associated Precipitation and Its Extremes Using Observational Data Sets and Climate Model Short-Range Hindcasts

Heavy precipitation, often associated with weather phenomena such as tropical cyclones, extratropical cyclones (ETCs), atmospheric rivers (ARs), and mesoscale convective systems (MCSs), can cause significant socio-economic loss. Here, in this study, we apply atmospheric feature trackers to quantify the contributions of these storm types in observational data sets and climate model short-range hindcasts. We generate a global hourly storm data set at 0.25° spatial resolution covering 2006–2020, based on the tracking results from TempestExtremes and Python FLEXible object TRacKeR. Our analyses show that these four storm types account for 67% of global annual mean precipitation and 82% of top 1% precipitation extremes, with MCSs mainly over the tropics, and ARs and ETCs over the midlatitudes. The percentage of precipitation contributions from these storms also show strong seasonality over many geographical locations. We further apply the tracking results to the Energy Exascale Earth System Model (E3SM) short-range hindcasts and evaluate how well these storms are simulated. The evaluation show that E3SM, with ∼1° resolution, significantly underestimates storm-associated precipitation totals and extremes, especially for MCSs in the tropics. Our analysis also suggests that model fails to capture the correct mean diurnal phases and amplitude of MCS precipitation. This phenomenon-based approach provides a better understanding of precipitation characteristics and can lead to enhanced model evaluation by revealing underlying problems in model physics related to precipitation processes associated with the heavy-precipitating storms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Investigating Resiliency of Transportation Network under Targeted and Potential Climate Change Disruptions

Ensuring robustness and resilience in intermodal transportation systems is essential for the continuity and reliability of global logistics. These systems are vulnerable to various disruptions, including natural disasters and technical failures. Despite significant research on freight transportation resilience, investigating the robustness of the system after targeted and climate-change-driven disruption remains a crucial challenge. Drawing on network science methodologies, this study models the interdependencies within the rail and water transport networks and simulates different disruption scenarios to evaluate system responses. Here, we use the data from the U.S. Department of Energy Volpe Center for network topology and tonnage projections. The proposed framework quantifies deliberate, stochastic, and climate-driven infrastructure failure, using higher resolution downscaled multiple Earth System Models’ simulations from Coupled Model Intercomparison Project Phase version 6. We show that the disruptions of a few nodes could have a larger impact on the total tonnage of freight transport than on network topology. For example, the removal of targeted 20 nodes can bring the total tonnage carrying capacity to 30% with about 75% of the rail freight network intact. This research advances the theoretical understanding of transportation resilience and provides practical applications for infrastructure managers and policymakers. By implementing these strategies, stakeholders and policymakers can better prepare for and respond to unexpected disruptions, ensuring sustained operational efficiency in transportation networks.

Climate Change Disruptions↗
DOE OSTI2022

CondiDiag1.0: a flexible online diagnostic tool for conditional sampling and budget analysis in the E3SM atmosphere model (EAM)

Abstract. Numerical models used in weather and climate prediction take into account a comprehensive set of atmospheric processes (i.e., phenomena) such as the resolved and unresolved fluid dynamics, radiative transfer, cloud and aerosol life cycles, and mass or energy exchanges with the Earth's surface. In order to identify model deficiencies and improve predictive skills, it is important to obtain process-level understanding of the interactions between different processes. Conditional sampling and budget analysis are powerful tools for process-oriented model evaluation, but they often require tedious ad hoc coding and large amounts of instantaneous model output, resulting in inefficient use of human and computing resources. This paper presents an online diagnostic tool that addresses this challenge by monitoring model variables in a generic manner as they evolve within the time integration cycle. The tool is convenient to use. It allows users to select sampling conditions and specify monitored variables at run time. Both the evolving values of the model variables and their increments caused by different atmospheric processes can be monitored and archived. Online calculation of vertical integrals is also supported. Multiple sampling conditions can be monitored in a single simulation in combination with unconditional sampling. The paper explains in detail the design and implementation of the tool in the Energy Exascale Earth System Model (E3SM) version 1. The usage is demonstrated through three examples: a global budget analysis of dust aerosol mass concentration, a composite analysis of sea salt emission and its dependency on surface wind speed, and a conditionally sampled relative humidity budget. The tool is expected to be easily portable to closely related atmospheric models that use the same or similar data structures and time integration methods.

58 GEOSCIENCES↗
DOE OSTI2023

CORPSE model with litter decomposition parameters derived from the LIDET dataset

This is a version of the CORPSE model (Carbon, Organisms, Rhizosphere and Protection in the Soil Environment, Sulman et al. 2014) that uses litter decomposition parameters derived from a modified Monte Carlo simulation using the LIDET litter decomposition dataset (Long-term Intersite Decomposition Experiment Team, Harmon 2013). The code also includes the Baseline parameters, and the eight other best parameter sets identified in a modified Monte Carlo simulation. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of files:The folder "Input Files" contains one folder for each LIDET site with data necessary to run the model. Note that "(site)" in the filenames below indicates where the LIDET site code appears (see Table 1 for site codes). Data streams include: CORPSE_full_spinup_litter.csv, CORPSE_full_spinup_rhizo.csv, CORPSE_full_spinup_bulk.csv, litterbag_init_100g_6spp.csv: initial C and N (kg C or N/m2) pool values for each soil layer, the litterbag_init_100_6spp.csv file is for the litterbag layer and is the same file for all sites. All initial C and N files have the same columns (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 soilT (site) DOY274start.csv: Average daily soil temperature (oC) interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). soilT (site) DOY274start.csv: Average daily soil volumetric water content (VWC) scalar interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). litter production.csv: Average daily litter production values for each site, data sources listed in Table S3 of related publication. litter (site) CN.csv: C:N ratio for each species from LIDET dataset (Table 2, Harmon 2013). (site).csv: Table indicating number of observations for each species decomposed at each site. Instructions: Save the model code ("CORPSE_LIDET.R") and "Input Files" folder in the same folder. Also make a folder for the model output (e.g., "results_Baseline") in the same folder. Set the working directory (setwd) in the model code to the folder with the files saved in step #1. Select the parameter set to use for the litter and litterbag compartments, comment out all other parameter sets. Run code. Output will be saved in the folder made in step 1. Output destination can be changed as necessary in code section called "Running the model." Table 1 LIDET sites and site codes used in model files. Site Code - Site AND - H.J. Andrews Experimental Forest BNZ - Bonanza Creek Experimental Forest BSF - Blodgett Research Forest CDR - Cedar Creek Natural History Area CPR - Central Plains Experimental Range HBR - Hubbard Brook Experimental Forest HFR - Harvard Forest JUN - Juneau KBS - Kellogg Biological Station KNZ - Konza Prairie Research Natural Area NWT - Niwot Ridge/Green Lakes Valley OLY - Olympic National Park OLY Conifer forest SEV - Sevilleta National Wildlife Refuge SMR - Santa Margarita Ecological Reserve UFL - University of Florida VCR - Virginia Coast Reserve Table 2 LIDET species and species codes used in model files (6 common species). Species - Species Code Sugar maple (Acer saccharum) - ACSA Drypetes (Drypetes glauca) - DRGL Red pine (Pinus resinosa) - PIRE Chestnut oak (Quercus prinus) - QUPR Western redcedar (Thuja plicata) - THPL Wheat (Triticum aestivum) - TRAE References:Bonan, G. B., Hartman, M. D., Parton, W. J., & Wieder, W. R. (2013). Evaluating litter decomposition in earth system models with long-term litterbag experiments: an example using the Community Land Model version 4 (CLM4). Global Change Biology, 19(3), 957-974. https://doi.org/https://doi.org/10.1111/gcb.12031 Harmon, M. (2013). LTER Intersite Fine Litter Decomposition Experiment (LIDET), 1990 to 2002. Long-Term Ecological Research. Forest Science Data Bank, Corvallis, OR. [Data set]. Accessed http://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023. https://doi.org/10.6073/pasta/f35f56bea52d78b6a1ecf1952b4889c5. Sulman, B. N., Phillips, R. P., Oishi, A. C., Shevliakova, E., & Pacala, S. W. (2014). Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO2. Nature Climate Change, 4, 1099 - 1102. https://doi.org/10.1038/nclimate2436

Juice, Stephanie↗
DOE OSTI2023

Properties and Controlling Processes of Aerosol and Cloud Condensation Nuclei in Marine Boundary Layer Over Eastern North Atlantic (Final Report)

Remote marine low cloud systems have large spatial and temporal coverages. Because of their relatively low optical thickness and background aerosol concentrations, marine low clouds are particularly susceptible to perturbations in aerosols associated with anthropogenic emissions. Recent studies show that a large fraction of the global aerosol indirect forcing can be attributed to changes in marine low clouds, despite their relatively long distance from most anthropogenic sources. Currently, there remain large uncertainties in the magnitude of the global aerosol indirect radiative forcing. One major contribution to this large uncertainty derives from poor understanding of the aerosols under natural conditions, the perturbation by anthropogenic emissions, and the processes that drive them. In this project, we took advantage of the comprehensive, multi-dimensional characterization of meteorological parameters, trace gases, aerosol, cloud, and drizzle fields during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) and the long-term measurements at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site on Graciosa Island. By synergistically combining in situ measurements, remote sensing retrievals, and model simulations, we quantified aerosol properties in the ENA, their vertical profiles, and variations with season, synoptic state, and airmass. We investigated key aerosol processes, including entrainment of free troposphere aerosol into marine boundary layer, new particle formation, growth of Aitken mode particles, and the influences of these processes on the cloud condensation nuclei (CCN) population in marine boundary layer. The representation of key aerosol processes in the Energy Exascale Earth System Model (E3SM) was evaluated using the observations. This project has led to 26 publications, which are listed at the end of this report.

54 ENVIRONMENTAL SCIENCES↗