Overview of the 1990–1995 eruption at Unzen Volcano

Following 198 years of dormancy, a small phreatic eruption started at the summit of Unzen Volcano (Mt. Fugen) in November 1990. A swarm of volcano-tectonic (VT) earthquakes had begun below the western flank of the volcano a year before this eruption, and isolated tremor occurred below the summit shortly before it. The focus of VT events had migrated eastward to the summit and became shallower. Following a period of phreatic activity, phreatomagmatic eruptions began in February 1991, became larger with time, and developed into a dacite dome eruption in May 1991 that lasted approximately 4 years. The emergence of the dome followed inflation, demagnetization and a swarm of high-frequency (HF) earthquakes in the crater area. After the dome appeared, activity of the VT earthquakes and the summit HF events was replaced largely by low-frequency (LF) earthquakes. Magma was discharged nearly continuously through the period of dome growth, and the rate decreased roughly with time. The lava dome grew in an unstable form on the shoulder of Mt. Fugen, with repeating partial collapses. The growth was exogenous when the lava effusion rate was high, and endogenous when low. A total of 13 lobes grew as a result of exogenous growth. Vigorous swarms of LF earthquakes occurred just prior to each lobe extrusion. Endogenous growth was accompanied by strong deformation of the crater floor and HF and LF earthquakes. By repeated exogenous and endogenous growth, a large dome was formed over the crater. Pyroclastic flows frequently descended to the northeast, east, and southeast, and their deposits extensively covered the eastern slope and flank of Mt. Fugen. Major pyroclastic flows took place when the lava effusion rate was high. Small vulcanian explosions were limited in the initial stage of dome growth. One of them occurred following collapse of the dome. The total volume of magma erupted was 2.1×10⁸ m³ (dense-rock-equivalent); about a half of this volume remained as a lava dome at the summit (1.2 km long, 0.8 km wide and 230–540 m high). The eruption finished with extrusion of a spine at the endogenous dome top. Several monitoring results convinced us that the eruption had come to an end: the minimal levels of both seismicity and rockfalls, no discharge of magma, the minimal SO₂ flux, and cessation of subsidence of the western flank of the volcano. The dome started slow deformation and cooling after the halt of magma effusion in February 1995.     

Theoretical Evaluation of the Mechanical Behavior of Unsaturated Soils

An evaluation method for the mechanical behavior of unsaturated soils is studied in this paper. Although the mechanical behavior of unsaturated soils is complicated, a simple modeling is preferable in practice. This is because the soil properties are not homogeneous and ground data is limited when structures are being designed. In addition, in order to evaluate the reliability of the design, the physical meanings of the parameters applied in the prediction model should be clear. Firstly, the authors study the relationship between compaction curves and compression indexes in the unsaturated state that is used in the proposed constitutive model. Based on the constitutive model, the stress paths for constant volume shear tests are formulated under a constant void ratio condition and the stress paths for undrained shear tests are calculated under a constant water content condition. In the case of unsaturated specimens, the volume of these specimens changes with the shear deformation and the stress paths depend on the initial degree of saturation. The results of the calculation qualitatively describe the test results by considering the changes in effective confining pressure in the undrained condition and the water retention curves.     

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.     

Analysis of the geographical accessibility of neurosurgical emergency hospitals in Sapporo city using GIS and AHP

To improve the geographical accessibility of neurosurgical emergency hospitals for elderly people, we developed several alternative site plans for a new neurosurgical emergency hospital in Sapporo, Japan. Hospitals, population data, routes, and the numerical information for the Analytic Hierarchy Process computations were input into a Geographical Information System. Pairwise comparison revealed the following weights which were assigned to each of the four criteria: 0.674 for availability of hospital beds; 0.169 for the maximum road distance of the shortest routes; 0.101 for the elderly population within a 3‐km radius; and 0.056 for the median road distance of the shortest routes. The alternative proposed could cover 4000 more elderly people in the 3‐km radius of the hospitals. The integration of Geographical Information Systems and the Analytic Hierarchy Process constitutes a powerful tool for analysing traffic conditions in mid‐sized cities and for suggesting city planning to improve prognosis of stroke.     

Correlation between dissolved He concentration and Cl in groundwater at Äspö, Sweden

Tunnel excavation at Äspö Island, Sweden, has caused severe groundwater disturbance, gradually extending deeper into the tunnel as present-day Baltic seawater intrudes through fractures connecting to the surface. However, the paleo-hydrogeochemical conditions have remained in the deep highly saline waters that have avoided mixing. A correlation has been observed between dissolved ⁴He concentration and Cl⁻ ion concentration, measured every two years from 1995 to 2001 at Äspö. Groundwater mixing conditions can be examined by the correlations between 1/Cl, ³⁶Cl/Cl, and ³H concentrations. Subsurface production is responsible for the majority of the ³⁶Cl and excess dissolved ⁴He of interstitial groundwater in fractures. The secular equilibrium ratio of ³⁶Cl/Cl in rock was theoretically estimated to be (5.05 ± 0.82) × 10⁻¹⁴ based on the neutron flux intensity, a value comparable to the measured ³⁶Cl/Cl ratio in rock and groundwater. The degassing crustal ⁴He flux was estimated to be 2.9 × 10⁻⁸  1.3 × 10⁻⁶ (ccSTP/cm²a) using the HTO diffusion coefficient for the Äspö diorite. The ⁴He accumulation rate ranges from 6.8 × 10⁻¹⁰ (for the in situ accumulation rate) to 7.0 × 10⁻⁹ (ccSTP/(g_water · a) considering both ⁴He in situ production and the degassing flux, assuming ⁴He is accumulated constantly in groundwater. By comparing the subsurface ³⁶Cl increase with ⁴He concentrations in groundwater, the ⁴He accumulation rate was determined from data for groundwater arriving at the secular equilibrium of ³⁶Cl/Cl. The ⁴He accumulation rate was found to be (1.83 ± 0.72) × 10⁻⁸ ccSTP/(g_water · a) without determining the magnitude of degassing ⁴He flux.     

An integrated model simulating the initiation and motion of earthquake and rain induced rapid landslides and its application to the 2006 Leyte landslide

A gigantic rapid landslide claiming over 1,000 fatalities was triggered by rainfalls and a small nearby earthquake in the Leyte Island, Philippines in 2006. The disaster presented the necessity of a new modeling technology for disaster risk preparedness which simulates initiation and motion. This paper presents a new computer simulation integrating the initiation process triggered by rainfalls and/or earthquakes and the development process to a rapid motion due to strength reduction and the entrainment of deposits in the runout path. This simulation model LS-RAPID was developed from the geotechnical model for the motion of landslides (Sassa 1988) and its improved simulation model (Sassa et al. 2004b) and new knowledge obtained from a new dynamic loading ring shear apparatus (Sassa et al. 2004a). The examination of performance of each process in a simple imaginary slope addressed that the simulation model well simulated the process of progressive failure, and development to a rapid landslide. The initiation process was compared to conventional limit equilibrium stability analyses by changing pore pressure ratio. The simulation model started to move in a smaller pore pressure ratio than the limit equilibrium stability analyses because of progressive failure. However, when a larger shear deformation is set as the threshold for the start of strength reduction, the onset of landslide motion by the simulation agrees with the cases where the factor of safety estimated by the limit equilibrium stability analyses equals to a unity. The field investigation and the undrained dynamic loading ring shear tests on the 2006 Leyte landslide suggested that this landslide was triggered by the combined effect of pore water pressure due to rains and a very small earthquake. The application of this simulation model could well reproduce the initiation and the rapid long runout motion of the Leyte landslide.