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.     

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.     

Hydrochemical dynamics of the “lake–spring” system in the crater of El Chichón volcano (Chiapas,Mexico)

El Chichón volcano (Chiapas, Mexico) erupted violently in March–April 1982, breaching through the former volcano–hydrothermal system. Since then, the 1982 crater has hosted a shallow (1–3.3 m, acidic (pH ∼ 2.2) and warm (∼ 30 °C) crater lake with a strongly varying chemistry (Cl/SO₄ = 0–79 molar ratio). The changes in crater lake chemistry and volume are not systematically related to the seasonal variation of rainfall, but rather to the activity of near-neutral geyser-like springs in the crater (Soap Pool). These Soap Pool springs are the only sources of Cl for the lake. Their geyser-like behaviour with a long-term (months to years) periodicity is due to a specific geometry of the shallow boiling aquifer beneath the lake, which is the remnant of the 1983 Cl-rich (24,000 mg/l) crater lake water. The Soap Pool springs decreased in Cl content over time. The zero-time extrapolation (1982, year of the eruption) approaches the Cl content in the initial crater lake, meanwhile the extrapolation towards the future indicates a zero-Cl content by 2009 ± 1. This particular situation offers the opportunity to calculate mass balance and Cl budget to quantify the lake–spring system in the El Chichón crater. These calculations show that the water balance without the input of SP springs is negative, implying that the lake should disappear during the dry season. The isotopic composition of lake waters (δD and δ¹⁸O) coincide with this crater lake-SP dynamics, reflecting evaporation processes and mixing with SP geyser and meteoric water. Future dome growth, not observed yet in the post-1982 El Chichón crater, may be anticipated by changes in lake chemistry and dynamics.     

Primary igneous anhydrite: Progress since its recognition in the 1982 El Chichón trachyandesite

Primary igneous anhydrite was first identified in 1982 El Chichón pumices. Analysis of the sulfur budget for the eruption provided compelling evidence that the pre-eruptive magma contained a significant gas phase at ∼ 7 km depth in order to account for the “excess gas release” of ∼ 5–9 million tons of SO₂ to the stratosphere by the eruption. Primary igneous anhydrite and a larger “excess gas release” of ∼ 20 million tons of SO₂ were noted for the significantly larger eruption of Mount Pinatubo in 1991, for which a separate gas phase at ∼ 7–9 km depth was also required by the sulfur budget. Pumices from both eruptions have mineral assemblages dominated by plagioclase and hornblende, with minor biotite, and show evidence for co-nucleation and mutual inclusions of anhydrite and apatite. Both magmas were also very water-rich and highly oxidized, with oxygen fugacities $1 log unit above the synthetic Ni–NiO buffer. Furthering the similarities between these two eruptions, ion-microprobe analyses of sulfur isotopic compositions of anhydrites in pumices from El Chichón and Mount Pinatubo both showed that individual crystals are isotopically homogeneous, but inter-crystalline variations in δ³⁴S are well beyond analytical error.     

Constraining P-wave velocity variations in the upper mantle beneath Southeast Asia

We have produced a P-wave model of the upper mantle beneath Southeast (SE) Asia from reprocessed short period International Seismological Centre (ISC) P and pP data, short period P data of the Annual Bulletin of Chinese Earthquakes (ABCE), and long period PP-P data. We used 3D sensitivity kernels to combine the datasets, and mantle structure was parameterized with an irregular grid. In the best-sampled region our data resolve structure on scale lengths less than 150 km. The smearing of crustal anomalies to larger depths is reduced by a crustal correction using an a priori 3D model. Our tomographic inversions reveal high-velocity roots beneath the Archean Ordos Plateau, the Sichuan Basin, and other continental blocks in SE Asia. Beneath the Himalayan Block we detect high seismic velocities, which we associate with subduction of Indian lithospheric mantle. This structure is visible above the 410 km discontinuity and may not connect to the remnant of the Neo-Tethys oceanic slab in the lower mantle. Our images suggest that only the southwestern part of the Tibetan plateau is underlain by Indian lithosphere and, thus, that the upper mantle beneath northeastern Tibet is primarily of Asian origin. Our imaging also reveals a large-scale high-velocity structure in the transition zone beneath the Yangtze Craton, which could have been produced in multiple subduction episodes. The low P-wave velocities beneath the Hainan Island are most prominent in the upper mantle and transition zone; they may represent counter flow from the surrounding subduction zones, and may not be unrelated to processes beneath eastern Tibet.     

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).     

Long term AVHRR observations of surface radiative flux from El Chichón crater lake (1996–2006)

It has been shown that due to the small surface of crater lakes, temperature surveillance is a problem using meteorological satellites. This is particularly true for El Chichón surface lake because it's about one tenth of an AVHRR pixel at nadir. In order to guarantee at least one unmixed pixel in AVHRR data, it is necessary to use only AVHRR data from NOAA satellite passes as close as possible to the nadir for the period 1996–2006, therefore AVHRR data of El Chichón's crater lake were only used it they were cloudless and had scan angles close to nadir. The analysis of the time series data shows that lake surface temperature had annual maximum values (> 35 °C) during 1996 and 1997 then surface temperature decay with a negative exponential trend reaching a steady state of about 30 °C in the last years (2004–2006). A seasonal temperature variation between the dry (December to May) and the wet (June to November) seasons is also observed. Differences between nocturnal and midday temperatures indicate the influence of lake energy emission (including reflectance) at midday under a strong short-wave solar radiation. Water surface radiative flux under these conditions reaches an average of 77.8 W m^− 2 and a maximum of 187.1 W m^− 2. Whereas nocturnal heat output from El Chichón crater lake has an average surface radiative flux of 20.4 W m^− 2 and a maximum of 74.3 W m^− 2.     

Flat-Slab Thermal Structure and Evolution Beneath Central Mexico

Recent seismic and magnetotelluric experiments, aimed at better characterizing the shape and state of the subducting slab and continental crust beneath Central Mexico, exposed significant differences with conclusions of previous studies. A new slab geometry is revealed in which the subducting Cocos slab is perfectly flat between 120 to 290?km from the trench, after which it plunges into the asthenosphere at a dip angle of ~65°, in sharp contrast with the previously proposed ~20° dip angle. Seismic tomography studies show negative P-wave velocity anomalies (?2 to ?4%) in the mantle wedge beneath the Mexican Volcanic Belt, and positive anomalies (+2 to +3%) for the subducted Cocos slab. Magnetotelluric experiments exposed a very low-resistivity area (1?C10? ??m) located within the continental crust just below the Mexican Volcanic Arc. Finally, several spots of non-volcanic tremors (NVTs) have been recorded inside the continental crust above the flat-slab segment. While all these experiments provide a better picture of the subduction system beneath Central Mexico, several key processes need further investigation. In this study, we take advantage of these new observations to better constrain the thermal structure beneath Central Mexico. Two different thermal models are computed for a mantle potential temperature (T _p) of 1,350 and 1,450°C, respectively. The new thermal structures are then converted into P-wave velocity anomalies and compared with the observed V _p anomalies. We found that a T _p of 1,450°C produced larger V _p anomalies that do not fit the observations. However, using a T _p of only 1,350°C, our predicted V _p anomalies are positive (+2 to +3%) for the cold slab and negative (?2 to ?4%) in the mantle wedge. These V _p estimates are consistent with the observed seismic tomography from P-wave arrivals, and therefore we conclude that a T _p of 1,350°C is a better estimate for the mantle potential temperature beneath Central Mexico. The new thermal model, in conjunction with phase diagrams for sediments, hydrated basalt and lithospheric mantle, have been used to estimate the amount and location of fluids released from the subducting Cocos slab. Several dehydration pulses have been identified along the slab interface where most of the fluids stored in sediments and oceanic crust are released into the overlying continental crust above the flat-slab. We found a good correlation between the pattern of these dehydration pulses and the location of NVTs, suggesting that slab dehydration is responsible for triggering the tremors. We suggest that NVT bursts localized above the flat slab segment represent the manifestation of ongoing continental crust hydration and weakening, a process that has been going on since 15?Ma ago when the Cocos slab entered into a flat-slab regime. Such continuous weakening would have reduced the suction forces that kept the slab in a flat regime in the last 15?Ma, allowing the slab to easily roll back. The continuous low-resistivity region recorded beneath the volcanic front in Central Mexico might represent the evidence of slab dehydration and crust weakening over time.     

Three-dimensional P-wave velocity structure derived from local earthquakes at the Katmai group of volcanoes,Alaska

The three-dimensional P-wave velocity structure beneath the Katmai group of volcanoes is determined by inversion of more than 10,000 rays from over 1000 earthquakes recorded on a local 18 station short-period network between September 1996 and May 2001. The inversion is well constrained from sea level to about 6 km below sea level and encompasses all of the Katmai volcanoes; Martin, Mageik, Trident, Griggs, Novarupta, Snowy, and Katmai caldera. The inversion reduced the average RMS travel-time error from 0.22 s for locations from the standard one-dimensional model to 0.13 s for the best three-dimensional model. The final model, from the 6th inversion step, reveals a prominent low velocity zone (3.6–5.0 km/s) centered at Katmai Pass and extending from Mageik to Trident volcanoes. The anomaly has values about 20–25% slower than velocities outboard of the region (5.0–6.5 km/s). Moderately low velocities (4.5–6.0 km/s) are observed along the volcanic axis between Martin and Katmai Caldera. Griggs volcano, located about 10 km behind (northwest of) the volcanic axis, has unremarkable velocities (5.0–5.7 km/s) compared to non-volcanic regions. The highest velocities are observed between Snowy and Griggs volcanoes (5.5–6.5 km/s). Relocated hypocenters for the best 3-D model are shifted significantly relative to the standard model with clusters of seismicity at Martin volcano shifting systematically deeper by about 1 km to depths of 0 to 4 km below sea level. Hypocenters for the Katmai Caldera are more tightly clustered, relocating beneath the 1912 scarp walls. The relocated hypocenters allow us to compare spatial frequency-size distributions (b-values) using one-dimensional and three-dimensional models. We find that the distribution of b is significantly changed for Martin volcano, which was characterized by variable values (0.8 < b < 2.0) with standard locations and more uniform values (0.8 < b < 1.2) after relocation. Other seismic clusters at Mageik (1.2 < b < 2.2), Trident (0.5 < b < 1.5) and Katmai Caldera (0.8 < b < 1.8) had stable b-values indicating the robustness of the observations. The strong high b-value region at Mageik volcano is mainly associated with an earthquake swarm in October, 1996 that possibly indicates a shallow intrusion or influx of gas. The new velocity and spatial b-value results, in conjunction with prior gravity (Bouguer anomalies up to − 40 mgal) and interferometry (several cm uplift) data, provide strong evidence in favor of partially molten rock at shallow depths beneath the Mageik–Katmai–Novarupta region. Moderately low velocities beneath Martin and Katmai suggest that old, mostly solidified intrusions exist beneath these volcanoes. Higher relative velocities beneath the Griggs and Snowy vents suggest that no magma is resident in the shallow crust beneath these volcanoes.     

Seismic Activity Around and Under Krakatau Volcano,Sunda Arc: Constraints to the Source Region of Island Arc Volcanics

There is general agreement that calc-alkaline volcanic rocks at convergent plate margins are genetically related to the process of subduction (Ringwood, 1974; Maaloe and Petersen, 1981; Hawkesworth et al., 1997). However, opinions on the mode and site of generation of primary magma for island arc volcanism differ substantially. The site of generation of calc-alkaline magma is thought to be either in the mantle wedge (Plank and Langmuir, 1988; McCulloch and Gamble, 1991) or in the subducting slab (White and Dupré, 1986; Defant and Drummond, 1990; Edwards et al., 1993; Ryan and Langmuir, 1993). We present seismological evidence in favour of the latter concept. A distinctive seismicity pattern around and under the Krakatau volcano was identified during systematic studies of the SE Asian convergent plate margins by means of global seismological data. A column-like cluster of events, probably associated with the dynamics of the volcano, is clearly separated from the events in the Wadati-Benioff zone. The accuracy of hypocentral determinations of the events of the cluster does not differ from the accuracy of the events belonging to the subducting slab. The depths of the cluster events vary from very shallow to about 100 km without any apparent discontinuity. On the other hand, there is a pronounced aseismic gap in the Wadati-Benioff zone directly beneath the volcano at depths between 100-150 km. The Krakatau cluster connects this aseismic gap to the volcano at the surface. The pervasive occurrence of earthquakes in the continental wedge between the subducting slab and the Earth surface bears witness to the brittle character of the continental lithosphere and casts doubt on the existence of large-scale melting of mantle material. The aseismic gap (Hanu and Vank, 1985), interpreted by us as a partially melted domain occurring in subducted slabs in practically all active subduction zones that reach depths greater than 100 km, is here used as evidence for the location of the primary source region of island arc volcanics in the subducting plate.     

Fine S wave velocity structure beneath Iwate volcano, northeastern Japan, as derived from receiver functions and travel times

A genetic algorithm inversion of receiver functions derived from a dense seismic network around Iwate volcano, northeastern Japan, provides the fine S wave velocity structure of the crust and uppermost mantle. Since receiver functions are insensitive to an absolute velocity, travel times of P and S waves propagating vertically from earthquakes in the subducting slab beneath the volcano are involved in the inversion. The distribution of velocity perturbations in relation to the hypocenters of the low-frequency (LF) earthquakes helps our understanding of deep magmatism beneath Iwate volcano. A high-velocity region (dV_S/V_S=10%) exists around the volcano at depths of 2–15 km, with the bottom depth decreasing to 11 km beneath the volcano’s summit. Just beneath the thinning high-velocity region, a low-velocity region (dV_S/V_S=−10%) exists at depths of 11–20 km. Intermediate-depth LF (ILF) events are distributed vertically in the high-velocity region down to the top of the low-velocity region. This distribution suggests that a magma reservoir situated in the low-velocity region supplies magma to a narrow conduit that is detectable by the hypocenters of LF earthquakes. Another broad low-velocity region (dV_S/V_S=−5 to −10%) occurs at depths of 17–35 km. Additional clusters of deep LF (DLF) events exist at depths of 32–37 km in the broad low-velocity zone. The DLF and ILF events are the manifestations of magma movement near the Moho discontinuity and in the conduit just beneath the volcano, respectively.     

The transport of water in subduction zones

The transport of water from subducting crust into the mantle is mainly dictated by the stability of hydrous minerals in subduction zones. The thermal structure of subduction zones is a key to dehydration of the subducting crust at different depths. Oceanic subduction zones show a large variation in the geotherm, but seismicity and arc volcanism are only prominent in cold subduction zones where geothermal gradients are low. In contrast, continental subduction zones have low geothermal gradients, resulting in metamorphism in cold subduction zones and the absence of arc volcanism during subduction. In very cold subduction zone where the geothermal gradient is very low(?5?C/km), lawsonite may carry water into great depths of ?300 km. In the hot subduction zone where the geothermal gradient is high(25?C/km), the subducting crust dehydrates significantly at shallow depths and may partially melt at depths of 80 km to form felsic melts, into which water is highly dissolved. In this case, only a minor amount of water can be transported into great depths. A number of intermediate modes are present between these two end-member dehydration modes, making subduction-zone dehydration various. Low-T/low-P hydrous minerals are not stable in warm subduction zones with increasing subduction depths and thus break down at forearc depths of ?60–80 km to release large amounts of water. In contrast, the low-T/low-P hydrous minerals are replaced by low-T/high-P hydrous minerals in cold subduction zones with increasing subduction depths, allowing the water to be transported to subarc depths of 80–160 km. In either case, dehydration reactions not only trigger seismicity in the subducting crust but also cause hydration of the mantle wedge. Nevertheless, there are still minor amounts of water to be transported by ultrahigh-pressure hydrous minerals and nominally anhydrous minerals into the deeper mantle. The mantle wedge overlying the subducting slab does not partially melt upon water influx for volcanic arc magmatism, but it is hydrated at first with the lowest temperature at the slab-mantle interface, several hundreds of degree lower than the wet solidus of hydrated peridotites. The hydrated peridotites may undergo partial melting upon heating at a later time. Therefore, the water flux from the subducting crust into the overlying mantle wedge does not trigger the volcanic arc magmatism immediately.     

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.     

Comparison between the Chile and Mexico triple junction areas substantiates slab window development beneath northwestern Mexico during the past 12-10 Myr

When combined with the Miocene-Recent volcanic record of Baja California, a parallel drawn between the Chile and Mexico triple junction areas substantiates slab window development beneath northwestern Mexico during the past 12-10 Myr. The slab-free zone manifestations challenge the notion that ridge subduction has not occurred beneath the southern Baja California peninsula. The geochemically distinctive rocks from the Santa Clara volcanic field of west-central Baja California, including coeval adakites and niobium-enriched basalt, are commonly inferred to signal partial melting of the subducting plate at shallow depths and relatively high temperatures, before slab dehydration occurs. Such PT conditions for slab melting have only been observed in association with spreading-ridge subduction. We propose that slab window development beneath southern Baja California and mainland Mexico (30° to 18°N) resulted from subduction of the East Pacific rise.     

Eruptive history and tectonic setting of Medicine Lake Volcano,a large rear-arc volcano in the southern Cascades

Medicine Lake Volcano (MLV), located in the southern Cascades ∼ 55 km east-northeast of contemporaneous Mount Shasta, has been found by exploratory geothermal drilling to have a surprisingly silicic core mantled by mafic lavas. This unexpected result is very different from the long-held view derived from previous mapping of exposed geology that MLV is a dominantly basaltic shield volcano. Detailed mapping shows that < 6% of the ∼ 2000 km² of mapped MLV lavas on this southern Cascade Range shield-shaped edifice are rhyolitic and dacitic, but drill holes on the edifice penetrated more than 30% silicic lava. Argon dating yields ages in the range ∼ 475 to 300 ka for early rhyolites. Dates on the stratigraphically lowest mafic lavas at MLV fall into this time frame as well, indicating that volcanism at MLV began about half a million years ago. Mafic compositions apparently did not dominate until ∼ 300 ka. Rhyolite eruptions were scarce post-300 ka until late Holocene time. However, a dacite episode at ∼ 200 to ∼ 180 ka included the volcano's only ash-flow tuff, which was erupted from within the summit caldera. At ∼ 100 ka, compositionally distinctive high-Na andesite and minor dacite built most of the present caldera rim. Eruption of these lavas was followed soon after by several large basalt flows, such that the combined area covered by eruptions between 100 ka and postglacial time amounts to nearly two-thirds of the volcano's area. Postglacial eruptive activity was strongly episodic and also covered a disproportionate amount of area. The volcano has erupted 9 times in the past 5200 years, one of the highest rates of late Holocene eruptive activity in the Cascades. Estimated volume of MLV is ∼ 600 km³, giving an overall effusion rate of ∼ 1.2 km³ per thousand years, although the rate for the past 100 kyr may be only half that. During much of the volcano's history, both dry HAOT (high-alumina olivine tholeiite) and hydrous calcalkaline basalts erupted together in close temporal and spatial proximity. Petrologic studies indicate that the HAOT magmas were derived by dry melting of spinel peridotite mantle near the crust mantle boundary. Subduction-derived H₂O-rich fluids played an important role in the generation of calcalkaline magmas. Petrology, geochemistry and proximity indicate that MLV is part of the Cascades magmatic arc and not a Basin and Range volcano, although Basin and Range extension impinges on the volcano and strongly influences its eruptive style. MLV may be analogous to Mount Adams in southern Washington, but not, as sometimes proposed, to the older distributed back-arc Simcoe Mountains volcanic field.     

Phase relations in hydrous MORB at 18-28 GPa: implications for heterogeneity of the lower mantle

The phase relations in hydrous and anhydrous mid-ocean ridge basalt were determined at pressures of 18-28 GPa. Liquidus phase relations in hydrous and anhydrous MORB are different. Garnet is the liquidus phase at pressures below 21 GPa, Ca-Al (CAS) phase and stishovite are the liquidus phases at pressures of 22-27 GPa, and stishovite and Ca-perovskite are the liquidus phases above 27 GPa, whereas Ca-perovskite is a liquidus phase of anhydrous MORB at pressures above 23 GPa. Under subsolidus conditions, we have found that in the hydrous MORB system the stability fields of Al-bearing perovskite and Na-Al (NAL) phase might shift to lower pressure by about 1.5 GPa compared to the dry MORB system. This shift could be explained by oxidation of a garnet-bearing assemblage by hydrous fluid and formation of Fe³⁺-bearing aluminous perovskite at lower pressures relative to the anhydrous system and/or differences in water solubility of the phases existing in perovskite-bearing assemblages. Our data indicate that hydrous basaltic crust remains denser than peridotite along the geotherm of a subducting slab, i.e. there is no density crossover between peridotite and basalt. Therefore, in slabs going through the 660 km discontinuity, basalt would gravitationally sink into the lower mantle under relatively hydrous conditions. The delamination of former basaltic crust near the 660 km discontinuity might be possible under relatively dry conditions of subduction. There are no stable highly hydrous phases in MORB above 10 GPa even at lower temperatures corresponding to subducting slabs. Therefore, MORB cannot be an important carrier of water to the deep Earth interior. However, it can be constantly supplied by water-bearing fluid from the underlying peridotite part of the descending slab. Thus, it is plausible that water can control subduction of the oceanic crust into the lower mantle.     

Receiver functions in northeast China - implications for slab penetration into the lower mantle in northwest Pacific subduction zone

Seismic studies of the subducting lithosphere and the upper mantle discontinuities in the northwest Pacific subduction zone beneath Japan and northeast China have suggested contrary subduction scenarios. There was little consensus on the issue whether the subducting slab penetrates the upper mantle discontinuities into the lower mantle or it is deflected atop of the 660-km discontinuity over several hundred kilometers. We calculate receiver functions from a recent seismic broadband station network located in northeast China and find topographic variations of the upper mantle discontinuities. A deeper-than-normal 660-km discontinuity is observed over an area of 400 km and it coincides with the stagnant slab imaged by seismic tomography. The 660-km discontinuity is locally depressed by more than 35 km and the transition zone is thickened by more than 20 km in the east of the region where it encounters the slab. These observations provide evidence of the slab accumulating in the mantle transition zone and locally penetrating into the lower mantle.     

Ukinrek Maars, Alaska, I. April 1977 eruption sequence, petrology and tectonic setting

During ten days of phreatomagmatic activity in early April 1977, two maars formed 13 km behind the Aleutian arc near Peulik volcano on the Alaska Peninsula. They have been named “Ukinrek Maars”, meaning “two holes in the ground” in Yupik Eskimo. The western maar formed at the northwestern end of a low ridge within the first three days and is up to 170 m in diameter and 35 m in depth. The eastern maar formed during the next seven days 600 m east of West Maar at a lower elevation in a shallow saddle on the same ridge and is more circular, up to 300 m in diameter and 70 m in depth. The maars formed in terrain that was heavily glaciated in Pleistocene times. The groundwater contained in the underlying till and silicic volcanics from nearby Peulik volcano controlled the dominantly phreatomagmatic course of the eruption.During the eruptions, steam and ash clouds reached maximum heights of about 6 km and a thin blanket of fine ash was deposited north and east of the vents up to a distance of at least 160 km. Magma started to pool on the floor of East Maar after four days of intense phreatomagmatic activity.The new melt is a weakly undersaturated alkali olivine basalt (Ne = 1.2%) showing some transitional character toward high-alumina basalts. The chemistry, an anomaly in the tholeitic basalt-andesite-dominated Aleutian arc, suggests that the new melt is primitive, generated at a depth of 80 km or greater by a low degree of partial melting of garnet peridotite mantle with little subsequent fractionization during transport.The Pacific plate subduction zone lies at a depth of 150 km beneath the maars. Their position appears to be tectonically controlled by a major regional fault, the Bruin Bay fault, and its intersection with cross-arc structural features. We favor a model for the emplacement of the Ukinrek Maars that does not link the Ukinrek conduit to the plumbing system of nearby Peulik volcano. The Ukinrek eruptions probably represent a genetically distinct magma pulse originating at asthenospheric depths beneath the continental lithosphere.     

Geological,tomographic, kinematic and geodynamic constraints on the dynamics of sinking slabs

We use geodynamic models with imposed plate velocities to test the forward-modeled history of subduction based on a particular plate motion model against alternative seismic tomography models. We utilize three alternative published reference frames: a hybrid moving hotspot-palaeomagnetic, a hybrid moving hotspot-true polar wander corrected-palaeomagnetic, and a Subduction Reference Frame, a plate model including longitudinal shifts of subduction zones by matching subduction volumes imaged by P-wave tomography, to assess which model best predicts present day mantle structure compared with seismic tomography and volumetrically derived subduction history. Geodynamic modeling suggests paleo-longitudinal corrections applied to the Subduction Reference Frame result in lower mantle slab material beneath North America and East Asia accumulating up to 10–15° westward of that imaged by tomography, whereas the hybrid models develop material offset by 2–9°. However, the Subduction Reference Frame geodynamic model produces slab material beneath the Tethyan Domain coinciding with slab volumes imaged by tomography, whereas the hybrid reference frame models do not, suggesting regional paleo-longitudinal corrections are required to constrain slab locations. We use our models to test inferred slab sinking rates in the mantle focusing on well-constrained regions. We derive a globally averaged slab-sinking rate of 13 ± 3 mm/yr by combining the ages of onset and cessation of subduction from geological data and kinematic reconstructions with images of subducted slabs in the mantle. Our global average slab-sinking rate overlaps with the 15–20 mm/yr rate implied by mantle convection models using a lower mantle viscosity 100 times higher than the upper mantle.     

中国东北地区远震P波走时层析成像研究

利用中国东北流动和固定台网的234个宽频带地震仪记录的远震波形数据,采用波形相关方法拾取了57251个有效相对走时残差数据,进一步采用FMTT(Fast Marching Teleseismic Tomography)层析成像的方法,反演获取了研究区下方深达800 km的P波速度结构.结果显示:在长白山下方发现有一个高速异常结构,这可能就是俯冲到欧亚大陆板块下方的太平洋板块,由于板块的部分下沉,使得板块的形状并没有呈现出明显的板片状.长白山、阿尔山、五大连池火山下方都有低速异常体,长白山和阿尔山下的低速异常向下延伸至地幔转换带,可能与其上部的火山形成有关.五大连池火山下方的低速异常向下延伸至200 km左右,不同埋深的低速异常结构可能意味着五大连池与长白山和阿尔山有着不同的成因.松辽盆地呈现以高速异常为主导高低速异常混合分布的特性,暗示松辽盆地可能有岩石圈拆沉的过程,盆地南部下方的低速异常与长白山和阿尔山下的低速异常有连通性,可能是下地幔热物质上涌的一个通道.