Analytical Model for Permeability Evolution in Microcracking Rock

— Analytical expressions to predict the enhancement of permeability due to stress-induced microcracking in initially low porosity rock are presented. A fracture mechanical microcrack model is employed to derive integrated effective hydraulic variables as a function of stress, which are then used to calculate the evolution of permeability using the statistically-based Dienes model. The model enables determination of permeability enhancement as a function of two loading parameters and three material parameters. Results are in reasonable agreement with experimental measurements and indicate that appreciable increases in permeability can be anticipated during brittle failure. The analytical nature of the model makes it easily incorporatable into numerical models that require quantification of the permeability evolution as a function of stress, for which there is currently no law.     

Deep submarine pyroclastic eruptions: theory and predicted landforms and deposits

Submarine pyroclastic eruptions at depths greater than a few hundred meters are generally considered to be rare or absent because the pressure of the overlying water column is sufficient to suppress juvenile gas exsolution so that magmatic disruption and pyroclastic activity do not occur. Consideration of detailed models of the ascent and eruption of magma in a range of sea floor environments shows, however, that significant pyroclastic activity can occur even at depths in excess of 3000 m. In order to document and illustrate the full range of submarine eruption styles, we model several possible scenarios for the ascent and eruption of magma feeding submarine eruptions: (1) no gas exsolution; (2) gas exsolution but no magma disruption; (3) gas exsolution, magma disruption, and hawaiian-style fountaining; (4) volatile content builds up in the magma reservoir leading to hawaiian eruptions resulting from foam collapse; (5) magma volatile content insufficient to cause fragmentation normally but low rise speed results in strombolian activity; and (6) volatile content builds up in the top of a dike leading to vulcanian eruptions. We also examine the role of bulk-interaction steam explosivity and contact-surface steam explosivity as processes contributing to volcaniclastic formation in these environments. We concur with most earlier workers that for magma compositions typical of spreading centers and their vicinities, the most likely circumstance is the quiet effusion of magma with minor gas exsolution, and the production of somewhat vesicular pillow lavas or sheet flows, depending on effusion rate. The amounts by which magma would overshoot the vent in these types of eruptions would be insufficient to cause any magma disruption. The most likely mechanism of production of pyroclastic deposits in this environment is strombolian activity, due to the localized concentration of volatiles in magma that has a low rise rate; magmatic gas collects by bubble coalescence, and ascends in large isolated bubbles which disrupt the magma surface in the vent, producing localized blocks, bombs, and pyroclastic deposits. Another possible mode of occurrence of pyroclastic deposits results from vulcanian eruptions; these deposits, being characterized by the dominance of angular blocks of country rocks deposited in the vicinity of a crater, should be easily distinguishable from strombolian and hawaiian eruptions. However, we stress that a special case of the hawaiian eruption style is likely to occur in the submarine environment if magmatic gas buildup occurs in a magma reservoir by the upward drift of gas bubbles. In this case, a layer of foam will build up at the top of the reservoir in a sufficient concentration to exceed the volatile content necessary for disruption and hawaiian-style activity; the deposits and landforms are predicted to be somewhat different from those of a typical primary magmatic volatile-induced hawaiian eruption. Specifically, typical pyroclast sizes might be smaller; fountain heights may exceed those expected for the purely magmatic hawaiian case; cooling of descending pyroclasts would be more efficient, leading to different types of proximal deposits; and runout distances for density flows would be greater, potentially leading to submarine pyroclastic deposits surrounding vents out to distances of tens of meters to a kilometer. In addition, flows emerging after the evacuation of the foam layer would tend to be very depleted in volatiles, and thus extremely poor in vesicles relative to typical flows associated with hawaiian-style eruptions in the primary magmatic gas case. We examine several cases of reported submarine volcaniclastic deposits found at depths as great as 3000 m and conclude that submarine hawaiian and strombolian eruptions are much more common than previously suspected at mid-ocean ridges. Furthermore, the latter stages of development of volcanic edifices (seamounts) formed in submarine environments are excellent candidates for a wide range of submarine pyroclastic activity due not just to the effects of decreasing water depth, but also to: (1) the presence of a summit magma reservoir, which favors the buildup of magmatic foams (enhancing hawaiian-style activity) and episodic dike emplacement (which favors strombolian-style eruptions); and (2) the common occurrence of alkalic basalts, the CO₂ contents of which favor submarine explosive eruptions at depths greater than tholeiitic basalts. These models and predictions can be tested with future sampling and analysis programs and we provide a checklist of key observations to help distinguish among the eruption styles.     

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.     

Tensor time domain electromagnetic resistivity measurements at Ngatamariki geothermal field, New Zealand

Experimental measurements in the Ngatamariki geothermal field, North Island, New Zealand were made to test the applicability of the time domain electromagnetic method for detailed investigation of the resistivity structure within a geothermal field. Low-frequency square wave signals were transmitted through three grounded bipole current sources sited about 8 km from the measurement lines. Despite high levels of electrical noise, transient electric field vectors could be determined reliably for times between 0.02 and 3.3 s after each step in the source current. Instantaneous apparent resistivity tensors were then calculated. Apparent resistivity pseudosections along the two measurement lines show smooth variations of resistivity from site to site. Over most of the field the images consistently show a three-layer resistivity structure with a conductive middle layer (3–10 Ωm) representing the conductive upper part of the thermal reservoir. A deep-seated region of low resistivity in the northwest of the field may indicate a conductive structure at about 1 km associated with a deeper diorite intrusion. Measurements sited closer than about 100 m to drillholes appear to have been disturbed by metallic casing in the holes. A change in resistivity structure in the east of the field may indicate a major geological or hydrothermal boundary.     

Predictive instantaneous optimal control of elastic structures during earthquakes

A predictive instantaneous optimal control (PIOC) algorithm is proposed for controlling the seismic responses of elastic structures. This algorithm compensates for the time delay that happens in practical control applications by predicting the structural response over a period that equals the time delay, and by substituting the predicted response in the instantaneous optimal control (IOC) algorithm. The unique feature of this proposed PIOC algorithm is that it is simple and at the same time compensates for the time delay very effectively. Numerical examples of single degree of freedom structures are presented to compare the performance of PIOC and IOC systems for various time delay magnitudes. Results show that a time delay always causes degradation of control efficiency, but PIOC can greatly reduce this degradation compared to IOC. The effects of the structure's natural periods and the choice of control gains on the degradation induced by the time delay are also analyzed. Results show that shorter natural periods and larger control gains are both more sensitive and more serious to the degradation of control efficiency. Finally, a practical application of PIOC is performed on a six‐story moment‐resisting steel frame. It is demonstrated that PIOC contributes significantly to maintain stability in multiple degree of freedom structures, and at the same time PIOC has a satisfactory control performance. Copyright © 2003 John Wiley & Sons, Ltd.     

Lateral tidal mixing in the Antarctic marginal seas

Tidal mixing plays an important role in the modification of dense water masses around the Antarctic continent. In addition to the vertical (diapycnal) mixing in the near-bottom layers, lateral mixing can also be of relevance in some areas. A numerical tide simulation shows that lateral tidal mixing is not uniformly distributed along the shelf break. In particular, strong mixing occurs all along the Ross Sea and Southern Weddell Sea shelf breaks, while other regions (e.g., the western Weddell Sea) are relatively quiet. The latter regions correspond surprisingly well to areas where indications for cross-shelf exchange of dense water masses have been found. The results suggest that lateral tidal mixing may account for the relatively small contribution of Ross Sea dense water masses to Antarctic Bottom Water.     

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.     

In Situ Monitoring of Vapor Phase TCE Using a Chemiresistor Microchemical Sensor

A chemiresistor microchemical sensor has been developed to detect and monitor volatile organic compounds in unsaturated and saturated subsurface environments. A controlled study was conducted at the HAZMAT Spill Center at the Nevada Test Site, where the sensor was tested under a range of temperature, moisture, and trichloroethylene (TCE) concentrations. The sensor responded rapidly when exposed to TCE placed in sand, and it also responded to decreases in TCE vapor concentration when clean air was vented through the system. Variations in temperature and water vapor concentration impacted baseline chemiresistor signals, but at high TCE concentrations the sensor response was dominated by the TCE exposure. Test results showed that the detection limit of the chemiresistor to TCE vapor in the presence of fluctuating environmental variables (i.e., temperature and water vapor concentration) was on the order of 1000 parts per million by volume, which is about an order of magnitude higher than values obtained in controlled laboratory environments. Automated temperature control and preconcentration is recommended to improve the stability and sensitivity of the chemiresistor sensor.     

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.