Double-difference Relocation of the Aftershocks of the Tecom��n, Colima, Mexico Earthquake of 22 January 2003

On 22 January 2003, the M _w?=?7.6 Tecomán earthquake struck offshore of the state of Colima, Mexico, near the diffuse triple junction between the Cocos, Rivera, and North American plates. Three-hundred and fifty aftershocks of the Tecomán earthquake with magnitudes between 2.6 and 5.8, each recorded by at least 7 stations, are relocated using the double difference method. Initial locations are determined using P and S readings from the Red Sismológica Telemétrica del Estado de Colima (RESCO) and a 1-D velocity model. Because only eight RESCO stations were operating immediately following the Tecomán earthquake, uncertainties in the initial locations and depths are fairly large, with average uncertainties of 8.0?km in depth and 1.4?km in the north?Csouth and east?Cwest directions. Events occurring between 24 January and 31 January were located using not only RESCO phase readings but also additional P and S readings from 11 temporary stations. Average uncertainties decrease to 0.8?km in depth, 0.3?km in the east?Cwest direction, and 0.7?km in the north?Csouth direction for events occurring while the temporary stations were deployed. While some preliminary studies of the early aftershocks suggested that they were dominated by shallow events above the plate interface, our results place the majority of aftershocks along the plate interface, for a slab dipping between approximately 20° and 30°. This is consistent with the slab positions inferred from geodetic studies. We do see some upper plate aftershocks that may correspond to forearc fault zones, and faults inland in the upper plate, particularly among events occurring more than 3?months after the mainshock.     

Using an Enhanced Dataset for Reassessing the Source Region of the 2003 Armer��a, Mexico Earthquake

We present a fresh look at the source region of the 22 January 2003 M _w 7.4 Armería earthquake, which occurred off the Pacific coast of the state of Colima, Mexico, near the town of Armería. The effects of this earthquake in the neighboring states of Colima and Jalisco were different and stronger than those of previous recent major earthquakes in the region. This earthquake and its aftershocks were recorded by two local telemetered seismograph networks (RESCO and RESJAL). From 22 January to 24 January 2003, no important seismicity was located on the plates interface, or within the Rivera Plate, and most epicenters were located west of the Armería River, which is the western border of the Colima Graben, and is located outside of the Colima Gap region. From 24 January to 31 January, the seismicity recorded by both networks showed a migration in depth, with an almost vertical offshore distribution between 4 and 24?km in depth. For this period, a seven-station portable digital seismograph network, equipped with three-component seismometers, was deployed in the epicentral area to study the aftershock sequence in detail. With this denser network more than 200 M _L?>?2.0 aftershocks were recorded. The aftershock foci were deeper than those recorded during the early period and most of them locate on a hypothetical 12° dipping interface between the Rivera and North American Plates. Composite focal mechanism solutions for the aftershocks located during both periods indicate a reverse fault character that changes with time. Analysis of the new dataset still indicates that the earthquake was a shallow intraplate event.     

Properties of the Aftershock Sequences of the 2003 Bingöl, M D = 6.4, (Turkey) Earthquake

Aftershock sequences of the magnitude M _{W =}6.4 Bingöl earthquake of 1 May, 2003 (Turkey) are studied to analyze the spatial and temporal variability of seismicity parameters of the b value of the frequency-magnitude distribution and the p value describing the temporal decay rate of aftershocks. The catalog taken from the KOERI contains 516 events and one month’s time interval. The b value is found as 1.49 ± 0.07 with Mc =3.2. Considering the error limits, b value is very close to the maximum b value stated in the literature. This larger value may be caused by the paucity of the larger aftershocks with magnitude M _D ≥ 5.0. Also, the aftershock area is divided into four parts in order to detect the differences in b value and the changes illustrate the heterogeneity of the aftershock region. The p value is calculated as 0.86 ± 0.11, relatively small. This small p value may be a result of the slow decay rate of the aftershock activity and the small number of aftershocks. For the fitting of a suitable model and estimation of correct values of decay parameters, the sequence is also modeled as a background seismicty rate model. Constant background activity does not appear to be important during the first month of the Bingöl aftershock sequences and this result is coherent with an average estimation of pre-existing seismicity. The results show that usage of simple modified Omori law is reasonable for the analysis. The spatial variability in b value is between 1.2 and 1.8 and p value varies from 0.6 to 1.2. Although the physical interpretation of the spatial variability of these seismicity parameters is not straightforward, the variation of b and p values can be related to the stress and slip distribution after the mainshock, respectively. The lower b values are observed in the high stress regions and to a certain extent, the largest b values are related to Holocene alluvium. The larger p values are found in some part of the aftershock area although no slip occurred after the main shock and it is interpreted that this situation may be caused by the alluvium structure of the region. These results indicate that the spatial distribution in b and p values are generally related to the rupture mechanism and material properties of an aftershock area.     

Focal Fault and Stress Field of the Yao‘ an M 6.5 Earthquake on January15,2001

Using the 78 focal mechanism solutions of the fore shocks, main shock and after shocks of the earthquake sequence for the Yao'an earthquake, the characteristics of the focal faults and stress field for the earthquake sequence are analyzed. The results show that the main rupture plane of the Yao'an earthquake sequence is a tectonic fault with N50°W strike and steep dip and all the main shock, the fore shocks and the vast majority of after shocks occurred on the main rupture plane. A tectonic fracture with NNE-NE strike also participated in development process of the sequence. The focal stress field of the sequence dominated by principal compressional stress with nearly horizontal orientation SSE is consistent with the regional tectonic stress field. In the sequence development, the stress field in the focal region was complex with multi-azimuths and multi-action models and the focal rupture showed complex features with multi-directions and multi-patterns.     

Stress Release and Source Scaling of the 2010–2011 Canterbury,New Zealand Earthquake Sequence from Spectral Inversion of Ground Motion Data

The Canterbury earthquake sequence beginning with the 2010 M _W 7.2 Darfield earthquake is one of the most notable and well-recorded crustal earthquake sequences in a low-strain-rate region worldwide and as such provides a unique opportunity to better understand earthquake source physics and ground motion generation in such a tectonic setting. Ground motions during this sequence ranged up to extreme values of 2.2 g, recorded during the February 2011 M _W 6.2 event beneath the city of Christchurch. A better understanding of the seismic source signature of this sequence, in particular the stress release and its scaling with earthquake size, is crucial for future ground motion prediction and hazard assessment in Canterbury, but also of high interest for other low-to-moderate seismicity regions where high-quality records of large earthquakes are lacking. Here we present a source parameter study of more than 200 events of the Canterbury sequence, covering the magnitude range M _W 3–7.2. Source spectra were derived using a generalized spectral inversion technique and found to be well characterized by the ω ^?2 source model. We find that stress drops range between 1 and 20 MPa with a median value of 5 MPa, which is a factor of 5 larger than the median stress drop previously estimated with the same method for crustal earthquakes in much more seismically active Japan. Stress drop scaling with earthquake size is nearly self-similar, and we identify lateral variations throughout Canterbury, in particular high stress drops at the fault edges of the two major events, the M _W 7.2 Darfield and M _W 6.2 Christchurch earthquakes.     

Geochemistry of the volcano-hydrothermal system of El Chichón Volcano, Chiapas, Mexico

 The 1982 eruption of El Chichón volcano ejected more than 1 km³ of anhydrite-bearing trachyandesite pyroclastic material to form a new 1-km-wide and 300-m-deep crater and uncovered the upper 500 m of an active volcano-hydrothermal system. Instead of the weak boiling-point temperature fumaroles of the former lava dome, a vigorously boiling crater spring now discharges  / 20 kg/s of Cl-rich (∼15 000 mg/kg) and sulphur-poor ( / 200 mg/kg of SO₄), almost neutral (pH up to 6.7) water with an isotopic composition close to that of subduction-type magmatic water (δD=–15‰, δ¹⁸O=+6.5‰). This spring, as well as numerous Cl-free boiling springs discharging a mixture of meteoric water with fumarolic condensates, feed the crater lake, which, compared with values in 1983, is now much more diluted (∼3000 mg/kg of Cl vs 24 030 mg/kg), less acidic (pH=2.6 vs 0.56) and contains much lower amounts of S ( / 200 mg/kg of SO₄, vs 3550 mg/kg) with δ³⁴S=0.5–4.2‰ (+17‰ in 1983). Agua Caliente thermal waters, on the southeast slope of the volcano, have an outflow rate of approximately 100 kg/s of 71  °C Na–Ca–Cl water and are five times more concentrated than before the eruption (B. R. Molina, unpublished data). Relative N₂, Ar and He gas concentrations suggest extensional tectonics for the El Chichón volcanic centre. The ³He/⁴He and ⁴He/²⁰Ne ratios in gases from the crater fumaroles (7.3R_a, 2560) and Agua Caliente hot springs (5.3R_a, 44) indicate a strong magmatic contribution. However, relative concentrations of reactive species are typical of equilibrium in a two-phase boiling aquifer. Sulphur and C isotopic data indicate highly reducing conditions within the system, probably associated with the presence of buried vegetation resulting from the 1982 eruption. All Cl-rich waters at El Chichón have a common source. This water has the appearence of a "partially matured" magmatic fluid: condensed magmatic vapour neutralized by interaction with fresh volcaniclastic deposits and depleted in S due to anhydrite precipitation. Shallow ground waters emerging around the volcano from the thick cover of fresh pumice deposits (Red waters) are Ca–SO₄–rich and have a negative oxygen isotopic shift, probably due to ongoing formation of clay at low temperatures. Received: 21 July 1997 / Accepted: 4 December 1997     

Flow Characteristics of Large-Scale Liquefaction-Slip of the Loess Strata in Shibei Tableland, Guyuan City, Induced by the 1920 Haiyuan M8 1/2 Earthquake

Based on the dynamic triaxial liquefaction test of the loess samples which are taken from Shibei tableland, Guyuan City, Ningxia, China, the characteristics of dynamic strain, dynamic stress and pore water pressure are studied under cyclic loading. Triaxial shear test is conducted immediately after the sample reaches liquefaction point. During the test, the property of the liquefied soil is analyzed through fluid mechanics method, whereby the fluidity of the liquefied soil is represented by apparent viscosity.The results show that the fluidity of liquefied loess changes from "shear thickening" to "shear thinning" as the shear force continues, and the fluidity of liquefied loess is closely related to its structure. In addition, in the process of forming a new stable state, the apparent viscosity and deviant stress change with axial strain in a similar approach. When the sample reaches its stable state, it meanwhile shows a relatively stable apparent viscosity. According to the fluid mechanics and the law of conservation of energy, the slip distance of the liquefied soil is estimated, and the results are in good agreement with the field investigation results.     

Insights into eruption dynamics from textural analysis: the case of the May, 2008, Chaitén eruption

The May, 2008, Chaitén (southern Chile) eruption was characterized by several explosive events, each associated with plumes which reached up to about 19?km above sea level on May 6. A study of the textural and physical features of the juvenile clasts erupted during the climactic phase of the 2008 eruption of Chaitén is presented. Pumice clasts show unimodal density distribution (main mode at 600?kg/m³), average vesicularity of about 69?%, a glassy groundmass with no microcrystals, and vesicles with dimension between ～1?μm and ～2?mm. They also show a unimodal vesicle size distribution with most frequent vesicle size in the range 0.05–0.08?mm and an estimated vesicle number density of 1.3?±?0.5?×?10⁵?mm^?3 related to a rapid nucleation event produced during the late phases of magma rise. This is confirmed by the absence of microcrystals that could otherwise have delayed vesicle formation and allowed the magma to maintain a low viscosity and a supersaturation in volatiles. Vesiculation and fragmentation were triggered by a sudden decompression of the melt associated with the opening of the volcanic conduit (～10?MPa?s^?1).     

Relocation of the October 9, 1505 Earthquake with M63/4

There have been 12 possible locations of the earthquakes occurring in the South Yellow Sea since 1505.In this paper,the location of the earthquake that occurred in 1505 has been determined by the collection of more historical data of the influenced field,referring the isoseismal data of earthquakes with the epicentral intensityⅨ and combining geophysical field data with tectonic condition,due to its great influence on seismic safety assessment of some significant engineering.     

Pyroclastic flow deposits of the 1991 eruption of Volcán de Colima,Mexico

The April 16, 1991, eruption of Volcán de Colima represents a classical example of partial dome collapse with the generation of progressively longer-runout, Merapi-type pyroclastic flows that traveled up to 4 km along the El Cordoban gullies (East, Central and West). The flows filled the gullies with block-and-ash flow deposits up to 10 m thick, of which, after 7 years of erosion, only remnants remained in the El Cordoban West and East gullies. The El Cordoban Central gully, however, provided a well-preserved and incised longitudinal section of the 1991 deposits. The deposits were emplaced as proximal and distal facies, separated by a change in slope angle from >30° to <20°. The proximal facies consists of massive, clast-supported flow units (up to 1 m thick) with andesite blocks locally supported by a matrix of coarse ash and devoid of segregation structures or grading. The distal facies consists of a massive, matrix-supported deposit up to 8 m thick, which contains dispersed andesite blocks in a fine ash matrix. In the distal facies, a train of blocks marks flow-unit upper boundaries and, although sorting is poor, some grading is present. Thin, finely stratified, or dune-bedded layers of fine ash material are locally present above or below units of both facies. Sedimentologic parameters show that the size or fraction of large pyroclasts (larger than –1 ) decreases from proximal to distal facies, as the percentage of matrix (0 to 4 ) increases, especially immediately beyond the break in slope. We propose that the propagation of the Colima pyroclastic flows is critically dependent on local slope angle, the presence of erodible slope debris, and the decrease in grain size with distance from the vent. The progressive fining is probably caused by some combination of erosion, clast breakup and deposition of larger pyroclasts, and is itself influenced by the slope angle. In the proximal region, the flows moved as granular avalanches, in which interacting grains ground each other and erosion occurred to produce an overriding dilute ash cloud. The maximum runout distance of the avalanches was controlled by the angle of repose of the material, and the volume and grain size of source and eroded material. Because the slope angle is close to the repose angle for this debris, granular avalanches were not able to propagate far beyond the change in slope. If, however, an avalanche had enough mass in finer grain size fractions, at least part of the flow continued beyond the break in slope and across the volcano apron, propagating in a turbulent state and depositing surge layers, or in an otherwise settling-modified state and depositing block-and-ash flow layers.Editorial responsibility: T Druitt     

The Characteristic Analysis and Seismic Triggering Study of the M6.2 and M6.1 Dayao Earthquake Sequences in 2003

INTRODUCTIONTheM6·2 andM6·1 earthquakes occurred successivelyin Dayao,Yunnan Province on July 21and October 16 ,2003 and bothinduced human casualty and economic loss to some extent .Adequatehistorical information reveals that ,earthquakes bigger thanM6·0 in Yunnan are distributed mainlyalong some active tectonic faults and seismic belts . Dayao and Yao an are in the Dianzhong tectonicblock andthere are noapparentfaultsthere(Su Youjing,2004) .Seismic activitylevel inthisregioniswea…     

Interpretation of the Spatial Distribution Characteristics of the Co-seismic Landslides Induced by the 1920 Haiyuan M8 1/2 Earthquake Using Remote Sensing Images

Analyzing the spatial distribution characteristics of earthquake-induced secondary disasters based on advanced techniques is significantly important, especially in understanding the process of strong earthquakes in the Loess Pateau. Using ArcGIS, this study interprets multi-temporal high-resolution satellite images, field investigation data, and historical seismic records. Major conclusions are obtained as follows:① Landslides induced by the Haiyuan earthquake are mainly distributed in the intersection area of the end of the Haiyuan fault and Liupanshan fault, as indicated by multiple dense distribution centers; ② The landslide distribution of the Haiyuan Earthquake is determined by the distance to the fault, topographic relief, slope, lithology, and other factors. In detail, the closer the distance to the fault, the greater the density of the landslide. The greater the slope and relief of the terrain, the greater the density and the smaller the average area of a landslide. Compared with tertiary strata, Quaternary strata has a larger average area, and the density of the landslides is smaller; ③ The density curve of the death toll in the Haiyuan earthquake can be used as a reference for the distribution of co-seismic landslides. Several Haiyuan co-seismic landslides are distributed in the Tongwei landslide area; however, the major landslides here are induced by the 1718 Tongwei earthquake rather than the 1920 Haiyuan earthquake; ④ The co-seismic landslides of the Haiyuan earthquake exhibits the "slope effect" in the south-west plate of Haiyuan fault, presenting the dominant sliding direction towards the fault and epicenter; however, the "slope effect" is not evident in the northeast plate of the fault.     

Source location of long period seismicity at Volcàn de Colima,México

This paper presents an analysis of seismicity associated with the volcanic activity of Volcàn de Colima (México) and recorded in the period November 2005–April 2006 during a field survey by the Istituto Nazionale di Geofisica e Vulcanologia (INGV)–Osservatorio Vesuviano, the Observatorio Vulcanologico de Colima of Colima University and the Instituto Andaluz de Geofisica, University of Granada. Three different types of volcanic earthquakes have been identified on the basis of their spectral properties: Type A (0.3–1 Hz), Type B (1–5 Hz) and Type C (3–4 Hz). Results of polarization analysis applied to Type A events show a predominance of radial motion, indicating that the wavefield comprises compressional waves (P) and shear waves polarized in the vertical plane (SV), while the signal always begins with a negative polarity. Type A, B and C earthquakes have been located using both a flat layered model and a 3D model including topography. Hypocentre distributions indicate that the source of Type A signals is very shallow and confined to a small volume lying about 1 km below the crater. In contrast, the source of Type B and C events is significantly deeper, with most hypocentres located in a volume of about 1 km³ centred at 2.5–3 km depth. A cluster analysis based on the cross-correlation among the waveforms of different events recorded at the same station was applied to Type A earthquakes. Only two clusters, which include only a small percentage of events were found, indicating that earthquake families were uncommon during the period of our survey.     

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     

Petrology, magnetostratigraphy and geochronology of the Miocene volcaniclastic Tepoztlán Formation: implications for the initiation of the Transmexican Volcanic Belt (Central Mexico)

The volcaniclastic Tepoztlán Formation (TF) represents an important rock record to unravel the early evolution of the Transmexican Volcanic Belt (TMVB). Here, a depositional model together with a chronostratigraphy of this Formation is presented, based on detailed field observations together with new geochronological, paleomagnetic, and petrological data. The TF consists predominantly of deposits from pyroclastic density currents and extensive epiclastic products such as tuffaceous sandstones, conglomerates and breccias, originating from fluvial and mass flow processes, respectively. Within these sediments fall deposits and lavas are sparsely intercalated. The clastic material is almost exclusively of volcanic origin, ranging in composition from andesite to rhyolite. Thick gravity-driven deposits and large-scale alluvial fan environments document the buildup of steep volcanic edifices. K-Ar and Ar-Ar dates, in addition to eight magnetostratigraphic sections and lithological correlations served to construct a chronostratigraphy for the entire Tepoztlán Formation. Correlation of the 577 m composite magnetostratigraphic section with the Cande and Kent (1995) Geomagnetic Polarity Time Scale (GPTS) suggests that this section represents the time intervall 22.8–18.8 Ma (6Bn.1n-5Er; Aquitanian-Burdigalian, Lower Miocene). This correlation implies a deposition of the TF predating the extensive effusive activity in the TMVB at 12 Ma and is therefore interpreted to represent its initial phase with predominantly explosive activity. Additionally, three subdivisions of the TF were established, according to the dominant mode of deposition: (1) the fluvial dominated Malinalco Member (22.8–22.2 Ma), (2) the volcanic dominated San Andrés Member (22.2–21.3 Ma) and (3) the mass flow dominated Tepozteco Member (21.3–18.8 Ma).