Effect of agricultural land use on the chemistry of groundwater from basaltic aquifers, Jeju Island, South Korea

Geochemical processes were identified as controlling factors of groundwater chemistry, including chemical weathering, salinization from seawater and dry sea-salt deposition, nitrate contamination, and rainfall recharge. These geochemical processes were identified using principal component analysis of major element chemistry of groundwater from basaltic aquifers in Jeju Island, South Korea, a volcanic island with intense agricultural activities. The contribution of the geochemical processes to groundwater chemistry was quantified by a simple mass-balance approach. The geochemical effects due to seawater were considered based on Cl contributions, whereas the effects due to natural chemical weathering were based on alkalinity. Nitrogenous fertilizers, and especially the associated nitrification processes, appear to significantly affect groundwater chemistry. A strong correlation was observed between Na, Mg, Ca, SO₄ and Cl, and nitrate concentrations in groundwater. Correspondingly, the total major cations, Cl, and SO₄ in groundwater were assessed to estimate relative effect of N-fertilizer use on groundwater chemistry. Cl originates more from nitrate sources than from seawater, whereas SO₄ originates mostly from rainwater. N-fertilizer use has shown the greatest effect on groundwater chemistry, particularly when nitrate concentrations exceed 6–7 mg/L NO₃–N. Nitrate contamination significantly affects groundwater quality and 18% of groundwater samples have contamination-dominated chemistry.     

Inverse groundwater modelling in the Willunga Basin,South Australia

A new inverse technique for modelling groundwater flow, based on a functional minimization technique, has been used to calibrate a groundwater flow model of a subregion of the Port Willunga aquifer within the Willunga Basin in South Australia. The Willunga Basin is the location of extensive viticulture, irrigated primarily by groundwater, the levels and quality of which have declined significantly over the last 40 years. The new method is able to generate estimates of transmissivity, storativity and groundwater recharge over the whole subregion as a time-varying continuous surface; previous methods estimate local discrete parameter values at specific times. The new method has also been shown to produce accurate head values for the subregion and very good estimates of groundwater recharge. Its ultimate goal will be to provide a new and invaluable tool for significantly improved groundwater resource management. Supported in part by US National Science Foundation grants, DMS-0107492 and DMS-0079478.     

The stiffness degradation and damping ratio evolution of Taipei Silty Clay under cyclic straining

The non-linear behavior of Taipei Silty Clay under cyclic strain loading was investigated through a series of undrained cyclic strain-controlled tests. The Ramberg–Osgood equation was used with our proposed stiffness degradation model to calculate degraded secant moduli. The proposed degradation model is simple in that it has only one more component than Idriss's model, the modulus ratio for the first cycle, which reflects the effects of the previous cyclic strain history and the current level of the cyclic strain amplitude, and can be used to describe softening and hardening behavior under irregular cyclic straining. It was found that the Ramberg–Osgood equation successfully predicts the damping ratio for small to medium strains. However, it overestimates the damping ratio for larger strains, so we suggest it can be corrected with a damping ratio index. In addition, the proposed equation for describing the evolution of the damping ratio provides the means to assess the variation for Taipei Silty Clay in the measured damping ratio with both the number of cycles and the strain amplitude.     

Determination of vertical displacements over the coastal area of Korea due to the ocean tide loading using GPS observations

This paper describes the GPS applicability for detecting the vertical displacements of ground stations caused by ocean tide loading effects. An experiment was carried out using 12 permanent GPS stations located in the coastal area of Korea using data in the period 1 July until 26 August 2003. The relative height differences were calculated from hourly DGPS data processing based on the carrier-phase observation. The power spectra of the M₂ and N₂ constituents of ocean tide loading were derived using the CLEAN algorithm. The differential vertical displacements generated by the ocean tide loading effect are typically 3–25 mm in coastal area of the Korea. We compared the results from GPS with those of the ocean tide models, NAO.99Jb regional model and GOT00.2, FES99 global models. The M₂ (N₂) amplitude differences of vertical displacements between GPS and GOT00.2 is 1.22 ± 3.61 mm (1.01 ± 1.48 mm), and that of the M₂ (N₂) amplitude difference between GPS and FES99 is 0.04 ± 4.64 mm (0.64 ± 1.75 mm), whereas the M₂ (N₂) amplitude difference between GPS and NAO.99Jb are 0.05 ± 1.03 mm (0.86 ± 1.18 mm). The highest vertical displacements at the PALM station are found for 24.5 ± 0.7 mm from GPS observation, and 22.9 mm from the regional model NAO.99Jb and 13.17 and 10.00 mm from the global models GOT00.2 and FES99, respectively. These values show that the vertical displacements derived from GPS are in good agreement with those of the regional model NAO.99Jb around Korea, more than with the global models. This result indicated that GPS is an effective tool to measure the vertical displacement caused by the ocean tide loading effect in the coastal area, and we need to use the NAO.99Jb ocean tide model rather than the global ocean tide models in and around the Korean peninsula for position determination with permanent GPS installations. This work demonstrates that vertical displacement caused by the M₂ and N₂ constituents of ocean tide loading can be measured by carrier-phase DGPS.     

Surface current variability in the Keum River Estuary (South Korea) during summer 2002 as observed by high-frequency radar and coastal monitoring buoy

High-frequency (HF) radar observations of surface currents were conducted for 3 months during summer 2002 in the Keum River estuary. A comparison between HF radar-derived currents and directly measured ones form a buoy showed that the regression slope is close to 1 and the correlation coefficient greater than 0.86, with an RMS difference less than 13 cm/s which is less than 17% of the tidal current. This fairly good agreement allows us to use HF radar observation in investigating the surface flow and circulation in this tidal-current-dominant coastal-plume area. To examine the spatial variation in tidal current characteristics, as well as currents associated with non-tidal forcing, the HF radar-derived currents were separated into tidal and sub-tidal frequency currents. The overall pattern of M₂-current ellipse distribution in the study area showed a counterclockwise rotation, with the offshore maximum current direction to the northeast. Eccentricity, the direction of maximum current, and the phase of net motion of the ellipse changed near the estuary mouth and near the gap of the Saemangeum reclamation tide dyke due to the complex coastal geometry and the out-flowing jet during the ebb period.     

Ekman motion in shallow open sea in the presence of time-harmonic variation of water depth

The time variation of the wind-induced flow in a homogeneous unbounded sea region has been analytically investigated. The time-dependent Ekman solution in a homogeneous, shallow open sea has been further extended by taking into account the time variation in water depth which might be caused by either tidal motion or change in the mixed layer thickness. The solution approach taken in this study is based on the Galerkin-eigenfunction method originally developed by Heaps [1972. On the numerical solution of three-dimensional hydrodynamical equations for tides and storm surges. Memoires de la Societe Royale des Sciences de Liege Serie 6(2), 143–180]. A series of calculations have been made with emphasis on the influence of the time variation in water depth upon the build-up of Ekman spirals in the presence of oscillatory variations in water depth. It has been found that two oscillations contribute to the wind drift current, the inertial oscillation and the depth-variation-induced oscillation; the inertial oscillation decays with time, but the depth-variation-induced oscillation remains undamped despite the presence of bottom friction. The presence of time-harmonic variation produces peculiar forms of hodograph with a curled or circular pattern according to the angular frequency of the water depth variation.     

Developing a bivariate spatial association measure: An integration of Pearson's r and Moran's I

 This research is concerned with developing a bivariate spatial association measure or spatial correlation coefficient, which is intended to capture spatial association among observations in terms of their point-to-point relationships across two spatial patterns. The need for parameterization of the bivariate spatial dependence is precipitated by the realization that aspatial bivariate association measures, such as Pearson's correlation coefficient, do not recognize spatial distributional aspects of data sets. This study devises an L statistic by integrating Pearson's r as an aspatial bivariate association measure and Moran's I as a univariate spatial association measure. The concept of a spatial smoothing scalar (SSS) plays a pivotal role in this task. Received: 07 November 2000 / Accepted: 02 August 2001