A model for the generation of intermediate and silicic igneous rocks is presented, based on experimental data and numerical modelling. The model is directed at subduction-related magmatism, but has general applicability to magmas generated in other plate tectonic settings, including continental rift zones. In the model mantlederived hydrous basalts emplaced as a succession of sills into the lower crust generate a deep crustal hot zone. Numerical modelling of the hot zone shows that melts are generated from two distinct sources; partial crystallization of basalt sills to produce residual H2O-rich melts; and partial melting of pre-existing crustal rocks. Incubation times between the injection of the first sill and generation of residual melts from basalt crystallization are controlled by the initial geotherm, the magma input rate and the emplacement depth. After this incubation period, the melt fraction and composition of residual melts are controlled by the temperature of the crust into which the basalt is intruded. Heat and H2O transfer from the crystallizing basalt promote partial melting of the surrounding crust, which can include meta-sedimentary and meta-igneous basement rocks and earlier basalt intrusions. Mixing of residual and crustal partial melts leads to diversity in isotope and trace element chemistry. Hot zone melts are H2O-rich. Consequently, they have low viscosity and density, and can readily detach from their source and ascend rapidly. In the case of adiabatic ascent the magma attains a super-liquidus state, because of the relative slopes of the adiabat and the liquidus. This leads to resorption of any entrained crystals or country rock xenoliths. Crystallization begins only when the ascending magma intersects its H2O-saturated liquidus at shallow depths. Decompression and degassing are the driving forces behind crystallization, which takes place at shallow depth on timescales of decades or less. Degassing and crystallization at shallow depth lead to large increases in viscosity and stalling of the magma to form volcano-feeding magma chambers and shallow plutons. It is proposed that chemical diversity in arc magmas is largely acquired in the lower crust, whereas textural diversity is related to shallow-level crystallization.

/pdf/the-genesis-of-intermediate-and-silicic-magmas-in-deep-3g95p0zv3p.pdf

The Genesis of Intermediate and Silicic Magmas in Deep Crustal Hot Zones

Subduction zones: subduction zones

[1] New thermal-petrologic models of subduction zones are used to test the hypothesis that intermediate-depth intraslab earthquakes are linked to metamorphic dehydration reactions in the subducting oceanic crust and mantle. We show that there is a correlation between the patterns of intermediate-depth seismicity and the locations of predicted hydrous minerals: Earthquakes occur in subducting slabs where dehydration is expected, and they are absent from parts of slabs predicted to be anhydrous. We propose that a subductingoceanicplatecanconsistoffourpetrologicallyandseismicallydistinctlayers:(1) hydrated, fine-grained basaltic upper crust dehydrating under equilibrium conditions and producing earthquakes facilitated by dehydration embrittlement; (2) coarse-grained, locally hydrated gabbroic lower crust that produces some earthquakes during dehydration but transformschieflyaseismicallytoeclogiteatdepthsbeyondequilibrium;(3)locallyhydrated uppermost mantle dehydrating under equilibrium conditions and producing earthquakes; and (4) anhydrous mantle lithosphere transforming sluggishly and aseismically to denser minerals. Fluid generated through dehydration reactions can move via at least three distinct flowpaths:percolationthroughlocal,transient,reaction-generatedhigh-permeabilityzones; flow through mode I cracks produced by the local stress state; and postseismic flow through fault zones. INDEX TERMS: 7218 Seismology: Lithosphere and upper mantle; 7230 Seismology: Seismicity and seismotectonics; 8123 Tectonophysics: Dynamics, seismotectonics; 8135 Tectonophysics: Evolution of the Earth: Hydrothermalsystems (8424); 3660 Mineralogyand Petrology: Metamorphicpetrology;

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2001JB001129

Subduction factory 2. Are intermediate-depth earthquakes in subducting slabs linked to metamorphic dehydration reactions?

We have created a digital age grid of the ocean floor with a grid node interval of 6 arc min using a self-consistent set of global isochrons and associated plate reconstruction poles. The age at each grid node was determined by linear interpolation between adjacent isochrons in the direction of spreading. Ages for ocean floor between the oldest identified magnetic anomalies and continental crust were interpolated by estimating the ages of passive continental margin segments from geological data and published plate models. We have constructed an age grid with error estimates for each grid cell as a function of (1) the error of ocean floor ages identified from magnetic anomalies along ship tracks and the age of the corresponding grid cells in our age grid, (2) the distance of a given grid cell to the nearest magnetic anomaly identification, and (3) the gradient of the age grid: i.e., larger errors are associated with high age gradients at fracture zones or other age discontinuities. Future applications of this digital grid include studies of the thermal and elastic structure of the lithosphere, the heat loss of the Earth, ridge-push forces through time, asymmetry of spreading, and providing constraints for seismic tomography and mantle convection models.

Digital isochrons of the world's ocean floor

[1] High-pressure metamorphic rocks provide evidence that in subduction zones material can return from depths of more than 100 km to the surface. The pressure-temperature paths recorded by these rocks are variable, mostly revealing cooling during decompression, while the time constraints are generally narrow and indicate that the exhumation rates can be on the order of plate velocities. As such, subduction cannot be considered as a single pass process; instead, return flow of a considerable portion of crustal and upper mantle material must be accounted for. Our numerical simulations provide insight into the self-organizing large-scale flow patterns and temperature field of subduction zones, primarily controlled by rheology, phase transformations, fluid budget, and heat transfer, which are all interrelated. They show the development of a subduction channel with forced return flow of low-viscosity material and progressive widening by hydration of the mantle wedge. The large-scale structures and the array of pressure-temperature paths obtained by these simulations favorably compare to the record of natural rocks and the structure of high-pressure metamorphic areas.

https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2002TC001406

Exhumation of high-pressure metamorphic rocks in a subduction channel: A numerical simulation

The geometry of the subducted Rivera and Cocos plates beneath the North American plate in southern Mexico was determined based on the accurately located hypocenters of local and teleseismic earthquakes. The hypocenters of the teleseisms were relocated, and the focal depths of 21 events were constrained using a body wave inversion scheme. The suduction in southern Mexico may be approximated as a subhorizontal slab bounded at the edges by the steep subduction geometry of the Cocos plate beneath the Caribbean plate to the east and of the Rivera plate beneath North America to the west. The dip of the interplate contact geometry is constant to a depth of 30 km, and lateral changes in the dip of the subducted plate are only observed once it is decoupled from the overriding plate. On the basis of the seismicity, the focal mechanisms, and the geometry of the downgoing slab, southern Mexico may be segmented into four regions : (1) the Jalisco region to the west, where the Rivera plate subducts at a steep angle that resembles the geometry of the Cocos plate beneath the Caribbean plate in Central America ; (2) the Michoacan region, where the dip angle of the Cocos plate decreases gradually toward the southeast, (3) the Guerrero-Oaxaca region, bounded approximately by the onshore projection of the Orozco and O'Gorman fracture zones, where the subducted slab is almost subhorizontal and underplates the upper continental plate for about 250 km, and (4) the southern Oaxaca and Chiapas region, in southeastern Mexico, where the dip of the subduction gradually increases to a steeper subduction in Central America. These drastic changes in dip do not appear to take place on tear faults, suggesting that smooth contortions accommodate these changes in geometry. The inferred 80 and 100 km depth contours of the subducted slab lie beneath the southern front of the Trans-Mexican Volcanic Belt, suggesting it is directly related to the subduction. Thus the observed nonparallelism with the Middle American Trench is apparently due to the changing geometry of the Rivera and Cocos plates beneath the North American plate in southern Mexico, and not to zones of weakness in the crust of the North American plate as some authors have suggested.

/pdf/shape-of-the-subducted-rivera-and-cocos-plates-in-southern-2f6c8m4sb3.pdf

Shape of the subducted Rivera and Cocos plates in southern Mexico: Seismic and tectonic implications

The Central Andes are the Earth's highest mountain belt formed by ocean-continent collision. Most of this uplift is thought to have occurred in the past 20 Myr, owing mainly to thickening of the continental crust, dominated by tectonic shortening. Here we use P-to-S (compressional-to-shear) converted teleseismic waves observed on several temporary networks in the Central Andes to image the deep structure associated with these tectonic processes. We find that the Moho (the Mohorovicic discontinuity--generally thought to separate crust from mantle) ranges from a depth of 75 km under the Altiplano plateau to 50 km beneath the 4-km-high Puna plateau. This relatively thin crust below such a high-elevation region indicates that thinning of the lithospheric mantle may have contributed to the uplift of the Puna plateau. We have also imaged the subducted crust of the Nazca oceanic plate down to 120 km depth, where it becomes invisible to converted teleseismic waves, probably owing to completion of the gabbro-eclogite transformation; this is direct evidence for the presence of kinetically delayed metamorphic reactions in subducting plates. Most of the intermediate-depth seismicity in the subducting plate stops at 120 km depth as well, suggesting a relation with this transformation. We see an intracrustal low-velocity zone, 10-20 km thick, below the entire Altiplano and Puna plateaux, which we interpret as a zone of continuing metamorphism and partial melting that decouples upper-crustal imbrication from lower-crustal thickening.

Subduction and collision processes in the Central Andes constrained by converted seismic phases

A critical reappraisal of great historical interplate earthquakes in the occidental margin of South America, including southern Peru and northern Chile, is carried out. A spacetime distribution of the earthquakes associated to the seismotectonics regions defined by the rupture zones of the greatest events (1868, Mw = 8.8 and 1877, Mw = 8.8) is obtained. Both regions are seismic gaps that are in the maturity state of their respective earthquake cycles. The region associated to the 1868 earthquake presents a notable seismic quiescence in the present century.

Reappraisal of great historical earthquakes in the Northern Chile and Southern Peru seismic gaps

A great earthquake (surface-wave magnitude, 7.8) occurred along the coast of central Chile on 3 March 1985, causing heavy damage to coastal towns. Intense foreshock activity near the epicenter of the main shock occurred for 11 days before the earthquake. The aftershocks of the 1985 earthquake define a rupture area of 170 by 110 square kilometers. The earthquake was forecast on the basis of the nearly constant repeat time (83 ± 9 years) of great earthquakes in this region. An analysis of previous earthquakes suggests that the rupture lengths of great shocks in the region vary by a factor of about 3. The nearly constant repeat time and variable rupture lengths cannot be reconciled with time- or slip-predictable models of earthquake recurrence. The great earthquakes in the region seem to involve a variable rupture mode and yet, for unknown reasons, remain periodic. Historical data suggest that the region south of the 1985 rupture zone should now be considered a gap of high seismic potential that may rupture in a great earthquake in the next few tens of years.

https://science.sciencemag.org/content/sci/233/4762/449.full.pdf

Mario Pardo

Papers

Shape of the subducted Rivera and Cocos plates in southern Mexico: Seismic and tectonic implications

Subduction and collision processes in the Central Andes constrained by converted seismic phases

Reappraisal of great historical earthquakes in the Northern Chile and Southern Peru seismic gaps

The 1985 central chile earthquake: a repeat of previous great earthquakes in the region?

Seismotectonic and stress distribution in the central Chile subduction zone