Gradiometria magnética e GPR aplicados à arqueologia em sítio estuarino em Penalva - MA

Esta pesquisa apresenta os resultados da aplicação de métodos geofísicos não invasivos, visando a indicação de possíveis locais de deposição de materiais, no levantamento arqueológicos em um sítio estuarino localizado no município de Penalva, no estado do Maranhão. Os métodos utilizados foram Magnetometria e Radar de Penetração no Solo (GPR). As medidas magnéticas foram utilizadas para indicar os locais para escavação. O GPR foi utilizado para confirmar as anomalias geradas pelo método magnético, proporcionando a diminuição de erros quando as anomalias magnéticas não são provenientes de feições arqueológicas. Foram realizados 14 perfis de GPR, utilizando uma antena de 400 MHz e 1664 medidas magnéticas, utilizando o magnetômetro. Durante o levantamento arqueológico realizado, baseando-se nas anomalias fornecidas pela geofísica, a equipe liderada pelos arqueólogos Deusdédit Carneiro Leite Filho e Fernando Luiz Tavares Marques realizou cinco escavações no sítio, o que lhes permitiu o resgate de várias peças cerâmicas. Os resultados apresentados nesta pesquisa permitiram verificar e avaliar a potencialidade da prospecção geofísica como método auxiliar em estudos arqueológicos.

Investigação geofísica ambiental e forense nos cemitérios do Benguí e do Tapanã (Belém - PA)

A Geofísica é uma ferramenta útil no estudo da contaminação de água provocada por cemitérios, assim como na investigação criminalística. A contaminação provocada por cemitérios é um tema ainda pouco estudado. Por outro lado, Geofísica forense é quase desconhecida no Brasil. Para estudar esses dois temas – cemitérios e alvos forenses – foram eleitos dois cemitérios públicos: o cemitério do Bengui, já desativado, e o cemitério do Tapanã. A investigação geofísica nos dois Cemitérios foi realizada por meio dos métodos eletromagnéticos ground penetrating radar (GPR) e Slingram de baixo número de indução (low induction number – LIN). No Cemitério do Bengui, testou-se o GPR na detecção do nível hidrostático, visto que, na implantação de um Cemitério, a informação mais importante é a profundidade do aquífero. No Cemitério do Tapanã, onde foi implantado o Campo de Testes controlados para Geofísica Forense, Ambiental e de Resgate, denominado FORAMB, utilizou-se GPR e Slingram em levantamentos mensais a fim de detectar o nível hidrostático, monitorar o desenvolvimento da contaminação e, especialmente, encontrar alvos forenses (um cadáver, uma caixa simulando um túnel e uma caixa com ferros simulando armas). O GPR mostrou-se eficaz em todos os estudos, especialmente na delineação dos alvos forenses. O Slingram não forneceu informações úteis.

Investigação com os métodos GPR e magnético de urnas funerárias do Sítio Nossa Senhora do Perpétuo Socorro (Pará)

Este trabalho apresenta os resultados da aplicação de métodos geofísicos na eleição de áreas potenciais para salvamento arqueológico no Sítio Nossa Senhora do Perpétuo Socorro no estado do Pará. O Sítio possui urnas funerárias indígenas aflorantes e sub-aflorantes de bastante interesse arqueológico. Os métodos utilizados foram Magnetometria e Radar de Penetração no Solo (GPR). As medidas magnéticas foram utilizadas para indicar os locais mais apropriados para escavação reduzindo bastante as áreas potenciais. O GPR foi usado para confirmar as anomalias detectadas pelas medidas magnéticas e proporcionar uma melhor avaliação espacial tanto horizontal quanto verticalmente, diminuindo os erros cometidos quando se identificam anomalias magnéticas que não são causadas por feições arqueológicas. A metodologia de estudo baseou-se na comparação de anomalias obtidas em terrenos desconhecidos com terrenos onde comprovadamente havia a existência de elementos mapeáveis como urnas aflorantes e raízes, caracterizando e agrupando essas anomalias.

O método GPR aplicado à localização de tubulações utilizadas no abastecimento de água na região urbana do município de Belém - Pará

Este trabalho consistiu em avaliar a viabilidade do método GPR (Ground Penetrating Radar) na localização de tubulações de aço, cimento-amianto, PVC e ferro fundido, de diversos diâmetros, utilizadas pela Companhia de Saneamento do Pará (COSANPA), no abastecimento de água à população da região urbana do município de Belém-Pará. Para o alcance deste objetivo foram realizados 14 perfis em 6 locais, distintos da área urbana de Belém, com o intuito de identificar 13 alvos, utilizando-se antenas de 200 MHz e 400 MHz, geralmente usadas neste tipo de levantamento. Os radargramas obtidos com o GPR permitiram identificar quase todos os alvos pesquisados. Apenas 3 dos 13 alvos, não foram plenamente identificados, ficando de fora, apenas os alvos de PVC, de pequena bitola (50 mm e 75 mm). Após o estudo concluímos que o método GPR é muito eficaz neste tipo de aplicação e que sua utilização tem grande valia na localização das redes de abastecimento de água. O que deve facilitar enormemente as obras de expansão da rede, assim como a sua manutenção.

Metodologia geofísica aplicada à pesquisa de salvamento arqueológico nos sítios Bittencourt, Jambuaçu e Jaburu, Estado do Pará

Este trabalho apresenta os resultados da aplicação de métodos geofísicos no salvamento arqueológico em três sítios no estado do Pará: Bittencourt, localizado em Abaetetuba; Jambuaçu, localizado em Moju; e Jaburu, localizado em Almeirim. Os métodos utilizados foram Magnetometria, Cintilometria e Radar de Penetração no Solo (GPR). As medidas magnéticas foram utilizadas para indicar os locais para escavação. O GPR foi usado para confirmar as anomalias detectadas pelas medidas magnéticas e proporcionou a diminuição de erros cometidos quando se identificam anomalias magnéticas que não são causadas por feições arqueológicas. As medidas de cintilometria foram usadas para tentar mapear os solos de Terra Preta Arqueológica (TPA) através da radiação gama natural do solo. Os locais indicados para escavação pela magnetometria e GPR possibilitaram a descoberta de materiais arqueológicos como cerâmica, louça, lítico e ferro. Os resultados obtidos com a cintilometria foram satisfatórios, mostrando haver relação entre as medidas mais baixas de cintilometria e as áreas com maior espessura de TPA. A metodologia geofísica utilizada demonstrou eficiência e rapidez sem causar impactos destrutivos nos sítios.

Radarfacies and sedimentological analysis: study of sedimentary substrate from an archaeological site (shell mound), southern Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Neotectonics in the northern equatorial Brazilian margin

An increasing volume of publications has addressed the role of tectonics in inland areas of northern Brazil during the Neogene and Quaternary, despite its location in a passive margin. Hence, northern South America plate in this time interval might have not been as passive as usually regarded. This proposal needs further support, particularly including field data. In this work, we applied an integrated approach to reveal tectonic structures in Miocene and late Quaternary strata in a coastal area of the Amazonas lowland. The investigation, undertaken in Marajo Island, mouth of the Amazonas River, consisted of shallow sub-surface geophysical data including vertical electric sounding and ground penetrating radar. These methods were combined with morphostructural analysis and sedimentological/stratigraphic data from shallow cores and a few outcrops. The results revealed two stratigraphic units, a lower one with Miocene age, and an upper one of Late Pleistocene-Holocene age. An abundance of faults and folds were recorded in the Miocene deposits and, to a minor extent, in overlying Late Pleistocene-Holocene strata. In addition to characterize these structures, we discuss their origin, considering three potential mechanisms: Andean tectonics, gravity tectonics related to sediment loading in the Amazon Fan, and rifting at the continental margin. Amongst these hypotheses, the most likely is that the faults and folds recorded in Marajo Island reflect tectonics associated with the history of continental rifting that gave rise to the South Atlantic Ocean. This study supports sediment deposition influenced by transpression and transtension associated with strike-slip divergence along the northern Equatorial Brazilian margin in the Miocene and Late Pleistocene-Holocene. This work records tectonic evidence only for the uppermost few ten of meters of this sedimentary succession. However, available geological data indicate a thickness of up to 6 km, which is remarkably thick for an area regarded as a passive margin. (C) 2012 Elsevier Ltd. All rights reserved.

Dynamics of the cold surface layer of polythermal Storglaciären, Sweden

Polythermal glaciers, i.e. glaciers with a combination of ice at and below the freezing point, are widespread in arctic and subarctic environments. The polythermal structure has major implications for glacier hydrology, ice flow and glacial erosion. However, the interplay of factors governing its spatial and temporal variations such as net mass balance, ice advection and water content in the ice is poorly investigated and as yet not fully understood. This study deals with a thorough investigation of the polythermal regime on Storglaciären, northern Sweden, a small valley glacier with a cold surface layer in the ablation area. Extensive field work was performed including mapping of the cold surface layer using ground-penetrating radar, ice temperature measurements, mass balance and ice velocity measurements. Analyses of these data combined with numerical modelling were used specifically to investigate the spatial and temporal variability of the cold surface layer, the spatial distribution of the water content just below the cold surface layer transition, the effect of radar frequency on the detection of the surface layer, and the sensitivity of the cold surface layer to changes in forcing. A comparison between direct temperature measurements in boreholes and ground-penetrating surveys shows that the radar-inferred cold-temperate transition depth is within ±1 m from the melting point of ice at frequencies above ~300 MHz. At frequencies below ~155 MHz, the accuracy degrades because of reduced scattering efficiency that occurs when the scatterers become much smaller compared to the wavelength. The mapped spatial pattern of the englacial cold-temperate transition boundary is complex. This pattern reflects the observed spatial variation in net loss of ice at the surface by ablation and vertical advection of ice, which is suggested to provide the predominant forcing of the cold surface layer thickness pattern. This is further supported by thermomechanical modeling of the cold surface layer, which indicates high sensitivity of the cold surface layer thickness to changes in vertical advection rates. The water content is the least investigated quantity that is relevant for the thermal regime of glaciers, but also the most difficult to assess. Spatial variability of absolute water content in the temperate ice immediately below the cold surface layer on Storglaciären was determined by combining relative estimates of water content from ground-penetrating radar data with absolute determination from temperature measurements and the thermal boundary condition at the freezing front. These measurements indicate large-scale spatial variability in the water content, which seems to arise from variations in entrapment of water at the firn-ice transition. However, this variability cannot alone explain the spatial pattern in the thermal regime on Storglaciären. Repeated surveys of the cold surface layer show a 22% average thinning of the cold surface layer on Storglaciären between 1989 and 2001. Transient thermomechanical modeling results suggest that the cold surface layer adapts to new equilibrium conditions in only a few decades after a perturbation in the forcing is introduced. An increased winter air temperature since mid-1980s seems to be the cause of the observed thinning of the cold surface layer. Over the last decades, mass balance measurements indicate that the glacier has been close to a steady state. The quasi-steady state situation is also reflected in the vertical advection, which shows no significant changes during the last decades. Increased winter temperatures at the ice surface would result in a slow-down of the formation of cold ice at the base of the cold surface layer and lead to a larger imbalance between net loss of ice at the surface and freezing of temperate ice at the cold-temperate transition.

Applicazione ed integrazione di metodi di prospezione geofisica: studio idrogeofisico della Zona Insatura in condizioni idriche controllate

Il lavoro è incentrato sull’applicazione ed integrazione di differenti tecniche di indagine geofisica in campo ambientale e ingegneristico/archeologico. Alcuni esempi sono stati descritti al fine di dimostrare l’utilità delle metodologie geofisiche nella risoluzione di svariate problematiche. Nello specifico l’attenzione è stata rivolta all’utilizzo delle tecniche del Ground Penetrating Radar e del Time Domain Reflectometry in misure condotte su un corpo sabbioso simulante una Zona Insatura. L’esperimento è stato realizzato all’interno di un’area test costruita presso l’azienda agricola dell’Università La Tuscia di Viterbo. Hanno partecipato al progetto le Università di Roma Tre, Roma La Sapienza, La Tuscia, con il supporto tecnico della Sensore&Software. Nello studio è stato condotto un approccio definito idrogeofisico al fine di ottenere informazioni da misure dei parametri fisici relativi alla Zona Insatura simulata nell’area test. Il confronto e l’integrazione delle due differenti tecniche di indagine ha offerto la possibilità di estendere la profondità di indagine all’interno del corpo sabbioso e di verificare l’utilità della tecnica GPR nello studio degli effetti legati alle variazioni del contenuto d’acqua nel suolo, oltre a determinare la posizione della superficie piezometrica per i differenti scenari di saturazione. Uno specifico studio è stato realizzato sul segnale radar al fine di stabilire i fattori di influenza sulla sua propagazione all’interno del suolo. Il comportamento dei parametri dielettrici nelle condizioni di drenaggio e di imbibizione del corpo sabbioso è stato riprodotto attraverso una modellizzazione delle proprietà dielettriche ed idrologiche sulla base della dimensione, forma e distribuzione dei granuli di roccia e pori, nonché sulla base della storia relativa alla distribuzione dei fluidi di saturazione all’interno del mezzo. La modellizzazione è stata operata sulle basi concettuali del Differential Effective Medium Approximation.

Ancient Greek harbours used as geo-archives for palaeotsunami research : case studies from Krane (Cefalonia), Lechaion (Gulf of Corinth) and Kyllini (Peloponnese)

Since historical times, coastal areas throughout the eastern Mediterranean are exposed to tsunami hazard. For many decades the knowledge about palaeotsunamis was solely based on historical accounts. However, results from timeline analyses reveal different characteristics affecting the quality of the dataset (i.e. distribution of data, temporal thinning backward of events, local periodization phenomena) that emphasize the fragmentary character of the historical data. As an increasing number of geo-scientific studies give convincing examples of well dated tsunami signatures not reported in catalogues, the non-existing record is a major problem to palaeotsunami research. While the compilation of historical data allows a first approach in the identification of areas vulnerable to tsunamis, it must not be regarded as reliable for hazard assessment. Considering the increasing economic significance of coastal regions (e.g. for mass tourism) and the constantly growing coastal population, our knowledge on the local, regional and supraregional tsunami hazard along Mediterranean coasts has to be improved. For setting up a reliable tsunami risk assessment and developing risk mitigation strategies, it is of major importance (i) to identify areas under risk and (ii) to estimate the intensity and frequency of potential events. This approach is most promising when based on the analysis of palaeotsunami research seeking to detect areas of high palaeotsunami hazard, to calculate recurrence intervals and to document palaeotsunami destructiveness in terms of wave run-up, inundation and long-term coastal change. Within the past few years, geo-scientific studies on palaeotsunami events provided convincing evidence that throughout the Mediterranean ancient harbours were subject to strong tsunami-related disturbance or destruction. Constructed to protect ships from storm and wave activity, harbours provide especially sheltered and quiescent environments and thus turned out to be valuable geo-archives for tsunamigenic high-energy impacts on coastal areas. Directly exposed to the Hellenic Trench and extensive local fault systems, coastal areas in the Ionian Sea and the Gulf of Corinth hold a considerably high risk for tsunami events, respectively.Geo-scientific and geoarcheaological studies carried out in the environs of the ancient harbours of Krane (Cefalonia Island), Lechaion (Corinth, Gulf of Corinth) and Kyllini (western Peloponnese) comprised on-shore and near-shore vibracoring and subsequent sedimentological, geochemical and microfossil analyses of the recovered sediments. Geophysical methods like electrical resistivity tomography and ground penetrating radar were applied in order to detect subsurface structures and to verify stratigraphical patterns derived from vibracores over long distances. The overall geochronological framework of each study area is based on radiocarbon dating of biogenic material and age determination of diagnostic ceramic fragments. Results presented within this study provide distinct evidence of multiple palaeotsunami landfalls for the investigated areas. Tsunami signatures encountered in the environs of Krane, Lechaion and Kyllini include (i) coarse-grained allochthonous marine sediments intersecting silt-dominated quiescent harbour deposits and/or shallow marine environments, (ii) disturbed microfaunal assemblages and/or (iii) distinct geochemical fingerprints as well as (iv) geo-archaeological destruction layers and (v) extensive units of beachrock-type calcarenitic tsunamites. For Krane, geochronological data yielded termini ad or post quem (maximum ages) for tsunami event generations dated to 4150 ± 60 cal BC, ~ 3200 ± 110 cal BC, ~ 650 ± 110 cal BC, and ~ 930 ± 40 cal AD, respectively. Results for Lechaion suggest that the harbour was hit by strong tsunami impacts in the 8th-6th century BC, the 1st-2nd century AD and in the 6th century AD. At Kyllini, the harbour site was affected by tsunami impact in between the late 7th and early 4th cent. BC and between the 4th and 6th cent. AD. In case of Lechaion and Kyllini, the final destruction of the harbour facilities also seems to be related to the tsunami impact. Comparing the tsunami signals obtained for each study areas with geo-scientific data from palaeotsunami events from other sites indicates that the investigated harbour sites represent excellent geo-archives for supra-regional mega-tsunamis.

Geophysical Analysis of Seasonal Montandon Gravel Ridge Water Table Fluctuation and Moisture Gradient Variation due to Storm Events

The purpose of this research project is to continue exploring the Montandon Long-Term Hydrologic Research Site(LTHR) by using multiple geophysical methods to obtain more accurate and precise information regarding subsurface hydrologic properties of a local gravel ridge,which are important to both the health of surrounding ecosystems and local agriculture. Through using non-invasive geophysical methods such as seismic refraction, Direct Current resistivity and ground penetrating radar (GPR) instead of invasive methods such as boreholedrilling which displace sediment and may alter water flow, data collection is less likely to bias the data itself. In addition to imaging the gravel ridge subsurface, another important researchpurpose is to observe how both water table elevation and the moisture gradient (moisture content of the unsaturated zone) change over a seasonal time period and directly after storm events. The combination of three types of data collection allows the strengths of each method combine together and provide a relatively strongly supported conclusions compared to previous research. Precipitation and geophysical data suggest that an overall increase in precipitation during the summer months causes a sharp decrease in subsurface resistivity within the unsaturated zone. GPR velocity data indicate significant immediate increase in moisture content within the shallow vadose zone (< 1m), suggesting that rain water was infiltrating into the shallow subsurface. Furthermore, the combination of resistivity and GPR results suggest that the decreased resistivity within the shallow layers is due to increased ion content within groundwater. This is unexpected as rainwater is assumed to have a DC resistivity value of 3.33*105 ohm-m. These results may suggest that ions within the sediment must beincorporated into the infiltrating water.

Monitoring Stormwater Redistribution Into Low Resistivity Soils Using Noninvasive Geophysical Techniques

Water held in the unsaturated zone is important for agriculture and construction and is replenished by infiltrating rainwater. Monitoring the soil water content of clay soils using ground-penetrating radar (GPR) has not been researched, as clay soils cause attenuation of GPR signal. In this study, GPR common-midpoint soundings (CMPs) are used in the clayey soils of the Miller Run floodplain to monitor changes in the soil water content (SWC) before and after rainfall events. GPR accomplishes this task because increases in water content will increase the dielectric constant of the subsurface material, and decrease the velocity of the GPR wave. Using an empirical relationship between dielectric constant and SWC, the Topp relation, we are able to calculate a SWC from these velocity measurements. Non-invasive electromagnetics, resistivity, and seismic were performed, and from these surveys, the layering at the field site was delineated. EM characterized the horizontal variation of the soil, allowing us to target the most clay rich area. At the CMP location, resistivity indicates the vertical structure of the subsurface consists of a 40 cm thick layer with a resistivity of 100 ohm*m. Between 40 cm and 1.5 m is a layer with a resistivity of 40 ohm*m. The thickness estimates were confirmed with invasive auger and trenching methods away from the CMP location. GPR CMPs were collected relative to a July 2013 and September 2013 storm. The velocity observations from the CMPs had a precision of +/- 0.001 m/ns as assessed by repeat analysis. In the case of both storms, the GPR data showed the expected relationship between the rainstorms and calculated SWC, with the SWC increasing sharply after the rainstorm and decreasing as time passed. We compared these data to auger core samples collected at the same time as the CMPs were taken, and the volumetric analysis of the cores confirmed the trend seen in the GPR, with SWC values between 3 and 5 percent lower than the GPR estimates. Our data shows that we can, with good precision, monitor changes in the SWC of conductive soils in response to rainfall events, despite the attenuation induced by the clay.

Investigating Peatland Stratigraphy and Hydrogeology Using Integrated Electrical Geophysics

Hydrology has been suggested as the mechanism controlling vegetation and related surficial pore-water chemistry in large peatlands. Peatland hydrology influences the carbon dynamics within these large carbon reservoirs and will influence their response to global warming. A geophysical survey was completed in Caribou Bog, a large peatland in Maine, to evaluate peatland stratigraphy and hydrology. Geophysical measurements were integrated with direct measurements of peat stratigraphy from probing, fluid chemistry, and vegetation patterns in the peatland. Consistent with previous field studies, ground-penetrating radar (GPR) was an excellent method for delineating peatland stratigraphy. Prominent reflectors from the peat-lake sediment and lake sediment-mineral soil contacts were precisely recorded up to 8 m deep. Two-dimensional resistivity and induced polarization imaging were used to investigate stratigraphy beneath the mineral soil, beyond the range of GPR. We observe that the peat is chargeable, and that IP imaging is an alternative method for defining peat thickness. The chargeability of peat is attributed to the high surface-charge density on partially decomposed organic matter. The electrical conductivity imaging resolved glaciomarine sediment thickness (a confining layer) and its variability across the basin. Comparison of the bulk conductivity images with peatland vegetation revealed a correlation between confining layer thickness and dominant vegetation type, suggesting that stratigraphy exerts a control on hydrogeology and vegetation distribution within this peatland. Terrain conductivity measured with a Geonics EM31 meter correlated with confining glaciomarine sediment thickness and was an effective method for estimating variability in glaciomarine sediment thickness over approximately 18 km(2). Our understanding of the hydrogeology, stratigraphy, and controls on vegetation growth in this peatland was much enhanced from the geophysical study.

GPR Reflection Profiles of Clark and Commonwealth Glaciers, Dry Valleys, Antarctica

Englacial horizons deeper than 100 m are absent within 100 MHz ground-penetrating radar (GPR) surface profiles we recorded on Clark and Commonwealth Glaciers in the Antarctic Dry Valleys region. Both glaciers show continuous bottom horizons to 280 m, with bottom signal-to-noise ratios near 30 dB. Density horizons should fade below 50 m depth because impermeable ice occurred by 36 m. Folding within Commonwealth Glacier could preclude radar strata beneath about 80 m depth, but there is no significant folding within Clark Glacier. Strong sulfate concentrations and contrasts exist in our shallow ice core. However, it appears that high background concentration levels, and possible decreased concentration contrasts with depth placed their corresponding reflection coefficients at the limit of, or below, our system sensitivity by about 77 m depth. Further verification of this conclusion awaits processing of our deep-core chemistry profiles.

Variability in Accumulation Rates from GPR Profiling on the West Antarctic Plateau

lsochronal layers in firn detected with ground-penetrating radar (GPR) and dated using results from ice-core analyses are used to calculate accumulation rates along a 100 km across-flow profile in West Antarctica. Accumulation rates are shown to be highly variable over short distances. Elevation measurements from global positioning system surveys show that accumulation rates derived from shallow horizons correlate well with surface undulations, which implies that wind redistribution of snow is the leading cause of this variability. Temporal changes in accumulation rate over 25-185 year intervals are smoothed to along-track length scales comparable to surface undulations in order to identify trends in accumulation that are likely related to changes in climate. Results show that accumulation rates along this profile have decreased in recent decades, which is consistent with core-derived time series of annual accumulation rates measured at the two ends of the radar profile. These results suggest that temporal variability observed in accumulation-rate records from ice cores and GPR profiles can be obscured by spatial influences, although it is possible to resolve temporal signals if the effects of local topography and ice flow are quantified and removed.