sizing up crystals and their melt inclusions a new approach to crystallization studies a correlation between melt inclusion composition and host crystal size has been discovered in a basaltic ash deposit from a volcanic arc setting least evolved melt inclusions featuring high h2o co2 saturation pressures are found in large olivine crystals melt inclusion compositions become more evolved and volatile saturation pressures decrease as host crystal size decreases this trend is interpreted as reflecting crystal nucleation and growth during decompression of a magma saturated with h2o and co2 in general melt inclusion compositions show large incompatible element enrichment and modest major element variation consistent with in situ crystallization crystal size is used to temporally constrain the formation age of melt inclusions and the decompression path of the magma an estimated crystal growth rate of 10 8 cm s 1 yields an initial average magma ascent rate of 4 8 km yr 1 at mid crustal depths at shallower depths the ascent rate was roughly an order of magnitude greater

seismic and geodetic progression of the 2018 summit caldera collapse of kīlauea volcano abstract the 2018 eruption of k lauea volcano hawai i resulted in a major collapse of the summit caldera along with an effusive eruption in the lower east rift zone the caldera collapse comprised 62 highly similar collapse cycles of strong ground deformation and earthquake swarms that ended with a magnitude 5 collapse event and one partial cycle that did not end with a collapse event we analyzed geodetic and seismic data to better understand how the caldera collapse progressed over 3 months of activity focusing on the cyclical activity we identified 3 main phases of collapse initial ring fault activation and small explosions phase 1 an eastward shift in activity and freeing of the central piston phase 2 and a recoupling of the piston to the reservoir followed by relatively steady behavior until the eruption s end phase 3 additionally we observed geodetic evidence of tangential motion from the localization of the main ring fault phase 2 and the formation of a major peripheral ring fault on the eastern side of the collapse caldera during phase 3 both geodetic and seismic parameters suggest that the collapse may have had an eastward component of motion after the ring fault system had formed the cyclical seismic and geodetic parameters show no obvious signs that the collapse was coming to an end with the only notable change being a significant increase in the ratio of cyclical displacement to co collapse displacement observed during the last complete cycle on gnss stations outside the caldera region

basaltic calderas collapse dynamics edifice deformation and variations of magma withdrawal 1 the incremental caldera collapses of fernandina 1968 miyakejima 2000 and piton de la fournaise 2007 are analyzed in order to understand the collapse dynamics in basaltic setting and the associated edifice deformation for each caldera the collapse dynamics is assessed through the evolution of the 1 time interval t between two successive collapse increments 2 amount of vertical displacement during each collapse increment and 3 magma outflow rate during the whole collapse caldera process we show from the evolution of t that piton de la fournaise and fernandina were characterized by a similar collapse dynamics despite large differences in the caldera geometry and the duration of the whole collapse caldera process this evolution significantly differs from that of miyakejima where t strongly fluctuated throughout the whole collapse process quantification of the piston vertical displacements enables us to determine the magma outflow rates between each collapse increment displacement data tiltmeter and or gps for piton de la fournaise and miyakejima are used to constrain the edifice overall deformation and the edifice deformation rates these data reveal that both volcanoes experienced edifice inflation once the piston collapsed into the magma chamber such a deformation which lasts during the first collapse increments only is interpreted as the result of larger volume of piston intruded in the magma chamber than magma withdrawn before each collapse increment once the effect of the collapsing rock column vanishes edifice deflates we also determine for each caldera the critical amount of magma evacuated before collapse initiation and compare it to analog models the significant differences between models and nature are explained by the occurrence of preexisting weak zones in nature i e the ring faults that are not taken into account in analog models finally we show that t at piton de la fournaise and fernandina was equally controlled by the frictional resistance along the ring faults and the magma outflow rate in addition to these two parameters the collapse dynamics of miyakejima was also influenced by variations of the magma bulk modulus which changed after the influx of deep gas rich magma into the collapse related magma chamber altogether our results show that the dynamics of caldera collapse in basaltic volcanoes proceeds in two phases phase 1 starting with the first collapse is characterized by the largest collapse amplitude an incremental edifice inflation and a step by step increase of the rate of magma outflow phase 2 shows a rapid decrease of the magma discharge rate to a low level concomitant with the continuous edifice deflation if deep magma is injected into the magma chamber as at miyakejima an additional phase occurs phase 3

constraining spatiotemporal characteristics of magma migration at piton de la fournaise volcano from pre eruptive seismicity volcano tectonic seismicity has been recorded for decades on various volcanoes and is linked with the magma transport and reservoir pressurization yet earthquakes often appear broadly distributed such that magma movement is difficult to infer from its analysis we explore the seismicity that occurred before eruptions at piton de la fournaise in the last 5 years using template matching and relocation techniques we produce a refined image of the summit seismicity demonstrating that most earthquakes are located on a ring structure however only a portion of this structure is activated before each eruption which provides an indication as to the direction of the future eruptive site furthermore we show that the delay between the pre eruptive swarm and the eruption onset is proportional to the distance of the eruptive fissures relative to the summit cone this reveals that the beginning of the intrusion already bears information regarding the future eruption location

intermittent inflations recorded by broadband seismometers prior to caldera formation at miyake jima volcano in 2000 abstract very long period vlp pulses with widths of 20 s on velocity seismograms were observed during volcanic activity at miyake jima volcano japan in 2000 the vlp events occurred repeatedly during a few days prior to caldera formation and essentially vanished following the onset of caldera collapse waveform inversions of the pulse like signals point to a source offset 3 5 km beneath and 1 km south of the summit a candidate for the source mechanism is the inflation of an elliptical cylinder with axis tilted 20 30 from vertical and major axis of the elliptical cross section oriented northeast southwest the inferred mechanism appears consistent with a step like pressurization of a magma reservoir impacted by a falling rock mass in response to gravitational instability the repeated occurrences of the rock collapses lead to the caldera formation at miyake jima

three dimensional shear velocity anisotropic model of piton de la fournaise volcano la reunion island from ambient seismic noise we cross correlate 4 years of seismic noise from the seismic network of piton de la fournaise volcano la reunion island to measure the group velocity dispersion curves of rayleigh and love waves we average measurements from vertical and radial components to obtain 577 rayleigh wave dispersion curves the transverse components provided 395 love wave dispersion curves we regionalize the group velocities measurements into 2 d velocity maps between 0 4 and 8 s finally we locally inverted these maps for a pseudo 3 d anisotropic shear velocity model down to 3 km below the sea level using a neighborhood algorithm the 3 d isotropic shear wave model shows three distinct high velocity anomalies surrounded by a low velocity ring the anomaly located below the present plaine des sables could be related to an old intrusive body at the location of the former volcanic center before it migrated toward its present location the second high velocity body located below the summit of the volcano likely corresponds to the actual preferential dyke intrusion zone as highlighted by the seismicity the third high velocity anomaly located below the grandes pentes and the grand brule areas and is an imprint of the solidified magma chamber of the dismantled les alizes volcano radial anisotropy shows two main anomalies positive anisotropy above sea level highlighting the recent edifice of piton de la fournaise with an accumulation of horizontal lava flows and the second one below the sea level with a negative anisotropy corresponding to the ancient edifice of piton de la fournaise dominated by intrusions of vertical dykes

the collapse of barðarbunga caldera iceland lying below vatnajokull ice cap in iceland bar arbunga stratovolcano began experiencing wholesale caldera collapse in 2014 august 16 one of the largest such events recorded in the modern instrumental era simultaneous with this collapse is the initiation of a plate boundary rifting episode north of the caldera observations using the international constellation of radar satellites indicate rapid 50 cm d 1 subsidence of the glacier surface overlying the collapsing caldera and metre scale crustal deformation in the active rift zone anomalous earthquakes around the rim of the caldera with highly nondouble couple focal mechanisms provide a mechanical link to the dynamics of the collapsing magma chamber a model of the collapse consistent with available geodetic and seismic observations suggests that the majority of the observed subsidence occurs aseismically via a deflating sill like magma chamber

using centroid time delays to characterize source durations and identify earthquakes with unique characteristics the relationship between m 0 and the rupture duration is often difficult to establish this is particularly true for large earthquakes for which the moment rate functions mrf generally have complicated shapes and the estimated durations can vary considerably depending on the methodology used to evaluate the mrf in this work we show that the centroid time delay c provides an alternative estimate of the source duration inverted mrfs often end gradually making the end of coseismic rupture difficult to detect in such cases when the rupture duration is not well defined the time delay c is a useful quantity to represent the first order temporal characteristics of the rupture process variations in stress parameter can be investigated by assuming a standard scaling relationship between the seismic moment m0m0 and c this simple scaling relationship can also be used to identify unusual earthquakes with unique source properties such as events involving complicated rupture processes or earthquakes characterized by unusual rupture velocities stress drops or aspect ratios

seismic source mechanism of degassing bursts at kilauea volcano hawaii results from waveform inversion in the 10 50 s band the current march 2008 to february 2009 summit eruptive activity at kilauea volcano is characterized by explosive degassing bursts accompanied by very long period vlp seismic signals we model the source mechanisms of vlp signals in the 10 50 s band using data recorded for 15 bursts with a 10 station broadband network deployed in the summit caldera to determine the source centroid location and source mechanism we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account the vlp signals associated with the bursts originate in a source region 1 km below the eastern perimeter of halemaumau pit crater the observed waveforms are well explained by the combination of a volumetric component and a vertical single force component for the volumetric component several source geometries are obtained which equally explain the observed waveforms these geometries include 1 a pipe dipping 64 to the northeast 2 two intersecting cracks including an east striking crack dike dipping 80 to the north intersecting a north striking crack another dike dipping 65 to the east 3 a pipe dipping 58 to the northeast intersecting a crack dipping 48 to the west southwest and 4 a pipe dipping 57 to the northeast intersecting a pipe dipping 58 to the west southwest using the dual crack model as reference the largest volume change obtained among the 15 bursts is 24 400 m3 and the maximum amplitude peak to peak of the force is 20 gn each burst is marked by a similar sequence of deflation and inflation trailed by decaying oscillations of the volumetric source the vertical force is initially upward synchronous with source deflation then downward synchronous with source reinflation followed by oscillations with polarity opposite to the volumetric oscillations this combination of force and volume change is attributed to pressure and momentum changes induced during a fluid dynamic source mechanism involving the ascent expansion and burst of a large slug of gas within the upper 150 m of the magma conduit as the slug expands upon approach to the surface and more liquid becomes wall supported by viscous shear forces the pressure below the slug decreases inducing conduit deflation and an upward force on the earth the final rapid slug expansion and burst stimulate vlp and lp oscillations of the conduit system which slowly decay due to viscous dissipation and elastic radiation consideration of the fluid dynamic arguments leads us to prefer the dual crack vlp source model as it is the only candidate model capable of producing plausible values of length scales and pressure changes the magnitudes of the vertical forces observed in the 15 bursts appear consistent with slug masses of 104 to 106 kg

crystal size distribution csd in rocks and the kinetics and dynamics of crystallization crystal size distributions csds measured in metamorphic rocks yield quantitative information about crystal nucleation and growth rates growth times and the degree of overstepping t of reactions during metamorphism csds are described through use of a population density function n dn dl where n is the cumulative number of crystals per unit volume and l is a linear crystal size plots of ln n vs l for olivine pyroxene and magnetite in high temperature 1000 c basalt hornfelses from the isle of skye define linear arrays indicating continuous nucleation and growth of crystals during metamorphism using the slope and intercept of these linear plots in conjunction with growth rate estimates we infer minimum mineral growth times of less than 100 years at t 10 c and nucleation rates between 10 4 and 10 1 cm3 s garnet and magnetite in regionally metamorphosed pelitic schists from south central maine have csds which are bell shaped we interpret this form to be the result of two processes 1 initial continuous nucleation and growth of crystals and 2 later loss of small crystals due to annealing the large crystals in regional metamorphic rocks retain the original size frequency distribution and may be used to obtain quantitative information on the original conditions of crystal nucleation and growth the extent of annealing increases with increasing metamorphic grade and could be used to estimate the duration of annealing conditions if the value of a rate constant were known finally the different forms of crystal size distributions directly reflect differences in the thermal histories of regional vs contact metamorphosed rocks because contact metamorphism involves high temperatures for short durations resulting csds are linear and unaffected by annealing similar to those produced by crystallization from a melt because regional metamorphism involves prolonged cooling from high temperatures primary linear csds are later modified by annealing to bell shapes

edifice strength and magma transfer modulation at piton de la fournaise volcano 1 from 2003 to 2007 eruptive activity at piton de la fournaise was shown to follow cycles comprising many summit proximal eruptions and finishing by a distal eruption gps measurements evidenced striking asymmetric deformation between its western and eastern flanks horizontal displacements recorded during interdistal periods showed a characteristic amplitude at the top of the eastern flank displacements recorded at the base of the summit cone showed a bimodal distribution with low amplitudes during interdistal periods and large ones during distal eruptions to account for displacement asymmetry characteristic amplitude and large flank displacement we modeled the volcanic edifice using a drucker prager elastoplastic rheology friction angles of 15 and 30 were needed to model the displacements respectively during distal eruptions and interdistal periods this change shows that strain weakening occurred during distal events large plastic displacement that occurred in the eastern flank during distal eruptions relaxed the horizontal elastic stress accumulated during interdistal periods it triggered summit deflation horizontal magma transfer and distal flank eruption and reset the eruptive cycle our elastoplastic models also show that simple source geometries may induce large eastern flank displacements that would be explained by a complex geometry in a linear elastic edifice magma supply is often thought to control volcano s eruptive activity with surface deformation reflecting changes in magma supply rate the volcano s response being linear our results bring some evidences that on piton de la fournaise time space discretization of magma transfer may be the result of the edifice s nonlinear response rather than changes in magma supply

olivine growth rates in a tholeiitic basalt an experimental study of melt inclusions in plagioclase abstract basaltic glass inclusions trapped in plagioclase phenocrysts an 84 are remnant of their parent magmatic liquid they can be used as natural reactors for the experimental investigation of olivine growth rate as a function of temperature the growth of one olivine nucleus can be observed at constant temperature supercooling from 15 to 150 c have been investigated growth habits vary from equant to feather in qualitative agreement with previous studies growth rates vary from 10 m s 1 to 6 10 7 m s 1 they vary with the direction the growth process planar or dendritic and the degree of supercooling chemical analysis of crystal overgrowth and the remaining liquid composition enables a mass balance calculation which confirms the rates determined optically the small number of results obtained so far does not permit to draw inferences on the growth mechanisms

characterization of westward propagating signals in the south atlantic from altimeter and radiometer records abstract radar altimeter data from topex poseidon and jason 1 and microwave radiometer data from trmm tmi are used to investigate the large scale variability between 10 5 s and 35 5 s in the south atlantic ocean the proposed method for the analysis of the longitude time diagrams of the cross correlation between ssh and sst anomalies shows that the variability in mid latitudes is a blend of first mode baroclinic rossby waves and propagating mesoscale eddy like structures the estimated phase speed of the wave cp and propagation speed of the eddies cv are similar in 70 of the cases the absolute difference between cp and cv is less than 11 in 40 of the cases this difference is less than 5 statistical results indicate that in the case of eddies as the thermocline deepens the sea surface temperature rises and vice versa however planetary waves show more complex yet self consistent results in lower latitudes 10 5 s 15 5 s the shallower thermocline and the weak thermal gradients impose a zero phase lag between temperature and height similar to eddies poleward of those latitudes sea surface temperature and height are in quadrature of phase this indicates that geostrophic advection of the relatively stronger thermal gradient is performed by rossby waves

size and volume evaluation of the caldera collapse on piton de la fournaise volcano during the april 2007 eruption using aster stereo imagery 1 the scale of summit collapse during the april 2007 eruption of piton de la fournaise volcano reunion island western indian ocean was evaluated using before and after differential digital elevation models dems derived from nadir and backward looking images from the advanced spaceborne emission and reflection radiometer aster instrument the dimensions of horizontal length and width volume and depth of the depression were estimated as 1 100 800 m 9 6 10 7 m 3 and 320 m respectively these estimates show excellent agreement with field survey data a ring shaped thermal anomaly diameter of about 500 m was found at dolomieu crater during or just after the collapse on nighttime aster swir image it may correspond to the high temperature areas associated with the cross section of a hydrothermal zone at a constant altitude inside the crater our analysis of aster orbital data documents topographic and related thermal changes as a result of the 2007 eruption and demonstrates the power of aster as a volcanological tool citation urai m n geshi and t staudacher 2007 size and volume evaluation of the caldera collapse on piton de la fournaise volcano during the april 2007 eruption using aster stereo imagery geophys

centroid single force inversion of seismic waves generated by landslides the centroid moment tensor cmt waveform inversion method of dziewonski et al 1981 is modified to analyze long period seismic waves generated by single force events such as the gigantic landslides associated with the 1980 eruption of mount st helens we refer to the method as centroid single force csf inversion as the result of the inversion we obtain the spatio temporal centroid and three components of a vector we call the csf vector which is the seismic point source representation of a landslide and whose direction is the opposite of the direction of sliding the sealer magnitude of the csf vector measures the overall size of the landslide and is given by mcsf md where m is the mass of the sliding object and d is the sliding distance we apply csf inversion to the long period surface waves generated by the 1980 mount st helens event the 1975 kalapana earthquake and the 1974 mantaro landslide using seismograms from the gdsn ida and hglp networks the csf solution for the st helens event is remarkably consistent with the actual geological observations and proves the efficacy of csf inversion the csf solution for the kalapana event does not fit the overall waveforms better than the cmt solution does however because the cmt solution does not explain the observed love wave radiation pattern it may be necessary to invoke a combination of both types of mechanism for this event although the geometry of the csf solution the direction of the force for the mantaro event is consistent with the actual landslide mcsf is about 5 times smaller than that expected from the mass and travel distance of the landslide estimated by geological observations this discrepancy may suggest either the relatively aseismic nature of this landslide or that the total volume of the slide was overestimated by analogy to single station cmt inversion single station csf inversion also appears to be stable and useful

gravity structure of piton de la fournaise volcano and inferred mass transfer during the 2007 crisis abstract the structure and mass transfers at piton de la fournaise volcano are investigated using new on land and offshore measurements the structures are defined using 2 and 3d modelling short wavelength positive anomalies to the w and nw of the central area are interpreted as shallow piles of thick dense lava flows which have filled volcano tectonic depressions several hundreds of meters deep this finding provides a new insight into the evolution of the volcano conversely the currently active central zone is associated with a negative short wavelength anomaly explained by the presence of thin highly vesiculated and fractured lava flows building the central cone low density formations underlying the rivieres des remparts and langevin areas are tentatively interpreted as breccias related to erosion and or to southward directed mass wasting topographic offshore highs in the continuation of the se and ne rift zones of piton de la fournaise are characterized by negative gravity anomalies and reversed magnetic anomalies their interpretation as hyaloclastites and pillow lavas provides a coherent explanation for both types of signals if we admit that they were primarily constructed during the alizes stage two main deeper dense structures are also found the grand brule complex and a complex beneath the plaine des sables and the enclos fouque both are identified as hypovolcanic intrusive complexes the grand brule complex 10 18 km in plan view shows a marked ns elongated geometry unlike the other comparable complexes on the island the plaine des sables enclos fouque complex a more or less circular structure 5 km in diameter is attributed to the ancient shield stage it is centered on the enclos fault suggesting that the paleo center of the ancient shield was located about 1 5 km eastward of the previously proposed location an apparent contradiction was noted between the seismic and gravity images of the internal structure in particular a high velocity cylinder beneath the summit was not detected in the gravity models this contradiction was resolved by considering the relative sensitivity of each method finally the gravity dataset has allowed us to explore the gravity changes that accompanied the major eruptive and volcano tectonic crisis of april 2007 350 m collapse of the summit and emission of 120 106 m3 of lava the change in the mass distribution of the edifice was large enough to be detected in the survey data collected in the central area before and after the crisis however the accuracy was too low for the precise assessment of the mass transfers

analysis of long period seismic waves excited by the may 18 1980 eruption of mount st helens a terrestrial monopole long period 100 to 260 s love and rayleigh waves excited by the eruption of mount st helens on may 18 1980 and recorded by ida sro and asro stations were analyzed to determine the mechanism of the eruption the amplitude radiation patterns of both rayleigh and love waves are two lobed with a nodal direction in e5 s for rayleigh waves and in n5 e for love waves these radiation patterns preclude any double couple mechanism the radiation pattern the initial phase the relatively large amplitude ratio of love to rayleigh waves and the existence of clear nodes in the radiation patterns of fundamental mode and higher mode rayleigh waves suggest that the source is represented by an almost horizontal less than 15 from the horizontal single force pointed toward s5 w the surface wave spectra fall off very rapidly at periods shorter than 75 s suggesting a very slow source process although the details of the source time history could not be determined a smooth bell shaped time function f o s t 1 2 f o l cos t for 0 t 2 and f o s t 0 for t 2 with 75 s is considered appropriate on the basis of comparison between synthetic and observed seismograms and of the shape of the source spectrum the peak value of the force f 0 is about 10 18 dynes the tailing end of the source time function could not be resolved and some overshoot may be added the magnitude and the time history of the force can be explained by a northward landslide followed by a lateral blast observed at the time of the eruption two distinct events about 110 s apart can be identified on body wave and short period surface wave records the first event may correspond to the earthquake which triggered the landslide and the lateral blast the second event appears to correspond to a second large earthquake and explosion which took place about 2 minutes after the first earthquake

cyclic magma storages and transfers at piton de la fournaise volcano la reunion hotspot inferred from deformation and geochemical data abstract piton de la fournaise is in a period of intense volcanic activity since 1998 to constrain the magma dynamics responsible for this activity we combined gps ground deformation monitoring interpreted through numerical modelling and geochemistry two cycles of continuous volcano inflation are evidenced for the may 2004 december 2005 period with a rest from march to october 2005 these inflations are consistent with two cycles of compatible major element enrichment in the emitted lavas numerical models indicate that the pressurization of a single magma reservoir may be responsible for the observed pre eruptive inflations of the volcano the reservoir located at 2300 m depth has a radius of 500 m at the beginning of each cycle dykes propagate from the roof of the reservoir and yield eruptions of differentiated basalt near the summit at the end of the cycle dykes propagate from the eastern sidewall of the reservoir and yield distal eruptions of primitive magmas away from the summit the volumes of magma emitted during the primitive eruptions seem too large to explain the surface deformations and therefore suggest some refill of the reservoir by deeper magmas our results may be used to predict the location and lava volume of future eruptions at piton de la fournaise volcano depending on the timing of these eruptions within a cycle of volcanic activity

relation of summit deformation to east rift zone eruptions on kilauea volcano hawaii an inverse relationship exists between the summit deflation of kilauea as recorded by summit tilt and the elevation of associated eruptive vents on the east rift zone this relationship implies that east rift eruptions drain the summit magma reservior to pressure levels that are dependent on the elevation of the eruptive vents

a refined technique for measuring crystal size distributions in thin section log normal size distributions are of the form n noe l where n number density l crystal length and is a constant a method for measuring three dimensional log normal crystal or grain size distributions csds in thin section has been deduced from computer experiments in which 2d sections were cut through assemblages of 3d solids the size ranges and distributions studied were appropriate for igneous microphenocryst to megacryst populations conversion from 2d to 3d is based on an exact correction for spheres of uniform diameter cumulate numbers of polygons with length l n2d are converted to n3d by the equation ln n3d ln n2d l s ln l s the number density is then obtained as n dn dl the parameters s and correct the measured lengths and no no number density at l 0 respectively and are functions of crystal shape the parameter is a weak function of the degree of spatial orientation of the crystals highly symmetrical shapes such as cubes octahedra and elongated prisms can be accurately measured when randomly oriented however rectangular solids with a b c cannot be accurately measured because they produce bimodal length distributions in cross section strongly oriented textures trachytic or lineated can be accurately measured regardless of crystal shape new csd data from alkaline rocks and a kimberlite give examples of csds modified by megacryst retention xenocryst addition phenocryst accumulation and groundmass nucleation

the april 2007 eruption and the dolomieu crater collapse two major events at piton de la fournaise la reunion island indian ocean abstract after 10 years of high activity and 24 eruptions piton de la fournaise pdf produced on april 2 to may 1 2007 one of its most voluminous and intense eruption since at least one century the eruption focused at 590 m elevation in the grand brule on the south east flank of piton de la fournaise 7 km away from summit it was located close to the southern cliff of the grand brule volcano tectonic collapse the rempart du tremblet and to the nearby village called le tremblet and it highly rattled and incommoded the inhabitants there eleven hours after the beginning of the eruption the lava flows reached the sea and created 0 45 km 2 of new land until the end of eruption in the night of april 5 the rock column beneath dolomieu crater started to collapse with a piston like mechanism into the magma chamber forming within less than 24 h the most recent and well documented caldera we report in this paper the chronology of these two related major events in the recent activity of pdf seen by the scientific networks of the volcanological observatory ovpf ipgp and field observations which allowed to back up the scenario of the caldera formation and to quantify these two exceptional episodes

very long period seismic signals and caldera formation at miyake island japan over a period of roughly 40 days starting on 8 july 2000 a caldera structure 1 7 kilometers in diameter developed by means of gradual depression and expansion of the summit crater at miyake island japan at the same time very long period vlp seismic signals were observed once or twice a day source mechanism analyses of the vlp signals show that the moment tensor solutions are smooth step functions over a time scale of 50 seconds with dominant volumetric change components we developed a model to explain the caldera and the vlp signals in which a vertical piston of solid materials in the conduit is intermittently sucked into the magma chamber by lateral magma outflow this model offers potential for making quantitative estimations of the characteristic physical properties of magma systems

magma dynamics and collapse mechanisms during four historic caldera forming events 1 four historic caldera forming events were studied to understand the relationship of magma withdrawal processes and caldera subsidence mechanisms two calderas are silicic katmai in 1912 and pinatubo in 1991 and two are basaltic fernandina in 1968 and miyakejima in 2000 all events have sufficient geophysical geologic and petrologic data with which to examine and model magma withdrawal and caldera collapse the data reveal that the magmas erupted at katmai and pinatubo were in a bubbly state in the reservoir immediately before and during caldera collapse the bubbly magma allowed for its efficient extraction from the reservoir causing significant underpressures to develop rapidly particularly in the case of katmai where the erupted rhyolite was voluminous and nearly aphyric and has very low viscosity the rapidly developing underpressures at katmai and pinatubo caused sudden en masse caldera collapse halfway through the climactic eruptions thereby liberating large amounts of seismic energy at fernandina and miyakejima by contrast caldera collapse was initiated early and continued for an extended period of time from weeks to months consisting of a series of discrete subsidence events manifested by large earthquakes at fernandina and by very long period vlp signals at miyakejima systematic changes in earthquake magnitudes and quiescent intervals at both volcanoes reveal changes in friction as collapse took place during extended time intervals

source mechanism of long period events at kusatsu shirane volcano japan inferred from waveform inversion of the effective excitation functions abstract we investigate the source mechanism of long period lp events observed at kusatsu shirane volcano japan based on waveform inversions of their effective excitation functions the effective excitation function which represents the apparent excitation observed at individual receivers is estimated by applying an autoregressive filter to the lp waveform assuming a point source we apply this method to seven lp events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1 3 hz the results of the waveform inversions show dominant volumetric change components accompanied by single force components common to all the events analyzed and suggesting a repeated activation of a sub horizontal crack located 300 m beneath the summit crater lakes based on these results we propose a model of the source process of lp seismicity in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack the rapid discharge of fluid causes the collapse of the fluid filled crack and excites acoustic oscillations of the crack which produce the characteristic waveforms observed in the lp events the presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid filled conduit linking the crack and the base of the crater lake

simple scaling of catastrophic landslide dynamics catastrophic landslides involve the acceleration and deceleration of millions of tons of rock and debris in response to the forces of gravity and dissipation their unpredictability and frequent location in remote areas have made observations of their dynamics rare through real time detection and inverse modeling of teleseismic data we show that landslide dynamics are primarily determined by the length scale of the source mass when combined with geometric constraints from satellite imagery the seismically determined landslide force histories yield estimates of landslide duration momenta potential energy loss mass and runout trajectory measurements of these dynamical properties for 29 teleseismogenic landslides are consistent with a simple acceleration model in which height drop and rupture depth scale with the length of the failing slope

crystal size distribution systematics and the determination of magma storage times the 1959 eruption of kilauea volcano hawaii abstract the theory of crystal size distribution csd theory is based on a steady state population balance that monitors the flux of crystals growing into and out of specificed size categories in precipitating solutions the conservation equation describing this balance permits crystal growth and nucleation rates to be determined directly from crystal size distribution data in this investigation it is shown that csd analysis can alternatively be used to calculate the residence time of crystals in the system or in volcanologic terms the magma storage time prior to eruption if the characteristic crystal growth rate can be independently determined or estimated the crystal size data needed for storage time determinations are easily obtained from thin sections of glassy eruption samples and the mathematical calculations are relatively simple analysis of the errors inherent in this new technique predicts storage estimates accurate to within an order of magnitude

field observations of the 2008 summit eruption at piton de la fournaise ile de la reunion and implications for the 2007 dolomieu collapse abstract more than one year after the most important eruption ever observed at piton de la fournaise in april 2007 a new eruption started in 2008 at piton de la fournaise with three different eruptive phases in september november and december they were located within the 340 m deep dolomieu crater due to the steep dip angle 45 80 and unstable walls of the new crater formed in april 5 2007 no measurements in situ of the 2008 eruption had been possible only observations from the dolomieu crater rim and from helicopter were performed using an infrared camera a portable rangefinder and numerous photographs these field observations allowed precise monitoring of the setting up of the lava flows with time the total erupted volume of lava of the three phases was 2 2 106 m3 and the average flow rates ranged between 0 3 and 1 m3 s 1 lava temperatures of up to 1150 c have been measured by an infrared camera overall infrared images of the dolomieu crater illustrate the control of the eruptive vents by the structure of the april 2007 dolomieu collapse

crystal size distribution csd in rocks and the kinetics and dynamics of crystallization ii makaopuhi lava lake crystal size distribution csd theory has been applied to drill core samples from makaopuhi lava lake kilauea volcano hawaii plagioclase and fe ti oxide size distribution spectra were measured and population densities n were calculated and analyzed using a steady state crystal population balance equation n n0 exp l g slopes on ln n versus crystal size l plots determine the parameter g a product of average crystal growth rate g and average crystal growth time the intercept is j g where j is nucleation rate known temperature depth distributions for the lava lake provide an estimate of effective growth time allowing nucleation and growth rates to be determined that are independent of any kinetic model plagioclase growth rates decrease with increasing crystallinity 9 9 5 4 10 11 cm s as do plagioclase nucleation rates 33 9 1 6 10 3 cm3 s ilmenite growth and nucleation rates also decrease with increasing crystallinity 4 9 3 4 10 10 cm s and 15 2 2 10 3 cm3 s respectively magnetite growth and nucleation rates are also estimated from the one sample collected below the magnetite liquidus g 2 9 10 10 cm s j 7 6 10 2 cm3 s moments of the population density function were used to examine the change in crystallization rates with time preliminary results suggest that total crystal volume increases approximately linearly with time after 50 crystallization a more complete set of samples is needed for material with 50 crystals to define the entire crystallization history comparisons of calculated crystallization rates with experimental data suggests that crystallization in the lava lake occurred at very small values of undercooling this interpretation is consistent with proposed thermal models of magmatic cooling where heat loss is balanced by latent heat production to maintain equilibrium cooling

the 2015 gorkha earthquake a large event illuminating the main himalayan thrust fault the 2015 gorkha earthquake sequence provides an outstanding opportunity to better characterize the geometry of the main himalayan thrust mht to overcome limitations due to unaccounted lateral heterogeneities we perform centroid moment tensor inversions in a 3 d earth model for the main shock and largest aftershocks in parallel we recompute s top and p to s receiver functions from the hi climb data set inverted centroid locations fall within a low velocity zone at 10 15 km depth and corresponding to the subhorizontal portion of the mht that ruptured during the gorkha earthquake north of the main shock hypocenter receiver functions indicate a north dipping feature that likely corresponds to the midcrustal ramp connecting the flat portion to the deep part of the mht our analysis of the main shock indicates that long period energy emanated updip of high frequency radiation sources previously inferred this frequency dependent rupture process might be explained by different factors such as fault geometry and the presence of fluids

april 2007 collapse of piton de la fournaise a new example of caldera formation collapse calderas are frequent in the evolution of volcanic systems but very few have formed during historical times piton de la fournaise is one of the world s most active basaltic shield volcanoes the caldera collapse which occurred during the april 2007 lateral eruption is one of the few large documented collapse events on this volcano it helps to understand the mode and origin of caldera collapses in basaltic volcanoes field observations gps and seismic data show that the collapse occurred at an early stage of the eruption the cyclic seismic signal suggests a step by step collapse that directly influenced the lateral eruption rate likely the caldera results from the combined effect of i the progressive collapse of the plumbing system above the magma chamber since 2000 and ii the large amount of magma withdrawal during the early stage of the eruption by both a significant intrusion within the edifice and an important emission rate

explosive activity of the summit cone of piton de la fournaise volcano la reunion island a historical and geological review summit explosive activity and collapses that form pit craters and calderas represent major volcanic hazards on a dominantly effusive frequently active volcano like piton de la fournaise only three summit collapse events 1986 2002 2007 have been recorded since the foundation of the piton de la fournaise volcano observatory ovpf in 1979 and two of them 1986 and 2007 were associated with weak phreatic activity at piton de la fournaise the normal explosive activity consists of short lived and mild 20 m high lava fountains which quickly evolve into strombolian activity during the eruptions based on comprehensive literature review and high resolution image analysis of surface outcrops and summit caldera walls we reconstructed the time distribution of recent explosive events phreatomagmatic phreatic and their link with summit collapses and lateral flank effusive eruptions in historical time post 1640 ce we recognise two main clusters of explosive events frequent and violent phreatomagmatic to phreatic explosions occurred during the oldest cluster 1708 1878 and alternated with long lasting periods years to decades of summit effusive activity in contrast scarce and on average weak explosions occurred during the youngest cluster 1897 2012 when discrete and short lived 6 months effusive eruptions represent the main eruptive dynamics historical summit collapses pit craters and caldera all localised at the top of the summit cone were related to voluminous lateral eruptions and were followed by a significant decrease in eruptive rate however magma draining during lateral eruptions was not systematically associated with summit collapses or explosions the long lasting occurrence of magma at very shallow depth below the volcano summit followed by a rapid lateral drainage apparently represents a critical condition favouring magma groundwater interaction to produce explosive activity the prehistoric growth of the piton de la fournaise summit cone results chiefly from long lasting to continuous activity centred below its western side bory crater containing lava lakes high lava fountains long lasting effusive activity lava lakes ash plumes and block ejections were common types of eruptive dynamics in the historical past between 1640 and 1878 in this perspective short lived small volume eruptions and long pauses up to six years during the last century of activity of piton de la fournaise can be considered as a lull despite the high frequency of eruption 1 eruption 9 months on average temporal and spatial variations in recurrence rate and eruptive dynamics of basaltic volcanism such as those recognised at piton de la fournaise should be considered in the formulation of hazard assessments and in the interpretation of precursory patterns

the 2007 eruptions and caldera collapse of the piton de la fournaise volcano la reunion island from tilt analysis at a single very broadband seismic station seismic records from la reunion island very broadband geoscope station are investigated to constrain the link between the 2007 eruptive sequence and the related caldera collapse of the piton de la fournaise volcano tilt estimated from seismic records reveals that the three 2007 eruptions belong to a single inflation deflation cycle tilt trend indicates that the small volume summit eruption of 18 february occurred during a phase of continuous inflation that started in january 2007 inflation decelerated 24 days before a second short lived small volume eruption on 30 march almost simultaneous with a sudden large scale deflation of the volcano deflation rate which had stabilized at relatively low level increased anew on 1 april while no magma was erupted followed on 2 april by a major distal eruption and on 5 april by a summit caldera collapse long term tilt variation suggests that the 2007 eruptive succession was triggered by a deep magma input

water co2 cl and f in melt inclusions in phenocrysts from three holocene explosive eruptions crater lake oregon rare melt inclusions 100 m in diameter trapped near the boundaries of corroded patchy zones in plagioclase phenocrysts from plinian pumice of three holocene eruptions were analyzed by ir spectroscopy for molecular h 2o oh groups and co 2 and by electron microprobe for cl and f the three rhyodacitic eruptions each of which began with aplinian phase occurred over 200 yr the llao rock and cleetwood eruptions ended with degassed lava flows and the subsequent climactic eruption with voluminous ignimbrite groundmass glass and melt inclusion compositions anhydrous are similar inclusions with total h 2o concentrations of 1 2 3 4 wt are adjacent to fractures or are hourglass inclusions suggesting partial degassing caused by depressurization melt inclusions in phenocrysts in climactic ignimbrite may have vesiculated for the same reason intact inclusions in plinian pumices have total h 2o concentrations believed to represent magmatic h 2o contents wt llao rock 5 3 5 3 cleetwood 3 8 4 7 climactic 3 9 0 2 l n 6 ratios of or molecular h 2o indicate closure temperatures of 200 500oc that reflect syn or posteruptive cooling co 2 concentrations are 25 ppm h 2o and co 2 concentrations indicate saturation pressures of 1 0 1 8 kbar or depths 5 km six inclusions from climactic pumice average 400 40 l ppm f 1880 70 ppm cl f and cl concentrations of the other samples are similar and not well correlated with each other or with total h 2o location of melt inclusions near boundaries of patchy zones which are mantled by oscillatory zoned overgrowths suggests that their h 2o concentrations represent magmatic values significantly before eruption although mean h 2o concentrations of analyzed melt inclusions decrease for the three successive eruptions it is not certain that this indicates a corresponding change in magmatic h 2o content during the interval between eruptions

source mechanism of seismic waves excited by pyroclastic flows observed at unzen volcano japan a number of pyroclastic flows have occurred at unzen volcano unzendake in kyushu japan since the end of may 1991 seismic waves which are excited by the pyroclastic flows contain low frequency components in their initial phases they were recorded by two long period seismographs in the seismic network deployed around the volcano this paper is an attempt to clarify the source mechanism of the seismic waves triggered by the pyroclastic flows visual data show that the pyroclastic flow which excites these seismic waves occurs in three successive steps first fall of a lava dome front dome collapse then its collision with the mountain slope and its fragmentation into small blocks and finally descent of the blocks along the slope we can build a simple mechanical model which gives theoretical time histories of vertical and horizontal components of a force acting on the spot of dome collapse on the mountain slope the model of a mechanism generating these seismic waves was constructed through analysis of the low frequency components with the aid of inversion methods which can estimate source time function of the force acting on the source the estimated time histories are close to those predicted by the simple mechanical model we conclude that the three steps mentioned above are essential in exciting the low frequency seismic signals this idea is supported by the fact that the volume of fallen material estimated from the visual data is proportional to the maximum vertical amplitude of seismic waves triggered by the pyroclastic flow this relation holds for the two pyroclastic flows of june 22 1991 and may 28 1992 which were studied in this paper these findings seismically confirm a visual observation that the pyroclastic flow is a result of dome collapse

the role of volatiles in magma chamber dynamics many andesitic volcanoes exhibit effusive eruption activity1 with magma volumes as large as 107 109 m3 erupted at rates of 1 10 m3 s 1 over periods of years or decades during such eruptions many complex cycles in eruption rates have been observed with periods ranging from hours to years2 3 4 5 6 7 longer term trends have also been observed and are thought to be associated with the continuing recharge of magma from deep in the crust and with waning of overpressure in the magma reservoir here we present a model which incorporates effects due to compressibility of gas in magma we show that the eruption duration and volume of erupted magma may increase by up to two orders of magnitude if the stored internal energy associated with dissolved volatiles can be released into the magma chamber this mechanism would be favoured in shallow chambers or volatile rich magmas and the cooling of magma by country rock may enhance this release of energy leading to substantial increases in eruption rate and duration

structures and evolution of the plumbing system of piton de la fournaise volcano inferred from clustering of 2007 eruptive cycle seismicity analysis of seismic activity associated with the eruptions of 2007 which led to the collapse of the dolomieu crater on april 5th reveals the link between the seismicity and the magma transfers at piton de la fournaise three eruptive phases occurred on february 18th march 30th and april 2nd 2007 at the summit 2 000 m and 600 m high on the south east flank respectively illustrating the three types of eruptions defined for the current piton de la fournaise activity we use cross correlation of seismic waveforms and clustering to improve the earthquake locations and determine the best constrained focal mechanisms with an average of 78 p phase polarities the pre eruptive seismicity of the february and march eruptions is composed of time extended clusters that also preceded other distal eruptions from 2000 to 2007 our analysis shows that the seismic swarm prior to the february eruption initiated the intrusion that led to the april 2007 eruption the seismicity preceding the dolomieu crater collapse consists of numerous but time limited clusters and specific seismic activity that accompanied the dolomieu crater collapse from april 1st to april 5th earthquakes with clvd mechanisms combined with normal faulting sources occurred between 0 8 and 0 km asl until the complete rupture of the shallow magma storage roof this collapse induced the propagation of a de pressurization front and triggered a migration of seismicity from 0 to 8 km along a very narrow path

triggering basaltic volcanic eruptions by bubble melt separation understanding the processes by which volcanic eruptions are triggered is crucial for volcanic hazard prediction and assessment in intermediate and silicic magmas where water is the dominant volatile phase1 2 second boiling provides an effective eruption trigger1 as the magma cools and crystallizes volatiles are increasingly concentrated in the residual melt where they eventually become saturated and are exsolved thereby raising the magma pressure to the point of eruption1 2 in contrast in basaltic magma chambers in which carbon dioxide is the dominant volatile phase the effect of second boiling is to decrease the chamber pressure but basaltic magmas are also relatively inviscid so volatile bubbles can separate from the turbulently convecting melt to produce a foam layer at the chamber roof3 4 decompression of the rising bubbles can increase the chamber pressure and so trigger an eruption5 6 here we model the complex competition between the processes of magma cooling and bubble ascent we find that in basalt of very low viscosity 1 10pas bubble separation dominates and can increase the reservoir pressure by as much as 10 mpa sufficient to trigger an eruption7 on timescales of 1 10 years comparable to the duration of eruptive episodes at kilauea volcano hawaii4 8

the 2007 caldera collapse of piton de la fournaise volcano source process from very long period seismic signals abstract in april 2007 piton de la fournaise volcano experienced its largest caldera collapse in at least 300 yr this event consisted of a series of 48 subsidence increments characterized by very long period vlp seismic signals equivalent to m w 4 4 to 5 4 source analysis of vlp events suggests a piston like mechanism with a collapsing crack source corresponding to the contraction of the magma reservoir and a single force representing the collapse of the above rock column we show that the collapse dynamics is primarily controlled by magma withdrawal from the reservoir which was likely in a bubbly state at the time of the event similar to the 2018 kilauea collapse we observe a reduction of the time interval between successive subsidence increments which results from the acceleration of magma withdrawal and a progressive weakening of the edifice at the beginning of the sequence

analysis of seismic body waves excited by the mount st helens eruption of may 18 1980 seismic body waves which were excited by the may 18 1980 eruption of mount st helens and recorded by the global digital seismographic network stations are analyzed to determine the nature and the time sequence of the events associated with the eruption the polarity of teleseismic p waves period 20s is identical at six stations which are distributed over a wide azimuthal range this observation together with a very small s to p amplitude ratio at 20 s suggests that the source is a nearly vertical single force a simple model shows that for seismic radiation a volcanic eruption can be represented by a single force applied in the direction opposite to the blast direction the time history of the vertical force suggests two distinct groups of events about 2 min apart each consisting of several subevents with a duration of about 25 s the magnitude of the force is approximately 2 6 10 12 n this vertical force is in contrast with the long period 150 s southward horizontal single force which has been determined by a previous study and interpreted to be due to the massive landslide an m s 5 2 earthquake initiated the eruption sequence the direction of the p wave first motion of this event observed at two nearby stations is consistent with the radiation pattern expected for the landslide and suggests that it represents the onset of the landslide the ground motions observed at station lon 67 km are dominated by rayleigh waves i e lamb pulse and provide tight constraints on the time sequence of the events

toward continuous quantification of lava extrusion rate results from the multidisciplinary analysis of the 2 january 2010 eruption of piton de la fournaise volcano la reunion the dynamics of the 2 12 january 2010 effusive eruption at piton de la fournaise volcano were examined through seismic and infrasound records time lapse photography so2 flux measurements deformation data and direct observations digital elevation models were constructed for four periods of the eruption thus providing an assessment of the temporal evolution of the morphology the volume and the extrusion rate of the lava flow these data were compared to the continuous recording of the seismic and infrasonic waves and a linear relationship was found between the seismic energy of the tremor and the lava extrusion rate this relationship is supported by data from three other summit eruptions of piton de la fournaise and gives total volume and average lava extrusion rate in good agreement with previous studies we can therefore provide an estimate of the lava extrusion rate for the january 2010 eruption with a very high temporal resolution we found an average lava extrusion rate of 2 4 m3s 1 with a peak of 106 6 m3s 1 during the initial lava fountaining phase we use the inferred average lava extrusion rate during the lava fountaining phase 30 23 m3s 1 to estimate the value of the initial overpressure in the magma reservoir which we found to range from 3 7 106 pa to 5 9 106 pa finally based on the estimated initial overpressure the volume of magma expelled during the lava fountaining phase and geodetic data we inferred the volume of the magma reservoir using a simple mogi model between 0 25 km3 and 0 54 km3 which is in good agreement with previous studies

real time monitoring of relative velocity changes using ambient seismic noise at the piton de la fournaise volcano la reunion from january 2006 to june 2007 we present the results of a real time method based on coda wave interferometry from seismic noise cross correlation functions for relative seismic velocity variations monitoring on a volcanic edifice the ambient seismic noise at the piton de la fournaise volcano on la reunion island is analyzed from january 2006 to june 2007 during this period five eruptions occurred showing a great diversity in eruption duration intensity and eruptive fissure location two different methods are used to compute the velocity variations in order to compare their stability in quasi real time routine we compare the obtained velocity variations with the surface deformation observed by gps and extensometers networks this allows us to identify and quantify three major processes at the origin of seismic wave velocity variations in the edifice firstly the observed pre eruptive summit inflation is accompanied by a decrease in seismic velocity secondly the edifice deflation following the opening of an eruptive fissure is characterized by an increase of the velocity finally the summit caldera collapse generates a strong velocity drop coda wave interferometry from seismic noise cross correlation functions in quasi real time may allow us to forecast eruption and constrain the processes taking place in the volcanic plumbing system

formation of the april 2007 caldera collapse at piton de la fournaise volcano insights from gps data abstract on april 5th 2007 a major collapse of the dolomieu summit crater occurred at piton de la fournaise volcano la reunion island to constrain the origin and the dynamics of the collapse we modelled continuous gps displacements preceding and accompanying this event our results reveal that the dolomieu collapse formed by a piston like mechanism in which the fast draining of the shallow magma reservoir at the beginning of the march 30th may 1st 2007 eruption caused a gravity driven downward displacement of its roof during the five days preceding the collapse summit deflation was caused by volume change through the fractured rock column above the shallow magma chamber due to the closure of voids or small scale cavities or the upward stoping and migration of an underground cavity since 2000 weakening of the rock column by repeated refilling and draining of the shallow magma chamber favoured its destabilization on april 5th 7 h before the onset of the collapse the deflation source was located at 300 m depth revealing the initiation of the collapse at very shallow depth the sudden collapse on april 5th was followed by continuous cyclic subsidence until april 6th each cycle was characterized by a progressive inward deflation of the summit zone and ended by sudden summit outward displacements caused by stress relaxation following a collapse event subsidence of the stoped column in the magma reservoir during the collapse of the crater on april 5th and 6th acted as a piston thus increasing the eruption rate the collapse occurred in only 24 h by successive events and increased the depth of the dolomieu crater by 340 m after april 6th only minor readjustments of the summit cone occurred with slight summit deflation and small landslides of the caldera walls