skip to main content


Title: Evaluating Ensemble Kalman Filter Analyses of Severe Hailstorms on 8 May 2017 in Colorado: Effects of State Variable Updating and Multimoment Microphysics Schemes on State Variable Cross Covariances

Ensemble Kalman filter (EnKF) analyses of the storms associated with the 8 May 2017 Colorado severe hail event using either the Milbrandt and Yau (MY) or the NSSL double-moment bulk microphysics scheme in the forecast model are evaluated. With each scheme, two experiments are conducted in which the reflectivity ( Z) observations update in addition to dynamic and thermodynamic variables: 1) only the hydrometeor mixing ratios or 2) all microphysical variables. With fewer microphysical variables directly constrained by the Z observations, only updating hydrometeor mixing ratios causes the forecast error covariance structure to become unreliable, and results in larger errors in the analysis. Experiments that update all microphysical variables produce analyses with the lowest Z root-mean-square innovations; however, comparing the estimated hail size against hydrometeor classification algorithm output suggests that further constraint from observations is needed to more accurately estimate surface hail size. Ensemble correlation analyses are performed to determine the impact of hail growth assumptions in the MY and NSSL schemes on the forecast error covariance between microphysical and thermodynamic variables. In the MY scheme, Z is negatively correlated with updraft intensity because the strong updrafts produce abundant small hail aloft. The NSSL scheme predicts the growth of large hail aloft; consequently, Z is positively correlated with storm updraft intensity and hail state variables. Hail production processes are also shown to alter the background error covariance for liquid and frozen hydrometeor species. Results in this study suggest that EnKF analyses are sensitive to the choice of MP scheme (e.g., the treatment of hail growth processes).

 
more » « less
NSF-PAR ID:
10155284
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Monthly Weather Review
Volume:
148
Issue:
6
ISSN:
0027-0644
Page Range / eLocation ID:
p. 2365-2389
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Hail forecast evaluations provide important insight into microphysical treatment of rimed ice. In this study we evaluate explicit 0–90-min EnKF-based storm-scale (500-m horizontal grid spacing) hail forecasts for a severe weather event that occurred in Oklahoma on 19 May 2013. Forecast ensembles are run using three different bulk microphysics (MP) schemes: the Milbrandt–Yau double-moment scheme (MY2), the Milbrandt–Yau triple-moment scheme (MY3), and the NSSL variable density-rimed ice double-moment scheme (NSSL). Output from a hydrometeor classification algorithm is used to verify surface hail size forecasts. All three schemes produce forecasts that predict the coverage of severe surface hail with moderate to high skill, but exhibit less skill at predicting significant severe hail coverage. A microphysical budget analysis is conducted to better understand hail growth processes in all three schemes. The NSSL scheme uses two-variable density-rimed ice categories to create large hailstones from dense, wet growth graupel particles; however, it is noted the scheme underestimates the coverage of significant severe hail. Both the MY2 and MY3 schemes produce many small hailstones aloft from unrimed, frozen raindrops; in the melting layer, hailstones become much larger than observations because of the excessive accretion of water. The results of this work highlight the importance of using a MP scheme that realistically models microphysical processes.

     
    more » « less
  2. Abstract

    Day-ahead (20–22 h) 3-km grid spacing convection-allowing model forecasts are performed for a severe hail event that occurred in Denver, Colorado, on 8 May 2017 using six different multimoment microphysics (MP) schemes including: the Milbrandt–Yau double-moment (MY2), Thompson (THO), NSSL double-moment (NSSL), Morrison double-moment graupel (MOR-G) and hail (MOR-H), and Predicted Particle Properties (P3) schemes. Hail size forecasts diagnosed using the Thompson hail algorithm and storm surrogates predict hail coverage. For this case hail forecasts predict the coverage of hail with a high level of skill but underpredict hail size. The storm surrogate updraft helicity predicts the coverage of severe hail with the most skill for this case. Model data are analyzed to assess the effects of microphysical treatments related to rimed ice. THO uses diagnostic equations to increase the size of graupel within the hail core. MOR-G and MOR-H predict small rimed ice aloft; excessive size sorting and increased fall speeds cause MOR-H to predict more and larger surface hail than MOR-G. The MY2 and NSSL schemes predict large, dense rimed ice particles because both schemes predict separate hail and graupel categories. The NSSL scheme predicts relatively little hail for this case; however, the hail size forecast qualitatively improves when the maximum size of both hail and graupel is considered. The single ice category P3 scheme only predicts dense hail near the surface while above the melting layer large concentrations of low-density ice dominate.

     
    more » « less
  3. Abstract

    This study uses a new, unique dataset created by combining multi-Doppler radar wind and reflectivity analysis, diabatic Lagrangian analysis (DLA) retrievals of temperature and water substance, and a complex hail trajectory model to create millions of numerically simulated hail trajectories in the Kingfisher, Oklahoma, supercell on 29 May 2012. The DLA output variables are used to obtain a realistic, 4D depiction of the storm’s thermal and hydrometeor structure as required input to the detailed hail growth trajectory model. Hail embryos are initialized in the hail growth module every 3 min of the radar analysis period (2251–0000 UTC) to produce over 2.7 million hail trajectories. A spatial integration technique considering all trajectories is used to identify locations within the supercell where melted particles and subsevere and severe hailstones reside in their lowest and highest concentrations. It is found that hailstones are more likely to reside for longer periods closer to the downshear updraft within the midlevel mesocyclone in a region of decelerated midlevel mesocyclonic horizontal flow, termed the downshear deceleration zone (DDZ). Additionally, clusters of trajectories are analyzed using a trajectory clustering method. Trajectory clusters show there are many trajectory pathways that result in hailstones ≥ 4.5 cm, including trajectories that begin upshear of the updraft away from ideal growth conditions and trajectories that grow within the DDZ. There are also trajectory clusters with similar shapes that experience widely different environmental and hailstone characteristics along the trajectory.

    Significance Statement

    The purpose of this study is to understand how hail grew in a thunderstorm that was observed by numerous instruments. The observations were input into a hail trajectory model to simulate hail growth. We found a part of the storm near the updraft where hailstones could remain aloft longer and therefore grow larger. Most modeled severe hailstones were found in the storm in this region. However, we also found that there are many different pathways hailstones can take to become large, although there are still some common characteristics among the pathways.

     
    more » « less
  4. Abstract

    Assimilating radar reflectivity into convective-scale NWP models remains a challenging topic in radar data assimilation. A primary reason is that the reflectivity forward observation operator is highly nonlinear. To address this challenge, a power transformation function is applied to the WRF Model’s hydrometeor and water vapor mixing ratio variables in this study. Three 3D variational data assimilation experiments are performed and compared for five high-impact weather events that occurred in 2019: (i) a control experiment that assimilates reflectivity using the original hydrometeor mixing ratios as control variables, (ii) an experiment that assimilates reflectivity using power-transformed hydrometeor mixing ratios as control variables, and (iii) an experiment that assimilates reflectivity and retrieved pseudo–water vapor observations using power-transformed hydrometeor and water vapor mixing ratios (qυ) as control variables. Both qualitative and quantitative evaluations are performed for 0–3-h forecasts from the five cases. The analysis and forecast performance in the two experiments with power-transformed mixing ratios is better than the control experiment. Notably, the assimilation of pseudo–water vapor with power-transformedqυas an additional control variable is found to improve the performance of the analysis and short-term forecasts for all cases. In addition, the convergence rate of the cost function minimization for the two experiments that use the power transformation is faster than that of the control experiments.

    Significance Statement

    The effective use of radar reflectivity observations in any data assimilation scheme remains an important research topic because reflectivity observations explicitly include information about hydrometeors and also implicitly include information about the distribution of moisture within storms. However, it is difficult to assimilate reflectivity because the reflectivity forward observation operator is highly nonlinear. This study seeks to identify a more effective way to assimilate reflectivity into a convective-scale NWP model to improve the accuracy of predictions of high-impact weather events.

     
    more » « less
  5. Abstract

    For data assimilation to provide faithful state estimates for dynamical models, specifications of observation uncertainty need to be as accurate as possible. Innovation-based methods based on Desroziers diagnostics, are commonly used to estimate observation uncertainty, but such methods can depend greatly on the prescribed background uncertainty. For ensemble data assimilation, this uncertainty comes from statistics calculated from ensemble forecasts, which require inflation and localization to address under sampling. In this work, we use an ensemble Kalman filter (EnKF) with a low-dimensional Lorenz model to investigate the interplay between the Desroziers method and inflation. Two inflation techniques are used for this purpose: 1) a rigorously tuned fixed multiplicative scheme and 2) an adaptive state-space scheme. We document how inaccuracies in observation uncertainty affect errors in EnKF posteriors and study the combined impacts of misspecified initial observation uncertainty, sampling error, and model error on Desroziers estimates. We find that whether observation uncertainty is over- or underestimated greatly affects the stability of data assimilation and the accuracy of Desroziers estimates and that preference should be given to initial overestimates. Inline estimates of Desroziers tend to remove the dependence between ensemble spread–skill and the initially prescribed observation error. In addition, we find that the inclusion of model error introduces spurious correlations in observation uncertainty estimates. Further, we note that the adaptive inflation scheme is less robust than fixed inflation at mitigating multiple sources of error. Last, sampling error strongly exacerbates existing sources of error and greatly degrades EnKF estimates, which translates into biased Desroziers estimates of observation error covariance.

    Significance Statement

    To generate accurate predictions of various components of the Earth system, numerical models require an accurate specification of state variables at our current time. This step adopts a probabilistic consideration of our current state estimate versus information provided from environmental measurements of the true state. Various strategies exist for estimating uncertainty in observations within this framework, but are sensitive to a host of assumptions, which are investigated in this study.

     
    more » « less