Understanding Necrosol pedogenetical processes in post
HomeHome > Blog > Understanding Necrosol pedogenetical processes in post

Understanding Necrosol pedogenetical processes in post

Nov 05, 2023

Scientific Reports volume 12, Article number: 10619 (2022) Cite this article

1149 Accesses

3 Citations

6 Altmetric

Metrics details

In Archaeology much emphasis is dedicated to bone preservation, but less attention is paid to the burial soil (i.e., Necrosol), despite its crucial role in governing the geochemical environment. The interaction between human remains and sediments starts after inhumation, leading to bidirectional physico-chemical changes. To approach these complex, bidirectional processes, we sampled at high resolution (n = 46) two post-Roman wooden coffin burials (one single and another double), and the coeval paleosol (n = 20; nearby pedo-sedimentary sequence). The samples were analysed for physical (grain size, colour) and chemical (pH; LOI; elemental composition: FTIR-ATR, XRF, C, N) properties. Principal component analysis enabled to identify five main pedogenetical processes: decalcification, melanization, acidification, neoformation of secondary minerals (i.e., clays) and enrichment in phosphorus. Melanization, acidification and phosphorous enrichment seem to be convergent processes in Necrosols—irrespective of the parent material. Decalcification may be restricted to carbonate containing soil/sediments. Despite not mentioned in previous research, clay formation might also be an overall process. Compared to the local, coeval paleosol, pedogenesis in the studied burial soils was low (double burial) to moderate (single burial). Our results also emphasize the need to study the finer soil fractions, as they provide clues both on soil formation and bone diagenesis.

Necrosol is a valuable archive of pre- and post-mortem information. This term coined by Graf1 in 1986 refers to cemetery and burial soils. In the second half of the twentieth century, studies of Necrosol began to be developed, but it was not until 2004 that it was first described as: Soils formed by special human activity in cemeteries and burial ground with specific soil horizons, specific physical, chemical and biological properties (p. 110)2. Necrosol formation results from the interaction of the soil with human remains and other materials associated with the burial, the presence of a human body and skeletal remains that decompose/alter in-situ being key for this soil to be named. Changes taking place in the sediment happen both at short and long-term scales, concurrently with the body's taphonomic processes3, turning it, most of the times, into a soil of rapid formation. After burial, flesh decomposition produces chemical compounds and physicochemical reactions that modify the surrounding soil/sediment. Once the body is skeletonized the alteration persists due to the direct contact between skeleton and soil/sediment, involving the diagenesis of bones and the pedogenesis of the Necrosol. While bone diagenesis is a well-stablished topic in the research agenda, especially in archaeological sciences4,5,6,7,8, the Necrosol has been hardly addressed.

Since bones incorporate elements by absorption and release them through chemical alteration, the pedogenetical/geochemical environment of the local soil/sediment influences bone preservation7. The importance of identifying pre-mortem acquisition from post-mortem change promoted the first studies on soils from archaeological burials. In 1988, Pate and Hutton9 analysed the exchange of chemical elements between inhumated bone and its associated sediments. A year later, Pate et al.10 emphasized the importance of the geochemical properties of burial sites and proposed a protocol for soil sampling during excavation. They also recommended taking samples from excavation profiles as to compare the general soil chemistry of the site with localised conditions in areas more adjacent to skeletons. From the 1990s onwards, more studies about Necrosol physicochemical properties have been published, mainly focused on the chemical properties11,12,13,14 and organic content15,16. Studies on the inorganic chemical composition only deal with a few elements (see17). In Archaeology, the increase of phosphorous content in soil was traditionally researched as a signal of skeletal remains or to identify a burial site18,19,20,21,22,23.

Today, Necrosol research is focused on present-day cemeteries24,25,26,27 or war mass graves22,28,29,30. Metal contamination is of major concern due to the ecological impact of large inhumation cemeteries25,28,31. Another growing research line concerns forensic cases32. Recent experiments simulate burials to analyse soil properties and animal tissue decomposition under controlled conditions33,34. Current studies are getting more complex and include a larger number of chemical elements 30,35,36, combining them with the study of soil micromorphology15,37,38 and complement bone diagenesis research with soil analysis39,40.

Although Necrosol has been described more than thirty years ago, to the best of our knowledge, no specific studies deal with its pedogenetic processes. As Lambert et al.41 stated, there are numerous investigations in soil chemical composition but not too many in soils associated to inhumated bones. For archaeological contexts to date, we have neither a complete description of Necrosol characteristic properties nor a standard approach for its characterization. Thus, present research applies different methodologies, analysing only specific characteristics or discussing single properties, most of the time as complement to the study of human remains. But understanding Necrosol formation has an important archaeological meaning. As Pickering et al.16 point out, Necrosol is a valuable archive of pre- and post-mortem information. Given that bone diagenesis acts differently depending on the geochemical environment (see a summary in7), it is necessary to reach a good comprehension of the Necrosol environment to be able to understand the archaeological context of inhumations and cremations.

It may be important to approach the signal of Necrosol formation as part of the heritage they left, in particular in areas where the predominant acidic soils do not allow the preservation of skeletons. To approach these aspects, the present study aims to describe the characteristic composition and pedogenesis of Necrosol through a multi-sampling study of two post-Roman burials from NW Spain. To this end, we collected and analysed 46 samples from two burials: one containing and individual and another with two individuals; both originally in wooden coffins. Sampling was done in two transects in every individual, transvers and longitudinal; samples outsides the burials were also collected. The three skeletons were found in supine position, burials were West–East oriented (Fig. 1). We also sampled the lower soil cycle (n = 20) of a reference pedo-sedimentary sequence located 10 m from the burials. The sequence was composed of 3 main stratigraphical layers, each corresponding to a dune deposit and, thus, dominated by sands with biogenic carbonates. The two upper layers showed weak soil pedogenetical transformation and the bottom layer showed a moderate pedogenesis. The latter corresponded to the soil surface coeval of the burials and was composed of an upper buried epipedon (Ab horizon) and a lower layer of non-altered sands (C horizon)—in which the burials were made. The sampling design is described in detail in the "Material and methods" section. We also discuss to which extent the study of Necrosol can help to get insights into aspects of past societies and serve as a complementary source of information to the study of the human skeletal remains.

(A) and (B) Aerial view of A Lanzada with the burial area marked by a square and the location of the pedo-sedimentary sequence by a star in (B) (A modified from42, https://bit.ly/3FwpZrE; B modified from43, https://bit.ly/3BBqxKy). (C) Sampling design of T1 and T5 burials (image by Aerorec S.L. & Deputación de Pontevedra modified from44). (D) The pedo-sedimentary sequence and soil sampling (yellow dots; image by Aerorec S.L. & Deputación de Pontevedra modified from44). Radiocarbon datings are included in (C) and (B).

Samples collected outside the burials were lighter and had higher hue (dune sands, L*: 71.7 ± 3.2; C*: 12.7 ± 1.8; h: 80.9 ± 1.9) than samples collected inside the burials (Necrosol, L*: 61.8 ± 4.4; C*: 14.3 ± 0.9; h: 79.1 ± 0.9) (SI_Fig. 1). However, both had a larger yellow than red component (yellow, b*: 13.7 ± 1.2; red, a*: 2.5 ± 0.5). Results from grain size analysis indicates a predominance of sand fractions, being ~ 85% in Necrosol and more than 95% in samples outside the burial. In addition, Necrosol has higher content of gravel (2.76 ± 1.7%), fine sands (22.22 ± 1.5%) and silt + clay (9.79 ± 2.81%) compared to samples outside the burials (gravel:1.46 ± 2.14%; FS: 19.62 ± 3.81%; SC: 2.32 ± 1.62%) (SI_Fig. 1). pH values reflect alkaline conditions, although samples from dune sands had larger values (9.3 ± 0.1 in water, 8.8 ± 0.2 in KCl) than those of the Necrosol (8.9 ± 0.2 in water, 8.3 ± 0.2 in KCl) (SI_Fig. 1).

Both soil horizons of the paleosol (Ab-C) had different characteristics. The colour of the buried epipedon (Ab) was dark brown; being a*, b* and chromaticity higher (a*: 3.3 ± 0.8; b*: 15.7 ± 1.7; C*: 16.0 ± 1.9) than in the dune sands (C horizon) (a*: 1.7 ± 0.4; b*: 11.2 ± 1.3; C*:11.4 ± 1.3). Whereas luminosity and hue were lower in the Ab (L: 63.5 ± 1.9 and h: 78.0 ± 1.4) compared to the C horizon (L: 75.8 ± 4.0 and h: 81.4 ± 1.12) (SI_Fig. 1). Grain size showed a predominance of sands in both horizons (more than 70%), although medium sands were more abundant in the C horizon, while coarse and fine sands were more abundant in the Ab. Gravel content was three times higher in Ab (19.6 ± 6.3%) than in C (6.0 ± 9.7). The difference is much larger for the silt + clay: the content in Ab (6.4 ± 3.3%) is about ten times higher than in C (0.6 ± 0.8%) (SI_Fig. 1). pH results indicate higher alkalinity in C (9.4 ± 0.1) than in Ab (8.9 ± 0.1), being pH very homogeneous within each soil horizon.

The chemical composition is represented in SI_Fig. 1. Necrosol's LOI (0.68 ± 0.15%) and N content (0.03 ± 0.01%) were higher than those of dune sands (LOI: 0.37 ± 0.07%; N: 0.014 ± 0.006%) (SI_Fig. 1), while C showed the opposite (Necrosol: 3.36 ± 0.23%; outside burials: 3.96 ± 0.51%). Regarding the other elements, some of them (i.e., S and Si) showed high concentrations in punctual samples, while other (i.e., Al, K, or Cr) presented more homogeneous concentrations among samples. Silicon, Ca, Rb, Sr, Zr and U concentrations were higher in samples outside the burials. In contrast, P, Cu, Zn and Br concentrations were more elevated in Necrosol samples.

In the paleosol, LOI values were higher in the Ab (2.68 ± 0.57%) than in the C horizon (1.09 ± 0.37%). Nitrogen showed a similar distribution (Ab: 0.06 ± 0.02%; C: 0.022 ± 0.015%) (SI_Fig. 1) but carbon was lower in the Ab (1.7 ± 0.11%) than in the C horizon (4.19 ± 0.66%). The Ab also presented higher concentrations for Fe, Ti, Ga, Rb, Y, Pb, Th and Br, while the C horizon had higher S, Ca and Sr content. Phosphorus content was higher and Mn content lower in the Necrosol than in the Ab of the paleosol.

The average and standard deviation spectra, as well as the average spectrum of the second derivative are represented in Fig. 2. Six main absorbance areas can be observed: 3700–3400 cm−1, 2520–2510 cm−1, 1560–1300 cm−1, 1220–620 cm−1, and < 550 cm−1. The standard deviation spectra shows that variability between samples is largest in the 1560–1300 cm−1 region, very high in the regions 1050–850 cm−1 and 600–500 cm−1, and moderate to low in the regions 1200–1050 cm−1 and 3700–3600 cm−1. The second derivative spectrum (Fig. 2) enables to identified characteristic absorbances of soil components: quartz (1165, 1094, 1080, 798, 777, 693, 460 cm−1), K-feldspar (647, 535, 417 cm−1), carbonates (both calcite and aragonite, 2514, 1478–1411, 874, 859, and 712 cm−1), and clay minerals (i.e., kaolinite, 3694, 3668, 3647 and 3621, 1030, 1005, 911 cm−1); small amounts of mica (1005, 960, 527 cm−1) are also probable45,46,47. Very low absorbances around 3000–2800 and 1700–1600 cm−1 may correspond with low amounts of soil organic matter (SOM)45,48,49. While most of the vibrations in the region 1200–1050 cm−1 correspond to absorbances of silicate minerals, the standard deviation spectrum enables to identify a shoulder of moderate variability at 1200–1100 cm−1 that can be associated to vibrations of biogenic silica50,51,52,53,54.

Top to bottom: average, standard deviation and average (reversed) second derivative spectra.

We computed difference spectra55 by subtracting the average spectrum of the dune sands samples (C horizon of the paleosol) to the average spectra of the Ab samples, the samples collected outside the burials and those from T1 and T5 Necrosol (Fig. 3). Negative differences are observed in the regions 2520–2510 (peaking at 2514 cm−1), 1560–1300 (peaking at 1411 cm−1), 900–850 (peaking at 874 and 859 cm−1) and 720–710 cm−1 (peaking at 712 cm−1). Positive differences are found in the regions 3700–3400, 1400–900, 600–500 and 480–420 cm−1. Negative differences correspond to carbonates’ vibrations, while positive differences correspond to quartz, clay and other silicates—and possibly also to SOM (OH vibrations around 3400 cm−1). Samples outside the burials do not show positive differences in the clay region (Fig. 3) and only low negative values for the carbonates’ region. The other samples showed an increasing trend in both negative values and positive values following the sequence: C horizon → outside burials → inside T5 → inside T1 → Ab horizon (Fig. 3).

(A) Average spectra of the paleosol (black, Ab; yellow, C horizon), inside T1 (orange), inside T5 (pink), samples outside burials (dark brown). (B) Difference spectra compared to the dune sands, in blue the sections where the difference is positive and in grey where it is negative.

We performed a PCA using all analytical data (colour parameters, grain size, soil reaction, elemental composition, and selected absorbances of the IR data corresponding to soil components) obtained for the samples. The loadings of the 74 individual variables are in SI_Table 1. The first 5 components (73% of the total variance) contained at least a significant proportion of the variance of more than one variable and are the ones described here.

The first component, Cp1, accounted for 48.2% of the variance and showed large positive (> 0.7) loadings for kaolinite (3694, 3651, 3619, 1029, 1005, 911, 693, 647, 606, 585, 531, 423 cm-1) and OH (3424 and 3215 cm−1) vibrations, chromaticity and colour components (a* and b*), total SOM indicators (LOI, N), organically-bound elements (Br), and metal elements (Fe, Pb, Th, Ti) (SI_Table 1). Moderate (0.3–0.7) positive loadings were also shown by the silt + clay fraction, silicate IR absorbances (431 cm−1), SOM absorbances (aliphatic SOM: 2922, 2879, 2853, and 2842 cm−1), and some major (Al) and metallic (Rb, Y, Ga, Cu, and Zn) elements (SI_Table 1). Variables with large (< − 0.7) negative loadings include carbonates’ absorbances (1478, 1448, 1411, 874, 859, 712 cm−1), biogenic silica absorbances (1247, 1204, 1165, 1142, 1114 cm−1), hue (h) and luminosity (L*), soil reaction (pHw, pHk), total C, Ca, S and Sr, and medium and fine sands (SI_Table 1). Moderate (-0.43 to -0.65) negative loadings were also found for quartz absorbances (1081 and 1094 cm−1) and U.

Cp2 accounted for 9.9% of the total variance, coarse sand and Zr showed a large positive loading, while P had a large negative loading (SI_Table 1). Moderate positive loadings were also found for many elements (Y, Rb, Mn, Si, Fe, Nb, N, Br, Pb), carbonates (i.e., aragonite, 859 and 712 cm−1), L*, soil reaction (pHw), and aliphatic SOM (2922 cm−1). Moderate negative loadings were found for fine soil components (silt + clay), silicates (quartz and kaolinite, 532, 449, 431 and 423 cm−1), Ca and Cu (SI_Table 1).

Cp3 (9.3% of the total variance) is dominated by large to moderate positive loadings of quartz (1094, 1081, 798, 777, 462, 431 cm−1), aliphatic SOM (2922, 2879, 2853, 2842 cm−1), biogenic silica (1165, 1142, 1114 cm−1), K and Rb. Carbonates (i.e., calcite and aragonite, 1411, 874, 859, 712 cm−1), clays (i.e., kaolinite, 3694, 911, 647, 606, 585, 531 cm−1) and S, have moderate negative loadings (SI_Table 1).

Cp4 (3.5% of total variance) shows no large loadings. Moderate positive loadings were obtained for Al, Ti, K, Si, U and colour components (a*, b* and chromaticity) and moderate negative loadings for Cr, Zn and hue (SI_Table_1). Cp5 (2% of total variance) is dominated by the anti-covariation of Cr, Si and U against Zn and Fe (SI_Table 1).

Figure 4 shows the components’ scores of the samples of the burials (outside and inside), and the paleosol. Positive Cp1 scores are found for the Ab and most of the samples from inside T1. The C horizon and samples collected outside the burials show negative loadings. Regarding Cp2, all samples from the paleosol have positive scores while almost all samples collected inside and outside the burials show negative scores, exception made of a few samples outside the burials (Fig. 4). Cp3 scores show an increase in value from the C horizon to the buried epipedon, to decrease again at the top of this horizon; samples from Necrosol show no clear trend, although those collected in the skulls have large negative scores while the rest have positive or slightly negative scores. No pattern is observed for Cp4 and Cp5 scores, both in the paleosol (Ab and C horizon) and burial samples (Fig. 4).

PCA scores of the analysed samples (burials, Ab-soil, and dune-C horizon), for the five main principal components. The burial samples have been grouped according to the place they were taken (out outside the burial), the transect (lg longitudinal) and the anatomical region (sk skull, sh shoulders, pel pelvis) with which they were related.

The extracted principal components represent the five major geochemical signals that characterize the the paleosol (Ab-C horizons; Arenosol), the Necrosol and the samples outside the burials (i.e., corresponding to the C horizon of the Arenosol) (see Fig. 5). The results indicate that the first three components are the ones that capture the main differences between the burials’ Necrosol and the contemporaneous soil, while the other two seem to reflect small-scale heterogeneity in the distribution of certain soil components (i.e., mica and metals).

Diagram of the Necrosol pedogenesis compared with that of the buried epipedon of the Arenosol.

The first principal component, Cp1, is the one that accounts for a larger number of soil properties. Properties with positive loadings are essentially related to the enrichment in secondary minerals (kaolinite and silt + clay fractions), SOM (total SOM, N, aliphatic SOM) (SI_Table 1) and changes in soil reaction (i.e., pH). Properties with negative loadings are related to the abundance of biogenic carbonates (calcite, aragonite, C, Ca and Sr) and biogenic silica. Thus, the component is reflecting the main pedogenetic processes occurring in the palesol and the Necrosol (Fig. 5): decalcification (i.e., weathering of biogenic carbonates), melanization (SOM accumulation), and neoformation of secondary minerals (i.e., clay formation) (see56). The progression of pedogenesis is also accompanied by changes in other soil properties, as a decrease in soil luminosity and hue, and an increase in chromaticity (and both redness and yellowness), increased soil acidity (i.e., lower pH), and increased concentrations of metals and organically bound elements (Al, Fe, Ti, Rb, Cu, Ga, Y, Pb, Zn, Br) (SI_Table 1). Carbon, usually correlated to N content in soil due to both being basic constituents of the SOM, is here decoupled to SOM content because the inorganic C of the biogenic carbonates (e.g., from shells) dominates the C pool. In a similar way, S seems to be more dependent on the biogenic carbonates that on SOM.

The scores of Cp1 can be taken as a measure of the degree of pedogenesis: negative loadings indicating no or only slight pedogenetic evolution and positive loadings indicating more intense pedogenetic evolution. The dune sands (C-horizon of the paleosol and samples collected outside the burials) have a very similar composition and the largest negative scores. This indicates no or weak pedogenesis, a fact that means that both are parent material of the Ab horizon and the Necrosol respectively. The paleosol shows the typical evolution in soil profiles of increasing pedogenetic transformation from the parent material to the epipedon (i.e., the Ab horizon, in our case) (Fig. 5). The lower values at the top two samples of the paleosol are due to the transition to the overlying dune cycle that buried the epipedon (Fig. 4).

Cp1 scores show much lower variation in Necrosol, suggesting a lower degree of pedogenesis than that found for the Ab; but clear differences can be observed (Fig. 4). In T1 values are almost all positive, while in T5 values are almost all slightly negative. The intensity of pedogenesis in T1 is comparable to that of the lower half of the Ab, and in T5 it is comparable to that observed in the transition between the C horizon and the Ab (Fig. 4). In T1, samples collected in the shoulders transect show the weakest pedogenesis. No clear differences were found in T5. We should keep in mind that while pedogenetic transformations can be similar between the Arenosol and the Necrosol, the origin of both soils is not alike. In the first case, it was a surface layer, while the Necrosol is always a buried layer that contains a body in decomposition/alteration (Fig. 5). Our results point to pedogenetic convergence given the same parental material.

Comparing the main pedogenetic processes occurring in the Necrosol studied by us with others around the world, all share similar processes. Although ours is the first study of the Necrosol that used quantitative colour (with a colorimeter), melanization has been also found in previous studies. Fiedler et al.22 and Majgier and Rahmonov12 recorded darker colours in the Necrosol using Munsell scale (by human-eye). Other authors only indicate that the Necrosol is dark-coloured16,36. The dark colour of the Necrosol can be observed in the rectangular shape around the skeletons, and the finding of iron nails in the borders, pointing to the use of wooden coffins.

The alkalization of Necrosols has been suggested as a key factor for skeletal remains preservation5. However, all samples from the burials (T1 and T5) are less alkaline than those of the parent material (C-horizon), pointing to some degree of acidification -as it happens in the Ab of the paleosol. Although bone chemical weathering may contribute to soil alkalinization, the acidity generated by the enrichment in organic matter may counter-balance the effect through intense chemical weathering and leaching of the biogenic carbonate (i.e., decarbonation). This process seems to have been slightly more intense in the thoracic area of individual L006, which could be linked to his more confined decomposition (see "Material and methods" section). Here, the organic matter may have caused a greater acidification of the soil in contact with this skeleton. Acidification of the Necrosol has been observed by other researchers11,12,13,30,57 in soils originally alkaline, and was also linked to the increase in organic matter and the decrease in calcium content (decalcification). To approach the influence of soil pH on bone conservation, it must be also determined in the pre-burial soil or sampling surrounding areas, since the pH of the Necrosol is highly influenced by the soil-skeleton exchange over time.

Metal enrichment has been described in other investigations28,35,36,57, which found an increase in Al, Fe, Pb, Zn and Rb in the Necrosols. Keeley et al.58, Amuno & Amuno28 and Charzyński et al.57 found higher Cu concentrations inside burials. However, Charzyński et al.57 and Pankowská et al.35, in investigations of Necrosols with cremated remains from Nazi concentration camps, suggested that Cu concentrations could respond to other post-depositional processes.

In our study, the enrichment in metals and elements associated with organic matter (i.e., N and Br) is accompanied by an increase of fine (i.e., silt and clay) soil fractions. Most metal elements are enriched in the fine fractions of soils and sediments59,60,61. To date, not much attention has been paid to the differences in the content of silt and clay fractions inside and outside the burials. To our knowledge, only our previous investigation in Álvarez-Fernández, et al.62 analysed the silt and clay fraction apart from the fine earth. In this study, we encouraged the analysis of fine fractions, when studying sandy soil/sediments in particular, as they are the most reactive and have the largest potential to contain information about the interaction between buried bodies and burial environment. However, this is the first study in which an enrichment in silt and clay was observed inside the burials.

The second principal component is dominated by the inverse distribution of coarse sands and Zr (probably reflecting the content in zircon) regarding P content. This component accounts for large differences between the Arenosol and the Necrosol (Fig. 4). The Necrosol is enriched in P and shows a predominance of smaller grain sizes than the Arenosol. Thus, phosphorous enrichment seems to be a specific process occurring in the Necrosol (Fig. 5), which is certainly related to the chemical weathering of the human remains (both soft tissues and bone). Differences are also observed between burials, since T1 shows a larger enrichment than T5. As it was the case for the intensity of pedogenesis, in T1 the lowest P enrichment is observed in the shoulders transect, while distribution is more homogeneous in T5. It is somewhat surprising that the burial containing two individuals, and thus more initial body mass for P release, has lower average P concentrations than the burial with only one individual. Two complementary aspects should account for this result: (1) T1 is a more confined environment and a larger proportion of the products of body decomposition may have accumulated within the burial area; (2) the individual buried in T1 was a senile woman (> 60 years), affected by age-related osteoporosis, which may have resulted in enhanced bone chemical weathering.

Phosphorus is an element traditionally used in Environmental Archaeology to detect the presence of human activity63,64,65,66 and to identify burials19. Although there are other elements that humans transfer to the soils, phosphorus is the least susceptible to change and leaching, and therefore lasts over time67,68. For burials, phosphorous and calcium are the major elements of the bone mineral component (i.e., hydroxyapatite) that are released to the surrounding soil due to bone alteration. The process underlying Cp2 is probably the result of this release as shown by the correlation between phosphorus and calcium and silt + clay. The correlation with calcium could be explained because when phosphorus is released into the soil it binds to calcium, iron or aluminium to form stable inorganic compounds67. However, we only see correlation with calcium and not with iron and aluminium, most likely because of the alkaline nature of the parent material66 and the wide availability of calcium. These compounds help to keep phosphorus in the soil, but it has also been documented that in sandy soils, such as the Necrosol studied here, leaching does occur. What is different about the Necrosol that allows phosphorus to persist? The most likely component responsible for that is the finer soil fraction (i.e., clay), that is the most reactive. Numerous studies describe the enrichment of phosphorus in Necrosols12,14,22,28,35,36 and it seems an important specific feature of Necrosol pedogenesis.

Component Cp3 captures a pedogenetic signal related to the accumulation of resistant minerals and resistant SOM compounds that is more intense in the Necrosol than in the Arenosol. In the paleosol, a clear trend (Fig. 4) of increasing accumulation of resistant minerals and SOM is observed from the C horizon to the upper part of the Ab—exception made of the upper two samples, for reasons already commented in the degree of pedogenesis (Fig. 4). This trend is consistent with an increase in the degree of pedogenesis. But, while the Ab shows low positive values, Cp3 scores in the burials’ soils are larger both outside and inside the burial, except for the samples taken in and around the skulls (T1 and T5) and the pelvis transect (of T1). This secondary enrichment in aliphatic SOM, which is more resistant to biological degradation than other organic compounds (as proteins and polysaccharides), is probably related to the decomposition of the body soft tissues and of the wooden coffins. At this stage, we do not have a proper explanation for the secondary enrichment in resistant minerals in the burials. But we cannot rule out that it derives from the digging and burying operations in post-Roman time—mixing of sand layers, deposition of external materials in the dune sands while the burials were open, etc.

Components Cp4 and Cp5 are related to the content in some minerals (probably micas) and metals (Zn, Fe, Cr, U). Although scores tend to be negative inside the burials, reflecting a somewhat lower content of mica, Cr and U, but higher content in Zn and Fe, distribution patterns are not clear and may reflect micro-scale heterogeneity in the composition of the soil/sediment and in the geochemical conditions driving the release, mobility and accumulation of the metal elements (see for example7).

This investigation aimed to characterize Necrosol pedogenesis and composition in two burials from post-Roman times (AD sixth century). Our results provide information on the main five processes that took place: decalcification, acidification, melanization, neoformation of secondary minerals and enrichment in phosphorus. As a result of these processes the Necrosol acquires some characteristics that describe it as having darker colour, lower alkalinity and higher content of fine particles, organic matter and phosphorus. Translated into archaeological information, by comparing the chemical composition of the Necrosol and the surrounding soil we determined that these two inhumations from A Lanzada were dug in the dune and the same material was used to cover the wood coffins; therefore, there was no reconditioning of the burial site. The decomposition of wood coffins and the individuals’ bodies together with bone diagenesis triggered the formation of the Necrosol.

Although many studies indicate the increase in phosphorus content as a key aspect of necropolises’ soils, our study reveals that other pedogenetic processes must also be considered. Melaninization and acidification are very characteristic of the Necrosol from the studied burials but also in the buried epipedon (Ab), which points to a common soil process related to the increase in SOM. Higher phosphorous content and increase in metals seem to be related to the diagenesis of the skeletons—in consequence, we would expect the same processes to occur when soil is in contact with the bone. Each pedogenetic process separately can be observed in other soil types, but their combination together with accumulation of resistant minerals and resistant SOM compounds are, in our opinion, the key for a Necrosol in sandy parent materials to be described. Although our results are promising and provide a suitable explanation for other observations made in the literature, more studies are needed considering different chronological frames and parent materials to better understand Necrosol formation. Our findings suggest that the study of the Necrosol complements the characterization of human remains, being both relevant to understand the burial ritual.

The site of "A Lanzada" (Noalla, Sanxenxo) is located in the province of Pontevedra, NW Spain (UTM 51023814X; 4697448.14Y) (Fig. 1). It has been the subject of numerous archaeological campaigns from the 1950s, the last one in 2016–201744. Archaeological remains point to a wide chronology of occupation from the eighth century BC (Late Bronze Age) to AD tenth century (Middle Ages)69,70,71,72. One of the most remarkable archaeological features of A Lanzada site is its necropolis (Fig. 1), with two well-defined funerary areas dated from Roman and post-Roman times72. The last archaeological campaign (2016–2017) was focused on the East area of the site (Fig. 1) in which, adjacent to a monumental size house (a possible religious building), two burials were excavated (by OLC and AMC): T1, with one individual (L01), and T5, with two individuals (L06 and L07) (Fig. 1). L01 was a senile age female (L01, > 60 years-old), while L06 was a male adolescent (13–20 years-old) and L07 a mature male (40–60 years-old). The three skeletons were buried in dune sands and a clear brownish colour, rectangular shape pattern was observed embedding them (Fig. 1). This fact, together with the presence of iron-made nails and slight displacements at join articulations, have been archaeologically interpreted as they were buried in wooden coffins. As it can be seen in Fig. 1, the skeletons were in supine position with stretched legs. L001 is West–East oriented with arms stretched along the body. L006 was also West–East oriented with arms stretched and hands over the pelvic area. L007 was East–West oriented, and his arms were crossed over the abdomen. Open space decomposition is noticeable in L001 and L007 due to the rotation of the innominate bones and the head of the femurs. L006 skeleton is more confined against the lateral of the coffin and most of his joins did not rotate.

We collected 46 samples from within both burials, using a multisampling design organised in two transects—longitudinal and transverse—along every individual. Additional samples (5) were also taken from inside the crania and the innominate bone area (Fig. 1). To contextualize the evolution of the site, a pedo-sedimentary sequence (SQ1) (Fig. 1), located 10 m away from the burial area, was also collected. For this study we selected 20 soil samples (5 cm in thickness) corresponding to the soil (i.e., paleosol) contemporaneous of the burials. The selected samples comprise two horizons: (1) a buried epipedon (Ab) that represents the surface of the soil that was in use when the burials were excavated, (2) the dune sands or underlying parent material (C horizon), that represent the layer where inhumations were placed. In a previous work7, we have classified the soils in the archaeological site as Haplic Arenosol (calcaric)73.

To determine the chronology of the paleosol and the burials, 14C dating was carried out. One soil sample (silt + clay fraction) of the sequence (SQ1.S27) and two samples of left ribs collagen of the individuals L01 and L06 were analysed. The results are presented in Fig. 1 and SI_Table 2 and indicate that the individuals were buried by the end of the sixth century AD and the buried epipedon may have started its formation in the first century AD.

All samples (fine earth, < 2 mm) were analysed for physical (grain size, colour) and chemical (pH, LOI, elemental composition: C and N, XRF, FTIR-ATR) properties. Grain size analysis was performed using a set of sieves to determine the percentage of coarse sand (2–0.5 mm), medium sand (0.5–0.2 mm), fine sand (0.2–0.05 mm) and silt + clay (< 0.05 mm). Soil/sediment reaction (pH) was measured in water (pHw) and KCl (pHKCl) suspensions (ratio 1:2.5) with a pHmeter74. Colour was measured using a CR-5 Konica Minolta colourimeter, using the CIELab colour space that provides five parameters: L* (luminosity), a* (green–red component), b* (blue–yellow component), C* (chromaticity), and h (hue). Loss on ignition (LOI) was obtained to estimate soil organic matter content, by heating the samples at 105 °C for 24 h and then ashing at 550 °C for 5 h in a muffle furnace. Carbon and nitrogen content was measured using a LECO-TruSpec CHNS analyser; while the concentrations of P, S, Si, Al, Fe, Ti, K, Ca, Ga, Rb, Sr, Y, Zr, Nb, Cr, Mn, Ni, Cu, Zn, Pb, Th, U and Br were determined by X-ray fluorescence. The XRF equipment used were calibrated with standard reference materials and both are hosted at the RIAIDT facility of the Universidade of Santiago de Compostela, Spain. Quantification limits were 0.01% for major elements (Si, Al, Fe, Ti, K and Ca), 100 μg g−1 for P, S and Mn, 1 μg g−1 for Ga, Rb, Sr, Y, Zr, Nb, Cu, Zn, Th, Ni, Cr, U and Br; for Pb quantification limit was 0.5 μg g−1.

Spectra were acquired in finely milled samples, in the mid-infrared (MIR) region (4000–400 cm−1) by attenuated total reflectance (ATR), using an Agilent Technologies Cary 630 spectrometer. Resolution was set to 4 cm−1 and each spectrum is the average of 200 scans. The equipment was thoroughly cleaned before each measurement and a background was collected before every sample. Average, standard deviation and second derivative spectra and peak identification (based on the second derivative spectrum) were obtained with {andurinha} R package75. Assignment of compounds related to vibrations and classes is based on literature (see references in "Secondary enrichment in resistant minerals and SOM" section), considering the limitations imposed on IR interpretation of complex samples49,76,77,78.

Principal components analysis (PCA) was performed on the set of properties analysed (74 variables in 66 samples), in correlation mode and using a non-rotated solution. One factor ANOVA was run on the components scores, grouping according to the different groups that represent soils/horizons (inside and outside the burials, buried epipedon and dune sands). Since grain size data and elemental concentrations are a case of close data79, we applied a centred log-ratio (clr) transformation prior to statistical analyses80. Statistical analyses were made using software R81 (package {andurinha}75) and SPSS Statistics 23.

All data generated or analysed during this study are included in this published article in the supplementary files.

Graf, A. Flora and vegetation der Friedhölfe in Berlin (West). Verh. Ber. Bot. Ver. 5, 1–210 (1986).

Google Scholar

Sobocka, J. Necrosol as a new antropogenic soil type. In Proceedings, Soil Anthropization VII (ed. Sobocka, J.) 107–112 (2004).

Grupe, G. & Harbeck, M. Thaphonomic and diagenetic processes. In Handbook of Palaeoanthropology (eds Henke, W. & Tattersall, I.) 417–439 (Springer, 2015). https://doi.org/10.1007/978-3-540-33761-4_7.

Chapter Google Scholar

Nielsen-Marsh, C. M. & Hedges, R. E. M. Patterns of diagenesis in bone I: The effects of site environment. J. Archaeol. Sci. 27(12), 1139–1150. https://doi.org/10.1006/jasc.1999.0537 (2000).

Article Google Scholar

Nielsen-Marsh, C. M. et al. Bone diagenesis in the European Holocene II: Taphonomic and environmental considerations. J. Archaeol. Sci. 34, 1523–1531. https://doi.org/10.1016/j.jas.2006.11.012 (2007).

Article Google Scholar

Smith, C. I., Nielsen-Marsh, C. M., Jans, M. M. E. & Collins, M. J. Bone diagenesis in the European Holocene I: Patterns and mechanisms. J. Archaeol. Sci. 34(9), 1485–1493. https://doi.org/10.1016/j.jas.2006.11.006 (2007).

Article Google Scholar

López-Costas, O., Lantes-Suárez, O. & Martínez Cortizas, A. Chemical compositional changes in archaeological human bones due to diagenesis: Type of bone vs soil environment. J. Archaeol. Sci. 67, 43–51. https://doi.org/10.1016/j.jas.2016.02.001 (2016).

Article CAS Google Scholar

Kendall, C., Eriksen, A. M. H., Kontopoulos, I., Collins, M. J. & Turner-Walker, G. Diagenesis of archaeological bone and tooth. Palaeogeogr. Palaeoclimatol. Palaeoecol. 491, 21–37. https://doi.org/10.1016/j.palaeo.2017.11.041 (2018).

Article Google Scholar

Pate, F. D. & Hutton, J. T. The use of soil chemistry data to address post-mortem diagenesis in bone mineral. J. Archaeol. Sci. 15(6), 729–739. https://doi.org/10.1016/0305-4403(88)90062-3 (1988).

Article Google Scholar

Pate, F. D., Hutton, J. T. & Norrish, K. Ionic exchange between soil solution and bone: Toward a predictive model. J. Appl. Geochem. 4(3), 303–316. https://doi.org/10.1016/0883-2927(89)90034-6 (1989).

Article CAS ADS Google Scholar

Charzyński, P., Bednarek, R., Switoniak, M. & Zolnowska, B. Ekranic tecnosols and Urbic technosols of Torun necropolis. Geologija 53(4), 179–185. https://doi.org/10.6001/GEOLOGIJA.V53I4.1905 (2011).

Article Google Scholar

Majgier, L. & Rahmorov, O. Selected chemical propierties of necrosols from the abandoned cemeteries Sławobo and Szymonka (Great Mazurian Lakes District). Bull. Geogr. Phys. Geogr. Ser. 5(1), 43–55. https://doi.org/10.2478/v10250-012-0003-8 (2012).

Article Google Scholar

Majgier, L. & Rahmorov, O. Necrosols of cemeteries in Masurian Lakeland. In Technogenic Soils of Poland (eds Charzyński, P. et al.) 95–109 (Polish Society of Soil Science, 2013).

Google Scholar

Vélez, S. et al. Caracterización de Necrosoles en el bosque seco tropical en Colombia: Un acercamiento desde la geología forense. Rev. Mex. Cienc. Geol. 37(2), 146–156. https://doi.org/10.22201/cgeo.20072902e.2020.2.1539 (2020).

Article Google Scholar

Burns, A. et al. Micromorphological and chemical investigation of late-Viking age grave fills at Hofstaðir, Iceland. Geoderma 306, 183–194. https://doi.org/10.1016/j.geoderma.2017.06.021 (2017).

Article CAS ADS Google Scholar

Pickering, M. D. et al. Signatures of degraded body tissues and environmental conditions in grave soils from a Roman and an Anglo-Scandinavian age burial from Hungate, York. J. Archaeol. Sci. 99, 87–98. https://doi.org/10.1016/j.jas.2018.08.007 (2018).

Article CAS Google Scholar

Clinton, C. K., Duncan, C. M., Shaw, R. K., Jackson, L. & Jackson, F. L. C. Identification of trace metals and potential anthropogenic influences on the historic New York African BurialGround population: A pXRF technology approach. Sci. Rep. 9, 18976. https://doi.org/10.1038/s41598-019-55125-7 (2019).

Article CAS PubMed PubMed Central ADS Google Scholar

Solecki, R. Notes on soil analysis and archaeology. Am. Antiq. 16(3), 254–256. https://doi.org/10.2307/276788 (1951).

Article Google Scholar

Dietz, E. F. Phosphorus accumulation in soil of an Indian habitation site. Am. Antiq. 22, 405–409. https://doi.org/10.2307/276142 (1957).

Article Google Scholar

Provan, D. M. J. Soil phosphate analysis as a tool in archaeology. Nor. Archaeol. Rev. 4(1), 37–50. https://doi.org/10.1080/00293652.1971.9965134 (1971).

Article Google Scholar

Farswan, Y. S. & Nautiyal, V. Investigation of phosphorus enrichment in the burial soil of Kauman, Mid-central Himalaya, India. J. Archaeol. Sci. 24(3), 251–258. https://doi.org/10.1006/jasc.1996.0108 (1997).

Article Google Scholar

Fiedler, S., Berger, J., Stahr, K. & Graw, M. Localisation of a mass grave from the Nazi Era: A case study. In Criminal and Environmental Soil Forensics (eds Ritz, K. et al.) 303–314 (Springer, 2009).

Chapter Google Scholar

Tallón Armada, R., López-Costas, O. & Martínez Cortizas, A. Análisis del contenido en fósforo en los suelos y sedimentos del yacimiento de Ventosiños (Coeses) como una alternativa al hallazgo de restos óseos. In Un yacimiento ceremonial en la transición del Bronce al Hierro: Coeses (Lugo) (eds Cano Pan, J. A. et al.) (Arqueoloxía do Noroeste SLU, 2015).

Google Scholar

Spongberg, A. L. & Becks, P. M. Organic contamination in soils associated with cemeteries. J. Soil Contam. 9(2), 87–97. https://doi.org/10.1080/10588330008984177 (2000).

Article CAS Google Scholar

Jonker, C. & Olivier, J. Mineral contamination from cemetery soils: Case study of Zandfontein Cemetery, South Africa. Int. J. Environ. Res. Public Health 9(2), 511–520. https://doi.org/10.3390/ijerph9020511 (2012).

Article CAS PubMed PubMed Central Google Scholar

Prestes da Silva, R. B. et al. Contamination of heavy metals in soils under cemetery occupation in Amazonas, Brazil. Soil Sediment. Contam. 29(2), 192–208. https://doi.org/10.1080/19320383.2019.1696280 (2020).

Article Google Scholar

Mohammed, M. A. & Abudeif, A. M. Use of the geophysical approaches for studying the environmental impact assessment of the human burying techniques to the soil and groundwater: A case study of Geheina cemeteries, Sohag, Egypt. J. Afr. Earth Sci. 172, 104010. https://doi.org/10.1016/j.jafrearsci.2020.104010 (2020).

Article Google Scholar

Amuno, S. A. & Amuno, M. M. Geochemical assessment of two excavated mass graves in Rwanda: A pilot study. Soil Sediment. Contam. 23, 144–165. https://doi.org/10.1080/15320383.2013.786021 (2014).

Article CAS Google Scholar

Żychowski, J. Geological aspects of decomposition of corpses in mass graves from WW1 and 2, located in SE Poland. Environ. Earth Sci. 64, 437–448. https://doi.org/10.1007/s12665-010-0867-x (2011).

Article CAS Google Scholar

Źychowski, J. Selected elements in the soils covering mass graves from World Wars I and II in Southeastern Poland. Minerals 11(3), 275. https://doi.org/10.3390/min11030275 (2021).

Article CAS ADS Google Scholar

Amuno, S. A. Potential ecological risk of heavy metal distribution in cemetery soils. Water Air Soil Pollut. 224, 1435. https://doi.org/10.1007/s11270-013-1435-2 (2013).

Article CAS ADS Google Scholar

Santillana, E., Cordero, J. C. & Alamilla, F. Forensic soil analysis: Case study of looting at a Roman-Visigothic Burial Vault. In Soil in Criminal and Environmental Forensics. Soil Forensics (eds Kars, H. & van den Eijkel, L.) 45–60 (Springer, 2016). https://doi.org/10.1007/978-3-319-33115-7_4.

Chapter Google Scholar

Carter, D. O., Yellowlees, D. & Tibbett, M. Cadaver decomposition in terrestrial ecosystems. Naturwissenschaften 94, 12–24. https://doi.org/10.1007/s00114-006-0159-1 (2007).

Article CAS PubMed ADS Google Scholar

Szelecz, I., Koenig, I., Seppey, C. V. W., Le Bayon, R. C. & Mitchell, E. A. D. Soil chemistry changes beneath decomposing cadavers over a one-year period. Forensic Sci. Int. 286, 155–165. https://doi.org/10.1016/j.forsciint.2018.02.031 (2018).

Article CAS PubMed Google Scholar

Pankowská, A., Moník, M. & Nechvátal, M. Reading the silhouettes of burnt dead: Using elemental analysis (pXRF) to identify late bronze and early iron age urn cenotaphs. Anthropologie 56(1), 39–52. https://doi.org/10.26720/anthro.17.08.28.1 (2018).

Article Google Scholar

Asare, M. O. et al. Human burials can affect soil elemental composition for millennia—Analysis of necrosols from the Corded Ware Culture graveyard in the Czech Republic. Archaeol. Anthropol. Sci. 12, 255. https://doi.org/10.1007/s12520-020-01211-1 (2020).

Article Google Scholar

Angelucci, D. E. Geoarchaeological insights from a Roman age incineration feature (ustrinum) at Enconsta de Sant’Ana (Lisbon, Portugal). J. Archaeol. Sci. 35(9), 2624–2633. https://doi.org/10.1016/j.jas.2008.04.020 (2008).

Article Google Scholar

Lang, C. The Hidden Archive of Historical Human Inhumations Locked Within Burial Soils. Ph.D. thesis (University of York, 2014).

Cascant, M. M. et al. Prediction of alkaline earth elements in bone remains by near infrared spectroscopy. Talanta 162, 428–434. https://doi.org/10.1016/j.talanta.2016.10.071 (2017).

Article CAS PubMed Google Scholar

Viani, A., Machová, D., Mácová, P., Mali, G. & Velemínský, P. Bone diagenesis in the medieval cemetery of Vratislavs’ Palace in Prague. Archaeol. Anthropol. Sci. 13, 39. https://doi.org/10.1007/s12520-021-01286-4 (2021).

Article Google Scholar

Lambert, J. B., Simpson, S. V., Buikstra, J. E. & Charles, D. K. Analysis of soil associated with woodland burials. Adv. Chem. 205, 97–113. https://doi.org/10.1021/ba-1984-0205.ch006 (1984).

Article CAS Google Scholar

Punta A Lanzada, O Grove (Galicia, Spain) 42° 25′ 44.61″ N 8° 52′ 29.31″ W elev 16 m eye alt 585 m. Google Earth. July 18, 2020. March 20, 2021. https://bit.ly/3FwpZrE.

A Lanzada site (Galicia, Spain) 42° 25′ 44.64″ N 8° 52″ 29.42″ W elev 16 m eye alt 549 m. Google Earth. July 18, 2020. October 12, 2021. https://bit.ly/3BBqxKy.

Rodríguez Martínez, R. M. Informe valorativo arqueológico para la recuperación patrimonial do xacemento de A Lanzada (Sanxenxo, Pontevedra). Fase II. (Diputación provincial de Pontevedra, 2017).

Tinti, A., Tugnoli, V., Bonora, S. & Francioso, O. Recent applications of vibrational mid-infrared (IR) spectroscopy for studying soil components: A review. J. Cent. Eur. Agric. 16(1), 1–22. https://doi.org/10.5513/JCEA01/16.1.1535 (2015).

Article Google Scholar

Müller, C. M. et al. Infrared attenuated total reflectance spectroscopy: Aninnovative strategy for analyzing mineral components in energy relevant systems. Sci. Rep. 4, 6764. https://doi.org/10.1038/srep06764 (2014).

Article CAS PubMed PubMed Central Google Scholar

Martínez Cortizas, A., López-Merino, L., Silva-Sánchez, N., Sjöström, J. K. & Kylander, M. E. Investigating the mineral composition of peat by combining FTIR-ATR and multivariate analysis. Minerals 11, 1084. https://doi.org/10.3390/min11101084 (2021).

Article CAS ADS Google Scholar

Artz, R. R. G. et al. FTIR spectroscopy can be used as a screening tool for organic matter quality in regenerating cut over peatlands. Soil Biol. Biochem. 40(2), 515–527. https://doi.org/10.1016/j.soilbio.2007.09.019 (2008).

Article CAS Google Scholar

Simonescu, C. M. Applications of FTIR spectroscopy in environmental studies. In Advanced Aspects of Spectroscopy (ed. Farrukh, M. A.) 49–84 (IntechOpen, 2012). https://doi.org/10.5772/48331.

Chapter Google Scholar

Liu, X., Colman, S. M., Brown, E. T., Minor, E. C. & Li, H. Estimation of carbonate, total organic carbon, and biogenic silica content by FTIR and XRF techniques in lacustrine sediments. J. Paleolimnol. 50, 387–398. https://doi.org/10.1007/s10933-013-9733-7 (2013).

Article CAS ADS Google Scholar

Meyer-Jacob, C. et al. Independent measurement of biogenic silica in sediments by FTIR spectroscopy and PLS regression. J. Paleolimnol. 52, 245–255. https://doi.org/10.1007/s10933-014-9791-5 (2014).

Article ADS Google Scholar

Vogel, H., Meyer-Jacob, C., Thöle, L., Lippold, J. A. & Jaccard, S. L. Quantification of biogenic silica by means of Fourier transform infrared spectroscopy (FTIRS) in marine sediments. Limnol. Oceanogr. Methods 12(14), 828–838. https://doi.org/10.1002/lom3.10129 (2016).

Article Google Scholar

Melucci, D., Zappi, A., Morozzi, P., Giglio, F. & Tositti, L. ATR-FTIR spectroscopy, a new non-destructive approach for the quantitative determination of biogenic silica in marine sediments. Molecules 24(21), 3927. https://doi.org/10.3390/molecules24213927 (2019).

Article CAS PubMed Central Google Scholar

Leiva-Dueñas, C. et al. Long-term dynamics of production in western Mediterranean seagrass meadows: Trade-offs and legacies of past disturbances. Sci. Total Environ. 754, 142117. https://doi.org/10.1016/j.scitotenv.2020.142117 (2021).

Article CAS PubMed ADS Google Scholar

Margenet, A. J., Calderón, F. J. & Parikh, S. J. Limitations and potential of spectral subtractions in Fourier-Transform Infrared Spectroscopy of soil samples. Soil Sci. Soc. Am. J. 80(1), 10–26. https://doi.org/10.2136/sssaj2015.06.0228 (2015).

Article CAS Google Scholar

Bockheim, J. G. & Gennadiyev, A. N. The role of soil-forming processes in the definition of taxa in soil taxonomy and the world soil reference base. Geoderma 95(1–2), 53–72. https://doi.org/10.1016/S0016-7061(99)00083-X (2000).

Article ADS Google Scholar

Charzyński, P., Markiewicz, M., Majorek, M. & Bednarek, R. Geochemical assessment of soils in the German Nazi concentration camp in Stutthof (Northern Poland). Soil Sci. Plant Nutr. 61, 47–54. https://doi.org/10.1080/00380768.2014.1000232 (2015).

Article CAS Google Scholar

Keely, H. C. M., Hudson, G. E. & Evans, J. Trace elements content of human bones in various states of preservation. The soil silhouette. J. Archaeol. Sci. 4, 19–24 (1977).

Article Google Scholar

Fuge, R. Sources of halogens in the environment, influences on human and animal health. Environ. Geochem. Health 10, 51–61 (1988).

Article CAS Google Scholar

Schuetz, L. Atmospheric mineral dust—Properties and source markers. In Paleoclimatology and Paleometeorology: Modern and Past Patterns of Global Atmospheric Transport (eds. Leinen, M. & Sarnthein M.). NATO ASI Series (Series C: Mathematical and Physical Sciences), vol 282, 359–384 (1989).

Taboada, T., Martínez-Cortizas, A., García, C. & García-Rodeja, E. Particle-size fractionation of titanium and zirconium during weathering and pedogenesis of granitic rocks in NW Spain. Geoderma 131, 218–236. https://doi.org/10.1016/j.geoderma.2005.03.025 (2006).

Article CAS ADS Google Scholar

Álvarez-Fernández, N., Martínez Cortizas, A., García-López, Z. & López-Costas, O. Approaching mercury distribution in burial environment using PLS-R modelling. Sci. Rep. 11, 21231. https://doi.org/10.1038/s41598-021-00768-8 (2021).

Article PubMed PubMed Central ADS Google Scholar

Sánchez, A., Cañabate, M. L. & Lizcano, R. Phosphorus analysis at archaeological sites: An optimization of the method and interpretation of the results. Archaeometry 38(1), 151–164. https://doi.org/10.1111/j.1475-4754.1996.tb00768.x (1996).

Article Google Scholar

Scudder, S. J., Foss, J. E. & Collins, M. E. Soil science and archaeology. Adv. Agron. 57, 1–76. https://doi.org/10.1016/s0065-2113(08)60922-0 (1996).

Article CAS Google Scholar

Holliday, V. T. & Gartner, W. G. Methods of soil P analysis in archaeology. J. Archaeol. Sci. 34(2), 301–333. https://doi.org/10.1016/j.jas.2006.05.004 (2007).

Article Google Scholar

Kolb, M. F. Analysis of carbon, nitrogen, pH, phosphorus, and carbonates as tools in geoarchaeological research. In Encyclopedia of Geoarchaeology. Encyclopedia of Earth Sciences Series (ed. Gilbert, A. S.) (Springer, 2017). https://doi.org/10.1007/978-1-4020-4409-0_11.

Chapter Google Scholar

Proudfoot, B. The analysis and interpretation of soil phosphorus inarchaeological contexts. In Geoarchaeology (eds Davidson, D. A. & Shackley, M. L.) 93–113 (Westview Press, 1976).

Google Scholar

Eidt, R. C. Detection and examination of anthrosols by phosphate analysis. Science 197, 1327–1333 (1977).

Article CAS ADS Google Scholar

Blanco Freijeiro, A., Fusté Ara, M. E. & Garcia Alen, A. La necrópolis galaico-romana de La Lanzada (Noalla, Pontevedra), I. Cuad. Estud. Gallegos 16, 141–158 (1961).

Google Scholar

Blanco Freijeiro, A., Fusté Ara, M. E. & Garcia Alen, A. La necrópolis galaico-romana de La Lanzada (Noalla, Pontevedra), II. Cuad. Estud. Gallegos 22, 5–23 (1967).

Google Scholar

Filgueira Valverde, J. E. & Blanco Freijeiro, A. Excavaciones de La Lanzada. Noticiario Arqueológico Hispánico V 1956–1961, 137–152 (1962).

Google Scholar

López-Costas, O. Taphonomy and burial context of the Roman/post-Roman funerary areas (2nd to 6th centuries AD) of A Lanzada, NW Spain. Estud. Quat. APEQ 12, 55–67. https://doi.org/10.30893/eq.v0i12.111 (2015).

Article Google Scholar

IUSS Working Group, W.R.B. World Reference Base for Soil Resources 2006, First Update 2007. World Resources Reports Food and Agriculture Organization of the United Nations. FAO, No. 13 (2007).

Hendershot, W. H., Lalande, H. & Duquette, M. Soil reaction and exchangeable acidity. In Soils Sampling and Methods of Analysis, 2nd edn, (eds Carter, M. S. & Gregorich, E. G.) 173–178. https://doi.org/10.1201/9781420005271.ch16 (Taylor & Francis Group, 2008).

Chapter Google Scholar

Álvarez Fernández, N. & Martínez Cortizas, A. Andurinha: Make Spectroscopic Data Processing Easier. R package version 0.0.2. https://CRAN.R-project.org/package_andurinha. (2020) .

Coates, J. Interpretation of infrares spectra, a practical approach. In Encyclopedia of Analytical Chemistry (ed. Meyers, R. A.) 10815–10837 (Wiley, 2000).

Google Scholar

Socrates, G. Infrared and Raman Characteristic Group Frequencies, Tables and Charts 362 (Wiley, 2001).

Google Scholar

Larkin, P. J. IR and Raman Spectroscopy, Principles and Spectral Interpretation 228 (Elsevier, 2011).

Google Scholar

Reimann, C., Filzmoser, P., Garret, R. & Dutter, R. Improving Data Behaviour for Statistical Analysis: Ranking and Transformations. In Statistical Data Analysis Explained: Applied Environmental Statistics with R (eds Reimann, C. et al.) 167–180 (Wiley, 2008).

Chapter Google Scholar

Egozcue, J. J., Pawlowsky-Glahn, V., Mateu-Figueras, G. & Barceló-Vidal, C. Isometric Logratio transformations for compositional data analysis. Math. Geol. 35, 279–300. https://doi.org/10.1023/A:1023818214614 (2003).

Article MathSciNet MATH Google Scholar

R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2021).

Google Scholar

Download references

Present research was funded by "Estudo de esqueletos humanos e de secuencias edafosedimentarias do xacemento de A Lanzada (2017-CP035)" by Deputación Provincial de Pontevedra. We thank Deputación de Pontevedra, Museum of Pontevedra and the Dirección Xeral de Patrimonio da Xunta de Galicia for providing access to the archaeological soil samples. Special thanks go to the director or archaeological campaign Rafael Rodríguez Martinez for his support in all studies related to A Lanzada. This project is funded by Grupos de Referencia Competitiva (ED431C 2021/32) by Xunta de Galicia. Authors would like to thank the use of RIAIDT-USC analytical facilities. ZGL is funded by Deputación provincial A Coruña. NAF is funded by the project Fomentar a actividade investigadora do persoal investigador finalista nas convocatorias de axudas da ERC do H2020 by GAIN (2021-CP052). OLC is funded by Ramón y Cajal 2020 (RYC2020-030531-I), JIN project (PID2019-111683RJ-I00) Spanish Ministerio de Ciencia e Innovacion and Beca Leonardo a Investigadores y Creadores Culturales 2020 (2020-PO048) de la Fundación BBVA.

CRETUS, EcoPast (GI-1553), Universidade de Santiago de Compostela, 15782, Santiago de Compostela, Spain

Zaira García-López, Antonio Martínez Cortizas & Noemi Álvarez-Fernández

CRETUS, EcoPast (GI-1553), Area of Archaeology, Department of History, Universidade de Santiago de Compostela, 15782, Santiago de Compostela, Spain

Olalla López-Costas

Archaeological Research Laboratory, Stockholm University, Wallenberglaboratoriet, 10691, Stockholm, Sweden

Olalla López-Costas

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Z.G.L., A.M.C., O.L.C. contributed to the conception and design of the work. Z.G.L., A.M.C., N.A.F. and O.L.C. contributed to the acquisition of the data. Z.G.L., A.M.C., O.L.C. contributed to the analysis and interpretation of the data. Z.G.L., A.M.C., N.A.F. and O.L.C. contributed to the draft and revised it.

Correspondence to Zaira García-López.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

García-López, Z., Martínez Cortizas, A., Álvarez-Fernández, N. et al. Understanding Necrosol pedogenetical processes in post-Roman burials developed on dunes sands. Sci Rep 12, 10619 (2022). https://doi.org/10.1038/s41598-022-14750-5

Download citation

Received: 08 March 2022

Accepted: 13 June 2022

Published: 23 June 2022

DOI: https://doi.org/10.1038/s41598-022-14750-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.