Search NASA⌕ Search

SEARCH · Search NASA

Results for “default”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Numerically Relevant Timescales in the MG2 Microphysics Model

Climate models rely on parameterizations of a variety of processes in the atmospheric physics, but a common concern is that the temporal resolution is too coarse to consistently resolve the behavior that individual parameterizations are designed to capture. This study examines timescales numerically derived from the Morrison-Gettelman (MG2) microphysics as implemented within the Energy Exascale Earth System Model, version 1 (E3SMv1). Numerically relevant timescales in MG2 are derived by computing the eigenspectrum of its Jacobian. These timescales are found to often be smaller than the default 5 min time step used for MG2. The fast timescales are then heuristically connected to individual microphysics processes. By substepping a few particular rain processes within MG2, the time discretization error for those processes was considerably reduced with minimal additional expense to the overall microphysics. While this improvement has a substantial effect on the target processes and on the vertical distribution of stratiform-derived rain within E3SMv1, the overall model climate is found to not be sensitive to MG2 time step. Here, we hypothesize that this is because the surface climate does not depend strongly on certain process rates, especially MG2's rain evaporation rate.

58 GEOSCIENCES↗
DOE OSTI2020

Understanding cloud droplet spectral dispersion effect using empirical and semi-analytical parameterizations in NCAR CAM5.3

Five parameterizations of cloud droplet spectral shape are implemented in a global climate model to investigate the dispersion effect and aerosol indirect effect (AIE). We design a series of experiments by modifying the microphysical cloud scheme of NCAR CAM5.3 (Community Atmospheric Model Version 5.3). We employ four empirical (Martin94, RLiu03, PengL03, Liu08) and one semi-analytical (LiuLi15) expressions for cloud droplet spectral shape parameters. Analysis focuses on the instantaneous differences in the simulated cloud microphysical properties and the comparison between model output and satellite data. The results show that RLiu03, PengL03 and LiuLi15 produce wider droplet spectrum and faster autoconversion rate, but Liu08 has a narrower droplet spectrum and slower autoconversion rate than the default parameterization (Martin94) in CAM5.3. Global dispersion effects caused by the five parameterizations modify the aerosol indirect effect by -10% (counteract) to 13% (strengthen). The simulated AIEs and dispersion effects exhibit noticeably spatial inhomogeneity. In the sensitive regions of AIE (Southeast Asia, North Pacific and west coast of South America), we decompose the response of shortwave cloud forcing to the change in droplet number for analysis. The varying dispersion effects can be explained by different responses of cloud properties in different spectral parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulated Precipitation Diurnal Variation With a Deep Convective Closure Subject to Shallow Convection in Community Atmosphere Model Version 5 Coupled With CLUBB

In order to improve the physical consistency between shallow and deep convection, we modify the deep convective closure in the Community Atmosphere Model version 5 (CAM5) coupled with a third-order turbulence closure parameterization (i.e., Cloud Layers Unified by Binormals; CLUBB). The revised closure reserves a portion of the total convective available potential energy for shallow convection via utilizing the heating and moistening profiles from CLUBB to distribute moisture and energy between shallow and deep convection. Simulations at two resolutions (i.e. 2° and 0.5°, respectively) are conducted to investigate the impacts of convective closure on the simulated precipitation diurnal variations. Results from low-resolution simulations show that the revised closure suppresses deep convection until the lower troposphere is sufficiently moistened by shallow convection, which improves the precipitation diurnal variation simulations compared to the default closure, with the precipitation diurnal peak over tropical lands delayed from 12LST to 19LST. The revised closure better simulates the diurnal variations for precipitation over the Asian monsoon region, such as the delayed precipitation onset, but still fails to well capture the nighttime peak for precipitation there. This deficiency is alleviated to some extent when applying the revised closure in highresolution simulations, but nighttime precipitation is still underestimated probably because key processes responsible for nighttime convection are missing. Altogether, our results indicate that establishing the consistency between shallow and deep convection is critical for the precipitation diurnal cycle simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Integrating Arctic Plant Functional Types in a Land Surface Model Using Above‐ and Belowground Field Observations

Abstract Accurate simulations of high‐latitude ecosystems are critical for confident Earth system model (ESM) projections of carbon cycle feedbacks to global climate change. Land surface model components of ESMs, including the E3SM Land Model (ELM), simulate vegetation growth and ecosystem responses to changing climate and atmospheric CO 2 concentrations by grouping heterogeneous vegetation into like sets of plant functional types (PFTs). Many such models represent high‐latitude vegetation using only two PFTs (shrub and grass), thereby missing the diversity of vegetation growth forms and functional traits in the Arctic. Here, we use field observations of biomass and leaf traits across a gradient of plant communities on the Seward Peninsula in northwest Alaska to replace the original ELM configuration for the first time with nine Arctic‐specific PFTs. The newly developed PFTs include: (1) nonvascular mosses and lichens, (2) deciduous and evergreen shrubs of various height classes, including an alder PFT, (3) graminoids, and (4) forbs. Improvements relative to the original model configuration included greater belowground biomass allocation, persistent fine roots and rhizomes of nonwoody plants, and better representation of variability in total plant biomass across sites with varying plant communities and depth to bedrock. Simulations through 2100 using the RCP8.5 climate scenario and constant PFT fractional areas showed alder‐dominated plant communities gaining more biomass and lichen‐dominated communities gaining less biomass compared to default PFTs. Our results highlight how representing the diversity of arctic vegetation and confronting models with measurements from varied plant communities improves the representation of arctic vegetation in terrestrial ecosystem models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Effect of the Grell-Freitas Deep Convection Scheme in Quasi-Uniform and Variable-Resolution Aquaplanet CAM Simulations

The Grell-Freitas deep convection scheme (GF) has been implemented in the Community Atmosphere Model (CAM) with the nonhydrostatic Model for Prediction Across Scales (MPAS) dynamical core. Aquaplanet simulations have been performed with a global quasi-uniform resolution mesh with ~120 km grid spacing and a variable resolution mesh employing a circular refined region with ~30 km grid spacing centering at the equator and ~120 km elsewhere. GF produces overall patterns of the general circulations comparable to those of other aquaplanet experiments, including those produced using the default Zhang-McFarlane convection scheme (ZM) in CAM Version 5.4. GF alleviates the overestimated occurrence of weak intensity precipitation characteristic of ZM due to GF’s different closure methods and deep-convection trigger conditions. Compared to ZM, the GF characteristics of the equatorial Kelvin waves are closer to observations. GF also better simulates precipitation at variable resolution, with no unrealistic wave distortion by the refinement. Overall, GF shows less sensitivity of convective precipitation and heating/moistening tendencies to changes in horizontal resolution compared to ZM. However, vertically-overlapped total cloud fraction simulated with GF in the tropics is more sensitive to the resolution than ZM due to the larger resolution-sensitivity of high-level cloud fraction in the tropics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Optimizing Carbon Cycle Parameters Drastically Improves Terrestrial Biosphere Model Underestimates of Dryland Mean Net CO 2 Flux and its Inter-Annual Variability

Drylands occupy ~40% of the land surface and are thought to dominate global carbon (C) cycle inter-annual variability (IAV). Therefore, it is imperative that global terrestrial biosphere models (TBMs), which form the land component of IPCC earth system models, are able to accurately simulate dryland vegetation and biogeochemical processes. However, compared to more mesic ecosystems, TBMs have not been widely tested or optimized using in situ dryland CO 2 fluxes. Here, for this work, we address this gap using a Bayesian data assimilation system and 89 site-years of daily net ecosystem exchange (NEE) data from 12 southwest US Ameriflux sites to optimize the C cycle parameters of the ORCHIDEE TBM. The sites span high elevation forest ecosystems, which are a mean sink of C, and low elevation shrub and grass ecosystems that are either a mean C sink or “pivot” between an annual C sink and source. We find that using the default (prior) model parameters drastically underestimates both the mean annual NEE at the forested mean C sink sites and the NEE IAV across all sites. Our analysis demonstrated that optimizing phenology parameters are particularly useful in improving the model's ability to capture both the magnitude and sign of the NEE IAV. At the forest sites, optimizing C allocation, respiration, and biomass and soil C turnover parameters reduces the underestimate in simulated mean annual NEE. Our study demonstrates that all TBMs need to be calibrated for dryland ecosystems before they are used to determine dryland contributions to global C cycle variability and long-term carbon-climate feedbacks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Model Evaluation and Intercomparison of Marine Warm Low Cloud Fractions With Neural Network Ensembles

Abstract Low cloud fractions (LCFs) and meteorological factors (MFs) over an oceanic region containing multiple cloud regimes are examined for three data sets: one Energy Exascale Earth System Model (E3SM) simulation with the default 72‐layer vertical grid (E3SM72), another one with 8‐times vertical resolution via the Framework for Improvement by Vertical Enhancement (E3SM 8), and one with MFs from ERA5 reanalysis and LCFs from the CERES SSF product (ERA5‐SSF). Neural networks (NNs) are trained to capture the relationship between MFs and LCF and to select the best‐performing MF subsets for predicting LCF. NN ensembles are used to (a) confirm the performance of selected MF subsets, (b) to serve as proxy models for each data set to predict LCFs for MFs from all data sets, and (c) to classify MFs into those in shared and uniquely occupied MF subspaces. Overall, E3SM72 and E3SM 8 have large fractions of MFs in shared MF subspace, but less so near the Californian and Peruvian stratocumulus decks. E3SM 8 and ERA5 have small fractions of MFs in shared MF subspace but greater than E3SM72 and ERA5, especially in the Southeast Pacific. The differences in LCFs between three pairs of data sets are decomposed into those associated with the differences in the LCF‐MF relationship and those involving different MFs. Given the same MFs, LCFs produced by E3SM 8 are greater than those produced by E3SM72 but are still different from those in ERA5‐SSF. In general, the shift in MFs dominates the difference in the LCFs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Impact of a New Cloud Microphysics Parameterization on the Simulations of Mesoscale Convective Systems in E3SM

Abstract Mesoscale convective systems (MCSs) are one of the most climatically significant forms of convection because of their large role in water and energy cycles. The mesoscale features associated with MCS are difficult to represent in climate models because the relevant dynamics and physics are absent or poorly represented with coarse model resolution (∼100 km). Using a regionally refined model (RRM) with 0.25° grid spacing embedded in the Energy Exascale Earth System Model (E3SM), we explore the impact of cloud microphysics parameterizations on the simulation of precipitation, particularly MCS precipitation over the contiguous United States. The Predicted Particle Properties (P3) cloud microphysics scheme has been modified and implemented into E3SM to overcome the limitations of the default Morrison and Gettelman (MG2) scheme in which rimed precipitating ice particles (graupel/hail) are absent and frozen particles are artificially partitioned into cloud ice and snow. We show that P3 improves the simulation of precipitation statistics including frequency distribution compared with MG2 with a limited effect on the diurnal cycle. P3 predicts higher hourly rain rates, resulting in 20% more MCSs and a higher total MCS precipitation (4.4%) compared to MG2, agreeing better with observations. The improvements with P3 mainly result from improved representations of ice microphysics, which not only produces higher rain rates through melting but also leads to a stronger large‐scale ascending motion by releasing more latent heating. This study suggests that improving microphysics parameterization is important for simulating MCS precipitation as future climate model resolutions continue to increase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Implementation and Evaluation of a Unified Turbulence Parameterization Throughout the Canopy and Roughness Sublayer in Noah‐MP Snow Simulations

Abstract The Noah‐MP land surface model (LSM) relies on the Monin‐Obukhov (M‐O) Similarity Theory (MOST) to calculate land‐atmosphere exchanges of water, energy, and momentum fluxes. However, MOST flux‐profile relationships neglect canopy‐induced turbulence in the roughness sublayer (RSL) and parameterize within‐canopy turbulence in an ad hoc manner. We implement a new physics scheme (M‐O‐RSL) into Noah‐MP that explicitly parameterizes turbulence in RSL. We compare Noah‐MP simulations employing the M‐O‐RSL scheme (M‐O‐RSL simulations) and the default M‐O scheme (M‐O simulations) against observations obtained from 647 Snow Telemetry (SNOTEL) stations and two AmeriFlux stations in the western United States. M‐O‐RSL simulations of snow water equivalent (SWE) outperform M‐O simulations over 64% and 69% of SNOTEL sites in terms of root‐mean‐square‐error (RMSE) and correlation, respectively. The largest improvements in skill for M‐O‐RSL occur over closed shrubland sites, and the largest degradations in skill occur over deciduous broadleaf forest sites. Differences between M‐O and M‐O‐RSL simulated snowpack are primarily attributable to differences in aerodynamic conductance for heat underneath the canopy top, which modulates sensible heat flux. Differences between M‐O and M‐O‐RSL within‐canopy and below‐canopy sensible heat fluxes affect the amount of heat transported into snowpack and hence change snowmelt when temperatures are close to or above the melting point. The surface energy budget analysis over two AmeriFlux stations shows that differences between M‐O and M‐O‐RSL simulations can be smaller than other model biases (e.g., surface albedo). We intend for the M‐O‐RSL physics scheme to improve performance and uncertainty estimates in weather and hydrological applications that rely on Noah‐MP.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

The Global Atmosphere‐aerosol Model ICON‐A‐HAM2.3–Initial Model Evaluation and Effects of Radiation Balance Tuning on Aerosol Optical Thickness

Abstract The Hamburg Aerosol Module version 2.3 (HAM2.3) from the ECHAM6.3‐HAM2.3 global atmosphere‐aerosol model is coupled to the recently developed icosahedral nonhydrostatic ICON‐A (icon‐aes‐1.3.00) global atmosphere model to yield the new ICON‐A‐HAM2.3 atmosphere‐aerosol model. The ICON‐A and ECHAM6.3 host models use different dynamical cores, parameterizations of vertical mixing due to sub‐grid scale turbulence, and parameter settings for radiation balance tuning. Here, we study the role of the different host models for simulated aerosol optical thickness (AOT) and evaluate impacts of using HAM2.3 and the ECHAM6‐HAM2.3 two‐moment cloud microphysics scheme on several meteorological variables. Sensitivity runs show that a positive AOT bias over the subtropical oceans is remedied in ICON‐A‐HAM2.3 because of a different default setting of a parameter in the moist convection parameterization of the host models. The global mean AOT is biased low compared to MODIS satellite instrument retrievals in ICON‐A‐HAM2.3 and ECHAM6.3‐HAM2.3, but the bias is larger in ICON‐A‐HAM2.3 because negative AOT biases over the Amazon, the African rain forest, and the northern Indian Ocean are no longer compensated by high biases over the sub‐tropical oceans. ICON‐A‐HAM2.3 shows a moderate improvement with respect to AOT observations at AERONET sites. A multivariable bias score combining biases of several meteorological variables into a single number is larger in ICON‐A‐HAM2.3 compared to standard ICON‐A and standard ECHAM6.3. In the tropics, this multivariable bias is of similar magnitude in ICON‐A‐HAM2.3 and in ECHAM6.3‐HAM2.3. In the extra‐tropics, a smaller multivariable bias is found for ICON‐A‐HAM2.3 than for ECHAM6.3‐HAM2.3.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Development of Land‐River Two‐Way Hydrologic Coupling for Floodplain Inundation in the Energy Exascale Earth System Model

Abstract Floodplain inundation links river and land systems through significant water, sediment, and nutrient exchanges. However, these two‐way interactions between land and river are currently missing in most Earth System Models. In this study, we introduced the two‐way hydrological coupling between the land component, E3SM Land Model, and the river component, Model for Scale Adaptive River Transport, in Energy Exascale Earth System Model (E3SM) to study the impacts of floodplain inundation on land and river processes. We calibrated the river channel geometry and developed a new data‐driven inundation scheme to improve the simulation of inundation dynamics in E3SM. The new inundation scheme captures 96% of the spatial variation of inundation area in a satellite inundation product at global scale, in contrast with 7% when the default inundation scheme of E3SM was used. Global simulations including the new inundation scheme performed at resolution with and without land‐river two‐way coupling was used to quantify the impact of coupling. Comparisons show that two‐way coupling modifies the water and energy cycle in 20% of the global land cells. Specifically, riverine inundation is reduced by two‐way coupling, but inland inundation is intensified. Wetter periods are more impacted by the two‐way coupling at the global scale, while regions with different climates exhibit different sensitivities. The two‐way exchange of water between the land and river components of E3SM provides the foundation for enabling two‐way coupling of land‐river sediment and biogeochemical fluxes. These capabilities will be used to improve understanding of the interactions between water and biogeochemical cycles and their response to human perturbations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improved Dust Representation and Impacts on Dust Transport and Radiative Effect in CAM5

Dust transport and spatial distribution are poorly represented in current global climate models (GCMs) including the Community Atmosphere Model version 5 (CAM5). Particularly, models lack explicit representation of super-coarse dust, which may have important implications for dust radiative forcing and impacts on biogeochemistry. A nine-mode version of the modal aerosol model (MAM9) has been developed to address these issues. In this new aerosol scheme, four dust modes have been designed to treat dust particles of sizes up to 20 μm. The MAM9-simulated results are compared with those from the default four-mode version of MAM (MAM4) and also with the in situ surface measurements of dust concentration and deposition flux, satellite-retrieved dust extinction profile, and in situ vertical measurements of dust concentrations from the NASA Atmosphere Tomography Mission (ATom). Overall, MAM9 improves the dust representation in remote regions while maintaining reasonably good results near the dust source regions. In addition, MAM9 reduces the fine dust burden and increases the coarse dust burden globally. The increased coarse dust burden has slightly increased the dust direct radiative effect by 0.01 W m –2 while it enhanced dust indirect radiative effect by 0.36 W m –2 , globally.

58 GEOSCIENCES↗
DOE OSTI2022

Global Dust Cycle and Direct Radiative Effect in E3SM Version 1: Impact of Increasing Model Resolution

Quantification of dust aerosols in Earth System Models (ESMs) has important implications for water cycle and biogeochemistry studies. This study examines the global life cycle and direct radiative effects (DREs) of dust in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1), and the impact of increasing model resolution both horizontally and vertically. The default 1° E3SMv1 captures the spatial and temporal variability in the observed dust aerosol optical depth (DAOD) reasonably well, but overpredicts dust absorption in the shortwave (SW). Simulations underestimate the dust vertical and long-range transport, compared with the satellite dust extinction profiles. After updating dust refractive indices and correcting for a bias in partitioning size-segregated emissions, both SW cooling and longwave (LW) warming of dust simulated by E3SMv1 are increased and agree better with other recent studies. The estimated net dust DRE of –0.42 Wm –2 represents a stronger cooling effect than the observationally based estimate –0.2 Wm –2 (–0.48 to +0.2), due to a smaller LW warming. Constrained by a global mean DAOD, model sensitivity studies of increasing horizontal and vertical resolution show strong influences on the simulated global dust burden and lifetime primarily through the change of dust dry deposition rate; there are also remarkable differences in simulated spatial distributions of DAOD, DRE, and deposition fluxes. Thus, constraining the global DAOD is insufficient for accurate representation of dust climate effects, especially in transitioning to higher- or variable-resolution ESMs. Better observational constraints of dust vertical profiles, dry deposition, size, and LW properties are needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

On the Links Between Ice Nucleation, Cloud Phase, and Climate Sensitivity in CESM2

Ice nucleation in mixed-phase clouds has been identified as a critical factor in projections of future climate. Here we explore how this process influences climate sensitivity using the Community Earth System Model 2 (CESM2). We find that ice nucleation affects simulated cloud feedbacks over most regions and levels of the troposphere, not just extratropical low clouds. However, with present-day global mean cloud phase adjusted to replicate satellite retrievals, similar total cloud feedback is attained whether ice nucleation is simulated as aerosol-sensitive, insensitive, or absent. These model experiments all result in a strongly positive total cloud feedback, as in the default CESM2. A microphysics update from CESM1 to CESM2 had substantially weakened ice nucleation, due partly to a model issue. Our findings indicate that this update reduced global cloud phase bias, with CESM2's high climate sensitivity reflecting improved mixed-phase cloud representation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Effect of Physically Based Ice Radiative Processes on Greenland Ice Sheet Albedo and Surface Mass Balance in E3SM

Abstract A significant portion of surface melt on the Greenland Ice Sheet (GrIS) is due to dark ice regions in the ablation zone, where solar absorption is influenced by the physical properties of the ice, light absorbing constituents (LACs), and the overlying crustal surface or melt ponds. Earth system models (ESMs) typically prescribe the albedo of ice surfaces as a constant value in the visible and near‐infrared spectral regions. This work advances ESM ice radiative transfer modeling by (a) incorporating a physically based radiative transfer model (SNow, ICe and Aerosol Radiation model Adding‐Doubling Version 4; SNICAR‐ADv4) into the Energy Exascale Earth System Model (E3SM), (b) determining spatially and temporally varying bare ice physical properties over the GrIS ablation zone from satellite observations to inform SNICAR‐ADv4, and (c) assessing the impacts on simulated GrIS albedo and surface mass balance associated with modeling of more realistic bare ice albedo. GrIS‐wide bare ice albedo in E3SMv2 is overestimated by ∼4% in the visible and ∼7% in the near‐infrared wavelengths compared to the Moderate Resolution Imaging Spectroradiometer. Our bare ice physical property retrieval method found that LACs, ice crustal surfaces, and melt ponds reduce visible albedo by 30% in the bare ice region of the GrIS ablation zone. The realistic bare ice albedo reduces surface mass balance by ∼145 Gt, or 0.4 mm of sea‐level equivalent between 2000 and 2021 compared to the default E3SM. This work highlights the importance of simulating bare ice albedo accurately and realistically to improve our ability to quantify changes in the GrIS surface mass and radiative energy budgets.

58 GEOSCIENCES↗
DOE OSTI2024

Improving Stratocumulus Cloud Amounts in a 200‐m Resolution Multi‐Scale Modeling Framework Through Tuning of Its Interior Physics

Abstract High‐Resolution Multi‐scale Modeling Frameworks (HR)—global climate models that embed separate, convection‐resolving models with high enough resolution to resolve boundary layer eddies—have exciting potential for investigating low cloud feedback dynamics due to reduced parameterization and ability for multidecadal throughput on modern computing hardware. However low clouds in past HR have suffered a stubborn problem of over‐entrainment due to an uncontrolled source of mixing across the marine subtropical inversion manifesting as stratocumulus dim biases in present‐day climate, limiting their scientific utility. We report new results showing that this over‐entrainment can be partly offset by using hyperviscosity and cloud droplet sedimentation. Hyperviscosity damps small‐scale momentum fluctuations associated with the formulation of the momentum solver of the embedded large eddy simulation. By considering the sedimentation process adjacent to default one‐moment microphysics in HR, condensed phase particles can be removed from the entrainment zone, which further reduces entrainment efficiency. The result is an HR that can produce more low clouds with a higher liquid water path and a reduced stratocumulus dim bias. Associated improvements in the explicitly simulated sub‐cloud eddy spectrum are observed. We report these sensitivities in multi‐week tests and then explore their operational potential alongside microphysical retuning in decadal simulations at operational 1.5° exterior resolution. The result is a new HR having desired improvements in the baseline present‐day low cloud climatology, and a reduced global mean bias and root mean squared error of absorbed shortwave radiation. We suggest it should be promising for examining low cloud feedbacks with minimal approximation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Incorporating an Interactive Fire Plume-Rise Model in the DOE's Energy Exascale Earth System Model Version 1 (E3SMv1) and Examining Aerosol Radiative Effect

The vertical distribution of biomass burning aerosol (BBA) is important in regulating their impacts on weather and climate. The plume-rise process affects the injection height of BBA and interacts with the air parcel lifting and cloud processes. However, these processes are not represented in most global climate models. In this study, we replaced the fixed vertical profiles of monthly BBA emissions in the Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) with an interactive fire plume-rise model. The vertical distribution of BBA emissions was calculated as a function of ambient thermodynamic conditions from the host E3SMv1, with distributions of fire sizes and sensible heat fluxes derived from the observations. The maximum fire radiative power (FRP) technique was used to determine the fire size. Scaling-FRP technique is used to calculate the wildfire heat release. Daily BBA emission, superimposed with a fire diurnal cycle retrieved from the satellite observation, was included in model simulations. The model shows improved agreement with satellite retrievals and in situ measurement during the National Oceanic and Atmospheric Administration Wildfire Experiment for Cloud chemistry, Aerosol absorption, and Nitrogen campaign. The model-observation comparison demonstrates the importance of the plume-rise model and fire diurnal cycle assumption in determining the BBA fields. We also find that E3SMv1 with new features produces a larger carbonaceous aerosol burden, leading to 0.13 W m –2 warming at the top of atmosphere compared to the default E3SMv1. This highlights the importance of accurately representing the BBA injection height and indicates a no-linear nature in the BBA-induced radiative effect.

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↗