Factors controlling groundwater chemistry in an agricultural area with complex topographic and land use patterns in mid‐western South Korea

Shallow and bedrock groundwater from granitic aquifers were investigated for the hydrogeochemistry of major and minor constituents in an agricultural area. Nitrate concentrations were observed up to 49 mg/l as NO₃‐N, with 22% of samples exceeding the drinking water standard, which could pose a significant threat because most residents rely on groundwater as their drinking water source. Principal component analysis revealed three principal components (PCs): (1) nitrate contamination, contributed by major cations, Cl^?, SO and NO , (2) reduction processes positively involving Fe, Mn and B, and negatively involving dissolved oxygen and NO and (3) natural mineralization, involving HCO and F^?. Cluster analysis, performed on the PC scores, resulted in seven sample groups, which were successfully identified by total depth, elevation and land use. The nitrate‐contaminated groups had mixed land uses, with locally concentrated residential areas. Uncontaminated groundwater groups were found in the natural environment, including high‐altitude spring water and bedrock groundwater with a higher degree of natural mineralization. Shallow groundwater groups in paddy fields in lowlands were affected by reducing environments, of which one group was characterized by high Fe, Mn and B, and negligible nitrate. Groundwater with intermediate nitrate and lower Cl^? and SO was found primarily in hilly terrains with orchards and vegetable gardens, indicating lower contaminant loadings than lowland areas. Higher concentrations of F^? and nitrate were observed in the nitrate‐contaminated water, which seemed unlikely to be explained by groundwater mixing. The strong acidity generated from nitrification may infiltrate deeper into the aquifer, induce accelerated weathering of bedrock and result in the coexistence of F^? and nitrate, which may be an evidence of intense nitrate loading, leading to soil acidification. Multivariate statistical analysis successfully delineated hydrochemical characteristics of groundwater attained by natural and anthropogenic processes in an agriculturally stressed area with complex topographic land use patterns. Copyright © 2009 John Wiley & Sons, Ltd.     

Common Evolution of Mechanical and Transport Properties in Thermally Cracked Westerly Granite at Elevated Hydrostatic Pressure

Increasing the damage and crack porosity in crustal rocks can result in significant changes to various key physical properties, including mechanical strength, elastic and mechanical anisotropy, and the enhancement of transport properties. Using a Non-Interactive Crack Effective Medium (NIC) theory as a fundamental tool, we show that elastic wave dispersion can be inverted to evaluate crack density as a function of temperature and is compared with optically determined crack density. Further, we show how the existence of embedded microcrack fabrics in rocks also significantly influences the fracture toughness (K_IC) of rocks as measured via a suite of tensile failure experiments (chevron cracked notch Brazilian disk). Finally, we include fluid flow in our analysis via the Guéguen and Dienes crack porosity-permeability model. Using the crack density and aspect ratio recovered from the elastic-wave velocity inversion, we successfully compare permeability evolution with pressure with the laboratory measurements of permeability.     

Fracture Toughness and Fracture Roughness Interrelationship in Thermally treated Westerly Granite

This paper presents an experimental work aimed at assessing the correlation between fracture toughness (K _IC ) and fracture roughness for a series of Westerly granite specimens thermally treated up to 850°C. Mode I fracture toughness as a function of thermal treatment is determined using Cracked Chevron Notched Brazilian Disc specimens. The degree of roughness of the resultant fracture surfaces is analyzed with the aid of a high accuracy, high precision stereo-topometric measurement system. Roughness and toughness values display a negative correlation as a function of temperature. Fracture toughness decreases with increasing temperature due to the gradual opening of grain-grain boundaries in response to thermal stresses. Mode I fractures preferentially follow these weakened grain-grain boundaries, which in addition to the thermal expansion of individual grains, result in rougher failure profiles with increasing temperature. At low temperature, a distinct anisotropy in roughness was observed in all fracture surfaces with higher roughness values perpendicular to the direction of fracture propagation. However, higher treatment temperatures resulted in the homogenization of fracture roughness in all directions. These results confirm the important link among petrofabric analysis, fracture toughness, and fracture roughness in response to thermal treatment.     

火山地震活动的实验室模拟

火山体内产生地震活动的物理过程非常复杂,目前仍然没有被完全理解。我们在实验室内用来自埃特纳火山(意大利)的玄武岩做了变形和破碎实验并报道了实验结果。该实验用一组围绕样品的传感器监测,以完成全波形记录及微震事件定位和分析。充满液体的孔隙和损伤区快速破裂后减压激发了许多低频事件,类似于长周期火山地震活动。低频事件与孔隙流体的减压有关,位于破裂样本的损伤区;这些事件都具有一个弱剪切滑动分量(双力偶),与发生在活火山下的流体驱动事件相一致。     

An englacial image and water pathways of the Fourcade glacier on King George Island, Antarctic Peninsula, inferred from ground-penetrating radar

The distribution of small fractures and water content of the Fourcade glacier on King George Island, Antarctica, was investigated in November 2006 and December 2007 by two ground-based (470-and 490-m-long profiles) and one helicopter-borne (470-m-long profile) ground-penetrating radar (GPR) surveys using 50-, 100-, and 500-MHz antennas. Radar images in the pre-migrated GPR sections are characterized by a smooth ice surface and irregular bed topography, numerous diffraction hyperbolas in the ice and at the g...     

2D acoustic‐elastic coupled waveform inversion in the Laplace domain

Although waveform inversion has been intensively studied in an effort to properly delineate the Earth's structures since the early 1980s, most of the time‐ and frequency‐domain waveform inversion algorithms still have critical limitations in their applications to field data. This may be attributed to the highly non‐linear objective function and the unreliable low‐frequency components. To overcome the weaknesses of conventional waveform inversion algorithms, the acoustic Laplace‐domain waveform inversion has been proposed. The Laplace‐domain waveform inversion has been known to provide a long‐wavelength velocity model even for field data, which may be because it employs the zero‐frequency component of the damped wavefield and a well‐behaved logarithmic objective function. However, its applications have been confined to 2D acoustic media. We extend the Laplace‐domain waveform inversion algorithm to a 2D acoustic‐elastic coupled medium, which is encountered in marine exploration environments. In 2D acoustic‐elastic coupled media, the Laplace‐domain pressures behave differently from those of 2D acoustic media, although the overall features are similar to each other. The main differences are that the pressure wavefields for acoustic‐elastic coupled media show negative values even for simple geological structures unlike in acoustic media, when the Laplace damping constant is small and the water depth is shallow. The negative values may result from more complicated wave propagation in elastic media and at fluid‐solid interfaces. Our Laplace‐domain waveform inversion algorithm is also based on the finite‐element method and logarithmic wavefields. To compute gradient direction, we apply the back‐propagation technique. Under the assumption that density is fixed, P‐ and S‐wave velocity models are inverted from the pressure data. We applied our inversion algorithm to the SEG/EAGE salt model and the numerical results showed that the Laplace‐domain waveform inversion successfully recovers the long‐wavelength structures of the P‐ and S‐wave velocity models from the noise‐free data. The models inverted by the Laplace‐domain waveform inversion were able to be successfully used as initial models in the subsequent frequency‐domain waveform inversion, which is performed to describe the short‐wavelength structures of the true models.     

Effects of spilled oil on microbial communities in a tidal flat

Effects of spilled oil on microbial communities in tidal flats were examined by use of a simulator for a tidal flat ecosystem. The simulator is composed of a wave generator, a tide control device, and a tidal flat. Sediment for the tidal flat was obtained at a natural tidal flat in Hiroshima Bay, Japan. After stabilizing the benthic organisms, fuel oil C was added to the surface of the flat at 1 lm(-2). Although the total number of micro-organisms remained at 1.5-3.5 x 10(9) cells g(-1) dry sediment irrespective of the addition of oil, bacterial communities which were analyzed based on the 16S rDNA showed clear changes after the addition of fuel oil C and after a subsequent recovery period. Bacterial colonies were randomly isolated from the oil-supplemented sediment during the experiments, and the isolates were examined for susceptibility to hydrocarbons in order to screen the oil-susceptible bacteria. The proportion of oil-susceptible bacteria in the isolates decreased with the addition of the oil. Oil-susceptible bacteria showed an inability to assimilate petroleum compounds as well as an inhibition of growth. The possibility of using oil-susceptible bacteria as an indicator of bioremediation in tidal flats was discussed.     

The Stratification of Jupiter's Troposphere at the Galileo Probe Entry Site

Galileo Probe Atmospheric Structure Investigation (ASI) pressure and temperature sensor data acquired during the parachute descent phase have been used to derive the static stability structure of Jupiter's troposphere at pressure levels of 0.5-22 bars using three techniques. The first approach utilizes both the p-sensor and T-sensor data, but since the p-sensor's zero offset was significantly affected by the thermal anomaly in the probe, two other approaches using only T-sensor data have also been developed. By making the physically reasonable assumptions of equilibrium descent for the probe and hydrostatic balance of the atmosphere, an algorithm for deriving the background static stability from T-sensor measurements alone is developed. Regions with static stability 0.1-0.2 K km⁻¹ are found at 0.5-1.7 bars, 3-8.5 bars, and 14-20 bars. Between these layers, regions of weaker static stability are present. Mean molecular weight gradients due to the vertical variation of water vapor abundance near the 11-bar pressure level appear to stabilize the atmosphere at this level. Oscillatory structures with vertical wavelength ∼15-30 km and amplitude ∼0.1-0.2 K are observed in the T-sensor data. For pressure <2 bars, these eddies are well above the noise level of the measurements and are consistent with the predictions of linear gravity wave theory for a wave with horizontal phase speed c_x=160 m s⁻¹ with respect to System III propagating through the static stability derived from the T-sensor data alone. They provide quantitative confirmation of the static stability derived from T-sensor data in the troposphere where p<2 bars. The observed static stability structure shows an inverse correlation with the regions of wind shear observed by the Doppler Wind Experiment: regions of highest shear in the horizontal wind appear to be associated with regions of lowest static stability. The particulate population detected by other experiments on the probe shows some correlation with the uppermost layer of static stability, suggesting enhanced solar energy deposition at these levels may play a role in producing the positive static stability.     

Standard and Non-Standard Solar Models

We summarize the physical input and assumptions commonly adopted in modern standard solar models that also produce good agreement with solar oscillation frequencies. We discuss two motivations for considering non-standard models: the solar neutrino problem and surface lithium abundance problem. We begin to explore the potential for mixed core models to solve the neutrino problem, and compare the structure, neutrino flux, and oscillation frequency predictions for several models in which the inner 25% of the radius is homogenized, taking into account the effects of non-local equilibrium abundances of ³He. The results for the neutrino flux and helioseismic predictions are far from satisfactory, but such models have the potential to reduce the predicted ⁷Be/⁸B neutrino flux ratio, and further studies are warranted. Finally, we discuss how much the neutrino problem can be alleviated in the framework of the standard solar model by using reaction rates, abundances and neutrino capture cross-sections at the limits of their uncertainties, while still satisfying the constraints of helioseismology.     

Photometric observations of NY Cephei

NewUBV photoelectric observations of NY Cep were made from 1982 to 1989 with GOTO 61-cm telescope at the Yonsei University Observatory. Seven times of minimum light were obtained. An improved orbital period of 15.27566 day was determined. The light curves were analyzed using the Wilson-Devinney method. It is confirmed that NY Cep is a detached system and the absolute dimensions are also obtained. It is found that the secondary eclipse does not occur due to large eccentricity and inclination.