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.     

System identification of the Vincent Thomas suspension bridge using earthquake records

The Vincent Thomas Bridge in the Los Angeles metropolitan area, is a critical artery for commercial traffic flow in and out of the Los Angeles Harbor, and is at risk in the seismically active Southern California region, particularly because it straddles the Palos Verdes fault zone. A combination of linear and non‐linear system identification techniques is employed to obtain a complete reduced‐order, multi‐input–multi‐output (MIMO) dynamic model of the Vincent Thomas Bridge based on the dynamic response of the structure to the 1987 Whittier and 1994 Northridge earthquakes. Starting with the available acceleration measurements (which consists of 15 accelerometers on the bridge structure and 10 accelerometers at various locations on its base), an efficient least‐squares‐based time‐domain identification procedure is applied to the data set to develop a reduced‐order, equivalent linear, multi‐degree‐of‐freedom model. Although not the main focus of this study, the linear system identification method is also combined with a non‐parametric identification technique, to generate a reduced‐order non‐linear mathematical model suitable for use in subsequent studies to predict, with good fidelity, the total response of the bridge under arbitrary dynamic environments. Results of this study yield measurements of the equivalent linear modal properties (frequencies, mode shapes and non‐proportional damping) as well as quantitative measures of the extent and nature of non‐linear interaction forces arising from strong ground shaking. It is shown that, for the particular subset of observations used in the identification procedure, the apparent non‐linearities in the system restoring forces are quite significant, and they contribute substantially to the improved fidelity of the model. Also shown is the potential of the identification technique under discussion to detect slight changes in the structure's influence coefficients, which may be indicators of damage and degradation in the structure being monitored. Difficulties associated with accurately estimating damping for lightly damped long‐span structures from their earthquake response are discussed. The technical issues raised in this paper indicate the need for added spatial resolution in sensor instrumentation to obtain identified mathematical models of structural systems with the broadest range of validity. Copyright © 2003 John Wiley & Sons, Ltd.     

Observing artificially induced strain: implications for subglacial deformation

By analysing a series of four successive thin‐sections from a ceramic clay that was subjected to uniaxial compression, we were able to monitor the development of microstructures in a fine‐grained sediment. The artificially induced microstructures, such as unidirectional clay reorientations and linear and circular grain arrangements, are identical to features that have been observed in thin‐sections of subglacially deformed tills, and therefore may be used as representative analogues. We argue that the structures, reflecting slip, planar shear displacements as well as rotational movements, can be explained by assuming a Coulomb‐plastic response to imposed shear. We conclude that sediments subjected to subglacial deformation behave as Coulomb materials, at least during the final stages of the deformation. The present study bridges the gap between field studies, experimental studies and theoretical modelling. The microscopic observations assist in visualising inferred subglacial processes and facilitate up‐ and downscaling between diverse methodological approaches. Copyright © 2003 John Wiley & Sons, Ltd.