Role of subduction dynamics on the unevenly distributed volcanism at the Middle American subduction system

A typical subduction of an oceanic plate beneath a continent is expected to be accompanied by arc volcanoes along the convergent margin. However, subduction of the Cocos plate at the Middle American subduction system has resulted in an uneven distribution of magmatism/volcanism along strike. Here we construct a new three-dimensional shear-wave velocity model of the entire Middle American subduction system, using full-wave ambient noise tomography. Our model reveals significant variations of the oceanic plates along strike and down dip, in correspondence with either weakened or broken slabs after subduction. The northern and southern segments of the Cocos plate, including the Mexican flat slab subduction, are well imaged as high-velocity features, where a low-velocity mantle wedge exists and demonstrate a strong correlation with the arc volcanoes. Subduction of the central Cocos plate encounters a thick high-velocity feature beneath North America, which hinders the formation of a typical low-velocity mantle wedge and arc volcanoes. We suggest that the presence of slab tearing at both edges of the Mexican flat slab has been modifying the mantle flows, resulting in the unusual arc volcanism.

Natural Time Analysis of Seismicity within the Mexican Flat Slab before the M7.1 Earthquake on 19 September 2017

One of the most important subduction zones in the world is located in the Mexican Pacific Coast, where the Cocos plate inserts beneath the North American plate. One part of it is located in the Mexican Pacific Coast, where the Cocos plate inserts beneath the North American plate with different dip angles, showing important seismicity. Under the central Mexican area, such a dip angle becomes practically horizontal and such an area is known as flat slab. An earthquake of magnitude M7.1 occurred on 19 September 2017, the epicenter of which was located in this flat slab. It caused important human and material losses of urban communities including a large area of Mexico City. The seismicity recorded in the flat slab region is analyzed here in natural time from 1995 until the occurrence of this M7.1 earthquake in 2017 by studying the entropy change under time reversal and the variability β of the order parameter of seismicity as well as characterize the risk of an impending earthquake by applying the nowcasting method. The entropy change ΔS under time reversal minimizes on 21 June 2017 that is almost one week after the observation of such a minimum in the Chiapas region where a magnitude M8.2 earthquake took place on 7 September 2017 being Mexico’s largest quake in more than a century. A minimum of β was also observed during the period February–March 2017. Moreover, we show that, after the minimum of ΔS, the order parameter of seismicity starts diminishing, thus approaching gradually the critical value 0.070 around the end of August and the beginning of September 2017, which signals that a strong earthquake is anticipated shortly in the flat slab.

Multivariate analysis, mass balance techniques, and statistical tests as tools in igneous petrology: application to the Sierra de las Cruces volcanic range (Mexican Volcanic Belt)

Magmatic processes have usually been identified and evaluated using qualitative or semiquantitative geochemical or isotopic tools based on a restricted number of variables. However, a more complete and quantitative view could be reached applying multivariate analysis, mass balance techniques, and statistical tests. As an example, in this work a statistical and quantitative scheme is applied to analyze the geochemical features for the Sierra de las Cruces (SC) volcanic range (Mexican Volcanic Belt). In this locality, the volcanic activity (3.7 to 0.5 Ma) was dominantly dacitic, but the presence of spheroidal andesitic enclaves and/or diverse disequilibrium features in majority of lavas confirms the operation of magma mixing/mingling. New discriminant-function-based multidimensional diagrams were used to discriminate tectonic setting. Statistical tests of discordancy and significance were applied to evaluate the influence of the subducting Cocos plate, which seems to be rather negligible for the SC magmas in relation to several major and trace elements. A cluster analysis following Ward's linkage rule was carried out to classify the SC volcanic rocks geochemical groups. Finally, two mass-balance schemes were applied for the quantitative evaluation of the proportion of the end-member components (dacitic and andesitic magmas) in the comingled lavas (binary mixtures).