Numerical models of the magmatic processes induced by slab breakoff

After the onset of continental collision, magmatism often persists for tens of millions of years, albeit with a different composition, in reduced volumes, and with a more episodic nature and more widespread spatial distribution, compared to normal arc magmatism. Kinematic modelling studies have suggested that slab breakoff can account for this post-collisional magmatism through the formation of a slab window and subsequent heating of the overriding plate and decompression melting of upwelling asthenosphere, particularly if breakoff occurs at depths shallower than the overriding plate. To constrain the nature of any melting and the geodynamic conditions required, we numerically model the collision of two continental plates following a period of oceanic subduction. A thermodynamic database is used to determine the (de)hydration reactions and occurrence of melt throughout this process. We investigate melting conditions within a parameter space designed to generate a wide range of breakoff depths, timings and collisional styles. Under most circumstances, slab breakoff occurs deeper than the depth extent of the overriding plate; too deep to generate any decompressional melting of dry upwelling asthenosphere or thermal perturbation within the overriding plate. Even if slab breakoff is very shallow, the hot mantle inflow into the slab window is not sustained long enough to sufficiently heat the hydrated overriding plate to cause significant magmatism. Instead, for relatively fast, shallow breakoff we observe melting of asthenosphere above the detached slab through the release of water from the tip of the heating detached slab. Melting of the subducted continental crust during necking and breakoff is a more common feature and may be a more reliable indicator of the occurrence of breakoff. We suggest that magmatism from slab breakoff alone is unable to explain several of the characteristics of post-collisional magmatism, and that additional geodynamical processes need to be considered when interpreting magmatic observations

Overriding Plate Thickness as a Controlling Factor for Trench Retreat Rates in Narrow Subduction Zones

Slab width is a significant factor in controlling subduction zone dynamics, particularly the retreat velocities, which tend to decrease with wider slabs. However, observations of natural narrow subduction zones reveal no correlation between slab width and trench velocities. This suggests that other factors may exert a greater influence. In this study, we employ 3D numerical subduction models to systematically assess the impact of slab width, strength of slab coupling to the lateral plate (LP), and overriding plate (OP) thickness on trench kinematics and geometry. Our models focus on narrow slabs (400–1,200 km), and the results demonstrate that, in the case of narrow subduction zones, the slab width has little effect on trench migration rates and the viscous coupling at the lateral slab edge is only important for very narrow subduction zones (≤800 km). Conversely, the OP thickness emerges as a crucial factor, with increasing plate thickness leading to a strong decrease in trench velocities. These findings provide an explanation for the observed trench velocities in natural narrow subduction zones, where an inverse relationship with OP thickness is evident. Furthermore, our study reveals that not only slab width, but also the OP thickness and the slab coupling to the LP, significantly influence trench geometry. Strong lateral coupling promotes the formation of concave trench geometries, while thick overriding plates favor the development of “w”‐shaped geometries. Overall, a comprehensive understanding of subduction processes necessitates considering the interplay between slab width, OP thickness, and slab coupling to the LP

The role of pre-existing structures during rifting, continental breakup and transform system development, offshore West Greenland

Continental breakup between Greenland and North America produced the small oceanic basins of the Labrador Sea and Baffin Bay, which are connected via the Davis Strait, a region mostly comprised of continental crust. This study contributes to the debate regarding the role of pre-existing structures on rift development in this region using seismic reflection data from the Davis Strait data to produce a series of seismic surfaces, isochrons and a new offshore fault map from which three normal fault sets were identified as (i) NE-SW, (ii) NNW-SSE and (iii) NW-SE. These results were then integrated with plate reconstructions and onshore structural data allowing us to build a two-stage conceptual model for the offshore fault evolution in which basin formation was primarily controlled by rejuvenation of various types of pre-existing structures. During the first phase of rifting between at least Chron 27 (ca. 62 Ma; Palaeocene), but potentially earlier, and Chron 24 (ca. 54 Ma; Eocene) faulting was primarily controlled by pre-existing structures with oblique normal reactivation of both the NE-SW and NW-SE structural sets in addition to possible normal reactivation of the NNW-SSE structural set. In the second rifting stage between Chron 24 (ca. 54 Ma; Eocene) and Chron 13 (ca. 35 Ma; Oligocene), the sinistral Ungava transform fault system developed due to the lateral offset between the Labrador Sea and Baffin Bay. This lateral offset was established in the first rift stage possibly due to the presence of the Nagssugtoqidian and Torngat terranes being less susceptible to rift propagation. Without the influence of pre-existing structures the manifestation of deformation cannot be easily explained during either of the rifting phases. Although basement control diminished into the post-rift, the syn-rift basins from both rift stages continued to influence the location of sedimentation possibly due to differential compaction effects. Variable lithospheric strength through the rifting cycle may provide an explanation for the observed diminishing role of basement structures through time

GEMSToolbox: A novel modelling tool for rapid screening of mines for geothermal heat extraction

Mine water geothermal heat extraction is a promising technology to provide long-term, zero carbon heating in former coal mining regions. However to allow the technology to develop further, modelling tools which are faster than 3D coupled-process simulators and more site specific than analytical solutions need to be developed. A novel modelling tool, GEMSToolbox, was designed for assessing the feasibility of disused and flooded mine workings for such geothermal heat extraction. The basis of GEMSToolbox is built upon tested, computationally efficient modelling methods, and combines efficient solvers for water flow through the mine workings with fast solution methods for heat exchange between the water in the mine galleries and the surrounding rock mass. It expands and improves on these methods by (1) allowing for arbitrarily complex multi-seam mine geometries, and (2) addressing potential thermal interaction between nearby mine galleries using a novel geometric weighting technique. Fast calculation allows for wide parameter investigation studies that is required given the often uncertain state of disused mine systems. This makes the tool ideally suitable for the feasibility stage of a project for which site-specific yet computationally efficient alternative tools are currently lacking. The tool is demonstrated on a case study using plans from a real mine system and shows how it can be used to evaluate the long-term thermal performance of a mine water geothermal heat scheme

The thinning of subcontinental lithosphere: The roles of plume impact and metasomatic weakening

Geologically rapid (tens of Myr) partial removal of thick continental lithosphere is evident beneath Precambrian terranes, such as North China Craton, southern Africa, and the North Atlantic Craton,and has been linked with thermomechanical erosion by mantle plumes. We performed numerical experiments with realistic viscosities to test this hypothesis and constrain the most important parameters that influence cratonic lithosphere erosion. Our models indicate that the thermomechanical erosion by a plume impact on typical Archean lithospheric mantle is unlikely to be more effective than long-term erosion from normal plate-mantle interaction. Therefore, unmodified cratonic roots that have been stable for billions of years will not be significantly disrupted by the erosion of a plume event. However, the buoyancy and strength of highly depleted continental roots can be modified by fluid-melt metasomatism, and our models show that this is essential for the thinning of originally stable continental roots. The long-term but punctuated history of metasomatic enrichment beneath ancient continents makes this mode of weakening very likely. The effect of the plume impact is to speed up the erosion significantly and help the removal of the lithospheric root to occur within tens of Myr if affected by metasomatic weakening

Variable water input controls evolution of the Lesser Antilles volcanic arc

Oceanic lithosphere carries volatiles, notably water, into the mantle through subduction at convergent plate boundaries. This subducted water exercises control on the production of magma, earthquakes, formation of continental crust and mineral resources. Identifying different potential fluid sources (sediments, crust and mantle lithosphere) and tracing fluids from their release to the surface has proved challenging1. Atlantic subduction zones are a valuable endmember when studying this deep water cycle because hydration in Atlantic lithosphere, produced by slow spreading, is expected to be highly non-uniform2. Here, as part of a multi-disciplinary project in the Lesser Antilles volcanic arc3, we studied boron trace element and isotopic fingerprints of melt inclusions. These reveal that serpentine—that is, hydrated mantle rather than crust or sediments—is a dominant supplier of subducted water to the central arc. This serpentine is most likely to reside in a set of major fracture zones subducted beneath the central arc over approximately the past ten million years. The current dehydration of these fracture zones coincides with the current locations of the highest rates of earthquakes and prominent low shear velocities, whereas the preceding history of dehydration is consistent with the locations of higher volcanic productivity and thicker arc crust. These combined geochemical and geophysical data indicate that the structure and hydration of the subducted plate are directly connected to the evolution of the arc and its associated seismic and volcanic hazards