Introduction
Endorheic basins, also known as closed or internally drained basins, are subject to intensive research, because they provide important insights into the evolution of sedimentary basins, climate changes, or tectonic processes (Mather, Reference Mather2000; García Castellanos et al., Reference García-Castellanos, Vergés, Gaspar-Escribano and Cloetingh2003; Sobel et al., Reference Sobel, Hilley and Strecker2003; García Castellanos and Cruz Larrasoaña, Reference García-Castellanos and Cruz Larrasoaña2015; Heidarzadeh et al., Reference Heidarzadeh, Ballato, Hassanzadeh, Ghassemi and Strecker2017; Bridgland et al., Reference Bridgland, Westaway and Hu2020; among many others). There are only a few examples of outcropping Quaternary continental records of endorheic basins (e.g., Silva et al., Reference Silva, Roquero, López-Recio, Huerta and Martínez-Graña2017; Stokes et al., Reference Stokes, Mather, Rodes, Kearsey and Lewin2018; Bridgland et al., Reference Bridgland, Westaway and Hu2020). This lack of outcropping endorheic basins occurs because, apart from drilling and geophysical surveys, such information is revealed only once the sedimentary record of the basin has been exposed, which implies the dissection of the basins by river networks during a subsequent exorheic phase. Therefore, the capture of an endorheic basin is a major milestone in its landscape evolution (e.g., Merritts et al., Reference Merritts, Vincent and Wohl1994; García Castellanos et al., Reference García-Castellanos, Vergés, Gaspar-Escribano and Cloetingh2003; Arboleya et al., Reference Arboleya, Babaut, Owen, Teixell and Finkel2008; Struth et al., Reference Struth, Garcia-Castellanos, Viaplana-Muzas and Vergés2019). Consequently, the detailed characterization of these basins is crucial to fully understand the contribution and interplay of tectonic and climate processes in their sedimentary filling, as well as the erosion dynamics during their later exorheic stages.
The Guadix-Baza Basin (GBB) is an example of a Plio-Quaternary endorheic basin that was captured and subsequently dissected by river incision and slope erosion during the Quaternary period (Fig. 1). Extensive sedimentation during the endorheic stage produced a continuous continental sedimentary succession that was several hundred meters thick. This succession represents the most well-preserved and accessible continental Plio-Quaternary stratigraphic record in Europe. This record includes more than 150 paleontological sites of vertebrates (Agustí, Reference Agustí1986; Alberdi and Bonadonna, Reference Alberdi and Bonadonna1989; Agustí et al., Reference Agustí, Lozano-Fernández, Oms, Piñero, Furió, Blain, López-García and Martínez-Navarro2015; Maldonado-Garrido et al., Reference Maldonado-Garrido, Piñero and Agustí2017; among many others). Moreover, the paleontological record of the GBB also comprises the oldest human fossils of western Europe (Toro-Moyano et al., Reference Toro-Moyano, Martínez-Navarro, Agustí, Souday, Bermúdez de Castro, Martinón-Torres and Fajardo2013) and some of the oldest stone tool industries in Europe (also older than 1 Ma; Toro-Moyano et al., Reference Toro-Moyano, Barska, Cauche, Celiberti, Grégoire, Lebegue, Moncel, Lumley and de2011). After the transition from endorheic to exorheic conditions (Calvache and Viseras, Reference Calvache and Viseras1997), extensive deposition in the basin ended and was replaced by intense incisions of the fluvial network into the flat top surface of the basin (glacis surface; sensu García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b). The new base level (Atlantic Ocean) triggered an intense fluvial network incision in the basin that reached a maximum of 300 m below the glacis surface (Fig. 2). However, although the capture of the GBB is a major event that configured most of the critical features of the basin, the age of this process remains elusive.
Geologic map of the Betic Cordillera showing the location of the Guadix-Baza Basin (GBB). The Guadiana Menor River (the main river of the GBB) and the Guadalquivir River are also depicted.
Oblique panoramic views of the badlands landscape of the Guadix-Baza Basin (GBB). The flat elevated surface is the glacis. (a) Western sector; (b) eastern sector.
In this work, we provide the oldest numerical ages obtained thus far for exorheic deposits in the GBB. We integrate these results with previously reported data for the youngest endorheic deposits to discuss the age of the capture. Our results not only provide crucial information for the geologic history of the GBB but also permit a more general approach that can be applied to other internally drained basins and that will help in understanding how these basins were captured and subsequently reshaped. Furthermore, our conclusions enhance the understanding of both geologic and geomorphic settings of some of the most relevant archaeological and paleontological sites in Europe.
Geologic Setting
The GBB is an intramontane basin more than 4000 km2 in size located in the Central Betic Cordillera (SE Spain; Figs. 1 and 3). This region undergoes active tectonic deformation related to the 5 mm/yr, NNW‒SSE convergence between the Nubian and Eurasian plates (DeMets et al., Reference DeMets, Gordon, Argus and Stein1994; Serpelloni et al., Reference Serpelloni, Vannucci, Pondrelli, Argnani, Casula, Anzidei, Baldi and Gasperini2007; Nocquet, Reference Nocquet2012). As a result, the Central Betic Cordillera experiences NNW‒SSE shortening (Sanz de Galdeano, Reference Sanz de Galdeano1983) and ENE‒WSW orthogonal extension (Galindo-Zaldívar et al., Reference Galindo-Zaldívar, Gil, Sanz de Galdeano, Lacy, García-Armenteros, Ruano, Ruiz, Martínez-Martos and Alfaro2015; Martin Rojas et al., Reference Martín Rojas, Alfaro, Galindo-Zaldívar, Borque, García-Tortosa, Sanz de Galdeano and Avilés2023). This deformation is responsible for the formation of several active faults and folds and for the high mean altitude of the GBB (ca. 1000 m above sea level).
Geologic map of the Guadix-Baza Basin (GBB) showing the locations of the dated deposits and paleontological sites mentioned in this study. FP1, Fonelas P-1 paleontological site; CB1, Cúllar-Baza 1 paleontological site; H1, Huéscar 1 paleontological site; SZ, Solana del Zamborino paleontological site; AT, Alicún travertines; ZT, Zújar travertines; PA, Puente Arriba fluvial terrace. The black traces represent active faults, including the Baza Fault.
Active structures deform the Pliocene and Quaternary infill of the GBB (Alfaro et al., Reference Alfaro, Delgado, Sanz de Galdeano, Galindo Zaldívar, García-Tortosa, López Garrido, López Casado, Marín, Gil and Borque2008, Reference Alfaro, Sánchez-Alzola, Martin-Rojas, García-Tortosa, Galindo-Zaldívar, Avilés and López Garrido2021; García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b, Reference García Tortosa, Alfaro, Galindo Zaldívar and Sanz de Galdeano2011; Fernández-Ibáñez et al., Reference Fernández-Ibáñez, Pérez-Peña, Azor, Soto and Azañón2010; Sanz de Galdeano et al., Reference Sanz de Galdeano, García-Tortosa, Peláez, Alfaro, Azañón, Galindo-Zaldívar, López Casado, López Garrido, Rodríguez-Fernández and Ruano2012; Castro et al., Reference Castro, Martin-Rojas, Medina-Cascales, García-Tortosa, Alfaro and Insua-Arévalo2018; Medina-Cascales et al., Reference Medina-Cascales, Martin-Rojas, García-Tortosa, Peláez and Alfaro2020; Fig. 3), leading to a heterogeneous basin architecture characterized by tectonically uplifted and lowered areas (García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar and Sanz de Galdeano2011). The most significant active structure is the east-dipping, normal Baza Fault (Fig. 3), whose total cumulative displacement has exceeded 2000 m since the late Miocene (Alfaro et al., Reference Alfaro, Delgado, Sanz de Galdeano, Galindo Zaldívar, García-Tortosa, López Garrido, López Casado, Marín, Gil and Borque2008). The slip rates reported for this active fault range between 0.3 ± 0.3 and 1.3 ± 0.4 mm/yr (Alfaro et al., Reference Alfaro, Delgado, Sanz de Galdeano, Galindo Zaldívar, García-Tortosa, López Garrido, López Casado, Marín, Gil and Borque2008, Reference Alfaro, Sánchez-Alzola, Martin-Rojas, García-Tortosa, Galindo-Zaldívar, Avilés and López Garrido2021; García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b, Reference García Tortosa, Alfaro, Galindo Zaldívar and Sanz de Galdeano2011). The regional ENE-WSW extension accommodated by these faults has characterized the geodynamic setting of the GBB during the Pliocene and Quaternary (Martin Rojas et al., Reference Martín Rojas, Alfaro, Galindo-Zaldívar, Borque, García-Tortosa, Sanz de Galdeano and Avilés2023). In the case of the Baza Fault, no significant change in the slip rate has been described during the Quaternary.
The displacement of the Baza Fault is responsible for some of the major geologic and geomorphic features of the GBB (García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b, Reference García Tortosa, Alfaro, Galindo Zaldívar and Sanz de Galdeano2011), the most remarkable being the division of the GBB into two sectors: the eastern sector (Baza Subbasin), located on the fault downthrown block, and the western sector (Guadix Subbasin), located on the fault upthrown block.
The endorheic–exorheic transition of the GBB
The Plio-Quaternary geologic, sedimentological, and geomorphological evolution of the GBB is divided into two main stages (Fig. 4): a first endorheic stage and a subsequent exorheic stage. The transition between the endorheic and exorheic stages took place when the GBB was captured by an outer fluvial system.
Sketches illustrating the Plio-Quaternary evolution of the Guadix-Baza Basin (GBB). During the initial endorheic stage, glacis developed throughout the entire basin. After the capture of the basin, the GBB became exorheic, and erosion has prevailed since that moment. pFR, Fardes paleo-River; FR, Fardes River; GQR, Guadalquivir River: GMR, Guadiana Menor River.
The endorheic stage and subsequent continental sedimentation in the GBB began at the end of the late Miocene, when the basin was disconnected from the Atlantic Ocean and the Mediterranean Sea (Soria et al., Reference Soria, Fernádez and Viseras1999). During this endorheic stage, the Baza Fault conditioned the sedimentary environments of the basin. In the downthrown eastern sector, a lake (Baza paleo-Lake) developed (Alfaro et al., Reference Alfaro, Delgado, Sanz de Galdeano, Galindo Zaldívar, García-Tortosa, López Garrido, López Casado, Marín, Gil and Borque2008; García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b). This lake was gradually filled by a thick lacustrine sedimentary succession (Vera, Reference Vera1970; Peña, Reference Peña1985; Gibert et al., Reference Gibert, Ortí and Rosell2007a). Moreover, the upthrown western sector was dominated by fluvial systems and, thus, detrital sedimentation (Vera, Reference Vera1970; Viseras, Reference Viseras1991). The main river of the western fluvial system was the Fardes paleo-River (Calvache and Viseras, Reference Calvache and Viseras1997; Fig. 4). This river drained from W to E toward paleo-Lake Baza, the main depocenter of the basin (Fig. 4). The Fardes paleo-River outlet into the Baza paleo-Lake was located to the north of Jabalcon Mountain (Fig. 4), where an alternation of fluvial and lacustrine sediments is observed. This interdigitation is located in the transition zone between the western fluvial sector and the eastern lacustrine sector.
A top basin glacis developed in the GBB at the end of its endorheic stage (García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b, Reference García Tortosa, Alfaro, Galindo Zaldívar and Sanz de Galdeano2011), extending almost the entire basin (Figs. 3 and 4). This glacis presents a mixed depositional/erosive nature (García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b), with erosion dominating in the outer parts close to the surrounding mountains and deposition prevailing in the inner sectors of the basin. Due to this “mixed” characteristic, the endorheic deposits just below the glacis present different ages along the GBB. This geomorphological surface has been used as a marker to estimate tectonic deformation rates (García Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b, Reference García Tortosa, Alfaro, Galindo Zaldívar and Sanz de Galdeano2011; Fernández Ibáñez et al., Reference Fernández-Ibáñez, Pérez-Peña, Azor, Soto and Azañón2010; Sanz de Galdeano et al., Reference Sanz de Galdeano, García-Tortosa, Peláez, Alfaro, Azañón, Galindo-Zaldívar, López Casado, López Garrido, Rodríguez-Fernández and Ruano2012) and fluvial incision rates in the GBB (Pérez-Peña et al., Reference Pérez-Peña, Azañón, Azor, Tuccimei, Della Seta and Soligo2009).
The glacis surface represents the youngest remaining feature of the endorheic stage. At some point during the middle Pleistocene), the exorheic stage began when the GBB was captured by the Guadiana Menor River (GMR), a tributary of the Guadalquivir River (Fig. 4) (Calvache and Viseras, Reference Calvache and Viseras1997). Previous works focused on this capture process propose, for instance, that the area from which the basin was captured was controlled by tectonics (Calvache and Viseras, Reference Calvache and Viseras1997; García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b; Moral and Balanyá, Reference Moral and Balanyá2020). However, no further details about the development of the capture process have been provided.
Once the basin was captured, its internal drainage was opened toward the Atlantic Ocean. The capture implied a major drop (more than 500 m) in the base level of the GBB drainage system, triggering intense headward erosion and fluvial incision processes that have dominated the basin since that moment (García-Tortosa et al., Reference García Tortosa, Sanz de Galdeano, Alfaro, Jiménez-Espinosa, Jiménez-Millán and Lorite Herrera2008a, Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b). During the exorheic stage, sedimentation was restricted to small alluvial systems around the basin borders and to valley bottoms, resulting in the formation of several fluvial terraces at different elevations. In addition, various travertine systems, such as the Alicún travertines (Díaz-Hernández and Juliá, Reference Díaz-Hernández and Juliá2006) and the Zújar travertines, precipitated and covered some of these exorheic fluvial terraces.
The controversial age of the capture
The age of the GBB's capture must be constrained between the youngest endorheic deposit and the older exorheic deposit. To date, several works have focused on one or another of the abovementioned constraints (Peña, Reference Peña1985; Calvache and Viseras, Reference Calvache and Viseras1997; Ortiz et al., Reference Ortiz, Torres, Llamas, Canoira, García-Alonso, García de la Morena and Lucini2000; Díaz-Hernández and Juliá, Reference Díaz-Hernández and Juliá2006; Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b; García-Tortosa et al., Reference García Tortosa, Sanz de Galdeano, Alfaro, Jiménez-Espinosa, Jiménez-Millán and Lorite Herrera2008a, Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b; Scott and Gibert, Reference Scott and Gibert2009). These different approaches led to controversial results, as the proposed ages for the capture range between 17 ka (Calvache and Viseras, Reference Calvache and Viseras1997) and 600 ka (Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b; Scott and Gibert, Reference Scott and Gibert2009). Therefore, the age of the GBB's capture remains an open scientific debate. In this section, we describe the works focused on this debate.
Most of the research regarding the timing of capture primarily relies on lower constraints. In this way, the first group of studies focused on dating the most recent endorheic deposits or the age of the glacis surface. These ages are determined by studying several paleontological sites located at the uppermost part of the endorheic sedimentary succession, very close to the glacis.
The most used paleontological site in relation to capture has been the Solana del Zamborino site (Botella, Reference Botella1975; Botella et al., Reference Botella, Vera and Porta1976; Casas et al., Reference Casas, Ruano and Torres1976; Fig. 3). Pioneer works dated this site as ca. 100 ka (Botella, Reference Botella1975; Botella et al., Reference Botella, Vera and Porta1976) because of the presence of the Acheulian stone tool industry. This age was widely used as a lower constraint for the GBB's capture (Peña, Reference Peña1985; Vera et al., Reference Vera, Rodríguez, Guerra and Viseras1994). For instance, Calvache and Viseras (Reference Calvache and Viseras1997) proposed that capture occurred between 100 and 17 ka, which corresponds to the age of several exorheic deposits (Jiménez de Cisneros, Reference Jiménez de Cisneros1994). More recently, Scott and Gibert (Reference Scott and Gibert2009) postulated an age of 750–770 ka for the Solana del Zamborino site based on the position of the Matuyama-Brunhes paleomagnetic reversal (ca. 781 ka), found just below this site, and sedimentation rates. In addition, Scott and Gibert (Reference Scott and Gibert2009), using their sedimentation rates, proposed an approximate age of 600 ka for the glacis surface and thus for the basin's capture. Using a similar approach but different sedimentation rates, Álvarez Posada et al. (Reference Álvarez Posada, Parés, Sala, Viseras and Pla-Pueyo2017) proposed an age of 480–300 ka for the Solana del Zamborino site.
Another paleontological site used to constrain the age of the capture was Cúllar-Baza 1 (CB1), which is in the eastern sector of the GGB (Fig. 3). This site was initially dated as younger than ca. 750 ka (Ruiz Bustos, Reference Ruiz Bustos1976; Alberdi and Bonadonna, Reference Alberdi and Bonadonna1989; Agustí et al., Reference Agustí, Oms and Parés1999; Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b; among others). This age agrees with later magnetostratigraphic analyses that placed the Matuyama-Brunhes reversal (ca. 781 ka) below the CB1 site (Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b). Furthermore, Gibert et al. (Reference Gibert, Scott, Martin and Gibert2007b) propose an age of ca. 600 ka for the glacis in the CB1 sector using the location of the Matuyama-Brunhes reversal ca. 19 m below this surface.
Further research regarding the age of the endorheic deposits was carried out by Azañón et al. (Reference Azañón, Tuccimei, Azor, Sánchez-Almazo, Alonso-Zarza, Soligo, Pérez-Peña, Alonso-Zarza and Tanner2006). These authors propose an estimated age of 43 ka for the capture event, based on 234U/230Th dating of a calcrete paleosoil located at the top surface of the basin in the western sector of the GBB, assuming that this calcrete was formed during the endorheic stage of the basin.
A second group of studies relied on dating the oldest exorheic deposits to establish a minimum age of the capture. Ortiz et al. (Reference Ortiz, Torres, Llamas, Canoira, García-Alonso, García de la Morena and Lucini2000) dated an exorheic fluvial terrace with amino acid racemization. These authors postulated that GBB capture occurred earlier than 239 ka. Díaz-Hernández and Juliá (Reference Díaz-Hernández and Juliá2006) did not directly propose an age for fluvial capture but estimated a time span for the development of the glacis and the river incision. For this purpose, they conducted 234U/230Th dating on several travertine platforms and calcretes in the Alicún travertines, found in the western sector of the GBB (Fig. 3, AT). In the case of these exorheic travertines, they obtained an age of ca. 220–190 ka for a platform located in a valley 190 ± 10 m below the glacis surface. For the calcrete, formed at the glacis level, the oldest sample provided an age of ca. 350 ka. From these data, they propose a timing between 350 and 205 ka for glacis development and between 115 and 48 ka for valley incision and the formation of erosive landforms.
The age of the capture has also been quantified using indirect criteria (García-Tortosa et al., Reference García Tortosa, Sanz de Galdeano, Alfaro, Jiménez-Espinosa, Jiménez-Millán and Lorite Herrera2008a, Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b). Using the fault slip rate of the Baza Fault and the glacis offset induced by this structure, García-Tortosa et al. (Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b) postulate a minimum age of 400 ka for the GBB capture.
New Ages of Exhoreic Deposits of The GBB
Dating methods
We numerically dated exorheic deposits located along the course of the Guadiana Menor, the main river in the GBB, to better constrain the age of the capture. The studied exorheic deposits were the Zújar travertines and the Puente Arriba fluvial terrace (Figs. 3 and 5–8).
(a) Geologic map of the Zújar travertine platforms. (b) Geologic cross sections along the Zújar travertines (location in a). (c) Topographic profile showing the position of the Zújar travertines related to the glacis and the present thalweg. GBB, Guadix-Baza Basin.
(a) Panoramic view showing the stepped arrangement of the Zújar travertine platforms formed in the Guadiana Menor River valley. (b) Detail of a travertine platform deposited over the exorheic detrital sediments of a previous fluvial terrace.
(a) Geologic map of the Puente Arriba fluvial terraces. (b) Geologic cross section along the Puente Arriba fluvial terraces (location in a). (c) Topographic profile showing the position of the dated terrace (T2) in relation to the glacis and the present thalweg.
(a) Panoramic view of the Puente Arriba fluvial terrace. (b) Detail of the sampling site (upper level of fine sediments).
We dated the Zújar travertines using the uranium-series disintegration (234U/230Th) method (Fig. 5, Table 1). For this purpose, three samples of the travertine platforms were collected and radiometrically dated at the Geochronology Laboratory of Geosciences Barcelona (GEO3BCN–CSIC) (Table 1). The radiometric ages were obtained through alpha-spectrometry using an ORTEC OCTETE PLUS spectrometer equipped with eight BR-024-450-100 detectors. The chemical separation of the radioisotopes and purification from travertine samples (~20 g) were conducted following the procedure described by Bischoff et al. (Reference Bischoff, Juliá and Mora1988), and isotope electrodeposition was performed according to the method of Talvitie (Reference Talvitie1972), modified by Hallstadius (Reference Hallstadius1984). Absolute ages were obtained employing the software designed by Rosenbauer (Reference Rosenbauer1991).
List of samples and numerical ages obtained for the Zújar travertine platforms and the Puente Arriba Terrace using U/Th and optically simulated luminescence (OSL) methods, respectively.
In the Puente Arriba Terrace (PAT), we collected one sample 1.5 m below the top of this terrace for optically stimulated luminescence (OSL) dating (Fig. 7, Table 1). This sample was dated at the Laboratory of Radioisotopes at the University of Seville. Dose rates were based on the average radionuclide activities of bulk material from each sample. High-resolution gamma spectrometry was used to measure the concentrations of 238U, 232Th, and 40K. Appropriate conversion factors (Adamiec and Aitken, Reference Adamiec and Aitken1998) were then used to derive the dose rates. A water content of 5 ± 2% was considered representative of the burial time. This value was used to calculate the attenuation of the dose. The contribution of cosmic radiation to the total dose rate was calculated as a function of latitude, altitude, burial depth, and average overburden density based on data by Prescott and Hutton (Reference Prescott and Hutton1994). Equivalent dose (D e) values were derived from the OSL measurements of quartz grain sizes ranging from 180 to 250 μm extracted from the sample. We measured 24 to 48 multigrain aliquots (~30 grains/aliquot) by applying the SAR blue OSL protocol (Murray and Wintle, Reference Murray and Wintle2000).
Numerical ages of the Zújar travertines
The Zújar travertines (Figs. 5 and 6) are in the western sector of the GBB, very close to the border with the eastern sector (Fig. 3). These travertines are located ca. 35 km upstream from the area where the GMR captured the former fluvial network of the basin. This group of travertine platforms was deposited on the left bank of the GMR and appear to be related to a hydrothermal spring, currently located 1 km to the west. This spring is characterized by a water output temperature up to 38°C and a discharge rate of 180 L/s (Cruz-Sanjulián and García-Rosell, Reference Cruz-Sanjulián and García-Rossell1972). At this position, the river valley is steeply incised into endorheic deposits. The travertine structure is made up of 10 carbonate platforms in a stepped arrangement (platform travertines P0 to P9 in Fig. 5a) deposited on a slight slope toward the river valley and positioned between ca. 300 and 250 m below the glacis surface. Some of the travertine bodies seem to partially overlap, although others remain individualized (Fig. 5b). Most of these platforms lie unconformably over Plio-Pleistocene endorheic sediments. However, some of them precipitated over previous exorheic deposits consisting of fluvial terraces of the GMR valley after its fluvial incision (Fig. 5a).
According to the derived 234U/230Th ages, the oldest travertine is platform P8 (sample ZT-3), which is located 250 ± 10 m below the glacis surface and 70 ± 1 m above the present thalweg of the GMR (Fig. 5). The P8 travertine was dated to 240.8 +29.27/−24.04 ka.
Platform P2 (sample ZT-2) is located 275 ± 10 m below the glacis surface and 45 ± 1 m above the present thalweg (Fig. 5). We obtained a radiometric date of 109.04 +6.77/−6.43 ka for P2.
The lowest travertine platform (P1, sample ZT-1) is located 300 ± 10 m below the glacis surface and 20 ± 1 m above the present thalweg (Fig. 5). Samples from this lower terrace yielded an age of 78.87 +4.53/−4.37 ka.
Numerical age of the PAT
The PAT (Figs. 7 and 8) is a fluvial terrace found 65 km upstream from the capture area and located in the eastern sector of the GBB (Fig. 3) on the NW side of the GMR valley. It forms part of a set of four fluvial terraces (T1 to T4, from older to younger) deposited at different heights with respect to the present thalweg of the GMR (Fig. 7). Unfortunately, T1 facies were not appropriate for OSL dating. The dated PAT (T2) is located 186 ± 10 m below the glacis surface and 30 ± 1 m above the thalweg (Fig. 7). It is an unpaired fill fluvial terrace that reaches a thickness of up to 10 m. The terrace is partially eroded and composed of clast-supported conglomerate and gravel deposits. Locally, sandy levels can be observed. The exorheic deposit of the PAT unconformably overlies Lower Pleistocene endorheic deposits. OSL dating of the PAT provided an age of 89.5 ± 5.3 ka (Table 1).
When Was The GBB Captured? Insights From New Exhoreic Ages
In this section, we discuss the meaning of our new numerical dates in terms of the age of the GBB capture. The time frame in which an internally drained basin is captured is constrained between the age of the youngest endorheic deposits and the oldest exorheic ones (Figs. 9–11). Therefore, to accurately date capture processes, effort must be put into dating these constraining horizons.
Numerical age and position of the exorheic deposits dated or included in the discussion of this study. They are depicted according to their depths below the glacis surface and the distances to the capture area, which is approximately the confluence between the Guadiana Menor and Fardes rivers (blue and purple longitudinal profiles). AT, Alicún travertines; ZT, Zújar travertines; PAT, Puente Arriba Terrace.
Sediment thickness between the glacis and the Matuyama-Brunhes reversal in different stratigraphic successions of the endorheic infilling of the Guadix-Baza Basin. The positions of the paleontological sites within these successions are indicated, along with the authors who identified the Matuyama-Brunhes reversal in each site (Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b; Scott and Gibert, Reference Scott and Gibert2009; Pla-Pueyo et al., Reference Pla, Viseras, Soria, Tent-Manclús and Arribas2011; Álvarez-Posada et al., Reference Álvarez Posada, Parés, Sala, Viseras and Pla-Pueyo2017).
Chronological table of the Guadix-Baza Basin (GBB) capture event. The left side presents the different age proposals from previous works (Peña, Reference Peña1985; Vera et al., Reference Vera, Rodríguez, Guerra and Viseras1994; Calvache and Viseras, Reference Calvache and Viseras1997; Ortiz et al., Reference Ortiz, Torres, Llamas, Canoira, García-Alonso, García de la Morena and Lucini2000; Díaz-Hernández and Juliá, Reference Díaz-Hernández and Juliá2006; Azañón et al., Reference Azañón, Tuccimei, Azor, Sánchez-Almazo, Alonso-Zarza, Soligo, Pérez-Peña, Alonso-Zarza and Tanner2006; Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b; García-Tortosa et al., Reference García Tortosa, Alfaro, Galindo Zaldívar, Gibert, López Garrido, Sanz de Galdeano and Ureña2008b; Scott and Gibert, Reference Scott and Gibert2009). The right side illustrates the time range we propose in this study for the capture process. The upper constraint of this range is the age of the oldest dated exorheic deposits presented in this work, with the additional time estimate for the drainage network to dissect the valley to the position of these deposits (ca. 250 m). The lower constraint is the Matuyama-Brunhes horizon plus the time required for the sedimentation of the thickness of endorheic deposits between the paleomagnetic reversal and the glacis. The ages of the youngest endorheic deposits are also supported by paleontological data (green bar).
New numerical ages acquired in this work fell in a time range between 24.08 +29.26/−24.04 and 78.87 + 4,53/−4,37 ka (Table 1). We thus provided the oldest numerical age to date for an exorheic deposit within the GBB. As described earlier, platform P8 is an exorheic deposit cropping out within the valley of the GMR (Figs. 5 and 9). Therefore, this travertine body was formed when the valley was already dissected approximately 250 m below the glacis surface. This age agreed with the ages of other exorheic travertines of the GBB such as the Alicún travertines (Díaz-Hernández and Juliá, Reference Díaz-Hernández and Juliá2006). The uppermost Alicún travertine platform is located 190 ± 10 m below the glacis and was dated as ca. 220–190 ka (Fig. 9).
Hence, the age of the Zújar travertine platform P8 represents an upper constraint for the fluvial capture of the GBB. This implies that the capture should be older than 240.8 +29.26/−24.04 ka. However, to estimate the age of the fluvial capture, it is necessary to add to this age the time span required by the drainage network to dissect 250 m until the position of the travertine platform (Figs. 5, 9, and 11). Calculating this time span using incision rates may, however, be problematic, as they are highly sensitive to tectonics, climate changes, and local base-level variations. Extensive literature has proven that these forcings do not act linearly through time, which in turn implies a lack of linearity of incision rates along time and space, especially when dealing with time spans of more than hundreds of thousands of years (Pazzaglia et al., Reference Pazzaglia, Gardner and Merritts1998; Whipple, Reference Whipple2001; Faust and Wolf, Reference Faust and Wolf2017). Table 2 shows incision rates estimated from the samples’ position above the present thalweg, ranging between 0.2 and 0.6 mm/yr. If we use these incision rates to estimate the age of the capture assuming that they are constant, we obtain ages older than 780 ka. These ages of the capture assuming constant incision rates are not in agreement with the sedimentary record of the basin. Therefore, we think that the incision rates in the GBB are not linear but varied in time. Consequently, due to the high uncertainty derived from the use of incision rates, we considered that these rates were not suitable to estimate the age of the glacis and, therefore, of the basin capture. In any case, the Zújar travertines may have indicated an early capture event. This early capture was supported by the age of the PAT (89.5 ± 5.3 ka). This age implied that by this time, the GMR had already dissected 186 ± 10 m below the glacis surface at a distance of more than 60 km from the capture area (Fig. 9).
Incision rates estimated from the age of the samples and their position with respect to the current thalweg.
a From the obtained incision rates, we estimate the time span required by the drainage network to dissect vertically to the position of the samples.
b This time span plus the age of the sample would give us an alleged age for the capture.
On the other hand, the age of the most recent endorheic deposits can provide a lower constraint on the age of the capture. Several controversial ages have been proposed for the younger horizons of the GBB. However, there is one datum consistently accepted for authors who worked in the GBB: the presence of the Matuyama-Brunhes reversal dated at ca. 781 ka. This paleomagnetic boundary is identified in several stratigraphic sections of the basin that also contain major paleontological sites. The Fonelas P-1 site (Arribas et al., Reference Arribas, Riquelme, Palmqvist, Garrido, Hernández, Laplana and Soria2001, Reference Arribas, Garrido, Viseras, Soria, Pla-Pueyo, Solano, Garcés, Beamud and Carrión2009; Viseras et al., Reference Viseras, Soria, Durán, Pla, Garrido, García-García and Arribas2006; Pla-Pueyo et al., Reference Pla, Viseras, Soria, Tent-Manclús and Arribas2011) and the Solana del Zamborino site (Botella, Reference Botella1975; Botella et al., Reference Botella, Vera and Porta1976; Casas et al., Reference Casas, Ruano and Torres1976; Martín Penela, Reference Martín Penela1988; Scott and Gibert, Reference Scott and Gibert2009; Álvarez-Posada et al., Reference Álvarez Posada, Parés, Sala, Viseras and Pla-Pueyo2017) are both in the western sector, while the CB1 and Huéscar 1 sites are both in the eastern sector of the basin (Ruiz-Bustos, Reference Ruiz Bustos1976, Reference Ruiz-Bustos1984; Mazo et al., Reference Mazo, Sesé, Ruiz-Bustos and Peña1985; Alberdi and Bonadonna, Reference Alberdi and Bonadonna1989; Agustí et al., Reference Agustí, Oms and Parés1999; Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b; among others; Figs. 3 and 10). Figure 10 depicts the thickness of endorheic sediments over the Matuyama-Brunhes reversal and below the glacis surface. This thickness was ca. 7 m in Fonelas P-1, ca. 15 m in Solana del Zamborino, ca. 19 m in CB1, and ca. 15 m in Huéscar-1. Therefore, to estimate the age of the youngest endorheic deposits, we needed to add to the 781 ka of the paleomagnetic reversal time span necessary to deposit between 7 and 19 m of endorheic deposits (Fig. 10). Unfortunately, there is no consensus on sedimentation rates in the GBB. An example of this last statement arose in the Solana del Zamborino paleontological site. Scott and Gibert (Reference Scott and Gibert2009) conducted a magnetostratigraphic analysis in this stratigraphic succession, identifying the Matuyama-Brunhes 15 m below the glacis. They used a sedimentation rate of 10 cm/ka to calculate the time span between the polarity reversal and the stratigraphic position of the site, obtaining an age of 770–750 ka for the Solana del Zamborino site. This sedimentation rate is obtained from a paleomagnetic and stratigraphic study in Cúllar (Fig. 3), in the easternmost part of the basin (Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b). Using a similar approach (paleomagnetism and sedimentation rates), Álvarez Posada et al. (Reference Álvarez Posada, Parés, Sala, Viseras and Pla-Pueyo2017) proposed an age of 480–300 ka for the Solana del Zamborino site based on a sedimentation rate of ca. 2 cm/ka. However, it has to be considered that this rate was obtained next to the Fonelas P1 site (Pla-Pueyo et al., Reference Pla, Viseras, Soria, Tent-Manclús and Arribas2011), where sedimentary facies (lacustrine facies) are different from those found in the Solana del Zamborino site (fluvial conglomerate facies). The use of sedimentation rates for calculating this time span is an unreliable approach, because it mainly depends on which number is selected and employed. This large disparity in the two proposed ages (Scott and Gibert, Reference Scott and Gibert2009; Álvarez Posada et al., Reference Álvarez Posada, Parés, Sala, Viseras and Pla-Pueyo2017) prevents the use of this site to constrain the age of the GBB's capture.
Other data that may contribute to better constrain the age of the recent endorheic deposits of the GBB are the overall ages of the paleontological sites existing in the basin. According to the faunal assemblages collected within the endorheic sediments, only 3 out of more than 150 sites in the basin are either close to or younger than the Matuyama-Brunhes reversal: the Caniles, CB1, and Solana del Zamborino sites (Ruíz-Bustos, Reference Ruiz-Bustos1984; Guerra Merchán and Ruiz Bustos, Reference Guerra Merchán and Ruiz Bustos1992; Scott and Gibert, Reference Scott and Gibert2009; Álvarez-Posada et al., Reference Álvarez Posada, Parés, Sala, Viseras and Pla-Pueyo2017). The other paleontological sites have ages much older than this paleomagnetic reversal. In addition, Demuro et al. (Reference Demuro, Arnold, Parés and Sala2015) reported ages of 570–420 ka for Huéscar 1 based on OSL analyses. However, other studies focusing on this site suggest an older age of 781 ka, supported by paleontological (Mazo et al., Reference Mazo, Sesé, Ruiz-Bustos and Peña1985; Alberdi and Bonadonna, Reference Alberdi and Bonadonna1989) and magnetostratigraphic data (Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b). Therefore, the Caniles, CB1, and Solana del Zamborino sites are considered the most recent sites in the basin (Maldonado-Garrido et al., Reference Maldonado-Garrido, Piñero and Agustí2017).
The Caniles site has an assigned age of ca. 781 ka based on its vertebrate faunal content (Guerra Merchán and Ruiz Bustos, Reference Guerra Merchán and Ruiz Bustos1992). The age of the CB1 site has been assigned as middle Pleistocene in different studies. Some authors propose an age between 500 and 750 ka based on its faunal assemblage (e.g., Ruiz Bustos, Reference Ruiz Bustos1976; Alberdi and Bonadonna, Reference Alberdi and Bonadonna1989), while others propose an age of ca. 781 ka using paleomagnetism (Gibert et al., Reference Gibert, Scott, Martin and Gibert2007b). However, considering the stratigraphic proximity of the site to the Matuyama-Brunhes reversal, it is more likely that its age is closer to that proposed by Gibert et al. (Reference Gibert, Scott, Martin and Gibert2007b).
In relation to calcrete ages of 350 ka (Díaz-Hernández and Juliá, Reference Díaz-Hernández and Juliá2006) and 42 ka (Azañón et al., Reference Azañón, Tuccimei, Azor, Sánchez-Almazo, Alonso-Zarza, Soligo, Pérez-Peña, Alonso-Zarza and Tanner2006), it is necessary to consider that calcretes could have formed while the glacis was active or after it was abandoned. Consequently, we consider these ages are not reliable to estimate the age of the capture.
In conclusion, we consider that there are only two data sets robust enough to quantitatively constrain the age of the GBB's capture. The upper quantitative constraint is our age of 240.8 +29.27/−24.04 ka for the Zújar travertine platform P8 (the oldest dated exorheic deposit). The lower quantitative constraint is 781 ka, owing to the polarity reversal present in the upper part of the endorheic sedimentary succession.
The time span between ca. 240 and 781 ka could be refined by adding (1) the time necessary to dissect 250 m below the glacis until the position of the P8 travertine (Fig. 9) and (2) the time span necessary to sediment 7 to 19 m of endorheic deposits (sediment thickness between the Brunhes-Matuyama reversal and the glacis) (Fig. 10). As discussed earlier, we consider that incision rates are not suitable to quantitatively estimate the age of the basin capture. We also discussed that the controversy related to the sedimentation rates hinders a quantitative approach to calculate the age of the glacis. Further sedimentological and geomorphological analyses would be necessary to overcome these limitations. Therefore, we consider that a further refinement of the 781–240 ka time span using incision and sedimentation rates can only be addressed qualitatively.
Additional qualitative data for the lower constraint are the absence of paleontological sites younger than ca. 750 ka, except for the controversial age of the Solana de Zamborino site. Based on this, we postulate that the capture could have occurred close to the lower constraint, that is, close to 781 ka, well before the oldest proposed age of 600 ka by Gibert et al. (Reference Gibert, Scott, Martin and Gibert2007b).
An earlier capture would not contradict, however, the age of 600 ka proposed by Gibert et al. (Reference Gibert, Scott, Martin and Gibert2007b) for the glacis in the eastern sector. At the end of the endorheic stage, the headward erosion of a tributary of the Guadalquivir River reached the divide between the GBB and the Guadalquivir Basin. At that moment, erosion started in the GBB, leading to a new river, the Guadiana Menor. At this early stage, erosion in the GBB was initially constrained to a small area around the capture area. The first phase of the capture process was initiated when the abovementioned tributary of the Guadalquivir River reached a first river of the GBB fluvial network. As the headward erosion of the GMR proceeded, it eventually reached the main river of the western sector of the GBB, that is, the Fardes paleo-River. This moment was a milestone in the evolution of the basin, as it implied the capture of most of the drainage network of the western sector.
We hypothesize that the capture of the GBB was not a simple event because of the basin configuration related to the presence of the Baza Fault. Most likely, endorheic conditions persisted for a longer period in the eastern sector of the basin, located in the downthrown block of the Baza Fault. This fault could have kept the eastern sector of the GBB uncaptured, allowing continuous sedimentation, while the western sector of the basin was already captured. Subsequently, the eastward-migrating erosion along the western sector reached the Baza paleo-Lake. Consequently, the rivers previously draining toward the lake were incorporated into the new fluvial network, which drained to the west through the GMR. At this moment, the entire GBB became exorheic. An approximation to the evolution of the capture process in the basin could be deduced by comparing the available ages of exorheic deposits with their distance to the capture area and their position below the glacis surface (Fig. 9). However, we consider that relying solely on three dated exorheic deposits is insufficient for conducting this type of discussion. A further analysis, focused on dating a greater number of exorheic deposits throughout the entire GBB, would be necessary in the future to establish a much more precise evolutionary model along time and space.
Conclusions
In this study, we provide new numerical data to constrain the age of Guadix-Baza Basin capture. We dated a series of exorheic deposits that crop out within the valley of the Guadiana Menor River: three platforms of the Zújar travertines (P1, P2, and P8) and the Puente Arriba fluvial terrace.
The oldest Zújar travertine platform is dated to 240.8 +29.27/−24.04 ka. It corresponds to a platform at a position of 250 ± 10 m below the glacis and 70 ± 2 m above the present thalweg. On the other hand, the OSL dating of the Puente Arriba terrace provided an age of 89.5 ± 5.3 ka.
We thus provide the oldest age recorded to date for an exorheic deposit in the basin (ca. 240 ka), establishing a new upper constraint for Guadix-Baza Basin capture. We infer that basin capture took place before ca. 240 ka plus the additional time necessary for the Guadiana Menor River to incise 250 m down to the position of the Zújar travertine platform. Furthermore, our new dating, together with previous data, support the possibility that the capture occurred earlier than previously suggested. First, the presence of the Matuyama-Brunhes reversal, dated ca. 781 ka, close to the top of the endorheic succession of the basin represents a robust quantitative lower constraint. Additionally, the ages of the paleontological sites throughout the basin consistently fall within the range of approximately 750 ka or older (except for some ages proposed for the Solana del Zamborino site). Therefore, although we do not provide a precise age of the capture, several findings suggest that the capture event likely predates even the oldest proposed ages of 600 ka.
Funding
This research was funded by the Spanish Ministry of Science, Innovation and University (Research Projects RTI2018-100737-B-I00 and PID2021-127967NB-I00), Junta de Andalucía regional government (RNM325 research group), Generalitat Valenciana Research Project AICO/2021/196, and the University of Alicante (Research Project VIGROB053).
Competing Interests
The authors declare no competing interests.
Open access
