Search

In mountainous areas, accurately estimating the long‐term climatology of seasonal precipitations is challenging due to the lack of high‐altitude rain gauges and the complexity of the topography. This study addresses these challenges by interpolating seasonal precipitation data from 3189 rain gauges across France over the 1982–2018 period, using geographical coordinates, and altitude. In this study, an additional predictor is provided from simulations of a Convection‐Permitting Regional Climate Model (CP‐RCM). The simulations are averaged to obtain seasonal precipitation climatology, which helps capture the relationship between topography and long‐term seasonal precipitation. Geostatistical and machine learning models are evaluated within a cross‐validation framework to determine the most appropriate approach to generate seasonal precipitation reference fields. Results indicate that the best model uses a machine learning approach to interpolate the ratio between long‐term seasonal precipitation from observations and CP‐RCM simulations. This method successfully reproduces both the mean and variance of observed data, and slightly outperforms the best geostatistical model. Moreover, incorporating the CP‐RCM outputs as an explanatory variable significantly improves interpolation accuracy and altitude extrapolation, especially when the rain gauge density is low. These results imply that the commonly used altitude‐precipitation relationship may be insufficient to derive seasonal precipitation fields. The CP‐RCM simulations, increasingly available worldwide, present an opportunity for improving precipitation interpolation, especially in sparse and complex topographical regions. [ABSTRACT FROM AUTHOR]

Seasonal precipitation estimation in ungauged mountainous areas is essential for understanding and modeling a physical variable of interest in many environmental applications (hydrology, ecology, and cryospheric studies). Precipitation lapse rates (PLRs), defined as the increasing or decreasing rate of precipitation amounts with the elevation, play a decisive role in high-altitude precipitation estimation. However, the documentation of PLR in mountainous regions remains weak even though their utilization in environmental applications is frequent. This article intends to assess the spatial variability and the spatial-scale dependence of seasonal PLRs in a varied and complex topographical region. At the regional scale (10 000 km 2), seven different precipitation products are compared in their ability to reproduce the altitude dependence of the annual/seasonal precipitation of 1836 stations located in France. The convection-permitting regional climate model (CP-RCM) AROME is the best in this regard, despite severe precipitation overestimation in high altitudes. The fine resolution of AROME allows for a precise assessment of the influence of altitude on winter and summer precipitation on 23 massifs at the sub-regional scale (∼ 1000 km 2) and 2748 small catchments (∼ 100 km 2) through linear regressions. With AROME, PLRs are often higher in winter at the catchment scale. The variability in the PLR is higher in high-altitude regions such as the French Alps, with higher PLRs at the border than inside the massifs. This study emphasizes the interest of conducting a PLR investigation at a fine scale to reduce spatial heterogeneity in the seasonal precipitation–altitude relationships. [ABSTRACT FROM AUTHOR]

Intensity‐Duration‐Area‐Frequency (IDAF) models provide the mathematical link between precipitation intensities (I), durations (D), areas (A), and frequency of occurrence (F). They play a critical role in hydrological design, areal rainfall hazard quantification, storm characterization, and early warning system development. IDAF models extend the conventional Intensity‐Duration‐Frequency models by accounting for the spatial extent of precipitation (i.e., the area). In this study, we develop IDAF models using the entire non‐zero precipitation intensities, not only the extremes. We use the extended generalized Pareto distribution (EGPD) to model the precipitation intensities. To build the IDAF models, we adopt a data‐driven approach that allows the linkage of EGPD parameters with duration and area, based on empirically determined parametric relationships. The inference of model parameters is done using a global maximum likelihood estimation, and uncertainties are assessed by the bootstrap method. The study area is Switzerland, a topographically complex region of 42,000 km2 with regional precipitation variability and clear seasonality. The study utilizes 17 years of data from CombiPrecip, a radar‐reanalysis product developed by geostatistically merging radar and rain gauge data in an operational setting. We build the IDAF models for the spatiotemporal range of 1–72 hr and 1 to 1,089 km2 at each pixel in the study area. To the best of our knowledge, our study is the first attempt to use the EGPD in IDAF curve modeling. It discusses the use and limitations of CombiPrecip in extreme value analysis and highlights the challenges of modeling areal precipitation in a complex topographical environment. Key Points: We develop seasonal Intensity Duration Area Frequency models at every pixel location in SwitzerlandWe use all the non‐zero precipitation data and model the intensities using the extended generalized Pareto distributionWe highlight the complexity of modeling areal precipitation in mountainous regions [ABSTRACT FROM AUTHOR]

Seasonal precipitation estimation in ungauged mountainous areas is essential for understanding and modeling a physical variable of interest in many environmental applications (hydrology, ecology, cryospheric studies). Precipitation Lapse Rates (PLRs), defined as the increasing or decreasing rate of precipitation amounts with the elevation, play a decisive role in high-altitude precipitation estimation. However, the documentation of PLR in mountainous regions remains weak even though their utilization in environmental applications is frequent. This article intends to assess the spatial variability and the spatial-scale dependence of seasonal PLRs in a varied and complex topography region. At the regional scale (10,000 km2), seven different precipitation products are compared in their ability to reproduce the altitude dependence of the annual/seasonal precipitations of 1,836 stations located in France. The Convection-Permitting Regional Climate Model (CP-RCM) AROME is the best in this regard, despite severe precipitation overestimation in high altitudes. The fine resolution of AROME allows for a precise assessment of the influence of altitude on winter and summer precipitations on 23 massifs at the sub-regional scale (∼ 1,000 km2) and 2,748 small catchments (∼ 100 km2) through linear regressions. With AROME, PLRs are mostly higher in winter at the catchment scale. The variability of PLR is higher in high-altitude regions such as the French Alps, with higher PLRs at the border than inside the massifs. This study emphasizes the interest of conducting PLR investigation at a fine scale to reduce spatial heterogeneity in the seasonal precipitations–altitude relationships. [ABSTRACT FROM AUTHOR]

Intensity‐duration‐frequency (IDF) curves are useful in water resources engineering for the planning and design of hydrological structures. As opposed to the common use of only extreme data to build IDF curves, here, we use all the non‐zero rainfall intensities, thereby making efficient use of the available information. We use the extended generalized Pareto distribution to model the distribution of the non‐zero rainfall intensities. We consider three commonly used approaches for building IDF curves. The first approach is based on the scale‐invariance property of rainfall, the second relies on the general IDF formulation of Koutsoyiannis et al. (1998, https://doi.org/10.1016/S0022-1694(98)00097-3), and the last approach is purely data‐driven(Overeem et al., 2008, https://doi.org/10.1016/j.jhydrol.2007.09.044). Using these three approaches, and some extensions around them, we build a total of 10 models for the IDF curves, and then we compare them in a split‐sampling cross‐validation framework. We consider a total of 81 stations at 10 min resolution in Switzerland. Due to the marked seasonality of rainfall in the study area, we performed a seasonal‐based analysis. The results reveal the model based on the data‐driven approach as the best model. It is able to correctly model the observed intensities across duration while being reliable and robust. It is also able to reproduce the space and time variability of extreme rainfall across Switzerland. Plain Language Summary: Intensity‐duration‐frequency (IDF) curves are used in water resources engineering for the safe design of hydraulic structures. They summarize the relationship between rainfall intensity, duration, and frequency of occurrence. IDF curves are commonly built using probability distributions that use only extreme data, which can be a drawback in terms of the efficient use of already scarce data. In this study, we use a distribution, the extended generalized Pareto distribution, that uses all the non‐zero rainfall intensities to model IDF curves. We build a total of 10 models, which are based on three commonly used approaches and some extensions around them. We compare them in terms of their reliability and robustness. The best IDF model is able to reproduce the known variability of extreme rainfall across Switzerland. Key Points: We build intensity duration frequency (IDF) curves using all the non‐zero precipitation data and model the intensities with the extended generalized Pareto distributionWe consider three approaches to building IDF curves: scale invariance, IDF formulation of Koutsoyiannis et al. (1998), https://doi.org/10.1016/S0022-1694(98)00097-3 and data‐driven methodWe compare and select the best model in a cross‐validation framework [ABSTRACT FROM AUTHOR]

In this article, we compare the performance of three regionalization approaches in improving the at-site estimates of daily precipitation. The first method is built on the idea of conventional RFA (regional frequency analysis) but is based on a fast algorithm that defines distinct homogeneous regions relying on their upper-tail similarity. It uses only the precipitation data at hand without the need for any additional covariate. The second is based on the region-of-influence (ROI) approach in which neighborhoods, containing similar sites, are defined for each station. The third is a spatial method that adopts generalized additive model (GAM) forms for the model parameters. In line with our goal of modeling the whole range of positive precipitation, the chosen marginal distribution model is the extended generalized Pareto distribution (EGPD) to which we apply the three methods. We consider a dense network composed of 1176 daily stations located within Switzerland and in neighboring countries. We compute different criteria to assess the models' performance in the bulk of the distribution and the upper tail. The results show that all the regional methods offered improved robustness over the local EGPD model. While the GAM method is more robust and reliable in the upper tail, the ROI method is better in the bulk of the distribution. [ABSTRACT FROM AUTHOR]

Estimates for rare to very rare floods are limited by the relatively short streamflow records available. Often, pragmatic conversion factors are used to quantify such events based on extrapolated observations, or simplifying assumptions are made about extreme precipitation and resulting flood peaks. Continuous simulation (CS) is an alternative approach that better links flood estimation with physical processes and avoids assumptions about antecedent conditions. However, long-term CS has hardly been implemented to estimate rare floods (i.e. return periods considerably larger than 100 years) at multiple sites in a large river basin to date. Here we explore the feasibility and reliability of the CS approach for 19 sites in the Aare River basin in Switzerland (area: 17 700 km 2) with exceedingly long simulations in a hydrometeorological model chain. The chain starts with a multi-site stochastic weather generator used to generate 30 realizations of hourly precipitation and temperature scenarios of 10 000 years each. These realizations were then run through a bucket-type hydrological model for 80 sub-catchments and finally routed downstream with a simplified representation of main river channels, major lakes and relevant floodplains in a hydrologic routing system. Comprehensive evaluation over different temporal and spatial scales showed that the main features of the meteorological and hydrological observations are well represented and that meaningful information on low-probability floods can be inferred. Although uncertainties are still considerable, the explicit consideration of important processes of flood generation and routing (snow accumulation, snowmelt, soil moisture storage, bank overflow, lake and floodplain retention) is a substantial advantage. The approach allows for comprehensively exploring possible but unobserved spatial and temporal patterns of hydrometeorological behaviour. This is of particular value in a large river basin where the complex interaction of flows from individual tributaries and lake regulations are typically not well represented in the streamflow observations. The framework is also suitable for estimating more frequent floods, as often required in engineering and hazard mapping. [ABSTRACT FROM AUTHOR]

Resume : En Amazonie, les petits bassins versants ne disposent pas de donnees de debit, ce qui empeche l'etude de ces bassins, par exemple, pour des fins de production d'energie [...]

In this article, we compare the performances of three regionalization approaches in improving the at-site estimates of daily precipitation. The first method is built on the idea of conventional RFA (Regional Frequency Analysis) but is based on a fast algorithm that defines distinct homogeneous regions relying on their upper tail similarity. It uses only the precipitation data at hand without the need for any additional covariate. The second is based on the region-of-influence (ROI) approach in which neighborhoods, containing similar sites, are defined for each station. The third is a spatial method that adopts Generalized Additive Model (GAM) forms for the model parameters. In line with our goal of modeling the whole range of positive precipitation, the chosen marginal distribution model is the Extended Generalized Pareto Distribution (EGPD) on which we apply the three methods. We consider a dense network composed of 1176 daily stations located within Switzerland and in neighboring countries. We compute different criteria to assess the models' performances both in the bulk of the distribution as well as in the upper tail. The results show that all the regional methods offered improved robustness over the local EGPD model. While the GAM method is more robust and reliable in the upper tail, the ROI method is better in the bulk of the distribution. [ABSTRACT FROM AUTHOR]

Natural risk studies such as flood risk assessments require long series of weather variables. As an alternative to observed series, which have a limited length, these data can be provided by weather generators. Among the large variety of existing ones, resampling methods based on analogues have the advantage of guaranteeing the physical consistency between local weather variables at each time step. However, they cannot generate values of predictands exceeding the range of observed values. Moreover, the length of the simulated series is typically limited to the length of the synoptic meteorological records used to characterize the large-scale atmospheric configuration of the generation day. To overcome these limitations, the stochastic weather generator proposed in this study combines two sampling approaches based on atmospheric analogues: (1) a synoptic weather generator in a first step, which recombines days of the 20th century to generate a 1000-year sequence of new atmospheric trajectories, and (2) a stochastic downscaling model in a second step applied to these atmospheric trajectories, in order to simulate long time series of daily regional precipitation and temperature. The method is applied to daily time series of mean areal precipitation and temperature in Switzerland. It is shown that the climatological characteristics of observed precipitation and temperature are adequately reproduced. It also improves the reproduction of extreme precipitation values, overcoming previous limitations of standard analogue-based weather generators. [ABSTRACT FROM AUTHOR]

Over the last years, several stationarity tests have been proposed. One of these methods uses time-frequency representations and stationarized replicas of the signal (known as surrogates) for testing wide-sense stationarity. In this letter, we propose a procedure to improve the original surrogate test. The proposed methodology can be seen as a guideline on how the surrogate test can be improved. We show mathematically that the modified test should exhibit improved classification performance. Numerical simulations on synthetic and real-world signals are carried out to evaluate the modified test against competing ones. [ABSTRACT FROM AUTHOR]

Reliable flood estimates are needed for designing safe and cost‐effective flood protection structures. Classical flood estimation methods applied for deriving such estimates focus on peak discharge and neglect other important flood characteristics such as flood volume and the interdependence among different flood characteristics. Furthermore, they do not account for potential nonstationarities in hydrological time series due to climate change. The consideration of both the interdependence between peak discharge and flood volume and its nonstationarity might help us derive more reliable flood estimates. A few studies have looked at changes in the general dependence of peak discharge and flood volume for small sets of catchments and explored ways of modeling such changes. However, spatial analyses of trends in this dependence or in their climatological drivers have not been carried out. The aim of this study was to help close this knowledge gap by first quantifying trends in the general dependence between peak discharge and flood volume as described by Kendall's tau on a spatially comprehensive data set of 307 catchments in Switzerland. Second, potential climatological drivers for changes in the dependence between peak discharge and flood volume were identified. Our results show that the dependence between peak discharge and flood volume and its trends are spatially heterogeneous. This pattern cannot be explained by one driver only but by an interplay of changes in precipitation, snowmelt, and soil moisture. Both the trends and the links between drivers and trends depend on the climate model chain considered and are therefore uncertain. Key Points: The dependence of peak discharge and flood volume and its trends are spatially heterogeneousSpatial differences in trends in the peak‐volume dependence are partly explained by differences in hydro‐climatic drivers, i.e. precipitation, snowmelt, and soil moistureTrends in the peak‐volume dependence are dependent on the model chain and therefore uncertain [ABSTRACT FROM AUTHOR]

Floods often affect not only a single location, but also a whole region. Flood frequency analysis should therefore be undertaken at a regional scale which requires the considerations of the dependence of events at different locations. This dependence is often neglected even though its consideration is essential to derive reliable flood estimates. A model used in regional multivariate frequency analysis should ideally consider the dependence of events at multiple sites which might show dependence in the lower and/or upper tail of the distribution. We here seek to propose a simple model that on the one hand considers this dependence with respect to the network structure of the region and on the other hand allows for the simulation of stochastic event sets at both gauged and ungauged locations. The new Fisher copula model is used for representing the spatial dependence of flood events in the nested Thur catchment in Switzerland. Flood event samples generated for the gauged stations using the Fisher copula are compared to samples generated by other dependence models allowing for modeling of multivariate data including elliptical copulas, R-vine copulas, and max-stable models. The comparison of the dependence structures of the generated samples shows that the Fisher copula is a suitable model for capturing the spatial dependence in the data. We therefore use the copula in a way such that it can be used in an interpolation context to simulate event sets comprising gauged and ungauged locations. The spatial event sets generated using the Fisher copula well capture the general dependence structure in the data and the upper tail dependence, which is of particular interest when looking at extreme flood events and when extrapolating to higher return periods. The Fisher copula was for a medium-sized catchment found to be a suitable model for the stochastic simulation of flood event sets at multiple gauged and ungauged locations. [ABSTRACT FROM AUTHOR]

Design flood estimates for a given return period are required in both gauged and ungauged catchments for hydraulic design and risk assessments. Contrary to classical design estimates, synthetic design hydrographs provide not only information on the peak magnitude of events but also on the corresponding hydrograph volumes together with the hydrograph shapes. In this study, we tested different regionalization approaches to transfer parameters of synthetic design hydrographs from gauged to ungauged catchments. These approaches include classical regionalization methods such as linear regression techniques, spatial methods, and methods based on the formation of homogeneous regions. In addition to these classical approaches, we tested nonlinear regression models not commonly used in hydrological regionalization studies, such as random forest, bagging, and boosting. We found that parameters related to the magnitude of the design event can be regionalized well using both linear and nonlinear regression techniques using catchment area, length of the main channel, maximum precipitation intensity, and relief energy as explanatory variables. The hydrograph shape, however, was found to be more difficult to regionalize due to its high variability within a catchment. Such variability might be better represented by looking at flood-type specific synthetic design hydrographs. [ABSTRACT FROM AUTHOR]

Reanalyses, based on numerical weather prediction methods assimilating past observations, provide continuous precipitation datasets and represent interesting options for assessing the climatology of regions with sparse station networks (e.g., northern Canada). However, reanalysis series cannot be used directly because of possible biases and mismatch between their spatial and temporal resolutions with that needed for local applications. To address these issues, a Stochastic Model Output Statistics (SMOS) approach was selected to post-process precipitation series simulated by the Climate Forecast System Reanalysis (CFSR) across Canada. This approach uses CFSR precipitation as a covariate and is based on two regression models: the first one is a logistic regression that deals with precipitation occurrence, and the second is a vector generalized linear model for precipitation intensity. At-site post-processed daily precipitation series are randomly generated using the SMOS approach, and selected climate indicators from the Expert Team on Climate Change Detection and Indices, which is jointly sponsored by the Commission for Climatology of the World Meteorological Organization's (WMO) World Climate Data and Monitoring Programme, the Climate Variability and Predictability Programme of the World Climate Research Programme, and the Joint WMO-IOC Technical Commission for Oceanography and Marine Meteorology (CCI/CLIVAR/JCOMM) are estimated and compared with corresponding observed and CFSR values. The two models in the SMOS approach, in addition to adequately correcting systematic biases, produced better predictions than the climatology of the wet and dry and intensity sequences. Additionally, the SMOS generally yields consistent climate indices when compared with those from CFSR without post-processing, though there is still room for improvement for specific indices (e.g., annual maximum of cumulative wet days). [ABSTRACT FROM AUTHOR]

Floods do often not only affect a single location but a whole region. Flood frequency analysis should therefore be undertaken at a regional scale which requires the considerations of the dependence of events at different locations. This dependence is often neglected even though its consideration is essential to derive reliable flood estimates. A model used in regional multivariate frequency analysis should ideally consider the dependence of events at multiple sites which might show dependence in the lower and/or upper tail of the distribution. We here seek at proposing a simple model that on the one hand considers this dependence with respect to the network structure of the region and on the other hand, allows for the simulation of stochastic event sets at both gauged and ungauged locations. The new Fisher copula model is used for representing the spatial dependence of flood events in the nested Thur catchment in Switzerland. Flood event samples generated for the gauged stations using the Fisher copula are compared to samples generated by other dependence models allowing for modeling multivariate data including elliptical copulas, R-vine copulas, and max-stable models. The comparison of the dependence structures of the generated samples shows that the Fisher copula is a suitable model for capturing the spatial dependence in the data. We therefore use the copula in a way such that it can be used in an interpolation context to simulate event sets comprising gauged and ungauged locations. The spatial event sets generated using the Fisher copula well capture the general dependence structure in the data and the upper tail dependence, which is of particular interest when looking at extreme flood events and when extrapolating to higher return periods. The Fisher copula is therefore a suitable model for the stochastic simulation of flood event sets at multiple gauged and ungauged locations. [ABSTRACT FROM AUTHOR]

Abstract: Flood hydrograph shapes contain valuable information on the flood‐generation mechanisms of a catchment. To make good use of this information, we express flood hydrograph shapes as continuous functions using a functional data approach. We propose a clustering approach based on functional data for flood hydrograph shapes to identify a set of representative hydrograph shapes on a catchment scale and use these catchment‐specific sets of representative hydrographs to establish regions of catchments with similar flood reactivity on a regional scale. We applied this approach to flood samples of 163 medium‐size Swiss catchments. The results indicate that three representative hydrograph shapes sufficiently describe the hydrograph shape variability within a catchment and therefore can be used as a proxy for the flood behavior of a catchment. These catchment‐specific sets of three hydrographs were used to group the catchments into three reactivity regions of similar flood behavior. These regions were not only characterized by similar hydrograph shapes and reactivity but also by event magnitudes and triggering event conditions. We envision these regions to be useful in regionalization studies, regional flood frequency analyses, and to allow for the construction of synthetic design hydrographs in ungauged catchments. The clustering approach based on functional data which establish these regions is very flexible and has the potential to be extended to other geographical regions or toward the use in climate impact studies. [ABSTRACT FROM AUTHOR]

Abstract: Design hydrographs described by peak discharge, hydrograph volume, and hydrograph shape are essential for engineering tasks involving storage. Such design hydrographs are inherently uncertain as are classical flood estimates focusing on peak discharge only. Various sources of uncertainty contribute to the total uncertainty of synthetic design hydrographs for gauged and ungauged catchments. These comprise model uncertainties, sampling uncertainty, and uncertainty due to the choice of a regionalization method. A quantification of the uncertainties associated with flood estimates is essential for reliable decision making and allows for the identification of important uncertainty sources. We therefore propose an uncertainty assessment framework for the quantification of the uncertainty associated with synthetic design hydrographs. The framework is based on bootstrap simulations and consists of three levels of complexity. On the first level, we assess the uncertainty due to individual uncertainty sources. On the second level, we quantify the total uncertainty of design hydrographs for gauged catchments and the total uncertainty of regionalizing them to ungauged catchments but independently from the construction uncertainty. On the third level, we assess the coupled uncertainty of synthetic design hydrographs in ungauged catchments, jointly considering construction and regionalization uncertainty. We find that the most important sources of uncertainty in design hydrograph construction are the record length and the choice of the flood sampling strategy. The total uncertainty of design hydrographs in ungauged catchments depends on the catchment properties and is not negligible in our case. [ABSTRACT FROM AUTHOR]

Statistical downscaling models (SDMs) are often used to produce local weather scenarios from large-scale atmospheric information. SDMs include transfer functions which are based on a statistical link identified from observations between local weather and a set of large-scale predictors. As physical processes driving surface weather vary in time, the most relevant predictors and the regression link are likely to vary in time too. This is well known for precipitation for instance and the link is thus often estimated after some seasonal stratification of the data. In this study, we present a two-stage analog/regression model where the regression link is estimated from atmospheric analogs of the current prediction day. Atmospheric analogs are identified from fields of geopotential heights at 1000 and 500 hPa. For the regression stage, two generalized linear models are further used to model the probability of precipitation occurrence and the distribution of non-zero precipitation amounts, respectively. The two-stage model is evaluated for the probabilistic prediction of small-scale precipitation over France. It noticeably improves the skill of the prediction for both precipitation occurrence and amount. As the analog days vary from one prediction day to another, the atmospheric predictors selected in the regression stage and the value of the corresponding regression coefficients can vary from one prediction day to another. The model allows thus for a day-to-day adaptive and tailored downscaling. It can also reveal specific predictors for peculiar and non-frequent weather configurations. [ABSTRACT FROM AUTHOR]

Many multi-site stochastic models have been proposed for the generation of daily precipitation, but they generally focus on the reproduction of low to high precipitation amounts at the stations concerned. This paper proposes significant extensions to the multi-site daily precipitation model introduced by Wilks, with the aim of reproducing the statistical features of extremely rare events (in terms of frequency and magnitude) at different temporal and spatial scales. In particular, the first extended version integrates heavy-tailed distributions, spatial tail dependence, and temporal dependence in order to obtain a robust and appropriate representation of the most extreme precipitation fields. A second version enhances the first version using a disaggregation method. The performance of these models is compared at different temporal and spatial scales on a large region covering approximately half of Switzerland. While daily extremes are adequately reproduced at the stations by all models, including the benchmark Wilks version, extreme precipitation amounts at larger temporal scales (e.g., 3-day amounts) are clearly underestimated when temporal dependence is ignored. [ABSTRACT FROM AUTHOR]

Rainfall is an important driver of erosion processes. The mean rainfall rate is often used to account for the erosive impact of a particular climate. However, for some erosion processes, erosion rate is a nonlinear function of rainfall, e.g., due to a threshold for erosion. When this is the case, it is important to take into account the full distribution of rainfall, instead of just the mean. In light of this, we have characterized the variability of daily rainfall over the Himalayan orogen using high spatial and temporal resolution rainfall data sets. We find significant variations in rainfall variability over the Himalayan orogen, with increasing rainfall variability to the west and north of the orogen. By taking into account variability of rainfall in addition to mean rainfall rate, we find a pattern of rainfall that, from a geomorphological perspective, is significantly different from mean rainfall rate alone. Using these findings, we argue that short-term rainfall variability may help explain observed short and long-term erosion rates in the Himalayan orogen. [ABSTRACT FROM AUTHOR]

Statistical Downscaling Methods (SDMs) are often used to produce local weather scenarios from large scale atmospheric information. SDMs include transfer functions which are based on a statistical link identified from observations between local weather and a set of large scale predictors. As physical processes generating surface weather vary in time, the most relevant predictors and the regression link are likely to also vary in time. This is well known for precipitation for instance and the link is thus often estimated after some seasonal stratification of the data. In this study, we present a hybrid model where the regression link is estimated from atmospheric analogs of the current prediction day. Atmospheric analogs are first identified from geopotential fields at 1000 and 500 hPa. For the regression stage, two Generalized Linear Models are further used to model the probability of precipitation occurrence and the distribution of non-zero precipitation amounts respectively. The hybrid model is evaluated for the probabilistic prediction of local precipitation over France. It noticeably improves the skill of the prediction for both precipitation occurrence and quantity. As the analog days vary from one prediction day to another, the atmospheric predictors selected in the regression stage and the value of the corresponding regression coefficients vary from one prediction day to another. The hybrid approach allows thus for a day-to-day adaptive and tailored downscaling. It can also reveal specific predictors for peculiar and non-frequent weather configurations. [ABSTRACT FROM AUTHOR]

Accurate estimates of flood peaks, corresponding volumes, and hydrographs are required to design safe and cost-effective hydraulic structures. In this paper, we propose a statistical approach for the estimation of the design variables peak and volume by constructing synthetic design hydrographs for different flood types such as flash-floods, short-rain floods, long-rain floods, and rain-on-snow floods. Our approach relies on the fitting of probability density functions to observed flood hydrographs of a certain flood type and accounts for the dependence between peak discharge and flood volume. It makes use of the statistical information contained in the data and retains the process information of the flood type. The method was tested based on data from 39 mesoscale catchments in Switzerland and provides catchment specific and flood type specific synthetic design hydrographs for all of these catchments. We demonstrate that flood type specific synthetic design hydrographs are meaningful in flood-risk management when combined with knowledge on the seasonality and the frequency of different flood types. [ABSTRACT FROM AUTHOR]

Estimates of flood event magnitudes with a certain return period are required for the design of hydraulic structures. While the return period is clearly defined in a univariate context, its definition is more challenging when the problem at hand requires considering the dependence between two or more variables in a multivariate framework. Several ways of defining a multivariate return period have been proposed in the literature, which all rely on different probability concepts. Definitions use the conditional probability, the joint probability, or can be based on the Kendall's distribution or survival function. In this study, we give a comprehensive overview on the tools that are available to define a return period in a multivariate context. We especially address engineers, practitioners, and people who are new to the topic and provide them with an accessible introduction to the topic. We outline the theoretical background that is needed when one is in a multivariate setting and present the reader with different definitions for a bivariate return period. Here, we focus on flood events and the different probability concepts are explained with a pedagogical, illustrative example of a flood event characterized by the two variables peak discharge and flood volume. The choice of the return period has an important effect on the magnitude of the design variable quantiles, which is illustrated with a case study in Switzerland. However, this choice is not arbitrary and depends on the problem at hand. WIREs Water 2016, 3:819-833. doi: 10.1002/wat2.1173 For further resources related to this article, please visit the . [ABSTRACT FROM AUTHOR]

River discharge is one of the most important quantities in hydrology. It provides fundamental records for water resources management and climate change monitoring. Even very short data-gaps in this information can cause extremely different analysis outputs. Therefore, reconstructing missing data of incomplete data sets is an important step regarding the performance of the environmental models, engineering, and research applications, thus it presents a great challenge. The objective of this paper is to introduce an effective technique for reconstructing missing daily discharge data when one has access to only daily streamflow data. The proposed procedure uses a combination of regression and autoregressive integrated moving average models (ARIMA) called dynamic regression model. This model uses the linear relationship between neighbor and correlated stations and then adjusts the residual term by fitting an ARIMA structure. Application of the model to eight daily streamflow data for the Durance river watershed showed that the model yields reliable estimates for the missing data in the time series. Simulation studies were also conducted to evaluate the performance of the procedure. [ABSTRACT FROM AUTHOR]

Despite many recent improvements, ensemble forecast systems for streamflow often produce under-dispersed predictive distributions. This situation is problematic for their operational use in water resources management. Many options exist for post-processing of raw forecasts. However, most of these have been developed using meteorological variables such as temperature, which displays characteristics very different from streamflow. In addition, streamflow data series are often very short or contain numerous gaps, thus compromising the estimation of post-processing statistical parameters. For operational use, a post-processing method has to be effective while remaining as simple as possible. We compared existing post-processing methods using normally distributed and gamma-distributed synthetic datasets. To reflect situations encountered with ensemble forecasts of daily streamflow, four normal distribution parameterizations and six gamma distribution parameterizations were used. Three kernel-based approaches were tested, namely, the 'best member' method and two improvements thereof, and one regression-based approach. Additional tests were performed to assess the ability of post-processing methods to cope with short calibration series, missing values or small numbers of ensemble members. We thus found that over-dispersion is best corrected by the regression method, while under-dispersion is best corrected by kernel-based methods. This work also shows key limitations associated with short data series, missing values, asymmetry and bias. One of the improved best member methods required longer series for the estimation of post-processing parameters, but if provided with adequate information, yielded the best improvement of the continuous ranked probability score. These results suggest guidelines for future studies involving real operational datasets. Copyright © 2014 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]