Search NASA⌕ Search

SEARCH · Search NASA

Results for “Precipitation parameterization”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A riming‐dependent parameterization of scattering by snowflakes using the self‐similar Rayleigh–Gans approximation

Abstract Riming is a key process of precipitation formation in ice‐containing clouds, but quantifying riming from observations is challenging, limiting our ability to evaluate the riming process in numerical weather models. One challenge for radar observations is that riming changes both the physical properties (mass, area cross‐section) and scattering properties of ice particles. These changes need to be implemented consistently as a function of riming in radar forward operators, which are required for retrievals and model evaluation in observation space. In this study, mass–size, cross‐section area–size, and backscattering cross‐section relations are developed as a function of the normalized rime mass for aggregates composed of various monomer types (columns, dendrites, needles, plates, and rosettes). The proposed framework allows us to simulate scattering properties of aggregated ice particles consistently as a function of riming in retrievals and radar forward operators. The parameterizations are developed from a large data set of simulated rimed aggregates of different sizes and monomer crystal types. The backscattering cross‐section parameterization (the “riming‐dependent parameterization”) is evaluated for radar frequencies of 35.6 and 94.0 GHz and is based on the Self‐Similar Rayleigh–Gans approximation (SSRGA), which is increasingly used to calculate microwave scattering of ice crystals and snowflakes. Compared with parameterizations from the literature that do not consider riming, the riming‐dependent parameterization leads to significantly smaller biases in terms of backscattering cross‐section. When using the particle masses and scattering properties of the individual particles simulated by the aggregation and riming model as a reference, the bias of our parameterization is below 1 dB when integrating over an exponential particle size distribution with sizes from 0.1–10 mm.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mesoscale Convective Systems in a Superparameterized E3SM Simulation at High Resolution

Accurately representing mesoscale convective systems (MCSs) is crucial to simulating the energy and water cycles in global climate models. Using a novel MCS identification and tracking algorithm applied to observations and model simulations, we evaluate how well the Energy Exascale Earth System Model simulates MCSs over the central US. Simulations performed by E3SM at 25 km grid spacing with (SP-E3SM) and without (E3SM) superparameterization are compared and evaluated against observations using multiple metrics highlighting important MCS characteristics. Compared to E3SM, SP-E3SM better simulates MCS number and MCS precipitation amount, diurnal cycle, propagation, and the probability distribution of precipitation rate in both spring and summer. The improvement from SP is partly contributed by improvement in simulating the large-scale environments, featuring enhanced atmospheric low-level moisture and larger moisture transport to the central US relative to E3SM. However, SP-E3SM still underestimates MCS precipitation amount, particularly in summer. This underestimation is closely related to the negative bias in MCS precipitation intensity, although the drier environments simulated during summer also contributes to underestimation of MCS frequency. Without SP, the larger bias in MCS precipitation amount is closely related to the negative bias in MCS frequency, while the substantial dry bias in the large-scale environments also contributes to underestimation of MCS intensity. Our results suggest that SP improves MCS simulation by improving modeling of the large-scale environments and convection initiation, which are both major limiting factors in E3SM even at 25 km grid spacing where deep convection is represented by a cumulus parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Tropical Cyclone Wind Shear-Relative Asymmetry in Reanalyses

Abstract While tropical cyclones (TCs) are axisymmetric vortices to the first order, they often exhibit noteworthy structural asymmetries. These often result from environmental vertical wind shear, which tilts the vortex and induces a wavenumber 1 pattern in the circulation and precipitation fields. Reanalyses and climate models have improved in representing the TC structure and climatology, but their relatively coarse resolution and dependence on parameterized physics cast doubt on their ability to capture the asymmetric TC structure. We perform the most comprehensive process-oriented assessment of TC asymmetry to date in reanalyses. Specifically, we analyze the composite shear-relative TC structure in ERA5 and Climate Forecast System Reanalysis (CFSR), which vary in their resolutions, physical parameterization suites, and data assimilation techniques. These structures are compared with aircraft reconnaissance radar observations. In agreement with the observations, the strongest tangential winds are usually found left-of-shear, while inner core rainfall, ascent, vortex tilt, and low-level inflow are favored directly downshear or in the downshear-left quadrant. Outer rainband convection generally peaks in the downshear-right quadrant. Thermodynamic asymmetries are also apparent, with anomalous low-level moisture right-of-shear, midlevel warmth in the upshear-right quadrant (uptilt), and cloud properties suggestive of a realistic precipitation life cycle from growth to fallout. We also decompose rainfall contributions from the convective parameterization and large-scale cloud schemes and highlight the roles of vorticity advection, buoyancy advection, and diabatic processes in driving asymmetric vertical motions in the inner core and outer rainband regions. Our results suggest that process-level studies of TC asymmetry and TC–wind shear interaction under future warming are viable using climate models. Significance Statement Asymmetries are common in tropical cyclones (TCs), influencing their intensity, track, and hazards. Vertical wind shear often plays a leading-order role in causing these asymmetries. It is uncertain how well asymmetric structures and processes are captured in reanalyses and global climate models (GCMs) with grid spacings of 0.25° and coarser. In this study, we first evaluate TC asymmetry in reanalyses, which have the benefit of being forced by observations. This helps to assess whether the resolutions associated with GCMs sufficiently capture asymmetric structures and processes and motivates upcoming work with free-running GCMs to study how TC asymmetry may change in a warming climate.

Carstens, Jacob D.↗
DOE OSTI2025

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Parameterizing Convective Organization Effects With a Moisture-PDF Approach in Climate Models: Concept and a Regional Case Simulation

We propose a parameterization scheme of convective organization effects based on a moisture-distribution approach, which can reflect aggregation of convective cells within a model grid as well as the interaction between convection and spatial heterogeneity in free-troposphere moisture. With this concept, convective cells in an aggregated state are surrounded by air that is moister than the grid-mean condition, which provides a “shielding” effect favorable for deeper convective updrafts. Such effects are represented quantitatively via utilizing a schematic diagram about the geometric interpretation of convective cluster, dry area, and their contact area, in which both the convective cluster size and dry area fraction increase during the development of convection. Our sensitivity analyses indicate that the new scheme performs well in capturing the mean precipitation features. Excluding the convective organization effects in the model leads to a considerable reduction in the simulated precipitation magnitude. Notably, there exists an inverse relationship between background mean moisture condition and sub-grid moisture variability for a given rain rate, suggesting that with a high spatial heterogeneity in free-troposphere moisture, convection can still develop and maintain its strength under relatively dry background conditions, consistent with the convection-permitting model simulation. Overall, our new parameterization of convective organization effects can successfully reproduce the relationship between precipitation and sub-grid moisture variability, which is a missing element in traditional convection parameterization schemes but important for the simulations of precipitation variability at various scales in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mechanisms in Regulating the Quasi‐Biennial Oscillation in Exascale Earth System Model Version 2

Abstract The update to the deep convection parameterization in the energy exascale Earth system model version 2 (E3SMv2) makes the simulated quasi‐biennial oscillation (QBO) have a shorter period than without the update when other tunable parameters the same since the update makes convection more intense but less frequent while leaving the time‐mean convective heating tendency almost unchanged. Amplitudes of momentum fluxes of parameterized gravity waves (GWs) are intensified since they are determined by the square of convective heating tendency. In addition, stronger planetary waves are simulated in E3SMv2 with the convection scheme update, partially due to the tropospheric precipitation change. Furthermore, there is evidence that planetary waves are intensified in the stratosphere due to the dissipation of enhanced parameterized GWs. These factors are found to be responsible for modulating the QBO simulated in E3SMv2.

Li, Yuanpu [Climate and Global Dynamics Division U↗
DOE OSTI2021

Transforming ESM Physical Parameterization Development Using Machine Learning Trained on Global Cloud-Resolving Models and Process Observations

ESMs robustly predict that 21st century greenhouse warming will slowly increase global mean precipitation, rapidly increase extreme precipitation, and increase subtropical drought. ESMs agree less about precipitation trends and extremes over particular land regions critical to human societies, e. g. in semi-arid regions such as California or the Sahel, or in wetter climates prone to monsoonal rainfall (e. g. southeast Asia) or to tropical cyclones and flooding from mesoscale convective systems (e. g. the southeastern U.S.) Deep convective parameterizations and poor representation of orography and complex vegetated land surfaces contribute to this inter-model spread; clouds, aerosols and sea-surface temperature biases are also key. Reducing regional precipitation projection uncertainty has enormous planning value for water supplies, land use, wildfire, hydropower, flood control, etc. IPCC-class ESMs are making painfully slow progress on this.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Remote Influence of Andean Convection on Amazonian Rainfall and Its Mechanisms

Models from Coupled Model Intercomparison Project Phase 6 produce too much precipitation over the Andes but too little over the Amazon or the Wet Andes-Dry Amazon (WADA) bias pattern. Unlike the conventional view that convection parameterization and land model deficiencies can contribute to Amazonian rainfall biases, we approach this long-standing biased model behavior through the lens of Andean convection. Using Community Earth System Model v1.1 and focusing on the wet season, our mechanism-denial experiments demonstrate that Andean convection notably reduces precipitation over the Amazon during austral summer. The Andean forced Amazonian response operates on weather timescale. Furthermore, the reduction of Amazonian rainfall is detectable within a few hours after initial Andean forcing. The precipitation response is primarily driven by variations in the moisture budget and is moderated by changes in convective available potential energy over the Amazon. Changes in the total advection of moisture over the Amazon are dominated by the vertical advection term and can be attributed to discrepancies in the dynamic omega field. In the experiments, the Andean east flank region is scrutinized where the vertical velocity and moisture fields play an intermediary role for the Andean driven WADA connection. The Andean forcing induces descending anomalies on the Andean east flank. The disturbances of wind and geopotential fields over the Andean east flank propagate eastward via Kelvin waves. Over the Amazon, descending anomalies and advective drying lead to reduction of mid-to-high level cloud, increase of shortwave cloud forcing and surface net radiation, and enhancement of themodynamic stability and rainfall reduction.

58 GEOSCIENCES↗
DOE OSTI2024

Implementation and Exploration of Parameterizations of Large-Scale Dynamics in NCAR's Single Column Atmosphere Model SCAM6

A single column model with parameterized large-scale (LS) dynamics is used to better understand the response of steady-state tropical precipitation to relative sea surface temperature under various representations of radiation, convection, and circulation. The large-scale dynamics are parametrized via the weak temperature gradient (WTG), damped gravity wave (DGW), and spectral weak temperature gradient (Spectral WTG) method in NCAR's Single Column Atmosphere Model (SCAM6). Radiative cooling is either specified or interactive, and the convective parameterization is run using two different values of a parameter that controls the degree of convective inhibition. Results are interpreted in the context of the Global Atmospheric System Studies -Weak Temperature Gradient (GASS-WTG) Intercomparison project. Using the same parameter settings and simulation configuration as in the GASS-WTG Intercomparison project, SCAM6 under the WTG and DGW methods produces erratic results, suggestive of numerical instability. However, when key parameters are changed to weaken the large-scale circulation's damping of tropospheric temperature variations, SCAM6 performs comparably to single column models in the GASS-WTG Intercomparison project. The Spectral WTG method is less sensitive to changes in convection and radiation than are the other two methods, performing qualitatively similarly across all configurations considered. Under all three methods, circulation strength, represented in 1D by grid-scale vertical velocity, is decreased when barriers to convection are reduced. This effect is most extreme under specified radiative cooling, and is shown to come from increased static stability in the column's reference radiative-convective equilibrium profile. This argument can be extended to interactive radiation cases as well, though perhaps less conclusively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Precipitation and nitrogen application stimulate soil nitrous oxide emission

Precipitation and nitrogen (N) fertilization are the two most important drivers for soil nitrous oxide (N 2 O) emission. However, the effects of changes in N fertilization and precipitation patterns (i.e., precipitation intensity and frequency) on N 2 O emissions in agricultural fields are still unclear. In this study, we simulated soil N 2 O emission under different precipitation patterns (6 precipitation intensities, and 12 precipitation frequencies by either merging or splitting precipitation events) and N fertilization rates (low, typical, and high N fertilization) in a cornfield using the DeNitrification-DeComposition model. The model was parameterized and validated using meteorological data and N experimental measurements in Nashville, Tennessee, USA. Results showed that soil water filled pore space (WFPS) and simulated soil N 2 O emission increased as precipitation intensity increased. Less frequent but high intensity precipitation treatments reduced the soil WFPS by 25.2% and stimulated soil N 2 O emission by 45.3%, while more frequent but low intensity precipitation treatments increased soil WFPS by 9.0% and reduced soil N 2 O emission by 23.9%. Compared to typical N fertilization, the sensitivity of soil N 2 O emission to precipitation was higher under high N than low N fertilization treatments, and the response ratios were 50.0% and 40.1%, respectively. There was significant interactive effect of precipitation intensity and N fertilization on soil N 2 O emission. Furthermore, these findings improved our understanding of precipitation and N impacts on soil N 2 O emissions and provided useful knowledge for irrigation and N fertilizer management in agriculture to mitigate greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Future changes in extreme precipitation over the San Francisco Bay Area: Dependence on atmospheric river and extratropical cyclone events

Extreme precipitation poses a major challenge for local governments, including the City and County of San Francisco, California, as flooding can damage and destroy infrastructure and property. As the climate continues to warm, reliable future precipitation projections are needed to provide the best possible information to decision makers. However, future changes in the magnitude of extreme precipitation are uncertain, as current state-of-the-art global climate models are typically run at relatively coarse horizontal resolutions that require the use of convective parameterization and have difficulty simulating observed extreme rainfall rates. Here, we performed ensembles of convection-permitting regional climate model simulations to investigate how five historically impactful extreme precipitation events over the San Francisco Bay Area could change if similar events occurred in future climates. We found that changes in storm-total precipitation depend strongly on storm type. Precipitation associated with an atmospheric river (AR) accompanied by an extratropical cyclone (ETC) is projected to increase at a rate exceeding (by up to 1.5 times) the theoretical Clausius Clapeyron scaling of 6–7% per °C warming. On the other hand, future precipitation changes are weak or negative for events characterized by an AR only, despite increases in precipitable water and integrated vapor transport that are similar to those of the co-occurring AR and ETC events. The differences in the sign of future precipitation change between AR-only events and co-occurring AR and ETC events is instead linked with changes in mid-tropospheric vertical velocity. Given that the majority of observed ARs are associated with an ETC, this research has important implications for future precipitation impacts over the Bay Area, as it indicates that storm-total precipitation associated with the most common type of storm event may increase by up to 26–37% in 2100 relative to historical.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Forecasting Remote Atmospheric Responses to Decadal Kuroshio Stability Transitions

Based on observational estimates and global ocean eddy-resolving coupled retrospective initialized predictions, we show that Kuroshio Extension variability affects rainfall variability along the west coast of North America. We show that the teleconnection between the current undulations and downstream rainfall can lead to improved subseasonal to seasonal predictions of precipitation over California, and we demonstrate that capturing these teleconnections requires coupled systems with sufficient ocean resolution (i.e., eddy-resolving), especially over time scales longer than one season. The improved forecast skill is diagnosed in terms of 35 years of retrospective initialized ensemble forecasts with an ocean eddy-resolving and an ocean eddy-parameterized coupled model. Not only does the ocean eddy-resolving model show sensitivity to Kuroshio Extension variability in terms of western North America precipitation, but the ocean eddy-resolving forecasts also show improved forecast skill compared to the ocean eddy-parameterized model. Additionally, the ocean eddy-parameterized coupled model shows no sensitivity to Kuroshio Extension variability. We also find near-decadal variability associated with a progression of a lower-tropospheric height dipole around the North Pacific and how these height anomalies lead to wind-driven Rossby waves that affect the eddy activity in the Kuroshio Extension with a time lag on the order of four years. This decadal-scale variability (~10 years) opens the possibility of multiyear predictability of western North American rainfall.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Assessing Raindrop Breakup Parameterizations Using Disdrometer Observations

In this work, an intercomparison of raindrop mean diameter frequency distribution (RDFD) is performed for numerical simulations of precipitating cloud systems using an array of models and microphysics schemes. This includes results from the Regional Atmospheric Modeling System (RAMS) double-moment microphysics, the Hebrew University Cloud Model bin microphysics (HUCM) interfaced to the RAMS parent model, and the Weather Research and Forecasting (WRF) Model with the Thompson, Morrison, double-moment 6-class (WDM6), and National Severe Storms Laboratory (NSSL) double-moment schemes. Simulations are examined with respect to the raindrop size distribution (DSD) volume-number mean diameter (D m ) and intercept parameter (N w ). When compared to a suite of disdrometer observations, the RDFD resulting from each microphysics scheme exhibits varying degrees of mean drop size constraints and peaks in the frequency distribution of D m . A more detailed investigation of the peaked RDFD from the RAMS simulations suggests that the parameterization of raindrop collisional breakup can impose strong limitations on the evolution of simulated drop growth. As such, a summary and comparison of the drop breakup parameterizations among the aforementioned microphysics schemes is presented. While some drop breakup parameterizations are adjusted toward the observations by modifying the threshold diameter for the onset of breakup, this study explores the use of a modified maximum breakup efficiency. This method permits the parameterization to retain its threshold breakup diameter, while limiting the strength of drop breakup and permitting a broader range of drop sizes. As a result, the simulated mean drop sizes are in better agreement with observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Convection–Kelvin Wave Coupling in a Global Convection-Permitting Model

A convectively coupled equatorial Kelvin wave (CCKW) was observed over the equatorial Indian Ocean in early November 2011 during the DYNAMO field campaign. This study examines the structure of the CCKW event using two simulations made using the MPAS model: one with 3-km grid spacing without convective parameterization and another with a 15-km grid and parameterized convection. Both simulations qualitatively capture the observed structure of the CCKW, including its vertical tilt and progression of cloud/precipitation structures. The two simulations, however, differ substantially in the amplitude of the CCKW-associated precipitation. While the 3-km run realistically captures the observed modulation of precipitation by the CCKW, the 15-km simulation severely underestimates its magnitude. To understand the difference between the two MPAS simulations regarding wave–convection coupling within the CCKW, the relationship of precipitation with convective inhibition, saturation fraction, and surface turbulent fluxes is investigated. Results show that the 15-km simulation underestimates the magnitude of the CCKW precipitation peak in association with its unrealistically linear relationship between moisture and precipitation. Precipitation, both in observations and the 3-km run, is predominantly controlled by saturation fraction and this relationship is exponential. In contrast, the parameterized convection in the 15-km run is overly sensitive to convective inhibition and not sensitive enough to environmental moisture. The implications of these results on CCKW theories are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Aerosol-Cloud Interactions From Aviation Soot Emissions

Current models estimate global aviation contributes approximately 5% to the total anthropogenic climate forcing, with aerosol-cloud interactions having the greatest effect. However, radiative forcing estimates from aviation aerosol-cloud interactions remain undetermined. There is an expected significant increase in aircraft emissions with aviation demand expected to rise by over 4% per year. Soot may play an important role in the ice nucleation of aircraft-induced cirrus formation due to a high emission rate, but the ice nucleating properties are poorly constrained. Understanding the microphysical processes leading to atmospheric ice crystal formation is crucial for the reliable parameterization of aerosol-cloud interactions in climate models due to their impact on precipitation and cloud radiative properties. Ice nucleation of aircraft-emitted soot is potentially affected by particle morphology with condensation of supercooled water occurring in pores followed by ice nucleation. However, soot has heterogeneous properties and undergoes atmospheric aging and oxidation that could change surface properties and contribute to complex ice nucleation processes. Further, this review synthesizes current knowledge of ice nucleation catalyzed by aviation in the cirrus regime and its effects on global radiative forcing. Further research is required to determine the ice nucleation and microphysical processes of cirrus cloud formation from aviation emissions in both controlled laboratory and field investigations to inform models for more accurate climate predictions and to provide efficient mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Decision support for United States—Canada energy integration is impaired by fragmentary environmental and electricity system modeling capacity

The renewable energy transition is leading to increased electricity trade between the United States and Canada, with Canadian hydropower providing firm lower-carbon power and buffering variability of wind and solar generation in the U.S. However, long-term power purchase agreements and transborder transmission projects are controversial, with two of four proposed transmission lines between Quebec, Canada and the northeast U.S. cancelled since 2018. Here, we argue that controversies are exacerbated by a lack of open-source data and tools to understand tradeoffs of new hydropower generation and transmission infrastructure in comparison to alternatives. This gap includes impacts that incremental transmission and generation projects have on the economics of the entire system, for example, how new transmission projects affect exports to existing markets or incentivize new generation. We identify priority areas for data synthesis and model development, such as integrating linked hydropower and hydrologic interactions in energy system models and openly releasing (by utilities) or back-calculating (by researchers) hydropower generation and operational parameters. Publicly available environmental (e.g. streamflow, precipitation) and techno-economic (e.g. costs, reservoir size,) data can be used to parameterize freely usable and extensible models. Existing models have been calibrated with operational data from Canadian utilities that are not publicly available, limiting the range of scientific and commercial questions these tools have been used to answer and the range of parties that have been involved. Studies conducted using highly resolved, national-scale public data exist in other countries, notably, the United States, and demonstrate how greater transparency and extensibility can drive industry action. Improved data availability in Canada could facilitate approaches that (1) increase participation in decarbonization planning by a broader range of actors; (2) allow independent characterizations of environmental, health, and economic outcomes of interest to the public; and (3) identify decarbonization pathways consistent with community values.

13 HYDRO ENERGY↗
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 OSTI2022

Subgrid-scale horizontal and vertical variation of cloud water in stratocumulus clouds: a case study based on LES and comparisons with in situ observations

Abstract. Stratocumulus clouds in the marine boundary layer cover a large fraction of ocean surface and play an important role in the radiative energy balance of the Earth system. Simulating these clouds in Earth system models (ESMs) has proven to be extremely challenging, in part because cloud microphysical processes such as the autoconversion of cloud water into precipitation occur at scales much smaller than typical ESM grid sizes. An accurate autoconversion parameterization needs to account for not only the local microphysical process (e.g., the dependence on cloud water content qc and cloud droplet number concentration Nc) but also the subgrid-scale variability of the cloud properties that determine the process rate. Accounting for subgrid-scale variability is often achieved by the introduction of a so-called enhancement factor E. Previous studies of E for autoconversion have focused more on its dependence on cloud regime and ESM grid size, but they have largely overlooked the vertical dependence of E within the cloud. In this study, we use a large-eddy simulation (LES) model, initialized and constrained with in situ and surface-based measurements from a recent airborne field campaign, to characterize the vertical dependence of the horizontal variation of qc in stratocumulus clouds and the implications for E. Similar to our recent observational study (Zhang et al., 2021), we found that the inverse relative variance of qc, an index of horizontal homogeneity, generally increases from cloud base upward through the lower two-thirds of the cloud and then decreases in the uppermost one-third of the cloud. As a result, E decreases from cloud base upward and then increases towards the cloud top. We apply a decomposition analysis to the LES cloud water field to understand the relative roles of the mean and variances of qc in determining the vertical dependence of E. Our analysis reveals that the vertical dependence of the horizontal qc variability and enhancement factor E is a combined result of condensational growth throughout the lower portion of the cloud and entrainment mixing at cloud top. The findings of this study indicate that a vertically dependent E should be used in ESM autoconversion parameterizations.

54 ENVIRONMENTAL SCIENCES↗