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.     

Induced Microearthquake Patterns in Hydrocarbon and Geothermal Reservoirs: Six Case Studies

—?The injection or production of fluids can induce microseismic events in hydrocarbon and geothermal reservoirs. By deploying sensors downhole, data sets have been collected that consist of a few hundred to well over 10,000 induced events. We find that most induced events cluster into well-defined geometrical patterns. In many cases, we must apply high-precision, relative location techniques to observe these patterns. At three sedimentary sites, thin horizontal strands of activity are commonly found within the location patterns. We believe this reflects fracture containment between stratigraphic layers of differing mechanical properties or states of stress. At a massive carbonate and two crystalline sites, combinations of linear and planar features indicate networks of intersecting fractures and allow us to infer positions of aseismic fractures through their influence on the location patterns. In addition, the fine-scale seismicity patterns often evolve systematically with time. At sedimentary sites, migration of seismicity toward the injection point has been observed and may result from slip-induced stress along fractures that initially have little resolved shear. In such cases, triggering events may be critical to generate high levels of seismic activity. At one crystalline site, the early occurrence of linear features that traverse planes of activity indicate permeable zones and possible flow paths within fractures. We hope the continued development of microseismic techniques and refinement of conceptual models will further increase our understanding of fluid behavior and lead to improved resource management in fractured reservoirs.     

Multiparametric Traveltime Inversion

In conventional seismic processing, the classical algorithm of Hubral and Krey is routinely applied to extract an initial macrovelocity model that consists of a stack of homogeneous layers bounded by curved interfaces. Input for the algorithm are identified primary reflections together with normal moveout (NMO) velocities, as derived from a previous velocity analysis conducted on common midpoint (CMP) data. This work presents a modified version of the Hubral and Krey algorithm that is designed to extend the original version in two ways, namely (a) it makes an advantageous use of previously obtained common-reflection-surface (CRS) attributes as its input and (b) it also allows for gradient layer velocities in depth. A new strategy to recover interfaces as optimized cubic splines is also proposed. Some synthetic examples are provided to illustrate and explain the implementation of the method.     

Reduction of the Trophic Status of a Water Body with the Use of Large-Size Zooplankton

It is shown that large-size zooplankton can be used as a natural filter in flow-through plants for reducing the trophic status of water bodies and improve water transparency. Experiments in Chistye Prudy pond in the period of blue-green algae blooming show large-size zooplankton to reduce phytoplankton biomass in the pond throughout the period of phytoplankton growth. Daphnia magna specimens (with a biomass of 92 to 3450 mg/l) in lothic ecosystems are found to reduce the biomass of dominating Anabaena spirodes blue-green algae when this biomass attains its maximum.