Issue
BSGF - Earth Sci. Bull.
Volume 192, 2021
Special Issue Orogen lifecycle: learnings and perspectives from Pyrenees, Western Mediterranean and analogues
Article Number 45
Number of page(s) 13
DOI https://doi.org/10.1051/bsgf/2021036
Published online 18 October 2021

© L. Struth et al., Published by EDP Sciences 2021

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction to drainage opening and transient topography

The study of fluvial network rearrangements provides a key to understand past and future landscape evolution. River profiles respond to tectonic or climatic perturbations on water supply, slope or base level changes (e.g. Whipple and Tucker, 1999). However, the large drainage-area changes triggered by fluvial captures (e.g., García-Castellanos et al., 2003; Willett et al., 2014; Viaplana-Muzas et al., 2018) impose even more abrupt changes in the incision and transport regime. An extreme scenario for such large drainage area changes is provided by the opening of internally-drained basins (endorheic-exorheic transition) by either river piracy, basin overfilling, basin overspilling or a combination of the three mechanisms (García-Castellanos, 2006).

Internally-drained (endorheic) basins have no connection to the oceanic base level, evaporating the collected water runoff in central lake systems. Accordingly, they trap all the products of erosion within the continent, often leading to the formation of low relief at high elevation (Han et al., 2019). Their transit to exorheism is driven by either capture, overfilling, or overspilling through a basin outlet located at a topographic sill, the lowest point along the basin’s drainage divide (García-Castellanos, 2006). This causes a sudden lowering of the geomorphological base level and an increase in the stream power near the spillway, where the water so far evaporated in endorheic lakes meets the steep slope at the divide with the river basin that takes the role of basin outlet. As the lake’s outlet is incised, the lake’s level is reduced and, due to its shallow waters, its area becomes rapidly reduced. This is independent on whether the opening is dominated by basin overfilling or by fluvial capture (piracy). The spillway and the just extended drainage basin then tend to adjust to the new base level (Howard et al., 1994) resulting in an upstream propagation of an erosion wave (e.g., Loget and Van Den Driessche, 2009; Struth et al., 2019; Rodríguez-Rodríguez et al., 2020a, 2020b). This retrogressive incision wave erodes the elevated low relief formed during the endorheic conditions. The time-scale of the landscape response to these perturbations is poorly understood (e.g. Whipple, 2001; Whittaker and Boulton, 2012). Previous works analyzed the effect of base level changes (e.g. Bishop, 1995; Burbank and Anderson, 2001), erodibility (e.g. DiBiase and Whipple, 2011), climate (e.g. Schumm, 1979), or elastic thickness (e.g. Chang and Lijun, 2019) in topographic evolution.

Here, we make use of a numerical model integrating all these mechanisms applied to a scenario inspired by the endorheic-exorheic transition of the Duero (Douro in Portuguese) and Ebro drainage basins, covering most of northern Iberia. Remarkably, these two basins have a similar catchment area (∼98,000 km2 and ∼86,000 km2, respectively) and geological history but a completely different present topography. During the Late Miocene, both the Duero and Ebro basins became overfilled with sediment, initiating their incision by a new drainage network. This is based on the mammal fauna (Santisteban et al., 1997) and on magnetostratigraphic dating of the top of the Páramo Formation of the northeastern most Duero Basin (9.7–9.6 Ma according to Krijgsman et al., 1996; ∼9.1 Ma according to Beamud et al., 2006). The same technique combined with flexural isostatic modeling at the Ebro Basin yields an age of opening of 7.5 to 12 Ma (García-Castellanos and Larrasoaña, 2015; García-Castellanos et al., 2003). However, the onset of drainage entrenchment through the Iberian Massif and the Catalan Coastal Range (CCR) barriers (Fig. 1A) seems much slower or delayed by several million years (Cunha et al., 2019a) and their present-day topography and drainage networks are very distinct, with a more eroded, concave-up, and equilibrated Ebro River profile relative to the convex-up Duero River profile (Fig. 1B). The paradox presented by the Duero Basin is that the end of lacustrine sedimentation (9.6 Ma; Krijgsman et al., 1996) in the central areas largely precedes the estimated onset of fluvial incision of the basin (3.7 to 1 Ma; e.g., Silva et al., 2017; Cunha et al., 2019a; Rodríguez-Rodríguez et al., 2020b) and the last minor deposits in the margins of the sedimentary basin (Quaternary alluvial fans and depositional glacis or rañas; e.g., de Vicente et al., 2018; Cunha et al., 2019b). These latest sedimentary infill is considered to be roughly coeval with the onset of basin incision and terrace formation. Recently, however, based on terrace cosmogenic dating, Rodríguez-Rodríguez et al. (2020b) have argued for an earlier onset of basin incision. The time gap between lacustrine sedimentation and the onset of exorheism and the excavation of the sedimentary basin is intriguing because, by definition, internal drainage is a trap for the erosion products of the drained areas and therefore lake sedimentation cannot end sooner than exorheism.

To investigate this seeming conflict, the nature and timing of this transient landscape and the key mechanisms controlling the basin drainage evolution, we use numerical modelling of isostasy and river incision constrained with published geological data, knickpoint analysis, and erosion and sedimentation rates and volumes (e.g. Vacherat et al., 2018; Struth et al., 2019; Cunha et al., 2019a) (Tab. 1). We evaluate specifically the effects of geographical (barrier width, basin elevation), lithospheric (elastic thickness; Te) and lithological (barrier erodibility) parameters on the evolution of the topography under the situation of a high-altitude spill point.

thumbnail Fig. 1

(A) Current topography of the Ebro and Duero basins in the northern part of Iberia. Black dashed line represents the current position of the bedrock incision knickpoint wave in the Duero Basin (Struth et al., 2019). (B) Swath profile of both basins indicating the maximum (red line), mean (green line) and minimum elevation (blue line; river profile). Zmax indicates the maximum infill elevation. CCR: Catalan Coastal Range. The x-axis indicates the age of the bedrock and sediment outcropping along the profile.

Table 1

Input parameters for the reference setup of the TISC surface-processes model.

2 Modelling methodology

We model the topographic evolution using the code TISC, which calculates lithospheric flexural isostasy (elastic thin plate) in response to surface mass transport. This transport is in turn calculated via a fluvial erosion-limited incision and stream-power sediment transport (García-Castellanos, 2002; García-Castellanos and Jiménez-Munt, 2015). The model therefore disregards a significant role of other sources of relative topographic changes such as eustasy of dynamic topography, but isostasy seems to capture the most of the vertical motions in the region (García-Castellanos and Larrasoaña, 2015). The initial model setup is an idealization of the Duero and Ebro basin geology (Figs. 13A). The input parameters are summarized in Table 1. We adopt a precipitation rate (actually the runoff) linearly increasing with elevation to roughly capture the present-day runoff distribution in the two drainage basins, neglecting past changes in precipitation. The values adopted (Tab. 1) correspond to the average runoff in the region and a conservative evaporation rate that allows for the lakes to cover an area similar to the lacustrine sediment outcrops. In TISC, water flows from cell to cell in a rectangular grid following the steepest descent among eight neighboring cells, forming lakes that fill topographic minima and finding their outlet (if after evaporation the lake keeps extra water, therefore exorheic) or the extent of the endorheic lake (the lake’s area is reduced until finding a water balance between its inflows and the lake surface evaporation). The erosion model follows a hybrid detachment- and transport-limited approach in which incision rates are a power law of the basal shear stress of the river flow: dzdt=ke(ττc)a(τ>τc),(1) where z is the elevation, t is time, ke is the erodability (García-Castellanos and O’Connor, 2018), τc is the shear stress at the base of the water flow, and τc is the critical shear stress. For simplicity, here we adopt τc = 0 and a= 1.5. This leads to a stream power law for erosion with the form (see García-Castellanos and Jiménez-Munt, 2015): dzdt=K, Qm,Sn,(2)with: K=ke(pg)a(nkw)sas,m=3a8,n=13a16,(3)where ρ is the density of water, g the gravitational acceleration, n the roughness coefficient, kw the proportionality of channel width to discharge. The sediment load q acquired by the river (i.e., the erosion rate in Eq. (2)) is limited by a transport capacity qeq, adopted proportional to stream power: qeq=KcapQS.(4)

Specifically, TISC multiplies the erosion rate in Equation (2) by (qeq-q)/qeq in order to cancel erosion and trigger aggradation when the river exceeds its transport capacity. When the river transports a sediment load q larger than qeq then dz/dt becomes positive, implying fluvial aggradation in areas of low channel slope, or delta sedimentation in lakes or in the ocean, where Kcap is assumed null.

As an innovation relative to the modeling methodology in García-Castellanos and Jiménez-Munt (2015), we have incorporated here the flow of underground water along the sedimentary units. For this, we solve the Darcy’s law in finite difference assuming an isotropic porous medium and using a permeability of 10−14 m2 for the sedimentary cover. Remarkably, underground flow multiplies by a factor 2 to 5 the migration velocity of the drainage divide between both sedimentary basins, allowing to reproduce the migration of tens of kilometers observed at the Ebro-Duero divide (Mikeš, 2009; Vacherat et al., 2018).

thumbnail Fig. 2

Four stages of the evolution of the reference model. Each stage shows the calculated elevation and drainage network (left) and the incision rate distribution (right). The horizontal red lines locate the sections in Figure 3. A video animation of this model can be accessed at https://youtu.be/nF0obRQWL60.

thumbnail Fig. 3

Reference model results. Panels (A) and (B) show the initial setup at t = −10 Myr. (C) Topography obtained at 0 Myr (present). (D) cross-section showing isostatic subsidence and rebound (dashed line and arrows). (E) Cumulative uplift (dashed) and subsidence (bold lines) and sediment thickness. (F) Sedimentation (−)/erosion (+) rates at t = 0. River profiles with knickpoint location for the Duero (G) and Ebro (H) river networks. The points displayed in these boxes indicate the location of bedrock (yellow, B_KP) and cover (red, C_KP) knickpoints (compare with Struth et al., 2019). Arrows indicate the distance traveled by the knickpoints from the initial location (yellow circle at x = 50 km).

3 Model setup

The model starts at t = −10 Myr with both basins set very close to overspilling (based on Krijgsman et al., 1996, for the Duero and García-Castellanos et al., 2003, for the Ebro), and ends at 0 Myr (present). We choose our preferred set of parameters (hereafter named reference model) by fitting the knickpoint wave position and retreat rate (Struth et al., 2019; Rodríguez-Rodríguez et al., 2020b), delta sediment volume and basin erosion rates (see Tab. 2). We set the preserved sedimentary infill elevation in the reference model at 800 m above sea level, a compromise between the 535–750 m derived for the Ebro Basin (García-Castellanos and Larrasoaña, 2015) and the 750–850 m at the Duero Basin (Antón et al., 2012, Antón et al., 2014).

This reference scenario adopts an equivalent elastic thickness for the lithosphere of Te= 15 km (representative for the values obtained by Ruiz et al., 2006, Kaban et al., 2018; Gaspar-Escribano et al., 2004), a width of 175 and 25 km for the Duero and Ebro barriers, respectively, with an erodibility of 2.6 10−8 m yr−1 Pa−1.5 for both bedrock barriers and 4.0e−7 m yr−1 Pa−1.5 for the sediment infill. The former erodibility values, needed to fit the measured erosion rates in the Duero Arribes outlet area, are the lowest ever derived from long-term river erosion, to our knowledge, consistent with the solid, unaltered granite lithology forming most of the highest reaches of the outlet of the Duero sedimentary basin (García-Castellanos and O’Connor, 2018). We calculate mean erosion and sedimentation rates for the basin and the delta area. Starting from this reference model, we test the effect of varying the lithospheric elastic thickness (Te), the barrier width and the barrier erodibility.

Table 2

Results from the reference model and comparison with published observations and results.

4 Modeling results

The results from the reference setup (Figs. 2 and 3; see animation in https://youtu.be/nF0obRQWL60) show that, as the basin’s drainage opens finding an outlet towards the ocean, incision gradually propagates into the continent first causing the extinction of the lakes, then the reintegration of the drainage, and then the incision of the sedimentary basins. For large values of elastic thickness Te, the lithosphere behaves too rigidly to induce distinct isostatic vertical motions in response to erosion (Fig. 4A). However, for values lower than 27.5 km (encompassing most Te estimations in Iberia) the isostatic rebound in response to erosion of the barrier and the subsidence related to offshore sedimentation significant cause vertical motions (Fig. 3E). These, in the case of the wider barrier of the western (Duero) basin, further delays in the integration and incision of the lake in the new drainage network. Wider barriers (Fig. 4B) and higher erodibility (Fig. 4C) both extend the time needed to fully capture and extinguish the lake and to entrench and incise the new drainage network. The earliest penetration of the incision wave into the sedimentary basin is thus obtained for higher effective elastic thickness, lower barrier width and higher erodibility.

In spite of the similar initial setting and overspill adopted, the times of maximum basin erosion and delta sedimentation rates (Fig. 5) differ by 7 Myr between the left and right basins of the reference model. The thin layer of sediment prescribed on the Duero barrier and the corresponding spillway allows a fast capture of the basin (a rapid end of endorheism) while the isostatic rebound imposes a slow incision propagation through the less erodible barrier bedrock. Narrow topographic barriers and smaller Te values imply that erosion and sedimentation in the margins induce more isostatic vertical uplift at the edge of the continental basin. The large sediment erosion in the Ebro Basin (unloading of ca. 40,000 km3, in agreement with estimations in García-Castellanos et al., 2003) results in a similar sediment accumulation in the delta (22% of the sediment is lost across the boundaries of the model) and its corresponding flexural subsidence (Fig. 3). This surface mass transport causes the isostatic rebound of the basin area by up to 600 m (Fig. 3D). This rebound is enhanced by the lesser contribution of the forebulge uplift in response to the sediment loading in the delta (Fig. 3D, E). Interestingly, this mechanism provides an explanation for the exhumation of older sediments in the eastern part of the Ebro basin (see the age of the outcropping sedimentary infill in Fig. 1B; also Lewis et al., 2000). Wider barriers inhibit the forebulge isostatic effect by separating the sediment source from the sink, as apparent in the Duero side in Figures 3D, E and 6.

The barrier width and erodibility also determine the erosion pattern and headwards propagation. On one hand, a wider and less erodible barriers delay the erosive wave retreat towards the basin. The difference in barrier width (175 km for the Duero barrier against 25 km for the Ebro barrier), results in a time lag between the erosion peak at both basins of about 7 Myr (see Fig. 4A). After reaching that peak, erosion rate decreases towards the steady-state in both basins. The results show with two different knickpoint retreat velocities for the bedrock and the cover at 0.02 to 0.03 m yr−1 and 0.10 to 0.20 m yr−1, respectively, close to what is derived from knickpoint analysis by Struth et al. (2019) and cosmogenic dating by Rodríguez-Rodríguez et al. (2020b). This has been tuned by finding the best-fitting ratio of basin-sediment erodability relative to bedrock (a factor 20; Tab. 1). The time needed for the faster cover-knickpoint wave to reach the central divide (x = 0) is 1.4 Myr in the case of the Ebro and between 7.4–6.0 Myr for the Duero. The sooner arrival of the Ebro Basin is due to the narrower topographic barrier, which also explains that both the cover and bedrock knickpoints propagate undistinguishable at the same speed.

Importantly, the adopted erodability values are also compatible with the average incision rates obtained for both basins, ensuring that the surface process model is capturing the first order response of both vertical incision and headwards erosion retreat. What ensures that this reference model is eroding at reasonable rates is the present-day (t = 0) average erosion rates of 21.2 m/Myr (Duero) and 70.1 m/Myr (Ebro), close to the 26 m/Myr measured by Schaller et al. (2016) in the Esla R. tributary of the Duero R. In contrast, the vertical incision along the main river in the reference model (Fig. 3F) is up to 144 m/Myr, in agreement with the 122–250 m/Ma estimated from cosmogenic dating (Rodríguez-Rodríguez et al., 2020b) for the central Duero River. As for the Ebro Basin, the erosion rates in Figure 3F satisfactorily reproduce the 82–127 m/Myr range compiled by Regard et al. (2021). These satisfactory results support the appropriateness of the analytical expression in Equation (3) (details in García-Castellanos and Jiménez-Munt, 2015), and particularly the exponents m = 0.56 and n = 1.22 derived for a = 1.5.

The 4–6 Myr delay (depending mainly on the erodibility of bedrock) between the onset of exorheism and the widespread incision of the basin predicted for the basin on the left side (Fig. 2) is a consequence of the wider hard-rock barrier separating this basin from the ocean. This mechanism also leads to the coeval incision of the center of the basin while small sedimentation rates are still occurring at the basin margins (best visible in the video). This suggests an explanation for the deposition of the Quaternary raña glacis (e.g. Alonso-Zarza et al., 2002) long after the end of bulk sedimentation in most of the Duero Basin.

thumbnail Fig. 4

Model sensitivity test. Age of lake extinction and basin incision onset for varying elastic thickness (A), barrier width (B), and erodibility (C). The grey band indicates the value used for the reference model. t = 0 Myr indicates present day.

thumbnail Fig. 5

Model results for the Ebro (black lines) and Duero (red lines) basins from −10 to +20 Myr. (A) Mean erosion rate Er of the two sedimentary basins, showing a trend currently decreasing in the Ebro and increasing in the Duero. (B) Mean sedimentation rates Sr in the delta areas. The inset shows the onset of sedimentation at the Duero delta and synchronous initiation of the incision of the basin. (C) Volume of sediment in the deltas, showing a current stabilization in the Ebro delta contrasting with the increasing volume of the Duero delta (diamonds indicate published data; Table DR-2).

thumbnail Fig. 6

Distribution of the cumulative sedimentation and erosion (black line) and the isostatic rebound (red dashed line) through the transect displayed in red in Figure 3.

5 Discussion on the timing of transient landscape evolution

Internally-drained basins form often amid uplifting orogens that create topographic lows while blocking them from the income of humid air (e.g. Carroll et al., 2010; Han et al., 2019). The dry conditions allow cancelling the outlet of these basins and prolonging their life (García-Castellanos, 2006). Our modeling results suggest that basins without tectonically active boundaries (or long after the ending of tectonic deformation) can also have their opening delayed by low values of lithosphere rigidity, implying a more local flexural response to the unloading of the basin divides (e.g. Heller et al., 2011). Eroding the boundaries of a closed basin in a low lithospheric rigidity region (normally associated with a young tectono-thermal age of the lithosphere) leads to flexural uplift of the basin flanks, and then to a delay in erosion and reintegration of the internal basin. This mechanism was first reported in the Ebro Basin (García-Castellanos et al., 2003) and the Great Divide Basin (Wyoming, USA) by Heller et al. (2011) and our results show that it is also consistent with the evolution of the Duero Basin and partially explains its different evolution relative to the neighboring Ebro Basin.

The timescales of transient landscape evolution are hard to determine because conditions of paleoclimate and tectonics are often not quantified, and because the identification of a steady-state topography in the field is highly speculative (Willett and Brandon, 2002; Goren, 2016). According to Beeson et al. (2017), the time required to reach the steady-state after a tectonic or climatic perturbation is related to the timescale of river steepness adjustment at individual basins, to balance erosion and rock uplift, defining a transient landscape time. Previous studies (e.g. Paola et al., 1992; Beaumont et al., 2000; Whipple, 2001) suggest that each geomorphic system has an intrinsic response time (or equilibrium time) to approach a balance between the relief-generating mechanisms (e.g., tectonics or isostasy) and the gradual levelling of the surface carried out by surface transport processes. This relates to the time-scales needed after a climatic perturbation, a change in drainage connectivity, or a change in tectonic deformation rates for the average elevation to asymptotically reach stability (Castelltort and Van Den Driessche, 2003) and for sedimentation/erosion rates to remain constant. This response time is often estimated in the order of millions of years (Pazzaglia, 2003, Whipple, 2001, Whipple and Tucker, 1999; García-Castellanos et al., 2003). Based on an analytical formulation, Whipple (2001) estimated a range of 0.25 to 2.5 Myr to adjust to a change in rock uplift rate or base level fall, depending on the non-linearity of the incision rule and the magnitude and type of perturbation. Whittaker and Boulton (2012) calculated a fluvial response time of 1–3 Myr for drainage areas between 10 and 70 km2 in Hatay Graben (Turkey) and the Apennines. Our results indicate that, even under a hypothetically simultaneous drainage opening with a similar base level fall of 800 m at ∼10 Ma for the two basins, the Ebro-like side reaches the stabilization of erosion after ∼9 Myr (time needed to recover 50% of the erosion rate peak; Fig. 5A) while the Duero-like barrier leads to a much slower and delayed network incision (see stabilization in terms of sediment volume in Fig. 5C). This highlights that the different initial topographic and lithological configuration may explain the differential timing of incision. To see the sensitivity to each parameter, and to allow the portability of our results to other scenarios, we also tested the effect of the drainage area and the initial basin elevation. Doubling the drainage area or precipitation rate leads to a 2.0 and 3.2 Myr earlier entrenchment of the left (Duero) basin in the model, respectively, due to the larger amount of water discharge collected by the main river. Figure 7 shows the result of doubling the base level fall from 800 m (initial elevation of the sedimentary basins in the reference model) to 1600 m, obtaining a nearly double value for the erosion/sediment rate and for the isostatic rebound (Fig. 7). Rates of erosion and sedimentation reach a maximum value peak before adapting to the new base level (Figs. 3B and 4A) and then decrease towards a new steady state. Topographic adjustment to the new base level is performed through the propagation of two knickpoint families, a slow one for the bedrock and a fast one along the sedimentary cover. The wider topographic barrier separating the Duero Basin from the ocean implies a time of transitory landscape (time needed to reach a concave-up river profile over the entire basin) of about 30 Myr, 2 to 3 times longer than in the Ebro Basin.

These results support a timescale of landscape transition significantly longer than obtained by the aforementioned studies, up to 30 Myr in our Duero Basin scenario, due to the hard lithology of Variscan granite forming its barrier. In absence of tectonic activity, the transient landscape evolution during the opening of internally-drained basins is mainly controlled by the width and erodibility of the topographic barrier, the size and elevation of the basin above the base level, and the rigidity of the lithosphere. Because erodibility varies by more than eight orders of magnitude (García-Castellanos and O’Connor, 2018), the bedrock lithology along the spillway can exert far more influence on the topographic evolution than the geographical configuration (basin area and elevation). The isostatic vertical motions due to the mass transport during the basin capture are less important, but low values of lithospheric rigidity (typical of regions characterized by young tectonothermal ages such as those affected by the Alpine orogeny) imply a higher amplitude of flexural rebound due to erosion, uplifting the basin margins and delaying the arrival of fluvial incision to the basin, sometimes by millions of years (Fig. 4A), as documented in our study region.

thumbnail Fig. 7

Results from an identical setup as in Figure 3A but for double elevation of the basins (double base level fall: 1600 m). (A) Topography; (B) cross section; (C) isostatic rebound and subsidence in response to surface mass transport (contours) and sediment thickness; sedimentation/erosion rates (D) and river profiles with knickpoint location for the Duero (E) and Ebro (F) river networks. Note the more pronounced erosion, integration, and uplift of both basins, relative to the results for the reference model in Figure 3.

6 Implications for the drainage evolution of North Iberia

The contrasting preservation of the top of the Duero Basin sedimentary infill versus the advanced erosion of the Ebro Basin can be explained by the wider and harder granite barrier of the former, and does not necessarily imply a much older transition to exorheism of the latter. Our modeling results show that the two basins, in spite of their similar area, base level fall, and age of their sedimentary overfilling, have undergone very different transient landscape evolution due to the wider and less erodible barrier that separates the Duero Basin from the Atlantic Ocean. The time lag between the end of the main Cenozoic sedimentation including lacustrine deposits in the Duero Basin (9.6 Ma; Krijgsman et al., 1996) and the onset of basin incision (possibly starting as early as 5 Ma but surely accelerating after 2.5 Ma; Cunha et al., 2019a; Rodríguez-Rodríguez et al., 2020b) can be explained by the slow propagation of the main knickpoint along the wide, resistant granitic barrier of the Duero Basin (Arribes gorge area). As exemplified by the reference model, a hydrological overspill of the Duero towards the Atlantic Ocean around 9.6 Ma may have terminated the shallow, playa-lake, endorheic system while the granite lithology along the outlet could then prevent the basin incision for several million years. Only when the main knickpoint fully traverses the granitic barrier and reaches the Duero sedimentary basin upstream from Miranda do Douro and Zamora does the incision fully propagate into the center of the basin. This implies that, under conditions of resistant lithology of a wide and flat topographic barrier bounding the endorheic basin, exorheism does not immediate lead to basin incision as often assumed (García-Castellanos, 2006; Cunha et al., 2019a), but in contrast a multi-million year delay can be expected. In other words, while the bulk of incision may have started only recently in geological time-scales (3.7 to 1 Ma), this does not exclude that exorheism may have taken place much earlier, explaining the lack of widespread sedimentary infill in the Duero Basin after 9.6 Ma (except for marginal glacis).

The reference model, constrained by measured river incision rates in both basins, also provides a process-based explanation for the two families of knickpoints (bedrock and sediment cover) retreating at rates of ∼0.02 and ∼0.20 m/yr described by Struth et al. (2019) and Rodríguez-Rodríguez et al. (2020b). In this scenario, the marginal depositional piedmont or glacis (also known as rañas) and fluvial fans next to the ranges can be regarded as the result of pediment formation at the foot of the peri-basinal ranges, in areas where the erosional wave had not yet reached during the Pliocene-Quaternary transition. The Duero Basin, in contrast with other dated basin openings/captures such as the Ebro Basin, provides a geological setting where basin overfilling and drainage overspilling and opening ca. 9.6 Ma does not rapidly translate into widespread basin erosion and drainage entrenchment. The Duero basin, due to its wide western topographic barrier made of very hard granitic lithology, underwent a rather delayed acceleration of incision rates only after 2.5 Ma (Rodríguez-Rodríguez, 2019b) or after 1.8 Ma (Cunha et al., 2019b). This scenario can explain the long time elapsed since the end of the main sedimentation in the center of the basin at 9.6 Ma (Krijgsman et al., 1996) until the younger development of incision and terrace formation coeval with pediment or glacis formation in the periphery of the sedimentary basin. It is also consistent with the late transition of the terrace source materials, from recycled basin materials to a bedrock (mountain range) origin, about 1–2 Ma (Rodríguez-Rodríguez et al., 2020b).

Acknowledgments

This work is supported by MITE(CGL2014-59516), SUBTETIS (PIE-CSIC-201830E039) and ALPIMED (PIE-CSIC-201530E082) projects, funded by the Spanish Government. L.R.-R. was funded by the Marie Curie-Clarín COFUND program financed jointly by the 7th WP of the European Union-Marie Curie Actions and Gobierno del Principado de Asturias (Reference ACA-17-19). The work by A.J.-D. was supported by a Juan de la Cierva-Formación postdoctoral contract (ref. FJCI-2016-28878) from the Spanish Ministry of Science, Innovation and Universities. Pedro Cunha, Laurent Husson, and Ciro Cerrone are kindly acknowledged for their detailed and constructive reviews of an earlier version of this manuscript.

References

  • Alonso-Zarza AM, Armenteros I, Braga JC, Munoz A, Pujalte V, Ramos E, et al. 2002. Tertiary. In: Gibbons W, Moreno T, eds. The Geology of Spain. Geological Society, pp. 293–334. [CrossRef] [Google Scholar]
  • Antón L, Rodés A, de Vicente G, Pallàs R, García-Castellanos D, Stuart FM, et al. 2012. Quantification of fluvial incision in the Duero Basin (NW Iberia) from longitudinal profile analysis and terrestrial cosmogenic nuclide concentrations. Geomorphology 165–166: 50–61. https://doi.org/10.1016/j.geomorph.2011.12.036. [CrossRef] [Google Scholar]
  • Antón L, De Vicente G, Muñoz-Martín A, Stokes M. 2014. Using river long profiles and geomorphic indices to evaluate the geomorphological signature of continental scale drainage capture, Duero basin (NW Iberia). Geomorphology 206: 250–261. [CrossRef] [Google Scholar]
  • Antón L, Muñoz-Martín A, de Vicente G, Finnegan NJ. 2018. Deciphering fluvial-capture-indiced erosional patterns at the continental scale on the Iberian Peninsula. In: AGU Fall Meeting Abstracts, pp. EP33A–1905. [Google Scholar]
  • Arche A, Evans G, Clavell E. 2010. Some considerations on the initiation of the present SE Ebro river drainage system: Post- or pre-Messinian? Journal of Iberian Geology 36: 73–85. [Google Scholar]
  • Bartrina MT, Cabrera L, Jurado MJ, Guimerà JA, Roca E. 1992. Evolution of the central Catalan margin of the Valencia Trough (western Mediterranean). Tectonophysics 203: 219–247. [CrossRef] [Google Scholar]
  • Beamud E, Garcés M, Montes M, Nozal F, Calvo JP, López-Olmedio F, et al. 2006. Magnetoestratigrafía del Mioceno de las cuencas del Tajo y del Duero. In: Proceedings of the IV Simposio de Paleomagnetismo Ibérico (MAGIBER). 4. pp. 9–12. [Google Scholar]
  • Beaumont C, Kooi H, Willett S. 2000. Coupled tectonic-surface process models with applications to rifted margins and collisional orogens. In: Summerfield MA, ed. Geomorphology and Global Tectonics. Chirchester: John Wiley and Sons, pp. 28–55. [Google Scholar]
  • Beeson HW, McCoy SW, Keen-Zebert A. 2017. Geometric disequilibrium of river basins produces long-lived transient landscapes. Earth and Planetary Science Letters 475: 34–43. https://doi.org/10.1016/j.epsl.2017.07.010. [CrossRef] [Google Scholar]
  • Bishop P. 1995. Drainage rearrangement by river capture, beheading and diversion. Progress in Physical Geography 19: 449–473. [CrossRef] [Google Scholar]
  • Burbank DW, Anderson RS. 2001. Tectonic Geomorphology. Massachusetts: Blackwell Science, 274 p. [Google Scholar]
  • Carroll AR, Graham SA, Smith ME. 2010. Walled sedimentary basins of China. Basin Research 22: 17–32. [CrossRef] [Google Scholar]
  • Castelltort S, Van Den Driessche J. 2003. How plausible are high-frequency sediment supply-driven cycles in the stratigraphic record? Sedimentary Geology 157: 3–13. https://doi.org/10.1016/S0037-0738(03)00066-6. [CrossRef] [Google Scholar]
  • Chang C, Lijun L. 2019. Distinct responses of intraplate sedimentation to different subsidence mechanisms. Insights from forward landscape evolution simulations. Journal of Geophysical Reasearch: Earth Surface 124: 1139–1159. https://doi.org/10.1029/2018JF004905. [CrossRef] [Google Scholar]
  • Cunha PP, Martins AA, Gomes A, Stokes M, Cabral J, Lopes FC, et al. 2019a. Mechanisms and age estimates of continental scale endorheic to exorheic drainage transition: Douro River, western Iberia. Global and Planetary Change 181: 102985. https://doi.org/10.1016/j.gloplacha.2019.102985. [CrossRef] [Google Scholar]
  • Cunha PP, de Vicente G, Martín-González F. 2019b. Cenozoic sedimentation Along the piedmonts of thrust related basement ranges and Strike-slip deformation belts of the Iberian Variscan Massif. In: Quesada C, Oliveira JT, eds. The Geology of Iberia: a Geodynamic Approach. Regional Geology Review Series 4. Springer Nature Switzerland AG. https://doi.org/10.1007/978-3-030-11190-8_5 (chapter 5, book ID: 473223. ISBN: 978-3-030-11189-2). [Google Scholar]
  • de Vicente G, Cunha PP, Muñoz-Martín A, Cloetingh SAPL, Olaiz A, Vegas R. 2018. The Spanish-Portuguese Central System: an example of intense intraplate deformation and strain partitioning. Tectonics 37: 4444–4469. [CrossRef] [Google Scholar]
  • DiBiase RA, Whipple KX. 2011. The influence of erosion thresholds and runoff variability on the relationships among topography, climate and erosion rate. Journal of Geophysical Reasearch 116: F04036. https://doi.org/10.1028/2011JF002095. [Google Scholar]
  • Evans G, Arche A. 2002. The flux of siliciclastic sediment from the Iberian Peninsula, with particular reference to the Ebro. In: Jones SJ, Frostick LE, eds. Sediment Flux to Basins: Causes, Controls and Consequences. Geological Society, London, Special Publication 191: 199–208. [CrossRef] [Google Scholar]
  • Gaspar-Escribano JM, García-Castellanos D, Roca E, Cloetingh S. 2004. Cenozoic vertical motions of the Catalan Coastal Ranges (NE Spain): the role of tectonics, isostasy, and surface transport. Tectonics 23. https://doi.org/10.1029/2003TC001511 (pdf). [Google Scholar]
  • García-Castellanos D. 2002. Interplay between lithospheric flexure and river transport in foreland basins. Basin Research 14: 89–104. [CrossRef] [Google Scholar]
  • García-Castellanos D. 2006. Long-term evolution of tectonic lakes: climatic controls on the development of internally drained basins. In: Willett SD, Hovius N, Brandon MT, Fisher DM, eds. Tectonics, climate, and landscape evolution. GSA Special Paper 398: 283–294. https://doi.org/10.1130/2006.2398(17). [Google Scholar]
  • García-Castellanos D, Jiménez-Munt I. 2015. Topographic evolution and climate aridification during continental collision: insights from computer simulations. PLoS One. https://doi.org/10.1371/journal.pone.0132252. [Google Scholar]
  • García-Castellanos D, Larrasoaña JC. 2015. Quantifying the post-tectonic topographic evolution of closed basins: the Ebro basin (northeast Iberia). Geology 43: 663–666. https://doi.org/10.1130/G36673.1. [CrossRef] [Google Scholar]
  • García-Castellanos D, Vergés J, Gaspar-Escribano J, Cloetingh S. 2003. Interplay between tectonics, climate, and fluvial transport during the Cenozoic evolution of the Ebro Basin (NE Iberia). Journal of Geophysical Research: Solid Earth 108(B7): 2347. https://doi.org/10.1029/2002JB002073. [Google Scholar]
  • García-Castellanos D, O’Connor J. 2018. Outburst floods provide erodability estimates consistent with long-term landscape evolution. Scientific Reports (Nature Pub) 8: 10573. https://doi.org/10.1038/s41598-018-28981-y. [CrossRef] [Google Scholar]
  • Goren L. 2016. A theoretical model for fluvial channel response time during time-dependent climatic and tectonic forcing and its inverse applications. Geophysical Research Letters 43(20): 10,753–10,763. https://doi.org/10.1002/2016GL070451. [CrossRef] [Google Scholar]
  • Han Z, Sinclair HD, Li Y, Wang C, Tao Z, Qian X, et al. 2019. Internal drainage has sustained low-relief Tibetan landscapes since the early Miocene. Geophysical Research Letters 46: 8741–8752. https://doi.org/10.1029/2019GL083019. [CrossRef] [Google Scholar]
  • Heller PL, McMillan ME, Humphrey N. 2011. Climate-induced formation of a closed basin: Grat Divide Basin, Wyoming. Geological Society of America Bulletin 123(1–2): 150–157. https://doi.org/10.1130/B30113.1. [CrossRef] [Google Scholar]
  • Howard AD, Dietrich WE, Seidl MA. 1994. Modeling fluvial erosion on regional to continental scales. Journal of Geophysical Research: Solid Earth 99(B7): 13971–13986. [CrossRef] [Google Scholar]
  • Kaban MK, Chen B, Tesauro M, Petrunin AG, El Khrepy S, Al-Arifi N. 2018. Reconsidering effective elastic thickness estimates by incorporating the effect of sediments: a case study for Europe. Geophysical Research Letters 45(18): 9523–9532. https://doi.org/10.1029/2018GL079732. [CrossRef] [Google Scholar]
  • Krijgsman W, Garcés M, Langereis CG, Daams R, van Dam J, van der Meulen AJ, et al. 1996. A new chronology for the middle to late Miocene continental record in Spain. Earth and Planetary Science Letters 142: 367–380. [CrossRef] [Google Scholar]
  • Lewis CJ, Vergés J, Marzo M. 2000. High mountains in a zone of extended crust: insights into the Neogene-Quaternary topographic development of northeastern Iberia. Tectonics 19: 86–102. [CrossRef] [Google Scholar]
  • Loget N, Van Den Driessche J. 2009. Wave train model for knickpoint migration. Geomorphology 106: 376–382. [CrossRef] [Google Scholar]
  • Maillard A, Mauffret A, Watts AB, Torné M, Pascal G, Buhl P, et al. 1992. Tertiary sedimentary history and structure of the Valencia Trough, (Western Mediterranean). Tectonophysics 204: 57–75. [CrossRef] [Google Scholar]
  • Mikeš D. 2009. Forward modelling of the Bureba escarpment (Rioja Trough, N-Spain). Central European Journal of Geosciences 1: 463–471. https://doi.org/10.2478/v10085-009-0037-0. [Google Scholar]
  • Nelson CH, Maldonado A. 1990. Factors controlling late Cenozoic continental margin growth from the Ebro Delta to the western Mediterranean Deep Sea. Marine Geology 95: 419–440. [CrossRef] [Google Scholar]
  • Paola C, Heller PL, Angevine CL. 1992. The large-scale dynamics of grain size variation in alluvial basins. 1, Theory. Basin Research 4: 73–90. [CrossRef] [Google Scholar]
  • Pazzaglia FK. 2003. Landscape evolution models. Developments in Quaternary Science 1: 247–274. https://doi.org/10.1016/S1571-0866(03)01012-1. [CrossRef] [Google Scholar]
  • Regard V, Vacherat A, Bonnet S, Mouthereau F, Nørgaard J, Knudsen MF. 2021. Late Pliocene-Pleistocene incision in the Ebro Basin (North Spain). Bulletin de la Société Géologique de France 192: 30. https://doi.org/10.1051/bsgf/2021020. [CrossRef] [EDP Sciences] [Google Scholar]
  • Rodríguez-Rodríguez L, Antón L, Pallàs R, García-Castellanos D, Jiménez-Munt I, Pastor-Martín C. 2020a. A GIS method to identify flat surfaces and restore relict fluvial long-profiles from terrace remnants gives new clues on how large basins respond to endorheic-exorheic transitions (Duero basin, sIberian Peninsula). Earth Surface Processes and Landforms 45(4): 1013–1027. https://doi.org/10.1002/esp.4795. [CrossRef] [Google Scholar]
  • Rodríguez-Rodríguez L, Antón L, Rodes Á, Pallàs R, García-Castellanos D, Jiménez-Munt I, et al. 2020b. Dates and rates of endo-exorheic drainage development: insights from fluvial terraces (Duero River, Iberian Peninsula). Global and Planetary Change 193: 103271. https://doi.org/10.1016/j.gloplacha.2020.103271. [CrossRef] [Google Scholar]
  • Ruiz J, Gomez-Ortiz D, Tejero R. 2006. Effective elastic thicknesses of the lithosphere in the Central Iberian Peninsula from heat flow: Implications for the rheology of the continental lithospheric mantle. Journal of Geodynamics 41(5): 500–509. ISSN: 0264-3707. https://doi.org/10.1016/j.jog.2006.01.005. [CrossRef] [Google Scholar]
  • Santisteban JI, Alcalá L, Mediavilla RM, Alberdi MT, Luque L, Mazo A, et al. 1997. El yacimiento de Tariego de Cerrato: El inicio de la red fluvial actual en el sector central de la Cuenca del Duero. Cuadernos de geología ibérica (22): 431–446. ISSN 0378-102X. [Google Scholar]
  • Schumm SA. 1979. Geomorphic thresholds: the concept and its applications. Transactions of the Institute of British Geographers 4: 485–515. [CrossRef] [Google Scholar]
  • Schaller M, Ehlers TA, Stor T, Torrent J, Lobato L, Christl M, Vockenhuber C. 2016. Timing of European fluvial terrace formation and incision rates constrained by cosmogenic nuclide dating. Earth and Planetary Science Letters 451: 221–231. https://doi.org/10.1016/j.epsl.2016.0.022. [CrossRef] [Google Scholar]
  • Silva PG, Roquero E, López-Recio M, Huerta P, Martínez-Graña AM. 2016. Chronology of fluvial terrace sequences for large Atlantic rivers in the Iberian Peninsula (Upper Tagus and Duero drainage basins, Central Spain). Quaternary Science Reviews 166: 188–203. https://doi.org/10.1016/j.quascirev.2016.05.027. [Google Scholar]
  • Silva PG, Roquero E, López-Recio M, Huerta P, Martínez-Graña AM. 2017. Chronology of fluvial terrace sequences for large Atlantic rivers in the Iberian Peninsula (Upper Tagus and Duero drainage basins, Central Spain). Quaternary Science Reviews 166: 188–203. [CrossRef] [Google Scholar]
  • Struth L, García-Castellanos D, Viaplana-Muzas M, Vergés V. 2019. Drainage network dynamics and knickpoint evolution in the Ebro and Duero basins: from endorheism to exorheism. Geomorphology 327: 554–571. https://doi.org/10.1016/j.geomorph.2018.11.033. [CrossRef] [Google Scholar]
  • Vacherat A, Bonnet S, Mouthereau F. 2018. Drainage reorganization and divide migration induced by the excavation of the Ebro basin (NE Spain). Earth Surface Dynamics 6: 369–387. https://doi.org/10.5194/esurf-6-369-2018. [CrossRef] [Google Scholar]
  • Viaplana-Muzas M, Babault J, Dominguez S, Van Den Driessche J, Legrand X. 2018. Modeling of drainage dynamics influence on sediment routing system in a fold-and-thrust belt. Basin Research 31(2): 290–310. https://doi.org/1.1111/bre.12321. [Google Scholar]
  • Whipple KX. 2001. Fluvial landscape response time: how plausible is steady-state denudation? American Journal of Science 301 (4–5): 313–325. [CrossRef] [Google Scholar]
  • Whipple KX, Tucker GE. 1999. Dynamics of the stream power river incision model: implications for height limits of mountain ranges, landscape response time scales and research needs. Journal of Geophysical Research 104: 17661–17674. https://doi.org/10.1029/1999JB900120. [CrossRef] [Google Scholar]
  • Whittaker AC, Boulton SJ. 2012. Tectonic and climatic controls on knickpoint retreat rates and landscape response times. Journal of Geophysical Research 117: F02024. https://doi.org/10.1029/2011JF002157. [Google Scholar]
  • Willett SD, Brandon MT. 2002. On steady states in mountain belts. Geology 30: 175–178. [CrossRef] [Google Scholar]
  • Willett SD, McCoy SW, Perron JT, Goren L, Chen CY. 2014. Dynamic re-organization of river basins. Science 343: 1248765. https://doi.org/10.1126/science.1248765. [CrossRef] [Google Scholar]

Cite this article as: Struth L, García-Castellanos D, Rodríguez-Rodríguez L, Viaplana-Muzas M, Vergés J, Jiménez-Díaz A. 2021. Topographic, lithospheric and lithologic controls on the transient landscape evolution after the opening of internally-drained basins. Modelling the North Iberian Neogene drainage, BSGF - Earth Sciences Bulletin 192: 45.

All Tables

Table 1

Input parameters for the reference setup of the TISC surface-processes model.

Table 2

Results from the reference model and comparison with published observations and results.

All Figures

thumbnail Fig. 1

(A) Current topography of the Ebro and Duero basins in the northern part of Iberia. Black dashed line represents the current position of the bedrock incision knickpoint wave in the Duero Basin (Struth et al., 2019). (B) Swath profile of both basins indicating the maximum (red line), mean (green line) and minimum elevation (blue line; river profile). Zmax indicates the maximum infill elevation. CCR: Catalan Coastal Range. The x-axis indicates the age of the bedrock and sediment outcropping along the profile.

In the text
thumbnail Fig. 2

Four stages of the evolution of the reference model. Each stage shows the calculated elevation and drainage network (left) and the incision rate distribution (right). The horizontal red lines locate the sections in Figure 3. A video animation of this model can be accessed at https://youtu.be/nF0obRQWL60.

In the text
thumbnail Fig. 3

Reference model results. Panels (A) and (B) show the initial setup at t = −10 Myr. (C) Topography obtained at 0 Myr (present). (D) cross-section showing isostatic subsidence and rebound (dashed line and arrows). (E) Cumulative uplift (dashed) and subsidence (bold lines) and sediment thickness. (F) Sedimentation (−)/erosion (+) rates at t = 0. River profiles with knickpoint location for the Duero (G) and Ebro (H) river networks. The points displayed in these boxes indicate the location of bedrock (yellow, B_KP) and cover (red, C_KP) knickpoints (compare with Struth et al., 2019). Arrows indicate the distance traveled by the knickpoints from the initial location (yellow circle at x = 50 km).

In the text
thumbnail Fig. 4

Model sensitivity test. Age of lake extinction and basin incision onset for varying elastic thickness (A), barrier width (B), and erodibility (C). The grey band indicates the value used for the reference model. t = 0 Myr indicates present day.

In the text
thumbnail Fig. 5

Model results for the Ebro (black lines) and Duero (red lines) basins from −10 to +20 Myr. (A) Mean erosion rate Er of the two sedimentary basins, showing a trend currently decreasing in the Ebro and increasing in the Duero. (B) Mean sedimentation rates Sr in the delta areas. The inset shows the onset of sedimentation at the Duero delta and synchronous initiation of the incision of the basin. (C) Volume of sediment in the deltas, showing a current stabilization in the Ebro delta contrasting with the increasing volume of the Duero delta (diamonds indicate published data; Table DR-2).

In the text
thumbnail Fig. 6

Distribution of the cumulative sedimentation and erosion (black line) and the isostatic rebound (red dashed line) through the transect displayed in red in Figure 3.

In the text
thumbnail Fig. 7

Results from an identical setup as in Figure 3A but for double elevation of the basins (double base level fall: 1600 m). (A) Topography; (B) cross section; (C) isostatic rebound and subsidence in response to surface mass transport (contours) and sediment thickness; sedimentation/erosion rates (D) and river profiles with knickpoint location for the Duero (E) and Ebro (F) river networks. Note the more pronounced erosion, integration, and uplift of both basins, relative to the results for the reference model in Figure 3.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.