Groundwater Impacts from the M5.8 Earthquake in Korea as Determined by Integrated Monitoring Systems

This paper describes the impacts of the M5.8(5.1) Gyeongju earthquakes on groundwater levels using data obtained from a unique coastal monitoring well. The monitoring strategy integrates conventional water level monitoring with periodic, continuous measurements of temperature and electrical conductivity (EC) within the water column of the well. Another important component of the monitoring system is a new instrument, the InterfacEGG, which is capable of dynamically tracking the freshwater-saltwater interface. Although the system was set up to monitor seawater intrusion related to over-pumping, as well as rainfall and tidal effects, it recorded impacts associated with a large earthquake and aftershocks approximately 241 km away. Seismic energies associated with the M5.8(5.1) Gyeongju earthquakes induced groundwater flows to the monitoring well through fractures and joints in the crystalline basement rocks. Temperature and EC logging data showed that the EC vertical profile declined from an average of approximately 5300 to 4800 μS/cm following the earthquakes. The temperature profile showed a trend toward lower temperatures as the depth increased, a feature not commonly observed in previous studies. Data from the InterfacEGG suggested that the rise in EC was not due to the saltwater intrusion, but from the tendency for brackish water entering the borehole to induce convective mixing at deeper depths as the seismic waves travel through the well-aquifer system. The increase in groundwater levels was caused by pulse of colder, less brackish water flowing into the well because of the earthquake. This behavior reflects an enhancement in rock permeability by removing precipitates and colloidal particles from clogged fractures, which improve the hydraulic connection with a nearby unit with a higher hydraulic head. This study suggests there is value added with a more aggressive monitoring strategy.     

A generic approach to modelling flexible confined boundary conditions in SPH and its application

In this paper, a new approach to applying confining stress to flexible boundaries in the smoothed particle hydrodynamics (SPH) method is developed to facilitate its applications in geomechanics. Unlike the conventional SPH methods that impose confining boundary conditions by creating extra boundary particles, the proposed approach makes use of kernel truncation properties of SPH approximations that occur naturally at free-surface boundaries. Therefore, it does not require extra boundary particles and, as a consequence, can be utilised to apply confining stresses onto any boundary with arbitrary geometry without the need for tracking the curvature change during the computation. This enables more complicated problems that involve moving confining boundaries, such as confining triaxial tests, to be simulated in SPH without difficulties. To further enhance SPH applications in elasto-plastic computations of geomaterials, a robust numerical procedure to implement Mohr-Coulomb plasticity model in SPH is presented for the first time to avoid difficulties associated with corner singularities in Mohr-Coulomb model. The proposed approach was first validated against two-dimensional finite element (FE) solutions for confining biaxial compression tests to demonstrate its predictive capability at small deformation range when FE solutions are still valid. It is then further extended to three-dimensional conditions and utilised to simulate triaxial compression experiments. Simulation results predicted by SPH show good agreement with experiments, FE solutions, and other numerical results available in the literature. This suggests that the proposed approach of imposing confining stress boundaries is promising and can handle complex problems that involve moving confining boundary conditions.     

Implementation of turbulent mixing over a stratocumulus-topped boundary layer and its impact in a GCM

The effect of a vertical diffusion scheme over a stratocumulus topped boundary layer (STBL) was investigated using the YONU AGCM (Yonsei University Atmospheric General Circulation Model). To consider the impact of clouds on the turbulence production, the turbulence mixing term, driven by radiative cooling at the cloud top, is implemented as an extended non-local diffiusion scheme. In the model with this new scheme, the STBL parameterization significantly influences the lower atmosphere over the tropical and...     

Interdecadal changes in the storm track activity over the North Pacific and North Atlantic

Analysis of NCEP-NCAR I reanalysis data of 1948–2009 and ECMWF ERA-40 reanalysis data of 1958–2001 reveals several significant interdecadal changes in the storm track activity and mean flow-transient eddy interaction in the extratropics of Northern Hemisphere. First, the most remarkable transition in the North Pacific storm track (PST) and the North Atlantic storm track (AST) activities during the boreal cold season (from November to March) occurred around early-to-mid 1970s with the characteristics of global intensification that has been noticed in previous studies. Second, the PST activity in midwinter underwent decadal change from a weak regime in the early 1980s to a strong regime in the late 1980s. Third, during recent decade, the PST intensity has been enhanced in early spring whereas the AST intensity has been weakened in midwinter. Finally, interdecadal change has been also noted in the relationship between the PST and AST activities and between the storm track activity and climate indices. The variability of storm track activity is well correlated with the Pacific Decadal Oscillation and North Atlantic Oscillation prior to the early 1980s, but this relationship has disappeared afterward and a significant linkage between the PST and AST activity has also been decoupled. For a better understanding of the mid-1970s’ shift in storm track activity and mean flow-transient eddy interaction, further investigation is made by analyzing local barotropic and baroclinic energetics. The intensification of global storm track activity after the mid-1970s is mainly associated with the enhancement of mean meridional temperature gradient resulting in favorable condition for baroclinic eddy growth. Consistent with the change in storm track activity, the baroclinic energy conversion is significantly increased in the North Pacific and North Atlantic. The intensification of the PST and AST activity, in turn, helps to reinforce the changes in the middle-to-upper tropospheric circulation but acts to interfere with the changes in the low-tropospheric temperature field.     

Changes in climate classification and extreme climate indices from a high-resolution future projection in Korea

We investigate the future changes in the climate zone and six extreme temperature indices in Korea, using the 20-km high-resolution atmospheric general circulation model (MRI-AGCM3.1S). The Trewartha and K?ppen climate classification schemes are applied, and four summer-based extreme temperature indices (i.e., summer days, tropical nights, growing degree days, and cooling degree days (CDD) and two winter-based indices (frost days and heating degree days (HDD) are analyzed. To represent significantly the change in threshold indices, the monthly mean bias is corrected in model. The model result reasonably captures the temporal and spatial distribution of the present-day extreme temperatures associated with topography. It was found that in the future climate, the area of the subtropical climate zone in Korea expands northward and increases by 21% under the Trewartha classification scheme and by 35% under the K?ppen classification scheme. The spatial change in extreme climate indices is significantly modulated by geographical characteristics in relation to land-ocean thermal inertia and topographical effects. The change is manifested more in coastal regions than in inland regions, except for that in summer days and HDD. Regions with higher indices in the present climate tend to reveal a larger increase in the future climate. The summer-based indices display an increasing trend, while the winter-based indices show a decreasing trend. The most significant increase is in tropical nights (+452%), whereas the most significant decrease is in HDD (?25%). As an important indicator of energy-saving applications, the changes in HDD and CDD are compared in terms of the frequency and intensity. The future changes in CDD reveal a higher frequency but a lower temperature than those in HDD. The more frequent changes in CDD may be due to a higher and less dispersed occurrence probability of extreme temperatures during the warm season. The greater increase in extreme temperature events during the summer season remains an important implication of projecting future changes in extreme climate events.     

Arsenic levels in human hair,Kandal Province,Cambodia: The influences of groundwater arsenic,consumption period,age and gender

This study focused on the analysis of As levels in human hair samples collected from six villages in the Kandal Province of Cambodia. Of interest were the influence of, and interactions among, certain factors affecting As intake into the human body: As concentrations in groundwater, period of groundwater consumption, age and gender. The results revealed As levels in human hair ranging from 0.06 to 30 μg g⁻¹ with median and arithmetic mean values of 0.61 and 3.20 μg g⁻¹ (n = 68), respectively. Furthermore, a linear relationship was found between As concentrations in human hair and in the local groundwater. Arsenic (III) is the dominant species in Kandal groundwater, constituting in most cases at least 60% of the total As. Arsenic concentration ranged from 5 to 1543 μg L⁻¹, with the median value 348 μg L⁻¹ and arithmetic mean 454 μg L⁻¹. In large rural, poor areas holding most of Kandal’s 1.1 million people, up to 2 in 1000 people are believed to be at risk of cancer through the As-enriched water they drink. A toxicity risk assessment provides a hazard quotient (HQ) equaling 5.12, also a clear indication of non-carcinogenic exposure risk. On the authors’ visit to Kampong Kong commune, Kandal Province, cases of arsenicosis were diagnosed in patients as a result of drinking As-enriched groundwater.     

Evaluation of multiple regression models using spatial variables to predict nitrate concentrations in volcanic aquifers

Multiple linear regression of spatial variables including land use, soil type, and topography was applied to predict nitrate concentration and evaluate major factors affecting nitrate occurrence in springs and wells in the southern and northern areas of Jeju volcanic island, Korea. Three types of contributing area surrogates (CAS), namely circle, semicircle, and wedge, were employed to calculate the spatial variables. The regression results showed R² of 0.81–0.84 for springs and 0.74–0.77 for wells; R² values for wedge and semicircular CAS were more than 10% higher than those for circular CAS. The R² of spring models was significantly affected by both the shape and size of CAS, with optimal radii of 150–250 m and 300–400 m in the southern and northern areas, respectively, corresponding to thinner upper basaltic aquifers, and implying shorter flow paths in the southern area. The most influential variables in springs were orchards and soil types related to agriculture including silty loam and silty clay loam, indicating that nitrate levels are strongly affected by N fertilization in cultivated areas. In contrast, wells showed much less sensitivity to both shapes and sizes of CAS, with less contribution of land use and soil type to the regression, which could be attributed to a mix of multiple aquifer zones and widely different factors in the installation and operation of wells. Field parameters of electrical conductivity (EC) and pH increased the R² up to 10%, suggesting that these can be useful when regression with spatial variables yields a lower R². The optimal spatial scales for prediction of nitrate concentration and spatial variables that significantly contribute to nitrate contamination can provide relevant criteria for establishing groundwater management policies, considering the increasing anthropogenic land‐use trends on the island, where groundwater is highly sensitive to changes in spatial variables. Copyright © 2015 John Wiley & Sons, Ltd.     

A comparison of climatological subseasonal variations in the wintertime storm track activity between the North Pacific and Atlantic: local energetics and moisture effect

Distinct differences of the storm track?Cjet relationship over the North Pacific and North Atlantic are investigated in terms of barotropic and baroclinic energetics using NCEP-2 reanalysis data for the period of 1979?C2008. From fall to midwinter the Pacific storm track (PST) activity weakens following the southward shift of the Pacific jet, whereas the Atlantic storm track (AST) activity remains steady in position and intensifies regardless of the slight southward shift of the Atlantic jet. This study is devoted to seeking for the factors that can contribute to this conspicuous difference between the two storm tracks on climatological subseasonal variation by analyzing eddy properties and local energetics. Different eddy properties over the two oceans lead to different contribution of barotropic energy conversion to the initiation of storm tracks. In the North Atlantic, meridionally elongated eddies gain kinetic energy efficiently from stretching deformation of the mean flow in the jet entrance. On the other hand, the term associated with shearing deformation is important for the initiation of PST. Analysis of baroclinic energetics reveals that the intensification of the AST activity in midwinter is mainly attributed to coincidence between location of maximum poleward and upward eddy heat fluxes and that of the largest meridional temperature gradient over slight upstream of the AST. The relatively large amount of precipitable water and meridional eddy moisture flux along baroclinic energy conversion axis likely provides a more favorable environment for baroclinic eddy growth over the North Atlantic than over the North Pacific. In the meantime, the midwinter minimum of the PST activity is attributable to the southward shift of the Pacific jet stream that leads to discrepancy between core region of poleward and upward heat fluxes and that of meridional thermal gradient. Weakening of eddy-mean flow interaction due to eddy shape and reduction of moist effect are also responsible for the weakening of storm track activities in midwinter when the strongest baroclinicity exists over the North Pacific.     

Assessing sustainability of groundwater resources on Jeju Island,South Korea,under climate change,drought, and increased usage

Numerical groundwater models were used to assess groundwater sustainability on Jeju Island, South Korea, for various climate and groundwater withdrawal scenarios. Sustainability criteria included groundwater-level elevation, spring flows, and salinity. The latter was studied for the eastern sector of the island where saltwater intrusion is significant. Model results suggest that there is a need to revise the current estimate of sustainable yield of 1.77?×?10⁶ m³/day. At the maximum extraction of 84  % of the sustainable yield, a 10-year drought scenario would decrease spring flows by 28 %, dry up 27 % of springs, and decrease hydraulic head by an island-wide average of 7 m. Head values are particularly sensitive to changes in recharge in the western parts of the island, due to the relatively low hydraulic conductivity of fractured volcanic aquifers and increased groundwater extraction for irrigation. Increases in salinity are highest under drought conditions around the current 2-m head contour line, with an estimated increase of up to 9 g/L under 100 % sustainable-yield use. The study lists recommendations towards improving the island’s management of potable groundwater resources. However, results should be treated with caution given the available data limitations and the simplifying assumptions of the numerical modeling approaches.     

Study on the hydraulic parameter differences between pumping and recharge in a confined aquifer of a soft soil area

Artificial recharge is an effective remediation measure for controlling groundwater level and subsidence in many coastal cities in China. Hydraulic parameters estimated by pumping tests are often used in the design of both pumping and recharge systems. However, the hydraulic parameters in the recharge process have been found to differ from those in the pumping process and should be studied in greater detail. Eight single-well pumping and recharge tests were conducted within a confined aquifer in a soft soil area in the city of Tianjin, and the differences in wellbore storage influences and well losses between the recharge and pumping processes were examined. Furthermore, based on the Hantush and Jacob model, an algorithm combining the Levenberg–Marquardt algorithm (LMA) and genetic algorithm (GA) was employed for estimation of the hydraulic parameters. The results illustrated that the combined algorithm eliminating wellbore storage influences could provide hydraulic parameters from which the groundwater level variation could be accurately simulated. The hydraulic conductivity and specific storage values obtained in the pumping tests were higher than those obtained in the recharge tests. In addition to slight plugging of the recharge well, the specific storage differences could be explained by the compression and rebound deformation characteristics of sand in the confined aquifer. The specific storage estimated by pumping tests should be adjusted when applied in groundwater recharge calculation.