Frequency-domain electromagnetic induction (FDEMI) methods are frequently used in non-invasive, area-wise mapping of the subsurface electromagnetic soil properties. A crucial part of data analysis is the geophysical inversion of the data, resulting in either conductivity and/or magnetic susceptibility subsurface distributions. We present a novel 1D stochastic optimization approach that combines dimension-adapting reversible jump Markov chain Monte Carlo (MCMC) with artificial bee colony (ABC) optimization for geophysical inversion, with specific application to frequency-domain electromagnetic induction (FDEMI) data. Several solution models of simplified model geometry and a variable number of model knots, which are found by the inversion method, are used to create re-sampled resulting average models. We present synthetic test inversions using conductivity models based on 14 direct-push (DP) EC logs from Greece, Italy, and Germany, as well as field data applications using multi-coil FDEMI devices from three sites in Azerbaijan and Germany. These examples show that the method can effectively lead to solutions that resemble the known DP input models or image reasonable stratigraphic and archaeological features in the field data. Neighboring 1D solutions on field data examples show high coherence along profiles even though each 1D inversion is independently handled. The computational effort for one 1D inversion is less than 120,000 forward calculations, which is much less than usually needed in MCMC inversions, whereas the resulting models show more plausible solutions due to the dimension-adapting properties of the inversion method.

We present a case study of a bog showing how an integrated approach of multi-method geophysical sounding and local soil sampling can be used to identify, differentiate, and map organic sediments. Our study is based on ground-penetrating radar (GPR), electrical resistivity tomography (ERT) and shear-wave seismic (SH seismic) profiling applied to sediments of the former Lake Duvensee (northern Germany), nowadays a bog. This is a well-known locality for remains from the Mesolithic hunter-gatherers’ occupation that has been attracting archaeological and geoarchaeological research for100 years. The bog is embedded in low conductive glacial sand and is characterized by layers of different gyttja sediments (detritus and calcareous). The present study was conducted in order to identify the bog morphology and the thickness of the peat body and lake sediments, in order to understand the basin evolution. To validate the geophysical results, derived from surface measurements, drilling, soil analyses as well as borehole guided wave analysis of electromagnetic waves and Direct-Push (DP-EC) have been carried out and used for comparison. It turned out that each method can distinguish between sediments that differ in grain size, particularly between peat, lake sediments (gyttjas and mud) and basal glacial sand deposits. GPR is even able to separate between strongly and weakly decomposed peat layers, which is also clear considering resistivity variations in the ERT computation. From the association between geophysical properties and sediment analysis (e.g., water content and organic matter) different gyttjas were distinguished (coarse and fine) and seismic velocity was correlated to bulk density. Moreover, GPR and SH-wave seismics present different resolutions, confirming that the latter allows measurements, which are more focused on determining the extension of basal sand deposits, the depth of which is difficult to reach with GPR. Representative values of electrical resistivity, dielectric permittivity, and shear wave velocity have been determined for each sediment type and are therefore available to complete the investigation of wetland environments. Fine grained lake sediments were difficult to differentiate by the applied methods. This could be a result of high ionic concentration within the permanent groundwater body, partly masking the sediment properties.

We investigate the landscape development of the early Mesolithic hunter-gatherer sites of Duvensee (10000–6500 cal. BCE). Based on ground-penetrating radar (GPR) and geoarchaeological drillings, we present for the first time a three-dimensional (3D) reconstruction of the palaeoenvironment of 63 ha covering subarea of the former lake during the Mesolithic. The archaeological aims were (1) to detect the location of former islands possibly hosting hunter-gatherer settlements and (2) to reconstruct the ancient landscape development for understanding prehistoric land use. The research in Duvensee lasts almost 100 years, providing vivid illustrations of early Mesolithic life. Clusters of Mesolithic camps have been found located on small sand hills that formed islands in the prehistoric lake. For this environment, we present depth maps of the three most important sedimentary facies interfaces of the ancient Lake Duvensee. Interface1 represents the transition between coarse organic sediments (peat and coarse detritus gyttja) and fine-grained organic sediments (fine detritus gyttja, calcareous gyttja), Interface2 represents the transition to the underlying clayish-loamy sediments, and Interface3 marks the top of the basal sand deposits at the lake bottom. From Interface3, we identified the location and extent of five former islands with Mesolithic camps. Stratigraphic information from the corings enabled us to create a 3D model of the spatio-temporal development of the Duvensee bog. The locations of the islands and their estimated dive-up times agree with the spatio-temporal pattern of the previous archaeological finds. The model shows where hunter-gatherers could settle and move from one island to another following the shorelines of the overgrowing lake. The 3D stratigraphic model provides growth and shrinking rates of the island and lake areas in the Mesolithic, and volumes of organic and non-organic deposited lake sediments. Besides, it provides a basis for a sustainable groundwater management needed for heritage preservation.

The shift to the early Holocene in northern Europe is strongly associated with major environmental and climatic changes that influenced hunter-gatherers’ activities and occupation during the Mesolithic period. The ancient lake Duvensee (10,000–6500 cal. BCE) has been studied for almost a century, providing archaeological sites consisting of bark mats and hazelnut-roasting hearths situated on small sand banks deposited by the glacier. No method is yet available to locate these features before excavation. Therefore, a key method for understanding the living conditions of hunter-gatherer groups is to reconstruct the paleoenvironment with a focus on the identification of areas that could possibly host Mesolithic camps and well-preserved archaeological artefacts. We performed a 16-channel MALÅ Imaging Radar Array (MIRA) system survey aimed at understanding the landscape surrounding the find spot Duvensee WP10, located in a hitherto uninvestigated part of the bog. Using an integrated approach of high-resolution ground radar mapping and targeted excavations enabled us to derive a 3D spatio-temporal landscape reconstruction of the investigated sector, including paleo-bathymetry, stratigraphy, and shorelines around the Mesolithic camps. Additionally, we detected previously unknown islands as potential areas for yet unknown dwelling sites. We found that the growth rates of the islands were in the order of approximately 0.3 m²/yr to 0.7 m²/yr between the late Preboreal and the Subboreal stages. The ground-penetrating radar surveying performed excellently in all aspects of near-surface landscape reconstruction as well as in identifying potential dwellings; however, the direct identification of small-scale artefacts, such as fireplaces, was not successful because of their similarity to natural structures.

Understanding the landscape evolution and human-environmental interaction within it is one of the key tasks of early Holocene research. As mobile hunter–gatherers leave few traces of structural organization, understanding their habitats is relevant for comprehending these people. Rarely does the spatial distribution of artifacts correspond to the real pattern of past human activity, but rather shows the pattern of identified artifacts. Geophysical investigations try to fill this gap and have been applied increasingly in archaeological prospection delivering landscape reconstruction, which are verified and fine-tuned using corings and excavations. Despite promising 3D models, a tool to predict the location of undiscovered former human presence and the conditions which influenced people to move across the landscape is not well developed. The primary goal of this paper is to present a methodology for connecting spatial patterns of past human activity based on archaeological and geophysical data. We discuss different GPR (ground-penetrating radar) facies classified at the shoreline of the former Lake Duvensee and geomorphological variables, which leads to the possibility of understanding where and why people chose preferred areas to settle on former islands. We also demonstrate that Mesolithic hunter–gatherer groups preferred dry areas with access to open water for short-term campsites and flatter and more protected areas for specialized and repeatedly occupied campsites. The cardinal orientation of a campsite seems to be secondary to the local peat over-growing process and access to water.