• The Korean Journal of Applied Statistics
• /
• v.35 no.1
• /
• pp.119-129
• /
• 2022

In this paper, we predict the earthquake magnitudes which were recently occurred in Gyeongju and Pohang, using statistical methods based on historical data. For this purpose, we use the five-year block maximum data of 1392~1771 period, which has a relatively high annual density, among the historical earthquake magnitude data of the Chosun Dynasty. Then, we present the prediction and analysis of earthquake magnitudes for the return level over return period in the Chosun Dynasty using the extreme value theory based on the distribution of generalized extreme values (GEV). We use maximum likelihood estimation (MLE) and L-moments estimation for parameters of GEV distribution. In particular, this study also demonstrates via the goodness-of-fit tests that the GEV distribution can be an appropriate analytical model for these historical earthquake magnitude data.

• Proceedings of the Korean Society of Agricultural Engineers Conference
• /
• 2005.10a
• /
• pp.284-287
• /
• 2005

The objective of this study is to check into variation trends of design rainfall according to change of the number of years for observed data. To make comparative study of the relation between design rainfall and recorded year, this study was used maximum rainfall for 24-hr consecutive duration at Gangneung, Seoul, Incheon, Chupungnyeong, Pohang, Daegu, Jeonju, Ulsan, Gwangju, Busan, Mokpo and Yeosu rainfall stations. The tests for Independence, Homogeneity and detection of outliers were used Wald-Wolfowitz's test, Mann-Whitney's test and Grubbs and Beck test respectively. To select appopriate distribution, the distribution of genaralized pareto(GPA), generalized extreme value(GEV), generalized logistic(GLO), lognormal and pearson type 3 distribution is judged by L-moment ratio diagram and Kolmogorov-Smirnov (K-S) test. Design rainfall was estimated by at-site frequency analysis using L-moments and Generalized extreme value(GEV) distribution according to change of the number of years for observed data. Through the comparative analysis for design rainfall induced by L-moments and GEV distribution, relationship between design rainfall and recorded year is provided.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.5
• /
• pp.482-490
• /
• 2008

For a coastal or harbor structure design, one of the most important environmental factors is the appropriate extreme highest tide level condition. Especially, the information of extreme highest tide level distribution is essential for reliability design. In this paper, 23 set of extreme highest tide level data obtained from National Oceanographic Research Institute(NORI) were analyzed for extreme highest tide levels. The probability distributions considered in this research were Generalized Extreme Value(GEV), Gumbel, and Weibull distribution. For each of these distributions, three parameter estimation methods, i.e. the method of moments, maximum likelihood and probability weighted moments, were applied. Chi-square and Kolmogorov-Smirnov goodness-offit tests were performed, and the assumed distribution was accepted at the confidence level 95%. Gumbel distribution which best fits to the 22 tidal station was selected as the most probable parent distribution, and optimally estimated parameters and extreme highest tide level with various return periods were presented. The extreme values of Incheon, Cheju, Yeosu, Pusan, and Mukho, which estimated by Shim et al.(1992) are lower than that of this result.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.21 no.3
• /
• pp.253-258
• /
• 2015

The purpose of this study is to find suitable probability distribution function of complex distribution data like multimodal. Normal distribution is broadly used to assume probability distribution function. However, complex distribution data like multimodal are very hard to be estimated by using normal distribution function only, and there might be errors when other distribution functions including normal distribution function are used. In this study, we experimented to find fit probability distribution function in multimodal area, by using AIS(Automatic Identification System) observation data gathered in Mokpo port for a year of 2013. By using chi-squared statistic, gaussian mixture model(GMM) is the fittest model rather than other distribution functions, such as extreme value, generalized extreme value, logistic, and normal distribution. GMM was found to the fit model regard to multimodal data of maritime traffic flow distribution. Probability density function for collision probability and traffic flow distribution will be calculated much precisely in the future.

• Journal of Korean Society of Water and Wastewater
• /
• v.30 no.1
• /
• pp.41-49
• /
• 2016

In this study, rainfall characteristics with stationary and non-stationary perspectives were analyzed using generalized extreme value (GEV) distribution and Gumbel distribution models with rainfall data collected in major cities of Korea to reevaluate the return period of sewer flooding in those cities. As a result, the probable rainfall for GEV and Gumbel distribution in non-stationary state both increased with time(t), compared to the stationary probable rainfall. Considering the reliability of ${\xi}_1$, a variable reflecting the increase of storm events due to climate change, the reliability of the rainfall duration for Seoul, Daegu, and Gwangju in the GEV distribution was over 90%, indicating that the probability of rainfall increase was high. As for the Gumbel distribution, Wonju, Daegu, and Gwangju showed the higher reliability while Daejeon showed the lower reliability than the other cities. In addition, application of the maximum annual rainfall change rate (${\xi}_1{\cdot}t$) to the location parameter made possible the prediction of return period by time, therefore leading to the evaluation of design recurrence interval.

• Proceedings of the Korean Society of Agricultural Engineers Conference
• /
• 2002.10a
• /
• pp.229-232
• /
• 2002

This study was conducted to estimate the design flood by the determination of best fitting order of LH-moments of the annual maximum series at six and nine watersheds in Korea and Australia, respectively. Adequacy for flood flow data was confirmed by the tests of independence, homogeneity, and outliers. Gumbel (GUM), Generalized Extreme Value (GEV), Generalized Pareto (GPA), and Generalized Logistic (GLO) distributions were applied to get the best fitting frequency distribution for flood flow data. Theoretical bases of L, L1, L2, L3 and L4-moments were derived to estimate the parameters of 4 distributions. L, L1, L2, L3 and L4-moment ratio diagrams (LH-moments ratio diagram) were developed in this study.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.259-267
• /
• 2009

Generalized Extreme Value (GEV) distribution is recommended for flood frequency and extreme rainfall distribution in many country. L-moment method is the most common estimation procedure for the GEV distribution. In this study, the relationships between the cross-site correlations between extreme events and the cross-correlation of estimators of L-moment ratios (L-moment Coefficient of Variation (L-CV) and L-moment Coefficient of Skewness (L-CS)) for data generated from GEV distribution were derived by Monte Carlo simulation. Those relationships were fit to the simple power function. In this Monte Carlo simulation, GEV+ distribution were employed wherein unrealistic negative values were excluded. The simple power models provide accurate description of the relationships between cross-correlation of data and cross-correlation of L-moment ratios. Estimated parameters and accuracies of the power functions were reported for different GEV distribution parameters combinations. Moreover, this study provided a description about regional regression approach using Generalized Least Square (GLS) regression method which require the cross-site correlation among L-moment estimators. The relationships derived in this study allow regional GLS regression analyses of both L-CV and L-CS estimators that correctly incorporate the cross-correlation among GEV L-moment estimators.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.54 no.6
• /
• pp.133-142
• /
• 2012

The objective of this study was to correct the bias of the Representative Concentration Pathways (RCP)-based future precipitation data using a quantile mapping method. This method was adopted to correct extreme values because it was designed to adjust simulated data using probability distribution function. The Generalized Extreme Value (GEV) distribution was used to fit distribution for precipitation data obtained from the Korea Meteorological Administration (KMA). The resolutions of precipitation data was 12.5 km in space and 3-hour in time. As the results of bias correction over the past 30 years (1976~2005), the annual precipitation was increased 16.3 % overall. And the results for 90 years (divided into 2011~2040, 2041~2070, 2071~2100) were that the future annual precipitation were increased 8.8 %, 9.6 %, 11.3 % respectively. It also had stronger correction effects on high value than low value. It was concluded that a quantile mapping appeared a good method of correcting extreme value.

• Communications for Statistical Applications and Methods
• /
• v.25 no.1
• /
• pp.15-27
• /
• 2018

Parametric method of flood frequency analysis involves fitting of a probability distribution to observed flood data. When record length at a given site is relatively shorter and hard to apply the asymptotic theory, an alternative distribution to the generalized extreme value (GEV) distribution is often used. In this study, we consider the beta-P distribution (BPD) as an alternative to the GEV and other well-known distributions for modeling extreme events of small or moderate samples as well as highly skewed or heavy tailed data. The L-moments ratio diagram shows that special cases of the BPD include the generalized logistic, three-parameter log-normal, and GEV distributions. To estimate the parameters in the distribution, the method of moments, L-moments, and maximum likelihood estimation methods are considered. A Monte-Carlo study is then conducted to compare these three estimation methods. Our result suggests that the L-moments estimator works better than the other estimators for this model of small or moderate samples. Two applications to the annual maximum stream flow of Colorado and the rainfall data from cloud seeding experiments in Southern Florida are reported to show the usefulness of the BPD for modeling hydrologic events. In these examples, BPD turns out to work better than $beta-{\kappa}$, Gumbel, and GEV distributions.

• Wind and Structures
• /
• v.26 no.3
• /
• pp.129-146
• /
• 2018

Extreme wind speed analysis has been carried out conventionally by assuming the extreme series data is stationary. However, time-varying trends of the extreme wind speed series could be detected at many surface meteorological stations in China. Two main reasons, exposure change and climate change, were provided to explain the temporal trends of daily maximum wind speed and annual maximum wind speed series data, recorded at Hangzhou (China) meteorological station. After making a correction on wind speed series for time varying exposure, it is necessary to perform non-stationary statistical modeling on the corrected extreme wind speed data series in addition to the classical extreme value analysis. The generalized extreme value (GEV) distribution with time-dependent location and scale parameters was selected as a non-stationary model to describe the corrected extreme wind speed series. The obtained non-stationary extreme value models were then used to estimate the non-stationary extreme wind speed quantiles with various mean recurrence intervals (MRIs) considering changing climate, and compared to the corresponding stationary ones with various MRIs for the Hangzhou area in China. The results indicate that the non-stationary property or dependence of extreme wind speed data should be carefully evaluated and reflected in the determination of design wind speeds.