A multi-technique approach to characterization: the Sant Martí de Tous chert as a prehistoric resource for the NE of the Iberian Peninsula

The Sant Genís Formation is located in the NE of the Iberian Peninsula (Catalonia, Spain) and is dated to the Priabonian (upper Eocene), being part of the evaporitic formations of the margin of the Ebro Basin. It is formed by a succession of sandy lutites, occasional limestone layers, marls, and local stratified gypsum and cherts, including the Sant Martí de Tous chert. The Sant Martí de Tous chert type is confirmed by its abundance at specific locations within the territory (NE Iberian Peninsula). This is an important raw material procurement area, as evidenced by the presence of this chert in the main prehistoric sites of the region (e.g., Abric Romaní) and the constant discovery of new sites in the area around the Sant Genís Formation, especially from the Neolithic period onwards (e.g., Cal Sitjo, La Guinardera Nord workshop). All these features, together with the great heterogeneity of the Cenozoic evaporitic cherts, prompted us to carry out a multi-proxy characterization for obtaining a valid criterion by which to identify this chert in the archaeological record. This very heterogeneity hampers the macroscopic characterization of archaeological cherts (e.g., at Abric Romaní), but through petrographic analysis we have been able to identify their origins and, albeit to a lesser extent, their lithostratigraphic unit of provenance within the formation. The complementary mineralogical and geochemical techniques applied in this study show that, although it is difficult to establish an exact origin, through intensive sampling and criteria such as the differential presence of Fe and Ca, some differences among the siliceous varieties within a formation can be detected.


Introduction
The Sant Martí de Tous chert (SMT) was widely exploited as a raw material in prehistoric times due to its quality for knapping and its high abundance within a well-delimited area. Indeed, this area became a key location for raw material procurement and stone-tool configuration. The exploitation of this territory has been documented through the continuous use of this chert type. From the Middle Paleolithic (e.g., Abric Romaní) through the Upper Paleolithic (e.g., Esquerda de les Alzines), the Mesolithic (e.g., Abric Agut, Cal Sitjo), the Neolithic (e.g., Vilars de Tous, Cal Sitjo), and even during modern times, when it was employed for the configuration of gunflints (e.g., La Guinardera Nord) (Clop et al. 2005;Rey-Solé et al. 2014;Gómez de Soler et al. 2020a, 2021. The SMT chert, an evaporitic lithotype from the Sant Genís Formation (Priabonian, upper Eocene), has been studied in recent years in what is considered a prominent chert outcropping area, concentrated in a small region. The area acted as a territorial organizer through the procurement of lithic raw materials among the prehistoric hunter-gatherers and the first farmers who inhabited the northeastern Iberian Peninsula. Located between inland and coastline, north and south, the area connects different ecotones (transitional space between two ecosystems), suggesting that it served as a strategic transit area. This is evidenced by the high number of sites located in the region from the Upper Pleistocene forward and the exponential increase in occupations from the Neolithic onwards. Relevant examples include the Middle Paleolithic sites of Mediona (Estévez et al. 1993), La Noguera (Sant Quintí de Mediona) (Freixas and Mir 1987) and Teixoneres cave (Moià) (Picin et al. 2020), the archaeological sequences of Cova Gran (Collbató) and Simanya Gran (Sant Llorenç del Munt) , Cova Foradada and Balma de la Griera (Calafell), Guineu (Font-Rubí) and Abric de la Consagració (Capellades) (Oms et al. 2016;Morales et al. 2019Morales et al. , 2022Fullola et al. 2020;Vaquero & Picin pers. comm.), the Upper Paleolithic sites of Can Manel and Pinyons (Capellades) (Vaquero et al. 2013), the Mesolithic site of Abric Agut (Capellades) (Vaquero et al. 2006), and the Neolithic sites of Can Sadurní and Guixeres del Penedès (Vilobí del Penedès) (Edo et al. 2012;Oms et al. 2021). Despite the diachrony, all these sites supplied chert from this formation, becoming a pole of attraction and a vertebra of the territorial structure over time. Obtaining a multi-technique dataset on SMT chert characterization will fulfil an important need, allowing different researchers to identify this chert type frequently found in the archaeological sites of the NE Iberian Peninsula.
In this work, we develop the first multi-scalar approach to the chert of the Sant Genís Formation (Priabonian, Cenozoic) of evaporitic origin with a view to obtaining a detailed characterization of the different lithological units yielding it. From a petrographic perspective, Cenozoic cherts of evaporitic origin are characterized by great heterogeneity, the absence of microfossils, and a high presence of post-genetic alterations such as the development of white patina, which frequently make it difficult to undertake a correct archaeological ascription.
SMT chert (Sant Martí de Tous, Barcelona) is found within evaporite lithofacies (sabkha-like), defining a depositional environment comprising shallow saline lakes with low-concentration brines of calcium sulfates, where massive gypsums, gypsiferous marls, carbonate sedimentation, gypsarenites, and red and gray lutites are deposited, depending on the horizontal location (lake center or margin) and stage (dilution, exposure) (Ortí et al. 2007). SMT chert outcrops are located within lutites, gypsums, marls, and calcarenites. This heterogeneous lithological context motivated our multimethod analysis of their petrographic and geochemical traits, our aim being to explore the homogeneity or heterogeneity among cherts from the same geological formation but different enclosing lithologies. Heterogeneity in evaporitic cherts poses a particular problem when it comes to their macroscopic characterization and ascription to archaeological raw material groups. Determining their diagenetic variability will make it possible to define these groups more precisely and reliably, so that the following analyses (thin section, chemical characterization, etc.) will show the reality of the SMT type chert as a whole rather than sample specificities that can lead to confusion.
We petrographically analyze 46 samples (enclosing rocks and cherts of geological and archaeological origin), of which six chert types representing the main lithologies of the formation are tested by multiple material analyses. Techniques including X-ray diffraction (XRD), magnetic property analysis, portable X-ray fluorescence (p-EDXRF), inductively coupled plasma optical emission spectroscopy (ICP-OES), and neutron activation analysis (NAA) are applied to explore the mineralogical and geochemical variability of the samples. The data obtained are compared with three archaeological chert samples (AR'01-M-Q47-49; AR'02-M-U47-220; and AR'02-M-P53-537) from level M (between 51.8 ± 1.4 and 55.8 ± 2.3 ka BP) of Abric Romaní, which were previously assigned to this formation (Gómez de Soler et al. 2020a), to test the ascription and confirm the validity of this multi-method approach.

Geological setting
The Sant Genís Formation (Colldeforns et al. 1994) is a Cenozoic unit that forms part of the sedimentary infill of the Ebro Basin, bordered by the Pyrenees and the Basque-Cantabrian Ranges to the north, by the Iberian Range to the south, and by the Catalan Coastal Range to the east. The depression drained mainly by the Ebro River is a Tertiary sedimentary basin with Paleogene deposits formed by erosive processes affecting the bordering ranges (Solé Sabarís 1964).
The Sant Genís Formation is of Priabonian age (upper Eocene) and is part of La Noguera lacustrine system (Anadón et al. 1989), falling within the Upper Gray Lutite unit (Saez 1987). At present, it extends a length of 12 km, with a maximum width of 1.5 km, showing a NE-SW orientation from the locality of Sant Genís (Jorba) through the locality of Fillol (Sant Martí de Tous), and covering a surface area of 5.6 km 2 (Fig. 1). It comprises a succession of 400 m of red sandy lutites, with horizontal lamination and ripples, including occasional limestone layers with charophytes and gastropods (Colldeforns et al. 1994). In geological maps, it is presented as two formations; one formed by red sandstones and marls, and the other by gray sandstones and marls, with local stratified gypsum (IGME 1975).
Silicification outcrops are associated with the unit of red sandy lutites and occasional tabular red sandstones. At the base of the formation, alternating layers of nodular and laminated gypsums are described, and to the SW of the locality of SMT, this changes to an alternation of massive gypsums with nodular cherts (Ortí et al. 2007).
At its base (NE), the geological section of the region shows the Artés Formation (lower Priabonian), formed by an alternation of marls and red sandstones. Above this lies the Sant Genís Formation, represented at its base by the 20-m-thick SMT gypsum, and towards the top by an 80 m deposit of upper Priabonian marls, red sandstones, and sporadic gypsum. On top (SE of the section), the Intermediate Lacustrine Units (Colldeforns et al. 1994) are located, represented by the Copons Formation (Sannoisian, lower Oligocene), consisting of gray marls, sandstones, and gypsums (Fig. 2).

Materials and methods
The chert outcrops of the Sant Genís Formation are located to the SW of the town of Sant Martí de Tous. Two sectors with a considerable chert concentration can be differentiated, the sectors of Guinardera and Fillol, which include a total of 14 chert outcrops (eight associated with the former sector and six with the latter). The names of these outcrops relate to the nearest toponyms (see Fig. 1).

Survey and sampling
The cartographic and geological data available for the area were analyzed to identify the different chert-bearing units. Once located, the chert outcrops were registered in a specifically designed database, where their main characteristics (location, accessibility, geological context, etc.) were recorded. In addition, the outcrop locations were imported in a vector format onto a digitized map with geological and cartographic data (scale 1: 50,000) to create a specific cartography of the chert availability in the area.

Petrographic analysis
The morphoscopic analysis consisted of naked-eye and stereoscopic microscope (ZEISS Stemi DV4 8-32x) observations. We analyzed external aspects such as the type of support (e.g., block, pebble, preform, core, flake) and the cortical area (presence/absence, type of cortex, morphology, lithological origin), and internal aspects such as the color (using the Munsell Rock Color Book, Munsell Color 2013) and its distribution within the surface, transparency, texture, the presence of sedimentary structures, fissures and fractures, post-depositional alterations, impurities, and the microfossil content (see Soto et al. 2014).
Petrological thin sections were analyzed by polarized light microscopy (Nikon Eclipse E400 POL 40-400x) to determine the mineral composition, textural characteristics, and impurities content, thus analyzing the siliceous matrix as well as noting the presence or absence of relict and authigenic minerals and possible alterations or compactions.
A total of 46 petrographic thin sections were prepared and analyzed at the IPHES-CERCA Geoarchaeology Laboratory; these comprised 22 geological chert samples, 13 enclosing rocks, and 11 archaeological samples.

X-ray diffraction (XRD)
XRD analyses were performed at the Scientific and Technical Resources Service of the Universitat Rovira i Virgili (URV) in Tarragona (Spain), using a Siemens D5000 diffractometer (Bragg-Brentano parafocusing geometry and vertical θ-θ goniometer) fitted with a curved graphite diffractedbeam monochromator, an incident-beam Soller slit, a 0.06° receiving slit, and a scintillation counter as a detector. The angular 2θ diffraction range was between 5° and 70°. The data were collected with an angular step of 0.03° at 6 s per step and sample rotation. Cuk α radiation was obtained from a copper X-ray tube operated at 40 kV and 30 mA.
Identification of the minerals was achieved by comparison of the XRD diffractogram with the ICDD database (release 2007) using Diffrac plus Evaluation software (Bruker-AXS 2009).
Quantitative phase analysis was derived from XRD powder patterns using the Rietveld (Rietveld 1969) software TOPAS (Bruker-AXS 2009). The instrumental contribution was introduced with the fundamental parameter approach (Cheary et al. 2004). The preferred orientation effect was corrected using the March- Dollase model (March, 1932) (Dollase 1986). For the moganite and calcite phases, the [011] and [104] directions were introduced respectively; for gypsum, the [010] direction was employed. The background for each diffraction pattern was fitted as a 5th-degree Chebyshev polynomial. For every profile, the zero shift and overall scale factor, cell parameters, and Lorentzian contribution to the peak width produced by the crystallite size were refined.
The quantitative phase analysis using the Rietveld method consists in fitting a calculated X-ray diffractogram to the experimental one and minimizing, by the least-square method, the difference between them. All the algorithms that are used in the Rietveld method try to minimize this difference between the intensity at each point of the observed profile (I i ) and the calculated intensity at the same point (I ci ) of the diffractogram. The quantity to be minimized during the refinement is: where i = 1/I i . The sum is performed for all points in the diffractogram. The calculated X-ray diffractogram is obtained from the crystal structure of each phase plus geometrical considerations of the diffractometer.
The following table shows the crystal structures employed for each phase ( Table 1).
The weight fraction of a w i phase in a system of n crystalline phases is given by the expression: where ρ i is the crystal density, V i is the unit cell volume, and S i is the refinable scale factor for phase i. The previous expression is applicable when all phases considered in the sample are crystalline so that:

Magnetic properties
A bulk sample (~ 500 mg) from each chert sample was extracted to analyze its magnetic properties. The analyses were carried out using a variable field translation balance (VFTB) and comprised the following sequence: (i) progressive isothermal remanent magnetization (IRM) acquisition curves, (ii) hysteresis loops (± 1 T), (iii) backfield coercivity curves, and (iv) thermomagnetic curves up to 700 °C in air. All these analyses were performed at the laboratory of paleomagnetism of the Universidad de Burgos (Spain).
After correcting the hysteresis cycles for their dia/paramagnetic fraction, various rock-magnetic parameters were determined using the Rock Mag Analyzer software (Leonhardt 2006). These included Ms (saturation magnetization), Mrs (remanent saturation magnetization), and Bc (coercive field). The remanent coercive field (Bcr) was determined separately from the backfield coercivity curves. In combination with other analyses such as thermomagnetic curves, these parameters are used as indicators of ferromagnetic (sensu lato, s.l.)strongly magnetic-mineral concentration and composition in the sample (Evans and Heller 2003). The Curie/Néel temperatures were determined from the thermomagnetic curves following the two-tangent method of Grommé et al. (1969). Another proxy considered was the S-ratio (Bloemendal et al. 1992), which is defined here as IRM@0.3 T/IRM@1 T, taking 1 T as the maximum available applied field. The S-ratio is generally used to constrain the ratio between high-coercivity (e.g., hematite/goethite) and low-coercivity (e.g., magnetite/ maghemite) minerals.

Portable energy-dispersive x-ray fluorescence (p-EDXRF)
Chert samples were analyzed at the Instituto de Ciencia de los Materials (ICMUV) of the Universitat de València (Spain), using a portable EDXRF spectrometer with an X-ray tube (Moxtek, Inc.) with a silver anode operating in transmission mode at a voltage of 30 kV, a current of 10 A and an acquisition time of 240 s. The spectrometer has a thermoelectrically cooled Si-PIN detector (Amptek, Inc.), with an active area of 6 mm 2 , a 12.5 m beryllium window, and an energy resolution of 165 eV (FWHM @ 5.9 keV). An aluminum collimator collimates the X-ray beam, and the beam diameter on the surface of the sample is about 5 mm.
These elements were fitted together on a mechanical device that enabled us to place the system in front of the analysis point and maintain the same geometry, ensuring the reproducibility of the measurements. The X-ray beam impinged perpendicularly upon the sample, and the tube and detector were placed with an excitation-detection geometry of 45° and a 2 cm sample-detector distance.
The selected samples were left untreated and analyzed at several discrete locations, using relatively flat areas of each chert sample with no cortex or inclusions. Net areas of the K-fluorescence lines generated at each location were calculated by means of the PyMCA analysis software (Solé et al. 2007).
The R environment (R Core Team 2021) and the "ggplot2" package (Wickham 2016) were used to create boxplots of the p-EDXRF counts for major elements in the samples. A principal component analysis (PCA) of the p-EDXRF readings was then performed using the R packages "FactoMineR" (Le et al. 2008) and "factoextra" (Kassambara & Mundt 2017).

Inductively coupled plasma optical emission spectroscopy (ICP-OES)
This technique involves vaporizing, dissociating, ionizing, and exciting the different chemical components of a sample within an argon plasma. When atoms are de-excited, electromagnetic radiation is emitted within the visible ultraviolet spectrum with characteristic wavelengths for each component. These radiations are separated according to their wavelength, and finally their intensity is measured.
ICP-OES analyses were performed at the Scientific and Technical Resources Service of the Universitat Rovira i Virgili (URV) in Tarragona (Spain), using an inductively coupled plasma optical emission spectrometer (ARCOS FHS16). These analyses consisted of two phases. The first phase was to dissolve the sample and perform the digestion with an Anton Paar microwave (Space Cube M627). The process was based on a 500 mg ground sample of chert with 2 mL of HNO 3 and 0.5 mL of HF and consisted of 5 min in the microwave. After various steps of five and 10 min with a total of approximately 40 min in the microwave adding 0.5 mL of HF in two stages, we ended with a total of 500 mg of chert, with 2 mL of HNO 3 and 1.5 mL of HF. Afterwards, the sample was fully diluted. The second phase consisted of direct measurement with an ICP spectrometer Spectro ARCOS FHS16. The following parameters were used for element detection: nebulization system cross-flow, plasma power up to 1400 W, pump pitch to 30, refrigerant flow at 13.00 (l/min), auxiliary flow at 1.00 (l/ min), nebulization flow at 0.90 (l/min), torch Z at 10.0 (mm), and measurement time of 28.0 (s). Line plots of the trace element values for each geological chert sample were created in R with "ggplot2," and ternary plots of Mg, Al, Fe and Al, Ca, Fe with the "ggtern" package (Hamilton & Ferry 2018).

Neutron activation analysis (NAA)
Using a corundum drill bit, about 80 mg of chert powder was obtained from each sample. The sample powder was then compressed to form pills, adding pure cellulose powder as a binder. The whole set of samples, together with six pills of the Bonn standard, was sent to the research reactor of the Reactor Instituut Delft (Netherlands) and was irradiated there for 10 h with a flux of 5 × 10 12 neutrons/(cm 2 s). The emitted radiation from each sample was measured three times during the following 4 weeks applying the routine NAA procedure (Mommsen et al. 1991) at the laboratory of the Helmholtz-Institut für Strahlen-und Kernphysik of the Universität Bonn (Germany). Altogether 30 elemental concentrations were determined. Due to the hardness of the chert samples, a certain number of rubbed-off parts of the corundum drill bit (pure Al oxide) were introduced into the samples, forcing the authors to interpret concentration ratios, instead of absolute concentrations.
Boxplots of the major elements (NA, K, Ca, and Fe) were created in R using "ggplot2." The measurements of the rareearth elements in the geological samples were included in a PCA ("FactoMineR" and "factoextra" R packages), and the results from the archaeological samples were included in the PCA space as supplementary individuals.

Survey and sampling
The stratigraphic sequence of the Sant Genís Formation allows the vertical succession of the outcrops and the horizontal relationship among the lithological units through the Guinardera and Fillol sectors to be identified (Fig. 3).
The Cal Perdut section is an 80 m sequence, presenting at its base 20 m of white marls with decimetric nodular cherts (PER-01), displaying a 20-30% silicification ratio. The overlying layer is formed of 2 m of secondary gypsum, where meganodular cherts are present (PER-02). Metric-sized erratic chert blocks, subangular in shape, confer a 100% chert ratio in this layer. Further up the section is a 15 m layer of calcarenitic rocks with discontinuous tabular chert at the top (PER-03), exhibiting a 10% chert ratio. At the top, a layer made up of alternating marls, green-red lutites, and calcarenitic rocks is found, with chert in both lutite units (the PER-04 and MDT outcrops). Chert appears in decimetric formats with a representation ratio of 10% (Fig. 4).
The Guinardera section is a 40 m sequence. At the base, corresponding to the top of Cal Perdut, 10 m of red lutites with sporadic red sandstones are described, yielding discontinuous tabular chert with a 5% silicification ratio (GUIN-01). This unit is followed by 30 m of gray lutites with chert at its base and top (GUIN-02 and GUIN-03), described as decimetric black and gray chert nodules with a 15% silicification ratio in each case (Fig. 5).
In the Fillol sector (Ortí et al. 2007), the Fillol and Cal Tomàs sections are described. The lower unit of the Fillol section is formed of 20 m of gypsiferous marls and secondary gypsum alternating with red and gray lutites and cherts (Cal Calaf (CAL) 01 and 02) with banded stratifications and nodular and meganodular structures. Chert is present in irregular shapes, with decimetric formats and an average chert ratio of 50-60%. Further up, an intermediate unit comprising 10 m of secondary massive gypsums with nodular and meganodular structures yields another chert outcrop (Fillol (FILL) 01). The cherts exhibit elongated nodular shapes and centimetric formats, with a 5-15% representation ratio (Fig. 6). At the top, a 40-m-thick succession of laminated calcarenites alternating with thin units of secondary gypsums without chert is described, covering the whole surface where the village of Fillol is located. This lithological succession, showing a progressive decrease in gypsum with an increase in calcarenites, indicates a shallowing of the sabkha basin, a displacement of its depocenter, and the location of a tidal plane in the Fillol section.
Applying the CAR method (Soto et al. 2018) to calculate the siliceous abundance in the territory, a volume of 0.0052 km 3 of chert in an exposed outcrop volume of 0.0207 km 3 was estimated, meaning a chert density in volume of 1.14% for the total exposed volume of the Sant Genís Formation (for more information see Gómez de Soler et al. 2020a).
Among the red lutite facies, GUIN-01 displays singular traits. Macroscopically, these cherts present a very fine texture, with a homogeneous red color (5R3/4 Dusky Red) and opaque diaphaneity and occasional fissures. Microscopically, they are formed of crypto-and microcrystalline quartz (55%) and abundant clay minerals (up to 40%). As an accessory mineral, iron oxides in the channel porosities (5%) are also present as massive hematite. All these characteristics allow us to define this chert type as a jasper (Fig. 12).

XRD
The geological chert samples analyzed by X-ray diffraction show contamination by tungsten carbide (WC), < 0.5 wt%, a chemical compound associated with grinding tools to powder samples when these are very hard (e.g., chert) (Fig. SI1). The results of the quantitative phase analysis for each sample are shown in Table 2, together with the conventional Rietveld agreement factor (R wp ): The percentages of the quartz phase are very homogeneous (88.6-97.5%). FILL-01 and MDT chert types stand out due to the slight increase in moganite content they show, and PER-03 due to the presence of calcite. These values with quartz and moganite as the only phases are typical of the cherts with fibrous silica (chalcedony) (Graetsch & Grünberg 2012) common among early-diagenetic cherts characteristic of evaporitic formations (Chowns & Elkins 1974). The presence of calcite in PER-03 may be due to a higher subaerial exposure, thus suggesting a more marginal environment, as its enclosing rock indicates.

Magnetic properties
The application of mineral-magnetic parameters to study the magnetic properties of geoarchaeological materials refers to the identification of iron oxides, hydroxides and some iron sulfides, and minerals generally occurring in trace amounts in many rocks and soils, as well as their variations in concentration and granulometry (Evans & Heller 2003). As expected, the studied cherts are characterized by dominant diamagnetic behavior. Diamagnetism is a property of all matter characterized by very weak magnetization antiparallel to the applied field. When the applied field disappears, diamagnetic materials such as quartz or calcite (typically dominant in cherts) do not retain a remanent or permanent magnetization (Dunlop & Özdemir 1997;Dekkers 2007). Nevertheless, despite being present in very small concentrations, some differences among the studied cherts can be distinguished by virtue of variations in the composition, concentration, and granulometry of their ferromagnetic mineralogy as we outline below. Table 3 compiles the main magnetic parameters calculated for this sample set. Ms values vary between 3.08 × 10 −3 and 8.08 × 10 −4 Am 2 kg −1 whereas Mrs values range between 7.59 × 10 −4 and 1.11 × 10 −4 Am 2 kg −1 . Bc varies between 134.24 and 9.91 mT, and Bcr between 322.74 and 22.64 mT. Like the Ms and Mrs values obtained from the hysteresis loops, the variation in the intensity of magnetization can also be inferred from the thermomagnetic curves. The magnetization value of the heating cycle at room temperature (J 30 ) ranges between 1.57 × 10 −3 and 3.64 × 10 −4 Am 2 kg −1 . In essence, the Ms and J 30 values are rather similar, and under the proviso that the magnetic mineralogy is relatively homogeneous as explained below, they can be considered indicators of the magnetic composition and concentration. As far as domain-state (1) Under PPL, samples show vacuolar porosity (vugs), associated with gypsiferous dissolution structures. Under XPL, the anhedral character of the nonuniform extinction of gypsum crystals (3a and 3b) and subhedral crystals and microcrystalline aggregates is observed; (2) shows typical microstructures of the alabastrine variety. B (1) Under XPL, intercrystalline porosity is observed, as well as calcite crystals infilling fissures and bioturbation structures (3a); (3b) detail of calcite crystals and porphyroblastic (2) and vein gypsum (5) microstructures under XPL, in addition to anhedral habits of non-uniform extinction, in a zone where a major accumulation of clay is observed (4) variations (granulometry) are concerned, we do not consider it realistic to represent the biparametric ratio (Mrs/Ms vs Bcr/Bc) in the so-called Day plot (Day et al. 1977). This is for two reasons: first, its interpretative limitations have been demonstrated (Roberts et al. 1995); second, as detailed below, some cherts also contain variable high-coercivity contributions, and the "Day plot" is theoretically only valid for (titano) magnetite-bearing assemblages.
A diagnostic parameter that qualitatively provides a first approximation to the type of ferromagnetic mineralogy present is the S-ratio ( Fig. 13 and Table 3), allowing highand low-coercivity minerals to be differentiated. The range of variation in this parameter is between 0 and 1 (or − 1 depending on the definition considered; see Bloemendal et al. 1992;Kruiver & Passier 2011). Values near unity indicate a predominance of soft or low-coercivity ferromagnetic minerals (e.g., magnetite), whereas lower values point to a greater contribution from hard or high-coercivity phases (e.g., hematite or goethite). In the present case, the lowest values correspond to FILL-01 (between 0.55 and 0.48) and the PER-04 sets (between 0.43 and 0.65). The other four chert types display values from 0.85 to 1.00, indicating the predominant contribution of a low-coercivity mineral to magnetization. Although S-ratio values close to 1 are often interpreted as representing only low-coercivity minerals, this does not necessary mean that hematite is totally absent (Bloemendal et al. 1992). Therefore, the interpretation of the S-ratio is not definitive, and the support of additional rockmagnetic analyses is required, as provided here. In any case, the S-ratio is a valid proxy as a first approximation for differentiating between soft and hard ferromagnetic minerals.
Although the magnetization is mostly dominated by a soft mineral, the occurrence of a high-coercivity phase identified in certain types of chert in the S-ratio is also confirmed by other experiments. The IRM acquisition curves of the FILL-01 (Fig. 14a) and PER-04 groups (Fig. 14d) do not reach saturation up to 1 T. That is, the magnetization continues to increase despite the increase in the applied field. The hysteresis loops of these types of chert are clearly wasp-waisted, indicative of the coexistence of low and high-coercivity   Fig. 14b and e). As expected, the highest Bc and Bcr values coincide with the groups (FILL-01 and PER-04) where the presence of highcoercivity minerals is most evident (Table 3). Among the generally weak magnetism of the studied collection, the other cherts are dominated by low-coercivity minerals (e.g., Fig. 14g-i).
These results are consistent with those observed in the thermomagnetic curves. With a view to obtaining smoother diagrams more amenable to interpretation, a three to five-point running mean was applied to the thermomagnetic curves. The thermomagnetic curves for most samples exhibit a drop in the heating cycle between 550 and 580 °C, consistent with the Curie temperature (Tc) of magnetite or of magnetite with a small amount of isomorphous substitution (Fig. 14i). Exceptionally, the FILL-01 samples exhibit an increase in their heating cycles from 400 to 420 °C, indicating the creation of magnetite during the laboratory heating (e.g., Fig. 14c). During the cooling a major increase in magnetization is found in these samples between 580 and 500 °C, indicating that additional secondary magnetite is created during the experiment. The high-coercivity phase observed in the IRM curves, hysteresis loops, and the S-ratio of the FILL-01 and PER-04 samples is hematite, according to the T N of around 660-675 °C detected in their thermomagnetic curves ( Fig. 14c and f). Table 4 shows the areas of the K-fluorescence lines of the elements detected in the spectra. The presence of W and Ni in the XRF spectra is not taken into account since both elements are present due to the milling process. Statistically, according to the intensity of the chemical elements analyzed, the boxplot shows that the six geological samples have homogeneous values. Exceptions are the value of Ca for PER-03, coherent with its enclosing rock consisting of calcarenites (as X-ray diffraction revealed), and the higher presence of Fe in samples FILL-01 and PER-04, as already shown by magnetic property analysis (Fig. 15). Also, the silica content shows that cherts from non-marginal positions within the sabkha (gypsiferous marls and gypsums) have a lower Si content, probably associated with  the process of silica precipitation and replacement (Ortí et al. 2007).

p-EDXRF
The PCA (components 1-2), representing 68.1% of the variance (Fig. 16a), highlights the major iron content in PER-03 and FILL-01. Also notable is that PER-03 and PER-04 are grouped according to their calcarenitic and lutitic enclosings. Likewise, MDT and CAL-01 are grouped due to their content in V. Lastly, PER-01 is revealed to be the most homogeneous in chemical composition. According to the PCA (56.8% of the variance) (Fig. 16b), the positive axes of components 1 and 3 are consistent with the previous PCA (components 1 and 2) (Fig. 16a). The negative axes again group the CAL-01, PER-01, and MDT samples according to their V content and the homogeneity in their chemical composition.

ICP-OES
A total of 24 trace elements were identified by means of ICP-OES. However, the Cr and Cd values were below the detection limits, and Ba was present only in samples CAL-01 and MDT. The P, V, Ag, Pb, Ga, and In values were obtained as a range (Table 5).
The values of Fe and Ca also show the similarities and differences among chert types revealed by this analytical technique. The Ca content is high for PER-03, as is the Fe content for the FILL-01 and PER-04 samples (Fig. 17). The Ca values for PER-04 and CAL-01 are probably due to the context of their formation (lake margin), which may have undergone several subaerial exposure episodes or post-depositional alterations. The higher values for Na in PER-01 and FILL-01 are coherent with the enclosing lithofacies, gypsiferous marls, and secondary gypsums, respectively.
The ternary Mg-Al-Fe graph (Fig. 18a) shows the PER-01 sample (within gypsiferous marls) as separated from the rest, which could indicate a greater burial depth. The rest of the enclosing rocks are either shallower (gypsum) or marginal (calcarenites and lutites). The ternary plot for Al-Ca-Fe content groups cherts according to their Fe content (FILL-01 and PER-04) and is consistent with the results observed from the magnetic properties (Fig. 18b).

NAA
A total of 30 element concentrations were measured by NAA, including trace, major, and rare-earth elements (Table SI1). Determining the content of diagenetically immobile rare-earth elements (REE) in chert provides information about the formation environment of the chert deposits (Murray et al. 1990;Murray 1994 and references therein). Important are especially the shale-normalized concentrations of Ce n , La n , and Nd n , all measurable with NAA. The so-called Ce-anomaly measures the shale-normalized concentrations of Ce against La and Nd [Ce n /Ce n * = 3 Ce n / (2 La n + Nd n ), sometimes also logarithmically (Piper & Bau 2013)]. Also in use is the correlative ratio La n /Ce n . Values of Ce n /Ce n * < ~ 0.8 and La n /Ce n > ~ 1 point to Ce depletion and indicate pelagic, open-ocean depositional environments, whereas values of Ce n /Ce n * > ~ 0.8 and La n /Ce n < ~ 1 are found for continental margin environments (Malyk-Selivanova et al. 1998). If an archaeological artifact differs in these ratios from a geological deposit, it certainly has a different origin. In Fig. 19, the Ce-anomaly is shown as a function of the ratio La n /Ce n . As shale normalization, the concentration values of the world shale average (WSA) as given in Piper & Bau 2013 are taken. The large 1σ uncertainties reflect the large measurement uncertainties (errors). As specified above, the high Ce-anomaly values as well as the low La n / Ce n ratios indicate a continental margin environment for the formation of the geological chert samples.
Plotting the major elements obtained by NAA highlights the following (Fig. 20): (1) a higher Na content is found in PER-01 compared with the rest of the cherts, emphasizing the absence of Na in the FILL-01 chert type, despite this having secondary gypsum as a bedding rock; (2) the PER-01, PER-03, and MDT cherts show a high K content; (3) the high Ca content of PER-03 is justified by its lower silicification, reflecting the presence of Ca from its host (calcarenites) and consistent with the presence of Ca phases according to X-ray diffraction, as also observed in the ICP-OES and p-EDXRF analyses; (4) as regards Fe content, differences are observed in PER-04 and FILL-01 with respect to the rest of the samples, as consistent with magnetic property analysis, ICP-OES and p-EDXRF.
NAA-REE plotting shows a great dispersion of values, the values being higher in PER-04 than in the rest of the samples (Fig. 21). The PCA (components 1-2; 92.6% of  the total variance) of the REE differentiates the PER-04 and FILL cherts from the rest of the samples according to their higher REE content (Fig. 22). Nevertheless, it seems that what we are observing represents different stages in the silicification process. At this point, it suggests that the composition of the bedrock does not play much of a role. Considering the geological map of the area (IGME 1975), the presence of REE corresponds to the Triassic and Paleozoic materials in the region. The Cenozoic materials seem to be riding on a reverse fault and to have been placed on top of the Paleozoic and Mesozoic materials, the former probably being of granitic origin.

Discussion
From a macroscopic point of view, it is difficult to associate the different SMT chert varieties with each of the lithofacies in the Sant Genís Formation. First of all, it is important to perform an exhaustive and systematic field work surveys. Secondly, the elaboration of columns and stratigraphic profiles, and the calculation of the chert ratio for each outcrop and formation. In this first stage, we were able partially to identify, quantify, and delimit the macroscopic variability of this key SMT chert type. In addition, the petrographic analysis of the geological samples allowed us microscopically to differentiate chert varieties of the SMT chert type by their quartz features, mineral components, and preserved sedimentary structures in accordance with the embedding lithology. With this background, we have been able to propose a multi-technique study to assess whether these varieties of the SMT chert type were indeed distinguishable through a variety of techniques of characterization and materials science. Considering the challenge posed by the problems associated with equating the macroscopic scale with the microscopic scale, the idea is to test whether the use of these techniques will lead to a greater and easier characterization of archaeological assemblages. On the other hand, the comparison of these techniques will allow us to evaluate their viability and which ones can be the most useful for archaeological purposes.

Multi-technique comparison
The nature of the enclosing deposits in the saline lakes of the Sant Genís Formation makes it difficult to differentiate among silicified materials. Our multi-technique approach provides with the insights to relate the geochemical composition with the enclosing lithofacies, their stratigraphic positions, and associated paleoenvironments within the geological formation.
Nevertheless, this work has identified similarities and differences among some of the chert varieties of the SMT chert type relating to the concentration of Fe and Ca. The XRD results clearly differentiate those chert varieties rich in calcite (PER-03) from others. These results are consistent with the p-EDXRF and with the nature of the enclosing rocks. The presence of a higher Ca concentration has to do with the preservation of the original embedding trait, calcarenites from the marginal environments of the evaporite basin, as pointed out by several authors for chert formations (Knauth 1994;Ortí et al. 2007Ortí et al. , 2014. The analysis of the magnetic properties has provided interesting data but being a relatively unfamiliar technique in sourcing lithics, a brief background information and its relevance is added to help the reader. Mineral-magnetic analyses mainly focus on the characterization of the composition (a specific type of mineral), concentration, and granulometry of ferromagnetic mineralogy. There are several types of ferromagnetism: ferro-, ferri-, and antiferromagnetism (Evans & Heller 2003). For the scope of this study, we will refer here to ferromagnetism (sensu lato, s.l.) in reference to materials with very large magnetic moments in applied fields (Dekkers 2007). Ferromagnetism is a property of a few minerals such as iron oxides (e.g., magnetite or hematite), iron hydroxides (e.g., goethite), and iron sulfides (e.g., greigite or pyrrhotite) capable of retaining a memory or remanence of the direction and/or intensity of a magnetic field even when this is no longer present. That is the basis of paleomagnetism. In other substances (diamagnetic and paramagnetic), a very low-magnetic moment is measured in applied magnetic fields and no remanent magnetic moment exists. Paramagnetic minerals such as siderite or pyrite yield small but positive magnetic moment that quickly outweighs diamagnetic moments. Diamagnetism is a property of all matter and exhibits small and negative magnetic moments under an applied field (Dekkers 2007). Minerals such as quartz or calcite, dominant in cherts including those studied here, are diamagnetic (Evans and Heller 2003;Dunlop and Özdemir 1997). Diamagnetism is hundreds of times smaller than paramagnetism and thousands of times smaller than ferromagnetism (s.l.). For that reason, subtle variations in the Fig. 21 Plot of the REE elements from the NAA ferromagnetic mineralogy of these cherts dominate the magnetic signal and such differences in their magnetic properties can be interpreted in terms of different catchment areas, diagenesis, or subaerial exposition, for instance.
The results of the magnetic property analysis obtained here are consistent with similar works, highlighting in general the weak magnetization of cherts (e.g., Borradaile et al. 1998;Thacker & Ellwood 2002;Larrasoaña et al. 2016;Ortega et al. 2018). However, among the geological samples studied, certain differences in the ferromagnetic mineral composition and concentration can be observed: magnetite is the dominant ferromagnetic mineral, but variable contributions of hematite are distinguished depending on the type of chert studied. A higher concentration of hematite is particularly visible in the FILL-01 and PER-04 chert types, whereas in the other cherts (CAL-01, MDT, PER-01, and PER-03), its contribution is very small or extremely difficult to distinguish by its magnetic properties alone. One of the main advantages of studying magnetic properties compared to other techniques is that instead of generically identifying iron (Fe) in the sample's composition, it allows detecting which specific ferromagnetic phase dominates the magnetization. This is the case of FILL-01 and PER-04 cherts where hematite was identified as the main iron oxide (Table 3 and Fig. 14), whereas, broadly speaking, p-EDXRF or ICP-OES only detects the element (Fe) per se. It is also interesting that these two varieties of chert (FILL-01 and PER-04) are precisely the ones showing the highest concentration of iron in p-EDXRF and ICP-OES analyses (Tables 4 and 5,  respectively). Basically, a good testimony of the complementarity of information that a multi-method comparative approach like this yield.
The homogeneous and high Fe values that grouped each two samples (FILL-01 and PER-04) are also consistent with the fieldwork data, as they are associated to gypsiferous and gray/red lutite lithofacies with a greater presence of oxide product of frequent subaerial exposures, and the results of thin-section (facies) analysis (see "Petrographic analysis"), p-EDXRF (see "p-EDXRF") and ICP-OES (see "ICP-OES").
The presence of Al detected for the FILL-01 and PER-04 samples by ICP-OES also confirms the occurrence of subaerial exposure episodes or post-genetic alterations, as Al is a major mineral component of exposure processes of gypsum and hydrated aluminum filosilicates conforming the enclosing rocks. When comparing the Mg-Al-Fe content (see Fig. 18a), the PER-01 and PER-03 cherts are separated from the rest, probably due to the deeper nature of its depositional environment and enclosing lithofacies (gypsiferous marls and calcarenites lithofacies), suffering lesser exposure processes. The rest of the grouped chert types share shallow (gypsums) or marginal (lutites) enclosing lithofacies as depositional environments of the Paleogene saline lake that produced the Sant Genís Formation (Ortí et al. 2007). Comparing the Al-Ca-Fe content, the cherts are grouped (FILL-01 and PER-04) according to their measured Fe, a result which is also consistent with the results observed on the basis of ferromagnetic properties and confirmed by the other applied techniques, showing higher Fe content. PER-03 is also separated showing a higher content of Ca due to their calcareous origin and deepest stratigraphical position with no subaerial exposure processes.
Relative fractionations of rare-earth elements (REE) through NAA analyses have proved useful environmental indicators in chert of all geological ages, determining tectonic history and the sediment's depositional position within ocean basins from the relative contribution of endmember terrigenous (non-depleted in Ce) and seawater (Cedepleted) sources (Murray 1994). In this work, REE analysis was applied to establish whether the so-called Ce-anomaly is valid for chert sources from inner platforms and continental environments. The results show that the Ce-anomaly in evaporitic cherts is always higher than 1 and the La n / Ce n ratio lower than 1, agreeing with the values from the continental margin environments noted by Murray (1994). Murray (1994) also pointed out that the presence of Ca, Mn, Mg, P, K, Na, Sr, Ba, Ge, Co, B, and V is in each case significantly affected by diagenetic processes and no longer reflects the depositional environment of the host sediment. PER-04 and FILL-01 are the only samples that present a major element marker with this technique, confirming also a higher Fe concentration (particularly hematite as shown by magnetic properties) ( Table SI1).
As mentioned in the results ("NAA"), and considering the geological data of the area (IGME 1975), the origin of the REE in cherts most probably is related to Triassic (Buntsandstein, Muschelkalk, and Keuper) and Paleozoic materials. These REE would be related to the different enclosing sediments and silicifications that have been incorporated throughout the time and diagenetic processes. The variability or similarities among samples likely depended on the composition of the older surrounding materials, as well as the sedimentary environment in which they were formed (Anadón et al. 1979;Gibbons & Moreno 2002).
As pointed out by Malyk-Selivanova et al. (1998) immobile elements detected by NAA analysis are excellent indicators of provenance and depositional environment because they can discriminate between different stratigraphic formations and units. As demonstrated, mobile elements, which can be detected using other techniques (such as EDXRF and ICP-OES), trace diagenetic environments by showing local variations. Therefore, the identification of mobile element concentrations (such as Fe and Ca) contributes to differentiating chert varieties within a geological formation. In our study, the continental evaporitic depositional environments of the Sant Genís cherts forced us to consider mobile elements as markers of the diagenetic processes (Bustillo et al. 2012). Thus, the presence of Fe and/or Ca in samples may suggest differences among the enclosing lithological units and, as a result, distinctions in the particular chert outcrop.

Consistency and challenges of the multi-technique approach
Three archaeological samples from the Abric Romaní site were chosen to assess the consistency/reliability of the geochemical analysis performed (Gómez de Soler et al. 2020a). Petrographically, these archaeological samples were assigned to various lithofacies from the analyzed formation: AR'02-M-U47-220 to PER-04; AR'02-M-P53-537 to CAL-01; and AR'01-M-Q47-49 to PER-01.
First, the Ce-anomaly obtained from these archaeological samples, as with the geological ones, shows a correlative ratio La n /Ce n where Ce n /Ce n * > ~ 0.8 and La n /Ce n < ~ 1, indicating similarities in their depositional environment (see Fig. 19).
XRD diffractograms indicate homogeneity among the geological and archaeological samples, with quartz as the main compositional phase, followed by moganite. However, variability in the results is detected by the p-EDXRF and in the content of trace elements and REE.
The PCA of the p-EDXRF (68.1% variance; components 1 and 2) (Fig. 23a) indicates similarities among the archaeological samples and the PER-01 chert type and MDT (56.8% variance; components 1 and 3) (Fig. 23b). The PCA of the NAA-REE results differs from the former petrographic ascriptions, suggesting the difficulty of ascribing earlydiagenetic chert types to formations such as our study case, i.e., Cenozoic evaporites with so much vertical and lateral variability (Fig. 24). These results highlight the challenge of accurately differentiating chert varieties formed within similar physicochemical conditions using various analytical techniques, underscoring the need for larger sample sizes to encompass the full geochemical variability of chert. At this point, depending on the scale of analysis, it may be possible to ascribe the archaeological samples to the geological formation, but not to specific lithofacies of that formation. One point in favor is that the archaeological application to establish the catchment areas and the mobility routes that these communities had is not affected. All the outcrops are part of the same geological formation and are delimited in a specific area to the south of the municipality of Sant Martí de Tous, being considered varieties of the SMT chert type. However, this would not be valid if this siliceous Fig. 23 a PCA of the p-EDXRF of the geological and archaeological samples with 68.1% variance (components 1 and 2); b PCA of the geological and archaeological samples with 56.8% variance (components 1 and 3). heterogeneity were to occur in a geological formation widely spread over a territory.

Insights of the multi-technique approach
What are the most adequate current techniques for the characterization of Cenozoic evaporitic chert? XRD is an easy, accessible, and economical technique, but as it identifies crystalline phases, for cherts, it generally represents quartz as the predominant crystalline phase. The non-destructive version of this technique, µ-XRD, only analyzes small surface areas (around 1mm 2 ) of the sample. This area can be increased by applying XY oscillation during the analysis. Magnetic property analysis is a fast and minimally invasive technique with more success in materials rich in ferromagnetic (s.l.) minerals like obsidian (Frahm et al. 2016), but it also has diagnostic potential for others lithologies including cherts (e.g., Larrasoaña et al. 2016;Ortega et al. 2018). Even though cherts are diamagnetic rocks (do not preserve remanent magnetization), the analyses carried out in this study have successfully detected compositional and concentration differences in the Fe-bearing ferromagnetic mineralogy. Specifically, hematite has been identified and quantified as the most present iron oxide in FILL-01 and PER-04 chert varieties along with minor contributions of magnetite also in the others. Identifying which specific ferromagnetic mineral is present is an information only provided in this study by the magnetic properties, regardless of other techniques that can also detect "iron" (broadly speaking). P-EDXRF has proved to be an excellent solution for characterization due to its non-destructive, accessible, and economical features (Gauthier et al. 2012;Roldán et al. 2015;Sánchez de la Torre et al. 2017b, 2020Eixea et al. 2021, among others). ICP-OES is a useful technique, as has been demonstrated by other studies (Ramacciotti et al. 2019(Ramacciotti et al. , 2022Agam 2020). However, it is a destructive technology, and difficulties in the digestion processes call its suitability into question, especially when ICP-MS detects a greater number of different elements and does so more precisely (Navazo et al. 2008;Finkel et al. 2022, among others). Yet this has the same problems as ICP-OES, being a destructive technique with difficulties in the process of digestion. Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) might be a more appropriate technique, as various works have previously shown (Andreeva et al. 2014;Speer 2014;Gurova et al. 2016;Sánchez de la Torre et al. 2017a, 2020Brandl et al. 2018;among others). NAA is a costly and destructive technique and not very versatile, as very few centers are equipped with the required facilities and technologies (Luedtke 1978; Malyk-Selivanova  Huckell et al. 2011;Speer 2016;Prudêncio et al. 2016;Eixea et al. 2022, among others). Beyond the most commonly used techniques for rock characterization, there are other complementary techniques used in different works that have given good results and can be explored in the future. Examples are visible/near-infrared reflectance (VNIR) spectroscopy (Hubbard et al. 2004;Parish 2011Parish , 2018, infrared reflectance microspectroscopy (IRRMS) (Hassler et al. 2013), Raman spectroscopy (Olivares et al. 2009;Capel Ferrón et al. 2015), trace element mapping with PIXE (Sánchez de la Torre et al. 2017a), and scanning electron microscopy with energy-dispersive X-ray spectroscopy (Abrunhosa et al. 2020).
Cenozoic evaporitic cherts are heterogeneous in nature, showing a high diversity of macroscopic, textural, geochemical, magnetic, and compositional traits that make it difficult to ascribe them precisely to different depositional environments and lithostratigraphic units (Colombo 1986;Ortí et al. 1997Ortí et al. , 2007Ortí et al. , 2014Soto et al. 2014Soto et al. , 2023Soto 2015;García-Simón & Domingo 2016;Gómez de Soler 2016;Gómez de Soler et al. 2020a). Facies analysis, systematic silicification samplings and samplings of enclosing rocks, and petrographic thin sections are shown to be traditional but appropriate techniques for differentiating chert varieties and different enclosing lithostratigraphic units.
Our analysis has shown that the individual application of modern compositional techniques may well fail to differentiate early-diagenetic cherts from a single geological formation, even from different lithological units. A combination of several techniques is to date the most adequate procedure for characterizing cherts and encompassing their variability. Exhaustive geoarchaeological prospections, facies recognition techniques, systematic sampling studies, macroscopic identification, petrographic thin-section analysis, and the application of several geochemical characterization techniques such as EDXRF and LC-ICP-MS would ensure proper characterization and geological ascriptions.
This study has highlighted the strengths and weaknesses of the different tests, evidencing the need for a variety of mineralogical and geochemical characterization techniques, as well as for a deep geological knowledge and exhaustive and systematic field recognition techniques to understand the full variability of prehistoric lithic raw materials. The differences and similarities shown by this study also bring to light the need for a systematic sampling of other regional formations yielding chert in order to establish a comparative framework. This would provide a starting point from which to consider the whole range of variability and the causes for the presence of the different trace elements and REE, allowing the different raw materials available to be quantitatively distinguished (Luedtke 1978;Malyk-Selivanova et al. 1998).

Implications of the SMT chert in the archaeological assemblages
Due to its high concentration in a relatively small area (calculated using the CAR method) and its heterogeneity (in terms of size, shape, color, texture, and knappability), the SMT chert seems to have been a lithic resource that was long exploited by Pleistocene and early Holocene groups inhabiting an area ranging up to more than 50 km from the Sant Genís Formation. Archaeological examples of its exploitation are the Middle Paleolithic occupations of the Abric Romaní site (Gómez de Soler et al. 2020a), the Teixoneres cave with the evaporitic cherts documented there (currently under study), the many Mesolithic and Neolithic sites documented around the town of Sant Martí de Tous and along the Anoia Basin (Enrich 1981;Clop et al. 2005;Cámara 2017;Pujol 2016;Gómez de Soler et al. 2021), and the archaeological sites of the Cinglera del Capelló (e.g., Can Manel, Abric Agut) in Capellades (Vaquero et al. 2013) (15 km east of Sant Martí de Tous), attesting the continuous use of this chert type in the central part of the Catalan region (NE Iberian Peninsula).
This region is delimited by other areas with abundant chert sources. These include the Gaià and Francolí Basins, with the evaporitic cherts of the Valldeperes Formation (Lutetian) and the Bosc d'en Borrás Formation (Bartonian) (Gómez de Soler 2016), and the lacustrine cherts of the Castelltallat Formation in the Pre-Pyrenees (Roy et al. 2013;Sánchez de la Torre et al. 2017a, b;Ortega et al. 2018). To the northeast and east, a territory almost empty of chert is detected until beyond the Pyrenees in the Languedoc-Roussillon region (France). The only presence of chert, albeit scarce, is restricted to the Jurassic unit of the Illes Medes, the Cretaceous Montgrí Formation (Ortega et al. 2016), and the Neogene opaline formation of the Camp dels Ninots (Gómez de Soler et al. 2012) and Puig de les Moleres (Vehí 2001). In studies of lithic procurement in these regions, no chert from these formations has ever been documented (Ortega & Maroto 2001;Ortega 2002).
The near-absence of chert in the northeast and east of the NE Iberian Peninsula and the abundance and high quality of the SMT chert compared to other available lithologies justify the diachronic selection of SMT as a key chert source for central Catalonia. Further systematization of the mineralogical and geochemical characterization of the archaeological assemblages will enable SMT chert exploitation to be identified in prehistoric settlements.

Conclusions
The SMT chert type is confirmed by its abundance at specific locations within a delimited territory in the NE Iberian Peninsula. The region is shown to be an important raw material procurement area, as evidenced by the presence of the SMT chert type in the main prehistoric sites of the region (e.g., Abric Romaní, Cova de les Texioneres) and the constant discovery of new sites in the region around the Sant Genís Formation, especially from the Neolithic period onwards (e.g., Cal Sitjo, La Guinardera workshop, La Guinardera Nord workshop).
Through the characterization of cherts from the Sant Genís Formation, a Cenozoic evaporite formation, systematically sampled vertically and horizontally, this study confirms the great heterogeneity of Cenozoic cherts from the Sant Genís Formation (SMT chert type), widely exploited during the Pleistocene and Holocene of the Iberian Peninsula. This heterogeneity hampers the macroscopic characterization of the archaeological cherts, but petrographic analysis has enabled us to identify their origins and even, to a lesser extent, their lithostratigraphic unit of provenance within the formation. Complementary mineralogical and geochemical techniques applied in this study show that the application of specific criteria such as the concentration of Fe and Ca can assist in the identification of these cherts in one of the chert varieties of the SMT chert type, and consequently ascribe the procurement of these cherts to one of the outcrops of the Sant Genís Formation (SMT chert type).
As archaeologists, we are used to facing the sempiternal problem of correlating macroscopic and microscopic analytical scales, but the simultaneous use of multiple techniques can minimize the frequent inconsistencies of qualitative analysis. Intensive geoarchaeological survey programs, the identification of possible differences among the chert outcrops within a specific geological formation, exhaustive facies analysis, and sampling for reference collections, followed by macroscopic and petrographic thinsection characterization, combined with mineralogical and geochemical techniques, stand out as the ideal and most adequate practices for encompassing the variability of geological chert and thus for the correct and precise characterization of its provenance.