Search

Reducing errors in streamflow simulations is one of the main issues for a reliable forecast system aimed to manage floods and water resources. Data assimilation is a powerful tool to reduce model errors. Unfortunately, its use in operational chains with distributed and physically based models is a challenging issue since many methodologies require computational times that are hardly compatible with operational needs. The implemented methodology corrects modelled water level in channels and root-zone soil moisture using real-time water level gauge stations. Model's variables are corrected locally, then the updates are propagated upstream with a simple approach that accounts for sub-basins' contributions. The overfitting issue, which arises when updating a spatially distributed model with sparse streamflow data, is hence here addressed in the context of a large-scale operational implementation working in real time thanks to the simplicity of the strategy. To test the method, a hindcast of daily simulations covering 18 months was performed on the Italian Tevere basin, and the modelling results with and without assimilation were compared. The setup was that currently in place in the operational framework in both cases. The analysis evidences a clear overall benefit of applying the proposed method even out of the assimilation time window. [ABSTRACT FROM AUTHOR]

The typical complex orography of the Mediterranean coastal areas support the formation of the so-called back-building mesoscale convective systems (MCS) producing torrential rainfall often resulting in flash floods. As these events are usually very small-scaled and localized, they are hardly predictable from a hydrometeorological standpoint, frequently causing a significant amount of fatalities and socioeconomic damage. Liguria, a northwestern Italian region, is characterized by small catchments with very short hydrological response time and is thus extremely prone to the impacts of back-building MCSs. Indeed, Liguria has been hit by three intense back-building MCSs between 2011 and 2014, causing a total death toll of 20 people and several hundred millions of euros of damages. Consequently, it is necessary to use hydrometeorological forecasting frameworks coupling the finescale numerical weather prediction (NWP) outputs with rainfall-runoff models to provide timely and accurate streamflow forecasts. Concerning the aforementioned back-buildingMCS episodes that recently occurred in Liguria, this work assesses the predictive capability of a hydrometeorological forecasting framework composed by a kilometer-scale cloud-resolving NWP model (WRF), including a 6-h cycling 3DVAR assimilation of radar reflectivity and conventional weather stations data, a rainfall downscaling model [Rainfall Filtered Autoregressive Model (RainFARM)], and a fully distributed hydrological model (Continuum). A rich portfolio of WRF 3DVAR direct and indirect reflectivity operators has been explored to drive the meteorological component of the proposed forecasting framework. The results confirm the importance of rapidly refreshing and data intensive 3DVAR for improving the quantitative precipitation forecast, and, subsequently, the flash flood prediction in cases of back-building MCS events. [ABSTRACT FROM AUTHOR]

The autumn of 2014 was characterized by a number of severe weather episodes over Liguria (northern Italy) associated with floods and remarkable damage. This period is selected as a test bed to evaluate the performance of a rainfall assimilation scheme based on the nudging of humidity profiles and applied to a convection-permitting meteorological model at high resolution. The impact of the scheme is assessed in terms of quantitative precipitation forecast (QPF) applying an object-oriented verification methodology that evaluates the structure, amplitude, and location (SAL) of the precipitation field, but also in terms of hydrological discharge prediction. To attain this aim, the meteorological model is coupled with the operational hydrological forecasting chain of the Ligurian Hydrometeorological Functional Centre, and the whole system is implemented taking operational requirements into account. The impact of rainfall data assimilation is large during the assimilation period and still relevant in the following 3 h of the free forecasts, but hardly lasts more than 6 h. However, this can improve the hydrological predictions. Moreover, the impact of the assimilation is dependent on the environment characteristics, being more effective when nonequilibrium convection dominates, and thus an accurate prediction of the local triggering for the development of the precipitation system is required. [ABSTRACT FROM AUTHOR]

Coupling meteorological and hydrological models is a common and standard practice in the field of flood forecasting. In this study, a numerical weather prediction (NWP) chain based on the BOLogna Limited Area Model (BOLAM) and the MOdello LOCale in Hybrid coordinates (MOLOCH) was coupled with the operational hydrological forecasting chain of the Ligurian Hydro-Meteorological Functional Centre to simulate two major floods that occurred during autumn 2011 in northern Italy. Different atmospheric simulations were performed by varying the grid spacing (between 1.0 and 3.0 km) of the high-resolution meteorological model and the set of initial/boundary conditions driving the NWP chain. The aim was to investigate the impact of these parameters not only from a meteorological perspective, but also in terms of discharge predictions for the two flood events. The operational flood forecasting system was thus used as a tool to validate in a more pragmatic sense the quantitative precipitation forecast obtained from different configurations of the NWP system. The results showed an improvement in flood prediction when a high-resolution grid was employed for atmospheric simulations. In turn, a better description of the evolution of the precipitating convective systems was beneficial for the hydrological prediction. Although the simulations underestimated the severity of both floods, the higher-resolution model chain would have provided useful information to the decision-makers in charge of protecting citizens. [ABSTRACT FROM AUTHOR]

Summary One of the main difficulties that flood forecasters are faced with is evaluating how errors and uncertainties in forecasted precipitation propagate into streamflow forecast. These errors, must be combined with the effects of different initial soil moisture conditions that generally have a significant impact on the final results of a flood forecast. This is further complicated by the fact that a probabilistic approach is needed, especially when small and medium size basins are considered (the variability of the streamflow scenarios is in fact strongly influenced by the aforementioned factors). Moreover, the ensemble size is a degree of freedom when a precipitation downscaling algorithm is part of the forecast chain. In fact, a change of ensemble size could lead to different final results once the other inputs and parameters are fixed. In this work, a series of synthetic experiments have been designed and implemented to test an operational probabilistic flood forecast system in order to augment the knowledge of how streamflow forecasts can be affected by errors and uncertainties associated with the three aforementioned elements: forecasted rainfall, soil moisture initial conditions, and ensemble size. [ABSTRACT FROM AUTHOR]

The forecast of rainfall-driven floods is one of the main themes of analysis in hydrometeorology and a critical issue for civil protection systems. This work describes a complete hydrometeorological forecast system for small- and medium-sized basins and has been designed for operational applications. In this case, because of the size of the target catchments and to properly account for uncertainty sources in the prediction chain, the authors apply a probabilistic framework. This approach allows for delivering a prediction of streamflow that is valuable for decision makers and that uses as input quantitative precipitation forecasts (QPF) issued by a regional center that is in charge of hydrometeorological predictions in the Liguria region of Italy. This kind of forecast is derived from different meteorological models and from the experience of meteorologists. Single-catchment and multicatchment approaches have been operationally implemented and studied. The hydrometeorological forecasting chain has been applied to a series of case studies with encouraging results. The implemented system makes effective use of the quantitative information content of rainfall forecasts issued by expert meteorologists for flood-alert purposes. [ABSTRACT FROM AUTHOR]

Polarimetric radars provide measurements that describe the shape and dimensions of hydrometeors and are unaffected by calibration, attenuation, and the presence of ice. These measurements can potentially lead to a more detailed description of hydrometeors and to an improvement in quantitative rainfall rate estimation. The authors present an algorithm that exploits polarimetric measures for rain-rate estimation and investigate its application in a real-time framework by using measurements from the C-band polarimetric radar at Monte Settepani in Savona, Italy. It is based on a flowchart decision tree that allows the use of the best rain-rate retrieval algorithm, depending on the value of polarimetric variables. The methodology was applied to a real-time framework for more than a year, and the results were presented for all the significant events observed during the test period. To evaluate the performance of the algorithm, a comparison was made with rain gauge observation from a dense regional network. The performances of the algorithm were compared with those obtained by standard operational Z–R formulations to evaluate the benefit of this approach for operational applications. [ABSTRACT FROM AUTHOR]

In recent years, the Liguria Region of Italy has been affected by several flash floods that have caused significant losses in terms of damage to infrastructure and human lives. The forecast of this kind of events is a critical issue for civil protection. The suitable approach to deal with this problem starts from the most accurate forecast of the rainfall field to be the input of the hydrological model. For this purpose, it is possible to take advantage of the output of the nowcasting, that is a forecast at short lead time (0 to 6 hours), but also of the forecast of the Numerical Weather Prediction Systems. The realization of an integrated hydrological nowcasting chain coupling the combination of different rainfall fields and the hydrological model can improve the discharge forecast. Within this study this coupling is performed joining in a unique tool the nowcasting technique PhaSt, a spectral-based nowcasting procedure, the NWPS MOLOCH corrected with data assimilation and the hydrological model Continuum, a continuous distributed hydrological model. The connection between the output of the forecasting models is performed through a technique called blending, that tries to linearly combine the rainfall fields according to their reliability function of the lead time. The blending has been approached in a novel way in this study exploiting the information retrieved from the NWPS both in terms of indication regarding the variation of total volume on the domain considered and in terms of location of the rainfall structures. An evaluation of the performances of the chain has been performed on some case events of autumn 2014 all over Liguria Region, northern Italy. A distributed analysis, executed on the entire domain on which the hydrological model is computed, has been possible obtaining representative results for all the basins involved in the events (more than 600 hydrological sections have been considered). According to the results obtained with the case studies it remains an open challenge the definition of a blending function that takes the best from every forecasting method. [ABSTRACT FROM AUTHOR]

In Italy, rainfall represents the most common triggering factor for landslides; thus, many Italian Regional Departments of Civil Protection are setting up warning systems based on rainfall recordings. Common methods are mainly based on empirical relationships that provide the rainfall thresholds above which the occurrence of landslide phenomena is likely to be expected. In recent years, the use of machine learning approaches has gained popularity in landslide susceptibility analysis and prediction. To support the operational early warning system of Liguria Civil Protection Department for landslides hazard, we propose the implementation of a polynomial Kernel regularized least squares regression (KRLS) algorithm, for predicting the daily occurrence of shallow landslides in the five Alert Zones in Liguria (North Western Italy). The model provides an estimate of the number of landslides associated with the set of three different hydrological features, also used for the Hydrological Assessment procedure: the soil moisture, the accumulated precipitation over 12 h and the precipitation peak over 3 h. Results of the model are converted to an Alert Scenario of landslide occurrence, based on the magnitude of the expected event and identified according to the National and Regional legislation (Regional Civil Protection guidelines D.G.R. n. 1116, 23/12/2020). The performance of the predictive model (e.g. accuracy of 93%) is deemed satisfactory and the methodology is considered a valuable support to the operational early warning system of Liguria Civil Protection Department. The choice of predictive variables allows, in future development, the values obtained from historical data to be replaced by those obtained from meteorological forecast models, introducing the use of the developed model in the operational forecasting chain. [ABSTRACT FROM AUTHOR]

We present IT-SNOW, a serially complete and multi-year snow reanalysis for Italy (∼ 301 × 10 3 km2) – a transitional continental-to-Mediterranean region where snow plays an important but still poorly constrained societal and ecological role. IT-SNOW provides ∼ 500 m daily maps of snow water equivalent (SWE), snow depth, bulk snow density, and liquid water content for the initial period 1 September 2010–31 August 2021, with future updates envisaged on a regular basis. As the output of an operational chain employed in real-world civil protection applications (S3M Italy), IT-SNOW ingests input data from thousands of automatic weather stations, snow-covered-area maps from Sentinel-2, MODIS (Moderate Resolution Imaging Spectroradiometer), and H SAF products, as well as maps of snow depth from the spatialization of over 350 on-the-ground snow depth sensors. Validation using Sentinel-1-based maps of snow depth and a variety of independent, in situ snow data from three focus regions (Aosta Valley, Lombardy, and Molise) show little to no mean bias compared to the former, and root mean square errors are of the typical order of 30–60 cm and 90–300 mm for in situ, measured snow depth and snow water equivalent, respectively. Estimates of peak SWE by IT-SNOW are also well correlated with annual streamflow at the closure section of 102 basins across Italy (0.87), with ratios between peak water volume in snow and annual streamflow that are in line with expectations for this mixed rain–snow region (22 % on average and 12 % median). Examples of use allowed us to estimate 13.70 ± 4.9 Gm3 of water volume stored in snow across the Italian landscape at peak accumulation, which on average occurs on 4 March ± 10 d. Nearly 52 % of the mean seasonal SWE is accumulated across the Po river basin, followed by the Adige river (23 %), and central Apennines (5 %). IT-SNOW is freely available at 10.5281/zenodo.7034956 and can contribute to better constraining the role of snow for seasonal to annual water resources – a crucial endeavor in a warming and drier climate. [ABSTRACT FROM AUTHOR]

By shifting winter precipitation into summer freshet, the cryosphere supports life across the world. The sensitivity of this mechanism to climate and the role played by the cryosphere in the Earth's energy budget have motivated the development of a broad spectrum of predictive models. Such models represent seasonal snow and glaciers with various complexities and generally are not integrated with hydrologic models describing the fate of meltwater through the hydrologic budget. We present Snow Multidata Mapping and Modeling (S3M) v5.1, a spatially explicit and hydrology-oriented cryospheric model that simulates seasonal snow and glacier evolution through time and that can be natively coupled with distributed hydrologic models. Model physics include precipitation-phase partitioning, snow and glacier mass balances, snow rheology and hydraulics, a hybrid temperature-index and radiation-driven melt parametrization, and a data-assimilation protocol. Comparatively novel aspects of S3M are an explicit representation of the spatial patterns of snow liquid-water content, the implementation of the Δh parametrization for distributed ice-thickness change, and the inclusion of a distributed debris-driven melt factor. Focusing on its operational implementation in the northwestern Italian Alps, we show that S3M provides robust predictions of the snow and glacier mass balances at multiple scales, thus delivering the necessary information to support real-world hydrologic operations. S3M is well suited for both operational flood forecasting and basic research, including future scenarios of the fate of the cryosphere and water supply in a warming climate. The model is open source, and the paper comprises a user manual as well as resources to prepare input data and set up computational environments and libraries. [ABSTRACT FROM AUTHOR]

We present IT-SNOW, a serially complete and multi-year snow reanalysis for Italy (300k+ km²) covering a transitional continental-to-Mediterranean region where snow plays an important, but still poorly constrained societal and ecological role. IT-SNOW provides ∼500-m, daily maps of Snow Water Equivalent (SWE), snow depth, bulk-snow density, and liquid water content for the 5 period 01/09/2010 - 31/08/2021, with future updates envisaged on a regular basis. As the output of an operational chain employed in real-world civil-protection applications (S3M Italy), IT-SNOW ingests input data from thousands of automatic weather stations, snow-covered-area maps from Sentinel 2, MODIS, and H-SAF products, and maps of snow depth from the spazialization of 350+ on-the-ground snow-depth sensors. Validation using Sentinel-1-based maps of snow depth and a variety of independent, in-situ snow data from three focus regions (Aosta Valley, Lombardia, and Molise) shows little to none mean bias compared to the former, and Root Mean Square Errors on the order of 30 to 60 cm and 90 to 300 mm for in-situ, measured snow depth and Snow Water Equivalent, respectively. Estimates of peak SWE by IT-SNOW are also well correlated with annual streamflow at the closure section of 102 basins across Italy (0.87), with ratios between peak SWE and annual streamflow that are in line with expectations for this mixed rain-snow region (22% on average). Examples of use allowed us to estimate 13.70 ± 4.9 Gm3 of SWE across the Italian landscape at peak accumulation, which on average occurs on the 4th of March. Nearly 52% of mean seasonal SWE is accumulated across the Po river basin, followed by the Adige river (23%), and central Apennines (5%). IT-SNOW is freely available with the following DOI: https://doi.org/10.5281/zenodo.7034956 (Avanzi et al., 2022b) and can contribute to better constraining the role of snow for seasonal to annual water resources - a crucial endevor in a warming and drier climate. [ABSTRACT FROM AUTHOR]

The Hydro-Meteorological Centre (CMI) of the Environmental Protection Agency of Liguria Region, Italy, is in charge of the hydrometeorological forecast and the in-event monitoring for the region. This region counts numerous small and very small basins, known for their high sensitivity to intense storm events, characterised by low predictability. Therefore, at the CMI, a radar-based nowcasting modelling chain called the Small Basins Model Chain, tailored to such basins, is employed as a monitoring tool for civil protection purposes. The aim of this study is to evaluate the performance of this model chain, in terms of: (1) correct forecast, false alarm and missed alarm rates, based on both observed and simulated discharge threshold exceedances and observed impacts of rainfall events encountered in the region; (2) warning times respect to discharge threshold exceedances. The Small Basins Model Chain is proven to be an effective tool for flood nowcasting and helpful for civil protection operators during the monitoring phase of hydrometeorological events, detecting with good accuracy the location of intense storms, thanks to the radar technology, and the occurrence of flash floods. [ABSTRACT FROM AUTHOR]

Hydrological monitoring systems relying on radar data and distributed hydrological models are now feasible at large-scale and represent effective early warning systems for flash floods. Here we describe a system that allows hydrological occurrences in terms of streamflow at a national scale to be monitored. We then evaluate its operational application in Italy, a country characterized by various climatic conditions and topographic features. The proposed system exploits a modified conditional merging (MCM) algorithm to generate rainfall estimates by blending data from national radar and rain-gauge networks. Then, we use the merged rainfall fields as input for the distributed and continuous hydrological model, Continuum, to obtain real-time streamflow predictions. We assess its performance in terms of rainfall estimates from MCM, using cross-validation and comparison with a conditional merging technique at an event-scale. We also assess its performance against rainfall fields from ground-based data at catchment-scale. We further evaluate the performance of the hydrological system in terms of streamflow against observed data (relative error on high flows less than 25% and Nash–Sutcliffe Efficiency greater than 0.5 for 72% and 46% of the calibrated study sections, respectively). These results, therefore, confirm the suitability of such an approach, even at national scale, over a wide range of catchment types, climates, and hydrometeorological regimes, and for operational purposes. [ABSTRACT FROM AUTHOR]

By shifting winter precipitation into summer freshet, the cryosphere supports life across the world. The sensitivity of this shifting mechanism to climate, as well as the role played by the cryosphere in the Earth energy budget, has motivated the development of a broad spectrum of predictive models. Such models rarely combine a high degree of physical realism in both the seasonal snow and glaciers, and generally are not integrated with hydrologic models describing the fate of meltwater through the hydrologic budget. We present S3M v5.1, a spatially explicit and hydrology-oriented cryospheric model that successfully reconstructs seasonal snow and glacier evolution through time and that can be natively coupled with distributed hydrologic models. Model physics include precipitation-phase partitioning, snow and glacier energy and mass balances, snow rheology and hydraulics, and a data-assimilation protocol. Comparatively novel aspects of S3M with respect to the existing literature are an explicit representation of the spatial patterns of snow liquid-water content, an hybrid approach to snowmelt that decouples the radiation- and temperature-driven contributions, the implementation of the ∆h parametrization for distributed ice-thickness change, and the inclusion of a distributed debris-driven melt factor. Focusing on its operational implementation in the Italian north-western Alps, we show that S3M provides robust predictions of the snow and glacier mass balances at multiple scales, thus delivering the necessary information to support real-world hydrologic operations. S3M is well suited for both operational flood forecasting and basic research, including future scenarios of the fate of the cryosphere and water supply in a warming climate. The model is open source, and the paper comprises an user manual as well as resources to prepare input data and set up computational environments and libraries. [ABSTRACT FROM AUTHOR]

• Two flash flood early warning systems are compared analysing three extreme events. • Hazard is assessed in terms of return period for the whole Liguria region. • Weather radar rainfall is used as input for the warning systems. The anticipation of flash flood events is crucial to issue warnings to mitigate their impact. This work presents a comparison of two early warning systems for real-time flash flood hazard forecasting at regional scale. The two systems are based in a gridded drainage network and they use weather radar precipitation inputs to assess the hazard level in different points of the study area, considering the return period (in years) as the indicator of the flash flood hazard. The essential difference between the systems is that one is a rainfall-based system (ERICHA), using the upstream basin-aggregated rainfall as the variable to determine the hazard level, while the other (Flood-PROOFS) is a system based on a distributed rainfall-runoff model to compute the streamflows at pixel scale. The comparison has been done for three rainfall events in the autumn of 2014 that resulted in severe flooding in the Liguria region (Northwest of Italy). The results obtained by the two systems show many similarities, particularly for larger catchments and for large return periods (extreme floods). [ABSTRACT FROM AUTHOR]