Use of ground penetrating radar to map subsurface archaeological features in an urban area

The ground penetrating radar (GPR) technique was used to investigate the subsurface in an urban area located in Mesagne (Italy) to obtain a map of the archaeological features in the ground. The GPR survey was undertaken at selected locations placed near (about 50 m) to a necropolis dating from the Messapian to the Roman imperial age, using a GSSI Sir System 2 incorporating 200 and 500 MHz centre frequency antennae. The selected areas (A and B) were surveyed along parallel 1 m spaced profiles using a 200 MHz antenna in area A and along parallel 0.5 m spaced profiles using a 500 MHz antenna in area B. For the selected areas the processed data were visualised in 3D space not only by means of the standard time slice technique, but also by means of a recently proposed approach, namely by iso-amplitude surfaces of the complex trace amplitude. The immediacy in revealing the spatial positioning of highly reflecting bodies, such as the anomaly interpreted as an old hypogeum room in area A, makes 3D visualisation techniques very attractive in archaeological applications of GPR. Their sensitivity to the signal/noise ratio is, on the other hand, highlighted by the quite poor performance in area B, where the only reliable result provided by all the techniques was the ancient living surface reflection, whereas none of them could effectively enhance the visibility of weak hyperbola reflections noted on 2D sections and probably related to the walls located on the ancient living surface. The performance of the various techniques in these two different situations allowed insights into their main advantages and drawbacks to be gained.     

Integrated seismic tomography and ground-penetrating radar (GPR) for the high-resolution study of burial mounds (tumuli)

The study primarily aims at providing adequate imaging resolution of large and prominent targets of archaeological interest, such as pyramids and tumuli, at all depth levels. We implemented an integrated seismic tomography and georadar (STG) technique to perform high-resolution imaging and characterization of tumuli (burial mounds). We tested the proposed technique on a preserved late Bronze Age burial mound in northern Italy, for which STG succeeded in performing an accurate 3-D reconstruction of the structure and stratigraphy as proved by later archaeological excavations. We completed two transmission seismic tomography measurements, at present ground level and at 1.5 m elevation, with a 24-channel seismograph and 15° angular separation between geophones. The ground-penetrating radar (GPR) dataset encompasses 12 250 MHz radial profiles and 12 common mid point gathers for velocity analysis. Shallow layers of the mound are successfully imaged by GPR, whilst the structure of the deep central part of the tomb is reconstructed from seismic traveltime inversion. In particular, GPR images lenses and layers of sediments forming the external part of the tumulus, evidences of a looting attempt, peripheral structures associated with later exploitation of the mound (furnaces) and, in the external sector of the tumulus, the top of the deep layer of silty sediments covering the funeral chamber. Tomographic results reveal seismic velocity anomalies of potential archaeological interest at ground level, which were successively validated by archaeological excavations. The integration of GPR and tomographic datasets is an effective strategy to overcome the imaging and interpretation problems related to the structure of such peculiar funeral monuments. STG can be applied to a virtually unlimited dimensional range and requires a limited data acquisition, processing and inversion effort. The results of the study allowed the identification of the funeral chamber and a detailed imaging of layering and structural details.     

GPR and inductive electromagnetic surveys applied in three coastal sambaqui (shell mounds) archaeological sites in Santa Catarina state, South Brazil

In this article, GPR – Ground Penetrating Radar and inductive electromagnetic (EM38 equipment) results are presented for three archaeological sites of shell mounds located between Laguna and Jaguaruna cities, mid-southern coast of Santa Catarina state, Brazil. The Jabuticabeira II, Santa Marta IV, and Encantada III sites were studied. These sites consist of carbonate shell mounds built by societies that inhabited this region during the pre-colonial period (7.5–1.3 ky BP). Geophysical responses obtained in these archaeological sites had allowed the identification of anomalous regions related to archaeological point targets as well as continuous targets, correlated to the shell mound's structure and buried geologic features. GPR numerical modelling studies were done in order to evaluate the performance of the GPR method in relation to materials found in that environment, besides helping interpretations of real results. The synthetic GPR model was a good match with real data obtained in the field. So, the real and synthetic results serve as a reference guide to direct excavation activities in these archaeological and geological structures. In the Jabuticabeira II shell mound, lithic materials and burials were found amidst archaeological layers characterized by sediment rich in organic material, ferruginous concretions, and loose shells. In the Santa Marta IV archaeological site, a concentration of ceramic materials, burials and fireplaces were found. In the Encantada III site, a strong anomaly characterized by a hyperbolic reflection was related to the presence of a silicified tree root; a great example of ambiguity in shallow geophysics interpretation in archaeological applications. Results functioned as a guide to archaeological excavations, allowing cost and time reduction as well as contributing to improve the knowledge of these archaeological sites situated in southern Brazil.     

Integrated geophysical surveys at the archaeological site of Tindari (Sicily,Italy)

The results of geophysical survey carried out at the archaeological site of Tindari, located 70 km west-north-west of Messina (Sicily, Italy), are presented and discussed. The site is one of the most important archaeological sites, about 90 ha large, on the north-eastern side of Sicily and it was one of the last Greek settlements in Sicily. Keeping in mind the vastness of the site, according to the archaeologists an area, of about 1200 m² was selected, for a geophysical survey. Electrical resistivity tomography (ERT), including induced polarization (IP) measurements, and seismic refraction tomography were applied in the course of geophysical exploration in the area. The objective of this preliminary geophysical investigation was to verify the effectiveness of the technique for a spatial definition of the buried archaeological structures (mostly walls, columns, etc.) to determine their characteristics and to study the presence of collapsed columns within the zone of archaeological interest, in view of their application in the rest of the site, that, for obvious reasons, could not be entirely excavated in brief times. The geophysical data, visualized in 3D space, revealed a distribution of low-contrast shallow anomalies that indicate the presence of different types of buried structures in the surveyed area. Also, the probable accumulations of collapsed columns could be determined. These results confirm the hypothesis of a large northern extension of the archaeological site and provide useful information to design a more efficient excavation plan.     

Three-dimensional image of seismic refraction tomography and electrical resistivity tomography survey in the castle of Occhiolà (Sicily,Italy)

A geophysical survey was carried out at the ruins of Occhiolà Castle (Sicily, Italy), a medieval village located at the north-western part of a hill named “Terravecchia” at 491 m asl. In order to map the archaeological structures (such as walls and burrows), Electrical resistivity tomography (ERT) and seismic-refraction tomography methods were used. The resistivity variations are known to correlate quite well with the lithology, thus providing important information for identifying the buried archaeological remains. On the other hand, seismic velocity variations provide information about the geometrical features of the remains. The two geophysical methods are used in an enjoined way to better aid the interpretation and evaluate the significance and reliability of the results obtained with each single method. The electrical and seismic data are displayed in three dimensions by using the iso-resistivity and iso-velocity surfaces, respectively. This allows us to effectively define the location of the anomalies present in the single electrical and seismic sections. The results obtained in the survey highlight the presence of structures of regular shape, probably due to features of archaeological interest.     

Cesium magnetometer surveys at a Pithouse Site near Silver City,New Mexico

The La Gila Encantada Site is located on an isolated ridge top north of Silver City, New Mexico in an area defined culturally as the Mimbres Mogollon region. The 180 m × 80 m (14400 m²) sized site was originally recorded as a dense scatter of ceramics, lithics, and ground stone along with a number of surface depressions that appear to represent pithouses. Cesium magnetometer surveys were conducted to identify hearths, pithouse boundaries, and activity areas outside of pithouses in support of archaeological investigations, and to test this instrument's ability to image these features. This paper presents a characterization of the magnetic signature of a pithouse as a magnetic high in the center caused by the central hearth, low magnetic variability along the floor of the house, and then increased magnetic variability at the pithouse boundary. This characterization was successfully confirmed for three pithouses using cesium magnetometery and archaeological excavation.     

New radiometric ages for the Early Upper Palaeolithic type locality of Brno-Bohunice (Czech Republic): comparison of OSL,IRSL, TL and C dating results

New radiometric data are reported from the recent excavation of the type locality of the Early Upper Palaeolithic entity of the Bohunician. Recently obtained radiocarbon (¹⁴C) data on charcoal are compared with new Optically Stimulated Luminescence (OSL) dating of sediment. OSL ages were determined on sediment from the archaeological occupation at Brno-Bohunice, as well as from the over- and underlying loessic sediments. Multiple techniques were applied, which all gave congruent results. While a dual protocol (post IR-OSL) failed the quality criteria tests, ages were obtained by Multiple-Aliquot-Additive-Dose (MAAD) on polymineral material and Single-Aliquot-Regeneration (SAR) on fine grain quartz extract as well as on polymineral material. Fading tests show significant loss of Infrared Stimulated Luminescence (IRSL) after storage for 3 and 12 months for one sample, but little or no fading for others. The resulting (uncorrected) age estimates are smaller than those on quartz by OSL methods. The latter are considered to be more reliable estimates of the sedimentation age of these deposits. The measured OSL doses do not show a simple distribution and the lowest 5% was used for age calculation to represent the most likely sedimentation age. The quartz from the loess overlying the archaeological layer is OSL dated to 30.9 ± 3.1 ka, while the sediment for the paleosol which contains the archaeological layer gave an age of 58.7 ± 5.8 ka. The attribution of this paleosol to the Hengelo interstadial is therefore questionable. However, if the Hengelo interstadial is correlated with the Dansgaard/Oeschger (D/O) event 12, statistical agreement within 2-σ is achieved. The OSL result for the archaeological layer is in accordance with a weighted average TL date on heated flint artifacts of 48.2 ± 1.9 ka from this layer as well as calibrated radiocarbon data (CalPal Hulu 2007) from nearby locations. However, radiocarbon data on charcoal samples obtained during excavation at Brno-Bohunice 2002 provide age estimates between 30 and 40 ka ¹⁴C-years, which translate to approximately (33) 35–44 ka on the calendric time scale according to the Hulu 2007 model. For the underlying loess a depositional age of 104.3 ± 10.6 ka was obtained by OSL. The presented OSL ages indicate that a simple correlation of soil sequences between sites within a region has to be verified by chronometric dating.     

Magnetic imaging and electrical resistivity tomography studies in a Roman military installation found in Satala archaeological site,northeastern Anatolia,Turkey

Satala is one of the last great military centers in the Roman East available for archaeological and historical investigations. This archaeological site is situated on the crossing of two singularly important routes in North-East Asia Minor. Only little archaeological fieldwork in and around Satala has so far been carried out, thus delimiting our knowledge of the site and its role within the Roman military structure of the East. In order to provide further data on the site, we carried out geophysical surveys including the application of magnetic and electrical resistivity techniques in an area north of the Sadak village. The geophysical surveys within this area (supposedly part of the Roman military camp) proceeded in two stages. The first stage saw magnetic gradiometer imaging studies being carried out on 2.1 ha, measured by a Geoscan FM-36 fluxgate gradiometer by using 0.5 × 1 m grid intervals. For the second stage a different geophysical technique was used—resistivity tomography. Resistivity data were collected using a number of combined 2D resistivity pseudosections in the eastern part of the area which contained very regular magnetic anomalies. The 3D data were obtained by the combination of all survey lines collected from 2D data sets, and thereafter the arranged data were processed by using the 3D robust inversion modified from the smoothness-constrained algorithm. Electrical resistivity tomography investigations revealed that the buried archaeological structures might be located near the surface, except for some structures found in the middle of the studied area. The archaeological structures were furthermore determined by realistic model sections and volumetric representations. Magnetic imaging and electrical resistivity tomography surveys show that the combined usage of these techniques advances the understanding of archaeological structures beneath the surface.     

Magnetic detection of ship ballast deposits and anchorage sites in King Herod's Roman harbour,Caesarea Maritima,Israel

Ballast stone deposits are a common feature of sediments in ancient harbour basins but are often overlooked as a potential source of archaeological information. Recent geophysical investigations at Caesarea Maritima in Israel have discovered a thick, laterally extensive ballast layer in the area seaward of the 1st c. BC Roman harbour. The ballast deposits were identified by low-relief mounds on the seabed with elevated magnetic intensities. Jet probing and excavation of magnetic anomalies at several locations revealed a 20–60 cm thick rubble layer containing large quantities of Late Roman and Byzantine pottery, local sedimentary boulders (kurkar sandstone, limestone cobbles) and foreign igneous and metamorphic boulders (granite, schist, volcanics; ca. 50%). The foreign boulders and pottery identify the rubble layer as ballast and ships refuse jettisoned by merchant ships outside the harbour. The strong magnetic contrast between the ballast deposits and the natural seabed sediments is attributed to the high magnetic susceptibility (>10⁻³ SI) of crystalline boulders and pottery materials within the ballast rubble.     

Can palynomorphs occur in burnt ancient potsherds? An experimental proof

Potsherds from archaeological sites are usually found barren of palynomorphs. Recently half burnt potsherds from an excavation site dated to 3320 ± 400 years BP to 2080 ± 80 years BP in Gangetic West Bengal, India yielded some palynotaxa. That palynomorphs can remain preserved in half burnt potsherds was experimentally proved with pots made from mud mixed with known pollen grains and burning for variable duration. Pots burnt for shorter duration showed recovery, while prolonged firing showed poor/no recovery. Present study suggests the use of pollen grains from under-fired ancient potsherds as a reliable tool in deciphering vegetation, climate, agriculture and related parameters of archaeological research.     

Joint geophysical measurements to investigate the Rossano of Vaglio archaeological site (Basilicata Region,Southern Italy)

A high resolution geophysical survey was carried out in the archaeological site of Rossano di Vaglio (Basilicata Region, Southern Italy), where an important ancient sanctuary is located. It was built during the IV century B.C. and devoted to the goddess Mephitis. The sanctuary rises in an area affected by a multiple and retrogressive rototraslational landslide, historically and presently subject to reactivation. The main objective of this work was the identification of buried structures of archaeological interest in an area designated by the Archaeological Superintendence of the Basilicata Region. The study was performed by means of the use of high resolution geophysical surveys. In particular, we made use of the joint application of three highly sensitive and non-invasive geophysical techniques, namely the Geoelectrical, the Magnetic and the Ground Probing Radar (GPR) methodologies. In such a way, we obtained two important results: first, we provided the archaeologists with information about the limits of the areas to be excavated; second, we could verify in real time the reliability of the geophysical results. The experimental results showed four main magnetic anomalies in the area of study, in agreement with the GPR results obtained for the same target. Finally, a partial excavation test of the investigated area revealed a buried building structure, located in correspondence of an anomaly identified by means of the geophysical prospecting.     

Structured light scanning for high-resolution documentation of in situ archaeological finds

Archaeologists strive to document the process of excavation and discovery as completely as possible. Over the past several decades archaeologists have incorporated a growing number of computerized techniques for documenting archaeological finds. Scanning is one such technique. There are a number of technologies that now allow archaeologists to scan structures, excavation surfaces and in situ artifacts to create high-resolution, 3D data sets. We report here on a trial application of one of these, a structured-light scanner, to create 3D representations of excavated surfaces and associated artifacts at two Middle Paleolithic sites in southwest France. In each instance, surfaces of approximately 2.5 m² were scanned in approximately 1 day. The resulting data sets are very good representations of the originals in terms of colors and spatial details, and as such provided an important piece of archaeological documentation. To use this equipment successfully in the field, however, required solving a number of logistical issues, and the amount of time required to learn to use this equipment was significant. Once these issues are addressed, this technology is appropriate for documenting extraordinary, unique finds where time and costs are offset by the importance of good documentation.     

Ground Penetrating Radar in the Northern Great Lakes: A Trial Survey of a Contact Period Occupation in Marquette County,Michigan

ABSTRACT

Timing, inclement weather, and limited budgets can obstruct archaeological survey efforts. Here, we ameliorate some of these with use of ground penetrating radar (GPR) at the Goose Lake Outlet #3 (GLO#3) site (20MQ140). GPR surveys to guide survey and excavation efforts in these well-drained sandy soils are limited. GPR imagery exhibited false positives; however, shovel probes, subsequent site excavation, and artifact analysis led to numerous discoveries at this protohistoric site—including glass trade beads dating to the 1630s. These discoveries solidify evidence for some of the earliest European-made trade items in the region and provide further confirmation for placement of an Indian trail and ancient travel corridor between Lake Superior and Lake Michigan. Two ¹⁴C samples obtained from moose-bone in features with trade goods fall within the expected occupation range of 1633–1668 and 1615–1647. Several anomolies, possible hearth clusters were seen on the GPR imagery; however, many features remained invisible.     

3D geometrical modeling of excavation units at the archaeological site of Tell ‘Acharneh (Syria)

3D geometric modeling consists of representing geometric and spatial relationships of volumetric objects. We think it could be helpful in the context of archaeological excavation units representation and analysis. This article presents a procedure developed to generate 3D models from GPS positions taken at the top and the bottom of the excavation units boundaries on the archaeological site of Tell ‘Acharneh (Syria). It shows and discusses two geometrical modeling approaches (voxel and tetrahedral) used in the Gocad 3D modeling tool. Once excavation units are geometrically modeled, it is possible to refer them within a trench or the entire archaeological site, to handle them in various ways (zoom, rotation, translation), to perform on them 3D spatial analysis such as volumetric calculus or intersection computation, to make various kinds of queries such as to find out excavation units that have a certain number of artefacts, to generate sections anywhere in the 3D model, and finally to publish it with VRML (Virtual Reality Modular Language). As well as improving data analysis techniques, we think that if this 3D modeling operation can be done during the excavation, it could greatly help archaeologists to plan more efficiently their daily excavation strategy.     

Marine electrical resistivity tomography for shipwreck detection in very shallow water: a case study from Agropoli (Salerno,southern Italy)

A marine electric resistivity survey was carried out over a submerged beach along the Agropoli shore (Salerno, Italy) to detect buried objects of archaeological interest below the sandy seabed. We found a shipwreck, a military vessel that probably sunk during the Salerno landing operations of the allied forces in the Second World War. Resistivity data provide information on the vertical and horizontal extension of the shipwreck, which is characterized by very low calculated resistivity values (about 2–5 ohm m). Such values differ significantly from the sand and the bedrock values (5–40 ohm m). Although the presence of the shipwreck is clearly visible from geoelectric data, the joint application of electric, magnetic and multibeam bathymetric techniques reduces the ambiguities inherent in each method. As shown in the Electrical Resistivity Tomography (ERT) and confirmed by the Digital Elevation Model (obtained from the processing of bathymetric data), the shipwreck extends more than 30 m in NE–SW direction and it is about 13 m wide. The global extension of the relic is consistent with the magnetic data, that are characterized by a magnetic anomaly with an amplitude of about 1800 nT and similar dimension, as inferred from the estimation of source boundaries obtained from the computation of the analytic signal. The results of our survey encourage the use of marine geoelectrical methods for the detection of buried archaeological targets, particularly in locations where the use of seismic prospecting is not effective (e.g. very shallow water with sandy sea-bottoms). The integration of different geophysical methods allows to better define the extension, depth and thickness of buried objects, suggesting that such an approach is the most effective for underwater archaeological investigations.     

Quantifying tomb geometries in resistivity images using watershed algorithms

Quantifying the geometries (defined here as width, height and depth of burial) of archeological structures within resistivity models produced as a result of the regularization constraints used in most inversion algorithms is difficult, especially when structures are closely spaced. Here we apply the watershed by simulated immersion method of boundary detection to smooth 2D resistivity images generated for synthetic and field data over 3D targets. The synthetic studies include a single cavity model, a model for two widely spaced cavities (spacing ? unit electrode spacing) and a model for two closely spaced cavities (spacing < unit electrode spacing). We also examine a single-cavity model where a relatively resistive overburden, common at archaeological sites in Egypt, is included. In the case of the single cavity models, the maximum error for any geometries are 18% for the model without the resistive overburden and 10% for the model where the overburden is included, whereas it increases to 24% for the widely spaced model and 40% for the closely spaced model. Despite, the higher errors in the closely spaced cavity model, application of the algorithm confirms the presence of two features, which is not ascertainable from the smooth resistivity images. Boundaries detected with the watershed algorithm are subsequently used to define a disconnect in the regularization, resulting in a markedly improved estimate of the resistivity structure (particularly for the closely spaced cavity model) in a second inversion step using the model obtained from the smoothness constraint inversion as the starting model. This revised resistivity model also results in a lower root-mean-square (rms) misfit between measured and theoretical data, and between synthetic and inverted models. We demonstrate how the method can be applied on images from the archaeological site at Qurnet Murai, Luxor City Egypt.     

Non-invasive archaeological exploration in stratigraphically complex rural settings: an example from Ferento (Viterbo, Italy)

This paper presents the results of a ground penetrating radar (GPR) survey conducted in the rural area of the ancient city of Ferento near Viterbo (Northern Latium, Italy), where the visible ruins testify an uninterrupted history from the Etruscan to the Medieval period. The soil covering the ruins has been intensively ploughed, cultivated and used for pasture over the centuries, and is characterised by a high content of clay minerals. In this scenario, the use of GPR is particularly difficult due to both possible ‘false positive anomalies’ and strong signal attenuation. In order to reconstruct the geometrical shape and size of the geophysical anomalies related to a medieval dwelling, multi-profile GPR data were collected in two different areas, using radar equipped with 250-MHz antennas. 2D and pseudo-3D migrated maps were created and signal velocity analysis was performed to estimate the target depths. Despite the clayey-calcareous material, the results show a good signal penetration and allow a clear reconstruction of the shape of the buried ruins, confirming that the anomalies are due to medieval buildings located at a depth in agreement with what was observed on the archaeological stratigraphy. Furthermore, our study suggests that in radar data interpretation, caution should be used if the surveyed area was subjected to prolonged and intensive agriculture.     

Methods for the extraction of archaeological features from very high-resolution Ikonos-2 remote sensing imagery,Hisar (southwest Turkey)

Archaeological research in the territory of Sagalassos is a multidisciplinary project covering scientific disciplines traditionally linked to archaeology and also new technologies such as very high-resolution remote sensing with sufficient radiometric and spatial resolution (<2.5 m). This paper focuses on the evaluation of GIS-, pixel- and object-based techniques for automatic extraction of archaeological features from Ikonos-2 satellite imagery, which are then compared to a visual interpretation of ancient structures. The study was carried out on the unexcavated archaeological site of Hisar (southwest Turkey). Although all techniques are able to detect archaeological structures from Ikonos-2 imagery, none of them succeed in extracting features in a unique spectral class. Various landscape elements, including archaeological remains, can be automatically classified when their spectral characteristics are different. However, major difficulties arise when extracting and classifying archaeological features such as wall remnants, which are composed of the same material as the surrounding substrate. Additionally, archaeological structures do not have unique shape or colour characteristics, which can make the extraction more straightforward. In contrast to automatic extraction methods, a simple visual interpretation performs rather well. The methods presented in this paper can be applied with variable success to archaeological structures composed of the same material as the surrounding substrate, which is often the case.     

Reduction of topography effect in inductive electromagnetic profiles: application on coastal sambaqui (shell mound) archaeological site in Santa Catarina state, Brazil

We present results of electrical conductivity profiles obtained with inductive electromagnetic geophysical method in Santa Marta archaeological site, Santa Catarina State, southern Brazil. This site is a sambaqui (shell mound), in which several human occupation remains are found during pre-colonial period such as buried lithic and bone artifacts, fire-place, etc. Most of these mounds include as well many human burials, which, in many cases, point to funerary ritual as a main agency for mound building. A set of profiles of apparent electrical conductivity and magnetic susceptibility was acquired in two sites aiming the identification of geophysical anomalies with potential interest for excavation. To enhance conductivity data, we applied an effective procedure to remove topographical effects in the apparent conductivity measurements, which are rather conditioned by the presence of a variable water table depth or conductive sediment layer. A linear dependence among conductivity values and the site elevation provided a simple linear model to remove the influence of topography. Corrected electric conductivity maps substantially improved the definition of anomalies, many of them rather subtle in raw data images. Corrected maps also show a better adherence with magnetic susceptibility maps, both of them identifying archaeological structures of interest: a well-structured fire-place and a concentration of ceramic fragments.     

Identification of fragmentary bone from the Pacific

We explore bone microstructure for taxonomic identification of archaeological bones too fragmentary to permit secure identification on morphological grounds. Backscattered electron (BSE) imaging is used to observe bone tissue types and the arrangement of vascular canals, and to facilitate quantification of osteonal canal dimensions. Examination of known examples of relevant taxa (humans [n = 8], pigs [n = 4] and dogs [n = 4]) shows significant differences among them. When the results of this examination are applied to a blind test of modern and archaeological specimens (humans [n = 8], pigs [n = 2]), 100% of specimens are identified correctly. The approach is applied to 13 morphologically unidentifiable fragments from Hawai’i and Fiji to evaluate its potential for identifying bone tools and to increase the number of samples available for dietary analysis. Potential applications of the approach for other contexts are discussed.