Evaluating soil moisture estimation from ground-penetrating radar hyperbola fitting with respect to a systematic time-domain reflectometry data collection in a boreal podzolic agricultural field

The upper 30 cm of the soil profile, which hosts the majority of the root biomass, can be considered as the shallow agricultural root zone of most temperate crops. The electromagnetic wave velocity in the soil obtained from reflection hyperbolas in ground-penetrating radar (GPR) data can be used to estimate soil moisture (SM). Finding shallow hyperbolas in a radargram and minimizing the subjective error associated with the hyperbola fitting are the main challenges in this approach. Nevertheless, we were motivated by the recent improvements of hyperbola fitting algorithms, which can reduce the subjective error and processing time. To overcome the difficulty of finding very shallow hyperbolas, we applied the hyperbola fitting method to reflections ranging from 27- to 50-cm depth using a 500-MHz centre-frequency GPR and compared the estimated moisture with vertically installed, 30-cm-long time-domain reflectometry (TDR) sensors. We also compared TDR and GPR sample areas in a 2-D plane using different GPR survey types and different hyperbola depths. SM measured with TDR and GPR were not significantly different according to Mann–Whitney's test. Our analyses showed that a root mean square error of 0.03 m³ m⁻³ was found between the two methods. In conclusion, the proposed method might be suitable to estimate SM with an acceptable accuracy within the root zone if the soil profile is fairly uniform within the application depth range.     

Determining the number of manual measurements required to improve peat thickness estimations by ground penetrating radar

Ground penetrating radar (GPR) is often used for investigating peat thickness. The quality of GPR measurements depends on electromagnetic wave (EMW) velocity estimates. The objective of this study is to determine the number of manual measurements required to minimize EMW velocity error in peatlands. A total of 175 manually measured peat thicknesses are used with a depth‐to‐target method to assess EMW velocity in two southern Quebec peatlands. Mean measured EMW velocities are 0·040 and 0·039 m ns^–1 with standard deviations of 0·013 and 0·008 m ns^–1. Statistical analyses show that at least 30 calibration points are required to minimize the EMW velocity error, regardless of the geological setting. Copyright © 2008 John Wiley & Sons, Ltd.     

Determination of electromagnetic wave velocity in horizontally layered sedimentary target: A ground-penetrating radar study from Silurian limestones,Estonia

Ground-penetrating radar (GPR) is a non-destructive geophysical technique to obtain information about shallow subsurface by transmitting electromagnetic waves into the ground and registering signals reflected from objects or layers with different dielectric properties. The present GPR study was conducted in Võhmuta limestone quarry in Estonia in order to describe the relationship between GRP responses to the variations in petrophysical properties. Sub-horizontally oriented cores for petrophysical measurements were drilled from the side wall of the quarry. The GPR profiles were run at the sloped trench floor and on the top of side wall in order to correlate traceable reflections with physical properties. Based on three techniques: (i) hyperbola fitting, (ii) wide angle reflection and refraction (WARR), and (iii) topographic, a mean electromagnetic wave velocity value of 9.25 cm ns^?1 (corresponding to relative dielectric permittivity of 10.5) was found to describe the sequence and was used for time-to-depth conversion. Examination of radar images against petrophysical properties revealed that major reflections appear in levels where the changes in porosity occur.     

Measuring snow ablation rates in alpine terrain with a mobile multioffset ground‐penetrating radar system

Ground‐penetrating radar (GPR) has become a promising technique in the field of snow hydrological research. It is commonly used to measure snow depth, density, and water equivalent over large distances or along gridded snow courses. Having built and tested a mobile lightweight set‐up, we demonstrate that GPR is capable of accurately measuring snow ablation rates in complex alpine terrain. Our set‐up was optimized for efficient measurements and consisted of a multioffset radar with four pairs of antennas mounted to a plastic sled, which was small enough to permit safe and convenient operations. Repeated measurements at intervals of 2 to 7 days were taken during the 2014/2015 winter season along 10 profiles of 50 to 200 m length within two valleys located in the eastern Swiss Alps. Resulting GPR‐based data of snow depth, density, and water equivalent, as well as their respective change over time, were in good agreement with concurrent manual measurements, in particular if accurate alignment between repeated overpasses could be achieved. Corresponding root‐mean‐square error (RMSE) values amounted to 4.2 cm for snow depth, 17 mm for snow water equivalent, and 22 kg/m³ for snow density, with similar RMSE values for corresponding differential data. With this performance, the presented radar set‐up has the potential to provide exciting new and extensive datasets to validate snowmelt models or to complement lidar‐based snow surveys.     

Acoustic and electromagnetic properties of soils saturated with salt water and NAPL

Electrical, seismic, and electromagnetic methods can be used for noninvasive determination of subsurface physical and chemical properties. In particular, we consider the evaluation of water salinity and the detection of surface contaminants. Most of the relevant properties are represented by electric conductivity, P-wave velocity, and dielectric permittivity. Hence, it is important to obtain relationships between these measurable physical quantities and soil composition, saturation, and frequency. Conductivity in the geoelectric frequency range is obtained with Pride's model for a porous rock. (The model considers salinity and permeability.) White's model of patchy saturation is used to calculate the P-wave velocity and attenuation. Four cases are considered: light nonaqueous phase liquid (LNAPL) pockets in water, dense nonaqueous phase liquid (DNAPL) pockets in water, LNAPL pockets in air, and DNAPL pockets in air. The size of the pockets (or pools), with respect to the signal wavelength, is modeled by the theory. The electromagnetic properties in the GPR frequency range are obtained by using the Hanai–Bruggeman equation for two solids (sand and clay grains) and two fluids (LNAPL or DNAPL in water or air). The Hanai–Bruggeman exponent (1/3 for spherical particles) is used as a fitting parameter and evaluated for a sand/clay mixture saturated with water.Pride's model predicts increasing conductivity for increasing salinity and decreasing permeability. The best-fit exponent of the Hanai–Bruggeman equation for a sand/clay mixture saturated with water is 0.61, indicating that the shape of the grains has a significant influence on the electromagnetic properties. At radar frequencies, it is possible to distinguish between a water-saturated medium and a NAPL-saturated medium, but LNAPL- and DNAPL-saturated media have very similar electromagnetic properties. The type of contaminant can be better distinguished from the acoustic properties. P-wave velocity increases with frequency, and has dissimilar behaviour for wet and dry soils.     

Modelling of GPR waves for lossy media obeying a complex power law of frequency for dielectric permittivity

The attenuation of ground‐penetrating radar (GPR) energy in the subsurface decreases and shifts the amplitude spectrum of the radar pulse to lower frequencies (absorption) with increasing traveltime and causes also a distortion of wavelet phase (dispersion). The attenuation is often expressed by the quality factor Q. For GPR studies, Q can be estimated from the ratio of the real part to the imaginary part of the dielectric permittivity. We consider a complex power function of frequency for the dielectric permittivity, and show that this dielectric response corresponds to a frequency‐independent‐Q or simply a constant‐Q model. The phase velocity (dispersion relationship) and the absorption coefficient of electromagnetic waves also obey a frequency power law. This approach is easy to use in the frequency domain and the wave propagation can be described by two parameters only, for example Q and the phase velocity at an arbitrary reference frequency. This simplicity makes it practical for any inversion technique. Furthermore, by using the Hilbert transform relating the velocity and the absorption coefficient (which obeys a frequency power law), we find the same dispersion relationship for the phase velocity. Both approaches are valid for a constant value of Q over a restricted frequency‐bandwidth, and are applicable in a material that is assumed to have no instantaneous dielectric response. Many GPR profiles acquired in a dry aeolian environment have shown a strong reflectivity inside dunes. Changes in water content are believed to be the origin of this reflectivity. We model the radar reflections from the bottom of a dry aeolian dune using the 1D wavelet modelling method. We discuss the choice of the reference wavelet in this modelling approach. A trial‐and‐error match of modelled and observed data was performed to estimate the optimum set of parameters characterizing the materials composing the site. Additionally, by combining the complex refractive index method (CRIM) and/or Topp equations for the bulk permittivity (dielectric constant) of moist sandy soils with a frequency power law for the dielectric response, we introduce them into the expression for the reflection coefficient. Using this method, we can estimate the water content and explain its effect on the reflection coefficient and on wavelet modelling.     

中尺度模拟土壤油污染区探地雷达探测效果分析

为使物理模拟实验效果与实际探测情形更为接近,在室外自然条件下建立中尺度土壤石油污染实验模型,油污染区扩展深度超过1 m,采用实地探测中常用的500 MHz雷达天线进行长期定时探测.通过实测雷达图像特征、土壤含水量含油量分析,并对比前人开展的小尺度室内模拟试验结果,综合评价探地雷达对油污染区的探测效果.研究表明探地雷达探测图像异常特征与污染区扩散阶段密切相关:包气带内油污染区会引起振幅增强;毛细带的油污染区则表现为水位面反射轴附近清晰可辨的高幅异常区,且水位面反射轴呈下凹状;随扩散过程持续进行,异常区下移与水位面反射轴相交,并产生水平扩张.当污染土含油饱和度大于20%时,可通过雷达图像异常区圈定污染范围;当污染土含油饱和度大于15%时,可通过频谱图出现低频响应的位置圈定污染区水平范围.中尺度实验结果与室内小尺度模拟结果具有一致性,可作为油污染区雷达图像异常的解译依据.     

Seasonal changes in ground‐penetrating radar signature observed at a polythermal glacier,Bylot Island,Canada

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 m³ 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.     

A conceptual model of supra‐glacial lake formation on debris‐covered glaciers based on GPR facies analysis

Supra‐glacial lakes and ponds can create hotspots of mass loss on debris‐covered glaciers. While much research has been directed at understanding lateral lake expansion, little is known about the rates or processes governing lake deepening. To a large degree, this knowledge gap persists due to sparse observations of lake beds. Here we report on the novel use of ground penetrating radar (GPR) surveys to simultaneously collect supra‐glacial lake bathymetry and bottom composition data from Spillway Lake (surface area of 2.4 × 10⁵ m²; volume of 9.5 × 10⁴ m³), which is located in the terminus region of the Ngozumpa Glacier in the Khumbu region of the Nepal Himalaya. We identified two GPR bottom signals corresponding to two sedimentary facies of (1) sub‐horizontal layered fine sediment drape and (2) coarse blocky diamict. We provide an understanding of the changes in subaqueous debris distribution that occur through stages of lake expansion by combining the GPR results with in situ observations of shoreline deposits matching the interpreted facies. From this, we present an updated conceptual model of supra‐glacial lake evolution, with the addition of data on the evolving debris environment, showing how dominant depositional processes can change as lakes evolve from perched lakes to multi‐basin base‐level lakes and finally onto large moraine‐dammed lakes. Throughout lake evolution, processes such as shoreline steepening, lakebed collapse into voids and conduit interception, subaerial and subaqueous calving and rapid areal expansion alter the spatial distribution and makeup of lakebed debris and sediments forcing a number of positive and negative feedbacks on lake expansion. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of overlying velocity,particle size,and biofilm growth on stream–subsurface exchange of particles

Stream–subsurface exchange strongly influences the transport of contaminants, fine particles, and other ecologically relevant substances in streams. We used a recirculating laboratory flume (220 cm long and 20 cm wide) to study the effects of particle size, overlying velocity, and biofilm formation on stream–subsurface exchange of particles. Sodium chloride was used as a non‐reactive dissolved tracer and 1‐ and 5‐µm fluorescent microspheres were used as particulate tracers. Surface–subsurface exchange was observed with a clean sand bed and a bed colonized by an autotrophic–heterotrophic biofilm under two different overlying velocities, 0·9 and 5 cm s^?1. Hydrodynamic interactions between the overlying flow and sand bed resulted in a reduction of solute and particle concentrations in the water column, and a corresponding accumulation of particles in both the sediments and in the biofilm. Increasing overlying velocity and particle size resulted in faster removal from the overlying water due to enhanced mass transfer to the bed. The presence of the biofilm did not affect solute exchange under any flow condition tested. The presence of the biofilm significantly increased the deposition of particles under an overlying velocity of 5 cm s^?1, and produced a small but statistically insignificant increase at 0·9 cm^?1. The particles preferentially deposited within the biofilm matrix relative to the underlying sand. These results demonstrate that hydrodynamic transport conditions, particle size, and biofilm formation play a key role in the transport of suspended particles, such as inorganic sediments, particulate organic matter, and pathogenic microorganisms in freshwater ecosystems, and should be taken into consideration when predicting the fate and transport of particles and contaminants in the environment. Copyright © 2009 John Wiley & Sons, Ltd.     

Hydrological connectivity of soil pipes determined by ground‐penetrating radar tracer detection

Soil pipes are common and important features of many catchments, particularly in semi‐arid and humid areas, and can contribute a large proportion of runoff to river systems. They may also signi?cantly in?uence catchment sediment and solute yield. However, there are often problems in ?nding and de?ning soil pipe networks which are located deep below the surface. Ground‐penetrating radar (GPR) has been used for non‐destructive identi?cation and mapping of soil pipes in blanket peat catchments. While GPR can identify subsurface cavities, it cannot alone determine hydrological connectivity between one cavity and another. This paper presents results from an experiment to test the ability of GPR to establish hydrological connectivity between pipes through use of a tracer solution. Sodium chloride was injected into pipe cavities previously detected by the radar. The GPR was placed downslope of the injection points and positioned on the ground directly above detected soil pipes. The resultant radargrams showed signi?cant changes in re?ectance from some cavities and no change from others. Pipe waters were sampled in order to check the radar results. Changes in electrical conductivity of the pipe water could be detected by the GPR, without data post‐processing, when background levels were increased by more than approximately twofold. It was thus possible to rapidly determine hydrological connectivity of soil pipes within dense pipe networks across hillslopes without ground disturbance. It was also possible to remotely measure travel times through pipe systems; the passing of the salt wave below the GPR produced an easily detectable signal on the radargram which required no post‐processing. The technique should allow remote sensing of water sources and sinks for soil pipes below the surface. The improved understanding of ?owpath connectivity will be important for understanding water delivery, solutional and particulate denudation, and hydrological and geomorphological model development. Copyright © 2004 John Wiley & Sons, Ltd.     

GPR reflection characteristics and depositional models of mud volcanic sediments—Wushanting mud volcano field, southwestern Taiwan

Ground penetrating radar (GPR) survey was conducted in the Wushanting mud volcano field (Yanchao, Kaohsiung) using a 500 MHz antennae, which allowed high-resolution imaging of subsurface structures. Seven GPR reflection characteristics are recognized. Sigmoid GPR reflection patterns resulted from a recent mud lobe deposited above an underlying older mud lobe front. Contorted GPR facies resulted from recent mud flow which encountered obstacles. Subparallel reflections resulted from mud volcano deposits of limited flowability, low velocity and gentle gradient. Hummocky reflection patterns are formed by interfingering of recent mud lobes building onto low land. Disrupted GPR facies were due to lateral breaks of continuity from mud cracks, which, according to field observation, can provide channels for erosion and form deeper erosion gullies. GPR time slices of different depths are rendered as a three-dimensional model. Approximately orbicular GPR reflection characteristics can indicate arcuate stacked mud lobe fronts of different periods. Some depositional models to explain GPR reflection characteristics can be founded upon observations of recent sedimentary phenomena. The models of this study may be applied to paleoenvironments and the depositional evolution of mud volcanoes in similar geological settings.     

A large landslide event in a post‐glacial landscape: rethinking glacial legacy

Threlkeld Knotts (c. 500 m above sea level) in the English Lake District has hitherto been considered to be a glacially‐modified intrusion of microgranite. However, its surface features are incompatible with glacial modification; neither can these nor the subsurface structures revealed by ground‐penetrating radar (GPR) be explained by post‐glacial subaerial processes acting on a glacially‐modified microgranite intrusion. Here we re‐interpret Threlkeld Knotts as a very large post‐glacial landslide involving the microgranite, with an estimated volume of about 4 × 10⁷ m³. This interpretation is tested against published and recent information on the geology of the site, the glacial geomorphic history of the area and newly‐acquired GPR data. More than 60 large post‐Last Glacial Maximum (LGM) rock–slope failures have significantly modified the glaciated landscape of the Lake District; this is one of the largest. Recognition of this major landslide deposit in such a well‐studied environment highlights the need to continuously re‐examine landscapes in the light of increasing knowledge of geomorphic processes and with available technology in currently active or de‐glaciating environments. Copyright © 2012 John Wiley & Sons, Ltd.     

Active microwave remote sensing for soil moisture measurement: a field evaluation using ERS‐2

Active microwave remote sensing observations of backscattering, such as C‐band vertically polarized synthetic aperture radar (SAR) observations from the second European remote sensing (ERS‐2) satellite, have the potential to measure moisture content in a near‐surface layer of soil. However, SAR backscattering observations are highly dependent on topography, soil texture, surface roughness and soil moisture, meaning that soil moisture inversion from single frequency and polarization SAR observations is difficult. In this paper, the potential for measuring near‐surface soil moisture with the ERS‐2 satellite is explored by comparing model estimates of backscattering with ERS‐2 SAR observations. This comparison was made for two ERS‐2 overpasses coincident with near‐surface soil moisture measurements in a 6 ha catchment using 15‐cm time domain reflectometry probes on a 20 m grid. In addition, 1‐cm soil moisture data were obtained from a calibrated soil moisture model. Using state‐of‐the‐art theoretical, semi‐empirical and empirical backscattering models, it was found that using measured soil moisture and roughness data there were root mean square (RMS) errors from 3·5 to 8·5 dB and r² values from 0·00 to 0·25, depending on the backscattering model and degree of filtering. Using model soil moisture in place of measured soil moisture reduced RMS errors slightly (0·5 to 2 dB) but did not improve r² values. Likewise, using the first day of ERS‐2 backscattering and soil moisture data to solve for RMS surface roughness reduced RMS errors in backscattering for the second day to between 0·9 and 2·8 dB, but did not improve r² values. Moreover, RMS differences were as large as 3·7 dB and r² values as low as 0·53 between the various backscattering models, even when using the same data as input. These results suggest that more research is required to improve the agreement between backscattering models, and that ERS‐2 SAR data may be useful for estimating fields‐scale average soil moisture but not variations at the hillslope scale. Copyright © 2004 John Wiley & Sons, Ltd.     

地质雷达在乌拉山山前断裂探测中的应用

地质雷达是利用电磁波对地下不同电性介质进行探测的地球物理仪器，其探测速率快、分辨率高，可弥补探槽和其他地球物理方法存在探测盲区的缺陷，正在越来越多地应用于活动断层探测领域。本文以乌拉山山前断裂为例开展地质雷达探测工作，使用无人机正射影像技术对测线进行地形校正，获得断层浅部地质雷达图像。研究结果表明，本文研究方法能有效反映探槽揭露的地层单元和断层分布。本次探测中，雷达波形图像特征为:浅地表的土壤层反射波总体较弱；粗粒沉积为主的砾石层反射波总体较强，同相轴连续性好；细粒沉积为主的砂层反射波弱于砾石层，波形以中、高频为主，同相轴具有弱连续性；对于洪冲积地区，地质雷达能分辨具有一定特征的地层单元，这为剖面图像的断层识别提供了标志；通过无人机正射影像技术对地质雷达测线进行地形校正，有利于获得更为准确的探测结果。     

Evaluating some factors that affect feasility of using ground penetrating radar for landmine detection

Ground penetrating radar （GPR） is one of the promising technologies that can be used to detect landmines. Many factors may affect the ability of GPR to detect landmines. Among those factors are： 1） the type of landmine material （metallic or plastic）, 2） conditions of the host soil （soil texture and soil moisture）, and 3） the radar frequency utilized. The impact of these factors on the ability of GPR to detect landmines is investigated by studying their effect on the dielectric permittivity contrast between the landmine and the host soil, as well as on the attenuation of the radar waves. The impact of each factor was theoretically reviewed and modeled using the Matlab and Mathcad software packages. Results of the computer modeling were correlated with GPR data acquired for metallic and plastic landmine types. It was found that the ability of GPR to detect landmines depends to a great extent on the landmine type, water content of the host soil, utilized radar frequency, and soil texture. The landmines are much easier to detect than plastic landmines for any soil conditions and any radar frequency. Increasing the soil＇s moisture content, regardless of soil texture, eases the detection of the plastic landmine and worsens the detection of the metallic mines. Increasing the percentage of clay in the soil causes the same effect as the moisture content. However, higher radar frequency delivers better results for landmine detection as long as the percentage of clay and the moisture content in the soil remains low. The results of this study are expected to help in selecting optimum radar antennae and data acquisition parameters depending on the landmine type and environmental conditions.     

GPR reflection characteristics and depositional models of mud volcanic sediments—Wushanting mud volcano field,southwestern Taiwan

Ground penetrating radar (GPR) survey was conducted in the Wushanting mud volcano field (Yanchao, Kaohsiung) using a 500 MHz antennae, which allowed high-resolution imaging of subsurface structures. Seven GPR reflection characteristics are recognized. Sigmoid GPR reflection patterns resulted from a recent mud lobe deposited above an underlying older mud lobe front. Contorted GPR facies resulted from recent mud flow which encountered obstacles. Subparallel reflections resulted from mud volcano deposits of limited flowability, low velocity and gentle gradient. Hummocky reflection patterns are formed by interfingering of recent mud lobes building onto low land. Disrupted GPR facies were due to lateral breaks of continuity from mud cracks, which, according to field observation, can provide channels for erosion and form deeper erosion gullies. GPR time slices of different depths are rendered as a three-dimensional model. Approximately orbicular GPR reflection characteristics can indicate arcuate stacked mud lobe fronts of different periods. Some depositional models to explain GPR reflection characteristics can be founded upon observations of recent sedimentary phenomena. The models of this study may be applied to paleoenvironments and the depositional evolution of mud volcanoes in similar geological settings.     

基于全波形反演的探地雷达数据逆时偏移成像

逆时偏移成像(RTM)常用来处理复杂速度模型,包括陡倾角及横向速度变化剧烈的模型.与常规偏移成像方法(如Kirchhoff偏移)相比,逆时偏移成像能提供更好的偏移成像结果,近些年逆时偏移成像越来越广泛地应用到勘探地震中,它逐渐成为石油地震勘探中的一种行业标准.电磁波和弹性波在动力学和运动学上存在相似性,故本文开发了基于麦克斯韦方程组的电磁波逆时偏移成像算法,并将其应用到探地雷达数据处理中.时间域有限差分(FDTD)用于模拟电磁波正向和逆向传播过程,互相关成像条件用于获得最终偏移结果.逆时偏移成像算法中,偏移成像结果受初始模型影响较大,而其中决定电磁波传播速度的介电常数的影响尤为重要.本文基于时间域全波形反演(FWI)算法反演获得了更为精确的地下介电常数模型,并将其反演结果作为逆时偏移成像的初始介电常数模型.为了验证此算法的有效性,首先构建了一个复杂地质结构模型,合成了共偏移距及共炮点探地雷达数据,分别应用常规Kirchhoff偏移算法及逆时偏移成像算法进行偏移处理,成像结果显示由逆时偏移成像算法得到的偏移结果与实际模型具有较高的一致性;此外本文在室内沙槽中进行了相关的物理模拟实验,采集了共偏移距及共炮点探地雷达数据,分别应用Kirchhoff和叠前逆时偏移成像算法进行处理,结果表明叠前逆时偏移成像在实际应用中能获得更好的成像效果.     

Measuring thaw depth beneath peat‐lined arctic streams using ground‐penetrating radar

In arctic streams, depth of thaw beneath the stream channel is likely a significant parameter controlling hyporheic zone hydrology and biogeochemical cycling. As part of an interdisciplinary study of this system, we conducted a field investigation to test the effectiveness of imaging substream permafrost using ground‐penetrating radar (GPR). We investigated three sites characterized by low‐energy water flow, organic material lining the streambeds, and water depths ranging from 0·2 to 2 m. We acquired data using a 200 MHz pulsed radar system with the antennas mounted in the bottom of a small rubber boat that was pulled across the stream while triggering the radar at a constant rate. We achieved excellent results at all three sites, with a clear continuous image of the permafrost boundary both peripheral to and beneath the stream. Our results demonstrate that GPR can be an effective tool for measuring substream thaw depth. Copyright © 2005 John Wiley & Sons, Ltd.