Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subgrid-scale processes”

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.

35 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Sensitivities of subgrid-scale physics schemes, meteorological forcing, and topographic radiation in atmosphere-through-bedrock integrated process models: a case study in the Upper Colorado River basin

Abstract. Mountain hydrology is controlled by interacting processes extending from the atmosphere through the bedrock. Integrated process models (IPMs), one of the main tools needed to interpret observations and refine conceptual models of the mountainous water cycle, require meteorological forcing that simulates the atmospheric process to predict hydroclimate then subsequently impacts surface–subsurface hydrology. Complex terrain and extreme spatial heterogeneity in mountainous environments drive uncertainty in several key considerations in IPM configurations and require further quantification and sensitivity analyses. Here, we present an IPM using the Weather Research and Forecasting (WRF) model which forces an integrated hydrologic model, ParFlow-CLM, implemented over a domain centered over the East River watershed (ERW), located in the Upper Colorado River basin (UCRB). The ERW is a heavily instrumented 300 km2 region in the headwaters of the UCRB near Crested Butte, CO, with a growing atmosphere-through-bedrock observation network. Through a series of experiments in the water year 2019 (WY19), we use four meteorological forcings derived from commonly used reanalysis datasets, three subgrid-scale physics scheme configurations in WRF, and two terrain shading options within WRF to test the relative importance of these experimental design choices for key hydrometeorological metrics including precipitation and snowpack, as well as evapotranspiration, groundwater storage, and discharge simulated by the ParFlow-CLM. Our hypothesis is that uncertainty from synoptic-scale forcings produces a much larger spread in surface–subsurface hydrologic fields than subgrid-scale physics scheme choice. Results reveal that the WRF subgrid-scale physics configuration leads to larger spatiotemporal variance in simulated hydrometeorological conditions, whereas variance across meteorological forcing with common subgrid-scale physics configurations is more spatiotemporally constrained. Despite reasonably simulating precipitation, a delay in simulated discharge peak is due to a systematic cold bias across WRF simulations, suggesting the need for bias correction. Discharge shows greater variance in response to the WRF simulations across subgrid-scale physics schemes (26 %) rather than meteorological forcing (6 %). The topographic radiation option has minor effects on the watershed-average hydrometeorological processes but adds profound spatial heterogeneity to local energy budgets (±30 W m−2 in shortwave radiation and 1 K air temperature differences in late summer). This is the first presentation of sensitivity analyses that provide support to help guide the scientific community to develop observational constraints on atmosphere-through-bedrock processes and their interactions.

54 ENVIRONMENTAL SCIENCES↗
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↗
DOE OSTI2022

Limitations of Bin and Bulk Microphysics in Reproducing the Observed Spatial Structure of Light Precipitation

Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios ( q c and q r , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulated q c closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulated q r properties, though, with bulk simulated q r consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produce q r that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Implicit learning of convective organization explains precipitation stochasticity

Accurate prediction of precipitation intensity is crucial for both human and natural systems, especially in a warming climate more prone to extreme precipitation. Yet, climate models fail to accurately predict precipitation intensity, particularly extremes. One missing piece of information in traditional climate model parameterizations is subgrid-scale cloud structure and organization, which affects precipitation intensity and stochasticity at coarse resolution. Here, using global storm-resolving simulations and machine learning, we show that, by implicitly learning subgrid organization, we can accurately predict precipitation variability and stochasticity with a low-dimensional set of latent variables. Using a neural network to parameterize coarse-grained precipitation, we find that the overall behavior of precipitation is reasonably predictable using large-scale quantities only; however, the neural network cannot predict the variability of precipitation (R 2 ~ 0.45) and underestimates precipitation extremes. The performance is significantly improved when the network is informed by our organization metric, correctly predicting precipitation extremes and spatial variability (R 2 ~ 0.9). The organization metric is implicitly learned by training the algorithm on a high-resolution precipitable water field, encoding the degree of subgrid organization. The organization metric shows large hysteresis, emphasizing the role of memory created by subgrid-scale structures. We demonstrate that this organization metric can be predicted as a simple memory process from information available at the previous time steps. These findings stress the role of organization and memory in accurate prediction of precipitation intensity and extremes and the necessity of parameterizing subgrid-scale convective organization in climate models to better project future changes of water cycle and extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Using aircraft measurements to characterize subgrid-scale variability of aerosol properties near the Atmospheric Radiation Measurement Southern Great Plains site

Abstract. Complex distributions of aerosol properties evolve in space and time as a function of emissions, new particle formation, coagulation, condensational growth, chemical transformation, phase changes, turbulent mixing and transport, removal processes, and ambient meteorological conditions. The ability of chemical transport models to represent the multi-scale processes affecting the life cycle of aerosols depends on their spatial resolution since aerosol properties are assumed to be constant within a grid cell. Subgrid-scale-dependent processes that affect aerosol populations could have a significant impact on the formation of particles, their growth to cloud condensation nuclei (CCN) sizes, aerosol–cloud interactions, dry deposition and rainout and hence their burdens, lifetimes, and radiative forcing. To address this issue, we characterize subgrid-scale variability in terms of measured aerosol number, size, composition, hygroscopicity, and CCN concentrations made by repeated aircraft flight paths over the Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystem (HI-SCALE) campaign. Subgrid variability is quantified in terms of both normalized frequency distributions and percentage difference percentiles using grid spacings of 3, 9, 27, and 81 km that represent those typically used by cloud-system-resolving models as well as the current and next-generation climate models. Even though the SGP site is a rural location, surprisingly large horizontal gradients in aerosol properties were frequently observed. For example, 90 % of the 3, 9, and 27 km cell mean organic matter concentrations differed from the 81 km cell around the SGP site by as much as ∼ 46 %, large spatial variability in aerosol number concentrations and size distributions were found during new particle formation events, and consequently 90 % of the 3, 9, and 27 km cell mean CCN number concentrations differed from the 81 km cell mean by as much as ∼ 38 %. The spatial variability varied seasonally for some aerosol properties, with some having larger spatial variability during the spring and others having larger variability during the late summer. While measurements at a single surface site cannot reflect the surrounding variability of aerosol properties at a given time, aircraft measurements that are averaged within an 81 km cell were found to be similar to many, but not all, aerosol properties measured at the ground SGP site. This analysis suggests that it is reasonable to directly compare most ground SGP site aerosol measurements with coarse global climate model predictions. In addition, the variability quantified by the aircraft can be used as an uncertainty range when comparing the surface point measurements with model predictions that use coarse grid spacings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Vertical dependence of horizontal variation of cloud microphysics: observations from the ACE-ENA field campaign and implications for warm-rain simulation in climate models

Abstract. In the current global climate models (GCMs), the nonlinearity effect of subgrid cloud variations on the parameterization of warm-rain process, e.g., the autoconversion rate, is often treated by multiplying the resolved-scale warm-rain process rates by a so-called enhancement factor (EF). In this study, we investigate the subgrid-scale horizontal variations and covariation of cloud water content (qc) and cloud droplet number concentration (Nc) in marine boundary layer (MBL) clouds based on the in situ measurements from a recent field campaign and study the implications for the autoconversion rate EF in GCMs. Based on a few carefully selected cases from the field campaign, we found that in contrast to the enhancing effect of qc and Nc variations that tends to make EF > 1, the strong positive correlation between qc and Nc results in a suppressing effect that tends to make EF < 1. This effect is especially strong at cloud top, where the qc and Nc correlation can be as high as 0.95. We also found that the physically complete EF that accounts for the covariation of qc and Nc is significantly smaller than its counterpart that accounts only for the subgrid variation of qc, especially at cloud top. Although this study is based on limited cases, it suggests that the subgrid variations of Nc and its correlation with qc both need to be considered for an accurate simulation of the autoconversion process in GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Updraft dynamics and microphysics: on the added value of the cumulus thermal reference frame in simulations of aerosol–deep convection interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Understanding drivers and uncertainty in projected African precipitation

We investigate the drivers of projected summer precipitation changes and their uncertainties across Africa in the second half of the 21st century under the SSP2-4.5 and SSP5-8.5 scenarios using CMIP6 models. Our results reveal distinct regional precipitation changes, particularly under SSP5-8.5, with robust increases of 75%, 24%, and 17% over the Sahara, South Eastern Africa, and South Central Africa, respectively, and a decline of up to 5% over West Southern Africa (WSAF). In most regions, precipitation increases are driven by enhanced vertical thermodynamic processes associated with temperature-induced moisture increases and enhanced moisture convergence. In contrast, the WSAF decrease is associated with vertical dynamic processes driven by a weakening of the Hadley circulation’s ascending branch. Model uncertainty accounts for over 85% of total projection uncertainty across all regions and is largely due to subgrid-scale parameterizations. Overall, this study enhances our understanding of climate change impacts on African precipitation.

Taguela, Thierry N↗
DOE OSTI2023

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Unified Entrainment and Detrainment Closures for Extended Eddy-Diffusivity Mass-Flux Schemes

The extended Eddy Diffusivity Mass Flux (EDMF) scheme is used as a unied parameterization that simulates subgrid-scale dynamics across a range of regimes.This work presents unied entrainment and detrainment closures that allow this scheme to simulate the dry convective boundary layer, shallow convection and deep convection , while adjacent works present other necessary closures for this unied scheme and test its performance for stable, stratocumulus topped and convective boundary layers. Entrainment and detrainment are modeled as a combination of turbulent and dynamical processes. The turbulent component is taken as downgradient diffusion between subdomains. The dynamical component is based on a ratio of buoyancy and vertical velocity that ensures consistent treatment of the ratio of mass-flux at vanishing area fraction. Dynamical detrainment is controlled by the difference in relative humidity between subdomains when liquid water is present. These closures are rst evaluated with entrainment and detrainment diagnosed from Large Eddy Simulations where decaying tracers are used to identify the different subdomains and mass exchanges between them. When implemented in a single column model (SCM), these closures are tested in simulations of several observationally motivated cases of dry convective boundary layer, shallow convection and deep convection. Results show good agreement between the SCM and Large Eddy Simulations in key properties for climate models and profiles of first, second and third moments. The development of diurnal cycles and the onset of and precipitation (with simplified microphysics) are well predicted by the SCM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗
DOE OSTI2020

Performance and Accuracy Implications of Parallel Split Physics-Dynamics Coupling in the Energy Exascale Earth System Atmosphere Model

Simultaneous calculation of atmospheric processes is faster than calculating processes one at a time. This type of parallelism is beneficial or perhaps even necessary to provide good performance on modern supercomputers, which achieve faster performance through increased processor count rather than improved clock speed. The scalability of the Energy Exascale Earth System Model (E3SM) Atmosphere Model (EAM) is limited by the fluid dynamics which scales up to the number of mesh cells in the global mesh. In contrast, the suite of physics parameterizations in EAM is scalable up to the total number of physics columns, which is an order of magnitude greater than the number of mesh cells. A proposed solution to unlocking the greater potential performance from the physics suite is to solve the physics and dynamics in parallel. This work represents a first attempt at parallel splitting of the grid-scale fluid dynamics model and the subgrid-scale physics parameterizations in a global atmosphere model. We will demonstrate that switching to parallel physics-dynamics coupling extends the scalability of the EAM to up to 3 times the previous peak scalability limit and is up to 20% faster than the sequentially split coupling at the highest core counts and the same time step. Decadal simulations of both coupling approaches show very little impact to the model climate. This improved performance does not come without drawbacks, however. Parallel splitting requires a shorter time step and other modifications which largely offset performance gains. A mass fixer is required for conservation. Techniques for mitigating these issues are also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Bimodality in simulated precipitation frequency distributions and its relationship with convective parameterizations

Abstract Bimodality in precipitation frequency distributions is often evident in atmospheric models, but rarely in observations. This study i) proposes a metric to objectively quantify the bimodality in precipitation distributions, ii) evaluates model simulations contributed to the Coupled Model Intercomparison Project (CMIP) phase 5 (CMIP5), phase 6 (CMIP6), and the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project by comparing them to satellite-based and reanalysis precipitation products, and iii) investigates possible origins of bimodal precipitation distributions. Our results reveal that about 83% (20 out of 24) of CMIP5 and 70% (21 out of 30) of CMIP6 models used in this study exhibit bimodal distributions. The few DYAMOND models that use a deep convective parameterization also show bimodal distributions, while most DYAMOND models do not. Predictably, the bimodality originates from the separation of precipitation process between resolved grid-scale and parameterized subgrid-scale. However, in a larger number of models bimodality arises from the parameterized subgrid-scale convective precipitation alone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Robustness of the Stochastic Parameterization of Subgrid-Scale Wind Variability in Sea Surface Fluxes

Abstract High-resolution numerical models have been used to develop statistical models of the enhancement of sea surface fluxes resulting from spatial variability of sea surface wind. In particular, studies have shown that flux enhancement is not a deterministic function of the resolved state. Previous studies focused on single geographical areas or used a single high-resolution numerical model. This study extends the development of such statistical models by considering six different high-resolution models, four different geographical regions, and three different 10-day periods, allowing for a systematic investigation of the robustness of both the deterministic and stochastic parts of the data-driven parameterization. Results indicate that the deterministic part, based on regressing the unresolved normalized flux onto resolved-scale normalized flux and precipitation, is broadly robust across different models, regions, and time periods. The statistical features of the stochastic part of the model (spatial and temporal autocorrelation and parameters of a Gaussian process fit to the regression residual) are also found to be robust and not strongly sensitive to the underlying model, modeled geographical region, or time period studied. Best-fit Gaussian process parameters display robust spatial heterogeneity across models, indicating potential for improvements to the statistical model. These results illustrate the potential for the development of a generic, explicitly stochastic parameterization of sea surface flux enhancements dependent on wind variability.

Endo, Kota↗
DOE OSTI2022

Evaluation of Model Summertime Boundary Layer Cloud Development over Complex Terrain in New York State

Weather forecasting over complex terrain with diverse land cover is challenging. Utilizing the high-resolution observations from New York State Mesonet (NYSM), we are able to evaluate the surface processes of the Weather Research Forecast (WRF) Model in a detailed, scale-dependent manner. Here, in the study, possible impacts of land–atmosphere interaction on surface meteorology and boundary layer cloud development are investigated with different model resolutions, land surface models (LSMs), and planetary boundary layer (PBL) physical parameterizations. The High-Resolution Rapid Refresh, version 3 (HRRR), forecasting model is used as a reference for the sensitivity evaluation. Results show that over complex terrain, the high-resolution simulations (1 km × 60 vertical levels) generally perform better compared to low-resolution (3 km × 50 levels) in both surface meteorology and cloud fields. LSMs play a more important role in surface meteorology compared to PBL schemes. The NoahMP land surface model exhibits daytime warmer and drier biases compared to the Rapid Update Cycle (RUC) due to better prediction of the Bowen ratio in RUC. The PBL schemes would affect the convective strength in the boundary layer. The Shin–Hong (SH) scale-aware scheme tends to produce the strongest convective strength in the PBL, while the ACM2 PBL scheme rarely resolved convection even at 1-km resolution. By considering the radiation effect of subgrid-scale (SGS) clouds, the Mellor–Yamada–Nakanishi–Niino eddy diffusivity mass flux (MYNN-EDMF) predicted the highest cloud coverage and lowest surface solar radiation bias. The configuration of SGS clouds in MYNN-EDMF would not only significantly reduce shortwave radiation bias, but also affect the convection behaviors through land surface–cloud–radiation interaction.

54 ENVIRONMENTAL SCIENCES↗