The 1963–1964 eruption of Agung volcano (Bali, Indonesia)

The February 1963 to January 1964 eruption of Gunung Agung, Indonesia’s largest and most devastating eruption of the twentieth century, was a multi-phase explosive and effusive event that produced both basaltic andesite tephra and andesite lava. A rather unusual eruption sequence with an early lava flow followed by two explosive phases, and the presence of two related but distinctly different magma types, is best explained by successive magma injections and mixing in the conduit or high level magma chamber. The 7.5-km-long blocky-surfaced andesite lava flow of ～0.1?km³ volume was emplaced in the first 26?days of activity beginning on 19 February. On 17 March 1963, a major moderate intensity (～4?×?10⁷?kg?s^?1) explosive phase occurred with an ～3.5-h-long climax. This phase produced an eruption column estimated to have reached heights of 19 to 26?km above sea level and deposited a scoria lapilli to fine ash fall unit up to ～0.2?km³ (dense rock equivalent—DRE) in volume, with Plinian dispersal characteristics, and small but devastating scoria-and-ash flow deposits. On 16 May, a second intense 4-h-long explosive phase (2.3?×?10⁷?kg?s^?1) occurred that produced an ～20-km-high eruption column and deposited up to ～0.1?km³ (DRE) volume of similar ash fall and pyroclastic flow deposits, the latter of which were more widespread than in the March phase. The two magma types, porphyritic basaltic andesite and andesite, are found as distinct juvenile scoria populations. This indicates magma mixing prior to the onset of the 1963 eruption, and successive injections of the more mafic magma may have modulated the pulsatory style of the eruption sequence. Even though a total of only ～0.4?km³ (DRE volume) of lava, scoria and ash fall, and scoria-and-ash pyroclastic flow deposits were produced by the 1963 eruption, there was considerable local damage caused mainly by a combination of pyroclastic flows and lahars that formed from the flow deposits in the saturated drainages around Agung. Minor explosive activity and lahar generation by rainfall persisted into early 1964. The climactic events of 17 March and 16 May 1963 managed to inject ash and sulfur-rich gases into the tropical stratosphere.     

Studies of mantle structure of U.S.S.R. territory on long-range seismic profiles

Long-range seismic sounding carried out during the last few years on the territory of the U.S.S.R. has shown a basic inhomogeneity of the uppermost mantle, as well as evidence of regularities in the distribution of its seismic parameters. The following data were used: times and apparent velocities of P- and S-waves for investigation of mantle velocities, converted waves for seismic discontinuity model studies and wave attenuation for Q-factor estimation. Strong regularities were distinguished in the distribution of average seismic velocities for the uppermost mantle, in their dependence on the age and type of geostructure and on their position relative to the central part of the continent. Old platforms and the inner part of the continent are marked by velocities under the Mohorovi?i? discontinuity of more than 8.2–8.3 km s^?1, young platforms and outer parts of the continent by 8.0–8.2 km s^?1, and orogenic and rift zones by 7.8–8.0 km s^?1. The difference becomes more pronounced at a depth of about 100–200 km: for the old platform mantle velocities of 8.5–8.6 km s^?1 are typical; beneath the orogenic and rift areas, inversion zones with velocities less than 7.8 km s^?1 are observed.The converted waves show fine inhomogeneities of the crust and uppermost mantle, the presence of many discontinuities with positive and negative changes of velocity, and anisotropy of seismic waves in some of the layers. Wave attenuation allowed the determination of the Q-factor in the mantle. It varied from one region to another but a close relation between Q and P-wave velocity is the main cause of its variation.     

Insights into eruption dynamics from textural analysis: the case of the May, 2008, Chaitén eruption

The May, 2008, Chaitén (southern Chile) eruption was characterized by several explosive events, each associated with plumes which reached up to about 19?km above sea level on May 6. A study of the textural and physical features of the juvenile clasts erupted during the climactic phase of the 2008 eruption of Chaitén is presented. Pumice clasts show unimodal density distribution (main mode at 600?kg/m³), average vesicularity of about 69?%, a glassy groundmass with no microcrystals, and vesicles with dimension between ～1?μm and ～2?mm. They also show a unimodal vesicle size distribution with most frequent vesicle size in the range 0.05–0.08?mm and an estimated vesicle number density of 1.3?±?0.5?×?10⁵?mm^?3 related to a rapid nucleation event produced during the late phases of magma rise. This is confirmed by the absence of microcrystals that could otherwise have delayed vesicle formation and allowed the magma to maintain a low viscosity and a supersaturation in volatiles. Vesiculation and fragmentation were triggered by a sudden decompression of the melt associated with the opening of the volcanic conduit (～10?MPa?s^?1).     

Volumetric characteristics of lava flows from interferometric radar and multispectral satellite data: the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos

We have used a suite of remotely sensed data, numerical lava flow modeling, and field observations to determine quantitative characteristics of the 1995 Fernandina and 1998 Cerro Azul eruptions in the western Galápagos Islands. Flank lava flow areas, volumes, instantaneous effusion rates, and average effusion rates were all determined for these two eruptions, for which only limited syn-eruptive field observations are available. Using data from SPOT, TOPSAR, ERS-1, and ERS-2, we determined that the 1995 Fernandina flow covers a subaerial area of 6.5×10⁶ m² and has a subaerial dense rock equivalent (DRE) volume of 42×10⁶ m³. Field observations, ATSR satellite data, and the FLOWGO numerical model allow us to determine that the effusion rate declined exponentially from a high of ~60–200 m³ s^-1 during the first few hours to <5 m³ s^-1 prior to ceasing after 73 days, with a mean effusion rate of 4–16 m³ s^-1. Integrating the ATSR-derived, exponentially declining effusion rate over the eruption duration produces a total (subaerial + submarine) DRE volume of between 27 and 100×10⁶ m³, the range in values being due to differing assumptions about heat loss characteristics; only values in the higher part of this range are consistent with the independently derived subaerial volume. Using SPOT, TOPSAR, ERS-1, and ERS-2 data, we determine that the 1998 Cerro Azul flow is 16 km long, covers 16 km², and has a DRE volume of 54×10⁶ m³. FLOWGO produces at-vent velocity and effusion rate values of 11 m s^-1 and ~600 m³ s^-1, respectively. The velocity value agrees well with the 12 m s^-1 estimated in the field. The mean effusion rate (total DRE volume/duration) was 7–47 m³ s^-1. Dike dimensions, fissure lengths, and pressure gradients along the conduit based on magma chamber depth estimates of 3–5 km produce mean effusion rates for the two eruptions that range over nearly four orders of magnitude, the range being due to uncertainty in the magma viscosity, dike dimensions, and pressure gradient between magma chamber and vent. Although somewhat consistent with mean effusion rates from other techniques, their wide range makes them less useful. The exponentially declining effusion rates during both eruptions are consistent with release of elastic strain being the driving mechanism of the eruptions. Our results provide independent input parameters for previously published theoretical relationships between magma chamber pressurization and eruption rates that constrain chamber volumes and increases in volume prior to eruption, as well as time constants of exponential decay during the eruption. The results and theoretical relationships combine to indicate that at both volcanoes probably 25–30% of the volumetric increase in the magma chamber erupted as lava onto the surface. In both eruptions the lava flow volumes are less than 1% of the magma chamber volume.     

The thermal boundary-layer interpretation of D″ and its role as a plume source

The “anomalous” layer in the lowermost mantle, identified as D″ in the notation of K.E. Bullen, appears in the PREM Earth model as a 150 km-thick zone in which the gradient of incompressibility with pressure, $\text{d}\text{K}\text{d}\text{P}$, is almost 1.6, instead of 3.2 as in the overlying mantle. Since PREM shows no accompanying change in density or density gradient, we identify D″ as a thermal boundary layer and not as a chemically distinct zone. The anomaly in $\text{d}\text{K}\text{d}\text{P}$ is related to the temperature gradient by the temperature dependence of K_s, for which we present a thermodynamic identity in terms of accessible quantities. This gives the numerical result (?K_s/?T)_P=?1.6×10⁷ Pa K^?1 for D″ material. The corresponding temperature increment over the D″ range is 840 K. Such a layer cannot be a static feature, but must be maintained by a downward motion of the lower mantle toward the core-mantle boundary with a strong horizontal flow near the base of D″. Assuming a core heat flux of 1.6 × 10¹² W, the downward speed is 0.07 mm y^?1 and the temperature profile in D″, scaled to match PREM data, is approximately exponential with a scale height of 73 km. The inferred thermal conductivity is 1.2 W m^?1 K^?1. Using these values we develop a new analytical model of D″ which is dynamically and thermally consistent. In this model, the lower-mantle material is heated and softened as it moves down into D″ where the strong temperature dependence of viscosity concentrates the horizontal flow in a layer ～ 12 km thick and similarly ensures its removal via narrow plumes.     

It goes both ways: measurements of simultaneous evapotranspiration and fog droplet deposition at a montane cloud forest

Fluxes of latent heat, sensible heat, and water vapor, including turbulent deposition of fog droplets, were measured for two months in autumn 2005 within a subtropical montane cypress forest in Taiwan. The goal of the study was to determine whether significant evapotranspiration can occur during foggy conditions. Water vapor fluxes, Q_W, as determined with the Bowen Ratio method, were compared to those simultaneously measured with the eddy covariance method. The median Bowen Ratio was 1.06, and the median Q_W flux was 5 · 2 × 10^?5 kg m^?2 s^?1. The vertical gradients of temperature and specific humidity over the forest, ΔT and Δq, peaked around noon during days without fog, and were reduced during foggy conditions. For 66% of the data points, ΔT and Δq were negative, corresponding to positive (upward) fluxes of sensible heat Q_H and latent heat Q_E. A Monte Carlo simulation proved that statistically significant evapotranspiration rates, i.e., upward water vapor fluxes, occurred during fog. At the same time, deposition fluxes of fog droplets occurred. Our results show that even during fog events, significant evapotranspiration may occur. Copyright © 2008 John Wiley & Sons, Ltd.     

Submarine pyroclastic deposits formed during the 20th May 2006 dome collapse of the Soufrière Hills Volcano, Montserrat

The 20th May 2006 lava dome collapse of the Soufrière Hills Volcano, Montserrat, had a total non-dense rock equivalent (non-DRE) collapse volume of approximately 115?×?10⁶?m³. The majority of this volume was deposited into the ocean. The collapse was rapid, 85% of the mobilized volume being removed in just 35?min, giving peak pyroclastic flow flux of 66?×?10³?m³?s^?1. Channel and levee facies on the submarine flanks of the volcano and formation of a thick, steep-sided ridge, suggest that the largest and most dense blocks were transported proximally as a high concentration granular flow. Of the submerged volume, 30% was deposited from the base of this granular flow, forming a linear, high-relief ridge that extends 7?km from shore. The remaining 70% of the submerged volume comprises the finer grain sizes, which were transported at least 40?km by turbidity currents on gradients of <2°. At several localities, the May 2006 distal turbidity currents ran up 200?m of topography and eroded up to 20?cm of underlying substrate. Multiple turbidites are preserved, representing current reflection from the graben margins and deflection around topography. The high energy of the May 2006 collapse resulted in longer submarine run out than the larger (210?×?10⁶?m³) Soufrière Hills dome collapse in July 2003.     

The mechanical and electrical properties of earth's asthenosphere

An adequate theory of continental drift can be based on heat transfer theory, but it does demand the acceptance of a large downward revision of traditional estimates of average upper mantle temperatures and a consistent understanding of lithosphere and asthenosphere in terms of a difference in rheological behaviour under prolonged non-hydrostatic stress. The recognition that an extremely viscous average state of the upper mantle is self regulating both requires and permits an explanation of magma generation at a strictly limited rate (when averaged for the whole planet over a few years) in terms of unsteady and local deformational heating.The activity of water as a reducer of silicate creep resistance is used to develop the hypothesis that water produced by an amphibole dehydration has been effectively trapped in the Earth and is the underlying cause of a low seismic Q ～ 50 and an electrical conductivity 10^?2 ?10^?1 ohm^?1 m^?1, at depths of ～ 100 km. At the predicted low horizontally-avera temperatures, the conductivity contrast of rock and aqueous solutions is very large, and mantle electrical conductivity studies now look best-suited to explore this trapping process, and the distinctly recognisable possibility that the uptake of ocean water in the subduction process exceeds the rate of loss that can be explained purely through magmatic activity.     

Modelling the dynamics and thermodynamics of volcanic degassing

 The rates of passive degassing from volcanoes are investigated by modelling the convective overturn of dense degassed and less dense gas-rich magmas in a vertical conduit linking a shallow degassing zone with a deep magma chamber. Laboratory experiments are used to constrain our theoretical model of the overturn rate and to elaborate on the model of this process presented by Kazahaya et al. (1994). We also introduce the effects of a CO₂–saturated deep chamber and adiabatic cooling of ascending magma. We find that overturn occurs by concentric flow of the magmas along the conduit, although the details of the flow depend on the magmas' viscosity ratio. Where convective overturn limits the supply of gas-rich magma, then the gas emission rate is proportional to the flow rate of the overturning magmas (proportional to the density difference driving convection, the conduit radius to the fourth power, and inversely proportional to the degassed magma viscosity) and the mass fraction of water that is degassed. Efficient degassing enhances the density difference but increases the magma viscosity, and this dampens convection. Two degassing volcanoes were modelled. At Stromboli, assuming a 2 km deep, 30% crystalline basaltic chamber, containing 0.5 wt.% dissolved water, the ∼700 kg s^–1 magmatic water flux can be modelled with a 4–10 m radius conduit, degassing 20–100% of the available water and all of the 1 to 4 vol.% CO₂ chamber gas. At Mount St. Helens in June 1980, assuming a 7 km deep, 39% crystalline dacitic chamber, containing 4.6 wt.% dissolved water, the ∼500 kg s^–1 magmatic water flux can be modelled with a 22–60 m radius conduit, degassing ∼2–90% of the available water and all of the 0.1 to 3 vol.% CO₂ chamber gas. The range of these results is consistent with previous models and observations. Convection driven by degassing provides a plausible mechanism for transferring volatiles from deep magma chambers to the atmosphere, and it can explain the gas fluxes measured at many persistently active volcanoes. Received: 26 September 1997 / Accepted: 11 July 1998     

Application of age and Sr isotope data to the petrogenesis of the Marangudzi ring complex,Rhodesia

The Marangudzi ring complex, Rhodesia, consists essentially of a gabbro mass intruded by ring dykes of quartz syenite and cone sheets of nepheline syenite. Five intrusive units (gabbro, two quartz syenite and two nepheline syenite units) have been studied using Rb-Sr and K-Ar methods. A total of 24 whole rock samples define a Rb-Sr isochron which gives an age of 186 ± 3m.y. and an initial (⁸⁷Sr/⁸⁶Sr)₀ ratio of 0.70769 ± 0.00006 (±2sigma; based on λ = 1.42 × 10^?11yr^?1). K-Ar and Rb-Sr analyses on biotite and hastingsite separates are consistent with this age assignment. Whole rock Rb-Sr isochrons for the different units treated individually agree with the above age and initial Sr ratio within analytical uncertainties. This is believed to indicate that the different rock types are comagmatic forming by fractional crystallization of a parental, mantle-derived, K₂O-rich basaltic magma, having an initial Sr ratio of 0.7077, without appreciable assimilation of the Precambrian country rock. The entire differentiation, emplacement and crystallization processes took place over a rather short time span.     

Stratigraphy, chemistry, and eruptive dynamics of the 12.4 ka plinian eruption of Apoyeque volcano, Managua, Nicaragua

Apoyeque volcano, located 9 km northwest of Managua city, erupted explosively at 12.4 ka. The Plinian eruption deposited a widespread pumice fall deposit known as the Upper Apoyeque Tephra (UAq). The UAq is massive, reversely graded, and consists of white juvenile pumice (~78 vol.%), a variety of cognate lithics and accidental altered lithics. The whole-rock pumice composition is rhyodacitic (SiO₂?=?66.9–68.5 wt.%) with a mineral paragenesis of plagioclase, orthopyroxene, clinopyroxene, amphibole, titanomagnetite, and ilmenite in a rhyolitic glass groundmass (SiO₂?=?74.4?±?0.6 wt.%). The deposit’s dispersal axis is to the south, with the deposit covering a minimum area of 877 km² within the 50 cm isopach and has a total volume of 3 km³ (dense rock equivalent, 1.15 km³). The eruption column reached a maximum height of ca.28 km. The eruption ejected a total mass of 3?×?10¹² kg at an average rate of 2?×?10⁸ kg/s, and based on available models, we infer duration of almost 4 h. Petrographic and geochemical characteristics suggest that the eruption was triggered by magma mixing.     

Determination of the lifetime of molecules O 2 + and estimation of the flux of penetrating electrons generating a cloud of excess concentration of O 2 + in the ionosphere at altitudes above 220 km

The occurrence of anomalous (nonthermal) profiles of green emission of oxygen atoms detected with a Fabry-Perot spectrometer in auroras with the effect of a rapid decrease in the intensity of the wings of their dissociative component has been investigated. Based on an analysis of these measured profiles, it has been found that the characteristic time of recombination of a molecular oxygen ion at altitudes of 200–400 km is about 5–7 s. It appears that these molecular ions occur in a horizontally limited region of the auroral ionosphere as a result of ionization by a space localized flux of soft electrons with energies of 0.2–0.4 keV penetrating up to altitudes of 200 km. The estimation of the electron flux produces a value of 10¹⁰–10¹³ electrons cm^?2 s^?1. They generate the excess concentration n(O ₂ ⁺ ) ~ 5.6 × 10⁵ cm^?3.     

Geoid anomalies across fracture zones and the thickness of the lithosphere

Recent advances in the measurement and interpretation of geoid height anomalies provide a new way to estimate the thickness of the oceanic lithosphere as a function of crustal age. GEOS-III satellite altimetry measurements show abrupt changes in sea level across fracture zones which separate areas of lithosphere with different ages. These changes have the correct location, amplitude, and wavelength to be caused by the combined gravitational attraction of the relief across the fracture zone and the isostatic support of this relief. Eight profiles of geoid height and bathymetry across the Mendocino fracture zone are inverted to determine the depth of the isostatic compensation, assuming that the compensation occurs in a single layer. These depths are then interpreted with a thermal boundary layer model of lithospheric growth. To explain satisfactorily the geoid measurements, the thermal diffusivity of the upper mantle must be 3.3 × 10^?3 cm² s^?1 and the thickness of the lithosphere, defined as the depth at which the geotherm reaches 95% of its maximum value, must be9.1km m.y.^?1/2 × t^1/2, where t is lithospheric age.     

Constraining tephra dispersion and deposition from three subplinian explosions in 2011 at Shinmoedake volcano,Kyushu, Japan

Constraining physical parameters of tephra dispersion and deposition from explosive volcanic eruptions is a significant challenge, because of both the complexity of the relationship between tephra distribution and distance from the vent and the difficulties associated with direct and comprehensive real-time observations. Three andesitic subplinian explosions in January 2011 at Shinmoedake volcano, Japan, are used as a case study to validate selected empirical and theoretical models using observations and field data. Tephra volumes are estimated using relationships between dispersal area and tephra thickness or mass/area. A new cubic B-spline interpolation method is also examined. Magma discharge rate is estimated using theoretical plume models incorporating the effect of wind. Results are consistent with observed plume heights (6.4–7.3 km above the vent) and eruption durations. Estimated tephra volumes were 15–34?×?10⁶ m³ for explosions on the afternoon of 26 January and morning of 27 January, and 5.0–7.6?×?10⁶ m³ for the afternoon of 27 January; magma discharge rates were in the range 1–2?×?10⁶ kg/s for all three explosions. Clast dispersal models estimated plume height at 7.1?±?1 km above the vent for each explosion. The three subplinian explosions occurred with approximately 12-h reposes and had similar mass discharge rates and plume heights but decreasing erupted magma volumes and durations.     

Amplitude-dependent effects of longitudinal seismic wave propagation in the interhole space

Field investigations of the amplitude dependence of the P wave velocity in dry and water-saturated rocks are carried out in the space between two shallow boreholes. The seismic wave velocity nonlinearly varies with the strain amplitude in the range ～(4–50) × 10^?8. The pattern of the velocity variation with amplitude depends on the pulse propagation direction. In dry and partially water-saturated rocks, the wave velocity decreases by 1.5% with the amplitude increasing within the range mentioned above and increases by 0.4% in completely water-saturated rocks (with an accuracy of up to 0.1%). Amplitude variations within a closed cycle (A _min … → A _max … → A _min) lead to hysteresis in the V _p (A _min-max-min) dependence (i.e., the ascending and descending branches of the curve do not coincide). If the hysteretic loop is not closed, the residual velocity component ΔV _p (A) is present. This effect is observed in dry and weakly saturated rocks. In a completely saturated rock, hysteresis of the velocity dependence is absent; the ascending and descending amplitude branches coincide. It is suggested that the amplitude characteristics and their hysteresis can be used in the future as an additional criterion for the differentiation of rocks by their fluid saturation.     

Dynamics of ash-dominated eruptions at Vesuvius: the post-512 AD AS1a event

Recent stratigraphic studies at Vesuvius have revealed that, during the past 4,000 years, long lasting, moderate to low-intensity eruptions, associated with continuous or pulsating ash emission, have repeatedly occurred. The present work focuses on the AS1a eruption, the first of a series of ash-dominated explosive episodes which characterized the period between the two Subplinian eruptions of 472 AD and 1631 AD. The deposits of this eruption consist of an alternation of massive and thinly laminated ash layers and minor well sorted lapilli beds, reflecting the pulsatory injection into the atmosphere of variably concentrated ash-plumes alternating with Violent Strombolian stages. Despite its nearly constant chemical composition, the juvenile material shows variable external clast morphologies and groundmass textures, reflecting the fragmentation of a magma body with lateral and/or vertical gradients in both vesicularity and crystal content. Glass compositions and mineralogical assemblages indicate that the eruption was fed by rather homogeneous phonotephritic magma batches rising from a reservoir located at ~ 4 km (100 MPa) depth, with fluctuations between magma delivery and magma discharge. Using crystal size distribution (CSD) analyses of plagioclase and leucite microlites, we estimate that the transit time of the magma in the conduit was on the order of ~ 2 days, corresponding to an ascent rate of around 2 × 10⁻² ms⁻¹. Accordingly, assuming a typical conduit diameter for this type of eruption, the minimum duration of the AS1a event is between about 1.5 and 6 years. Magma fragmentation occurred in an inertially driven regime that, in a magma with low viscosity and surface tension, can act also under conditions of slow ascent.     

Three-dimensional numerical analysis of magma transport through a pre-existing fracture in the crust

Magmas are transported through pre-existing fractures in many repeatedly erupting volcanoes. The study of this special process of magma transport is fundamentally important to understand the mechanisms and conditions of volcanic eruptions. In this paper, we numerically simulate the magma propagation process through a pre-existing vertical fracture in the crust by using the combined finite difference method (FDM), finite element method (FEM) and discontinuous deformation analysis (DDA) approach. FDM is used to analyze magma flow in the pre-existing fracture, FEM is used to calculate the opening of the fracture during magma intrusion, and DDA is used to deal with the contact of the closed fracture surfaces. Both two-dimensional (2D) and three-dimensional (3D) examples are presented. Parametric studies are carried out to investigate the influence of various physical and geometric parameters on the magma transport in the pre-existing fracture. We have considered magma chamber depth ranging from 7 km to 10 km under the crust surface, magma viscosity ranging from 2 × 10⁻² to 2 × 10⁻⁷ MPa s, and the density difference between the magma and host rock ranging from 300 to 700 kg/m³. The numerical results indicate that (1) the fluid pressure p varies gradually along the depth, (2) the shape of the magma body during propagation is like a torch bar and its width ranges from 2 m to 4 m approximately in the 3D case and 10 m to 50 m in the 2D case for the same physical parameters used, (3) the crust surface around the pre-existing fracture begins to increase on both sides of the fracture, forms a trough between them, then gradually uplifts during the transport of the magma, and finally takes the shape of a crater when the magma reaches the surface. We have also examined the influence of physical and geometric parameters on the minimum overpressure for magma transport in the 3D case. The numerical results show that our numerical technique presented in this paper is an effective tool for simulating magma transport process through pre-existing fractures in the crust.     

The magma budget of Volcan Arenal, Costa Rica from 1968 to 1980

The volume of magma emitted by Volcan Arenal from July 1968 to March 1980 has been calculated to be 304 × 10⁶ m³ (dense rock equivalent). Most of this magma has been emplaced as block lava flows on the western flanks of the volcano following the initial explosive eruptions in 1968. From 1968 to 1973 the volumetric discharge rate of magma decreased from about 3-2 m³ s⁻¹ to about 1 m³ s⁻¹. During a break in activity in late 1973 the site of effusion moved from Crater A to Crater C about 400 m higher. Subsequent effusion was at a lower rate (0.3 m³ s⁻¹) which remained constant for the next six years. Comparison of dry-tilt measurements during this latter period of steady-state effusion with numerical finite-element models of Arenal's elastic response to the evacuation of magma from an underlying reservoir favor a very shallow reservoir (< 2 km depth) to explain the data. However, the constraints imposed by the measured volumes of magma are not compatible with such a reservoir. Instead, it is argued that the steady downward tilting of the volcano's summit was caused by the loading of the western side of the volcano by about 19 × 10⁶ m³ of lava. Surface loading by lava flows may be an important deformational effect at other volcanoes. A system of magma supply involving open conduits (pipes) for the uppermost one kilometer and transitory conduits (cracks) to a crustal reservoir is proposed. This crustal reservoir initially contained a compositionally graded magma which was evacuated from 1968 to 1973. The subsequent abrupt decrease in effusion rate is compatible with the increased magmatic head required to reach Crater C. The constancy of magma composition and effusion rate from 1974 to 1980 implies a homogeneous magma reservoir.     

Theoretical and Observed Depth Correction for M s

—?Modal summation technique is used to generate 5000, three-component theoretical seismograms of Love and Rayleigh waves, assuming modified PREM (PREM-C) and AK135F global earth models. The focal depth h and the geometrical fault parameters are randomly chosen so as to uniformly cover possible source mechanisms and obtain uniform distribution of log h in the interval 1?h?h?M _s of the form:¶ΔM _s (h)=0 forh< 20km, ΔM _s (h)=0.314log(h)-0.409 for 20≠h< 60km, ΔM _s (h)=1.351log(h)-2.253 for 60≠h< 100km, ΔM _s (h)=0.400log(h)-0.350 for 100≠h< 600km .¶After applying the above correction, the relationship between the surface wave magnitude and the scalar seismic moment for the observational data set significantly improves, and becomes independent of the source depth. In relation to CTBT, no depth correction is needed for M _S when the m _b ???M _S discriminant is computed, because the proposed correction is zero for earthquakes with foci above 20?km.