首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
We propose a Bayesian fusion approach to integrate multiple geophysical datasets with different coverage and sensitivity. The fusion strategy is based on the capability of various geophysical methods to provide enough resolution to identify either subsurface material parameters or subsurface structure, or both. We focus on electrical resistivity as the target material parameter and electrical resistivity tomography (ERT), electromagnetic induction (EMI), and ground penetrating radar (GPR) as the set of geophysical methods. However, extending the approach to different sets of geophysical parameters and methods is straightforward. Different geophysical datasets are entered into a trans-dimensional Markov chain Monte Carlo (McMC) search-based joint inversion algorithm. The trans-dimensional property of the McMC algorithm allows dynamic parameterisation of the model space, which in turn helps to avoid bias of the post-inversion results towards a particular model. Given that we are attempting to develop an approach that has practical potential, we discretize the subsurface into an array of one-dimensional earth-models. Accordingly, the ERT data that are collected by using two-dimensional acquisition geometry are re-casted to a set of equivalent vertical electric soundings. Different data are inverted either individually or jointly to estimate one-dimensional subsurface models at discrete locations. We use Shannon's information measure to quantify the information obtained from the inversion of different combinations of geophysical datasets. Information from multiple methods is brought together via introducing joint likelihood function and/or constraining the prior information. A Bayesian maximum entropy approach is used for spatial fusion of spatially dispersed estimated one-dimensional models and mapping of the target parameter. We illustrate the approach with a synthetic dataset and then apply it to a field dataset. We show that the proposed fusion strategy is successful not only in enhancing the subsurface information but also as a survey design tool to identify the appropriate combination of the geophysical tools and show whether application of an individual method for further investigation of a specific site is beneficial.  相似文献   

2.
Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten–Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.  相似文献   

3.
Methods for estimating the parameter distributions necessary for modeling fluid flow and contaminant transport in the shallow subsurface are in great demand. Soil properties such as permeability, porosity, and water retention are typically estimated through the inversion of hydrological data (e.g., measurements of capillary pressure and water saturation). However, ill-posedness and non-uniqueness commonly arise in such non-linear inverse problems making their solutions elusive. Incorporating additional types of data, such as from geophysical methods, may greatly improve the success of inverse modeling. In particular, ground-penetrating radar (GPR) methods have proven sensitive to subsurface fluid flow processes and appear promising for such applications. In the present work, an inverse technique is presented which allows for the estimation of flow parameter distributions and the prediction of flow phenomena using GPR and hydrological measurements collected during a transient flow experiment. Specifically, concepts from the pilot point method were implemented in a maximum a posteriori (MAP) framework to allow for the generation of permeability distributions that are conditional to permeability point measurements, that maintain specified patterns of spatial correlation, and that are consistent with geophysical and hydrological data. The current implementation of the approach allows for additional flow parameters to be estimated concurrently if they are assumed uniform and uncorrelated with the permeability distribution. (The method itself allows for heterogeneity in these parameters to be considered, and it allows for parameters of the petrophysical and semivariogram models to be estimated as well.) Through a synthetic example, performance of the method is evaluated under various conditions, and some conclusions are made regarding the joint use of transient GPR and hydrological measurements in estimating fluid flow parameters in the vadose zone.  相似文献   

4.
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.  相似文献   

5.
A high resolution geophysical survey was carried out on the Pont De Coq, a medieval stone arch bridge located in Normandy (France) in 2011 and 2012. Two complementary methods are used: Electrical Resistivity Tomography (ERT) and Ground PenetratingRadar (GPR). They allow to evaluate the structural state of the bridge and to characterize the subsurface around and beneath the bridge. An excellent correlation is obtained between the geophysical methods and the geological data obtained around the bridge. In order to improve the restitution of the geophysical data, an advanced photogrammetric method is performed, providing a high resolution 3D Digital Terrain Model (DTM) of the Pont de Coq. The advanced photogrammetry enhances the presentation of the GPR and ERT data. This approach is an easy-to-use, rapid and cost-effective tool for stakeholders. Finally, it is a promising and original method for improved interpretations of future geophysical surveys.  相似文献   

6.
The complex ecohydrological processes of rangelands can be studied through the framework of ecological sites (ESs) or hillslope‐scale soil–vegetation complexes. High‐quality hydrologic field investigations are needed to quantitatively link ES characteristics to hydrologic function. Geophysical tools are useful in this context because they provide valuable information about the subsurface at appropriate spatial scales. We conducted 20 field experiments in which we deployed time‐lapse electrical resistivity tomography (ERT), variable intensity rainfall simulation, ground‐penetrating radar (GPR), and seismic refraction, on hillslope plots at five different ESs within the Upper Crow Creek Watershed in south‐east Wyoming. Surface runoff was measured using a precalibrated flume. Infiltration data from the rainfall simulations, coupled with site‐specific resistivity–water content relationships and ERT datasets, were used to spatially and temporally track the progression of the wetting front. First‐order constraints on subsurface structure were made at each ES using the geophysical methods. Sites ranged from infiltrating 100% of applied rainfall to infiltrating less than 60%. Analysis of covariance results indicated significant differences in the rate of wetting front progression, ranging from 0.346 m min?1/2 for sites with a subsurface dominated by saprolitic material to 0.156 m min?1/2 for sites with a well‐developed soil profile. There was broad agreement in subsurface structure between the geophysical methods with GPR typically providing the most detail. Joint interpretation of the geophysics showed that subsurface features such as soil layer thickness and the location of subsurface obstructions such as granite corestones and material boundaries had a large effect on the rate of infiltration and subsurface flow processes. These features identified through the geophysics varied significantly by ES. By linking surface hydrologic information from the rainfall simulations with subsurface information provided by the geophysics, we can characterize the ES‐specific hydrologic response. Both surface and subsurface flow processes differed among sites and are directly linked to measured characteristics.  相似文献   

7.
综合地球物理技术在采空区的探测中发挥了重要作用.目前通常采用单方法反演、仅对不同方法反演结果进行对比解释的综合勘探方式,单方法反演的多解性严重降低了其探测精度.如何提高采空区的探测精度,对采空区进行有效探测一直被认为是地球物理技术面临的首要难题.为了提高地震与电法技术的探测精度,基于交叉梯度联合反演理论,设计了地震初至折射走时数据和高密度电法数据的联合反演算法流程,对采空区理论模型和野外实际数据进行了联合反演处理.结果发现通过两者的联合反演,不仅可以提高采空区电阻率反演模型的成像效果,而且能够获得地震单方法反演难以成像的采空区低速异常体,从而提高了地震与电法技术对采空区的探测精度.表明地震与电法探测数据联合反演是一种提高采空区探测精度的有效方法.  相似文献   

8.
Herbaceous vegetation in the Sahel grows almost exclusively on sandy soils which preferentially retain water through infiltration and storage. The hydrological functioning of these sandy soils during rain cycles is unknown. One way to tackle this issue is to spatialize variations in water content but these are difficult to measure in the vadose zone. We investigated the use of Electrical Resistivity Tomography (ERT) as a technique for spatializing resistivity in a non-destructive manner in order to improve our knowledge of relevant hydrological processes. To achieve this, two approaches were examined. First, we focused on a possible link between water tension (which is much easier to measure in the field by point measurements than water content), and resistivity (spatialized with ERT). Second, because ERT is affected by solution non-uniqueness and reconstruction smoothing, we improved the accuracy of ERT inversion by comparing calculated solutions with in-situ resistivity measurements. We studied a natural microdune during a controlled field experiment with artificial sprinkling which reproduced typical rainfall cycles. We recorded temperature, water tension and resistivity within the microdune and applied surface ERT before and after the 3 rainfall cycles. Soil samples were collected after the experiment to determine soil physical characteristics. An experimental relationship between water tension and water content was also investigated. Our results showed that the raw relationship between calculated ERT resistivity and water tension measurements in sand is highly scattered because of significant spatial variations in porosity. An improved correlation was achieved by using resistivity ratio and water tension differences. The slope of the relationship depends on the soil solution conductivity, as predicted by Archie's law when salted water was used for the rain simulation. We found that determining the variations in electrical resistivity is a sensitive method for spatializing the differences in water tension which are directly linked with infiltration and evaporation/drainage processes in the vadose zone. However, three factors complicate the use of this approach. Firstly, the relation between water tension and water content is generally non-linear and dependent on the water content range. This could limit the use of our site-specific relations for spatializing water content with ERT through tension. Secondly, to achieve the necessary optimization of ERT inversion, we used destructive resistivity measurements in the soil, which renders ERT less attractive. Thirdly, we found that the calculated resistivity is not always accurate because of the smoothing involved in surface ERT data inversion. We conclude that further developments are needed into ERT image reconstruction before water tension (and water content) can be spatialized in heterogeneous sandy soils with the accuracy needed to routinely study their hydrological functioning.  相似文献   

9.
The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.  相似文献   

10.
Crosshole ground penetrating radar (GPR) tomography has been widely used and has the potential to improve the obtained subsurface models due to its high spatial resolution compared to other methods. Recent advances in full-waveform inversion of crosshole GPR data show that higher resolution images can be obtained compared to conventional ray-based GPR inversion because it can exploit all information present in the observed data. Since the first application of full-waveform inversion on synthetic and experimental GPR data, the algorithm has been significantly improved by extending the scalar to a vectorial approach, and changing the stepped permittivity and conductivity update into a simultaneous update. Here, we introduce new normalized gradients that do not depend on the number of sources and receivers which enable a comparison of the gradients and step lengths for different crosshole survey layouts. An experimental data set acquired at the Boise Hydrogeophysics Research Site is inverted using different source–receiver setups and the obtained permittivity and conductivity images, remaining gradients and final misfits are compared for the different versions of the full-waveform inversion. Moreover, different versions of the full-waveform inversion are applied to obtain an overview of all improvements. Most improvements result in a reducing final misfit between the measured and synthetic data and a reducing remaining gradient at the final iteration. Regions with relatively high remaining gradient amplitudes indicate less reliable inversion results. Comparison of the final full-waveform inversion results with Neutron–Neutron porosity log data and capacitive resistivity log data show considerably higher spatial frequencies for the logging data compared to the full-waveform inversion results. To enable a better comparison, we estimated a simple wavenumber filter and the full-waveform inversion results show an improved fit with the logging data. This work shows the potential of full-waveform inversion as an advanced method that can provide high resolution images to improve hydrological models.  相似文献   

11.
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 m2), 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.  相似文献   

12.
The results of geophysical studies conducted with selected electrical and electromagnetic methods in the Kondratowa Valley in the Tatra Mountains (the Carpathian Mountains, Poland) are presented in the article. The surveys were performed with the following methods: electrical resistivity tomography (ERT), georadar (GPR) and conductivity meter (CM). The objective of the noninvasive geophysical measurements was to determine the thickness of the Quaternary postglacial sediments that fill the bottom of the valley and to designate the accumulation of boulders deposited on Quaternary sediments. The results of ERT surveys conducted along the axis of the valley allowed to determine the changeability of the thickness of the postglacial sediments and allowed to designate a few areas of occurrence of boulders. The ERT, GPR and CM surveys conducted across the valley allowed to designate with high accuracy the thickness of the accumulation of boulders sliding down the valley bottom from the couloirs surrounding the valley.  相似文献   

13.
Geodynamic activity in the area of West Bohemia is typified by the occurrence of earthquake swarms, Quaternary volcanism and high flux of mantle-derived CO2. The highest swarm activity occurs beneath the eastern edge of the Cheb basin, which is delineated by the NW-SE trending morphologically pronounced Mariánské Lázn?? Fault (MLF) controlling the formation of the basin. The previous trenching survey across the MLF zone has identified several fault strands with possible Quaternary activity. In this paper we present the results of the geophysical survey focused to trace the faults signatures in geophysical sections and to build an image of near surface tectonics. The method of electric resistivity tomography (ERT) along two profiles parallel to the trench identified a strong resistivity contrast between the bodies of sandy gravels in the middle and conductive clayey sands to the west and weathered crystalline basement to the east. The 2-D ground penetration radar (GPR) sections show direct correlation of reflections with lithological boundaries identified in the trench. As expected, the GPR signal amplitudes increase with the resistivities found in the ERT sections. Two of the four faults identified in the trench are indicated in the resistivity and GPR sections. A 3-D GPR measurement has identified a spot of high amplitudes elongated parallel to the MLF trend, which coincides with the high resistivity body. To improve the signal-to-noise ratio of the time slices we stacked the GPR time slices within vertically homogeneous blocks. This provided a contrast image of the sand-gravel body including its boundaries in three dimensions. The detailed analysis of the 3-D GPR cube revealed additional fault that limits the highly reflective sands and appears to be offset by another younger fault. Our results suggest a complex fault pattern in the studied area, which deserves a further study.  相似文献   

14.
The Carlo's V Castle, located in Crotone Town, on the Ionian coast of the Calabria Region (Italy), date back to the 13th century d.C. (Fig. 1). During its long life, the building changed several owners and sustained the damages and the consequent reconstructions due to the innumerable naval battles. Moreover, the castle suffered the action of the earthquakes which always afflict the region.With the principal aim of detecting the location, depth and geometry of the rests of destroyed structures, a systematic Ground Penetrating Radar (GPR) survey was carried out in the area inside the boundary walls. The results are sixty-two one-meter-spaced, filtered and migrated radargrams arranged in four 3D data-sets. From each data-set, the most significant time-slice was extracted.To reduce the ambiguity in the GPR data interpretation, additional geophysical techniques, such as Magnetic (M), and Electrical Resistivity Tomography (ERT), were carried out with a partial superimposition with the GPR data. A comparison and a joint interpretation amongst different geophysical data pointed out some very remarkable features associated to buried remains and possible buried cannonballs.With the secondary aim to check the presence of an old military walkway linking two bastions a GPR profile was carried out on the sea side boundary wall. The GPR results are in agreement with an ERT survey carried out on the same profile and consistent with the presence of an underground passage.  相似文献   

15.
A geophysical study that involved different techniques was carried out with the aim to improve the knowledge of the archaeological site where the Basilica of Maxentius was founded and to discern individual covered structures (foundations).Vertical Electrical Sounding (VES), seismic refraction and Electrical Resistivity Tomography (ERT) studies were performed in the archaeological site. VES and seismic refraction allowed to characterise the main geological formations of the hill where the Basilica was built and to distinguish the concrete floor and backfilling. Electrical data were processed using different algorithms; their results were compared to appraise the inverted models' robustness.ERT inversion algorithms were used to delineate shape and size of a much more complex structure, that were originally expected from archaeological excavation plan. The results of the commercial program were used as a posteriori information to include them in the algorithm proposed by the authors; the sequential use of the programs defined a processing procedure.The integrated use of different geophysical techniques reduced a great deal the intrinsic ambiguities of each method. Direct explorations (boreholes and archaeological excavations) confirmed the geophysical results.  相似文献   

16.
Geotechnical projects usually rely on traditional sounding and drilling investigations. Drilling only provides point information and the geology needs to be interpolated between these points. Near surface geophysical methods can provide information to fill those gaps. Norwegian case studies are presented to illustrate how two-dimensional electrical resistivity tomography (ERT) can be used to accurately map the extent of quick clay deposits. Quick clay may be described as highly sensitive marine clay that changes from a relatively stiff condition to a liquid mass when disturbed. Quick clay slides present a geo-hazard and therefore layers of sensitive clay need to be mapped in detail. They are usually characterized by higher resistivity than non-sensitive clay and ERT is therefore a suitable approach to identify their occurrence. However, our experience shows that ERT cannot resolve this small resistivity contrast near large anomalies such as a bedrock interface. For this reason, a constrained inversion of ERT data was applied to delineate quick clay extent both vertically and laterally. As compared to the conventional unconstrained inversions, the constrained inversion models exhibit sharper resistivity contrasts and their resistivity values agree better with in situ measurements.  相似文献   

17.
The objective of data assimilation is to provide physically consistent estimates of spatially distributed environmental variables. In this study a relatively simple data assimilation method has been implemented in a relatively complex hydrological model. The data assimilation technique is Newtonian relaxation or nudging, in which model variables are driven towards observations by a forcing term added to the model equations. The forcing term is proportional to the difference between simulation and observation (relaxation component) and contains four-dimensional weighting functions that can incorporate prior knowledge about the spatial and temporal variability and characteristic scales of the state variable(s) being assimilated. The numerical model couples a three-dimensional finite element Richards equation solver for variably saturated porous media and a finite difference diffusion wave approximation based on digital elevation data for surface water dynamics. We describe the implementation of the data assimilation algorithm for the coupled model and report on the numerical and hydrological performance of the resulting assimilation scheme. Nudging is shown to be successful in improving the hydrological simulation results, and it introduces little computational cost, in terms of CPU and other numerical aspects of the model’s behavior, in some cases even improving numerical performance compared to model runs without nudging. We also examine the sensitivity of the model to nudging term parameters including the spatio-temporal influence coefficients in the weighting functions. Overall the nudging algorithm is quite flexible, for instance in dealing with concurrent observation datasets, gridded or scattered data, and different state variables, and the implementation presented here can be readily extended to any of these features not already incorporated. Moreover the nudging code and tests can serve as a basis for implementation of more sophisticated data assimilation techniques in a Richards equation-based hydrological model.  相似文献   

18.
The history of subsidence, fires, flooding and other kinds of environmental hazards related to shallow coal workings in India goes back to colonial times some 300 years ago. As coal production accelerated in modern times, so did the environmental and socio-economic drawbacks related to exploitation. In the mid-1980s, a hydropneumatic sand-stowing method was developed to fill in abandoned galleries but their exact location had to be known. Unfortunately, most of these old workings are uncharted and consequently large tracts of land cannot be stabilized. A research programme making use of integrated surface, borehole and cross-hole geophysical methods was undertaken over a five-year span to try to solve this problem. Surface geophysical methods, being cheaper and faster than their cross- and downhole counterparts, were used to cover larger areas on an exploratory basis, while cross-hole methods were employed to locate more accurately one or a network of galleries to be perforated by drillhole(s) and used as a conduit for sand stowing. We report the results of one of the cross-hole geophysical methods: electrical resistivity tomography (ERT). A pole–dipole configuration is used and both cross-hole and surface-borehole methodologies are tested. Forward modelling and inversion of synthetic data making use of downhole and surface physical and geometrical parameters are presented first. This phase is followed by the inversion of real data. It is concluded that ERT is not applicable for the detection of dry voids, but is effective in a waterlogged environment which is estimated to represent 85–90% of the cases. In waterlogged galleries, ERT is applicable in both cross-hole and surface-downhole modes, the latter allowing a larger surface coverage at low cost. ERT is thus a reliable geophysical tool to image water-filled voids and an adequate technique to address environmental and geotechnical problems.  相似文献   

19.
The Peinan archaeological site is the most intact Neolithic village with slate coffin burial complexes in Taiwan. However, the area that potentially contains significant ancient remains is covered by dense vegetation. No reliable data show the distribution of the ancient village, and no geophysical investigation has been performed at this site. To evaluate various geophysical methods under the geological setting and surface condition of the site, the physical properties of the remains were measured and four geophysical methods involving magnetic, electromagnetic (EM), electrical resistivity tomography (ERT), and ground-penetrating radar (GPR) were tested along three parallel profiles. The results imply that the EM and magnetic methods are much cost-effective and suitable for investigating the entire area. GPR and ERT methods can provide high resolution subsurface image, which are much suitable for subsequently detail investigation.The EM and magnetic surveys were thus conducted over the entire Peinan Cultural Park to understand the distribution of the ancient building remains at the Peinan site. The results of this study were verified by subsequent excavations, which indicate that the EM survey was successful in delineating the majority of the ancient village because the basements of building are highly resistive in comparison to the background sediment. The results of this investigation suggest that the ancient village was broadly distributed over the eastern part of the Peinan Culture Park and extended to the southeast.  相似文献   

20.
In recent years, ground‐penetrating radar (GPR) has been increasingly used for characterization of subglacial and englacial environments at polythermal glaciers. The geophysical method is able to exploit the dielectric difference between water, air, sediment and ice, allowing delineation of subsurface hydrological, thermal and structural conditions. More recent GPR research has endeavoured to examine temporal change in glaciers, in particular the distribution of the cold ice zone at polythermal glaciers. However, the exact nature of temporal change that can be identified using GPR has not been fully examined. This research presents the results of three GPR surveys conducted over the course of a summer ablation season at a polythermal glacier in the Canadian Arctic. A total of approximately 30 km of GPR profiles were collected in 2002 repeatedly covering the lower 2 km of Stagnation Glacier, Bylot Island (72°58′ N 78°22′ W). Comparison between profiles indicated changes in the radar signature, including increased noise, appearance and disappearance of englacial reflections, and signal attenuation in the latter survey. Further, an area of chaotic returns in up‐glacier locations, which was interpreted to be a wet temperate ice zone, showed marked recession over the course of the ablation season. Combining all the temporal changes that were detected by GPR, results indicate that a polythermal glacier may exhibit strongly seasonal changes in hydrological and thermal characteristics throughout the ice body, including the drainage of 17 000 m3 of temporarily stored intra‐glacial meltwater. It is also proposed that the liquid water content in the temperate ice zone of polythermal glaciers can be described as a fraction of a specific retention capacity. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号