Introduction

Abandoned medieval church sites in Norway tend to lead a fairly anonymous existence. However, in 2014 the Norwegian Directorate for Cultural Heritage and Hedmark County Council received alarming information on the state of the church site at Furulund, some 15 km north of the town of Kongsvinger (Figure 1). Human skeletal remains had started surfacing as a result of ploughing, and there were concerns that the graveyard was rapidly being decimated. A subsequent site visit confirmed that the site had seen substantial erosion and that remedial action was urgently needed. The Directorate therefore sought advice from the Norwegian Institute for Cultural Heritage Research (NIKU) on how to map the church site by way of non-intrusive methods. Two methods were proposed and ultimately employed; preliminary geochemical sampling and analysis using portable Xray fluorescence (pXRF) followed by highresolution groundpenetrating radar (GPR) surveys. The aim of this paper is to illustrate how the complementary non-intrusive techniques can record archaeological features, place them in their environmental setting, assess the damage and provide sufficient information for protection strategies. To have wider relevance to archaeology, methods, approaches and results have to be contextually evaluated, and therefore this article also sets the medieval cemetery site in what is known of its archaeological and historical context, and its current position in cultural heritage management.

Figure 1 

Furulund is situated on the western banks of the River Glomma, between the town of Kongsvinger and the village of Kirkenær in Hedmark County.

Map source: The Norwegian Mapping Authority, Geovekst and Municipalities, 2018.

Reactive archaeology is often a reality for reasons of resources, need and capacity. However, when the aim is preservation in-situ over preservation by record via excavation, such as in the case of the large number of medieval graveyards in Norway, there is a need for a cost-effective, non-destructive approach. Unmarked medieval graveyards constitute a considerable challenge for the governmental bodies tasked with their protection and management, and the researchers for whom these sites constitute a valuable source of information. In Norway, up to 2,000 graveyards are believed to have existed in the Middle Ages. Of these, 647 are still in use, while a further 614 are securely attested in historical sources, but are now recorded as abandoned. The remainder are simply hinted at in toponymical and folkloric sources (Brendalsmo & Eriksson 2015). The unmarked graveyards are often found in presently cultivated fields, and more often than not, accurate locational information is lacking. Although protected by the Norwegian Cultural Heritage Act, this renders these vulnerable sites inadequately maintained and threatened by continual natural processes, agricultural activities or intentional and unintentional acts of destruction. The need to accurately locate and delimit such sites is therefore imperative both from a research and management point of view. Given the large number of potential sites, however, their mapping by way of intrusive methods is deemed unfeasible on account of the costs involved. In light of this, geophysical and geochemical methods present themselves as advantageous, both in terms of cost-efficiency and their non-destructive nature.

The Church Site and its Setting

Natural setting

The site is surrounded by arable fields and forests at the base of an alluvial valley, some 120 m from the current course of the River Glomma, a large, mature lowland river. The site is situated on top of a gentle slope caused by an older natural levee on the river’s western bank, 8 m above present water levels (Figure 2). The prevailing soils and sediments consist of alluvial silts (NGU 2017), which have formed fluvic soils, with a typical fine to coarse silt texture and a low proportion of fine sands within the laminated structure. The soils are classified as Endostagnic Fluvic Cambisols on the site and the immediate surroundings, changing to Fluvic Cambisols further afield (NIBIO 2017). Common in fluvial sediments where new flood events regularly deposit new sediment, Cambisols are young, poorly developed soils with weak structural development. The poorly drained property (endostagnic) is most apparent immediately east of the church site, where the levee slopes fairly sharply down to a small, linear depression. This was noticeably wetter, with a far higher organic content in the plough soil, suggesting consistently poorer drainage and colluvial inputs. As the river course comprises frequent meanders, oxbow lakes and palaeochannels, the site could easily sit beside such a feature, with peaty, poorly draining soils. Additionally, its close proximity to the river will cause a high water table, especially in periods of high discharge, such as spring.

Figure 2 

ALS derived DTM of the survey site and its surroundings. The survey site is located on the western bank of the River Glomma. The tributary Agnåa can be seen to the west, with numerous palaeochannels and terraces beside the current river beds.

Map source: The Norwegian Mapping Authority, Geovekst and Municipalities, 2018.

Historical background

Because of their location, rural medieval graveyards in Norway are seldom subject to archaeological investigations and, indeed, none have been investigated in their entirety. The character of their layout and of the graves themselves must therefore draw on evidence from studies of graveyards still in use, on excavations in urban contexts, medieval law texts and recorded tradition. Graveyards were always located next to a church, in an area delineated by a ditch or a simple fence, and usually measuring c. 30 × 30 m (Brendalsmo 2014: 78). Within this area shallow, rectangular graves were dug, usually disregarding the presence of older graves. The graves themselves were orientated east – west according to Christian convention and the deceased may have been interred in rectangular or trapezoidal wooden coffins, or in a simple linen shroud (Brendalsmo 1989: 27). Divisions within the graveyard were strictly governed by social standing and sex. More affluent individuals were buried closer to the church building than the less privileged, and women were buried to the south of the church and men to the north. Divisions based on geographic affiliation, where individual hamlets were assigned discrete plots have also been noted (Brendalsmo 2014: 81).

Little is known with certainty about the church that once served the small community at Furulund. Judging from a surviving limestone baptismal font, it predates the year 1200, and written sources inform us that it was in use until 1648, when it was considered too dilapidated for further upkeep, and consequently torn down. Furulund is briefly mentioned in a handful of written sources between the 14th and 17th centuries, which provide clues as to the status and, arguably, the physical size of the church. The medieval cadastres and urbariums in which the church figures, record that only a nominal tax was paid every year to the bishop and the Diocese (Brendalsmo 2015: 7). Later sources indicate that service was only held once per month. From the few written sources available, therefore, it seems the building at Furulund constituted a small rural church of modest, but perhaps not uncommon standing.

No systematic archaeological investigations have been carried out in order to locate the church or its graveyard. Map regression studies indicate that the areas surrounding the site were forested at least until the early 19th century, before they were converted to the arable land we see today. At this point the deterioration of the site must have started to accelerate, as intensive ploughing aerates and erodes topsoil and also incorporates subsoil, with the consequence of truncating and degrading buried archaeology. Unsupervised investigations by amateur metal detectorists in the 1980s recovered nails and tacks believed to stem from the church. Their distribution pattern formed a rough square which was thought to represent the original extent of the building, and this particular area has been demarcated by the landowner by retaining a small section of unploughed land in the field, measuring approximately 10 × 10 m. Rather than protecting the site, however, this has had the unfortunate effect of increasing the damage to the surrounding soil as turning the plough repeatedly along the island edges has caused intensified erosion of both topsoil and subsoil. According to local informants, human bones were exposed during the excavation of a trench for a water pipe near the assumed location of the church. This was dug without archaeological supervision and so the exact location of the disturbed graves is unknown. Representatives from the Directorate for Cultural Heritage and Hedmark County Council have previously defined the approximate extent of the site based on available evidence, and concluded that the graveyard was roughly circular in shape, measuring some 33 m in width, and 36 m in length. Despite this approximate delineation, and the fact that the site is scheduled, it currently sees annual ploughing and therefore continual truncation.

Methods

The choice of GPR over other geophysical methods was based on previous experiences where GPR had successfully detected inhumation graves in cemeteries, whereas the application of pXRF for cursory soil analysis was more experimental. Portable XRF has been successfully applied to a variety of archaeological settlement and production sites in Norway (Cannell 2017), but had never before been used on a mortuary site as a prospection tool. Geochemistry using pXRF was chosen on the assumption that the systematic mapping of certain elements across the church site would yield relatively enhanced values that would map differential land use, and thus delimit the graveyard. The combination of geophysics and geochemical prospection has been successful at delimiting settlement and activity areas in recent research (Fleisher & Sulas 2015; Booth et al. 2017). In their case study, Fleisher and Sulas found the geochemical data added activity details to geophysical results, and evidenced activities in areas void of geophysical responses, thus the complementing methods provided information on the social and practical use of space.

The use of geophysical methods to detect and map graves, clandestine or otherwise, has a well-established history, and a considerable body of literature exists on the subject (e.g., Vaughan 1986; Bevan 1991; Davenport 2001; Cheetham 2005; Jones 2008). Due to its comparatively high spatial resolution and its capability to resolve smaller targets, whilst simultaneously providing depth information, GPR is generally considered the most suitable solution for mapping inhumation burials (Conyers 2006; Jones 2008; Moffat 2015). Alternative geophysical techniques such as magnetometer surveys (e.g., Linford 2004), and electromagnetic induction (e.g., Bigman 2012) have also seen some success, particularly when used in combination (Nobes 1999; Dalan, De Vore & Clay 2010; Bigman 2014). Geophysical techniques have only recently been employed as tools for locating and mapping funerary monuments in Norway. This can largely be attributed to the introduction of the concept of high-resolution surveys, that is, surveys with a crossline resolution of 25 cm or less. Its use has generally been confined to prehistoric grave mounds (Trinks, Gansum & Hinterleitner 2010; Gustavsen 2014: 262–264; Schneidhofer et al. 2017), although medieval and later graveyards have also been investigated (e.g., Solli & Stamnes 2013; Davis et al. 2000).

Early published studies on mortuary sites using multi-element geochemical approaches, such as Keeley, Hudson and Evans (1977) and Bethell and Smith (1989), sampled graves in order to test the feasibility of identifying an inorganic chemical signature for a body. Bethell and Smith concluded that, whilst the local enhancement of certain elements could be directly associated with the presence of a body, the results were environmentally specific. This conclusion is not unique to the application of geochemistry on mortuary sites. The interplay between soil processes and the major and minor elemental enhancement from past human land use should be seen as site specific, and therefore it is the relative enhancement or depletion of certain elements that form the basis of interpretation (Oonk, Slomp & Huisman 2009; Wilson, Davidson & Cresser 2009; Vyncke et al. 2011). In Norway, single and multi-element archaeological geochemistry has largely seen use on settlement sites, where it has been used to delimit areas of activity and to determine settlement morphology (Provan 1971; Provan 1973; Bakkevig 1980; Prøsch-Danielsen 1996; Prøsch-Danielsen 2005; Cannell 2017). The application to mortuary sites, however, is more limited. Within a rescue archaeology setting, the E18 Project conducted excavations of an extensive burial mound cemetery at Gulli in Vestfold County in 2005. A total of 1850 phosphate samples were taken from the topsoil over selected areas of the cemetery, but upon excavation the data from the phosphate mapping was difficult to ascribe to the documented archaeological features (Gjerpe & Samdal 2005). Similar conclusions were reached at the Ormen Lange Project in Møre og Romsdal County in 2003, where phosphate samples were extracted from a Neolithic grave (Åstveit 2008: 272–274).

X-ray fluorescence has been applied to archaeological research for over 50 years (e.g., Shackley 2011; Charlton 2013; Frahm & Doonan 2013), although more recent technological innovation has allowed for a reduction in instrument size, to the extent where field portable and handheld instruments are widely available. The use of portable XRF on soils is relatively new, especially in archaeology. It has been used for some years for environmental monitoring (e.g., Kalnicky & Singhvi 2001; Berger, Zou & Schliecher 2009; Photos-Jones & Hall 2011; Ulmanu et al. 2011; Parsons et al. 2012), but it is still experimental to some degree and the instrument methodologically contested (e.g., Speakman & Shackley 2013). That said, it is rapidly becoming established as an efficient and effective means of analysing for multiple elements to parts per million (ppm) accuracy with minimal sample preparation. This is corroborated by a rapidly increasing volume of published studies using the instrument on archaeological soils and sediments (e.g., Donais & George 2012; Gauss et al. 2013; Hayes 2013; Grattan et al. 2014; Bissett & Claassen 2016; Booth et al. 2017).

Geochemical sampling

Geochemical sampling had to be achieved in a minimally intrusive, cost-effective manner. Using an Eijkelkamp windowed soil auger, 3 cm in diameter, the base of the plough soil was sampled, taking samples at regular intervals in transects through the area. The auger was pushed 25–30 cm into the plough soil (the average topsoil depth), and removed. A sample from the base of the plough soil was then extracted using single-use plastic equipment to reduce cross-contamination. The auger was then cleaned before the next sample was taken. Using transects allowed a greater areal coverage, with a focus upon the potential limits of the site rather than the centre.

A total of five transects were initially established over the area, with samples spaced 5 m along the transect lines, the large distance between samples necessary to keep to a restricted budget. Further sample points were subsequently added, in order to verify trends observed in the initial analyses, bringing the total number of samples covering the assumed location of the graveyard to 61 (Figure 3). The sampling was carried out in November 2014 in cold, dry conditions.

Figure 3 

GPR survey area, GPR profile lines and pXRF sample points at Furulund. The GPR area encompasses approx. 1.8 hectares. The pXRF samples were taken in five main transects across the site.

Map source: The Norwegian Mapping Authority, Geovekst and Municipalities, 2018.

Sampling occurred in the topsoil only, a layer that is regularly disturbed by ploughing and therefore no longer considered in-situ archaeological remains. The sampling, therefore does not directly disturb archaeology, although this could be deemed minimally intrusive over non-intrusive. More intrusively, two additional, deeper cores were taken approximately 100 m north and 50 m west of the main survey area, well away from any protected archaeological remains. Here the auger was hammered 1 m into the subsoil in order to determine soil structure and development, and to sample for background elemental values. Samples were taken from all identified master soil horizons (Ap, B, C).

In-situ analysis is of course possible with pXRF, however for lighter elements including P, this can produce unrepresentative data due to moisture effects, issues of surface geometry, and variable distances between the sample and the instrument exacerbating issues of signal attenuation. Such analysis is purely qualitative and cannot be reproduced. Analysis was therefore conducted on processed samples in the laboratory to produce more interpretable data. Analysis typically requires 10–20 g of dried soil, so the initial samples from the site were sub-sampled. In order to remove the negative effects of moisture on sample results (Bastos, Melquiades & Biasi 2012; Coronel et al. 2014; Schneider et al. 2016), all samples were dried at 105°C for 24 hours. As the soils at Furulund were well-sorted silts, sieving was not deemed necessary. After homogenising the sample, they were then placed in clean plastic sample cups with a 6µm polypropylene window. Analysis was undertaken using a Thermo Scientific Niton™ XLt3 GOLDD portable XRF in mining mode, resting on a stand to ensure constant distance between the sample and the instrument. The instrument can measure the elements Mg to Pb (and above, but not relevant for this study), with a resolution down to 10 ppm, depending on the element’s atomic weight. The analytical time totalled 300 seconds between all filters (50 seconds for the main filter, 100 seconds for the high, 50 for the low and 100 for the light). These settings are considered optimal for the lighter elements of interest, without taking the extra step of using helium in the pXRF (Cannell 2017). Prior to each sample run and at regular intervals, a blank sample, followed by purchased, matrix matched international standards were analysed. These were NIST 2711a, NIST 1646a, Sigma-Aldrich trace metals clay 2 and Sigma-Aldrich trace metals loamy clay 2. These are suitably similar in elemental range and matrix to allow the calibration of results and between them covered all of the relevant elements. The instrument limits of detection (LOD) were set to 2σ. The calibrated values for the selected elements were imported into a geographical information system where interpolated and gridded surfaces representing trends in the values were generated using ordinary kriging. These were then combined with other data sources for further analysis.

Ground-penetrating radar

The GPR survey was conducted in October 2015 following two weeks of increasingly cold (c. 12-0°C), and unusually dry (0–0.2 mm precipitation) weather. An area of 1.8 hectares was surveyed using a motorized 16-channel, 400 MHz MALÅ Imaging Radar Array (MIRA) from MALÅ Geoscience, covering the assumed churchyard and considerable areas to the east, west and north of this. The channel spacing was 10.5 cm and the measurements were GPS-triggered at a rate of 50 Hz, resulting in an inline resolution of approximately 5 cm, depending on the driving speed. Ideally, elongated targets such as inhumation graves should be surveyed so that the survey lines cross the targets perpendicularly, as this will increase the likelihood of imaging and defining the targets (Moffat 2015: 46). Surveys of Christian cemeteries and graveyards, therefore, should be carried out in a north-south direction. At Furulund this was not an option as, although compacted and reasonably levelled, the surface still bore signs of the previous spring’s ploughing in the form of narrow, east-west running ridges and troughs.

Once collected, the data were processed using the ApRadar software, developed by ZAMG ArcheoProspections®/LBI ArchPro. Here, trace interpolation, time-zero corrections, band-pass frequency filtering, spike removal, de-wow filters, average-trace-removal, amplitude-gain corrections, amplitude balancing and Hilbert transformations were applied. Time-to-depth conversion was based on hyperbola fitting and set to a velocity of 10 cm/ns for the upper parts of the dataset, down to 10 ns, whilst decreasing to 8 cm/ns at 20 ns and beyond. The data were then resampled to a resolution of 8 × 8 cm, and subsequently interpolated into a 3D data block from which georeferenced depth slices of 5 cm thickness were generated. Visualisation, analyses and interpretation of the depth slices was performed in ArcGIS, and subsequently combined with the pXRF and additional data sources.

Results

Ground-penetrating radar

Extending well beyond the assumed boundaries of the church site, the GPR data enables a tentative characterisation of the site’s former and present natural setting. To the east of the church site, on the slope towards the river, the dataset is characterised by substantial and sharply defined interconnected areas with reflecting and absorbing properties. The areas follow the slope closely, and provide a good indication of where the plough causes mechanical erosion of the site, as it cuts into and truncates the natural stratigraphy. In the western part of the site an extensive feature containing interleaving absorbing and reflecting material, can be observed. Crossing the site diagonally from the southwest towards the northeast, this feature likely represents an infilled palaeochannel, comparable to those observed in the vicinity in aerial LiDAR data (Figure 4).

Figure 4 

A: GPR profile 1, running across the area north of the graveyard and showing the natural stratigraphy of the area. B: GPR profile 2, running across the graveyard, indicating that the graves are cutting the natural stratigraphy – see Figure 3 for the profiles’ position in relation to the graveyard.

Scattered in the area west of the graveyard are a series of anomalies of varying sizes and shapes, but invariably decreasing in extent with depth. Their fills appear homogeneous with absorbing properties, and they are clearly defined against the more reflective subsoil background. Because of their diverging appearance, and the fact that they do not form any clear geometric patterns, the anomalies have been interpreted as tree throws formed by the mechanical removal of vegetation in the area, as it was converted to arable land. The eastern part of the site is characterised by a considerable number of drainage ditches and modern pipes. None of these appear to have disturbed the graveyard significantly, bar a single ditch diagonally crossing the area in the south-eastern corner. A north-south running ditch and pipe skirts the eastern extent of the graveyard, narrowly missing the burials.

In the GPR data a cluster of anomalies interpreted as graves has been identified immediately north of the present villa and gardens (Figure 5). The anomalies have attenuating properties and intermittently form a strong contrast to the surrounding, highly reflective subsoil. Covering this is a fairly homogeneous deposit with attenuating properties, where no archaeological features can be observed. The features interpreted as graves start appearing in the datasets as a more or less continuous, rectangular cluster at a depth of c. 60–80 cm below the surface, and individual features only become identifiable at a depth of c. 100 cm, before disappearing at c. 160 cm below the surface. The cluster is poorly defined against the surrounding background to the east, on the brink of the slope towards the river. Furthermore, due to their tight spacing or even intercutting, a secure delineation of individual graves is difficult. Those features that can be positively identified as graves are rectangular to sub-rectangular in plan. They measure between 80–250 cm in length, and 35–80 cm in width, and are orientated E–W. Because of the level of intercutting, it is difficult to determine the graves’ individual depths, but for the singular graves around the edges of the graveyard, it is estimated to c. 40 cm. All told, a total of 130 individual graves have been identified, 84 of which have been classified as ‘certain’, the remaining 46 classified as ‘possible’. The features form a rectangular area measuring c. 32 × 22 m. Its borders are clearly defined to the west and north, whereas the eastern limits are less pronounced due to erosion and the southern limits obscured by the present day gardens.

Figure 5 

A: Aggregated amplitude map of the graveyard, showing depth slices between 80 and 100 cm. B: Interpretation of the anomalies observed in the depth range 60–160 cm. C: Interpretation map of the area surrounding the graveyard.

Map source: The Norwegian Mapping Authority, Geovekst and Municipalities, 2018.

Geochemical Prospection

Of the 34 elements measured, 17 were selected for calibration. These elements are those that, from previous published studies, have produced interpretable data relating to previous human activity (Entwistle, Abrahams & Dodgshon 1998; Middleton 2004; Wilson, Davidson & Cresser 2009; Vyncke et al. 2011). Whilst some elements are often interpreted as specific to an activity, all relative enhancement or depletion are dependent upon local soil conditions (Cook & Heizer 1965; Bethell & Smith 1989; Wilson, Davidson & Cresser 2009; Cannell 2017), therefore of these 17, it was not assumed all would be interpreted as related to the graveyard location. In fact, only four elements could be directly related to the graveyard location; P and Ca via the enhancement of the local soils from skeletal remains and associated human activity, Fe from the changes in soil processes in the affected area, and Cu (Table 1). Other elements, when subject to interpolation or ordinary kriging, either spatially determined other, more recent activity, such as raised Zn by a trackway, or had no distinct spatial patterning in the sampled area. Each element (Ca, P, Fe and Cu) will be briefly considered in turn below, after a short evaluation of the elements (Pb and Sr) that were expected to be related to the graveyard, but did not produce interpretable spatial patterning.

Table 1

Descriptive statistics for each element for all 61 sample points.

Cu Fe Ca P

Maximum 204.130 30510.740 7287.790 1778.630
Minimum 0.000 22333.910 5361.330 647.800
Mean 14.418 25854.613 6100.688 1175.210
Median 0.000 25866.190 6064.720 1122.090
Standard Deviation 32.174 1767.643 430.554 258.104
Kurtosis 19.949 –0.064 –0.580 –0.123
Skewness 3.896 0.311 0.247 0.517
Range 204.130 8176.830 1926.460 1130.830
Count 61.000 61.000 61.000 61.000

In the Middle Ages, exposure to lead via utensils, medication, pigments and personal ornaments was common, however the amount varied greatly dependent upon the location and status of the individual. Lead is primarily retained in bone; over 90% of absorbed lead is stored within the bone structure (Drasch 1982). The measured values for lead (Pb) were low; up to 35.8 ppm. Whilst there was a slight elevation over the graveyard area, due to the inconsistency of enhancement over the area where bones were present on the surface, the data was not deemed to represent the graves. Strontium values are environmentally dependent, but as a small amount of Sr replaces Ca in the bone structure (Price, Burton & Bentley 2002), it was considered possible that Sr could be slightly enhanced by the presence of human bone. This was not the case, perhaps because the amounts taken up by humans via groundwater are very low, and environmentally dependent thus the local population reflect the local conditions without creating a reservoir effect.

The data for Fe shows a normal distribution, with values ranging from 22,333 to 30,510 ppm within the sampled area (Figure 6). The mean of 25,272 ppm is somewhat lower than background values, which have a mean of 26,664.08 ppm. The background values show considerable variation, from 25,525.57 to 38,723.81 ppm, which reflects the fact that iron is both soluble in anaerobic (i.e. waterlogged) conditions, and highly variable in alluvial sediments. The distribution of Fe indicates a generally lower iron concentration in the graveyard area, near respecting the graves located via GPR (see below). However it is difficult to directly assign this to past activity; it potentially a result of changing drainage caused by the nearby slope to the east of the graveyard, and the impact of graves themselves on the local drainage and Fe distribution. Hence it is not a product of the graves, but perhaps an effect caused by their presence. That said, in certain soil types, relative iron depletion due to loading of the soils with organic material has been documented (Oonk et al. 2009), and this effect can be combined with the changing drainage.

Figure 6 

Interpolated distribution of Fe (ordinary kriging), with the sample points as graduated circles overlaid over the GPR interpretations.

Copper has repeatedly been associated with past human activity, usually in conjunction with metal working (Aston, Martin & Jackson 1998; Cook et al. 2009; Carey et al. 2014), or as traces in certain ashes (Canti 2003). It is present in trace amounts in organic materials, including bones, and as a direct product of copper metal use. The distribution of copper at Furulund is concentrated in a well-defined 25 × 25 m area, fairly central to the study area (Figure 7). Background values for Cu are negligible, with one exception of 22.96 ppm. Therefore, the concentration of copper centrally within the graveyard is highly likely to be from human activity. Although the values are not high – up to 204.23 ppm – the instrument has a 10 ppm or higher resolution for Cu. The additional samples taken over the area repeat and confirm the pattern of enhancement over the area of the graveyard where the density of graves reduces, and slightly into the denser area. Because the samples are taken from plough soils, the exact location of enhancement cannot be assumed, as some degree of movement has occurred. The effects of ploughing are assumed to be spatially limited, as discussed in Ammerman (1985) and Gustavsen et al. (2017). The exact source of the copper is unknown; however, the elevated values do not correlate with the area where bones are being ploughed to the surface. It is tentatively suggested that this is the location of the small timber church, which is discussed further below.

Figure 7 

Interpolated distribution of Cu (ordinary kriging), with the sample points as graduated circles overlaid over the GPR interpretations.

One automatically associates calcium with bones; over 70% of human bone is composed of calcium phosphate crystals (Branca 1997). The measured values range from soil samples was 5,370.63 to 7,287.79 ppm (Figure 8), slightly lower than the background value range of 5,872.68–7,931.46 ppm. These have a similar range, suggesting that there is considerable natural variation in the sediments. The distribution of values in the graveyard area, however, suggests some input from the inhumations. Where the bones are visible on the surface, in the eastern part of the sampled area, the calcium values are consistently higher, reflecting the fact that weathering inputs have contributed to the soil composition. Where the burials remain below the plough depth, to the north and west, the calcium values are consistently lower. Whilst this does not define the graveyard area itself, it does indicate a large proportion of the area has burials disturbed by ploughing and plough erosion.

Figure 8 

Interpolated distribution of Ca (ordinary kriging), with the sample points as graduated circles overlaid over the GPR interpretations.

Phosphate has a long history of association with archaeology as an indicator of past human activity. Besides being present in all organic material, it also is a major component of bone. The background values vary from 422.09–1,239.61 ppm, with the topsoil values being greatly enhanced compared to the subsoil samples. The measured values in the sampled area range from 647.80 to 1,497.08 ppm (Figure 9), with the highest values being concentrated in the central and eastern areas.

Figure 9 

Interpolated distribution of P (ordinary kriging), with the sample points as graduated circles overlaid over the GPR interpretations.

The distribution of values shows two central trends: firstly that the central area of the graveyard is enhanced; an area that corresponds to the elevated Cu values. Secondly, there is a ‘halo’ of lower values immediately outside the grave area as identified via GPR. To the east this is less clear, however this area has considerable plough erosion. This ring of lower values is potentially archaeologically significant, and is discussed further below.

Discussion and Conclusions

Interpreting the GPR dataset

As is evident from this survey and others elsewhere (e.g., Schneidhofer et al. 2016), GPR data can provide valuable information on environmental setting, both past and present, and allows for a broader contextual interpretation of the site. At Furulund, the data shows that the church site is located atop a series of now backfilled palaeochannels and flood deposits in an area that was once more densely forested than it appears today. In addition, information on modern intrusions such as drainage ditches can be gleaned, as can the extent and degree of erosion on the slope to the east of the graveyard (Figure 10).

Figure 10 

The GPR data provides information on the past and present setting of the church site.

Map source: The Norwegian Mapping Authority, Geovekst and Municipalities, 2018.

The successful mapping of features and deposits in GPR data is to a large degree dependent on the relative proportions of water retained in the soil. In soils where the difference in relative dielectric permittivity (RDP) between layers is low, the radar signal will pass through with little or no reflections generated. Conversely, where the difference in water content, and therefore RDP, is high, and occurs over a short distance, high amplitude reflections will be generated at the interface (Conyers 2012: 34–40). At Furulund, the upper deposits overlying the graveyard appear to form a homogenous, attenuating surface where no archaeological features can be identified. This conflicts with the human remains observed in the topsoil, which again indicates that graves close to the surface are being truncated by ploughing. The first indications of a presence of a graveyard only appear at a depth of c. 60–80 cm, where the surrounding deposits generate relatively high amplitude reflections, while the anomalies representing graves retain low amplitude properties. A working hypothesis for Furulund, therefore, is that the deposits covering the site are fairly homogenous and well-draining, and compositionally identical to the backfill of the graves, hampering a positive identification in the upper parts of the dataset. Below these deposits, the fluvial deposits into which the graves have been cut must be structurally different, presumably retaining more water. This change is fairly abrupt, as the deposits clearly reflect considerably more of the energy in the radar signal. Constructing the original graves has disrupted this reflective surface, and must have altered its water retaining properties. The vertical grave shafts have formed free-draining pockets of soil and, as there is little structural difference between these and the overlying deposits, relatively little energy is reflected compared to the surrounding soil matrix.

The successful detectability of graves in geophysical datasets is inherently linked to the grave’s original form and its present state. This, in turn, is governed by cultural factors such as the status of the individual, mortuary traditions or depth of burial, as well as natural conditions such as local pedology, geomorphology and geology (Dick et al. 2017). Inhumation graves should, in theory, be detectable using GPR, particularly at the contact zone between grave shaft and backfill (Conyers 2006: 66), and graves containing surviving coffins should readily lend themselves to detection on account of the presence of voids, or solid remnants of the casket (Conyers 2012). Graves containing bodies simply wrapped in shrouds are naturally harder to detect. Here, the backfill needs to be significantly altered texturally or the grave shaft needs to disrupt the natural stratigraphy in order to create a distinct contrast to the surrounding soil matrix (Conyers 2006: 70). In literature regarding the detection of graves by means of radar, these are almost invariably recognised and presented as hyperbolic point-source reflections, generated from the top or sides of surviving caskets, from voids created by partially collapsed coffins within the grave (e.g., Conyers 2006: 69), or even by cavities formed by skeletal remains (Damiata et al. 2013). As the graves mapped at Furulund have backfills with low reflection properties, it is perhaps reasonable to assume that the graves were originally without caskets or with simple wrappings such as tree bark. Any such remains will have, due to the present soil conditions, decomposed to such a degree that they are now undetectable. Regardless of the original character of the graves, it is clear that their construction must have altered the natural stratigraphy of the site, thereby altering the local drainage properties.

Although poorly defined as individual features, the graves form a distinct cluster with a relatively clear outline and delineation, making an adjustment of the original interpretation possible. According to the earlier delineation, the graveyard was positioned directly north of the present villa and gardens and, to a degree, encompassing these. It was believed to be semi-circular in shape, measuring approximately 36 m in length (N–S) and 33 m in width (E–W). According to our interpretations of the GPR data, however, the graveyard is actually smaller and has a markedly different layout. It measures approximately 32 m in length, and can be securely traced to a width of 22 m. This tallies well with proportions observed in previously excavated medieval graveyards in Scandinavia (Reitan 2006: 256), although the layout of the Furulund graveyard differs somewhat. The surveys indicate that the graveyard had a rectangular shape, with relatively sharp delineations to the north and west. The eastern boundary is less clear and is damaged by the down-slope erosion, and the southern edge is obscured by modern gardens. It is, however, reasonable to assume that the linear brink of the slope has originally formed the eastern boundary of the graveyard. As for its position, the GPR survey shows that this also differs from the previous interpretation, with a centre point further north than what was previously thought. The GPR survey yielded no trace of solid building remains, neither in terms of foundations or a graveyard demarcation, presumably as a result of the intense cultivation of the area since being converted to agricultural land. Small rural churches in early medieval Norway were exclusively of timber construction supported by posts set in stone lined postholes, or later featuring sill beams resting on stone foundations (Brendalsmo 2015: 105). Neither method leaves much in terms of physical remains upon removal, and would be difficult to detect using geophysical methods following centuries of intense ploughing. If, as speculated here, that the enhanced Cu in the soil is associated with the church location, then the building can be tentatively located fairly centrally to the graveyard area.

Interpreting the geochemical dataset

In the case of Furulund, the application of archaeological geochemistry was judged to have the potential to detect the long-lasting use of the site for both a church and burials; activity which over several hundred years will have imported objects and human remains to that location. The contrast to the surrounding area, which was probably cultivated or forested but not used for consistent burials, is likely to be marked and measurable from elemental enhancement. The human body contains substantial amounts of calcium, phosphate, iron and other elements in lesser amounts (Bethell & Smith 1989), and coffins, shrouds or clothing, with or without personal ornamentation, all adds elemental enhancement to the soil within a specified, culturally, ritually and perhaps physically delineated area. The bones, and therefore potentially the enhanced soils, are now beginning to surface, as each time the plough turns the soil, it scrapes the subsoil to a greater or lesser degree, and the elemental enhancement is becoming incorporated into the plough soil. Studies have shown that the movement of objects in the soil, and by implication the soil itself by the plough, is often less than 5 m from its point of origin (Gustavsen et al. 2017), suggesting that any enhancement will be broadly representative of the graveyard area.

Early sources for Christian law in Norway such as Frostatingsloven, which in part dates to the late eleventh century, state that burial in a church graveyard was both a law and a right, with exceptions such as those who were excommunicated. As consecrated land, graveyard demarcation was necessary and stipulated by law although, as Brendalsmo notes, there are discrepancies between early written laws and practice as revealed by archaeological investigations (Brendalsmo 2000). As for the nature of the graveyard demarcation, this could have been an insubstantial construction such as a wooden fence or a simple ditch. According to the early Christian laws of Norway, rural graveyards were traditionally enclosed by wooden fences, the congregation being responsible for its upkeep (Jensenius 2000: 62), and previous excavations of rural medieval church sites have also revealed insubstantial features demarcating the churchyard (Reitan 2006: 256–258, with refs.). If such an enclosure was present at Furulund, its detection by GPR would have depended on its state of preservation and the character of the backfill. Given the soil conditions on site, and the evident truncation of other archaeological features by modern activities, however, it is unlikely to have survived as anything but an area of differentiation in use and thus invisible to GPR. As suggested in the geochemical data for the northern limit of the graveyard in the area of lower P values, particularly to the north and west, we have a possible location of the graveyard boundary, thus the geochemical data adds detail to the GPR results.

That the enhanced P correlates with the elevated Cu could suggest that this was the focus of activity in the churchyard, and likely the church itself. No trace of the church is visible in the GPR data, but here the geochemical data may offer insight. Explaining the presence of Cu is challenging as there are few comparative datasets available, however the church was at the fore of introducing and expanding monetary use in early medieval Norway (1050–1250 A.D.), and encouraged offerings in the church in coin form, a practice that became firmly established by the 1200’s. Large numbers of coins have been found distributed directly under medieval churches, as has been attested by excavation at several medieval churches (Gullbekk 2015: 237–242). The lack of other metals associated with coinage, such as silver, at Furulund, creates interpretative problems although, as this was a meagre, local church, the wealth of the congregation was perhaps limited to offering lower value coins. It should also be noted that evidence of fine smithing (bronze ingots, crucible fragments, globules) has been observed on at least one near contemporary church site (Reitan 2006). Here, fragments of bronze, tentatively interpreted as book fittings, were also found within the limits of the church itself. Given the proliferation of coins found on medieval church sites elsewhere it seems reasonable that the increased Cu values at Furulund represent coin deposition rather than the, so far, rarely observed smithing activity.

Methodological assessment

Methodologically, the surveys at Furulund have yielded valuable information on the feasibility of using GPR and soil chemical techniques for the delineation of graveyards containing inhumation graves on cultivated land.

In the process of sampling, mapping, selecting and interpreting, it must always be acknowledged the inherent bias in the process. These are caused by a plethora of issues, some more openly or easily addressed than others. On the technical side, we have the instrument parameters, detection limits, internal and empirical calibration performed, and also the sampling limitations. For example, sampling the base of the topsoil assumes a constant and even scouring of the plough into the subsoil, and even amalgamation of this material into the base of the topsoil. This is obviously not always the case, nor the assumption that any artificial or natural fertilisers used over the site have a universal and even effect. Element selection is in part based on our reliance on previous studies for reference and how the data is handled via interpolative techniques. Of course, the success of the technique is firmly established via published material, but this does not mean that each individual case study should ignore the ever present bias in the data produced. A further complication in geochemical analysis is the local environmental conditions which cannot, and should not, be compared wholesale to other case studies (Oonk, Slomp & Huisman 2009; Bethell & Smith 1989). The use of pXRF as an analytical instrument is also a source of error due to issues of sample heterogeneity, surface geometry, instrument resolution that varies per element, and issues of signal interference and peak overlap (Schneider et al. 2016; Cannell 2017). Whilst sample pre-treatment and drying can mitigate some of these issues, for pXRF the final data set is a proportional representation of the actual amounts in the sample, on a per-element basis. The data does not represent absolute values. As interpretation is always based on intra-site relative enhancement and considered contextually, this was deemed as suitable to answer the aims of the project. Future application of the method to similar sites will test if these elements can be established as more universal indicators of medieval settlement and burial. Despite these limitations in the use of pXRF, the data has produced valuable and interpretable data that adds further dimensions to the site use and extent when combined with the GPR data.

The GPR survey has demonstrated that it is possible to identify individual graves and clusters of graves, given that the present soil conditions are favourable. As observed at Furulund, the topmost 60 cm of the dataset covering the graveyard are all but void of archaeological features, whilst human remains are continually being brought to the surface by the plough. This indicates that, while graves are probably present fairly close to the surface, the geophysical contrast between these and the upper deposits is too low to produce identifiable anomalies. The cause of this lack of contrast is at present not known, but could be related to the current ploughing regime at the site, where the graves are decimated to such a degree that they are unlikely to be identified even using traditional methods such as topsoil stripping. It could, however, also be related to the original construction of the graves where they have been backfilled using material compositionally identical to that of the surrounding soils. The successful identification of graves at Furulund only occurs at depth where the deposits are more likely to retain water and where the construction of graves has disturbed the natural stratigraphy. This has formed pockets of well-drained soil which contrasts with the surrounding soil matrix. Given the particular soil conditions at Furulund, therefore, it is not possible to evaluate the potential of GPR on its own on sites of similar date and character but with set in areas with completely different soil conditions.

Overall, therefore, the combination of techniques can be used to mitigate individual weaknesses, and enhance the overall interpretation. Here, the GPR response was poor for the upper depth range due to environmental factors, whereas the geochemical data is a surface technique dependent on the site truncation by ploughing. GPR provided an insight into the composition, density and layout of the graveyard, whereas the geochemical data provided information on the location of the graveyard boundary, and the presence and potential activities within the church. The damage to the east of the site is present in both datasets, but the geochemical data, using Ca from bone, suggests the extent of it. Together, therefore, they create a more nuanced picture of the previous site use and current state, in a minimally/non-intrusive manner.