The Lowermost Tejo River Terrace at Foz do Enxarrique, Portugal: A Palaeoenvironmental Archive from c. 60–35 ka and Its Implications for the Last Neanderthals in Westernmost Iberia

: Reconstruction of Pleistocene environments and processes in the sensitive geographical location of westernmost Iberia, facing the North Atlantic Ocean, is crucial for understanding impacts on early human communities. We provide a characterization of the lowest terrace (T6) of the Lower Tejo River, at Vila Velha de R ó d ã o (eastern central Portugal). This terrace comprises a lower gravel bed and an upper division consisting of ﬁne to very ﬁne sands and coarse silts. We have used a multidisciplinary and poor preservation of phytoliths. We provide updated OSL ages for the lower two Tejo terraces, obtained by post infra-red stimulated luminescence: (i) T5 is c. 140 to 70 ka; (ii) T6 is c. 60 to 35 ka. The single archaeological and fossiliferous level located at the base of the T6 upper division, recording the last regional occurrence of megafauna (elephant and rhinoceros) and Mousterian artefacts, is now dated at 44 ± 3 ka. With reference to the arrival of Neanderthals in the region, probably by way of the Tejo valley (from central Iberia), new dating suggests a probable age of 200–170 ka for the earliest Mousterian industry located in the topmost deposits of T4.


Introduction
The Vila Velha do Ródão and Arneiro depressions are located in the furthest upstream reach of the Lower Tejo River, about 20 km from the Spanish border (this is Portuguese Reach I as defined by reference [1]) ( Figure 1). From the sensitive location of the study area close to the North Atlantic, we can assume strong interaction between marine and terrestrial processes and environmental conditions. This highlights the relevance of this area for research on past climate and environmental change. Previous research in the study area has mainly focused on characterizing the geological setting, geomorphic genesis (e.g., of the terrace staircase), active tectonics, archaeological background and landscape evolution in general (e.g., see review in reference [2]). Climatic and environmental changes in the Lower Tejo Basin during the last glacial cycle and especially during Marine Isotope Stage (MIS) 3, remain poorly documented.
In this contribution, we focus on the sedimentary record from the lowest terrace (T6) of the Tejo at Vila Velha de Rodão, which represents an important terrestrial archive of relevance to the possible relation between palaeoenvironmental conditions and early human occupation dynamics in the region. Herein, we reconstruct and discuss the climate and environmental conditions during the last glacial cycle, through integration of evidence from optically stimulated luminescence (OSL) dating, grain-size distribution, rock magnetic properties (low field magnetic susceptibility, frequency-dependent magnetic susceptibility and isothermal remanent magnetization curves), sediment mineralogy, phytolites and palynology, as well as reviewing published paleontological and archaeological data.

Geological and Geomorphological Setting
The Tejo is one of the largest systems of Western Europe and flows E-W across almost the whole of Iberia; it is an ancient river (c. 3.7 Ma) with an important sedimentary record [23,24]. In the uppermost Portuguese reach of the Lower Tejo, the river flows through two quartzite ridges by way of the Ródão gorge (named "Portas de Ródão"), which separates the Ródão (upstream) and Arneiro (downstream) depressions.
The oldest bedrock comprises the Neoproterozoic and lower Cambrian schists and metagreywackes of the Beiras Group and the Ordovician Armorican Quartzite Formation. The latter is dominated by resistant ridges that topographically dominate (by c. 150 m) the extensive adjacent planation surface developed on phylites and metagreywackes. The Cenozoic is represented by the Cabeço do Infante Formation, the Silveirinha dos Figos Formation and the Murracha Group. The first two are dominated by soft sandstones and gravels, while the Murracha Group consists of gravels interbedded with fine sediments [23,25].
In Lower Tejo Reach I, below a culminant sedimentary unit (the Falagueira Formation, at c. +260 m-above the river bed) corresponding to the ancestral Tejo River before drainage network entrenchment, the Pleistocene to Holocene record is summarized as follows [1,2,26] (Figure 1): (i) T1, with the surface at +111 m, without artefacts; (ii) T2, at +83 m, without artefacts; (iii) T3, at +61 m, without artefacts; (iv) T4, at +34 m, with Acheulean in the basal and middle levels and Mousterian in the uppermost levels; (v) T5, at +18 m, with Mousterian industries through the entire fluvial sequence; (vi) T6, at +10 m, with Mousterian industries at the lower deposits; (vii) Carregueira Formation (aeolian sands) (32 to 12 ka), with Upper Palaeolithic to Epi-Palaeolithic industries; (viii) Alluvial plain and a cover unit of aeolian sands (Holocene), with Mesolithic and more recent industries. Immediately upstream of the Ródão gorge, the modern river bed is at c. 72 m above sea level (a.s.l.). In this area, geomorphological evidence for late Cenozoic tectonics arises from interpretation of valley asymetry and drainage patterns, fault scarps, tectonic lineaments, fracture-controlled valleys, and vertical displacement of planation surfaces and terraces [26]. The sedimentary controls for the formation of the Lower Tejo terrace staircase (mainly glacio-eustasy and differential uplift) are different from those affecting the Middle and Upper Tejo, because of separation (between the Middle and Lower Tejo) by a lengthy knick zone through hard basement [27,28]. Quaternary 2019, 3 3 of 29 is dominated by resistant ridges that topographically dominate (by c. 150 m) the extensive adjacent planation surface developed on phylites and metagreywackes. The Cenozoic is represented by the Cabeço do Infante Formation, the Silveirinha dos Figos Formation and the Murracha Group. The first two are dominated by soft sandstones and gravels, while the Murracha Group consists of gravels interbedded with fine sediments [23,25]. In Lower Tejo Reach I, below a culminant sedimentary unit (the Falagueira Formation, at c. +260 m-above the river bed) corresponding to the ancestral Tejo River before drainage network entrenchment, the Pleistocene to Holocene record is summarized as follows [1,2,26] (Figure 1): (i) T1, with the surface at +111 m, without artefacts; (ii) T2, at +83 m, without artefacts; (iii) T3, at +61 m, without artefacts; (iv) T4, at +34 m, with Acheulean in the basal and middle levels and Mousterian in the uppermost levels; (v) T5, at +18 m, with Mousterian industries through the entire fluvial sequence; (vi) T6, at +10 m, with Mousterian industries at the lower deposits; (vii) Carregueira Formation (aeolian sands) (32 to 12 ka), with Upper Palaeolithic to Epi-Palaeolithic industries; (viii) Alluvial plain and a cover unit of aeolian sands (Holocene), with Mesolithic and more recent industries. Immediately upstream of the Ródão gorge, the modern river bed is at c. 72 m above sea level (a.s.l.). In this area, geomorphological evidence for late Cenozoic tectonics arises from interpretation of valley asymetry and drainage patterns, fault scarps, tectonic lineaments, fracture-controlled valleys, and vertical displacement of planation surfaces and terraces [26]. The sedimentary controls for the formation of the Lower Tejo terrace staircase (mainly glacio-eustasy and differential uplift) are different from those affecting the Middle and Upper Tejo, because of separation (between the Middle and Lower Tejo) by a lengthy knick zone through hard basement [27,28].   T1; 3-T1;  4-T2; 5-T3; 6-erosion level correlative with T3; 7-T4; 8-erosion level correlative with T4; 9-T5; 10-erosion level correlative with T5; 11-T6; 12-alluvial plain; 13-colluvium; 14 Previous dating of the Pleistocene litostratigraphic divisions was undertaken using Uranium-series, Thermoluminescence (TL), Optically stimulated luminescence on quartz (Quartz-OSL) and infra-red stimulated luminescence (IRSL) on K-feldspar. Quartz-OSL was used to date: the Carregueira Formation (32 to 12 ka; 3 ka) [2]; the alluvial deposits at the Azinhal archaeological site, linked to T6 (61 ± 2 ka: GLL code 050302); the topmost deposits of T4 at the Pegos do Tejo-2 archaeological site (minimum age of 135 ± 21 ka: GLL code 050301) [29]. In Reach I of the Tejo, IRSL dating provided a first temporal framework for the Lower Tejo terraces [26], but at that stage of research T5 and T6 were not yet separated as distinct terraces. Recently, T4 in Lower Tejo Reach IV was dated to c. 340 to 155 ka, using post infra-red stimulated luminescence (pIRIR) [28].

Methods
Geomorphological, stratigraphical, sedimentological and chronological data were obtained using standard methodology (e.g., [30]): (1) geomorphological study, complemented by local detailed investigations and the production of a detailed map using Geographic Information System (GIS), (2) field descriptions of the sedimentary units, (3) sedimentological characterization of the deposits and (4) luminescence dating.

Field Work
The T6 deposits at Foz do Enxarrique were studied in detail in order to improve our understanding of the local stratigraphy and sedimentology. Exposures of T5 and T4 in the study area were also revisited. Fieldwork included stratigraphic logging and sedimentological characterization of the sedimentary deposits in order to obtain data on the depositional facies, including sediment colour, texture, maximum particle size, clast lithology, fossil content, bedding and depositional architecture.
At the studied stratigraphic section of T6 at Foz do Enxarrique, continuous sediment sampling was undertaken manually, every 1 cm, to a depth of 5.00 m. Samples were labelled as follows: e.g., for "T6FE0.21", T6 identifies the terrace code, FE the site and 0.21 the sample depth. We also collected a present-day sediment sample from the Foz do Enxarrique stream bed (FE-modern). Each sampled horizon (1 cm) was characterised according to its colour (using the Munsell system), texture and the relative abundance of carbonate concretions. Phytolith analyses were undertaken on carbonate concretions from five levels within the sequence (T6FE0. 72 From T4, samples for OSL dating were obtained from the Vilas Ruivas, Rodense Bolaria (Vila Velha de Ródão) and Pegos do Tejo-2 (Arneiro) sites ( Figure 1).

Optically Stimulated Luminescence Dating
From Reach I of the Lower Tejo, samples for OSL dating were previously collected from the sedimentary successions of the T4, T5 and T6 terraces and dated by IRSL, with a correction for the anomalous fading effect [26], because the quartz-OSL signal was found to be in saturation. However, it was later documented that the fading correction used was inappropriate, leading to age underestimation. In order to improve the chronology of the terrace sequences and their associated lithic industries, we now use the pIRIR protocol. K-feldspar grains from the samples dated in 2008 (IRSL) were measured (Equivalent doses) in 2013 by pIRIR (the most up-to-date protocol); these new results are now presented. The storage of K-feldspar grains does not affect the luminescence properties or the resulting ages. In summary, several samples were selected for OSL dating (pIRIR protocol): three from T6 at Foz do Enxarrique (upper unit), five from T5 (collected at Foz do Enxarrique and Vilas Ruivas; Figure 1) and four from T4 (collected at Pegos do Tejo-2, Vilas Ruivas and Rodense Bolaria/V.V.Ródão; Figure 1).
OSL is an absolute dating technique that measures the time elapsed since sedimentary quartz or feldspar grains were last exposed to daylight [31]. Exposure to daylight during sediment transport removes the latent luminescence signal from those minerals. After burial, the luminescence signal (trapped charge) starts to accumulate in the mineral grains due to ionising radiation. The annual dose of a sediment sample is related to the decay of 238 U, 232 Th and 40 K present in the sediment itself, to cosmic ray bombardment and to the water content of the sediment. In the laboratory, the equivalent dose (D e , assumed to be the dose absorbed since the last exposure to light, i.e., the burial dose, expressed in Grays -Gy) is determined by comparing the natural luminescence signal resulting from charge trapped during burial with that trapped during a laboratory irradiation. In this study, the radionuclide concentrations were measured by high-resolution gamma spectrometry [32]. These concentrations were then converted to environmental dose rates using the specified conversion factors [33]. For the calculation of the dose rate of sand-sized K-feldspar grains, an internal K content of 12.5 ± 0.5% was assumed [34]. Dividing the D e by the environmental dose rate (in Gy/ka) gives the luminescence age of the sediment.
Sample preparation for luminescence analyses was carried out in darkroom conditions, at the Department of Earth Sciences of the University of Coimbra. Samples were wet-sieved to separate the 180-250 µm grain-size fraction, followed by HCl (10%) and H 2 O 2 (10%) treatments to remove carbonates and organic matter, respectively. The K-feldspar-rich fraction was floated off using a heavy liquid solution of sodium polytungstate (ρ = 2.58 g/cm 3 ). The K-feldspar fraction was treated with 10% HF for 40 min to remove the outer alpha-irradiated layer and to clean the grains. After etching, the fraction was treated with HCl (10%) to dissolve any remaining fluorides.
At the Nordic Laboratory for Luminescence Dating (NLL), OSL were conducted using automated luminescence Risø TL/OSL-20 readers (Roskilde, Denmark), each containing a calibrated beta source. Small (2 mm) aliquots of K-feldspar were mounted in stainless steel cups. The K-feldspar equivalent doses (D e ) were measured with a pIRIR SAR protocol using a blue filter combination [35,36]. Preheating was at 320 • C for 60 s and the cut-heat 310 • C for 60 s. After preheating the aliquots were IR bleached at 50 • C for 200 s (IR 50 signal) and subsequently stimulated again with IR at 290 • C for 200 s (pIRIR 290 signal). It has been shown [36] that the post-IR IRSL signal measured at 290 • C can give accurate results without the need to correct for signal instability. For all IR 50 and pIRIR 290 calculations, the initial 2 s of the luminescence decay curve less a background derived from the last 50 s was used.

Grain-Size Measurements
Grain-size analyses of uncemented sediment samples was carried out using a Beckman Coulter LS230 laser granulometer (Brea, CA, USA), with a measurement range of 0.04 to 2000 µm and a Quaternary 2019, 2, 3 6 of 30 relative error less than 2%. Visual inspection of grain-size distribution curves allowed the identification and interpretation of unimodal or multimodal subpopulations. The T6FE sediment samples of the 5.00-3.20 m depth interval were analysed at a 5 cm spacing; the 3.18-2.60 m and 2.50-0.30 m depth intervals were analysed at a 1 cm spacing in order to provide a better distinction between fluvial and aeolian deposition.

Mineral Composition
Analyses of sediment composition were based on binocular microscope observation and X-ray diffraction (Department of Earth Sciences-University of Coimbra), as well as Scanning Electron Microscopy (SEM) and Energy Dispersive Spectometry (EDS) of selected carbonate concretions (UNESP Laboratory, at the University of Rio Claro-Brazil). A Philips PW 3710 X-ray diffractometer (Virginia, USA), with a Cu tube, at 40 kV and 20 nA was used for mineralogical identification within carbonate concretions and for clay mineralogy. The mineralogical composition of the <2 µm fraction was obtained in oriented samples before and after ethylene glycol treatment and heating up to 550 • C. The percentages of the clay minerals in each sample were determined through the peak areas of the mineral present, with the use of specific correction parameters.

Rock Magnetism
For magnetic susceptibility measurement, samples were dried at 40 • C and transferred into plastic bags for subsequent analysis. Rock magnetic properties were measured in the Instituto Dom Luis, University of Lisboa and in the Department of Earth Sciences of the University of Coimbra, and consisted of low field mass specific magnetic susceptibility (χ in m 3 /kg), frequency-dependent magnetic susceptibility (Kfd in %) and isothermal remnant magnetization (IRM). Magnetic susceptibility measures the ability of a material to be magnetized and includes contributions (in proportion to their abundance) from all diamagnetic (calcite), paramagnetic (clays), and ferromagnetic (magnetite) minerals present in the sediment. Low-field magnetic susceptibility was measured with a MFK1 (AGICO Inc, Brno, Czech Republic) apparatus operating with magnetic field intensity of 200 A/m and frequency of 978 Hz. Data were reported as mass-normalized values (m 3 /kg). Frequency-dependent magnetic susceptibility is an indicator of the presence of superparamagnetic particles (SP), generally produced during pedogenic processes. Low (0.47 kHz) and high (4.7 kHZ) frequency-dependent magnetic susceptibility was measured with a Bartington Instruments magnetic susceptibility meter coupled to a MS2B sensor and reported in percentage as follows: Kfd (%) = 100*(Klf-Khf)/Klf. After cleaning by alternating field demagnetization up to 100 mT, samples were subsequently submitted to stepwise isothermal remanent magnetization (IRM) acquisition with an impulse magnetizer (model IM-10-30). We applied maximum fields of 1.2T following approximately 30 steps. Remanence was measured with a JR-6A (AGICO Inc, Brno, Czech Republic) magnetometer.
Data were analysed using a cumulative log-Gaussian (CLG) function with software developed for the purpose [37]. The S-ratio was calculated with the formula −IRM −0.3T /IRM 1T .

Phytoliths
For phytolith analyses, samples with a volume of 1 cm 3 were placed in an Erlenmeyer flask and dissolved in 20 mL of HNO 3 and H 2 SO 4 solution at 1:4. The material was heated for 3 h at 90 • C on a hot plate. After cooling at~25 • C, 10 ml of H 2 O 2 was added, before washing in distilled water, centrifuging (1500 rotations per minute up to neutralization (pH~7.0), and washing with alcohol. For slide preparation, 50 µL of material was extracted by pipette, placed on slides and dried on hot plates. Coverslips were fixed using Entelan ® resin (Hatfield, UK). Phytoliths were analyzed through optical microscopy (×160 and ×640), identified with reference to literature [38][39][40] and named according to the International Code of Phytolith Nomenclature [41].

Palynology
For palynological studies, thirteen sediment samples were selected: two from lower and middle levels within T5; three from the T6 archaeological level and eight samples from the upper division of T6 (silty very fine sands and sandy silts). These samples were subjected to a physical and chemical pollen concentration pre-treatment. The pollen residue was isolated with a standard palynological preparation methodology [42], with some modifications: omitting acetolysis and sieving, the latter in an attempt to increase pollen concentration. The pollen residue was assembled on thin glass slides to allow its identification and counting. It was embedded in glycerine and sealed with histolaque, to permit movement of the grains for more complete observation of the morphological features of pollen and non-pollen-palynomorphs. Grains were identified (based on references [43][44][45][46]) and counted using an optical transmitted-light microscope.

Geochemical Analyses
Geochemical analyses on sediment samples collected from the upper unit of T6 were performed, at the laboratory, with a X-ray fluorescence spectrometer (Niton XL3t Ultra Analyser-Thermo Fisher Scientific; Waltham, MA, USA).

Geomorphology, Lithostratigraphy and Sedimentology
In the confluence area of the Enxarrique stream with the Tejo, several geomorphic units are T6 (silty very fine sands and sandy silts). These samples were subjected to a physical and chemical pollen concentration pre-treatment. The pollen residue was isolated with a standard palynological preparation methodology [42], with some modifications: omitting acetolysis and sieving, the latter in an attempt to increase pollen concentration. The pollen residue was assembled on thin glass slides to allow its identification and counting. It was embedded in glycerine and sealed with histolaque, to permit movement of the grains for more complete observation of the morphological features of pollen and non-pollen-palynomorphs. Grains were identified (based on references [43][44][45][46]) and counted using an optical transmitted-light microscope.

Geochemical Analyses
Geochemical analyses on sediment samples collected from the upper unit of T6 were performed, at the laboratory, with a X-ray fluorescence spectrometer (Niton XL3t Ultra Analyser -Thermo Fisher Scientific; Waltham, MA, USA).

Geomorphology, Lithostratigraphy and Sedimentology
In the confluence area of the Enxarrique stream with the Tejo, several geomorphic units are  In the Vila Velha de Ródão -Arneiro area, T4 comprises a lower boulder gravel (up to 2 m-thick) and an upper division of gravelly coarse sands (up to 2 m thick). The main exposures are at Pegos do Tejo-2 and Vilas Ruivas, both located in the Arneiro depression ( Figure 1). The massive clast- Tejo-2 and Vilas Ruivas, both located in the Arneiro depression ( Figure 1). The massive clast-supported boulder gravel comprises clasts that are subrounded, with MPS (mean diameter of the 10 largest clasts) = 32 cm, of quartzite (75%), white quartz and rare slates/metagreywackes. The T5 terrace usually has a basal gravelly pavement (c. 0.2 m thick) overlain by very fine sands (up to 2.3 m thick), with some very thin (<1 cm) levels containing calcium carbonate concretions. The main outcrops are provided by the sites at Foz do Enxarrique (up to c. 2.5 m thick) and Vilas Ruivas (0.5 m thick).
The sedimentary deposits of T6 are well exposed at the Foz do Enxarrique site (Figures 2-4) and were studied in detail during the work reported here. The sedimentary infill of T6 comprises a lower gravel, up to 0.40 m thick, and an upper unit, c. 5.60 m thick, dominated by very fine sands and coarse silts.  The lower artificial pavement is placed approximately at the top of the bed containing artefacts and fossil bones (the archaeological level); the lower gravel bed is located just below this. A metal staircase and a platform has been built in front of the exposed section and provides access to the top as well some protection for the exposure. Note one of the several explanatory panels provided at the site. The Tejo River can be seen in the right background.
The upper unit of T6, generally lacking bedding but with rare lamination, can be divided into three layers: (i) from the surface to a depth of 4.55 m, an upper bed comprising sandy silts; (ii) from 4.55 to c. 5.40 m, a middle division consisting of micaceous very fine to fine sands with some interbedded thin gravel stringers; (iii) from c. 5.40-5.60 m, a bed comprising Mousterian artefacts and fossil bones in a matrix of micaceous fine sands. The predominant colour of the T6 upper unit deposits is yellowish brown to bright brown. The thin (1-3 cm thick) levels with calcium carbonate concretions and rhizoliths ( Figure 5), 1-2 cm wide, are intercalated between thick intervals of uncemented sediment. At the depths of 4.55-4.57 and 4.73-4.82 m, quartzite and quartz clasts (<0.5 cm and <10 cm in size, respectively) and rolled concretions were observed dispersed in micaceous very fine to fine sands. A massive calcium carbonate level is present at a depth of 5.33-5.36 m, just above laminated fine sand. The levels containing calcium carbonate concretions do not show evidence of any erosive surface and they dip (up to 5 degrees) toward the Tejo River, progressively increasing in thickness.
The lower bed T6 consists of clast-supported boulder-pebble gravels, c. 0.4 m thick. The maximum clast diameter is 40 cm, with MPS = 31 cm. The gravel clast lithologies are quartzite (58%), milky quartz (27%) and metagreywakes/phylites (15%). The clasts are sub-rounded to angular. At about 76 m a.s.l., the T6 deposits overlie metamorphic basement, by an unconformity. The lower artificial pavement is placed approximately at the top of the bed containing artefacts and fossil bones (the archaeological level); the lower gravel bed is located just below this. A metal staircase and a platform has been built in front of the exposed section and provides access to the top as well some protection for the exposure. Note one of the several explanatory panels provided at the site. The Tejo River can be seen in the right background.
The upper unit of T6, generally lacking bedding but with rare lamination, can be divided into three layers: (i) from the surface to a depth of 4.55 m, an upper bed comprising sandy silts; (ii) from 4.55 to c. 5.40 m, a middle division consisting of micaceous very fine to fine sands with some interbedded thin gravel stringers; (iii) from c. 5.40-5.60 m, a bed comprising Mousterian artefacts and fossil bones in a matrix of micaceous fine sands. The predominant colour of the T6 upper unit deposits is yellowish brown to bright brown. The thin (1-3 cm thick) levels with calcium carbonate concretions and rhizoliths ( Figure 5), 1-2 cm wide, are intercalated between thick intervals of uncemented sediment. At the depths of 4.55-4.57 and 4.73-4.82 m, quartzite and quartz clasts (<0.5 cm and <10 cm in size, respectively) and rolled concretions were observed dispersed in micaceous very fine to fine sands. A massive calcium carbonate level is present at a depth of 5.33-5.36 m, just above laminated fine sand. The levels containing calcium carbonate concretions do not show evidence of any erosive surface and they dip (up to 5 degrees) toward the Tejo River, progressively increasing in thickness.

Luminescence Dating Results
The luminescence dating results of the samples collected from T6, T5 and T4 are summarized in Tables 1 and 2. The geographical and stratigraphic locations of each sample are provided by Table 2. Four out of the twelve samples have D e values lying above 2xD 0 (corresponding to 86% of luminescence saturation). In view of the possibility of small systematic errors in the measurement of the luminescence signal, we feel it prudent to adopt the strategy of only presenting minimum doses for these samples (equivalent to 2xD 0 ) and the derived minimum ages [47]. Two of the four samples collected from T4 were collected from the lower gravel (052208 and 052246) and provided age estimates of c. 280 ka by IRSL and minimum ages of >220 ka by pIRIR. Samples 062202 and 062203 were collected from the base and middle of the T4 upper division which mainly consists of soft sandstones, and gave underestimated ages of 209 ± 11 ka and 129 ± 11 ka by IRSL and minimum ages of >190 ka and >160 ka by pIRIR, respectively. Sample 062203 was collected below the Middle Palaeolithic site of Pegos do Tejo-2, which is located c. 1 m bellow the T4 surface. As the uppermost deposits of T4 were recently dated at 155 ka ( [28] (Cunha et al., 2017)) and sample 062203 provided a minimum age of 190 ka by pIRIR, it is very probable that the Middle Palaeolithic in the region had started by 200-170 ka.
From T5, which is onlaped by T6 at the Foz do Enxarrique site, ages of 145 ± 10 ka (at the terrace base) and 135 ± 9 ka (1.5 m below the terrace surface) were obtained by pIRIR dating. Samples collected from the T5 sequence at Vilas Ruivas were also dated, giving ages of 117 ± 7 ka, 99 ± 7 ka and 71 ± 4 ka by pIRIR.
From the T6 upper unit at Foz do Enxarrique, pIRIR dating provided the following three ages: 44 ± 3 ka, sample 052201 collected from the archaeological and fossiliferous bed (at a depth of 5.50-5.40 m); 43 ± 4 ka, sample 052202 collected from the micaceous fine sands located 20 cm higher (5.30-5.20 m); 37 ± 2 ka, sample 062201 collected c. 90 cm below the terrace surface (0.93-0.89) m. Previously, the archaeological bed was dated by Uranium series on equid and bovid teeth, providing an average age of 33.6 ± 5 ka ( [48] (Raposo, 1995)), and by fading-corrected IRSL, giving an age of 38.5 ± 1.6 ka ( [26] (Cunha et al., 2008)). It appears that both the U-series and conventional IRSL ages are underestimates, compared with the pIRIR ages. As previously noted, the IRSL dating used a correction for the anomalous fading effect which is not fully appropriate, leading to age underestimation. It should also be noted that U-series on teeth can easily be impaired as a result of the uptake of young uranium. The pIRIR protocol constitutes the best approach to date the type of material in consideration and should be considered to provide the best estimates for burial ages and for dating the associated lithic industries.

Mineral Composition
The sand and silt fractions of the T6 upper unit mainly comprise detrital grains of quartz (predominant), K-feldspar, plagioclase and muscovite, locally cemented by calcite. SEM and EDS analysis of 17 selected carbonate concretions also confirmed that these consist of calcite (micrite) developed in detrital silt-fine sand consisting of quartz (predominant), microcline, plagioclase and mica grains.
Exoscopy (SEM) (Figures 6 and 7) document well preserved microcline and quartz grains. The quartz grains are sub-angular and show conchoidal fractures, but with no evidences of abrasion or dissolution. Biotite grains show some dissolution features.
analysis of 17 selected carbonate concretions also confirmed that these consist of calcite (micrite) developed in detrital silt-fine sand consisting of quartz (predominant), microcline, plagioclase and mica grains.
Exoscopy (SEM) (Figures 6 and 7) document well preserved microcline and quartz grains. The quartz grains are sub-angular and show conchoidal fractures, but with no evidences of abrasion or dissolution. Biotite grains show some dissolution features.  The mineralogical composition of the <2 μm fraction (Table 3) of 8 samples collected from the T6 upper unit does not indicate any vertical change and consists of smectite 35-52%, illite 22-35%, kaolinite 20-31% and vermiculite 17% (only present in sample T6FE 3.46-3.50). For comparison, two samples from T5 were also studied and contained smectite 53-54%, illite 21-29% and kaolinite 18-25%. Thus both terraces have a similar clay-mineral association, although T6 has less smectite. The modern sample also has vermiculite and chlorite.  The mineralogical composition of the <2 µm fraction (Table 3) of 8 samples collected from the T6 upper unit does not indicate any vertical change and consists of smectite 35-52%, illite 22-35%, kaolinite 20-31% and vermiculite 17% (only present in sample T6FE 3.46-3.50). For comparison, two samples from T5 were also studied and contained smectite 53-54%, illite 21-29% and kaolinite 18-25%. Thus both terraces have a similar clay-mineral association, although T6 has less smectite. The modern sample also has vermiculite and chlorite.

Grain-Size Analysis
The sediment samples collected from the upper unit of T6 and analyzed by laser granulometry provided the results presented in Figure 4 and Supplementary Material. This unit is clearly dominated by coarse silt, but sampling at the 1 cm scale has documented significant grain-size oscillations as well as major grain size cycles.
-For the 4.33-5.00 m depth, the sediment is increasingly coarser, with a mean grain-size between 32 and 110 µm. In average, the sediment consists of 57% sand, 36% silt and 7% clay. Between depths of 4.73 and 4.82 m, rolled calcium carbonate concretions occur. Quaternary 2019, 3 14 of 29

Grain-Size Analysis
The sediment samples collected from the upper unit of T6 and analyzed by laser granulometry provided the results presented in Figure 4 and Supplementary Material. This unit is clearly dominated by coarse silt, but sampling at the 1 cm scale has documented significant grain-size oscillations as well as major grain size cycles.
-For the 4.33-5.00 m depth, the sediment is increasingly coarser, with a mean grain-size between 32 and 110 μm. In average, the sediment consists of 57% sand, 36% silt and 7% clay. Between depths of 4.73 and 4.82 m, rolled calcium carbonate concretions occur.

Magnetic Susceptibility
Mass specific magnetic susceptibility of samples from the Foz do Enxarrique T6 profile varies from c. 8 × 10 −8 m 3 /kg to c. 6 × 10 −7 m 3 /kg (Table S1- Supplementary Material; Figures 4, 9  and 10). These values are comparable with those of suspended sediments from other fluvial systems, such as the Jackmoor Brook catchment in South West England [49]. The lowest values of magnetic susceptibility (c. 7 × 10 −8 m 3 /kg) are found at the base of the profile, from c. 5.00 up to 4.58 m, and corresponding to sand-rich sediment (Figure 4). Magnetic susceptibility is primary controlled by lithological changes, but can also reflect pedogenic and/or post-depositional weathering driven by environmental and climatic forcing in the case of continental sediments [50]. Quartz is diamagnetic (negative and very low magnetic susceptibility), explaining the very low values of magnetic susceptibility observed at the base of the profile (Figure 4). Phyllosilicates, including clays, are paramagnetic (positive and weak magnetic susceptibility), giving significant higher values of magnetic susceptibility in the silt interval from 4.32 m to 0.40 m. The short-term cyclical variations observed in this part of the profile do not seem to be controlled by the lithology, which is very homogenous in the field, but may reflect changes in the magnetic mineralogy (ex. proportion of magnetite and hematite), including magnetic enhancement by weathering and/or pedogenic processes. The saw-tooth oscillation observed from 4.57 m to 4.32 m probably reflects the mixture of gravel with different origin and composition. In contrast, the increase of magnetic susceptibility in the topsoil horizon, which corresponds in the field to a black sediment, rich in organic material, probably reflects magnetic enhancement resulting from burning, pedogenic processes and/or anthropogenic pollution [51][52][53].

Frequency-Dependence Magnetic Susceptibility
Frequency-dependence magnetic susceptibility varies between 4 and 10 along the entire profile, indicating a significant contribution from superparamagnetic particles. These superparamagnetic particles are generally ultra-fine magnetite and/maghemite produced during pedogenic processes, as a by-product of the metabolism of iron-reducing bacteria [54,55]. The presence of pedogenic magnetic particles in the T6 profile agrees with the ubiquitous occurrence of rhizoliths (root moulds and concretions) observed in the field. However, these pedogenic particles can also have a detrital origin, from the incorporation of soil sediments or clays.

Isothermal Remanent Magnetization Curves
Isothermal Remanent Magnetization (IRM) measurements were conducted on eight samples from the T6 profile in order to investigate the nature and origin of the ferromagnetic particles present in the sediment and their contribution in the short-term cyclical oscillations observed in the magnetic susceptibility curve. In addition, a sample of modern sediments was analyzed for comparison. Results are shown in Table 3 and Figure 9. After treatment of the raw data by the cumulative Log-Gaussian function [37], two main components are obtained in all samples. Component 1 has B 1/2 values of 21-24 mT and DP (dispersion parameter) of 0.25-0.30, corresponding to the coercivity range of detrital magnetite [56]. Magnetite contributes to 77-83% of the total remanence. Component 2 has B 1/2 values of 200-234 mT and DP of 0.50-0.55, typical of hematite [56]. Hematite contributes to 17-22% of the total remanence. Goethite has not been identified although the data are noisy above imparted fields of 1 T. The S-ratio varies between 0.86 and 0.92 and reflects the dominance of magnetite in the total remanence. Except for the slight variations of the S-ratio, no significant changes in the magnetic mineralogy (coercivity and grain-size) was noted within the eight analyzed samples, nor with modern sediments (Figure 9). However, there is a clear correlation between magnetic susceptibility and saturation isothermal remnant magnetization (SIRM), suggesting that magnetic susceptibility is mostly controlled by the concentration of ferromagnetic iron oxides (magnetite and hematite) ( Figure 10). In addition, samples having the lowest magnetic susceptibility values also have lower S-ratio ( Table 4), suggesting that the short-term cyclical variations observed in the magnetic susceptibility curve reflect variation in the relative proportion of magnetite and hematite.
Quaternary 2019, 3 16 of 29 Except for the slight variations of the S-ratio, no significant changes in the magnetic mineralogy (coercivity and grain-size) was noted within the eight analyzed samples, nor with modern sediments (Figure 9). However, there is a clear correlation between magnetic susceptibility and saturation isothermal remnant magnetization (SIRM), suggesting that magnetic susceptibility is mostly controlled by the concentration of ferromagnetic iron oxides (magnetite and hematite) ( Figure 10). In addition, samples having the lowest magnetic susceptibility values also have lower S-ratio ( Table 4), suggesting that the short-term cyclical variations observed in the magnetic susceptibility curve reflect variation in the relative proportion of magnetite and hematite.    Table 4. Rock magnetic parameters calculated from the analysis of the IRM curves (see Figure 9). The % column represents the percentage of the contribution of the magnetic phase; SIRM is the IRM at saturation is proportional to the concentration of the magnetic phase; B1/2 corresponds to the coercivity and DP is the dispersion parameter of the Gaussian curves.   Table 4. Rock magnetic parameters calculated from the analysis of the IRM curves (see Figure 9). The % column represents the percentage of the contribution of the magnetic phase; SIRM is the IRM at saturation is proportional to the concentration of the magnetic phase; B1/2 corresponds to the coercivity and DP is the dispersion parameter of the Gaussian curves.

Phytolith Analysis
Phytoliths were identified in those carbonate concretions corresponding to samples with codes T6FE0.72, T6FE2.00, T6FE2.90, T6FE3.48 and T6FE4.12, but constituted fewer than 10 units per observed concretion, which limits their potential for environmental interpretation. The identified morphotypes are predominantly larger than 20 µm, of Bulliform and Cuneiform types ( Figure 11) and rare Elongates and Trapeziforms. The Bulliform and Cuneiform types are usually produced by grasses undergoing hydric stress. The poor preservation of the studied phytoliths is probably due to the alkalinity of the concretions, although there is some evidence of aeolian abrasion. Quaternary 2019, 3 18 of 29

Phytolith Analysis
Phytoliths were identified in those carbonate concretions corresponding to samples with codes T6FE0.72, T6FE2.00, T6FE2.90, T6FE3.48 and T6FE4.12, but constituted fewer than 10 units per observed concretion, which limits their potential for environmental interpretation. The identified morphotypes are predominantly larger than 20 μm, of Bulliform and Cuneiform types ( Figure 11) and rare Elongates and Trapeziforms. The Bulliform and Cuneiform types are usually produced by grasses undergoing hydric stress. The poor preservation of the studied phytoliths is probably due to the alkalinity of the concretions, although there is some evidence of aeolian abrasion.

Palynology
From a general point of view, the prepared residues from both the T6 and T5 terraces are very poor in organic matter and pollen. A small number of highly degraded pollen grains were observed, of which some are indeterminable because of their damaged state. The small number of identified pollen grains and, in some cases, the sterility of the samples, precludes the construction of pollen diagrams. The results are presented in Table 5, organized by sample; taxa for each sample are divided into arboreal types, followed by shrubs and herbaceous and, finally, indeterminate pollen grains.

Palynology
From a general point of view, the prepared residues from both the T6 and T5 terraces are very poor in organic matter and pollen. A small number of highly degraded pollen grains were observed, of which some are indeterminable because of their damaged state. The small number of identified pollen grains and, in some cases, the sterility of the samples, precludes the construction of pollen diagrams. The results are presented in Table 5, organized by sample; taxa for each sample are divided into arboreal types, followed by shrubs and herbaceous and, finally, indeterminate pollen grains.  The two samples from the lower part of the T5 sequence from outcrops at the Foz do Enxarrique site (T5FE015-0.23 and T5FE0.92-1.00) only allowed the identification of a pollen grain of Pinus sp. For the three samples collected from the basal layer (archaeological strata) of the T6 upper member only Chenopodiacea could be identified. In the samples collected in the upper part of T6, the few taxa identified indicate the presence of Pinus sp., the only tree species, Ericaceae (shrubs) and, in the herbaceous group, Poaceae, Artimisia, Chenopodiaceae and Asteraceae tubuliflora. In conclusion, the results obtained from the palynological study of the T6 and T5, from oxidized sediments, revealed very few palynomorphs and do not allow any palaeoenvironmental interpretation.

Geochemical Analyses
The results of geochemical analyses performed with the Thermal Scientific Analyzer Nikon XL3t spectrometer, on samples collected from the upper unit of T6, indicate that Ca has high values on the calcium carbonate concretions/rizoliths and very low values on uncemented siliciclastic sediment (coarse silt to fine sand) (Appendix A, Figure A1).

Palaeolithic Sites and Fauna Associated
No evidence has been found to date of industries prior to the Acheulean in this area. However, considering the presence of Oldowan material at several sites in Iberia, some dating back to c. 1.4 Ma, this absence must be seen as a problem of archaeological visibility [57].
The oldest archaeological evidence in this region appears in the basal deposits of T4, which are dated as older than 280 ka and perhaps as old as c. 340 ka [26,28]. A large Acheulean quartzite assemblage, with handaxes, cleavers, very large flakes and respective cores, was recovered from Monte do Famaco (Vila Velha de Ródão), in a colluvium resulting from the erosion of basal T4 deposits [58]. Test pits were excavated at this site and a large assemblage was also collected from the surface, but no thorough excavation was undertaken, nor has the assemblage been published in detail. Other evidence of Acheulean (e.g., identification of single or small numbers of handaxes and cleavers, sometimes along with quartzite flakes) has been reported from T4 deposits at Vilas Ruivas (Arneiro), but also at the surface of Pegos do Tejo-1, Monte das Nove Oliveiras, Monte do Arneiro and Monte da Cabeça Gorda; all these sites at Arneiro are still to be excavated [59,60].
The earliest Mousterian (Middle Paleolithic) industries occur on both sides of the Tejo in the uppermost deposits of T4. At Pegos do Tejo-2 (probably dating from c. 200-170 ka), the identification of a ring of cobbles was interpreted as a possible hearth (fire pit). These occur along with Levallois, discoidal, opportunistic and bipolar cores, flakes (often thick), denticulates, notches and sidescrapers [61,62]. A similar age (c. 165-155 ka) was attributed to Cobrinhos (Vila Velha de Ródão), a dense lithic assemblage spread over 1300 m 2 and dominated by Levallois, discoidal, radial and opportunistic cores and Levallois and Kombewa flakes (and also some points and blades), pseudo-Levallois points, sidescrapers, denticulates and notches [22].
At Monte da Revelada (Vila Velha de Ródão), also at the top of T4, the lowest archaeological context has discoidal and Levallois cores, flakes with dihedral and faceted platforms, along with radial dorsal patterns, including some Levallois flakes and pseudo-Levallois points that can be associated with a Mousterian occupation. These artefacts occur along with oval stone features that that again may correspond to hearths. At the nearby site of Alto da Revelada, Mousterian artefacts occur scattered on the surface, on the Cabeço do Infante Formation (Paleogene), in the western area where it has not been covered by aeolian sands [63].
In the T5 deposits at Vilas Ruivas (Arneiro), now dated to c. 117-71 ka, it was again possible to recognize hearths, windshields and post-holes [59,60], associated with a small but well-preserved Mousterian assemblage. This was composed of Levallois, discoidal and centripetal cores, Levallois flakes, pseudo-Levallois points, some blades and points, denticulates, notches and sidescrapers, mostly quartzite, but without faunal remains.
At the Foz do Enxarrique site, several archaeological campaigns were undertaken between 1982 and 2001. A single archaeological level was identified, at the base of the T6 upper unit, with a diverse Pleistocene faunal assemblage associated with a rich Mousterian industry represented by more than 10,000 artefacts [48,64]. Among them are Levallois, discoidal and centripetal cores, Levallois flakes and pseudo-Levallois points, plus blades and points, denticulates, notches and sidescrapers, most of them produced from quartzite (c. 66%), quartz (c. 25%) and flint (c. 9%). A complete sequence of debitage is present. The explotation of raw materials was interpreted as resulting from local and opportunistic procurement [48].
The Foz do Enxarrique faunal assemblage from the campaigns carried out between 1982 and 1999, a total of 785 bone specimens, was studied by Brugal and consists of red deer (Cervus elaphus-54.8%), horse (Equus caballus-36.5%), aurochs (Bos primigenius-2.2%), elephant (1.2%), rhinoceros (0.5%), rabbit (4.2%), and a very few carnivores such as fox (0.5%), hyena (0.2%) and lynx (0.2%) [65][66][67][68]. The fauna collected during later campaigns (2000 and 2001) has the same characteristics: a clear predominance of deer, horse and aurochs; the predominance of these herbivores and the weak occurrence of carnivores points to an anthropogenic accumulation [69]. This seems to be reinforced by the distribution of the anatomical units of the skeleton identified from the most commonly represented species (Cervus and Equus). The cranial specimen is the most represented (47.3% for the deer and 75.3% for equidae), the axial and appendicular skeleton representing only 52.7% and 24.7%, respectively [66]. This suggested the idea that the site served as a hunting zone, because the edible parts (associated to the appendicular and axial skeleton) have been removed. Such interpretation is more difficult for the megafauna assemblage as, for instance, in the case of the elephant remains (classified as Elephas antiquus (=Paleoloxodon antiquus)) that are represented by a lamella of a superior molar and four bone fragments [66][67][68]. Worthy of note is the occurrence of birds, fish and mollusc remains [66]; however, in a subsequent study no bird bones have been identified [69]. Amongst the studied remains, just over 50% are small fragments of unidentified bones; there are 42% of identified remains, which represents a high rate compared to other Paleolithic sites, where the fauna remains appear very fragmented (10 to 20% of total remains: TNR) [66]. About 13 specimens were identified with marks from the action of carnivores and 15 with marks from human butchery (nine with cutmarks and six with evidence of burning) [66]. Results from use-wear study on the lithic assemblage performed on 110 quartzite artefacts showed a dominance of wear traces associated to butchering, cutting of soft-animal tissue and scraping on bone activities [70]. Based on this, the authors conclude that the site was primarily a hunting zone, in which the animals were killed, quartered and then mainly taken to other places.
The Carregueira Formation (the cover unit of aeolian sans) contains Late Palaeolithic to Neolithic industries. At Monte da Revelada, the upper bed is composed of aeolian sands, disturbed by ploughing, and contains a mixture of Epipaleolithic and Neolithic material. A similar context, but here with backed bladelets, was found within the aeolian sands covering the eastern area of Alto da Revelada. At Vilas Ruivas, the aeolian unit, again disturbed by ploughing, has a Late Upper Palaeolithic or Epipaleolithic industry, in which flint occurs in greater quantity and the lithic assemblage has prismatic cores, blades and bladelets, small flakes, burins and endscrapers. Therefore, it is possible that the MIS 2 aeolian sands of the Carregueira Formation may preserve the first modern human occupations in the region.

Discussion
Regarding the environmental interpretation of the T6 record in the study area, the lower boulder-pebble gravel, 0.4 m thick and overlying a strath cut in metamorphic basement, corresponds to the coarse river-bed sedimentation near the margin of the energetic Tejo palaeochannel, probably during the interval c. 60-45 ka (MIS 3). No lithic artefacts or fossils, which could have helped in the interpretation, were found in this bed.
The T6 upper unit, mainly consisting of fine to very fine sands, grading upwards to coarse silt, is attributed overbank sedimentation. The detailed environmental interpretation of the various stratigraphic subunits is discussed in the following paragraphs.
The c. 20 cm lower bed (at 5.60-5.40 m depth), comprising Mousterian artefacts and fossil bones in a matrix of gravelly micaceous fine sands, is interpreted as overbank deposits close to the channel margin that also record hominin activities of hunting and butchery [48,64,66,70]. During this period, the Tejo channel moved laterally towards the west, preserving this record of human occupation. Thus, this confined place at the confluence of the Enxarrique and Açafal streams was used by animals for drinking and they were easily hunted. As this thin layer, now dated to 44 ± 3 ka, represents the last regional evidence of Mousterian industries and the megafauna, it may correspond to a cold and dry period that negatively impacted the animals and the Neanderthals. It is possible the Neanderthals relied heavily on some specific biotic resource that may have been reduced during cold climatic conditions, so that they faced difficulties in adapting. The remaining Neanderthals could have been induced to move toward the better climate conditions of SW Iberia and were later absorbed into modern human populations. Also relevant to this discussion is the fact that the ages of the stratigraphic levels with the youngest Mousterian industries in westernmost Iberia are progressively younger toward the SSW: the Cantabrian region, c. 48 [58,[72][73][74][75]; Murcia, 37 ka [11]; Gibraltar, between c. 33 and 24 ka [9].
The sediments of the 5.40-4.55 m depth interval within T6, consisting of micaceous fine sands with some interbedded thin gravel stringers, are also attributed to fluvial overbank deposition close to the channel margin. No fossils or artefacts were found here or in the upper deposits of T6. The onset of a new sedimentation phase (overbank fine-grained sediments) without artefacts is not necessarily evidence for cultural breakdown. However, in other areas of the Lower Tejo where T6 is preserved, no Mousterian artefacts were found in younger stratigraphic levels.
Regarding the T6 upper unit, the coarse silts from 4.55 cm depth to the surface are attributed to overbank sedimentation, but some characteristics point to the possibility of short-distance transport by wind, namely the lack of lamination and of erosion surfaces, the absence of dispersed coarser grains, the low clay content, a mean grain size in the coarse silt range, the fine skewed distributions and evidence of aeolian abrasion provided by phytolith analysis. The literature shows grain-size-distribution curves of loess deposits to be very similar to these silty T6FE deposits (e.g., [76][77][78][79][80]. Possible short-distance aeolian transport of exposed overbank fines could have been promoted by strong winds coming from the west and penetrating through the Ródão gorge ( Figure 1). However, if there had there been significant aeolian transport, the resultant sediments should also cover higher terrace levels, and this is not evidenced by field observation.
Sediment magnetic properties have been widely applied to fluvial sediments and loess and may provide useful information about fluvial activity, climate and environmental changes, as well as pedogenesis [50,51,[81][82][83]. In fluvial sediments from Beijing, for example, high magnetic susceptibility values generally reflect warm-climate conditions, whereas lower values match colder periods [84]. In wind-blown sediments and buried soils from southern Siberia [81], colder high-wind periods that are associated with an absence of soil formation show low values of frequency-dependence of magnetic susceptibility, whereas higher values are observed in episodes with less wind. In the classic loess-palaeosol sequence of the central Chinese Loess Plateau, there is a striking correlation between magnetic susceptibility and grain size [85], which are good indicators of summer and winter monsoon intensity respectively [86,87]. In general, ferromagnetic crystals in soils derive from both primary (detrital) and secondary (enhanced) iron minerals. The latter are most often of stable single-domain size or less and associated with the clay fraction, whereas the former are usually associated with sand and coarse silt-size fractions [49]. Regarding the studied T6 sediments, no significant correlation (R 2 = 0.0261) is observed between sedimentary grain size and magnetic susceptibility, suggesting that the latter is not controlled by mineralogy (ex. paramagnetic clay and phyllosilicates versus diamagnetic quartz) ( Figure 5). In contrast, a striking correlation between SIRM and magnetic susceptibility (R 2 = 0.9554 for magnetite; R 2 = 0.769 for hematite) indicates that magnetic susceptibility is dominantly controlled by the iron-oxide content, in the form of magnetite and hematite ( Figure 10), probably of a detrital origin. The proportion and contribution of magnetite and hematite in the bulk remanence is mostly similar in all samples, as illustrated by narrow S-ratio values of 0.86-0.92. In particular, the IRM curves of the T6 samples are very similar to those of modern sediments (Figure 9), suggesting a common source for the entire sedimentary profile during the last millennium. Superparamagnetic particles, generally interpreted as a product of pedogenic processes, are present in all samples and may have been formed in situ or transported from the surrounding soils. The presence of numerous rhizoliths observed in the field suggests that pedogenic magnetic particles may have precipitated in-situ, during soil development. However, the poor correlation (R 2 = 0.0145) between mass specific magnetic susceptibility and frequency-dependent magnetic susceptibility indicates that pedogenic processes alone cannot explain the short-term cyclical variations observed in the magnetic susceptibility curve of the T6 profile ( Figure 4). Conversely, a slight but significant correlation (R 2 = 0.5023) between the S-ratio, i.e., the relative proportion of magnetite versus hematite, and magnetic susceptibility imply that the short-term cycles may correspond to changes in weathering regime and climate. More exactly, warmer/drier periods would enhance the oxidation of magnetite (or maghemite) and promote precipitation of hematite, and the reverse. Because the magnetic susceptibility of hematite is lower than magnetite, this provides a potential explanation for the cyclical oscillations observed in the magnetic susceptibility curve.
The 0.60-5.40 m depth interval contains c. 30 thin levels of calcium carbonate concretions and rizoliths intercalated in uncemented coarse silt. The characteristics of these levels, which have no evidence of any erosive surface and dip toward the palaeochannel, progressively increasing in thickness, point to a secondary origin for the carbonate concretions. A relatively stable surface and a certain amount of rainfall during the represented periods are indicated. The most probable source of calcium carbonate for these pedogenetic concretions is the dissolution of dolomitic and calcium carbonates that occur at the base of the Cabeço do Infante Formation (Paleogene), which crops out at a short distance from the site. Upstream sources of calcium carbonates are located at least at 200 km away (in the Madrid Cenozoic basin).
The results obtained from palynological study of the T6 and T5 deposits do not allow palaeoenvironmental interpretation. Regarding the interval represented by the T6 terrace, the palynological record of the MD95-2039 ocean core points to an open landscape with steppe vegetation and low values of tree pollen, suggesting a severely cold and dry climate, during 61-59 ka; during the interval 57-31 ka, there are fluctuations in the expansion and contraction of arboreal pollen and Ericaceae related to the alternation of warmer and humid conditions during the interstadials and the cold and dry stadial minima of the last glacial cycle [92]. In the Puente Pino sequence (Toledo, Spain), for the 42-34 ka interval, declining woodland and the increasing herbaceous pollen taxa are observed, related to adverse climate conditions of cold and dry character [97]. The palaeoenvironmental data provided by these two adjacent contexts support the possibility of an open landscape in the region of the study area and could explain the several levels with pedogenic calcium carbonate rizoliths and concretions in the T6 upper unit.
The mineralogical composition of the <2 µm fraction of the samples collected from T6 and T5 are similar, consisting of smectite, illite and kaolinite, although T6 has less smectite. Even if a significant part of the clay minerals could have been sourced from erosion of the local Paleogene and Miocene formations, very rich in smectite [23,25], this clay mineral association seems compatible with the regional climatic conditions during MIS 5-3 [91,92].
By 32 ka, the climate had changed to cold and dry conditions and aeolian deposition dominated the valley landscape, preserving Upper Palaeolithic industries.

Conclusions
Updated ages from the three Lower Tejo terrace sequences, containing Mousterian industries, were obtained by pIRIR, as follows: (i) OSL dating of the oldest Mousterian industry, stratigraphically situated in the uppermost T4 deposits, suggests a probable age of c. 200-170 ka for the arrival of the Neanderthals in this region, probably by way of the Tejo River valley from central Iberia; (ii) T5 dates from c. 140 ka at the base and 70 ka at the top; (iii) T6 dates from c. 60 ka at the base and c. 35 ka near the top; (iv) The new date of 44 ± 3 ka for a level located at the base of the T6 upper unit records the last regional occurrence of a Mousterian industry and of the megafauna.
The T6 upper unit (so far without Mousterian or Early Upper Palaeoleolithic industries), consists of fine sands to coarse silts interpreted as overbank sediments. It has a large number of intercalated thin levels of carbonate concretions and rizoliths, suggesting episodic evaporation and development of paleosols in a seasonal dry period, in agreement with the occurrence of phytoliths.

Supplementary Materials:
The following is available online at http://www.mdpi.com/2571-550X/2/1/3/s1. Table S1-Grain-size statistical parameters and mass specific magnetic susceptibility of samples collected from of the T6 upper unit at Foz do Enxarrique site.