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 Northern Peloponnese is not only home of a series of ancient poleis that are being studied by archaeologists, but it is also located on the southern shoulder of the most active extensional crustal structure in the world; the Corinthian rift. This rift has shaped the Northern Peloponnese as we now see it today since the Pliocene. Normal faulting, the tectonic uplift of syn-rift sediments and sea level changes, has shaped a landscape of steps rising from the coast to the ridges in the hinterland that provides challenging conditions to a geophysical survey. Where we can find coarse grained slope and delta deposits of conglomerate on top of banks of marl on ridges and slopes, the lower marine terraces and the coastal plain as well as valleys show the protective caprock eroded and the marl covered by young alluvial deposits. These materials show only a small contrast in their magnetic properties, which reduces the importance of magnetic mapping for the archaeological prospection in this region. The human utilization of the coastal plain and the urban areas pose additional challenges. These challenges have been overcome through various approaches that are shown in exemplary case studies from Aigeira and Sikyon. Whereas a combination of magnetic mapping and ground-penetrating radar (GPR) works very well on the ridges and along the slopes where we find coarser sediments in addition to the magnetic mapping, it is not suitable in the coastal plain due to the attenuating properties of the alluvial sediment. Here, electrical resistivity tomography (ERT) proved to be very successful in mapping entire parts of a settlement in great detail. Seismic soundings were also sucessfully applied in determining the bedrock depth, the detection of walls and in the question of locating the harbor basin. In the presented six exemplary case studies, the following findings were made: (1) A fortification wall and building foundations at a depth of 0.4–1.2 m on a plateau northwest of the acropolis of Aigeira was found by 400 MHz GPR. (2) A honeycomb-shaped pattern of magnetic anomalies that suggested cavities could be identified as a weathering pattern of conglomerate rocks. (3) A rock basement 2.3 m deep and remains of an enclosing wall of the Aigeira theater area were found by shear wave refraction measurements. (4) Extensive ERT surveys detected several building remains in Sikyon like a potential building and grave monuments as well as several small houses. (5) A silted-up depression in the sediments of the coastal plane located through Love wave measurements, could be taken as evidence for either a silted harbor or a navigable riverbed.

We show an extensive multimethod geophysical study of focusing on some enigmatic subsurface structures found at Ancient Aigeira (N Peloponnes, Greece) that could be interpreted either as prehistoric chamber tombs or complex weathering patterns of the local marl–conglomerate rock sequences. It turns out that the nonseismic methods do not allow to distinguish between an archaeological and a geological origin of the observed patterns with certainty. In contrast, we demonstrate how shear‐wave seismics and full‐waveform inversion (FWI) can be used in archaeological prospection for distinguishing between these alternative essentially different interpretational models that are not distinguishable through nonseismic prospection data. The example site Aigeira is strategically well located on a hill on the Northern Peloponnese overlooking the Corinthian Gulf and has been inhabited with occupational gaps since Middle Neolithic times until the 12th to early 14th century ce. Magnetics, ground‐penetrating radar (GPR) and electrical resistivity tomography (ERT) reveal a honeycomb‐shaped anomaly pattern that could have been interpreted as a system of prehistoric chamber tombs. The time‐domain SH‐FWI strategy based on a sequential inversion of low‐pass and band‐pass filtered data results in subsurface models for shear‐wave velocity and density that accurately fits the complicated seismic data set. A highly heterogeneous subsurface is revealed that is characterized by linear cracks on a decimetre scale. The seismic FWI results are compared in detail with GPR, ERT and among each other. It turns out that the FWI result is consistent with each of these other geophysical methods but provides a more comprehensive subsurface characterization that it is supported by corings in addition. With the help of the seismic survey, we can reject the interpretation hypothesis of a prehistoric cemetery with chamber tombs and confirm that the enigmatic geophysical patterns represent a geological weathering structure that could be addressed as a reincised fan delta draped by reddish palaeosols.

Mycenaean chamber tombs are composed of air-filled burial cavities and, therefore, can be detected by non-invasive geophysical methods. In this study, an electrical resistivity tomography (ERT)-based approach was used to detect graves at three different sites in the area surrounding ancient Olympia. We collected ERT data at the necropolis of Mageiras-Kioupia, where a number of chamber tombs had already been located and partly excavated in order to ground-truth the ERT response for chamber tombs in different states of preservation. Near the village of Epitalio, a preceding archaeological survey discovered remains of a one chamber tomb in mid-slope position of a larger hill. There, we applied the ERT-approach, calibrated at Mageiras-Kioupia, to the tomb and detected ten more tomb candidates at this site. In addition to ERT, we carried out ground penetrating radar (GPR) studies at the Epitalio site. In general GPR, data are highly consistent with ERT results. We were thus able to localize an unknown, most probably Mycenaean cemetery where no systematic archaeological information previously existed. In the third area of investigation, the Dartisa plateau, no potential chamber tomb comparable to the ones excavated and measured at Mageiras-Kioupia and Epitalio were identified.

Highlights:• Horizontal-to-Vertical Spectral Ratio is used for mapping the thickness of the alluvial fill.• Recording of ambient noise with mini arrays is used to map the shear-wave velocity.• Imaged depression filled with silty material might be connected to former harbor.