Regional frequency analysis of daily rainfall extremes using L-moments approach

Daily extreme precipitation values are among environmental events with the most disastrous consequences for human society. Information on the magnitudes and frequencies of extreme precipitations is essential for sustainable water resources management, planning for weather-related emergencies, and design of hydraulic structures. In the present study, regional frequency analysis of maximum daily rainfalls was investigated for Golestan province located in the northeastern Iran. This study aimed to find appropriate regional frequency distributions for maximum daily rainfalls and predict the return values of extreme rainfall events (design rainfall depths) for the future. L-moment regionalization procedures coupled with an index rainfall method were applied to maximum rainfall records of 47 stations across the study area. Due to complex geographic and hydro-climatological characteristics of the region, an important research issue focused on breaking down the large area into homogeneous and coherent sub-regions. The study area was divided into five homogeneous regions, based on the cluster analysis of site characteristics and tests for the regional homogeneity. The goodness-of-fit results indicated that the best fitting distribution is different for individual homogeneous regions. The difference may be a result of the distinctive climatic and geographic conditions. The estimated regional quantiles and their accuracy measures produced by Monte Carlo simulations demonstrate that the estimation uncertainty as measured by the RMSE values and 90% error bounds is relatively low when return periods are less than 100 years. But, for higher return periods, rainfall estimates should be treated with caution. More station years, either from longer records or more stations in the regions, would be required for rainfall estimates above T=100 years. It was found from the analyses that, the index rainfall (at-site average maximum rainfall) can be estimated reasonably well as a function of mean annual precipitation in Golestan province. Index rainfalls combined with the regional growth curves, can be used to estimate design rainfalls at ungauged sites. Overall, it was found that cluster analysis together with the L-moments based regional frequency analysis technique could be applied successfully in deriving design rainfall estimates for northeastern Iran. The approach utilized in this study and the findings are of great scientific and practical merit, particularly for the purpose of planning for weather-related emergencies and design of hydraulic engineering structures.     

Rainfall-triggered landslides in an unsaturated soil: a laboratory flume study

Extreme and/or prolonged rainfall events frequently cause landslides in many parts of the world. In this study, infiltration of rainfall into an unsaturated soil slope and triggering of landslides is studied through laboratory model (flume) tests, with the goal of obtaining the triggering rainfall intensity–duration (I–D) threshold. Flume tests with fine sand at two different relative densities (34 and 48%) and at slope angle of 56.5° are prepared, and rainfall (intensity in the range of 18 to 64 mm/h) is applied via a mist sprinkler system to trigger landslides. Soil water characteristic curve and hydraulic conductivity function of the fine sand are also presented. In flume tests, suction in the soil is measured with tensiometers, the progress of wetting front with time and deformations in the soil are also measured. Some of the findings of this study are: for the fine sand used in this study (a) the failure mechanism is infinite-slope type (mostly translational), and the failure surface is generally coincident with the wetting front or is in its vicinity, (b) the deformations leading to a landslide occurred abruptly, (c) both relatively high-intensity–short-duration rainfalls and relatively low-intensity–long duration rainfalls triggered landslides, (d) the shape of the I–D threshold is demonstrated to be a bilinear relation in log intensity–log duration plot, (e) below a certain rainfall intensity landslides are not triggered, (f) the effect of relative density of the soil on the I–D threshold is demonstrated by physical laboratory tests (as the relative density of the soil increases, the triggering rainfall intensity–duration threshold moves to larger rainfall events). The results of this study could be useful for accurate numerical modeling of rainfall-triggered landslides.     

Analysis of rainfall preceding debris flows on the Smědavská hora Mt., Jizerské hory Mts., Czech Republic

In August 2010, extreme rainfall affected the north of the Czech Republic and caused regional floods and landslides. Three torrential debris flows originated in the Jizerské hory Mts., close to Bílý Potok on the north slope of the Smědavská hora Mt. The rainfall situation which triggered the debris flow was analyzed and compared with the rainfall situation in 1958 when a debris flow occurred in the same area. The rainfall data were obtained from rain gauges of the Czech Hydrometeorological Institute. Four rain gauges were chosen close to the Smědavská hora Mt. with data of daily amounts from 1983 to 2013 and 10-min intensity or hourly amounts from the specific period. The data from 1958 were available from three different rain gauges (only daily amounts). The data series were not complete so linear regression was applied to interpolate them. A number of analyses were carried out including daily rainfall, 2-day/3-day moving values, antecedent precipitation index (API) of 5/10/30 days, 10-min intensity, and hourly amounts, and the trigger factor of the debris flow in the study area was also investigated. It was determined that for the triggering of debris flows, both high API values as well as high-intensity short-duration rainfall is needed. It was documented that in cases of solely high API indices or high-intensity short-duration rainfalls, no debris flows were initiated.     

Quantifying uncertainties associated with depth duration frequency curves

Uncertainty in depth–duration–frequency (DDF) curves is usually disregarded in the view of difficulties associated in assigning a value to it. In central Iran, precipitation duration is often long and characterized with low intensity leading to a considerable uncertainty in the parameters of the probabilistic distributions describing rainfall depth. In this paper, the daily rainfall depths from 4 stations in the Zayanderood basin, Iran, were analysed, and a generalized extreme value distribution was fitted to the maximum yearly rainfall for durations of 1, 2, 3, 4 and 5 days. DDF curves were described as a function of rainfall duration (D) and return period (T). Uncertainties of the rainfall depth in the DDF curves were estimated with the bootstrap sampling method and were described by a normal probability density function. Standard deviations were modeled as a function of rainfall duration and rainfall depth using 10⁴ bootstrap samples for all the durations and return periods considered for each rainfall station.     

Regional Rainfall Depth–Duration–Frequency Equations for an Alpine Region

In this paper the scaling hypotheses are applied to annual maximum series of rainfall depth for different rainfall duration to derive the depth–duration–frequency (DDF) curve. It is shown that, based on the empirically observed scaling properties of rainfall and some general assumptions about the cumulative distribution function for the annual maximum of the rainfall depth, it is possible to derive a simple DDF relationship. This general framework provides a basis for the generation of maps that can be used to infer DDF curves at any point of a particular area. Data from a dense raingauge network in a mountainous region in north-eastern Italy (the Trentino Province) are used to clarify the methodology for the construction and regionalization of the DDF relationship. The geographical variation of short-duration (i.e., less than 60 xmin) rainfall extremes is also evaluated by using the same framework. It is found that depth–duration ratios, defined as the ratios of the t-min to the 60-min rainfall depths of the same return period, may be considered independent of return period and geographical location for any storm duration less than 60 min.     

Rainfall characteristics and thresholds for periglacial debris flows in the Parlung Zangbo Basin,southeast Tibetan Plateau

The Parlung Zangbo Basin in the southeastern Tibet Plateau is affected by the summer monsoon from the Indian Ocean, which produces large rainfall gradients in the basin. Rainfall data during 2012–2015 from five new meteorological stations are used to analyse the rainfall characteristics. The daily rainfall, rainfall duration, mean rainfall intensity, and peak rainfall intensity are consistent, but sometimes contrasting. For example, these values decrease with increasing altitude, and the gradient is large downstream and small upstream, respectively. Moreover, the rainfall intensity peaks between 01:00 and 06:00 and increases during the afternoon. Based on the analysis of 14 debris flow cases in the basin, differences in the rainfall threshold differ depending on the location as sediment varieties. The sediment in the middle portions of the basin is wet and well structured; thus, long-duration, high-intensity rainfall is required to generate debris flows. Ravels in the upstream area are arid and not well structured, and short-duration rainfall is required to trigger debris flows. Between the above two locations, either long-duration, low-intensity rainfall or short-duration, high-intensity rainfall could provoke debris flows. Clearly, differences in rainfall characteristics and rainfall thresholds that are associated with the location must be considered in debris flow monitoring and warnings.     

Examination of influences of rainfall patterns on shallow landslides due to dissipation of matric suction

The influences of rainfall patterns on shallow landslides due to the dissipation of matric suction are examined in this study. Four representative rainfall patterns including the uniform, advanced, intermediated, and delayed rainfalls are adopted. The results show that not only the occurrence of shallow landslides but also the failure depth and the time of failure are affected by the rainfall pattern. The different rainfall patterns seem to have the same minimum landslide-triggering rainfall amount. There is a rainfall duration threshold for landslide occurrence for a rainfall event with larger than the minimum landslide-triggering rainfall amount. For each rainfall pattern, the rainfall duration threshold for landslide occurrence decreases to constant with the increase of rainfall amount. The uniform rainfall has the least rainfall duration threshold for landslide occurrence, followed by the advanced rainfall, and then the intermediated rainfall. For each rainfall pattern, the failure depths and the times of failure from the same amount of rainfall with different durations could be largely different. In addition, the differences of the failure depths and the times of failure between various rainfall patterns with the same amount and duration of rainfall could be also significant. The failure depth and the time of failure, as compared with the occurrence of shallow landslides, are more sensitive to the rainfall condition. In other words, in comparison with the evaluation of the occurrence of shallow landslides, it needs more accurate rainfall prediction to achieve reliable estimations of the failure depth and the time of failure.     

世界实测与调查最大点雨量及其外包线公式

介绍了世界范围内大量实测与调查的点雨量极值,据此求得了世界最大点雨量外包线的新公式R=491D0.452,并对世界最大点雨量的地区分布特点、暴雨成因以及点暴雨极值在可能最大降水估算中的运用作了简要的说明.     

The rainfall intensity–duration control of shallow landslides and debris flows: an update

A global database of 2,626 rainfall events that have resulted in shallow landslides and debris flows was compiled through a thorough literature search. The rainfall and landslide information was used to update the dependency of the minimum level of rainfall duration and intensity likely to result in shallow landslides and debris flows established by Nel Caine in 1980. The rainfall intensity–duration (ID) values were plotted in logarithmic coordinates, and it was established that with increased rainfall duration, the minimum average intensity likely to trigger shallow slope failures decreases linearly, in the range of durations from 10 min to 35 days. The minimum ID for the possible initiation of shallow landslides and debris flows was determined. The threshold curve was obtained from the rainfall data using an objective statistical technique. To cope with differences in the intensity and duration of rainfall likely to result in shallow slope failures in different climatic regions, the rainfall information was normalized to the mean annual precipitation and the rainy-day normal. Climate information was obtained from the global climate dataset compiled by the Climate Research Unit of the East Anglia University. The obtained global ID thresholds are significantly lower than the threshold proposed by Caine (Geogr Ann A 62:23–27, 1980), and lower than other global thresholds proposed in the literature. The new global ID thresholds can be used in a worldwide operational landslide warning system based on global precipitation measurements where local and regional thresholds are not available..     

暴雨和降水偏态系数分析

采用统一的统计参数计算方法分析计算了全国0.51万~1.36万个雨量站历时为10 min~3 d的暴雨资料、210个站历时为1 d~90 d和1 a的降水资料的偏态系数和变差系数比值R_sv。探讨了皮尔逊Ⅲ型频率分布曲线的R_sv特性;分析了R_sv的地域分布、R_sv与雨量历时的关系,并研究了对R_sv有影响的多个因子。研究表明,在目前资料条件下,不宜对各个地区各种历时的暴雨频率分析均采用单一的R_sv值(3.5)。     

Setting up evaluate indicators for slope control engineering based on spatial clustering analysis

Rainfall-induced landslides (RILs) have been a source of social and economic disruption in the mountainous Baguio area in northern Philippines. Prolonged heavy rainfall usually happens during tropical cyclone and southwest monsoon activity. A pragmatic approach to RIL mitigation is to develop rainfall-based early warning. We implemented a modified regression method to derive the empirical minimum intensity (I)–duration (D) threshold I = 6.46 D ^?0.28 and a normalized ID threshold NI = 0.002 D ^?0.28 for rainfall duration ranging between 24 and 264 h. Using a separate data set to evaluate the applicability of the threshold, 93% of the landslide-triggering rainfall events fell above the derived threshold. RILs also occurred when 24-h rainfall was 0.02–28% of the mean annual precipitation or after accumulating at least 500 mm of rainfall from the onset of the rainy season. The thresholds may be further refined as more landslide data become available in the future.     

北京山区基于历史资料的泥石流临界雨量研究

通过对北京山区历史泥石流降雨资料的分析,根据有雨量站和无雨量站的泥石流易发区的资料情况,分别建立了不同前期雨量的预报模型,并根据预报模型初步建立了不同前期雨量前提条件下的黄,橙,红三色预报模式。     

Rainfall patterns for shallow landsliding in perialpine Slovenia

This paper presents two types of analysis: an antecedent rainfall analysis based on daily rainfall and an intensity-duration analysis of rainfall events based on hourly data in perialpine Slovenia in the ?kofjelo?ko Cerkljansko hills. For this purpose, eight rainfall events that are known to have caused landslides in the period from 1990 to 2010 were studied. Over the observed period, approximately 400 records of landslides were collected. Rainfall data were obtained from three rain gauges. The daily rainfall from the 30 days before landslide events was investigated based on the type of landslides and their geo-environmental setting, the dates of confirmed landslide activity and different consecutive rainfall periods. The analysis revealed that the rainfall events triggering slope failure can be divided into two groups according to the different antecedent periods. The first group of landslides typically occurred after short-duration rainstorms with high intensity, when the daily rainfall exceeded the antecedent rainfall. The second group comprises the rainfall events with a longer antecedent period of at least 7 days. A comparison of the plotted peak and mean intensities indicates that the rainfall patterns that govern slope failure are similar but do not necessarily reflect the rainfall intensity at the time of shallow landslides in the Dav?a or Poljane areas, where the majority of the landslides occurred. Because of several limitations, the suggested threshold cannot be compared and evaluated with other thresholds.     

Rainfall thresholds for prediction of shallow landslides around Chamoli-Joshimath region,Garhwal Himalayas,India

Majority of landslides in the Indian sub-continent are triggered by rainfall. Several attempts in the global scenario have been made to establish rainfall thresholds in terms of intensity-duration and antecedent rainfall models on global, regional and local scales for the occurrence of landslides. However, in the context of the Indian Himalayas, the rainfall thresholds for landslide occurrences are not yet understood fully. Neither on regional scale nor on local scale, establishing such rainfall thresholds for landslide occurrences in Indian Himalayas has yet been attempted. This paper presents an attempt towards deriving local rainfall thresholds for landslides based on daily rainfall data in and around Chamoli-Joshimath region of the Garhwal Himalayas, India. Around 128 landslides taken place in last 4 years from 2009 to 2012 have been studied to derive rainfall thresholds. Out of 128 landslides, however, rainfall events pertaining to 81 landslides were analysed to yield an empirical intensity–duration threshold for landslide occurrences. The rainfall threshold relationship fitted to the lower boundary of the landslide triggering rainfall events is I?=?1.82 D ^?0.23 (I?=?rainfall intensity in millimeters per hour and D?=?duration in hours). It is revealed that for rainfall events of shorter duration (≤24 h) with a rainfall intensity of 0.87 mm/h, the risk of landslide occurrence in this part of the terrain is expected to be high. Also, the role of antecedent rainfall in causing landslides was analysed by considering daily rainfall at failure and different period cumulative rainfall prior to failure considering all 128 landslides. It is observed that a minimum 10-day antecedent rainfall of 55 mm and a 20-day antecedent rainfall of 185 mm are required for the initiation of landslides in this area. These rainfall thresholds presented in this paper may be improved with the hourly rainfall data vis-à-vis landslide occurrences and also data of later years. However, these thresholds may be used in landslide warning systems for this particular region of the Garhwal Himalayas to guide the traffic and provide safety to the tourists travelling along this pilgrim route during monsoon seasons.     

Application of Monte Carlo simulation technique for flood estimation for two catchments in New South Wales,Australia

The currently adopted rainfall-based design flood estimation method in Australia, known as design event approach (DEA), has a flaw that is widely criticized by the hydrologists. The DEA is based on the assumption that a rainfall depth of a certain frequency can be transformed to a flood peak of the same frequency by adopting the ‘representative values’ of other model input variables, such as temporal patterns and losses. To overcome the limitation associated with the DEA, this paper develops stochastic model inputs to apply Monte Carlo simulation technique (MCST) for design flood estimation. This uses data from 86 pluviograph stations and six catchments from eastern New South Wales (NSW), Australia, to regionalize the distributions of various input variables (e.g., rainfall duration, inter-event duration, intensity and temporal patterns and loss and routing characteristics) to simulate thousands of flood hydrographs using a nonlinear runoff routing model. The regionalized stochastic inputs are then applied with the MCST to two catchments in eastern NSW. The results indicate that the developed MCST provide more accurate flood quantile estimates than the DEA for the two test catchments. The particular advantage of the new MCST and stochastic design input variables is that it reduces the subjectivity in the selection of model input values in flood modeling. The developed MCST can be adapted to other parts of Australia and other countries.     

Comparison of statistical methods and multi-time validation for the determination of the shallow landslide rainfall thresholds

Shallow landslides are unforeseeable phenomena often resulting in critical conditions in terms of people’s safety and damage. The main purpose of this paper is the comparison of different statistical methods used to determine the rainfall thresholds for the shallow landslide occurrence. Rainfall data over a 46-year period were collected for one rain gauge located in a test area of northwest Italy (Riviera Spezzina; RS). In the RS, intense rainfalls often induce shallow landslides causing damage and sometimes casualties. The rainfall events occurred in the 1967–2006 period were classified as events inducing shallow landslides (SLEs_1967–2006) and events that did not trigger shallow landslides (NSLEs_1967–2006). Thresholds for various percentiles of SLEs_1967–2006 were computed by identifying the lower limit above which shallow landslides occurred. Another set of thresholds, corresponding to different probabilities of occurrence, was determined using SLEs_1967–2006 and NSLEs_1967–2006. The least-squares linear fit (LSF) and the quantile regression (QR) techniques were employed in the former approach, while the logistic regression (LR) was applied in the latter. The thresholds were validated with the same data used for their definition and with the data recorded in the 2008–2014 period. Contingency tables were created and contingencies and skill scores were computed. The 10% probability threshold obtained using the LR method is characterized by the best values of at least two skill scores for both periods considered; therefore, it may be considered the “best” threshold for the RS. The results of this work can help the choice of the best statistical method to determine the shallow landslide rainfall thresholds.     

Flood frequency analysis based on simulated peak discharges

Flood frequency approaches vary from statistical methods, directly applied on the observed annual maximum flood series, to adopting rainfall–runoff simulation models that transform design rainfalls to flood discharges. Reliance on statistical flood frequency analysis depends on several factors such as the selected probability distribution function, estimation of the function parameters, possible outliers, and length of the observed flood series. Through adopting the simulation approach in this paper, watershed-average rainfalls of various occurrence probabilities were transformed into the corresponding peak discharges using a calibrated hydrological model. A Monte Carlo scheme was employed to consider the uncertainties involved in rainfall spatial patterns and antecedent soil moisture condition (AMC). For any given rainfall depth, realizations of rainfall spatial distribution and AMC conditions were entered as inputs to the model. Then, floods of different return periods were simulated by transforming rainfall to runoff. The approach was applied to Tangrah watershed in northeastern Iran. It was deduced that the spatial rainfall distribution and the AMCs exerted a varying influence on the peak discharge of different return periods. Comparing the results of the simulation approach with those of the statistical frequency analysis revealed that, for a given return period, flood quantiles based on the observed series were greater than the corresponding simulated discharges. It is also worthy to note that existence of outliers and the selection of the statistical distribution function has a major effect in increasing the differences between the results of the two approaches.     

中国极端降水事件的频数和强度特征

使用1951-2004年中国738个测站逐日降水资料,采用百分位的方法定义极端降水事件的阈值,分析了不同持续时间的极端降水事件的时空分布及变化趋势特征。结果表明,极端降水事件多发于35°N以南,特别是在长江中下游和江南地区以及高原东南部,且在这些地区极端降水事件持续时间也较长。季节分布上,主要出现在夏季,以低持续性事件为主。在中国东部地区,持续时间越长的极端降水其强度往往越强。趋势分析表明,全国持续1d极端事件的相对频数具有上升趋势而平均强度具有下降趋势,其空间上均表现为全国大部分上升、华北和西南等地下降的趋势。持续2d以上极端事件在长江中下游流域、江南地区和高原东部等地区有显著增多和增强的趋势,而在华北和西南地区有减少和减弱趋势,但全国平均的趋势不显著。     

甘肃省舟曲县"8.8"特大山洪泥石流灾害与气象成因分析及其应对建议

通过调查舟曲地质灾害历史资料.分析降雨与地质灾害的相关性,发现地质灾害与灾害前1～3 h累积降水量、月平均降水量和月暴雨频数明显正相关.然后分析舟曲"8.8"特大山洪泥石流的气象和降水成因.在此基础上,提出舟曲地质灾害的防治措施.     

变化环境对城市暴雨及排水系统影响研究进展

近年来,变化环境（气候变化和城镇化）导致城市暴雨的发生频率或强度增加,加剧了城市暴雨洪涝问题。因此,分析变化环境对城市暴雨及排水系统的影响对城市水资源规划管理、市政规划设计和城市防灾减灾有着重要意义。总结了未来高精度降雨预估技术、变化环境下排水系统设计标准等关键问题的研究进展,综述了变化环境下城市短历时暴雨演变规律及变化环境对城市排水系统影响的研究现状,探讨了变化环境对城市短历时暴雨的影响机理,并归纳了当前研究变化环境对城市暴雨及排水系统影响的主要方法。指出今后应重点加强的研究:①加强气候变化和城镇化对城市短历时暴雨影响机制的研究;②提升区域气候模式对城市区域下垫面和大气相互作用的描述能力,并加强公里尺度对流可解析模型在城市气候变化影响研究中的应用;③加强对变化环境下排水系统设计标准的研究;④综合评估气候变化和城镇化对城市排水系统超载、污染物转移和水安全等的影响。