Abstract
Surface processes and magmatism condition the structural evolution of continental rifts and passive margins through mechanical and thermal effects on the lithosphere rheology. However, their inter-relationships in extensional settings are largely unknown. Here, I use coupled thermo-mechanical geodynamic and landscape evolution numerical modeling to assess the links between erosion of rift shoulders, sedimentation within the rift basin and extensional rock melting. Results suggest that, when the crust is thinner than ~40 km, the extension rate is slower than ~2 cm/yr and the mantle potential temperature is below ~1230 °C, efficient surface processes may double crustal melting by Moho lowering and inhibit mantle decompression melting by ~50% through sediment loading within the rift basin. It is thus likely that surface processes significantly influenced the magmatic activity of a number of extensional settings worldwide – e.g. the Mediterranean, the Gulf of California, the Iberia-Newfoundland margin, and the South China Sea. Because magmatism and surface processes affect jointly the geological carbon cycle, the surface processes forcing on extensional rock melting investigated here involves an additional means of linkage between plate tectonics and climate changes.
Similar content being viewed by others
Introduction
Continental extension, rifting and breakup involve the formation of thick sedimentary basins and escarpments at rift flanks that rise even a few kilometres above sea level. Progressive thinning and cooling of a heated lithosphere, strain localization, small-scale mantle convection and flexural forces are classically considered as some of the most important tectonic processes controlling the evolution of these prominent topographic features1,2,3. In turn, acting at similar rates to those of tectonic processes and generating stresses in the order of hundreds of MPa, sedimentation into the rift basin of the material eroded from the uplifted flanks affects the rheological response of the lithosphere to the tectonic forcing4,5,6,7. Enhanced basin subsidence by outward flow of the lower viscous crust due to sediment loading and enhanced uplift of the rift shoulders by surface unloading due to erosion are amongst the most relevant rheological implications of surface processes during rifting4 (Fig. 1).
Magmatism is a distinctive feature of continental lithospheric stretching8,9. Depending on the volumes of magma produced, two end-member types of rifted margins can be defined: volcanic or magma-poor10. At volcanic rifted margins, voluminous magmas erupt or intrude crustal units over short times, usually during continental breakup11,12, although the mechanisms responsible for the production of large volumes of magmas and the active contribution of mantle plumes to rifting are controversial13,14. Magma-poor rifted margins are characterized by wide domains of extended lithosphere with local magmatic bodies intruding syn-rift and especially post-breakup units10,15,16. In both cases, because the stresses required to rupture a typical continental lithosphere are estimated to be higher than those generated by extensional tectonics17,18, lithospheric heating/weakening and magma production likely provide a substantial contribution to lithospheric rupturing19,20,21. For instance, the magma migration path is expected to follow fractures along forming fault zones22,23, and fluid supply will then increase the pore fluid pressure, thereby lowering the plastic yield strength of fractured rocks and further localising the strain along weakening fault zones24,25,26.
While the relationships between lithospheric strain and surface processes and between lithospheric strain and magmatism are classically the subject of intense research4,27,28,29,30,31,32, the links between surface processes and rock melting across extensional settings are poorly constrained (Fig. 1). Yet, peaks of igneous activity due to enhanced mantle decompression melting have been ascribed to surface unloading by ongoing deglaciation33,34,35 and associated erosion36 or sea level lowering37,38. Modifications of the surface topography during the structural evolution of a magmatic province may affect the stress pattern within the elastic upper crust, with associated effects on the propagation of dykes39, the exsolution of magmatic gasses, and the pressure/overpressure field40. These factors co-determine the probability of magmatic, volcanic, and degassing events to occur as well as the amount of volatiles released into the ocean and atmosphere41. Because carbon degassing from extensional settings influenced climate changes at geological timescales42,43 and climate exerts a primary control on surface processes27, assessing likely relationships between surface processes and rock melting in rifting contexts is timely and key to our understanding of the geological carbon cycle and the surface-deep Earth processes coupling44.
In this study, I use coupled thermo-mechanical and landscape evolution numerical modeling to investigate the role of erosion and sedimentation in affecting the magmatic activity during continental extension and rifting. Given a set of experimentally/petrologically-determined rheological and partial rock melting relationships and parameters, one may simulate a range of plausible continental rift histories accounting for variable prescribed extension rates, crustal thickness, mantle potential temperature (sensu ref. 45) and erosion/sedimentation rates. The comparison between models enables assessing under which conditions surface processes can affect extensional rock melting as well as the mechanisms and feedbacks involved.
The numerical model
The main question addressed here is: can the surface mass redistribution by erosion and sedimentation affect partial melting of a stretching continental lithosphere? The general model includes three main components: (1) a surface process model, (2) a rheological lithospheric model, and (3) a model for partial rock melting. The evolution of the surface loads due to erosion and sedimentation is based on empirically determined laws of the surface transport of geological materials46,47. Rheological and partial melting laws/parameters are based on experimental rock mechanics and petrology48,49. The different components of the model are described in more detail hereafter.
Surface processes
The simplest representation of surface processes consists in the removal of rocks from uplifted regions (erosion), transport of the eroded material to the nearby basins and deposition within them (sedimentation). When integrated over long timescales (>105 yr) and large length scales (>105 m), these processes can be jointly reproduced by a linear downslope diffusion equation applied to the evolving topography in response to the tectonic strain4,7,47,50. The sediment flux at the surface, \({{\rm{q}}}_{s}\), is related to the local slope, \(\nabla z\), by \({{\rm{q}}}_{s}=k\nabla z\). The effective diffusivity, \(k\), given by the velocity of transport of the eroded material multiplied by its thickness, is a measure of the efficiency of surface processes. The assumption of mass conservation leads to the linear diffusion equation for erosion and sedimentation,
where \(t\) is time. Equation 1 can be solved numerically through integration on a discrete topography with constant-elevation boundary conditions4,50 (Fig. 2), thereby simulating the smoothing of the landscape by the cumulative effect of erosion and sediment deposition46,47. At each time step, the modelled landscape is updated for the effects of erosion/sedimentation and the associated surface load changes are computed.
Rheology
This model component accounts for the visco-elasto-plastic rheology of the lithosphere based on rock mechanics data for a quartz-dominated crust and an olivine-dominated mantle48. This implies a power law stress and exponential temperature dependence of the strain rate within the ductile part of the lithosphere. The elastic behaviour is given by linearly related stress and strain assuming typical Young’s modulus, \(\mu \), and Poisson’s ratio, \(\nu \), while a linear dependence between the brittle rock strength and pressure is assumed for plastic deformation (Table 1).
The structural changes of an evolving lithosphere are governed by the mass, momentum and energy conservation equations,
where \(\rho \) is the density, \({\boldsymbol{v}}\) is the velocity tensor, \({\boldsymbol{\sigma }}\) is the stress tensor, \(g\) is the acceleration due to gravity, \({C}_{p}\) is the specific heat capacity, \(T\) is temperature, \(c\) is the thermal conductivity, and \({H}_{r}\), \({H}_{s}\), \({H}_{a}\), and \({H}_{l}\) are the radiogenic, shear, adiabatic and latent heat production per unit volume, respectively51. Values for \({H}_{r}\) are listed in Table 1. \({H}_{s}\), related to dissipation of the mechanical energy during irreversible deformation, is calculated as \({H}_{s}={\sigma }_{ij}^{\text{'}}\dot{{\varepsilon }_{ij}^{\text{'}}}\), where \({\boldsymbol{\sigma }}{\boldsymbol{{\prime} }}\) is the deviatoric stress tensor, \(\dot{{\boldsymbol{\varepsilon }}{\boldsymbol{{\prime} }}}\) is the deviatoric strain rate tensor, i and j are coordinate indices (x, y) and repeated ij indices denotes summation. The adiabatic heat production, related to changes in pressure, is calculated as \({H}_{a}=T\alpha \frac{DP}{Dt}\), where \(\,\alpha \) is the thermal expansion. \({H}_{l}\), due to rock melting and crystallisation, is accounted for by computing the effective heat capacity, \(C{p}_{eff}\), and thermal expansion, \({\alpha }_{eff}\), of partially molten rocks (0 < \(\xi \) ≤ 1, see following) as52:
where \({Q}_{l}\) is the lithology-related latent heating (Table 1) and \(\xi \) is the volumetric fraction of melt computed as described hereafter.
Partial rock melting
The numerical model allows for partial melting and crystallisation of magma in the pressure-temperature domain between the wet solidus and dry liquidus of corresponding rocks (Table 1)49,53,54. At constant pressure, \(\xi \) is assumed to increase linearly with temperature according to the relations52.
where \({T}_{s}\) and \({T}_{l}\) are the solidus and liquidus of the considered rock, respectively. The effective density, \({\rho }_{eff}\), of partially molten rocks is then calculated as
where \({\rho }_{s}^{0}\) and \({\rho }_{l}^{0}\) are the standard densities of solid and molten rocks, respectively, and \({\rho }_{s}={\rho }_{0}[1+\beta (P-{P}_{0})]\times [1-\alpha (T-{T}_{0})]\) is the density of solid rocks at given P-T conditions (where \(\beta \) is the compressibility, P is pressure, and \({\rho }_{0}\), \({P}_{0}\) and \({T}_{0}\) are the density, pressure and temperature of rocks at surface conditions) (Table 1). The effective viscosity, \(\eta \), of partially molten rocks with \(\xi \) > 0.1 is assigned a low constant value of 1016 Pa s, which is an oversimplification, but it allows obviating numerical hurdles due to too high viscosity jumps between different materials55.
Integrated model components, reference setup and boundary conditions
Equations 1, 2 and 5 are coupled via stress, temperature and velocity continuity conditions, adopting the finite-differences approximation scheme and a fully staggered bi-dimensional grid, following the approach described in refs. 25 and55. The initial domain (Fig. 3) measures 400 × 300 km in the x and y dimensions, resolved by 161 × 61 grid points respectively, distributed on an irregular Eulerian grid that accounts for a resolution of 2 km along both directions in the central-upper part of the model. 400x300 Lagrangian markers are randomly distributed in the x and y dimensions and used for advecting the material properties. The material properties carried by Lagrangian markers are then interpolated onto the Eulerian grid via a 4th order Runge-Kutta interpolation scheme. The reference model includes a 100 km thick continental lithosphere and a 35 km thick continental crust, consistently with Moho depths commonly observed at continental margins in global crustal models56. The velocity boundary conditions are free slip at all boundaries (x = 0 and x = 400 km; y = 0 and y = 300 km). The left and right boundaries (x = 0 and x = 400 km) also account for x-parallel velocities, which define the extension rate within the model being equally distributed on the two boundaries so that the total extension rate measures up to a few cm/yr (Table 2), in agreement with common plate velocity values57. The lower boundary (y = 300 km) also accounts for y-parallel velocity to compensate for horizontal extension and ensure global mass conservation. The top surface of the lithosphere is calculated dynamically as an internal free surface through a 10 km thick layer of “sticky air”55. The initial temperature gradient in the asthenospheric mantle is 0.4 °C/km (adiabatic)51. The thermal boundary conditions are 0 °C for the upper boundary, nil horizontal heat flux across the vertical boundaries, and temperatures between 1320–1420 °C at the lower model boundary in order to account for different mantle potential temperatures (see Table 2 and a discussion about plausible mantle potential temperature in section 4). Except for the initial thermal state, there are no temperature conditions imposed within the model domain throughout the simulations. A seed of week material with 2 km radius is imposed at the centre of the model domain at the Moho to initiate the lithospheric rupture. At each time step, which is limited by the Courant criteria58, vertical topographic load changes are computed based on Eq. 1.
The parametric study focuses on the imposed extensional velocity, crustal thickness and mantle potential temperature (Table 2), which are arguably amongst the most important parameters controlling the structural and magmatic evolution of continental rifts4,20,21. Note that, because compositional changes between the lithospheric and asthenospheric mantle are not accounted for (Table 1), the parametric study on the mantle potential temperature also provides insights about the effects of changing the lithospheric thickness, inherently set by the depth of the 1300 °C isotherm59. For each combination of imposed conditions, one simulation with efficient and one with inefficient surface processes (controlled by the effective diffusivity, \(k\), in Eq. 1) are performed. A total of 36 numerical experiments were performed and analysed, results are summarized in the following.
Results
The evolution of the reference model is shown in Fig. 4. The details and timing of the model evolution depend upon the assumed thermal and mechanical parameters (see following), but the overall sequence of events, coherent with that of “type II” continental rifts following the classification of ref. 60 is robust and includes: (1) lithospheric stretching, early breakup of the mantle lithosphere and asthenospheric upwelling to Moho depths, (2) prolonged crustal stretching to continental breakup and oceanization between the two newly formed continental passive margins. Faster lithospheric mantle than crustal thinning leads to heating of the Moho and the production of migmatitic core complexes and/or andesitic/dacitic effusive rocks by basal crustal melting in the early stages of the model evolution. Similar early extensional products were modeled before61 and observed, for instance, in the Rhodopes, the Vøring Plateau and the Basin and Range62,63,64,65. Crustal melting in latter two cases is driven by heating from basaltic intrusions due to active upwelling of hot asthenosphere63,65. Here, crustal melting rather occurs by syn-extensional diffusive heating from the asthenosphere, more similar to, for instance, the Eocene Rhodope migmatitic core complexes62. Regardless of the source of heat for extensional crustal melting, which may be multiple66, a first order consistency between modeled and observed magmatic products is the essential requirement for assessing a possible surface processes forcing. Later decompression melting of the asthenosphere begins prior to continental breakup and continues during syn- and post-breakup phases. Overall, the topography produced by the model is consistent with that of rift flanks and basins observed in natural case studies, where the rift basin is filled at the expenses of the eroded rift shoulders (Figs. 4 and 5).
The results of the parametric study are summarized in Fig. 6. The estimated volume of melts generated during rifting at 1 cm/yr with the Moho at 35 km depth and mantle potential temperature of 1200 °C (i.e. reference model) is subject to strong variations depending on the redistribution of the surface masses (Fig. 6a). Early crustal melting at Moho levels is enhanced by efficient sediment delivery to the rift basin, which dampens crustal thinning, increases subsidence and facilitates basal heating due to asthenospheric upwelling through the thinning mantle lithosphere (Figs. 5a and 7a). In turn, a greater crustal thickness and higher Moho temperature due to efficient surface processes promote ductile straining of the crust with respect to brittle rupturing, leading to delayed continental breakup (Fig. 6a). During later stages of the model evolution, mantle melting is inhibited by efficient sediment delivery to the rift basin (Fig. 6a) because an increasing sedimentary load dampens decompression of the upwelling asthenosphere, thereby reducing the portion of the geotherm above the mantle solidus/liquidus (Fig. 7b,c).
When the reference model setup (Fig. 6a) is modified to account for a higher mantle potential temperature (e.g. 1250 °C, Fig. 6b), faster extensional velocity (e.g. 2 cm/a, Fig. 6c), and thicker crust (e.g. 45 km, Fig. 6d), the maximum amount of melts produced increases by ~ five-, three- and six-fold respectively. Enhanced partial melting of both the crust and asthenosphere when a higher mantle potential temperature with respect to the reference model is accounted for (Fig. 6b) is expected, because this implies a warmer geotherm across the model domain. Doubling the extensional velocity with respect to the reference model leads to faster lithospheric rupturing (breakup after only ~8 Myr compared to ~17 Myr in the reference model), which in turn enhances mantle melt production by faster upwelling and decompression rates (Fig. 6c). Increased crustal thickness with respect to the reference model implies crustal material to greater depths and higher temperatures and, thus, enhanced crustal melt production (Fig. 6d). Similar trends in the modulation by surface processes to the amount of melts produced are observed when the mantle potential temperature is 1250 °C, the extension velocities are 2 cm/yr or the initial crustal thickness is 45 km (Fig. 6b–d). However, these trends appear smoothened in the former two cases (Fig. 6b,c) and significantly dampened in the latter case (Fig. 6d), when crustal melting dominates on mantle decompression melting and particularly voluminous melts are produced.
Discussion
Coupled thermo-mechanical and landscape evolution numerical modeling allows reproducing histories of continental rifting and partial melting consistent with those of previous studies and observations from natural settings4,60,67,68,69. It is thus possible to draw general relationships from the numerical results and evaluate their applicability to and implications for natural rifts. Without feedbacks between surface and subsurface processes, the topographic anomalies generated during continental rifting would be rapidly reset4. To a first approximation, the rift flanks and the underlying Moho rise by approximately the ~5/6th of the crustal thickness lost by erosion within a given time window and, by the same principle, the crust bends downward as the rift basin is filled, which generates more accommodation space for sediments4,70. This study, however, shows that the feedbacks between surface and deep Earth processes in rifting settings are not limited to topographic features, but also involve partial rock melting (Figs. 4 and 6–8). If thermal blanketing due to efficient rift basin filling may rise lower crustal temperatures by ~50–100 °C (ref. 4), thereby contributing to enhance crustal partial melting, sediment loading in the rift basin lowers the geotherm with respect to the crustal solidus, which dampens or even overrides the effect of thermal blanketing (Fig. 7). The competing effects of thermal blanketing and solidus lowering on crustal melting, however, is outpaced by flexural bending of the Moho in a stretching and warming lithosphere, which supplies additional crustal material at near/above-solidus conditions, ultimately enhancing the amount of crustal melts (Figs. 6 and 7). Lithospheric stretching at rates of a few cm/yr forces the asthenosphere to upwell and decompress at a similar rate (Figs. 4 and 5). Decompression partial melting of the asthenosphere, however, is dampened by surface loading due to efficient filling of the rift basin (Fig. 6). The dampening is a function of the rate of basin deepening/filling, the sediment density, and the surface-to-depth stress change transfer of the rift system. For the reference model setup, the crustal melt production is roughly doubled and the mantle melt production is reduced by ~50% when surface processes redistribute the surface masses so efficiently to nearly reset the topography through time (i.e. the rate of erosion/deposition is similar that of tectonically-controlled surface uplift/subsidence) (Figs. 5, 6 and 8). This upper bound estimate may be used as a reference to infer possible modulations by surface processes to rock melting in natural rifting settings, which, for a given erosion/deposition rate, appears to be inversely correlated to the extensional velocity, mantle potential temperature and initial Moho depth (Fig. 6). Particularly noteworthy is that, while faulting and hence brittle-plastic deformation plays an important role in setting the rate and amount of melt percolation along fracture zones22,71, melt percolation reduces the brittle-plastic strength of rocks through an increase in pore fluid pressure72. This feedback, which further localises deformation along weakening fault zones and facilitates the ascent of partially molten rocks through the lithosphere, is not taken into account in the numerical models presented here, which has two main implications. First, the reference upper bound estimate regarding the possible modulation of partial melting by surface processes (Fig. 8) is likely underestimated to some extent. Second, thermal blanketing and flexural stresses due to abundant sediment supply to the rift basin weaken the lithosphere and favour prolonged viscous stretching upon abrupt brittle/plastic rupturing4, thereby explaining the retardation of continental breakup (Fig. 6). However, since brittle-plastic strain localisation due to melts-rocks interactions is not taken into account, the influence of surface processes on the timing of lithospheric rupturing cannot be properly assessed. I anticipate that future studies will investigate the relative contributions of thermal blanketing and melts-rocks interactions modulated by surface processes in affecting the timing of continental breakup in extensional settings.
In terms of driving mechanisms, the surface processes forcing on rock melting quantified here is similar to the proposed deglacial or eustatic control on extensional rock melting in, for instance, Iceland, the Mediterranean, or the Australian-Antarctic ridge33,36,37,38. However, the improvement with respect to previous studies is that the parametric investigation performed with fully coupled landscape evolution and thermo-mechanical geodynamic models allows constraining a range of conditions in which the surface processes forcing on extensional rock melting is likely conspicuous. That is, when the crust is thinner than ~40 km, the extension rate is slower than ~2 cm/yr and the mantle potential temperature is below ~1230 °C (Fig. 6). While the first two conditions are commonly observed56,57, the global mean mantle potential temperature at mid-oceanic ridges was estimated as ~ 1360 °C (refs. 73,74,75). However, long wavelength temperature variations in the sub-lithospheric convective upper mantle are at least \(\pm \) ~ 200 °C (refs. 76,77,78). Therefore, plausible mantle potential temperatures range between ~ 1160–1560 °C (ref. 13,45). According to ref. 76,78 negative anomalies with respect to the mean value may reach up to ~ −250 °C and ~ −150 °C at “cold” mantle regions such as non-ridge domains or nearby subduction zones, respectively. Underneath the continents and far from mantle plumes, the mantle potential temperature is estimated as ~1200 °C, that is at least 100 °C cooler than that beneath the oceans79. The mantle geotherm below the Corinth Rift and nearby the Gulf of Lion rift system, for instance, may be cooler than the 1100 °C adiabat80, possibly due to the proximal Hellenic and Alpine-Apennine slabs. Estimated mantle temperatures for eastern North America and offshore of Iberia suggest mantle geotherms significantly below the 1280 °C mantle adiabat prior and during rifting, which is consistent with the magma-poor structure of the Iberia-Newfoundland margin79,80 ref. 81 estimate the modern average geotherm for the western United States as asymptotic to the 1300 °C adiabat due to the influx of oceanic asthenosphere from the Gulf of California. This scenario implies that the mantle potential temperature at the base of the lithosphere prior to oceanization was closer to 1200 °C than to 1300 °C79. These lines of evidence suggest syn-extensional mantle potential temperatures as low as ~ 1230 °C, or even cooler, in the Mediterranean region, the Iberia-Newfoundland margin and the Gulf of California. I propose that in these and possibly other extensional settings such as the South China Sea, surface processes may have affected partial melting by up to doubling the amount of melts produced (Figs. 6a and 8)82. I further speculate that the mass redistribution due to surface processes across these settings may have contributed to the rift geometry, migration, segmentation or failure, for instance through modulation of crustal or mantle melting and associated magma-assisted lithospheric rupturing19,72,83 or the modulation of inherited structures84,85. If correct, an important implication is that the role of climate-controlled surface processes in conditioning the surface expressions of plate tectonics are not limited to the modulation of the topography27,86, but also extend to continental drifting and, thus, the current global plate configuration.
Recent investigations show a correlation between the worldwide rifts or mid oceanic ridges length over the last ~200 Ma and atmospheric CO2 proxies42,43, in turn suggesting that extensional rock melting contributes to setting long term global climate trends. However, climate is the main driver of surface processes. Thus, the surface processes forcing on extensional rock melting investigated here may involve a previously unknown feedback mechanism between the Solid Earth and the Surface Earth, which likely conditioned the long term (multi-Myr) evolution of the Earth system. A decrease in solidus temperatures of peridotites in presence of CO287,88,89 would enhance this newly-recognized surface processes feedback on climate via extensional melting. Since the magma migration velocities and signal propagation rates in sedimentary systems commonly reach up to the tens of cm per year51,90, it is likely that the long term surface processes forcing on extensional rock melting addressed here also finds expression at shorter (e.g. centennial to multi-millenial) timescales. Milankovitch frequencies in tephra or bathymetric records worldwide37,91 corroborate this speculation. However, in order to validate these assertions, more advanced erosion-deposition physical laws including, for instance, fluvial and dispersed sediment transport, respectively dominant in channels and hillslopes92, are required. Accounting for variable rock erodibility, precipitation rates and river network geometries in space and time93 would further improve the description of the modifications of the surface load, in turn allowing to better constrain the surface processes forcing on extensional rock melting. Because unravelling the interactions between plate tectonics and climate through the geological carbon cycle is undoubtedly one of the major challenges in the Earth sciences44, future works that will address the many possible nuances in the first order relationships between surface processes and extensional rock melting identified here are warranted.
Conclusions
The numerical models presented here show that the surface mass redistribution by erosion and sediment deposition conditions rock melting across extensional settings. Prominent surface processes may double the amount of magma generated by crustal melting and reduce by ~50% the amount of magma generated by mantle decompression melting, if continental stretching is slower than ~2 cm/yr, the crust is thinner than ~40 km, and the mantle potential temperature is below ~1230 °C. The limitations the numerical model is subject to, particularly those regarding the influence of fluids on the plastic deformation of rocks, are such that these inferences may be underestimated and that the range of conditions in which surface processes significantly affect extensional rock melting may be wider. The implications of these findings are, at least, twofold. First, climate-controlled erosion, sediment transport and deposition, condition the structural evolution of extensional settings not only through effects on the topography and subsurface stress and strain field4, but also through effects on partial rock melting and associated magmatic and rheological changes. Second, a previously unrecognized mutual feedback between climate and tectonics exists and takes place through a control of surface processes on extensional rock melting and of rift magmatism on surface processes via climatic effects.
Data availability
No datasets were generated or analysed during the current study, but the codes used are available from the corresponding author on request.
References
Royden, L. & Keen, C. E. Rifting process and thermal evolution of the continental margin of Eastern Canada determined from subsidence curves. Earth Planet. Sci. Lett. 51, 343–361 (1980).
Watts, A. B. & Tome, M. Crustal structure and the mechanical properties of extended continental lithosphere in the Valencia through (western Mediterranean). J. Geol. Sot. London 149, 813–827 (1992).
Cloetingh, S. & Burov, E. B. Thermomechanical stmcture of European lithosphere: constraints from rheological profiles and EET estimates. Geophys. J. Int. 124, 695–723 (1996).
Burov, E. B. & Cloetingh, S. Erosion and rift dynamics: new thermomechanical aspects of post-rift evolution of extensional basins. Earth and Planetary Science Letters 150, 7–26 (1997).
England, P. & Richardson, S. W. The influence of erosion upon the mineral facies of rocks from different metamorphic environments. J. Geol. Sot. London 134, 201–213 (1977).
Erdős, Z., Huismans, R. S., van der Beek, P. & Thieulot, C. Extensional inheritance and surface processes as controlling factors of mountain belt structure. Journal of Geophysical Research: Solid Earth 119(112), 9042–9061 (2014).
Burov, E. & Poliakov, A. Erosion and rheology controls on synrift and postrift evolution: Verifying old and new ideas using a fully coupled numerical model. Journal of Geophysical Research: Solid Earth 106(1B8), 16461–16481 (2001).
Bailey, D. K. Lithosphere control of continental rift magmatism. Journal of the Geological Society 133(11), 103–106 (1977).
White, Rs et al. Magmatism at rifted continental margins. Nature 330(16147), 439 (1987).
Franke, D. Rifting, lithosphere breakup and volcanism: Comparison of magma-poor and volcanic rifted margins. Marine and Petroleum Geology, 43, 63–87 (2013).
R. S. White, Initiation of the Iceland Plume and Opening of the North Atlantic, In: Extensional Tectonics and Stratigraphy of the North Atlantic Margins, pp. 149-154, 1989.
Coffin, M. F. & Eldholm, O. Large igneous provinces: crustal structure, dimensions, and external consequences. Reviews of Geophysics 32, 1–36 (1994).
White, R. & McKenzie, D. Magmatism at rift zones: the generation of volcanic continental margins and flood basalts. Journal of Geophysical Research 94, 7685–7729 (1989).
Foulger, G. R. & Natland, J. H. Is “Hotspot” Volcanism a consequence of plate tectonics? Science 300, 921–922 (2003).
Whitmarsh, R. B., Manatschal, G. & Minshull, T. A. Evolution of magma-poor continental margins from rifting to seafloor spreading. Nature 413(16852), 150 (2001).
Manatschal, G. New models for evolution of magma-poor rifted margins based on a review of data and concepts from West Iberia and the Alps. International Journal of Earth Sciences 93(13), 432–466 (2004).
W. R. Buck, Consequences of Asthenospheric Variability in Continental Rifting, in Rheology and Deformation of the Lithosphere at Continental Margins (eds. Karner, G. D., Taylor, B., Driscoll, N. W. & Kohlstedt, D. L.) 1–30 (Columbia Univ. Press, New York)., 2004.
Bott, M. H. Sublithospheric loading and plate-boundary forces. Phil. Trans. R. Soc. Lond. 337, 83–93 (1991).
Kendall, J. M., Stuart, G. W., Ebinger, C. J., Bastow, I. D. & Keir, D. Magma-assisted rifting in Ethiopia. Nature 433, 146–148 (2005).
Lavecchia, A., Beekman, F., Clark, S. R. & Cloetingh, S. A. Thermo-rheological aspects of crustal evolution during continental breakup and melt intrusion: The Main Ethiopian Rift, East Africa. Tectonophysics 686, 51–62 (2016).
Lavecchia, A., Clark, S. R., Beekman, F., Cloetingh, S. A. & Burov, E. Thermal perturbation, mineral assemblages, and rheology variations induced by dyke emplacement in the crust. Tectonics 35(15), 1137–1152 (2016).
Connolly, J. A. & Podladchikov, Y. Y. Compaction-driven fluid flow in viscoelastic rock. Geodinamica Acta 11, 55–84 (1998).
Spence, D. A., Sharp, P. W. & Turcotte, D. L. Buoyancy-driven crack propagation: a mechanism for magma migration. Journal of Fluid Mechanics 174, 135–153 (1987).
Katz, R. F. Magma dynamics with the enthalpy method: benchmark solutions and magmatic focusing at mid-ocean ridges. Journal of Petrology 49, 2099–2121 (2008).
Gerya, T. & Yuen, D. Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems. Phys. Earth Planet. Inter. 163, 83–105 (2007).
Turcotte, D. L. Magma migration. Annual Review of Earth and Planetary Sciences 10(11), 397–408 (1982).
W. F. Ruddiman, Tectonic Uplift and Climate Change, Springer, 1997.
B. B. Smith, Uplift, erosion and stability: perspectives on long-term landscape development, Geological Society of London, 1999.
Avouac, J. P. & Burov, E. B. Erosion as a driving mechanism of intracontinental mountain growth. Journal of Geophysical Research - all series 101, 17–747 (1996).
Cloetingh, S. A., Ziegler, P. A. & Bogaard, P. J. TOPO-EUROPE: The geoscience of coupled deep Earth-surface processes. Global and Planetary Change 58, 1–118 (2007).
Schmeling, H. Dynamic models of continental rifting with melt generation. Tectonophysics 480(11-4), 33–47 (2010).
Corti, G. et al. Analogue modelling of continental extension: a review focused on the relations between the patterns of deformation and the presence of magma. Earth-Science Reviews 63(13-4), 169–247 (2003).
Jull, M. & McKenzie, D. The effect of deglaciation on mantle melting beneath Iceland. Journal of Geophysical Research 101(1B10), 21815–21828 (1996).
Huybers, P. & Langmuir, C. Feedback between deglaciation, volcanism, and atmospheric CO2. Earth and Planetary Science Letters 286, 479–491 (2009).
Schmidt, P. et al. Effects of present‐day deglaciation in Iceland on mantle melt production rates. Journal of Geophysical Research: Solid Earth 118(17), 3366–3379 (2013).
P. Sternai, L. Caricchi, S. Castelltort & J.-. D. Champagnac, Deglaciation and glacial erosion: a joint control on magma productivity by continental unloading, Geophysical Research Letters, vol. https://doi.org/10.1002/2015GL067285, (2016).
Crowley, J. W., Katz, R. F., Huybers, P., Langmuir, C. & Park, S.-H. Glacial cycles drive variations in the production of oceanic crust. Science 347(16227), 1237–1240 (2015).
Sternai, P. et al. Magmatic pulse driven by sea-level changes associated with the Messinian salinity crisis. Nature Geoscience https://doi.org/10.1038/NGEO3032, (2017).
Maccaferri, F., Bonafede, M. & Rivalta, E. A quantitative study of the mechanisms governing dike propagation, dike arrest and sill formation. Journal of Volcanology and Geothermal Research 208(11-2), 39–50 (2011).
Jellinek, A. M. & De Paolo, D. J. A model for the origin of large silicic magma chambers: precursors of caldera-forming eruptions. Bulletin of Volcanology 65(15), 363–381 (2003).
Albino, F., Pinel, V. & Sigmundsson, F. Influence of surface load variations on eruption likelihood: application to two Icelandic subglacial volcanoes, Grímsvötn and Katla. Geophysical journal international 181(13), 1510–1524 (2010).
Brune, S., Williams, S. E. & Müller, R. D. Potential links between continental rifting, CO2 degassing and climate change through time. Nature Geoscience 10, 941–948 (2017).
Van Der Meer, D. G. et al. Plate tectonic controls on atmospheric CO2 levels since the Triassic. PNAS 111(112), 4380–4385 (2014).
Sternai, P. et al. Magmatic Forcing of Cenozoic Climate? Journal of Geophysical Research - Solid Earth https://doi.org/10.1029/2018JB016460, (2019).
McKenzie, D. & Bickle, M. J. The volume and composition of melt generated by extension of the lithosphere. Journal of petrology 29(13), 625–679 (1988).
M. A. Carson & M. J. Kirkby, Hillslope Form and Processes, Cambridge Univ. Press, Cambridge, p. 475, 1972.
Gossman, H. Slope modelling with changing boundary conditions - effects of climate and lithology. Z. Geomorphol. N.F., Suppl. Bd. 25, 72–88 (1976).
G. Ranalli, Rheology of the Earth, Deformation and Flow Processes in Geophysics and Geodynamics, Chapman & Hall, 1995.
Hirschmann, M. M. Mantle solidus: experimental constraints and the effects of peridotite composition. Geochemistry Geophysics Geosystems 1, 1–21 (2000).
Gerya, T. V. & Yuen, D. A. Characteristics-based marker-in-cell method with conservative finite-differences schemes for modeling geological flows with strongly variable transport properties. Physics of the Earth and Planetary Interiors 140(14), 293–318 (2003).
D. L. Turcotte & G. Schubert, Geodynamics, Cambridge: Cambridge University Press, 2002.
Burg, J. P. & Gerya, T. V. The role of viscous heating in Barrovian metamorphism of collisional orogens: thermomechanical models and application to the Lepontine Dome in the Central Alps. Journal of Metamorphic Geology 23(12), 75–95 (2005).
P. C. Hess, Origin of Igneous Rocks, Harward University Press, London, UK, 1989.
W. Johannes The significance of experimental studies for the formation of migmatites, In: Migmatites (ed. Ashworth, J. R.),Blackie, Glasgow, UK, p. 36–85, 1985.
T. Gerya, Introduction to Numerical Geodynamic Modelling, Cambridge University Press, 2010.
Mooney, W. D., Laske, G. & Masters, T. G. Mooney, Laske and Masters, CRUST 5.1: a global crustal model at 5x5 degrees. Journal of Geophysical Research 103, 727–747 (1998).
R. D. Müller, et al., Ocean basin evolution and global-scale plate reorganization events since Pangea breakup, Annual Reviews of Earth and Planetary Science, vol. 44, 2016.
Courant, R. On the partial difference equations of mathematical physics. Mathematische Annalen 100, 32–74 (1928).
Clauser, C. & Huenges, E. Thermal conductivity of rocks and minerals. Rock physics and phase relations: a handbook of physical constants 3, 105–126 (1995).
Huismans, R. & Beaumont, C. Depth-dependent extension, two-stage breakup and cratonic underplating at rifted margins. Nature 473, 74–79 (2011).
Schenker, F. L., Gerya, T. & Burg, J. P. Bimodal behavior of extended continental lithosphere: Modeling insight and application to thermal history of migmatitic core complexes. Tectonophysics 579, 88–103 (2012).
Jones, C. E., Tarney, J., Baker, J. H. & Gerouki, F. Tertiary granitoids of Rhodope, northern Greece: Magmatism related to extensional collapse of the Hellenic Orogen? Tectonophysics 210(13-4), 295–314 (1992).
Eldholm, O. Magmatic tectonic evolution of a volcanic rifted margin. Marine Geology 102(11-4), 43–61 (1991).
Viereck, L. G. et al. Origin of the Palaeogene Vøring Plateau volcanic sequence. Geological Society, London, Special Publications 39(11), 69–83 (1988).
P. B. Gans, Synextensional magmatism in the Basin and Range province: A case study from the eastern Great Basin, Geological Society of America., vol. 233, 1989.
D. L. Whitney, C. Teyssier & O. Vanderhaeghe, Gneiss domes and crustal flow, Geological Society of America, vol. Special Paper 380, pp. 1-19, 2004.
Brune, S., Heine, C., Clift, P. D. & Pérez-Gussinyé, M. Rifted margin architecture and crustal rheology: Reviewing Iberia-Newfoundland, Central South Atlantic, and South China Sea. Marine and Petroleum Geology 79, 257–281 (2017).
G. Corti et al. Transition from continental break‐up to punctiform seafloor spreading: How fast, symmetric and magmatic, Geophysical Research Letters, vol. 30, 112, 2003.
Ziegler, P. A. & Cloetingh, S. Dynamic processes controlling evolution of rifted basins. Earth-Science Reviews 64(11-2), 1–50 (2004).
A. B. Watts, Isostasy and flexure of the lithosphere, Cambridge: Cambridge University Press, 2009.
Roche, V. et al. Emplacement of metamorphic core complexes and associated geothermal systems controlled by slab rollback. Earth and Planetary Science Letters 498, 322–333 (2018).
Clemens, J. & Mawer, C. K. Granitic magma transport by fracture propagation. Tectonophysics 204(13-4), 339–360 (1992).
C. Herzberg & P. D. Asimow, Petrology of some oceanic island basalts: PRIMELT2. XLS software for primary magma calculation, Geochemistry, Geophysics, Geosystems, vol. 9, 19, 2008.
Lee, C. T., Luffi, P., Plank, T., Dalton, H. & Leeman, W. P. Constraints on the depths and temperatures of basaltic magma generation on Earth and other terrestrial planets using new thermobarometers for mafic magmas. Earth and Planetary Science Letters 279(11-2), 20–33 (2009).
Putirka, K. Excess temperatures at ocean islands: Implications for mantle layering and convection. Geology 36(14), 283–286 (2008).
Anderson, D. L. The scale of mantle convection. Tectonophysics 284, 1–17 (1998).
Anderson, D. L. The thermal state of the upper mantle; no role for mantle plumes. Geophysical Research Letters 27(122), 3623–3626 (2000).
Kaula, W. Minimal upper mantle temperature variations consistent with observed heat flow and plate velocities. Journal of Geophysical Research: Solid Earth 88(1B12), 10323–10332 (1983).
Reston, T. J. & Morgan, J. P. Continental geotherm and the evolution of rifted margins. Geology 32(12), 133–136 (2004).
Goes, S., Govers, R. & Vacher, A. P. Shallow mantle temperatures under Europe from P and S wave tomography. Journal of Geophysical Research: Solid Earth 105(1B5), 11153–11169 (2000).
S. Goes & S. van der Lee, Thermal structure of the North American uppermost mantle inferred from seismic tomography, Journal of Geophysical Research, vol. 107, 1B3, https://doi.org/10.1029/2000JB000049, 2002.
Lei, C. et al. Structure and sediment budget of Yinggehai–Song Hong basin, South China Sea: Implications for Cenozoic tectonics and river basin reorganization in Southeast Asia. Tectonophysics 655, 177–190 (2015).
Corti, G. et al. Magma-induced strain localization in centrifuge models of transfer zones. Tectonophysics 348(14), 205–218 (2002).
G. Corti, J. van Wijk, S. Cloetingh & C. K. Morley, Tectonic inheritance and continental rift architecture: Numerical and analogue models of the East African Rift system, Tectonics, vol. 26, 16, 2007.
Brune, S., Corti, G. & Ranalli, G. Controls of inherited lithospheric heterogeneity on rift linkage: Numerical and analog models of interaction between the Kenyan and Ethiopian rifts across the Turkana depression. Tectonics 36(19), 1767–1786 (2017).
Molnar, P. & England, P. Late Cenozoic uplift of mountain ranges and global climate change: chicken or egg? Nature 346(16279), 29–34 (1990).
Wyllie, P. J. Effects of H2O and CO2 on magma generation in the crust and mantle. Journal of the Geological Society 134, 215–234 (1977).
Dasgupta, R., Hirschmann, M. M. & Smith, N. D. Water follows carbon: CO2 incites deep silicate melting and dehydration beneath mid-ocean ridges. Geology 35(12), 135–138 (2007).
Gudfinnsson, G. H. & Presnall, D. C. Continuous gradations among primary carbonatitic, kimberlitic, melilititic, basaltic, picritic, and komatiitic melts in equilibrium with garnet lherzolite at 3-8 GPa. Journal of Petrology 46, 1645–1659 (2005).
B. W. Romans, S. Castelltort, J. A. Covault, A. Fildani & J. P. Walsh, Environmental signal propagation in sedimentary systems across timescales, Earth-Science Reviews, 2015.
S. Kutterolf, J. C. Schindlbeck, M. Jegen, A. Freundt & S. M. Straub, Milankovitch frequencies in tephra records at volcanic arcs: The relation of kyr-scale cyclic variations in volcanism to global climate changes, Quaternary Science Reviews, pp. 1-16, 2019.
Simpson, G. A dynamic model to investigate coupling between erosion, deposition, and three-dimensional (thin-plate) deformation. Journal of Geophysical Research 109, 1–12 (2004).
Ueda, K., Willett, S. D., Gerya, T. & Ruh, J. Geomorphological–thermo-mechanical modeling: Application to orogenic wedge dynamics. Tectonophysics 659, 12–30 (2015).
Acknowledgements
Pietro Sternai was supported by the Italian MIUR (Rita Levi Montalcini grant, DM 694–26/2017). This study is also part of the project MIUR – Dipartimenti di Eccellenza 2018–2022, Department of Earth and Environmental Sciences, University of Milano-Bicocca. Nadia Malaspina provided suggestions regarding the link between the surface processes forcing on extensional rock melting and climate. I thank three anonymous reviewers for providing careful revisions.
Author information
Authors and Affiliations
Contributions
P.S. conceived the scientific question, performed the study and wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sternai, P. Surface processes forcing on extensional rock melting. Sci Rep 10, 7711 (2020). https://doi.org/10.1038/s41598-020-63920-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-63920-w
This article is cited by
-
Towards a process-based understanding of rifted continental margins
Nature Reviews Earth & Environment (2023)
-
Magmatic evolution of the migrating central Urumieh–Dokhtar arc, Iran: implications for magma production
International Journal of Earth Sciences (2023)
-
3D geodynamic-geomorphologic modelling of deformation and exhumation at curved plate boundaries: Implications for the southern Alaskan plate corner
Scientific Reports (2022)
-
Feedbacks between sea-floor spreading, trade winds and precipitation in the Southern Red Sea
Nature Communications (2022)
-
Climatic control on the location of continental volcanic arcs
Scientific Reports (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.