S-wave shadows in the Krafla Caldera in NE-Iceland,evidence for a magma chamber in the crust

During the present tectonic activity in the volcanic rift zone in NE-Iceland it has become apparent that the attenuation of seismic waves is highly variable in the central region of the Krafla volcano. Earthquakes associated with the inflation of the volcano have been used to delineate two regions of high attenuation of S-waves within the caldera. These areas are located near the center of inflation have horizontal dimensions of 1–2 km and are interpreted as the expression of a magma chamber. The top of the chamber is constrained by hypocentral locations and ray paths to be at about 3 km depth. Small pockets of magma may exist at shallower levels. The bottom of the chamber is not well constrained, but appears to be above 7 km depth. Generally S-waves propagate without any anomalous aftenuation through laver 3 (v_p=0.5 km sec^?1) across the volcanic rift zone in NE-Iceland. The rift zone therefore does not appear to be underlain by an estensive magma chamber at crustal levels. The Krafla magma chamber is a localized feature of the Krafla central volcano.     

Seismic constraints on magma chambers at Hekla and Torfajökull volcanoes,Iceland

Hekla and Torfajökull are active volcanoes at a rift–transform junction in south Iceland. Despite their location next to each other they are physically and geologically very different. Hekla is an elongate stratovolcano, built mainly of basaltic andesite. Torfajökull is a prominent rhyolitic centre with a 12-km-diameter caldera and extensive geothermal activity. The scope of this study is to examine the propagation of body waves of local earthquakes across the Hekla–Torfajökull area and look for volumes of anomalous S-wave attenuation, which can be evidence of magma chambers. So far the magma chamber under Hekla has been modelled with various geophysical means, and its depth has been estimated to be 5–9 km. A data set of 118 local earthquakes, providing 663 seismic rays scanning Hekla and Torfajökull, was used in this study. The major part, 650 seismograms, did not show evidence for S-wave attenuation under these volcanoes. Only six seismograms had clear signs of S-wave attenuation and seven seismograms were uncertain cases. The data set samples Hekla well at depths of 8–14 km, and south part of it also at 4–8 km and 14–16 km. Western Torfajökull is sampled well at depths of 4–14 km, eastern and southern Torfajökull at 6–12 km. Conclusions cannot be drawn regarding the existence of magma beyond these depth ranges. Also, magma volumes of smaller dimensions than about 800 m cannot be detected with this method. If a considerable molten volume exists under Hekla, it must be located either above 4 km or below 14 km. The former possibility seems unlikely, because Hekla lacks geothermal activity and persistent seismicity, usually taken as expressions of a shallow magma chamber. An aseismic volume with a diameter of 4 km at the depth of 8 km in the west part of Torfajökull has been inferred in earlier studies and interpreted as evidence for a cooling magma chamber. Our results indicate that this volume cannot be molten to a great extent because S-waves travelling through it are not attenuated. Intense geothermal activity and low-frequency earthquakes are possibly signs of magma in the south part of Torfajökull, but a magma chamber was not detected there in the areas sampled by this study.Editorial responsibility: T. Druitt     

Modelling of surface ground deformations in the Phlegraean Fields volcanic area,Italy

A preliminary finite elements model of the ground deformations observed at Phlegraean Fields is proposed. The model assumes an oblate-spheroid magma chamber at the depth of 5.4 km with major semiaxis of 1.5 km and minor semiaxis of 0.75 km. The dimensions of the magma chamber have been evaluated by using a thermal model based on the assumptions that a progressively cooling huge magmatic body is responsible for the volcanic activity at Phlegraean Fields in the last 35,000 years. Surface deformations caused by an over-pressure of 30 MPa in the magma chamber have been calculated. Constant, and temperature-dependent elastic parameters of the surrounding medium have been considered. Vertical displacements of the order of those presently observed at Phlegraean Fields can be obtained only with temperature-dependent elastic properties of the medium.     

Magma contamination by direct wall rock interaction: constraints from xenoliths from the walls of a carbonate-hosted magma chamber (Vesuvius 1944 eruption)

During the 1944 eruption of Vesuvius different types of xenoliths were ejected. They represent fragments of the walls of a low volume (<0.5 km³) shallow (3–4 km depth) magma chamber. The study of these xenoliths enables us to estimate the amount of contamination occurring at the boundary of a high-T alkaline magma chamber hosted in carbonate rocks. The process of contamination of the magma by carbonates can be modelled, using isotopic and chemical data, as a mixing between magma and marbles. Mass exchanges occur at the boundary between the crystallizing magma and marble wall rocks, where endoskarn forms. The contamination of the solidification front of the chamber is very limited. The solidification front and the skarn shell effectively isolate the interior of the magma chamber from new inputs of contaminants from the carbonate wall rocks. Therefore, the main volume of magma, hosted in the magma chamber, did not undergo any significant mass exchange with the wall rocks.     

A geomechanical interpretation of the local seismicity related to eruptions and renewed activity on Tolbachik,Koryakskii, and Avacha Volcanoes,Kamchatka, in 2008–2012

The local seismicity during the 2012–2013 eruption of Tolbachik Volcano and the 2008–2009 steam–gas eruption of Koryakskii Volcano is here considered as resulting from injections of magma that produced dikes, sills, and renewed activity at preexisting faults. We identified plane-oriented earthquake clusters in order to reveal the above zones using earthquake catalogs made at the Kamchatka Branch of the Geophysical Service of the Russian Academy of Sciences (KB GS RAS). Subsequent space–time analysis of these observations lends itself to the following interpretation. The November 27, 2012 Tolbachik lava eruption was preceded by an injection of magma resulting in a series of dikes trending west-northwestward in the range of absolute depths between–4 and +3 km in a zone situated southeast of the Ploskii Tolbachik Volcano edifice. The dikes penetrated into a nearly horizontal permeable zone at an absolute depth of approximately zero, producing sills and emplacing a magma-conducting dike along the top of the zone of cinder cones (the dip angle is 50° toward the azimuth 300°) 5.5 km from the epicenter of the initial magma injection. The summit steam–gas eruption of Koryakskii Volcano in 2008–2009 was preceded by magma filling a crustal chamber (the top of the chamber is at–3 km absolute depth; the chamber is 2.5 km across) close to the southwestern base of Koryakskii. Further, magma injection in a nearly north–south zone (7.5 by 2.5 km), the absolute depth between–2 and–5 km) in the north sector of Koryakskii Volcano was occurring concurrently with the summit steam–gas eruption. The injection of magma into the cone of Avacha Volcano (2010) produced sills (at altitudes between +1600 and +1900 m) and dikes (mostly striking northwest).     

长白山天池火山岩浆系统分析

针对外媒报道长白山天池火山在近2年内有可能喷发的言论,在长白山天池火山区布测了一条长度约为103 km的二维大地电磁测深观测剖面,对火山区深部电性结构进行探测研究.由于研究区内不明来源的电磁干扰非常强,对数据采用了远参考处理、Robust处理、Rhoplus分析、张量阻抗分解和基于一维层状介质电阻率与相位互算方法等先进处理技术,获取到一批在强干扰区质量较为可靠的电磁数据,利用数据计算分析了长白山天池火山区二维构造走向和感应矢量特征,采用NLCG二维反演技术对资料进行了二维反演解释,并将反演结果与前人探测结果进行了对比分析.探测结果表明:在天池火山口下方存在明显的直立型岩浆通道,岩浆通道在下方约5~8 km位置形成关闭;在火山口下方往北方向附近,在埋深位置约7 km深处存在一个明显的低阻异常体,电阻率小于10 Ωm,且与岩浆通道对接,推测其可能是地表浅部发育的岩浆囊;在长白山山门附近C07-C09号测点之间和C04-C05号测点之间,在埋深约7~17 km深处发现近直立型低阻带,低阻带与下方低阻体直接相连,推测低阻带内赋存有活动的岩浆;随着埋深的增加,从天池火山口南部约20 km位置往北方向,在埋深13~30 km之间壳内广泛发育明显的低阻异常体,推测其可能是活动的岩浆囊.反演结果与前人探测结果整体电性特征相似,但又局部不同.     

A model for the origin of large silicic magma chambers: precursors of caldera-forming eruptions

The relatively low rates of magma production in island arcs and continental extensional settings require that the volume of silicic magma involved in large catastrophic caldera-forming (CCF) eruptions must accumulate over periods of 10 ⁵ to 10 ⁶ years. We address the question of why buoyant and otherwise eruptible high-silica magma should accumulate for long times in shallow chambers rather than erupt more continuously as magma is supplied from greater depths. Our hypothesis is that the viscoelastic behavior of magma chamber wall rocks may prevent an accumulation of overpressure sufficient to generate rhyolite dikes that can propagate to the surface and cause an eruption. The critical overpressure required for eruption is based on the model of Rubin (1995a). An approximate analytical model is used to evaluate the controls on magma overpressure for a continuously or episodically replenished spherical magma chamber contained in wall rocks with a Maxwell viscoelastic rheology. The governing parameters are the long-term magma supply, the magma chamber volume, and the effective viscosity of the wall rocks. The long-term magma supply, a parameter that is not typically incorporated into dike formation models, can be constrained from observations and melt generation models. For effective wall-rock viscosities in the range 10 ¹⁸ to 10 ²⁰ Pa s ^–1, dynamical regimes are identified that lead to the suppression of dikes capable of propagating to the surface. Frequent small eruptions that relieve magma chamber overpressure are favored when the chamber volume is small relative to the magma supply and when the wall rocks are cool. Magma storage, leading to conditions suitable for a CCF eruption, is favored for larger magma chambers (>10 ² km ³) with warm wall rocks that have a low effective viscosity. Magma storage is further enhanced by regional tectonic extension, high magma crystal contents, and if the effective wall-rock viscosity is lowered by microfracturing, fluid infiltration, or metamorphic reactions. The long-term magma supply rate and chamber volume are important controls on eruption frequency for all magma chamber sizes. The model can explain certain aspects of the frequency, volume, and spatial distribution of small-volume silicic eruptions in caldera systems, and helps account for the large size of granitic plutons, their association with extensional settings and high thermal gradients, and the fact that they usually post-date associated volcanic deposits.     

Plinian-type eruptions in the medicine lake highland, california, and the nature of the underlying magma

Contemporaneous Plinian eruptions of rhyolite pumice from Glass Mountain and Little Glass Mountain during the last 1100 years B.P. were followed by extrusion of lava flows. 1.2 km of material was erupted and 10% by volume is tephra. All of the tephra deposits consist of very poorly sorted coarse ash and lapilli that are mostly pumice pyroclasts.Eruptive sequences, chemical composition and petrographic character of the rhyolites at Little Glass Mountain and Glass Mountain suggest that they came from the same magma body. The 1:9 ratio of tephra to lavas is typical of small silicic magma chambers. Eruption from a small chamber, 4–6 km deep, at vents 15 km apart is possible if magma rose along cone sheets with dips of 45–60°. The caldera rim and arcuate lines of vents near it may represent the surface expression of several concentric cone sheets.Pumice pyroclasts erupted at Glass Mountain and Little Glass Mountain may have formed in the following manner: (1) vesicle growth and coalescence beginning at 1–2 km depths; (2) elongation of the vesicles by flow within the cone sheets; (3) disruption of the vesiculated magma when it reached the surface by an expansion wave passing down through it; and (4) eruption of comminution products as pumice pyroclasts. Plinian activity at Little Glass Mountain and Glass Mountain continued until the volatile-rich top of the magma chamber had been depleted.     

Multiple volcano deformation sources in a post-rifting period: 1989–2005 behaviour of Krafla,Iceland constrained by levelling,tilt and GPS observations

The Krafla rifting episode, which occurred in North Iceland in 1975–1984, was followed by inflation of a shallow magma chamber until 1989. At that time, gradual subsidence began above the magma chamber and has continued to the present at a declining rate. Pressure decrease in a shallow magma chamber is not the only source of deformation at Krafla, as other deformation processes are driven by exploitation of two geothermal fields, together with plate spreading. In addition, deep-seated magma accumulation appears to take place, with its centre ∼ 10 km north of the Krafla caldera. The relative strength of these sources has varied with time. New results from a levelling survey and GPS measurements in 2005 allow an updated view on the deformation field. Deformation rates spanning 2000–2005 are the lowest recorded in the 30-year history of geodetic studies at the volcano. The inferred rate of 2000–2005 subsidence related to processes in the shallow magma chamber is less than 0.3 cm/yr whereas it was ∼ 5 cm/yr in 1989–1992. Currently, the highest rate of subsidence takes place in the Leirbotnar area, within the Krafla caldera, and appears to be a result of geothermal exploitation.     

地壳中岩浆沿已有断层运移的数值模拟——在长白山天池火山的应用

本文利用有限元及有限差分的方法,对壳内岩浆房或岩浆囊中的岩浆在构造应力及由于围岩与岩浆的密度差产生的浮力作用下,沿已有断层向上运移的动力学过程进行了数值模拟.在岩浆囊顶部与上覆岩层接触处,沿着已有微小破裂,岩浆在一定超压力条件下使已有断层张开并继续向上延伸,从而形成岩浆向上运移的通道.研究了岩浆黏度、密度差、模型深度对最小超压力（岩浆运移到地表所需的最小岩浆房超压力）的影响.在10 km深度的地壳中,若岩浆黏度为0.1~10³ Pa·s,当超压力达到17~20 MPa时,岩浆压力可以驱动岩浆运移到地表层;同时,岩浆动力黏度越大,使岩脉运移到地表需要的超压力就越大.当密度差为300~700 kg·m^-3,其变化对超压力的影响比较小.本文亦对比了三维应力条件和二维平面应变条件下不同结果,比较了不同条件下岩浆运移造成的地表垂直位移变化.结合长白山天池火山地区的区域地质环境,对长白山天池火山岩浆运移条件进行了参数试验性计算分析,估算了在给定长白山天池火山模型条件下地下可能存在的岩浆囊的大小,其结果对认识长白山天池火山地区岩浆活动及相关的预测和监控有参考意义.     

Identification of a magma chamber in the Ghoubbet-Asal rift (Djibouti) from a magnetotelluric experiment

A very broad band (10^?3 to 10⁴ Hz) magnetotelluric investigation of the axial zone of the Ghoubbet-Asal rift (Djibouti) has revealed a shallow (2–4 km) magma chamber which can be mapped in some detail. The suggested roof of the chamber is shallowest very close to the Ardoukoba volcano which was built during the November 1978 rifting episode.     

The magmatic feeding system of El Reventador volcano (Sub-Andean zone,Ecuador) constrained by texture,mineralogy and thermobarometry of the 2002 erupted products

After a 26 years long quiescence El Reventador, an active volcano of the rear-arc zone of Ecuador, entered a new eruptive cycle which lasted from 3 November to mid December 2002. The initial sub-Plinian activity (VEI 4 with andesite pyroclastic falls and flows) shifted on 6 and 21 November to an effusive stage characterized by the emission of two lava flows (andesite to low-silica andesite Lava-1 and basaltic andesite Lava-2) containing abundant gabbro cumulates. The erupted products are medium to high-K calc-alkaline and were investigated with respect to major element oxides, mineral chemistry, texture and thermobarometry. Inferred pre-eruptive magmatic processes are dominated by the intrusion of a high-T mafic magma (possibly up to 1165 ± 15 °C) into an andesite reservoir, acting as magma mixing and trigger for the eruption. Before this refilling, the andesite magma chamber was characterized by water content of 5.3 ± 1.0%, high oxygen fugacity (> NNO + 2) and temperatures, in the upper and lower part of the reservoir, of 850 and 952 ± 65 °C respectively. Accurate amphibole-based barometry constrains the magma chamber depth between 8.2 and 11.3 km (± 2.2 km). The 6 October 2002 seismic swarm (hypocenters from 10 to 11 km) preceding El Reventador eruption, supports the intrusion of magmas at these depths. The widespread occurrence of disequilibrium features in most of the andesites (e.g. complex mineral zoning and phase overgrowths) indicates that convective self-mixing have been operating together with fractional crystallization (inferred from the cognate gabbro cumulates) before the injection of the basic magma which then gave rise to basaltic andesite and low-silica andesite hybrid layers. Magma mixing in the shallow chamber is inferred from the anomalous SiO₂–Al₂O₃ whole-rock pattern and strong olivine disequilibria. Both lavas show three types of amphibole breakdown rims mainly due to heating (mixing processes) and/or relatively slow syn-eruptive ascent rate (decompression) of the magmas. The lack of any disequilibrium textures in the pumices of the 3 November fall deposit suggest that pre-eruptive mixing did not occur in the roof zone of the chamber. A model of the subvolcanic feeding system of El Reventador, consistent with the intrusion of a low-Al₂O₃ crystal-rich basic magma into an already self-mixed andesite shallow reservoir, is here proposed. It is also inferred that before entering the shallow chamber the “basaltic” magma underwent a polybaric crystallization at deeper crustal levels.     

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     

夏威夷一维深部电性构造研究

利用加拿大凤凰公司在夏威夷大岛地区布设的长期大地电磁观测台站,选取2018年5月至2019年4月的观测数据,反演并获得了该台站下方大岛深部的电性构造及其随时间变化的特征,并将反演结果与前人探测结果进行对比。研究结果表明：在观测点下方2—3 km处存在一个明显的低阻薄层,电阻率大约为20Ω·m,推测其可能是地表浅部发育的岩浆囊;100 km深处电阻率变化范围在20—70Ω·m,推测该深度下可能为活跃的复杂岩浆储层,反演结果与前人在邻近地区的探测结果具有可比性。     

Magma convection and mixing dynamics as a source of Ultra-Long-Period oscillations

Many volcanic eruptions are shortly preceded by injection of new magma into a pre-existing, shallow (<10 km) magma chamber, causing convection and mixing between the incoming and resident magmas. These processes may trigger dyke propagation and further magma rise, inducing long-term (days to months) volcano deformation, seismic swarms, gravity anomalies, and changes in the composition of volcanic plumes and fumaroles, eventually culminating in an eruption. Although new magma injection into shallow magma chambers can lead to hazardous event, such injection is still not systematically detected and recognized. Here, we present the results of numerical simulations of magma convection and mixing in geometrically complex magmatic systems, and describe the multiparametric dynamics associated with buoyant magma injection. Our results reveal unexpected pressure trends and pressure oscillations in the Ultra-Long-Period (ULP) range of minutes, related to the generation of discrete plumes of rising magma. Very long pressure oscillation wavelengths translate into comparably ULP ground displacements with amplitudes of order 10⁻⁴–10⁻² m. Thus, new magma injection into magma chambers beneath volcanoes can be revealed by ULP ground displacement measured at the surface.     

Stress field around the Coloumbo magma chamber,southern Aegean: Its significance for assessing volcanic and seismic hazard in Santorini

Coloumbo submarine volcano lies 6.5 km offshore the NE part of the Santorini island complex and exhibits high seismicity along with vigorous hydrothermal activity. This study models the local stress field around Coloumbo's magma chamber and investigates its influence on intrusion emplacement and geometry. The two components of the stress field, hoop and radial stress, are calculated using analytical formulas that take into account the depth and radius of the magma chamber as these are determined from seismological and other observations. These calculations indicate that hoop stress at the chamber walls is maximum at an angle of 74° thus favouring flank intrusions, while the radial stress switches from tensile to compressive at a critical distance of 5.7 km from the center of the magma chamber. Such estimates agree well with neotectonic and seismological observations that describe the local/regional stress field in the area. We analyse in detail the case where a flank intrusion reaches the surface very near the NE coast of Thera as this is the worst-case eruption scenario. The geometrical features of such a feeder dyke point to an average volumetric flow rate of 9.93 m³ s⁻¹ which corresponds to a Volcanic Explosivity Index of 3 if a future eruption lasts about 70 days. Hazards associated with such an eruption include ashfall, ballistic ejecta and base surges due to explosive mixing of magma with seawater. Previous studies have shown that areas near erupting vents are also foci of moderate to large earthquakes that precede or accompany an eruption. Our calculations show that a shallow event (3–5 km) of moment magnitude 5.9 near the eruptive vent may cause Peak Ground Acceleration in the range 122–177 cm s⁻² at different locations around Santorini. These values indicate that seismic hazard even due to a moderate earthquake near Coloumbo, is not trivial and may have a significant impact especially on older buildings at Thera island.     

The volcano-magnetic effect

Summary A geometrically simple volcano is considered, havig a spherical magma chamber of 2.5 km radius centred at 10 km depth. The Curie point isotherm is assumed to be a plane at 20 km depth, except for the spherical volume which is also non-magnetic. The stress pattern in the vicinity of the spherical chamber, due to regional stress of sufficient intensity to cause an eruptions, is used to calculate the change in magnetization which results from the piezomagnetic effect through the volume of solid rock. The consequent magnetic field anomaly at the surface is then obtaied by numerical integration of the dipole law of force over the stressed volume. For rocks of the type found on the volcanic island of St. Vincent (West Indies), this model gives a maximum local volcano-magnet c effect of about 7 gammas.     

2-D crustal Poisson's ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson's ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson's ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson's ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson's ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson's ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson's ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region.     

Tectonic and structural evolution of Gough volcano: A volcanological model

Morphostructural, stratigraphic and tectonic data indicate that the evolution of Gough volcano is similar to other oceanic intraplate volcanoes, is older than 1 Ma, and is related to a transform fault. At least six evolutionary stages can be distinguished within two major magmatostructural periods dominated by basaltic and trachytic magmas, respectively.The basaltic shield volcano is characterized by a curved, elongated shape in plan and a rift zone with a high density of dykes, combined with a radial intrusive system. The latter is interpreted as being fed by a magma chamber some 4 km below the surface. The activity of the volcano became more centralized at the end of the basaltic period and its slopes became steeper. This corresponds to the development of a shallower and narrower central conduit in the edifice. The basaltic period was terminated by formation of a shield caldera related to the 4 km deep magma chamber. The term “shield caldera” is used for a collapse structure that is postmagmatic, large in comparison with the diameter of the volcano, and delimited by normal faults that do not show a closed circular pattern but rather a series of arcs. In contrast, summit calderas are defined as smaller, circular-shaped, centrally situated, synmagmatic features, related to a central shallow column. During the basaltic period, landslides were generated on the flanks of the edifice as a result of slope stability factors which are not easy to determine at present, and dynamic factors among which the intrusion of magma along a curved zone certainly played a major role.The trachytic period is characterized by comparatively rare pyroclastic deposits and a large volume of thick flows extruded from domes. These extrusions, as well as plugs, formed from vertical cylindrical columns of magma rising from shallow individual magma pockets fed by the main reservoir.     

Net gravity decrease at Askja volcano,Iceland: constraints on processes responsible for continuous caldera deflation, 1988–2003

Askja caldera in northeast Iceland has been in a state of unrest for decades. Ground-deformation surveys show that the rate of deformation, i.e., deflation, is much higher then observed at any other dormant volcano in Iceland. This work presents the results from microgravity and deformation studies at Askja from 1988 to 2003. The deflation reaches a maximum of −0.46 m in the centre of the caldera, relative to a station outside the caldera, during the study period. The source of deformation is inferred to be at ∼3 km depth and a recent study infers a second deeper source at ∼16 km depth. The deflation is consistent with a subsurface volume change of −0.018 km³. We find a net microgravity decrease of 115 μGal in the centre of the caldera relative to the same station. This corresponds to a subsurface mass decrease of 1.6×10¹¹ kg between 1988 and 2003 based on the use of a point source model. A combination of magma drainage and cooling and contraction of the shallow magma reservoir at 3 km depth is our favoured model, consistent with the integrated observations. We suggest that extensional tectonic forces generate space in the ductile part of the crust to accommodate ongoing magma drainage from the shallow magma chamber.