Temporal evolution of long-period seismicity at Etna Volcano, Italy, and its relationships with the 2004–2005 eruption

Between December 2004 and August 2005, more than 50,000 long-period events (LP) accompanied by very-long period pulses (VLP) were recorded at Mt. Etna, encompassing the effusive eruption which started in September 2004. The observed activity can be explained by the injection of a gas slug formed within the magmatic column into an overlying cavity filled by either magmatic or hydrothermal fluids, thus triggering cavity resonance. Although a large number of LP events exhibit similar waveforms before the eruption, they change significantly during and after the eruption. We study the temporal evolution of the LP-VLP activity in terms of the source movement, change of the waveforms, temporal evolution of the dominant resonance frequencies and the source Q factor and changes in the polarization of the signal. The LP source locations before and after the eruption, respectively, do not move significantly, while a slight movement of the VLP source is found. The intensity of the LP events increases after the eruption as well as their dominant frequency and Q factor, while the polarization of the signals changes from predominantly transversal to pure radial motion. Although in previous studies a link between the observed LP activity and the eruption was not found, these observations suggest that such a link was established at the latter end of the eruptive sequence, most likely as a consequence of a reestablishment of the pressure balance in the plumbing system, after it was undermined due to the discharge of large amounts of resident magma during the eruption. Based on the polarization properties of the signal and geological setting of the area, a fluid-filled crack is proposed as the most likely source geometry. The spectral analysis based on the autoregressive-models (SOMPI) is applied to the signals in order to analyse the resonance frequencies and the source Q-factors. The results suggest water and basalt at low gas volume fraction as the most likely fluids involved in the source process. Using theoretical relations for the “slow waves” radiated from the fluid-filled crack, we also estimate the crack size for both fluids, respectively

Moment tensor inversion of Explosive Long Period events recorded on Arenal Volcano, Costa Rica, constrained by synthetic tests

International audienceIn order to constrain the moment tensor solution of an explosive seismic event recorded on Arenal volcano, Costa Rica, we perform tests using synthetic data. These data are generated using a 3D model including the topography of the volcano and the best estimation of the velocity model available for Arenal. Solutions for (i) the moment tensor components, and (ii) the moment tensor plus single forces, are analyzed. When noisy data and mislocated sources are used in the inversion, spurious single forces are easily generated in the solution for the moment tensor components plus single forces. Forces also appear when the inversion is performed using an explosive event recorded on Arenal in 2005. Synthetic tests indicate that these forces might be spurious. However the mechanism is correctly retrieved by the inversion in both solutions. The ability to recover the explosive mechanism for the 2005 event combined with the interpretative aids from the synthetics tests will enable us to invert for the large variation in events observed on Arenal

Temporal evolution of Long-Period seismicity at Etna Volcano, Italy, and its relationships with the 2004-2005 eruption

Between December 2004 and August 2005, more than 50,000 long-period events (LP) accompanied by very-long period pulses (VLP) were recorded at Mt. Etna, encompassing the effusive eruption which started in September 2004. The observed activity can be explained by the injection of a gas slug formed within the magmatic column into an overlying cavity filled by either magmatic or hydrothermal fluids, thus triggering cavity resonance. Although a large number of LP events exhibit similar waveforms before the eruption, they change significantly during and after the eruption. We study the temporal evolution of the LP-VLP activity in terms of the source movement, change of the waveforms, temporal evolution of the dominant resonance frequencies and the source Q factor and changes in the polarization of the signal. The LP source locations before and after the eruption, respectively, do not move significantly, while a slight movement of the VLP source is found. The intensity of the LP events increases after the eruption as well as their dominant frequency and Q factor, while the polarization of the signals changes from predominantly transversal to pure radial motion. Although in previous studies a link between the observed LP activity and the eruption was not found, these observations suggest that such a link was established at the latter end of the eruptive sequence, most likely as a consequence of a reestablishment of the pressure balance in the plumbing system, after it was undermined due to discharge of large amounts of resident magma during the eruption. Based on the polarization properties of the signal and geological setting of the area, a fluid- filled crack is proposed as the most likely source geometry. The spectral analysis based on the autoregressive-models (SOMPI) is applied to the signals in order to analyse the resonance frequencies and the source Q-factors. The results suggest water and basalt with the low gas volume fraction as the most likely fluids involved in the source process. Using theoretical relations for the “slow waves” radiated from the fluid-filled crack, we also estimate the crack size for both fluids, respectively

Analysis of sustained long-period activity at Etna Volcano, Italy

Following the installation of a broadband network on Mt. Etna, sustained Long- Period (LP) activity was recorded accompanying a period of total quiescence and the subsequent onset of the 2004–2005 effusive episode. From c. about 56000 events detected by an automatic classification procedure, we analyse a subset of about 3000 signals spanning the December 17th, 2003–September 25th, 2004, time interval. LP spectra are characterised by several, unevenly-spaced narrow peaks spanning the 0.5–10 Hz frequency band. These peaks are common to all the recording sites of the network, and different from those associated with tremor signals. Throughout the analysed time interval, LP spectra and waveforms maintain significant similarity, thus indicating the involvement of a non-destructive source process that we interpret in terms of the resonance of a fluid-filled buried cavity. Polarisation analysis indicates radiation from a non-isotropic source involving large amounts of shear. Concurrently with LP signals, recordings from the summit station also depict Very-Long-Period (VLP) pulses whose rectilinear motion points to a region located beneath the summit craters at depths ranging between 800 and 1100 m beneath the surface. Based on a refined repicking of similar waveforms, we obtain robust locations for a selected subset of the most energetic LP events from probabilistic inversion of travel-times calculated for a 3D heterogenous structure. LP sources cluster in a narrow volume located beneath the summit craters, and extending to a maximum depth of ≈800 m beneath the surface. No causal relationships are observed between LP, VLP and tremor activities and the onset of the 2004–2005 lava effusions, thus indicating that magmatic overpressure played a limited role in triggering this eruption. These data represent the very first observation of LP and VLP activity at Etna during non-eruptive periods, and open the way to the quantitative modelling of the geometry and dynamics of the shallow plumbing system

Seismic source mechanisms of tremor recorded on Arenal volcano, Costa Rica, retrieved by waveform inversion

International audienceMoment tensor inversions of volcanic tremor are synthetically tested and subsequently applied to a dataset recorded on Arenal volcano in 2005. We selected 15 harmonic tremor events showing an emergent but very clear onset with a fundamental frequency range of 0.8-1.5 Hz. These events constitute an excellent opportunity to broaden our knowledge of tremor generation models. The inversions were performed for a common source location, the position of which was retrieved through the evaluation of the joint probability density function of the misfit values obtained by the moment tensor inversion of all the events and all predefined source locations (grid search). Events are located beneath the summit crater, in a shallow position. The inversion procedure was tested through the use of synthetic tremor, generated using full wavefield numerical simulations. The maximum amplitude decomposition method is used to analyse the solutions of the synthetic tests. The results highlight the stability of the moment tensor solution when the whole length of the signal is used in the inversion procedure. Hence the whole length of the tremor is utilised to retrieve the source mechanism generating the 15 tremor events. A sliding window is utilised to investigate the time variability of our solution. A crack dipping 20° to the North-Northeast is reconcilable with all 15 tremor events. This mechanism is found to be constant for the entire length of each tremor and for different events

Hemolymph microbiome of Pacific oysters in response to temperature, temperature stress and infection

Microbiota provide their hosts with a range of beneficial services, including defense from external pathogens. However, host-associated microbial communities themselves can act as a source of opportunistic pathogens depending on the environment. Marine poikilotherms and their microbiota are strongly influenced by temperature, but experimental studies exploring how temperature affects the interactions between both parties are rare. To assess the effects of temperature, temperature stress and infection on diversity, composition and dynamics of the hemolymph microbiota of Pacific oysters (Crassostrea gigas), we conducted an experiment in a fully-crossed, three-factorial design, in which the temperature acclimated oysters (8 or 22 °C) were exposed to temperature stress and to experimental challenge with a virulent Vibrio sp. Strain. We monitored oyster survival and repeatedly collected hemolymph of dead and alive animals to determine the microbiome composition by 16s rRNA gene amplicon pyrosequencing. We found that the microbial dynamics and composition of communities in healthy animals (including infection survivors) were significantly affected by temperature and temperature stress, but not by infection. The response was mediated by changes in the incidence and abundance of operational taxonomic units (OTUs) and accompanied by little change at higher taxonomic levels, indicating dynamic stability of the hemolymph microbiome. Dead and moribund oysters, on the contrary, displayed signs of community structure disruption, characterized by very low diversity and proliferation of few OTUs. We can therefore link short-term responses of host-associated microbial communities to abiotic and biotic factors and assess the potential feedback between microbiota dynamics and host survival during disease

Shallow Tectonic Stress Magnitudes at the Hikurangi Subduction Margin, New Zealand

AbstractQuantifying tectonic stress magnitudes is crucial in understanding crustal deformation processes, fault geomechanics, and variable plate interface slip behaviors in subduction zones. The Hikurangi Subduction Margin (HSM), New Zealand, is characterized by along‐strike variation in interface slip behavior, which may be linked to tectonic stress variations within the overriding plate. This study constrains in situ stress magnitudes of the shallow (&lt;3 km) overriding plate of the HSM to better understand its tectonics and how they relate to larger scale subduction dynamics. Results reveal σ3: Sv ratios of 0.6–1 at depths above 650–700 m TVD and 0.92–1 below this depth interval along the entire HSM. Additionally, for depths below 650–700 m TVD, SHmax: Sv ratios of 0.95–1.81 in the central HSM and 0.95–2.3 in the southern HSM are estimated. These stress ratios suggest a prevalent thrust to strike‐slip (σ1 = SHmax) faulting regime across the central and southern HSM. In the central HSM, the presence of NNE‐NE striking reverse faults co‐existing with a modern σ1 (SHmax) aligned ENE‐WSW suggests that overtime the stress state here evolved from a contractional to a strike‐slip state, where the compressional direction changes from perpendicular (NW‐SE) to oblique (ENE‐WSW) to the Hikurangi margin. This temporal change in stress state may be explained by forearc rotation, likely combined with the development of upper plate overpressures. In the southern HSM, the modern WNW‐ESE/NW‐SE σ1 (SHmax) and pre‐existing NNE‐NE striking reverse faults indicate that stress state remains contractional and perpendicular (NW‐SE) to the Hikurangi margin overtime.</jats:p

Amplification of strong ground motion in the city of Zagreb, Croatia, estimated by computation of synthetic seismograms. Soil Dynamics and Earthquake Engineering

A hybrid technique consisting of modal summation and subsequent \uaenite differences modelling is applied for the computation of synthetic accelerograms along a pro\uaele crossing the city of Zagreb, the capital of Croatia. Assuming the source geometry is known, the ampli\uaecation properties of the underlying soil may be determined by comparison of synthetics and their response spectra computed for a bedrock model with the ones obtained under the assumption of a realistic laterally varying local model. The peak ground acceleration is larger by a factor of up to 3.5 than the value obtained for the bedrock model. The ampli\uaecation of the response spectra is most prominent for frequencies below 2 Hz, and increases sharply to the SW from the mapped fault running through the centre of the city. q 2002 Published by Elsevier Science Ltd

Moment-tensor inversion of LP events recorded on Etna in 2004 using constraints obtained from wave simulation tests

The persistent occurrence of long period (LP) events at Mt Etna became apparent with the installation of the first fixed broad-band seismic network in late 2003. Repeating similar LP events from Nov. ‘03 to Sept. ‘04 indicate a non-destructive source process. We perform moment tensor (MT) inversions on a stacked high S/N ratio representative LP signal, conducting a grid search for the source geometry and L2-inversion for the source time function. Results indicate a NNW-SSE oriented resonating sub-vertical crack as the most probable source. This result is consistent with deformation and GPS observations. Crucial to this result are constraints imposed by detailed 3D full waveform numerical simulations in a heterogeneous tomographic model with topography, and in particular a detailed assessment of the influence of very near surface velocity structure on LP signals. Pulsating gas injection is hypothesised as the most likely LP trigger

Long Period (LP) Events on Mt. Etna volcano (Italy): the influence of velocity structures on moment tensor inversion

International audienceSince a few decades volcanic Long Period (LP) events have been recorded on many active volcanoes and their study has been recognized as an important tool to characterize volcanic activity. LP event analyses through moment tensor (MT) inversions have led to kinematic descriptions of various source mechanisms. The main challenge in these inversions is to " strip out " the propagation effect in order to isolate the source; hence the velocity model used controls the accuracy of the retrieved source mechanism. We first carry out several synthetic tests of inversions on Mt. Etna volcano (Italy). Four geological models with topography are considered with increasing complexity: the most complex model is used to generate synthetic data, while the other three models are used to calculate the Greens' functions for inversions. The retrieved solutions from the three velocity models are similar. The MT solutions for a deeper source are well retrieved, while a shallower source test suffers from high uncertainties and strong misinterpretation of the source orientation. The homogeneous model gives the lowest misfit value, but source location and mechanism decomposition are inaccurate. When a complex model different from the true one is used, a high misfit value and a wrong solution is obtained. We then incorporate our findings into the MT inversion of an LP event recorded on Mt Etna in 2008. We obtain very different solutions among the three models in terms of source location and mechanism decomposition. The overall shape of the retrieved source time functions are similar, but some amplitude differences arise, especially for the homogeneous model. Our work highlights the importance of including the unconsolidated surface materials in the computation of Green’s functions especially when dealing with shallow sources