Because of their limited spatial resolution, numerical weather prediction and climate models have to rely on parameterizations to represent atmospheric turbulence and convection. Historically, largely independent approaches have been used to represent boundary layer turbulence and convection, neglecting important interactions at the subgrid scale. Here we build on an eddy‐diffusivity mass‐flux (EDMF) scheme that represents all subgrid‐scale mixing in a unified manner, partitioning subgrid‐scale fluctuations into contributions from local diffusive mixing and coherent advective structures and allowing them to interact within a single framework. The EDMF scheme requires closures for the interaction between the turbulent environment and the plumes and for local mixing. A second‐order equation for turbulence kinetic energy (TKE) provides one ingredient for the diffusive local mixing closure, leaving a mixing length to be parameterized. Here, we propose a new mixing length formulation, based on constraints derived from the TKE balance. It expresses local mixing in terms of the same physical processes in all regimes of boundary layer flow. The formulation is tested at a range of resolutions and across a wide range of boundary layer regimes, including a stably stratified boundary layer, a stratocumulus‐topped marine boundary layer, and dry convection. Comparison with large eddy simulations (LES) shows that the EDMF scheme with this diffusive mixing parameterization accurately captures the structure of the boundary layer and clouds in all cases considered.
This content will become publicly available on June 1, 2024
“Supermodeling” climate by allowing different models to assimilate data from one another in run time has been shown to give results superior to those of any one model and superior to any weighted average of model outputs. The only free parameters, connection strengths between corresponding variables in each pair of models, are determined using some form of machine learning. It is demonstrated that supermodeling succeeds because near critical states, interscale interactions are important but unresolved processes cannot be effectively represented diagnostically in any single parameterization scheme. In two examples, a pair of toy quasigeostrophic (QG) channel models of the midlatitudes and a pair of ECHAM5 models of the tropical Pacific atmosphere with a common ocean, supermodels dynamically combine parameterization schemes so as to capture criticality, associated critical structures, and the supporting scale interactions. The QG supermodeling scheme extends a previous configuration in which two such models synchronize with intermodel connections only between medium-scale components of the flow; here the connections are trained against a third “real” model. Intermittent blocking patterns characterize the critical behavior thus obtained, even where such patterns are missing in the constituent models. In the ECHAM-based climate supermodel, the corresponding critical structure is the single ITCZ pattern, a pattern that occurs in neither of the constituent models. For supermodels of both types, power spectra indicate enhanced interscale interactions in frequency or energy ranges of physical interest, in agreement with observed data, and supporting a generalized form of the self-organized criticality hypothesis.
In a “supermodel” of Earth’s climate, alternative models (climate simulations), which differ in the way they represent processes on the smallest scales, are trained to exchange information as they run, adjusting to one another much as weather prediction models adjust to new observations. They form a consensus, capturing atmospheric behaviors that have eluded all the separate models. We demonstrate that simplified supermodels succeed, where no single approach can, by correctly representing
- Award ID(s):
- 2015618
- NSF-PAR ID:
- 10497228
- Publisher / Repository:
- American Meteorological Society
- Date Published:
- Journal Name:
- Journal of the Atmospheric Sciences
- Volume:
- 80
- Issue:
- 6
- ISSN:
- 0022-4928
- Page Range / eLocation ID:
- 1565 to 1584
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Turbulence parameterization plays a critical role in the simulation of many weather regimes. For challenging cases such as the stratocumulus-capped boundary layer (SCBL), traditional schemes can produce unrealistic results even when a fine large-eddy-simulation (LES) resolution is used. Here we present an implicit generalized linear algebraic subfilter-scale model (iGLASS) to better represent unresolved turbulence in the simulation of the atmospheric boundary layer, at both standard LES and so-called terra incognita (TI) resolutions. The latter refers to a range of model resolutions where turbulent eddies are only partially resolved, and therefore the simulated processes are sensitive to the representation of unresolved turbulence. iGLASS is based on the truncated conservation equations of subfilter-scale (SFS) fluxes, but it integrates the full equations of the SFS turbulence kinetic energy and potential energy to retain “memory” of the SFS turbulence. Our evaluations suggest iGLASS can perform significantly better than traditional eddy-diffusivity models and exhibit skills comparable to the dynamic reconstruction model (DRM). For a neutral boundary layer case run at LES resolution, the simulation using iGLASS exhibits a wind profile that reasonably matches the similarity-theory solution. For an SCBL case with 5-m vertical resolution, iGLASS maintains more realistic cloud water profiles and boundary layer structure than traditional schemes. The SCBL case is also tested at TI resolution, and iGLASS also exhibits superior performance. iGLASS permits significant backscatter, whereas traditional models allow forward scatter (diffusion) only. As a physics-based approach, iGLASS appears to be a viable alternative for turbulence parameterization.
-
more » « less
Abstract Mesoscale climate models provide indispensable tools to understand land‐atmosphere interactions over urban regions. However, uncertainties in urban canopy parameters (UCPs) and parameterization schemes lead to degraded representation of the drag effect in complex built terrains. In particular, for the widely applied single‐layer urban canopy model (SLUCM) coupled with the Weather Research and Forecasting (WRF) model, near‐surface horizontal wind speed is known to be overestimated systematically. In this study, idealized large eddy simulations (LES) and WRF‐SLUCM simulations are conducted to study the separate effect of UCPs and aerodynamic parameterization on atmospheric boundary layer processes and rainfall variabilities in Phoenix, Arizona. For LES that explicitly resolves surface geometry, significant differences between three‐dimensional (3D) versus two‐dimensional (2D) representation of urban morphology are found in the surface layer and above. When surface drag is parameterized following SLUCM, surface morphologies have little impacts on the mean momentum transfer. WRF‐SLUCM simulation results, incorporated with 3D urban morphology data, indicate that simply refining the frontal area index will reduce the surface drag, which further amplifies the systematic positive bias of SLUCM in predicting horizontal wind speed. Replacing the drag parameterization in SLUCM by LES‐based aerodynamic parameters has evident impacts on near‐surface wind speed. The impact of urban roughness representation becomes the most evident during rainfall periods, due to the important role of surface drag in dictating moisture convergence. Our study underlines that apart from intensive efforts in obtaining detailed UCPs, it is also critical to enhance the urban momentum exchange parameterization schemes.
-
more » « less
The Weather Research and Forecasting (WRF) Model has been extensively used for wind energy applications, and current releases include a scheme that can be applied to examine the effects of wind turbine arrays on the atmospheric flow and electricity generation from wind turbines. Herein we present a high-resolution simulation using two different wind farm parameterizations: 1) the “Fitch” parameterization that is included in WRF releases and 2) the recently developed Explicit Wake Parameterization (EWP) scheme. We compare the schemes using a single yearlong simulation for a domain centered on the highest density of current turbine deployments in the contiguous United States (Iowa). Pairwise analyses are applied to diagnose the downstream wake effects and impact of wind turbine arrays on near-surface climate conditions. On average, use of the EWP scheme results in small-magnitude wake effects within wind farm arrays and faster recovery of full WT array wakes. This in turn leads to smaller impacts on near-surface climate variables and reduced array–array interactions, which at a systemwide scale lead to summertime capacity factors (i.e., the electrical power produced relative to nameplate installed capacity) that are 2%–3% higher than those from the more commonly applied Fitch parameterization. It is currently not possible to make recommendations with regard to which wind farm parameterization exhibits higher fidelity or to draw inferences with regard to whether the relative performance may vary with prevailing climate conditions and/or wind turbine deployment configuration. However, the sensitivities documented herein to the wind farm parameterization are of sufficient magnitude to potentially influence wind turbine array siting decisions. Thus, our research findings imply high value in undertaking combined long-term high-fidelity observational studies in support of model validation and verification.
-
more » « less
Abstract A field campaign at Siple Dome in West Antarctica during the austral summer 2019/20 offers an opportunity to evaluate climate model performance, particularly cloud microphysical simulation. Over Antarctic ice sheets and ice shelves, clouds are a major regulator of the surface energy balance, and in the warm season their presence occasionally induces surface melt that can gradually weaken an ice shelf structure. This dataset from Siple Dome, obtained using transportable and solar-powered equipment, includes surface energy balance measurements, meteorology, and cloud remote sensing. To demonstrate how these data can be used to evaluate model performance, comparisons are made with meteorological reanalysis known to give generally good performance over Antarctica (ERA5). Surface albedo measurements show expected variability with observed cloud amount, and can be used to evaluate a model’s snowpack parameterization. One case study discussed involves a squall with northerly winds, during which ERA5 fails to produce cloud cover throughout one of the days. A second case study illustrates how shortwave spectroradiometer measurements that encompass the 1.6-
μ m atmospheric window reveal cloud phase transitions associated with cloud life cycle. Here, continuously precipitating mixed-phase clouds become mainly liquid water clouds from local morning through the afternoon, not reproduced by ERA5. We challenge researchers to run their various regional or global models in a manner that has the large-scale meteorology follow the conditions of this field campaign, compare cloud and radiation simulations with this Siple Dome dataset, and potentially investigate why cloud microphysical simulations or other model components might produce discrepancies with these observations.significance statement Antarctica is a critical region for understanding climate change and sea level rise, as the great ice sheets and the ice shelves are subject to increasing risk as global climate warms. Climate models have difficulties over Antarctica, particularly with simulation of cloud properties that regulate snow surface melting or refreezing. Atmospheric and climate-related field work has significant challenges in the Antarctic, due to the small number of research stations that can support state-of-the-art equipment. Here we present new data from a suite of transportable and solar-powered instruments that can be deployed to remote Antarctic sites, including regions where ice shelves are most at risk, and we demonstrate how key components of climate model simulations can be evaluated against these data.
