Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atmospheric Model Parameterizations”

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

Improving Predictability of Mixed-Phase Clouds and Aerosol Interactions in the Community Earth System Model (CESM) with ARM Measurements

The objective of this project is to improve the simulation and predictability of mixed-phase clouds and aerosol interactions in the Community Earth System Model (CESM) through comparisons with the DOE ARM observations. There are four major goals of the proposed study: (1) Improve the representation of ice microphysical processes in mixed-phase clouds; (2) Develop a long-term multi-sensor mixed-phase cloud dataset; (3) Test the performance of improved ice microphysics in CESM-CAM5 with the ARM data; and (4) Examine mixed-phase cloud microphysics-aerosol-dynamics-radiation interactions in CESM-CAM5. In this project, we have (1) Improved the representation of ice microphysical processes in mixed-phase clouds that include ice nucleation, ice depositional growth through the Wegener–Bergeron–Findeisen (WBF) process, feedbacks between cloud microphysics, dynamics and radiation, and provides links to aerosol properties in the Community Earth System Model (CESM)-Community Atmosphere Model version 5 (CAM5); (2) Developed a long-term multi-sensor mixed-phase cloud dataset. The mixed-phase cloud retrieval algorithms are improved by refining the phase determinations and ice number concentration retrievals. A multi-year mixed-phase cloud observation dataset (1999-2004; October 2013-February 2017) is regenerated at the ARM Barrow site with the improved algorithms. The dataset has been used in the global model evaluation to improve the model parameterizations; (3) Tested the performance of improved ice microphysics in CESM-CAM5 with the ARM data. CESM-CAM5 is run in the DOE Cloud-Associated Parameterizations Testbed (CAPT) and in the single-column model (SCM) mode to facilitate comparison with ARM observations. Multi-year CAPT simulations of seasonal variations of cloud microphysical properties, cloud liquid water path (LWP) and ice water path (IWP), cloud longwave and shortwave forcings, cloud occurrence frequency are evaluated against the ARM multi-sensor mixed-phase cloud retrievals; and (4) Examined mixed-phase cloud microphysics-aerosol-dynamics-radiation interactions in CESM-CAM5 that include the impacts of different ice nucleation mechanisms, sensitivity to dust IN concentration, treatment of the WBF process on mixed-phase cloud properties. Aerosol effect on mixed-phase clouds through the liquid phase (droplet activation) and ice phase processes (e.g., the glaciation indirect effect) is investigated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Data-driven emulation of modal aerosol microphysics via neural operator-based modeling

The complexity and the small characteristic scales of aerosol microphysical processes pose a big challenge for accurate and efficient Earth system simulations at regional and global scales. In this work, we construct and evaluate a surrogate model: the aerosol deep operator network (ADON), a physics-inspired dual-net architecture for emulating the aerosol microphysics parameterization suite in the version 2 of the Energy Earth System Model (E3SMv2). The current version of the surrogate model is trained on a dataset comprising 9.8 million samples obtained from a global E3SMv2 simulation with the horizontal resolution of about one degree under cloud-free conditions. Incorporating domain spatial and temporal coordinates, as well as principle components extracted from training data, the dual-net surrogate model effectively captures the intricate representations of aerosol and the relationship with atmospheric state variables, achieving an R-squared score over $$95.7\%$$ for all the lognormal aerosol modes in the extrapolated regime. The validated model provides feature importance of input variables and their impact on the predictive capacity of the surrogate model in relation to the E3SM. The computational cost of online inference time deployed on CPUs and GPUs with lower precisions highlights ADON’s efficiency and potential in robust predictive modeling for large-scale Earth system computations.

Bai, Zhe↗
DOE OSTI2020

Influence of sun zenith angle on canopy clumping and the resulting impacts on photosynthesis

Addressing the impact of vegetation architecture on shortwave radiation transfer in land surface models is important for accurate weather forecasting, carbon budget estimates, and climate predictions. This paper investigates to what extent it is possible to retrieve structural parameters of two different parameterization schemes from direct transmittance derived from digital hemispherical photography and 3D radiative transfer modeling for two study sites with different vegetation canopy architectures. Neglecting the representation of 3D canopy structure in radiative transfer schemes leads to significant errors in shortwave radiation partitioning (up to 3.5 times more direct transmittance in the 3D model). Structural parameters, referred to as whole canopy ‘clumping indices’, were obtained in order to evaluate the impact of angular variation in clumping on shortwave radiation transfer. Impacts on photosynthesis were evaluated at site level with the UKESM land surface model, JULES. A comparison between flux tower derived and modeled photosynthesis indicates that considering zenith angular variations of structural parameters in the radiative transfer scheme of the UKESM land surface model significantly improves photosynthesis prediction in light limited ecosystems (from RMSE = 2.91 μmol CO 2 .m -2 .s -1 to RMSE = 1.51 μmol CO 2 .m -2 .s -1 , 48% smaller), typically with enhanced photosynthesis from bottom layers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hybridizing Machine Learning and Physically-based Earth System Models to Improve Prediction of Multivariate Extreme Events (AI Exploration of Wildland Fire Prediction)

Focal Areas: This project responds to two focal areas identified in the DOE Call for AI4ESP White Papers: 1) Predictive modeling through the use of artificial intelligence (AI) techniques, and 2) insights gleaned from complex data using explainable AI and big data analytics. Science Challenge: Large wildland fires (hereafter wildfires) appearing as high-impact compound climate extreme events are closely related to hydroclimate and water cycle extremes that modulate surface fuel supply and combustibility. These compound events have multivariate climatic features (e.g., temperature, precipitation, relative humidity, wind, lightning) and societal drivers (e.g., forest management, land use change, human caused ignitions). Meanwhile, they induce strong feedbacks to the coupled atmosphere, biosphere, and hydrosphere by perturbing regional and global radiation budget as well as ecological, biogeochemical, and water cycles across multiple spatiotemporal scales. The nonlinear interactions between these natural and anthropogenic components of the Earth system are too complex to be completely and adequately represented in today’s Earth system models (ESMs). The inherent stochastic nature of fire activity at all scales further increases the difficulty of its prediction using ESMs that are usually developed from deterministic equations and parameterizations. Besides, concurrence of long-term (decadal to interdecadal) global climate change and fire regime shifts overlapping with short-term (intraseasonal to interannual) variations of regional fire weather and burning activity confound predictability of these compound extreme events. We propose to address the above scientific challenges by using machine learning (ML)-based data-driven modeling techniques to integrate observations and physically-based ESMs’ simulations in a computationally efficient hybrid prediction system. This prediction system is supposed to characterize the wildfire’s sensitivity to climate and exogenous drivers at high resolution (~ 0.25°) on subseasonal to seasonal (S2S) timescales providing improved predictability and explainability. We will use the system to help identify: (1) What are the computational elements of a hybrid system needed to predict compound climate extreme events such as global wildfires? (2) What are the key drivers (either natural or anthropogenic) that modulate short-term variations of multivariate fire weather and burning activity over different regions? How can one take advantage of those driver-response relationships to improve the predictability of large wildfires on S2S time scales? (3) What are the underlying physical mechanisms and sources of improved predictability? Which ML techniques are optimal in revealing and adapting these mechanisms?

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Polar primary aerosols across the ocean-sea ice-snow-atmosphere interface: From sources to impacts

Primary aerosols play a critical role in polar climate systems, influencing cloud formation, precipitation, radiative balance, and surface energy budgets. This paper provides a comprehensive synthesis of primary aerosol sources, transformation and removal processes, and broader atmospheric impacts in polar regions, emphasizing their links to ocean and sea ice biogeochemistry. These aerosols (including sea salt, primary organic aerosol, and primary biological aerosol particles) originate from marine and cryospheric environments and are emitted through physical processes, such as wave breaking, bubble bursting, and blowing snow. Emission sources include seawater, sea ice, snow, and freshwater from river discharge and glacial runoff. Once airborne, these particles can serve as a chemical reservoir, influencing atmospheric composition and reactivity, and as seeds for cloud droplet and ice crystal formation, influencing cloud microphysics and polar climate. Despite their importance, many of the processes governing primary aerosol emissions and transformations remain poorly constrained. The most pressing knowledge gaps pertain to emission processes, limited spatiotemporal observational coverage, instrumentation constraints, parameterization development, and the integration of interdisciplinary expertise. To improve our understanding of primary aerosol drivers and their response to climate, future research efforts should prioritize strategically coordinated and cross-disciplinary process studies, advancements in measurement technologies and coverage, and close collaboration between modelers and observational scientists to inform and refine model parameterizations. As polar regions continue to undergo profound changes marked by increased precipitation, reduced sea and land ice, freshening oceans, and shifting ecosystem dynamics, characterizing present-day primary aerosol populations is vital. Improved understanding will be essential for anticipating future changes in aerosol-radiation and aerosol-cloud interactions and their implications for polar and global climate systems.

Aerosol-cloud↗
DOE OSTI2023

Quantifying the Impact of Vertical Resolution on the Representation of Marine Boundary Layer Physics for Global-Scale Models

Abstract While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement. Significance Statement In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Cosmological constraints on dark matter interactions with ordinary matter

Dark matter interactions with electrons or protons during the early Universe leave imprints on the cosmic microwave background and the matter power spectrum, and can be probed through cosmological and astrophysical observations. These interactions lead to momentum and heat exchange between the ordinary and dark matter components, which in turn results in a transfer of pressure from the ordinary to the dark matter. Here, we explore these interactions using a diverse suite of data: cosmic microwave background anisotropies, baryon acoustic oscillations, the Lyman-α forest, and the abundance of Milky-Way subhalos. We derive constraints using model-independent parameterizations of the dark matter–electron and dark matter–proton interaction cross sections and map these constraints onto concrete dark matter models. Our constraints are complementary to other probes of dark matter interactions with ordinary matter, such as direct detection, big bang nucleosynthesis, various astrophysical systems, and accelerator-based experiments. They exclude sufficiently large cross sections for a large range of dark matter masses, which cannot be accessed by direct-detection experiments due to the overburden from the Earth’s atmosphere or crust.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Soil respiration spatial and temporal variability in China between 1961 and 2014

Soil respiration (Rs) plays an important role in terrestrial-atmospheric carbon exchange but remains one of the least studied components of the carbon cycle. How environmental changes influence Rs, and in turn, how Rs influences terrestrial carbon storage in China is unclear. Here, we estimated spatial patterns and temporal trends in Rs from 1961 to 2014 to determine the influence of the recent warming hiatus on the Rs temporal variability in China. We evaluated the relationship between Rs and a set of environmental factors and developed biome-specific Rs models that were parameterized based on 2,862 Rs records and 18 continuous variables (e.g., air temperature, precipitation and leaf area index). Based on the biome-specific Rs models and environmental information, we estimated the annual Rs as well as the change trend of Rs in China during the study period. The mean annual Rs from 1961 to 2014 was 3.50 (confidence interval +/- 0.87) Pg center dot C center dot year(-1). Rs in China increased at a rate of 0.008 Pg center dot C center dot year(-2) between 1982 and 1998 (warming period), but was static between 1999 and 2014 (warming hiatus). This reduction in Rs increase rate may have an important role in enhancing terrestrial carbon sequestration. These findings confirm that temporal heterogeneity in Rs could be influenced by regional environmental changes, and further help in improving the understanding of soil carbon dynamics in China.

Jian, Jinshi↗
DOE OSTI2023

Description and Evaluation of a New Deep Convective Cloud Model Considering In-cloud Inhomogeneity

Convections still need to be parameterized in general circulation models (GCMs) in the coming decades. Performances of GCMs are significantly influenced by the convection schemes used. In contrast to most conventional cloud models that ignore in-cloud inhomogeneities, a new convective cloud model explicitly considering in-cloud inhomogeneities is developed. The new model adopts a single bulk plume approach, but divides the plume into a series of interacting sub-plumes in order to mimic the transition from the plume core to its edges. We implemented the new cloud model in NCAR Community Atmosphere Model version 5.0 (CAM5) and evaluated its performance. Single-column tests show a stronger mass flux profile and low-level detrainment with the new model. In global simulations, the long-lasting wet bias in the tropical free troposphere of CAM5 is alleviated. Spatial pattern and intensity distribution of precipitation is improved, but the global hydrological cycle is slightly over-enhanced. The diurnal cycle of tropical precipitation is also better captured though the precipitation peak still occurs earlier than the observations and the amplitude of the cycle tends to weaken. MJO simulation is slightly improved with the new scheme. In addition, the new cloud model generated more tropical cyclones, in better agreement with observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Coupled Ocean Wave-Atmosphere Models for Offshore Wind Energy Applications

Coupling ocean wave models to mesoscale atmospheric models is necessary to represent the effect of waves on wind turbine hub-height winds. In this report we provide a review of the most widely used ocean waves models and the phased-averaged spectral wave modeling paradigm that they are based on. Methodologies used to couple these wave models to mesoscale atmospheric models are described along with details of existing coupled modeling systems. We summarize impacts on offshore wind resource assessments that have been investigated with such coupled modeling systems to-date. Specifically, coupling in the North and Baltic Seas was shown to have a small negative effect on the offshore wind energy resource at weak-to-moderate wind speeds, with little impact at higher wind speeds. Finally, limitations of these existing coupled modeling systems and impact assessments are discussed, including examples of and the potential for improved parameterizations of the air-sea fluxes, and the use of fine-scaled simulations using phased-resolved wave models.

17 WIND ENERGY↗
DOE OSTI2025

Improving the QBO Forcing by Resolved Waves With Vertical Grid Refinement in E3SMv2

Abstract The quasi‐biennial oscillation (QBO) is the dominate mode of variability in the tropical stratosphere and plays an important role in stratospheric dynamics and chemistry. The QBO is notably deficient in many climate models, including the Energy Exascale Earth System Model (E3SM) developed by the US Department of Energy. In this work, we refine the lower stratospheric vertical grid spacing from roughly 1 km to 500 m to facilitate more realistic equatorial wave activity in the lower stratosphere in E3SM version 2. The refinement results in a simulated QBO with a reasonable amplitude and easterly‐westerly transition in both directions, but still has a longer period than observed, slower easterly downward propagation speed, and shallower vertical depth. Similar refinement in the multi‐scale modeling framework configuration of E3SM yields similar improvements. By analyzing the forcing contributions from different wave types, we find that most of the QBO forcing still comes from parameterized gravity wave drag from convection. The improved QBO forcing contributions from resolved waves, especially equatorial Kelvin waves and resolved small scale waves, can be attributed to the grid refinement.

E3SM↗
DOE OSTI2023

Long-term variability in immersion-mode marine ice-nucleating particles from climate model simulations and observations

Abstract. Ice-nucleating particles (INPs) in the Southern Ocean (SO) atmosphere have significant impacts on cloud radiative and microphysical properties. Yet, INP prediction skill in climate models remains poorly understood, in part because of the lack of long-term measurements. Here we show, for the first time, how model-simulated INP concentrations compare with year-round INP measurements during the Macquarie Island Cloud Radiation Experiment (MICRE) campaign from 2017–2018. We simulate immersion-mode INP concentrations using the Energy Exascale Earth System Model version 1 (E3SMv1) by combining simulated aerosols with recently developed deterministic INP parameterizations and the native classical nucleation theory (CNT) for mineral dust in E3SMv1. Because MICRE did not collect aerosol measurements of super-micron particles, which are more effective ice nucleators, we evaluate the model's aerosol fields at other high-latitude sites using long-term in situ observations of dust and sea spray aerosol. We find that the model underestimates dust and overestimates sea spray aerosol concentrations by 1 to 2 orders of magnitude for most of the high-latitude sites in the Southern Hemisphere. We next compare predicted INP concentrations with concentrations of INPs collected on filter samples (typically for 2 or 3 d) and processed offline using the Colorado State University ice spectrometer (IS) in immersion freezing mode. We find that when deterministic parameterizations for both dust and sea spray INPs are used, simulated INPs are within a factor of 10 of observed INPs more than 60 % of the time during summer. Our results also indicate that the E3SM's current treatment of mineral dust immersion freezing in the SO is impacted by compensating biases – an underprediction of dust amount was compensated by an overprediction of its effectiveness as INPs. We also perform idealized droplet freezing experiments to quantify the implications of the time-dependent behavior assumed by the E3SM's CNT-parameterization and compare with the ice spectrometer observations. We find that the E3SM CNT 10 s diagnostic used in this study is a reasonable approximation of the exact formulation of CNT, when applied to ice spectrometer measurements in low-INP conditions similar to Macquarie Island. However, the linearized 10 s diagnostic underestimates the exact formula by an order of magnitude or more in places with high-INP conditions like the Sahara. Overall, our findings suggest that it is important to correct the biases in E3SM's simulated dust life cycle and update E3SM's INP parameterizations. INP prediction errors of 2 to 3 orders of magnitude can have considerable impacts on the simulated cloud and radiative properties in global climate models. On comparing INP concentrations during MICRE against ship-based campaigns, Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) and Antarctic Circumnavigation Expedition (ACE), we find that INPs from the latter are significantly higher only in regions closer to Macquarie Island. This alludes to the fact that physical, chemical and biological processes affecting INP concentrations as stimulated by the island could be partly responsible for the high INP concentrations observed at Macquarie Island during the MICRE campaign. Therefore, improvements to both aerosol simulation and INP parameterizations are required to adequately simulate INPs and their cloud impacts in E3SM. It will be helpful to include a parallel measurement of the size-resolved aerosol composition and explore opportunities for long-term measurement platforms in future field campaigns studying INP sources in remote marine regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗
DOE OSTI2021

A Bimodal Diagnostic Cloud Fraction Parameterization. Part II: Evaluation and Resolution Sensitivity

A wide range of approaches exists to account for subgrid cloud variability in regional simulations of the atmosphere. This paper addresses the following questions: 1) Is there still benefit in representing subgrid variability of cloud in convection-permitting simulations? 2) What is the sensitivity to the cloud fraction parameterization complexity? 3) Are current cloud fraction parameterizations scale-aware across convection-permitting resolutions? These questions are addressed for regional simulations of a 6-week observation campaign in the U.S. southern Great Plains. Particular attention is given to a new diagnostic cloud fraction scheme with a bimodal subgrid saturation-departure PDF, described in Part I. The model evaluation is performed using ground-based remote sensing synergies, satellite-based retrievals, and surface observations. It is shown that not using a cloud fraction parameterization results in underestimated cloud frequency and water content, even for stratocumulus. The use of a cloud fraction parameterization does not guarantee improved cloud property simulations, however. Diagnostic and prognostic cloud schemes with a symmetric subgrid saturation-departure PDF underestimate cloud fraction and cloud optical thickness, and hence overestimate surface shortwave radiation. These schemes require empirical bias-correction techniques to improve the cloud cover. The new cloud fraction parameterization, introduced in Part I, improves cloud cover, liquid water content, cloud-base height, optical thickness, and surface radiation compared to schemes reliant on a symmetric PDF. Furthermore, cloud parameterizations using turbulence-based, rather than prescribed constant subgrid variances, are shown to be more scale-aware across convection-permitting resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Tropical dry forest response to nutrient fertilization: a model validation and sensitivity analysis

Abstract. Soil nutrients, especially nitrogen (N) and phosphorus (P), regulate plant growth and hence influence carbon fluxes between the land surface and atmosphere. However, how forests adjust biomass partitioning to leaves, wood, and fine roots in response to N and/or P fertilization remains puzzling. Recent work in tropical forests suggests that trees increase fine root production under P fertilization, but it is unclear whether mechanistic models can reproduce this dynamic. In order to better understand mechanisms governing nutrient effects on plant allocation and improve models, we used the nutrient-enabled ED2 model to simulate a fertilization experiment being conducted in a secondary tropical dry forest in Costa Rica. We evaluated how different allocation parameterizations affected model performance. These parameterizations prescribed a linear relationship between relative allocation to fine roots and soil P concentrations. The slope of the linear relationship was allowed to be positive, negative, or zero. Some parameterizations realistically simulated leaf, wood, and fine root production, and these parameterizations all assumed a positive relationship between relative allocation to fine roots and soil P concentration. Model simulations of a 30-year timeframe indicated strong sensitivity to parameterization and fertilization treatment. Without P fertilization, the simulated aboveground biomass (AGB) accumulation was insensitive to the parameterization. With P fertilization, the model was highly sensitive to the parameterization and the greatest AGB accumulation occurred when relative allocation to fine roots was independent of soil P. Our study demonstrates the need for simultaneous measurements of leaf, wood, and fine root production in nutrient fertilization experiments and for longer-term experiments. Models that do not accurately represent allocation to fine roots may be highly biased in their simulations of AGB, especially on multi-decadal timescales.

Environmental Sciences & Ecology↗
DOE OSTI2023

Tracking Aerosol Convection Interactions Experiment (TRACER) Field Campaign Report

Convective clouds serve a critical role in the Earth’s energy and water cycles through their transport of heat, moisture, momentum, and chemical species through the troposphere driving the global circulation (e.g., Hartmann et al. 1984, Del Genio et al. 2012, Su et al. 2014). On more local scales, convective clouds impact the atmospheric heating profile through diabatic heating effects, removal of water from the atmospheric column through precipitation, and conditioning of the local environment impacting further development of clouds (e.g., Sullivan and Voigt 2021). These critical roles underscore the importance of realistic representation of convective processes across scales of models from large-eddy simulation (LES), to convection-permitting models (CPM; e.g., Kendon et al. 2020, Marinescu et al. 2021), to numerical weather prediction (NWP) models used for operational weather forecasting, to Earth system models used to predict climate sensitivity (Sanderson et al. 2011, Sherwood et al. 2014, Tomassini et al. 2014, Zhao et al. 2016, Cronin et al. 2017). A key component of improving model representation of convective clouds is better quantification and parameterization of updraft microphysics and dynamics, including their interactions with the surrounding environment and storm organization (Bony et al. 2015, Hagos and Houze 2016, Donner et al. 2016, Morrison et al. 2020). Aerosol is an important environmental factor that could affect convective clouds and precipitation since cloud droplet and ice formation processes are initiated by it. Andrae et al. (2004) hypothesized that aerosols associated with increased biomass burning particles acting as cloud condensation nuclei (CCN) result in smaller and more monodisperse cloud droplets leading to suppression of warm rain formation, ultimately leading to more cloud water being lofted above the freezing level based on observations in the Amazon region. The subsequent increase in latent heat release increases the buoyancy of rising convective parcels invigorating the deep convection. This work was followed by a description of the theoretical basis for this “cold-phase invigoration” by Rosenfeld et al. (2008), who argued that it could have a significant effect for deep convective clouds with warm cloud-bases. Several modeling studies (e.g., Khain et al. 2005, 2009, van den Heever et al. 2006, Fan et al. 2007, 2009, 2012, Lee et al. 2008, Storer et al. 2010, Lebo et al. 2012, Storer and van den Heever 2013, Chen et al. 2020, Dagan et al. 2022) have investigated these aerosol-convection interactions and the environmental factors that influence their relative importance and magnitude. More recently, several studies have indicated that “warm-phase invigoration”, the enhancement of convection through condensational heating, also appears to play a role in enhancing both shallow cumuli (Seiki and Nakajima 2014, Saleeby et al 2015) and deeper tropical convection (Lebo and Seinfeld 2011, Khain et al. 2012, Sheffield et al 2015, Fan et al. 2018, Igel and van den Heever 2021), as well as Houston thunderstorms (Fan et al. 2007, 2020). However, still other studies have provided additional evidence of systematic biases in simulated convective outflow ice size distribution properties, which are consistent with a lack of poorly understood secondary ice production within convective updrafts (e.g., Fridlind et al. 2017). To help address these critical gaps in our understanding of cloud processes, aerosol processes and aerosol-cloud interactions, the Tracking Aerosol Convection Interactions Experiment was designed building upon efforts by the Aerosol, Cloud, Precipitation and Climate (ACPC) Initiative (http://acpcintiative.org/), a joint effort of the International Geosphere-Biosphere Programme (IGBP) and the World Climate Research Program (WCRP) that focused on resolving uncertainties in the interactions between aerosol and clouds towards better understanding the role that these interactions play in the climate system. The TRACER campaign was motivated by recommendations from a number of pilot studies undertaken by ACPC (van den Heever et al. 2017, Fridlind et al. 2019, Hu et al. 2019, Fan et al. 2020, Marinescu et al. 2021, Hernandez-Deckers et al. 2022) that pointed towards the southeastern Texas region as a locale where aerosol-convection interactions could be studied owing to the copious occurrence of isolated convection during the summer months accompanied by diverse and significant sources of aerosols from both anthropogenic and natural sources. The TRACER campaign began on 01 October 2021 and extended through 30 September 2022 with an intensive operational period (IOP) during June-September 2022. Three main sites (Table 1) were managed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

QUIC-URB and QUIC-fire extension to complex terrain: Development of a terrain-following coordinate system

Ensemble-based approaches to prescribed fire planning cannot be supported by CFD-based models like FIRETEC and WFDS because they are too computationally expensive and cannot leverage LES approaches like CAWFE and WRF-SFIRE because too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, typical for instance of the Western United States. In this paper, we describe the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. In particular, the paper presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captured the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. We report the model agreed well with measurements on the upwind side of the peak, but overestimated speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field was obtained in 6 s, making the solver appealing for coupled fire–atmosphere simulations. The Socorro Mountain was highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect was addressed in this work). Results highlight that future work should focus on the implementation of parameterizations of wake-eddies, similar to QUIC-URB’s building parameterizations, and on thermodynamic-driven flow.

54 ENVIRONMENTAL SCIENCES↗