Snow water equivalent modeling components in NewAge-JGrass

This paper presents a package of modified temperature-index-based snow water equivalent models as part of the hydrological modeling system NewAge-JGrass. Three temperature-based snow models are integrated into the NewAge-JGrass modeling system and use many of its components such as those for radiation balance (short wave radiation balance, SWRB), kriging (KRIGING), automatic calibration algorithms (particle swarm optimization) and tests of goodness of fit (NewAge-V), to build suitable modeling solutions (MS). Similarly to all the NewAge-JGrass components, the models can be executed both in raster and in vector mode. The simulation time step can be daily, hourly or sub-hourly, depending on user needs and availability of input data. The MS are applied on the Cache la Poudre River basin (CO, USA) using three test applications. First, daily snow water equivalent is simulated for three different measurement stations for two snow model formulations. Second, hourly snow water equivalent is simulated using all the three different snow model formulae. Finally, a raster mode application is performed to compute snow water equivalent maps for the whole Cache la Poudre Basin

The JGrass-NewAge system for forecasting and managing the hydrological budgets at the basin scale: models of flow generation and propagation/routing

Abstract. This paper presents a discussion of the predictive capacity of the implementation of the semi-distributed hydrological modeling system JGrass-NewAge. This model focuses on the hydrological budgets of medium scale to large scale basins as the product of the processes at the hillslope scale with the interplay of the river network. The part of the modeling system presented here deals with the: (i) estimation of the space-time structure of precipitation, (ii) estimation of runoff production; (iii) aggregation and propagation of flows in channel; (v) estimation of evapotranspiration; (vi) automatic calibration of the discharge with the method of particle swarming. The system is based on a hillslope-link geometrical partition of the landscape, combining raster and vectorial treatment of hillslope data with vector based tracking of flow in channels. Measured precipitation are spatially interpolated with the use of kriging. Runoff production at each channel link is estimated through a peculiar application of the Hymod model. Routing in channels uses an integrated flow equation and produces discharges at any link end, for any link in the river network. Evapotranspiration is estimated with an implementation of the Priestley-Taylor equation. The model system assembly is calibrated using the particle swarming algorithm. A two year simulation of hourly discharge of the Little Washita (OK, USA) basin is presented and discussed with the support of some classical indices of goodness of fit, and analysis of the residuals. A novelty with respect to traditional hydrological modeling is that each of the elements above, including the preprocessing and the analysis tools, is implemented as a software component, built upon Object Modelling System v3 and jgrasstools prescriptions, that can be cleanly switched in and out at run-time, rather than at compiling time. The possibility of creating different modeling products by the connection of modules with or without the calibration tool, as for instance the case of the present modeling chain, reduces redundancy in programming, promotes collaborative work, enhances the productivity of researchers, and facilitates the search for the optimal modeling solution

Hydrological control of soil thickness spatial variability on the initiation of rainfall-induced shallow landslides using a three-dimensional model

Thickness and stratigraphic settings of soils covering slopes potentially control susceptibility to initiation of rainfall-induced shallow landslides due to their local effect on slope hydrological response. Notwithstanding the relevance of the assessment of hazard to shallow landsliding at a distributed scale by approaches based on a coupled modelling of slope hydrological response and slope stability, the spatial variability of soil thickness and stratigraphic settings are factors poorly considered in the literature. Under these premises, this paper advances the well-known case study of rainfall-induced shallow landslides involving ash-fall pyroclastic soils covering the peri-Vesuvian mountains (Campania, southern Italy). In such a unique geomorphological setting, the soil covering is formed by alternating loose ash-fall pyroclastic deposits and paleosols, with high contrasts in hydraulic conductivity and total thickness decreasing as the slope angle increases, thus leading to the establishment of lateral flow and an increase of pore water pressure in localised sectors of the slope where soil horizon thickness is less. In particular, we investigate the effects, on hillslope hydrological regime and slope stability, of irregular bedrock topography, spatial variability of soil thickness and vertical hydraulic heterogeneity of soil horizons, by using a coupled three-dimensional hydrological and a probabilistic infinite slope stability model. The modelling is applied on a sample mountain catchment, located on Sarno Mountains (Campania, southern Italy), and calibrated using physics-based rainfall thresholds derived from the literature. The results obtained under five simulated constant rainfall intensities (2.5, 5, 10, 20 and 40 mm h−1) show an increase of soil pressure head and major failure probability corresponding to stratigraphic and morphological discontinuities, where a soil thickness reduction occurs. The outcomes obtained from modelling match the hypothesis of the formation of lateral throughflow due to the effect of intense rainfall, which leads to the increase of soil water pressure head and water content, up to values of near-saturation, in narrow zones of the slope, such as those of downslope reduction of total soil thickness and pinching out of soil horizons. The approach proposed can be conceived as a further advance in the comprehension of slope hydrological processes at a detailed scale and their effects on slope stability under given rainfall and antecedent soil hydrological conditions, therefore in predicting the most susceptible areas to initiation of rainfall-induced shallow landslides and the related I-D rainfall thresholds. Results obtained demonstrate the occurrence of a slope hydrological response depending on the spatial variability of soil thickness and leading to focus slope instability in specific slope sectors. The approach proposed is conceived to be potentially exportable to other slope environments for which a spatial modelling of soil thickness would be possible

Trends in heat and cold wave risks for the Italian Trentino-Alto Adige region from 1980 to 2018

Heat waves (HWs) and cold waves (CWs) can have considerable impact on people. Mapping risks of extreme temperature at local scale, accounting for the interactions between hazard, exposure, and vulnerability, remains a challenging task. In this study, we quantify risks from HWs and CWs for the Trentino-Alto Adige region of Italy from 1980 to 2018 at high spatial resolution. We use the Heat Wave Magnitude Index daily (HWMId) and the Cold Wave Magnitude Index daily (CWMId) as the hazard indicators. To obtain HWs and CW risk maps we combined the following: (i) occurrence probability maps of the hazard obtained using the zero-inflated Tweedie distribution (accounting directly for the absence of events for certain years), (ii) normalized population density maps, and (iii) normalized vulnerability maps based on eight socioeconomic indicators. The methodology allowed us to disentangle the contributions of each component of the risk relative to total change in risk. We find a statistically significant increase in HW hazard and exposure, while CW hazard remained stagnant in the analyzed area over the study period. A decrease in vulnerability to extreme temperature spells is observed through the region except in the larger cities where vulnerability increased. HW risk increased in 40 % of the region, with the increase being greatest in highly populated areas. Stagnant CW hazard and declining vulnerability result in reduced CW risk levels overall, except for the four main cities where increased vulnerability and exposure increased risk levels. These findings can help to steer investments in local risk mitigation, and this method can potentially be applied to other regions where there are sufficient detailed data.</p

Modeling the water budget of the Upper Blue Nile basin using the JGrass-NewAge model system and satellite data

The Upper Blue Nile basin is one of the most data-scarce regions in developing countries, and hence the hydrological information required for informed decision making in water resource management is limited. The hydrological complexity of the basin, tied with the lack of hydrometeorological data, means that most hydrological studies in the region are either restricted to small subbasins where there are relatively better hydrometeorological data available, or on the whole-basin scale but at very coarse timescales and spatial resolutions. In this study we develop a methodology that can improve the state of the art by using available, but sparse, hydrometeorological data and satellite products to obtain the estimates of all the components of the hydrological cycle (precipitation, evapotranspiration, discharge, and storage). To obtain the water-budget closure, we use the JGrass-NewAge system and various remote sensing products. The satellite product SM2R-CCI is used for obtaining the rainfall inputs, SAF EUMETSAT for cloud cover fraction for proper net radiation estimation, GLEAM for comparison with NewAge-estimated evapotranspiration, and GRACE gravimetry data for comparison of the total water storage amounts available in the whole basin. Results are obtained at daily time steps for the period 1994–2009 (16 years), and they can be used as a reference for any water resource development activities in the region. The overall water-budget analysis shows that precipitation of the basin is 1360 ± 230 mm per year. Evapotranspiration accounts for 56 % of the annual water budget, runoff is 33 %, storage varies from −10 to +17 % of the water budget

Future hot-spots for hydro-hazards in Great Britain: a probabilistic assessment

In an increasing hydro-climatic risk context as a result of climate change, this work aims to identify future hydro-hazard hot-spots as a result of climate change across Great Britain. First, flood and drought hazards were defined and selected in a consistent and parallel approach with a threshold method. Then, a nation-wide systematic and robust statistical framework was developed to quantify changes in frequency, magnitude, and duration, and assess time of year for both droughts and floods, and the uncertainty associated with climate model projections. This approach was applied to a spatially coherent statistical database of daily river flows (Future Flows Hydrology) across Great Britain to assess changes between the baseline (1961–1990) and the 2080s (2069–2098). The results showed that hydro-hazard hot-spots are likely to develop along the western coast of England and Wales and across north-eastern Scotland, mainly during the winter (floods) and autumn (droughts) seasons, with a higher increase in drought hazard in terms of magnitude and duration. These results suggest a need for adapting water management policies in light of climate change impact, not only on the magnitude, but also on the timing of hydro-hazard events, and future policy should account for both extremes together, alongside their potential future evolution.</p

Evaluating performance of simplified physically based models for shallow landslide susceptibility

Rainfall-induced shallow landslides can lead to loss of life and significant damage to private and public properties, transportation systems, etc. Predicting locations that might be susceptible to shallow landslides is a complex task and involves many disciplines: hydrology, geotechnical science, geology, hydrogeology, geomorphology, and statistics. Two main approaches are commonly used: statistical or physically based models. Reliable model applications involve automatic parameter calibration, objective quantification of the quality of susceptibility maps, and model sensitivity analyses. This paper presents a methodology to systemically and objectively calibrate, verify, and compare different models and model performance indicators in order to identify and select the models whose behavior is the most reliable for particular case studies.The procedure was implemented in a package of models for landslide susceptibility analysis and integrated in the NewAge-JGrass hydrological model. The package includes three simplified physically based models for landslide susceptibility analysis (M1, M2, and M3) and a component for model verification. It computes eight goodness-of-fit indices by comparing pixel-by-pixel model results and measurement data. The integration of the package in NewAge-JGrass uses other components, such as geographic information system tools, to manage input–output processes, and automatic calibration algorithms to estimate model parameters. The system was applied for a case study in Calabria (Italy) along the Salerno–Reggio Calabria highway, between Cosenza and Altilia. The area is extensively subject to rainfall-induced shallow landslides mainly because of its complex geology and climatology. The analysis was carried out considering all the combinations of the eight optimized indices and the three models. Parameter calibration, verification, and model performance assessment were performed by a comparison with a detailed landslide inventory map for the area. The results showed that the index distance to perfect classification in the receiver operating characteristic plane (D2PC) coupled with the model M3 is the best modeling solution for our test case

Estimating the index flood with continuous hydrological models: an application in Great Britain

Estimating peak river discharge, a critical issue in engineering hydrology, is essential for designing and managing hydraulic infrastructure such as dams and bridges. In the UK, practitioners typically apply the Flood Estimation Handbook (FEH) statistical method which estimates the design flood as the product of a relatively frequent flow estimate (the index flood, IF) and a regional growth factor. For gauged catchments the IF is estimated from observations. For ungauged catchments it is computed through a multiple regression model. While the FEH IF method provides peak flow estimates that are statistically robust, it does not readily take into account catchment heterogeneity or effect of environmental change on river flows. This study presents a new methodology to estimate the IF at national scale using continuous simulation from a physically based hydrological model (Grid-to-Grid). The methodology is tested across Great Britain and compares well with IF estimates at 550 gauging stations (R2 = 0.91). The promising results for Great Britain support the aspiration that continuous simulation from large-scale hydrological models coupled with increasing availability of global weather and climate products, could be used to estimate design floods in regions with limited gauge data or affected by environmental change

Differential orographic impact on sub-hourly, hourly, and daily extreme precipitation

Extreme precipitation of multiple durations is responsible for major natural hazards in mountainous regions, such as flash floods and debris flows. Understanding the orographic impact on the statistics of precipitation extremes is therefore crucial for improving hydrological design and risk management strategies. Here, we use a novel statistical approach for the analysis of extremes based on ordinary events to improve our understanding of the orographic impact on extreme precipitation of durations ranging between 5 min and 24 h. We focus on Trentino, a rough orographic region in the eastern Italia Alps, and use data from 78 quality-controlled rain gauges with 5-minute resolution. We show that our framework well reproduces the statistical properties of the observed annual maxima (Nash-Sutcliffe 0.82–0.95, Bias from -4% to 7%) as well as their relation with orography. We then exploit the reduced uncertainty of this approach to quantify the orographic impact on precipitation right-tail statistics and on extreme return levels using a regression analysis. We identify two main modes of orographic relationship: a reverse orographic effect for hourly and sub-hourly durations (10–20% decrease per 1000 m elevation) and an orographic enhancement for durations of ∼8 h or longer (7.5–10% increase per 1000 m elevation). We observe that these two modes result from three main precipitation regimes, which show different proportion between extreme and very-extreme events and which emerge at very short durations (∼20 min or shorter), mid durations (∼30 min-1 hour) and long durations (∼2 h or longer). These findings are of interest for risk management applications and climate change impact studies

Modeling shortwave solar radiation using the JGrass-NewAge system

This paper presents two new modeling components based on the object modeling system v3 (OMS3) for the calculation of the shortwave incident radiation (Rsw&darr;) on complex topography settings, and the implementation of several ancillary tools. The first component, NewAGE-SwRB, accounts for elevation slope, aspect, shadow of the sites, and uses suitable parameterization for obtaining the cloudless irradiance. A second component, NewAGE-DEC-MOD's is implemented to estimate the irradiance reduction due to the presence of clouds according to three parameterizations. To obtain a working modeling composition that is comparable with ground data at measurement stations the two components are connected to a kriging component. With the help of an additional component, NewAGE-V (verification package), the performance of modeled (Rsw&darr;) is quantitatively evaluated. The two components (and the various parameterizations they contain) are tested using the data from three basins, and some simple verification tests were carried out to assess the goodness of the methods used. Moreover, a raster mode test is performed in order to show the capability of the system in providing solar radiation raster maps. The components are part of a larger system, JGrass-NewAGE, their input and outputs are geometrical objects immediately displayed in a geographical information system (GIS). They can be used seamlessly with the various modeling solutions available in JGrass-NewAGE for the estimation of long wave radiation, evapotranspiration, and snow melting, as well as standalone components to just estimate shortwave radiation for various uses. The modularity of the approach leads to more accurate physical-statistical studies aimed to assess in depth the components' performances and extends their results spatially, without the necessity of recoding any part of the component