Comparison between 3D-Var and 4D-Var data assimilation methods for the simulation of a heavy rainfall case in central Italy

This work aims to provide a comparison between three dimensional and four dimensional variational data assimilation methods (3D-Var and 4D-Var) for a heavy rainfall case in central Italy. To evaluate the impact of the assimilation of reflectivity and radial velocity acquired from Monte Midia Doppler radar into the Weather Research Forecasting (WRF) model, the quantitative precipitation forecast (QPF) is used. The two methods are compared for a heavy rainfall event that occurred in central Italy on 14 September 2012 during the first Special Observation Period (SOP1) of the HyMeX (HYdrological cycle in Mediterranean EXperiment) campaign. This event, characterized by a deep low pressure system over the Tyrrhenian Sea, produced flash floods over the Marche and Abruzzo regions, where rainfall maxima reached more than 150 mm 24 h−1. To identify the best QPF, nine experiments are performed using 3D-Var and 4D-Var data assimilation techniques. All simulations are compared in terms of rainfall forecast and precipitation measured by the gauges through three statistical indicators: probability of detection (POD), critical success index (CSI) and false alarm ratio (FAR). The assimilation of conventional observations with 4D-Var method improves the QPF compared to 3D-Var. In addition, the use of radar measurements in 4D-Var simulations enhances the performances of statistical scores for higher rainfall thresholds.