Causes of microgravity change at Poás volcano,Costa Rica: an active but non-erupting system

Microgravity changes with time, not always consistent with the Bouguer-corrected free air gradient, have been recorded and associated with cyles of eruptive activity at Krafla, Kilauea, Pacaya and Etna volcanoes. In contrast, over the non-erupting yet active fumarolic vents at Poás (Costa Rica), real-time gravity observations over three periods during the years 1983–1985 have identified ca. 140µGal amplitude, cyclic gravity variations. Their decrease in amplitude with distance from the active crater, coupled with a static sub-surface structural model, have allowed the effects of a variety of possible causative dynamic phenomena to be evaluated. It is concluded that cyclic changes in the average density (by ca. 0.03 Mg m^–3) in the magma pipe at depths below 500 m, rather than variation in magma chamber or water table geometry, are responsible for the observed gravity variations. Specifically, anaverage of 1 % fluctuation in the volume of gas in crystal-free magma, a process driven by thermal convection cycles, probably accounts for the density/gravity change.     

The evolution of andesite volcano structures: new evidence from gravity studies in Costa Rica

Published gravity data on active volcanoes generally reflecteither the low density scoriaceous/pumiceous deposits that are localized within ring-fracture collapse depressions, such as the calderas of mature silicic volcanoes,or the high density frozen magma conduits that occur beneath basaltic shields and cones. The intensive gravity surveys reported here over three complex andesite volcanoes reveal features of both types. Their multi-component gravity fields have crater-centred positive anomalies (1–2 km diameter) surrounded by broader zones of negative gravity with similar amplitudes but greater width (5–10 km). The former are thought to reflect sub-crater magma pipes ofnormal density (ca. 2.5–2.6 Mg m⁻³) surrounded by pyroclastic scoria, ashes and occasional lava flows of muchlower net density (1.8–2.4 Mg m⁻³) which, in turn, account for the negative anomalous zones because the deeper, more consolidated and older parts of these andesite volcano edifices have more normal densities (2.3–2.6 Mg m⁻³).The low density materials are particularly interesting because they appear to have filled topographic depressions to depths of several hundred metres, especially where old caldera-like structures have been postulated from the steep gravity gradients over perimeter ring faults. A model is developed whereby short periods of caldera collapse, associated with intermittent, large high level magma bodies, are interspersed by normal crater-like activity with narrow sub-surface magma pipes. Dominantly pyroclastic activity from summit craters generates the materials that gradually fill earlier-formed topographic depressions. This study demonstrates the unique value of detailed gravity surveys, combined with surface geological information, for modelling and understanding the evolution of active volcano summit regions.     

Magnetic monitoring at Merapi volcano, Indonesia

Merapi volcano, located 30 km north of the heavily populated city of Yogjakarta, Java, is one of the most active of the 129 volcanoes in Indonesia. About every 2 years a new phase of activity is observed. Depending on the past activity the unrest gives rise either to an endogenous dome which partly collapses in the southwest direction or to pyroclastic flows which travel as far as 15 km. The 1990–1997 period has involved a plume emission on 30 August 1990, an extrusion on 20 January 1992, and a pyroclastic eruption on 22 November 1994. The intensity of the Earth magnetic field has been measured simultaneously and digitally recorded at four stations since 1990. Two Overhauser magnetometers with resolution of 0.01 nT have been installed in the summit area to strengthen the volcano monitoring. Outstanding magnetic changes appear to correlate with volcanic activity. Three types of volcanomagnetic signals can be identified: long-term trends up to 15 nT with period >10 years; medium-term cyclic variations, at most 3 nT in amplitude and with 1–2 years period; and small events, reaching 1.5 nT, lasting a few months, and associated with any remarkable volcanic activity. Merapi volcano began a new cycle of activity in 1995 leading to a dome growth in July 1996, and accompanied by 27 nuées ardentes in August. The comparison between magnetic data, seismicity, and surface phenomena suggests that some long-term trends of decade periods could be of thermomagnetic origin, while mid-term volcanomagnetic variations associated with the cycles of Merapi activity could be of piezomagnetic origin. Short-term variations of a few weeks duration, less than 1.5 nT, are well correlated with the 1995–1996 seismic activity.     

Post-eruptive volcanic dome evolution as revealed by deformation and microgravity observations at Usu volcano (Hokkaido, Japan)

Usu volcano (Hokkaido, Japan) is a dacitic volcano, known for its high production rate of lava domes and crypto-domes. It is thus a good target to study processes of volcanic dome evolution (upheaval and/or relaxation). We carried out repeated GPS and microgravity surveys on the three most recent domes of Mt. Usu (1910: Meiji Shinzan; 1943–1945: Showa-Shinzan and 1977–1982: Usu-Shinzan). The repeat period was 1 to 2 months and extended from October 1996 to June 1997. We also compare new data with results from former studies. More than 20 years after the start of Usu-Shinzan dome growth, there is still subsidence at a maximum rate of about 7 to 8 cm/year. The reasons for this subsidence are discussed. Repeated gravity surveys revealed an increase of gravity on the domes (about 60±10 microgal/year for Usu-Shinzan, about 15 microgal at Showa-Shinzan and 10 to 20 microgal for Meiji-shinzan); this gravity increase exceeds that expected due to subsidence. We discuss and interpret the excess gravity change in terms of a density increase in the edifice, caused by a combination of processes (contraction of the edifice, water level change, devesiculisation, cooling and magma intrusion). Quantification of these processes at Usu volcano may help to understand the processes of evolution at domes on other volcanoes such as Merapi (Indonesia), Unzen (Japan) or Montserrat (West Indies).     

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.     

Lahars at Merapi volcano, Central Java: an overview

Merapi volcano, in Central Java, is one of the most active volcanoes in the world. At least 23 of the 61 reported eruptions since the mid-1500s have produced source deposits for lahars. The combined lahar deposits cover about 286 km² on the flanks and the surrounding piedmonts of the volcano. At Merapi, lahars are commonly rain-triggered by rainfalls having an average intensity of about 40 mm in 2 h. Most occur during the rainy season from November to April, and have average velocities of 5–7 m/s at 1000 m in elevation. A wide range of facies may be generated from a single flow, which may transform downvalley from debris flow to hyperconcentrated streamflow.Because of the high frequency and magnitude of the lahar events, lahar-related hazards are high below about 450–600 m elevation in each of the 13 rivers which drain the volcano. Hazard-zone maps for lahar were produced by Pardyanto et al. (Volcanic hazard map, Merapi volcano, Central Java (1/100,000). Geol. Surv. of Indonesia, Bandung, II, 4, 1978) and the Japanese–Indonesian Cooperation Agency (Master plan for land conservation and volcanic debris control in the area of Mt Merapi, Jakarta, 1980), but these maps are of a very small scale to meet modern zoning requirements. More recently, a few large-scale maps (1/10,000- and 1/2000-scale) and risk assessments have been completed for a few critical river systems.     

On the interpretation of vertical gravity gradients produced by magmatic intrusions

A method for the computation and interpretation of gravity and height changes and vertical gravity gradients produced by magmatic intrusions in a layered elastic–gravitational medium is presented. The methodology assumes a planar medium geometry, which consists of welded layers overlying a half-space. The medium is elastic and gravitating. The intrusion is treated as a point source and can be located at any depth inside the medium. The theoretical elastic–gravitational model allows the computation of the so-called geometric and orthometric vertical displacements as well as gravity changes of different types. The corresponding vertical gravity gradients can also be computed. We present several examples of theoretical computations and study different types of these geometric and orthometric gravity gradients and the information we can get from the application of the methodology described. The results presented show that the use of these gradients is a useful tool to obtain information on the dynamics of the injection processes, including the detection of new magma recharge. We show that using the elastic–gravitational deformation model we can explain non-linear gravity–height relationships that appear in volcanic areas. We also present the application of the methodology to Mayon volcano, Philippines, delineating the intrusion of new magma, consistent with the produced eruptions after the observation period.     

Combined discrete and continuous gravity observations at Mount Etna

Systematic investigation of discrete gravity measurements has continued at Mount Etna since 1986. The network now covers an area of 400 km² with about 70 stations 0.5–3 km apart. Mass redistributions occurring at depths ranging between about 8 km below sea level and a few hundred metres below the surface (magma level changes within the shallower parts of the feeding conduits) have been identified from these data. Conventional (discrete) microgravity monitoring on a network of stations furnishes only instantaneous states of the mass distribution at continuously active systems. In order to obtain information on the rate at which the volcanic processes (and thus mass transfers) occur, three stations for continuously recording gravity where installed on Mount Etna in 1998. A 16-month long sequence from one of the continuously running stations (PDN, located 2 km from the active northeast crater at the summit of Etna volcano) is presented. After removing the effects of Earth Tide and tilt, the correlation of the residual gravity sequence with simultaneous recordings of meteorological parameters acquired at the same station was analysed. Once the meteorological effects have also been removed, continuous gravity changes are within 10 μGal of gravity changes measured using conventional microgravity observations at sites very close to the continuous station. This example shows how discrete and continuous gravity observations can be used together at active volcanoes to get a fuller and more accurate picture of the spatial and temporal characteristics of volcanic processes.     

Microgravity and height changes caused by volcanic activity: Four

Microgravity measurements and levelling surveys on volcanoes are not always easy to make, but are useful for studying volcanic processes quantitatively. Gravity changes associated with volcanic activity are not always significant. Precision of microgravity measurements depend critically on the procedures adopted, and those applied in the present paper are described. Levelling technique is now orthodox, and some empirical laws relating ground deformation to volcanic activity are deduced from the accumulated data. Gravity changes occur at the same time and places as ground deformations. The relationship between microgravity and height changes are discussed from the standpoint of analyzing the data obtained on volcanoes. The observational results obtained on four volcanoes in Japan are separately analyzed because each volcano exhibits different patterns of gravity changes and deformations. During the 1977–1982 activity of Usu volcano, deformation was accompanied by microgravity changes frequently observed at a particular benchmark at the base of the volcano for about five years. The gravity changes prove to be not a direct effect of magma movements but to be caused by the deformations of ground strata and aquifers around the benchmark. The 1983 eruption of Miyakejima volcano was associated with local gravity changes around the eruptive fissures due to magma intrusion which was approximately modelled. Similarly the 1986 eruption of Ooshima volcano caused gravity changes on the volcano, but these were poorly correlated with elevation changes and their origins were not uniquely interpreted. To detect gravity changes associated with the activity of Sakurajima volcano, an equigravity point was selected at the north of the volcano besides the gravity points on and around the volcano itself. The probable gradual accumulation of magmas beneath the volcano for eight years is substantiated by observed microgravity and elevation changes.     

Interpretation of 1992–1994 Gravity Changes around Mayon Volcano,Philippines, Using Point Sources

Significant gravity changes observed around the Mayon Volcano (Philippines) between 1992 and 1994 at 26 stations are interpreted in terms of an increase of mass and pressure changes at several point sources modelled using a fast inversion process. This inversion approach attempts to fit gravity and elevation changes by combining a random search for the positions of the sources and a linear least-squares fit for the incremental mass, pressure and possible common regional values for gravity or elevation changes. Some stabilizer terms are included in the misfit function. Models with one and two sources were tested against the observed changes at Mayon. Models with only one-source give a best fit for a shallow source with a positive mass increment, horizontally displaced far from the summit. The study using two sources gives a best fit that is similar to the one-source model, but in addition indicates anomalous behavior at stations in the SW. Neglecting the stations located southward from a local fracture, the best-fitting model suggests one central positive mass change source, which is likely to be an intrusion of about 0.5 MU with a depth of about 5 km beneath the volcano. Standard deviation for the residuals ranges from 7–8 μGal for one-source models to 6–7 μGal for models with two sources. Both of the cases are below the error value of 9.4 μGal estimated for the gravity data, so that it is not possible to discriminate between both possible interpretations without additional information.     

Temporal variations in magma composition at Merapi Volcano (Central Java, Indonesia): magmatic cycles during the past 2000 years of explosive activity

Merapi Volcano (Central Java, Indonesia) has been frequently active during Middle to Late Holocene time producing basalts and basaltic andesites of medium-K composition in earlier stages of activity and high-K magmas from 1900 ¹⁴C yr BP to the present. Radiocarbon dating of pyroclastic deposits indicates an almost continuous activity with periods of high eruption rates alternating with shorter time spans of distinctly reduced eruptive frequency since the first appearance of high-K volcanic rocks. Geochemical data of 28 well-dated, prehistoric pyroclastic flows of the Merapi high-K series indicate systematic cyclic variations. These medium-term compositional variations result from a complex interplay of several magmatic processes, which ultimately control the periodicity and frequency of eruptions at Merapi. Low eruption rates and the absence of new influxes of primitive magma from depth allow the generation of basaltic andesite magma (56–57 wt% SiO₂) in a small-volume magma reservoir through fractional crystallisation from parental mafic magma (52–53 wt% SiO₂) in periods of low eruptive frequency. Magmas of intermediate composition erupted during these stages provide evidence for periodic withdrawal of magma from a steadily fractionating magma chamber. Subsequent periods are characterised by high eruption rates that coincide with shifts of whole-rock compositions from basaltic andesite to basalt. This compositional variation is interpreted to originate from influxes of primitive magma into a continuously active magma chamber, triggering the eruption of evolved magma after periods of low eruptive frequency. Batches of primitive magma eventually mix with residual magma in the magmatic reservoir to decrease whole-rock SiO₂ contents. Supply of primitive magma at Merapi appears to be sufficiently frequent that andesites or more differentiated rock types were not generated during the past 2000 years of activity. Cyclic variations also occurred during the recent eruptive period since AD 1883. The most recent eruptive episode of Merapi is characterised by essentially uniform magma compositions that may imply the existence of a continuously active magma reservoir, maintained in a quasi-steady state by magma recharge. The whole-rock compositions at the upper limit of the total SiO₂ range of the Merapi suite could also indicate the beginning of another period of high eruption rates and shifts towards more mafic compositions.     

New Results at Mayon,Philippines, from a Joint Inversion of Gravity and Deformation Measurements

In this paper, we detail the combination of the genetic algorithm (GA) inversion technique with the elastic-gravitational model originally developed by Rundle and subsequently refined by Fernández and others. A sensitivity analysis is performed for the joint inversion of deformation and gravity to each of the model parameters, illustrating the importance of proper identification of both the strengths and limitations of any source model inversion, and this technique in particular. There is a practical comparison of the theoretical results with the inversion of geodetic data observed at the Mayon volcano in the Philippines, where there are gravity changes without significant deformation, after the 1993 eruption.     

10,000 Years of explosive eruptions of Merapi Volcano, Central Java: archaeological and modern implications

Stratigraphy and radiocarbon dating of pyroclastic deposits at Merapi Volcano, Central Java, reveals 10,000 years of explosive eruptions. Highlights include:(1) Construction of an Old Merapi stratovolcano to the height of the present cone or slightly higher. Our oldest age for an explosive eruption is 9630±60 ¹⁴C y B.P.; construction of Old Merapi certainly began earlier.(2) Collapse(s) of Old Merapi that left a somma rim high on its eastern slope and sent one or more debris avalanche(s) down its southern and western flanks. Impoundment of Kali Progo to form an early Lake Borobudur at 3400 ¹⁴C y B.P. hints at a possible early collapse of Merapi. The latest somma-forming collapse occurred 1900 ¹⁴C y B.P. The current cone, New Merapi, began to grow soon thereafter.(3) Several large and many small Buddhist and Hindu temples were constructed in Central Java between 732 and 900 A.D. (roughly, 1400–1000 ¹⁴C y B.P.). Explosive Merapi eruptions occurred before, during and after temple construction. Some temples were destroyed and (or) buried soon after their construction, and we suspect that this destruction contributed to an abrupt shift of power and organized society to East Java in 928 A.D. Other temples sites, though, were occupied by “caretakers” for several centuries longer.(4) A partial collapse of New Merapi occurred <1130±50 ¹⁴C y B.P. Eruptions 700–800 ¹⁴C y B.P. (12–14th century A.D.) deposited ash on the floors of (still-occupied?) Candi Sambisari and Candi Kedulan. We speculate but cannot prove that these eruptions were triggered by (the same?) partial collapse of New Merapi, and that the eruptions, in turn, ended “caretaker” occupation at Candi Sambisari and Candi Kedulan. A new or raised Lake Borobudur also existed during part or all of the 12–14th centuries, probably impounded by deposits from Merapi.(5) Relatively benign lava-dome extrusion and dome-collapse pyroclastic flows have dominated activity of the 20th century, but explosive eruptions much larger than any of this century have occurred many times during Merapi's history, most recently during the 19th century.Are the relatively small eruptions of the 20th century a new style of open-vent, less hazardous activity that will persist for the foreseeable future? Or, alternatively, are they merely low-level “background” activity that could be interrupted upon relatively short notice by much larger explosive eruptions? The geologic record suggests the latter, which would place several hundred thousand people at risk. We know of no reliable method to forecast when an explosive eruption will interrupt the present interval of low-level activity. This conclusion has important implications for hazard evaluation.     

Rate of sediment yield following small‐scale volcanic eruptions: a quantitative assessment at the Merapi and Semeru stratovolcanoes,Java, Indonesia

Sediment yields were calculated on the ?anks of Merapi and Semeru volcanoes in Java, Indonesia, using two different methods. During the ?rst year following the 22 November 1994 eruption of Merapi, a sediment yield in excess of 1·5 × 10⁵ m³ km^?2 yr^?1 was calculated in the Boyong River drainage basin, based on the volumes of sediment that were trapped by ?ve check dams. At Semeru, sediment discharges were assessed in the Curah Lengkong River from direct measurements on the lahars in motion and on the most signi?cant stream?ows. The calculated rate of sediment yield during one year of data in 2000 was 2·7 × 10⁵ m³ km^?2 yr^?1. Sediment yields are dominated by rain‐triggered lahars, which occur every rainy season in several drainage basins of Merapi and Semeru volcanoes, mostly during the rainy season extending from October to April. The return period of lahars carrying sediment in excess of 5 × 10⁵ m³ is about one year in the Curah Lengkong River at Semeru. At Merapi, the volume of sediments transported by a lahar did not exceed 2·8 × 10⁵ m³ in the Boyong River during the rainy season 1994–95. On both volcanoes, the sediments are derived from similar sources: pyroclastic‐?ow/surges deposits, rockfalls from the lava domes, and old material from the riverbed and banks. However, daily explosions of vulcanian type at Semeru provide a more continuous sediment supply than at Merapi. Therefore, sediment yields are larger at Semeru. Copyright © 2004 John Wiley & Sons, Ltd.     

Temporal gravity and elevation changes at Pacaya volcano, Guatemala

Repetitive gravity surveys at Pacaya Volcano from 1975 to 1979 revealed time-dependent changes in the gravity field, which although related to volcanic activity, could not be uniquely attributed to elevation changes or mass redistributions because elevation control was lacking. Elevation control was established in July 1979 using precision leveling. Relative elevation and gravity measurements in June and July of 1979, January 1980 and June 1980 indicate concurrent gravity and elevation changes contemporaneous with variations in eruptive activity. From June 1979 to January 1980, while fumarolic activity was dominant, relative to the most remote station, the volcano deflated by at least 195 mm and the gravity field increased by up to 221 μgal. From January 1980 to June 1980, preceding a Strombolian eruption beginning about June 1980, the volcano inflated by at least 19 mm and the gravity field decreased by up to 231 μgal. Gravity change maps for the intervals of January 1978 to June 1979, June 1979 to January 1980, and January 1980 to June 1980 show areas subject to repeated positive and negative gravity change. Some of those areas coincide with areas of maximum elevation change observed in the June 1979–January 1980 and January 1980–June 1980 intervals; however, gravity changes were observed in areas lacking elevation changes. Adjusting observed gravity changes for elevation changes using a free-air value of −3.086 μgal/cm does not substantially alter the pattern, position, or amplitude of the gravity changes. The relationship between gravity changes, elevation changes, and volcanic activity requires a mechanism producing gravity decreases with little inflation during times of increased eruptive activity, and producing gravity increases with subsidence during times of declining eruptive activity. Such a pattern of changes could be produced by a near-surface magma body in which high-density degassed magma is displaced volume for volume by low-density vesiculated magma during time of increased eruptive activity, and in which loss of gasses by fumarolic activity produces a density increase and a reduction in volume of the magma body during periods of declining eruptive activity. Such a pattern of changes could also be induced by a low-density, vesiculated magma body moving upward in the volcanic pile by piecemeal stoping where the high-density rocks of the volcano are replaced on a volume for volume basis by low-density magma during periods of increasing eruptive activity; and by later density increases and magma body volume reductions accompanying devolatilization and devesiculation during periods of declining eruptive activity. Simple density change and density contrast models involving shallow magma bodies at depths of 100 to 200 m indicate density changes or contrasts of about 0.4 g/cm³ could produce the gravity changes.     

Behavioral patterns of Fuego volcano, Guatemala

The times of activity at Fuego (one of the most active volcanoes in the world) since 1800 correlate with the activity of other Central American volcanoes. Approximately 0.7 km³ of olivine-bearing, high-Al₂O₃ basalt has been erupted since 1932, and about 1.7 km³ has been produced during 450 years of historic records. A minimum of 13,000 years and a maximum of 100,000 years were required to build Fuego's cone of 50 km³. Within the recent cluster of activity since 1932, rates of magma production have increased to 0.5 m³/s and the trend has been toward more eruptions (shorter reposes) of progressively more mafic basalt. 47% of the eruptions occurred within 2 days of the fortnightly tidal maximum and 56% occurred within 2 hours of the semi-diurnal minimum of the vertical tidal gravity acceleration. Thus the maximum compressional component of the tidal cycles can trigger an eruption at Fuego. Eruptions with higher effusion rates produce larger volumes of materials, although they only last a few hours. The 20–70 year clusters of activity beginning at 80–170-year intervals are interpreted as reflecting the ascent of primary batches of magma. A deeper (8–16 km), larger (> 1 km³) primary chamber and a shallower (2–5 km), smaller (0.1 km³), dike-like secondary chamber best explain Fuego's behavioral pattern.     

Magma movements in Etna volcano associated with the major 1991–1993 lava eruption: evidence from gravity and deformation

The 1991–1993 eruption was probably the largest on Mt. Etna for 300 years. Since then the volcano has entered an unusually quiescent period. A comprehensive record of gravity and ground deformation changes presented here bracket this eruption and give valuable insight into magma movements before, during and after the eruption. The gravity and deformation changes observed before the eruption (1990–1991) record the intrusion of magma into the summit feeder and the SSE-trending fracture system which had recently been active in 1978, 1979, 1983 and 1989, creating the feeder dyke for the 1991–1993 eruption. In the summit region gravity changes between 1992 and 1993 (spanning the end of the eruption) reflect the withdrawal of magma from the conduit followed more recently (1993–1994) by the re-filling of magma in the conduit up to pre-eruption levels. In contrast, in the vicinity of the fracture zone, gravity has remained at the 1991–1992 level, indicating that no withdrawal has occurred here. Rather, magma has solidified in the fracture system and sealed it such that the 1993–1994 increase in magma level in the conduit was not accompanied by further intrusion into the flanks. Mass calculations suggest that a volume of at least 10⁷ m³ of magma has solidified within the southeastern flank of the volcano.     

Ground deformation and gravity changes accompanying the March 1981 eruption of Mount Etna

We present reults from simultaneous precise levelling and gravity surveys on Mount Etna covering the period August 1980–August 1981. The flank eruption of March 1981 erupted 18–35 × 10⁵m³ of lava. Following it, upward movements of more than 17 cm were observed close to the new fissure and a broad, apparently independent, uplift of 5 cm was observed 4 km to the west. A zone of about 2 cm depression to the east of the fissure is insufficient to account for the volume of magma erupted. Gravity results show positive changes of up to 63 microgal, and display good positive correlation with elevation changes. Both sets of measurements appear to be due to new intrusion of magma rather than subsurface magma drainage. Ground deformation close to the new fissure is well modelled by intrusion of a dyke in the zone 100–500 m below the surface, striking along the fissure and of dip between 75–90°. The gravity changes are modelled as due to a deeper intrusion of magma, along the same line but some 1500 m below the surface. The changes were not present immediately after the eruption but occurred during the ensuing 5 months. It is proposed that this introduction of matter occurred by crack propagation along the fissure in the aftermath of the eruption. Towards the west of the fissure, and some 4 km west of the summit, ground deformation is modelled by intrusion of a dyke in the zone 300–1500 m below the surface and dipping at 80–85°. Again, gravity changes appear to be due to magma intrusion at greater depth, close to sea level. In this case gravity changes are interpreted as due to magma density changes, as a result of pressure increase in a larger scale fissure zone. This same pressure increase may be forcing the new intrusion close to the surface, and makes this part of the volcano a region of especially high risk.     

Modeling of Stress Changes at Mayon Volcano,Philippines

In this paper, we analyze the stress change associated with the inverted volcanic source models at Mayon volcano, Philippines, where there are gravity changes without significant deformation after the 1993 eruption. We detail the applicable data and the associated inversion techniques and models prior to calculating the appropriate stress changes. It is determined that the stress change associated with the central magmatic source produces compressional stress changes at a secondary source to the northwest, prompting a change in the local water storage in the underlying fractured rock.     

The first absolute gravity measurements in Indonesia

For the purposes of the calibration of the superconducting gravimeter (SG) in Bandung and the establishment of the absolute gravity (AG) points, we carried out AG measurements for the first time in Indonesia in November 2002. The measurements in Bandung were conducted between November 15th and 20th by means of a FG5 (#210), and 14,520 effective drops were obtained. The gravity value newly determined at the AG point in Bandung is 977976701.2 μgal (1 μgal = 10⁻⁸ ms⁻²) and the scale factor for the SG is −52.22 μgal/V. We also established another AG point in Yogyakarta near Merapi volcano and carried out AG measurements in Yogyakarta between November 22nd and 26th. The gravity value determined for this station is 978203093.5 μgal.