HS7.8
Hydrometeorologic stochastics for hydrologic applications: extremes, scales, probabilities

HS7.8

EDI
Hydrometeorologic stochastics for hydrologic applications: extremes, scales, probabilities
Convener: Hannes Müller-Thomy | Co-conveners: Alberto Viglione, Jose Luis Salinas Illarena, Auguste Gires, Gaby Gründemann
Presentations
| Tue, 24 May, 11:05–11:45 (CEST)
 
Room L2

Presentations: Tue, 24 May | Room L2

Chairpersons: Hannes Müller-Thomy, Jose Luis Salinas Illarena, Alberto Viglione
11:05–11:10
|
EGU22-549
|
Presentation form not yet defined
|
Giuseppe Formetta, Francesco Marra, Eleonora Dallan, Mattia Zaramella, and Marco Borga

Extreme precipitation in mountainous regions is the main trigger of hydrological hazards such as flash floods and debris flows, among the most dangerous natural-hazards worldwide, both for social and for economic losses. Mountains significantly influence weather and climate, including altered distribution of precipitation and of its extremes, with different impacts at different durations. 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, which are defined as the finite independent samples of the analyzed stochastic process (e.g. Marani and Ignaccolo, 2015), to improve our understanding of the orographic impact on extreme precipitation of durations ranging between 5 minutes and 24 hours. 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 validated our statistical framework against statistical properties of the observed annual maxima (Nash-Sutcliffe and Bias) 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 and an orographic enhancement for durations of ~8 hours or longer. 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 mid durations and long durations. These findings are of interest for risk management applications and climate change impact studies.

References

Marani, M., & Ignaccolo, M. (2015). A metastatistical approach to rainfall extremes. Advances in Water Resources, 79, 121-126.

 

How to cite: Formetta, G., Marra, F., Dallan, E., Zaramella, M., and Borga, M.: Differential orographic impact on sub-hourly, hourly, and daily extreme precipitation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-549, https://doi.org/10.5194/egusphere-egu22-549, 2022.

11:10–11:15
|
EGU22-570
|
ECS
|
Presentation form not yet defined
Andrea Magnini, Michele Lombardi, Elena Valtancoli, and Attilio Castellarin

Due to the limited length of locally available sequences of precipitation extremes, estimates of design rainstorms at a given location (i.e. point rainfall depth associated with given durations and non-exceedance probabilities) are traditionally obtained from regional frequency analysis. Several statistical regionalization methods proposed in the literature enable one to exploit sequences of precipitation extremes observed at a number of sites that supposedly share the same frequency regime of rainfall extremes with the site of interest (herein also referred to as a homogeneous pooling group of sites). Homogeneous pooling groups of sites can be identified by looking at specific climatic descriptors; for instance, some reliable authors successfully utilize Mean Annual Precipitation (MAP) as the sole proxy for locally characterizing the frequency regime of sub-daily rainfall extremes and for grouping sequences of rainfall extremes records. We aim at advancing this traditional approach (1) by relaxing the hypothesis of the existence of a homogeneous pooling group of sites characterized by a unique regional parent distribution and (2) by incorporating additional morphological and climatic information in the regional model. Our research focuses on a large study area in Northern Italy, counting more than 2350 Annual Maximum Series of rainfall depth for different time-aggregation intervals between 1 and 24 hours, that have been collected between 1928 and 2011 in the Italian Rainfall Extreme Dataset (I2-RED).  We refer to local MAP value as well as to several other morphologic descriptors (e.g. minimum distance to the coast, elevation of orographic barriers, aspect, terrain slope, etc.) for characterizing the frequency regime of sub-daily rainfall extremes. We train a probabilistic neural network that uses the descriptors cited above as input layers for modeling the local frequency regime of observed rainfall annual maxima. We resort to a Generalized Extreme Value (GEV) distribution whose parameters are data-driven functions of the local morphoclimatic descriptors as well as the time-aggregation interval. We then perform a series of cross-validation experiments targeted at assessing the accuracy of the developed data-driven regional frequency model relative to a simpler regional model in which GEV parameters are functions of MAP and time aggregation intervals.

Our results address the following research problems: (a) identification of the most descriptive morphological proxies for representing the frequency regime of sub-daily rainfall extremes, (b) assessment of potential, limitations, and robustness of data-driven multivariate regional frequency models of sub-daily rainfall extremes relative to simpler and more traditional regionalization schemes.

How to cite: Magnini, A., Lombardi, M., Valtancoli, E., and Castellarin, A.: Regional predictions of sub-daily rainfall extremes through data-driven blends of morphoclimatic descriptors, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-570, https://doi.org/10.5194/egusphere-egu22-570, 2022.

11:15–11:20
|
EGU22-3805
|
On-site presentation
|
Francesco Marra and Simon M. Papalexiou

The parent distribution of daily precipitation is usually not known, and the exceedance probability of extremes is described using a Generalized Extreme Value distribution (GEV) fitting the annual maxima. However, knowing the parent distribution would allow us to use ordinary statistics to describe extremes, with two advantages: (i) a decreased parameter estimation uncertainty; (ii) the possibility to establish direct relations between ordinary and extreme events. Recent studies suggest that daily precipitation could have Weibull tails, meaning that the probability of exceeding large values decrease as a stretched exponential. Here, we exploit a global dataset of long and quality-controlled continuous rain gauge records (~8,000 stations, ≥50 complete years) to investigate this question.

We find that the observed annual maxima are likely samples from Weibull tails in ~88% of the stations worldwide. On average, ~36% of the wet days belong to these tails. We find a strong climatic dependence in their definition, with smaller portions of data in the Weibull tails in central Europe, US east coast and southern Australia. We then generate synthetic records with the same characteristics (yearly number of wet days, Weibull tails with the same shape parameter) and increasing lengths (10-110 years); we estimate the corresponding GEV shape parameters and contrast them with the ones obtained from very long annual maxima records (~15,000 stations, 40-163 years; Papalexiou and Koutsoyiannis, 2013, https://doi.org/10.1029/2012WR012557). We show that parent distributions with Weibull tails well explain the properties of the observed GEV shape parameters. These GEV tails (type-III, Frechet) are heavier than the limiting GEV for Weibull parent distributions (type-I, Gumbel); this implies a pre-asymptotic behavior: the average yearly number of wet days (globally, n=100±50) is not large enough to fulfill the asymptotic assumption (n~∞) of extreme value theory. Contrasting our results with generalized Pareto tails, as predicted by extreme value theory for high threshold exceedances, we find that the two models are equivalent within the observational uncertainties; the Weibull model, however, describes a portion of data which is, on average, 7 times larger.

How to cite: Marra, F. and Papalexiou, S. M.: Daily precipitation with stretched-exponential tails could explain the statistics of observed annual maxima, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-3805, https://doi.org/10.5194/egusphere-egu22-3805, 2022.

11:20–11:25
|
EGU22-2738
|
ECS
|
On-site presentation
Golbarg Goshtsasbpour, Uwe Haberlandt, Abbas El Hachem, Jochen Seidel, and András Bárdossy

Precipitation extremes are space-time phenomena and traditionally the statistical analyses on such occurrences have treated them merely as point events. Many of the consequences of such events like floods are related to the water volume, hence the spatial aspect of them cannot be neglected. This work aims to bring the areal aspect of the extreme rainfall into play by introducing the area into the Extreme Value Analysis (EVA) and providing Area-Depth-Duration-Frequency (ADDF) curves. For this purpose, different spatial rainfall products have been used and compared with each other. Processed raw radar data, a product of conditional merging of the radar and station data as well as the RADKLIM data (a product of the German Weather Service designed for climate research) have been used for the EVA. Unexpected patterns have been observed in the ADDF curves based on the processed radar data which were not in agreement with the assumptions of the classical approach of areal reduction factor.  Usually, it is assumed that areal precipitation extremes increase with decreasing area, so in practice reduction factors are used to estimate areal precipitation extreme values from point observations. This behavior was observed as expected mostly for durations shorter than 2 hours in all the study locations whereas the opposite was present for longer durations, where the precipitation is increasing with increasing area, so that the ADDF curves, representing different areas, show crossings at these durations. Different hypotheses about the reason for the crossings like seasonality and spatial non-stationarity have been tested and did not explain the crossings. On the other hand, the ADDF curves of the merged rainfall product hardly showed such patterns and followed the classical assumptions. Therefore, the appearance of such crossings in the ADDF curves of a spatial rain product might be an indicator of artifacts in the radar rainfall product. It has to be investigated in further tests if these results hold and if these crossings could be used as an indicator for unplausible radar data.  

How to cite: Goshtsasbpour, G., Haberlandt, U., El Hachem, A., Seidel, J., and Bárdossy, A.: Statistical Analysis of Space-times Dynamics of Extreme Precipitation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2738, https://doi.org/10.5194/egusphere-egu22-2738, 2022.

11:25–11:30
|
EGU22-987
|
ECS
|
On-site presentation
Bora Shehu and Uwe Haberlandt

Rainfall intensity-duration-frequency (IDF) curves are required for the design of several water systems and protection works. Typically, long (more than 40 years) station data are employed first to generate annual extremes (AMS) for different durations and then to fit a GEV probability distribution. Since station data are only point measurements, regionalization techniques are applied to estimate IDF curves at ungauged locations.  Prior results revealed that the best way to obtain IDF maps for Germany was kriging interpolation of parameters from very long stations with the parameters of the short stations acting as an external drift. However, how certain the obtained IDFs values are, and how to derive the uncertainty range at each location, remain still unanswered. Therefore, it is the objective of this study, to investigate the propagation of uncertainty in the regionalization of the IDF curves for Germany.

For this purpose, the available station data from the German Weather Service (DWD) for whole Germany are employed, which includes; 1100 sub-hourly (5min) recordings with observations period shorter than 20 years, and 89 sub-hourly (5min) recordings with 60-70 years of observations. Annual extremes are extracted at each location for different durations (from 5mins up to 7days), and local IDF curves are estimated according to Koutsoyiannis et al. (1998). The parameters of the obtained IDF functions are then interpolated using external drift kriging. Finally, quantiles are derived for each location, duration and given return period (Ta=2, 10, 20, 50 and 100 years). Through a non-parametric bootstrap, the uncertainty is estimated for three different components of the regionalization: i) local estimation of parameters, ii) variogram estimation and iii) spatial sampling distribution.  Simulated annealing is employed to ensure that the spatial resampling of locations represents the obtained variograms. The final uncertainty range is then considered as the 95% confidence interval of the obtained IDF curve for each location, duration and return period.      

The results of this study reveal how the uncertainty of annual rainfall extremes propagates from local estimation to the regionalization of IDF curves based on kriging theory. The comparison of the three components will shed light to the following questions: Which is the contribution of each component to the final uncertainty range of IDFs curve? How is the uncertainty range changing based on different durations and return periods? Are there any spatial trends in Germany regarding the uncertainty range of IDFs curves? 

How to cite: Shehu, B. and Haberlandt, U.: Uncertainty analysis of regionalized intensity-duration-frequency curves in Germany, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-987, https://doi.org/10.5194/egusphere-egu22-987, 2022.

11:30–11:35
|
EGU22-7181
|
Virtual presentation
|
Thordis Thorarinsdottir, Thea Roksvåg, Julia Lutz, Lars Grinde, and Anita Dyrrdal

As a warming climate leads to more frequent heavy rainfall, the importance of accurate rainfall statistics is increasing. Rainfall statistics are often presented as intensity-duration-frequency (IDF) curves showing the rainfall intensity (return level) that can be expected at a location for a duration, and the frequency of this intensity (return period). IDF curves are commonly constructed by fitting generalized extreme value (GEV) distributions to observed annual maximum rainfall for several target durations, where the available observation data sources may vary for the different durations. As the estimation is performed independently across durations, the resulting IDF curves may be inconsistent across durations and return periods. We discuss how consistent estimates across the different durations may be derived by post-processing independently obtained Bayesian posterior distributions for each duration. The proposed methods are evaluated for simulated data and for Norwegian rainfall data from 83 locations, for 16 durations between 1 minute and 24 hours, where the post-processing yields consistent and accurate estimates.

How to cite: Thorarinsdottir, T., Roksvåg, T., Lutz, J., Grinde, L., and Dyrrdal, A.: A Bayesian framework to derive consistent intensity-duration-frequency curves from multiple data sources, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7181, https://doi.org/10.5194/egusphere-egu22-7181, 2022.

11:35–11:40
|
EGU22-2423
|
ECS
|
On-site presentation
|
Felix Fauer, Jana Ulrich, and Henning Rust

Extreme precipitation is currently the biggest climate-change-related thread in middle Europe with flooding events leading to high death tolls and huge existential and financial losses. Evaluating the probabilities of these extremes can help preventing casualties and reducing impact consequences. Our analysis is based on Intensity-Duration-Frequency (IDF) curves which describe the major statistical characteristics of extreme precipitation events. They provide information on the probability of exceedance of certain precipitation intensities for a range of durations and can help to visualize how extreme the event for different durations is.  A popular way of evaluating return periods of extremes is to model the underlying distribution of block maxima with the Generalized Extreme Value (GEV) distribution. A core problem when modeling extremes is the scarce availability of data. This can be addressed by using the available data more efficiently. Therefore, including maxima from different measurement durations is useful for (1) gathering more information from the data and (2) estimating return periods for different time scales with a consistent modeling approach. Duration-dependence is implemented directly into the parameter estimation (Koutsoyiannis et al., 1998) and enables a consistent model, i.e. without quantile-crossing. We present a new parameterization of the duration-dependent GEV (d-GEV) that is more flexible with respect to long ranges of durations and is considering the different time scales on which extremes occur in winter and summer (Fauer et al., 2021). Applying the new model to the extreme rain event on 14 July 2021 in Ahrtal, Germany reveals that the event was most extreme on a time scale of 20-30 hours.

Investigating the impact of large-scale atmospheric flow on extremes will help to learn how extremes changed in the past and make projections about their change in the future. Large scale variables are incorporated into the model as covariates in generalized linear models for the d-GEV parameters. An ongoing study tests for the inclusion of NAO, a blocking index, monthly mean temperature, etc., as predictor variables. First results show a significant correlation (5%-level) between monthly precipitation maxima and NAO/blocking for some durations and some seasons. It will be analyzed whether this connection can be useful for modelling d-GEV parameters with large-scale variables.

How to cite: Fauer, F., Ulrich, J., and Rust, H.: Efficient usage of information for modeling precipitation extremes and large scale influence, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-2423, https://doi.org/10.5194/egusphere-egu22-2423, 2022.

11:40–11:45
|
EGU22-6720
|
On-site presentation
Ross Pidoto and Uwe Haberlandt

Long term time series of meteorologic variables are generally lacking and is especially the case at sub-daily temporal resolutions. These time series are needed for applications such as hydrological modelling of catchments and derived flood frequency analyses. Stochastic weather generators are one such solution and are able to generate long synthetic time series of arbitrary length.

This study explores the coupling of an hourly space-time rainfall model with a non-parametric K-NN resampling approach for non-rainfall climate variables such as temperature, humidity and global radiation. A daily gridded observational climate dataset is used for the resampling. As a last step, a simple disaggregation technique is applied to the resampled non-rainfall climate variables to achieve an hourly timestep.

To study the effectiveness and performance of the hybrid weather generator, synthetic time series for 400 mesoscale catchments within Germany consisting of 700 rainfall stations were generated and compared to observations. Results show that the hybrid weather generator adequately reproduces observed statistics for rainfall and the non-rainfall climate variables in addition to maintaining cross correlations between the climate variables.

How to cite: Pidoto, R. and Haberlandt, U.: A multi-scale space-time hybrid weather generator, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-6720, https://doi.org/10.5194/egusphere-egu22-6720, 2022.