Multi-Analytical Techniques to Deﬁne the Mineralogical and Petrophysical Characteristics and Provenance of Siliceous Lithic Findings: The Case Study of La Calvera Rock Shelter (Cantabria, Spain)

: This archaeometric study aims at characterizing the archaeological ﬁnds belonging to the lithic industry from La Calvera rock shelter (Camaleño, Cantabria) and at hypothesizing the possible provenance of each material. The site, located in the mountainous area of Picos de Europa National Park (more than 1000 m a.s.l.) close to the megalithic complex of Peña Oviedo, is characterized by the presence of hearths and charcoal remains, ancient pottery, and a rich lithic assemblage composed of siliceous rocks dating back to >8000 BP and linked to the ﬁrst Holocene occupations of the Cantabrian Mountains. For the study of the rock shelter’s lithic assemblage, a multi-analytical approach was used: SEM-EDS and XRD analyses were performed to deﬁne the microtextural characteristics of samples and to identify the amorphous/crystalline phases; physical and mechanical tests were conducted to deﬁne the petrophysical properties (density, porosity, imbibition and saturation indexes, mechanical strength) of the samples. Preliminary data of chemical compositions, obtained by portable XRF, are also presented. The results show that the different siliceous materials (quartzite, cherts, hyaline quartz) can be distinguished through the basic analytical techniques used here. In addition, most of the archaeological samples have mineralogical and petrographic features similar to the natural samples coming from nearby outcrops, corroborating the hypothesis of a local material supply. The presence of local sources of useful raw materials could have favoured the site’s occupation. Finally, the diverse compositional and textural feature of the analysed materials result in different mechanical properties (porosity, density, hardness, workability), so they likely had different uses and technical functions.


Introduction and Aims
A rock shelter is an opening of modest size and extension that is generally formed by weathering and erosion processes (e.g., induced by water run-off) of a rock that is less resistant than the surrounding rocks.Shelters are horizontally shallow, unlike caves generated by karstification, which are much deeper.Rock shelters can have an archaeological importance since they were often used by humans as refuge from the weather [1][2][3][4].Prehistoric people frequented such shelters as a place to live, leaving behind tools and other artefacts that now assume a high archaeological significance [5][6][7][8].
La Calvera rock shelter is located in Picos de Europa National Park, 1180 m above sea level, close to Camaleño (Cantabria; Figure 1), in the area of the Peña de Oviedo megalithic complex.From the geologic point of view, the Picos de Europa is an imposing mountain  During the excavations in the La Calvera rock shelter, the archaeologists documented the presence of hearths and charcoal, historical pottery in the upper levels, and the presence of a lithic assemblage composed of different siliceous rocks (chert, rock crystal, quartzite) Minerals 2023, 13, 666 3 of 22 dating back to >8000 BP and linked to the first Holocene occupations of the Cantabrian Mountains [9].More than 500 lithic fragments were discovered with some retouched pieces (mainly cores, scrapers and small blades).
Chert is a very fine-grained siliceous rock composed almost exclusively of microcrystalline quartz and chalcedony.It is normally found in the form of layers or as nodules and lenses within carbonate rocks or alternated with clay-rich shale levels.Stratified cherts are essentially of biogenic origin and owe their formation to the accumulation of the siliceous shell or skeleton organisms (e.g., radiolarians, diatoms, silico-flagellates and sponges) [10], even if recent studies suggest a significant role of diagenesis in silica redistribution in bedded cherts [11].Nodular cherts generally have a diagenetic origin and would form as a result of silicification processes inside the host rocks [10].The contribution of hydrothermal fluids interacting with seawater in the depositional environment has been reported locally by several authors [12,13].Given their high hardness, chemical and physical resistance, conchoidal fracturing and its use as fire starter, this material had an important use in the ancient lithic industry, especially in prehistoric times.On the basis of the different chipping techniques of cherts, subperiods and, thus, the working material culture of civilizations can be identified.
Contrary to chert, the term quartzite is more ambiguous since historically it has been used with different meanings to describe a wide variety of rocks of both metamorphic and sedimentary origin [14].Moreover, distinguishing among different types of quartzite is often not easy due to the similar mineralogy and textures that different quartzites can share.As recently pointed out by Prieto et al. [15] the full characterization of this material should be achieved through a detailed petrographic analysis integrated with digital imaging to recognize and quantify the textural features of quartz grains; however, this kind of approach is necessarily destructive and, thus, cannot be applied to highly valuable archaeological finds.In this paper, the term quartzite is used to refer to silicified quartzarenites (also known as orthoquartzite) almost completely composed of quartz.The discrimination between quartzite and cherts has been performed through macroscopic and microscopic observation (under reflected light), following the criteria proposed by [14] and concerning the fracturing, the luster, the grain size of the samples.
The provenance and lithic supply in Cultural Heritage studies [16,17] are fundamental to define the dynamics involving human populations and their surrounding environment, mobility, and possible relations among different settlements; however, matching prehistoric artefacts and raw materials from potential quarries just on the basis of naked-eye examination is often controversial.Only a multidisciplinary approach, based on the use of different analytical techniques, allows for the full characterization of the archaeological geomaterials and their natural counterparts [17][18][19][20][21][22], leading to more objective and robust results [23][24][25].
The study of the lithic assemblage from La Calvera rock shelter, has been performed through a two-stage, multi-analytical approach.Firstly, a large set of samples was analysed using non-destructive to micro-destructive techniques, based on smartphone imaging and molecular spectroscopy techniques [26], to understand if it is possible to distinguish siliceous materials coming from different source areas and finally to identify the supply area of archaeological finds.Secondly, according to the first analyses, a subset of samples was selected to carry out a petrographic, mineralogical and petrophysical characterization (this study).This second stage is aimed at further characterizing the selected samples and to test the hypothesis postulated after the former analyses.To achieve this goal, the microstructural and textural features of both the amorphous and crystalline phases within the samples were observed through a scanning electron microscope (SEM-EDS) and analysed by powder X-ray diffraction (XRD) to detect minor phases and the degree of crystallinity.Finally, the physical and mechanical characteristics were determined to understand the technical properties of the artefacts in relation to their possible use as tools.A summary of the chemical composition of the different lithotypes is also reported.

Geological Setting and Natural Chert Occurrences
The Cordillera Cantabrica extends for about 480 km along the northern Spain coast, from Galicia (to the west) to Basque Country (to the east).The Picos de Europa area, where the La Calvera rock shelter is located (close to the Camaleño area, Figure 1), is in the central part of the Cordillera.The current relief of the Cantabrian zone is the result of the superposition of the Variscan and Alpine orogeneses during the Palaeozoic and Cainozoic Eras, respectively [27,28].The orogenic cycles led to the stacking of different tectonic units that mostly share the same sequence of lithostratigraphic formations, even if a local variability is observed.Indeed, from one unit to another, the same formation could have, or could even lack, different features and different thickness.Additionally, the name of the same formation can change depending on the unit or on the locality in which it occurs [29].
The area of Peña de Oviedo, where the La Calvera shelter rests, is located on the south-eastern side of the Picos de Europa park, close to the Camaleño village (Figure 1).In this area, only the Picos de Europa Unit (northwest) and the Pisuerga-Carrión Unit (south-east) were cropped out, separated by a regional thrust.The former is characterized by a predominance of upper Carboniferous limestone [30] belonging to the formations of Picos de Europa and Barcaliente (previously known as Caliza de montaña).Small slices of limestones belonging to the Las Portillas and Alba (or Genicera) Fm., also occur in the southern side of the Barcaliente Fm.The Pisuerga-Carrión Unit consists of several formations (Barcena, Viorna, Narova, Remoña and Campollo Fms.), mainly siliciclastic (slates, shales, sandstone and conglomerates) with limestone intercalation, of the Carboniferous age [31].A sequence of allochthonous Devonian to Carboniferous rocks covers large areas of the Pisuerga-Carrión Unit; this sequence, referred to as Palentine Nappes or Palentian domain [32] consists of alternating siliciclastic and calcareous deposits grouped in different formations.From the bottom to the top, they are the Gustalapiedra Cardaño Fm., Murcia Fm, Vidrieros Fm., Vegamian Fm., Alba Fm., and finally the Carboniferous sequence of the Potes group.
Quaternary covers consist of moraines, glacial/fluvio-glacial deposits, cemented landslides and boulder flows from the Pleistocene, and torrential cones, slope deposits and alluvium from the Holocene [33,34].It was during this period, dominated by important glacial systems of the Cantabrian Mountains and interglacial stages [35][36][37], that the Picos de Europa area assumed its current physiognomy and the typical rock shelters developed, mainly in the calcareous lithologies.
Several authors investigated the outcropping formations in the Cantabrian region to identify the potential sources of the lithic artefacts found in prehistoric human settlements [29,38,39].In particular, Herrero-Alonso et al. [29] compiled an updated and complete inventory of the chert-bearing formations in a wide area of the Cantabrian region that includes the study area of this paper.According to these authors, in the surroundings of Peña de Oviedo, cherts, radiolarites and/or quartzite of knappable quality can be found in the Las Portillas, Barcaliente and Picos de Europa formations (all belonging to the Picos de Europa unit) and in the Vegamián formation (Palentine Nappe).The Alba (or Genicera) formation, also hosting knappable cherts, is found in both units forming decametre-thick slices.Quartzite occurrences in the surroundings of La Calvera have been reported in the Murcia Fm. [40] and in the Vidrieros Fm. [41], associated with chert nodules.None of the cited quartzites have a metamorphic origin, rather, they are the result of sedimentary/diagenetic processes.Chert-bearing formations are documented also in some areas of the Pisuerga-Carrión Unit, but they do not occur in the study area.The features of these cherts and the comparison with the archaeological finds of La Calvera will be outlined in the Section 5.

Materials and Methods
Based on the spectroscopic analyses and colour features on a large sample set, ten representative samples were selected and analysed in this work (Table 1, Figure 2).One cobble of macrocrystal quartz or rock crystal (AR10), three samples of local grey chert (flakes: AR12.1 and AR18.2, chunk: RM03.4), one sample of ochre chert (AR21), two samples of quartzite (AR29 and AR34) and one sample of yellowish chert (AR37).Finally, for comparison with macroscopically similar chert, two samples of Domeño chert collected from the geological outcrop of Andilla (Valencia, Spain) were also analysed.Each sample was cleaned with deionized water and a brush in order to remove superficial contaminations and incrustations prior to the analyses.
A small chip of each sample was powdered to carry out XRD analyses, which were performed at the Department of Chemical and Geological Sciences (University of Cagliari).XRD patterns were acquired by a PANalytical X'Pert Pro diffractometer (Malvern PANalytical, Almelo, The Netherlands) that works with theta-theta geometry using Ni-filtered Cu Kα1 radiation (λ = 1.540598Å) and equipped with a X'Celerator detector.Operative conditions were acquisition range 5-70 • , step size 0.008 • , 0.19 s per step, voltage 40 kV and current 40 mA.Data were processed by X'Pert HighScore Plus (TM) 2.1.2software using the PDF2 database (released in 2010 by ICDD, Newtown Square, PA, USA).
SEM analyses were performed on ten samples using a Quanta Fei 200 equipped with a ThermoFischer Ultradry EDS detector, at the CeSAR laboratories (University of Cagliari, Italy).Raw samples were put into the sample chamber without conductive coating to preserve them for further analyses, thus, low-vacuum conditions (0.3 to 0.5 torr) were used to dissipate electrons from incident beam.Variable spot sizes of 4-5 (in arbitrary units given by the Quanta Fei equipment) and an accelerating voltage of 15-25 kV were adopted during the analytical sessions.
Petrophysical and mechanical tests on the samples were performed on ten specimens of the most significant samples.Physical tests were carried out according to Buosi et al. [42] and Columbu et al. [43,44] (see the Supplementary SM1).
Multielement analysis was performed on a larger sample set, including ten chunks of grey chert used as raw materials, twenty-four fragments of the same grey chert, six cherts of different types, fifteen fragments of quartzite, six rock crystal fragments and eleven fragments of Domeño chert.Elemental concentrations of Al, K, Ca, Fe, Ti and Zr were detected using a S1 Titan portable energy dispersive X-ray fluorescence spectrometer (pXRF) equipped with a Rh X-ray tube (50 kV) and X-Flash ® silicon drift detector (Bruker, Billerica, MA, USA).Internal calibration Geochem-trace was used.Each sample was analysed between two and up to five spots, and the results were then averaged.
Based on the spectroscopic analyses and colour features on a large sample set, ten representative samples were selected and analysed in this work (Table 1, Figure 2).One cobble of macrocrystal quartz or rock crystal (AR10), three samples of local grey chert (flakes: AR12.1 and AR18.2, chunk: RM03.4), one sample of ochre chert (AR21), two samples of quartzite (AR29 and AR34) and one sample of yellowish chert (AR37).Finally, for comparison with macroscopically similar chert, two samples of Domeño chert collected from the geological outcrop of Andilla (Valencia, Spain) were also analysed.

Compositional Characteristics
The characterization of mineralogical and amorphous (or with low crystallinity) phases and microtextural features was performed on 10 major samples considered most significant of the facies studied.After studying the macro-mesoscopic characteristics by optical microscopy analysis in reflected light (Table 2), the samples were analysed using scanning electron microscopy (SEM-EDS) and X-ray diffractometry (XRD).The results, reported in the next two paragraphs, allow for a preliminary comparison between the archaeological and geological samples in terms of characteristic compositional aspects aimed at the study of material provenance from the rock outcrops.

SEM-EDS Analysis
The results are reported in Table 3.Back scattered electrons (BSE) imaging at low magnifications (about 100×) shows that most samples look similar, with a dark-grey surface disseminated with light-grey spots, particles, and/or undefined stains (Figure 3a,b).Higher magnifications (400-800×) reveal that the surface is commonly quite irregular due to the presence of micrometre-sized cracks, holes and euhedral microcrystals (Figure 3c,d).Exceptions to this general appearance are the samples AR10 (rock crystal, Figure 3e), and AR34, whose surface is more compact and homogeneous.fragments and REE-phosphates were locally observed (Figures 3a,d and 4a).Characteristic snowflake-like patinas of NaCl have been observed only in sample AR10 (Figure 4f).Looking at Table 3, which reports a summary of the phases detected by the SEM analyses, the discrimination between the different categories of lithic fragments appears quite hard.Generally, a higher abundance of metal oxides, especially Fe, Ti and Zn, has Energy dispersive spectroscopy (EDS) analyses, even if just qualitative, provide evidence of the different phases in each sample.The most common phases (other than quartz) are represented by Si-Al-bearing (±K, Ca, Na, Fe, Mg) minerals in form of finegrained alteration patinas (Figure 3c,d) or having a phyllosilicate-like, tabular habitus (Figure 3f).Fe and/or Ti oxides are also widespread and were found in all samples in variable amounts and occurrences (isolated crystals, aggregates, oriented trails, patinas, etc.) (Figure 3a,d and 4a,b).Other minor phases are calcite (Figure 4c,d), whose former presence is testified also by euhedral holes (Figure 4d,e), gypsum (AR10, AR18.2, not shown) and pyrite in framboidal aggregates (AN1.2) (Figure 4c).Traces of Ni-Cr oxides, zircon fragments and REE-phosphates were locally observed (Figure 3a,d and Figure 4a).Characteristic snowflake-like patinas of NaCl have been observed only in sample AR10 (Figure 4f).been observed in greyish cherts and in RM3.4 whose surface was largely covered by r dish oxides encrustation (Figure 4b).Oxides seems to be less abundant in quartzite sa ples, even in the yellowish samples.The two samples from Andilla have very diffuse cite crystals likely linked to the geology of the surrounding rocks.Despite these dif ences, samples belonging to different categories can host the same phases, and sample the same type can have different phases; in other words, there are no minerals diagno of a given type.This is likely due to the samples' heterogeneity but also to the fact t SEM imaging investigates small areas, so it cannot be considered statistically represen tive of the whole sample.Looking at Table 3, which reports a summary of the phases detected by the SEM analyses, the discrimination between the different categories of lithic fragments appears quite hard.Generally, a higher abundance of metal oxides, especially Fe, Ti and Zn, has been observed in greyish cherts and in RM3.4 whose surface was largely covered by reddish oxides encrustation (Figure 4b).Oxides seems to be less abundant in quartzite samples, even in the yellowish samples.The two samples from Andilla have very diffuse calcite crystals likely linked to the geology of the surrounding rocks.Despite these differences, samples belonging to different categories can host the same phases, and samples of the same type can have different phases; in other words, there are no minerals diagnostic of a given type.This is likely due to the samples' heterogeneity but also to the fact that SEM imaging investigates small areas, so it cannot be considered statistically representative of the whole sample.

XRD Analysis
The mineralogical composition and a synoptic scheme of all XRD patterns are reported in Table 4 and in Figure 5, respectively.To highlight the minor peaks, the Y-axes (counts) have been reported in square root.Furthermore, to facilitate the pattern comparison, all Y-axes were set to a maximum value of 70,000 counts.All samples exhibit the typical patterns of well crystalline quartz that, using the "search and match" tool of X'Pert HighScore Plus software, provide a score of 75%-85% with synthetic quartz (Ref.code 00-046-1045 in PDF database).Despite the similarity of all samples, a difference can be observed: samples AR10, AR29 and AR34 show a higher crystallinity compared to the other samples.Indeed, they show higher intensities and lower FWHM (full width at half maximum) of all peaks, and those that are more evident are positioned at about 20.87 • , 50.13 • and 59.95 • 2θ angles.This means that quartzite samples are almost as crystalline as the pure rock crystal.All samples exhibit the typical patterns of well crystalline quartz that, using the "search and match" tool of X'Pert HighScore Plus software, provide a score of 75%-85% with synthetic quartz (Ref.code 00-046-1045 in PDF database).Despite the similarity of all samples, a difference can be observed: samples AR10, AR29 and AR34 show a higher crystallinity compared to the other samples.Indeed, they show higher intensities and lower The measurement of the crystallinity degree (reported in Table 4) has been performed comparing the background signals of each sample with the background of a reference material of known crystallinity.The AR10 crystal rock was chosen as a reference material of almost perfect crystallinity, assigning it the value of 100%.The results of this calculation method should be carefully considered since it works when comparing similar samples (i) whose only difference is in the crystallinity degree and (ii) that are measured using the same operative conditions.In this case, the second requirement is met, whereas the first is only partially met since samples slightly differ from each other due to the presence of minor phases.To lessen the effect produced by other phases, only the region of the XRD pattern ranging from 19 • to 23 • 2θ, where the {100} quartz peak occurs at 20.86 • 2θ, was considered.This calculation confirms the higher crystallinity of quartzite samples AR29 and AR34 (99% and 95%, respectively, Table 4).Greyish AR12.1 and AR18.2 and yellowish AR37 cherts share a similar crystallinity of 84%, 88% and 80%, respectively.Andilla cherts have crystallinity of 84%, so they are highly comparable with local and archaeological cherts.Significantly lower values are found in the grey-ochre chert AR21 = 62% (Table 4).
The matrix crystallinity shows a general positive correlation (R 2 = 0.68, Figure 6) with the particle grain size obtained by SEM analysis.We can observe three populations of samples with different behaviours due to the diverse petrogenesis of these siliceous rocks, especially regarding the crystal rock and the quartzite samples with respect to the cherts.3).
The identification of minor phases is quite difficult since the intense si peaks hides the contribution of minor phases, and also because the latter very low amounts.Thus, in many cases, only the main peak of the minor observed, and it is not possible to assess the presence of a phase with certain or two peaks.However, by comparing the patterns of all samples and by c pieces of information with those of the SEM analyses, it is possible to identif with good reliability, whereas some other can be just hypothesized.Calcite tified in sample AN1.2 by the detection of its six higher peaks; the most i 29.40° 2θ, has been also found in AN5.1, AR21 and RM3.4 and, even if it is ascribable to calcite in these samples, its presence is reliable since it is sup SEM-EDS analyses.The presence of a ferroan dolomite/ankerite can be AR18.2 and RM3.4, which show a peak at 30.94° 2θ angle, but this is the on this mineral; therefore, it is too weak to clearly assess its presence.Pyrite h  3).
The identification of minor phases is quite difficult since the intense signal of quartz peaks hides the contribution of minor phases, and also because the latter are present in very low amounts.Thus, in many cases, only the main peak of the minor phases can be observed, and it is not possible to assess the presence of a phase with certainty just by one or two peaks.However, by comparing the patterns of all samples and by coupling these pieces of information with those of the SEM analyses, it is possible to identify some phases with good reliability, whereas some other can be just hypothesized.Calcite is clearly identified in sample AN1.2 by the detection of its six higher peaks; the most intense one, at 29.40 • 2θ, has been also found in AN5.1, AR21 and RM3.4 and, even if it is the only peak ascribable to calcite in these samples, its presence is reliable since it is supported by the SEM-EDS analyses.The presence of a ferroan dolomite/ankerite can be supposed in AR18.2 and RM3.4, which show a peak at 30.94 • 2θ angle, but this is the only evidence of this mineral; therefore, it is too weak to clearly assess its presence.Pyrite has been identified by five peaks in AR12.1, and three of them have been found in AR18.2 and RM3.4.Interestingly, the only sample where pyrite was identified by SEM analysis (AN1.2) does not show peaks referrable to it.An undefined phyllosilicate phase (illite/muscovite) has been recognized in AR12.1, AR18.2 and AR34 through the identification of very small peaks at about 8.8 • , 19.8 • and 34.9 • 2θ angles; this is a tentative association since the maximum peak of these phyllosilicates (at about 26.6 • ) is hidden by the strongest quartz peak; however, the finding of tabular Si-Al-phases by SEM-EDS suggests the correctness of this assumption.Another plausible phase, identified by XRD analyses, is a Ti oxide that could be either anatase or rutile; these minerals share a similar position of the strongest peak at about 25.3 • 2θ but this is the only visible peak (in AN5.1, AR12.1 and AR18.2) that can be assigned to this phase.However, the presence of a Ti-oxide, as well as the phyllosilicate phase, is supported by SEM-EDS observation and analyses (Figures 4 and 5; Table 3).The sample AR34 does not show the anatase/brookite peak but is the only sample that has an unassigned peak at 27.4 • 2θ and that is compatible with the strongest peak of rutile (another Ti oxide polymorph).Different unassigned peaks are found in all patterns the most noticeable of which are those at 25.3 • and 29.0 • 2θ that recur in almost all patterns.Other unassigned peaks are found locally at 25.8 • 2θ (in four samples) and at 56.3 • 2θ only in AR34.Despite its ubiquitous presence, Fe oxides were not detected by XRD, probably due to their low crystallinity.Likewise, the other accessory minerals identified by SEM (zircon, P-REE, gypsum, etc.) were not detected due to their very low content.

Petrophysical Properties
For a complete characterization of the analysed archaeological and geological samples, and to understand the technical properties of the artefacts, the following main petrophysical properties were determined: real density, bulk density, open porosity to helium, open and closed porosity to water, imbibition coefficient (expressed in weight), saturation index, punching resistance index, and the compressive and tensile strengths calculated indirectly (Table 5).The physical and mechanical properties are variously affected by different compositional features of the samples.Figure 7 shows the graphs that plot the data (Table 5) of the main physical properties that have been determined.In Figure 7a, which shows the open porosity to helium against the bulk density, a high variability (Table 5) of the open porosity (i.e., interconnected) is observed, which is mainly induced by the compositional heterogeneity and only subordinately by the decay of the samples.index, punching resistance index, and the compressive and tensile strengths calculated indirectly (Table 5).The physical and mechanical properties are variously affected by different compositional features of the samples.Figure 7 shows the graphs that plot the data (Table 5) of the main physical properties that have been determined.In Figure 7a, which shows the open porosity to helium against the bulk density, a high variability (Table 5) of the open porosity (i.e., interconnected) is observed, which is mainly induced by the compositional heterogeneity and only subordinately by the decay of the samples.The real density, related to the characteristics of the solid phases (mineralogical and amorphous), reflects the compositional differences.The samples of Domeño cherts (AN01.2,AN05.1 from the Andilla area) and grey-ochre chert (AR21), which have similar values ranging from 2.52 to 2.56 g/cm 3 for real density and from 2.38 to 2.41 g/cm 3 for bulk density, belong to the first homogeneous group (Figure 7a), while the remaining samples (i.e., AR12.1, AR18.2, AR37, AR10, AR29, AR34, RM03.4) constitute a second homogeneous group, in which sample AR37 also falls, which, however, shows intermediate behaviour between the two groups (Figure 7a).Given that the real density of the amorphous phases is lower (about 2.3-2.55 g/cm 3 ) than the crystalline phases of α-quartz (2.65 g/cm 3 ), a light positive correlation (R 2 = 0.52) between the real density and the degree of matrix crystallinity is observed (Figure 7b).
The open porosity to helium, ranging from 5.1% to 5.9% in Group 1 and from 1.2 to 3.2% in Group 2 (Table 5), affects the bulk density, as shown in Figure 7c, with a high correlation coefficient (R 2 = 0.96).The porosity of these samples is mainly represented by intraphase pores with planar geometry; thus, it shows a positive correlation with the size of the amorphous/crystalline phases (Figure 7d) with a high correlation coefficient (R 2 = 0.98).The water open porosity shows a trend similar to that of helium open porosity.Additionally, it shows lower values ranging from 2.5% to 5.6% in Group 1 and from 1.0 to 3.1% in Group 2 (Table 5, Figure 7e).The saturation index (graphically shown in Figure 7e) shows values that are always lower than 100% and generally range from 85% to 98% (Table 5), except two samples (i.e., Domeño chert AN01.2 and raw material RM03.4) belonging to the first and second sample groups, respectively, which have lower values (i.e., 49% and 53%, Table 5).
Regarding the physical and mechanical properties, the samples show a great dispersion of data with high standard deviations (Table 5).The first group shows values of the PLT punching strength index from 8 to 15.2 MPa, while the Group 2 shows values ranging from 4.9 to 15.3 MPa (Table 5, Figure 7d).
Contrary to what is commonly reported in the literature, the strength index does not show a clear negative correlation with porosity (Figure 7f), probably due to the very low porosity values with small differences between samples.Only within Group 1 does a clear link between the two inversely proportional properties become evident.
In addition to porosity, the mechanical strength is also variously and subordinately influenced by the grain size of the phases (see Figure 8a), the degree of crystallinity of the matrix (see Figure 8b), and in some cases, by the presence of anisotropy (i.e., samples AR12.1, RN3.4) or crystal (i.e., sample AR10) planes that weaken the microstructure from a physical and mechanical point of view.Contrary to what is commonly reported in the literature, the strength index does show a clear negative correlation with porosity (Figure 7f), probably due to the very porosity values with small differences between samples.Only within Group 1 does a c link between the two inversely proportional properties become evident.
In addition to porosity, the mechanical strength is also variously and subordina influenced by the grain size of the phases (see Figure 8a), the degree of crystallinity of matrix (see Figure 8b), and in some cases, by the presence of anisotropy (i.e., samp AR12.1, RN3.4) or crystal (i.e., sample AR10) planes that weaken the microstructure fr a physical and mechanical point of view.

Preliminary Results of Multielement Analysis
The average elemental concentrations and standard deviation for the different cla of samples are shown in Table 6.

Preliminary Results of Multielement Analysis
The average elemental concentrations and standard deviation for the different classes of samples are shown in Table 6.As can be observed in Figure 9, grey chert flakes and chunks have similar levels for all the analysed elements, suggesting that the former were probably obtained from the latter as raw material.On the other hand, the other chert types show lower concentrations for Al, K, Ti, Fe and Zr.Domeño chert has the highest amounts of Ca.Average concentrations of quartzite are higher than those of cherts for most elements.Quartzite samples also have a higher variance.

Discussion
The investigations performed on the archaeological and geological samples yielded significant results regarding the sourcing, use and compositional characteristics of siliceous findings and raw materials coming from the area surrounding the La Calvera rock shelter archaeological site.
First, the petrographic and mineralogical investigations have allowed for the classification of the studied archaeological samples, mainly as cherts and subordinately as quartzites and "rock-crystals".The mineralogical analyses by XRD on the lithic finds have revealed both higher crystallinity and grain size in the quartzite and crystal rock samples

Discussion
The investigations performed on the archaeological and geological samples yielded significant results regarding the sourcing, use and compositional characteristics of siliceous findings and raw materials coming from the area surrounding the La Calvera rock shelter archaeological site.
First, the petrographic and mineralogical investigations have allowed for the classification of the studied archaeological samples, mainly as cherts and subordinately as quartzites and "rock-crystals".The mineralogical analyses by XRD on the lithic finds have revealed both higher crystallinity and grain size in the quartzite and crystal rock samples in comparison with cherts.The different crystallinity and grain size have also been confirmed by the petrophysical analysis.In detail: the real density to which it is positively correlated, while the porosity of rocks shows a negative correlation, confirms the presence of intraparticle planar pores among the particles (grains).The physical tests highlight the presence of two main different behaviours of the samples: the first population consisting of the archaeological cherts, quartzites and chunks of raw materials coming from the Calvera rock shelter, and a second subordinate sample population consisting of Domeño cherts from Andilla and the unclassified grey-ochre chert.Mechanical strength is high and comparable in the different samples analysed, due to the low overall porosity of the materials, which does not exceed 5%; however, the lower grain size and crystallinity of the chert microstructure, characterised by a more "brittle" physical and mechanical behaviour with pseudo-conchoidal microfracturing, probably facilitates the processing and production of sharper edges than quartzites and crystal rocks.The lithic assemblage of La Calvera shows that chert artefacts are commonly smaller than quartzite ones, likely due to different uses, which, in turn, are influenced by different petrophysical properties.The higher flakeability of cherts resulted in its suitability for arrowheads and cutting tools, whereas the coarser, less porous (<3%), and denser (>2.6 g/cm3) quartzites produced thicker and heavier tools, such as pestles or scrapers.
The measurement of the physical properties performed in this study was also aimed at testing the effectiveness of this approach in distinguishing among the different materials and, in future studies, using this information to identify potential sources.The data obtained seem to be inconclusive, but it must be considered that they are just preliminary results and that the method needs a wider sample set to be tested and statistically validated.From the available data, it can be observed that the foreign samples of Domeño, introduced to test the method, show a distinct behaviour in all diagrams of Figure 7 and can be grouped together with ochre cherts.Yellowish cherts have physical features between those of the Domeño and La Calvera materials.Considering that the analysed samples display detectable differences in their petrophysical properties, even if similar from a mineralogical point of view, we believe this approach is worth being tested further.
SEM-EDS allowed for the definition of the microstructural and textural characteristics of samples, highlighting specific discriminant compositional information on single samples by the identification of minor or accessory phases, e.g., the relatively high content of calcite, distinctive of Andilla Domeño cherts and also confirmed by XRD results.However, a clear signature of each specific source (or material) cannot be found by SEM-EDS and/or XRD analyses since it would require detailed partially destructive investigation on every lithic fragment.
Based on the above-described analyses, some provenance hypothesis can be made.
The excavation campaigns at La Calvera rock shelter discovered abundant chunks and cores of grey chert (samples AN01.2,AN05.1, A12.1, AR18.2) from different archaeological levels, making this lithic material the most representative of the site.Considering these abundant findings and the proximity of La Calvera to chert-bearing rocks [29,41], a local supply of this grey chert can be reasonably hypothesized.
Grey and black cherts are also commonly found within the Cantabrian Range but only the Alba Fm. (Lower Carboniferous), Barcaliente Fm. (previously known as Caliza de Montaña, Mountain Limestone, Upper Carboniferous) and Picos de Europa Fm. (Upper Carboniferous) are locally present [29].In addition, a report of greyish chert associated with dark quartzites is described by Castillo-Diez [41] within the Vidrieros Fm., in close proximity to the shelter of La Calvera, and in contact with the quartzites of the Murcia Fm.Although it is difficult to distinguish between the aforementioned cherts based on of the obtained results, the presence of some mineralogical peculiarities detected by XRD and SEM-EDS (i.e., presence of phyllosilicates, massive texture, absence of fossils, negligible amounts of Ca and Mg), would lean toward the Barcaliente Fm. cherts [29] or to the Vidrieros/Murcia Fms.The latter can be identified just by its proximity and by the macroscopic description since no analytical data can be found in the literature.On the other hand, we cannot exclude that further outcrops of the same chert, nor that cherts and quartzite boulders of secondary origin (i.e., from colluvial or alluvial deposits), were exploited [40].
As regards the yellow chert (AR37), petrographic and mineralogical data corroborate the first identification with a Piloña chert outcropping, based on chemical and macroscopic features, in eastern Asturias about 50 km NW of the studied site [38,45].This chert circulated in the area and is present in other Mesolithic sites of the area [46].The performed analyses pointed out the difference between the grey-ochre chert artefact (AR21) and local chert, suggesting a different and possibly non-local provenance.Concerning quartzite samples, little information is available in the literature ([40] and references therein), but a local source is the most plausible assumption since dark quartzite occurs in the Vidrieros Fm. [41] and in the Murcia Fm [40].Rock crystal, as well, is possibly local [46], but a precise identification of the source area cannot be established.

Conclusions
The multidisciplinary research gave interesting results for the characterisation of siliceous rock artefacts and the raw materials provenance study, confirming the support of the non-destructive SEM analyses.
The petrographic and mineralogical characterisation carried out on a subset of samples previously analysed by colour analysis and other spectroscopic techniques allowed for the definition of the provenance of the different siliceous rocks, especially for cherts.According to the analytical results, most of the archaeological chert samples have chemical characteristics compatible with natural ones outcropping in the same area, confirming the close supply, which could have possibly favoured the occupation of the site.As regards the origin of the quartzite from the territory, the previous studies do not provide enough information, and samples from geological outcrops should be added to obtain more robust data.Rock crystal samples probably have a local raw material supply, because quartz crystals are commonly found in Mesolithic contexts in the southern areas of the Picos de Europa Fm. and in eastern Asturias.
The study of the physical and mechanical properties has been proved to be a new and very fruitful approach in the characterization of cherts because it can provide useful information on the different mechanical behaviours of siliceous samples, which certainly conditioned the workability and, thus, the uses, technical functions and production of tools of antiquity.The lower crystallinity of cherts affects its flakeability and the typical conchoidal fracturing, resulting in sharp tools suitable for arrowheads and small cutting tools.Quartzite, which is coarser, less porous and denser, was used to produce larger tools such as beating masses and scrapers.

Figure 1 .
Figure 1.Geological map of the Peña de Oviedo sector (west territory of Camaleño village) at scale 1:50.000from Sheet 81-Potes modified, Instituto Geològico y Minero de España web site at link: http://info.igme.es/cartografiadigital/geologica/Magna50Hoja.aspx?language=es&id=81 (accessed on 15 March 2023).The inlets show the position of the study area in Spain (red rectangle) and a sketch map of the tectonic units in which the area is subdivided.The numbers of the formations are taken from the original map (sheet 81-Potes); those in bold indicate the formation hosting cherts and/or quartzites.

Figure 1 .
Figure 1.Geological map of the Peña de Oviedo sector (west territory of Camaleño village) at scale 1:50.000from Sheet 81-Potes modified, Instituto Geològico y Minero de España web site at link: http://info.igme.es/cartografiadigital/geologica/Magna50Hoja.aspx?language=es&id=81 (accessed on 15 March 2023).The inlets show the position of the study area in Spain (red rectangle) and a sketch map of the tectonic units in which the area is subdivided.The numbers of the formations are taken from the original map (sheet 81-Potes); those in bold indicate the formation hosting cherts and/or quartzites.

Figure 3 .
Figure 3. BSE images representative of analysed samples: (a) AR12.1 (grey chert) at low magnification, showing an apparent compact texture with cracks, holes (sometimes filled by phyllosilicates) and bright spots of zircon and metal oxides; (b) millimetres-long trail of Fe and Ti oxides within the quartz matrix (AN5.1,Grey-beigechert from Domeño); (c) AR37 (yellow chert) at moderate magnification that shows the rough surface covered by small crystals and plagues of light-grey phyllosilicates; (d) quartz matrix of AR29 (dark quartzite) hosting a REE-phosphate grain and Ti oxide; (e) compact and homogeneous surface of AR10 (rock crystal) samples with dark-grey stains due to carbon (organic?)impurities; (f) phyllosilicate in quartz matrix (AR34).

Figure 3 .
Figure 3. BSE images representative of analysed samples: (a) AR12.1 (grey chert) at low magnification, showing an apparent compact texture with cracks, holes (sometimes filled by phyllosilicates) and bright spots of zircon and metal oxides; (b) millimetres-long trail of Fe and Ti oxides within the quartz matrix (AN5.1,Grey-beigechert from Domeño); (c) AR37 (yellow chert) at moderate magnification that shows the rough surface covered by small crystals and plagues of light-grey phyllosilicates; (d) quartz matrix of AR29 (dark quartzite) hosting a REE-phosphate grain and Ti oxide; (e) compact and homogeneous surface of AR10 (rock crystal) samples with dark-grey stains due to carbon (organic?)impurities; (f) phyllosilicate in quartz matrix (AR34).

Minerals 2023 ,
13, x FOR PEER REVIEW 9 o

inerals 2023 ,Figure 6 .
Figure 6.Positive correlation between the matrix crystallinity determined by compa ground signals of each sample with the background of a reference known crystallin crystal rock and the grain size means determined by SEM analysis (see Table3).

Figure 6 .
Figure 6.Positive correlation between the matrix crystallinity determined by comparing of the background signals of each sample with the background of a reference known crystallinity of the AR10 crystal rock and the grain size means determined by SEM analysis (see Table3).

Figure 7 .
Figure 7. Petrophysical and mechanical properties of archaeological and geological samples: (a) real density vs. bulk density; (b) real density vs. matrix crystallinity; (c) open helium porosity vs. bulk

Figure 7 .
Figure 7. Petrophysical and mechanical properties of archaeological and geological samples: (a) real density vs. bulk density; (b) real density vs. matrix crystallinity; (c) open helium porosity vs. bulk density; (d) helium open porosity vs. phase size mean; (e) helium open porosity vs. water open porosity; (f) open helium porosity vs. PLT punching strength index.

Figure 8 .
Figure 8. Petrophysical and mechanical properties of archaeological and geological samples phase size mean vs. PLT punching strength index; (b) matrix crystallinity vs. PLT punching stren index.

Figure 8 .
Figure 8. Petrophysical and mechanical properties of archaeological and geological samples: (a) phase size mean vs. PLT punching strength index; (b) matrix crystallinity vs. PLT punching strength index.
Minerals 2023, 13, x FOR PEER REVIEW 17 of 21 latter as raw material.On the other hand, the other chert types show lower concentrations for Al, K, Ti, Fe and Zr.Domeño chert has the highest amounts of Ca.Average concentrations of quartzite are higher than those of cherts for most elements.Quartzite samples also have a higher variance.

Table 1 .
Analysed samples (artefacts and chunk raw materials) from La Calvera rock shelter and geological samples from the Domeño area in Andilla (Valencia, Spain).

Table 2 .
Macroscopic characteristics of archaeological finds and geological samples by optical microscopy in reflect light.
AR21Grey-ochre chert Grey-beige, compact with conchoid-type fracturing and with porous appearanceAR29 Dark quartzitePresence of dark-to-light colour gradient in thickness: dark side more porous and altered, light side shiny as formed by microcrystals (similar to sample AN5.1) AR34 Dark quartzite Dark in colour, porous and altered (similar to the dark side of sample ARCH29 but shinier) AR37 Yellowish chert Yellow-brown colour, with varied surface appearance: from smooth and firm in some places, to extremely porous in others RM03.4 Grey chert Black colour, shiny, almost obsidian-like appearance, with conchoid fractures, rich in alterations and also in diffuse patina (similar to sample 12.1) AN1.2 Blackish chert Blackish colour, compact with conchoid fracturing, low porosity, clean surface (no soil residue on the surface) AN5.1 Grey-beige chert Grey-beige surface with no soil residue on the surface, porous but compact appearance with conchoid-type fracturing

Table 3 .
Microscopic and mineralogical characteristics of archaeological finds and geological samples using SEM-EDS analyses.

Table 4 .
Mineralogy of archaeological finds and geological samples by XRD analysis.

Table 5 .
Petrophysical and mechanical properties of archaeological and geological samples.
Legend: ρ R = real density; ρ B = bulk density; Φ O He = helium open porosity; Φ O H 2 O = water open porosity; IC W = imbibition coefficient; SI = saturation index; Is (50) = point load test strength index; R C = indirect calculated compression strength; R T = indirect calculated tensile strength.

Table 6 .
Average elemental concentrations of the different classes of samples obtained by pXRF

Table 6 .
Average elemental concentrations of the different classes of samples obtained by pXRF.Note: Elemental concentrations are expressed as mass percentage, except for Zr, which are expressed as mg/kg.<LD: below the limits of detections; n = number of samples; SD: standard deviation.