Stochastic Kinetic Energy Backscatter Scheme Research Articles

Enhanced regional ocean ensemble data assimilation through atmospheric coupling in the SKRIPS model

We investigate the impact of ocean data assimilation using the Ensemble Adjustment Kalman Filter (EAKF) from the Data Assimilation Research Testbed (DART) on the oceanic and atmospheric states of the Red Sea. Our study extends the ocean data assimilation experiment performed by Sanikommu et al. (2020) by utilizing the SKRIPS model coupling the MITgcm ocean model and the Weather Research and Forecasting (WRF) atmosphere model. Using a 50-member ensemble, we assimilate satellite-derived sea surface temperature and height and in situ temperature and salinity profiles every three days for one year, starting January 01 2011. Atmospheric data are not assimilated in the experiments. To improve the ensemble realism, perturbations are added to the WRF model using several physics options and the stochastic kinetic energy backscatter (SKEB) scheme. Compared with the control experiments using uncoupled MITgcm with ECMWF ensemble forcing, the EAKF ensemble mean oceanic states from the coupled model are better or insignificantly worse (root-mean-square errors are 23% to −1.3% smaller), especially when the atmospheric model uncertainties are accounted for with stochastic perturbations. We hypothesize that the ensemble spreads of the air–sea fluxes are better represented in the downscaled WRF ensembles when uncertainties are well accounted for, leading to improved representation of the ensemble oceanic states from the new experiments with the coupled model. This indicates the ocean model assimilation will be improved with coupled models and may relax the need for operational centers to provide atmospheric ensembles to drive ocean forecasts. Although the feedback from ocean to atmosphere is included in this two-way regional coupled configuration, we find no significant effect of ocean data assimilation on the ensemble mean latent heat flux and 10-m wind speed over the Red Sea. This suggests that the improved skill using the coupled model is not from the two-way coupling, but from downscaling the ensemble atmospheric forcings (one-way coupled) to drive the ocean model.

Testing Stochastic and Perturbed Parameter Methods in an Experimental 1-km Warn-on-Forecast System Using NSSL’s Phased-Array Radar Observations

Abstract The success of the National Severe Storms Laboratory’s (NSSL) experimental Warn-on-Forecast System (WoFS) to provide useful probabilistic guidance of severe and hazardous weather is mostly due to the frequent assimilation of observations, especially radar observations. Phased-array radar (PAR) technology, which is a potential candidate to replace the current U.S. operational radar network, would allow for even more rapid assimilation of radar observations by providing full-volumetric scans of the atmosphere every ∼1 min. Based on previous studies, more frequent PAR data assimilation can lead to improved forecasts, but it can also lead to ensemble underdispersion and suboptimal observation assimilation. The use of stochastic and perturbed parameter methods to increase ensemble spread is a potential solution to this problem. In this study, four stochastic and perturbed parameter methods are assessed using a 1-km-scale version of the WoFS and include the stochastic kinetic energy backscatter (SKEB) scheme, the physically based stochastic perturbation (PSP) scheme, a fixed perturbed parameters (FPP) method, and a novel surface-model scheme blending (SMSB) method. Using NSSL PAR observations from the 9 May 2016 tornado outbreak, experiments are conducted to assess the impact of the methods individually, in different combinations, and with different cycling intervals. The results from these experiments reveal the potential benefits of stochastic and perturbed parameter methods for future versions of the WoFS. Stochastic and perturbed parameter methods can lead to more skillful forecasts during periods of storm development. Moreover, a combination of multiple methods can result in more skillful forecasts than using a single method. Significance Statement Phased-array radar technology allows for more frequent assimilation of radar observations into ensemble forecast systems like the experimental Warn-on-Forecast System. However, more frequent radar data assimilation can eventually cause issues for prediction systems due to the lack of ensemble spread. Thus, the purpose of this study is to explore the use of four stochastic and perturbed parameter methods in a next-generation Warn-on-Forecast System to generate ensemble spread and help prevent the issues from frequent radar data assimilation. Results from this study indicate the stochastic and perturbed parameter methods can improve forecasts of storms, especially during storm development.

The Impact of Stochastic Perturbations in Physics Variables for Predicting Surface Solar Irradiance

We present a probabilistic framework tailored for solar energy applications referred to as the Weather Research and Forecasting-Solar ensemble prediction system (WRF-Solar EPS). WRF-Solar EPS has been developed by introducing stochastic perturbations into the most relevant physical variables for solar irradiance predictions. In this study, we comprehensively discuss the impact of the stochastic perturbations of WRF-Solar EPS on solar irradiance forecasting compared to a deterministic WRF-Solar prediction (WRF-Solar DET), a stochastic ensemble using the stochastic kinetic energy backscatter scheme (SKEBS), and a WRF-Solar multi-physics ensemble (WRF-Solar PHYS). The performances of the four forecasts are evaluated using irradiance retrievals from the National Solar Radiation Database (NSRDB) over the contiguous United States. We focus on the predictability of the day-ahead solar irradiance forecasts during the year of 2018. The results show that the ensemble forecasts improve the quality of the forecasts, compared to the deterministic prediction system, by accounting for the uncertainty derived by the ensemble members. However, the three ensemble systems are under-dispersive, producing unreliable and overconfident forecasts due to a lack of calibration. In particular, WRF-Solar EPS produces less optically thick clouds than the other forecasts, which explains the larger positive bias in WRF-Solar EPS (31.7 W/m2) than in the other models (22.7–23.6 W/m2). This study confirms that the WRF-Solar EPS reduced the forecast error by 7.5% in terms of the mean absolute error (MAE) compared to WRF-Solar DET, and provides in-depth comparisons of forecast abilities with the conventional scientific probabilistic approaches (i.e., SKEBS and a multi-physics ensemble). Guidelines for improving the performance of WRF-Solar EPS in the future are provided.

Multivariate ensemble sensitivity analysis applied for an extreme rainfall over Indian subcontinent

Ensemble sensitivity analysis (ESA) uses sample statistics of ensemble forecasts to estimate relationships between forecast metrics and initial conditions. The ensemble sensitivity analysis is often considered as a simple univariate regression as it includes an approximation of analysis covariance matrix with corresponding diagonal elements. In this work, univariate ensemble sensitivity is extended to multivariate ensemble sensitivity that incorporates the contribution from the full covariance matrix in the sensitivity calculations. The performance of multivariate ensemble sensitivity over univariate is examined for meso- and convective scale ensemble forecasts of a heavy rainfall event that happened over the Chennai city in India in December 2015. The ensemble forecasts and analyses are generated using the Advanced Research - Weather Research and Forecasting (WRF) model Data Assimilation Research Testbed (DART) based Ensemble Kalman Filter (EnKF). Multivariate ensemble sensitivity shows organized sensitivity patterns, while the sensitivity values are found to be broadly distributed in univariate ensemble sensitivity. Both the methods are validated using a perturbed initial condition approach, and the results indicate that the multivariate ensemble sensitivity method is effective in predicting the forecast response closest to the actual model response compared to the univariate ensemble sensitivity. The impact of model error on sensitivity calculations is examined by generating a new set of ensembles that uses the Stochastic Kinetic Energy Backscatter Scheme (SKEBS). In the presence of added model error, the forecast response estimated by multivariate using SKEBS ensembles compares better with the actual response. It is found that the performance of the multivariate approach depends on the optimal choice of localization radius, and if insufficient localization is applied, the spurious long-distance correlation contaminates the performance of the multivariate ensemble sensitivity method. The impact of various forecast lead times on the univariate and multivariate ensemble sensitivity analysis indicates that responses using multivariate ensemble sensitivity are more accurate than univariate, especially at longer lead times when nonlinearity becomes significant. The performance of univariate and multivariate methods in convection-permitting scale is examined by using the high-resolution ensemble forecasts, and it is found that the multivariate sensitivity with localization substantially improves the estimates in finer scales.

The Development of the NCEP Global Ensemble Forecast System Version 12

Abstract The Global Ensemble Forecast System (GEFS) is upgraded to version 12, in which the legacy Global Spectral Model (GSM) is replaced by a model with a new dynamical core—the Finite Volume Cubed-Sphere Dynamical Core (FV3). Extensive tests were performed to determine the optimal model and ensemble configuration. The new GEFS has cubed-sphere grids with a horizontal resolution of about 25 km and an increased ensemble size from 20 to 30. It extends the forecast length from 16 to 35 days to support subseasonal forecasts. The stochastic total tendency perturbation (STTP) scheme is replaced by two model uncertainty schemes: the stochastically perturbed physics tendencies (SPPT) scheme and stochastic kinetic energy backscatter (SKEB) scheme. Forecast verification is performed on a period of more than two years of retrospective runs. The results show that the upgraded GEFS outperforms the operational-at-the-time version by all measures included in the GEFS verification package. The new system has a better ensemble error–spread relationship, significantly improved skills in large-scale environment forecasts, precipitation probability forecasts over CONUS, tropical cyclone track and intensity forecasts, and significantly reduced 2-m temperature biases over North America. GEFSv12 was implemented on 23 September 2020.

Potential of the Coupled WRF/WRF-Hydro Modeling System for Flood Forecasting in the Ouémé River (West Africa)

Since the beginning of the 2000s, most of the West-African countries, particularly Benin, have experienced an increased frequency of extreme flood events. In this study, we focus on the case of the Ouémé river basin in Benin. To investigate flood events in this basin for early warning, the coupled atmosphere–hydrology model system WRF-Hydro is used, and analyzed for the period 2008–2010. Such a coupled model allows exploration of the contribution of atmospheric components into the flood event, and its ability to simulate and predict accurate streamflow. The potential of WRF-Hydro to correctly simulate streamflow in the Ouémé river basin is assessed by forcing the model with operational analysis datasets from the European Centre for Medium-Range Weather Forecasts (ECMWF). Atmospheric and land surface processes are resolved at a spatial resolution of 5 km. The additional surface and subsurface water flow routing are computed at a resolution of 500 m. Key parameters of the hydrological module of WRF-Hydro are calibrated offline and tested online with the coupled WRF-Hydro. The uncertainty of atmospheric modeling on coupled results is assessed with the stochastic kinetic energy backscatter scheme (SKEBS). WRF-Hydro is able to simulate the discharge in the Ouémé river in offline and fully coupled modes with a Kling–Gupta efficiency (KGE) around 0.70 and 0.76, respectively. In the fully coupled mode, the model captures the flood event that occurred in 2010. A stochastic perturbation ensemble of ten members for three rain seasons shows that the coupled model performance in terms of KGE ranges from 0.14 to 0.79. Additionally, an assessment of the soil moisture has been developed. This ability to realistically reproduce observed discharge in the Ouémé river basin demonstrates the potential of the coupled WRF-Hydro modeling system for future flood forecasting applications.

Lateral terrestrial water flow contribution to summer precipitation at continental scale – A comparison between Europe and West Africa with WRF‐Hydro‐tag ensembles

AbstractIt is well accepted that summer precipitation can be altered by soil moisture condition. Coupled land surface – atmospheric models have been routinely used to quantify soil moisture – precipitation feedback processes. However, most of the land surface models (LSMs) assume a vertical soil water transport and neglect lateral terrestrial water flow at the surface and in the subsurface, which potentially reduces the realism of the simulated soil moisture – precipitation feedback. In this study, the contribution of lateral terrestrial water flow to summer precipitation is assessed in two different climatic regions, Europe and West Africa, for the period June–September 2008. A version of the coupled atmospheric‐hydrological model WRF‐Hydro with an option to tag and trace land surface evaporation in the modelled atmosphere, named WRF‐Hydro‐tag, is employed. An ensemble of 30 simulations with terrestrial routing and 30 simulations without terrestrial routing is generated with random realizations of turbulent energy with the stochastic kinetic energy backscatter scheme, for both Europe and West Africa. The ensemble size allows to extract random noise from continental‐scale averaged modelled precipitation. It is found that lateral terrestrial water flow increases the relative contribution of land surface evaporation to precipitation by 3.6% in Europe and 5.6% in West Africa, which enhances a positive soil moisture – precipitation feedback and generates more uncertainty in modelled precipitation, as diagnosed by a slight increase in normalized ensemble spread. This study demonstrates the small but non‐negligible contribution of lateral terrestrial water flow to precipitation at continental scale.

Assimilation of the Rain Gauge Measurements Using Particle Filter

AbstractThe well‐recognized constraint of nonlinear and non‐Gaussian distribution of rainfall observation limits its assimilation in the high‐dimensional numerical weather prediction (NWP) model. In this study, rainfall observed from Indian Meteorological Department (IMD) rain gauges over Indian landmass is assimilated in the Weather Research and Forecasting (WRF) model using particle filter. In the framework of imperfect weather models, particles (or ensembles) for rainfall predictions are created with various combinations of model physics (viz., cumulus parameterization, microphysics and planetary boundary layer schemes). The multiple hypotheses are used to determine the weights for different particles, and this is the step where IMD rainfall data are used for assimilation. Further, a resampling step is performed to generate new particles from high weight particles using stochastic kinetic‐energy backscatter scheme (SKEBS) method in which dynamical variables are perturbed into the model physics. Results, based on rainfall verification scores, suggest that the assimilation of the rainfall using particle filter could improve the prediction of rainfall over CNT runs (unweighted particles; without assimilation). Moreover, surface and vertical profile of temperature, water vapor mixing ratio (WVMR), and wind speed are also improved in 24‐hr forecasts.

Assessment of the Forecast Skill of Multiphysics and Multistochastic Methods within the GRAPES Regional Ensemble Prediction System in the East Asian Monsoon Region

Abstract To more comprehensively and accurately address model uncertainties in the East Asia monsoon region, a single-physics suite, where each ensemble member uses the same set of physics parameterizations as the control member in combination with multiple stochastic schemes, is developed to investigate if the multistochastic schemes that combine different stochastic schemes together can be an alternative to a multiphysics suite, where each ensemble member uses a different set of physics parameterizations (e.g., cumulus convection, boundary layer, surface layer, microphysics, and shortwave and longwave radiation). For this purpose, two experiments are performed for a summer monsoon month over China: one with a multiphysics suite and the other with a single-physics suite combined with multistochastic schemes. Three stochastic schemes are applied: the stochastically perturbed parameterizations (SPP) scheme, consisting of temporally and spatially varying perturbations of 18 parameters in the microphysics, convection, boundary layer, and surface layer parameterization schemes; the stochastically perturbed parameterization tendencies (SPPT) scheme; and the stochastic kinetic energy backscatter (SKEB) scheme. The combination of the three stochastic schemes is compared with the multiphysics suite in the Global and Regional Assimilation and Prediction Enhanced System–Regional Ensemble Prediction System with a horizontal grid spacing of 15 km. Verification results show that, overall, a single-physics suite that combines SPP, SPPT, and SKEB outperforms the multiphysics suite in precipitation verification and verification for upper-air weather variables, 10-m zonal wind, and 2-m temperature in the East Asian monsoon region. The indication is that a single-physics suite combining SPP, SPPT, and SKEB may be an appropriate alternative to a multiphysics suite. This finding lays a foundation for the development and design of future regional and global ensembles.

Sensitivity of a Bowing Mesoscale Convective System to Horizontal Grid Spacing in a Convection-Allowing Ensemble

The bow echo, a mesoscale convective system (MCS) responsible for much hail and wind damage across the United States, is associated with poor skill in convection-allowing numerical model forecasts. Given the decrease in convection-allowing grid spacings within many operational forecasting systems, we investigate the effect of finer resolution on the character of bowing-MCS development in a real-data numerical simulation. Two ensembles were generated: one with a single domain of 3-km horizontal grid spacing, and another nesting a 1-km domain with two-way feedback. Ensemble members were generated from their control member with a stochastic kinetic-energy backscatter scheme, with identical initial and lateral-boundary conditions. Results suggest that resolution reduces hindcast skill of this MCS, as measured with an adaptation of the object-based Structure–Amplitude–Location method. The nested 1-km ensemble produces a faster system than in both the 3-km ensemble and observations. The nested 1-km simulation also produced stronger cold pools, which could be enhanced by the increased (fractal) cloud surface area with higher resolution, allowing more entrainment of dry air and hence increased evaporative cooling.

Sensitivity Analysis of the Spatial Structure of Forecasts in Mesoscale Models: Noncontinuous Model Parameters

AbstractIn 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.

Representing Model Uncertainty by Multi-Stochastic Physics Approaches in the GRAPES Ensemble

To represent model uncertainties more comprehensively, a stochastically perturbed parameterization (SPP) scheme consisting of temporally and spatially varying perturbations of 18 parameters in the microphysics, convection, boundary layer, and surface layer parameterization schemes, as well as the stochastically perturbed parameterization tendencies (SPPT) scheme, and the stochastic kinetic energy backscatter (SKEB) scheme, is applied in the Global and Regional Assimilation and Prediction Enhanced System—Regional Ensemble Prediction System (GRAPES-REPS) to evaluate and compare the general performance of various combinations of multiple stochastic physics schemes. Six experiments are performed for a summer month (1–30 June 2015) over China and multiple verification metrics are used. The results show that: (1) All stochastic experiments outperform the control (CTL) experiment, and all combinations of stochastic parameterization schemes perform better than the single SPP scheme, indicating that stochastic methods can effectively improve the forecast skill, and combinations of multiple stochastic parameterization schemes can better represent model uncertainties; (2) The combination of all three stochastic physics schemes (SPP, SPPT, and SKEB) outperforms any other combination of two schemes in precipitation forecasting and surface and upper-air verification to better represent the model uncertainties and improve the forecast skill; (3) Combining SKEB with SPP and/or SPPT results in a notable increase in the spread and reduction in outliers for the upper-air wind speed. SKEB directly perturbs the wind field and therefore its addition will greatly impact the upper-air wind-speed fields, and it contributes most to the improvement in spread and outliers for wind; (4) The introduction of SPP has a positive added value, and does not lead to large changes in the evolution of the kinetic energy (KE) spectrum at any wavelength; (5) The introduction of SPPT and SKEB would cause a 5%–10% and 30%–80% change in the KE of mesoscale systems, and all three stochastic schemes (SPP, SPPT, and SKEB) mainly affect the KE of mesoscale systems. This study indicates the potential of combining multiple stochastic physics schemes and lays a foundation for the future development and design of regional and global ensembles.

Diurnal Effects of Regional Soil Moisture Anomalies on the Great Plains Low-Level Jet

AbstractThe Great Plains (GP) low-level jet (GPLLJ) contributes to GP warm season water resources (precipitation), wind resources, and severe weather outbreaks. Past research has shown that synoptic and local mesoscale physical mechanisms (Holton and Blackadar mechanisms) are required to explain GPLLJ variability. Although soil moisture–PBL interactions are central to local mechanistic theories, the diurnal effect of regional soil moisture anomalies on GPLLJ speed, northward penetration, and propensity for severe weather is not well known. In this study, two 31-member WRF-ARW stochastic kinetic energy backscatter scheme ensembles simulate a typical warm season GPLLJ case under CONUS-wide wet and dry soil moisture scenarios. In the GP (24°–48°N, 103°–90°W), ensemble mean differences in sensible heating and PBL height of 25–150 W m−2 and 100–700 m, respectively, at 2100 UTC (afternoon) culminate in GPLLJ 850-hPa wind speed differences of 1–4 m s−1 12 hours later (0900 UTC; early morning). Greater heat accumulation in the daytime PBL over dry soil impacts the east–west geopotential height gradient in the GP (synoptic conditions and Holton mechanism) resulting in a deeper thermal low in the northern GP, causing increases in the geostrophic wind. Enhanced daytime turbulent mixing over dry soil impacts the PBL structure (Blackadar mechanism), leading to increased ageostrophic wind. Overnight geostrophic and ageostrophic winds constructively interact, leading to a faster nocturnal GPLLJ over dry soil. Ensemble differences in CIN (~50–150 J kg−1) and CAPE (~500–1000 J kg−1) have implications for severe weather predictability.

Validation of Snow Multibands in the Comma Head of an Extratropical Cyclone Using a 40-Member Ensemble

Abstract This paper investigates the ability of the Weather Research and Forecasting (WRF) Model in simulating multiple small-scale precipitation bands (multibands) within the extratropical cyclone comma head using four winter storm cases from 2014 to 2017. Using the model output, some physical processes are explored to investigate band prediction. A 40-member WRF ensemble was constructed down to 2-km grid spacing over the Northeast United States using different physics, stochastic physics perturbations, different initial/boundary conditions from the first five perturbed members of the Global Forecast System (GFS) Ensemble Reforecast (GEFSR), and a stochastic kinetic energy backscatter scheme (SKEBS). It was found that 2-km grid spacing is adequate to resolve most snowbands. A feature-based verification is applied to hourly WRF reflectivity fields from each ensemble member and the WSR-88D radar reflectivity at 2-km height above sea level. The Method for Object-Based Diagnostic Evaluation (MODE) tool is used for identifying multibands, which are defined as two or more bands that are 5–20 km in width and that also exhibit a >2:1 aspect ratio. The WRF underpredicts the number of multibands and has a slight eastward position bias. There is no significant difference in frontogenetical forcing, vertical stability, moisture, and vertical shear between the banded versus nonbanded members. Underpredicted band members tend to have slightly stronger frontogenesis than observed, which may be consolidating the bands, but overall there is no clear linkage in ambient condition errors and band errors, thus leaving the source for the band underprediction motivation for future work.

A methodology for sensitivity analysis of spatial features in forecasts: the stochastic kinetic energy backscatter scheme

Stochastic kinetic energy backscatter schemes (SKEBSs) are introduced in numerical weather forecast models to represent uncertainties related to unresolved subgrid‐scale processes. These schemes are formulated using a set of parameters that must be determined using physical knowledge and/or to obtain a desired outcome. Here, a methodology is developed for assessing the effect of four factors on spatial features of forecasts simulated by the SKEBS‐enabled Weather Research and Forecasting model. The four factors include two physically motivated SKEBS parameters (the determining amplitude of perturbations applied to stream function and potential temperature tendencies), a purely stochastic element (a seed used in generating random perturbations) and a factor reflecting daily variability. A simple threshold‐based approach for identifying coherent objects within forecast fields is employed, and the effect of the four factors on object features (e.g. number, size and intensity) is assessed. Four object types are examined: upper‐air jet streaks, low‐level jets, precipitation areas and frontal boundaries. The proposed method consists of a set of standard techniques in experimental design, based on the analysis of variance, tailored to sensitivity analysis. More specifically, a Latin square design is employed to reduce the number of model simulations necessary for performing the sensitivity analysis. Fixed effects and random effects models are employed to assess the main effects and the percentage of the total variability explained by the four factors. It is found that the two SKEBS parameters do not have an appreciable and/or statistically significant effect on any of the examined object features.

Uncertainty Propagation in Coupled Atmosphere–Wave–Ocean Prediction System: A Study of Hurricane Earl (2010)

This study focuses on understanding the evolution of Hurricane Earl (2010) with respect to random perturbations in the storm’s initial strength, size, and asymmetry in wind distribution. We rely on the Unified Wave Interface-Coupled Model (UWIN-CM), a fully coupled atmosphere–wave–ocean system to generate a storm realization ensemble, and use polynomial chaos (PC) expansions to build surrogate models for time evolution of both the maximum wind speed and minimum sea level pressure in Earl. The resulting PC surrogate models provide statistical insights on probability distributions of model responses throughout the simulation time span. Statistical analysis of rapid intensification (RI) suggests that initial perturbations having intensified and counterclockwise-rotated winds are more likely to undergo RI. In addition, for the range of initial conditions considered RI seems mostly sensitive to azimuthally averaged maximum wind speed and asymmetry orientation, rather than storm size and asymmetry magnitude; this is consistent with global sensitivity analysis of PC surrogate models. Finally, we combine initial condition perturbations with a stochastic kinetic energy backscatter scheme (SKEBS) forcing in the UWIN-CM simulations and conclude that the storm tracks are substantially influenced by the SKEBS forcing perturbations, whereas the perturbations in initial conditions alone had only limited impact on the storm-track forecast.

Stochastically Perturbed Parameterizations in an HRRR-Based Ensemble

Abstract A stochastically perturbed parameterization (SPP) approach that spatially and temporally perturbs parameters and variables in the Mellor–Yamada–Nakanishi–Niino planetary boundary layer scheme (PBL) and introduces initialization perturbations to soil moisture in the Rapid Update Cycle land surface model was developed within the High-Resolution Rapid Refresh convection-allowing ensemble. This work is a follow-up study to a work performed using the Rapid Refresh (RAP)-based ensemble. In the present study, the SPP approach was used to target the performance of precipitation and low-level variables (e.g., 2-m temperature and dewpoint, and 10-m wind). The stochastic kinetic energy backscatter scheme and the stochastic perturbation of physics tendencies scheme were combined with the SPP approach and applied to the PBL to target upper-level variable performance (e.g., improved skill and reliability). The three stochastic experiments (SPP applied to PBL only, SPP applied to PBL combined with SKEB and SPPT, and stochastically perturbed soil moisture initial conditions) were compared to a mixed-physics ensemble. The results showed a positive impact from initial condition soil moisture perturbations on precipitation forecasts; however, it resulted in an increase in 2-m dewpoint RMSE. The experiment with perturbed parameters within the PBL showed an improvement in low-level wind forecasts for some verification metrics. The experiment that combined the three stochastic approaches together exhibited improved RMSE and spread for upper-level variables. Our study demonstrated that, by using the SPP approach, forecasts of specific variables can be improved. Also, the results showed that using a single-physics suite ensemble with stochastic methods is potentially an attractive alternative to using multiphysics for convection allowing ensembles.

Impact of Representing Model Error in a Hybrid Ensemble-Variational Data Assimilation System for Track Forecast of Tropical Cyclones over the Bay of Bengal

Uncertainties in the numerical weather prediction models are generally not well-represented in ensemble-based data assimilation (DA) systems. The performance of an ensemble-based DA system becomes suboptimal, if the sources of error are undersampled in the forecast system. The present study examines the effect of accounting for model error treatments in the hybrid ensemble transform Kalman filter—three-dimensional variational (3DVAR) DA system (hybrid) in the track forecast of two tropical cyclones viz. Hudhud and Thane, formed over the Bay of Bengal, using Advanced Research Weather Research and Forecasting (ARW-WRF) model. We investigated the effect of two types of model error treatment schemes and their combination on the hybrid DA system; (i) multiphysics approach, which uses different combination of cumulus, microphysics and planetary boundary layer schemes, (ii) stochastic kinetic energy backscatter (SKEB) scheme, which perturbs the horizontal wind and potential temperature tendencies, (iii) a combination of both multiphysics and SKEB scheme. Substantial improvements are noticed in the track positions of both the cyclones, when flow-dependent ensemble covariance is used in 3DVAR framework. Explicit model error representation is found to be beneficial in treating the underdispersive ensembles. Among the model error schemes used in this study, a combination of multiphysics and SKEB schemes has outperformed the other two schemes with improved track forecast for both the tropical cyclones.

The impact of stochastic parametrisations on the representation of the Asian summer monsoon

The impact of the stochastic schemes Stochastically Perturbed Parametrisation Tendencies (SPPT) and Stochastic Kinetic Energy Backscatter Scheme (SKEBS) on the representation of interannual variability in the Asian summer monsoon is examined in the coupled climate model CCSM4. The Webster–Yang index, measuring anomalies of a specified wind-shear index in the monsoon region, is used as a metric for monsoon strength, and is used to analyse the output of three model integrations: one deterministic, one with SPPT, and one with SKEBS. Both schemes show improved variability, which we trace back to improvements in the El Niño-Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD). SPPT improves the representation of ENSO and through teleconnections thereby the monsoon, supporting previous work on the benefits of this scheme on the model climate. SKEBS also improves monsoon variability by way of improving the representation of the IOD, in particular by breaking an overly strong coupling to ENSO.

The impact of stochastic physics on tropical rainfall variability in global climate models on daily to weekly time scales

AbstractMany global atmospheric models have too little precipitation variability in the tropics on daily to weekly time scales and also a poor representation of tropical precipitation extremes associated with intense convection. Stochastic parameterizations have the potential to mitigate this problem by representing unpredictable subgrid variability that is left out of deterministic models. We evaluate the impact on the statistics of tropical rainfall of two stochastic schemes: the stochastically perturbed parameterization tendency scheme (SPPT) and stochastic kinetic energy backscatter scheme (SKEBS), in three climate models: EC‐Earth, the Met Office Unified Model, and the Community Atmosphere Model, version 4. The schemes generally improve the statistics of simulated tropical rainfall variability, particularly by increasing the frequency of heavy rainfall events, reducing its persistence and increasing the high‐frequency component of its variability. There is a large range in the size of the impact between models, with EC‐Earth showing the largest improvements. The improvements are greater than those obtained by increasing horizontal resolution to ∼20 km. Stochastic physics also strongly affects projections of future changes in the frequency of extreme tropical rainfall in EC‐Earth. This indicates that small‐scale variability that is unresolved and unpredictable in these models has an important role in determining tropical climate variability statistics. Using these schemes, and improved schemes currently under development, is therefore likely to be important for producing good simulations of tropical variability and extremes in the present day and future.