Beta-�� distribution and its application to hydrologic events

The beta-κ distribution is a distinct case of the generalized beta distribution of the second kind. In previous studies, beta-p and beta-κ distributions have played important roles in representing extreme events, and thus, the present paper uses the beta-κ distribution. Further, this paper uses the method of moments and the method of L-moments to estimate the parameters from the beta-κ distribution, and to demonstrate the performance of the proposed model, the paper presents a simulation study using three estimation methods (including the maximum likelihood estimation method) and beta-κ and non beta-κ samples. In addition, this paper evaluates the performance of the beta-κ distribution by employing two widely used extreme value distributions (i.e., the GEV and Gumbel distributions) and two sets of actual data on extreme events.     

Long-term simulations of the sulfur concentrations over the China, Japan and Korea: A model comparison study

Three comprehensive acid deposition models were used to simulate the sulfur concentrations over northeast Asia over the period covering entire year of 2002, and discussed the aggregated uncertainties and discrepancies of the three models. The participating models are from the countries participating in the project of Longrange Transboundary Air Pollutants in Northeast Asia (LTP): China, Japan and Korea. The Eulerian Model-3/CMAQ (by China), Regional Air Quality Model (RAQM, by Japan), and Comprehensive Acid Deposition Model (CADM, by Korea) were employed by each country with common emissions data established by the administrative agencies of China, Japan and Korea. The episodic simulation results between 1 to 15, March 2002 are also presented, during which aircraft measurements were carried out over the Yellow sea. The episodic results show both a wide short-term variability in simulations against measurements, and maximum concentration differences of 3～5 times among the three models, requiring that further attention before confidence among the three models can be claimed for short-term simulations. However, the year-long cumulative simulations showed almost the same general features, with lower aggregated uncertainties between the three models, produced by the long term integration over northeast Asia.     

Representation of tropical storms in the northwestern pacific by the Modern-Era Retrospective analysis for research and applications

This study examines the tropical storms simulated in the Modern-Era Retrospective analysis for Research and Applications (MERRA) global atmospheric reanalysis for the recent 12 years (1998–2009), focusing on the tropical storm activity over the Northwestern Pacific. For validation, the International Best Track Archive for Climate Stewardship (IBTrACS) dataset is used as an observational counterpart. Climatological-mean features of the tropical storm genesis, tracks and their maximum intensity are the primary interests in this study. Regarding the genesis location of tropical storms, MERRA is reasonable in resolving major development regions over the South China Sea and the Northwestern Pacific close to the Philippines. The seasonal variation of the number of storms is also reproduced in a realistic way in MERRA, with peak values occurring from July to September. In addition, MERRA tends to reproduce the observed interannual variation of the number of tropical storms during the 12-years, though with a limited accuracy. The simulated paths toward higher latitudes are also reasonable in MERRA, where the reanalysis corresponds well with the observations in resolving frequent paths of westward moving storms and recurving storms toward the northeast. Regarding the intensity, MERRA captures the linear relationship between the minimum center pressure and the maximum wind speed near the surface at the maximum development. Some discrepancies from the observed features are found in the reanalysis, such as less frequent development of storms over the South China Sea and less frequent paths over this region. The reanalysis also does not attain the observed maximum intensity for the resolved tropical storms, particularly underestimating the center pressure. These deficiencies are likely related to limitations in the horizontal resolution and the parameterized physics of the data assimilation system.     

An algorithm for fast 3D inversion of surface electrical resistivity tomography data: application on imaging buried antiquities

In this work a new algorithm for the fast and efficient 3D inversion of conventional 2D surface electrical resistivity tomography lines is presented. The proposed approach lies on the assumption that for every surface measurement there is a large number of 3D parameters with very small absolute Jacobian matrix values, which can be excluded in advance from the Jacobian matrix calculation, as they do not contribute significant information in the inversion procedure. A sensitivity analysis for both homogeneous and inhomogeneous earth models showed that each measurement has a specific region of influence, which can be limited to parameters in a critical rectangular prism volume. Application of the proposed algorithm accelerated almost three times the Jacobian (sensitivity) matrix calculation for the data sets tested in this work. Moreover, application of the least squares regression iterative inversion technique, resulted in a new 3D resistivity inversion algorithm more than 2.7 times faster and with computer memory requirements less than half compared to the original algorithm. The efficiency and accuracy of the algorithm was verified using synthetic models representing typical archaeological structures, as well as field data collected from two archaeological sites in Greece, employing different electrode configurations. The applicability of the presented approach is demonstrated for archaeological investigations and the basic idea of the proposed algorithm can be easily extended for the inversion of other geophysical data.     

A DEM-based evaluation of potential flood risk to enhance decision support system for safe evacuation

A practical, DEM-based practical method is proposed to enhance flood risk management in fluvial areas by quantifying relative risk as a function of vulnerability to inland and evacuation difficulty. Both measures are based mainly on the topography of the region, so the method does not require detailed data on the physical characteristics of the land. First, we use the deterministic 8-node method on a digital elevation map (DEM) to trace storm waterways. Second, we repeat the process on a reversed DEM to trace evacuation routes that avoid the waterways and zones dangerously close to the rivers. Finally, on the basis of such two flow lines of evacuee and storm water, we proposed the protocol to evaluate the flood risk at every point on the map taking into account both the minimum time required for floodwater to arrive and duration of an evacuation from that location. The time that must be allocated for safe evacuation is defined as the potential flood risk of evacuation (PFRE). The method is demonstrated on a fluvial area of the Kaki River in Nagaoka city, Japan. In addition, we illustrated the application of the PFRE map to divide the region into areas of greater or lesser evacuation urgency.     

Effect of stress history on CPT and DMT results in sand

In order to investigate the effect of stress history on in-situ test results in granular sediments, a series of CPTs and DMTs are performed on Busan sand prepared in the calibration chamber. K_D is found to be the most sensitive to the stress history among CPT and DMT measurements. E_D and q_c are observed to be similarly affected by the stress history and, therefore, the E_D–q_c relation appears to be almost independent of the stress history. The K_D–D_R relation established without considering the stress history is likely to overestimate the relative density of OC sand. It is shown that the existence of the pre-stress of the granular sediment can be indirectly recognized by an estimation of the relative density larger than 100% when using the K_D–σ_v′–D_R relation suggested for NC sand. Although q_c/σ_v′–K_D/K₀ and E_D/σ_v′–K_D/K₀ relations are heavily influenced by the stress history, q_c/σ_m′–K_D/K₀ and E_D/σ_m′–K_D/K₀ relations are observed to be independent of the stress history. Based on these relations, charts to evaluate the K₀ value from q_c and/or DMT indices are developed for both NC and OC sands. The design chart based on E_D/σ_m′–K_D/K₀ and E_D/σ_v′–K_D/K₀ relations is expected to be practically useful as the usage of this chart requires only DMT indices. The developed design charts are applicable to Busan sand but different sets of equations and charts may be developed for other sands.     

A biogeochemical orientation survey in the Moisan gold-mineralized area, Haenam district in Korea

A biogeochemical orientation survey was carried out in the vicinity of an epithermal Au deposit in the Moisan Au–Ag mineralized area, Haenam district in Korea. The Au–Ag bearing quartz veins of the mine occur as narrow open-space fillings within Cretaceous silicic pyroclastics. The vein minerals consist mainly of quartz, sericite, pyrite, chalcopyrite, and galena, with some electrum and argentite. The main objectives of this study were to study the geochemical characteristics of rocks, soils and plants in this area, to investigate the spatial relationship between Au and associated elements in rock–soil–plant system, and to evaluate the applicability of biogeochemical prospecting for Au vein occurrences in Korea. Samples of rocks and soils, and leaves of three plant species (Japanese red pine — P. densiflora, oriental white oak — Q. aliena, Japanese mallotus — M. japonicus) were collected from the target mineralized area and control barren locations, and analyzed for trace elements by instrumental neutron activation analysis. Sampling lines were composed of one slope line which is almost parallel to the mineralized quartz-veins, and four transect lines spaced 100 m apart across the veins at 20 m sampling intervals. From the multi-element data of rock samples (n = 9), high values of Au (maximum 2030 ppb) are spatially related to Au–quartz veins. Soil samples (n = 61) collected from five sampling lines show higher values of Au (24–825 ppb) whereas soil samples from the control locations have lower values of Au (below 25 ppb). Many plant species collected from the vicinity of the veins have high Au contents compared with those at the control locations, but the ranges of Au values are variable among plant species. In a total of 128 samples of plant leaves, Q. aliena yielded Au values of 0.4 to 6.9 ppb, and M. japonicus 0.9 to 4.1 ppb. Gold contents in P. densiflora ranged from 0.1 to 5.6 ppb. Plant leaves from control areas show less than 1.6 ppb Au. The biological absorption coefficient (BAC) of Au in plants decreases in the order of Q. aliena > M. japonicus > P. densiflora. Based on the results of the study, Q. aliena appeared to be the best sampling media for biogeochemical prospecting of Au in the study area.     

Comment on ‘Clinoform quantification for assessing the effects of external forcing on continental margin development’

The application of geometric modeling to shelf‐margin stratigraphy has the potential to constrain interpretations of external forcings on margin development. Here we apply such a model to the Ebro margin in order to complement the analysis of Kertznus & Kneller (2009) . Our results suggest that increased mass wasting in the Pleistocene was unlikely to have been a factor in the observed long‐term shelf‐edge trajectory, and that the trajectory can be explained by the interaction of sediment flux, relative sea‐level rise, and basin shape.     

Numerical analysis and validation of weakly nonlinear ship motions and structural loads on a modern containership

This paper aims to validate a numerical seakeeping code based on a 3D Rankine panel method by comparing its results with experimental data. Particularly, the motion response and hull-girder loads on a real modern ship, a 6500 TEU containership, are considered in this validation study. The method of solution is a 3D Rankine panel method which adopts B-spline basis function in the time domain. The numerical code is based on the weakly nonlinear scheme which considers nonlinear Froude-Krylov and restoring forces. The main focus of this study is given to investigate the nonlinear characteristics of wave-induced loads, and to validate this present scheme for industrial use in the range of low Froude number. The comparisons show that the nonlinear motions and hull-girder loads, computed by the present numerical code, have good overall agreements with experimental results. It is found that, for the better accuracy of computational results, particularly at extreme waves in oblique seas, the careful treatment of soft-spring (or compatible) system is recommended to the control of non-restoring motions such as surge, sway, and yaw.