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.     

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.     

Predictive instantaneous optimal control of inelastic structures during earthquakes

A new inelastic structural control algorithm is proposed by incorporating the force analogy method (FAM) with the predictive instantaneous optimal control (PIOC) algorithm. While PIOC is very effective in compensating for the time delay for elastic structures, the FAM is highly efficient in performing the inelastic analysis. Unlike conventional inelastic analysis methods of changing stiffness, the FAM analyzes structures by varying the structural displacement field, and therefore the state transition matrix needs to be computed only once. This greatly simplifies the computation and makes inelastic analysis readily applicable to the PIOC algorithm. The proposed algorithm compensates for the time delay that happens in practical control systems by predicting the inelastic structural response over a period that equals the magnitude of the time delay. A one‐story frame with both strain‐hardening and strain‐softening inelastic characteristics is analyzed using this algorithm. Results show that the proposed control algorithm is feasibile for any inelastic structures. While the control efficiency deteriorates with the increase in magnitude of the time delay, the PIOC maintains acceptable performance within a wide range of time delay magnitudes. Finally, a computer model of a six‐story moment‐resisting steel frame is analyzed to show that PIOC has good control results for real inelastic structures. 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.     

Post-caldera activity in the Alban Hills volcanic district (Italy): 40Ar/39Ar geochronology and insights into magma evolution

We present 24 ⁴⁰Ar/³⁹Ar ages for the youngest volcanic products from the Alban Hills volcanic district (Rome). Combined with petrological data on these products, we have attempted to define the chronology of the most recent phase of activity and to investigate the magma evolution of this volcanic district. The early, mainly explosive activity of the Alban Hills spanned the interval from 561±1 to 351±3 ka. After approximately 50-kyr of dormancy, a mainly effusive phase of activity took place, accompanied by the strombolian activity of a small central edifice (Monte delle Faete). This second phase of activity spanned the interval 308±2 to 250±1 ka. After another dormancy period of approximately 50-kyr, a new, hydromagmatic phase of activity started at 200 ka at several centers located to the southwest of the Monte delle Faete edifice. After an initial recurrence period of approximately 50-kyr, which also characterized this new phase of activity, the longest dormancy period (approximately 80-kyr) in the history of the volcanic district preceded the start of the activity of the Albano and Giuturna centers at 70±1 ka. Results of our study suggest a quasi-continuous magmatic activity feeding hydromagmatic centers with a new acme of volcanism since around 70 ka. Based on data presented in this paper, we argue that the Alban Hills should not be considered an extinct volcanic district and a detailed re-assessment of the volcanic hazard for the area of Rome is in order. Published online: 4 April 2003 Editorial responsibility: J. Donnelly-Nolan     

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.