Fluid sources for the La Guitarra epithermal deposit (Temascaltepec district, Mexico): Volatile and helium isotope analyses in fluid inclusions

The La Guitarra deposit (Temascaltepec district, South-Central Mexico), belongs to the low/intermediate sulfidation epithermal type, has a polymetallic character although it is currently being mined for Ag and Au. The mineralization shows a polyphasic character and formed through several stages and sub-stages (named I, IIA, IIB, IIC, IID, and III). The previous structural, mineralogical, fluid inclusion and stable isotope studies were used to constrain the selection of samples for volatile and helium isotope analyses portrayed in this study. The N₂/Ar overall range obtained from analytical runs on fluid inclusion volatiles, by means of Quadrupole Mass Spectrometry (QMS), is 0 to 2526, and it ranges 0 to 2526 for stage I, 0 to 1264 for stage IIA, 0 to 1369 for stage IIB, 11 to 2401 for stage IIC, 19 to 324 for stage IID, and 0 to 2526 for stage III. These values, combined with the CO₂/CH₄ ratios, and N₂-He-Ar and N₂-CH₄-Ar relationships, suggest the occurrence of fluids from magmatic, crustal, and shallow meteoric sources in the forming epithermal vein deposit. The helium isotope analyses, obtained by means of Noble Gas Mass Spectrometry, display R/Ra average values between 0.5 and 2, pointing to the occurrence of mantle-derived helium that was relatively diluted or “contaminated” by crustal helium. These volatile analyses, when correlated with the stable isotope data from previous works and He isotope data, show the same distribution of data concerning sources for mineralizing fluids, especially those corresponding to magmatic and crustal sources. Thus, the overall geochemical data from mineralizing fluids are revealed as intrinsically consistent when compared to each other.The three main sources for mineralizing fluids (magmatic, and both deep and shallow meteoric fluids) are accountable at any scale, from stages of mineralization down to specific mineral associations. The volatile and helium isotope data obtained in this paper suggest that the precious metal-bearing mineral associations formed after hydrothermal pulses of predominantly oxidized magmatic fluids, and thus it is likely that precious metals were carried by fluids with such origin. Minerals from base-metal sulfide associations record both crustal and magmatic sources for mineralizing fluids, thus suggesting that base metals could be derived from deep leaching of crustal rocks. At the La Guitarra epithermal deposit there is no evidence for an evolution of mineralizing fluids towards any dominant source. Rather than that, volatile analyses in fluid inclusions suggest that this deposit formed as a pulsing hydrothermal system where each pulse or set of pulses accounts for different compositions of mineralizing fluids.The positive correlation between the relative content of magmatic fluids (high N₂/Ar ratios) and H₂S suggests that the necessary sulfur to carry mostly gold as bisulfide complexes came essentially from magmatic sources. Chlorine necessary to carry silver and base metals was found to be abundant in inclusion fluids and although there is no evidence about its source, it is plausible that it may come from magmatic sources as well.     

Aftershocks series monitoring of the September 18, 2004 M = 4.6 earthquake at the western Pyrenees: A case of reservoir-triggered seismicity?

On September 18, 2004, a 4.6 mbLg earthquake was widely felt in the region around Pamplona, at the western Pyrenees. Preliminary locations reported an epicenter less than 20 km ESE of Pamplona and close to the Itoiz reservoir, which started impounding in January 2004. The area apparently lacks of significant seismic activity in recent times. After the main shock, which was preceded by series of foreshocks reaching magnitudes of 3.3 mbLg, a dense temporal network of 13 seismic stations was deployed there to monitor the aftershocks series and to constrain the hypocentral pattern. Aftershock determinations obtained with a double-difference algorithm define a narrow epicentral zone of less than 10 km², ESE–WNW oriented. The events are mainly concentrated between 3 and 9 km depth. Focal solutions were computed for the main event and 12 aftershocks including the highest secondary one of 3.8 mbLg. They show mainly normal faulting with some strike-slip component and one of the nodal planes oriented NW–SE and dipping to the NE. Cross-correlation techniques applied to detect and associate events with similar waveforms, provided up to 33 families relating the 67% of the 326 relocated aftershocks. Families show event clusters grouped by periods and migrating from NW to SE. Interestingly, the narrow epicentral zone inferred here is located less than 4 km away from the 111-m high Itoiz dam. These hypocentral results, and the correlation observed between fluctuations of the reservoir water level and the seismic activity, favour the explanation of this foreshock–aftershock series as a rapid response case of reservoir-triggered seismicity, burst by the first impoundment of the Itoiz reservoir. The region is folded and affected by shallow dipping thrusts, and the Itoiz reservoir is located on the hangingwall of a low angle southward verging thrust, which might be a case sensible to water level fluctuations. However, continued seismic monitoring in the coming years is mandatory in this area to infer more reliable seismotectonic and hazard assessments.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Evidence for a seismic activity mainly constituted of hybrid events at Cayambe volcano,Ecuador. Interpretation in a iced-domes volcano context

The high activity level of Hybrid Events (HE) detected beneath the Cayambe volcano since 1989 has been more thoroughly investigated with data from a temporary array. The unusual HE spectral content allows separating a high-frequency signal riding on a low-frequency one, with a probable single source. HEs are interpreted as high frequency VT events, produced by the interaction between magmatic heat and an underground water system fed by thaw water from the summital glacier, which trigger simultaneous low-frequency fluid resonance in the highly fractured adjacent medium. Pure VTs are interpreted as ‘aborted’ HEs occurring probably in the oldest and coldest part of the volcano complex. To cite this article: B. Guillier, J.-L. Chatelain, C. R. Geoscience 338 (2006).     

Great earthquakes of variable magnitude at the Cascadia subduction zone

Comparison of histories of great earthquakes and accompanying tsunamis at eight coastal sites suggests plate-boundary ruptures of varying length, implying great earthquakes of variable magnitude at the Cascadia subduction zone. Inference of rupture length relies on degree of overlap on radiocarbon age ranges for earthquakes and tsunamis, and relative amounts of coseismic subsidence and heights of tsunamis. Written records of a tsunami in Japan provide the most conclusive evidence for rupture of much of the plate boundary during the earthquake of 26 January 1700. Cascadia stratigraphic evidence dating from about 1600 cal yr B.P., similar to that for the 1700 earthquake, implies a similarly long rupture with substantial subsidence and a high tsunami. Correlations are consistent with other long ruptures about 1350 cal yr B.P., 2500 cal yr B.P., 3400 cal yr B.P., 3800 cal yr B.P., 4400 cal yr B.P., and 4900 cal yr B.P. A rupture about 700-1100 cal yr B.P. was limited to the northern and central parts of the subduction zone, and a northern rupture about 2900 cal yr B.P. may have been similarly limited. Times of probable short ruptures in southern Cascadia include about 1100 cal yr B.P., 1700 cal yr B.P., 3200 cal yr B.P., 4200 cal yr B.P., 4600 cal yr B.P., and 4700 cal yr B.P. Rupture patterns suggest that the plate boundary in northern Cascadia usually breaks in long ruptures during the greatest earthquakes. Ruptures in southernmost Cascadia vary in length and recurrence intervals more than ruptures in northern Cascadia.     

New Approach for Estimation of Static and Seismic Active Earth Pressure

To estimate static and seismic active earth pressure (P_ad) on a rigid retaining wall, numerical analyses using different step sizes have been carried out in this paper, based on the modified Culmann line method by considering Coulomb’s planar rupture surface. Equivalent pseudo-static seismic forces are considered in the analysis. A new concept of modified unit weight by considering ground surcharge is introduced under static and seismic conditions. By numerical analysis, area of soil (A) has been estimated to obtain the ratio of A/A₀ where A₀ is θh², θ is the angle between retaining structure and ground surface and h is the vertical height of the wall. This ratio remains constant for a particular type of soil and has been used to estimate the maximum active earth pressure using force diagram. Results are provided in tabular form for easy calculation of the coefficient of static and seismic active earth pressure. Present results by considering the new technique, compares well with the results obtained by earlier researchers.     

Centrifuge modelling-based verification of response spectra design methods for some vulnerable structures

Residential RC framed structures suffered heavily during the 2001 Bhuj earthquake in Gujarat, India. These types of structures also saw severe damage in other earthquakes such as the 1999 Kocaeli earthquake in Turkey and 921 Ji-Ji earthquake in Taiwan. In this paper the seismic response of residential structures was investigated using physical modelling. Idealised soft storey and top heavy, two degrees of freedom (2DOF) portal frame structures were developed and tested on saturated and dry sand models at 25 g using the Schofield Centre 10-m Beam Centrifuge. It was possible to recreate observed field behaviour using these models. As observed in many of the recent earthquakes, soft storey structures were found to be particularly vulnerable to seismic loads. Elastic response spectra methods are often used in the design of simple portal frame structures. The seismic risk of these structures can be significantly increased due to modifications such as removal of a column or addition of heavy water tanks on the roof. The experimental data from the dynamic centrifuge tests on such soft storey or top-heavy models was used to evaluate the predictions obtained from the response spectra. Response spectra were able to predict seismic response during small to moderate intensity earthquakes, but became inaccurate during strong earthquakes and when soil structure interaction effects became important. Re-evaluation of seismic risk of such modified structures is required and time domain analyses suggested by building codes such as IBC, UBC or NEHRP may be more appropriate.     

用现今小震推断洪洞、临汾两次历史大震的震源断层

山西临汾地区是一个历史强震多发区，1303年和1695年发生了洪洞（M＝8）和临汾（M＝73／4）两次特大地震，这两次地震所在区域至今仍在持续不断的小地震活动，具有明显的大震区地震长期活动特征，我们对临汾无线传输地震台网记录的1987－1999年期间发生的1670次中，小地震重新进行了震源定位，根据对这些地震震源位置三维空间分布特征和震源机解制的分析，认为洪洞地震的震源断层应是长80km，埋深5－26km的NNE走向，高倾角的右旋走滑型断层，而临汾地震的震源断层是长70km，埋深5－22km的NWW走向，高倾角的左旋走滑型断层。这与洪洞，临汾两次大震极震区的等震线及该地区应力场的构造环境是吻合的。     

Moment-tensor statistics of the 1997 Long Valley microearthquake swarm

Abstract. A simple, fast, moment-tensor inversion method using bandpass-filtered P-amplitudes was used to study the moment-tensor statistics of Long Valley caldera microearthquakes. The events were recorded in the summer of 1997, during a swarm in the caldera. The swarm was associated with geodetic extension, uplift, and subsequent moderate earthquake activity. The moment tensor solutions for 1,993 events were calculated using the new method. The majority of the resulting focal mechanisms appear to be explained in terms of double couple mechanisms. Since some events did exhibit considerable deviation from double-couples, the moment data were studied for their statistical significance. The moments of the actual data were compared to the moments of synthetic data with varying degrees of random noise in their spectra. The results of this study suggested that unless data from more than 20 stations are used and the earthquake epicenter is located inside or very close to the network area, moment-tensor inversion does not correctly resolve the non-double-couple components of microearthquakes. Analysis of the inversion residuals shows that the average noise in the P-wave spectra was close to 20%. The fluctuations of the volumetric components of the moment-tensor are in good agreement with those of the synthetic pure double-couples with 20% of added noise. Thus the moment-tensor statistics suggests that little if any volume change is required to explain the observed seismic energy release in the swarm. However, the statistics do show that a significant compensated-linear-vector-dipole component maybe present in the bulk of the seismicity. Given the network used in this study, such a component could not be precisely resolved for individual earthquakes. This possibility deserves further investigation because of its bearing on the nature of fluid-fault-earthquake processes in swarms.