Archaeogeophysical-Based Approach for Inca Archaeology: Overview and one operational application

Even if, in recent decades, the use of remote sensing technologies (from satellite, aerial and ground) for archaeology is stepping into its golden age, in Southern America geophysics for preventive archaeology is more recent and less used than in Europe, Central America and Middle East. In this paper, we provide a brief overview and show the preliminary results obtained from the investigations conducted in Chachabamba (Peru). The archaeological area is located on a strategic terrace overlooking three Inca roads, which served the most important ceremonial centres (including Machu Picchu) of the Urubamba Valley also known as the Sacred Valley. In particular, Chachabamba investigations were conducted with two principal aims: (1) to give new impetus to archaeological research with targeted investigations aimed at improving and completing the site’s knowledge framework; (2) to experiment and validate an archaeogeophysical approach to be reapplied in other sites of the Urubamba valley, including Machu Picchu, having similar characteristics as those found in Chachabamba.     

Testing the accuracy of OSL and pIR IRSL dating of young geoarchaeological sediments in coastal Peru

Luminescence dating of late Quaternary sediments in Peru is challenging, especially on the Peruvian coast. Earlier studies have shown that quartz grains often exhibit a thermally unstable, medium signal that caused the underestimation of Optically Stimulated Luminescence (OSL) ages. InfraRed Stimulated Luminescence (IRSL) dating has shown to produce more reliable ages, depending, amongst other factors, on the age model (Central or Minimum Age Model), and the IRSL signal. IRSL dating of geoarchaeological sediments has, however, hardly been carried out, let along validated, against an independent age dataset. This dating approach is, nonetheless, the only promising way to date the geological substrate in which many of Peru's archaeological sites are buried. Peru contains some of the oldest and most important archaeological heritage sites, yet not much is known of the environmental context in which many of its early civilizations prospered. A better understanding of which luminescence method works best could therefore help in a better understanding of the geological-stratigraphical context of many of Peru's sites.To investigate this matter more fully, we compared the luminescence dating results of seven sediment samples from the top layer of the Lima alluvial fan and from geoarchaeological layers of the Maranga Complex (San Miguel, Lima), with an independent dataset of sixteen ¹⁴C ages. Our results showed that the quartz OSL ages always underestimated the expected ages due to a signal dominated by medium and slow components, and that the post-IR IRSL₂₂₅ (pIR IRSL₂₂₅) and IRSL₅₀ ages of K-feldspars, on basis of the Central Age Model (CAM), always overestimated the expected ages. The Minimum Age Model (MAM) on the other hand, correctly predicted the expected ages for the early Holocene, Lima alluvial fan sediments using the pIR IRSL₂₂₅ signal of K-feldspars, and the late Holocene, geoarchaeological ages using the IRSL₅₀ signal.     

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.     

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.     

Ground Penetrating Radar (GPR) attribute analysis for archaeological prospection

We evaluate the applicability and the effectiveness of the GPR attribute analysis for archaeological purposes and we test the attribute analysis on GPR data obtained in the river harbor area of the Aquileia Archaeological Park, NE Italy, where cultural heritage of the Roman imperial period is buried at different depths beneath a silty loam layer at an average depth not greater than 3–4 m. A multi-attribute approach is used to characterize the subsurface through several attribute categories, including instantaneous, section and volume attributes applied to quantities related not only to the reflection amplitude, but also to frequency and phase or other more complex parameters. The results allow the identification of different buried archaeological remains, from brick walls to isolated stone chunks also in a very complex environment characterized by several interconnected or superimposed archaeological remains.     

ESR dating in Song Terus cave (East Java, Indonesia)

Electron spin resonance (ESR) dating on teeth has been used to determine the antiquity of the human presence in the area of Gunung Sewu (East Java, Indonesia). The Song Terus cave is one of the archaeological sites located in this area, which has yielded an important collection of fossil bones and prehistoric lithic industry. ESR dating shows that this cave belongs to a karstic system which has been in place since Middle Pleistocene. The results of coupled ESR/Uranium-series dating indicates closed system for uranium in tooth tissues yielding ages of between 216 and 392 ka for the oldest archaeological levels, confirming the possibility of finding Middle Pleistocene human occupations in such karstic sites in Java.     

Seismic gaps and source zones of recent large earthquakes in coastal Peru

The earthquakes of central coastal Peru occur principally in two distinct zones of shallow earthquake activity that are inland of and parallel to the axis of the Peru Trench. The interface-thrust (IT) zone includes the great thrust-fault earthquakes of 17 October 1966 and 3 October 1974. The coastal-plate interior (CPI) zone includes the great earthquake of 31 May 1970, and is located about 50 km inland of and 30 km deeper than the interface thrust zone. The occurrence of a large earthquake in one zone may not relieve elastic strain in the adjoining zone, thus complicating the application of the seismic gap concept to central coastal Peru. However, recognition of two seismic zones may facilitate detection of seismicity precursory to a large earthquake in a given zone; removal of probable CPI-zone earthquakes from plots of seismicity prior to the 1974 main shock dramatically emphasizes the high seismic activity near the rupture zone of that earthquake in the five years preceding the main shock. Other conclusions on the seismicity of coastal Peru that affect the application of the seismic gap concept to this region are: (1) Aftershocks of the great earthquakes of 1966, 1970, and 1974 occurred in spatially separated clusters. Some clusters may represent distinct small source regions triggered by the main shock rather than delimiting the total extent of main-shock rupture. The uncertainty in the interpretation of aftershock clusters results in corresponding uncertainties in estimates of stress drop and estimates of the dimensions of the seismic gap that has been filled by a major earthquake. (2) Aftershocks of the great thrust-fault earthquakes of 1966 and 1974 generally did not extend seaward as far as the Peru Trench. (3) None of the three great earthquakes produced significant teleseismic activity in the following month in the source regions of the other two earthquakes. The earthquake hypocenters that form the basis of this study were relocated using station adjustments computed by the method of joint hypocenter determination.     

Geologic Setting, Field Survey and Modeling of the Chimbote, Northern Peru, Tsunami of 21 February 1996

—Whereas the coast of Peru south of 10°S is historically accustomed to tsunamigenic earthquakes, the subduction zone north of 10°S has been relatively quiet. On 21 February 1996 at 21:51 GMT (07:51 local time) a large, tsunamigenic earthquake (Harvard estimate M _w = 7.5) struck at 9.6°S, 79.6°W, approximately 130 km off the northern coast of Peru, north of the intersection of the Mendaña fracture zone with the Peru–Chile trench. The likely mechanism inferred from seismic data is a low-angle thrust consistent with subduction of the Nazca Plate beneath the South American plate, with relatively slow rupture characteristics. Approximately one hour after the main shock, a damaging tsunami reached the Peruvian coast, resulting in twelve deaths. We report survey measurements, from 7.7°S to 11°S, on maximum runup (2–5m, between 8 and 10°S), maximum inundation distances, which exceeded 500 m, and tsunami sediment deposition patterns. Observations and numerical simulations show that the hydrodynamic characteristics of this event resemble those of the 1992 Nicaragua tsunami. Differences in climate, vegetation and population make these two tsunamis seem more different than they were. This 1996 Chimbote event was the first large (M _w >7) subduction-zone (interplate) earthquake between about 8 and 10°S, in Peru, since the 17th century, and bears resemblance to the 1960 (M _w 7.6) event at 6.8°S. Together these two events are apparently the only large subduction-zone earthquakes in northern Peru since 1619 (est. latitude 8°S, est. M _w 7.8); these two tsunamis also each produced more fatalities than any other tsunami in Peru since the 18th century. We concur with Pelayo and Wiens (1990, 1992) that this subduction zone, in northern Peru, resembles others where the subduction zone is only weakly coupled, and convergence is largely aseismic. Subduction-zone earthquakes, when they occur, are slow, commonly shallow, and originate far from shore (near the tip of the wedge). Thus they are weakly felt, and the ensuing tsunamis are unanticipated by local populations. Although perhaps a borderline case, the Chimbote tsunami clearly is another wake-up example of a "tsunami earthquake."     

Geophysical prospection,archaeological excavation,and dating in two Hellenic pyramids

Geophysical prospection employing magnetometry and electromagnetic measurements has been applied in and around two small sized pyramids of Hellenikon and Ligourio in Argolid, Greece. The magnetic anomalies appropriately assessed were interpreted as possible archaeological targets. Subsequent test excavations revealed the presence of room foundations and parts of walls, as well as a plethora of ceramic ware. Typological study of the ceramics classified them to as early as the proto-Helladic period and to as late as the first centuries A.D. The earlier periods have been also confirmed by a novel application of thermoluminescence (TL) dating of ceramics and the megalithic stones themselves. The present results question earlier attempts classifying these pyramids at the Classical period and favour much earlier periods.     

The Arequipa (Peru) earthquake of June 23, 2001

The Arequipa earthquake of 23 June 2001 hasbeen the largest earthquake (Mw = 8.3)occurredin the last century in southern Peru with amaximum intensity of VIII (MM scale). Focalmechanisms of main shock and three largeraftershocks have been studied, showingthrusting solutions for main shock and twoaftershocks and normal motion for the eventof July, 5. The rupture area has beenobtained from distribution of aftershocks.The occurrence of the Arequipa earthquakeis related with the convergence processbetween the Nazca and South Americaplates.     

电容耦合式电阻率法在城市遗址考古勘探中的应用

城市密集柏油公路网和坚硬水泥地面给地下遗址考古勘探造成了很大困难,传统的高密度电阻率法需要插置电极,无法在城市环境中使用。电容耦合式电阻率法采用非侵入式的方式,通过特制电极与大地形成电容器产生电流进行电法勘探,特别适合城市硬地表介质的条件。为了检验电容耦合式电阻率法在城市遗址考古勘探中的应用效果,在地处杭州市主城区的吴越捍海塘遗址进行了电容耦合式电阻率法探测试验,采集多个收发距的电法数据,对数据进行反演处理,将反演结果与考古发掘结果进行对比分析。结果表明,反演的电阻率图像推断结果与考古发掘结果基本吻合,电容耦合式电阻率法是城市遗址考古勘探的有效方法,具有广泛的应用前景。      

Temporal changes of warm-season pan evaporation in a semi-arid basin in Western Turkey

Warm-season (from May to September) pan evaporation series recorded at nine sites in semi-arid Büyük Menderes Basin in Western Turkey were analyzed for evidence of long-term changes using non-parametric Mann–Kendall test. Only two sites showed significant linear trends, one in decreasing direction and the other in increasing direction. Area-averaged normalized anomalies of pan evaporation stayed unchanged over the period 1975–2006. Based on correlation and multiple regression analysis with visual inspection of the time series plots, sunshine duration which is used as a proxy for global radiation was found to be the major factor explaining the change and variability of pan evaporation in warm-season in semi-arid Büyük Menderes Basin. Stepwise multiple regression analysis showed no significant effect of rainfall on pan evaporation, which means that the effect of complementary relationship is rather weak and even non-existent.     

Archaeological evidence for a destructive earthquake in Patras,Greece

Oriented collapse of columns, large-scale destruction debris and temporary abandonment of the area deduced from an archaeological excavation provide evidence for a major (intensity IX) earthquake in Patras, Greece. This, and possibly a cluster of other earthquakes, can be derived from archaeological data. These earthquakes are not included in the historical seismicity catalogues, but can be used to put constraints to the seismic risk of this city. Patras was affected by a cluster of poorly documented earthquakes between 1714 and 1806. The city seems to be exposed to risks of progressive reactivation of a major strike-slip fault. A magnitude 6.4 earthquake in 2008 has been related to it. This fault has also been associated with a total of four events in the last 20 years, a situation reminiscent of the seismic hazard at the western edge of the North Anatolian Fault.     

GPR prospecting in a prehispanic village,NW Argentina

Palo Blanco is an approximately 1600 year old archaeological site located in Fiambalá Valley, in the Andean region of Argentina. Pioneer archaeological studies published in the 70s reported the existence of five residential units in this site. Also a small cemetery which included three circular tombs was discovered near to these buildings. Since that time, a profuse sedimentation covered the zone, so nowadays there are no evidences of most of these buildings on surface. Because of an imprecise location of the structures, most of the buildings became in fact missed. Then, in this work we aimed to re-localize two of the missed buildings, a residential unit and a tomb, by applying ground penetrating radar (GPR) methodology. We used fast fixed offset GPR configurations to investigate two areas in which these buildings could be expected. We used experimental and synthetic patterns to aid the identification of the characteristic signals due to the archaeological targets. We applied migration to the data to better define and resolve unclear anomalous signals. The employed methodology revealed the location of both buildings. Also a number of new non-reported structures were predicted and confirmed.     

Identification of modal quantities from two earthquake responses

The work deals with the identification of modal parameters of a structure from earthquake records when the input ground motion is unknown. This may occur, for example, owing to instrumental malfunctions. The procedure is based on the assumption that at least two responses are available and consists of two main steps. In the first one, modal frequencies are estimated by searching relative minima of a function that involves the ratio of the Fourier amplitudes of the two records, while the second phase is devoted to the identification of other modal quantities (i.e. effective participation factors and modal dampings). Once the identification process has been completed, an estimate of the unknown base input may be performed by means of the IFFT algorithm. The proposed approach has been checked against both finite element simulations of simple structures and field measurements on real buildings.     

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.     

Identification of coal seam strata from geophysical logs of borehole using Adaptive Neuro-Fuzzy Inference System

Different parameters obtained through well-logging geophysical sensors such as SP, resistivity, gamma–gamma, neutron, natural gamma and acoustic, help in identification of strata and estimation of the physical, electrical and acoustical properties of the subsurface lithology. Strong and conspicuous changes in some of the log parameters associated with any particular stratigraphy formation, are function of its composition, physical properties and help in classification. However some substrata show moderate values in respective log parameters and make difficult to identify or assess the type of strata, if we go by the standard variability ranges of any log parameters and visual inspection. The complexity increases further with more number of sensors involved.An attempt is made to identify the type of stratigraphy from borehole geophysical log data using a combined approach of neural networks and fuzzy logic, known as Adaptive Neuro-Fuzzy Inference System. A model is built based on a few data sets (geophysical logs) of known stratigraphy of in coal areas of Kothagudem, Godavari basin and further the network model is used as test model to infer the lithology of a borehole from their geophysical logs, not used in simulation. The results are very encouraging and the model is able to decipher even thin cola seams and other strata from borehole geophysical logs. The model can be further modified to assess the physical properties of the strata, if the corresponding ground truth is made available for simulation.     

Integrated Geophysical Investigations at the Greek Kamarina Site (Southern Sicily,Italy)

Kamarina, located in southern Sicily (Italy), was an important Greek colony since its foundation in the sixth century BC. Archaeological excavations, carried out since the twentieth century, uncovered only limited portions of the site so far. Despite the importance of the Greek colony, the presence of remarkable buildings that archaeologists expected to bring to light has not found fully correspondence in the archaeological excavations. Consequently, the integrated geophysical prospection carried out in the study area is aimed to support and address the future archaeological investigations. After the photographic and thermographic survey obtained by an unmanned aerial vehicle, we performed a systematic survey through ground magnetic and GPR methods over an area of 6200 m². The acquisition procedures have been optimized in order to get the best results combining high resolution and elevated speed of acquisition. The results derived from the three geophysical techniques have been conveniently combined by means of a cluster analysis, allowing us to clearly identify a series of buried archaeological features. Because of their geometrical characteristics, often in good agreement with the spatial arrangement of the archaeological remains at the surface, these buried archaeological features can be interpreted as roads, walls, or buildings foundations in which the various construction phases of the city can be clearly recognized. The integrated approach has proven to be essential for a robust interpretation of the archaeogeophysical investigation.     

Review of archaeointensity methods

Most of the traditional methods of determining the intensity of the ancient geomagnetic field from archaeological materials utilized thermal demagnetization of the natural remanent magnetization (NRM) and of the laboratory induced thermoremanent magnetization (TRM). When applied rigorously these methods are foolproof. They are, however, very time consuming and the number of samples with which they can be used is limited. Attempts to speed up these traditional methods have generally led to the use of subjective criteria in assessing the reliability of the results and archaeomagnetic research has recently been concentrated on extending the range of samples to which the method can be applied. Through the use of alternating field, rather than thermal, demagnetization of NRM and TRM it has become possible to apply corrections for alteration occurring during laboratory firing of the archaeological samples and develop objective criteria of reliability. Recent research has shown that it may be possible to determine archaeointensities the laboratory redeposition of lake sediments.