Multiple Bayesian modelling approaches to a suite of radiocarbon dates from ovens excavated at Ysgol yr Hendre,Caernarfon, North Wales

In recent years, the use of Bayesian statistical analysis for modelling radiocarbon dates has become much more commonplace within archaeology, however there remains very little literature to aid new practitioners with understanding the archaeological assumptions that underpin many of the modelling choices. Using a suite of radiocarbon dates from a sample of seven ovens, from the 18 excavated across the site at Ysgol yr Hendre, Caernarfon, North Wales, this paper explores the various methods of modelling this group of features. It stresses the underlying archaeological assumptions of each model, while also highlighting some of the tools that should especially be considered when developing very precise models. Finally it compares the results of the models.     

Archaeoseismology: Methodological issues and procedure

Archaeoseismic research contributes important data on past earthquakes. A limitation of the usefulness of archaeoseismology is due to the lack of continuous discussion about the methodology. The methodological issues are particularly important because archaeoseismological investigations of past earthquakes make use of a large variety of methods. Typical in situ investigations include: (1) reconstruction of the local archaeological stratigraphy aimed at defining the correct position and chronology of a destruction layer, presumably related to an earthquake; (2) analysis of the deformations potentially due to seismic shaking or secondary earthquake effects, detectable on walls; (3) analysis of the depositional characteristics of the collapsed material; (4) investigations of the local geology and geomorphology to define possible natural cause(s) of the destruction; (5) investigations of the local factors affecting the ground motion amplifications; and (6) estimation of the dynamic excitation, which affected the site under investigation. Subsequently, a ‘territorial’ approach testing evidence of synchronous destruction in a certain region may delineate the extent of the area struck by the earthquake. The most reliable results of an archaeoseismological investigation are obtained by application of modern geoarchaeological practice (archaeological stratigraphy plus geological–geomorphological data), with the addition of a geophysical-engineering quantitative approach and (if available) historical information. This gives a basic dataset necessary to perform quantitative analyses which, in turn, corroborate the archaeoseismic hypothesis. Since archaeoseismological investigations can reveal the possible natural causes of destruction at a site, they contribute to the wider field of environmental archaeology, that seeks to define the history of the relationship between humans and the environment. Finally, through the improvement of the knowledge on the past seismicity, these studies can contribute to the regional estimation of seismic hazard.     

Geophysical Survey and Archaeological Data at Masseria Grasso (Benevento,Italy)

The use of geophysical methods in metrology is a significant tool within the wide research topic of landscape archaeology context. Since 2011, the Ancient Appia Landscapes Project aims to recognize dynamics, shapes and layout of the ancient settlement located along the Appia road east of Benevento, and cyclical elements and human activities that influenced the choice of landscapes. The integration of geophysical data with an archaeological infra-site analysis allowed us to investigate the area of Masseria Grasso, about 6 km from Benevento (Campania region, Italy). In this framework, an archaeogeophysical approach (Geomagnetic and Ground Penetrating Radar) was adopted for detecting anomalies potentially correlated with buried archaeological evidences. The geomagnetic results have given a wide knowledge of buried features in a large survey highlighting significant anomalies associated with the presence of buildings, roads and open spaces. These geophysical results permitted us to define the first archaeological excavations and, successively, a detailed Ground Penetrating Radar approach has been provided highlighting the rooms and paved spaces. The overlap between archaeological dataset and geophysical surveys has also allowed recognizing the path of the ancient Appia road near the city of Benevento and hypothesize the settlement organization of the investigated area, which has been identified with the ancient Nuceriola.     

ESR chronology of the Somme River Terrace system and first human settlements in Northern France

The Somme Valley, Northern France, is famous for its archaeological sequence, where numerous rich Palaeolithic sites, such as Saint-Acheul, the type site of the Acheulian, have been discovered. The archaeological levels are often directly associated with fossil alluvial sediments of the River Somme or with slope deposits, including loess and palaeosols. In the middle reaches of the valley, near Amiens, the system of fossil-stepped fluvial terraces is particularly well developed and preserved, and occurs on 10 alluvial formations. These terraces, from +5 to +55 m above the present-day valley bedrock, allow the study of the environmental changes and the human settlement of this area through the Pleistocene.

Since 1988, ESR dating was systematically applied on bleached quartz extracted from the fossil fluvial deposits, in order to better describe the geological evolution of the stepped system. More recently, U-series/ESR dating has also been performed on teeth collected from the different terrace deposits. Here we present a synthetic review of the main ESR results, and propose an ESR chronology for the geological evolution of the Somme fluvial system and for the Middle Pleistocene human settlement of northern France.     

Multi-frequency Electromagnetic Induction Survey for Archaeological Prospection: Approach and Results in Han Hangu Pass and Xishan Yang in China

This study presents the potential of multi-frequency electromagnetic induction (EMI) in archaeology. EMI is currently less employed for archaeological prospection with respect to other geophysical techniques. It is capable of identifying shallow subsurface relics by simultaneously measuring the apparent electrical conductivity (EC_a) and apparent magnetic susceptibility (MS_a). Moreover, frequency sounding is able to quantify the depths and vertical shapes of buried structures. In this study, EMI surveys with five frequencies were performed at two heritage sites with different geological conditions: Han Hangu Pass characterized by cinnamon soil and Xishan Yang by sandy loams. In the first site, high EC_a values were observed with variations in depth correlated to archaeological remains. Moreover, electromagnetic anomalies related to an ancient road and five kiln caves were identified. In the second site, an ancient tomb, indicating extremely low EC_a and high MS_a, was discovered. Its electromagnetic properties are attributed to the cavity and ferroferric oxides.     

High resolution 3D ERT to help GPR data interpretation for researching archaeological items in a geologically complex subsurface

Muro Leccese (Lecce) contains one the most important Messapian archaeological sites in southern Italy.The archaeological interest of the site arises from the discovery of the remains of Messapian walls, tombs, roads, etc. (4th–2nd centuries BC) in the neighbourhood. The archaeological remains were found at about 0.3 m depth.At present the site belongs to the municipality, which intends to build a new sewer network through it. The risk of destroying potentially interesting ancient archaeological structures during the works prompted an archaeological survey of the area. The relatively large dimensions of the area (almost 10,000 m²), together with time and cost constraints, made it necessary to use geophysical investigations as a faster means to ascertain the presence of archaeological items. Since the most important targets were expected to be located at a soil depth of about 0.3 m, a ground-penetrating radar (GPR) survey was carried out in an area located near the archaeological excavations. Unfortunately the geological complexity did not allow an easy interpretation of the GPR data.Therefore a 3D electrical resistivity tomography (ERT) scan was conducted in order to resolve these interpretation problems.A three-way comparison of the results of the dense ERT measurements parallel to the x axis, the results of the measurements parallel to the y axis and the combined results was performed.Subsequently the synthetic model approach was used to provide a better characterization of the resistivity anomalies visible on the ERT field data.The 3D inversion results clearly illustrate the capability to resolve in view of quality 3D structures of archaeological interest. According to the presented data the inversion models along one direction (x or y) seems to be adequate in reconstructing the subsurface structures.Naturally field data produce good quality reconstructions of the archaeological features only if the x-line and y-line measurements are considered together. Despite the increased computational time required by the 3D acquisition and 3D inversion schemes, good quality results can be produced.     

Seismic hazard impact of the Lower Tagus Valley Fault Zone (SW Iberia)

The seismic hazard of SW Iberia is composed of two contributions: offshore, large to very large events on the plate boundary between Africa and Eurasia such as the Lisbon earthquake of 1755 or the Gorringe Bank earthquake of 1969; and onshore, moderate to strong intraplate earthquakes on inherited crustal fractures. One of these zones of crustal weakness is the Lower Tagus Valley (LTV) fault zone, which displays the highest level of seismic hazard in Western Iberia. In this paper we review the active tectonics and seismicity of the LTV, integrating previous geophysical data with recent results of paleoseismological investigations, and discuss its impact on the seismic hazard of SW Iberia. We conclude that the seismic zonation for hazard assessment currently in force in the building code is biased towards the scenario of distant offshore rupture, and does not take adequately into account the LTV seismic source.     

Corona Satellite Pictures for Archaeological Studies: A Review and Application to the Lost Forbidden City of the Han–Wei Dynasties

This paper provides a state of the art of the use of declassified satellite Corona for archaeological application and deals with the use of historical archives of Corona and aerial photographs for detecting ancient hidden cultural relics that may be completely lost in the modern landscape veiled by current land uses and land covers. Geospatial analysis is adopted to identify, extract and map the subtle archaeological features and cultural landscape of the lost Han–Wei Forbidden City in the Henan Province (China). It was one of the most important cities of ancient China, served as the capital for a long period, from the Eastern Zhou (2600 BP, Zhou Dynasty) to the Eastern Han and Wei Dynasties (1900–1500 BP). Outputs from our investigations clearly provided evidence of numerous buried structures of potential archaeological interest that were compared with historical documentation and excavations carried out in the area. The geospatial analysis, applied both to aerial photograph and to Corona pictures, enabled us to discover a number of unknown sub-palaces and structures, thus providing new insights for this important archaeological area.     

安徽固镇垓下遗址自然变形遗迹的发现及特征初探

古遗址、墓地发掘过程中自然变形的识别研究是反演或充实历史或史前地震变形资料的重要工作。地震部门与文物考古部门合作,对安徽固镇垓下遗址考古现场进行了专题发掘研究,发现了断层和裂缝,初步研究显示为同期2次层位变动事件的遗迹,发生时代约相当于大汶口文化晚期,该工作有助于推动中国东部人类文化层中地震遗迹的识别     

2-D and 3-D processing and interpretation of multi-fold ground penetrating radar data: a case history from an archaeological site

A 2.5-D and 3-D multi-fold GPR survey was carried out in the Archaeological Park of Aquileia (northern Italy). The primary objective of the study was the identification of targets of potential archaeological interest in an area designated by local archaeological authorities. The second geophysical objective was to test 2-D and 3-D multi-fold methods and to study localised targets of unknown shape and dimensions in hostile soil conditions. Several portions of the acquisition grid were processed in common offset (CO), common shot (CSG) and common mid point (CMP) geometry. An 8×8 m area was studied with orthogonal CMPs thus achieving a 3-D subsurface coverage with azimuthal range limited to two normal components. Coherent noise components were identified in the pre-stack domain and removed by means of FK filtering of CMP records. Stack velocities were obtained from conventional velocity analysis and azimuthal velocity analysis of 3-D pre-stack gathers. Two major discontinuities were identified in the area of study. The deeper one most probably coincides with the paleosol at the base of the layer associated with activities of man in the area in the last 2500 years. This interpretation is in agreement with the results obtained from nearby cores and excavations. The shallow discontinuity is observed in a part of the investigated area and it shows local interruptions with a linear distribution on the grid. Such interruptions may correspond to buried targets of archaeological interest. The prominent enhancement of the subsurface images obtained by means of multi-fold techniques, compared with the relatively poor quality of the conventional single-fold georadar sections, indicates that multi-fold methods are well suited for the application to high resolution studies in archaeology.     

OSL dating of fine and coarse quartz from a Palaeolithic sequence on the Bistrița Valley (Northeastern Romania)

Previous works focused on optically stimulated luminescence dating of quartz extracted from Romanian and Serbian loess reported significant discrepancies between ages obtained on fine (4–11 μm) and coarse (63–90 μm) quartz. The present study is directed at expanding these investigations. The SAR-OSL and double SAR-OSL protocols are applied on quartz of different grain sizes belonging to 9 samples extracted from a newly identified archaeological site at Bistricioara-Lutărie III on the Bistrița Valley (NE Romania). Radiocarbon ages are also obtained for the 3 uppermost cultural layers hosted in the loess-like deposit.Discrepant ages are obtained between fine (4–11 μm) and coarse (63–90 and 90–200 μm, respectively) quartz for equivalent doses higher than ∼80–100 Gy. However, a very good agreement is achieved for the youngest sample, with an age of ∼8 ka (with a D_e of 38 Gy for fine and 35 Gy for coarse grains, respectively). The comparison with independent control provided by radiocarbon dating suggests better agreement of ages calculated for coarse quartz. Our results are once again proof that concerns should be raised regarding the reliability of the equivalent doses obtained on quartz samples for which the laboratory dose response cannot be fitted by a single saturating exponential function. Further systematic investigations are required regarding the very different saturation characteristics of fine and coarse grained quartz.     

Mapping landscape evolution in 3D: Climate change,natural hazard and human settlements across the 1618 Piuro landslide in the Italian Central Alps

Natural and anthropogenic mountain landscapes coevolve responding on different temporal scales to climate changes and geodynamics by a series of increments that cause the dynamic association of morphological stabilization surfaces, stratigraphic units and landforms. Understanding the incremental history of palimpsest landscapes helps to recognize and forecast the effects of climate change on the sensitive mountain environments, contributes to archaeological and historical reconstruction and supports management strategies for natural risks prevention and mitigation. The Italian Bregaglia Valley provides an excellent site to unravel the recent/historical increments of evolution of landforms and human settlement, permitting to map the paleo-digital terrain models (DTMs) corresponding to the relevant landscape turning points. After the last de-glaciation, two large-scale landslides reshaped the valley floor, both predisposed by deep-seated gravitational slope deformations and one surely triggered by intense rainfalls. The most recent and impacting event buried in 1618 the rich border town of Piuro, the ancient one occurred in the same area at least 1.5 ka before. Combining stratigraphic, geomorphological, topographic, archaeological and historical data, we drew the paleo-DTMs of the pre- and post-1618 settings of Piuro, sketching the landscape evolution. Since two millennia, human settlements took advantage of the decadal to secular most stable surfaces, represented by the inactive lobes of debris-flow fans, the highest trunk river terraces and the top of humps formed by the ancient landslide body in the valley centre. Stratigraphic relationships, archaeological findings and age determinations show that both landslides diverted the trunk river and covered the existing fan lobes. On a secular timescale, fan progradation and trunk river terracing buried and reworked both the landslide bodies. The paleo-DTMs show their original areal extent and permit to compute their volume and to sketch the setting of the buried Piuro settlements, drawing the changes of the Mera trunk river course and the chronology of activity of the lateral debris-flow fan lobes.     

An algorithm for fast 3D inversion of surface electrical resistivity tomography data: application on imaging buried antiquities

In this work a new algorithm for the fast and efficient 3D inversion of conventional 2D surface electrical resistivity tomography lines is presented. The proposed approach lies on the assumption that for every surface measurement there is a large number of 3D parameters with very small absolute Jacobian matrix values, which can be excluded in advance from the Jacobian matrix calculation, as they do not contribute significant information in the inversion procedure. A sensitivity analysis for both homogeneous and inhomogeneous earth models showed that each measurement has a specific region of influence, which can be limited to parameters in a critical rectangular prism volume. Application of the proposed algorithm accelerated almost three times the Jacobian (sensitivity) matrix calculation for the data sets tested in this work. Moreover, application of the least squares regression iterative inversion technique, resulted in a new 3D resistivity inversion algorithm more than 2.7 times faster and with computer memory requirements less than half compared to the original algorithm. The efficiency and accuracy of the algorithm was verified using synthetic models representing typical archaeological structures, as well as field data collected from two archaeological sites in Greece, employing different electrode configurations. The applicability of the presented approach is demonstrated for archaeological investigations and the basic idea of the proposed algorithm can be easily extended for the inversion of other geophysical data.     

Measurement and calibration of weathering using the schmidt hammer,on wave washed moraines on the upper norrland coast,Sweden

The Schmidt hammer has for some years been used to measure the hardness of different rock surfaces and to date relatively boulders in moraine deposits. In this paper applications on isostatically raised boulder-beaches and rocky shores are described. The research area is the northern Swedish Bothnian coast where present isostatic uplift is between 0.8-1.0 cm year^?1. Elevations above sea-level can be converted to a timescale. Rebound values (R) of the Schmidt hammer are correlated with the altitude and the degree of rock weathering can be estimated for each site. Eight sites with boulder-beaches and nine sites with rocky shores have been analysed. For boulder-beaches correlation-coefficients of 0.82-0.97 for weathering values (R) and the elevations above sea-level were obtained. For the rocky shores these values were calculated to 0.75-0.97. This method has applications in both geomorphology and archaeology. For example, a deep layer of weathered dolerite in the Nordingrå region, which previously has been regarded as a post-Weichselian formation, can be shown to be of pre-Weichselian origin. Several archaeological applications are also mentioned in the paper.     

AMS dating of insect chitin – A discussion of new dates,problems and potential

Results from AMS dating applied to insect chitin from a variety of contexts and different preservation conditions and retrieval methods are presented. Secure contexts, which include other dated organic material from different geographic locations ranging from Egypt to Greenland and different chronological periods, from Lateglacial to Medieval, have been used. In addition, insect species with different dietary requirements have been selected for dating purposes in order to provide an understanding as to whether diet plays a role in the chitin dating results. Dates from each context/site are discussed separately in the context of their stratigraphy and/or archaeology. Our research concentrates on the results from pre-treatment methods which require small quantities of chitin as these could be applied in a variety of Quaternary and archaeological contexts. The dates from carbonised and desiccated remains where no chemicals had been involved in storage fell within the range of dates from other organics or the archaeology. Although some of the dates from waterlogged contexts were successful, problems were encountered and these have been linked with long term storage in various alcohols of uncertain provenance. Whilst short term immersion in paraffin (kerosene) and alcohols during processing probably has no impact, it is recommended that chitinous material for dating be stored in acidified distilled water. Our results demonstrate the potential of chitin as a dating medium and provide a basis for its wider application.     

Non-invasive Geophysical Surveys in Search of the Roman Temple of Augustus Under the Cathedral of Tarragona (Catalonia,Spain): A Case Study

An integrated geophysical survey has been conducted at the Tarragona’s Cathedral (Catalonia, NE Spain) with the aim to confirm the potential occurrence of archaeological remains of the Roman Temple dedicated to the Emperor Augustus. Many hypotheses have been proposed about its possible location, the last ones regarding the inner part of the Cathedral, which is one of the most renowned temples of Spain (twelfth century) evolving from Romanesque to Gothic styles. A geophysical project including electrical resistivity tomography (ERT) and ground probing radar (GPR) was planned over 1 year considering the administrative and logistic difficulties of such a project inside a cathedral of religious veneration. Finally, both ERT and GPR have been conducted during a week of intensive overnight surveys that provided detailed information on subsurface existing structures. The ERT method has been applied using different techniques and arrays, ranging from standard Wenner–Schlumberger 2D sections to full 3D electrical imaging with the advanced Maximum Yield Grid array. Electrical resistivity data were recorded extensively, making available many thousands of apparent resistivity data to obtain a complete 3D image after a full inversion. In conclusion, some significant buried structures have been revealed providing conclusive information for archaeologists. GPR results provided additional information about shallowest structures. The geophysical results were clear enough to persuade religious authorities and archaeologists to conduct selected excavations in the most promising areas that confirmed the interpretation of geophysical data. In conclusion, the significant buried structures revealed by geophysical methods under the cathedral were confirmed by archaeological digging as the basement of the impressive Roman Temple that headed the Provincial Forum of Tarraco, seat of the Concilium of Hispania Citerior Province.     

Deformation of “Villafranchian” lacustrine sediments in the Chisone Valley (Western Alps, Italy)

The Chisone Valley is located in the internal NW Alps, in the Pinerolese District, an area characterized by present low to medium seismicity. Fine-grained sediments (sand, silt and clay with interbedded gravel) crop out in the lower Chisone Valley: they were first interpreted as glaciolacustrine deposits, and then as a lacustrine infilling of the valley floor probably due to differential uplifting of the valley mouth. Review of this data, together with new field and palynological observations, lead us to refer the lacustrine deposits to approximately the Lower Pleistocene (Villafranchian). In many outcrops, the lacustrine deposits show strong soft-sediment deformation such as convolute laminations, water-escape structures and disrupted beds, some of them associated with folds and faults (cm to dm in size); only two sites show metric to decametric folds and faults trending E-W and N-S. Detailed structural analysis conducted along a recently exposed section (Rio Gran Dubbione site) shows several soft-sediment deformation features on the limbs of mesoscale folds. Because of their intimate structural association, the origin of these minor structures seems to be connected to synsedimentary activity on reverse and normal faults (m to dm in size) affecting the lacustrine deposits in the same locality. Soft-sediment deformation features can be interpreted as possible paleoseismites. If so, the present seismicity of the Pinerolese District, which is the major area of such activity in NW Italy, cannot be considered an isolated episode in the geological evolution of the region; even if there is no supporting evidence for continuous seismicity, the deformations in the lacustrine sediments of the Chisone Valley testify to Early Pleistocene seismic activity, probably related to the recent tectonic evolution of the internal side of the NW Alps.     

Chronometry of pedogenic and stratigraphic events from calcite produced by earthworms

Dating phases of pedogenesis, soil carbonate deposition or even the burial of whole soil profiles using ¹⁴C is a valuable goal in archaeology and pedology, but one that has been consistently hampered by the presence of old carbon skewing the measurements to produce apparent dates older than the true formation date. Calcite produced by earthworms could be a useful alternative source of datable carbon. Since earthworms both inhabit and ingest soils with an old carbon content, however, the granules could yield a ¹⁴C date older than the date of their formation. In this study, by examining granules from two sites of known-age stratigraphy, we show that the radiocarbon date derived from the granules' calcite closely reflects their true formation date, opening up the possibility of using the granules either individually or as distributions of dates to understand soil processes and date sealed archaeological layers.     

Evaluation and archaeological potential of geophysical investigations of the Roman fort and baths in Northern Muntenia,România

The Roman fort from Sfârleanca is one of the most representative archaeological sites dating from the Roman period (2nd century A.D.) in Northern Muntenia. The existing natural and anthropic features of the environment required the application of geophysical methods in order to outline the spatial pattern of the buried remains, to define the geometry of the anthropogenic settlements and to obtain detailed information about different archaeological materials without digging. During the survey, two different geophysical methods have been employed: total magnetic field measurements and electrical resistance mapping using Twin-probe array. The instrument consists of GSM19W Overhauser magnetometers with GPS, in base-rover system, and a twin-probe array LGM 4-Point light hp. The measurements were used to draw primary maps of the physical parameters (total magnetic field strength/intensity, electrical resistance), and also processed maps (filtering, derivative). The magnetic results obtained by interpreting the anomalies yielded information about the limits of the fort, about the internal organization of the military structure (its axial road, partially its secondary road, the localization and the shape of its constructions) and at the bath and heating installation. A previously unknown element is the possible water supply pipe made of ceramic material highlighted by the mathematical modelling of the data obtained by the magnetic investigations. Electrical resistance results provide complementary information to the magnetic survey concerning the limits of the baths and the remains of the fort structure. This paper brings to light geophysical investigations into this Roman fort and baths, extending the picture produced by previous archaeological excavations that only dealt with a small part of the site. It indicates the importance of using geophysical methods in preliminary archaeological research and the advantages of combining total magnetic field measurements and electrical resistance mapping when investigating an archaeological site characterized by a number of environmental difficulties.