Volcanic Markers of the Post-Subduction Evolution of Baja California and Sonora, Mexico: Slab Tearing Versus Lithospheric Rupture of the Gulf of California

The study of the geochemical compositions and K-Ar or Ar-Ar ages of ca. 350 Neogene and Quaternary lavas from Baja California, the Gulf of California and Sonora allows us to discuss the nature of their mantle or crustal sources, the conditions of their melting and the tectonic regime prevailing during their genesis and emplacement. Nine petrographic/geochemical groups are distinguished: ??regular?? calc-alkaline lavas; adakites; magnesian andesites and related basalts and basaltic andesites; niobium-enriched basalts; alkali basalts and trachybasalts; oceanic (MORB-type) basalts; tholeiitic/transitional basalts and basaltic andesites; peralkaline rhyolites (comendites); and icelandites. We show that the spatial and temporal distribution of these lava types provides constraints on their sources and the geodynamic setting controlling their partial melting. Three successive stages are distinguished. Between 23 and 13 Ma, calc-alkaline lavas linked to the subduction of the Pacific-Farallon plate formed the Comondú and central coast of the Sonora volcanic arc. In the extensional domain of western Sonora, lithospheric mantle-derived tholeiitic to transitional basalts and basaltic andesites were emplaced within the southern extension of the Basin and Range province. The end of the Farallon subduction was marked by the emplacement of much more complex Middle to Late Miocene volcanic associations, between 13 and 7 Ma. Calc-alkaline activity became sporadic and was replaced by unusual post-subduction magma types including adakites, niobium-enriched basalts, magnesian andesites, comendites and icelandites. The spatial and temporal distribution of these lavas is consistent with the development of a slab tear, evolving into a 200-km-wide slab window sub-parallel to the trench, and extending from the Pacific coast of Baja California to coastal Sonora. Tholeiitic, transitional and alkali basalts of subslab origin ascended through this window, and adakites derived from the partial melting of its upper lip, relatively close to the trench. Calc-alkaline lavas, magnesian andesites and niobium-enriched basalts formed from hydrous melting of the supraslab mantle triggered by the uprise of hot Pacific asthenosphere through the window. During the Plio-Quaternary, the ??no-slab?? regime following the sinking of the old part of the Farallon plate within the deep mantle allowed the emplacement of alkali and tholeiitic/transitional basalts of deep asthenospheric origin in Baja California and Sonora. The lithospheric rupture connected with the opening of the Gulf of California generated a high thermal regime associated to asthenospheric uprise and emplaced Quaternary depleted MORB-type tholeiites. This thermal regime also induced partial melting of the thinned lithospheric mantle of the Gulf area, generating calc-alkaline lavas as well as adakites derived from slivers of oceanic crust incorporated within this mantle.     

Heat flow,radioactive heat generation,and theoretical tectonics for northwestern Mexico

Nine heat flow values from Baja California and Sonora represent the first determinations of heat flux from Mexico and, in general, justify southerly extensions of the Basin and Range and Southern California thermal provinces of the southwestern United States. Heat flow in northern Baja California is low (0.84 hfu), consistent with known values from the Sierra Nevada Batholith, but a value of 1.23 hfu was recorded in the central area of the peninsula. Seven measurements from Sonora, averaging 2.5 hfu, are similar to Basin and Range values.Reduced heat flow values (q*) have been calculated for Sonora and are interpreted as indicative of a separate Basin and Range thermal subprovince in Mexico. A tectonic model developed for northwestern Mexico attributes the high heat flow of the Basin and Range to behind arc extension associated with the termination of Cenozoic subduction. Opening of the Gulf of California and sea-floor spreading therein has provided a vent for the release of heat accumulated under northwestern Mexico. Thus, the southern portion of the Basin and Range province (in Mexico) appears to be cooling and narrowing and a subprovince of high reduced heat flow is defined east and perhaps southwest of the Gulf of California.     

Thermal consequences of the formation of a slab window beneath the Mid-Cretaceous southwest Japan arc: A 2-D numerical analysis

Abstract   The development of voluminous granitic magmatism and widespread high-grade metamorphism in Mid-Cretaceous southwest Japan have been explained by the subduction of a spreading ridge (Kula–Pacific or Farallon–Izanagi plate boundaries) beneath the Eurasian continent and the formation of a slab window. In the present study, the thermal consequences of the formation of a slab window beneath a continental margin are evaluated through a 2-D numerical simulation. The model results are evaluated by comparison with the Mid-Cretaceous geology of southwest Japan. Of particular interest are the absence of an amphibolite- to granulite-facies metamorphic belt near the Wadati–Benioff plane, and significant melting of the lower crustal-mafic rocks sufficient to form a large amount of granitic magma. Because none of the model results simultaneously satisfied these two geological interpretations, it is suggested that subduction of plate boundaries in Mid-Cretaceous southwest Japan was not associated with the opening of a slab window. According to previous studies, and the results of the present study, two different tectonic scenarios could reasonably explain the geological interpretations for Mid-Cretaceous southwest Japan: (i) The spreading ridge did not subduct beneath the Eurasian continent, but was located off the continental margin, implying the continuous subduction of very young oceanic lithosphere; (ii) ridge subduction beneath the continental margin occurred after active spreading had ceased. Consequently, in both tectonic scenarios, the subduction of plate boundaries at the Mid-Cretaceous southwest Japan was not associated with a slab window, but very young (hot) oceanic lithosphere.     

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.     

Geochemistry of Cenezoic volcanic rocks, Baja California, Mexico: Implications for the petrogenesis of post-subduction magmas

Late Cenozoic volcanism in Baja California records the effects of cessation of subduction at a previously convergent, plate margin. Prior to 12.5 m.y., when subduction along the margin of Baja ceased, the predominant volcanic activity had a calc-alkaline signature, ranging in composition from basalt to rhyolite. Acidic pyroclastic activity was common, and possibly represented the westermost, distal edge of the Sierra Madre Occidental province. After 12.5 m.y., however, the style and composition of the magmatic products changed dramatically. The dominant rock type within the Jaraguay and San Borja volcanic fields is a magnesian andesite, with up to 8% MgO at 57% SiO₂, low Fe/Mg ratios, and high Na/K ratios. These rocks have unusual trace-element characteristics, with high abundances of Sr (up to 3000 ppm), low contents of Rb; K/Rb ratios are very high (usually over 1000, and up to 2500), and Rb/Sr ratios are low (less than 0.01). Furthermore, La_n/Yb_n ratios are high, consistent with derivation from a mantle source with fractionated REE patterns. ⁸⁷Sr/⁸⁶Sr ratios are less than 0.7048, and usually less than 0.7040, whereas the pre-12.5 m.y. lavas have ⁸⁷Sr/⁸⁶Sr ratios between 0.7038 and 0.7063. We have previously termed these rocks bajaites, in order to distinguish them from other magnesian andesites. Bajaites also occur in southernmost Chile and the Aleutian Islands, areas which also have histories of attempted or successful ridge subduction.It is proposed that the bajaite series is produced during the unusual physico-chemical conditions operating during the subduction of young oceanic lithosphere, or subduction of a spreading centre. During normal subduction, the oceanic crust dehydrates, releasing volatiles (water, Rb and other large-ion lithophile elements) into the overlying wedge. Subduction of younger crust will result in a progressive decrease, and eventual cessation of the transfer of volatiles when subduction stops. Thermal rebound of the mantle may cause the slab to melt, perhaps under eclogitestable conditions. The resulting melt will be heavy-REE-depleted, perhaps dacitic, but will otherwise inherit MORB-like Rb/Sr and K/Rb ratios. The ascending melt will react with the mantle to form the source of the bajaitic rocks. Furthermore, any amphibole in the mantle, stabilised during the higher PH₂O conditions of earlier subduction, will break down and contribute a high-K/Rb ratio component.The implications of this study are that firstly, the subducted slab does not contribute a highly fractionated REE component in most modern arcs (i.e. the slab does not melt); secondly, Rb has a very short residence time in the mantle, and its abundance in arc rocks is a direct reflection of the input from the dehydrating slab; and thirdly, bajaitelike rocks may provide recognition of attempted or successful ridge subduction in the geologic past.     

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.     

Adakitic lavas in the Central Luzon back‐arc region,Philippines: lower crust partial melting products?

Abstract Volcanism in the back-arc side region of Central Luzon, Philippines, with respect to the Manila Trench is characterized by fewer and smaller volume volcanic centers compared to the adjacent forearc side-main volcanic arc igneous rocks. The back-arc side volcanic rocks which include basalts, basaltic andesites, andesites and dacites also contain more hydrous minerals (ie, hornblende and biotite). Adakite-like geochemical characteristics of these back-arc lavas, including elevated Sr, depleted heavy rare earth elements and high Sr/Y ratios, are unlikely to have formed by slab melting, be related to incipient subduction, slab window magmatism or plagioclase accumulation. Field and geochemical evidence show that these adakitic lavas were most probably formed by the partial melting of a garnet-bearing amphibolitic lower crust. Adakitic lavas are not necessarily arc–trench gap region slab melts.     

Slab partial melts from the metasomatizing agent to adakite,Tafresh Eocene volcanic rocks,Iran

Geochemical and mineralogical characteristics of the Eocene volcanic succession in Tafresh area of the Urumieh–Dokhtar Magmatic Assemblage (UDMA) are unique in the 2000‐km‐length assemblage. Demonstrating rather steep rare earth element (REE) patterns and the widespread presence of amphibole (+biotite) phenocrysts are two distinct characters that dominate the Eocene volcanic succession of mainly andesitic composition. Coincidence of the geochemical and mineralogical characteristics of the whole volcanic succession with adakites, rather amphibole‐ (+biotite) rich dacitic (with 61–64 wt% SiO₂) stocks and dykes, is considered as the key in unraveling the role of ‘slab‐derived melt contribution’ in petrogenesis of the volcanic succession. Slab‐derived melting has been an ongoing process that metasomatized some parts of the mantle wedge from which hybrid rocks (andesites) are derived. Basalts with distinct signatures of slab melt metasomatism are yet another support for the occurrence of slab melting. Interlayering of normal, island‐arc‐type calc‐alkaline volcanic rocks with the slab‐melt metasomatized basalts and hybrid andesites suggests that the slab melting has been motivated by the subduction. Formation of the Tafresh Caldera, the likely consequence of an explosive eruption, is compatible with the volatile‐bearing nature of the adakitic volcanism in the study area. It is indicated by the ubiquitous presence of the hydrous minerals. Beneath the Tafresh area, in Eocene time, the subducting slab seems to have reached a critical high depth that is enough for the development of amphibolite–eclogite. The slab deformation, motivated by the geometry of subduction and/or the underlying mantle's steeper geotherms, is suggested to have resulted in the slab melting that helped develop a rock assemblage unique to the UDMA.     

Double-difference tomography of P- and S-wave velocity structure beneath the western part of Java,Indonesia

West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults. In this study, double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java. To infer the geometry of the structure beneath the study area, precise earthquake hypocenter determination was first performed before tomographic imaging. For this, earthquake waveform data were extracted from the regional Meteorological, Climatological, Geophysical Agency (BMKG) network of Indonesia from South Sumatra to Central Java. The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected, the key features being events of magnitude > 3, azimuthal gap < 210° and number of phases > 8. A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters. The hypocenter locations were then relocated using double-difference tomography (tomoDD). A significant reduction of travel-time (root mean square basis) and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java. Double-difference tomography was found to give a clear velocity structure, especially beneath the volcanic arc area, i.e., under Mt Anak Krakatau, Mt Salak and the mountains complex in the southern part of West Java. Low velocity anomalies for the P and S waves as well as the v_P/v_S ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.     

The petrology and geochemistry of volcanic rocks from the Northern Harşit river area, Pontid Volcanic Province, Northeast Turkey

Two volcanic cycles can be distinguished, in the Pontid magmatic arc. They comprise an Upper Cretaceous, Lower Volcanic Cycle of which only the waning stage contains abundant pyroclastic rocks. The latter show spatial association to the fault pattern and are closely related to mineralisation. The overlying, early Tertiary, Upper Volcanic Cycle shows evidence of explosive vulcanicity in the Upper Basic Series. Dacites and rhyodacites are only locally developed and again show spatial association with the faulting.Comparison of the major and minor element chemistries of the two cycles demonstrates the clear separation into a lower tholeiitic and an upper calc-alkaline cycle. The rocks show similar chemistry to volcanic suites from island arcs in other areas.The origin of the tholeiitic magma is ascribed to melting of “dry” amphibolite during early subduction of Tethyan ocean floor beneath “Pontian Land”. This resulted in low K abundances and K/Rb ratios, and some Fe enrichment in the tholeiitic basalts.The calc-alkaline magma is thought to be derived from a later stage in the subduction process when melting of amphibole was joined by melting of biotite or phlogopite. The Upper Volcanic Cycle is thus enriched in K and shows no Fe enrichment due to a probable higher water content. The higher Cr and Ni contents of the Upper Volcanic Cycle, together with K, may imply some melting of lherzolite overlying the subducted slab.     

马尼拉海沟俯冲带热结构的模拟研究

俯冲带热结构的数值模拟研究是对地表观测研究的重要补充,也是验证地球动力学模型的重要方法.本文沿马尼拉海沟俯冲带东火山链(EVC)和西火山链(WVC)各取一条剖面,依据地质、地球物理条件,进行了有限元热模拟计算.计算过程中,分析了摩擦和剪切热对俯冲带热结构的影响,模拟了EVC和WVC两条测线下俯冲带的热结构,并结合岩石学实验结果预测了俯冲板块发生脱水和部分熔融的位置.模拟结果表明,在100 km深度处,考虑摩擦和剪切热时,俯冲板块表面的温度约为865 ℃;而不考虑摩擦和剪切时,俯冲板块表面的温度仅为770 ℃,二者温差可达95 ℃.在相同深度处,考虑摩擦和剪切热时,在EVC和WVC测线下俯冲板块表面的温度分别为865 ℃和895 ℃,俯冲洋壳底部温度分别为560 ℃和605 ℃.俯冲板块表面少量矿物开始脱水的深度小于50 km,但大量脱水和部分熔融主要发生在深度100 km左右,这与地表观测的火山活动位置一致.     

Magmatic response to early aseismic ridge subduction: the Ecuadorian margin case (South America)

A geochemical and isotopic study of lavas from Pichincha, Antisana and Sumaco volcanoes in the Northern Volcanic Zone (NVZ) in Ecuador shows their magma genesis to be strongly influenced by slab melts. Pichincha lavas (in fore arc position) display all the characteristics of adakites (or slab melts) and were found in association with magnesian andesites. In the main arc, adakite-like lavas from Antisana volcano could be produced by the destabilization of pargasite in a garnet-rich mantle. In the back arc, high-niobium basalts found at Sumaco volcano could be produced in a phlogopite-rich mantle. The strikingly homogeneous isotopic signatures of all the lavas suggest that continental crust assimilation is limited and confirm that magmas from the three volcanic centers are closely related. The following magma genesis model is proposed in the NVZ in Ecuador: in fore arc position beneath Pichincha volcano, oceanic crust is able to melt and produces adakites. En route to the surface, part of these magmas metasomatize the mantle wedge inducing the crystallization of pargasite, phlogopite and garnet. In counterpart, they are enriched in magnesium and are placed at the surface as magnesian andesites. Dragged down by convection, the modified mantle undergoes a first partial melting event by the destabilization of pargasite and produces the adakite-like lavas from Antisana volcano. Lastly, dragged down deeper beneath the Sumaco volcano, the mantle melts a second time by the destabilization of phlogopite and produces high-niobium basalts. The obvious variation in spatial distribution (and geochemical characteristics) of the volcanism in the NVZ between Colombia and Ecuador clearly indicates that the subduction of the Carnegie Ridge beneath the Ecuadorian margin strongly influences the subduction-related volcanism. It is proposed that the flattening of the subducted slab induced by the recent subduction (<5 Ma?) of the Carnegie Ridge has permitted the progressive warming of the oceanic crust and its partial melting since ca. 1.5 Ma. Since then, the production of adakites in fore arc position has deeply transformed the magma genesis in the overall arc changing from ‘typical’ calc-alkaline magmatism induced by hydrous fluid metasomatism, to the space- and time-associated lithology adakite/high-Mg andesite/adakite-like andesite/high-Nb basalts characteristic of slab melt metasomatism.     

Imaging the Moho and Subducted Oceanic Crust at the Isthmus of Tehuantepec, Mexico, from Receiver Functions

Using teleseismic data recorded along a transect, which we call VEOX (for Veracruz-Oaxaca seismic line), of 46 broadband stations installed across the Isthmus of Tehuantepec in southern Mexico, we obtained receiver functions and stacked them to study the Moho topography and back projected them to visualize the subducted slab geometry beneath the isthmus. We observed a back-azimuth dependent Moho thickness across the transect, particularly beneath the Los Tuxtlas Volcanic Field. Also, we observed the Cocos plate which subducts with an angle of 26° between 140 and 310?km from the trench. Comparison with regional seismicity indicates that it occurs below the oceanic crust.     

An evaluation of the global variations in the major element chemistry of arc basalts

Arc volcanoes occur at convergent margins with a wide range in subduction parameters, and variations in these parameters might be expected to lead to variations in the chemistry of magmas parental to arcs. Major element analyses from approximately 100 volcanic centers within 30 arcs, normalized to 6% MgO to minimize the effects of crystal fractionation, display wide variations. Na₂O and CaO at 6% MgO (Na_6.0 and Ca_6.0) correlate remarkably well with the thickness of the overlying crust. These systematics are consistent with two possible models. In the first model, the crust behaves as a chemical filter; where the crust is thick, magmas crystallize at higher pressure and interact more extensively with the arc crust. Modeling of high-pressure crystallization and assimilation, however, does not reproduce the associated variations in Na_6.0 and Ca_6.0 without calling upon complicated combinations of fractionating phases and assimilants. In the second model, crustal thickness determines the height of the mantle column available for melting beneath arc volcanoes. If melting begins beneath arcs at similar depths, then the column of mantle that undergoes decompression melting is much shorter beneath the thickest arc crust. The shorter mantle column for arcs built on thick crust will lead to smaller extents of melting in the mantle, and hence higher Na_6.0 and lower Ca_6.0 in the parental magmas. Modeling shows that variations in the extent of melting in the mantle can easily account for the associated variations in Ca_6.0 and Na_6.0. The abundances of the other major elements at 6% MgO do not correlate well with crustal thickness, or any other subduction parameter. Co-variation of some of these other major elements (e.g., Si_6.0 and Fe_6.0) within individual arcs suggests that they are strongly influenced by local crustal level processes that obscure partial melting systematics. Correction for the crustal processes improves the relationship between Na_6.0 and Ca_6.0 that is so readily explained by partial melting. The extents of melting in the mantle beneath arc volcanoes estimated from the ranges in Na_6.0 and Ca_6.0 are remarkably similar to those estimated beneath mid-ocean ridges. This observation provides further evidence that the mantle wedge, and not the slab, melts beneath arc volcanic fronts.     

An Evaluation of Proposed Mechanisms of Slab Flattening in Central Mexico

Central Mexico is the site of an enigmatic zone of flat subduction. The general geometry of the subducting slab has been known for some time and is characterized by a horizontal zone bounded on either side by two moderately dipping sections. We systematically evaluate proposed hypotheses for shallow subduction in Mexico based on the spatial and temporal evidence, and we find no simple or obvious explanation for the shallow subduction in Mexico. We are unable to locate an oceanic lithosphere impactor, or the conjugate of an impactor, that is most often called upon to explain shallow subduction zones as in South America, Japan, and Laramide deformation in the US. The only bathymetric feature that is of the right age and in the correct position on the conjugate plate is a set of unnamed seamounts that are too small to have a significant effect on the buoyancy of the slab. The only candidate that we cannot dismiss is a change in the dynamics of subduction through a change in wedge viscosity, possibly caused by water brought in by the slab.     

Genesis of lavas of the Taupo Volcanic Zone, North Island, New Zealand

The Taupo Volcanic Zone forms part of the Taupo-Hikurangi subduction system, and comprises five volcanic centres: Tongariro, Taupo, Maroa, Okataina and Rotorua. Tongariro Volcanic Centre is formed almost entirely of andesite while the other four centres contain predominantly rhyolitic volcanics and later fissure eruptions of high-Al basalt. Estimated total volume of each lava type are as follows: 2 km³ of high-Al basalt (< 0.1%); 260 km³ of andesite (< 2.5%); 5 km³ of dacite (< 0.1%); > 10,000 km³ of rhyolite and ignimbrite (> 97.4%).The location of the andesites and vent alignments suggest a source from a subduction zone underlying the area. However, the lavas differ chemically from island-arc andesites such as those of Tonga; in particular by having higher contents of the alkali elements, light REE and Sr and Pb isotopes. This suggests some crustal contamination, and it is considered that this may occur beneath the wide accretionary prism of the subduction system. Amphibolite of the subduction zone will break down between 80 and 100 km and a partial melt will rise. A multi-stage process of magma genesis is then likely to occur. High-Al basalts are thought to be derived from partial melting of a garnet-free peridotite near the top of the mantle wedge overlying the subduction zone, locations of the vents controlled largely by faults within the crust. Rhyolites and ignimbrites were probably derived from partial melting of Mesozoic greywacke and argillite under the Taupo Volcanic Zone. Initial partial melting may have been due to hydration of the base of the crust; the “water” having come from dehydration of the downgoing slab. The partial melts would rise to form granodiorite plutons and final release of the magma to form rhyolites and ignimbrites was allowed because of extension within the Taupo graben.Dacites of the Bay of Plenty probably resulted from mixing of andesitic magma with small amounts of rhyolitic magma, but those on the eastern side of the Rotorua-Taupo area were more likely formed by a higher degree of partial melting of the Mesozoic greywacke-argillite basement. This may be due to intrusion of andesite magma on this side of the Taupo volcanic zone.     

Geology and K-Ar dating of the Tuxtla Volcanic Field,Veracruz, Mexico

The Tuxtla Volcanic Field (TVF) is located on the coast of the Gulf of Mexico in the southern part of the state of Veracruz, Mexico. Volcanism began about 7 my ago, in the Late Miocene, and continued to recent times with historical eruptions in ad 1664 and 1793. The oldest rocks occur as highly eroded remnants of lava flows in the area surrounding the historically active cone of San Martín Tuxtla. Between about 3 and 1 my ago, four large composite volcanoes were built in the eastern part of the area. Rocks from these structures are hydrothermally altered and covered with lateritic soils, and their northern slopes show extensive erosional dissection that has widened preexisting craters to form erosional calderas. The eastern volcanoes are composed of alkali basalts, hawaiites, mugearites, and benmoreites, with less common calc-alkaline basaltic andesites and andesites. In the western part of the area, San Martín Tuxtla Volcano and its over 250 satellite cinder cones and maars produced about 120 km³ of lava over the last 0.8 my. A ridge of flank cinder cones blocked drainage to the north to form Laguna Catemaco. Lavas erupted from San Martín and its flank vents are restricted to compositions between basanite and alkali basalt. The alignment of major volcanoes and flank vents along a N55°W trend suggests an extensional stress field in the crust with a minimum compressional stress orientation of N35° E. In total, about 800 km³ of lava has been erupted in the TVF in the last 7 my. This gives a magma output rate of about 0.1 km³/1000 year, a value smaller than most composite cones, but similar to cinder cone fields that occur in central Mexico. Individual eruptions over the last 5000 years had volumes on the order of 0.1km³, with average recurrence intervals of 600 years. The alkaline compositions of the TVF lavas contrast markedly with the calc-alkaline compositions erupted in the subduction-related Mexican Volcanic Belt to the west, leading previous workers to suggest that the TVF is not related to subduction. Trace-element signatures of TVF lavas indicate, however, that they are probably related to subduction. We suggest that the alkaline character of the TVF lavas is the result of low degrees of melting of a mantle source coupled with a stress regime that allows these small-volume melts to reach the surface in the TVF.     

圣安德烈斯断层门多西诺三联点迁移的构造热效应

由于圣安德烈斯断层门多西诺三联点向北迁移，或北美板块的相对向南运动，驱动角落流使地幔热物质向板片窗上涌，造成地幔温度升高和部分熔融．本文通过二维暂态有限单元法计算表明，当北美盖层厚度为20km时，计算预测的地热分布、部分熔融百分比、岩浆岩活动的时空分布特征等都与实际地热、地震和地质观测很好地吻合，确认了角落流机制和它引起的暂态热过程的重要性．     

Curie Point Depth Estimates and Correlation with Subduction in Mexico

We investigate the regional thermal structure of the crust in Mexico using Curie Point Depth (CPD) estimates. The top and bottom of the magnetized crust were calculated using the power-density spectra of the total magnetic field from the freely available ??Magnetic Anomaly Map of North America??. We applied this method to estimate the regional crustal thermal structure in overlapping square windows of 2°?×?2°. The CPD estimates range between 10 and 40?km and show several regions of relatively shallow and deep magnetic sources, with a general inverse correlation with measured heat flow. A deep CPD region (20?C30?km) is located in the fore-arc area where the subducting Cocos plate has a flat-slab geometry. This deep region is bound to the NW and SE by shallow CPD areas beneath the states of Michoacan (CPD?=?12?C16?km) and Oaxaca (CPD?=?~16?km), respectively. There is a good spatial correlation between this deep CPD area and two main fracture zones located on the incoming Cocos plate (Orozco and O??Gorman fracture zones), suggesting that subduction plays an important role in setting apart different CPD provinces along the Mexican coast. Another deep CPD (16?C32?km) area corresponds to the region where the Rivera plate subducts beneath Jalisco block. The Trans-Mexican Volcanic Belt is characterized by a decrease in Curie depths from west (16?C20?km) to east (10?C12?km). Finally, several deep CPD areas are situated in the back-arc region where old Mesozoic terrains are present. Our results suggest that the main control on the crust??s regional thermal structure in the fore-arc and volcanic arc regions is due to the subduction of the Cocos and Rivera plates beneath Mexico.     

智利三联点南部扩张洋脊俯冲区域岩石层热结构的数值模拟

智利三联点(CTJ)位于纳兹卡板块、南极洲板块与南美板块的交界处,由南极洲—纳兹卡板块之间的智利洋脊俯冲到智利海沟而形成.巴塔哥尼亚板片窗的发展是智利洋脊长期扩张俯冲的结果之一.随着纳兹卡板块的不断东向俯冲,纳兹卡板块范围逐渐变小,CTJ同时向北移动.本文采用数值模拟方法,建立了关于洋脊海沟碰撞的简单二维模型,来研究智利三联点南部扩张洋脊俯冲区域岩石层的热结构.模拟结果表明,洋脊的位置、板块相对汇聚速度及上覆大陆板块的存在均对俯冲区域海洋板块的温度结构有着很大影响,并且大陆板块下方海洋板块温度变化最大的位置距洋脊的水平距离与洋脊到板片窗范围的水平距离两者之间具有较好的一致性.同时,当存在两两板块间的相对汇聚时,洋脊右侧大陆板块下表面的温度升高,俯冲带内海洋板块温度接近于地幔温度.纳兹卡板块以7.8 cm·a~(-1)的速度急速俯冲于南美板块之下的过程中,同时伴随着智利洋脊的持续扩张俯冲,在智利三联点南部,南美板块之下纳兹卡板块的温度因而可以更快地达到地幔软流层的约1300℃温度,并最终消亡于地幔之中.