Borehole Breakouts in Berea Sandstone Reveal a New Fracture Mechanism

— Vertical drilling experiments in high-porosity (22% and 25%) Berea sandstone subjected to critical true triaxial far-field stresses, in which σ _H (maximum horizontal stress) >σ _v (vertical stress) >σ _h (least horizontal stress), revealed a new and non-dilatant failure mechanism that results in thin and very long tabular borehole breakouts that have the appearance of fractures, and which counterintuitively develop orthogonally to σ _H . These breakouts are fundamentally different from those induced in crystalline rocks, as well as limestones and medium-porosity Berea sandstone. Breakouts in these rocks are typically dog-eared in shape, a result of dilatant multi-cracking tangential to the hole and subparallel to the maximum far-field horizontal stress σ _H , followed by progressive buckling and shearing of detached rock flakes created by the cracks. In the high-porosity sandstone a narrow layer of grains compacted normal to σ _H is observed just ahead of the breakout tip. This layer is nearly identical to “compaction bands” observed in the field. It is suggested that when a critical tangential stress concentration is reached along the σ _h spring line at the borehole wall, grain bonding breaks down and a compaction band is formed normal to σ _H . Debonded loose grains are expelled into the borehole, assisted by the circulating drilling fluid. As the breakout tip advances, the stress concentration ahead of it persists or may even increase, extending the compaction band, which in turn leads to breakout lengthening.     

Earthen barriers to control lava flows in the 2001 eruption of Mt. Etna

Preceded by four days of intense seismicity and marked ground deformation, a new eruption of Mt. Etna started on 17 July and lasted until 9 August 2001. It produced lava emission and strombolian and phreatomagmatic activity from four different main vents located on a complex fracture system extending from the southeast summit cone for about 4.5 km southwards, from 3000 to 2100 m elevation (a.s.l.). The lava emitted from the lowest vent cut up an important road on the volcano and destroyed other rural roads and a few isolated country houses. Its front descended southwards to about 4 km distance from the villages of Nicolosi and Belpasso. A plan of intervention, including diversion and retaining barriers and possibly lava flow interruption, was prepared but not activated because the flow front stopped as a consequence of a decrease in the effusion rate. Extensive interventions were carried out in order to protect some important tourist facilities of the Sapienza and Mts. Silvestri zones (1900 m elevation) from being destroyed by the lava emitted from vents located at 2700 m and 2550 m elevation. Thirteen earthen barriers (with a maximum length of 370 m, height of 10–12 m, base width of 15 m and volume of 25 000 m³) were built to divert the lava flow away from the facilities towards a path implying considerably less damage. Most of the barriers were oriented diagonally (110–135°) to the direction of the flow. They were made of loose material excavated nearby and worked very nicely, resisting the thrust of the lava without any difficulty. After the interventions carried out on Mt. Etna in 1983 and in 1991–1992, those of 2001 confirm that earthen barriers can be very effective in controlling lava flows.     

Steel frames with bracing mechanism and hysteretic dampers

A new earthquake resistant structural system for multi‐storey frame structures, based on a dual function of its bracing components, is developed. This consists of a hysteretic damper device and a cross‐bracing mechanism with a kinetic closed circuit, working only in tension, so that cable members can be used for this purpose. Solutions are presented regarding the connections' design of three types of structural frame system, that are concerned throughout the study: braced moment free frame, braced moment resisting frame with moment free supports, and with moment resisting supports. The dynamic behaviour of the system is investigated on the basis of an SDOF model, and based on the response spectra method an approximate design approach of the controlled structures is shown. From the time history analysis of the structural systems for the El Centro earthquake the areas of appropriate stiffness relations of the frames to the hysteretic dampers and the cable braces are deduced, so that the energy dissipation of the system may be controlled by the damper‐cable bracing mechanism. Based on the results of these studies, a predesign approach is developed for the implementation of the control system in frame structures. Copyright © 2003 John Wiley & Sons, Ltd.     

An experimental study of nominally symmetric and asymmetric structures isolated with the FPS

This experimental investigation deals with the earthquake behaviour of a nominally symmetric and a mass‐asymmetric three‐storey structural model isolated with the frictional pendulum system (FPS). Both accidental and natural torsion are evaluated in the structure by using recorded accelerations in all building floors and measured deformations at the isolation level. A 3D‐shaking table was used to subject the model to five different ground motions, including impulsive as well as far‐field subduction‐zone type earthquakes. Results show that the analytical predictions of the earthquake behaviour of the isolated structure, as obtained from a physical model of the FPS, are in close agreement with the true complex inelastic measured behaviour of the FPS. Besides, experimental results also validate previous observations about the importance of accounting for the variability of the normal loads in modelling the earthquake behaviour of FPS isolators. Measured torsional deformation amplifications at the base of the building vary, in the mean, from 2.5% to 6% for the symmetric and asymmetric structural configurations, respectively. In relation to the fixed base structure, the reduction factors for the base shear of the isolated structure are, in the mean, about 3.9 for both configurations. Finally, it is concluded that the FPS is capable of controlling the lateral–torsional motions of mass‐asymmetric structures quite effectively by aligning the centre of mass of the superstructure with the centre of pendular and frictional resistance of the isolation system. Copyright © 2003 John Wiley & Sons, Ltd.     

Application of DARLAM to Regional Haze Modeling

— The CSIRO Division of Atmospheric Research limited area model (DARLAM) is applied to atmospheric transport modeling of haze in southeast Asia. The 1998 haze episode is simulated using an emission inventory derived from hotspot information and adopting removal processes based on SO₂.¶Results show that the model is able to simulate the transport of haze in the region. The model images closely resemble the plumes of NASA Total Ozone Mapping Spectrometer and Meteorological Service Singapore haze maps. Despite the limitation of input data, particularly for haze emissions, the three-month average pattern correlation obtained for the whole episode is 0.61. The model has also been able to reproduce the general features of transboundary air pollution over a long period of time. Predicted total particulate matter concentration also agrees reasonably well with observation.¶The difference in the model results from the satellite images may be attributed to the large uncertainties of emission, simplification of haze deposition and transformation mechanisms and the relatively coarse horizontal and vertical resolution adopted for this particular simulation.     

A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

Field Survey and Numerical Simulations: A Review of the 1998 Papua New Guinea Tsunami

— The Papua New Guinea (PNG) tsunami of 1998 is re-examined through a detailed review of the field survey as well as numerous numerical computations. The discussion of the field survey explores a number of possible misinterpretations of the recorded data. The survey data are then employed by a numerical model as a validation tool. A Boussinesq model and a nonlinear shallow water wave (NLSW) model are compared in order to quantify the effect of frequency dispersion on the landslide-generated tsunami. The numerical comparisons indicate that the NLSW model is a poor estimator of offshore wave heights. However, due to what appears to be depth-limited breaking seaward of Sissano spit, both numerical models are in agreement in the prediction of maximum water elevations at the overtopped spit. By comparing three different hot-start initial profiles of the tsunami wave, it is shown that the initial shape and orientation of the tsunami wave is secondary to the initial displaced water mass in regard to prediction of water elevations on the spit. These numerical results indicate that agreement between numerical prediction of runup values with field recorded values at PNG cannot be used to validate either a NLSW tsunami propagation model or a specific landslide tsunami hot-start initial condition. Finally, with the use of traditional tsunami codes, a new interpretation of the PNG runup measurements is presented.