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

No abstract

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

This work proposes a multi-parameter method for the detection of cloud-to-ground stroke rate (SR CG ) associated to convective cells, based on the measurements of a low-cost single-polarization X-band weather radar. To train and test our procedure, we built up a multi-year dataset, collecting 1575 radar reflectivity volumes that were acquired in the pilot study area of Naples metropolitan environment matched with the LIghtning NETwork (LINET) strokes and meteorological in-situ data. Three radar-based variables are extracted simultaneously for each rain cell and properly merged together, using "ad hoc" classification methods, to produce an estimation of the expected lightning activity for each rain cell. These variables, proxies of mixed-phase particles and ice amount into a convective cell, are combined into a single label to cluster the SR CG into two categories: SR CG = 0 (no production of strokes) or SR CG > 0 (stroke production), respectively. Overall, the main results are comparable with those that were obtained from more advanced radar systems, showing a Critical Success Index of 0.53, an Equitable Threat Score of 0.34, a Frequency Bias Index of 1.00, a Heidke Skill Score of 0.42, a Hanssen-Kuiper Skill Score of 0.42, and an area under the curve of probability of detection as a function of false alarm rate (usually referred as ROC curve) equal to 0.78. The developed technique, although with some limitations, outperforms those based on the use of single stroke proxy parameters. stroke discharges with a spatial resolution as high as up to 100 m and with a detection efficiency up to 95%, although these performances depend on the network density and the type of sensors.The real-time surveillance of stroke occurrence can also rely on weather radar measurements, which are able to track and characterize the three-dimensional (3D) structure of rain cells, thus allowing for identifying the developing cycle of cells and the areas much prone to stroke activity, even before the occurrence of the first lightning event (e.g., [5]). Therefore, the set-up of a reliable, affordable, and accessible radar-based stroke detection system, complementary to traditional ground-based stroke networks, can be very useful for risk prevention and for safety of human life, goods and services. In addition, a stand-alone radar-based stroke detection system could cover those areas where data from lightning networks are not freely accessible or where their detection efficiency levels are not constant over large domains, due to the irregular distribution of lightning sensors.This work is aimed at proposing a new algorithm for the radar-based detection of stroke activity based on a multi-variable approach.To explain our approach, it is useful to briefly summarize the atmosphere electrification mechanisms and the radar-based approaches so far proposed in the state of the art literature. Electrification mechanisms in thunderstorms are explained following the widely accepted Non-Inductive Charging (NIC) theory, whose evidences have been supported...

Abstract. The weather forecasts for precipitation have considerably improved in recent years thanks to the increase of computational power. This allows for the use of both a higher spatial resolution and the parameterization schemes specifically developed for representing sub-grid scale physical processes at high resolution. However, precipitation estimation is still affected by errors that can impact the response of hydrological models. To the aim of improving the hydrological forecast and the characterization of related uncertainties, a regional-scale meteorological–hydrological ensemble is presented. The uncertainties in the precipitation forecast and how they propagate in the hydrological model are also investigated. A meteorological–hydrological offline coupled ensemble is built to forecast events in a complex-orography terrain where catchments of different sizes are present. The Best Discharge-based Drainage (BDD; both deterministic and probabilistic) index, is defined with the aim of forecasting hydrological-stress conditions and related uncertainty. In this context, the meteorological–hydrological ensemble forecast is implemented and tested for a severe hydrological event which occurred over Central Italy on 15 November 2017, when a flood hit the Abruzzo region with precipitation reaching 200 mm (24 h)−1 and producing damages with a high impact on social and economic activities. The newly developed meteorological–hydrological ensemble is compared with a high-resolution deterministic forecast and with the observations (rain gauges and radar data) over the same area. The receiver operating characteristic (ROC) statistical indicator shows how skilful the ensemble precipitation forecast is with respect to both rain-gauge- and radar-retrieved precipitation. Moreover, both the deterministic and probabilistic configurations of the BDD index are compared with the alert map issued by Civil Protection Department for the event showing a very good agreement. Finally, the meteorological–hydrological ensemble allows for an estimation of both the predictability of the event a few days in advance and the uncertainty of the flood. Although the modelling framework is implemented on the basins of the Abruzzo region, it is portable and applicable to other areas.

Abstract. The weather forecasts for precipitation have considerably improved in recent years thanks to the increase of computational power. This allows to use both a higher spatial resolutions and the newly developed parameterization schemes for representing sub-grid scale physical processes. However, precipitation estimation is still affected by errors that can impact on the response of hydrological models. To the aim of considering the uncertainties in the precipitation forecast and how they propagate in the hydrological model, an ensemble approach is investigated. A meteo-hydro ensemble system is built to forecast events in a complex orography terrain where catchments of different size are present. In this context, the meteo-hydrological forecast system is implemented and tested for a severe hydrological event occurred over Central Italy on November 15, 2017. During this period, a flash flood hit the Abruzzo region causing precipitation up to 200 mm/24 hours and producing damages with a high impact on social and economic activities.The newly developed meteo-hydro ensemble system is compared with a high resolution deterministic forecast and with the observations over the same area, showing a very good response. In addition, the ensemble allows for an estimation of the predictability of the event a few days in advance and of the uncertainty of this flood. Although the modelling framework is implemented on the basins of Abruzzo region, it is portable and applicable to other areas.

Abstract. The precipitation forecast over the Mediterranean basin is still a challenge because of the complex orographic region which amplifies the need for local observation to correctly initialize the forecast. In this context the data assimilation techniques play a key role in improving the initial conditions and consequently the timing and position of precipitation pattern. For the first time, the ability of a cycling 4D-Var to reproduce a severe weather event in central Italy, as well as to provide a comparison with the largely used cycling 3D-Var, is evaluated in this study. The radar reflectivity measured by the Italian ground radar network is assimilated in the WRF model to simulate an event occurred on May 3, 2018 in central Italy. In order to evaluate the impact of data assimilation, several simulations are objectively compared by means of a Fraction Skill Score (FSS), which is calculated for several threshold values, and a Receiver Operating Characteristic (ROC) curve. The results suggest that both assimilation methods in cycling mode improve the 1, 3 and 6-hourly quantitative precipitation estimation. More specifically, the cycling 4D-Var with a warm start initialization shows the highest FSS values in the first hours of simulation both with light and heavy precipitation. Finally, the ROC curve confirms the benefit of 4D-Var: the area under the curve is 0.91 compared to the 0.88 of control experiment without data assimilation.

Abstract. This work exploits the potentiality of hail warning, based on single-polarization X-band weather radar measurements and tested on a large and well-documented data set of thunderstorm events in southern Italy near Naples. Even though X-band radars may suffer of two-way path attenuation especially at long ranges, due to their relatively low cost their use is rapidly 10 increasing for short-range applications such as urban environments. To identify hail through radar measurements, two different methodologies have been selected and adapted to X-band data within the study area: one uses the Waldvogel (WAL) approach, whereas the other one uses the Vertically-Integrated Liquid Density (VIL-Density) product. The study aims at developing a Probability-of-Hail (POH) index in order to support hail risk management at urban scales. In order to find the optimal threshold values to discriminate between hail and severe rain, an extensive intercomparison between outcomes of the two methodologies 15 and ground truth observations of hail has been performed, using a 2x2 contingency table and statistical scores.The results show that both methods are accurate for hail detection in the area of interest, although VIL-Density product is less satisfactory than WAL method in terms of false alarm ratio. The relationship between the output of these two methodologies and POH has been derived through a heuristic approach, using a third-order polynomial fitting curve. As an example, the POH indexes have been applied for the thunderstorm event occurred on 21 July 2014, proving to be reliable for hail core detection. 20