A 3100-year history of argon isotopic and compositional variation at El Chichón volcano

El Chichón volcano has produced at least 11 eruptive events in the past 8000 years, all of which display a relatively constant trachyandesitic composition. The constancy of the eruptive products attests to the tapping of a long-lived magma chamber and suggests a system held in steady state between the influx of basaltic magma and differentiation of resident magma. We have sampled eruptive products from eight of these events, subjecting their dominant phenocryst phases (plagioclase and hornblende) to argon isotopic analysis in order to investigate the dynamics of a steady state magmatic system. Plagioclase from the older eruptions, 1500 yr BP, 1600 yr BP, 2000 yr BP and 3100 yr BP, display variable enrichment of ⁴⁰Ar (excess argon), whereas hornblende from the same eruptions shows little or no enrichment. In contrast, both plagioclase and hornblende phenocrysts from the younger eruptions in 1982, 550 yr BP, 900 yr BP and 1250 yr BP have argon isotopic ratios near-atmospheric values. Isochron analysis of these mineral phases rules out xenocrystic contamination as the source of excess argon, while mafic enclaves and isotopic and compositional zoning evidence frequent recharge events, suggesting the source of this argon is most likely the same as the source of this basaltic magma; correlation with Cl points to melt/vapor inclusions as the primary host for Ar within the phenocrysts. Argon isotopic variations point to a disturbance of the system between 1500 yr BP and 1250 yr BP, while compositional evidence for a major mafic input is present in the 900 yr BP eruption, indicating a complex relationship between recharge and isotopic signature of eruptive products. The amount of excess argon within a plagioclase phenocryst are a function of variations in melt and vapor inclusion abundance, time elapsed between melt/vapor inclusion entrapment and eruption, variations in Ar abundance in melt (itself a function of vapor exsolution timing) and time variations in Ar isotopic composition of the melt; subdued behavior of hornblende is due to slower diffusion and minimal inclusions.     

Some regularity in the process of re-awakening of andesite and dacite volcanoes: specific features of the 1982 El Chichón volcano,México reactivation

El Chichón volcano is an andesite stratovolcano in southern México. It erupted in March 1982, after about 550 years of quiescence. The 1982 eruption of El Chichón has not been followed by the growth of a lava dome within the newly formed crater. This is rather anomalous since the construction of a new dome after the destruction of an old one is a common process during the eruptions at andesite and dacite volcanoes. To discuss this anomalous aspect of the El Chichón eruption, some regularity in the process of re-awakening of dormant (here defined as a period of quiescence of more than 100 years) andesite and dacite volcanoes are studied based on the seismic activity recorded at the volcanoes Bezymianny, Mount St. Helens, El Chichón, Unzen, Pinatubo and Soufrière Hills. Three stages were identified in the re-awakening activity of these volcanoes: (1) preliminary seismic activity, leading up to the first phreatic explosion; (2) activity between the first and the largest explosions; (3) post-explosion dome-building process. The eruptions were divided into two groups: low-VEI (Volcanic Explosivity Index) and the long duration stage-1 events (Unzen, 1991 and Soufrière Hills volcano, 1995) and high-VEI and the short duration stage-1 events (Bezymianny, 1956; Mount St. Helens, 1980; El Chichón, 1982 and Pinatubo, 1992). The comparative analysis of the seismo-eruptive activity of two eruptions of the second group, the 1980 of Mt. St. Helens and the 1982 of El Chichón, produced an explanation the absence of new dome building during the 1982 eruption of El Chichón volcano. It may be explained in terms of the unusually rapid emission of gas and water from the magmatic and hydrothermal system beneath the volcano during a relatively short sequence of large explosions that could have sharply increased the viscosity of the magma making impossible its exit to the surface.     

On the origin of El Chichón volcano and subduction of Tehuantepec Ridge: A geodynamical perspective

The origin of El Chichón volcano is poorly understood, and we attempt in this study to demonstrate that the Tehuantepec Ridge (TR), a major tectonic discontinuity on the Cocos plate, plays a key role in determining the location of the volcano by enhancing the slab dehydration budget beneath it. Using marine magnetic anomalies we show that the upper mantle beneath TR undergoes strong serpentinization, carrying significant amounts of water into subduction. Another key aspect of the magnetic anomaly over southern Mexico is a long-wavelength (∼ 150 km) high amplitude (∼ 500 nT) magnetic anomaly located between the trench and the coast. Using a 2D joint magnetic-gravity forward model, constrained by the subduction P–T structure, slab geometry and seismicity, we find a highly magnetic and low-density source located at 40–80 km depth that we interpret as a partially serpentinized mantle wedge formed by fluids expelled from the subducting Cocos plate. Using phase diagrams for sediments, basalt and peridotite, and the thermal structure of the subduction zone beneath El Chichón we find that ∼ 40% of sediments and basalt dehydrate at depths corresponding with the location of the serpentinized mantle wedge, whereas the serpentinized root beneath TR strongly dehydrates (∼90%) at depths of 180-200 km comparable with the slab depths beneath El Chichón (200-220 km). We conclude that this strong deserpentinization pulse of mantle lithosphere beneath TR at great depths is responsible for the unusual location, singularity and, probably, the geochemically distinct signature (adakitic-like) of El Chichón volcano.     

Deposition temperature of some PDC deposits from the 1982 eruption of El Chichón volcano (Chiapas,Mexico) inferred from rock-magnetic data

Thermal remanent magnetization (TRM) analyses were carried out on lithic fragments from two different typologies of pyroclastic density current (PDC) deposits of the 1982 eruption of El Chichón volcano, in order to estimate their equilibrium temperature (T_dep) after deposition. The estimated T_dep range is 360–400 °C, which overlaps the direct measurements of temperature carried out four days after the eruption on the PDC deposits. This overlap demonstrates the reliability of the TRM method to estimate the T_dep of pyroclastic deposits and to approximate their depositional temperature. These results also constraint the time needed for reaching thermal equilibrium within four days for the studied PDC deposits, in agreement with predictions of theoretical models.     

Estimating thermal inflow to El Chichón crater lake using the energy-budget,chemical and isotope balance approaches

El Chichón crater lake appeared immediately after the 1982 catastrophic eruption in a newly formed, 1-km wide, explosive crater. During the first 2 years after the eruption the lake transformed from hot and ultra-acidic caused by dissolution of magmatic gases, to a warm and less acidic lake due to a rapid “magmatic-to-hydrothermal transition” — input of hydrothermal fluids and oxidation of H₂S to sulfate. Chemical composition of the lake water and other thermal fluids discharging in the crater, stable isotope composition (δD and δ¹⁸O) of lake water, gas condensates and thermal waters collected in 1995–2006 were used for the mass-balance calculations (Cl, SO₄ and isotopic composition) of the thermal flux from the crater floor. The calculated fluxes of thermal fluid by different mass-balance approaches become of the same order of magnitude as those derived from the energy-budget model if values of 1.9 and 2 mmol/mol are taken for the catchment coefficient and the average H₂S concentration in the hydrothermal vapors, respectively. The total heat power from the crater is estimated to be between 35 and 60 MW and the CO₂ flux is not higher than 150 t/day or ~ 200 gm^− 2 day^− 1.     

Hazard map of El Chichón volcano,Chiapas, México: Constraints posed by eruptive history and computer simulations

El Chichón volcano consists of a 2-km wide Somma crater compound cone 0.2 Ma old with peripheral domes with a central crater reactivated several times during the Holocene. The most recent eruption at El Chichón occurred from March 28 to April 4, 1982, resulting in the worst volcanic disaster during historical times in Mexico, killing more than 2000 people and destroying nine towns and small communities. The volcanic hazard map of El Chichón is based on detailed field work that documented twelve eruptions during the last 8000 years, and computer simulations. To validate the results, computer simulations were first performed over pre-1982 topography mimicking the extent of the actual deposits produced and afterwards run over post-1982 topography. These eruptions have produced pyroclastic fall, surge, flow and lahar deposits. Pyroclastic flows have different volumes and Heim coefficients varying from 0.2 (pumice flows), to 0.15 (block-and-ash flows) and 0.10 (ash flows). Simulations using FLOW3D and TITAN2D indicate that pumice flows and block-and-ash flows can fill the moat area and follow main ravines up to distances of ca. 3 km from the crater, with no effect on populations around the volcano. On the other hand, more mobile ash flows related to column-collapse events can reach up to 4 km from the vent, but will always follow the same paths and still not affect surrounding populations. The energy-cone model was used to simulate the outflow of pyroclastic surges based on the 1982 event (H/L = 0.1 and 0.2), and shows that surges may reach some towns around the volcano.     

Major and trace element geochemistry of neutral and acidic thermal springs at El Chichón volcano,Mexico: Implications for monitoring of the volcanic activity

Four groups of thermal springs with temperatures from 50 to 80 °C are located on the S–SW–W slopes of El Chichón volcano, a composite dome-tephra edifice, which exploded in 1982 with a 1 km wide, 160 m deep crater left. Very dynamic thermal activity inside the crater (variations in chemistry and migration of pools and fumaroles, drastic changes in the crater lake volume and chemistry) contrasts with the stable behavior of the flank hot springs during the time of observations (1974–2005). All known groups of hot springs are located on the contact of the basement and volcanic edifice, and only on the W–SW–S slopes of the volcano at almost same elevations 600–650 m asl and less than 3 km of direct distance from the crater. Three groups of near-neutral (pH ≈ 6) springs at SW–S slopes have the total thermal water outflow rate higher than 300 l/s and are similar in composition. The fourth and farthest group on the western slope discharges acidic (pH ≈ 2) saline (10 g/kg of Cl) water with a much lower outflow rate (< 10 l/s).     

Morphology of ash aggregates from wet pyroclastic surges of the 1982 eruption of El Chichón Volcano, Mexico

The detailed stratigraphic study of the pyroclastic surge units S1, IU, and S3 produced during the most violent phases of the 1982 eruption of El Chichón volcano, contains a complex succession of hydromagmatic events triggered by the interaction of different proportions of magma and external water. Component analyses of the horizons within single units reveal that almost all wet and cohesive horizons contain ash aggregates. Based on their morphology and internal structure four different types of aggregates were distinguished: (a) accretionary lapilli, (b) armored lapilli, (c) irregular aggregates, and (d) cylindrical aggregates. The first three types have been described in the volcanological literature (field and experimental studies); cylindrical forms are reported here for the first time. These hollow cylindrical aggregates consist of concentric layers of crystals and glass fragments set in a finer-grained matrix. They formed around millimeter-size foliage fragments that are locally preserved in the interior of the aggregates as scorched or completely carbonized vestiges. SEM analyses suggest different mechanisms of formation for the four types of aggregates. Irregular aggregates and armored lapilli formed nearly instantaneously, whereas accretionary lapilli and cylindrical aggregates resulted from progressive aggregation of ash in different regions of the eruptive cloud.All types of ash aggregates contain fractured particles. This common feature suggests that particles ruptured during fragmentation prior to the growth of the aggregates. Broken clasts with cracks filled by a fine-grained matrix only occur inside the cylindrical ash aggregates and to a lesser degree in some types of accretionary lapilli. This suggests that small thermal contrasts at the contact of warm particles with the colder fine-grained matrix of the aggregate cause existing small fractures to propagate and open as the already weakened clasts deform slightly. The occurrence of all four types of aggregates in some horizons indicates that several mechanisms of aggregation occurred nearly simultaneously. The pyroclastic clouds therefore were not only stratified in terms of density but the content of fluid phases also were not uniform. A dark-red, Fe-rich amorphous film (locally rich in P and S) envelops the particles and fosters their preservation in the deposits by forming a hard shell. The composition of this cement reflects the abundance of these elements in acid fluids of hydrothermal systems that were intersected by the conduit during the eruption. In distal areas, fallout aggregates were incorporated by dissipating pyroclastic surges.     

Chemical composition of fumarolic gases and spring discharges from El Chichòn volcano, Mexico: causes and implications of the changes detected over the period 1998–2000

Since the March–April 1982 eruption of El Chichòn volcano, intense hydrothermal activity has characterised the 1-km-wide summit crater. This mainly consists of mud and boiling pools, fumaroles, which are mainly located in the northwestern bank of the crater lake. During the period 1998–2000, hot springs and fumaroles discharging inside the crater and from the southeastern outer flank (Agua Caliente) were collected for chemical analyses. The observed chemical fluctuations suggest that the physico-chemical boundary conditions regulating the thermodynamic equilibria of the deep rock/fluid interactions have changed with time. The chemical composition of the lake water, characterised in the period 1983–1997 by high Na⁺, Cl⁻, Ca²⁺ and SO₄²⁻ contents, experienced a dramatic change in 1998–1999, turning from a Na⁺–Cl⁻- to a Ca²⁺–SO₄²⁻-rich composition. In June 2000, a relatively sharp increase in Na⁺ and Cl⁻ contents was observed. At the same time, SO₂/H₂S ratios and H₂ and CO contents in most gas discharges increased with respect to the previous two years of observations, suggesting either a new input of deep-seated fluids or local variations of the more surficial hydrothermal system. Migration of gas manifestations, enhanced number of emission spots and variations in both gas discharge flux and outlet temperatures of the main fluid manifestations were also recorded. The magmatic-hydrothermal system of El Chichòn is probably related to interaction processes between a deep magmatic source and a surficial cold aquifer; an important role may also be played by the interaction of the deep fluids with the volcanic rocks and the sedimentary (limestone and evaporites) basement. The chemical and physical changes recorded in 1998–2000 were possibly due to variations in the permeability of the conduit system feeding the fluid discharges at surface, as testified by the migration of gas and water emanations. Two different scenarios can be put forward for the volcanic evolution of El Chichòn: (1) build-up of an infra-crater dome that may imply a future eruption in terms of tens to hundreds of years; (2) minor phreatic–phreatomagmatic events whose prediction and timing is more difficult to constrain. This suggests that, unlike the diminished volcanic activity at El Chichòn after the 1982 paroxistic event, the volcano-hydrothermal fluid discharges need to be more constantly monitored with regular and more frequent geochemical sampling and, at the same time, a permanent network of seismic stations should be installed.     

Eruption episodes and magma recharge events in andesitic systems: Mt Taranaki,New Zealand

Eruption episodes, where a series of eruption events are generically related, can include the eruption of a wide spectrum of volcanic activity over decadal periods. This paper concentrates on the opening phases of an eruption episode which occurred approximately 1800 yrs BP from Mt Taranaki, New Zealand. These events spanned the eruption of differing bulk compositions and styles from two distinct vent locations; an andesitic sub-plinian eruption from the summit vent and a scoria cone-building eruption of basaltic magma from a satellite vent. Compositional profiles and zoning textures of plagioclase, amphibole and clinopyroxene phenocrysts from the opening andesitic event show evidence of magma mixing and subsequent crystallisation just prior to the initiation of the eruption episode. Titanomagnetite grain morphology and Ti variation suggest that the magma mixing event occurred within a few days to weeks before the eruption acting as a trigger for it. We present a magmatic model which is constrained by the petrological observations and eruptions of the episode. In this model magma differentiation at depth causes its rise and recharging of a mid-crustal magma storage area at 5–7 km. Although the recharging magma differed slightly in oxygen fugacity and temperature, it was compositionally and physically similar enough to the residing andesitic magma to allow efficient mixing. The petrological characteristics described here can be readily observed and enable identification of mixing events in other recent eruption episodes.     

Phase equilibria and metamorphic evolution of garnet‐mica‐plagioclase gneisses from the Qiliping terrane in the western Dabie Orogen,central China

The extensive gneisses in the high‐pressure and ultrahigh‐pressure metamorphic terrane in the Dabie‐Sulu orogen usually show no evidence of eclogite‐facies metamorphism. The garnet‐mica‐plagioclase gneisses from the Qiliping region in the western Dabie Orogen, comprise garnet, phengite, biotite, plagioclase, quartz, rutile, ilmenite, chlorite, epidote, and hornblende. The garnet porphyroblasts, with inclusions of quartz, epidote, and rutile, exhibit slight compositional zonations, from core to mantle with an increase in pyrope and a decrease in spessartine, and from mantle to rim with a decrease in pyrope and grossular and an increase in spessartine. The high‐Si phengite indicates that the gneisses may be subjected to a high‐pressure metamorphism. By the P–T pseudosections calculated in a system NCKMnFMASHTO (Na₂O‐CaO‐K₂O‐MnO‐FeO‐MgO‐Al₂O₃‐SiO₂‐H₂O‐TiO₂‐O) for two representative samples, the metamorphic P–T path, reconstructed by the compositionally zoned garnet, shows that the prograde metamorphism is characterized by a temperature increase with a slight pressure increase from the conditions of 17.6 ± 1.5 kbar at 496 ± 15°C to the peak‐pressure ones of 21.8 ± 1.5–22.7 ± 1.5 kbar at 555 ± 15–561 ± 15°C; the early retrograde stage is dominated by decompression with a temperature increase to the maximum of 608 ± 15–611 ± 18°C at 10.3 ± 1.5–11.0 ± 1.5 kbar; and the late retrograde one is predominated by pressure and temperature decreases. The mineral assemblages in the prograde metamorphism are predicted to contain garnet, glaucophane, jadeite, lawsonite, phengite, quartz, rutile, and/or chlorite, which is different from those observed at present. Such high‐pressure metamorphism can partly be reconstructed by the P–T pseudosection in combination with the high‐Si phengite and garnet compositions in the core and mantle. This provides an important constraint on the subduction and exhumation of the terrane during the continent–continent collision between the Yangtze and Sino‐Korean cratons.     

Crustal structure in the Falcón Basin area, northwestern Venezuela, from seismic and gravimetric evidence

The Falcón Basin in northwestern Venezuela has a complex geological history driven by the interactions between the South American and Caribbean plates. Igneous intrusive bodies that outcrop along the axis of the basin have been associated with crustal thinning, and gravity modeling has shown evidence for a significantly thinned crust beneath the basin. In this study, crustal scale seismic refraction/wide-angle reflection data derived from onshore/offshore active seismic experiments are interpreted and forward-modeled to generate a P-wave velocity model for a 450 km long profile. The final model shows thinning of the crust beneath the Falcón Basin where depth to Moho decreases to 27 km from a value of 40 km about 100 km to the south. A deeper reflected phase on the offshore section is interpreted to be derived from the downgoing Caribbean slab. Velocity values were converted to density and the resulting gravimetric response was shown to be consistent with the regional gravity anomaly. The crustal thinning proposed here supports a rift origin for the Falcón Basin.     

Rate and historical change of direct recharge from precipitation constrained by unsaturated zone profiles of chloride and oxygen‐18 in dry river bed of North China Plain

Quantifying of direct recharge derived from precipitation is crucial for assessing sustainability of well‐irrigated agriculture. In the North China Plain, the land use is dominated by groundwater‐irrigated farmland where the direct recharge derived from precipitation and irrigation. To characterize the mean rate and historical variance of direct recharge derived from precipitation, unsaturated zone profiles of chloride and δ¹⁸O in the dry river bed of the Beiyishui River were employed. The results show that archival time scale of the profile covers the duration from 1980 to 2002 (corresponding to depths from 5 to 2 m) which is indicated by matching the δ¹⁸O peaks in the isotope profile with the aridity indexes gained by instrumental records of annual precipitation and annual potential evaporation. Using the chloride mass balance method, the mean rate of the direct recharge corresponding to the archival time scale is estimated to be 3·8 ± 0·8 mm year^?1, which accounts for about 0·7% of the long‐term average annual precipitation. Further, the direct recharge rates vary from 2·1 to 6·8 mm year^?1 since 1980. Despite the subhumid climate, the estimate of recharge rates is in line with other findings in semiarid regions. The low rate of direct recharge is considered as a result of the relative dry climate in recent decades. In dry river bed, unsaturated zone profiles of chloride and δ¹⁸O combined with instrumental records could offer valuable information about the direct recharge derived from precipitation during droughts. Copyright © 2011 John Wiley & Sons, Ltd.     

Groundwater recharge in the Akaki catchment,central Ethiopia: evidence from environmental isotopes (δ18O, δ2H and 3H) and chloride mass balance

Recharge patterns, possible flow paths and the relative age of groundwater in the Akaki catchment in central Ethiopia have been investigated using stable environmental isotopes δ¹⁸O and δ²H and radioactive tritium (³H) coupled with conservative chloride measurements. Stable isotopic signatures are encoded in the groundwater solely from summer rainfall. Thus, groundwater recharge occurs predominantly in the summer months from late June to early September during the major Ethiopian rainy season. Winter recharge is lost through high evaporation–evapotranspiration within the unsaturated zone after relatively long dry periods of high accumulated soil moisture deficits. Chloride mass balance coupled with the isotope results demonstrates the presence of both preferential and piston flow groundwater recharge mechanisms. The stable and radioactive isotope measurements further revealed that groundwater in the Akaki catchment is found to be compartmentalized into zones. Groundwater mixing following the flow paths and topography is complicated by the lithologic complexity. An uncommon, highly depleted stable isotope and zero‐³H groundwater, observed in a nearly east–west stretch through the central sector of the catchment, is coincident with the Filwoha Fault zone. Here, deep circulating meteoric water has lost its isotopic content through exchange reactions with CO₂ originating at deeper sources or it has been recharged with precipitation from a different rainfall regime with a depleted isotopic content. Copyright © 2006 John Wiley & Sons, Ltd.     

Suspended sediment yield and metal contamination in a river catchment affected by El Niño events and gold mining activities: the Puyango river basin,southern Ecuador

The suspended sediment yield and the transfer of polluted sediment are investigated for the Puyango river basin in southern Ecuador. This river system receives metal (Cd, Cu, Hg, Pb and Zn) and cyanide pollution generated by mining, and is associated with large‐scale hydrological variability, which is partly governed by El Niño events. Field sampling and statistical modelling methods are used to quantify the amount of mine tailings that is discharged into the basin. Annual suspended sediment yields are estimated using a novel combination of the suspended sediment rating method and Monte Carlo simulations, which allow for propagation of the uncertainties of the calculations that lead to final load estimates. Geochemical analysis of suspended and river bed sediment is used to assess the dispersion and long‐term fate of contaminated sediment within the river catchment. Knowledge of the inter‐ and intra‐annual variation in suspended sediment yield is shown to be crucial for judging the importance of mining discharges, and the extent to which the resultant pollution is diluted by river flows. In wet years, polluted sediments represent only a very small proportion of the yield estimates, but in dry years the proportion can be significant. Evidence shows that metal contaminated sediments are stored in the Puyango river bed during low flows. Large flood events flush this sediment periodically, both on an annual cycle associated with the rainy season, and also related to El Niño events. Therefore, environmental impacts of mining‐related discharges are more likely to be severe during dry years compared with wet years, and in the dry season rather than the wet season. The hydrological consequences of El Niño events are shown to depend upon the extent to which these events penetrate inland. It is, thus, shown that the general conclusion that El Niño events can significantly affect suspended sediment yields needs evaluation with respect to the particular way in which those events affect a given catchment. Copyright © 2003 John Wiley & Sons, Ltd.