U–Pb zircon chronostratigraphy of early-Pliocene ignimbrites from La Pacana, north Chile: implications for the formation of stratified magma chambers

High spatial resolution U–Pb dates of zircons from two consanguineous ignimbrites of contrasting composition, the high-silica rhyolitic Toconao and the overlying dacitic Atana ignimbrites, erupted from La Pacana caldera, north Chile, are presented in this study. Zircons from Atana and Toconao pumice clasts yield apparent ²³⁸U/²⁰⁶Pb ages of 4.11±0.20 Ma and 4.65±0.13 Ma (2σ), respectively. These data combined with previously published geochemical and stratigraphic data, reveal that the two ignimbrites were erupted from a stratified magma chamber. The Atana zircon U–Pb ages closely agree with the eruption age of Atana previously determined by K–Ar dating (4.0±0.1 Ma) and do not support long (>1 Ma) residence times. Xenocrystic zircons were found only in the Toconao bulk ignimbrite, which were probably entrained during eruption and transport. Apparent ²³⁸U/²⁰⁶Pb zircon ages of 13 Ma in these xenocrysts provide the first evidence that the onset of felsic magmatism within the Altiplano–Puna ignimbrite province occurred approximately 3 Myr earlier than previously documented.     

Exploration and discovery in Yellowstone Lake: results from high-resolution sonar imaging, seismic reflection profiling, and submersible studies

‘No portion of the American continent is perhaps so rich in wonders as the Yellow Stone’ (F.V. Hayden, September 2, 1874)Discoveries from multi-beam sonar mapping and seismic reflection surveys of the northern, central, and West Thumb basins of Yellowstone Lake provide new insight into the extent of post-collapse volcanism and active hydrothermal processes occurring in a large lake environment above a large magma chamber. Yellowstone Lake has an irregular bottom covered with dozens of features directly related to hydrothermal, tectonic, volcanic, and sedimentary processes. Detailed bathymetric, seismic reflection, and magnetic evidence reveals that rhyolitic lava flows underlie much of Yellowstone Lake and exert fundamental control on lake bathymetry and localization of hydrothermal activity. Many previously unknown features have been identified and include over 250 hydrothermal vents, several very large (>500 m diameter) hydrothermal explosion craters, many small hydrothermal vent craters (1–200 m diameter), domed lacustrine sediments related to hydrothermal activity, elongate fissures cutting post-glacial sediments, siliceous hydrothermal spire structures, sublacustrine landslide deposits, submerged former shorelines, and a recently active graben. Sampling and observations with a submersible remotely operated vehicle confirm and extend our understanding of the identified features. Faults, fissures, hydrothermally inflated domal structures, hydrothermal explosion craters, and sublacustrine landslides constitute potentially significant geologic hazards. Toxic elements derived from hydrothermal processes also may significantly affect the Yellowstone ecosystem.     

Asymmetric one‐storey elastic systems with non‐linear viscous and viscoelastic dampers: Simplified analysis and supplemental damping system design

Investigated is the accuracy in estimating the response of asymmetric one‐storey systems with non‐linear viscoelastic (VE) dampers by analysing the corresponding linear viscous system wherein all non‐linear VE dampers are replaced by their energy‐equivalent linear viscous dampers. The response of the corresponding linear viscous system is determined by response history analysis (RHA) and by response spectrum analysis (RSA) extended for non‐classically damped systems. The flexible and stiff edge deformations and plan rotation of the corresponding linear viscous system determined by the extended RSA procedure is shown to be sufficiently accurate for design applications with errors generally between 10 and 20%. Although similar accuracy is also shown for the ‘pseudo‐velocity’ of non‐linear VE dampers, the peak force of the non‐linear VE damper cannot be estimated directly from the peak damper force of the corresponding linear viscous system. A simple correction for damper force is proposed and shown to be accurate (with errors not exceeding 15%). For practical applications, an iterative linear analysis procedure is developed for determining the amplitude‐ and frequency‐dependent supplemental damping properties of the corresponding linear viscous system and for estimating the response of asymmetric one‐storey systems with non‐linear VE dampers from the earthquake design (or response) spectrum. Finally, a procedure is developed for designing non‐linear supplemental damping systems that satisfy given design criteria for a given design spectrum. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Characterisation of the 1997 Vulcanian explosions of Soufrière Hills Volcano,Montserrat, by video analysis

The activity of Convention at Montserrat Soufrière Hills Volcano, Montserrat, during the period 1995–1999 included numerous violent explosions. Two major cycles of Vulcanian explosions occurred in 1997: a first of 13 explosions between 4 and 12 August and a second of 75 between 22 September and 21 October. The explosions were short-lived events lasting a few tens of seconds during which partial fountain collapse generated pyroclastic surges and pyroclastic flows, and buoyant plumes ascended 3–15 km into the atmosphere. Each explosion discharged on average 3×10⁵ m³ (dense-rock equivalent, DRE) of magma, draining the conduit to depths of 1–2 km. The paper focuses on the first few seconds of three explosions of the 75 that occurred in September/October 1997: 6 October 1997 at 17:50, 7 October 1997 at 16:02 and 9 October 1997 at 12:32. Physical parameters such as exit velocities, magmatic water contents and magma pressures at fragmentation are estimated by following and modelling the ascent of individual momentum-dominated finger jets visible on videos during the initial stages of each explosion. The model treats each finger jet as an incompressible flow sustained by a steady flux of gas and particles during the few seconds of ascent, and produces results that compare favourably with those using a multiphase compressible code run using similar eruptive parameters. Each explosion reveals a progressive increase in eruptive intensity with time, jet exit velocities increasing from 40 m s^–1 at the beginning of the explosion up to 140 m s^–1 after a few seconds. Modelling suggests that the first magma to exit was largely degassed, whereas that discharged after a few seconds contained up to 2 wt% water. Magma overpressures up to ~10 MPa are estimated to have existed in the conduit immediately prior to each explosion. Progressive increases in jet exit velocity with time over the first few seconds of each explosion provide direct evidence for strong pre-eruptive gradients in water content and magma pressure in the upper reaches (probably 100–500 m) of the conduit. Fountain collapse occurred during the first 10–20 s of each explosion because the discharging jets had bulk densities up to 100 times that of the atmosphere and were unable to entrain enough air to become buoyant. Such high eruptive densities were due to the presence of partially degassed magma in the conduit.Editorial responsibility: A. Woods     

Indoor Radon Radioactivity at the University of Brunei Darussalam

— Indoor radon radioactivity in the rooms on the ground floor and first floor of the Physics Department, Faculty of Science, Universiti Brunei Darussalam was measured using a system that consists of an air filter pump, ZnS detector, photomultiplier tube and counter. Ground floor rooms' radon radioactivity was found to be about three times higher than that of the first floor. The maximum ground floor indoor radioactivity is only 0.39 Bqm^?3, a value relatively low and safe compared to the mean outdoor radon concentration of 1.41 Bqm^?3 measured (HU and TAN, 2000). The main source of radon emanation originates from the ground soil rather than the building materials.     

Source mechanism of the 22/11/1995 Gulf of Aqaba earthquake and its aftershock sequence

The Gulf of Aqaba earthquake occurred on 22/11/1995 with M _W = 7.2 and was the largest event to occur along the Dead Sea Transform in at least a century. We determined the source mechanism of the event based on the inversion of the P and SH waveforms observed by teleseismic stations. Our solution consists of 2 similar subevents, where the first exhibits a left-lateral strike-slip motion with a normal component and the second subevent shows an almost pure left-lateral strike-slip motion along the gulf major fault system. The total seismic moment is 7.7 × 10²⁶ dyne-cm, with source duration of 15 seconds. The aftershock sequence was recorded by the regional broad band stations BGIO and JER, of the Israel Seismic Network, and KEG, of the MEDNET regional network for earthquakeswith M _W > 4.0. Aftershock epicenter seems to cluster mainly in the Eilat basin and the Aragonese basin. The dominant mechanism indicates left-lateral strike slip motion along the gulf trend, similar to the mainshock.     

Optimum pilot sweep

The successful application of high-resolution seismic methods requires evaluating each element in the seismic system and ensuring that each part of the system contributes optimally to the success of the method. Unfortunately, unlike data processing, seismic signal generation is not carefully optimized. The purpose of our study was to optimize the source signal in order to better coordinate field operations with subsequent data processing to achieve their common objective. We developed an iterative method for a rational frequency distribution of the energy of a seismic source. The method allows the optimum amplitude spectrum of a source signal to be calculated, thus providing the best data quality at the end of the processing. We assume that the source signal is affected by a total transfer function, by the reflectivity function of a target interval, and by ambient noise, whose characteristics, if not known, can be estimated or measured in practice. The transfer function includes data processing other than the correlation stage and the final trace-optimizing filter. The variance of a reflectivity estimate is considered to be a measure of the data quality and improvement of the characteristic corresponds to a decrease in the variance. For this reason, a constrained Wiener deconvolution filter is used as the final trace-optimizing filter. It not only minimizes the variance of a reflectivity estimate but also ensures a specific signal-to-noise ratio. The method is made feasible by following the Vibroseis technique, primarily because of the versatility of the technique in controlling the signal spectrum. With the optimum amplitude spectrum obtained, the corresponding optimum pilot sweep can be readily calculated. Examples using synthetic data are presented to illustrate the method.