Abstract
Studies of past human-landscape interactions rely upon the integration of archaeological, biological and geological information within their geographical context. However, detecting the often ephemeral traces of human activities at a landscape scale remains difficult with conventional archaeological field survey. Geophysical methods offer a solution by bridging the gap between point finds and the surrounding landscape, but these surveys often solely target archaeological features. Here we show how simultaneous mapping of multiple physical soil properties with a high resolution multi-receiver electromagnetic induction (EMI) survey permits a reconstruction of the three-dimensional layout and pedological setting of a medieval reclaimed landscape in Flanders (Belgium). Combined with limited and directed excavations, the results offer a unique insight into the way such marginal landscapes were reclaimed and occupied during the Middle Ages. This approach provides a robust foundation for unravelling complex historical landscapes and will enhance our understanding of past human-landscape interactions.
Similar content being viewed by others
Introduction
In landscape archaeology, geophysical methods are increasingly being applied to conduct large area surveys1. Although these techniques allow obtaining high-resolution archaeological information at a landscape-scale2,3,4, they often neglect natural landscape variations. To fully understand the driving mechanisms behind human land-use, the integration of pedological and geomorphological information in these prospection stages is crucial.
In Europe, one of the most characteristic examples of past human-landscape interaction, is the reclamation of wetlands and forest that followed the urbanisation of the historical County of Flanders (Fig. 1b). During Medieval periods, these lowlands experienced a considerable population growth, making the County one of the most densely populated areas in Europe between the 11th and the 15th centuries5. To meet the demands of emerging cities, such as Ghent, Bruges and Ypres6, entire natural landscapes were reclaimed, transforming these into a landscape of dynamic exploitation7. The ruling Counts of Flanders spearheaded this evolution by endowing feudal lords with lands, while abbeys were deliberately installed on marginal land all over Flanders.
More recent, natural processes (e.g. flooding) and modernisation have altered the cultural landscape of medieval Flanders, leaving the nature and exact layout of these designed landscapes largely unknown8. Although social formation and economic exploitation of reclaimed landscapes during the High Middle Ages is, to some extent, documented in historical and cartographical sources, archaeological evidence remains scarce. One of the main reasons for this scarcity is that traditional investigations have focused on recovering architectural remains rather than situating the structures within their broader environmental contexts. By focusing on elite residences and monastic buildings, archaeological investigations have contributed to the biased image of medieval settlement landscapes that still prevails.
The situation was no different for the Cistercian abbey of Boudelo located in the north of the County of Flanders (Fig. 1b). In 1197, a small religious community settled in this area9, which at that time was an outback of the County, dominated by marshes and wetlands. Early historical accounts, such as those by monks from the Abbey of Clairvaux in Bourgogne (France), give testament to the harsh environmental conditions. They describe life at Boudelo as pauper and misserimus: poor and full of misery10. However, the community improved its property by cultivating the surrounding land. Between the 13th and 16th centuries, the monastic estate expanded to over 1000â ha, making the community one of the leading cultivators in the County. After religious as well as military struggles and successive floods, the monks were forced to abandon their grounds in Boudelo in 1578. They found refuge in Ghent and the monastic buildings were sold and dismantled.
Although the extensive reclamation strategy of the abbey is attested in historical records mentioning embankment, drainage, stockbreeding, extraction of peat and clay and the production of building ceramics, it remains hidden in the archaeological record. Excavations in the 1970s revealed the remnants of the cloister range and the abbey church, leading archaeologists to interpret the abbey grounds as limited to the coversand ridge11 (Fig. 1a). The border of the monastic precinct was believed to coincide with the edge of a Late Glacial palaeolake, which was and remains, a waterlogged environment dominated by peat and lime-rich lacustrine deposits (Fig. 1, Supplementary Fig. 1â2). However, in 2011 an electromagnetic induction (EMI) survey of this wetland area drastically altered this interpretation, as it unveiled traces of the abbey's outer court that once was part of the monastic precinct; a previously unknown designed landscape from which the monks directed their cultivation of the surrounding area.
The multi-receiver EMI sensor we used, allows simultaneously recording the apparent electrical conductivity (Ïa) and apparent magnetic susceptibility (κa) of four different soil volumes12. Whereas in non-saline environments Ïa can be directly related to soil texture (clay, silt, sand) and is influenced by soil organic matter and water content13,14, κa anomalies are often the result from ferrimagnetic soil enrichment15,16, the anthropogenic disturbance17 of top soils (e.g. ditches), or the heating of soil18 (e.g. hearths, bricks). This contrasts to other geophysical survey techniques, such as magnetometry19 and ground penetrating radar20, that each target only one specific variable and are inadequate to appreciate the full range of soil textural variability of surveyed sites. However, integrating electrical and magnetic soil properties in a single EMI survey allows targeting the pedological, geomorphological and archaeological variations jointly. By simultaneously mapping multiple soil volumes with a multi-receiver instrument, vertical discrimination potential is added to the electrical and magnetic data layers.
Results
The Ïa data revealed a complex environment with alternations of highly conductive lake deposits and peat accumulations (Fig. 2a). Apart from scarce sandy outcrops the measurements indicate a buried marshland influenced by long lasting, saturated conditions. Several ditches were revealed as highly conductive anomalies surrounding the entire survey area and enclosing two smaller zones. The high conductivity of the ditches is caused by high concentrations of fine-grained sediments (clay) and peat that contrast with the more resistive underlying sand. The extent and layout of the ditches and their persisting in the present-day drainage system indicates considerable medieval investment in reclaiming this wetland. The two smaller enclosed compounds appear to centre around sandy areas within the larger enclosure (S1 and S2 in Fig. 2). On the central platform within each compound, the κa measurements showed highly magnetic structures (Fig. 2b and d). At the easternmost site (S1), at least 18 features with a minimum diameter of 1â m, suggest the layout of a rectangular structure some 20â m by 30â m square (Fig. 2d). A central structure can be observed at the western site (S2) (Fig. 2b). Here, the elevated κa of the ditches indicates their infilling with brick and other magnetic debris, whereas a trail of high κa values leading north suggests displacing rubble towards the sand ridge.
To interpret and complement the geophysical data excavation trenches were laid out across the ditches and the central structure of each site (Supplementary Fig. 1â2). At both sites brick structures were found, as indicated by the κa measurements (Fig. 2bâe). These were identified as block foundations used to support a larger, most likely wooden, structure and to keep it away from the saturated soil. While a functional interpretation of S2 has been hampered by the poor conservation of the excavated building remains, the central structure at S1 has been interpreted as a large monastic barn. The morphology of the sandy outcrop, combined with elevated ground water levels caused instability of the foundations leading to the collapse of the building. The excavated material culture and construction materials indicate that both sites were occupied between the 13th and the early 14th centuries. This links their abandonment to a period of documented increased flooding in the region21.
Based on the Ïa and excavation data, we modelled the relative medieval topography of the site22,23. Although the heterogeneity of the ditch infillings, the peat layers and the targeted sandy layers limited the precision of this modelling procedure, an accurate model of the medieval surface was obtained. At 84 validation locations in the excavation trench at S2 a RMSEE of 0.37â cm and a Spearman correlation coefficient of 0.80 was obtained between the observed (Fig. 3a) and modelled (Fig. 3b) depths to the archaeological layers. Although the selection of the validation observations was not statistically random, the results confirmed the relative model correspondence to the medieval topography. A comparison of the modelled data to the observed surface along the validation transect (Fig. 3), shows the smoothing effect on the modelled data, but confirms the accurate representation of the targeted archaeological layers.
By integrating the κa data into the model, a comprehensive reconstruction of the reclaimed landscape was obtained (Fig. 4). Both S1 and S2 are located on prominently higher and sandier locations, which are connected to the bordering sand ridge in the north. However, while the barn (at S1) was built on an existing sandy outcrop, excavations showed that S2 was artificially raised (Supplementary Figure 2d) by removing sand from the ditches and the adjacent sand ridge. The design of both sites confirms a link to the nearby abbey buildings as their topographic position was chosen or altered in order to allow easy access to the sand ridge. The ditch enclosing the area was largest east of S1 with a width of up to 10â m and a depth of 1.5â m below the surface (Supplementary Figure 1d). Here, a larger drainage capacity was needed, as the ditch was part of the complex-wide drainage system. This large ditch also expressed a separation between the inside and outside worlds; a symbolic reminder of the separation between religious and secular life24. The ditches circling S2 and the remainder of S1 rarely surpassed a depth of 1â m and had a width ranging from 2â m to 8â m (Supplementary Fig. 1d and 2e). These wide ditches did not only supply additional drainage for the enclosed zones, but also formed a physical boundary that embodied aspects of identity and status25. Individual compounds, such as S1 and S2, with an (artificially) raised platform and enclosed by ditches, can be defined as late-medieval moated sites25. In addition, the surface model shows dykes neighbouring the ditches that add to the embankment and visibility of both moated sites. Remnants of these earthworks were attested in the excavation trenches, either still in situ or thrown into the moat ditches (Supplementary Figure 2e).
Discussion
Supported only by very limited invasive research, this reconstruction gives a comprehensive and unique insight into a designed medieval environment and shows, despite its transient nature, the amount of effort that was put into it. Our results also underline the discrepancies that can exist between historical information and uncovered archaeological realities that can contribute to understanding the response of past societies to social and environmental changes. The 3-D mapping of multiple soil properties, combined with limited and directed invasive research, provides a broad foundation for further geoarchaeological research. In the future, this methodology could be optimized by integrating additional physical soil properties, such as dielectric permittivity, with the Ïa and κa data. With the incorporation of the presented approach into different studies and environments, our knowledge of past human-landscape interactions can be significantly improved.
Methods
EMI survey
We used a motorized setup of a Dualem-21S EMI sensor12 to survey the study area (Fig. 1). This EMI sensor simultaneously measures Ïa and κa with four coil configurations, i.e. with receiver coils at 1â m, 1.1â m, 2 and 2.1â m from the transmitting coil, in two orientations, i.e. horizontal coplanar (HCP) and perpendicular (PRP) orientation, resulting in 1â m HCP, 2â m HCP, 1.1â m PRP and 2.1â m PRP. The survey was conducted by driving across the study area along parallel lines, 0.75â m apart whereby measurements were taken every 0.25â m. During surveying, soil temperature was recorded at a depth of 0.3â m to correct the Ïa data for temperature variations26 between survey days. The resulting measurements were drift corrected12. These data were then interpolated to a 0.1â m by 0.1â m grid using ordinary kriging27 with the software Surfer v. 10 (Golden Software). All Ïa datasets were used but only the κa measurements from the HCP configurations were considered, as the κa data of the PRP configurations had a poor signal to noise ratio.
Excavations
In August 2011 and August 2012, two excavation trenches were laid out across the most apparent Ïa and κa anomalies (Supplementary Figures 1,2). All profiles and unearthed brick structures were recorded using computer vision-based 3-D registration28.
Depth modelling
The four Ïa measurements, representative for soil volumes down to 0.5â m, 1.0â m, 1.5â m and 3.2â m below the surface, were combined to model the depth to the soil layers below the more conductive plough layer and the peaty to clayey infilling of the ditches (Fig. 3). Based on the excavation data, the morphology of the targeted soil layers was considered representative for the medieval topography. We followed the procedure from Saey et al.22 and De Smedt et al.29, based on the depth response functions of the Dualem-21S coil configurations as described by Wait30 and McNeill31. Model calibration23 was performed, using the observed depth to these layers obtained at 34 locations in the excavation trench at S1.
References
Kvamme, K. L. Geophysical surveys as landscape archaeology. American Antiquity 68, 435â457 (2003).
Gaffney, C. et al. The Stonehenge Hidden Landscapes Project. Archaeological Prospection 19, 147â155 (2012).
Keay, S. et al. The role of integrated geophysical survey methods in the assessment of archaeological landscapes: the case of Portus. Archaeological Prospection 16, 154â166 (2009).
Trinks, I. et al. Efficient, large-scale archaeological prospection using a true three-dimensional ground-penetrating Radar Array system. Archaeological Prospection 17, 175â186 (2010).
Pervenier, W. La démographie des villes du comté de Flandre aux XIVe et XVe siècles. Revue du Nord 65, 255â275 (1983).
Verhulst, A. The rise of cities in North-Western Europe. (Cambridge University Press, 1999).
Brenner, R. in From peasants to farmers? The transformation of the rural economy and society in the Low Countries (middle ages-beginning 20th century) in the light of the Brenner Debate (eds P. Hoppenbrouwers & J. L. Van Zanden) 275â338 (Brepols, 2001).
De Clercq, W. et al. Towards an Integrated Methodology for Assessing Rural Settlement Landscapes in the Belgian Lowlands. Archaeological Prospection 19, 141â145 (2012).
Abdij van Boudelo. Vol. record 8 sheet 80 r°âv° (Ghent State Archives).
Abdij van Boudelo. Vol. record 6, sub no. 50 (Ghent State Archives).
De Belie, A. De Boudelo abdij archeologisch onderzocht, Belsele. (De Belie, Alfons, 1997).
Simpson, D. et al. Evaluating the multiple coil configurations of the EM38DD and DUALEM-21S sensors to detect archaeological anomalies. Archaeological Prospection 16, 91â102 (2009).
Rhoades, J. D., Raats, P. A. C. & Prather, R. J. Effects of Liquid-phase Electrical Conductivity, Water Content and Surface Conductivity on Bulk Soil Electrical Conductivity. Soil Sci. Soc. Am. J. 40, 651â655 (1976).
Sudduth, K. A. et al. Relating apparent electrical conductivity to soil properties across the north-central USA. Computers and Electronics in Agriculture 46, 263â283 (2005).
Fassbinder, J. W. E., Stanjekt, H. & Vali, H. Occurrence of magnetic bacteria in soil. Nature 343, 161â163 (1990).
Le Borgne, E. Susceptibilité magnetique anormale du sol superficial. Annales de Géophysique 11, 399â419 (1955).
Gaffney, C. & Gater, J. Revealing the buried past. Geophysics for archaeologists. (Tempus, 2003).
Le Borgne, E. Influence du feu sur les propriétés magnétiques du sol et sur celles du schiste et du granite. Annales de Géophysique 16, 159â195 (1960).
Aspinall, A., Gaffney, C. & Schmidt, A. Magnetometry for archaeologists. (Altamira Press, 2008).
Jol, H. M. Ground penetrating radar: theory and applications. (Elsevier, 2009).
Buisman, J. Duizend jaar weer, wind en water in de Lage Landen. Vol. 1: tot 1300 (Van Wijnen, 1995).
Saey, T. et al. Reconstrucing the paleotopography beneath the loess cover with the aid of an electromagnetic induction sensor. Catena 74, 58â64 (2008).
De Smedt, P. et al. A multidisciplinary approach to reconstructing Late Glacial and Early Holocene landscapes. Journal of Archaeological Science 40, 1260â1267 (2013).
Cassidy-Welch, M. Monastic Spaces and their Meanings: Thirteenth-Century English Cistercian Monasteries. (Brepols, 2001).
Verhaeghe, F. in Medieval Moated Sites in North-West Europe. BAR International Series (eds F. A. Aberg & A. E. Brown) 127â172 (Archaeopress, 1981).
Sheets, K. R. & Hendrickx, J. M. H. Noninvasive Soil Water Content Measurement Using Electromagnetic Induction. Water Resour. Res. 31, 2401â2409 (1995).
Goovaerts, P. Geostatistics for Natural Resources Evaluation. (Oxford University Press, 1997).
De Reu, J. et al. Towards a three-dimensional cost-effective registration of the archaeological heritage. Journal of Archaeological Science 40, 1108â1121 (2013).
De Smedt, P. et al. Reconstructing palaeochannel morphology with a mobile multicoil electromagnetic induction sensor. Geomorphology 130, 136â141 (2011).
Wait, J. R. A note on the electromagnetic response of a stratified earth. Geophysics 27, 382â385 (1962).
McNeill, J. D. Electromagnetic terrain conductivity measurement at low induction numbers. Technical Note 6. 13 (Geonics Limited, Ontario, 1980).
AGIV. DHM-Vlaanderen (2001-2004): Digitaal Hoogtemodel Vlaanderen (CD-ROM). (Ondersteunend centrum GIS-Vlaanderen, 2004).
AGIV. Orthofoto's, middenschalig, kleur (2002), provincie Oost-Vlaanderen (CD-ROM). (Ondersteunend centrum GIS-Vlaanderen, 2003).
Acknowledgements
We thank Russell Palmer for his help with the text. We would also like to thank H. De Coene for granting unlimited access to his land and thank the volunteers who assisted in the excavation campaigns. Lastly, we thank Valentijn Van Parys for assisting in the geophysical surveys.
Author information
Authors and Affiliations
Contributions
P.D.S. coordinated the geophysical fieldwork and performed the geophysical data processing. P.D.S., E.M. and T.S. composed and validated the depth model. D.H. conducted the historical research. J.D.R. recorded 3-D excavation data and supplied the depth model calibration data. W.D.C. led the excavation campaign. M.V.M., P.C. and W.D.C. coordinated the research project. P.D.S., M.V.M., D.H. and W.D.C. wrote the paper. All authors commented on the results and the manuscript.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Supplementary Information
Supplementary information: figures
Rights and permissions
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivs 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-nd/3.0/
About this article
Cite this article
De Smedt, P., Van Meirvenne, M., Herremans, D. et al. The 3-D reconstruction of medieval wetland reclamation through electromagnetic induction survey. Sci Rep 3, 1517 (2013). https://doi.org/10.1038/srep01517
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep01517