Series Of Rainfall Events Research Articles

Screening of representative rainfall event series for long-term hydrological performance evaluation of grassed swales.

Evaluation of the hydrological performance of grassed swales usually needs long-term monitoring data. At present, suitable techniques for simulating the hydrological performance using limited monitoring data are not available. Therefore, current study aims to investigate the relationship between saturated hydraulic conductivity (Ks) fitting results and rainfall characteristics of various events series length. Data from a full-scale grassed swale (Enschede, the Netherlands) were utilized as long-term rainfall event series length (95 rainfall events) on the fitting outcomes. Short-term rainfall event series were extracted from these long-term series and used as input in fitting into a multivariate nonlinear model between Ks and its influencing rainfall indicators (antecedent dry days, temperature, rainfall, rainfall duration, total rainfall, and seasonal factor (spring, summer, autumn, and winter, herein refer as 1, 2, 3, and 4). Comparison of short-term and long-term rainfall event series fitting results allowed to obtain a representative short-term series that leads to similar results with those using long-term series. A cluster analysis was conducted based on the fitting results of the representative rainfall event series with their rainfall event characteristics using average values of influencing rainfall indicators. The seasonal index (average value of seasonal factors) was found to be the most representative short rainfall event series indicator. Furthermore, a Bayesian network was proposed in the current study to predict if a given short-term rainfall event series is representative. It was validated by a data series (58 rainfall events) from another full-scale grassed swale located in Utrecht, the Netherlands. Results revealed that it is quite promising and useful to evaluate the representativeness of short-term rainfall event series used for long-term hydrological performance evaluation of grassed swales.

Hydrological and erosion responses of steep spoil heaps to taproot and fibrous root grasses under simulated rainfalls

Spoil heaps, resulting from excavation and backfilling activities at construction sites, are one of the most important risk sources of soil erosion, threatening local ecosystems and even leading to geological disasters. However, after revegetation, spoil heaps can be transformed into available land resources such as cultivated land and forest land. Previous studies have reported the runoff and sediment reduction benefits of different vegetation measures on spoil heaps, while little information has been documented about what differences in the taproot and fibrous root grasses were transferred into the differences in hydrological and erosion processes. This study compared the effects of taproot (Artemisia gmelinii (A. gmelinii)) and fibrous root (Cynodon dactylon (C. dactylon)) grasses on soil erosion and runoff production processes and the reduction benefits of spoil heaps with silt loam soil. A series of rainfall events at three intensities (0.8, 1.2, and 1.8 mm min−1) were simulated on runoff plots (length of 3.0 m and width of 1.0 m) under treatments of A. gmelinii, C. dactylon, and bare soil (BS, as control) (each treatment was replicated), and the initial runoff time (IRT), saturated hydraulic conductivity, runoff rate (RR), flow velocity, and soil loss rate (SLR) were measured. The results showed that grasses delayed the IRTs on spoil heaps by 0.69–12.96 min compared BS, and the delay effects were more significant for fibrous root grass than for taproot grass. The soil loss reduction benefits (91.1 %) of the two kinds of grass were approximately 2.2 times greater than those of runoff volume reduction (42.1 %), and shear stress could be used to predict the soil loss rate following power functions (R2 = 0.72–0.85) for the grass-covered plots. The taproot grass exhibited a higher erosion reduction than fibrous root grass and is thus preferred to control the soil erosion of spoil heaps. Taproot grass has a greater stem diameter and can insert deeper into spoil heaps to consolidate the soil particles, which can improve roughness and hydraulic resistance, thereby producing a lower effective hydraulic shear stress that is directly exerted on soil particles. This study enhances the understanding of different grasses on the water erosion process and mechanism of spoil heaps and provides scientific guidance for selecting erosion control practices.

Short-Term Effects of Forest Fire on Water Quality along a Headwater Stream in the Immediate Post-Fire Period

Changes in water quality were examined during selected rainfall events in a headwater stream severely damaged by a forest fire on 21–23 February 2021. Seven water quality parameters were analyzed: pH, electrical conductivity, biochemical oxygen demand, chemical oxygen demand, total suspended solids (TSS), total nitrogen, and total phosphorous. First-flush effect and event mean concentration analyses were conducted in spring (dry season) and summer (wet season) immediately after the forest fire. In spring, the non-biodegradable organic matter concentrations increased along with pollutants related to ash-enhanced soil water repellency in water affected by first-flush effects. In summer, TSS and nutrient concentrations increased along with pollutants related to fire-induced soil surface disturbances after a series of rainfall events. First-flush analyses indicated that cumulative pollutant loads were greater in the spring than in the summer due to a higher concentration of forest fire ash in the early storm runoff immediately after the forest fire. The event mean concentrations revealed that pollutant loads were associated with both forest fire ash and storm events in the immediate post-fire period. Therefore, this study indicated that water quality along a headwater stream is affected by short-term effects of large and intense forest fires and rainfall seasonality.

Spatial and temporal profiles of enteric viruses in the coastal waters of Tokyo Bay during and after a series of rainfall events

Recreational activities in coastal waters that are polluted by enteric viruses can result in gastroenteritis etc. In this study, the pollution profiles of enteric viruses were examined in the coastal area of Tokyo Bay, Japan, by collecting 57 water samples from three different depths (0.5 m, 3.0 m, and 5.0 m) during and after a series of heavy rainfall events. Vertically spatial and temporal changes in the concentrations of NoV genogroup I (GI) and genogroup II (GII), pepper mild mottle virus (PMMoV), and Aichi virus (AiV) were determined using quantitative reverse transcription-polymerase chain reaction, while those of the bacterial indicator, Escherichia coli, and F-specific RNA bacteriophages (FRNA phages) were monitored using culture methods. PMMoV was highly abundant (1.4 × 104–6.8 × 106 genome copies/L), whereas the concentrations of the other enteric viruses were relatively low (AiV, 1.3 × 102–2.9 × 104; GI, 2.9 × 10–5.6 × 103; GII, 2.5 × 10–1.2 × 104 genome copies/L). All of the viruses showed lower fluctuations in concentration than E. coli, which increased up to 460-fold after the rainfall event and then decreased over the subsequent two weeks. The maximum vertical difference in E. coli concentration was observed immediately after the rainfall. The E. coli reached the surface and then gradually spread down, whereas the virus concentrations exhibited few fluctuations due to the remaining effects of the previous combined sewer overflows. These findings indicate that viruses have a relatively long retention period over fecal indicator bacteria in this coastal area.

Temporal and spatial evolution of soil surface roughness on stony plots

Soil surface roughness (SSR) is widely recognized as an important factor influencing water erosion processes. On semiarid hillslopes with stony soils, rock fragments accumulate as the result of preferential erosion of fine materials, often creating a rough, rocky surface. A series of rainfall events were simulated on a stony plot (2 × 6.1 m) at three slope gradients (5 %, 12 %, and 20 %) and rock cover was measured. Surface elevations were sampled by terrestrial LiDAR at high resolutions. Roughness indices, including random roughness (RR), fractal dimension, crossover length, and generalized fractal dimension, were calculated from LiDAR points directly. Results showed: 1) SSR displayed an increasing trend as the rainfall simulation proceeded for all three slope treatments; 2) the steeper slope developed greater surface roughness; and 3) both the increase of surficial exposed rocks and the formation of erosional features, e.g., rills and depressions, contributed to the spatiotemporal variations of SSR. Results also showed that the fractal dimension was not a good indicator of soil surface roughness, but rather was an index of the form of the surface. Crossover length was a measure of roughness at a scale of a few millimeters, while random roughness was a measure of elevation variations on the scale of the length of the transect measured, and thus encompassed larger morphological features including rills. We also established a new method for multifractal analysis that characterized the heterogeneity of soil surface roughness. These results improve our understanding of the evolution of semiarid stony hillslopes and the dynamic feedback mechanism between erosion, surface morphology and hydraulics.

Reconstruction of July 2006 Heavy Rainfall Event Series in the Chungju Dam Watershed using a Regional Climate Model

The approach that uses a regional climate model (RCM) is emphasized to prevent and reduce extreme rainfall events that are becoming more powerful due to climate change. This paper presented a preliminary study to produce Probable Maximum Precipitation (PMP) or mega-disaster scenarios using RCMs in Korea. The heavy rainfall event in July 2006, during which the highest average precipitation during the Changma season was recorded, was reconstructed using the Weather Research and Forecasting (WRF) model for the Chungju Dam watershed. Even when considering errors in the rainfall observation data, the WRF model could successfully reconstruct the July 2006 heavy rainfall event. It is expected that more reliable results will be obtained if the flow observation data is used in combination with a runoff model. Keywords: Extreme Rainfall Event, Regional Climate Model, Changma, WRF

Spatio-Temporal Synthesis of Continuous Precipitation Series Using Vine Copulas

Long and continuous series of precipitation in a high temporal resolution are required for several purposes, namely, urban hydrological applications, design of flash flood control structures, etc. As data of the temporally required resolution is often available for short period, it is advantageous to develop a precipitation model to allow for the generation of long synthetic series. A stochastic model is applied for this purpose, involving an alternating renewal process (ARP) describing a system consisting of spells that can take two possible states: wet or dry. Stochastic generation of rainfall time series using ARP models is straight forward for single site simulation. The aim of this work is to present an extension of the model to spatio-temporal simulations. The proposed methodology combines an occurrence model to define in which locations rainfall events occur simultaneously with a multivariate copula to generate synthetic events. Rainfall series registered in different regions of Germany are used to develop and test the methodology. Results are compared with an existing method in which long independent time series of rainfall events are transformed to spatially dependent ones by permutation of their order. The proposed model shows to perform as a satisfactory extension of the ARP model for multiple sites simulations.

Numerical Simulation of Heavy Rainfall in August 2014 over Japan and Analysis of Its Sensitivity to Sea Surface Temperature

This study evaluated the performance of the Weather Research and Forecasting (WRF) model version 3.7 for simulating a series of rainfall events in August 2014 over Japan and investigated the impact of uncertainty in sea surface temperature (SST) on simulated rainfall in the record-high precipitation period. WRF simulations for the heavy rainfall were conducted for six different cases. The heavy rainfall events caused by typhoons and rain fronts were similarly accurately reproduced by three cases: the TQW_5km case with grid nudging for air temperature, humidity, and wind and with a horizontal resolution of 5 km; W_5km with wind nudging and 5-km resolution; and W_2.5km with wind nudging and 2.5-km resolution. Because the nudging for air temperature and humidity in TQW_5km suppresses the influence of SST change, and because W_2.5km requires larger computational load, W_5km was selected as the baseline case for a sensitivity analysis of SST. In the sensitivity analysis, SST around Japan was homogeneously changed by 1 K from the original SST data. The analysis showed that the SST increase led to a larger amount of precipitation in the study period in Japan, with the mean increase rate of precipitation being 13 ± 8% K−1. In addition, 99 percentile precipitation (100 mm d−1 in the baseline case) increased by 13% K−1 of SST warming. These results also indicate that an uncertainty of approximately 13% in the simulated heavy rainfall corresponds to an uncertainty of 1 K in SST data around Japan in the study period.

Relationship between man-made environment and slope stability: the case of 2014 rainfall events in the terraced landscape of the Liguria region (northwestern Italy)

ABSTRACTIn the autumn of 2014, a series of rainfall events affected several sectors of the Liguria region, triggering many shallow landslides and causing three casualties and severe structural damages. The most intensely unstable area covered 385 km2, in which more than 1600 landslides have been identified. After these events, an airborne Light Detection and Ranging survey was carried out. The survey yielded a high-resolution digital terrain model (DTM) and aerial images that provided a means of identifying and mapping all the occurred landslides. The distribution analysis of slope instabilities highlighted the link with various human activities. In fact, the majority of the detected landslides occurred in man-modified areas. Geospatial and statistical analyses provided the identification of three main anthropic factors: terraces, their level of maintenance and road network. Moreover, they quantified their role in landslide triggering. These factors were not analysed as separate elements, but as a continuous process, overlapping in time, in man-made influence on landscape. The identification of such factors is a key element for a correct behaviour characterization of this landscape towards extreme flash floods events.

Comparison of annual maximum rainfall events of modern rain gauge data (1961–2010) and Chukwooki data (1777–1910) in Seoul, Korea

Abstract In this study, the annual maximum rainfall event series were constructed and compared for both the modern flip-bucket type rainfall data, collected since 1961 (the modern data), and the old Chukwooki rainfall data, collected from 1777 to 1910 (the Chukwooki data). First, independent rainfall events were derived, by applying the same rainfall threshold of 2 mm and data collection time interval of 2 hours, to both the Chukwooki and the modern data. Annual maximum rainfall event series were then constructed, by applying Freund's bivariate exponential distribution annually. Finally, bivariate frequency analysis was done for the annual maximum rainfall event series constructed, by applying the bivariate logistic model to evaluate and quantify their characteristics. The results are in summary: (1) characteristics of the Chukwooki rainfall events and modern rainfall events are very similar to each other; (2) the annual maximum rainfall events of modern data are slightly larger than those of the Chukwooki data. The total rainfall depth per rainfall event for any given return period is thus estimated to be a little higher for the modern data than that of the Chukwooki data. However, based on the findings in this study, it could not be concluded that the rainfall characteristics have significantly changed during the last 200 years.

Transport-distance specific SOC distribution: Does it skew erosion induced C fluxes?

The net effect of soil erosion by water, as a sink or source of atmospheric carbon dioxide (CO2), is determined by the spatial (re-)distribution and stability of eroded soil organic carbon (SOC), and the dynamic replacement of eroded C by the production of new photosynthate. The depositional position of eroded SOC is a function of the transport distances of soil fractions where the SOC is stored. In theory, the transport distances of soil fractions are related to their settling velocities under given flow conditions. Yet, very few field investigations have been conducted to examine the actual movement of eroded soil fractions along hillslopes, let alone the re-distribution pattern of SOC fractions. Eroding sandy soils and sediment were sampled after a series of rainfall events along a slope on a freshly seeded cropland in Jutland, Denmark. All the soil samples were fractionated into five settling classes using a settling tube apparatus. The spatial distribution of soil settling classes shows a coarsening effect immediately below the eroding slope, followed by a fining trend at the slope tail. These findings support the validity of the conceptual model proposed by Starr et al. (Land Degrad Dev 11:83–91, 2000) to predict SOC redistribution patterns along hillslopes. The δ13C values of soil fractions were more positive at the footslope than on the slope shoulder or at the slope tail, suggesting enhanced decomposition rate of fresh SOC input at the footslope during or after erosion-induced transport. Pronounced CO2 emission rates at the slope tail also suggest a higher potential for decomposition of the eroded SOC after deposition. Overall, our results illustrate that immediate deposition of fast settling soil fractions and the associated SOC at footslopes, and potential CO2 emissions during or immediately after transport, must be appropriately accounted for in attempts to quantify the role of soil erosion in terrestrial C sequestration. A SOC erodibility parameter based on actual settling velocity distribution of eroded fractions is needed to better calibrate soil erosion models.

Earthquake-altered flooding hazard induced by damage to storm water systems

Major earthquakes can cause extensive transformations to the land underlying cities, leading to decreased capacity in natural and built drainage systems and, as a consequence, to Increased Flooding Hazard (IFH). This phenomenon causes some areas, which previously were not exposed to flooding, to have the potential to flood, and already flood-prone areas to likely experience increased flood depth during the next rainfall events. This scenario occurred in Christchurch city, New Zealand, after the 2010–2011 Canterbury Earthquake Sequence (CES). The IFH was observed in many urban areas during a series of rainfall events occurred in the years after the CES. This paper proposes a method for analysing and assessing to what extent the earthquake-induced damage to storm water pipelines and the consequent impacts on the connectivity and capacity levels of the pipeline storm water network could contribute to the IFH. A probabilistic analysis, through a Monte Carlo simulation, is suggested for the proposed method so that the uncertainty affecting several key parameters can be accounted for. The proposed probabilistic method for IFH was implemented as an additional module within a recently developed open-source simulation tool, OOFIMS. Results from the added OOFIMS module are presented in terms of maps and cumulative distribution functions of increased flood height and flooded area, impact metrics that can be useful for emergency managers and infrastructure owners. The effectiveness of the proposed method to assess earthquake-altered flooding hazard and the relative OOFIMS-added module are tested using Christchurch as a case-study.

The 2011 flood event in the Assiniboine River Basin: causes, assessment and damages

In spring and summer of 2011, the Assiniboine River and its tributaries underwent a flood of unprecedented proportions. It was, by several measures, including the computed return period and duration of flooding, a more extreme event than the 1997 Red River “Flood of the Century,” and may have been the most severe flood experienced in the history of Canada. The 2011 flood was the largest recorded in the over 100 years that flow records have been kept on the Assiniboine River. As reported in the 2013 Manitoba 2011 Flood Review Task Force Report, antecedent conditions, winter snowpack and summer rains combined to produce one of the largest floods ever experienced in Manitoba. The spring of 2011 was preceded by a very wet fall in 2010 when precipitation averaged between 150 and 200% of normal over a wide area in Manitoba and Saskatchewan This produced extremely high regional antecedent moisture conditions whereby soil moisture levels were 100 to 250% of the long-term average. Added to this, the winter snowpack was relatively high, varying geographically in the range of 90 to 130% of normal snow water content and lower-than-normal winter temperatures, causing significant frost penetration. The combination of high soil moisture in the previous summer/fall and deep frost penetration reduced the soil’s ability to absorb spring meltwater and increased runoff volumes contributing to the magnitude of the spring peak. The severe winter was followed by a series of rainfall events that resulted in unprecedented rainfall volumes in the region throughout the months of May, June and July, when some areas received rainfall amounts that were 350% of normal and produced up to two additional flood peaks at some locations. Flooding on the Assiniboine River near Brandon lasted for 120 days. Large volumes of floodwater were diverted to Lake Manitoba due to the use of the Portage Diversion to mitigate flooding downstream of the City of Portage La Prairie. This large volume of inflow resulted in a wind-effect-eliminated Lake Manitoba level that peaked at 249.1 m in late July of 2011, which was close to 1 m higher than flood stage and almost 1.5 m higher than the top of the desirable range. Peak outflow from Lake Manitoba into the Fairford River was almost double historic maxima and resulted in Lake St. Martin levels that peaked in July of 2011 at 245.55 m, which was more than 0.6 m higher than the historic peak that occurred under natural conditions in 1955, and 1.4 m above flood stage. The damage estimate for this flood event has not been finalized but is likely in excess of CAD $1 billion. The majority of flood damages were incurred by communities and infrastructure along major tributaries like the Souris and Qu’Appelle Rivers, and the communities on Lake Manitoba and Lake St. Martin. The Lake St. Martin First Nations community was virtually destroyed, and as of 2014 some First Nations people from the communities along Lake St. Martin are still evacuated and waiting for their homes and communities to be rebuilt. Temporary flood control works were largely effective in preventing infrastructure damage to communities along the Assiniboine River, such as Brandon, although the agricultural damage and damage to individual properties were extensive.

Evaluation of the Concept of Critical Rainfall Duration by Bivariate Frequency Analysis of Annual Maximum Independent Rainfall Event Series in Seoul, Korea

AbstractThis study investigated the possible problems of using the concept of critical rainfall duration (CRD) for estimating the design flood. The investigation focused on the evaluation of the difference between the annual maximum independent rainfall event series (IRES) and the annual maximum short-duration rainfall series in Seoul, Korea. Annual maximum IRES was determined by applying Freund’s bivariate exponential distribution, which were also used for frequency analysis by applying the Gumbel’s bivariate logistic model. Both the annual maximum IRES and the frequency analysis results were then compared with those of the conventional annual maximum fixed-duration rainfall series (FDRS). The results are summarized as follows. First, the mean rainfall intensity of the annual maximum IRES was found to be little changed regardless of their durations. Thus, most of the annual maximum short-duration rainfall series would be selected in the annual maximum IRES. Second, results of the bivariate frequency anal...

Modelling the interaction between flooding events and economic growth

Abstract. Socio-hydrology describes the interaction between the socio-economy and water. Recent models analyze the interplay of community risk-coping culture, flooding damage and economic growth (Di Baldassarre et al., 2013; Viglione et al., 2014). These models descriptively explain the feedbacks between socio-economic development and natural disasters like floods. Contrary to these descriptive models, our approach develops an optimization model, where the intertemporal decision of an economic agent interacts with the hydrological system. In order to build this first economic growth model describing the interaction between the consumption and investment decisions of an economic agent and the occurrence of flooding events, we transform an existing descriptive stochastic model into an optimal deterministic model. The intermediate step is to formulate and simulate a descriptive deterministic model. We develop a periodic water function to approximate the former discrete stochastic time series of rainfall events. Due to the non-autonomous exogenous periodic rainfall function the long-term path of consumption and investment will be periodic.

Factors affecting the spatial pattern of bedrock groundwater recharge at the hillslope scale

AbstractThe spatial patterns of groundwater recharge on hillslopes with a thin soil mantle overlying bedrock are poorly known. Complex interactions between vertical percolation of water through the soil, permeability contrasts between soil and bedrock and lateral redistribution of water result in large spatial variability of water moving into the bedrock. Here, we combine new measurements of saturated hydraulic conductivity of soil mantle and bedrock of the well‐studied Panola Mountain experimental hillslope with previously collected (sub)surface topography and soil depth data to quantify the factors affecting the spatial pattern of bedrock groundwater recharge.We use geostatistical characteristics of the measured permeability to generate spatial fields of saturated hydraulic conductivity for the entire hillslope. We perform simulations with a new conceptual model with these random fields and evaluate the resulting spatial distribution of groundwater recharge during individual rainstorms and series of rainfall events. Our simulations show that unsaturated drainage from soil into bedrock is the prevailing recharge mechanism and accounts for 60% of annual groundwater recharge. Therefore, soil depth is a major control on the groundwater recharge pattern through available storage capacity and controlling the size of vertical flux. The other 40% of recharge occurs during storms that feature transient saturation at the soil‐bedrock interface. Under these conditions, locations that can sustain increased subsurface saturation because of their topographical characteristics or those with high bedrock permeability will act as hotspots of groundwater recharge when they receive lateral flow. Copyright © 2015 John Wiley & Sons, Ltd.

A novel conceptual framework for long‐term leaching of autochthonous soil particles during transient flow

SummaryTo date, most of the experimental and modelling research on the mobilization mechanisms of autochthonous colloidal‐sized soil particles has focused on single rainfall events. However, in the field, natural events interrupted by rainless periods of various durations follow one another. Some studies have shown that the amount of leached particles is much affected by the chronology of rainfall events. In this paper, we propose a model to compute the amount of mobilized particles during the transient flow regime of long series of successive rainfall events. Particles that can be mobilized are assumed to be located at the surface of preferential flow paths in contact with flowing water. As water passes, particles are mobilized, thus uncovering new particles and renewing the stock of mobilizable particles. We also report experimental data on the impact of long series of rainfall events on particle mobilization in undisturbed calcareous soil columns under controlled conditions. The model provides a framework to understand the variation of mobilization observed during the transient flow regime of these rainfall events. These variations resulted in part from the interplay between macropore water content at the onset of rain and the chronology of previous rainfall events. Additionally, the model provides a feedback loop between particle mobilization and minute macropore structure modifications. Once coupled with a particle transport model this feature may be useful to model soil structure changes during long series of successive rainfall events.

Bivariate Frequency Analysis of Annual Maximum Rainfall Event Series in Seoul, Korea

AbstractThe return period of a rainfall event is estimated by the frequency analysis for a given rainfall duration. Thus, it is possible to derive different return periods with different rainfall durations for a given rainfall event. The longest derived return period is generally cited to represent the rainfall event. However, it is not clear if the longest derived return period is a representative measure of the given rainfall event. In this study, as a solution for this problem, a bivariate frequency analysis was introduced. As a first step, annual maximum rainfall events were selected by applying a bivariate exponential distribution. As an application, a total of 1,534 rainfall events observed in Seoul, Korea, over the last 46 years were analyzed. The annual maximum rainfall event series were then analyzed by applying a bivariate logistic model. The results were also compared with those from a conventional univariate frequency analysis. The findings of this study are summarized as follows: (1) the biva...

Early warning thresholds for partially saturated slopes in volcanic ashes

Rainfall-induced landslides in steep soil slopes of volcanic origin are a major threat to human lives and infrastructure. In the context of constructing early warning systems in regions where extensive data on landslide occurrences and associated rainfall are inexistent, physically-based tools offer the possibility to establish thresholds for measurable field quantities. In this paper, a combined finite element infinite slope model is presented to study the transient hydraulic response of volcanic ash slopes to a series of rainfall events and to estimate seasonal safety factors. Furthermore, analytical considerations of partially saturated infinite slopes are made to define capillary stress thresholds for a landslide early warning system.

독립 호우사상에 대한 이변량 강우빈도해석 및 강우-유출해석

본 연구에서는 독립 호우사상을 이용하여 이변량 빈도해석 및 유출해석을 수행하고, 이를 기존 단변량 빈도해석 결과와 비교 평가하였다. 본 연구는 규모가 다른 중랑천, 청계천, 우이천유역 등 세 유역에 대해 수행되었다. 유출모형으로 Clark 모형을 이용하였고, 유효우량은 SCS 방법을 적용하여 계산하였다. 강우의 시간분포 모형으로 교호블록 방법 및 Huff 방법을 적용하여 그 결과가 비교될 수 있도록 하였다. 그 결과를 정리하면 다음과 같다. (1) 연최대치 독립 호우사상계열의 이변량 빈도해석 결과, 지속기간이 짧은 경우에는 단변량 빈도해석 결과와의 차이가 매우 크나 지속기간이 길어짊에 따라서 그 차이가 현저히 줄어드는 것으로 나타났다. 아울러 지속기간이 짧은 경우, 단변량 빈도해석 결과가 이변량 빈도해석 결과보다 더욱 크게 나타났으나 특정 지속기간 이상부터 그 결과가 역전되는 것으로 나타났다. (2) 교호블록 시간분포 방법을 적용하는 경우가 Huff 방법을 적용한 경우보다 더욱 큰 첨두유출량을 발생시키는 것으로 나타났다. 아울러 교호블록 방법을 적용하는 경우에는 강우 지속기간의 증가에 따라서 첨두유출량이 점차 증가하는 것으로 나타났으나, 강우 지속기간이 대략 24시간 정도 되었을 때 그 값이 거의 수렴하는 것으로 나타났다. (3) 중랑천 유역에 대해 Huff 방법을 적용하여 유출해석을 수행한 결과에서는 이변량 설계강우를 적용한 경우가 단변량 설계강우를 적용한 경우보다 더욱 큰 홍수량을 발생시키는 것으로 나타났다. 반면에 청계천 및 우이천 유역의 경우에는 이변량 설계강우를 적용한 경우보다 단변량 설계강우를 적용한 경우의 홍수량이 다소 큰 것으로 나타났다. 그러나 교호블록 방법을 적용한 경우에서는 모든 유역에 대해 이변량 설계강우를 적용한 경우가 단변량 설계강우를 적용한 경우보다 큰 홍수량을 발생시키는 것으로 나타났다. In this study, the bivariate frequency analysis of the independent annual rainfall event series was done to be used for the runoff analysis, whose results were also compared with those from the conventional univariate frequency analysis. This study was applied to three differently-sized basins such as the Joongryang Stream, Chunggye Stream, and Ooyi Stream. The Clark model was used as the runoff model, and the SCS method was applied for the calculation of the effective rainfall. The alternating block method and the Huff method were considered to be compared for the temporal distribution of rainfall event. Summarizing the results are as follows. (1) The difference between the univariate and bivariate frequency analysis results were large when the rainfall duration was short, but significantly decreased as the rainfall duration increased. The univariate frequency analysis results were bigger when the rainfall duration was short, but smaller in opposite case. (2) The peak flow derived by applying the alternating block method was bigger than that by the Huff method. Also, the peak flow when applying the alternating block method increased as the rainfall duration increased, but converged smoothly around the rainfall duration of 24 hours. (3) For the Joongryang Stream, when applying the Huff method, the peak flow derived for the bivariate frequency analysis was bigger than that for the univariate case, but for the other two basins, the results were opposite. When applying the alternating block method, the results were consistent for all three basins that the peak flow derived by applying the bivariate frequency analysis was bigger than those by the univariate frequency analysis.