A cellular model of Holocene upland river basin and alluvial fan evolution

The CAESAR (Cellular Automaton Evolutionary Slope And River) model is used to simulate the Holocene development of a small upland catchment (4·2 km²) and the alluvial fan at its base. The model operates at a 3 m grid scale and simulates every flood over the last 9200 years, using a rainfall record reconstructed from peat bog wetness indices and land cover history derived from palynological sources. Model results show that the simulated catchment sediment discharge above the alluvial fan closely follows the climate signal, but with an increase in the amplitude of response after deforestation. The important effects of sediment storage and remobilization are shown, and findings suggest that soil creep rates may be an important control on long term (>1000 years) temperate catchment sediment yield. The simulated alluvial fan shows a complex and episodic behaviour, with frequent avulsions across the fan surface. However, there appears to be no clear link between fan response and climate or land use changes suggesting that Holocene alluvial fan dynamics may be the result of phases of sediment storage and remobilization, or instabilities and thresholds within the fan itself. Copyright © 2002 John Wiley & Sons, Ltd.     

Effective stress analyses of liquefaction-induced deformation in river dikes

A series of effective stress analyses is carried out on the seismic performance of river dikes based on the case histories during the 1993 Hokkaido-Nansei-oki and 1995 Hyogoken-Nambu earthquakes in Japan. Seven case histories selected for the analyses involve a crest settlement ranging from none to 2.7 m in the dikes 3–6 m high with evidence of liquefaction at foundation soil. The effective stress model used is based on a multiple shear mechanism and was developed by one of the authors. The soil parameters are evaluated based on the site investigation and laboratory test results. The results of the analyses are basically consistent with the observed performance of the river dikes. In particular, the effective stress model shows a reasonable capability to reproduce the varying degree of settlements depending on the geotechnical conditions of foundation soils beneath the dikes. The analyses also indicate that the effect of a cohesive soil layer mixed with the liquefiable sand layers beneath the dikes can be a primary factor for reducing the liquefaction-induced deformation of dikes.     

Production of selected cosmogenic radionuclides by muons: 2. Capture of negative muons

We have determined the production yields for radionuclides in Al₂O₃, SiO₂, S, Ar, K₂SO₄, CaCO₃, Fe, Ni and Cu targets, which were irradiated with slow negative muons at the Paul Scherrer Institute in Villigen (Switzerland). The fluences of the stopped negative muons were determined by measuring the muonic X-rays. The concentrations of the long-lived and short-lived radionuclides were measured with accelerator mass spectrometry (AMS) and γ-spectroscopy, respectively. Special emphasis was put on the radionuclides ¹⁰Be, ¹⁴C and ²⁶Al produced in quartz targets, ²⁶Al in Al₂O₃ and S targets, ³⁶Cl in K₂SO₄ and CaCO₃ targets, and ⁵³Mn in Fe₂O₃ targets. These targets were selected because they are also the naturally occurring target minerals for cosmic ray interactions in typical rocks. We also present results of calculations for depth-dependent production rates of radionuclides produced after cosmic ray μ⁻ capture, as well as cosmic ray-induced production rates of geologically relevant radionuclides produced by the nucleonic component, by μ⁻ capture, by fast muons and by neutron capture.     

A rifted inside corner massif on the Mid-Atlantic Ridge at 5°S

The structure of the Mid-Atlantic Ridge at 5°S was investigated during a recent cruise with the FS Meteor. A major dextral transform fault (hereafter the 5°S FZ) offsets the ridge left-laterally by 80 km. Just south of the transform and to the west of the median valley, the inside corner (IC – the region bounded by the ridge and the active transform) is marked by a major massif, characterized by a corrugated upper surface. Fossil IC massifs can also be identified further to the west. Unusually, a massif almost as high as the IC massif also characterizes the outside corner (OC) south of the inactive fracture zone and to the east of the median valley. This OC massif has axis-parallel dimensions identical to the IC massif and both are bounded on their sides closest to the spreading axis by abrupt, steep slopes. An axial volcanic ridge is well developed in the median valley both south of the IC/OC massifs and in an abandoned rift valley to the east of the OC massif, but is absent along the new ridge-axis segment between the IC and OC massifs. Wide-angle seismic data show that between the massifs, the crust of the median valley thins markedly towards the FZ. These observations are consistent with the formation of the OC massif by the rifting of an IC core complex and the development of a new spreading centre between the IC and OC massifs. The split IC massif presents an opportunity to study the internal structure of the footwall of a detachment fault, from the corrugated fault surface to deeper beneath the fault, without recourse to drilling. Preliminary dredging recovered gabbros from the scarp slope of the rifted IC massif, and serpentinites and gabbros from the intersection of this scarp with the corrugated surface. This is compatible with a concentration of serpentinites along the detachment surface, even where the massif internally is largely plutonic in nature.     

Shallow seismicity, stress distribution and crustal deformation pattern in the Andaman-West Sunda arc and Andaman Sea, northeastern Indian Ocean

Shallow seismicity and available source mechanisms in the Andaman–westSunda arc and Andaman sea region suggest distinct variation in stressdistribution pattern both along and across the arc in the overriding plate.Seismotectonic regionalisation indicates that the region could be dividedinto eight broad seismogenic sources of relatively homogeneousdeformation. Crustal deformation rates have been determined for each oneof these sources based on the summation of moment tensors. The analysisshowed that the entire fore arc region is dominated by compressive stresseswith compression in a mean direction of N23^°, and the rates ofseismic deformation velocities in this belt decrease northward from 5.2± 0.65 mm/yr near Nias island off Sumatra and 1.12 ±0.13 mm/yr near Great Nicobar islands to as much as 0.4 ±0.04 mm/yr north of 8^°N along Andaman–Nicobar islandsregion. The deformation velocities indicate, extension of 0.83 ±0.05 mm/yr along N343^° and compression of 0.19 ±0.01 mm/yr along N73^° in the Andaman back arc spreadingregion, extension of 0.18 ± 0.01 mm/yr along N125^° andcompression of 0.16 ± 0.01 mm/yr along N35^° in NicobarDeep and west Andaman fault zone, compression of 0.84 ±0.12 mm/yr N341^° and extension of 0.77 ± 0.11 mm/yralong N72^° within the transverse tectonic zone in the Andamantrench, N-S compression of 3.19 ± 0.29 mm/yr and an E-Wextension of 1.24 ± 0.11 mm/yr in the Semangko fault zone ofnorth Sumatra. The vertical deformation suggests crustal thinning in theAndaman sea and crustal thickening in the fore arc and Semangko faultzones. The apparent stresses calculated for all major events range between0.1–10 bars and the values increase with increasing seismic moment.However, the apparent stress estimates neither indicate any significantvariation with faulting type nor display any variation across the arc, incontrast to the general observation that the fore arc thrust events showhigher stress levels in the shallow subduction zones. It is inferred that theoblique plate convergence, partial subduction of 90^°E Ridge innorth below the Andaman trench and the active back arc spreading are themain contributing factors for the observed stress field within the overridingplate in this region.     

The Influence of 3-D Structure on the Propagation of Seismic Waves Away from Earthquakes

-- The seismic records from significant earthquakes are profoundly affected by 3-D variations in crustal structure both in the source zone itself and in propagation to some distance. Even in structurally complex zones such as Japan and Mexico relatively coherent arrivals are found associated with different classes of propagation paths. The presence of strong lateral variations can disrupt the arrivals, and impose significant variations in propagation characteristics for different directions from the source as illustrated by observations for the 1995 Kobe and 2000 Tottori-ken Seibu earthquakes in western Japan. Such effects can be modelled in 3 dimensions using a hybrid scheme with a pseudospectral representation for horizontal coordinates and finite differences in depth. This arrangement improves parallel implementation by minimising communication costs. For a realistic 3-D model for the structure in western Japan the 3-D simulations to frequencies close to 1 Hz provide a good representation of the observations from subduction zones events such as the 1946 Nankai earthquake and the 2000 Tottori-ken Seibu earthquake. The model can therefore be used to investigate the pattern of ground motion expected for future events e.g., in current seismic gaps.     

Efficacy of Annual Bacteria Monitoring and Shock Chlorination in Wells Finished in a Floodplain Aquifer

Public health authorities generally recommend annual water-quality monitoring of rural water wells and shock chlorination if coliforms are detected. It is implicitly assumed that shock chlorination is effective in ridding most wells of bacteriological pathogens for months to years. Neither annual monitoring nor shock chlorination was effective in addressing coliform contamination of selected water wells in a small town developed on an alluvial aquifer where septic system effluents are impacting well water quality. Considerable temporal variation in total and fecal coliforms was observed in water wells monitored for a six-month period. Individual wells intermittently met and exceeded the drinking water criteria, indicating annual sampling was insufficient. Shock chlorination of three contaminated wells and their associated distribution systems proved ineffective because colonies apparently originated from outside the wells and reappeared over relatively short time periods (ranging from less than one week up to 21 weeks). The relatively fast and similar rate of recovery of total heterotrophic bacteria suggested they are related to biofilm formation in the wells and not to ground water contamination.     

Surface vibration due to a sequence of high speed moving harmonic rectangular loads

The transmission of vibrations over the surface of the ground, due to high-speed moving, vertical harmonic rectangular loads, is investigated theoretically. The problem is three-dimensional and the interior of the ground is modelled as an elastic half-space or a multilayered ground. The transformed solutions are obtained using the Fourier transform on the space variable. A new damping model in the spatial wavenumber domain, presented in Lefeuve-Mesgouez et al. [J. Sound. Vibr. 231 (2000) 1289] is used. Numerical results for the displacements on the surface are presented for loads moving with speeds up to and beyond the Rayleigh wave speed of the half-space.     

Optimized Seismic Threshold Monitoring – Part 2: Teleseismic Processing

v--vThis second paper (Part 2) pertaining to optimized site-specific threshold monitoring addresses the application of the method to regions covered by a teleseismic or a combined regional-teleseismic network. In the first paper (Part 1) we developed the method for the general case, and demonstrated its application to an area well-covered by a regional network (the Novaya Zemlya nuclear test site). In the present paper, we apply the method to the Indian and Pakistani nuclear test sites, and show results during the periods of nuclear testing by these two countries in May 1998. Since the coverage by regional stations in these areas is poor, an optimized approach requires the use of selected, high-quality stations at teleseismic distances.¶To optimize the threshold monitoring of these test sites, we use as calibration events either one of the nuclear explosions or a nearby earthquake. From analysis of the calibration events we derive values for array beamforming steering delays, filter bands, short-term averages (STA) lengths, phase travel times (P waves), and amplitude-magnitude relationships for each station. By applying these parameters, we obtain a monitoring capability of both test sites ranging from m_b 2.8-3.0 using teleseismic stations only. When including the nearby Nilore station to monitor the Indian tests, we show that the threshold can be reduced by about 0.4 magnitude units. In particular, we demonstrate that the Indian tests on 13 May, 1998, which were not detected by any known seismic station, must have corresponded to a magnitude (m_b) of less than 2.4.¶We also discuss the effect of a nearby aftershock sequence on the monitoring capability for the Pakistani test sites. Such an aftershock sequence occurred in fact on the day of the last Pakistani test (30 May, 1998), following a large (m_b 5.5) earthquake in Afghanistan located about 1100 km from the test site. We show that the threshold monitoring technique has sufficient resolution to suppress the signals from these interfering aftershocks without significantly affecting the true peak of the nuclear explosion on the threshold trace.     

3-D Viscoelastic FEM Modeling of Crustal Deformation in Northeast Japan

-- As a first step toward establishing a standard earthquake cycle model in Japan, we simulate the crustal deformation during the past 100 years in northeast Japan, using a 3-D FEM based on the kinematic model. Then, we compare the computed results with the observed long-term leveling data and the recent GPS data. On the whole, although the effect of the subducting PAC is dominant, coseismic deformation of the interplate earthquakes can be clearly seen in the inland. Moreover, the postseismic deformation of the earthquakes due to the viscoelastic upper mantle seriously affects the inland movements, and continues for a few decades. Our modeling, including the effects of the interplate earthquakes and the three-dimensional viscoelastic inhomogeneity, reasonably explains the observed movement. Finally, we stress that the viscoelastic effect should be taken into consideration in the analyses, even if no earthquakes occur in the analyzed period.