A New Automatic Subsurface Gas Monitoring System for Seismogeochemical Studies, Installed in Haruno Borehole, Shizuoka Prefecture, Central Japan

The results of subsurface gas monitoring by application of gas chromatography (GC) to the gas composition of bubbles associated with groundwater for seismogeochemical studies are reported. An automated gas monitoring system was used to determine gas compositions in a 500-m borehole at the Haruno Crustal Movement Observation Site (HOS), central Japan during period 1, from December 1999 to December 2000. The average ± two standard deviation (2SD) compositions of gases in this period were He = 82 ± 29 ppmV, H₂ = 170 ± 62 ppmV, Ar = 0.05 ± 0.07%, N₂ = 50 ± 8%, and CH₄ = 45 ± 6%. A new automated gas monitoring system equipped with a micro-GC was installed in the borehole at the HOS, and gas bubbles from the borehole were monitored during period 2, from December 2006 to March 2007. The average ± two standard deviation (2SD) compositions of gases in this period were He = 8 ± 7 ppmV, H₂ = 13 ± 15 ppmV, Ar = 0.6 ± 0.3%, N₂ = 66 ± 7%, and CH₄ = 14 ± 14%. The gas concentration ratios (He/Ar, H₂/Ar, N₂/Ar, and CH₄/Ar) fluctuated significantly over time and repeatedly showed abrupt spike-like increases during period 2. The gas compositions obtained in period 1 and 2 were markedly different. Over the period from 2006 to 2007, the gas bubbles were depleted in He, H₂, and CH₄ of deep origin, but enriched in Ar and N₂ of atmospheric origin. This difference can be interpreted as being due to an irreversible change of the aquifer/gas system. The present deep component in the HOS gas is estimated to have composition He = 63 ppmV, H₂ = 37 ppmV, Ar = 0.17%, N₂ = 63%, and CH₄ = 37%. The new monitoring system is able to analyze the gas composition using a smaller volume of sample gas and with greater precision than the previous system. During the 3-month monitoring period 2, the separation capacity of the capillary column of the micro-GC was sufficiently maintained to determine gas-chromatographic peak areas for the five gaseous species examined. This study confirms that the new monitoring system with micro-GC is promising for continuous subsurface gas monitoring for earthquake prediction studies.     

Appearance pattern of tooth germs with developmental process inMylopharyngodon piceus

The pharyngeal dental formula ofMylopharyngodon piceus is 4–5 as a rule, and the dentition is asymmetrical. It is difficult to identify each tooth in the larval dentition. In this paper the appearance pattern of tooth germ with development process in this fish is described in detail. The formation pattern of the left dentition is contrasted with that of the right one. In the developmental process, the left pharyngeal dentition lacks teeth at position An3. Thus the left dentition is D-type as designated by Nakajima (1984), while the right one is A-type. This project was funded by the International Cooperation of Japan-China and the National Natural Science Foundation of China.     

Specific affinity for phosphate uptake and specific alkaline phosphatase activity as diagnostic tools for detecting phosphorus-limited phytoplankton and bacteria

We analyzed responses of soluble reactive phosphorus (SRP), bioavailable phosphate (PO₄), particulate phosphorus, turnover time of orthophosphate (T_t), and alkaline phosphatase activity (APA) to varying degrees of nutrient stress. The nutrient stress was evoked by different treatments in concentration and combination of inorganic nitrogen (N) and phosphorus (P), and labile organic carbon (glucose) to mesocosms in experiments carried out in eutrophic southern (Odense Fjord, Denmark) and northern (Tvärminne Archipelago, Finland) coastal zones of the Baltic Sea. Despite seasonal and geographical differences, similar responses were observed in both experiments. Low SRP (<100 nmol l^?1), shortT _t (<10 h), and increased levels of APA were observed in both N+P balanced and P deficient treatments, while the opposite trend was observed in P replete treatments. The shortestT _t and the highest APA were found when glucose was combined with N treatment. Bioavailable PO₄ was estimated usingT _t and P uptake rates as derived from stoichiometric conversion of carbon based primary and bacterial production. With shorterT _t, the PO₄ pool declined to <1 nmol-P l^?1, whereas the SRP background pool (difference between SRP and PO₄) remained relatively constant (c. 50 nmol l^?1). APA was inversely related to PO₄ but not to SRP. Responses of specific APA and specific affinity for PO₄ uptake, which are APA and PO₄ uptake rates (inverse ofT _t), respectively, normalized to the summed P biomass of phytoplankton and bacteria, responded consistently to the pattern and magnitude of nutrient limitation evoked in our experiments. Our results, together with a literature survey, suggest that both parameters can be useful in examining PO₄ availability for the natural phytoplankton and bacteria community in P starved aquatic systems.     

A lower trophic ecosystem model including iron effects in the Okhotsk Sea

We applied a three-dimensional ecosystem-physical coupled model including iron the effect to the Okhotsk Sea. In order to clarify the sources of iron, four dissolved iron compartments, based on the sources of supply, were added to Kawamiya et al.'s [1995, An ecological-physical coupled model applied to Station Papa. Journal of Oceanography, 51, 635-664] model (KKYS) to create our ecosystem model (KKYS-Fe). We hypothesized that four processes supply iron to sea water: atmospheric loadings from Northeastern Asia, input from the Amur River, dissolution from sediments and regeneration by zooplankton and bacteria. We simulated one year, from 1 January 2001 to 31 December 2001, using both KKYS-Fe and KKYS. KKYS could not reproduce the surface nitrate distribution after the spring bloom, whereas KKYS-Fe agreed well with observations in the northwestern Pacific because it includes iron limitation of phytoplankton growth. During the spring bloom, the main source of iron at the sea surface is from the atmosphere. The contribution of riverine iron to the total iron utilized for primary production is small in the Okhotsk Sea. Atmospheric deposition, the iron flux from sediment and regeneration of iron in the water column play important roles in maintaining high primary production in the Okhotsk Sea.     

Estimation of the apertures of water-saturated fractures by nuclear magnetic resonance well logging

Fracture aperture is an important transport property in subsurface hydrology because it influences well productivity and the volume of the water resource. Nuclear magnetic resonance (NMR) well logging measures the hydrogen‐bearing fluid molecules in porous or fractured strata, and the NMR signal intensity increases with the amount of fluid in the sensed region of the NMR sonde. Fluid confined in a large fracture of >>0.2 mm in aperture has T2 (i.e. spin‐spin relaxation time) values as long as those of the bulk fluid. The bulk‐fluid porosity (i.e. porosity calculated using this long T2 component in a T2 histogram data) increases linearly with aperture. Therefore, NMR logging enables quantitative estimation of fracture apertures of >>0.2 mm using the bulk‐fluid porosity data if the calibration of the NMR sonde is performed adequately. We applied NMR logging to a borehole in a Holocene andesite lava at Sumikawa, Japan, to estimate the aperture of open fractures within the lava. A test well of 100 m depth and 20 cm diameter, filled with bentonite drilling mud, was scanned with an NMR sonde to obtain a profile of the porosity and the T2 histogram of the andesite. The bulk‐fluid porosity was calculated from the T2 histogram data, as the porosity at which the T2 value is larger than or equal to a threshold T2 of bulk bentonite mud. The bulk‐fluid porosity of a specific inclined fracture responsible for the total loss of circulation at 61.2 m depth during drilling was calculated assuming a threshold or T2 cut‐off of 33 ms, and again for a cut‐off of 100 ms. Calibration of the NMR sensor in a laboratory and measurement of the fracture dip angle by electrical microimaging logging enabled us to estimate the fracture aperture as 1.7 cm, assuming a T2 cut‐off of 33 ms, or 1.6 cm for a T2 cut‐off of 100 ms. The method of aperture determination described in this study is independent of fluid species and lithology, and is applicable to various hydrogen‐bearing borehole fluids (clean water, mud and oil) and geological settings.     

Application of Continuous 222Rn Monitor with Dual Loop System in a Small Lake

To estimate the spatial distribution of groundwater discharge from the bottom of a small lake of Kumamoto in Japan, we applied continuous radon measurements with a dual loop system and verified the results obtained using the radon method by visual diving surveys. Time‐shifting correction in the dual‐loop system is reasonable to obtain the true radon activity in water. Distributions of radon activity and water temperature in the study area reveal the effects on groundwater discharge and mixing situation of lake water. The estimated discharge zone ascertained using the radon method agrees with the groundwater discharge distribution observed through diving surveys. Although the data resolution of the radon method is much greater than the width of observed discharge zones, the general distribution of groundwater discharge is recognizable. The dual loop system of radon measurement is useful for smaller areas.     

Comparison of stable isotopes,ratios of Cl/Cl and I/I in brine and deep groundwater from the Pacific coastal region and the eastern margin of the Japan Sea

Fifty-three samples, including brines associated with oil and natural gas reservoirs and groundwater samples from deep boreholes, were collected from the Pacific and Japan Sea coastal regions in Japan. The ¹²⁹I/¹²⁷I and ³⁶Cl/Cl ratios, and stable isotopes (δD and δ¹⁸O) are compared to investigate differences related to the geotectonic settings of the two regions. The δD and δ¹⁸O data indicate that brine and groundwater from the Pacific coastal region reflect mixing of meteoric water with connate seawater in the pores of sedimentary rocks. On the other hand, brine and groundwater from the Japan Sea coastal region have been hydrothermally altered. In particular, brines associated with petroleum accumulations at Niigata and Akita showed the same isotopic characteristics as fluids found in the Kuroko deposits of the Green Tuff region in northeastern Japan. There is little difference in the ³⁶Cl/Cl ratios in brine and groundwater from the Pacific and Japan Sea coasts. Most brine and some deep groundwater, except those from the Pleistocene Kazusa Group, have already reached the average secular equilibrium ratio of 9.9 ± 2.7 × 10⁻¹⁵ for their mudstone and sandstone reservoirs. There was no correlation between the ³⁶Cl/Cl ratios and differences in geotectonic setting between the Pacific and the Japan Sea coast. The molar I/Br ratio suggests that the I in all of water samples was of biogenic origin. The average ¹²⁹I/¹²⁷I ratio was 290 ± 130 × 10⁻¹⁵ to 294 ± 105 × 10⁻¹⁵ in both regions, showing no relationship to the different geotectonic settings. The uncontaminated brine and groundwater samples are likely to have retained the original ¹²⁹I/¹²⁷I ratios of marine I released from the old organic matter stored in sedimentary rock.     

Assessment of rainfall-induced shallow landslides in Phetchabun and Krabi provinces,Thailand

Shallow landslides are a common type of rainfall-induced landslide, and various methods are currently used to predict their occurrence on a regional scale. Physically based models, such as the shallow landslide instability prediction (SLIP) model, have many advantages because these models can assess the hazards of shallow landslides dynamically, based on physical stability equations that consider rainfall as a triggering factor. The main objective of this research is to test the SLIP model’s potential to predict shallow landslide hazards in Thailand. To achieve this goal, the SLIP model was applied to two massive landslide events in Thailand. The results predicted by the SLIP model for the two study areas are outlined, and the model prediction capabilities were evaluated using the receiver operating characteristic plot. The Phetchabun results showed that the western part of the catchment had the lowest factor of safety (F _S) value, whereas the Krabi results showed that the slopes surrounding the peak of Khao Panom Mountain had the lowest F _S value, explaining the highest potentials for shallow landslides in each area. The SLIP model showed good performance: The global accuracies were 0.828 for the Phetchabun area and 0.824 for the Krabi area. The SLIP model predicted the daily time-varying percentage of unstable areas over the analyzed periods. The SLIP model simulated a negligible percentage of unstable areas over all considered periods, except for expected dates, suggesting that the prediction capability is reasonably accurate.     

A Simulation Study of the Plasma Wave Enhancements in the Earth’s Equatorial Plasmasphere

In the equatorial plasmasphere, plasma waves are frequently observed. To improve our understanding of the mechanism generating plasma waves from instabilities, a comparison of observations, linear growth-rate calculations, and simulation results is presented. To start the numerical experiments from realistic initial plasma conditions, we use the initial parameters inferred from observational data obtained around the plasma-wave generation region by the Akebono satellite. The linear growth rates of waves of different modes are calculated under resonance conditions, and compared with simulation results and observations. By employing numerical experiments by a particle code, we first show that upper hybrid-, Z-, and whistler-mode waves are excited through instabilities driven by a ring-type velocity distribution. The simulation results suggest a possibility that energetic electrons with energies of some tens of keV confined around the geomagnetic equator are responsible for the observed enhancements of Z- and whistler-mode waves. While the comparison between linear growth-rate calculations and observations shows the different tendency of wave amplitude of Z-mode and whistler-mode waves, the wave amplitude of these wave modes in the simulation results is consistent with the observation.     

Distribution and dispersal of the early larval stages of the Japanese spiny lobster Panulirus japonicus in the East China Sea

The annual distribution and dispersal of early‐stage phyllosoma of the Japanese spiny lobster Panulirus japonicus were examined in the East China Sea and the Sea of Japan off western Japan. Early‐stage larvae were sampled mainly in summer and relatively near the coast off western Kyusyu Island. Few larvae were found in the other three seasons in the East China Sea and the Sea of Japan. This finding suggests that P. japonicus larvae are retained in coastal areas for a few months after hatching until stage V and are transported offshore until autumn, before they grow to stage VI. Examination of archival drifter data off western Kyusyu Island during the spawning season of P. japonicus supported the relatively rapid dispersal of the larvae from the East China Sea near Japan to the Pacific Ocean or the Sea of Japan. Larval transport from the East China Sea to the Pacific Ocean, which is considered to be a main distribution area of middle‐ and late‐stage larvae, would occur in the south at approximately 32–33°N in the East China Sea near western Kyusyu Island.