Search

The knowledge of snowpack dynamics is of critical importance to several real-time applications especially in mountain basins, such as agricultural production, water resource management, flood prevention, hydropower generation. Since simulations are affected by model biases and forcing data uncertainty, an increasing interest focuses on the assimilation of snow-related observations with the purpose of enhancing predictions on snowpack state. The study aims at investigating the effectiveness of snow multivariable data assimilation (DA) at an Alpine site. The system consists of a snow energy-balance model strengthened by a multivariable DA system. An Ensemble Kalman Filter (EnKF) scheme allows assimilating ground-based and remotely sensed snow observations in order to improve the model simulations. This research aims to investigate and discuss: (1) the limitations and constraints in implementing a multivariate EnKF scheme in the framework of snow modelling, and (2) its performance in consistently updating the snowpack state. The performance of the multivariable DA is shown for the study case of Torgnon station (Aosta Valley, Italy) in the period June 2012 - December 2013. The results of several experiments are discussed with the aim of analyzing system sensitivity to the DA frequency, the ensemble size, and the impact of assimilating different observations.

Every year Africa is hit by extreme floods which, combined with high levels of vulnerability and increasing population exposure, often result in humanitarian crises and population displacement. Impact-based forecasting and early warning for natural hazards is recognized as a step forward in disaster risk reduction, thanks to its focus on people, livelihoods, and assets at risk. Yet, the majority of the African population is not covered by any sort of early warning system. This article describes the setup and the methodological approach of Flood-PROOFS East Africa, an impact-based riverine flood forecasting and early warning system for the Greater Horn of Africa (GHA), with a forecast range of 5 d. The system is based on a modeling cascade relying on distributed hydrological simulations forced by ensemble weather forecasts, link to inundation maps for specific return period, and application of a risk assessment framework to estimate population and assets exposed to upcoming floods. The system is operational and supports the African Union Commission and the Disaster Operation Center of the Intergovernmental Authority on Development (IGAD) in the daily monitoring and early warning from hydro-meteorological disasters in eastern Africa. Results show a first evaluation of the hydrological reanalysis at 78 river gauging stations and a semi-quantitative assessment of the impact forecasts for the catastrophic floods in Sudan and in the Nile River basin in summer 2020. More extensive quantitative evaluation of the system performance is envisaged to provide its users with information on the model reliability in forecasting extreme events and their impacts. [ABSTRACT FROM AUTHOR]

Climate change and water security are among the grand challenges of the 21st century, but literacy on these matters among high-school students is often unsystematic and/or detached from the real world. This study aims to introduce the educational objectives, methods, and early results of "Water and Us", a three-module initiative that can contribute to advancing water education in a warming climate by focusing on the natural and anthropogenic water cycle, climate change, and emerging water conflicts. The method of Water and Us revolves around storytelling to aid understanding and generate new knowledge, learning by doing, a flipped-classroom environment, and a constant link to examples from the real world (such as ongoing droughts across the world or seeds of conflict regarding transnational river basins). Water and Us was established in 2021–2022 and, during that school year, involved ≥200 students as part of a proof of concept to test the complete didactic approach using small-scale experiments. Results from ≥40 h of proof-of-concept events confirmed the effectiveness of this approach with respect to conveying the essential elements of the natural and anthropogenic water cycle, the most commonly recurring concepts related to climate change and water as well as the possible conflicts and solutions related to water scarcity in a warming climate. The Water and Us team remains interested in networking with colleagues and potential recipients to upscale and further develop this work. [ABSTRACT FROM AUTHOR]

State-of-the-art snow sensing technologies currently provide an unprecedented amount of data from both remote sensing and ground sensors, but their assimilation into dynamic models is bounded to data quality, which is often low – especially in mountain, high-elevation, and unattended regions where snow is the predominant land-cover feature. To maximize the value of snow-depth measurements, we developed a random forest classifier to automatize the quality assurance and quality control (QA/QC) procedure of near-surface snow-depth measurements collected through ultrasonic sensors, with particular reference to the differentiation of snow cover from grass or bare-ground data and to the detection of random errors (e.g., spikes). The model was trained and validated using a split-sample approach of an already manually classified dataset of 18 years of data from 43 sensors in Aosta Valley (northwestern Italian Alps) and then further validated using 3 years of data from 27 stations across the rest of Italy (with no further training or tuning). The F1 score was used as scoring metric, it being the most suited to describe the performances of a model in the case of a multiclass imbalanced classification problem. The model proved to be both robust and reliable in the classification of snow cover vs. grass/bare ground in Aosta Valley (F1 values above 90 %) yet less reliable in rare random-error detection, mostly due to the dataset imbalance (samples distribution: 46.46 % snow, 49.21 % grass/bare ground, 4.34 % error). No clear correlation with snow-season climatology was found in the training dataset, which further suggests the robustness of our approach. The application across the rest of Italy yielded F1 scores on the order of 90 % for snow and grass/bare ground, thus confirming results from the testing region and corroborating model robustness and reliability, with again a less skillful classification of random errors (values below 5 %). This machine learning algorithm of data quality assessment will provide more reliable snow data, enhancing their use in snow models. [ABSTRACT FROM AUTHOR]

Every year Africa is hit by extreme floods which, combined with high levels of vulnerability and increasing population exposure, often result in humanitarian crises and population displacement. Impact-based forecasting and early warning for natural hazards is recognized as a step forward in disaster risk reduction, thanks to its focus on people, livelihoods and assets at risk. Yet, the majority of the African population is not covered by any sort of early warning system. This article describes the setup of Flood-PROOFS East Africa, an impact-based riverine flood forecasting and early warning system for the Greater Horn of Africa (GHA), with a forecast range of 5 days. The system is based on a modeling cascade relying on distributed hydrological simulations forced by ensemble weather forecasts, link to inundation maps for specific return period, and application of a risk assessment framework to estimate population and assets exposed to upcoming floods. The system is operational and supports the African Union Commission and the IGAD Disaster Operation Center in the daily monitoring and early warning from hydro-meteorological disasters in Eastern Africa. Results show a first evaluation of the hydrological reanalysis at 78 river gauging stations and a semi-quantitative assessment of the impact forecasts for the catastrophic floods in Sudan and in the Nile River Basin in Summer 2020. [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]

Hydrological models often have issues in simulating streamflow (Q) during droughts, because of hard-to-capture feedback mechanisms across precipitation deficit, actual evapotranspiration (ET), and terrestrial water storage anomalies (TWSA). To gain more insights into these performance drops and move toward more robust hydrological models in the anthropogenic era, we evaluated Q, ET, and TWSA simulations during droughts of different severity and their sensitivity to the climatic conditions of the calibration period. We used the distributed hydrological model Continuum over the heavily human-affected Po river basin (northern Italy, period 2010-2022) and independent ground- and remote sensing-based datasets of Q, ET, and TWSA as benchmarks. Across the 38 study sub-catchments, Continuum simulated Q comparably well during wet years (2014 and 2020) and moderate droughts (2012 and 2017) with mean KGE = 0.59±0.32 during wet years and = 0.55±0.25 during moderate droughts. The model simulated well Q for the outlet section of the basin also for the severe 2022 drought (KGE = 0.82). However, performances for 2022 declined across the other sub-catchments (mean KGE = 0.18±0.69, meaning the model still preserved some skill over a climatological mean). The model properly represented seasonality of Q, ET, and TWSA over the basin, as well as a declining trend in TWSA. We explained the performance drops in 2022 with an increased uncertainty in ET anomalies, in particular in human-affected croplands. Calibrating during a moderate drought (2017) did not improve model performances during the severe 2022 drought (mean KGE = 0.18±0.63), pointing to the fairly unique conditions of this period in terms of hydrological processes and human interference on the hydrological cycle. By highlighting increased uncertainty of hydrological models specifically during severe droughts which are expected to increase in frequency, these findings provide relevant guidelines for assessments of model robustness in a changing climate and so for informing water management, disaster risk reduction, and climate change adaptation strategies. [ABSTRACT FROM AUTHOR]

Satellite-based Earth observations (EO) are an accurate and reliable data source for atmospheric and environmental science. Their increasing spatial and temporal resolutions, as well as the seamless availability over ungauged regions, make them appealing for hydrological modeling. This work shows recent advances in the use of high-resolution satellite-based EO data in hydrological modeling. In a set of six experiments, the distributed hydrological model Continuum is set up for the Po River basin (Italy) and forced, in turn, by satellite precipitation and evaporation, while satellite-derived soil moisture (SM) and snow depths are ingested into the model structure through a data-assimilation scheme. Further, satellite-based estimates of precipitation, evaporation, and river discharge are used for hydrological model calibration, and results are compared with those based on ground observations. Despite the high density of conventional ground measurements and the strong human influence in the focus region, all satellite products show strong potential for operational hydrological applications, with skillful estimates of river discharge throughout the model domain. Satellite-based evaporation and snow depths marginally improve (by 2 % and 4 %) the mean Kling–Gupta efficiency (KGE) at 27 river gauges, compared to a baseline simulation (KGE mean= 0.51) forced by high-quality conventional data. Precipitation has the largest impact on the model output, though the satellite data on average shows poorer skills compared to conventional data. Interestingly, a model calibration heavily relying on satellite data, as opposed to conventional data, provides a skillful reconstruction of river discharges, paving the way to fully satellite-driven hydrological applications. [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]

Climate change and water security are among the grand challenges of the 21st century, but literacy on these matters among high-school students is often unsystematic and/or far from the real world. To contribute advancing education in a warming climate and prepare next generations to play their role in future societies, we designed 'Water and Us', a three-module initiative focusing on the natural and anthropogenic water cycle, climate change, and conflicts. The method of Water and Us resolves around storytelling to aid understanding and generate new knowledge, learning by doing, a flipped classroom environment, and a constant link to the real world – such as the archetypal events of the California snow drought or the seeds of conflicts around transnational river basins. Water and Us was established in 2021, and since then has involved 200+ students in a proof of concept to test the didactic approach in small-scale experiments. Results from 40+ hours of events confirm that students are generally aware of climate change (90 %), but have sparse knowledge of the concrete actions that are in place to mitigate or adapt (up to 20 %). Understanding of the water cycle by students is often anchored to a naturalistic, but fictitious view where human interference is minimal. Our approach conveys key elements of the contemporary, natural/anthropogenic water cycle, how this cycle is challenged by warmer temperatures and declining snowpacks, and how education can contribute to avoiding maladaptation and conflicts. While this initiative is being channelled in awareness projects at various levels, the Water and Us team remains interested in networking with colleagues and potential recipients to scale up and further develop this work. [ABSTRACT FROM AUTHOR]

The vulnerability of flood-prone areas is determined by the susceptibility of the exposed assets to the hazard. It is a crucial component in risk assessment studies, both for climate change adaptation and disaster risk reduction. In this study, we analyse patterns of vulnerability for the residential sector in a frequently hit urban area of Milan, Italy. The conceptual foundation for a quantitative assessment of the structural dimensions of vulnerability is based on the modified source–pathway–receptor–consequence model. This conceptual model is used to improve the parameterization of the flood risk analysis, describing (i) hazard scenario definitions performed by hydraulic modelling based on past event data (source estimation) and morphological features and land-use evaluation (pathway estimation) and (ii) the exposure and vulnerability assessment which consists of recognizing elements potentially at risk (receptor estimation) and event losses (consequence estimation). We characterized flood hazard intensity on the basis of variability in water depth during a recent event and spatial exposure also as a function of a building's surroundings and buildings' intrinsic characteristics as a determinant vulnerability indicator of the elements at risk. In this sense the use of a geographic scale sufficient to depict spatial differences in vulnerability allowed us to identify structural vulnerability patterns to inform depth–damage curves and calculate potential losses from mesoscale (land-use level) to microscale (building level). Results produces accurate estimates of the flood characteristics, with mean error in flood depth estimation in the range 0.2–0.3 m and provide a basis to obtain site-specific damage curves and damage mapping. Findings show that the nature of flood pathways varies spatially, is influenced by landscape characteristics and alters vulnerability spatial distribution and hazard propagation. At the mesoscale, the "continuous urban fabric" Urban Atlas 2018 land-use class with the occurrence of at least 80 % of soil sealing shows higher absolute damage values. At microscale, evidence demonstrated that even events with moderate magnitude in terms of flood depth in a complex urbanized area may cause more damage than one would expect. [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]

In a warming climate, periods with lower than average precipitation will increase in frequency and intensity. During such periods, known as meteorological droughts, the decline in annual runoff may be proportionally larger than the corresponding decline in precipitation. Reasons behind this exacerbation of runoff deficit during dry periods remain largely unknown, and this challenges the predictability of when this exacerbation will occur in the future and how intense it will be. In this work, we tested the hypothesis that runoff deficit exacerbation during droughts is a common feature across climates, driven by evaporation enhancement. We relied on multidecadal records of streamflow and precipitation for more than 200 catchment areas across various European climates, which distinctively show the emergence of similar periods of exacerbated runoff deficit identified in previous studies, i.e. runoff deficit on the order of -20 % to -40 % less than what expected from precipitation deficits. The magnitude of this exacerbation is two to three times larger for basins located in dry regions than for basins in wet regions, and is qualitatively correlated with an increase in annual evaporation during droughts, in the order of +11 % and +33 % over basins characterized by energy-limited and water-limited evaporation regimes, respectively. Thus, enhanced atmospheric and vegetation demand for moisture during dry periods induces a nonlinear precipitation-runoff relationship for low-flow regimes, which results in an unexpectedly large decrease in runoff during periods of already low water availability. Forecasting onset, magnitude, and duration of these drops in runoff have paramount societal and ecological implications, especially in a warming climate, given their supporting role for safeguarding water, food, and energy. The outcome that water basins are prone to this exacerbation of runoff deficit for various climates and evaporation regimes makes further understanding of its patterns of predictability an urgent priority for water-resource planning and management in a warming and drier climate. [ABSTRACT FROM AUTHOR]

The vulnerability of flood-prone areas is determined by the susceptibility of the exposed assets to the hazard. It is a crucial component in risk assessment studies, both for climate change adaptation and disaster risk reduction. In this study, we analyse patterns of vulnerability for the residential sector in a frequently hit urban area of Milano, Italy. The conceptual foundation for a quantitative assessment of the structural dimensions of vulnerability is based on the modified Source-Pathway-Receptor-Consequence model. This conceptual model is used to improve the parameterization of the flood risk analysis describing: (i) hazard scenarios definition performed by hydraulic modelling based on past event data (Source estimation) and morphological features and land use evaluation (Pathway estimation); (ii) the exposure and vulnerability assessment which consists of recognizing elements potentially at risk (Receptor estimation) and event losses (Consequence estimation). The structural dimension of vulnerability is mapped at building level and used in loss estimation for the residential sector at meso and micro-scale. Results produces accurate estimates of the flood characteristics, with mean error in flood depths estimation in the range 0.2-0.3 m and provide a basis to obtain site-specific damage curves and damage mapping. Findings show that the nature of flood pathways varies spatially and is influenced by landscape characteristics and alters vulnerability spatial distribution and hazard propagation. At the mesoscale, the 'Continuous urban fabric' Urban Atlas 2018 land-use class with the occurrence of at least 80 % of soil sealing shows higher absolute damage values. At microscale, evidence demonstrated that even events with moderate magnitude in terms of flood depth in a complex urbanized area may cause more damage than it would expect. [ABSTRACT FROM AUTHOR]

Satellite Earth observations (EO) are an accurate and reliable data source for atmospheric and environmental science. Their increasing spatial and temporal resolution, as well as the seamless availability over ungauged regions, make them appealing for hydrological modeling. This work shows recent advances in the use of high- resolution satellite-based Earth observation data in hydrological modelling. In a set of experiments, the distributed hydrological model Continuum is set up for the Po River Basin (Italy) and forced, in turn, by satellite precipitation and evaporation, while satellite-derived soil moisture and snow depths are ingested into the model structure through a data-assimilation scheme. Further, satellite- based estimates of precipitation, evaporation and river discharge are used for hydrological model calibration, and results are compared with those based on ground observations. Despite the high density of conventional ground measurements and the strong human influence in the focus region, all satellite products show strong potential for operational hydrological applications, with skillful estimates of river discharge throughout the model domain. Satellite- based evaporation and snow depths marginally improve (by 2 % and 4 %) the mean Kling-Gupta efficiency (KGE) at 27 river gauges, compared to a baseline simulation (KGEmean = 0.51) forced by high-quality conventional data. Precipitation has the largest impact on the model output, though the satellite dataset on average shows poorer skills compared to conventional data. Interestingly, a model calibration heavily relying on satellite data, as opposed to conventional data, provides a skillful reconstruction of river discharges, paving the way to fully satellite-driven hydrological applications. [ABSTRACT FROM AUTHOR]

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]

In a warming climate, periods with below-than-average precipitation will increase in frequency and intensity. During such periods, known as meteorological droughts, sparse but consistent pieces of evidence show that the decline in annual runoff may be proportionally larger than the corresponding decline in precipitation (e.g., -40% vs. -20%). Reasons behind this exacerbation of runoff deficit during dry periods remain largely unknown, which challenges generalization at larger scales (i.e., beyond the single catchment), as well as the predictability of when this exacerbation will occur and how intense it will be. Here, we tested the hypothesis that runoff-deficit exacerbation during droughts is a common feature of droughts across climates and is driven by evapotranspiration enhancement.We support this hypothesis by relying on multidecadal records of streamflow and precipitation for more than 200 catchments across various European climates, which distinctively show the emergence of similar periods of exacerbated runoff deficit identified in previous studies, i.e., runoff deficit on the order of -20% to -40% less than what expected from precipitation deficit. The magnitude of this exacerbation is two to three times larger for basins located in dry regions than for basins in wet regions and is qualitatively correlated with an increase in annual evapotranspiration during droughts, on the order of 11% and 33% over basins characterized by energy- and water-limited evapotranspiration regimes, respectively. Thus, enhanced atmospheric and vegetation demand for moisture during dry periods induces a nonlinear and potentially hysteretic precipitation-runoff relationship for low-flow regimes, which results in an unexpectedly large decrease in runoff during periods of already low water availability. Forecasting onset, magnitude, and duration of these drops in runoff availability has paramount societal implications, especially in a warming climate, given their supporting role for water, food, and energy security. The outcome that water basins are prone to this exacerbation of runoff deficit for various climates and evapotranspiration regimes, compounded by the lack of specific parametrizations of this process in the majority of hydrological and land-surface models, make further understanding of its patterns of predictability an urgent priority for water-resource planning and management 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 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]

Precipitation orographic enhancement is the result of both synoptic circulation and topography. Since high-elevation headwaters are often sparsely instrumented, the magnitude and distribution of this enhancement, as well as how they affect precipitation lapse rates, remain poorly understood. Filling this knowledge gap would allow a significant step ahead for hydrologic forecasting procedures and water management in general. Here, we hypothesized that spatially distributed, manual measurements of snow depth (courses) could provide new insights into this process. We leveraged over 11 000 snow course data upstream of two reservoirs in the western European Alps (Aosta Valley, Italy) to estimate precipitation orographic enhancement in the form of lapse rates and, consequently, improve predictions of a snow hydrologic modeling chain (Flood-PROOFS). We found that snow water equivalent (SWE) above 3000 m a.s.l. (above sea level) was between 2 and 8.5 times higher than recorded cumulative seasonal precipitation below 1000 m a.s.l., with gradients up to 1000 mm w.e. km -1. Enhancement factors, estimated by blending precipitation gauge and snow course data, were consistent between the two hydropower headwaters (median values above 3000 m a.s.l. between 4.1 and 4.8). Including blended gauge course lapse rates in an iterative precipitation spatialization procedure allowed Flood-PROOFS to remedy underestimations both of SWE above 3000 m a.s.l. (up to 50 %) and – importantly – of precipitation vs. observed streamflow. Annual runoff coefficients based on blended lapse rates were also more consistent from year to year than those based on precipitation gauges alone (standard deviation of 0.06 and 0.19, respectively). Thus, snow courses bear a characteristic signature of orographic precipitation, which opens a window of opportunity for leveraging these data sets to improve our understanding of the mountain water budget. This is all the more important due to the essential role of high-elevation headwaters in supporting water security and ecosystem services worldwide. [ABSTRACT FROM AUTHOR]

Precipitation orographic enhancement depends on both synoptic circulation and topography. Since high-elevation headwaters are often sparsely instrumented, the magnitude and distribution of this enhancement remain poorly understood. Filling this knowledge gap would allow a significant step ahead for hydrologic-forecasting procedures and water management in general. Here, we hypothesized that spatially distributed, manual measurements of snow depth (courses) could provide new insights into this process. We leveraged 11,000+ snow-course data upstream two reservoirs in the Western European Alps (Aosta Valley, Italy) to estimate precipitation orographic enhancement in the form of lapse rates and consequently improve predictions of a snow-hydrologic modeling chain (Flood- PROOFS). We found that Snow Water Equivalent (SWE) above 3000 m ASL was between 2 and 8.5 times higher than recorded cumulative seasonal precipitation below 1000 m ASL, with gradients up to 1000 mm w.e. km-1. Enhancement factors estimated by blending precipitation-gauge and snow-course data were quite consistent between the two hydropower headwaters (median values above 3000 m ASL between 4.1 and 4.8). Including blended gauge-course lapse rates in an iterative precipitation-spatialization procedure allowed Flood-PROOFS to remedy underestimations of both SWE above 3000 m ASL (up to 50 %) and importantly precipitation vs. observed streamflow. Runoff coefficients based on blended lapse rates were also more consistent from year to year that those based on precipitation gauges alone (standard deviation of 0.06 and 0.19, respectively). Thus, snow courses bear a characteristic signature of orographic precipitation, which opens a window of opportunity for leveraging these data sets to improve our understanding of the mountain water budget. This is all the more important due to their essential role in supporting water security and ecosystem services worldwide. [ABSTRACT FROM AUTHOR]

Snow models are usually evaluated at sites providing high-quality meteorological data, so that the uncertainty in the meteorological input data can be neglected when assessing model performances. However, high-quality input data are rarely available in mountain areas and, in practical applications, the meteorological forcing used to drive snow models is typically derived from spatial interpolation of the available in situ data or from reanalyses, whose accuracy can be considerably lower. In order to fully characterize the performances of a snow model, the model sensitivity to errors in the input data should be quantified. In this study we test the ability of six snow models to reproduce snow water equivalent, snow density and snow depth when they are forced by meteorological input data with gradually lower accuracy. The SNOWPACK, GEOTOP, HTESSEL, UTOPIA, SMASH and S3M snow models are forced, first, with high-quality measurements performed at the experimental site of Torgnon, located at 2160 m a.s.l. in the Italian Alps (control run). Then, the models are forced by data at gradually lower temporal and/or spatial resolution, obtained by (i) sampling the original Torgnon 30 min time series at 3, 6, and 12 h, (ii) spatially interpolating neighbouring in situ station measurements and (iii) extracting information from GLDAS, ERA5 and ERA-Interim reanalyses at the grid point closest to the Torgnon site. Since the selected models are characterized by different degrees of complexity, from highly sophisticated multi-layer snow models to simple, empirical, single-layer snow schemes, we also discuss the results of these experiments in relation to the model complexity. The results show that, when forced by accurate 30 min resolution weather station data, the single-layer, intermediate-complexity snow models HTESSEL and UTOPIA provide similar skills to the more sophisticated multi-layer model SNOWPACK, and these three models show better agreement with observations and more robust performances over different seasons compared to the lower-complexity models SMASH and S3M. All models forced by 3-hourly data provide similar skills to the control run, while the use of 6- and 12-hourly temporal resolution forcings may lead to a reduction in model performances if the incoming shortwave radiation is not properly represented. The SMASH model generally shows low sensitivity to the temporal degradation of the input data. Spatially interpolated data from neighbouring stations and reanalyses are found to be adequate forcings, provided that temperature and precipitation variables are not affected by large biases over the considered period. However, a simple bias-adjustment technique applied to ERA-Interim temperatures allowed all models to achieve similar performances to the control run. Regardless of their complexity, all models show weaknesses in the representation of the snow density. [ABSTRACT FROM AUTHOR]

The accuracy of hydrological predictions in snow-dominated regions deeply depends on the quality of the snowpack simulations, with dynamics that strongly affect the local hydrological regime, especially during the melting period. With the aim of reducing the modelling uncertainty, data assimilation techniques are increasingly being implemented for operational purposes. This study aims to investigate the performance of a multivariate sequential importance resampling - particle filter scheme, designed to jointly assimilate several ground-based snow observations. The system, which relies on a multilayer energy-balance snow model, has been tested at three Alpine sites: Col de Porte (France), Torgnon (Italy), andWeissfluhjoch (Switzerland). The implementation of a multivariate data assimilation scheme faces several challenging issues, which are here addressed and extensively discussed: (1) the effectiveness of the perturbation of the meteorological forcing data in preventing the sample impoverishment; (2) the impact of the parameter perturbation on the filter updating of the snowpack state; the system sensitivity to (3) the frequency of the assimilated observations, and (4) the ensemble size. The perturbation of the meteorological forcing data generally turns out to be insufficient for preventing the sample impoverishment of the particle sample, which is highly limited when jointly perturbating key model parameters. However, the parameter perturbation sharpens the system sensitivity to the frequency of the assimilated observations, which can be successfully relaxed by introducing indirectly estimated information on snow-mass-related variables. The ensemble size is found not to greatly impact the filter performance in this point-scale application. [ABSTRACT FROM AUTHOR]

A reliable estimation of soil moisture conditions is fundamental for rivers’ discharge predictions, especially in small catchments where flash floods occur. In this context, microwave remote sensing can be exploited to estimate soil moisture at large scale. These estimates can be used to enhance the predictions of hydrological models using data assimilation techniques. Flash flood early warning systems can, thus, be improved. This study tested the effect of the assimilation of three different ASCAT-derived soil moisture products, processed and distributed within the EUMETSAT H-SAF framework (SM-OBS-1, SM-OBS-2, SM-DAS-2), into a distributed physically based hydrological model (Continuum). The study areas were three Italian catchments, representative of the typical Mediterranean small basins prone to flash floods. The products were first preprocessed in order to be comparable with the model soil moisture state estimate. Subsequently, they were assimilated using three Nudging-based techniques. Then, observed discharges were compared with the modeled one in order to understand the impact of the assimilation. The analysis was executed for a multiyear period ranging from July 2012 to June 2014 in order to test the assimilation algorithms for operational purposes in real-cases scenarios. Findings showed that the assimilation of H-SAF soil moisture products with simple preprocessing and assimilation techniques can enhance discharge predictions; the improvements significantly affect high flows. Although SM-OBS-2 and SM-DAS-1 are added-value products with respect to SM-OBS-1 (respectively, higher spatial and temporal resolution), they may not necessarily perform better. The impact of the assimilation strongly relies on the permanent catchment characteristics (e.g., topography, hydrography, land cover). [ABSTRACT FROM PUBLISHER]

During the autumn of 2011 two catastrophic, very intense rainfall events affected two different parts of the Liguria Region of Italy causing various flash floods. The first occurred in October and the second at the beginning of November. Both the events were characterized by very high rainfall intensities (> 100mm h-1) that persisted on a small portion of territory causing local huge rainfall accumulations (> 400mm6 h-1). Two main considerations were made in order to set up this work. The first consideration is that various studies demonstrated that the two events had a similar genesis and similar triggering elements. The second very evident and coarse concern is that two main elements are needed to have a flash flood: a very intense and localized rainfall event and a catchment (or a group of catchments) to be affected. Starting from these assumptions we did the exercise of mixing the two flash flood ingredients by putting the rainfall field of the first event on the main catchment struck by the second event, which has its mouth in the biggest city of the Liguria Region: Genoa. A complete framework was set up to quantitatively carry out a "what if" experiment with the aim of evaluating the possible damages associated with this event. A probabilistic rainfall downscaling model was used to generate possible rainfall scenarios maintaining the main characteristics of the observed rainfall fields while a hydrological model transformed these rainfall scenarios in streamflow scenarios. A subset of streamflow scenarios is then used as input to a 2- D hydraulic model to estimate the hazard maps, and finally a proper methodology is applied for damage estimation. This leads to the estimation of the potential economic losses and of the risk level for the people that stay in the affected area. The results are interesting, surprising and in a way worrying: a rare but not impossible event (it occurred about 50 km away from Genoa) would have caused huge damages estimated between 120 and EUR 230 million for the affected part of the city of Genoa, Italy, and more than 17 000 potentially affected people. [ABSTRACT FROM AUTHOR]

We explore the performance of three types of stochastic models used for spatial rainfall downscaling and assess their ability to reproduce the statistics of precipitation fields observed during the GATE radar experiment. We consider a bounded multifractal cascade, an autoregressive linear process passed through a nonlinear static filter (sometimes called a meta-Gaussian model), and a model based on the presence of individual rainfall cells with power law profile. As test statistics we use the low-order moments of the amplitude distribution, the distribution of generalized fractal dimensions, the generalized scaling exponents, the slope of the power spectrum, and the properties of the spatial autocorrelation. The results of the analysis indicate that all models provide, on average, a satisfactory representation of the statistical properties of the GATE rainfall fields (including the anomalous scaling behavior), with a slightly better performance of the model based on individual rainfall cells. All models, however, display large scatter in the field-to-field comparison with the data. These results indicate that data analysis alone does not allow, at the moment, for preferring one downscaling approach over another. [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]

Since the 1980s, vegetated lands have experienced widespread greening at the global scale. Numerous studies have focused on spatial patterns and mechanisms of this phenomenon, especially in the Arctic and sub-Arctic regions. Greening trends in the European Alps have received less attention, although this region has experienced strong climate and land-use changes during recent decades. We studied the rates and spatial patterns of greening in an inner-alpine region of the Western Alps. We used MODIS-derived normalized difference vegetation index (NDVI) at 8-day temporal and 250 m spatial resolution, for the period 2000–2018, and removed areas with disturbances in order to consider the trends of undisturbed vegetation. The objectives of this study were to (i) quantify trends of greening in a representative area of the Western Alps; and (ii) examine mechanisms and causes of spatial patterns of greening across different plant types. We show that 63% of vegetated areas experienced significant trends during the 2000–2018 period, of which only 8% were negative. We identify (i) a climatic control on spring and autumn phenology with contrasting effects depending on plant type and elevation, and (ii) land-use change dynamics, such as shrub encroachment on abandoned pastures and colonization of new surfaces at high elevation. Below 1500 m, warming temperatures promote incremental greening in the transition from spring to summer, but not in fall, suggesting either photoperiod or water limitation. In the alpine and sub-alpine belts (>1800 m asl), snow prevents vegetation development until late spring, despite favorable temperatures. Instead, at high elevation greening acts both in summer and autumn. However, photoperiod limitation likely prevents forested ecosystems from fully exploiting warmer autumn conditions. We furthermore illustrate two emblematic cases of prominent greening: recent colonization of previously glaciated/non vegetated areas, as well as shrub/tree encroachment due to the abandonment of agricultural practices. Our results demonstrate the interplay of climate and land-use change in controlling greening dynamics in the Western Alps. [ABSTRACT FROM AUTHOR]

Information on snow properties is of critical relevance for a wide range of scientific studies and operational applications, mainly for hydrological purposes. However, the ground-based monitoring of snow dynamics is a challenging task, especially over complex topography and under harsh environmental conditions. Remote sensing is a powerful resource providing snow observations at a large scale. This study addresses the potential of using Sentinel-2 high-resolution imagery to assess moderate-resolution snow products, namely H10—Snow detection (SN-OBS-1) and H12—Effective snow cover (SN-OBS-3) supplied by the Satellite Application Facility on Support to Operational Hydrology and Water Management (H-SAF) project of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). With the aim of investigating the reliability of reference data, the consistency of Sentinel-2 observations is evaluated against both in-situ snow measurements and webcam digital imagery. The study area encompasses three different regions, located in Finland, the Italian Alps and Turkey, to comprehensively analyze the selected satellite products over both mountainous and flat areas having different snow seasonality. The results over the winter seasons 2016/17 and 2017/18 show a satisfying agreement between Sentinel-2 data and ground-based observations, both in terms of snow extent and fractional snow cover. H-SAF products prove to be consistent with the high-resolution imagery, especially over flat areas. Indeed, while vegetation only slightly affects the detection of snow cover, the complex topography more strongly impacts product performances. [ABSTRACT FROM AUTHOR]