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.     

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.     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

Combined GPS and InSAR Models of Postseismic Deformation from the Northridge Earthquake

--Models of combined Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR) data collected in the region of the Northridge earthquake indicate that significant afterslip on the main fault occurred following the earthquake. Additional shallow deformation occurred to the west of the main rupture plane. Both data sets are consistent with logarithmic time-dependent behavior following the earthquake indicative of afterslip rather than postseismic relaxation. Aftershocks account for only about 10% of the postseismic motion. The two data sets are complimentary in determining the postseismic processes. Fault afterslip and shallow deformation dominate the deformation field in the two years following the earthquake. Lower crustal deformation may play an important role later in the earthquake cycle.     

Identification of Mining Blasts at Mid- to Far-regional Distances Using Low Frequency Seismic Signals

—?This paper reports results from two recent monitoring experiments in Wyoming. Broadband seismic recordings of kiloton class delay-fired cast blasts and instantaneous calibration shots in the Black Thunder coal mine were made at four azimuths at ranges from 1° to 2°. The primary focus of this experiment was to observe and to explain low-frequency signals that can be seen at all azimuths and should routinely propagate above noise to mid-regional distances where most events will be recorded by International Monitoring System (IMS) stations.¶The recordings clearly demonstrate that large millisecond delay-fired cast blasts routinely produce seismic signals that have significant spectral modulations below 10?Hz. These modulations are independent of time, the azimuth from the source and the orientation of the sensor. Low-frequency modulations below 5?Hz are seen beyond 9°. The modulations are not due to resonance as they are not produced by the calibration shots. Linear elastic modeling of the blasts that is guided by mine-blast reports fails to reproduce the fine detail of these modulations but clearly indicates that the enhanced “spectral roughness” is due to long interrow delays and source finiteness. The mismatch between the data and the synthetics is likely due to source processes, such as nonlinear interactions between shots, that are poorly understood and to other effects, such as variations of shot time and yield from planned values, that are known to be omnipresent but cannot be described accurately. A variant of the Automated Time-Frequency Discriminant (Hedlin, 1998b), which uses low-frequency spectral modulations, effectively separates these events from the calibration shots.¶The experiment also provided evidence that kiloton class cast blasts consistently yield energetic 2–10 second surface waves. The surface waves are strongly dependent on azimuth but are seen beyond 9°. Physical modeling of these events indicates that the surface waves are due mainly to the extended source duration and to a lesser extent to the slap-down of spalled material. The directionality is largely a path effect. A discriminant that is based on the partitioning of energy between surface and body waves routinely separates these events from the calibration shots.¶The Powder River Basin has essentially no natural seismic activity. How these mining events compare to earthquake observations remains to be determined.     

Particle recycling in volcanic plumes

We have developed a new theoretical model of an eruption column that accounts for the re-entrainment of particles as they fall out of the laterally spreading umbrella cloud. The model illustrates how the mass flux of particles in the plume may increase with height in the plume, by a factor as large as 2.5 because of this recycling. Three important consequences are that (1) the critical velocity required to generate a buoyant eruption column for a given mass flux increases, (2) the total height of rise of the column may decrease, and (3) we infer that in relatively wind-free environments, for eruption columns near the conditions for collapse, the recycling of particles may lead to an unsteady oscillating motion of the plume, which, in time, may lead to the formation of interleaved fall and flow deposits.     

Modelling Seismic Waves Around Underground Openings in Fractured Rock

—?The potential for large excavation-induced seismic events may be recognised, even if the timing of an event may be inherently unpredictable. In this case, modelling the wave propagation from a potential event could allow the dynamic motions around an excavation to be projected, and for areas of danger to be anticipated. However, the above and other potential applications require accurate models of wave interaction with the openings, as well as with the fractured rock which surrounds such excavations. This paper considers real recorded waveforms and how well these waveforms are modelled by explicit mechanical models of the source, the medium and the excavation. Models of experiments at three different scales of the problem are presented: small and large amplitude waveforms recorded around a deep-level mining tunnel in a synthetic rockburst experiment; waveforms from laboratory experiments of waves through plates of steel representing fractures; waveforms from active pulses in an acoustic emission experiment in a small volume of fractured rock at the surface of an underground excavation. The results show that elastic wave propagation around an excavation was a first approximation for small amplitude waves, but was less successful for modelling large amplitude waves and more fractured rock. Fractures in the models were represented explicitly with displacement discontinuities. Waveforms through known fracture geometries were particularly well-reproduced, and indicate the importance of fracture stiffness, the in situ stress state, and stress-dependence of the fractures in such models. Overall, the models are sufficiently successful at representing recorded behaviour, to be encouraging for the goal of representing accurate wave motions around excavations.     

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.     

"EFFECTS OF FLOW INTENSITY, OBSTACLE ALIGNMENT AND CROSS-SECTION GEOMETRY ON SCOUR AT BRIDGE ABUTMENTS" CARDOSO, A.H., SANTOS, J.S. and ROCA, M.

The authors present new measurements of maximum (equilibrium) local scour depth at bridge abutments, the new data being analysed in terms of the multiplying factors given in Melville (1992,1997) and Melville and Coleman (2000). The design method presented in these publications is referred to here as the University of Auckland (UoA) Method. The authors' data apply to the following factors:∏f (flow intensity), ∏θ (abutment alignment) and ∏g (approach channel geometry).     

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.