Reconstruction of the volcanic history of the Tacámbaro-Puruarán area (Michoacán, México) reveals high frequency of Holocene monogenetic eruptions

The 690?km² Tacámbaro-Puruarán area located at the arc-front part of the Michoácan-Guanajuato volcanic field in the Trans-Mexican Volcanic Belt (TMVB) records a protracted history of volcanism that culminated with intense monogenetic activity in the Holocene. Geologic mapping, ⁴⁰Ar/³⁹Ar and ¹⁴C radiometric dating, and whole-rock chemical analyses of volcanic products provide insights to that history. Eocene volcanics (55–40?Ma) exposed at uplifted blocks are related to a magmatic arc that preceded the TMVB. Early TMVB products are represented by poorly exposed Pliocene silicic domes (5–2?Ma). Quaternary (<2?Ma) volcanoes (114 mapped) are mainly scoria cones with lavas (49 vol.%), viscous lava flows (22 vol.%), and lava shields (22 vol.%). Erupted products are dominantly either basaltic andesites (37 vol. %), or andesites (17 vol.%), or span across both compositions (28 vol.%). Basalts (9 vol.%), dacites (4 vol.%), shoshonites (2 vol.%), and other alkali-rich rocks (<3 vol.%) occur subordinately. Early-Pleistocene volcanism was bimodal (dacites and basalts) and voluminous while since 1?Ma small-volume eruptions of intermediate magmas have dominated. Higher rates of lithospheric extension in the Quaternary may have allowed a larger number of small, poorly evolved dikes to reach the surface during this period. Eruptive centers as old as 1.7?Ma are aligned in a NE direction parallel to both, basement faults and the direction of regional compressive stress, implying structural control on volcanic activity. Data suggest that volcanism was strongly pulsatory and fed by localized low-degree partial melting of mantle sources. In the Holocene, at least 13 eruptions occurred (average recurrence interval of 800?years). These produced ~3.8?km³ of basaltic andesitic to andesitic magma and included four eruptions dated at ~1,000; 4,000; 8,000; and 11,000?years bc (calibrated ¹⁴C ages). To date, this is one of the highest monogenetic eruption frequencies detected within such a small area in a subduction-related arc-setting. These anomalous rates of monogenetic activity in an area with thick crust (>30?km) may be related to high rates of magma production at depth and a favorable tectonic setting.     

Correlations of three-component ground motions from the 23 June 2010 Val-des-Bois, Québec, earthquake

In this paper, we investigate the characteristics of three-component ground motions recorded during the $M_\text{w}\,5.2$ Val-des-Bois (Québec) earthquake, which occurred on the 23 June 2010. The earthquake is the largest recorded event in eastern Canada within the last decade. The records analyzed were provided by a strong motion monitoring network, comprising accelerometers located at sites with different soil conditions. The two orthogonal horizontal components and one vertical component at each recording station are uncorrelated to determine their principal directions, and the results obtained are used to characterize intensity ratios between the three uncorrelated components. A new hodograph representation is proposed to highlight the correlation between accelerations and displacement trajectories along various time increments at each recording station. The principal components are discussed in light of seismographic data, local site conditions, and trajectories. Time–frequency analyses of the uncorrelated records are also conducted to compare the distribution of spectral amplitudes and frequency content along the three principal components during the shaking. The results of this work shed more light on the characteristics of three-component ground motions from an important Eastern North America earthquake, and could be used to calibrate simulated multicomponent ground motions in this region.     

Predicted ground motion after the L��Aquila 2009 earthquake (Italy, M w 6.3): input spectra for seismic microzoning

After the April 6th 2009 L’Aquila earthquake (M _w 6.3), where 306 people died and a further 60,000 were displaced, seismic microzoning investigations have been carried out for towns affected by a macroseismic intensity equal to or greater than 7 MCS. Based upon seismotectonic data, historical seismicity and strong motion records, we defined input spectra to be used in the numerical simulations of seismic microzoning in four key municipalities, including the town of L’Aquila. We adopted two main approaches: uniform hazard response spectra are obtained by a probabilistic seismic hazard assessment introducing some time-dependency for individual faults on the study area; a deterministic design spectrum is computed from magnitude/distance pairs extracted by a stationary probabilistic analysis of historical intensities. The uniform hazard spectrum of the present Italian building code represents the third, less restrictive, response spectrum to be used for the numerical simulations in seismic microzoning. Strong motions recordings of the main shock of the L’Aquila sequence enlighten the critical role played by both the local response and distances metric for sites located above a seismogenic fault; however, these time-histories are compatible with the uncertainties of a deterministic utilization of ground motion predictive equations. As recordings at very near field are rare, they cannot be neglected while defining the seismic input. Disaggregation on the non-Possonian seismotectonic analysis and on the stationary site-intensity estimates reach very similar results in magnitude-distance pairs identification; we interpret this convergence as a validation of the geology-based model by historical observations.     

Self-sustained oscillations at Volcán de Colima (México) inferred by Independent Component Analysis

We have analyzed Long-Period (LP) seismic events of Volcán de Colima (Mexico) recorded by four three-component broad-band seismometers during January 2006. Frequency-domain analysis shows spectra mainly peaked in the range 0.3–2 Hz and characterized by a monotonic decreasing envelope, as those produced in a cylindrically-symmetric self-oscillating cavity interacting with a confined jet. Independent Component Analysis, a time-series decomposition method, extracts two or three nonlinear oscillation modes depending on the station, with the fundamental one peaked at 0.4–0.5 Hz. These decomposed waves have a well defined spectral content and are self-sustained oscillations with low dimensionality, as shown by dynamic system reconstruction methods. Moreover, they show radial polarization in near field and transverse polarization in far field in North-West South-East direction. Finally, we hypothesize a branched plumbing system with two principal conduits having alignments orthogonal to one another.     

Rupture process of the Lorca (southeast Spain) 11 May 2011 (M w = 5.1) earthquake

On 11 May 2011, a M _w ?=?5.1 earthquake shook the town of Lorca (SE Spain) causing a disproportionately large damage for its magnitude. In order to contribute to knowledge of the behavior of the active faults present in the region and define the parameters which control their motion, we made a detailed study of the rupture process of this earthquake from inversion of body waves at regional and teleseismic distances. Ground motion displacements obtained in this way are in agreement with near-field strong motion data and GPS observations recorded in Lorca. We have obtained a partly bilateral rupture propagating to WSW (238°, 54°, 59°) with 27 cm of maximum slip and shallow focus (4 km). The fault plane orientation corresponds to that of the Cejo de los Enamorados Fault located NE of the Lorca town and parallel to the Alhama de Murcia Fault. The distribution of slip on the fault plane can explain the lack of any observed surface rupture as we found that the rupture started at 4-km depth along a plane dipping at 54°, with motion propagating upward to stop at 1.5 km below the surface. The strong motion and GPS data recorded near the epicenter are in agreement with the maximum slip on the fault. Directivity effects and the extreme shallowness of the rupture could explain the considerable damage that the earthquake caused in the town of Lorca.     

Seismicity of the Baikal rift zone for the digital recording period of earthquake observation (2001–2006)

Some characteristics of the seismic regime of the Baikal rift zone for the analogue (1960–2000) and digital (2001–2006) registration modes are considered. It was ascertained that the general number of seismic events recorded by digital equipment increased by two times in comparison with the analogue observation period and the catalog completeness level was lowered to K _P = 7.0 as well. The focal mechanisms of 110 earth-quakes with K _P ≥ 8.0 were examined. The overwhelming majority of them during 2001–2006 were due to normal faulting and normal faulting with strike-slip components. Special attention was given to study of the strong events (K _P ≥ 13.5 and MPSP ≥ 4.7), whose focal solutions reveal the regional stress field.     

A Mw 6.3 earthquake scenario in the city of Nice (southeast France): ground motion simulations

The southern Alps–Ligurian basin junction is one of the most seismically active zone of the western Europe. A constant microseismicity and moderate size events (3.5 < M < 5) are regularly recorded. The last reported historical event took place in February 1887 and reached an estimated magnitude between 6 and 6.5, causing human losses and extensive damages (intensity X, Medvedev–Sponheuer–Karnik). Such an event, occurring nowadays, could have critical consequences given the high density of population living on the French and Italian Riviera. We study the case of an offshore Mw 6.3 earthquake located at the place where two moderate size events (Mw 4.5) occurred recently and where a morphotectonic feature has been detected by a bathymetric survey. We used a stochastic empirical Green’s functions (EGFs) summation method to produce a population of realistic accelerograms on rock and soil sites in the city of Nice. The ground motion simulations are calibrated on a rock site with a set of ground motion prediction equations (GMPEs) in order to estimate a reasonable stress-drop ratio between the February 25th, 2001, Mw 4.5, event taken as an EGF and the target earthquake. Our results show that the combination of the GMPEs and EGF techniques is an interesting tool for site-specific strong ground motion estimation.     

The characteristics of seismic signals produced by lahars and pyroclastic flows: Volcán de Colima,México

The general characteristics of seismic signals produced by pyroclastic flows (generated by either the collapse of a lava dome or an eruptive column) and lahars at Volcán de Colima, México are discussed. The paper concentrates on the 2004–2006 activity associated with and following the extrusion of andesitic block-lava in October–November 2004. It is shown that the duration of the broad-band seismic records of pyroclastic flows lasts a few minutes while the duration of seismic records of lahars continues for tens of minutes or hours. The spectra of seismic records produced by pyroclastic flows are characterized by lower peak frequencies (around 3–4 Hz) than for lahars (around 6–8 Hz). This difference in the frequency content together with the difference in the duration of seismic signals allows early diagnostic of the events in real time.     

Evaluation of site effects in the Aterno river valley (Central Italy) from aftershocks of the 2009 L��Aquila earthquake

A temporary network of 33 seismic stations was deployed in the area struck by the 6th April 2009, Mw 6.1 (Scognamiglio et?al. in Seism Res Lett 6/81, 2010), L??Aquila earthquake (central Italy), with the aim to investigate the site amplification within the Aterno river Valley. The seismograms of 18 earthquakes recorded by 14 of the 33 stations were used to evaluate the average horizontal to vertical spectral ratio (HVSR) for each site and the standard horizontal spectral ratio (SSR) between a site and a reference station. The obtained results have been compared to the geological and geophysical information in order to explain the resonance frequencies and the amplification levels with respect to surface geology of the valley. The results indicate that there is no uniform pattern of amplification, because of the complex geologic setting, as the thickness and degree of cementation of the deposits is highly variable.     

Seismic tomography and anisotropy in the source area of the 2008 Iwate–Miyagi earthquake (M 7.2)

In order to understand the generation mechanism of the Iwate–Miyagi earthquake (M 7.2) that occurred on 14 June 2008 in Northeast Japan, we determined high-resolution 3D seismic velocity (Vp, Vs) and Poisson's ratio (σ) structures as well as P-wave azimuthal anisotropy in and around the source area using 246,623 P-wave and 189,153 S-wave high-quality arrival times from 1941 aftershocks of the 2008 Iwate–Miyagi earthquake as well as 4878 shallow and intermediate-depth earthquakes in this area recorded by the dense High-Sensitivity Seismic Network (Hi-net) and the seismic stations of Tohoku University. Our results show that hypocenters of the mainshock and three large aftershocks are located in a boundary zone where both seismic velocity and Poisson's ratio change drastically in a short distance. A zone with pronounced low-velocity and high Poisson's ratio is revealed in the lower crust and uppermost mantle under the source area, which may reflect the arc magma and fluids ascending from the upper-mantle wedge. The P-wave azimuthal anisotropy is complex under the study area, which may also indicate the complex crustal structure there. These results indicate that the generation of the 2008 Iwate–Miyagi earthquake was influenced by the ascending arc magma and fluids associated with the dehydration reactions of the subducting Pacific slab under Northeast Japan.     

Process of source dynamics of the Jingtai earthquake (M=6. 2 )

ＰｒｏｃｅｓｓｏｆｓｏｕｒｃｅｄｙｎａｍｉｃｓｏｆｔｈｅＪｉｎｇｔａｉｅａｒｔｈｑｕａｋｅ（Ｍ＝６．２）Ｘｕｅ－ＢｉｎＤＵ（杜学彬）；Ｚｈｉ－ＤｅＹＡＮ（闫志德）ａｎｄＭｉｎｇ－ＷｕＺＯＵ（邹明武）（ＥａｒｔｈｑｕａｋｅＲｅｓｅａｒｃｈＩｎｓｔｉｔｕ...     

A long-term earthquake forecast for the Kuril-Kamchatka arc for the period from September 2011 to August 2016. The likely location,time, and evolution of the next great earthquake with M ≥ 7.7 in Kamchatka

We consider the results from the ongoing 2010–2011 work on long-term earthquake prediction for the Kuril-Kamchatka arc based on the pattern of seismic gaps and the seismic cycle. We develop a forecast for the next 5 years, from September 2011 to August 2016, for all segments of the Kuril-Kamchatka arc earthquake-generating zone. For 20 segments we predict the appropriate phases of the seismic cycle, the normalized rate of small earthquakes (A₁₀), the magnitudes of moderate earthquakes to be expected with probability 0.8, 0.5, and 0.15, and the maximum possible magnitudes and probability of occurrence for great (M ≥ 7.7) earthquakes. This study serves as another confirmation that it is entirely necessary to continue the work in seismic retrofitting in the area of Petropavlovsk-Kamchatskii.     

Performance of the L’Aquila (central Italy) gas distribution network in the 2009 (\mathrm{{M}}_\mathrm{{w}} 6.3) earthquake

Despite the seismic vulnerability of gas systems and the significance of the direct and indirect consequences that loss functionality might have on large communities, the analysis of the earthquake performance and of post-earthquake management for this kind of distribution networks appears under-represented in the international literature, with respect to other lifeline systems. To contribute on this matter, the study presented comprises an investigation of the impact of L’Aquila 2009 earthquake ( $\text{ M }_\mathrm{w}$ 6.3) on the performance of the local medium- and low-pressure gas distribution networks. The assessment of the physical impact of the earthquake to the buried components of network, namely pipes, valves, and demand nodes, was carried out when processing post-earthquake repair activity reports. Repair data, along with geometrical and constructive features, were collected in a geographic information system linked to the digitized maps of the network, and compared with the interpolated map of recorded transient ground motion, measured in terms of peak ground velocity (i.e., a $Shakemap^\mathrm{TM}$ ). The impact of permanent ground deformation was also investigated and found to be limited in the study area. The resulting observed repair rates (number of repairs per km), presented for different pipeline materials, were compared with repair ratio fragility functions available in literature, showing relatively agreement especially to those for steel pipes, likely also because of the uncertainties in the estimations. Finally, the management of the L’Aquila gas system in the emergency phase and the resilience (functionality recover versus time) of the system was discussed.     

Volcanomagnetic signals during the recent Popocatépetl (México) eruptions and their relation to eruptive activity

An interdisciplinary approach correlating magnetic anomalies with composition of the ejecta in each eruption, as well as with seismicity, was used to study the effect of magmatic activity on the local magnetic record at Popocatépetl Volcano located 65 km southeast of México City. Eruptions began on December, 1994, and have continued with dome growth and ash emissions since then. The Tlamacas (TLA) geomagnetic total field monitoring station, located 5 km away from Popocatépetl’s crater, was installed in December, 1997, in order to detect magnetic anomalies induced by this activity.Spatial correlation and weighted difference methods were applied to detect temporal geomagnetic anomalies using TLA’s record and the Teoloyucan Magnetic Observatory as a reference station. Weighted differences were applied to cancel the effects of non-vulcanogenic external field variations. Magnetic anomalies over a 2-year time span were classified into four types correlating them with geochemical, seismic and visual monitoring of the volcanic activity. Magnetic anomalies are believed to be caused by magma injection and gas pressure build-up, which is sensitive to vent morphology and clearing during eruption, although some anomalies appear to be thermally related, changes in the stress field are very important. Most magnetic anomalies are short time signals that reverse to baseline level. Decreasing anomalies (−0.5 to −6.8 nT) precede eruptions by 1–8 days.The presence of a mafic magmatic component was determined by mineral examination and silica and magnesium analyses on the ejecta from the 1997–1999 eruptions. Whole rock analyses ranged from dacitic (65% SiO₂) to andesitic (57% SiO₂) with 2–6.6% MgO. The higher MgO, lower silica samples contain forsteritic olivine (Fo90). SiO₂ does not increase and MgO does not increase with time, suggesting ascent of small magma pulses which are consistent with the magnetic data.     

Seismicity of an andesitic volcano during block-lava effusion: Volcán de Colima, México, November 1998–January 1999

The block-lava effusion at Volcán de Colima, México began on November 20, 1998, after 12 months of seismic activity, and ended about 80 days later. Three types of seismic events were observed during the lava effusion. Volcano—tectonic earthquakes occurred mainly at the very beginning and after the termination of lava effusion. Explosion earthquakes occurred frequently during the period of the maximum rate in lava effusion. The remainder of the seismic signals were associated with pyroclastic flows and rockfalls from the lava dome. These latter signals increased sharply in number at the onset of lava effusion. The rate of occurrence remained high when the lava discharge rate decreased but gradually decreased after the termination of lava effusion. Maximum daily durations of seismic signals are proportional to the daily volumetric output of lava, indicating the dependence of the number of pyroclastic flows on the rate of lava output. A log-log plot of seismic signal duration vs. number of events with this duration displays a linear relationship. The short-period seismic signals can be divided into three categories based on duration: short events with durations less than 100 s; intermediate events with durations between 100 and 250 s; and long events with durations longer than 250 s. We infer that long events correspond to pyroclastic flows with mean deposit volume 2×10⁵ m³, and intermediate events represent pyroclastic flows with mean deposit volume 1×10³ m³.Editorial responsibility: J McPhie     

Dissolution/crystallization kinetics recorded in the 2002–2003 lavas of Stromboli (Italy)

On 28 December 2002, new vents opened on the flanks of Stromboli, just below the summit craters, interrupting the persistent activity of the volcano with a 7-month-long effusive eruption. We here report on the plagioclase size distribution (PlgSD) in lava samples collected following the chronology of the 2002–2003 eruption. Data reveal a linear PlgSD similar to that found in samples of normal Stromboli activity, indicating that the switch from Strombolian explosive to effusive activity is not associated with changes in texture. Nevertheless, the crystal size distribution slopes and intercepts exhibit slight sinusoidal temporal variations that are here ascribed to a magma supply mechanism able to induce “resonance” in the crystal size distribution, with an amplitude that depends on the supply rate.