• Journal of Korea Water Resources Association
• /
• v.45 no.8
• /
• pp.827-837
• /
• 2012

The estimation of the rainfall quantile is of great importance in designing hydrologic structures. Conventionally, the rainfall quantile is estimated by univariate frequency analysis with an appropriate probability distribution. There is a limitation in which duration of rainfall is restrictive. To overcome this limitation, bivariate frequency analysis by using 3 copula models is performed in this study. Annual maximum rainfall events in 5 stations are used for frequency analysis and rainfall depth and duration are used as random variables. Gumbel (GUM), generalized logistic (GLO) distributions are applied for rainfall depth and generalized extreme value (GEV), GUM, GLO distributions are applied for rainfall duration. Copula models used in this study are Frank, Joe, and Gumbel-Hougaard models. Maximum pseudo-likelihood estimation method is used to estimate the parameter of copula, and the method of probability weighted moments is used to estimate the parameters of marginal distributions. Rainfall quantile from this procedure is compared with various marginal distributions and copula models. As a result, in change of marginal distribution, distribution of duration does not significantly affect on rainfall quantile. There are slight differences depending on the distribution of rainfall depth. In the case which the marginal distribution of rainfall depth is GUM, there is more significantly increasing along the return period than GLO. Comparing with rainfall quantiles from each copula model, Joe and Gumbel-Hougaard models show similar trend while Frank model shows rapidly increasing trend with increment of return period.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.46 no.5
• /
• pp.27-39
• /
• 2004

In the first part of this study, five homogeneous regions in view of topographical and geographically homogeneous aspects except Jeju and Ulreung islands in Korea were accomplished by K-means clustering method. A total of 57 rain gauges were used for the regional frequency analysis with minimum rainfall series for the consecutive durations. Generalized Extreme Value distribution was confirmed as an optimal one among applied distributions. Drought rainfalls following the return periods were estimated by at-site and regional frequency analysis using L-moments method. It was confirmed that the design drought rainfalls estimated by the regional frequency analysis were shown to be more appropriate than those by the at-site frequency analysis. In the second part of this study, LH-moment ratio diagram and the Kolmogorov-Smirnov test on the Gumbel (GUM), Generalized Extreme Value (GEV), Generalized Logistic (GLO) and Generalized Pareto (GPA) distributions were accomplished to get optimal probability distribution. Design drought rainfalls were estimated by both at-site and regional frequency analysis using LH-moments and GEV distribution, which was confirmed as an optimal one among applied distributions. Design rainfalls were estimated by at-site and regional frequency analysis using LH-moments, the observed and simulated data resulted from Monte Carlotechniques. Design drought rainfalls derived by regional frequency analysis using L1, L2, L3 and L4-moments (LH-moments) method have shown higher reliability than those of at-site frequency analysis in view of RRMSE (Relative Root-Mean-Square Error), RBIAS (Relative Bias) and RR (Relative Reduction) for the estimated design drought rainfalls. Relative efficiency were calculated for the judgment of relative merits and demerits for the design drought rainfalls derived by regional frequency analysis using L-moments and L1, L2, L3 and L4-moments applied in the first report and second report of this study, respectively. Consequently, design drought rainfalls derived by regional frequency analysis using L-moments were shown as more reliable than those using LH-moments. Finally, design drought rainfalls for the classified five homogeneous regions following the various consecutive durations were derived by regional frequency analysis using L-moments, which was confirmed as a more reliable method through this study. Maps for the design drought rainfalls for the classified five homogeneous regions following the various consecutive durations were accomplished by the method of inverse distance weight and Arc-View, which is one of GIS techniques.

• The Korean Journal of Applied Statistics
• /
• v.26 no.3
• /
• pp.483-494
• /
• 2013

The importance of financial risk management has been highlighted after several recent incidences of global financial crisis. One of the issues in financial risk management is how to measure the risk; currently, the most widely used risk measure is the Value at Risk(VaR). We can consider to estimate VaR using extreme value theory if the financial data have heavy tails as the recent market trend. In this paper, we study estimations of VaR using Peaks over Threshold(POT), which is a common method of modeling fat-tailed data using extreme value theory. To use POT, we first estimate parameters of the Generalized Pareto Distribution(GPD). Here, we compare three different methods of estimating parameters of GPD by comparing the performance of the estimated VaR based on KOSPI 5 minute-data. In addition, we simulate data from normal inverse Gaussian distributions and examine two parameter estimation methods of GPD. We find that the recent methods of parameter estimation of GPD work better than the maximum likelihood estimation when the kurtosis of the return distribution of KOSPI is very high and the simulation experiment shows similar results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.39 no.1
• /
• pp.165-174
• /
• 2019

As a nonstationarity is observed in hydrological data, various studies on nonstationary frequency analysis for hydraulic structure design have been actively conducted. Although the inherent diversity in the atmosphere-ocean system is known to be related to the nonstationary phenomena, a nonstationary frequency analysis is generally performed based on the linear trend. In this study, a nonstationary frequency analysis was performed using climate indices as covariates to consider the climate variability and the long-term trend of the extreme rainfall. For 11 weather stations where the trend was detected, the long-term trend within the annual maximum rainfall data was extracted using the ensemble empirical mode decomposition. Then the correlation between the extracted data and various climate indices was analyzed. As a result, autumn-averaged AMM, autumn-averaged AMO, and summer-averaged NINO4 in the previous year significantly influenced the long-term trend of the annual maximum rainfall data at almost all stations. The selected seasonal climate indices were applied to the generalized extreme value (GEV) model and the best model was selected using the AIC. Using the model diagnosis for the selected model and the nonstationary GEV model with the linear trend, we identified that the selected model could compensate the underestimation of the rainfall quantiles.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.37 no.6
• /
• pp.103-111
• /
• 1995

Drought index calculation based on the principal hydrological parameters, such as rainfall and reservoir storage, can estimate the duration and intensity of drought in irrigation reservoirs. It is difficult to build up a drought criteria since the conditions change variously by the reliability of rainfall. Because of the increasing water demands, it is urgent to prepare a generalized positive countermeasure to overcome drought. Water demands can at calculated but the estimation of drought characteristics, and the effective water management method can be established. The purpose of this study is to obtain a drought index and build up a data-base on the reservoir basins for establishing the fundamental hydrological data-base. This Index can observe the behavior of the WSI(Water Supply Index) and the component indices. The results summarized through this study are as follows. 1. WSI value of zero does not correspond to 100% in average due to the skewness in the probability distributions. 2. WSI is not a linear index; that is, given change in terms of water volume or percentage of average does not result in a proportional change on the WSI scale. 3. WSI is not always between the reservoir and the rainfall index in magnitude. This is only true if the component indices are of opposite sign. If they are of the same sign, the SWSI will often have a mangitude greater than either of the component indices. This is easily understood, because the concurrence of extreme values of the same sign for the two components is rarer than the occurrence of extreme values for either of the two components individually.

• The Korean Journal of Applied Statistics
• /
• v.23 no.3
• /
• pp.471-485
• /
• 2010

VaR is used extensively as a tool for risk management by financial institutions. For convenience, the normal distribution is usually assumed for the measurement of VaR, but recently the method using extreme value theory is attracted for more accurate VaR estimation. So far, GEV and GPD models are used for probability models of EVT for the VaR estimation. In this paper, the PP model is suggested for improved VaR estimation as compared to the traditonal EV models such as GEV and GPD models. In view of the stochastic process, the PP model is regarded as a generalized model which include GEV and GPD models. In the empirical analysis, the PP model is shown to be superior to GEV and GPD models for the performance of VaR estimation.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.162-162
• /
• 2017

최근 우리나라는 전 지구적인 기후변화로 인하여 집중호우 및 돌발 홍수와 같은 극치 사상들이 증가하고 있는 추세이며, 이에 대한 분석을 위해 극치 분포를 이용한 수문통계적 특성에 대한 접근이 주로 이루어지고 있다. 이를 위해서는 충분한 수의 자료가 필요하나 우리나라 강우자료는 지점별로 자료 보유 년 수가 비교적 많지 않기 때문에, 이러한 문제를 극복하기 위하여 하나의 지역, 즉 주어진 지점을 포함하여 수문학적으로 동일한 조건을 만족하는 주변 지점의 자료를 모두 포함하여 빈도해석을 실시하는 지역빈도해석이 필요하다. 따라서 본 연구에서는 지역빈도해석과 두 개의 형상매개변수를 포함하여 다양한 극치 수문통계특성을 나타낼 수 있다고 알려진 Burr XII 분포를 이용하여 우리나라 강우자료에 대한 그 적용성을 살펴보았다. 이를 위해 군집분석을 통한 강우지점의 지역화 과정을 거치고 분류된 지역을 L-moment ratio diagram에 도시하여, Burr XII 분포 영역 내 포함여부를 통해 Burr XII 분포의 적합도를 도시적으로 살펴보고, Hosking and Wallis (1997)이 제안한 적합성 척도($^{IST}$)를 통한 적합성 여부를 판별하였다. 또한 우리나라 강우자료에 비교적 적합하다고 알려진 분포인 generalized extreme value, generalized logistic, Gumbel 분포와의 비교를 위해, 전체 지역에 대하여 재현기간에 따른 상대편의 (relative bias)와 상대평균제곱근오차 (relative root mean square error)를 산정하여 Burr XII 분포형의 적용 가능성을 살펴보았다.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.6
• /
• pp.1056-1063
• /
• 2015

Waves can be expressed in terms of a spectrum; that is, the energy density distribution of a representative wave can be determined using statistical analysis. The JONSWAP, PM and BM spectra have been widely used for the specific target wave data set during storms. In this case, the extracted wave data are usually discontinuous and independent and cover a very short period of the total data-recording period. Previous studies on the continuous wave spectrum have focused on wave deformation in shallow water conditions and cannot be generalized for deep water conditions. In this study, the Generalized Extreme Value (GEV) function is proposed as a more-optimal function for the fitting of the continuous wave spectral shape based on long-term monitored point wave data in deep waters. The GEV function was found to be able to accurately reproduce the wave spectral shape, except for discontinuous waves of greater than 4 m in height.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.200-200
• /
• 2022

본 연구에서는 총 강우량과 강우강도을 고려한 이변수 분석으로 연최대 호우사상을 선별하고, 두 변수를 Copula 함수로 결합하여 최적의 모델조합을 찾는 확률호우사상 산정 방법론을 제시하였다. 국내 69개 관측소의 2020년까지의 관측 자료를 대상으로 1mm 이하의 강우는 제거한 뒤, IETD(Inter-Event Time Definition) 12시간을 기준으로 강우자료를 독립적인 호우사상으로 분리하였다. 호우사상의 여러 특성 중 양의 상관관계를 갖는 총 강우량과 강우강도를 변수로 선택해 이변수 지수분포에 대입하였고, 각 지점의 연최대 호우사상 시계열을 생성하였다. 2변수 지수분포의 매개변수는 전체 기간과 연도별로 나누어 추정해 본 결과 연도별 변동성이 큰 것을 확인해 연도별 추정 방식을 선택하였다. 연최대 강우사상 시계열의 총 강우량과 강우강도는 극한 강우에 적용하는 확률분포형 중 Lognarmal, Gamma, Gumbel, GEV(Generalized Extreme Value), GPD(Generalized Pareto Distribution) 5가지를 사용하여 각각 CDF(Cumulative distribution Function) 값을 추정하였다. 계산된 CDF 값은 3가지 Copula 모형으로 결합해 joint CDF 값을 산출하였다. 총 75개의 모델조합 중 최적 모델을 찾기 위해 CVM(Cramer-von-Mises) 적합도 검정을 시행하였다. CVM의 통계량 Sn 값이 가장 작은 모델조합을 해당 지점의 최적 모델조합으로 선정하였다.

• Wind and Structures
• /
• v.38 no.4
• /
• pp.309-323
• /
• 2024

The gust factor and turbulence intensity are two crucial parameters that characterize the properties of turbulence. In tropical cyclones (TCs), these parameters exhibit significant variability, yet there is a lack of established formulas to account for their probabilistic characteristics with consideration of their inherent connection. On this condition, a probabilistic analysis of gust factors and turbulence intensities of TCs is conducted based on fourteen sets of wind data collected at the Sutong Cable-stayed Bridge site. Initially, the turbulence intensities and gust factors of recorded data are computed, followed by an analysis of their probability densities across different ranges categorized by mean wind speed. The Gaussian, lognormal, and generalized extreme value (GEV) distributions are employed to fit the measured probability densities, with subsequent evaluation of their effectiveness. The Gumbel distribution, which is a specific instance of the GEV distribution, has been identified as an optimal choice for probabilistic characterizations of turbulence intensity and gust factor in TCs. The corresponding empirical models are then established through curve fitting. By utilizing the Gumbel distribution as a template, the nexus between the probability density functions of turbulence intensity and gust factor is built, leading to the development of a generalized probabilistic model that statistically describe turbulence intensity and gust factor in TCs. Finally, these empirical models are validated using measured data and compared with suggestions recommended by specifications.