Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

The possible contribution of circumferential fault intrusion to caldera resurgence

In the geological record, the intrusion of substantial amounts of magma into circumferential faults and ring fractures is commonly observed. Finite element modelling is used here to investigate the strain field that may be expected from such intrusive events. Two simple vertical scenarios are explored, one for a caldera with a central block of thickness to diameter ratio of ~1:1 (similar to Rabaul) and one with a ratio much less than 1:1 (similar to the Valles type). Surface deformation in both cases is similar with central uplift, the development of a moat (or trough) like feature just outside of the intersection of the azimuth of the intruded ring fault and the free surface, and broader scale tumescence at a scale several times larger than the calderas radius. The response of the block and sub-caldera magma chamber for the two scenarios, however, is different. The blocks are in effect squeezed; the high aspect ratio one deforms upwards at the surface and downwards at its base, whereas the low aspect ration one experiences up arching (or bending) of the central part of the caldera block. Central uplift still occurs when only a short arc of a ring fracture system or a circumferential fault is intruded. In both models, tumescence in the centre of the caldera from single ring dyke intrusion can only account for decimetres to metres of surface uplift. Repeated intrusions over tens to hundreds of thousands of years, however, may cause incremental up doming of the caldera block leading to larger scale resurgent features. The amount of uplift possible due to squeezing of a high aspect ratio block is limited. It is proposed, however, that where bending of plate-like blocks occur above a decompressible and/or malleable magma body, ring fault intrusion may be a significant contributor to resurgence. In the simple conceptual models shown here, the amount of ring dyke-induced central uplift will be >40–50% of the width of the ring complex. In the geological record the accumulation of intrusions into some ring fractures has led to annular or arcuate plutons of hundreds of meters to several kilometres in thickness. At certain calderas such intrusions may be a control on the marked concentration of uplift within the restricted area defined by the caldera faults. The complex nature of the horizontal displacements associated with the intrusion of ring and arcuate dykes is also explored. Intrusion into ring fracture zones will tend to take place into those sectors of the annular zone which are perpendicular to the least compressive stress vector. This may be a factor in the observed difference for caldera evolution in extensional and compressional areas. The unrest at several modern calderas is tentatively related to circumferential fault intrusion.Editorial responsibility: J. Stix     

Ground deformation and gravity changes accompanying the 1982 Pozzuoli uplift

During the summer of 1982 a continuous uplift began at Phlegraean Fields, an active volcanic area in southern Italy (Fig. 1), that persisted up to September 1984. The uplift, which reached a maximum value of about 160 cm in the central part of the Phlegraean caldera, was characterized by variable velocity and occurred within an area that extended about 7 km outward from the town of Pozzuoli (Fig. 1); the surface deformation performs a bell-shape pattern.The uplift was accompanied by horizontal displacements and gravity changes that closely correlate with the described elevation changes.The horizontal displacements displayed an anomalous pattern within a narrow belt about 1 km from the center of the uplift and approximately coincident with the area of maximum seismic activity. The change in gravity is attributed primarily to a free-air effect to which a small Bouguer effect must be added.Several models have been invoked in order to explain the observed phenomenon. The one which gives the best fit to the observed data is an increasing pressure source of radial simmetry, at a constant depth of about 3 km beneath the town of Pozzuoli, and having a diameter of several hundred meters. Migration of magma at depth is believed to be responsible for the observed activity in the Phlegraean Fields caldera.     

Simulation of the Seismic Response of Sedimentary Basins with Vertical Constant-Gradient Velocity for Incident <Emphasis Type="Italic">SH</Emphasis> Waves

The simulation of the seismic response of heterogeneous sedimentary basins under incident plane waves is computed using the Indirect Boundary Element Method (IBEM). To deal with these kinds of basins we used approximate analytical expressions for the two-dimensional Greens functions of a medium with constant-gradient wave propagation velocity. On the other hand, for the homogeneous half space underlying the sedimentary basin, the full space Greens functions were used. The response of semi-circular heterogeneous basins under incident SH waves is explored by means of the displacements in the frequency-space diagrams and synthetic seismograms. Moreover, we compared these results with those obtained for other homogeneous semi-circular models. The principal differences among them are pointed out. This simulation provided interesting results that displayed a complex amplification pattern in a rich spectrum of frequencies and locations. The maximum amplitudes levels were found around the edges of the heterogeneous sedimentary basins. In time domain some features characterize the seismic response of the basin which include enhanced generation and trapping of surface waves inside the sediments, and the reduced emission of seismic energy to the hard rock. In the heterogeneous models the lateral reflections of surface waves greatly influence the total displacements at the free surface in comparison with the homogeneous models where the displacements have a shared influence among both vertical and lateral reflections.     

Far-field Gravity and Tilt Signals by Large Earthquakes: Real or Instrumental Effects?

A wide set of dynamics phenomena (i.e., Geodynamics, Post Glacial Rebound, seismicity and volcanic activity) can produce time gravity changes, which spectrum varies from short to long (more than 1 year) periods. The amplitude of the gravity variations is generally in the order of consequently their detection requires instruments with high sensitivity and stability: then, high quality experimental data. Spring and superconducting gravimeters are intensively used with this target and they are frequently jointed with tiltmeters recording stations in order to measure the elasto-gravitational perturbation of the Earth. The far-field effects produced by large earthquakes on records collected by spring gravimeters and tiltmeters are investigated here. Gravity and tilt records were analyzed on time windows spanning the occurrence of large worldwide earthquakes; the gravity records have been collected on two stations approximately 600 km distant. The background noise level at the stations was characterized, in each season, in order to detect a possible seasonal dependence and the presence of spectral components which could hide or mask other geophysical signals, such as, for instance, the highest mode of the Seismic Free Oscillation (SFO) of the Earth. Some spectral components (6.5; 8; 9; 14, 20, 51) have been detected in gravity and tilt records on the occasion of large earthquakes and the effect of the SFO has been hypothesized. A quite different spectral content of the EW and NS tiltmeter components has been detected and interpreted as a consequence of the radiation pattern of the disturbances due to the earthquakes. Through the analysis of the instrumental sensitivity, instrumental effects have been detected for gravity meters at very low frequency.F.S.E. (Fondo Sociale Europeo -European Community -)     

A Local Least-Squares Modification of Stokes’ Formula

The combination of Stokes formula and an Earth Gravity Model (EGM) for geoid determination has become a standard procedure. However, the way of modifying Stokes formula vary from author to author, and numerous methods of modification exist. Most methods are deterministic, with the primary goal of reducing the truncation bias committed by limiting the area of Stokes integration around the computation point, but there are also some stochastic methods with the explicit goal to reduce the global mean square error of the geoid height estimator stemming from the truncation bias as well as the random errors of the EGM and the gravity data. The latter estimators are thus, at least from a theoretical point of view, optimal in a global mean sense, but in a local sense they may be far from optimality.Here we take advantage of the error variance-covariance matrices of the EGM and the terrestrial gravity data to derive the modification parameters of Stokes kernel in a local least-squares sense. The solution is given for the unbiased type of modification of Stokes formula of Sjöberg (1991).     

A note on the general concept of wave breaking for Rossby and gravity waves

A recently proposed general definition of wave breaking is further discussed, in order to deal with some points on which misunderstanding appears to have arisen. As with surface and internal gravity waves, the classification of Rossby waves into breaking and not breaking is a generic classification based on dynamical considerations, and not a statement about any unique signature or automatically recognizable shape. Nor is it a statement about passive tracers uncorrelated with potential vorticity on isentropic surfaces. A strong motivation for the definition is that proofs of the nonacceleration theorem of wave, mean-flow interaction theory rely, explicitly or implicitly, on a hypothesis that the waves do not break in the sense envisaged.The general definition refers to the qualitative behaviour of a certain set of material contours, namely those, and only those, which would undulate reversibly, with small slopes, under the influence of the waves' restoring mechanism, in those circumstances for which linearized, nondissipative wave theory is a self-consistent approximation to nonlinear reality. The waves' restoring mechanism depends upon the basic-state vertical potential density gradient in the case of gravity waves, and upon the basic-state isentropic gradient of potential vorticity in the case of Rossby waves. In the usual linearized theory of planetary scale Rossby waves on a zonal shear flow, the relevant material contours lie along latitude circles when undisturbed.     

Standing Alfvén waves with m ≫ 1 in an axisymmetric magnetosphere excited by a non-stationary source

As a continuation of our earlier paper, we consider here the case of the excitation of standing Alfvén waves by a source of the type of sudden impulse. It is shown that, following excitation by such a source, a given magnetic shell will exhibit oscillations with a variable frequency which increases from the shells poloidal to toroidal frequency. Simultaneously, the oscillations will also switch over from poloidally (radially) to toroidally (azimuthally) polarized. With a reasonably large attenuation, only the start of this process, the stage of poloidal oscillations, will be observed in the ionosphere.     

Participatory methods of incorporating scientific with traditional knowledge for volcanic hazard management on Ambae Island,Vanuatu

Ambae Island is the largest of Vanuatus active volcanoes. It is also one of the nations potentially most dangerous, with 60 million m³ of lake-water perched at over 1340 m in the summit caldera and over the active vent. In 1995, small phreatic explosions, earthquake swarms and heightened gas release led to calls for evacuation preparation and community volcanic hazard awareness programs for the ~9500 inhabitants. Differences in perspective or world-view between the island dwellers adhering to traditional beliefs (Kastom) and external scientists and emergency managers led to a climate of distrust following this crisis. In an attempt to address these issues, rebuild dialogue and respect between communities, outside scientists and administrators, and move forward in volcanic hazard education and planning for Ambae, we adapted and applied Participatory Rural Appraisal (PRA) approaches. Initial gender-segregated PRA exercises from two representative communities provided a mechanism for cataloguing local traditional viewpoints and hazard perceptions. Ultimately, by combining elements of these viewpoints and perceptions with science-based management structures, we derived volcanic hazard management guidelines, supported by an alert system and map that were more readily accepted by the test communities than the earlier top-down plans imposed by outside governmental and scientific agencies. The strength of PRA approaches is that they permit scientists to understand important local perspective issues, including visualisations of volcanic hazards, weaknesses in internal and external communication systems, and gender and hierarchy conflicts, all of which can hinder community emergency management. The approach we describe has much to offer both developing and industrialised communities that wish to improve their awareness programs and mitigative planning. This approach should also enhance communication and understanding between volcanologists and the communities they serve.     

Plinian fall deposits from Vulsini Volcanic District (Central Italy)

Stratigraphic investigations of the Vulsini Volcanic District indicate the existence of five prominent plinian fall deposits, deposited over a period of about 0.34 Ma. The oldest deposit (Basal pumices) crops out mainly in the peripheral areas and is one of the largest plinian events in the Vulsini District, with a volume of about 9 km³ and a calculated column height exceeding 30 km. Subcircular patterns of isopachs and isopleths around the Bolsena lake indicate emplacement in wind-free conditions and suggest a source vent location in the northernmost sector of the present Bolsena lake, where it is possible that an older central volcanic structure existed. The four younger plinian fall deposits are related to Bolsena activity. The oldest of these, named the Ponticello and Orvieto-Bagnoregio pumices, are mainly distributed on the north-eastern sector of the Vulsini District. Their volumes are an order of magnitude lower than that of the basal unit, and the estimated column heights do not exceed 20 km. The younger Ospedaletto pumice deposit has a NE-SW dispersal axis, whereas the youngest Casetta pumice deposit is found only in the north-eastern sector of the district. Their estimated volumes are respectively 1.2 and 0.1 km³, whereas the inferred vent areas appear to be located slightly to the east of the Ponticello and Orvieto-Bagnoregio pumice deposit source areas. The chronology of the Ospedaletto and Casetta pumices indicates that the final plinian activity from the Bolsena complex is contemporaneous with the Latera activity. The decreasing volumes from the oldest to the youngest units, together with the progressive shifting northeastwards of the source vents, may be related to the volcano-tectonic subsidence of the Bolsena area. The source vents for the post-Basal pumices events are mainly clustered just north of Bolsena village where the greatest displacement during subsidence occurred. Accordingly, source vent shifting appears to be related to the progressive opening of normal faults from the inner to the outer margin of the Bolsena depression.     

On the conditions of sheet injections and eruptions in stratovolcanoes

Nearly all eruptions in stratovolcanoes (composite volcanoes, central volcanoes) are supplied with magma through fractures. Consequently, a primary physical condition for an eruption to occur in a stratovolcano is that a magma-driven fracture is able to propagate to the surface. Magma-filled fractures, frozen or fluid, are referred to as sheet intrusions. More specifically, they are named dykes when subvertical, and inclined (or cone) sheets when inclined. Field observations indicate that most sheet intrusions do not reach the surface to feed eruptions but rather become arrested at various crustal depths. For this reason periods of volcanic unrest with sheet injections are much more common than volcanic eruptions. Whether a sheet intrusion becomes arrested or, alternatively, propagates to the surface depends primarily on the stress field in the stratovolcano. A stratovolcano normally consists of layers of contrasting mechanical properties, such as soft (low Youngs modulus) pyroclastic units and stiff (high Youngs modulus) lava flows. We present numerical models indicating that volcanoes composed of such layers commonly develop stress fields encouraging sheet and dyke arrest. The models indicate that a necessary condition for a sheet intrusion to reach the surface and feed a volcanic eruption is that the stress field along the sheet pathway becomes homogenised. We propose that much of the activity in a stratovolcano during a volcanic cycle encourages stress-field homogenisation. Field studies show that the sheet intrusions in individual stratovolcanoes have various dips: some are vertical dykes, others inclined sheets, and still others horizontal sills. Analytical models indicate that the dip of a sheet reaching the surface can have great effects on the magma transport during an eruption. This effect is normally greater for a flat volcano such as a collapse caldera than for a stratovolcano that forms a topographic high. We conclude that the shallower the dip of a sheet intrusion, the less will be its volumetric magma transport to the surface of a stratovolcano.Editorial responsibility: D Dingwell     

Modelling gravity variations consistent with ground deformation in the Campi Flegrei caldera (Italy)

In recent years (1970–72 and 1982–84) two inflation episodes took place in the Campi Flegrei caldera (Italy), characterized by significant ground uplift and gravity variations. An elastic half-space model with vertical density stratification is employed to compute the displacement field and the gravity variations produced by the deformation of buried layers, following the inflation of a spherically symmetric deformation source. Contributions to gravity variations are produced by dilation/contraction of the medium, by the displacements of density interfaces (the free surface and subsurface layers) and of source boundaries and, possibly, by new mass input from remote distances into the source volume. Three cases were examined in detail: In case I, the magma chamber is identified as the deformation source and volume and pressure increase in the magma chamber is due to input of new magma from remote distances; in case II deformation is due to magma differentiation within the magma chamber (deformation source with constant mass); in case III the geothermal system is identified as the deformation source and a pressure increase, possibly driven by the exsolution of high temperature and high pressure volatiles in the magma chamber, is assumed to play a dominant role. From the comparison between measured and computed gravity residuals (free-air-corrected gravity variations) we can assess that, in case I, an inflation source with constant density would predict gravity residuals compatible with observations, whereas an expansion at constant mass (case II) would predict gravity residuals much lower than observed. The resolving power of gravity data however prevents accurate assessment of the density of the emplaced material. In case III, the pervasive density increase of the geothermal fluids induced by pressure increase is assumed to be the main source of gravity variations. The average porosity value required for this model to match both the ground deformation and the gravity residuals is found to be ˜10%, a value which is compatible with measured porosity values at Campi Flegrei in deep wells. The subsidence phases following both inflation episodes and the gravity residuals during subsidence lead us to consider case III as more plausible, even if a suitable combination of cases I and III cannot be discarded.     

Seismic energy partitioning and scattering in laterally heterogeneous ocean crust

We present finite difference forward models of elastic wave propagation through laterally heterogeneous upper oceanic crust. The finite difference formulation is a 2-D solution to the elastic wave equation for heterogeneous media and implicitly calculatesP andSV propagation, compressional to shear conversion, interference effects and interface phenomena. Random velocity perturbations with Gaussian and self-similar autocorrelation functions and different correlation lengths (a) are presented which show different characteristics of secondary scattering. Heterogeneities scatter primary energy into secondary body waves and secondary Stoneley waves along the water-solid interface. The presence of a water-solid interface in the model allows for the existence of secondary Stoneley waves which account for much of the seafloor noise seen in the synthetic seismograms for the laterally heterogeneous models.Random incoherent secondary scattering generally increases aska (wavenumber,k, and correlation length,a) approaches one. Deterministic secondary scattering from larger heterogeneities is the dominant effect in the models aska increases above one. Secondary scattering also shows up as incoherence in the primary traces of the seisograms when compared to the laterally homogeneous case. Cross-correlation analysis of the initialP-diving wave arrival shows that, in general, the correlation between traces decreases aska approaches one. Also, because many different wave types exist for these marine models, the correlation between traces is range dependent, even for the laterally homogeneous case.     

Ages of calderas,large explosive craters and active volcanoes in the Kuril-Kamchatka region,Russia

The ages of most of calderas, large explosive craters and active volcanoes in the Kuril-Kamchatka region have been determined by extensive geological, geomorphological, tephrochronological and isotopic geochronological studies, including more than 600 ¹⁴C dates. Eight Krakatoa-type and three Hawaiian-type calderas and no less than three large explosive craters formed here during the Holocene. Most of the Late Pleistocene Krakatoa-type calderas were established around 30 000–40 000 years ago. The active volcanoes are geologically very young, with maximum ages of about 40 000–50 000 years. The overwhelming majority of recently active volcanic cones originated at the very end of the Late Pleistocene or in the Holocene. These studies show that all Holocene stratovolcanoes in Kamchatka were emplaced in the Holocene only in the Eastern volcanic belt. Periods of synchronous, intensified Holocene volcanic activity occurred within the time intervals of 7500–7800 and 1300–1800 ¹⁴C years BP.     

3-D Modelling of Campi Flegrei Ground Deformations: Role of Caldera Boundary Discontinuities

Campi Flegrei is a caldera complex located west of Naples, Italy. The last eruption occurred in 1538, although the volcano has produced unrest episodes since then, involving rapid and large ground movements (up to 2 m vertical in two years), accompanied by intense seismic activity. Surface ground displacements detected by various techniques (mainly InSAR and levelling) for the 1970 to 1996 period can be modelled by a shallow point source in an elastic half-space, however the source depth is not compatible with seismic and drill hole observations, which suggest a magma chamber just below 4 km depth. This apparent paradox has been explained by the presence of boundary fractures marking the caldera collapse. We present here the first full 3-D modelling for the unrest of 1982–1985 including the effect of caldera bordering fractures and the topography. To model the presence of topography and of the complex caldera rim discontinuities, we used a mixed boundary elements method. The a priori caldera geometry is determined initially from gravimetric modelling results and refined by inversion. The presence of the caldera discontinuities allows a fit to the 1982–1985 levelling data as good as, or better than, in the continuous half-space case, with quite a different source depth which fits the actual magma chamber position as seen from seismic waves. These results show the importance of volcanic structures, and mainly of caldera collapses, in ground deformation episodes.     

Note on fault-slip motion inferred from a study of micro-cataclastic particles from an underground shear rupture

Mining of a highly-stressed remnant in a deep South African gold mine was accompanied by considerable seismic activity and some significant rockbursts. The larger seismic events were registered some 60 km away at a WSSN station and several shear ruptures corresponding to these events were encountered during mining operations.A careful study based on detailed exploration of two of these ruptures proved them to be the source of two of the larger rockbursts.Certain striking features revealed by a scanning electron microscopic study of some of the fresh cataclastic rock-flour forming part of the comminuted filling of these ruptures provide strong evidence of violent shock rebound phenomena in the faulting process. This interpretation could provide useful insight into earthquake source mechanisms and also has practical significance in the understanding of mine rockbursts.     

An experimental study of the seismic ringing of the moon

Summary The seismic records obtained from the seismic station ALSEP (Apollo Lunar Surface Experiments Package) on the moon, due to the impact of the LM (Lunar Module) of the Apollo XII and the S-IVB rocket state of the Apollo XIII moon missions show the seismic signals decaying in amplitude slowly with time. The observed ringing phenomenon on these records along with the fact that the disturbances lasted for a long period were unexpected. Laboratory experiments on homogeneous spheres have produced this ringing. The physical explanation of this may be complicated, and it is quite possible that this observed ringing is due to the combined effect of the filtering, reflections and scattering of energy due to the curvature of the homogeneous (in the seismic sense) sphere. The explanation of the ringing of the moon may very well lie in such a study.     

A comparison of short-wave and long-wave velocities in layered media

Summary Ray theory based on Snell's law is strictly valid only when the variation of physical constants is small within a wavelength. In explosion seismology, Snell's law may apply to the arrival times of the first breaks, but not to those of the seismic reflections. For a stack of homogeneous, isotropic plane layers, the average velocity at its short-wavelength limit is greater than that at the long-wavelength limit for any given angle of incidence. Considering a few example of layered media, calculations are made for both short-wave and long-wave average velocities as a function of the angle of incidence. The former velocity shows a steady increase with the angle of incidence whereas the latter velocity may sometimes pass through a minimum before attaining its maximum value for horizontal propagation.