Search

The parameterization of orographic gravity wave drag (OGWD) is essential for accurate numerical weather prediction in regions of complex terrain. Current OGWD schemes assume hydrostatic orographic gravity waves (OGWs) but the parameterized OGWs in fine-resolution models with narrow subgrid-scale orography can be significantly affected by the nonhydrostatic effects (NHE). In our recent work, the OGWD scheme in the Model for Prediction Across Scales (MPAS) was revised by accounting for the NHE on the surface wave momentum flux of upward-propagating OGWs. Herein, the revised OGWD scheme is implemented in the Weather Research and Forecasting (WRF) Model to evaluate its performance in short-range weather forecast. Two sets of 36-h WRF simulations are conducted for nine Northeast China cold vortices (NECVs) that occurred in the warm season of 2011 using the original and revised OGWD schemes. Results show that the WRF Model tends to underestimate the intensity of the NECVs, producing too high geopotential height. When accounting for the NHE in the OGWD scheme, the NECV intensity biases are significantly reduced. Analyses reveal that the NHE act to weaken the lower-tropospheric OGWD by decreasing the surface wave momentum flux, which strengthens the NECV in the lower troposphere. Consequently, the strengthened low-level cyclonic circulation increases the posttrough cold advection to the southwest of the NECV which in turn enhances the NECV in the mid–upper troposphere with reduced geopotential height. The NHE are found to increase as the model horizontal resolution increases, suggesting greater importance of NHE in the OGWD parameterization of high-resolution numerical models. [ABSTRACT FROM AUTHOR]

The snow albedo is a vital component of land–atmosphere coupling models. It plays a critical role in regulating land surface energy exchange by controlling incoming solar radiation absorbed by the land surface and influencing the timing and rate of snowmelt. Accurate snow albedo simulation is essential to obtain surface energy balance and snow-cover estimates. Here, the simulation of albedo and snow cover using the Weather Research and Forecasting Model and an improved snow albedo scheme is verified against satellite-retrieved products during and immediately following eight snowfall events over the Tibetan Plateau. The improved model successfully characterizes the spatial pattern and inverted U-shaped temporal pattern of albedo over the entire Tibetan Plateau. This is attributed to the local optimization of snow-age parameters and explicit consideration of snow depth in the improved scheme. Compared with the previous model, the model proposed herein greatly decreases the overestimated albedo (by 0.13–0.27), yielding a bias range of ±0.08, mean relative bias decrease of 70%, and significant increase in the spatial correlation coefficient of 0.03–0.39 (mean: 0.13). The significant improvements of albedo estimates appear in deep snow-covered regions, largely attributed to parameter optimization related to snow albedo decay, while less improvements appear over the shallow snow-covered regions. Accurate reproduction of the spatiotemporal variation in albedo alleviated snow-cover overestimation by small amounts. For snow-cover estimates, the improved model consistently decreases the false-alarm rate by 0.03, and increases the overall accuracy and equitable threat score by 0.04 and 0.03, respectively. Moreover, the improved scheme shows an equivalent improvement of albedo estimates at both 1- and 5-km grid spacing over the eastern Tibetan Plateau; this is also true for snow-cover estimates. Significance Statement: Snow albedo schemes in widely used numerical weather prediction models show notable shortcomings in complex mountainous regions, hindering accurate surface energy balance and snow-cover prediction. The purpose of this study is to better understand the role of snow albedo on snow-cover estimates and reveal the application potential of an improved snow albedo scheme across the Tibetan Plateau. This is important because snow albedo influences the timing and rate of snowmelt, and in turn snow-cover estimates, through regulating the surface energy budget. Our results highlight the strong application potential of our improved scheme in reducing snow simulation errors, confirm the importance of snow depth on snow albedo, and provide a new perspective for improving the accuracy of snow forecast over the topographically high Tibetan Plateau. [ABSTRACT FROM AUTHOR]

In this paper, we present the implementation and evaluate the impact of the new roughness length configuration in the ALARO canonical model configuration of the ALADIN system at the edge of the orographic gravity wave drag gray zone. As an essential input for turbulence parameterization, the roughness length affects the near-surface turbulent fluxes and the screen-level interpolation of meteorological parameters. We utilize GMTED2010 and ECOCLIMAP-II databases to derive orographic and vegetation components of the effective roughness length and introduce tuning parameters enabling us to optimize predicted near-surface turbulent momentum fluxes and 10-m wind speed. Based on sensitivity tests, we (i) prove the necessity of tuning the roughness length fields, (ii) considerably reduce the RMSE of near-surface turbulent momentum fluxes (6%–7%) and 10-m wind speed for different groups of stations (3%–10%), and (iii) identify the tree height as the most influential input parameter in our computational domain. The RMSE decomposition indicates that the improvement of 10-m wind speed mostly comes from a decrease in the random error and bias of the mean. The variability is slightly underestimated, thus reducing the model accuracy for wind speeds above the 95th percentile but at an acceptable level. We explain that roughness length tuning also compensates for the missing roughness sublayer correction in our system. Finally, we show that, although the impact of the orographic gravity wave drag scheme at a horizontal mesh size of 1.8 km is generally small, it is still beneficial for capturing some finer features observed in atmospheric soundings. Significance Statement: Aiming to improve the 10-m wind speed forecast without sacrificing the accuracy of turbulent momentum fluxes in the kilometric resolution numerical weather prediction model, we derived new roughness length fields from high-resolution physiography databases. Therein, we proved the importance of tuning the input orography and vegetation fields and, depending on the time of day and year, reduced the root-mean-square error of 10-m wind speed by 3%–10%. Further, we demonstrated that orographic gravity wave drag parameterization is still needed to predict finer details seen in wind profiles from atmospheric soundings. Finally, we discussed the related simplifications in our model and their implications and proposed steps toward a more consistent and complete treatment of the near-surface turbulence. [ABSTRACT FROM AUTHOR]

This study investigates the influences of initial soil moisture over the Tibetan Plateau (TP) on precipitation simulation, and the respective effects of boundary layer vertical diffusion for heat (Kh) and vapor (Kq). Results indicate that the responses of boundary layer vertical diffusion to soil moisture are obvious mainly in the daytime. Wetter land surface corresponds to weaker vertical diffusion, which could strengthen thermal forcing and dynamic lifting in the lower atmosphere, and encourage water vapor saturation near the top of boundary layer to prevent the environmental dry air entrainment/invasion, which would be beneficial to more convection and precipitation. Wetter land surface over the TP could enhance the contrast between the cold in the northwestern TP and the warm in the southeastern TP, which would be conducive to the southeastward propagation of precipitation. The simulation of heat and moisture in the boundary layer could be improved by perturbing the relative intensity of Kh and Kq. From the perspective of heat and moisture, Kh affects atmospheric stability, while Kq affects moisture and its vertical transport in the boundary layer. The Kh and Kq have competitive effects on precipitation intensity by influencing the relative importance of moisture and atmospheric stability conditions in the boundary layer. Adjusting the relative intensity of Kh and Kq would deactivate the competitive effects. Stronger Kh but weaker Kq would alleviate the overestimated precipitation by inhibiting vertical transport of moisture to the top of boundary layer and attenuating convective instability in the boundary layer. Significance Statement: The purpose of this study is to better understand the effects of boundary layer vertical heat and moisture diffusion in the response of precipitation to soil moisture. This is important because boundary layer vertical diffusion is a crucial factor influencing the relation between soil moisture and precipitation. Our results reveal the competitive effects of boundary layer vertical diffusion for heat and vapor on the simulation of precipitation. These results point a potential way toward better understanding the response of precipitation to soil moisture. [ABSTRACT FROM AUTHOR]

While tropical cyclone (TC) risk is a global concern, high regional differences exist in the quality of available data. This paper introduces InCyc, a globally consistent database of high-resolution full-physics simulations of historical cyclones. InCyc is designed to facilitate analysis of TC wind risk across basins and is made available to research institutions. We illustrate the value of this database with a case study focused on key wind risk parameters, namely, the location and intensity of peak winds for the North Atlantic and western North Pacific basins. A novel approach based on random forest algorithms is introduced to predict the full distribution of these TC wind risk parameters. Based on a leave-one-storm-out evaluation, the analysis of the predictions shows how this innovative approach compares to other parametric models commonly used for wind risk assessment. We finally discuss why capturing the full distribution of variability is crucial as well as the broader use in the context of TC risk assessment systems (i.e., "catastrophe models"). [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

To simulate the large-scale impacts of wind farms, wind turbines are parameterized within mesoscale models in which grid sizes are typically much larger than turbine scales. Five wind-farm parameterizations were implemented in the Weather Research and Forecasting (WRF) Model v4.3.3 to simulate multiple operational wind farms in the North Sea, which were verified against a satellite image, airborne measurements, and the FINO-1 meteorological mast data on 14 October 2017. The parameterization by Volker et al. underestimated the turbulence and wind speed deficit compared to measurements and to the parameterization of Fitch et al., which is the default in WRF. The Abkar and Porté-Agel parameterization gave close predictions of wind speed to that of Fitch et al. with a lower magnitude of predicted turbulence, although the parameterization was sensitive to a tunable constant. The parameterization by Pan and Archer resulted in turbine-induced thrust and turbulence that were slightly less than that of Fitch et al., but resulted in a substantial drop in power generation due to the magnification of wind speed differences in the power calculation. The parameterization by Redfern et al. was not substantially different from Fitch et al. in the absence of conditions such as strong wind veer. The simulations indicated the need for a turbine-induced turbulence source within a wind-farm parameterization for improved prediction of near-surface wind speed, near-surface temperature, and turbulence. The induced turbulence was responsible for enhancing turbulent momentum flux near the surface, causing a local speed-up of near-surface wind speed inside a wind farm. Our findings highlighted that wakes from large offshore wind farms could extend 100 km downwind, reducing downwind power production as in the case of the 400-MW Bard Offshore 1 wind farm whose power output was reduced by the wakes of the 402-MW Veja Mate wind farm for this case study. Significance Statement: Because wind farms are smaller than the common grid spacing of numerical weather prediction models, the impacts of wind farms on the weather have to be indirectly incorporated through parameterizations. Several approaches to parameterization are available and the most appropriate scheme is not always clear. The absence of a turbulence source in a parameterization leads to substantial inaccuracies in predicting near-surface wind speed and turbulence over a wind farm. The impact of large clusters of offshore wind turbines in the wind field can exceed 100 km downwind, resulting in a substantial loss of power for downwind turbines. The prediction of this power loss can be sensitive to the chosen parameterization, contributing to uncertainty in wind-farm economic planning. [ABSTRACT FROM AUTHOR]

NOAA has been developing a fully coupled Earth system model under the Unified Forecast System framework that will be responsible for global (ensemble) predictions at lead times of 0–35 days. The development has involved several prototype runs consisting of bimonthly initializations over a 7-yr period for a total of 168 cases. This study leverages these existing (baseline) prototypes to isolate the impact of substituting (one-at-a-time) parameterizations for convection, microphysics, and planetary boundary layer on 35-day forecasts. Through these physics sensitivity experiments, it is found that no particular configuration of the subseasonal-length coupled model is uniformly better or worse, based on several metrics including mean-state biases and skill scores for the Madden–Julian oscillation, precipitation, and 2-m temperature. Importantly, the spatial patterns of many "first-order" biases (e.g., impact of convection on precipitation) are remarkably similar between the end of the first week and weeks 3–4, indicating that some subseasonal biases may be mitigated through tuning at shorter time scales. This result, while shown for the first time in the context of subseasonal prediction with different physics schemes, is consistent with findings in climate models that some mean-state biases evident in multiyear averages can manifest in only a few days. An additional convective parameterization test using a different baseline shows that attempting to generalize results between or within modeling systems may be misguided. The limitations of generalizing results when testing physics schemes are most acute in modeling systems that undergo rapid, intense development from myriad contributors—as is the case in (quasi) operational environments. [ABSTRACT FROM AUTHOR]

Obtaining a faithful probabilistic depiction of moist convection is complicated by unknown errors in subgrid-scale physical parameterization schemes, invalid assumptions made by data assimilation (DA) techniques, and high system dimensionality. As an initial step toward untangling sources of uncertainty in convective weather regimes, we evaluate a novel Bayesian data assimilation methodology based on particle filtering within a WRF ensemble analysis and forecasting system. Unlike most geophysical DA methods, the particle filter (PF) represents prior and posterior error distributions nonparametrically rather than assuming a Gaussian distribution and can accept any type of likelihood function. This approach is known to reduce bias introduced by Gaussian approximations in low-dimensional and idealized contexts. The form of PF used in this research adopts a dimension-reduction strategy, making it affordable for typical weather applications. The present study examines posterior ensemble members and forecasts for select severe weather events between 2019 and 2020, comparing results from the PF with those from an ensemble Kalman filter (EnKF). We find that assimilating with a PF produces posterior quantities for microphysical variables that are more consistent with model climatology than comparable quantities from an EnKF, which we attribute to a reduction in DA bias. These differences are significant enough to impact the dynamic evolution of convective systems via cold pool strength and propagation, with impacts to forecast verification scores depending on the particular microphysics scheme. Our findings have broad implications for future approaches to the selection of physical parameterization schemes and parameter estimation within preexisting data assimilation frameworks. [ABSTRACT FROM AUTHOR]

AbstractThe role of nonlocal transport on the development and maintenance of marine stratocumulus (Sc) clouds in coarse-resolution models is investigated, with a special emphasis on the downdraft contribution. A new parameterization of cloud-top-triggered downdrafts is proposed and validated against large-eddy simulation (LES) for two Sc cases. The applied nonlocal mass-flux scheme is part of the stochastic multiplume eddy-diffusivity/mass-flux (EDMF) framework decomposing the turbulent transport into local and nonlocal contributions. The complementary local turbulent transport is represented with the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme. This EDMF version has been implemented in the Weather Research and Forecasting (WRF) single-column model (SCM) and tested for three model versions: without mass flux, with updrafts only, and with both updrafts and downdrafts. In the LES, the downdraft and updraft contributions to the total heat and moisture transport are comparable and significant. The WRF SCM results show a good agreement between the parameterized downdraft turbulent transport and LES. While including updrafts greatly improves the modeling of Sc clouds over the simulation without mass flux, the addition of downdrafts is less significant, although it helps improve the moisture profile in the planetary boundary layer.

Numerical weather prediction models contain parameters that are inherently uncertain and cannot be determined exactly. It is thus desirable to have reliable objective approaches for estimation of optimal values and uncertainties of these parameters. Traditionally, the parameter tuning has been done manually, which can lead to the tuning process being a maze of subjective choices. In this paper we present how to optimize 20 key physical parameters in the atmospheric model Open Integrated Forecasting System (OpenIFS) that have a strong impact on forecast quality. The results show that simultaneous optimization of O(20) parameters is possible with O(100) algorithm steps using an ensemble of O(20) members; the results also show that the optimized parameters lead to substantial enhancement of predictive skill. The enhanced predictive skill can be attributed to reduced biases in low-level winds and upper-tropospheric humidity in the optimized model. We find that the optimization process is dependent on the starting values of the parameters that are optimized (starting from better-suited values results in a better model). The results show also that the applicability of the tuned parameter values across different model resolutions is somewhat limited because of resolution-dependent model biases, and we also found that the parameter covariances provided by the tuning algorithm seem to be uninformative. Significance Statement: The purpose of this work is to show how to use algorithmic methods to optimize a weather model in a computationally efficient manner. Traditional manual model tuning is an extremely laborious and time-consuming process, so algorithmic methods have strong potential for saving the model developers' time and accelerating development. This paper shows that algorithmic optimization is possible and that weather forecasts can be improved. However, potential issues related to the use of the optimized parameter values across different model resolutions are discussed as well as other shortcomings related to the tuning process. [ABSTRACT FROM AUTHOR]

The momentum roughness length (z0) significantly impacts wind predictions in weather and climate models. Nevertheless, the impacts of z0 parameterizations in different wind regimes and various model configurations on the hurricane size, intensity, and track simulations have not been thoroughly established. To bridge this knowledge gap, a comprehensive analysis of 310 simulations of 10 real hurricanes using the Weather Research and Forecasting (WRF) Model is conducted in comparison with observations. Our results show that the default z0 parameterizations in WRF perform well for weak (category 1–2) hurricanes; however, they underestimate the intensities of strong (category 3–5) hurricanes. This finding is independent of model resolution or boundary layer schemes. The default values of z0 in WRF agree with the observational estimates from dropsonde data in weak hurricanes while they are much larger than observations in strong hurricanes regime. Decreasing z0 close to the values of observational estimates and theoretical hurricane intensity models in high wind regimes (≳45 m s−1) led to significant improvements in the intensity forecasts of strong hurricanes. A momentum budget analysis dynamically explained why the reduction of z0 (decreased surface turbulent stresses) leads to stronger simulated storms. [ABSTRACT FROM AUTHOR]

Traditional atmospheric surface layer theory assumes homogeneous surface conditions. Regardless, nearly all surface layer parameterization schemes employed within numerical weather prediction models utilize the same techniques within highly heterogeneous coastal regimes as for homogeneous environments. We compare predicted surface weather and fluxes of momentum, heat, and moisture—focusing mainly on momentum—from regional simulations using the Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) atmospheric model to observations collected from offshore buoys, inland flux towers, and radiosonde profiles during the Coastal Land-Air-Sea Interaction (CLASI) project throughout the summer of 2021 around Monterey Bay, California. Results reveal that modeled cross-coastal surface flux gradients are spuriously discontinuous, leading to systematically overestimated fluxes and weak winds inland of the coastline during onshore flow periods. Additionally, contrary to observations, modeled surface exchange coefficients are insensitive to wind direction on both sides of the coast, which degrades predictive skill downstream from the coastline. Over the central bay, prediction degrades when near-surface wind directions deviate from the prevailing flow direction as the parameterized stress–wind relationship fails during these cases. Predictive skill over the bay is therefore linked to variations in wind direction. Offshore of the geographically complex peninsula, systematic biases are less clear; however, bifurcations in drag coefficients based on wind direction were measured here as well. Last, increasing the horizontal grid spacing from 333 m to 3 km does not significantly affect surface layer prediction. This work highlights the need to reevaluate surface layer parameterization methods for modeling within coastal regions. Significance Statement: Understanding surface layer weather is critical for many purposes, such as infrastructure design and weather forecasting. Within the context of numerical modeling and weather prediction, skillful forecasts of surface winds and temperature rely on accurate portrayal of the surface layer. By comparing observations collected during the Coastal Land-Air-Sea Interaction field program to numerical model solutions, we show that prediction of the surface layer fluxes of momentum, heat, and moisture break down near the coastline, which leads to biases in the predicted surface layer weather both inland and over the water. As surface layer parameterization methods across nearly all numerical models are rooted in the same practices, our results call into question the use of traditional methods near the coastline. [ABSTRACT FROM AUTHOR]

This study investigates mixed-phase cloud (MPC) processes along the warm conveyor belts (WCBs) of two extratropical cyclones observed during the North Atlantic Waveguide and Downstream Impact Experiment (NAWDEX). The aim is to investigate the effect of two radically distinct parameterizations for MPCs on the WCB and the ridge building downstream: the first one (REF) drastically limits the formation of liquid clouds, while the second one (T40) forces the liquid clouds to exist. REF exhibits a stronger heating below 6-km height and a more important cooling above 6-km height than T40. The stronger heating at lower levels is due to more important water vapor depositional processes while the larger cooling at upper levels is due to differences in radiative cooling. The consequence is a more efficient potential vorticity destruction in the WCB outflow region and a more rapid ridge building in REF than T40. A comparison with airborne remote sensing measurements is performed. REF does not form any MPCs whereas T40 does, in particular in regions detected by the radar–lidar platform like below the dry intrusion. Comparison of both ice water content and reflectivity shows there may be too much pristine ice and not enough snow in REF and not enough cold hydrometeors in general in T40. The lower ice-to-snow ratio in T40 likely explains its better distribution of hydrometeors with respect to height compared to REF. These results underline the influence of MPC processes on the upper-tropospheric circulation and the need for more MPC observations in midlatitudes. Significance Statement: The diabatic processes occurring in the warm conveyor belt (WCB) of extratropical cyclones impact the jet stream structure at midlatitudes. This study highlights some sensitivity of upper-level dynamics to mixed-phase-cloud-related processes. Comparisons of two different microphysical schemes for mixed-phase clouds shows that the ratio of liquid to solid clouds along the WCB ascents impacts the latent heat release and the radiation. Data from the NAWDEX campaign helps to determine room for improvement for both schemes and point out the need of a better understanding of these processes for an improved prediction of upper-level dynamics. [ABSTRACT FROM AUTHOR]

The satellite observational data assimilation community requires consistent hydrometer descriptions—including mass–size relation and particle size distribution—to be used in both the forecast model and observation operator. We develop a microphysics-scheme-consistent snow and graupel single-scattering property database to meet this requirement. In this database, snowflakes are modeled as a mixture of small column and large aggregated ice particles, the mixing ratios of which may be adjusted to satisfy a given mass–size relation. Snow single-scattering properties are computed for four different mass–size relations. Subsequently, the snow description in the Thompson microphysics scheme is used as an example to demonstrate how microphysics-scheme-consistent snow bulk optical properties are derived. The Thompson-scheme-consistent snow bulk optical properties are added to the Community Radiative Transfer Model (CRTM), version 2.4.0. With CloudSat Cloud Profiling Radar (CPR) snow and liquid precipitation retrievals as the inputs, CRTM simulations are performed over global oceans and compared with four collocated Global Precipitation Measurement (GPM) Microwave Imager (GMI) high-frequency channel observations. The CRTM simulated brightness temperatures show agreement with the GMI observed brightness temperatures in cases of light-to-moderate precipitation over extratropical and polar ice-free oceans, with root-mean-square errors of 4.3, 13.0, 1.8, and 3.3 K in the 166-GHz (vertical polarization), 166-GHz (horizontal polarization), 183 ± 3-GHz (vertical polarization), and 183 ± 7-GHz (vertical polarization) channels, respectively. The result demonstrates the potential of using the newly developed microphysics-scheme-consistent snow optical parameterization in data assimilation applications. [ABSTRACT FROM AUTHOR]

During the winter of 2020/21 an ensemble of FV3-LAM forecasts was produced over the contiguous United States for the Winter Weather Experiment using five physics suites. These forecasts are evaluated with the goal of optimizing physics parameterizations within the future operational Rapid Refresh Forecast System (RRFS) in the Unified Forecast System (UFS) realm and for selecting suitable physics suites for a multiphysics RRFS ensemble. The five physics suites have different combinations of land surface models (LSMs), planetary boundary layer (PBL) parameterizations, and surface layer schemes, chosen from those used in current and possible future operational systems and likely to be supported in the operational UFS. Full-season evaluation reveals a persistent near-surface cold bias in the U.S. Northeast from one suite and a nighttime warm bias in the southern Great Plains in another suite, while other suites have smaller biases. A representative case is chosen to diagnose the cause for each of these biases using sensitivity simulations with different physics combinations or modified parameters and verified with additional mesonet observations. The cold bias in the Northeast is attributed to aspects of the Noah-MP LSM over snow cover, where Noah-MP simulates lower soil water content, and thus lower thermal conductivity than other LSMs, leading to less upward ground heat flux during nighttime and consequently lower surface temperature. The nighttime warm bias found in the southern Great Plains is attributed to overestimation of vertical mixing in the K-profile-based eddy-diffusivity mass-flux (K-EDMF) PBL scheme and insufficient land–atmospheric coupling from the GFS surface layer scheme over short vegetation. A few key parameters driving these systematic biases are identified. [ABSTRACT FROM AUTHOR]

Recent work has shown that bulk-Richardson (Rib) parameterizations for friction velocity, sensible heat flux, and latent heat flux have similar, and in some instances better, performance than long-standing parameterizations from Monin–Obukhov similarity theory (MOST). In this work, we expanded upon new Rib parameterizations and developed parameterizations of turbulence statistics, i.e., standard deviations in the 30-min u (horizontal), υ (meridional), and w (vertical) wind components (i.e., σu, συ, and σw, respectively), which allowed us to derive Rib-based parameterizations of turbulent kinetic energy (e), and standard deviations in the 30-min temperature and moisture measurements (σθ and σq, respectively). We used datasets from three 10-m micrometeorological towers installed during the Land Atmosphere Feedback Experiment (LAFE) conducted in Oklahoma from 1 to 31 August 2017 and evaluated the new parameterizations by comparing them against parameterizations from MOST. We used the LAFE datasets and fully independent datasets obtained from two micrometeorological towers installed in Alabama between February 2016 and April 2017 to evaluate the performance of the parameterizations. Based on the slope of the relationship between the observed and parameterized turbulence statistics (mb) and the coefficient of correlation (r), we found that the Rib relationships generally performed better than MOST at parameterizing συ, σw, σθ, and σq, and the Rib relationships performed better at low wind speeds than at high wind speeds. These results, coupled with recent developments of Rib parameterizations for surface-layer momentum, heat, and moisture fluxes, provide further evidence to consider using Rib-based parameterizations in weather forecasting models. Significance Statement: Deficiencies in Monin–Obukhov similarity theory (MOST) are well known, yet MOST forms the basis in weather forecasting models for describing heat, moisture, and momentum transfer between the land surface and atmosphere. We expanded upon previous work suggesting a MOST alternative called the bulk-Richardson approach. We used data collected from meteorological towers installed in Oklahoma and compared the bulk-Richardson approach with MOST. We evaluated these two approaches using data from meteorological towers installed in Oklahoma and Alabama and found that, overall, the bulk-Richardson approach performed better than MOST in determining the 30-min variability in temperature, moisture, and wind. This result provides additional motivation to use a bulk-Richardson approach in weather forecasting models because doing so will likely yield improved forecasts. [ABSTRACT FROM AUTHOR]

Numerical weather, climate, or Earth system models involve the coupling of components. At a broad level, these components can be classified as the resolved fluid dynamics, unresolved fluid dynamical aspects (i.e., those represented by physical parameterizations such as subgrid-scale mixing), and nonfluid dynamical aspects such as radiation and microphysical processes. Typically, each component is developed, at least initially, independently.Once development ismature, the components are coupled to deliver a model of the required complexity. The implementation of the coupling can have a significant impact on the model.As the error associated with each component decreases, the errors introduced by the coupling will eventually dominate. Hence, any improvement in one of the components is unlikely to improve the performance of the overall system. The challenges associated with combining the components to create a coherentmodel are here termed physics-dynamics coupling. The issue goes beyond the coupling between the parameterizations and the resolved fluid dynamics. This paper highlights recent progress and some of the current challenges. It focuses on three objectives: to illustrate the phenomenology of the coupling problemwith references to examples in the literature, to show howthe problem can be analyzed, and to create awareness of the issue across the disciplines and specializations. The topics addressed are different ways of advancing full models in time, approaches to understanding the role of the coupling and evaluation of approaches, coupling ocean and atmosphere models, thermodynamic compatibility between model components, and emerging issues such as those that arise as model resolutions increase and/ormodels use variable resolutions.

Wind-farm parameterizations in weather models can be used to predict both the power output and farm effects on the flow; however, their correctness has not been thoroughly assessed. We evaluate the wind-farm parameterization of the Weather Research and Forecasting Model with large-eddy simulations (LES) of the wake performed with the same model. We study the impact on the velocity and turbulence kinetic energy (TKE) of inflow velocity, roughness, resolution, number of turbines (one or two), and inversion height and strength. We compare the mesoscale with the LES by spatially averaging the LES within areas correspondent to the mesoscale horizontal spacing: one covering the turbine area and two downwind. We find an excellent agreement of the velocity within the turbine area between the two types of simulations. However, within the same area, we find the largest TKE discrepancies because in mesoscale simulations, the turbine-added TKE has to be highest at the turbine position to be advected downwind. Within the downwind areas, differences between velocities increase as the wake recovers faster in the LES, whereas for the TKE both types of simulations show similar levels. From the various configurations, the impact of inversion height and strength is small for these heights and inversion levels. The highest impact for the one-turbine simulations appears under the low-speed case due to the higher thrust, whereas the impact of resolution is low for the large-eddy simulations but high for the mesoscale simulations. Our findings demonstrate that higher-fidelity simulations are needed to validate wind-farm parameterizations. [ABSTRACT FROM AUTHOR]

The ability of a stochastically perturbed parameterization (SPP) approach to represent uncertainties in the model component of the Canadian Global Ensemble Prediction System was demonstrated in Part I of this investigation. The goal of this second step in SPP evaluation is to determine whether the scheme represents a viable alternative to the current operational combination of a multiphysics configuration and stochastically perturbed parameterization tendencies (SPPT). An assessment of the impact of each model uncertainty estimate in isolation reveals that, although the multiphysics configuration is highly effective at generating ensemble spread, it is often the result of differing biases rather than a reflection of flow-dependent error growth. Moreover, some of the members of the multiphysics ensemble suffer from large errors on regional scales as a result of suboptimal configurations. The SPP scheme generates a greater diversity of member solutions than the SPPT scheme in isolation, and it has an impact on forecast performance that is similar to that of current operational uncertainty estimates. When the SPP framework is combined with recent upgrades to the model physics suite that are only applicable in the stochastic perturbation context, the quality of global ensemble guidance is significantly improved. Significance Statement: The stochastically perturbed parameterization (SPP) technique was introduced in Part I to represent model uncertainties in forecasts generated by an operational global ensemble prediction system. We focus here on the viability of this technique as a replacement for the system's current uncertainty estimates: multiphysics and stochastic perturbations of physics tendencies. Despite the practical success of this combination, it suffers from physical inconsistencies and poor conservation properties. The adoption of SPP allows the ensemble to benefit from a recent set of model updates that couple with this new representation of model uncertainty to yield significant improvements in the quality of forecasts generated by the system. [ABSTRACT FROM AUTHOR]

Accurately representing model-based sources of uncertainty is essential for the development of reliable ensemble prediction systems for NWP applications. Uncertainties in discretizations, algorithmic approximations, and diabatic and unresolved processes combine to influence forecast skill in a flow-dependent way. An emerging approach designed to provide a process-level representation of these potential error sources, stochastically perturbed parameterizations (SPP), is introduced into the Canadian operational Global Ensemble Prediction System. This implementation extends the SPP technique beyond its typical application to free parameters in the physics suite by sampling uncertainty both within the dynamical core and at the formulation level using "error models" when multiple physical closures are available. Because SPP perturbs components within the model, internal consistency is ensured and conservation properties are not affected. The full SPP scheme is shown to increase ensemble spread to keep pace with error growth on a global scale. The sensitivity of the ensemble to each independently perturbed "element" is then assessed, with those responsible for the bulk of the response analyzed in more detail. Perturbations to surface exchange coefficients and the turbulent mixing length have a leading impact on near-surface statistics. Aloft, a tropically focused error model representing uncertainty in the advection scheme is found to initiate growing perturbations on the subtropical jet that lead to forecast improvements at higher latitudes. The results of Part I suggest that SPP has the potential to serve as a reliable representation of model uncertainty for ensemble NWP applications. Significance Statement: Ensemble systems account for the negative impact that uncertainties in prediction models have on forecasts. Here, uncertain model parameters and algorithms are subjected to perturbations representing impact on forecast errors. By initiating error growth within the model calculations, the equally skillful members of the ensemble remain physically realistic and self-consistent, which is not guaranteed by other depictions of model error. This "stochastically perturbed parameterization" technique (SPP) comprises many small error sources, each analyzed in isolation. Each source is related to a limited set of processes, making it possible to determine how the individual perturbations affect the forecast. We conclude that SPP in the Canadian Global Ensemble Forecasting System produces realistic estimates of the impact of model uncertainties on forecast skill. [ABSTRACT FROM AUTHOR]

Tropical cyclone (TC) intensity has been shown to have limited predictability in numerical weather prediction models; therefore, ensemble forecasting may be critical. An ensemble prediction system (EPS) should ideally cover all sources of uncertainty; however, most meso- and convective-scale EPSs typically consider initial-condition uncertainty alone, with limited treatment of model uncertainty, even though the evolution of mesoscale features is highly dependent on uncertain parameterization schemes. The role of stochastic treatment of model error in the Hurricane Weather Research and Forecasting (HWRF) EPS is evaluated by applying independent stochastically perturbed parameterization (iSPPT) scheme to individual parameterization schemes for four TCs from 2017 to 2018. Experiments with Hurricane Irma (2017) indicate that TC intensity ensemble standard deviation is most sensitive to the amplitude of the stochastic perturbation field, with smaller impact from adjusting the decorrelation time scale and spatial length scale. Results from all four TC cases show that stochastic perturbations to the turbulent mixing scheme can increase the ensemble standard deviation in intensity metrics over a 72-h simulation without introducing significant differences in mean error or bias. By contrast, stochastic perturbations to the microphysics, radiation, and cumulus tendencies have negligible effects on intensity standard deviation. [ABSTRACT FROM AUTHOR]

Modern general circulation models continue to require parameterizations of subgrid transport due to planetary boundary layer (PBL) turbulence and convection. Some schemes that unify these processes rely on assumed joint probability distributions of vertical velocity and moist conserved thermodynamic variables to predict the subgrid-scale contribution to the mean state of the atmosphere. The multivariate double-Gaussian mixture has been proposed as an appropriate model for PBL turbulence and shallow convection, but it is unable to reproduce important features of shallow cumulus convection. In this study, a novel unified PBL turbulence–convection–cloud macrophysics scheme is presented based on the eddy-diffusivity/mass-flux framework. The new scheme augments the double-Gaussian representation of subgrid variability with multiple stochastic mass-flux plumes at minimal added computational cost. Improved results for steady-state maritime and transient continental shallow convection from a single-column model implementation of the new scheme are shown with respect to reference large-eddy simulations. Improvements are seen in the cloud layer due to mass-flux plumes occupying the extreme moist, low liquid-water potential temperature tail of the joint temperature–moisture distribution. Significance Statement: Computer models of the atmosphere used to predict future climate are unable to directly represent air motion at small spatial scales because it would take too long to run the model over the entire planet. Instead, models typically use coarse model grid spacing and a simplified statistical representation of the physical processes that cause small-scale motions. This paper improves a particular simplified representation by adding a mechanism to represent statistically rare events of strong small-scale air motion that coherently transport air from near the surface to higher in the atmosphere. This increased transport also improves the representation of clouds, a particularly difficult phenomenon to simulate in models. [ABSTRACT FROM AUTHOR]

A new one-dimensional 1.5-order planetary boundary layer (PBL) scheme, based on the K–ε turbulence closure applied to the Reynolds-averaged Navier–Stokes (RANS) equations, is developed and implemented within the Weather Research and Forecasting (WRF) Model. The new scheme includes an analytic solution of the coupled equations for turbulent kinetic energy and dissipation rate. Different versions of the PBL scheme are proposed, with increasing levels of complexity, including a model for the calculation of the Prandtl number, a correction to the dissipation rate equation, and a prognostic equation for the temperature variance. Five different idealized cases are tested: four of them explore convective conditions, and they differ in initial thermal stratification and terrain complexity, while one simulates the very stable boundary layer case known as GABLS. For each case study, an ensemble of different large-eddy simulations (LES) is taken as reference for the comparison with the novel PBL schemes and other state-of-the-art 1- and 1.5-order turbulence closures. Results show that the new PBL K–ε scheme brings improvements in all the cases tested in this study. Specifically, the more significant are obtained with the turbulence closure including a prognostic equation for the temperature variance. Moreover, the largest benefits are obtained for the idealized cases simulating a typical thermal circulation within a two-dimensional valley. This suggests that the use of prognostic equations for dissipation rate and temperature variance, which take into account their transport and history, is particularly important with the increasing complexity of PBL dynamics. [ABSTRACT FROM AUTHOR]

nerating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δx) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δx ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional meso- scale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. [ABSTRACT FROM AUTHOR]

The forecast sensitivity to observations (FSO) is embedded into a new optimization framework for improving the observation performance in atmospheric data assimilation. Key ingredients are introduced as follows: the innovation-weight parameterization of the analysis equation, the FSO-based evaluation of the forecast error gradient to parameters, a line search approach to optimization, and an efficient mechanism for step length specification. This methodology is tested in preliminary numerical experiments with the Naval Research Laboratory Atmospheric Variational Data Assimilation System-Accelerated Representer (NAVDAS-AR) and the U.S. Navy's Global Environmental Model (NAVGEM) at a T425L60 resolution. The experimental setup relies on a verification state produced by the European Centre for Medium-Range Weather Forecasts (ECMWF) to estimate the analysis and short-range forecast errors. Parameter tuning is implemented in a training stage valid for 1–14 April 2018 and aimed at improving the use of assimilated observations in reducing the initial-condition errors. Assessment is carried out for 15 April–31 May 2018 to investigate the performance of the weighted assimilation system in reducing the errors in analyses and 24-h model forecasts. In average, as compared with the control run and verified against the ECMWF analyses, the weighted approach provided 17.4% reduction in analysis errors and 3.1% reduction in 24-h forecast errors, measured in a dry total energy norm. Observation impacts are calculated to assess the use of various observation types in reducing the analysis and forecast errors. In particular, assimilation of satellite wind data is significantly improved through the innovation-weighting procedure. Significance Statement: A new methodology is introduced to improve the information content of observations in numerical weather prediction (NWP). The computational procedure relies on observation sensitivity tools developed at all major NWP centers, and therefore, it appeals to a large audience for implementation and testing in practical applications. Our approach retains all available observations and provides a judicious optimization-based guidance to identify system deficiencies and improve the weighting assigned to various observing system components. The practical ability to implement this methodology is demonstrated in a computational environment that features all elements necessary for NWP applications. Preliminary results show that proper specification of the innovation weights can significantly improve the observation performance in reducing both the analysis errors and short-range forecast errors. [ABSTRACT FROM AUTHOR]

Modeling atmospheric turbulence in the convective boundary layer is challenging at kilometer and subkilometer resolutions, as the horizontal grid spacing approaches the size of the most energetic turbulent eddies. In this range of resolutions, termed terra incognita or gray zone, partially resolved convective structures are grid dependent and neither traditional 1D mesoscale parameterizations nor 3D large-eddy simulations closures are theoretically appropriate. Leveraging on a new set of one-way nested, full-physics multiscale numerical experiments, we quantify the magnitude of the errors introduced at gray zone resolutions in a real-case application and we provide new perspectives on recently proposed modeling approaches. The new set of experiments is forced by real-time-varying boundary conditions, spans a wide range of scales, and includes traditional 1D schemes, 3D closures, scale-aware parameterizations, and strategies to suppress resolved convection at gray zone resolutions. The study area is Riyadh (Saudi Arabia), where deep CBLs develop owing to strong convective conditions. Detailed analyses of our experiments, including validation with radiosonde data, calculations of spectral features, and partitioning of turbulent fluxes between resolved and subgrid scales, show that (i) grid-dependent convective structures entail minor impacts on the first-order characteristics of the fully developed boundary layer due to some degree of implicit scale awareness of 1D parameterizations and (ii) 3D closures and scale-aware schemes outperform traditional 1D schemes especially in the surface layer, among other findings. The new suite of experiments provides a benchmark of real simulations that can be extended to assess how new turbulence closures perform at gray zone resolutions. Significance Statement: As recent advances in high-performance computing are leading to a new era in numerical simulations, regional atmospheric models can now increase their resolution to the widely unexplored kilometer and subkilometer range. While increasing the resolution of atmospheric models is desirable to (i) have more realistic weather and air quality predictions and (ii) better represent boundary conditions for microscale models, kilometer and subkilometer grid spacings pose some theoretical challenges that need to be addressed by the atmospheric modeling community. In this work we run a set of numerical experiments for a real case study that aim to offer new perspectives on recently developed modeling strategies and identify the most promising directions that should be investigated by follow-up studies. [ABSTRACT FROM AUTHOR]

Recent studies have demonstrated that high-resolution (~25 km) Earth System Models (ESMs) have the potential to skillfully predict tropical cyclone (TC) occurrence and intensity. However, biases in ESM TCs still exist, largely due to the need to parameterize processes such as boundary layer (PBL) turbulence. Building on past studies, we hypothesize that the depiction of the TC PBL in ESMs is sensitive to the configuration of the PBL parameterization scheme, and that the targeted perturbation of tunable parameters can reduce biases. The Morris one-at-a-time (MOAT) method is implemented to assess the sensitivity of the TC PBL to tunable parameters in the PBL scheme in an idealized configuration of the Community Atmosphere Model, version 6 (CAM6). The MOAT method objectively identifies several parameters in an experimental version of the Cloud Layers Unified by Binormals (CLUBB) scheme that appreciably influence the structure of the TC PBL. We then perturb the parameters identified by the MOAT method within a suite of CAM6 ensemble simulations and find a reduction in model biases compared to observations and a high-resolution, cloud-resolving model. We demonstrate that the high-sensitivity parameters are tied to PBL processes that reduce turbulent mixing and effective eddy diffusivity, and that in CAM6 these parameters alter the TC PBL in a manner consistent with past modeling studies. In this way, we provide an initial identification of process-based input parameters that, when altered, have the potential to improve TC predictions by ESMs. [ABSTRACT FROM AUTHOR]

The stochastically perturbed parameterizations scheme (SPP) is here implemented and tested in HarmonEPS-the convection-permitting limited area ensemble prediction system by the international research program High Resolution Limited Area Model (HIRLAM) group. SPP introduces stochastic perturbations to values of chosen closure parameters representing efficiencies or rates of change in parameterized atmospheric (sub)processes. The impact of SPP is compared to that of the stochastically perturbed parameterization tendencies scheme (SPPT). SPP in this first version in HarmonEPS perturbs 11 parameters, active in different atmospheric processes and under various weather conditions. The main motivation for this study is the lack of variability seen in cloud products in HarmonEPS, as reported by duty forecasters. SPP in this first version is able to increase variability in a range of weather variables, including the cloud products. However, for some weather variables the root-mean-squared error of the ensemble mean is increased and the mean bias is impacted, especially in winter. This indicates that (some) parameter perturbation distributions are not optimal in the current configuration, and a further sensitivity analysis is required. SPPT resulted in a smaller increase in variability in the ensemble, but the impact was almost completely masked out when combined with perturbations of the initial state, lateral boundaries, and surface properties. An in-depth investigation into this lack of impact from SPPT is here presented through examining, among other things, accumulated tendencies from the model physics. [ABSTRACT FROM AUTHOR]

Stochastic model error schemes, such as the stochastic perturbed parameterization tendencies (SPPT) and independent SPPT (iSPPT) schemes, have become an increasingly accepted method to represent model error associated with uncertain subgrid-scale processes in ensemble prediction systems (EPSs). While much of the current literature focuses on the effects of these schemes on forecast skill, this research examines the physical processes by which iSPPT perturbations to the microphysics parameterization scheme yield variability in ensemble rainfall forecasts. Members of three 120-member Weather Research and Forecasting (WRF) Model ensemble case studies, including two distinct heavy rain events over Taiwan and one over the northeastern United States, are ranked according to an area-averaged accumulated rainfall metric in order to highlight differences between high- and low-precipitation forecasts. In each case, high-precipitation members are characterized by a damping of the microphysics water vapor and temperature tendencies over the region of heaviest rainfall, while the opposite is true for low-precipitation members. Physically, the perturbations to microphysics tendencies have the greatest impact at the cloud level and act to modify precipitation efficiency. To this end, the damping of tendencies in high-precipitation forecasts suppresses both the loss of water vapor due to condensation and the corresponding latent heat release, leading to grid-scale supersaturation. Conversely, amplified tendencies in low-precipitation forecasts yield both drying and increased positive buoyancy within clouds. [ABSTRACT FROM AUTHOR]

Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, sudden stratospheric warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These optimal forecasts concretely measure the event's progress. Direct numerical simulation can estimate them in principle but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are short compared to the time scale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass and modified by stochastic forcing. While highly idealized, this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in time scales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology. [ABSTRACT FROM AUTHOR]

Monin–Obukhov similarity theory (MOST) has long been used to represent surface–atmosphere exchange in numerical weather prediction (NWP) models. However, recent work has shown that bulk Richardson (Rib) parameterizations, rather than traditional MOST formulations, better represent near-surface wind, temperature, and moisture gradients. So far, this work has only been applied to unstable atmospheric regimes. In this study, we extended Rib parameterizations to stable regimes and developed parameterizations for the friction velocity (u*), sensible heat flux (H), and latent heat flux (E) using datasets from the Land-Atmosphere Feedback Experiment (LAFE). We tested our new Rib parameterizations using datasets from the Verification of the Origins of Rotation in Tornadoes Experiment-Southeast (VORTEX-SE) and compared the new Rib parameterizations with traditional MOST parameterizations and MOST parameterizations obtained using the LAFE datasets. We found that fitting coefficients in the MOST parameterizations developed from LAFE datasets differed from the fitting coefficients in classical MOST parameterizations which we attributed to the land surface heterogeneity present in the LAFE domain. Regardless, the new Rib parameterizations performed just as well as, and in some instances better than, the classical MOST parameterizations and the MOST parameterizations developed from the LAFE datasets. The improvement was most evident for H, particularly for H under unstable conditions, which was based on a better 1:1 relationship between the parameterized and observed values. These findings provide motivation to transition away from MOST and to implement bulk Richardson parameterizations into NWP models to represent surface–atmosphere exchange. Significance Statement: Models used to forecast the weather represent physical processes occurring in Earth's atmosphere using complex mathematical relationships. Many weather models use an approach developed in the 1950s referred to as Monin–Obukhov similarity theory, or MOST, to determine how heat, moisture, and wind change vertically above the land surface. However, MOST limitations are well known within the scientific community. We implemented a method different from MOST called the bulk Richardson approach. We found the bulk Richardson approach generally yielded predictions of how heat, moisture, and wind change with height near the ground that were just as good as, and in some instances better than, predictions of these quantities obtained using MOST. This result motivates us to consider implementing the bulk Richardson approach into weather models, as doing so is expected to lead to their improvement. [ABSTRACT FROM AUTHOR]

Coarse-grained results from a large-eddy simulation (LES) using the Weather Research and Forecasting (WRF) Model were compared in this study with the WRF simulations at a typical convection-permitting horizontal grid spacing of 3 km for an idealized case of deep moist convection. The purpose of this comparison is to identify major differences at the subgrid process level between two widely used deep convection parameterization schemes in the WRF Model. It is shown that there are considerable differences in subgrid process representations between the two schemes due to different parameterization formulations and underlying assumptions. The two schemes not only differ in trigger function, subgrid cloud model, and closure assumptions but also disagree with the coarse-grained LES results in terms of vertical mass flux profiles. Thus, it is difficult to discern which scheme is more advantageous over the other at the subgrid process level. The conclusions from this study highlight the importance of establishing benchmarks using observations and LES to develop and evaluate convection parameterization schemes suitable for models at convection-permitting resolution. [ABSTRACT FROM AUTHOR]

Recently, large-eddy simulation (LES) has been increasingly employed in meteorological simulations because it is a promising method for turbulent parameterization. However, it is still difficult to affirm that the numerical accuracy required for a dynamical core is fully understood. In this study, we derived two theoretical criteria for the order of accuracy of the advection term in a typical situation of the atmospheric boundary layer, and demonstrate their validity by numerical experiments. In the targeted grid-spacing of O(10) m, we determined the required order of accuracy as follows: Based on the criterion of the numerical diffusion error, the upwind scheme must have at least seventh-order accuracy. The fourth-order central scheme is barely acceptable with fourth-order explicit diffusion, provided that its coefficient is one or two orders of magnitude smaller than the implicit diffusion coefficient of the third-order upwind scheme. Based on the criterion of numerical dispersion error, at minimum, the seventh or eighth order is required. The dispersion error was indirect for the energy spectra, although we expect it may affect the local turbulence mechanism. We also investigated the effects of temporal discretization for compressible models, and found that relatively lower-order time schemes are available up to the O(10) m grid spacing if the time step is sufficiently small due to sound wave limitations. The importance of the derived criteria is that the required order of accuracy increases as the grid spacing decreases. This suggests that considerable care should be taken regarding the numerical error problem for future high-resolution LES. [ABSTRACT FROM AUTHOR]

We present a detailed evaluation of the turbulence forecast product eddy dissipation parameter (EDP) used at the Deutscher Wetterdienst (DWD). It is based on the turbulence parameterization scheme TURBDIFF, which is operational within the Icosahedral Nonhydrostatic (ICON) numerical weather prediction model used operationally by DWD. For aviation purposes, the procedure provides the cubic root of the eddy dissipation rate ε1/3 as an overall turbulence index. This quantity is a widely used measure for turbulence intensity as experienced by aircraft. The scheme includes additional sources of turbulent kinetic energy with particular relevance to aviation, which are briefly introduced. These sources describe turbulence generation by the subgrid-scale action of wake eddies, mountain waves, and convection, as well as horizontal shear as found close to fronts or the jet stream. Furthermore, we introduce a postprocessing calibration to an empirical EDR distribution, and we demonstrate the potential as well as limitations of the final EDP-based turbulence forecast by considering several case studies of typical turbulence events. Finally, we reveal the forecasting capability of this product by verifying the model results against one year of aircraft in situ EDR measurements from commercial aircraft. We find that the forecasted EDP performs favorably when compared to the Ellrod index. In particular, the turbulence signal from deep convection, which is accounted for in the EDP product, is advantageous when spatial nonlocality is allowed in the verification procedure. [ABSTRACT FROM AUTHOR]

The simulated winds within the urban canopy of landfalling tropical cyclones are sensitive to the representation of the planetary boundary and urban canopy layers in numerical weather prediction models. To assess the subgrid-scale parameterizations of these layers, mesoscale model simulations were executed and evaluated against near-surface observations as the outer wind field of Hurricane Irma (2017) interacted with the built-up region from downtown Miami northward to West Palm Beach. Four model simulations were examined, comprising two different planetary boundary layer (PBL) parameterizations (a local closure scheme with turbulent kinetic energy prediction and a nonlocal closure scheme) and two different urban canopy models (UCMs) [a zeroth-order bulk scheme and a multilayer building effect parameterization (BEP) that mimics the three-dimensionality of buildings]. Overall, the simulated urban canopy winds were weakly sensitive to the PBL scheme and strongly sensitive to the UCM. The bulk simulations compared most favorably to an analyzed wind swath in the urban environment, while the BEP simulations had larger negative biases in the same region. There is uncertainty in magnitude of the urban environment biases due to the lack of many urban sheltered measurements in the wind swath analysis. Biases in the rural environment were similar among the bulk and BEP simulations. An improved comparison with the analyzed wind swath in the urban region was obtained by reducing the drag coefficient in BEP in one of the PBL schemes. The usefulness of BEP was demonstrated in its ability to predict realistic heterogeneous near-surface velocity patterns in urban regions. [ABSTRACT FROM AUTHOR]

Subgrid-scale turbulence in numerical weather prediction models is typically handled by a PBL parameterization. These schemes attempt to represent turbulent mixing processes occurring below the resolvable scale of the model grid in the vertical direction, and they act upon temperature, moisture, and momentum within the boundary layer. This study varies the PBL mixing strength within 4-km WRF simulations of a 26-29 January 2015 snowstorm to assess the sensitivity of baroclinic cyclones to eddy diffusivity intensity. The bulk critical Richardson number for unstable regimes is varied between 0.0 and 0.25 within the YSU PBL scheme as a way of directly altering the depth and magnitude of subgridscale turbulent mixing. Results suggest that varying the bulk critical Richardson number is similar to selecting a different PBL parameterization. Differences in boundary layer moisture availability, arising from reduced entrainment of dry, free tropospheric air, lead to variations in the magnitude of latent heat release above the warm frontal region, producing stronger upper-tropospheric downstream ridging in simulations with less PBL mixing. The more amplified flow pattern impedes the northeastward propagation of the surface cyclone and results in a westward shift of precipitation. In addition, trajectory analysis indicates that ascending parcels in the less-mixing simulations condense more water vapor and terminate at a higher potential temperature level than do ascending parcels in the more-mixing simulations, suggesting stronger latent heat release when PBL mixing is reduced. These results suggest that spread within ensemble forecast systems may be improved by perturbing PBL mixing parameters that are not well constrained. [ABSTRACT FROM AUTHOR]

Monthlong simulations targeting four Madden–Julian oscillation events made with several global model configurations are verified against observations to assess the roles of grid spacing and convective parameterization on the representation of tropical convection and midlatitude forecast skill. Specifically, the performance of a global convection-permitting model (CPM) configuration with a uniform 3-km mesh is compared to that of a global 15-km mesh with and without convective parameterization, and of a variable-resolution "channel" simulation using 3-km grid spacing only in the tropics with a scale-aware convection scheme. It is shown that global 3-km simulations produce realistic tropical precipitation statistics, except for an overall wet bias and delayed diurnal cycle. The channel simulation performs similarly, although with an unrealistically higher frequency of heavy rain. The 15-km simulations with and without cumulus schemes produce too much light and heavy tropical precipitation, respectively. Without convection parameterization, the 15-km global model produces unrealistically abundant, short-lived, and intense convection throughout the tropics. Only the global CPM configuration is able to capture eastward-propagating Madden–Julian oscillation events, and the 15-km runs favor stationary or westward-propagating convection organized at the planetary scale. The global 3-km CPM exhibits the highest extratropical forecast skill aloft and at the surface, particularly during week 3 of each hindcast. Although more cases are needed to confirm these results, this study highlights many potential benefits of using global CPMs for subseasonal forecasting. Furthermore, results show that alternatives to global convection-permitting resolution—using coarser or spatially variable resolution—feature compromises that may reduce their predictive performance. [ABSTRACT FROM AUTHOR]

As wind farms grow in number and size worldwide, it is important that their potential impacts on the environment are studied and understood. The Fitch parameterization implemented in the Weather Research and Forecasting (WRF) Model since version 3.3 is a widely used tool today to study such impacts. We identified two important issues related to the way the added turbulent kinetic energy (TKE) generated by a wind farm is treated in the WRF Model with the Fitch parameterization. The first issue is a simple "bug" in the WRF code, and the second issue is the excessive value of a coefficient, called CTKE, that relates TKE to the turbine electromechanical losses. These two issues directly affect the way that a wind farm wake evolves, and they impact properties like near-surface temperature and wind speed at the wind farm as well as behind it in the wake. We provide a bug fix and a revised value of CTKE that is one-quarter of the original value. This 0.25 correction factor is empirical; future studies should examine its dependence on parameters such as atmospheric stability, grid resolution, and wind farm layout. We present the results obtained with the Fitch parameterization in the WRF Model for a single turbine with and without the bug fix and the corrected CTKE and compare them with high-fidelity large-eddy simulations. These two issues have not been discovered before because they interact with one another in such a way that their combined effect is a somewhat realistic vertical TKE profile at the wind farm and a realistic wind speed deficit in the wake. All WRF simulations that used the Fitch wind farm parameterization are affected, and their conclusions may need to be revisited. [ABSTRACT FROM AUTHOR]

Numerical models that are unable to resolve moist convection in the atmosphere employ physical parameterizations to represent the effects of the associated processes on the resolved-scale state. Most of these schemes are designed to represent the dominant class of cumulus convection that is driven by latent heat release in a conditionally unstable profile with a surplus of convective available potential energy (CAPE). However, an important subset of events occurs in low-CAPE environments in which potential and symmetric instabilities can sustain moist convective motions. Convection schemes that are dependent on the presence of CAPE are unable to depict accurately the effects of cumulus convection in these cases. A mass-flux parameterization is developed to represent such events, with triggering and closure components that are specifically designed to depict subgrid-scale convection in low-CAPE profiles. Case studies show that the scheme eliminates the "bull's-eyes" in precipitation guidance that develop in the absence of parameterized convection, and that it can represent the initiation of elevated convection that organizes squall-line structure. The introduction of the parameterization leads to significant improvements in the quality of quantitative precipitation forecasts, including a large reduction in the frequency of spurious heavy-precipitation events predicted by the model. An evaluation of surface and upper-air guidance shows that the scheme systematically improves the model solution in the warm season, a result that suggests that the parameterization is capable of accurately representing the effects of moist convection in a range of low-CAPE environments. [ABSTRACT FROM AUTHOR]

Proper behavior of physics parameterizations in numerical models at grid sizes of order 1 km is a topic of current research. Modifications to parameterization schemes to accommodate varying grid sizes are termed ‘‘scale aware.’’ The general problem of grids on which a physical process is partially resolved is called the ‘‘gray zone’’ or ‘‘terra incognita.’’ Here we examine features of the Mellor–Yamada–Nakanishi–Niino (MYNN) boundary layer scheme with eddy diffusivity and mass flux (EDMF) that were intended to provide scale awareness, as implemented in WRF, version 4.1. Scale awareness is provided by reducing the intensity of nonlocal components of the vertical mixing in the scheme as the grid size decreases. However, we find that the scale-aware features cause poorer performance in our tests on a 600-m grid. The resolved circulations on the 600-m grid have different temporal and spatial scales than are found in large-eddy simulations of the same cases, for reasons that are well understood theoretically and are described in the literature. The circulations [model convectively induced secondary circulations (M-CISCs)] depend on the grid size and on details of the model numerics. We conclude that scale awareness should be based on effective resolution, and not on grid size, and that the gray-zone problem for boundary layer turbulence and shallow cumulus cannot be solved simply by reducing the intensity of the parameterization. Parameterizations with different characteristics may lead to different conclusions. [ABSTRACT FROM AUTHOR]

Subgrid-scale (SGS) parameterizations in atmosphere and ocean models are often defined independently in the horizontal and vertical directions because the grid spacing is not the same in these directions (anisotropic grids). In this paper, we introduce a new anisotropic SGS model in large-eddy simulations (LES) of stratified turbulence based on horizontal filtering of the equations of motion. Unlike the common horizontal SGS parameterizations in atmosphere and ocean models, the vertical derivatives of the horizontal SGS fluxes are included in our anisotropic SGS scheme, and therefore the horizontal and vertical SGS dissipation mechanisms are not disconnected in the newly developed model. Our model is tested with two vertical grid spacings and various horizontal resolutions, where the horizontal grid spacing is comparatively larger than that in the vertical. Our anisotropic LES model can successfully reproduce the results of direct numerical simulations, while the computational cost is significantly reduced in the LES. We suggest the new anisotropic SGS model as an alternative to current SGS parameterizations in atmosphere and ocean models, in which the schemes for horizontal and vertical scales are often decoupled. The new SGS scheme may improve the dissipative performance of atmosphere and ocean models without adding any backscatter or other energizing terms at small horizontal scales. [ABSTRACT FROM AUTHOR]

We develop a methodology for identification of candidate observables that best constrain the parameterization of physical processes in numerical models. This methodology consists of three steps: (i) identifying processes that significantly impact model results, (ii) identifying observables that best constrain the influential processes, and (iii) investigating the sensitivity of the model results to the measurement error and vertical resolution of the constraining observables. This new methodology is applied to the Jet Propulsion Laboratory stochastic multiplume Eddy-Diffusivity/Mass-Flux (JPL-EDMF) model for two case studies representing nonprecipitating marine stratocumulus and marine shallow convection. The uncertainty of physical processes is characterized with uncertainty of model parameters. We find that the most uncertain processes in the JPL-EDMF model are related to the representation of lateral entrainment for convective plumes and parameterization of mixing length scale for the eddy-diffusivity part of the model. The results show a strong interaction between these uncertain processes. Measurements of the water vapor profile for shallow convection and of the cloud fraction profile for the stratocumulus case are among those measurements that best constrain the uncertain JPL-EDMF processes. The interdependence of the required vertical resolution and error characteristics of the observational system are shown. If the observations are associated with larger error, their vertical resolution has to be finer and vice versa. We suggest that the methodology and results presented here provide an objective basis for defining requirements for future observing systems such as future satellite missions to observe clouds and the planetary boundary layer. [ABSTRACT FROM AUTHOR]

Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the It^o correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the It^o correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the It^o correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized It^ o correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized It^o correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation. [ABSTRACT FROM AUTHOR]

In a recent work, a sensitivity analysis methodology was described that allows for a visual display of forecast sensitivity, with respect to model parameters, across a gridded forecast field. In that approach, sensitivity was assessed with respect to model parameters that are continuous in nature. Here, the analogous methodology is developed for situations involving noncontinuous (discrete or categorical) model parameters. The method is variance based, and the variances are estimated via a random-effects model based on 2k−p fractional factorial designs and Graeco-Latin square designs. The development is guided by its application to model parameters in the stochastic kinetic energy backscatter scheme (SKEBS), which control perturbations at unresolved, subgrid scales. In addition to the SKEBS parameters, the effect of daily variability and replication (both, discrete factors) are also examined. The forecasts examined are for precipitation, temperature, and wind speed. In this particular application, it is found that the model parameters have a much weaker effect on the forecasts as compared to the effect of daily variability and replication, and that sensitivities, weak or strong, often have a distinctive spatial structure that reflects underlying topography and/or weather patterns. These findings caution against fine-tuning methods that disregard 1) sources of variability other than those due to model parameters, and 2) spatial structure in the forecasts. [ABSTRACT FROM AUTHOR]

To accurately simulate the atmospheric state within the planetary boundary layer (PBL) by PBL parameterization scheme in different regions with their dominant weather/climate regimes is important for global/regional atmospheric models. In this study, we introduce the turbulence kinetic energy (TKE) and TKE dissipation rate (ε) based 1.5-order closure PBL parameterization (E–ε, EEPS) in the Weather Research and Forecasting (WRF) Model. The performances of the newly implemented EEPS scheme and the existing Yonsei University (YSU) scheme, the University of Washington (UW) scheme, and Mellor–Yamada–Nakanishi–Niino (MYNN) scheme are evaluated over the stratocumulus dominated southeast Pacific (SEP) and over the Southern Great Plains (SGP) where strong PBL diurnal variation is common. The simulations by these PBL parameterizations are compared with various observations from two field campaigns: the Variability of American Monsoon Systems Project (VAMOS) Ocean–Cloud–Atmosphere–Land Study (VOCALS) in 2008 over the SEP and the Land–Atmosphere Feedback Experiment (LAFE) in 2017 over the SGP. Results show that the EEPS and YSU schemes perform comparably over both regions, while the MYNN scheme performs differently in many aspects, especially over the SEP. The EEPS (MYNN) scheme slightly (significantly) underestimates liquid water path over the SEP. Compared with observations, the UW scheme produces the best PBL height over the SEP. The MYNN produces too high PBL height over the western part of the SEP while both the YSU and EEPS schemes produce too low PBL and cloud-top heights. The differences among the PBL schemes in simulating the PBL features over the SGP are relatively small. [ABSTRACT FROM AUTHOR]

There is a longstanding challenge in numerical weather and climate prediction to accurately model tropical wave variability, including convectively coupled equatorial waves (CCEWs) and the Madden–Julian oscillation. For subseasonal prediction, the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS) has been shown to be superior to the NOAA Global Forecast System (GFS) in simulating tropical variability, suggesting that the ECMWF model is better at simulating the interaction between cumulus convection and the large-scale tropical circulation. In this study, we experiment with the cumulus convection scheme of the ECMWF IFS in a research version of the GFS to understand which aspects of the IFS cumulus convection scheme outperform those of the GFS convection scheme in the tropics. We show that the IFS cumulus convection scheme produces significantly different tropical moisture and temperature tendency profiles from those simulated by the GFS convection scheme when it is coupled with other physics schemes in the GFS physics package. We show that a consistent treatment of the interaction between parameterized convective plumes in the GFS planetary boundary layer (PBL) and the IFS convection scheme is required for the GFS to replicate the tropical temperature and moisture profiles simulated by the IFS model. The GFS model with the IFS convection scheme, and the consistent treatment between the convection and PBL schemes, produces much more organized convection in the tropics, and generates tropical waves that propagate more coherently than the GFS in its default configuration due to better simulated interaction between low-level convergence and precipitation. [ABSTRACT FROM AUTHOR]

The parameterization of deep moist convection as a subgrid-scale process in numerical models of the atmosphere is required at resolutions that extend well into the convective "gray zone," the range of grid spacings over which such convection is partially resolved. However, as model resolution approaches the gray zone, the assumptions upon which most existing convective parameterizations are based begin to break down. We focus here on one aspect of this problem that emerges as the temporal and spatial scales of the model become similar to those of deep convection itself. The common practice of static tendency application over a prescribed adjustment period leads to logical inconsistencies at resolutions approaching the gray zone, while more frequent refreshment of the convective calculations can lead to undesirable intermittent behavior. A proposed parcel-based treatment of convective initiation introduces memory into the system in a manner that is consistent with the underlying physical principles of convective triggering, thus reducing the prevalence of unrealistic gradients in convective activity in an operational model running with a 10 km grid spacing. The subsequent introduction of a framework that considers convective clouds as persistent objects, each possessing unique attributes that describe physically relevant cloud properties, appears to improve convective precipitation patterns by depicting realistic cloud memory, movement, and decay. Combined, this Lagrangian view of convection addresses one aspect of the convective gray zone problem and lays a foundation for more realistic treatments of the convective life cycle in parameterization schemes. [ABSTRACT FROM AUTHOR]

Kilometer-scale models allow for an explicit simulation of deep convective overturning but many subgrid processes that are crucial for convective initiation are still poorly represented. This leads to biases such as insufficient convection triggering and late peak of summertime convection. A physically based stochastic perturbation scheme (PSP) for subgrid processes has been proposed (Kober and Craig) that targets the coupling between subgrid turbulence and resolved convection. The first part of this study presents four modifications to this PSP scheme for subgrid turbulence: an autoregressive, continuously evolving random field; a limitation of the perturbations to the boundary layer that removes artificial convection at night; a mask that turns off perturbations in precipitating columns to retain coherent structures; and nondivergent wind perturbations that drastically increase the effectiveness of the vertical velocity perturbations. In a revised version, PSP2, the combined modifications retain the physically based coupling to the boundary layer scheme of the original scheme while removing undesirable side effects. This has the potential to improve predictions of convective initiation in kilometer-scale models while minimizing other biases. The second part of the study focuses on perturbations to account for convective initiation by subgrid orography. Here the mechanical lifting effect is modeled by introducing vertical and horizontal wind perturbations of an orographically induced gravity wave. The resulting perturbations lead to enhanced convective initiation over mountainous terrain. However, the total benefit of this scheme is unclear and we do not adopt the scheme in our revised configuration. [ABSTRACT FROM AUTHOR]

This study aims to investigate the impact of increased vertical resolution on global medium-range forecasts. For this purpose, the dependencies of simulated atmospheric temperature and specific humidity on the lowest model level height and vertical grid spacing are investigated. The reduced first model level increases vertical turbulent mixing by determining a higher planetary boundary layer (PBL) height and associated turbulent diffusivities and velocity scales. This contributes to warming/drying within the PBL and cooling/moistening above the PBL. Resulting dryness near the surface enhances an increase in surface moisture flux from the ocean, which results in the apparent moistening of the troposphere. Consequently, large-scale precipitation increases due to more humid atmospheric conditions, while convective precipitation decreases because of drier conditions near the starting level of convection. Meanwhile, the reduced vertical grid spacing resolves the overshooting layer well in the cumulus convection process, which decreases the detrained moisture at the convective cloud top. This leads to a noticeable downward shift of ice clouds in the upper troposphere, and further contributes to the enhancement of longwave cooling via cloud–radiation processes. In medium-range forecasts, the increased vertical resolution exerts a significant impact on the simulated features of tropospheric temperature and humidity, while changes in the prediction accuracy precipitation are negligible owing to compensation between convective and large-scale precipitation. Finally, the discussion of possible methods to minimize the sensitivities of model's physics to vertical resolution is presented. [ABSTRACT FROM AUTHOR]