Plio‑Pleistocene hydrothermal events of the Baza Basin (Betic Cordillera, SE Spain) and their paleoecological implications

The Lower Pliocene to upper Middle Pleistocene continental sedimentary infillings of the Baza Basin (Guadix-Baza Depression, SE Spain) are worldwide unique by their high stratigraphic completeness and exceptional preservation of their fossil record of terrestrial vertebrates. These sediments were deposited in fluvio-lacustrine environments and preserve huge assemblages of large mammals, including the oldest evidence of hominin presence in Western Europe at the late Lower Pleistocene sites of Barranco León and Fuente Nueva-3, dated to ∼ 1.4 Ma. Since latest Miocene times, the basin was subject to intense tectonic activity, showing sedimentary deposits originated from ancient hot springs. These hot springs were developed in a spatio-temporal network related to tectonic fractures, showing a higher concentration in the Orce sub-basin during Late Pliocene to Holocene times. Compositional analyses of cherts and other facies associated with these hot springs show geochemical and mineralogical markers linked to hydrothermal scenarios. The contribution of hot springs in the Baza Basin resulted in rich and productive biotopes, which supported a high diversity and biomass of terrestrial mammals. Similar geological, hydrochemical and ecological environments are found along hotspot lines in East Africa, the Levantine Cor-ridor and the circum-Mediterranean realm. We hypothesize here that the presence of these hot springs corridors favored the dispersal of hominins and other large mammals during the Plio-Pleistocene, which led to the first hominin arrival in Europe.


Geological and paleontological setting
According to the classification of Miall (1984), the Baza Basin is an intramountainous basin with tectonic subsidence and continental sedimentation, features inherited from an ancient marine basin formed during the Late Miocene in the interior of the Betic mountain range (Platt et al., 2013;Sanz de Galdeano & Vera, 1992).The transition from a marine to a continental sedimentary regime took place 7.5-7.0Ma ago, being the result of major regional tectonic movements that caused the isostatic uplifting of the basin (average uplift rate of ~ 200 m/Ma, as estimated from Upper Neogene coastal marine conglomerates and coral reefs; Braga et al., 2003), which resulted in its disconnection from the Mediterranean Sea (Soria, 1993;García-Aguilar, 1997;Playá et al., 2000;Guerra-Merchán et al., 2014).The Baza Basin is bounded to the north by reliefs of Mesozoic age, mostly carbonates (Thetys type), which belong to the External Zones of the Betic Cordillera.To the south, it is surrounded by reliefs composed of sedimentary and metamorphic rocks of Paleozoic-Triassic age, which belong to the Internal Zones.
To the west, it is connected with the Neogene Guadix Basin and to the southeast with the Neogene Almanzora Corridor Basin (Fig. 1B).
The first TSU, which is of Pliocene age, shows a sediment thickness of ~ 150 m.It is characterized by the alternation of marls/calcilutites/lutites and limestones deposited in lacustrine systems, which change towards the edges of the basin to sands and conglomerates of alluvial origin.Forest ecosystems were formed on the margins of the lake as inferred for the Ruscinian (~ 4.5 Ma) Baza-1 paleontological site (Ros-Montoya et al., 2017).Tectonic activity during this stage induced the formation of many soft-sediment deformation structures (i.e., paleo-seismites), which are mainly located at the marginal sites of a paleo-lake (Alfaro et al., 1997) and possibly resulted in deposits with hydrothermal affinities (García-Aguilar et al., 2014, 2015).Climate was warm and humid during this stage, which relates to an environment typical of modern subtropical or rainy tropical climates (Foucault & Mélières, 2000;Haywood & Valdés, 2004).
The upper section of this TSU outcrops mainly in the Orce area, a sub-basin in the NE sector of the Baza Basin, which has an extent of ~ 170 km 2 (Fig. 1B, C) and shows ~ 35 m thick extensive outcrops that correspond to the Lower to Middle Pleistocene.The Orce sub-basin preserves a unique fossil record of terrestrial vertebrates, including two archaeological sites dated to ~ 1.4 Ma (see review in Palmqvist et al., 2016) that mark the earliest presence of hominins in Western Europe: Barranco León (BL) and Fuente Nueva-3 (FN3).BL and FN3 preserve huge assemblages of skeletal remains of large mammals and Oldowan tools, which are linked through the presence of anthropogenic cut and percussion marks, as well as a human deciduous tooth in the case of BL, which is at this moment the oldest hominin fossil found in Western Europe (Espigares et al., 2013(Espigares et al., , 2019;;Martínez-Navarro et al., 1997, 2014;Oms et al., 2000aOms et al., , 2000b;;Palmqvist et al., 2023;Titton et al., 2021;Toro-Moyano et al., 2013;Yravedra et al., 2021Yravedra et al., , 2022)).Moreover, the Venta Micena site (VM, 1.6-1.5 Ma), which is slightly older than BL/FN3 and has been interpreted as a denning area of the giant hyena Pachycrocuta brevirostris (Arribas & Palmqvist, 1998;Palmqvist & Arribas, 2001;Palmqvist et al., 2022b), provides an excellent record of the large mammal community that preceded the hominin dispersal in the basin (Martínez-Navarro, 1991;Mendoza et al., 2005;Rodríguez-Gómez et al., 2017).The varied lithological composition of the Orce area, with alternations of fluvial red detrital facies and lacustrine carbonate-marly facies that locally show gypsum deposits (restricted to the center of the lacustrine system in this area) and cherts, points to a paleogeographic model of fluvio-lacustrine sedimentation highly variable in time and space.This resulted in open plains with forest patches in the emerged areas surrounding the fluviolacustrine system (García-Aguilar & Palmqvist, 2011;Mendoza et al., 2005;Palmqvist et al., 2022a;Scott et al., 2007).
The third TSU of the Baza Basin has a Middle Pleistocene to Upper Pleistocene age, showing a maximum thickness of 70 m.During this stage, the activity of major faults (Baza and Galera) resulted in some tectonic blocks that developed an effective subsidence, which allowed a continuous sedimentation.This TSU is composed of conglomerates and sands towards the basin edges, which change to lutites, marls and marly limestones deposited in shallow lacustrine environments towards the basin center.In the upper part of this unit and at the basin edges, alluvial facies (mainly red conglomerates and sands) prograde over lacustrine facies.This TSU is represented in several sectors of the basin perimeter (Fig. 1B, C) and the most significant paleontological sites are Huéscar-1 and Cúllar-Baza-1 (Ruiz-Bustos, 1976;Peña et al., 1977;Mazo et al., 1985;Alberdi & Bonadonna, 1989;Agustí et al., 2000;Torrente Casado, 2010;Ros-Montoya et al., 2018).
Finally, during the late Middle Pleistocene (~ 400 ka: García-Tortosa et al., 2008a, 2008b) or in Late Pleistocene times (100-17 ka: Viseras & Fernández, 1992;Calvache & Viseras, 1997;42 ka: Azañón et al., 2006), the hydrographic network of the sedimentary depression was captured by the Guadiana Menor River, a tributary of the Guadalquivir River.This resulted in the transit from an endorheic regime to an exorheic one, leading to a stage in which erosion dominated over sedimentation.As a result, a badlands landscape predominates now in the region, linked to a mean annual precipitation (MAP) of ~ 300 mm in the innermost part of the sedimentary depression (https:// es.clima te-data.org).However, climate was more humid during the late Early Pleistocene, with an estimated MAP of ~ 700 mm in BL/ FN3 and ~ 850 mm in VM (Palmqvist et al., 2022c).

Hot springs currently active in the Baza Basin
Tectonic dynamism resulted in hot springs, several of which are currently active in the Baza Basin and other nearby areas (Cruz-Sanjulián et al., 1972;Cruz-Sanjulián & García-Rosell, 1972, 1975;Diputación de Granada-ITGE, 1990; Fig. 1A, B, and Table 2).There are twelve hot springs in the Baza Basin and four in the Guadix Basin, all of them associated with regional faults.In the case of the hot springs of the Baza Basin, their water output temperatures range between 18 °C and 22 °C, except in the case of the 'Baños de Zújar' spring, with 38 °C.The hydrochemistry of the low temperature springs shows a predominance of calcium bicarbonate.In contrast, the high temperature springs mainly show silica, chlorides, and sulphates.This correspondence between high water output temperatures and siliceous-sulphate-chloride hydrochemistry appears also in the > 30 °C hot springs of the Guadix Basin.The current climatic scenario of the Baza Basin shows an average annual temperature of 13.5 °C (13.8 °C in the Baza station and 13.3 °C in the Orce station; https:// es.clima te-data.org), with strong annual oscillations.Specifically, average temperature ranges between 5 °C and 9 °C from November to March, averaging ~ 25 °C during summer.Likewise, the historical record of annual rainfall averages 335 mm at the Baza station and 366 mm at the Orce station (https:// es.clima te-data.org).In any case, the average water output temperature of the hot springs in the Baza Basin is 21 °C, which is 7.5 °C higher than the average annual atmospheric temperature in the basin.These data indicate that the water output temperatures of the springs are slightly lower than the average temperatures in summer, but much higher than in autumn and winter.This underlines the thermal character of these hot springs, whose flow rate varies between 1 and 400 l/s.

Tectonic dynamism
From a tectonic point of view, the Baza Basin lies in the interior of the alpine-type Betic mountain range, showing three main fracture systems, which are NE-SW, NNW-SSE and NNE-SSW oriented (Sanz de Galdeano, 1983Galdeano, , 1990Galdeano, , 2008;;Sanz de Galdeano & Vera, 1991, 1992;Sanz de Galdeano et al., 2007, 2012).These fractures played an important role in the tectono-sedimentary evolution of the Betic mountain range and generated a deep steplike geometry, which is indicative of epicortical extension processes (Platt et al., 2013;Sanz de Galdeano & Vera, 1992).Sanz de Galdeano et al. (2007) established two depocenters in the Baza Basin, induced by NNW-SSE and N70ºE faults that resulted in a regional ENE extension, as deduced by gravity maps of regional, residual and Bouguer anomalies: −150 mGal in Benamaurel (central zone of Baza Basin) and −16 mGal at the north of the Orce village (central zone of Orce sub-basin).Considering basement densities of 2.67 g/cm 3 and 2.32 g/cm 3 for the sedimentary infillings of both sectors, respectively (Sanz de Galdeano et al., 2007), these result in basement depths of ~ 1300 m and ~ 500 m, respectively.In this sense, Sanz de Galdeano (2008) highlighted the Cadiz-Alicante N70ºE fault as the boundary of the inferred depocenters in the central Baza and Orce-North sectors.Likewise, Sanz de Galdeano et al., (2012Galdeano et al., ( , 2020) ) described using geological criteria the main active fractures present in the Baza Basin, which have most likely induced the historical earthquakes recorded in this region: the Galera Fault (23 km long and > 5 km deep), the Baza Fault (37 km long and > 10 km deep), the  The position of hydrothermal upwellings is shown in Figure 1B.Average composition values (in mg/l), temperatures, water flows and conductivities according to Diputación de Granada-ITGE (1990) and Cruz-Sanjulian et al. (1972).Hot springs A, B, C, and D are located in the Guadix Basin.The other twelve hot springs are placed in the Baza Basin.GeoT geothermometers estimated using the Na/K calculation method (Fournier, 1977;Fournier & Truesdell, 1973;Giggenbach et al., 1988) Botardo-Alfahuara Fault (22 km long and > 5 km deep), the Benamaurel Fault (19 km long and > 10 km deep), and the Cúllar Fault (North and West, ~ 9 km long and 5-10 km deep, deduced in this case only by seismic criteria).Seismic phenomena were frequent during the Plio-Pleistocene, as evidenced by the presence of paleo-seismites (Alfaro et al., 1997;2000;2010;Gibert et al., 2005).Therefore, considering a standard value of the upper crustal geothermal gradient of 30º C/km (Gupta & Roy, 2007;Harrison et al., 1990), the depth calculated for these faults and the possible hydrothermal phenomena associated with them, we estimate temperatures of at least 150°-300 °C in their deepest zones.

Objectives
The main goal of this study is to establish a causal link between the presence of hydrothermal upwelling-related deposits in the continental Plio-Pleistocene infillings of the Baza Basin, their tectono-sedimentary origin, and their paleoecological consequences.For doing so, we have: (i) reviewed and documented the tectono-sedimentary context of the Baza Basin; (ii) performed a detailed field study, sampling selected stratigraphic sections; (iii) studied the sedimentology of these sections, analyzing the geochemical and mineralogical composition of the samples; and (iv) performed both local and high scale (Africa and Eurasia) correlations of the relationship between tectonic activity and thermal springs, on the one hand, and the presence of archeological and paleontological sites, on the other.

Materials and Methods
We have reviewed and documented the tectono-sedimentary context of the main archaeological and paleontological sites of the Baza Basin, as well as the sedimentary and compositional models of many modern and ancient thermal springs worldwide distributed together with the paleoanthropological sites linked to these springs.Field work in the Baza Basin consisted of geological mapping and observations at a mesoscopic scale of lithologies and their sedimentary structures, as well as of the main faults affecting the sedimentary infillings of the basin.Detailed stratigraphic profiles were obtained and the different sedimentary facies present in them were sampled (Fig. 2).Extensive sediment sampling was carried out along the stratigraphic record surveyed (Fig. 1B, C, Fig. 2) to perform mineral and chemical analyses for identifying those facies potentially associated with a hydrothermal origin.
In this study, we have analyzed 56 samples located throughout the Baza Basin, which correspond to eight types of lithologies: calcilutites, lutites, dark lutites, marls, travertines, cherts, gypsum, and native sulfur.All samples were obtained in clean field sections.For each sample, ∼200 g of sediments were stored in an airtight bag, which was marked with the corresponding acronym (abbreviations of samples collected from sections/sites are indicated in Fig. 2).These samples were subjected in the laboratory to grinding and mixing in an agate mortar.Then, ∼20 g were randomly selected from each sample for performing the analyses.The samples were analyzed geochemically by means of X-ray fluorescence (XRF).In addition, 14 of these samples were selected for analyzing their mineralogy by means of X-ray diffraction (XRD).All measurements were carried out at the Central Services of Support to Research ('Servicios Centrales de Apoyo a la Investigación', SCAI) of the University of Málaga (UMA).

Elemental analysis
Chemical analyses of the abundance of major and some trace elements were carried out using a Thermo Fisher ADVANT'XP wavelength dispersive atomic fluorescence spectrometer (XRF).This equipment consists of an autosampler with capacity for a total of 161 samples, a Rh tube and a goniometer with three collimators and four crystals that allow measuring from F to U. The measurements were obtained in a He atmosphere with the sample in bead.The results were analyzed with Thermofisher's UniQuant software, which allows the semi-quantification of unknown materials without the need of standards or reference materials.Measurement errors ranged between 0.3% and 1.5% for major elements, raising to 5% in the case of minor elements.The detection limit of these analyses was of 0.01%, which was quoted in the corresponding tables as "bld" (i.e., below the detection limit).

Mineralogical analysis
Fourteen samples, selected a priori as potentially indicative of being generated in a hydrothermal context, were analyzed by X-ray Diffraction (XRD) based on the XRF data.The measurements were obtained with an X'Pert PRO MPD diffractometer from PANalytical and measured in high resolution conditions.The X-ray tube has a Cu anode and the diffractometer in the incident beam consists of a Johansson-type primary monochromator with a Ge (111) crystal, which provides a strictly monochromatic Cu K α1 radiation.Random powder X-ray diffraction (XRD) of bulk-rocks was used for the identification of mineral phases, measuring the diffraction angles and the intensities of the diffracted beams.This allows to compare the resulting diffraction pattern with a reference database.The comparison was carried out with PDF-4 + 2021 (Gates-Rector & Blanton, 2019).When the acquired data were of sufficient quality (i.e., measured with a wide range of angles, good resolution and high intensity/noise ratio), they were used not only to identify mineral phases but also to measure them in a semi-quantitative fashion.To determine the proportion of phyllosilicate in the global samples, the phyllosilicates (4.46 Å)/quartz (3.343 Å) peak-area ratio was estimated from XRD patterns of randomly oriented powder.D-spacing and the full width at half-high maximum (FWHM) were determined on the (101) reflection peak of the opal specimen in order to establish the type of opal polymorphs (Elzea et al., 1994).For mineral abundance, the results of the analyses were expressed semiquantitatively; the error of the estimates ranged from 5 to 10%.For each sample, a relative range of abundance was indicated for each mineral.
The field study showed the presence of 17 levels of cherts, all located in the Orce sub-basin, mainly to the south of the town of Galera (Fig. 1C).Some paleontological sites (e.g., FN-1, GA-2 and NU-2) are situated just a few decimeters above these siliceous deposits (Fig. 2).Chert beds typically occur in irregular shaped strata, with a lateral extent ranging from 7 m to > 200 m, and a thickness of 8-80 cm (Fig. 3A,  C).The cartographic relationship between the appearance of cherts and the presence of tectonic fractures is remarkable (Fig. 1C; Fig. 3G, I, J).Lithologically, these cherts appear as white, dark gray or light brown massive facies.At an outcrop scale, they show abundant diaclases, micropores, changes of coloration, and the presence of conchoidal fractures.In some cases, the cherts appear with decimeter to centimeter-scale gypsum mineralizations within them (Fig. 3H).Another remarkable aspect is that these cherts appear in different sedimentary contexts, in stratigraphic contact with other six types of facies: limestones, marls, gypsum, (calci)lutites, and sandstones (Fig. 2).

Elemental analysis
Quantitative analyses of chemical elements were performed (Tables 3 and 4).The results were plotted in the Al-Fe-Mn diagram and the two following ratios were calculated: (Fe + Mn)/Ti and Al/(Al + Fe + Mn) (Fig. 5 and Table 5).A common feature of hydrothermal deposits is an Al/(Al + Fe + Mn) ratio of < 0.4, along with a (Fe + Mn)/Ti ratio of > 15 (Boström et al., 1973;Murray, 1994;Zhang et al., 2023).Forty-two out of the fifty-six samples analyzed did not show mineral and/or geochemical indicators that could a priori be linked to hydrothermal contexts, while 14 provided preliminary results that are compatible with a hydrothermal origin according to the Al-Fe-Mn diagram and the (Fe + Mn)/Ti − Al/(Al + Fe + Mn) ratios (Boström et al., 1973;Adachi et al., 1986;Yamamoto, 1987;Murray, 1994;Yamamoto et al., 2013, Fig. 5 and Table 5).Eight of these samples (AQ-2, GA-3A, GA-3B, GA-6, GA-7, SX-A, FN1A and CN-1) correspond to cherts, three to travertines (CCH-B, TCV and TE100), two to calcilutites (PL-1 and OR-1) and one to dark lutites (BZ1-A).Except for the dark lutites of the Pliocene site of Baza-1 (Figs. 1B and 2), all the others come from materials in the Orce sub-basin of Upper Pliocene to Upper Pleistocene age.The (Fe + Mn)/Ti and Al/(Al + Fe + Mn) ratios of these samples agree with conditions of hydrothermalism for the samples of travertines and calcilutites, but not for sample the dark lutites.The last sample provided results of 13.24 and 0.42, respectively.Although these values are very close to the limit of hydrothermal conditions (> 15 and < 0.4, respectively), they are placed outside the area of hydrothermalism.These geochemical ratios could not be applied to chert samples, because neither Al nor Ti were measured in them.
Chert samples showed SiO 2 concentrations in excess of 90%, with the exception of two samples, FN1A and GA-6, which exhibited higher concentrations of other compounds (e.g., 14.23% of CaO in FN1A and 5.72% in GA-6).Their variance in composition may be attributed to the presence of calcite in the case of FN1A and gypsum in the case of GA-6, in which a concentration of 6.6% of SO 3 was measured (Table 3).Other trace elements are present in proportions < 1% by weight.The geochemical composition of the cherts exhibits an uneven distribution of trace elements (Table 4).Al was absent in all cases, and this was also the case for Mn except for samples CN1 and FN1A, in which this element was recorded with abundances of < 0.03%.Fe, Ca, and Cl were consistently detected in all samples.Moreover, Mg, S, K, and Cu were present in every sample except CN1, but the proportions of these elements were lower than in the other samples analyzed in this study.
The five additional samples linked a priori to a hydrothermal origin include travertines, dark lutites and calcilutites.
The travertine levels (samples CCH-B, TCV and TE100, Fig. 1C) are located to the west of the town of Orce.They consist of brown stratiform layers (Fig. 3E) in sedimentary unconformity with the materials deposited below them, which form patchy structures with a lateral development limited to a few meters or, in the case of samples TCV and TE100, to ∼25 m.Internal annular carbonate growths developed over previous vegetal structures are visible in these travertines (Fig. 3F).Pisolitic gravels of chaotic aspect and centimeter-scale speleothems are also common.The other travertines present in the Baza Basin show stratiform levels of decimetric thickness in the Venta Micena sector (Fig. 2), where they are stratigraphically close to the paleontological sites of this area, and there are also travertine buildings of metric thickness in the vicinity of currently active hydrothermal springs in the Zújar sector (Fig. 1B).According to the dating of similar travertine deposits in the Alicún de las Torres current hot spring, which is located in the Guadix Basin (Crespo et al., 2014;Hernández et al., 2000), the age of these travertines is Middle to Late Pleistocene.
The dark lutites that show a geochemical profile compatible with a hydrothermal origin consist of a ~ 5 cm thick level that is located at 45 cm from the top of the Baza-1 paleontological site (Figs. 2 and 3D).From a lithological point of view, these lutites show dark gray tones and contain sub-millimeter-scale parallel laminations, root remains, small fragments of gastropod shells, oxidations, and micas randomly distributed in the sediment.The presence of dark lutites in the Baza Basin is commonly recorded in thin layers present in most of the stratigraphic series of Plio-Pleistocene age (Fig. 2).Mesoscopic scale sedimentary structures of traction origin and linked to detrital deposits (e.g., sandy microlevels, coarsening/finning upward sequences, imbrications, or scour-marks) are absent.Analyses of sample BZ1-A (Tables 3 and 4) showed a main chemical composition of Si > Fe > Al > Ca > K > Mg.Compared to other dark lutites studied, BZ1-A stands out for its notably higher concentration of Fe 2 O 3 , which reaches 11.19%.This iron content exceeds those of other samples of dark lutites.Additionally, this sample exhibits elevated levels of K, Mn, S, Cl, Ba, Cu, Rb, Pb, Sn, and Cs compared to the average contents of the samples of dark lutites.
Calcilutites are present in all stratigraphic profiles except Zújar (Fig. 2) and use to appear as light-colored levels of decimetric to metric thickness.From a lithological point of view, they typically exhibit massive facies or scarce internal arrangement, although sometimes there are internal shallowing-upward sequences with edaphogenic-type facies of reddish tones situated at their top.The PL-1 calcilutite sample shows a main geochemical composition of Ca > Si > Fe > Al > Mg (Tables 3 and 4).The rest of the elements detected appeared with concentrations < 0.6%.In comparison with other calcilutites of the basin, this sample has the highest CaO values and the lowest of SiO 2 , Al 2 O 3 , K 2 O, TiO 2 , Na 2 O, and SO 3 .Finally, the OR-1 calcilutite sample shows a main geochemical composition of Si > Ca > Mn > Al > Fe > Mg.The rest of the elements detected appear with concentrations < 0.7%.This sample exhibits the highest MnO value (~ 10 times more) compared to the average of the calcilutite samples.Additionally, geochemical analyses revealed the presence of some trace elements in the samples (Table 4).Of these, Ba stands out, being exclusively detected in four samples a priori linked to hydrothermal environments: BZ1-A (dark lutites), OR-1 (calcilutites), TCV (travertines), and TE100 (travertines).

Mineralogical analysis
The XRD data of cherts show that opal is the principal phase (Table 6, Fig. 4).The d (101) (> 4.07 Å) and the full width at half-high maximum (FWHM) reflection peak intensity of opal ( 101) indicate that CT is the polymorph of the opal (Elzea et al., 1994).Most of the samples taken in Galera showed the presence of gypsum in the XRD data, as evidenced by geochemistry (high CaO and SO 3 contents, Table 3) and field observations (Fig. 3H).The presence of silicates and carbonates in some samples (e.g., CN-1, GA-3B and SX-A) may come from the remains of the sedimentary materials above and/or below the levels where these samples were obtained, due to the irregular morphology of the chert strata.Travertine samples mainly consist of calcite (Table 6) and, to a lesser extent, pyrolusite.In addition, other mineral species such as gypsum and halite were found.The XRD data from these samples agree with the expectations from geochemical (Tables 3 and 4) and lithological data.The analyses carried out on the calcilutites of samples PL-1 and OR-1 (Table 6) showed a mineral composition of calcite, quartz and phyllosilicates (muscovite and paragonite).In Fig. 3 Field images of cherts and other facies linked to hydrothermal scenarios and associated to tectonic fractures.A Deposits of light colored cherts, which are located 3 km to the NW of the Orce town (sample SX-A, Fig. 1C); B detail of white cherts located close to the FN1 paleontological site (sample FN1A, Fig. 2); C thin level of grey to blue cherts (just at the top of the hammer), located at south of the Galera town (sample GA-6, Fig. 2); note the presence of deformed layers at the bottom of this image.D dark level located at the BZ-1 Pliocene site (sample BZ-1A, Fig. 2).E travertines (sample TE100, Fig. 1C) located west of Orce.F detail of these travertine facies.G NE-SW oriented tectonic fault present at the south of Galera (marked with an orange arrow), which affects the chert levels; the fault plane appears to be filled with gypsum ore.H mineral association between gypsum crystals (marked with an orange circle) and cherts, near the locality of Galera.I ~ 3 m wide chert dyke (the limits of the dyke are marked with orange arrows) with a NE-SW orientation, which coincides with a tectonic fracture near Galera (Fig. 1C).J detail of the N105°E tectonic contact between a chert layer and a deposit of calciarenites (marked with an orange arrow), which appears close to the paleontological site of FN1 (Fig. 1C)  The results are expressed as oxides weight %.LOI column indicates Loss of Ignition of each sample (%) relative with volatile components.Facies key: cherts (white); travertines (orange); sulphur (yellow); gypsum (light yellow); dark lutites (grey); lutites (blue); calcilutites (pink); marls (violet).Red colored cells indicate the maximum amounts for this element.Green colored cells indicate the minimum amount for this element, above limit detection.bld below limit detection (< 0.01%) addition, dolomite and plagioclase were both present as trace amounts.Finally, the geochemistry of the sample of dark lutites with hydrothermal markers (BZ1-A) consisted mainly of phyllosilicates (> 50%: muscovite, paragonite, kaolinite, and chlorite) and quartz (Table 6).Accessory minerals included feldspars, calcite, dolomite, lepidocrocite, halite, and gypsum.The high FeO 3 , CaO, SO 3 , and Cl (Tables 3 and  4) contents justified the presence of lepidocrocite, gypsum and halite.

Discussion
On a cartographic scale, 14 samples linked to hydrothermal environments were found in the vicinity of traces of the main tectonic fractures present in the basin (García-Tortosa et al., 2011;Sanz de Galdeano et al., 2007, 2012; Fig. 1B, C): samples AQ-2, GA-3A, GA-3B, GA-6, GA-7, and PL-1 are close to the NE-SW Galera Fault; sample CN1 is close to the NW-SE Alfahuara-Botardo Fault; sample BZ-1A is close to the NNW-SSE Baza Fault; and samples SX-A, FN1A, OR-1, CCH-B, TCV, and TE100 are all in the vicinity of WNW-ESE oriented fractures.This suggests the connection between materials with evidence of hydrothermal activity and their position close to tectonic fractures, as evidenced by the fact that the current thermal springs in the basin are found on or in the vicinity of traces of tectonic fractures (Fig. 1B).At an outcrop scale, the relationship between chert deposits and tectonic fractures is evident (Fig. 3G, I, J).

Origin of cherts and other facies potentially generated in hydrothermal contexts
Chert layers appear in the NE sector of the Baza Basin (Orce sub-basin) and are cartographically associated with fracture traces.In fact, a chert dyke (Fig. 3I) was recognized to the SE of the town of Galera, coincident with the trace of a NE-SW fault that shows synsedimentary features (Sánchez-Roa et al., 2016).In the case of Fuente Nueva-1 and Galera, the chert deposits are found near levels rich in paleontological remains (Fig. 2).Hydrothermal deposits often show high concentrations of iron and manganese, while aluminum and titanium are reliable indicators of sedimentary input (Adachi et al., 1986;Boström et al., 1973;Murray, 1994;Zhang et al., 2023).This is why many research studies use the ratios (Fe + Mn)/ Ti and Al/(Al + Fe + Mn) to identify hydrothermal influence (Table 5).The application of the Al-Fe-Mn diagram (Adachi et al., 1986;Yamamoto, 1987;Yamamoto et al., 2013) to the chert samples (Fig. 4) showed a geochemical The results are expressed as oxides weight %.Facies key: cherts (white); travertines (orange); sulphur (yellow); gypsum (light yellow); dark lutites (grey); lutites (blue); calcilutites (pink); marls (violet).Red colored cells indicate the maximum amounts for this element.Green colored cells indicate the minimum amount for this element, above limit detection.bld below limit detection (< 0.01%)   Mineral abbreviations according to Warr ( 2021 composition in agreement with a hydrothermal origin. Although Al-Fe-Mn diagrams have been applied to deposits of marine origin (Adachi et al., 1986;Yamamoto, 1987;Boncheva et al., 2023;El-Moghazy et al., 2023;Jurkowska & Swierczewska-Gladysz, 2024; among others), they have also been applied to continental environments (Madondo et al., 2021;Zhou et al., 2022).The results obtained in the analysis of the chert samples from the Baza Basin are compatible with a hydrothermal origin for them if we consider other evidence (see below) and also according to the Al-Fe-Mn contents and main trace elements present in the siliceous rocks of the North Pacific (Adachi et al., 1986).On the opposite, the (Fe + Mn)/Ti and Al/(Al + Fe + Mn) ratios have been used in sediments to determine their hydrothermal origin (Boström et al., 1973;Murray, 1994;Zhang et al., 2023).Therefore, the data obtained from the chert samples are expected to be more reliable than those derived from samples without chert contents.In any case, the use of both methods is complementary and the comparison of their results in these continental environments allows their evaluation as proxies in this sedimentary context.The main mineral composition of opal-CT in the chert samples indicates a silica mineral phase linked to hydrothermal conditions (Flörke et al., 1990;Fröhlich, 2020;Kano & Taguchi, 1982;Lynne et al., 2005).Silica-rich deposits linked to > 250ºC hydrothermal fluids have been associated with tectonic fractures in Permian cherts based on their geochemical analyses and Al-Fe-Mn diagrams (for review and references, see Zhou et al., 2022).
According to the Na/K calculation method (Fournier, 1977;Fournier & Truesdell, 1973;Giggenbach et al., 1988), the estimates obtained with geothermometers for the high-temperature hot springs that are currently active in the Guadix-Baza Depression are 88-89 °C.A similar situation occurs in Alhama de Granada, SE Spain (López-Chicano et al., 2001).Moreover, the hot springs of the Guadix-Baza Depression (Diputación de Granada-ITGE, 1990; Table 2) show a direct relationship between their output water temperature and silica contents (García-Aguilar et al., 2014), which can reach up to 54 mg/l at 44 °C in the Graena spring.According to the standard values of the upper crustal geothermal gradient, tectonic and geophysical gravimetry data, Sanz de Galdeano et al. ( 2007) suggested a temperature range in the deep traces of the main faults of the Baza Basin of 150-300 °C, which agrees with a hydrothermal scenario.These temperatures are compatible with those inferred for the formation of Permian chert nodules in the Junggar Basin (Zhou et al., 2022).
The trace-element contents of the cherts of the Baza Basin are very variable depending on the sample considered, because each of them shows a particular geochemical profile.The case of sample FN1A, which is located a few decimeters below the Fuente Nueva-1 site, stands out for its variety of trace elements, many of them not recorded in other deposits (e.g., Sr, As, Mo, and W; Table 4).In this sense, the presence of As has been associated with hydrothermal fluids near tectonic fractures west of the Orce sub-basin (Sánchez-Roa et al., 2016) and the presence of Sr has also been linked to hydrothermal contexts in the Plio-Pleistocene lacustrine environments of the Guadix-Baza Depression (García-Aguilar et al., 2014).On the other hand, the variety of geochemical profiles found in the chert samples is compatible with the variety of hydrochemical profiles (both qualitatively and quantitatively) and output temperatures present in the hot springs that are currently active in the basin (Table 2).These hydrochemical profiles must be linked to the geometry of the tectonic fractures associated with each hot spring and to the different lithological substrates that the hydrothermal fluids passed through on their way to the surface.This cause-effect relationship could be extrapolated to the ancient Plio-Pleistocene hydrothermal upwellings in the Baza Basin, whose salt contents, output temperatures and water flow would have varied depending on the tectonic, lithological, and structural features present in each case.
The hydrochemical composition of the hot springs that are currently active in the basin with output temperatures of > 35 °C (Cruz-Sanjulián et al., 1972;Cruz-Sanjulián & García-Rosell, 1972, 1975;Diputación de Granada-ITGE, 1990;Table 2) shows silica contents of > 40 mg/l, high sulfate values and the presence of chlorides.Similar conditions are widely documented in hydrothermally active regions, such as the Dead Sea-Jordan River Basin in Israel (Levitte & Eckstein, 1978), Yellowstone in USA (Nordstrom et al.,  5.The green area indicates a hydrothermal origin (Adachi et al., 1986;Yamamoto, 1987).The samples labeled outside the green area are those that show other geochemical evidence of hydrothermalism (see Table 5) 2005), Mount Tafta in Iran (Shakeri et al., 2008) or Denizli in Turkey (Özler, 2000), among others.On the other hand, White et al. (1956) and Fournier and Rowe (1966) established the hydrochemical conditions linked to the solubility and subsequent precipitation of silica in thermal waters.They defined a solubility of 315 mg/l for silica at 90 °C and 110 mg/l at 25 °C, demonstrating a correspondence between the output temperature of the hot springs and their silica contents.The rapid rise of hot waters to the surface results in silica oversaturation due to cooling and degasification.Dissolution by high-temperature waters of deep rocks and their subsequent ascent through tectonic fractures is common.In the case of cherts, however, these hot springs should have shown much higher output temperatures (> 35 °C) than in the atmosphere.Considering the climatic conditions inferred for the Early Pleistocene in the Baza Basin, these thermal pools would have originated small areas with a stable and warm microclimate throughout the year, suitable for the ecological development of a diverse mammalian community, especially during the cold periods of the Pleistocene (Lisiecki & Raymo, 2007;Paillard, 2017).
Under these conditions, cherts originated from the dissolution at depth conditions of rocks with elevated siliceous contents (e.g., quartzites and mica-schists of the Internal Zones of the Betic mountain range; Baena Pérez et al., 1979;Martín García et al., 1980;Villalobos et al., 2006) in hightemperature fluids.According to the most plausible calculated Na/K geothermometers (88-89 °C) and considering the standard geothermal gradient in the upper crust (25-30 °C/ km, Harrison et al., 1990;Gupta & Roy, 2007), the estimated depth in the Baza Basin for the dissolution of silica agrees with the basement depth calculated from seismic profiles, 2150 m (Alfaro et al., 2008).This depth is similar to the estimation for present-day hot springs in Algeria (Belhai et al., 2017) and indicates the presence of a deep fracture network affecting the basin substrate in the Guadix-Baza Depression, where the hot springs would dissolve components of the rocks bound to them.These depth values linked to the origin of thermal upwellings in the basin are lower than the depth range of the tectonic fractures in the basin (5-10 km), which has been estimated using geological and geophysical criteria (Sanz de Galdeano et al., 2012Galdeano et al., , 2020)).
Travertine levels of decimetric thickness and associated with a hydrothermal context have been detected in the southern bank of the Orce River, west to the town of Orce.The application of the Al-Fe-Mn diagram and the (Fe + Mn)/ Ti and Al/(Al + Fe + Mn) ratios to these samples showed a composition compatible with a hydrothermal origin (Fig. 5, Table 5).
Furthermore, the Al-Fe-Mn diagram suggests that sample BZ1-A of dark lutites may have had some hydrothermal influence.Moreover, the mineralogy of this sample (Table 6) indicates a primarily detrital deposit (quartz, paragonite and chlorite) with halite precipitation due to evaporation processes.It is noteworthy that the presence of Sr and Mo was detected in this sample (Table 4), as these elements have been found under a hydrothermal context in the Galera deposits (Sánchez-Roa et al., 2016).Despite having very high iron contents, the hydrothermal ratios (Fe + Mn)/Ti and Al/(Al + Fe + Mn)) suggest that this level has not a hydrothermal origin.Given the discrepancy between the Al-Fe-Mn diagram and these ratios, it is important to note that the diagram was designed for analyzing siliceous materials (Adachi et al., 1986) while the ratios are typically used for sediments (You et al., 2020;Zhang et al., 2023).As a result, it is plausible that the ratios provide a more accurate indication of hydrothermalism in these dark lutites.However, although the values in these ratios are close to those expected in hydrothermal environments, they are not within the hydrothermal zone.This casts doubts on whether BZ1-A was influenced by hydrothermalism.In any case, it is remarkable that BZ1-A is the only sample among the 17 samples of dark lutites analyzed that fits the diagram of hydrothermalism.
The calcilutite samples that provided geochemical indicators of hydrothermal activity correspond to two metric scale levels located near the localities of Huéscar (sample PL-1) and Orce (sample OR-1, Fig. 2).In both cases, they are placed in the vicinity of paleontological sites with macrofauna.The application of the (Fe + Mn)/Ti and Al/ (Al + Fe + Mn) ratios to these samples (Table 5) shows values that are compatible with a hydrothermal origin for them.In the case of sample PL-1, the Al-Fe-Mn diagram (Fig. 4) shows further evidence of its hydrothermal origin.From the mineralogical point of view, the composition of these samples shows quartz, calcite, dolomite, muscovite, paragonite, and plagioclase (Table 6).This association of minerals can be considered detrital-inherited by sedimentary processes, as evidenced by the lithological characteristics observed in the outcrops where the samples were taken.In addition to the previous arguments, the presence of significant amounts of Ba in samples BZ1-A (dark lutites), OR-1 (calcilutites), TCV (travertines), and TE100 (travertines) adds further evidence on the hydrothermal context associated with their deposition (Jankovic, 1982;Canet et al., 2005).
The hydrothermal markers in samples Pl-1, OR-1 and BZ1-A, whose lithologies contain inherited detrital elements, could be interpreted as due to minerals originating from source areas with metamorphic substrates.However, these hydrothermal markers do not appear in other deposits that are lithologically and sedimentologically similar to these three samples, which indicates that the high Fe and Mn contents should have a hydrothermal origin.Concerning the sedimentary context of the cherts and other deposits linked to hydrothermal conditions and their association with archaeological and paleontological sites, it is worth noting that these levels frequently appear on the top and/or bottom of the sites as well as in lateral connection with facies linked to a diversity of sedimentary scenarios.

Paleogeographical reconstruction of hydrothermal activity in the basin
The spatio-temporal correlation between the coordinates of the 14 samples with evidence of past hydrothermal activity and the position of the hot springs that are currently active in the Baza Basin allows us to propose a general model of paleogeographic evolution of the basin in relation to its hydrothermal activity (Fig. 6).This model makes it possible to specify the cause-effect relationships between the tectonic dynamism of the basin and the occurrence of hydrothermal phenomena linked to tectonic fractures (Fig. 1B, C) at specific places and times that generated in each case different geochemical and mineral markers in the sedimentary environments.These environments can be very diverse, since they do not show a direct genetic relationship with the position of the hydrothermal upwellings.In this way, the cherts, travertines, calcilutites, and dark lutites (i.e., the facies in which hydrothermal markers have been detected) correspond to different sedimentary environments, which include lacustrine zones with some depth (calcilutites), marshy environments (dark lutites), occasional upwellings in favor of tectonic fractures (cherts) and waterfalls or pools of thermal waters (travertines), according with the most plausible sedimentary models (Prothero & Schwab, 2004;Reading, 1978).Among these sedimentary environments, those linked to the cherts could be considered the only ones with a direct causal relationship between these facies and the presence of hydrothermal activity.
As shown in Fig. 6, hydrothermal phenomena appear discontinuously in the Baza Basin from the late Pliocene (Ruscinian) to the present.In this sense, there seems to be a gap between the upper Ruscinian and the middle Villafranchian, together with a gap in the lower upper Villafranchian.However, the presence of hydrothermal activity during these periods cannot be ruled out.Its demonstration would depend on a wide sampling of the different sedimentary facies linked to these periods, preferably located near tectonic fractures.Regarding the geographical distribution, it is evident that hydrothermal activity concentrates in the Galera-Orce-Cañada del Negro sector (Orce sub-basin).This points to an intense tectonic activity in this sub-basin during the Plio-Pleistocene, particularly in four intervals: late Ruscinian, middle Villafranchian, lattermost Villafranchian and Galerian to present, with the presence of hydrothermal phenomena that induced mineral and/or geochemical markers in the sedimentary environments associated to these places and times.Likewise, the Baza area seems to show a similar pattern to the Orce sub-basin for the Ruscinian, associated with an important tectonic fracture, the Baza Fault.It is also evident that each node of hydrothermal activity seems to demonstrate a particular physico-chemical characteristic, as shown by the compositional variety of the deposits linked to this activity.In this regard, factors such as the hydrochemical composition and physical parameters of the hydrothermal fluids, the sedimentary environments where they acted and (micro)biological factors would give rise to a wide range of situations.This pattern could have played an important role in the development of local ecosystems in the vicinity of the nodes of hydrothermal activity.Finally, it is evident how these hydrothermal phenomena coincide quite closely in their spatio-temporal coordinates with the position of important paleontological sites, as happens in the case of BZ-1, FN1, HU-1, and BO-7 (Fig. 2).If we further extend these spatio-temporal coordinates to a wider range, we see other important sites (e.g., BL, FN-3 and VM) that appear in a geographic and/or chronostratigraphic position close to the nodes of hydrothermal activity.This allows us to establish a possible paleoecological relationship between the presence of hydrothermal phenomena and the occurrence of paleontological sites in the Baza Basin.

Paleoecological implications of hydrothermal activity in the basin
A contribution of thermal waters would play a role in maintaining during the dry season a minimum level in the water table of the lacustrine environments of the Baza Basin.This would in turn allow the development of a luxurious vegetational cover in the surroundings of the hot springs, which would sustain a large vertebrate community.The hot springs of the Baza Basin led to the appearance of thermal pools of variable extent and length, with waters that differed from the main lacustrine systems in their composition and  7).This is reflected in some of the most important deposits of the basin, such as the Pliocene sites of Baza-1 and Fuente Nueva-1, which are both linked to a level of dark lutites and to a level of cherts, respectively.From a sedimentary point of view, this provides evidence on the role of hydrothermal events as hot spots for the concentration of terrestrial fauna.
The relationship between some ecological factors linked to contexts with thermal waters is shown in Fig. 8A, which depicts a direct, statistically significant relationship between prey biomass (kg/km 2 ) and annual rainfall (mm) values in a set of African Natural Parks and Game Reserves (data compiled from Bell (1981) and Hatton et al. (2015); see also Coe et al. (1976) for a similar relationship).However, there is considerable scatter around the least-squares regression line, which primarily depends on soil type and, secondarily, on the presence of thermal springs (Fig. 8B), which increase the net primary productivity thanks to their saline contents (García-Aguilar et al., 2014, 2015;Palmqvist et al., 2022a).Figure 8 also includes the mean estimates of carrying capacity for the Orce sites (3827.8kg/km 2 for VM and 3307.6 kg/km 2 for BL/FN-3), obtained using an approach based on the Weibull model (see Rodríguez-Gómez et al., 2022, 2024), and the estimates of annual rainfall, based on the composition of the herpetological assemblages in the case of BL/FN-3 (700 mm for BL and 738 mm for FN-3;Blain et al., 2016) and on the abundance of nitrogen isotopes in the bone collagen retrieved from the fossils of large mammals from VM (800 mm; Palmqvist et al., 2022cPalmqvist et al., , 2023)).As deduced for the Baza Basin, the presence of thermal waters provided a variety of salts and trace elements.These mineral nutrients allowed the development of a diverse plant community and the translation of this biodiversity to the local trophic pyramid, thus configuring a sort of ecological "Eden" on a limited spatio-temporal scale.Given that some species of mammals, for example the hippopotamuses, used these thermal pools and watered them, their salinity (at least where the fossil remains are found) must have been of < 5 g/l (Palmqvist et al., 2022a).In this environmental context, the ancient hominin settlers of the Baza Basin sought such places as preferential sources of water and a variety of trophic resources.
Figure 8B shows that although the productivity of the African National Parks and Game Reserves mainly depends on annual precipitation, as noted above, it also relates to soil type and to the presence of thermal waters, which contribute to increase the productivity of the terrestrial ecosystems and their herbivore biomass.Specifically, the areas with volcanic and Rift sediments have more positive residuals in the regression equation of Fig. 8A than those with other types of sediments (i.e., basement, floodplain, granite, Kalahari sands, and marine sediments; see also Bell, 1981).Similarly, the areas with active hot springs show greater values of prey biomass than those in which thermal activity is absent (Fig. 8B).
The network of spatio-temporal nodes of hydrothermal phenomena in the Baza Basin seems to be related with the presence of hominin settlements, as evidenced by four archaeological localities: the late Early Pleistocene sites of Barranco León and Fuente Nueva-3, the Early-Middle Pleistocene site of Huéscar-1, and the early Middle Pleistocene site of Cúllar-Baza-1.This possible association led us to the proposal of a new idea on the distribution of paleoanthropological sites on a global scale during the Plio-Pleistocene.
A number of African and Eurasian sites can be compared with those that characterized the Plio-Pleistocene environments of the Baza Basin.Throughout the East African Rift Valley, from Mozambique-Malawi to Ethiopia and Eritrea, following the Red Sea and the Levantine Corridor to Anatolia and Georgia, the main archaeological deposits are found in volcanic contexts and scenarios characterized by intense tectonic activity, which resulted in geothermal anomalies during the Neogene-Quaternary (Adamia, 2013;Glerum et al., 2020;Morley et al., 1992;Schlüter, 1997).More specifically, a relationship between hydrothermal activity and paleoanthropological deposits of 1.7 Ma has been evidenced by microbial markers in Olduvai Gorge by Sistiaga et al. (2020), who used a similar model to the one proposed here for the Baza Basin.Such correspondence also occurs in several European sites linked to recent hydrothermal and seismo-tectonic activity, for example the Acheulean site of Rodafnidia, Lesvos Island (Galanidou et al., 2013(Galanidou et al., , 2016)).This shows the common relationship between areas with intense geothermal activity and the presence of archaeological sites with hominin remains and fossils of large mammals.
Tectonic activity is a source of geomorphological variability (i.e., landscape roughness) and ecological diversity (Coppens, 1994;Foley, 1987;Gamble, 1993) and has been linked to the first hominins that evolved in the African Rift, because it allowed them to obtain stable resources: the tectonic and volcanic movements that generated the East African Rift System resulted in a broken, hilly rough country of complex topography where the early hominins could hunt and hide (Bailey et al., 2000(Bailey et al., , 2011)).Moreover, environmental roughness probably characterized the whole route of early hominin dispersal (King & Bailey, 2006).We suggest here that tectonic and hydrothermal activity are significantly linked to the evolution of the human lineage due to the consumption strategies of hominins and their risk of exposure to confrontation with carnivores.In fact, carnivores could have predated the own hominins (Treves & Palmqvist, 2007) and were also competitors in the access to meat, although they were also providers of carrion (Espigares et al., 2013;Palmqvist et al., 2022aPalmqvist et al., , 2023)).According to the evidence unearthed from the fossil sites of Orce, the hominins developed a scavenging behavior as a successful strategy to obtain animal resources in the Baza Basin (Espigares et al., 2019;Palmqvist et al., 2023).Analyses on resource partitioning among the members of the carnivore guild suggest that scavenging was the optimal strategy for these hominin populations (Rodríguez-Gómez et al., 2016).Anatomical evidence from late Early Pleistocene sites younger than those of Orce, such as level TD6 from Gran Dolina in Sierra de Atapuerca, suggests that these humans had lesser throwing skills than modern humans, which would hinder a hunting activity based on throwing cobbles against prey, although it would not have limited other hunting strategies (Bermúdez de Castro et al., 2020; however, for an opposite interpretation see García-Martínez et al., 2021).
A heterogeneous environment would be useful for exploiting different ecological niches in the case of an opportunistic species like Homo sp., relatively weak for a direct confrontation with the large carnivores.Moreover, these mosaic habitats would also provide protection from cursor predators in the open environments and from ambushing carnivores in the closed ones, as the hominins could inhabit the boundary zones between different environments, both for protecting themselves from predators and for exploiting the resources of these environments (see Treves & Palmqvist, 2007).In ecology, the transition zones between different ecosystems are called ecotones.We propose that rather than a generalist species, the early members of Homo would be specialized in exploiting the ecotones, thus showing adaptations to consume resources from different environments.According to an analysis of avifauna in Lower, Middle and Upper Paleolithic sites (Finlayson et al., 2011), there is a significant association between hominin presence and mosaic habitats (i.e., habitats with many ecotones).
On the other hand, hominin absence of Sierra de Atapuerca during the early Middle Pleistocene has been linked to a greater environmental homogeneity than during previous periods (Rodríguez et al., 2011;Rodríguez-Gómez et al., 2014).Today, the areas with highest ecosystem diversity are two regions with active tectonics, the African Rift and the Himalayan range, in accordance with the ecoregions defined by Olson et al. (2001).The richness of ecoregions per unit area observed in the African Rift is probably lower nowadays than during the Plio-Pleistocene, because tectonic activity was higher at that time (Bailey et al., 2000;King & Bailey, 2006).This diversity of ecoregions seems to relate to the diversity in the values of abiotic and biotic variables that can be observed in the region.For these reasons, the African Rift would provide a high diversity of habitats in which a generalist or ecotone-specialist species, such as early Homo, could have efficiently exploited these environments for resources and shelter.This model can be proposed for the ancient hominin populations found in the archaeological sites of the Baza Basin (for a review on the paleodemography of these populations, see Palmqvist et al., 2022a).
From a paleoecological point of view, the Ngorongoro Crater in the African Rift Valley is a model that can be correlated with the sites of the Baza Basin and other deposits found in geologically active areas.In Ngorongoro, the hydrothermal influence of the crater lakes has a seasonal pattern.As a result, the populations of large mammals are dispersed in the wet season, while during the dry season they concentrate near the pastures that surround the hot springs.These waters lose their ecological influence at ~ 700 m from their point of emergence, changing to alkaline (pH ~ 8) waters with high saline contents (DeoCampo, 2002).The list of species that live in the crater within a radius of ~ 100 m around the thermal springs includes buffaloes, gazelles, hippopotamuses, rhinoceroses, zebras, elephants, antelopes, wild boars, hyenas, foxes, cheetahs, jackals, and lions.This fauna is strikingly similar to the one that inhabited the Baza Basin during the Plio-Pleistocene in terms of their distribution among modes of locomotion, feeding types and body size classes (Mendoza et al., 2005).

Conclusions
Hot springs have been important for faunal and hominin survival during the Plio-Pleistocene, favoring the dispersal of mammalian species as documented in the Baza Basin, which preserves one of the continental sequences with greater stratigraphic completeness of Europe, including the earliest archaeological evidence of hominin presence in the northern Mediterranean latitudes during the late Early Pleistocene.In this basin, hot springs form a network of spatiotemporal nodes, which are particularly concentrated in the Orce sub-basin, showing permanent siliceous waters with temperatures > 35 °C.These hot springs played an important ecological role during the summer, when seasonal rainfed springs became dry, and also helped to maintain locally warm temperatures in winter, thus aiding the survival of water dependent vertebrates adapted to warm and moist climates, including the members of early Homo.
Nearly one hundred Plio-Pleistocene archaeological and paleontological localities are recorded in the perimeter of the Baza Basin.Most of these sites were generated in connection with tectono-sedimentary processes, as evidenced by the presence of paleo-seismites, ancient hot springs and faulting.The lithological, mineralogical and geochemical composition of the chert levels present in the Orce sub-basin suggests a high-temperature hydrothermal origin for these deposits.Other types of sedimentary facies present in the basin, such as travertines, dark lutites or calcilutites, were influenced by this network of spatio-temporal nodes of hot springs, showing diverse geochemical and mineralogical indicators associated with this hydrothermal origin.The mineral nutrients provided by the hot springs prompted the development of a luxurious vegetation, which in turn allowed the development of a diverse community of large mammals, including the first hominins that dispersed in Western Europe.

Fig. 2
Fig. 2 Stratigraphic profiles of the Baza Basin.The codes with the position of these series are shown in Fig. 1B, C

Fig. 5
Fig. 5 Al-Fe-Mn diagram applied to all samples analyzed in this study.The points represented correspond to the values in Table5.The green area indicates a hydrothermal origin(Adachi et al., 1986;Yamamoto, 1987).The samples labeled outside the green area are those that show other geochemical evidence of hydrothermalism (see Table5)

Fig. 6
Fig. 6 Spatio-temporal position of the samples with a hydrothermal origin and of those hot springs currently active in the Baza Basin.The X-axis indicates a SW-NE line along the Baza Basin, between the localities of Zújar and Cañada del Negro (NE of Orce), with a length of ~ 50 km.Key for localities: ZU-Zújar, BZ-Baza, BN-Benamaurel, GA-Galera, HU-Huéscar, OR-Orce, CN-Cañada del

Fig. 7
Fig. 7 Paleoecological reconstruction of the landscape in the vicinity of a lacustrine system of the Baza Basin linked to a thermal upwelling.A Highstand level in the lacustrine system.B Lowstand level in the lacustrine system.The trace of the fault positioned to the left indicates the position of a local hot spring.Drawing by Sandra Ramírez Cherbuy-InkScience

Table 2
Hot springs currently active in the Guadix-Baza Depression

Table 3̄
Quantitative elemental analysis for major elements using XRF and means and standard deviations for each lithology

Table 4
Quantitative elemental analysis for minor elements using XRF and means and standard deviations for each lithology

Table 5
Geochemical ratios of (Fe + Mn)/Ti and Al/(Al + Fe + Mn) of analyzed samples

Table 6
X-ray mineralogical data of the samples taken in the stratigraphic profiles of the Baza Basin with a marked hydrothermal context (see Fig.4): xxx and dark blue cells: dominant (> 50%); xx and (Ms, Pg, Chl?)