# Lanthanides Revealing Anthropogenic Impact within a Stratigraphic Sequence

Authors: Gianni Gallello; Agustín Pastor; Agustín Diez; Joan Bernabeu
Journal: Journal of Archaeology
 (2014)

---

## Abstract

Difficulties to differentiate between anthropogenic and natural processes in the formation of archaeological deposits are crucial for a correct interpretation not only of the actions involved in the development of archaeological sites, but also of their occupation-abandonment dynamics and the understanding of their spatial behaviors and relationship with the environment. We have carried out lanthanides (rare earth elements “REE”) analysis to distinguish anthropogenic from natural stratigraphic units in sediments using the advantage of the high sensibility, precision, and accuracy of ICP-MS measurements. In the Neolithic site of Mas d’Is (Alacant, Spain), we have applied REE analysis in a huge stratigraphic sequence called Pit 6, which was known to contain a large anthropogenic component. Randomly collected soil samples were sequentially taken in order to identify anthropogenic soil formations and to prove the proposed method blind testing has been used. In the specific case of Mas d’Is excavation a recurring question is whether paleosols are at the origin of the human occupation of the sites or it was the occupation of this areas which triggered the paleosols development. Our purpose was to distinguish the degree of human contribution to paleosols formation between samples sequentially taken at few centimeters of distances in a giant stratigraphic sequence (Pit 6) employing REE analysis.


---

## Body

## 1. Introduction

Neolithic site of Mas d’Is (Alacant, Spain) and its hinterland (Figure1) has led us to face some problems about the interaction of anthropogenic and natural processes [1]. At Mas d’Is we have discovered huge Neolithic stratigraphic sequences >5 meters that are not at all common in western Mediterranean for an open-air site [2] or elsewhere for the filling of some previously excavated features. New methodological proposals that can contribute to a better understanding of some sedimentological problems in this archaeological site have been tested. We have successfully applied lanthanides (rare earth elements “REE”) analysis to distinguish anthropogenic from natural stratigraphic units in predefined classes [1]. The potential of REE to distinguish anthropogenic versus natural sediments was useful in Mas d’Is case, where other well-known methodologies (sedimentology, microsedimentology, and petrography) fail in discriminating between human and natural soils [1]. In this study we have applied REE analysis in one of those huge stratigraphic sequences called Pit 6. Randomly collected soil samples were sequentially taken in order to identify anthropogenic soil formations and to prove the proposed method blind testing has been used.Figure 1
Location of the upper and middle Serpis Valley.Recently, many studies have concentrated their attention in different environments to know how human activity is recorded in soil formation. Major and trace elements chemical analysis has a long standing application in archaeology [3–6]. Multielemental soil analysis has been employed as a tool of geochemical prospecting [7–10] and to support interpretation of past workplaces and activity spaces around and within archaeological structures [11, 12]. In the last three decades REE soil analysis has started to be part of archaeometric studies. It is the strong partition of rear earth elements into the mineral particulate phase, their coherent behavior during weathering, erosion and fluvial transportation, and their high resistance to chemical mobilization [13–16] which underlie the use of REE in sediments as provenance indicators in geochemical studies for some time now [17–22]. Due to their peculiar characteristics, there have been several attempts to apply REE soil analysis to different archaeological scenarios [23–26].At the Mas d’Is, to test our methodological approach we have analyzed 15 soil samples. REE were measured using inductively coupled plasma mass spectrometry (ICP-MS). The multielement capabilities, high sensibility, precision, and accuracy of ICP-MS have been exploited to evaluate variations on REE between our samples. In the specific case of Mas d’Is excavation a recurring question is whether paleosols are at the origin of the human occupation of the sites or it was the occupation of this areas which triggered the paleosols development. Our purpose was to distinguish the degree of human contribution to paleosols formation between samples sequentially taken at few centimeters of distances in a giant stratigraphic sequence (Pit 6) employing REE analysis.


## 2. Material and Method

### 2.1. Geographical and Archaeological Background

The samples were collected on Mas d’Is archaeological site located in the mountainous inland area of the northernmost Alacant province (Figure1), within Penàguila headwaters, one of the main Serpis tributaries. The area where the site is located is an old platform where heavy erosion processes have led to a series of deep ravines. Holocene landscape development appears to have progressed from a period of stability to slope denudation with aggradation (stream infilling) followed by rapid incision which initiated sometime near the time of occupation-in fact. Some structures of Mas d’Is are eroded by ravines, which confirms that some incision process occurs after early Neolithic occupation (ca. 5600–4800 cal. BC). The result is a series of deep canyons cut by interfluves receiving the distinctive name of Les Puntes, one of which houses the Neolithic site of Mas d’Is. The current profile is the result of ongoing agricultural changes that ultimately have levelled the previous agricultural terracing topography masking and destroying much of the close-to-surface archaeological record.The sedimentological environment at Penàguila valley can be characterized as follows: on top of the bedrock, very homogenous beds of Miocene marls and almost 1000 meters thick a set of sediments ranging from weathered marls to alluvial deposits. Those alluvial deposits are framed by the development of important paleosols predating the Neolithic occupation of the site. Around ca. 5600 cal. BC, the first farmers appear at Mas d’Is, on the top of one of these paleosols.Archaeological investigations, beginning at the end of 19th century, have identified an important number of human sites dated to the Early and Middle Holocene. More recently, extensive surveys and excavations of Neolithic villages have provided more detailed information about human settlements in the Alcoi valley. Currently, we have documented more than 100 Neolithic open-air sites that almost without exception are related to paleosols.The excavated area at Mas d’Is provides different kinds of structures and archaeological materials that chronologically span between the very first Neolithic in the region (ca. 5600 cal. B.C.) and present times. The space occupied by the Neolithic site has been an agricultural plot since the Early Neolithic. In the sixth millennium cal BC the site could be divided in two main areas: the village, where the domestic features (houses, hearths, and mills) are found [27] and the area of the ditches, over 150 m east. Some ditches have an important sequence of infilling sediments, covering more than one millennium (from ca. 5600–4000 cal. BC). On the top, a disturbed layer of Bronze Age ends the archaeological occupation of this area [1]. Pit 6 is one of these features.


### 2.2. The Samples

Fifteen samples (M1–M15) were collected (Figures2-3) and were analyzed by ICP-MS (Table 1). Some samples (M1–M12) were taken from the button to the top in Pit 6 (Table 1 and Figures 2 and 3). Sample M1 was collected, 12 cms above the bad rock. This is apparently clayey silt with some organic content to judge by its color (dark gray). Sample M2 was taken 12 cms above M1. It is composed by dark-colored clayey silt with some gravel intrusions. Sample M3 was taken 10 cms above M2 and it is formed by mudstones intrusion with very dark color. Sample M4 was taken 25 cms above M3 and it is composed of a bundle of gray silt with gravel intrusions. Sample M5 was collected 15 cms above M4 and it is composed by a horizontal band above clayey silt with high organic content judging by its dark color. Sample M6 is located 10 cms above M5 and it is composed of a package of clear silt, with abundant intrusions of loam blocks, and appears to have been originated from a loosening of Pit 6 walls. The sample M7 was collected 45 cms above M6 and it is composed of a package of dark silt with somewhat of organic matter. Sample M8 is situated 70 cms above M7, a package of dark silts that may correspond to one of the paleosols documented filling this gap. Sample M9 was taken 40 cms above M8 and is composed of a package of silt with little intrusion.Table 1
A total of 15 samples (M1–M15) were analysed. Stratigraphic Unit (SU). Centimeters of distance between the samples (cms).
SAMPLE
ZONE
SU
CRONOLOGY
ANALYSIS
Cms
M1
PIT 6
—
—
REE
12
M2
PIT 6
—
—
REE
22
M3
PIT 6
—
—
REE
33
M4
PIT 6
—
—
REE
58
M5
PIT 6
—
—
REE
73
M6
PIT 6
—
—
REE
83
M7
PIT 6
—
—
REE
128
M8
PIT 6
—
—
REE
198
M9
PIT 6
—
—
REE
238
M10
PIT 6
104355
—
REE
263
M11
PIT 6
—
—
REE
338
M12
PIT 6
—
>5000 cal BC
REE
358
M13
PIT 6
Tertiary
REE
—
M14
C
Bedrock
Tertiary
REE
—
M15
BC105
—
Bronze
REE
—Figure 2
Profile of Foso 6 infilling layers at the Mas d’Is. Location of samples (M1–M7).Figure 3
Profile of Foso 6 infilling layers at the Mas d’Is. Location of samples (M8–M12). Sample M10 (UE104355).Sample M10 (UE104355) was located 25 cms above M9. This sample is composed of a darkish silt package that seems to correspond to a deposition period. In this stratigraphic unit (UE104355) a large core of flint was found indicating that this sample was collected in a unit where anthropogenic components are present. Sample M11 was collected 75 cms above sample M10 and is formed by a clear silt package that interposes a well-documented paleosol. Sample M12, 20 cms above the previous M11, is composed of dark silt package that is associated with a paleosol that fills an earlier erosion period dated after the ca. 5000 cal. BC. Sample M13 is marl in contact with the UE104355 (M10). The M14 (marl) has been directly taken from the parent rock. The M15 sample was sampled in the sector BC105 and it is composed by a dark area and silty sediment from the Bronze Age dated by the recovered ceramic fragments.


### 2.3. Chemical Analysis of Soil

The sampling was carried out collecting the sediments sequentially from Pit 6 stratigraphic wall (Figures2-3) and not sieved. The 15 samples have been sampled using microspoon spatula made of stainless steel, always cleaned before taking a new sample, and stored in 15 mL test tubes. Sample preparation and digestion have been carried out as described by Gallello et al. [1]. The process includes a thermal treatment at 450°C. Successively the samples were crushed and homogenized employing an agate mortar. A partial soil analysis has been carried out and we have developed a digestion method and a range of dilutions from the digested solution in order to provide reproducible and comparable results compatible with the sensitivity of the analytical method employed.The digestion method consisted of the addition of 1.5 mL HCl and 1.5 mL HNO3 to 0.5 g of sample in glass tubes placing them in a water bath at 100°C for 40 min. In addition, two blanks were prepared using the same acid volumes and digestion conditions of the samples. Subsequently, the digested solutions have been carefully poured into plastic tubes of 15 mL, bringing the volume to 15 mL with purified water. This concentrated solution (A) was diluted 1 : 10 obtaining a solution (B) to measure La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Sc, and Y. Dilution containing 10% concentrated HCl and 10% concentrated HNO3 has been maintained in all solutions. A multielemental stock solution containing REE, Sc, and Y at concentration of 100 μg/mL was prepared. The calibration standards were prepared employing 50 mL volumetric flasks and adding 5 mL HCl + 5 mL HNO3 plus the corresponding volume of standard solution and filling up to volume with pure water. Solutions were analyzed by ICP-MS with Perkin Elmer Elan DRCII (Concord, Ontario, Canada). To avoid the obstruction of the nebulizer system the sample was filtered employing filter paper (Whatman N.1 of 70 mm). Concentrations ranging between 0 and 0.6 μg/mL have been used for La, Ce, Pr, and Nd and concentrations ranging between 0 and 0.1 μg/mL for Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu. All standards were acquired from Sharlab S.L. (Barcelona). The standard error of readings during the analysis ranged from 3% to 9% for REE. Soil GBW07408 has been used as standard reference material for evaluating the analytical method. Rh was used as internal standard. REE concentrations were normalized to Post-Archaean Australian Shales (PAAS) concentrations using the values of Taylor and McLennan [28].


## 2.1. Geographical and Archaeological Background

The samples were collected on Mas d’Is archaeological site located in the mountainous inland area of the northernmost Alacant province (Figure1), within Penàguila headwaters, one of the main Serpis tributaries. The area where the site is located is an old platform where heavy erosion processes have led to a series of deep ravines. Holocene landscape development appears to have progressed from a period of stability to slope denudation with aggradation (stream infilling) followed by rapid incision which initiated sometime near the time of occupation-in fact. Some structures of Mas d’Is are eroded by ravines, which confirms that some incision process occurs after early Neolithic occupation (ca. 5600–4800 cal. BC). The result is a series of deep canyons cut by interfluves receiving the distinctive name of Les Puntes, one of which houses the Neolithic site of Mas d’Is. The current profile is the result of ongoing agricultural changes that ultimately have levelled the previous agricultural terracing topography masking and destroying much of the close-to-surface archaeological record.The sedimentological environment at Penàguila valley can be characterized as follows: on top of the bedrock, very homogenous beds of Miocene marls and almost 1000 meters thick a set of sediments ranging from weathered marls to alluvial deposits. Those alluvial deposits are framed by the development of important paleosols predating the Neolithic occupation of the site. Around ca. 5600 cal. BC, the first farmers appear at Mas d’Is, on the top of one of these paleosols.Archaeological investigations, beginning at the end of 19th century, have identified an important number of human sites dated to the Early and Middle Holocene. More recently, extensive surveys and excavations of Neolithic villages have provided more detailed information about human settlements in the Alcoi valley. Currently, we have documented more than 100 Neolithic open-air sites that almost without exception are related to paleosols.The excavated area at Mas d’Is provides different kinds of structures and archaeological materials that chronologically span between the very first Neolithic in the region (ca. 5600 cal. B.C.) and present times. The space occupied by the Neolithic site has been an agricultural plot since the Early Neolithic. In the sixth millennium cal BC the site could be divided in two main areas: the village, where the domestic features (houses, hearths, and mills) are found [27] and the area of the ditches, over 150 m east. Some ditches have an important sequence of infilling sediments, covering more than one millennium (from ca. 5600–4000 cal. BC). On the top, a disturbed layer of Bronze Age ends the archaeological occupation of this area [1]. Pit 6 is one of these features.


## 2.2. The Samples

Fifteen samples (M1–M15) were collected (Figures2-3) and were analyzed by ICP-MS (Table 1). Some samples (M1–M12) were taken from the button to the top in Pit 6 (Table 1 and Figures 2 and 3). Sample M1 was collected, 12 cms above the bad rock. This is apparently clayey silt with some organic content to judge by its color (dark gray). Sample M2 was taken 12 cms above M1. It is composed by dark-colored clayey silt with some gravel intrusions. Sample M3 was taken 10 cms above M2 and it is formed by mudstones intrusion with very dark color. Sample M4 was taken 25 cms above M3 and it is composed of a bundle of gray silt with gravel intrusions. Sample M5 was collected 15 cms above M4 and it is composed by a horizontal band above clayey silt with high organic content judging by its dark color. Sample M6 is located 10 cms above M5 and it is composed of a package of clear silt, with abundant intrusions of loam blocks, and appears to have been originated from a loosening of Pit 6 walls. The sample M7 was collected 45 cms above M6 and it is composed of a package of dark silt with somewhat of organic matter. Sample M8 is situated 70 cms above M7, a package of dark silts that may correspond to one of the paleosols documented filling this gap. Sample M9 was taken 40 cms above M8 and is composed of a package of silt with little intrusion.Table 1
A total of 15 samples (M1–M15) were analysed. Stratigraphic Unit (SU). Centimeters of distance between the samples (cms).
SAMPLE
ZONE
SU
CRONOLOGY
ANALYSIS
Cms
M1
PIT 6
—
—
REE
12
M2
PIT 6
—
—
REE
22
M3
PIT 6
—
—
REE
33
M4
PIT 6
—
—
REE
58
M5
PIT 6
—
—
REE
73
M6
PIT 6
—
—
REE
83
M7
PIT 6
—
—
REE
128
M8
PIT 6
—
—
REE
198
M9
PIT 6
—
—
REE
238
M10
PIT 6
104355
—
REE
263
M11
PIT 6
—
—
REE
338
M12
PIT 6
—
>5000 cal BC
REE
358
M13
PIT 6
Tertiary
REE
—
M14
C
Bedrock
Tertiary
REE
—
M15
BC105
—
Bronze
REE
—Figure 2
Profile of Foso 6 infilling layers at the Mas d’Is. Location of samples (M1–M7).Figure 3
Profile of Foso 6 infilling layers at the Mas d’Is. Location of samples (M8–M12). Sample M10 (UE104355).Sample M10 (UE104355) was located 25 cms above M9. This sample is composed of a darkish silt package that seems to correspond to a deposition period. In this stratigraphic unit (UE104355) a large core of flint was found indicating that this sample was collected in a unit where anthropogenic components are present. Sample M11 was collected 75 cms above sample M10 and is formed by a clear silt package that interposes a well-documented paleosol. Sample M12, 20 cms above the previous M11, is composed of dark silt package that is associated with a paleosol that fills an earlier erosion period dated after the ca. 5000 cal. BC. Sample M13 is marl in contact with the UE104355 (M10). The M14 (marl) has been directly taken from the parent rock. The M15 sample was sampled in the sector BC105 and it is composed by a dark area and silty sediment from the Bronze Age dated by the recovered ceramic fragments.


## 2.3. Chemical Analysis of Soil

The sampling was carried out collecting the sediments sequentially from Pit 6 stratigraphic wall (Figures2-3) and not sieved. The 15 samples have been sampled using microspoon spatula made of stainless steel, always cleaned before taking a new sample, and stored in 15 mL test tubes. Sample preparation and digestion have been carried out as described by Gallello et al. [1]. The process includes a thermal treatment at 450°C. Successively the samples were crushed and homogenized employing an agate mortar. A partial soil analysis has been carried out and we have developed a digestion method and a range of dilutions from the digested solution in order to provide reproducible and comparable results compatible with the sensitivity of the analytical method employed.The digestion method consisted of the addition of 1.5 mL HCl and 1.5 mL HNO3 to 0.5 g of sample in glass tubes placing them in a water bath at 100°C for 40 min. In addition, two blanks were prepared using the same acid volumes and digestion conditions of the samples. Subsequently, the digested solutions have been carefully poured into plastic tubes of 15 mL, bringing the volume to 15 mL with purified water. This concentrated solution (A) was diluted 1 : 10 obtaining a solution (B) to measure La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Sc, and Y. Dilution containing 10% concentrated HCl and 10% concentrated HNO3 has been maintained in all solutions. A multielemental stock solution containing REE, Sc, and Y at concentration of 100 μg/mL was prepared. The calibration standards were prepared employing 50 mL volumetric flasks and adding 5 mL HCl + 5 mL HNO3 plus the corresponding volume of standard solution and filling up to volume with pure water. Solutions were analyzed by ICP-MS with Perkin Elmer Elan DRCII (Concord, Ontario, Canada). To avoid the obstruction of the nebulizer system the sample was filtered employing filter paper (Whatman N.1 of 70 mm). Concentrations ranging between 0 and 0.6 μg/mL have been used for La, Ce, Pr, and Nd and concentrations ranging between 0 and 0.1 μg/mL for Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu. All standards were acquired from Sharlab S.L. (Barcelona). The standard error of readings during the analysis ranged from 3% to 9% for REE. Soil GBW07408 has been used as standard reference material for evaluating the analytical method. Rh was used as internal standard. REE concentrations were normalized to Post-Archaean Australian Shales (PAAS) concentrations using the values of Taylor and McLennan [28].


## 3. Results and Discussion

### 3.1. Composition of Samples

15 samples have been analyzed by ICP-MS to determine REE, Sc, and Y contents. The analyzed isotopes, instrumental LOD and LOQ, andR
2 are listed in Table 2. The results expressed in μg/g are shown in Table 3. The samples (from M1 to M12) have been sequentially collected from the bottom to the top in Pit 6 (Figures 2 and 3). The sample M13 is marl situated in contact to UE104355 where the sample M10 was taken. The sample M14 is marl belonging to another excavation sector. Also the sample M15 belongs to an archaeological stratigraphy located in a different excavation sector dated in bronze age (BC105). Observing the results (Table 3) M1 contains the lowest values of La, Ce, Nd, Eu, Tm, Yb, Sc, and Y between the samples. Gradually REE values arise sequentially between M2, M3, and M4. Samples M5, M6, M7, M8, and M9 maintain lanthanides concentrations similar to M2. Sample M10 belonging to UE104355 is, throughout the sequence, the richest unit in archaeological materials and shows the highest REE values. The sample M13 (marl) also shows REE values slightly lower than M10 probably due to the relation between this marl and UE104355 already mentioned. Samples M11 and M12 are the latest belonging to the sequence and show low REE concentrations in relation to the rest of samples except M1 which maintains the lowest levels. Sample M14 (marl) belongs to an excavation area far from Pit 6. REE values in this sample (M14) are in general similar to the sequence samples which contain the lowest values. Sample M15 dated in Bronze Age was taken in another excavation sector (BC105); this sample contains REE values that are similar to the average of the sequence (M1–M12). Finally the sequence sampled in Pit 6 reveals that the anthropogenic units (M10) and the directly related deposits (M13) show a REE enrichment.Table 2
ICP-MS analysis. Mass. Detection limits (LOD), quantification limits (LOQ), andR
2.
Elements
Mass
[Da]
LOD
LOQ
R
2
La
139
0.0004
0.0014
0.9997
Ce
140
0.0005
0.0018
0.9997
Pr
141
0.00010
0.0003
0.9997
Nd
142
0.0003
0.0001
0.9985
Sm
152
0.0003
0.0011
0.9999
Eu
151
5
E
-
05
0.00018
0.9998
Gd
158
0.00015
0.0005
0.9998
Tb
159
5
E
-
05
0.00017
0.9977
Dy
162
1.1
E
-
05
4
E
-
05
0.9998
Ho
165
3
E
-
05
0.00011
0.9983
Er
166
0.00013
0.0004
0.9999
Tm
169
1.6
E
-
05
5
E
-
05
0.9985
Yb
172
7
E
-
05
0.0002
0.9999
Lu
175
1.7
E
-
05
6
E
-
05
0.9991
Sc
45
0.013
0.04
0.9998
Y
89
0.0005
0.0016
0.9996
Rh*
103
Note: Instrumental LOD and LOQ. *Internal standard.Table 3
Rare Earth Elements sample concentrations and their standard deviations (±st/√
n).
M
Ml
M2
M3
M4
M5
M6
M7
M8
M9
M10
M11
M12
M13
M14
M15
La
10.4
13.7
17.0
16.6
13.9
14.6
12.2
14.3
13.8
18.4
13.1
12.2
17.8
12.4
13.5
±st/√
n
0.8
0.8
0.3
1.5
1.4
0.8
0.4
1.5
0.9
0.5
0.2
0.5
1.4
0.8
1.2
Ce
12.4
16.79
21.2
19.6
17.1
17.4
14.7
16.8
16.58
21.4
15.86
14.6
22.0
15.24
16.9
±st/√
n
0.9
1.09
0.2
1.7
1.9
0.9
0.6
1.8
1.09
0.4
0.14
0.6
1.6
0.96
1.4
Pr
2.2
2.9
3.42
3.4
2.9
3.0
2.54
3.0
2.8
3.71
2.58
2.56
3.6
2.60
2.9
±st/√
n
0.2
0.2
0.06
0.3
0.4
0.2
0.12
0.3
0.2
0.08
0.09
0.09
0.3
0.14
0.2
Nd
10.1
13.0
15.0
15.6
13.0
13.2
11.6
13.3
12.4
16.71
11.4
11.4
15.7
11.7
12.9
±st/√
n
0.5
0.8
0.4
1.3
1.2
0.6
0.7
1.6
0.9
0.05
0.4
0.2
1.5
0.7
1.2
Sm
1.79
2.25
2.49
2.7
2.2
2.26
2.05
2.3
2.08
2.84
2.06
1.98
2.6
2.09
2.3
±st/√
n
0.10
0.12
0.04
0.3
0.3
0.09
0.11
0.2
0.12
0.10
0.05
0.02
0.2
0.14
0.2
Eu
0.43
0.52
0.562
0.62
0.52
0.51
0.47
0.54
0.47
0.66
0.44
0.458
0.60
0.48
0.54
±st/√
n
0.03
0.05
0.006
0.07
0.07
0.02
0.02
0.05
0.03
0.02
0.04
0.010
0.04
0.04
0.04
Gd
1.83
2.22
2.40
2.69
2.21
2.20
2.04
2.30
2.03
2.88
1.88
1.94
2.42
2.07
2.22
±st/√
n
0.11
0.08
0.06
0.2
0.4
0.12
0.11
0.3
0.11
0.09
0.06
0.05
0.18
0.14
0.31
Tb
0.29
0.348
0.361
0.42
0.34
0.341
0.32
0.36
0.32
0.45
0.29
0.304
0.37
0.33
0.34
±st/√
n
0.02
0.013
0.010
0.04
0.05
0.013
0.02
0.05
0.02
0.02
0.02
0.005
0.04
0.03
0.04
Dy
1.38
1.64
1.660
2.00
1.59
1.59
1.51
1.71
1.48
2.12
1.36
1.42
1.68
1.52
1.62
±st/√
n
0.11
0.12
0.008
0.16
0.20
0.09
0.04
0.18
0.10
0.03
0.06
0.02
0.17
0.12
0.11
Ho
0.30
0.36
0.349
0.43
0.34
0.337
0.327
0.37
0.32
0.457
0.3
0.3059
0.36
0.33
0.35
±st/√
n
0.03
0.02
0.010
0.05
0.05
0.009
0.015
0.05
0.02
0.013
0.02
0.0011
0.03
0.03
0.04
Er
0.77
0.89
0.87
1.10
0.85
0.85
0.83
0.93
0.80
1.16
0.737
0.77
0.91
0.83
0.87
±st/√
n
0.04
0.06
0.06
0.11
0.14
0.05
0.06
0.11
0.03
0.04
0.009
0.02
0.08
0.05
0.10
Tm
0.107
0.1271
0.120
0.15
0.12
0.116
0.1163
0.128
0.109
0.162
0.103
0.108
0.128
0.117
0.123
±st/√
n
0.002
0.0013
0.007
0.02
0.02
0.007
0.0013
0.009
0.005
0.011
0.002
0.001
0.011
0.008
0.019
Yb
0.61
0.72
0.689
0.87
0.67
0.67
0.66
0.75
0.63
0.94
0.59
0.61
0.73
0.67
0.71
±st/√
n
0.02
0.03
0.003
0.06
0.11
0.03
0.03
0.09
0.04
0.07
0.03
0.02
0.06
0.03
0.11
Lu
0.096
0.110
0.105
0.135
0.105
0.103
0.101
0.117
0.096
0.146
0.093
0.094
0.114
0.1032
0.111
±st/√
n
0.003
0.008
0.004
0.011
0.018
0.002
0.005
0.010
0.006
0.006
0.002
0.004
0.011
0.0002
0.010
Sc
6.1
7.3
8.6
8.1
6.9
8.5
7.1
7.6
7.8
9.4
6.42
8.1
11.4
7.4
8.5
±st/√
n
0.4
0.6
0.2
0.5
0.4
0.6
0.5
0.9
1.3
0.4
0.08
1.4
0.7
0.6
0.9
Y
11.5
13.8
13.0
17.1
12.4
15.3
12.2
14.8
13.4
18.3
11.3
13.4
13.95
12.5
13.3
±st/√
n
0.8
0.8
0.3
1.2
1.4
1.3
0.7
1.7
1.2
0.8
0.2
1.2
0.98
0.9
1.2
Note: Value expressed inµg/g.


### 3.2. Relation of REE/RE
E
PAAS

Ratios La/YbPAAS, La/GdPAAS [15], La/SmPAAS, and Sm/YbPAAS [29] were used to characterize the REE/PAAS profiles in terms of relative enrichment or depletion of light REE (La-Nd), medium REE (Sm-Ho), and heavy REE (Er-Lu).Table4 shows differences in La/YbPAAS, La/GdPAAS, La/SmPAAS, and Sm/YbPAAS ratio between the samples similar to those observed in Table 3. M1 contains the lowest values of La/YbPAAS, La/GdPAAS, La/SmPAAS, and Sm/YbPAAS. Those values arise gradually trough the samples M2, M3, and M4. Samples M5, M6, M7, M8, and M9 show values similar to M2. The sample M10 (UE104355), rich in archaeological materials, presents again the highest La/YbPAAS, La/GdPAAS, La/SmPAAS, and Sm/YbPAAS values. Sample 13 (marl) related to M10 shows values slightly lower than M10. Samples M11 and M12 are the latest belonging to the sequence and show lower REE concentrations in relation to the rest of samples except M1 which maintains the lowest levels. Sample M14 (marl) values are in general similar to the sequence samples which contain the lowest values. Sample M15 (Bronze Age) contains again values similar to the average of the sequence (M1–M12).Table 4
La/YbPAAS, La/GdPAAS, La/SmPAAS, Sm/YbPAAS (REE values normalized to REEPAAS) and their standard deviations (±st/√
n).
SAMPLES
La/YbPAAS
±st/√
n
La/GdPAAS
±st/√
n
La/SmPAAS
±st/√
n
Sm/YbPAAS
±st/√
n
M1
3.7
0.3
2.2
0.2
1.9
0.14
0.63
0.04
M2
4.8
0.3
2.9
0.2
2.5
0.15
0.80
0.04
M3
6.031
0.10
3.65
0.06
3.06
0.05
0.88
0.02
M4
5.9
0.5
3.6
0.3
3.0
0.3
0.95
0.09
M5
4.9
0.5
3.0
0.3
2.5
0.3
0.80
0.09
M6
5.2
0.3
3.1
0.2
2.63
0.15
0.80
0.03
M7
4.3
0.2
2.62
0.09
2.20
0.08
0.73
0.04
M8
5.1
0.5
3.1
0.3
2.57
0.3
0.82
0.07
M9
4.9
0.3
2.96
0.18
2.48
0.15
0.74
0.04
M10
6.5
0.2
3.94
0.12
3.307
0.10
1.01
0.04
M11
4.66
0.09
2.82
0.05
2.37
0.04
0.73
0.02
M12
4.3
0.2
2.61
0.11
2.19
0.09
0.701
0.008
M13
6.3
0.5
3.8
0.3
3.2
0.2
0.92
0.09
M14
4.4
0.3
2.7
0.2
2.24
0.15
0.74
0.05
M15
4.8
0.4
2.9
0.3
2.4
0.2
0.81
0.08Summing up, M10 contains the highest values of La/YbPAAS, La/GdPAAS, La/SmPAAS, and Sm/YbPAAS and a mayor enrichment of the sequence (M1–M12). A moderate reduction in LREE, MREEs, and HREE can be observed for the intermediate samples of the sequence (M2-M3-M4-M5-M6-M7-M8-M9–M11). A mayor reduction in LREE, MREE, and HREE can be observed at the edge of the sequence, at the bottom (M1), and at the top (M12). Probably this happens because the depreciation of the filling of the Pit 6 was produced with materials similar to the baseline once the area ceased to be regularly occupied by Neolithic groups. Therefore the enrichment in REE seems to be related to stratigraphic units with a mayor anthropogenic activity. On the other hand a reduction in LREE, MREE, and HREE is detected in units with less human presence. The REE analysis of samples taken sequentially can be a useful tool to identify stratigraphic area with mayor anthropogenic impact and define the impact area of taphonomic environment. Furthermore those REE features are not clearly reproduced by mayor elements (Ca and Mg) and trace elements values in the studied samples (supplemental table, available online at http://dx.doi.org/10.1155/2014/767085), most likely due to their unintelligible, unknown, and unpredictable natural activities in the studied context.


### 3.3. Ce and Eu Anomalies

Ce
* and Eu
* anomalies are calculated as a ratio of normalized measured concentrations from the sample [30]:Ce
*
= (Ce/Ces)/(1/2 La/Las + 1/2 Pr/Prs) − 1,
Eu
* = (Eu/Eus)/(1/2 Sm/Sms + 1/2 Gd/Gds) − 1.The “*’’ and the s subscript denote the anomaly and the concentration of standard Post-Archean Australian Shale (PAAS).In this set of samplesEu
* and Ce
* values are less than one (Figure 4) and show a strong negative europium anomaly and cerium negative anomaly. That could point to more reducing conditions, without oxygen access in the archaeological site, proving stability in the studied stratigraphic units. The same Eu
* and Ce
* results have been previously obtained at the same site [1].Figure 4
Ce* and Eu* anomalies;normalized ratio<1 in samples (M1–M15) and their standard deviations.


## 3.1. Composition of Samples

15 samples have been analyzed by ICP-MS to determine REE, Sc, and Y contents. The analyzed isotopes, instrumental LOD and LOQ, andR
2 are listed in Table 2. The results expressed in μg/g are shown in Table 3. The samples (from M1 to M12) have been sequentially collected from the bottom to the top in Pit 6 (Figures 2 and 3). The sample M13 is marl situated in contact to UE104355 where the sample M10 was taken. The sample M14 is marl belonging to another excavation sector. Also the sample M15 belongs to an archaeological stratigraphy located in a different excavation sector dated in bronze age (BC105). Observing the results (Table 3) M1 contains the lowest values of La, Ce, Nd, Eu, Tm, Yb, Sc, and Y between the samples. Gradually REE values arise sequentially between M2, M3, and M4. Samples M5, M6, M7, M8, and M9 maintain lanthanides concentrations similar to M2. Sample M10 belonging to UE104355 is, throughout the sequence, the richest unit in archaeological materials and shows the highest REE values. The sample M13 (marl) also shows REE values slightly lower than M10 probably due to the relation between this marl and UE104355 already mentioned. Samples M11 and M12 are the latest belonging to the sequence and show low REE concentrations in relation to the rest of samples except M1 which maintains the lowest levels. Sample M14 (marl) belongs to an excavation area far from Pit 6. REE values in this sample (M14) are in general similar to the sequence samples which contain the lowest values. Sample M15 dated in Bronze Age was taken in another excavation sector (BC105); this sample contains REE values that are similar to the average of the sequence (M1–M12). Finally the sequence sampled in Pit 6 reveals that the anthropogenic units (M10) and the directly related deposits (M13) show a REE enrichment.Table 2
ICP-MS analysis. Mass. Detection limits (LOD), quantification limits (LOQ), andR
2.
Elements
Mass
[Da]
LOD
LOQ
R
2
La
139
0.0004
0.0014
0.9997
Ce
140
0.0005
0.0018
0.9997
Pr
141
0.00010
0.0003
0.9997
Nd
142
0.0003
0.0001
0.9985
Sm
152
0.0003
0.0011
0.9999
Eu
151
5
E
-
05
0.00018
0.9998
Gd
158
0.00015
0.0005
0.9998
Tb
159
5
E
-
05
0.00017
0.9977
Dy
162
1.1
E
-
05
4
E
-
05
0.9998
Ho
165
3
E
-
05
0.00011
0.9983
Er
166
0.00013
0.0004
0.9999
Tm
169
1.6
E
-
05
5
E
-
05
0.9985
Yb
172
7
E
-
05
0.0002
0.9999
Lu
175
1.7
E
-
05
6
E
-
05
0.9991
Sc
45
0.013
0.04
0.9998
Y
89
0.0005
0.0016
0.9996
Rh*
103
Note: Instrumental LOD and LOQ. *Internal standard.Table 3
Rare Earth Elements sample concentrations and their standard deviations (±st/√
n).
M
Ml
M2
M3
M4
M5
M6
M7
M8
M9
M10
M11
M12
M13
M14
M15
La
10.4
13.7
17.0
16.6
13.9
14.6
12.2
14.3
13.8
18.4
13.1
12.2
17.8
12.4
13.5
±st/√
n
0.8
0.8
0.3
1.5
1.4
0.8
0.4
1.5
0.9
0.5
0.2
0.5
1.4
0.8
1.2
Ce
12.4
16.79
21.2
19.6
17.1
17.4
14.7
16.8
16.58
21.4
15.86
14.6
22.0
15.24
16.9
±st/√
n
0.9
1.09
0.2
1.7
1.9
0.9
0.6
1.8
1.09
0.4
0.14
0.6
1.6
0.96
1.4
Pr
2.2
2.9
3.42
3.4
2.9
3.0
2.54
3.0
2.8
3.71
2.58
2.56
3.6
2.60
2.9
±st/√
n
0.2
0.2
0.06
0.3
0.4
0.2
0.12
0.3
0.2
0.08
0.09
0.09
0.3
0.14
0.2
Nd
10.1
13.0
15.0
15.6
13.0
13.2
11.6
13.3
12.4
16.71
11.4
11.4
15.7
11.7
12.9
±st/√
n
0.5
0.8
0.4
1.3
1.2
0.6
0.7
1.6
0.9
0.05
0.4
0.2
1.5
0.7
1.2
Sm
1.79
2.25
2.49
2.7
2.2
2.26
2.05
2.3
2.08
2.84
2.06
1.98
2.6
2.09
2.3
±st/√
n
0.10
0.12
0.04
0.3
0.3
0.09
0.11
0.2
0.12
0.10
0.05
0.02
0.2
0.14
0.2
Eu
0.43
0.52
0.562
0.62
0.52
0.51
0.47
0.54
0.47
0.66
0.44
0.458
0.60
0.48
0.54
±st/√
n
0.03
0.05
0.006
0.07
0.07
0.02
0.02
0.05
0.03
0.02
0.04
0.010
0.04
0.04
0.04
Gd
1.83
2.22
2.40
2.69
2.21
2.20
2.04
2.30
2.03
2.88
1.88
1.94
2.42
2.07
2.22
±st/√
n
0.11
0.08
0.06
0.2
0.4
0.12
0.11
0.3
0.11
0.09
0.06
0.05
0.18
0.14
0.31
Tb
0.29
0.348
0.361
0.42
0.34
0.341
0.32
0.36
0.32
0.45
0.29
0.304
0.37
0.33
0.34
±st/√
n
0.02
0.013
0.010
0.04
0.05
0.013
0.02
0.05
0.02
0.02
0.02
0.005
0.04
0.03
0.04
Dy
1.38
1.64
1.660
2.00
1.59
1.59
1.51
1.71
1.48
2.12
1.36
1.42
1.68
1.52
1.62
±st/√
n
0.11
0.12
0.008
0.16
0.20
0.09
0.04
0.18
0.10
0.03
0.06
0.02
0.17
0.12
0.11
Ho
0.30
0.36
0.349
0.43
0.34
0.337
0.327
0.37
0.32
0.457
0.3
0.3059
0.36
0.33
0.35
±st/√
n
0.03
0.02
0.010
0.05
0.05
0.009
0.015
0.05
0.02
0.013
0.02
0.0011
0.03
0.03
0.04
Er
0.77
0.89
0.87
1.10
0.85
0.85
0.83
0.93
0.80
1.16
0.737
0.77
0.91
0.83
0.87
±st/√
n
0.04
0.06
0.06
0.11
0.14
0.05
0.06
0.11
0.03
0.04
0.009
0.02
0.08
0.05
0.10
Tm
0.107
0.1271
0.120
0.15
0.12
0.116
0.1163
0.128
0.109
0.162
0.103
0.108
0.128
0.117
0.123
±st/√
n
0.002
0.0013
0.007
0.02
0.02
0.007
0.0013
0.009
0.005
0.011
0.002
0.001
0.011
0.008
0.019
Yb
0.61
0.72
0.689
0.87
0.67
0.67
0.66
0.75
0.63
0.94
0.59
0.61
0.73
0.67
0.71
±st/√
n
0.02
0.03
0.003
0.06
0.11
0.03
0.03
0.09
0.04
0.07
0.03
0.02
0.06
0.03
0.11
Lu
0.096
0.110
0.105
0.135
0.105
0.103
0.101
0.117
0.096
0.146
0.093
0.094
0.114
0.1032
0.111
±st/√
n
0.003
0.008
0.004
0.011
0.018
0.002
0.005
0.010
0.006
0.006
0.002
0.004
0.011
0.0002
0.010
Sc
6.1
7.3
8.6
8.1
6.9
8.5
7.1
7.6
7.8
9.4
6.42
8.1
11.4
7.4
8.5
±st/√
n
0.4
0.6
0.2
0.5
0.4
0.6
0.5
0.9
1.3
0.4
0.08
1.4
0.7
0.6
0.9
Y
11.5
13.8
13.0
17.1
12.4
15.3
12.2
14.8
13.4
18.3
11.3
13.4
13.95
12.5
13.3
±st/√
n
0.8
0.8
0.3
1.2
1.4
1.3
0.7
1.7
1.2
0.8
0.2
1.2
0.98
0.9
1.2
Note: Value expressed inµg/g.


## 3.2. Relation of REE/RE
E
PAAS

Ratios La/YbPAAS, La/GdPAAS [15], La/SmPAAS, and Sm/YbPAAS [29] were used to characterize the REE/PAAS profiles in terms of relative enrichment or depletion of light REE (La-Nd), medium REE (Sm-Ho), and heavy REE (Er-Lu).Table4 shows differences in La/YbPAAS, La/GdPAAS, La/SmPAAS, and Sm/YbPAAS ratio between the samples similar to those observed in Table 3. M1 contains the lowest values of La/YbPAAS, La/GdPAAS, La/SmPAAS, and Sm/YbPAAS. Those values arise gradually trough the samples M2, M3, and M4. Samples M5, M6, M7, M8, and M9 show values similar to M2. The sample M10 (UE104355), rich in archaeological materials, presents again the highest La/YbPAAS, La/GdPAAS, La/SmPAAS, and Sm/YbPAAS values. Sample 13 (marl) related to M10 shows values slightly lower than M10. Samples M11 and M12 are the latest belonging to the sequence and show lower REE concentrations in relation to the rest of samples except M1 which maintains the lowest levels. Sample M14 (marl) values are in general similar to the sequence samples which contain the lowest values. Sample M15 (Bronze Age) contains again values similar to the average of the sequence (M1–M12).Table 4
La/YbPAAS, La/GdPAAS, La/SmPAAS, Sm/YbPAAS (REE values normalized to REEPAAS) and their standard deviations (±st/√
n).
SAMPLES
La/YbPAAS
±st/√
n
La/GdPAAS
±st/√
n
La/SmPAAS
±st/√
n
Sm/YbPAAS
±st/√
n
M1
3.7
0.3
2.2
0.2
1.9
0.14
0.63
0.04
M2
4.8
0.3
2.9
0.2
2.5
0.15
0.80
0.04
M3
6.031
0.10
3.65
0.06
3.06
0.05
0.88
0.02
M4
5.9
0.5
3.6
0.3
3.0
0.3
0.95
0.09
M5
4.9
0.5
3.0
0.3
2.5
0.3
0.80
0.09
M6
5.2
0.3
3.1
0.2
2.63
0.15
0.80
0.03
M7
4.3
0.2
2.62
0.09
2.20
0.08
0.73
0.04
M8
5.1
0.5
3.1
0.3
2.57
0.3
0.82
0.07
M9
4.9
0.3
2.96
0.18
2.48
0.15
0.74
0.04
M10
6.5
0.2
3.94
0.12
3.307
0.10
1.01
0.04
M11
4.66
0.09
2.82
0.05
2.37
0.04
0.73
0.02
M12
4.3
0.2
2.61
0.11
2.19
0.09
0.701
0.008
M13
6.3
0.5
3.8
0.3
3.2
0.2
0.92
0.09
M14
4.4
0.3
2.7
0.2
2.24
0.15
0.74
0.05
M15
4.8
0.4
2.9
0.3
2.4
0.2
0.81
0.08Summing up, M10 contains the highest values of La/YbPAAS, La/GdPAAS, La/SmPAAS, and Sm/YbPAAS and a mayor enrichment of the sequence (M1–M12). A moderate reduction in LREE, MREEs, and HREE can be observed for the intermediate samples of the sequence (M2-M3-M4-M5-M6-M7-M8-M9–M11). A mayor reduction in LREE, MREE, and HREE can be observed at the edge of the sequence, at the bottom (M1), and at the top (M12). Probably this happens because the depreciation of the filling of the Pit 6 was produced with materials similar to the baseline once the area ceased to be regularly occupied by Neolithic groups. Therefore the enrichment in REE seems to be related to stratigraphic units with a mayor anthropogenic activity. On the other hand a reduction in LREE, MREE, and HREE is detected in units with less human presence. The REE analysis of samples taken sequentially can be a useful tool to identify stratigraphic area with mayor anthropogenic impact and define the impact area of taphonomic environment. Furthermore those REE features are not clearly reproduced by mayor elements (Ca and Mg) and trace elements values in the studied samples (supplemental table, available online at http://dx.doi.org/10.1155/2014/767085), most likely due to their unintelligible, unknown, and unpredictable natural activities in the studied context.


## 3.3. Ce and Eu Anomalies

Ce
* and Eu
* anomalies are calculated as a ratio of normalized measured concentrations from the sample [30]:Ce
*
= (Ce/Ces)/(1/2 La/Las + 1/2 Pr/Prs) − 1,
Eu
* = (Eu/Eus)/(1/2 Sm/Sms + 1/2 Gd/Gds) − 1.The “*’’ and the s subscript denote the anomaly and the concentration of standard Post-Archean Australian Shale (PAAS).In this set of samplesEu
* and Ce
* values are less than one (Figure 4) and show a strong negative europium anomaly and cerium negative anomaly. That could point to more reducing conditions, without oxygen access in the archaeological site, proving stability in the studied stratigraphic units. The same Eu
* and Ce
* results have been previously obtained at the same site [1].Figure 4
Ce* and Eu* anomalies;normalized ratio<1 in samples (M1–M15) and their standard deviations.


## 4. Conclusion

ICP-MS analysis has shown to be an effective technique to observe objective internal relation between samples to test our hypotheses. We consider that REE analyses are useful to identify anthropogenic patterns due to the particular chemical characteristics of those elements. Our data show that REE may be good discriminators when there are not only important, but also hidden differences between stratigraphic units. The results of samples sequentially taken in Pit 6 have revealed that the sample M10 (UE104355) contains the greatest enrichment in REE of the entire sequence (M1–M12) and confirm the anthropogenic origin of this unit. A moderate reduction in ​​ LREE, MREE, and HREE values can be observed in intermediate samples of the sequence (M2-M3-M4-M5-M6-M7-M8-M9–M11). A reduction of values on the edge of the sequence (bottom sample M1 and top sample M12) may reflect the natural origin of these units. Analysis of samples taken sequentially in a huge section (Pit 6) shows that rare earth elements may help to identify areas influenced by a higher human impact and to define the degree of anthropization on the environment also in a range of few centimetres.Poor preservation of ancient remains in archaeological stratigraphy could carry uncertainty to identify anthropogenic traces. REE analysis can indicate whether the formation of soil is due to human activities or natural processes. Many techniques may be able to shed light on ancient human group activities, but we are confident that the commonly employed methodologies can be enhanced by the use of REE. Our approach allows us to evaluate the degree of human contribution to the formation of paleosols, providing new ways to deal with one of the most classic archaeological problems.Certainly, more analyses of rare earths are needed to better characterize the composition of the parent rock and to gain more knowledge about the dynamics that take place to differentiate natural and anthropogenic soils, not only in the archaeological area of Mas d’Is but also at regional level.



---

*Source: 767085-2014-07-16.xml*