Abstract

Forty-two mortar samples, from two archaeological excavations located in Sagunto (Valencian Community, Spain), were analysed by both portable energy dispersive X-ray fluorescence spectroscopy (pED-XRF) and inductively coupled plasma mass spectrometry (ICP-MS) to determine major and minor elements and traces including rare earth elements (REEs). Collected data were crossed with those previously obtained from Sagunto Castle mortars, and principal component analysis (PCA) was applied to discriminate the construction phases of the unearthed buildings. REE permitted to ascribe most of the masonries to the Roman Imperial period. Moreover, a statistical model was built by employing partial least squares discriminant analysis (PLS-DA) in order to classify the mortars from Roman Imperial period and from Islamic period due to the problematic overlapping between these two phases. Results confirmed the effectiveness of the developed indirect chronology method, based on REE data, to discriminate among historic mortars from different construction periods on a wide scale including different Sagunto archaeological sites.

1. Introduction

Mortar is a building material composed essentially of binder and aggregate fractions and, in some cases, of additives of different types [13]. In particular, lime mortar had a key function in Roman architecture: Roman people made their walls and structures with this material adding reactive materials, like pozzolan materials (i.e., pulvis puteolanus and cocciopesto), to give a hydraulic character to the mortars [48]. The use of mortar in architecture is documented during the Middle Ages and in the following historical periods [613].

The chemical analysis by statistical approach, together with mineralogical and petrographic characterization of ancient mortars and polished stones, has shown to be a useful tool in the interpretation of the construction phases of several archaeological sites and historical complexes [1420]. Moreover, ancient mortars are subjected to decay phenomena which also must be detected and evaluated for conservation issues [2124].

This paper shows the results of the analyses carried out on ancient mortars collected from buildings discovered during two recent archaeological excavations at Sagunto, a town located in the Eastern Spain, ca. 30 km north of Valencia, close to the Costa del Azahar on the Mediterranean Sea. Sagunto is well known in the world for its complex history, and its area has been occupied since the Iberian Age. During the Roman period, Sagunto was interested in constructing important buildings such as the Circus and the Theatre. Thereafter, Sagunto was occupied by Islamic people, and during the Modern Ages, the city was involved in the Napoleonic Wars [25].

The collected mortar samples were analysed by both portable energy dispersive X-ray fluorescence spectroscopy (pED-XRF) and by an inductively coupled plasma mass spectrometer (ICP-MS) to determine their major elements, minor elements, and trace elements. The mineral element concentration and, in particular, the measured contents of rare earth elements (REE) of these samples were compared with those previously obtained from Sagunto Castle [26] by applying principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) to determine the construction periods.

2. Materials and Methods

2.1. Sampling

The studied mortars (Table 1) were sampled during two archaeological excavations located in the city of Sagunto (Figure 1). Twenty-nine samples were collected during the archaeological excavation of Los Huertos St., an area characterized by the presence of buildings and materials from the Roman Imperial phase and from the Islamic occupation in the Middle Ages. Twenty-three samples of lime mortars were collected from these buildings: nine samples come from the Circus (C), a building dating back to the Roman Imperial period, one from the cloaca (CLC), one from a jamb (JMN), two from two different walls (Eastern Wall, EW, and Northern Wall, NW), one from a noria (NR), and four from a pilaster (P). Six samples of earth mortar come from two rooms that the archaeologists interpreted as may belong to the Islamic occupation (Room 1, samples R1 and Room 2, samples R2), and in this study, they were analysed as control samples. Two samples of lime mortar were collected from Room 3 (R3), two from an unidentified building (UN), and one from a well (WLL). Thirteen samples of lime mortars were collected from different masonries during the emergency archaeological excavation of an area close to Sagunto’s Railroad Station (RS named samples). Data of twenty-three lime mortar samples from different buildings of the Castle of Sagunto, studied by Gallello et al. [26], were taken into account as the calibration set to perform data analysis and to date the structures of the above-quoted excavations. Part of the calibration samples come from masonries dating back to the Roman Imperial Period: five samples from the Theatre (TR) and one from the Curia (CUR); the other samples come from masonries dated back to the Islamic occupation phase: four samples from the first part of the Islamic Wall (MI), one of the wall which was considered modern before the study (MM), one from a Hermitage (ERM), two from two Islamic reworks of the Torre Central Estudiantes (TCE), and two from the tabernae of the Imperial Forum. Moreover, two samples from the Curia (CUR) and five samples from the Basilica of the Imperial Forum (FBI), whose dating was uncertain between the above-quoted two phases, were added to the dataset.

2.2. Major and Minor Elements Determination

All the analyses were carried out on each entire mortar sample previously pulverized (Dmax < 63 μm) and homogenized through Agatha mortar and pestle. Major and minor element concentrations were obtained by using a S1 Titan energy dispersive portable X-ray fluorescence spectrometer (pED-XRF) from Bruker (Kennewick, Washington DC, USA) equipped with an Rh X-ray tube and X-Flash® SDD. Geochem-trace calibration was used to perform the quantitative analyses, and S1 Sync software from Bruker was employed to measure Al2O3, SiO2, CaO, Ti, and Fe. The accuracy of the measurements was verified by using the following certified reference materials: soil NIM GBW07408 and limestone NCS DC60108a (Table 2). All the reading’s standard errors range between 1 and 5 wt.%, except for Al2O3 measured on NCS DC60108a, whose standard deviation increases up to 22 wt.% for concentration less than 0.5 wt.% probably due to the low sensitivity of the instrument for the determination of this element.

2.3. Trace Elements and REE Determination

Previously pulverized and homogenized samples of each entire mortar were prepared for the inductively coupled plasma mass spectrometry (ICP-MS) analysis. The mortars digestion method was developed by Gallello et al. [26] in order to provide reproducible and comparable results compatible with the sensitivity of the analytical method employed. The wet digestion consisted in the addition of 1.35 mL of HCl (37%) and 0.45 mL of HNO3 (69%) to ca. 0.15 g of each sample in glass tubes, which were heated in a boiling water bath for about 40 min. Then, the solutions were poured into plastic tubes and brought up to 25 mL with purified water. The concentration of the following elements: Ba, Bi, Cd, Cr, Co, Cu, Pb, Li, Mn, Mo, Ni, Sr, Tl, V, and Zn, REE (La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu), Sc, and Y were determined. Two multielement stock solutions for the ICP analysis in 5% HNO3, containing the abovequoted elements at a concentration of 1000 mg/L, were used as stock standards for calibration. 5 mL volumetric flasks were used in adding the corresponding volume of standard solutions 0.15 mL of HNO3, 0.45 mL of HCl, and the water necessary to reach the final volume. The concentration of trace elements ranges from 1 μg/L to 600 μg/L except for REE, Y, and Sc that ranges from 1 μg/L to 100 μg/L. The measurement accuracy was verified by using the certified reference materials soil NIM GBW07408 and limestone NCS DC60108a. As the internal standard, 50 μL of a 1000 mg/L Rh solution was added to each sample and to each calibration standard. The analyses were performed through an Agilent 7900 inductively coupled plasma mass spectrometer. The measurement conditions that are shown in Table 3 shows the ICP-MS parameters employed for the analyses, and Table 4 shows the main analytical features obtained for the measured mass of each considered isotope, including the instrumental detection and quantification limits (LOD and LOQ), and the coefficient of determination (R2) of the corresponding calibration lines.

2.4. Statistical Data Processing

The PCA models were built by using data obtained from a total of forty-two analysed samples, and the set of data obtained from samples dated to the Roman Imperial period and the Islamic occupation phase obtained previously by Gallello et al. [26]. Major and minor elements and trace elements including REE were employed as variables. This technique was used to explore large geochemical datasets by reducing the number of variables and providing a deep insight into the structure of the variance of the dataset. Data were processed through mean center and autoscale prior to modelling, and the obtained model was cross validated through leave-one-out method.

The PLS-DA model was built employing twenty-two samples as a calibration set: nine samples from the Circus (Los Huertos Street excavation) and samples studied by Gallello et al. [26] including five from the Roman Theatre, TR, built during the Roman Imperial Period and eight samples from Islamic masonries (two samples from the Imperial Forum tabernae, TFI, four samples from the first part of the Islamic Wall, MI, and two from the Torre Central Estudiantes, TCE). Five of these samples (C07, C08, C10b, and two TR samples) were randomly selected as an internal validation set to test the model (data not shown). Sixteen variables were employed (REE, Y, and Sc). Finally, thirty-seven mortar samples dated back to Roman Imperial period, Islamic occupation phase, and also some uncertain samples were included as the test set to predict their class. Data were preprocessed employing mean centering and autoscale, and “leave one out” cross validation method was employed for the evaluation of results. Data analysis was carried out using the PLS Toolbox 6.5 for Eigenvector Research Inc., (Wenatchee, WA, USA) running in Matlab R2014b from MathWorks Inc., (Natick, MA, USA).

3. Results and Discussion

3.1. Geochemical Results

The analytical results of major and minor elements as well as trace elements and REE data are reported in the Supplementary Materials (Annexes 1 and 2).

Lime mortar samples from the Railroad Station and from the Circus area buildings have comparable amounts of Al2O3 and SiO2 which range from 0.44 to 1.05 wt.% and from 11.30 to 26.97 wt.%, respectively, except for samples RS25A (Al2O3: 3.02 wt.%; SiO2: 34.86 wt.%) and RS11C (Al2O3: 5.17 wt.%; SiO2: 36.72 wt.%), which are both from the Railroad Station and show particularly high values of these two chemical components. On the contrary, Railroad Station’s samples have lower concentration of CaO (30 ± 4 wt.%) and higher one of Fe (1.3 ± 0.1 wt.%) and Ti (0.12 ± 0.01 wt.%) than Los Huertos Street mortars (CaO: 40 ± 3 wt.%; Fe: 0.99 ± 0.14 wt.%; Ti: 0.08 ± 0.02 wt.%). The samples of earth mortars from Room 1 (R1a, R1b, and R1c) and Room 2 (R2a, R2b, and R2d) show the lowest values of CaO (15.3 ± 1.1 wt.%) and the highest ones of Al2O3 (10.9 ± 1.1 wt.%), SiO2 (37.8 ± 1.7 wt.%), Fe (2.4 ± 0.3 wt.%), and Ti (0.25 ± 0.02 wt.%).

The lime mortar samples from Los Huertos Street and from the Railroad Station have comparable concentrations in almost all the measured trace elements. In particular, it is worth noticing that the two groups of samples show similar REE total amounts (∑REE) (Los Huertos Street: 20 ± 4 μg/g; Railroad Station: 20 ± 2 μg/g) and Sc concentrations (Los Huertos Street: 1.1 ± 0.3 μg/g; Railroad Station: 1.1 ± 0.1 μg/g), while samples from Los Huertos Street have a slightly higher concentration of Y than the ones from the Railroad Station excavation (4.1 ± 0.8 μg/g and 3.5 ± 0.4 μg/g resp.). Concerning the earth mortars from Room 1 and Room 2, they have comparable contents of REE, Y, and Sc. REE total amounts go from 37 μg/g to 53 μg/g, and Y and Sc range from 5 μg/g to 10 μg/g and from 2 μg/g to 3 μg/g, respectively.

To give a deep insight into the mechanisms that influence the amount of REE in the mortar samples, it is interesting to look at the correlation among the REE and analysed elements for each group of mortars. Table 5 shows the Pearson correlation coefficient () considering as variables REE total amounts and the indicated elements for each group of mortars considered separately (see Supplementary Materials Annex 2). Only statistically significant data were taken into account given . REE concentration of Los Huertos Street lime mortar samples show positive correlations, as well as for Sc (), which has a chemical behaviour similar to lanthanides, with Al2O3 (), Tl (), Co (), and Cr () and negative correlations with Ca () and Sr (), which suggests that the major contribution to the amount of REE comes from limestone clay impurities or from aluminosilicates and silicate rock clasts present in the aggregate fraction [27]. Concerning the lime mortar samples from the Railroad Station excavation, significant correlations were not detected among major elements and REE. However, significant positive correlations were found with Y (), Co () and V (). As regards earth mortars (R1 and R2) from Los Huertos Street, in this case, the REE total amount has high positive correlations with most of the trace elements (for Bi, Tl, Ba, Mo, Sr, Zn, Cu, Ni, Co, Mn, Cr, V, Sc, and Li, goes from 0.71 to 0.97) and, among the major elements, with Fe () and Ti (), suggesting that the major contribution in lanthanides comes probably from clay fraction, other silicate minerals, and silicate rock clasts.

3.2. Chemometrics for Construction Phase Discrimination

The identification of the construction phases of the buildings was conducted by comparing the mortars from the Railroad Station and Los Huertos Street to the ones from Sagunto Castle published by Gallello et al. [26], and whose construction periods were retraced in the same work. In particular, this previous study was focused on REE concentrations due to their proved effectiveness in archaeometric studies about provenance and raw materials of lithic and lithoid archaeological artifacts [2830].

To compare the samples from the excavations to the Imperial and Islamic mortars from Sagunto Castle [26], the principal component analysis (PCA) was carried using all variables (i.e., elements) (Figure 2). Figure 2(b) shows the contribution of the variables in PC1 (55.00% of samples’ variance) and PC2 (13.37% of samples’ variance). Figure 2(a) shows the position of each sample in a diagram where the x-axis corresponds to PC1 coordinates and y-axis to PC2. Most of the samples of Los Huertos Street and all Railroad Station are grouped together in the left side of the plot, while the samples dated back both to Imperial Roman period and Islamic period from de Castle are grouped in the right part of the graph. The samples of earth mortars from Room 1 (R1a, R1b, and R1c) and Room 2 (R2a, R2b, and R2d) of Los Huertos Street are plotted in the upper right side of the diagram, at the limit of the confidence interval or out of it.

In Figure 3(b), the contribution of the REE employed as variables is shown. PC1 explains 88.34% of the data variance and PC2 explains 5.40%. Lime mortar samples from Los Huertos Street and all Railroad Station excavations are grouped with the mortars from the Roman Theatre (TR) (Figure 3(a)), suggesting that the dating to the Roman Imperial period was also confirmed by the archaeological data while the earth mortar samples are plotted in PC1 and PC2 positive directions. Finally, a separation between mortars from Imperial Period (PC1 negative direction), and mortars from the Islamic Period, located in the left lower area of PC can be appreciated.

Since the discrimination of lime mortars from Roman Imperial and Islamic buildings seems to be the most problematic issue, the principal least squares discriminant analysis (PLS-DA) was employed to build a statistical model to classify uncertain samples. The calibration set and validation set were established by using samples from buildings and masonries that were previously classified by Gallello et al. [26] and by archaeological data as Roman Imperial (C and TR) or as Islamic (MII, TFI, and TCE) (Figure 4(a)). Figure 4(b) shows the contribution of each variable to the model. As it can be seen, the model clearly separates between mortars of the two different construction phases. Subsequently, the model was applied to a test set composed of uncertain Roman Imperial period or Islamic mortars (Figure 4(a)). The PLS-DA results confirmed that all the mortars from Los Huertos Street and Railroad Station are classified as Roman Imperial period materials except for sample WLL that has been classified as Islamic, which is coherent with the archaeological data that confirm the occupation of the area during the Islamic phase. As observed, on the contrary of Los Huertos Street, in Railroad Station samples, major elements are not correlated with REE (Table 5), and this could suggest that the mortars were made in a different way or this is an indicator of recycling building materials, and, according to the archaeological wall stratigraphy interpretation, the Railroad Station structures could be dated after the Imperial period. Thus, maybe they were built during the Late Roman Age. On the contrary, one mortar sample from the Basilica of the Imperial Forum (FBI SC) and one mortar from the Curia testify the presence of both Roman Imperial masonries and Islamic reworks in Roman buildings due to the heavy reemployment of the Imperial Forum area during the Islamic occupation [26].

The use of REE, Y, and Sc as variables has shown their high effectiveness in the discrimination of different kinds of mortars to retrace the construction history of archaeological sites and monuments. The comparison among mortars coming from different archaeological excavations and monuments of the same geographical area seems to be a very viable method to solve chronological issues related to the construction phases of both structures unearthed during archaeological excavations and historical complexes.

4. Conclusion

The study of the historic mortars from two archaeological excavations (i.e., Los Huertos Street and Railroad Station) located in the city of Sagunto permitted us to solve issues that had remained unanswered in the previous studies on these building materials.

The use of multivariate statistics employing REE as variables allowed us to classify the lime mortars belonging to the Imperial Roman period, in line with the archaeological data; however, the case of the Railroad Station structure needs deeper studies to confirm the chronology as Imperial Roman or Late Roman Ages. Also, it was possible to establish the chronological phase of some masonries of Sagunto Castle whose period of construction was uncertain and that can be attributed to the phase of the Islamic occupation or to the Roman Imperial period. The PLS-DA model was helpful to discriminate among mortars of uncertain attribution and that can be useful for possible future studies in Sagunto area.

The correlation among REE and major/minor and trace elements shows that the mechanisms that allowed REE to be discriminating parameters between Sagunto mortars belonging to different periods are related to chemical weathering processes involving limestone clay impurities, aluminosilicates, and silicate rock clasts present in the aggregate fraction. Therefore, further works need to be developed focusing the understanding of the chemical processes that are influencing the REE distributions.

Summarizing, the present study has confirmed the effectiveness of REE data in building materials as markers of different construction phases and has proved the usefulness of the application of the ancient mortar analysis on a geographical scale wider than a single archaeological excavation or monuments.

Data Availability

The data used to support the findings of this study are included within the article and supplementary materials.

Conflicts of Interest

The authors declare that there are no conflicts of interest with any institution or funding body.

Acknowledgments

The authors acknowledge also the analytical support of the Central Service for Experimental Research (SCSIE) of the University of Valencia, and of Dr. Laredo Ortíz and Dr. Soriano Vega. The authors would like to thank all the archaeologists working in Railroad Station and Los Huertos Street excavations, who have contributed the realization of this study and especially the Archaeologist at the Sagunto’s Council, Ms. Carla Flors. Authors acknowledge the financial support of Generalitat Valenciana (PROMETEO project II/2014/077) and Ministerio de Economia y Competitividad-Feder (Project CTQ 2014-52841-P and Project CTQ 2012-38635). Gianni Gallello acknowledges the financial support of the European Commission (Project H2020-MSCA-IF-2015-704709-MATRIX).

Supplementary Materials

Annex 1: concentrations of the measured major and minor elements obtained by pED-XRF. Annex 2: trace elements and REE results obtained by ICP-MS. (Supplementary Materials)