RADAR DETECTION OF BOREHOLES IN ADVANCE OF MINING*

Boreholes present a potential hazard to mining in rock salt mines. The only rock salt mine in England was almost lost when they mined into an unknown borehole open to an aquifer. Rock salt is a good transmitter of VHF radar waves, so a good method of detecting boreholes is by using radar probing. This paper investigates theoretically the detection of a borehole in advance of mining using two different radar systems, Bravo II and Charlie II. Bravo II has a peak power of 20 kW at 230 MHz and Charlie II has 3 W at 440 MHz. The borehole detectability depends critically on the radar wave attenuation in salt, which is characterized by its loss tangent, tan δ. The backscattering cross section σ of the borehole is also important. These parameters are related in the radar signal-to-noise equation which determines the signal-to-noise ratio (S/N) for a given range in salt to the borehole. Theoretically, in salt with a tan δ of 10^-3, Charlie II can detect (S/N= 3 dB) a 12.7 cm borehole 225 m in advance of mining while Bravo II, using only 100 W, can detect the same borehole 680 m ahead. These values are with the antenna E field parallel to the axis of the borehole. Theory shows this antenna orientation to be optimum for maximum detection probability. For salt with a tan δ of 2 × 10^-5, the smallest value measured for salt, Bravo II can detect the borehole up to 15.3 km using maximum power of 20 kW.     

Renewed uplift at the Yellowstone Caldera measured by leveling surveys and satellite radar interferometry

 A first-order leveling survey across the northeast part of the Yellowstone caldera in September 1998 showed that the central caldera floor near Le Hardy Rapids rose 24±5 mm relative to the caldera rim at Lake Butte since the previous survey in September 1995. Annual surveys along the same traverse from 1985 to 1995 tracked progressive subsidence near Le Hardy Rapids at an average rate of –19±1 mm/year. Earlier, less frequent surveys measured net uplift in the same area during 1923–1976 (14±1 mm/year) and 1976–1984 (22±1 mm/year). The resumption of uplift following a decade of subsidence was first detected by satellite synthetic aperture radar interferometry, which revealed approximately 15 mm of uplift in the vicinity of Le Hardy Rapids from July 1995 to June 1997. Radar interferograms show that the center of subsidence shifted from the Sour Creek resurgent dome in the northeast part of the caldera during August 1992 to June 1993 to the Mallard Lake resurgent dome in the southwest part during June 1993 to August 1995. Uplift began at the Sour Creek dome during August 1995 to September 1996 and spread to the Mallard Lake dome by June 1997. The rapidity of these changes and the spatial pattern of surface deformation suggest that ground movements are caused at least in part by accumulation and migration of fluids in two sill-like bodies at 5–10 km depth, near the interface between Yellowstone's magmatic and deep hydrothermal systems. Received: 30 November 1998 / Accepted: 16 April 1999     

Radar rainfall uncertainty modelling influenced by wind

Radar‐based estimates of rainfall are affected by many sources of uncertainties, which would propagate through the hydrological model when radar rainfall estimates are used as input or initial conditions. An elegant solution to quantify these uncertainties is to model the empirical relationship between radar measurements and rain gauge observations (as the ‘ground reference’). However, most current studies only use a fixed and uniform model to represent the uncertainty of radar rainfall, without consideration of its variation under different synoptic regimes. Wind is such a typical weather factor, as it not only induces error in rain gauge measurements but also causes the raindrops observed by weather radar to drift when they reach the ground. For this reason, as a first attempt, this study introduces the wind field into the uncertainty model and designs the radar rainfall uncertainty model under different wind conditions. We separate the original dataset into three subsamples according to wind speed, which are named as WDI (0–2 m/s), WDII (2–4 m/s) and WDIII (>4 m/s). The multivariate distributed ensemble generator is introduced and established for each subsample. Thirty typical events (10 at each wind range) are selected to explore the behaviours of uncertainty under different wind ranges. In each time step, 500 ensemble members are generated, and the values of 5th to 95th percentile values are used to produce the uncertainty bands. Two basic features of uncertainty bands, namely dispersion and ensemble bias, increase significantly with the growth of wind speed, demonstrating that wind speed plays a considerable role in influencing the behaviour of the uncertainty band. On the basis of these pieces of evidence, we conclude that the radar rainfall uncertainty model established under different wind conditions should be more realistic in representing the radar rainfall uncertainty. This study is only a start in incorporating synoptic regimes into rainfall uncertainty analysis, and a great deal of more effort is still needed to build a realistic and comprehensive uncertainty model for radar rainfall data. Copyright © 2014 John Wiley & Sons, Ltd.     

Impact of Doppler Radar Wind in Simulating the Intensity and Propagation of Rainbands Associated with Mesoscale Convective Complexes Using MM5-3DVAR System

Pre-monsoon rainfall around Kolkata (northeastern part of India) is mostly of convective origin as 80% of the seasonal rainfall is produced by Mesoscale Convective Systems (MCS). Accurate prediction of the intensity and structure of these convective cloud clusters becomes challenging, mostly because the convective clouds within these clusters are short lived and the inaccuracy in the models initial state to represent the mesoscale details of the true atmospheric state. Besides the role in observing the internal structure of the precipitating systems, Doppler Weather Radar (DWR) provides an important data source for mesoscale and microscale weather analysis and forecasting. An attempt has been made to initialize the storm-scale numerical model using retrieved wind fields from single Doppler radar. In the present study, Doppler wind velocities from the Kolkata Doppler weather radar are assimilated into a mesoscale model, MM5 model using the three-dimensional variational data assimilation (3DVAR) system for the prediction of intense convective events that occurred during 0600 UTC on 5 May and 0000 UTC on 7 May, 2005. In order to evaluate the impact of the DWR wind data in simulating these severe storms, three experiments were carried out. The results show that assimilation of Doppler radar wind data has a positive impact on the prediction of intensity, organization and propagation of rain bands associated with these mesoscale convective systems. The assimilation system has to be modified further to incorporate the radar reflectivity data so that simulation of the microphysical and thermodynamic structure of these convective storms can be improved.     

Challenges and solutions of remote sensing at offshore wind energy developments

Radar is becoming an important tool used to gather data on bird and bat activity at proposed and existing land-based wind energy sites. Radar will likely play an even more important role at the increasing development of wind energy offshore, given both the lack of knowledge about bird and bat activity offshore and the increased difficultly in obtaining offshore information. Most radar studies to date have used off-the-shelf or modified marine radars. However, there are several issues that continue to hinder the potential usefulness of radar at wind energy sites, with offshore sites providing a particular suite of challenges. We identify these challenges along with current or developing solutions.     

The role of visual saliency in the automation of seismic interpretation

In this paper, we propose a workflow based on SalSi for the detection and delineation of geological structures such as salt domes. SalSi is a seismic attribute designed based on the modelling of human visual system that detects the salient features and captures the spatial correlation within seismic volumes for delineating seismic structures. Using this attribute we cannot only highlight the neighbouring regions of salt domes to assist a seismic interpreter but also delineate such structures using a region growing method and post‐processing. The proposed delineation workflow detects the salt‐dome boundary with very good precision and accuracy. Experimental results show the effectiveness of the proposed workflow on a real seismic dataset acquired from the North Sea, F3 block. For the subjective evaluation of the results of different salt‐dome delineation algorithms, we have used a reference salt‐dome boundary interpreted by a geophysicist. For the objective evaluation of results, we have used five different metrics based on pixels, shape, and curvedness to establish the effectiveness of the proposed workflow. The proposed workflow is not only fast but also yields better results as compared with other salt‐dome delineation algorithms and shows a promising potential in seismic interpretation.     

Inversion of swell frequency from a 1-year HF radar dataset collected in Brittany (France)

This article presents long period ocean wave (swell) frequencies inverted from a 13-month dataset of high-frequency (HF) phased array radars and an assessment of these estimates by comparison with WAVEWATCH III model data. The method of swell frequency inversion from high-frequency radar sea echo Doppler spectra is described. Radar data were collected from a two-site HF Wellen Radar (WERA) radar system on the west coast of Brittany (France) operating at 12 MHz. A standard beam-forming processing technique has been used to obtain Doppler spectra of processed radar cells. Swell frequencies are obtained from the frequencies of particular spectral peaks of the second-order continuum in hourly averaged Doppler spectra. The data coverage of effective Doppler spectra considered for swell frequency estimates shows the influence of islands and shallow water effects. Swell estimates from both radar stations are in good agreement. The comparison of radar-derived results to WAVEWATCH III (WW3) estimates shows that radar measurements agree quite well with model results. The bias and standard deviation between two estimates are very small for swells with frequency less than 0.09 Hz (period >11 s), whereas radar estimates are generally lower than model estimates for shorter swells, along with higher standard deviation. Statistical analysis suggests that radar measurement uncertainty explains most of the difference between radar and model estimates. For each swell event, time series of frequency exhibits a quasi-linear frequency increase which is associated with the dispersive property of wave phase velocity. The use of swell frequency estimates from both radars on common radar cells only slightly increases the accuracy of swell frequency measurement.     

Assessing Groundwater Depletion and Dynamics Using GRACE and InSAR: Potential and Limitations

In the last decade, remote sensing of the temporal variation of ground level and gravity has improved our understanding of groundwater dynamics and storage. Mass changes are measured by GRACE (Gravity Recovery and Climate Experiment) satellites, whereas ground deformation is measured by processing synthetic aperture radar satellites data using the InSAR (Interferometry of Synthetic Aperture Radar) techniques. Both methods are complementary and offer different sensitivities to aquifer system processes. GRACE is sensitive to mass changes over large spatial scales (more than 100,000 km²). As such, it fails in providing groundwater storage change estimates at local or regional scales relevant to most aquifer systems, and at which most groundwater management schemes are applied. However, InSAR measures ground displacement due to aquifer response to fluid‐pressure changes. InSAR applications to groundwater depletion assessments are limited to aquifer systems susceptible to measurable deformation. Furthermore, the inversion of InSAR‐derived displacement maps into volume of depleted groundwater storage (both reversible and largely irreversible) is confounded by vertical and horizontal variability of sediment compressibility. During the last decade, both techniques have shown increasing interest in the scientific community to complement available in situ observations where they are insufficient. In this review, we present the theoretical and conceptual bases of each method, and present idealized scenarios to highlight the potential benefits and challenges of combining these techniques to remotely assess groundwater storage changes and other aspects of the dynamics of aquifer systems.     

Plasma convection at high latitudes using the EISCAT VHF and ESR incoherent scatter radars

The recent availability of substantial data sets taken by the EISCAT Svalbard Radar allows several important tests to be made on the determination of convection patterns from incoherent scatter radar results. During one 30-h period, the Svalbard Radar made 15 min scans combining local field aligned observations with two, low elevation positions selected to intersect the two beams of the Common Programme Four experiment being simultaneously conducted by the EISCAT VHF radar at Troms. The common volume results from the two radars are compared. The plasma convection velocities determined independently by the two radars are shown to agree very closely and the combined three-dimensional velocity data used to test the common assumption of negligible field-aligned flow in this regime.     

Characterization of an unstable rock mass based on borehole logs and diverse borehole radar data

Unstable rocky slopes are major hazards to the growing number of people that live and travel though mountainous regions. To construct effective barriers to falling rock, it is necessary to know the positions, dimensions and shapes of structures along which failure may occur. To investigate an unstable mountain slope distinguished by numerous open fracture zones, we have taken advantage of three moderately deep (51.0–120.8 m) boreholes to acquire geophysical logs and record single-hole radar, vertical radar profiling (VRP) and crosshole radar data. We observed spallation zones, displacements and borehole radar velocity and amplitude anomalies at 16 of the 46 discontinuities identified in the borehole optical televiewer images. The results of the VRP and crosshole experiments were disappointing at our study site; the source of only one VRP reflection was determined and the crosshole velocity and amplitude tomograms were remarkably featureless. In contrast, much useful structural information was provided by the single-hole radar experiments. Radar reflections were recorded from many surface and borehole fracture zones, demonstrating that the strong electrical property contrasts of these features extended some distance into the adjacent rock mass. The single-hole radar data suggested possible connections between 6 surface and 4 borehole fractures and led to the discovery of 5 additional near-surface fracture zones. Of particular importance, they supplied key details on the subsurface geometries and minimum subsurface lengths of 8 of the 10 previously known surface fracture zones and all of the newly discovered ones. The vast majority of surface fracture zones extended at least 40–60 m into the subsurface, demonstrating that their depth and surface dimensions are comparable.     

利用钻孔雷达探测地下含水裂缝

钻孔雷达是地质雷达的一种特殊模式.和通常的地质雷达比起来,它有好几个显著的特点.比如,利用钻孔接近地下的某个区域,天线可以相对接近所要测量的地质异常体或目标,从而导致在较深的范围对目标的精确测量.实验场地位于北京西部的一个石灰岩的小山上.该处有一组钻孔,并被大量的裂缝所切割.在这些钻孔中,我们进行了单孔反射测量和跨孔测量.测量的数据经过处理和解释表明,在该场地雷达的径向探测范围可达30米,地下裂缝的分布可以形成清晰的图像.我们可以看到很多裂缝,它们距钻孔的距离及倾角都可以确定下来.在有些情况下,确定裂缝分布的方位是可能的.     

Origin and karst geomorphological significance of the enigmatic Australian Nullarbor Plain ‘blowholes’

The Australian Nullarbor Plain, one of the world's largest limestone platforms (~200 000 km²), has few distinctive surface karst features for its size, but is known for its enigmatic ‘blowholes’, which can display strong barometric draughts. Thousands of these vertical tubes with decimetre–metre (dm–m) scale diameter puncture the largely featureless terrain. The cause and distribution of these has remained unclear, but they have been thought to originate from downward dissolution and/or salt weathering. To elucidate blowhole distribution and mode of formation we (i) correlated existing location data with Shuttle Radar Topography Mission (SRTM) data, which distinguishes the subtle undulations (< 10 m per km) of the landscape, (ii) surveyed blowhole morphology and (iii) determined their rock surface hardness. Over a sampled area of 4200 km², the distribution of 615 known blowholes is not correlated with present topography. Blowholes are often connected to small or, in some cases extensive, but typically shallow cavities, which exhibit numerous ‘cupolas’ (dome‐shaped pockets) in their ceilings. Statistical arguments suggest that cavities with cupolas are common, but in only a few cases do these puncture the surface. Hardness measurements indicate that salt weathering is not their main cause. Our observations suggest that blowholes do not develop downwards, but occur where a cupola breaks through the surface. Lowering of the land surface is suggested to be the main cause for this breakthrough. Although cupolas may undergo some modification under the current climate, they, as well as the shallow caves they are formed in, are likely to be palaeokarst features formed under a shallower water table and wetter conditions in the past. The findings presented have implications for theories of dissolutional forms development in caves worldwide. The environmental history of the Nullarbor platform allows testing of such theories, because many other factors, which complicate karst evolution elsewhere, have not interfered with landform evolution here. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrological model sensitivity to parameter and radar rainfall estimation uncertainty

Radar estimates of rainfall are being increasingly applied to flood forecasting applications. Errors are inherent both in the process of estimating rainfall from radar and in the modelling of the rainfall–runoff transformation. The study aims at building a framework for the assessment of uncertainty that is consistent with the limitations of the model and data available and that allows a direct quantitative comparison between model predictions obtained by using radar and raingauge rainfall inputs. The study uses radar data from a mountainous region in northern Italy where complex topography amplifies radar errors due to radar beam occlusion and variability of precipitation with height. These errors, together with other error sources, are adjusted by applying a radar rainfall estimation algorithm. Radar rainfall estimates, adjusted and not, are used as an input to TOPMODEL for flood simulation over the Posina catchment (116 km²). Hydrological model parameter uncertainty is explicitly accounted for by use of the GLUE (Generalized Likelihood Uncertainty Estimation). Statistics are proposed to evaluate both the wideness of the uncertainty limits and the percentage of observations which fall within the uncertainty bounds. Results show the critical importance of proper adjustment of radar estimates and the use of radar estimates as close to ground as possible. Uncertainties affecting runoff predictions from adjusted radar data are close to those obtained by using a dense raingauge network, at least for the lowest radar observations available. Copyright © 2004 John Wiley & Sons, Ltd.     

Quantifying snow water equivalent using terrestrial ground penetrating radar and unmanned aerial vehicle photogrammetry

This study demonstrates the potential value of a combined unmanned aerial vehicle (UAV) Photogrammetry and ground penetrating radar (GPR) approach to map snow water equivalent (SWE) over large scales. SWE estimation requires two different physical parameters (snow depth and density), which are currently difficult to measure with the spatial and temporal resolution desired for basin-wide studies. UAV photogrammetry can provide very high-resolution spatially continuous snow depths (SD) at the basin scale, but does not measure snow densities. GPR allows nondestructive quantitative snow investigation if the radar velocity is known. Using photogrammetric snow depths and GPR two-way travel times (TWT) of reflections at the snow-ground interface, radar velocities in snowpack can be determined. Snow density (RSN) is then estimated from the radar propagation velocity (which is related to electrical permittivity of snow) via empirical formulas. A Phantom-4 Pro UAV and a MALA GX450 HDR model GPR mounted on a ski mobile were used to determine snow parameters. A snow-free digital surface model (DSM) was obtained from the photogrammetric survey conducted in September 2017. Then, another survey in synchronization with a GPR survey was conducted in February 2019 whilst the snowpack was approximately at its maximum thickness. Spatially continuous snow depths were calculated by subtracting the snow-free DSM from the snow-covered DSM. Radar velocities in the snowpack along GPR survey lines were computed by using UAV-based snow depths and GPR reflections to obtain snow densities and SWEs. The root mean square error of the obtained SWEs (384 mm average) is 63 mm, indicating good agreement with independent SWE observations and the error lies within acceptable uncertainty limits.     

Extrusion and gravity current of a fluid: Implications for salt tectonics

Numerical models of the extrusion and gravity current of a viscous incompressible fluid are studied to determine the shapes of salt structures formed on the Earth’s surface and the velocities of rock salt extrusion and lateral spreading. Two main types of salt extrusion are examined. In the case of active extrusion, salt rises to the surface at a velocity of about 30–35 cm/yr, forming a salt dome about 550 m in height. In the case of passive extrusion, the velocity of salt extrusion from under a newly developing sedimentary minibasin is significantly lower. In the course of its evolution, the salt dome becomes topographically lower and transforms into a broad plateau. The extrusion velocity of salt controls the shape, size, and velocity of its gravity current. The shapes of salt domes modeled in this study agree well with observations. The gravity current velocities in the models vary from 3 m/yr to 60 cm/yr, depending on the proximity to the current orifice. Numerical modeling of salt extrusion and gravity current in various geodynamic settings can be used for a detailed analysis of the geological and geophysical evolution of structures containing salt layers and the related oil and gas fields.     

Radar tomography for NDT: comparison of techniques

Feasibility and potential of tomography by Ground Penetrating Radar are investigated through experiments on laboratory models. The aim is the development of radar tomography procedures for inspection of structures like walls or pillars in historical buildings. Two different approaches are explored to satisfy high-resolution requirements. The first approach improves the results of classical traveltime (TT) and amplitude tomography (AT) on thin straight or curved rays through a progressive reduction of the null space of the problem. TT is a quantitative tool based on the thin ray assumption that allows a good tradeoff between robustness and resolution. AT is as robust as TT, but its results have only qualitative contents, since the energy transferred to the medium is basically unknown and the scattering effects are not taken into account. In the second approach, GPR is considered as a diffracting source, so that migration (MIG) and diffraction tomography (DT) are applied to overcome the geometrical optic approximations. While DT is in principle the best tool to invert the scattered field and to achieve the maximum resolution, MIG can be a more robust solution that requires less preprocessing of the data. All these advantages and drawbacks of the different approaches are discussed with some examples on synthetic and real data.     

AVTIS: A novel millimetre-wave ground based instrument for volcano remote sensing

During May and June 2004, we tested a novel, millimetre-wave ground based, dual-mode radar and radiometer on Soufrière Hills Volcano, Montserrat. AVTIS (All-weather Volcano Topography Imaging Sensor) has an active (radar) mode designed to image the distance to lava domes on volcanoes whose summits are commonly obscured by persistent cloud, such as Soufrière Hills Volcano. The passive (radiometer) mode can be employed to measure the surface temperature of the imaged topography. In its current form, AVTIS can be deployed by two people and takes 50 min to acquire a 20° × 5° scene at 0.1° increments. During the fieldwork period the lava dome was not growing and only the radar mode was used. The data recorded indicate that the maximum distance imaged was about 3800 m. Combining datasets acquired from different viewpoints can potentially provide a full 3D topographic model. The accuracy and completeness of this reconstruction are reduced by two factors. Firstly, relatively small grazing angles of the ground-based line-of-sight rendered incised valleys invisible. Secondly, methods currently used to orient the instrument limit the accuracy of the resulting topographic information. Nevertheless, valuable information on new topographic surfaces was obtained in an area north of the lava dome where a valley has been infilled by deposits and in the amphitheatre created by the giant collapse event of July 2003.     

Coupled groundwater flow and transport: 1. Verification of variable density flow and transport models

This work examines variable density flow and corresponding solute transport in groundwater systems. Fluid dynamics of salty solutions with significant density variations are of increasing interest in many problems of subsurface hydrology. The mathematical model comprises a set of non-linear, coupled, partial differential equations to be solved for pressure/hydraulic head and mass fraction/concentration of the solute component. The governing equations and underlying assumptions are developed and discussed. The equation of solute mass conservation is formulated in terms of mass fraction and mass concentration. Different levels of the approximation of density variations in the mass balance equations are used for convection problems (e.g. the Boussinesq approximation and its extension, fully density approximation). The impact of these simplifications is studied by use of numerical modelling.Numerical models for nonlinear problems, such as density-driven convection, must be carefully verified in a particular series of tests. Standard benchmarks for proving variable density flow models are the Henry, Elder, and salt dome (HYDROCOIN level 1 case 5) problems. We studied these benchmarks using two finite element simulators - ROCKFLOW, which was developed at the Institute of Fluid Mechanics and Computer Applications in Civil Engineering and FEFLOW, which was developed at the Institute for Water Resources Planning and Systems Research Ltd. Although both simulators are based on the Galerkin finite element method, they differ in many approximation details such as temporal discretization (Crank-Nicolson vs predictor-corrector schemes), spatial discretization (triangular and quadrilateral elements), finite element basis functions (linear, bilinear, biquadratic), iteration schemes (Newton, Picard) and solvers (direct, iterative). The numerical analysis illustrates discretization effects and defects arising from the different levels of the density of approximation. We contribute new results for the salt dome problem, for which inconsistent findings exist in literature. Applications of the verified numerical models to more complex problems, such as thermohaline and three-dimensional convection systems, will be presented in the second part of this paper.     

Predicted signatures of pulsed reconnection in ESR data

Early in 1996, the latest of the European inco-herent-scatter (EISCAT) radars came into operation on the Svalbard islands. The EISCAT Svalbard Radar (ESR) has been built in order to study the ionosphere in the northern polar cap and in particular, the dayside cusp. Conditions in the upper atmosphere in the cusp region are complex, with magnetosheath plasma cascading freely into the atmosphere along open magnetic field lines as a result of magnetic reconnection at the dayside magnetopause. A model has been developed to predict the effects of pulsed reconnection and the subsequent cusp precipitation in the ionosphere. Using this model we have successfully recreated some of the major features seen in photometer and satellite data within the cusp. In this paper, the work is extended to predict the signatures of pulsed reconnection in ESR data when the radar is pointed along the magnetic field. It is expected that enhancements in both electron concentration and electron temperature will be observed. Whether these enhancements are continuous in time or occur as a series of separate events is shown to depend critically on where the open/closed field-line boundary is with respect to the radar. This is shown to be particularly true when reconnection pulses are superposed on a steady background rate.