Search

Abstract The UNESCO World Heritage Wadden Sea holds remains of a medieval cultural landscape shaped by interactions between man and natural forces. From the Netherlands to Denmark, human efforts of cultivating low-lying areas created a unique coastal landscape. Since the Middle Ages, storm floods widely drowned embanked cultural land and especially affected North Frisia (Germany), where once fertile marshland was permanently turned into tidal flats. One key region, the Edomsharde, was widely destroyed in 1362 AD. Medieval settlement remains still occur in the tidal flats around the island Hallig Südfall and are commonly associated with Edomsharde’s trading centre Rungholt—ever since a symbol for the region’s drowned landscapes and focus of this study. We present a first-time comprehensive reconstruction of this medieval settlement by means of new geophysical, geoarchaeological and archaeological data. Our results reveal remains of up to 64 newly found and rectified dwelling mounds, abundant drainage ditches, a seadike, and especially the discovery of Edomshardes’s main church as important landmark in this former cultural landscape. These finds together with the documented imported goods confirm a thriving society, involved in transregional trade and thereby close a significant gap in medieval history not only for North Frisia, but the entire Wadden Sea region.

Study Region: This study investigates the on- and offshore fresh to saline groundwater distribution as well as the submarine groundwater discharge within a coastal karstic aquifer system in Antalya (Türkiye). Study Focus: Groundwater acts as a critical link between on- and offshore environments, connecting freshwater systems to the saline oceans. With an increasing population along coastlines, understanding coastal groundwater interactions is essential. In this study we utilize on- and offshore hydrogeophysical techniques to image variations in subsurface porewater salinities. Offshore continuous resistivity profiling and hydroacoustic measurements are complemented by land based electromagnetic measurements. New hydrological insights for the region: Our findings show how groundwater discharges into the sea through coastal and submarine freshwater springs while simultaneously facing the risk of saltwater intrusions nearby. Diffuse groundwater flow dominates along the main beach area of Antalya, while channelized flow occurs primarily along the cliff coastline where submarine springs occur. The springs create distinct salinity anomalies in the water column and remove the sediment cover from the carbonate aquifer, exposing hard rock on the seafloor and creating different habitats. The multi-method geophysical approach provides a detailed picture of freshwater-saltwater interactions and highlights regions of varying groundwater discharge into the sea. Our findings contribute valuable insights for water management strategies in Antalya, with implications for safeguarding freshwater resources and mitigating the risk of saltwater intrusions.

Hyperbolae in radargrams are caused by a variety of small subsurface objects. The analysis of their curvature enables the determination of propagation velocity in the subsurface, which is important for exact time-to-depth conversion and migration and also yields information on the water content of the soil. Using deep learning methods and fitting (DLF) algorithms, it is possible to automatically detect and analyze large numbers of hyperbola in 3D Ground-Penetrating Radar (GPR) datasets. As a result, a 3D velocity model can be established. Combining the hyperbola locations and the 3D velocity model with reflection depth sections and timeslices leads to improved archaeological interpretation due to (1) correct time-to-depth conversion through migration with the 3D velocity model, (2) creation of depthslices following the topography, (3) evaluation of the spatial distribution of hyperbolae, and (4) derivation of a 3D water content model of the site. In an exemplary study, we applied DLF to a 3D GPR dataset from the multi-phased (2nd to 12th century CE) archaeological site of Goting on the island of Föhr, Northern Germany. Using RetinaNet, we detected 38,490 hyperbolae in an area of 1.76 ha and created a 3D velocity model. The velocities ranged from approximately 0.12 m/ns at the surface to 0.07 m/ns at approx. 3 m depth in the vertical direction; in the lateral direction, the maximum velocity variation was ±0.048 m/ns. The 2D-migrated radargrams and subsequently created depthslices revealed the remains of a longhouse, which was not known beforehand and had not been visible in the unmigrated timeslices. We found hyperbola apex points aligned along linear strong reflections. They can be interpreted as stones contained in ditch fills. The hyperbola points help to differentiate between ditches and processing artifacts that have a similar appearance as the ditches in time-/depthslices. From the derived 3D water content model, we could identify the thickness of the archaeologically relevant layer across the whole site. The layer contains a lot of humus and has a high water retention capability, leading to a higher water content compared to the underlying glacial moraine sand, which is well-drained.

An investigation of the Magdalenian occupation at Étiolles-Les Coudray (France) was conducted using geophysical methods. Based on ground-penetrating radar (GPR), electrical resistivity tomography (ERT), electromagnetic induction (EMI), and stratigraphic information, we present a reconstruction of the main sediment interfaces carrying the Magdalenian occupation. Étiolles-Les Coudray is one of the most important open-air campsites in the Paris Basin, where consecutive settlements distributed along the Hauldres stream were preserved by silts. The geoarchaeological goals were, in particular, the reconstruction of the ancient environment in which hunter–gatherers settled, providing spatialized known stratigraphies able to find an echo in the Seine Valley. Moreover, a focus on the capability of geophysical methods to detect archaeological features is also presented and discussed. We observed that the major reflections in the GPR records were generated from interfaces that have grain size variation: (1) the bottom of the Holocene colluvium and (2) the bottom of the upper Late Glacial silt. EMI and ERT show a very clear horizon associated with the upper Late Glacial silt, in some places even more clearly defined than with GPR. We confirmed the presence of a channel along the slope, placed under Locus 1, and a second channel of the same type globally following the paleotopography of Locus 2. We created a thickness map of the “beige sandy silt” and hypothesized a high probability of good preservation conditions of Magdalenian evidence. Finally, the detection of several localized diffraction hyperbolas in the GPR record offers the possibility to obtain the ground truth of the geophysical results in the near future and verify the nature (archaeological or geological) of these features.

The Lieth Moor area, located in the district of Pinneberg, Schleswig-Holstein, Germany, is a hotspot of Late Palaeolithic settlement activity. The exceptional abundance of archaeological sites is commonly attributed to the presence of a large palaeolake. However, in the Weichselian Late glacial, there were numerous large lakes in Schleswig-Holstein. Thus, a well-founded explanation for the find concentration in Lieth Moor is still lacking, and forming a research desideratum until today. To improve our understanding of this Late Pleistocene landscape and its appeal to hunter–gatherer groups of that time, we conducted a large-scale archaeogeophysical study focusing on a possible ford of the potential palaeolake. We employed Ground-Penetrating Radar and Electromagnetic Induction measurements, supplemented by existing legacy drill-probing data, to identify and map limnic gyttja (organic lake mud) sediments and their spatial distribution within the area. The findings of our study indicate that during the Late Pleistocene to Early Holocene, the Lieth Moor area comprised a cluster of small ponds rather than a continuous lake. These ponds likely interconnected during periods of increased water levels. The presence of dry islands within the region corresponds with archaeological evidence, suggesting that Late Palaeolithic communities visited some of these islands. The absence of the previously postulated palaeolake places the known findings within a completely new palaeoenvironmental context: instead of the previously suspected ford, we assume that the proximity to the Elbe Palaeovalley played a decisive role in the repeated habitation of Lieth Moor. This area, rich in fresh water and fish, along with the dune chain situated to the west, serving as both a vantage point and windbreak, presented an ideal location for awaiting animals migrating along the river Elbe and/or as a resting place within the settlement system of mobile hunter–fisher–gatherer groups.

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.

We performed geophysical and geoarchaeological investigations in the Wadden Sea off North Frisia (Schleswig-Holstein, Germany) to map the remains and to determine the state of preservation of the medieval settlement of Rungholt, especially its southern dyke segment, called the Niedam dyke. Based on archaeological finds and historical maps, Rungholt is assumed to be located in the wadden sea area around the island Hallig Südfall. During medieval and early modern times, extreme storm events caused major land losses, turning cultivated marshland into tidal flats. Especially the 1st Grote Mandrenke (or St. Marcellus' flood), an extreme storm surge event in 1362 AD, is addressed as the major event that flooded and destroyed most of the Rungholt cultural landscape. Cultural traces like remains of dykes, drainage ditches, tidal gates, dwelling mounds or even plough marks were randomly surveyed and mapped in the tidal flats by several authors at the beginning of the 20th century. Due to the tidal flat dynamics with frequently shifting tidal creeks and sand bars, the distribution of cultural remains visible at the surface is rapidly changing, making it hard to create a comprehensive map of the cultural landscape by surveying. Today, the Niedam dyke area is fully covered by tidal flat sediments, depriving any remains from further archaeological investigation. Since little is known about the precise location or state of preservation of these remains, our investigation aimed at the rediscovery of the medieval dyke system and associated structure with modern and accurate geophysical, geodetical and geoarchaeological methods. Magnetic gradiometry revealed a large part of the medieval dyke, confirming two tidal gates and several terps connected inland with the dyke, providing a detailed example of a Frisian medieval dyke system. Based on our results, the so far inaccurate and incomplete maps of this part of Rungholt can now be specified and completed. Beyond that, seismic reflection profiles give a first depth resolving insight in the remains of the dyke system, revealing a severe threat to the medieval remains by erosion. The site is exemplary for the entire North Frisian coast, that was influenced by multiple flood events in the middle ages to modern times.

In this study, we investigate GPR measurements in freshwater of less than 5 m at four different locations to derive rules of thumb in terms of depth penetration, resolution, and material contrasts of the method for 200 and 400 MHz antennas under field conditions. The objective is to improve the attractiveness of the method for archaeological issues in water, as there are hardly any studies on this subject so far. The depth penetration of 2–4 m is negligibly influenced by the choice of the 200 or 400 MHz antenna. Organic material in the water column also does not affect the water depth but offers new fields of applications for mapping and volume estimation of biomass in lakes with GPR. The horizontal resolution in the cm range in the direction of the profile and in the dm range across the profile could not be improved by the narrow antenna radiation pattern of 0.1 for materials of 45% porosity to water-saturated sand. Large reflection coefficients also show, for example, granite to sand and gyttja to wood. The water column has a considerable effect on the data quality of the 400 MHz antenna from a depth of 2 m due to the antenna ringing. Furthermore, multiples must be expected in a water column

Hyperbolic diffractions in Ground Penetrating Radar (GPR) data are caused by a variety of subsurface objects such as pipes, stones, or archaeological artifacts. Supplementary to their location, the propagation velocity of electromagnetic waves in the subsurface can be derived. In recent years, it was shown that deep learning tools can automatically detect hyperbola in radargrams using data measured over urban infrastructure, which are relatively clear. In contrast, in this study, we used an archaeological dataset with diverse underground structures. In the first step we used the deep learning network RetinaNet to detect hyperbola automatically and achieved an average precision of 0.58. In the next step, 10 different approaches for hyperbola fitting and thus velocity determination were applied. The derived information was validated with manually determined velocities and apex points. It was shown that hyperbola extraction by using a threshold and a column connection clustering (C3) algorithm followed by simple hyperbola fitting is the best method, which had a mean velocity error of 0.021 m/ns compared to manual determination. The average 1D velocity-depth distribution derived in 10 ns intervals was in shape comparable to the manually determined one, but had a systematic shift of about 0.01 m/ns towards higher velocities.

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 m2/yr to 0.7 m2/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.

Magnetic mapping is a common method for investigating archaeological sites. Typically, the magnetic field data are treated with basic signal improving processing followed by image interpretation to derive the location and outline of archaeological objects. However, the magnetic maps can yield more information; we present a two-step automatized interpretation scheme that enables us to infer the social structure from the magnetic map. First, we derive the magnetization distribution via inverse filtering by assuming a constant depth range for the building remains. Second, we quantify the building remains in terms of their total magnetic moment. In our field example, we consider this quantity as a proxy of household prosperity and its distribution as a social indicator. The inverse filtering approach is tested on synthetic data and cross-checked with a least-squares inversion. An extensive modeling study highlights the influence of depth and thickness of the layer for the filter construction. We deduce the rule of thumb that by choosing a rather too deep and too thick a layer, errors are smaller than for layers too low and too thin. The interpretation scheme is applied to the magnetic gradiometry map of the Chalcolithic Cucuteni–Tripolye site Maidanetske (Ukraine) that comprises the anomalies of about 2300 burned clay buildings. The buildings are arranged along concentric ellipses around an inner vacant space and a vacant ring. Buildings along this ring corridor have increased total magnetic moments. The total magnetic moment indicates the remaining building material and therefore the architecture. Lastly, architecture can reflect economic or social status. Consequently, the increased magnetic moment of buildings along the ring corridor indicates a higher economic or social status. The example of Maidanetske provides convincing evidence that the inversion of magnetic data and the quantification of buildings in terms of their magnetic moments enables the investigation of the social structure within sites.

This paper shows that local differences in house orientation in settlements from the Early Neolithic in Central Europe reflect a regular chronological trajectory based on Bayesian calibration of 14C-series. This can be used to extrapolate the dating of large-scale settlement plans derived from, among other methods, geophysical surveys. In the southwest Slovakian settlement of Vráble, we observed a progressive counter-clockwise rotation in house orientation from roughly 32° to 4° over a 300 year period. A survey of published and dated village plans from other LBK regions confirms that this counter-clockwise rotation per settlement is a wider Central European trend. We explain this observation as an unintentional, unconscious but systematic leftward deviation in the house builders' cardinal orientation, which has been termed "pseudoneglect" in studies of human perception. This means that whenever houses were intended to be oriented towards a specific direction and be parallel to each other, there was an error in perception causing slight counter-clockwise rotation. This observation is used as a basis to reconstruct dynamics of Early Neolithic settlement in the Slovakian Žitava valley, showing a rapid colonization, followed by increased agglomeration into large villages consisting of strongly autonomous farmsteads.

Water-layer multiples pose a major problem in shallow water seismic investigations as they interfere with primaries reflected from layer boundaries or archaeology buried only a few meters below the water bottom. In the present study we evaluate two model-driven approaches (“Prediction and Subtraction” and “RTM-Deco”) to attenuate water-layer multiple reflections in very shallow water using synthetic and field data. The tests comprise both multi- and constant-offset data. We compare the multiple removal efficiency of the evaluated methods with two traditional methods (Predictive Deconvolution and SRME). Both model-driven approaches yield satisfactory results concerning the enhancement of primary energy and the attenuation of multiple energy. For the synthetic test cases, the multiple energy is reduced by at least 80% for the Prediction and Subtraction approach, and by more than 60% for the RTM-Deco approach. The application to two field data sets shows a significant amplification of primaries formerly hidden by the first water-layer multiple, with a reduction of multiple energy of up to 50%. The waveforms obtained from FD modeling match the true waveforms of the field data well and small deviations in time and amplitude can be removed by a time shift of the traces as well as an amplitude adaption to the field data. The field data examples should be emphasized, where the tested Prediction and Subtraction approach works significantly better than the traditional methods: the multiples are effectively predicted and attenuated while primary signals are highlighted. In conclusion, this shows that this method is particularly suitable in shallow water applications. Both evaluated multiple attenuation approaches could be successfully transferred to two other 3D systems used in shallow water near surface investigations. Especially the Prediction and Subtraction approach is able to enhance the primaries for both tested 3D systems with the multiple energy being reduced by more than 50%.

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 applicability of offshore geoelectrical profiling in the littoral zone, e.g., for archaeological prospection, sediment classification and investigations on coastal ground water upwelling. We performed field measurements with a 20 m long multi-electrode streamer in inverse Schlumberger configuration, which we used to statistically evaluate measurement uncertainty and the reproducibility of offshore electric resistivity tomography. We compared floating and submerged electrodes, as well as stationary and towed measurements. We found out that apparent resistivity values can be determined with an accuracy of 1% to 5% (1σ) depending on the measurement setup under field conditions. Based on these values and focusing on typical meter-scale targets, we used synthetic resistivity models to theoretically investigate the tomographic resolution and depth penetration achievable near-beach underneath a column of brackish water of about 1 m depth. From the analysis, we conclude that offshore geoelectric sounding allows the mapping of archaeological stone settings. The material differentiation of low-porosity rock masses < 15% is critical. Submerged wooden objects show a significant resistivity contrast to sand and rocks. Distinguishing brine-saturated sandy sediments from cohesive silty-clayey sediments is difficult due to their equal or reversed resistivity contrasts. Submarine freshwater discharges in sandy aquifers can be localized well, though difficulties may occur if the seafloor encounters massive low-porosity rock masses. As to the measurement setups, submerged and floating electrodes differ in their spatial resolution. Whereas stone settings of 0.5 to 1 m can still be located with submerged electrodes within the uppermost 4 m underneath the seafloor, they have to be >2 m if floating electrodes are used. Therefore, we recommend using submerged electrodes, especially in archaeological prospection. Littoral geological and hydrogeological mapping is also feasible with floating electrodes in a more time-saving way.

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 investigate how suitable ground penetrating radar (GPR) and geoelectrics are to prospect the remains of submerged wooden archaeological constructions in the water column. For this purpose, we determined the contrasts of electric resistivity and dielectric permittivity from measurements on present-day wood samples, serving as simplified approximations of water saturated and undegraded archaeological wood. As common substitutes of hard and soft construction wood, we investigated oak and spruce wood. The electric resistivity and dielectric permittivity were determined with increasing moisture content from small-scale electric and GPR measurements using a Wenner alpha array and a 2 GHz Palm antenna in a watering experiment. In a water-saturated state, resistivity values of 20 were observed. The anisotropy effects and deviations of the wood species were seen to be up to 30%. On the basis of this, the relative material contrasts of wood with respect to fresh water, sand, and clay were calculatedand compared to values found in the literature for seismic contrasts. Geoelectric, GPR, and seismic measurements show contrasts of 0.3 to 0.8, −0.4 to 0.2, and −0.24 to 0.35, depending on the surrounding material and structural orientation of the wood. The highest contrasts were found for wood in fresh water, followed by clayey and sandy subsoils. On the basis of the determined contrasts, analytical calculations were performed showing that an object of 0.5 m diameter can be detected at depths between 0.5 m and 1.5 m with geoelectrics (Schlumberger) and at depths between 0.5 m and 3 m with ground penetrating radar measurements (400 MHz).

Throughout mankind’s history, the need to secure and protect the home settlement was an essential one. This holds especially true for the city of Ainos (modern Enez) in Turkish Thrace. Due to its continuous settlement history since the 7th/6th century BC, several different types of city walls were built—sometimes even on top of each other—several of which have been preserved over time. To decipher the construction style, the course and the age of a buried city wall segment in the southern part of the former city, a geoscientific multi-proxy approach including magnetic gradiometer (MG) and electrical resistivity tomography (ERT) measurements in combination with granulometrical, sedimentological and microfaunistical investigations on sediment cores was applied. We were able to (1) present reasonable arguments for its Hellenistic age; (2) reveal the course of this wall segment and extrapolate it further north into a less studied area; and (3) demonstrate that in this near-coastal area, the former swampy terrain had been consolidated for constructing the wall. Our multi-proxy approach serves as a valuable example for investigating buried structures in archaeological contexts, avoiding a less-economical, time-consuming, or even forbidden excavation.