Assimilation Of Reflectivity Data Research Articles

Development of Multiscale EnKF within GSI and Its Applications to Multiple Convective Storm Cases with Radar Reflectivity Data Assimilation Using the FV3 Limited-Area Model

Abstract To improve the representation of all relevant scales in initial conditions for large-domain convection-allowing models, a new multiscale ensemble Kalman filter (MEnKF) algorithm is developed and implemented within the Gridpoint Statistical Interpolation analysis system (GSI) data assimilation framework coupled with the Finite-Volume Cubed-Sphere Dynamical Core (FV3) limited-area model. The algorithm utilizes ensemble background error covariances filtered to match the observations assimilated. This is realized in a sequential manner. 1) When assimilating coarse-resolution observations such as radiosondes, ensemble background perturbations are filtered to remove scales smaller than those the observations can represent, along with relatively large horizontal localization radii to ensure low-noise and balanced analysis increments. 2) The resulting ensemble analyses from the first step then serve as the background to assimilate denser observations such as radar data with smaller localization radii. Several passes can be taken to assimilate all observations. In this paper, vertically increasing horizontal filter scales are used when assimilating rawinsonde and surface observations together, while radar data are assimilated in the second step. The algorithm is evaluated through six convective storm cases during May 2021, with cycled assimilation of either conventional data only or with additional radar reflectivity followed by 24-h ensemble forecasts. Overall, positive impacts of the MEnKF on forecasts are obtained regardless of reflectivity data; its advantage over the single-scale EnKF is most significant in surface humidity and temperature forecasts up to at least 12 h. More accurate hourly precipitation forecasts with MEnKF can last up to 24 h for light rain. Furthermore, MEnKF forecasts higher ensemble probabilities for the observed hazardous events.

Impact of Direct Radar Reflectivity Data Assimilation on the Simulation of Mesoscale Descending Inflow and Secondary Eyewall Formation in Hurricane Matthew (2016)

AbstractThe impact of assimilating ground‐based radar reflectivity on the rainband structure and secondary eyewall formation (SEF) of Hurricane Matthew (2016) is investigated within the framework of the Hurricane Weather Research and Forecasting model and its hybrid three‐dimensional ensemble‐variational data assimilation (DA) system. Compared to the control experiment (no radar reflectivity DA), the radar reflectivity DA experiment shows a clear signal of concentric eyewall and eyewall replacement cycle. Results demonstrate that radar reflectivity DA improves the stratiform rainband analysis, resulting in the mid‐level cooling associated with mesoscale descending inflow (MDI). The MDI further contributes to the low‐level acceleration maximum with boundary layer dynamics and triggers new convective updrafts in the SEF region. Momentum budget analysis also suggests that the mean vertical advection of absolute angular momentum plays an important role in the local momentum tendency in the SEF region in Hurricane Matthew (2016).

A Comparison of 3DEnVar and 4DEnVar for Convective-Scale Direct Radar Reflectivity Data Assimilation in the Context of a Filter and a Smoother

Abstract The Gridpoint Statistical Interpolation (GSI)-based four- and three-dimensional ensemble–variational (4DEnVar and 3DEnVar) methods are compared as a smoother and a filter, respectively, for rapidly changing storms using the convective-scale direct radar reflectivity data assimilation (DA) framework. Two sets of experiments with varying DA window lengths (WLs; 20, 40, 100, and 160 min) and radar observation intervals (RIs; 20 and 5 min) are conducted for the 5–6 May 2019 case. The RI determines the temporal resolution of ensemble perturbations for the smoother and the DA interval for the filter spanning the WL. For experiments with a 20-min RI, evaluations suggest that the filter and the smoother have comparable performance with a 20-min WL; however, extending the WL results in the outperformance of the filter over the smoother. Diagnostics reveal that the degradation of the smoother is attributed to the increased degree of nonlinearity and the issue of time-independent localization as the WL extends. Evaluations for experiments with different RIs under the same WL indicate that the outperformance of the filter over the smoother diminishes for most forecast hours at thresholds of 30 dBZ and above when shortening the RI. Diagnostics show that more frequent interruptions of the model introduce model imbalance for the filter, and the increased temporal resolution of ensemble perturbations enhances the degree of nonlinearity for the smoother. The impact of model imbalance on the filter overwhelms the enhanced nonlinearity on the smoother as the RI reduces. Significance Statement The background uncertainties of rapidly changing storms suffer from fast error growth and high degrees of nonlinearities during the data assimilation (DA) period. Two variants of the ensemble-based DA method can account for such temporal evolution. The smoother uses background ensemble from multiple observation times over an assimilation period to estimate the propagation of statistics. The filter frequently calculates the statistics at multiple observation times over the same period. Current comparisons of the smoother and the filter were mostly performed using simple models; however, unknowns remain for convection-allowing forecasts with additional complexities. This study compares the filter and the smoother for the convective-scale analysis and prediction using a real-data study and finds that the comparison varies with the assimilation period and the observation interval.

Impacts of adding hydrometeor control variables on the radar reflectivity data assimilation for the 6–8 August 2018 mesoscale convective system case

Hydrometeor information from radar reflectivity observations plays an important role in initializing convective-scale numerical weather forecast (NWP). The introduction of hydrometeor control variables facilitates rational division of the reflectivity information to obtain more accurate initial conditions of cloud and precipitation. Using the Weather Research and Forecasting (WRF) model, the effects of adding hydrometeor control variables on assimilation and forecasting are explored for a strong convective process occurring in Liaoning Province in this study. The results show that, compared with the experiment of no hydrometeor control variable, the presences of hydrometeor control variables are beneficial to the analyses of horizontal and vertical distribution of hydrometeors, which enables a better match of the analyzed structures with the reflectivity observations. Moreover, the reasonable hydrometeor distribution is able to reduce the generation of false echoes during assimilation. According to the results of the forecast experiments, using hydrometeor control variables tends to slow down the dissipation of hydrometeors, leading improved forecasting skills of the quantitative precipitation.

Impact of Radar Reflectivity Data Assimilation Frequency on Convection-Allowing Forecasts of Diverse Cases over the Continental United States

Abstract The sensitivity of convection-allowing forecasts over the continental United States to radar reflectivity data assimilation (DA) frequency is explored within the Gridpoint Statistical Interpolation (GSI)-based ensemble–variational (EnVar) system. Experiments with reflectivity DA intervals of 60, 20, and 5 min (RAIN60, RAIN20, and RAIN5, respectively) are conducted using 10 diverse cases. Quantitative verification indicates that the degree of sensitivity depends on storm features during the radar DA period. Five developing storms show high sensitivity, whereas five mature or decaying storms do not. The 20-min interval is the most reliable given its best overall performance compared to the 5- and 60-min intervals. Diagnostics suggest that the differences in analyzed cold pools (ACPs) among RAIN60, RAIN20, and RAIN5 vary by storm features during the radar DA period. Such ACP differences result in different forecast skills. In the case where RAIN20 outperforms RAIN60 and the case where RAIN5 outperforms RAIN20, assimilation of reflectivity with a higher frequency commonly produces enhanced and widespread ACPs, promoting broader storms that match better with reality than a lower frequency. In the case where RAIN5 performs worse than RAIN20, the model imbalance of RAIN5 overwhelms information gain associated with frequent assimilation, producing overestimated and spuriously fast-moving ACPs. In the cases where little sensitivity to the reflectivity DA frequency is found, similar ACPs are produced.

Comparison of the WRF-FDDA-Based Radar Reflectivity and Lightning Data Assimilation for Short-Term Precipitation and Lightning Forecasts of Severe Convection

This work evaluates and compares the performance of the radar reflectivity and lightning data assimilation schemes implemented in weather research and forecasting-four-dimensional data assimilation (WRF-FDDA) for short-term precipitation and lightning forecasts. All six mesoscale convective systems (MCSs) with a duration greater than seven hours that occurred in the Guangdong Province of China during June 2020 were included in the experiments. The results show that both the radar reflectivity data assimilation and lightning data assimilation improved the analyses and short-term forecasts of the precipitation and lightning of the MCSs. On average, for precipitation forecasts, the experiments with radar reflectivity data assimilation performed better than those with lightning data assimilation; however, for lightning forecasts, the experiments with lightning data assimilation performed better in the analysis period and 1 h forecast, and for some cases, the superiority lasted to three forecast hours. This highlights the potential of lightning data assimilation in short-term lightning forecasting.

Direct Assimilation of Radar Reflectivity Data Using Ensemble Kalman Filter Based on a Two-Moment Microphysics Scheme for the Analysis and Forecast of Typhoon Lekima (2019)

We investigate the impact of directly assimilating radar reflectivity data using an ensemble Kalman filter (EnKF) based on a double-moment (DM) microphysics parameterization (MP) scheme in GSI-EnKF data assimilation (DA) framework and WRF model for a landfall typhoon Lekima (2019). Observations from a single operational coastal Doppler are quality-controlled and assimilated. Compared with the baseline experiment initialized by GFS analysis, the reflectivity data assimilation experiment (Z-DA) resulted in an obvious improvement in both structural analysis and typhoon forecast skills in terms of intensity, precipitation, and track. Sensitivity experiments were conducted to evaluate the ability of EnKF to update certain state variables considering that the degree of freedom of analytical variables increased with a DM MP scheme. When either the total number concentration or other large-scale state variables that are not directly linked to reflectivity observations via the observation operator are not updated, the tendency of RMSIs and PS to be imbalanced is significantly increased during DA cycles compared to those of Z-DA with updating a full set of state variables, resulting in increased intensity and track forecast errors. These results indicate that the reliable ensemble covariance could handle the underconstraint issue associated with the DM scheme, and helps in obtaining more physically balanced analytical fields.

Assimilation of Radar Reflectivity Using a Time‐Lagged Ensemble Based Ensemble Kalman Filter With the “Cloud‐Dependent” Background Error Covariances

AbstractThis study proposes a method for direct assimilation of radar reflectivity in cycled deterministic forecast systems using a time‐lagged ensemble based ensemble Kalman filter with the “cloud‐dependent” background error covariances. With this method, forecasts from multiple deterministic forecast cycles that are initialized at different times but valid at the same time, are collected as a time‐lagged ensemble to compute the background error covariances. An algorithm to compute the specified “cloud‐dependent” background error covariances was devised, in which the vertical background error covariance is computed for each of the “cloud‐feature” bins which are divided according to the column maximum radar reflectivity and cloud top height. The single‐observation experiments and the observing system simulation experiments (OSSEs) indicate that the radar reflectivity data assimilation method was able to capture the main thermodynamic and microphysical features of convective clouds and was effective in reducing the spin‐up time of convection, and generating convective cells which were missed in the background. Meanwhile, assimilating “zero‐reflectivity” observations helped suppress spurious convection. The method provided overall more accurate thermodynamic analysis, radar reflectivity and precipitation analysis and forecast of the MCS, while did not significantly increase computational cost, compared to a cloud‐analysis based latent‐heat nudging approach.

Rapid Update with EnVar Direct Radar Reflectivity Data Assimilation for the NOAA Regional Convection-Allowing NMMB Model over the CONUS: System Description and Initial Experiment Results

This study first describes the extended Grid-Point Statistical Interpolation analysis system (GSI)-based ensemble-variational data assimilation (DA) system within the North American Mesoscale Rapid Refresh (NAMRR) system for the Nonhydrostatic Multiscale Model on the B grid (NMMB). Experiments were conducted to examine three critical aspects of data assimilation configuration in this system. Ten retrospective high-impact convective cases during the warm season of 2015–2016 were adopted for testing. A 10-member, 18 h ensemble forecast was launched for each experiment. Specifically, the experiment using horizontal (vertical) localization radii (Lr) of 300 km (0.55-scaled height measured in the nature log of pressure) overall had more skills than that of 500 km (1.1-scaled height) for conventional in-situ observation assimilation. Diagnostics suggest that the higher forecast skills could be attributed to applying smaller Lr in the boundary with large temperature and moisture gradients. For radar DA, the experiment was more skillful with horizontal (vertical) Lr of 15 km (1.1-scaled height) than that of 12 km (0.55-scaled height). Diagnostics suggest that the improved forecasts were achieved by using wider Lr to spread radar observations into unobserved areas more effectively. Slight forecast skill differences between the relaxation inflation factors of 95% and 65% are presented. The impact of varying inflation magnitudes primarily occurred in the upper-level spread.

Development of Convective-Scale Static Background Error Covariance within GSI-Based Hybrid EnVar System for Direct Radar Reflectivity Data Assimilation

AbstractA convective-scale static background-error covariance (BEC) matrix is further developed to include the capability of direct reflectivity assimilation and evaluated within the GSI-based three-dimensional variational (3DVar) and hybrid ensemble–variational (EnVar) methods. Specific developments are summarized as follows: 1) Control variables (CVs) are extended to include reflectivity, vertical velocity, and all hydrometeor types. Various horizontal momentum and moisture CV options are included. 2) Cross correlations between all CVs are established. 3) A storm intensity-dependent binning method is adopted to separately calculate static error matrices for clear-air and storms with varying intensities. The resultant static BEC matrices are simultaneously applied at proper locations guided by the observed reflectivity. 4) The EnVar is extended to adaptively incorporate static BECs based on the quality of ensemble covariances. Evaluation and examination of the new static BECs are first performed on the 8 May 2003 Oklahoma City supercell. Detailed diagnostics and 3DVar examinations suggest zonal/meridional winds and pseudo–relative humidity are selected as horizontal momentum and moisture CVs for direct reflectivity assimilation, respectively; inclusion of cross correlations favors spin up and maintains the analyzed storms; application of binning improves characteristics and persistence of the simulated storm. Relative to an experiment using the full ensemble BECs (Exp-PureEnVar), incorporating static BECs in hybrid EnVar reduces spinup time and better analyzes reflectivity distributions while the background ensemble is deficient in sampling errors. Compared to both pure 3DVar and Exp-PureEnVar, hybrid EnVar better predicts reflectivity distributions and better maintains a strong mesocyclone. Further examination through the 20 May 2013 Oklahoma supercells confirms these results and additionally demonstrates the effectiveness of adaptive hybridization.

Impact of Radar Reflectivity and Lightning Data Assimilation on the Rainfall Forecast and Predictability of a Summer Convective Thunderstorm in Southern Italy

Heavy and localized summer events are very hard to predict and, at the same time, potentially dangerous for people and properties. This paper focuses on an event occurred on 15 July 2020 in Palermo, the largest city of Sicily, causing about 120 mm of rainfall in 3 h. The aim is to investigate the event predictability and a potential way to improve the precipitation forecast. To reach this aim, lightning (LDA) and radar reflectivity data assimilation (RDA) was applied. LDA was able to trigger deep convection over Palermo, with high precision, whereas the RDA had a key role in the prediction of the amount of rainfall. The simultaneous assimilation of both data sources gave the best results. An alert for a moderate–intense forecast could have been issued one hour and a half before the storm developed over the city, even if predicting only half of the total rainfall. A satisfactory prediction of the amount of rainfall could have been issued at 14:30 UTC, when precipitation was already affecting the city. Although the study is centered on a single event, it highlights the need for rapidly updated forecast cycles with data assimilation at the local scale, for a better prediction of similar events.

Radar reflectivity data assimilation method based on background‐dependent hydrometeor retrieval: Comparison with direct assimilation for real cases

AbstractAssimilating radar reflectivity into NWP models is one of the keys to improve the accuracy of convective‐scale numerical weather prediction (NWP). There are generally two major branches of research in radar reflectivity assimilation: directly assimilating radar reflectivity, and indirectly assimilating reflectivity, i.e. assimilating hydrometeors retrieved from radar reflectivity. In this study, the indirect assimilation method based on background‐dependent hydrometeor retrievals is compared with the direct assimilation method using frequent data assimilation cycles for five real data cases. The retrieved hydrometeors in the indirect assimilation method are first verified against the hydrometeor types obtained from a polarimetric hydrometeor classification method. It is illustrated that the background‐dependent hydrometeor retrieval method can obtain reasonable model‐equivalent hydrometeors from radar reflectivity. The analysis increments for hydrometeors with both radar data assimilation methods show similar patterns but their magnitudes are different. Both quantitative and qualitative evaluations of forecasted composite reflectivities and accumulated precipitation indicate that the indirect assimilation method predicts the location and intensity of the simulated convection more accurately than the direct method. Furthermore, the indirect assimilation method is more efficient, which is valuable in real‐time applications by helping deliver forecasts quickly and thus helping forecasters make more timely warning decisions.

Improvement of mesoscale numerical weather prediction WRF-ARW in the Republic of Belarus by assimilation of radar wind and reflectivity data

The forecasts of severe weather events obtained with the WRF numerical mesoscale model with the adapted system for assimilation of reflectivity and radial velocity data from the network of Belarusian Doppler weather radars used in Belhydromet in 2019 are analysed. A description of the system for the echo quality control based on the radar dual-polarisation characteristics and the method for three-dimensional variational assimilation (3D-VAR) used to assimilate data in the WRF model are described. The results of case studies on the simulation of precipitation and strong wind for various circulation types in Belarus with and without radar data assimilation are given. The statistical and object-oriented verification of these forecasts is provided. The results of the comprehensive assessment reveal a decrease in the forecast error for 10-m wind speed for the early forecast hours (+6 h) by 1.34 m/s, as well as a more accurate forecast of the location, orientation of the cloud systems and precipitation zones, and a decrease in the number of false alarms in the version with assimilation. A preliminary conclusion on the possibility of using the forecast results in nowcasting systems is also made.

A radar reflectivity data assimilation method based on background-dependent hydrometeor retrieval: An observing system simulation experiment

Radar reflectivity contains information about hydrometeors and plays an important role in the initialization of convective-scale numerical weather prediction (NWP). In this study, a new background-dependent hydrometeor retrieval method is proposed and retrieved hydrometeors are assimilated into the Weather Research and Forecasting model (WRF), with the aim of improving short-term severe weather forecasts. Compared to traditional approaches that are mostly empirical and static, the retrieval parameters for hydrometeor identification and reflectivity partitioning in the new scheme are extracted in real-time based on the background hydrometeor fields and observed radar reflectivity. It was found that the contributions of hydrometeors to reflectivity change a lot in different reflectivity ranges and heights, indicating that adaptive parameters are necessary for reflectivity partitioning and hydrometeor retrieval. The accuracy of the background-dependent hydrometeor retrieval method and its impact on the subsequent assimilation and forecast were examined through observing system simulation experiments (OSSEs). Results show that by incorporating the background information, the retrieval accuracy was greatly improved, especially in mixed-hydrometeor regions. The assimilation of retrieved hydrometeors helped improve both the hydrometeor analyses and forecasts. With an hourly update cycling configuration, more accurate hydrometeor information was properly transferred to other model variables, such as temperature and humidity fields through the model integration, leading to an improvement of the short-term (0−3 h) precipitation forecasts.

Reflectivity and velocity radar data assimilation for two flash flood events in central Italy: A comparison between 3D and 4D variational methods

AbstractThe aim of this study is to provide an evaluation of the impact of two largely used data assimilation techniques, namely three‐ and four‐dimensional variational data assimilation systems (3D‐Var and 4D‐Var), on the forecasting of heavy precipitation events using the Weather Research and Forecasting (WRF) model. For this purpose, two flash flood events in central Italy are analysed. The first occurred on September 14, 2012 during an Intensive Observation Period of the Hydrological cycle in the Mediterranean experiment (HyMeX) campaign, while the other occurred on May 3, 2018. Radial velocity and reflectivity acquired by C‐band weather radars at Mt. Midia (central Italy) and San Pietro Capofiume (northern Italy), as well as conventional observations (SYNOP and TEMP), are assimilated into the WRF model to simulate these damaging flash flood events. In order to evaluate the impact of the 3D‐Var and 4D‐Var assimilation systems on the estimation of short‐term quantitative precipitation forecasts, several experiments are carried out using conventional observations with and without radar data. Rainfall evaluation is performed by means of point‐by‐point and filtering methodologies. The results point to a positive impact of the 4D‐Var technique compared to results without assimilation and with 3D‐Var experiments. More specifically, the 4D‐Var system produces an increase of up to 22% in terms of the Fractions Skill Score compared to 3D‐Var for the first flash flood event, while an increase of about 5% is achieved for the second event. The use of a warm start initialization results in a considerable reduction in the spin‐up time and a significant improvement in the rainfall forecast, suggesting that the initial precipitation spin‐up problem still occurs when using 4D‐Var.

Impact of airborne cloud radar reflectivity data assimilation on kilometre-scale numerical weather prediction analyses and forecasts of heavy precipitation events

Abstract. This article investigates the potential of W-band radar reflectivity to improve the quality of analyses and forecasts of heavy precipitation events in the Mediterranean area. The “1D+3DVar” assimilation method, operationally employed to assimilate ground-based precipitation radar data in the Météo-France kilometre-scale numerical weather prediction (NWP) model AROME, has been adapted to assimilate the W-band reflectivity measured by the airborne cloud radar RASTA (Radar Airborne System Tool for Atmosphere) during a 2-month period over the Mediterranean area. After applying a bias correction, vertical profiles of relative humidity are first derived via a 1-D Bayesian retrieval, and then used as relative humidity pseudo-observations in the 3DVar assimilation system of AROME. The efficiency of the 1-D Bayesian method in retrieving humidity fields is assessed using independent in-flight humidity measurements. To complement this study, the benefit brought by consistent thermodynamic and dynamic cloud conditions has been investigated by separately and jointly assimilating the W-band reflectivity and horizontal wind measurements collected by RASTA in the 3 h 3DVar assimilation system of AROME. The data assimilation experiments are conducted for a single heavy precipitation event and then also for 32 cases. Results indicate that the W-band reflectivity has a larger impact on the humidity, temperature and pressure fields in the analyses compared to the assimilation of RASTA wind data alone. Besides, the analyses get closer to independent humidity observations if the W-band reflectivity is assimilated alone or jointly with RASTA wind data. Nonetheless, the impact of the W-band reflectivity decreases more rapidly as the forecast range increases when compared to the assimilation of RASTA wind data alone. Generally, the joint assimilation of the W-band reflectivity with wind data results in the best improvement in the rainfall precipitation forecasts. Consequently, results of this study indicate that consistent thermodynamic and dynamic cloud conditions in the analysis leads to an improvement of both model initial conditions and forecasts. Even though to a lesser extent, the assimilation of the W-band reflectivity alone also results in a slight improvement of the rainfall precipitation forecasts.

Mesoscale Resolution Radar Data Assimilation Experiments with the Harmonie Model

This study presents a pre-processing approach adopted for the radar reflectivity data assimilation and results of simulations with the Harmonie numerical weather prediction model. The proposed method creates a 3D regular grid in which a horizontal size of meshes coincides with the horizontal model resolution. This minimizes the representative error associated with the discrepancy between resolutions of informational sources. After such preprocessing, horizontal structure functions and their gradients for radar reflectivity maintain the sizes and shapes of precipitation patterns similar to those of the original data. The method shows an improvement of precipitation prediction within the radar location area in both the rain rates and spatial pattern presentation. It redistributes precipitable water with smoothed values over the common domain since the control runs show, among several sub-domains with increased and decreased values, correspondingly. It also reproduces the mesoscale belts and cell patterns of sizes from a few to ten kilometers in precipitation fields. With the assimilation of radar data, the model simulates larger water content in the middle troposphere within the layer from 1 km to 6 km with major variations at 2.5 km to 3 km. It also reproduces the mesoscale belt and cell patterns of precipitation fields.

Statistics of the Z–R Relationship for Strong Convective Weather over the Yangtze–Huaihe River Basin and Its Application to Radar Reflectivity Data Assimilation for a Heavy Rain Event

The relationship between the radar reflectivity factor (Z) and the rainfall rate (R) is recalculated based on radar observations from 10 Doppler radars and hourly rainfall measurements at 6529 automatic weather stations over the Yangtze–Huaihe River basin. The data were collected by the National 973 Project from June to July 2013 for severe convective weather events. The Z–R relationship is combined with an empirical qr–R relationship to obtain a new Z–qr relationship, which is then used to correct the observational operator for radar reflectivity in the three-dimensional variational (3DVar) data assimilation system of the Weather Research and Forecasting (WRF) model to improve the analysis and prediction of severe convective weather over the Yangtze–Huaihe River basin. The performance of the corrected reflectivity operator used in the WRF 3DVar data assimilation system is tested with a heavy rain event that occurred over Jiangsu and Anhui provinces and the surrounding regions on 23 June 2013. It is noted that the observations for this event are not included in the calculation of the Z–R relationship. Three experiments are conducted with the WRF model and its 3DVar system, including a control run without the assimilation of reflectivity data and two assimilation experiments with the original and corrected reflectivity operators. The experimental results show that the assimilation of radar reflectivity data has a positive impact on the rainfall forecast within a few hours with either the original or corrected reflectivity operators, but the corrected reflectivity operator achieves a better performance on the rainfall forecast than the original operator. The corrected reflectivity operator extends the effective time of radar data assimilation for the prediction of strong reflectivity. The physical variables analyzed with the corrected reflectivity operator present more reasonable mesoscale structures than those obtained with the original reflectivity operator. This suggests that the new statistical Z–R relationship is more suitable for predicting severe convective weather over the Yangtze–Huaihe River basin than the Z–R relationships currently in use.

Assimilation of ZDR Columns for Improving the Spinup and Forecast of Convective Storms in Storm-Scale Models: Proof-of-Concept Experiments

Achieving accurate storm-scale analyses and reducing the spinup time of modeled convection is a primary motivation for the assimilation of radar reflectivity data. One common technique of reflectivity data assimilation is using a cloud analysis, which inserts temperature and moisture increments and hydrometeors deduced from radar reflectivity via empirical relations to induce and sustain updraft circulations. Polarimetric radar data have the ability to provide enhanced insight into the microphysical and dynamic structure of convection. Thus far, however, relatively little has been done to leverage these data for numerical weather prediction. In this study, the Advanced Regional Prediction System’s cloud analysis is modified from its original reflectivity-based formulation to provide moisture and latent heat adjustments based on the detection of differential reflectivity columns, which can serve as proxies for updrafts in deep moist convection and, subsequently, areas of saturation and latent heat release. Cycled model runs using both the original cloud analysis and above modifications are performed for two high-impact weather cases: the 19 May 2013 central Oklahoma tornadic supercells and the 25 May 2016 north-central Kansas tornadic supercell. The analyses and forecasts of convection qualitatively and quantitatively improve in both cases, including more coherent analyzed updrafts, more realistic forecast reflectivity structures, a better correspondence between forecast updraft helicity tracks and radar-derived rotation tracks, and improved frequency biases and equitable threat scores for reflectivity. Based on these encouraging results, further exploration of the assimilation of dual-polarization radar data into storm-scale models is warranted.

Impact of multiple radar reflectivity data assimilation on the numerical simulation of a flash flood event during the HyMeX campaign

Abstract. An analysis to evaluate the impact of multiple radar reflectivity data with a three-dimensional variational (3-D-Var) assimilation system on a heavy precipitation event is presented. The main goal is to build a regionally tuned numerical prediction model and a decision-support system for environmental civil protection services and demonstrate it in the central Italian regions, distinguishing which type of observations, conventional and not (or a combination of them), is more effective in improving the accuracy of the forecasted rainfall. In that respect, during the first special observation period (SOP1) of HyMeX (Hydrological cycle in the Mediterranean Experiment) campaign several intensive observing periods (IOPs) were launched and nine of which occurred in Italy. Among them, IOP4 is chosen for this study because of its low predictability regarding the exact location and amount of precipitation. This event hit central Italy on 14 September 2012 producing heavy precipitation and causing several cases of damage to buildings, infrastructure, and roads. Reflectivity data taken from three C-band Doppler radars running operationally during the event are assimilated using the 3-D-Var technique to improve high-resolution initial conditions. In order to evaluate the impact of the assimilation procedure at different horizontal resolutions and to assess the impact of assimilating reflectivity data from multiple radars, several experiments using the Weather Research and Forecasting (WRF) model are performed. Finally, traditional verification scores such as accuracy, equitable threat score, false alarm ratio, and frequency bias – interpreted by analysing their uncertainty through bootstrap confidence intervals (CIs) – are used to objectively compare the experiments, using rain gauge data as a benchmark.