Two‐dimensional joint inversions of cross‐hole resistivity data and resolution analysis of combined arrays

In this study, a new two‐dimensional inversion algorithm was developed for the inversion of cross‐hole direct current resistivity measurements. In the last decades, various array optimisation methods were suggested for resistivity tomography. However, researchers have still collected data by using classical electrode arrays in most cross‐hole applications. Therefore, we investigated the accuracy of both the individual and the joint inversion of the classical cross‐hole arrays by using both synthetic and field data with the developed algorithm. We showed that the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole electrode arrays gives inverse solutions that are closer to the real model than the individual inversions of the electrode array datasets for the synthetic data inversion. The model resolution matrix of the suggested arrays was used to analyse the inversion results. This model resolution analysis also showed the advantage of the joint inversion of bipole–bipole, pole–bipole, bipole–pole, and pole–tripole arrays. We also used sensitivity sections from each of the arrays and their superpositions to explain why joint inversion gives better resolution than the any individual inversion result.     

Cross‐hole resistivity tomography using different electrode configurations

This paper investigates the relative merits and effectiveness of cross‐hole resistivity tomography using different electrode configurations for four popular electrode arrays: pole–pole, pole–bipole, bipole–pole and bipole–bipole. By examination of two synthetic models (a dipping conductive strip and a dislocated fault), it is shown that besides the popular pole–pole array, some specified three‐ and four‐electrode configurations, such as pole–bipole AM–N, bipole–pole AM–B and bipole–bipole AM–BN with their multispacing cross‐hole profiling and scanning surveys, are useful for cross‐hole resistivity tomography. These configurations, compared with the pole–pole array, may reduce or eliminate the effect of remote electrodes (systematic error) and yield satisfactory images with 20% noise‐contaminated data. It is also shown that the configurations which have either both current electrodes or both potential electrodes in the same borehole, i.e. pole–bipole A–MN, bipole–pole AB–M and bipole–bipole AB–MN, have a singularity problem in data acquisition, namely low readings of the potential or potential difference in cross‐hole surveying, so that the data are easily obscured by background noise and yield images inferior to those from other configurations.     

Numerical modelling of resolution and sensitivity of ERT in horizontal boreholes

Resistivity in horizontal boreholes can give useful detailed information about the geological conditions for construction in rock, i.e. in front of a tunnel bore machine. This paper is an attempt to identify a suitable methodology for an effective measuring routine for this type of geophysical measurements under actual construction site conditions.Prior to any measurements numerical modelling was done in order to evaluate the resolution of different electrode arrays. Four different arrays were tested; dipole–pole, cross-hole dipole–dipole, cross-hole pole–tripole and multiple gradient array. Additionally the resolution of a combination of cross-hole dipole–dipole and multiple gradient was assessed. The 2D sensitivity patterns for various arrangements of the cross-hole dipole–dipole and multiple gradient array were examined. The sensitivity towards inaccurate borehole geometry and the influence of water in the boreholes were also investigated. Based on the model study the cross-hole dipole–dipole array, multiple gradient array and a combination of these were found to give the best result and therefore were used for test measurements in horizontal boreholes. The boreholes were 28.5 m long and drilled 6.5 m apart. Prototypes of semi-rigid borehole cables made it possible to insert multi electrode cables in an efficient way, allowing fast measurement routines. These measurements were then studied to determine their accuracy and applicability. The results showed a high resistivity rock mass at the site. A transition from high resistivity to slightly lower resistivity coincides well with a change in lithology from gneiss-granite to gneiss. It is likely that the shotcrete on the tunnel wall is seen as a low resistivity zone.The measurements are a valuable tool, but further development of the cables and streamlining of measuring routines have to be performed before the resistivity tomography can be used routinely in pilot holes during construction in rocks.     

A fracture study by using bipole–bipole cross-well logging

Electrical Resistivity Tomography (ERT) can provide images of subsurface electrical structure between two boreholes. Data quality control is a key issue before ERT inversion. However, there is no effective data quality control method on an ERT survey. In this paper, a method called common current gather for a bipole–bipole array (CCGbb) was proposed to check ERT data quality in a rapid way. Synthetic models were conducted to compare the response difference between pole–pole array and bipole–bipole array. A field work at granite area was tested to verify the applicability of the proposed CCGbb method. From the results of this study, we suggested that conducting CCGbb before ERT inversion and a cross-borehole tracer test for both field data quality control and possible water conducting fractures (WCFs) delineation.     

2D Optimized Electrode Arrays for Borehole Resistivity Tomography and CO2 Sequestration Modelling

Modern optimization approaches for electrode configurations can significantly improve the resolution of 2.5D resistivity imaging surveys. This study presents a brief review of the 2.5D optimization approach, particularly for borehole–borehole surveys with applications for mapping virtual CO₂ plumes sequestrated in deep saline reservoir formations. The applied algorithm searches for arrays that maximize the spatial resolution of the survey among the comprehensive dataset of best possible spatial resolution (i.e. least temporal resolution). A main goal of this study is to increase the temporal resolution of ERT borehole–borehole surveys by selecting optimized electrode configurations in order to minimise the required data acquisition time while sustaining a high spatial resolution. The optimized dataset starts with a base set and is iteratively increased based on the model resolution matrix (R ) until the required number of data points is achieved. Among four different optimization methods, the compare R (CR) method of the best resolution is applied to directly calculate R for each new array added to the optimized dataset. Small optimized datasets generated by this technique are only <5% of their comprehensive sets but of an average resolution ratio (R _r) of >0.95 (i.e. almost the same resolution). With increasing the size of the optimized dataset (during its generation), the algorithm progressively enhances R _r values in the central interwell region (of low sensitivities and low resolution) far higher than in the near borehole region (of high sensitivities). Also the inverted tomogram reliability increases by increasing the optimized data size. Briefly, the optimized arrays improve the resolution in the interwell region which is commonly low in borehole–borehole ERT studies. The inverted output model is evaluated quantitatively using the model difference relative to the input model. The results reflect the common smearing effects and artefacts of varying degrees that overpredict volumes, underpredict magnitudes and blur boundaries of the target anomalies. This input model is a synthetic resistivity model that was used to generate synthetic (forward solution) data used during the inversion. Applications on synthetic CO₂ models show that the mapping resolution for optimized datasets is better than that for other highly resolving arrays of the same number of data points. Problems of smeared boundaries and thin layers are less visible in the optimized array than in the other highly resolving arrays.     

A numerical comparison of 2D resistivity imaging with 10 electrode arrays

Numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays. The arrays analysed include pole‐pole (PP), pole‐dipole (PD), half‐Wenner (HW), Wenner‐α (WN), Schlumberger (SC), dipole‐dipole (DD), Wenner‐β (WB), γ‐array (GM), multiple or moving gradient array (GD) and midpoint‐potential‐referred measurement (MPR) arrays. Five synthetic geological models, simulating a buried channel, a narrow conductive dike, a narrow resistive dike, dipping blocks and covered waste ponds, were used to examine the surveying efficiency (anomaly effects, signal‐to‐noise ratios) and the imaging capabilities of these arrays. The responses to variations in the data density and noise sensitivities of these electrode configurations were also investigated using robust (L₁‐norm) inversion and smoothness‐constrained least‐squares (L₂‐norm) inversion for the five synthetic models. The results show the following. (i) GM and WN are less contaminated by noise than the other electrode arrays. (ii) The relative anomaly effects for the different arrays vary with the geological models. However, the relatively high anomaly effects of PP, GM and WB surveys do not always give a high‐resolution image. PD, DD and GD can yield better resolution images than GM, PP, WN and WB, although they are more susceptible to noise contamination. SC is also a strong candidate but is expected to give more edge effects. (iii) The imaging quality of these arrays is relatively robust with respect to reductions in the data density of a multi‐electrode layout within the tested ranges. (iv) The robust inversion generally gives better imaging results than the L₂‐norm inversion, especially with noisy data, except for the dipping block structure presented here. (v) GD and MPR are well suited to multichannel surveying and GD may produce images that are comparable to those obtained with DD and PD. Accordingly, the GD, PD, DD and SC arrays are strongly recommended for 2D resistivity imaging, where the final choice will be determined by the expected geology, the purpose of the survey and logistical considerations.     

Out‐of‐plane effects in 2D borehole‐to‐surface resistivity tomography and applications in mineral exploration

In this paper, we discuss the effects of anomalous out‐of‐plane bodies in two‐dimensional (2D) borehole‐to‐surface electrical resistivity tomography with numerical resistivity modelling and synthetic inversion tests. The results of the two groups of synthetic resistivity model tests illustrate that anomalous bodies out of the plane of interest have an effect on two‐dimensional inversion and that the degree of influence of out‐of‐plane body on inverted images varies. The different influences are derived from two cases. One case is different resistivity models with the same electrode array, and the other case is the same resistivity model with different electrode arrays. Qualitative interpretation based on the inversion tests shows that we cannot find a reasonable electrode array to determine the best inverse solution and reveal the subsurface resistivity distribution for all types of geoelectrical models. Because of the three‐dimensional effect arising from neighbouring anomalous bodies, the qualitative interpretation of inverted images from the two‐dimensional inversion of electrical resistivity tomography data without prior information can be misleading. Two‐dimensional inversion with drilling data can decrease the three‐dimensional effect. We employed two‐ and three‐dimensional borehole‐to‐surface electrical resistivity tomography methods with a pole–pole array and a bipole–bipole array for mineral exploration at Abag Banner and Hexigten Banner in Inner Mongolia, China. Different inverse schemes were carried out for different cases. The subsurface resistivity distribution obtained from the two‐dimensional inversion of the field electrical resistivity tomography data with sufficient prior information, such as drilling data and other non‐electrical data, can better describe the actual geological situation. When there is not enough prior information to carry out constrained two‐dimensional inversion, the three‐dimensional electrical resistivity tomography survey is the better choice.     

Electrical resistivity tomography survey for delineating uncharted mine galleries in West Bengal,India*

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.     

Electrical resistivity characterization and defect detection on a geosynthetic clay liner (GCL) on an experimental site

In this paper we analyze the onsite characterization of a geosynthetic clay liner (GCL) that serves to ensure the impermeability of a landfill cap by DC electrical methods. The imaging of the GCL geoelectrical properties is a challenging problem because it is a very thin (between 4 and 7 mm thick) and resistive layer (from 100,000 to 2,000,000 Ω·m) depending on meteorological conditions and aging. We compare results obtained using electrical resistivity tomography (ERT) using two different kinds of arrays (dipole–dipole DD and Wenner–Schlumberger) on an experimental site with engineered defects. To confirm these results and to find the real onsite GCL resistivity we have performed sampling of the posterior distribution of this parameter using vertical electrical sounding (VES) inversions. Different VES methods were extracted from ERT with DD array and converted into a Schlumberger array.As a main conclusion the dipole–dipole array provides a better resistivity resolution of the defects than the Wenner–Schlumberger array. On ERT images, the defect detection seems to be impossible if the GCL has very high resistivity, as it happened when it was put in place. Taking into account the equivalence rules, the inversions are in both cases (ERT and VES) compatible. The GCL resistivity estimated from PSO (particle swarm optimization) varies from 3.0 10⁵ to 1.10⁶ Ω·m depending on saturation conditions during the twenty first months of its placing. Then, the resistivity dropped to 4.10⁴–9.10⁴ Ω·m, indicating a probable chemical damage of the GCL due to aging. Finally the fact that the VES inversions are solved via PSO sampling allows for the detection of a very thin and resistive layer and opens the possibility of performing micro VES surveys along the landfill to detect possible GCL defects.     

Three‐dimensional VS profiling using microtremors in Kushiro,Japan

A practical method is presented for determining three‐dimensional S‐wave velocity (V_S) profile from microtremor measurements. Frequency–wave number (f–k) spectral analyses of microtremor array records are combined, for this purpose, with microtremor horizontal‐to‐vertical (H/V) spectral ratio techniques. To demonstrate the effectiveness of the proposed method, microtremor measurements using arrays of sensors were conducted at six sites in the city of Kushiro, Japan. The spectral analyses of the array records yield dispersion characteristics of Rayleigh waves and H/V spectra of surface waves, and joint inversion of these data results in V_S profiles down to bedrock at the sites. Conventional microtremor measurements were performed at 230 stations within Kushiro city, resulting in the H/V spectra within the city. Three‐dimensional V_S structure is then estimated from inversion of the H/V spectra with the V_S values determined from the microtremor array data. This reveals three‐dimensional V_S profile of Kushiro city, together with an unknown hidden valley that crosses the central part of the city. The estimated V_S profile is consistent with available velocity logs and results of subsequent borings, indicating the effectiveness of the proposed method. Copyright © 2008 John Wiley & Sons, Ltd.     

Pilot‐scale field validation of the long electrode electrical resistivity tomography method

A validation experiment, carried out in a scaled field setting, was attempted for the long electrode electrical resistivity tomography method in order to demonstrate the performance of the technique in imaging a simple buried target. The experiment was an approximately 1/17 scale mock‐up of a region encompassing a buried nuclear waste tank on the Hanford site. The target of focus was constructed by manually forming a simulated plume within the vadose zone using a tank waste simulant. The long electrode results were compared to results from conventional point electrodes on the surface and buried within the survey domain. Using a pole‐pole array, both point and long electrode imaging techniques identified the lateral extents of the pre‐formed plume with reasonable fidelity but the long electrode method was handicapped in reconstructing vertical boundaries. The pole‐dipole and dipole‐dipole arrays were also tested with the long electrode method and were shown to have the least favourable target properties, including the position of the reconstructed plume relative to the known plume and the intensity of false positive targets. The poor performance of the pole‐dipole and dipole‐dipole arrays was attributed to an inexhaustive and non‐optimal coverage of data at key electrodes, as well as an increased noise for electrode combinations with high geometric factors. However, when comparing the model resolution matrix among the different acquisition strategies, the pole‐dipole and dipole‐dipole arrays using long electrodes were shown to have significantly higher average and maximum values within the matrix than any pole‐pole array. The model resolution describes how well the inversion model resolves the subsurface. Given the model resolution performance of the pole‐dipole and dipole‐dipole arrays, it may be worth investing in tools to understand the optimum subset of randomly distributed electrode pairs to produce maximum performance from the inversion model.     

Precise Relative Location of 25-ton Chemical Explosions at Balapan Using IMS Stations

—?We test how well low-magnitude (m _bLg 1.8 to 2.6), 25-ton chemical explosions at Balapan, Kazakhstan, can be located using IMS stations and standard earth models, relying on precisely determined relative arrival times of nearly similar, regional and teleseismic waveforms. Three 1997 Balapan explosions were recorded by a number of currently reporting and surrogate IMS stations. Three regional stations and two teleseismic arrays yielded consistent waveforms appropriate for relative picking. Master-event locations based on the AK135 model and ground-truth information from the first, shallowest and best-recorded explosion, fell under 1 km from known locations, for depths constrained to that of the master event. The resulting 90% confidence ellipses covered 12–13?km² and contained the true locations; however, results for depth constrained to true depth were slightly less satisf actory. From predictions based on ground truth, we found a P _g -coda phase at Makanchi, Kazakhstan to be misidentified and poorly modeled. After accounting for this, 90% ellipses shrank to 2–3?km² and true-depth mislocation vectors became more consistent with confidence-ellipse orientations. These results suggest that a high level of precision could be provided by a tripartite array of calibration shots in cases where models are poorly known. We hope that the successful relocation of these small Balapan shots will support the role of calibration explosions in verification monitoring and special event studies, including on-site inspection.     

Geophysical investigation of tumuli by means of surface 3D Electrical Resistivity Tomography

Tumuli are artificially erected small hills that cover monumental tombs or graves. In this work, the surface three-dimensional (3D) Electrical Resistivity Tomography (ERT) method, composed of dense parallel two-dimensional (2D) tomographies, was used to investigate the properties of the tumuli filling material and to resolve buried archaeological structures inside the tumuli.The effectiveness of the method was investigated by numerical modeling and through 3D inversion of synthetic apparent resistivity data. A resistivity model that simulates the inhomogeneous tumulus material and the tombs that are buried inside the tumulus was assumed. The Dipole–Dipole (DD), Pole–Dipole (PD), Pole–Pole (PP), Gradient (GRAD), Midpoint-Potential-Referred (MPR) and Schlumberger Reciprocal (SCR) arrays, which are suitable for multichannel resistivity instruments, were tested. The tumulus topography (pyramid or capsized cup) was incorporated into the inversion procedure through a distorted finite element mesh. The inversion procedure was based on a smoothness constrained Gauss–Newton algorithm in which the Active Constraint Balancing (ACB) method was also applied in order to enhance the least-squares resolving power and stability.Synthetic modeling showed that the different tumulus layers and the horizontal contact of the artificial tumulus material with the natural background soil were reconstructed by all of the tested electrode arrays. Generally, PD and the GRAD arrays comprise the optimum choices to investigate the subsurface properties of a tumulus and locate buried tombs. The MPR model was inferior to the GRAD model, while the DD, PP and SCR models had the poorest resolution. It was also shown that the inversion models are practically independent from the survey direction and the topography shape of the tumulus.The real field data collected employing the PD array along a small tumulus from the archaeological site of Vergina in northern Greece enhanced the synthetic modeling results. The inversion model outlined a number of archaeological structures that exhibit a high possibility to correlate with graves. Overall, this work signifies that the surface 3D ERT method can provide a valuable tool in the non-destructive archaeological exploration of tumuli.     

A 3D ERT study of solute transport in a large experimental tank

A high resolution, cross-borehole, 3D electrical resistivity tomography (ERT) study of solute transport was conducted in a large experimental tank. ERT voxels comprising the time sequence of electrical images were converted into a 3D array of ERT estimated fluid conductivity breakthrough curves and compared with direct measurements of fluid conductivity breakthrough made in wells. The 3D ERT images of solute transport behaviour were also compared with predictions based on a 3D finite-element, coupled flow and transport model, accounting for gravity induced flow caused by concentration differences.The tank (dimensions 185×245×186 cm) was filled with medium sand, with a gravel channel and a fine sand layer installed. This heterogeneous system was designed to complicate solute transport behaviour relative to a homogeneous sand tank, and to thus provide a challenging but insightful analysis of the ability of 3D ERT to resolve transport phenomena. Four ERT arrays and 20 piezometers were installed during filling. A NaCl tracer (conductivity 1.34 S/m) was injected and intensively monitored with 3D ERT and direct sampling of fluid chemistry in piezometers.We converted the bulk conductivity estimate for 250 voxels in the ERT imaged volume into ERT estimated voxel fluid conductivity by assuming that matrix conduction in the tank is negligible. In general, the ERT voxel response is in reasonable agreement with the shape of fluid conductivity breakthrough observed in six wells in which direct measurements of fluid conductivity were made. However, discrepancies occur, particularly at early times, which we attribute to differences between the scale of the image voxels and the fluid conductivity measurement, measurement errors mapped into the electrical inversion and artificial image roughness resulting from the inversion.ERT images revealed the 3D tracer distribution at 15 times after tracer injection. The general pattern and timing of solute breakthrough observed with ERT agreed with that predicted from the flow/transport modelling. However, the ERT images indicate a vertical component of tracer transport and preferential flow paths in the medium sand. We attribute this to transient vertical gradients established during tracer injection, and heterogeneity caused by sorting of the sand resulting from the filling procedure. In this study, ERT provided a unique dataset of 250 voxel breakthrough curves in 1.04 m³. The use of 3D ERT to generate an array of densely sampled estimated fluid conductivity breakthrough curves is a potentially powerful tool for quantifying solute transport processes.     

GPR to constrain ERT data inversion in cavity searching: Theoretical and practical applications in archeology

I used theoretical forward models to show that a cavity embedded in a stratified sedimentary sequence can induce an equivalence problem in the ERT data inversion. Conductive top soil increases the misfit between the ground feature and the ERT model. The misfit depends on array and stratigraphy sequences. The latter induce an equivalence problem that manifests itself as wrong cavity depth positioning. The misfit is greater in the data acquired with Schlumberger array than with dipole–dipole.The ambiguity of ERT data inversion problems was tested in the detection of cavities linked to an 8th–6th century B.C. Sabine tomb, 3 m wide × 3 m long × 2 m high, excavated from a shaly gray volcanic ash (cinerite) layer covered by semi-lithoid tuff and top soil layers. In the real study I reduced the ambiguity in the inverse problem of ERT data using a priori information on geometry and resistivity of the cavity. The constrains were carried out from georadar data acquired with 80 and 200 MHz antenna. I demonstrate that this procedure has a practical application in cavity detection, and is a key to the reduction of the uncertainty inherent in the inversion process of ERT data.     

Increasing the effectiveness of electrical resistivity tomography using γ11n configurations

A new array type, i.e., the γ_11n arrays, is introduced in this paper, in which the sequence of the current (C) and potential (P) electrodes is CPCP, and the distance between the last two electrodes is n times the distance between the first two ones and that of the second one and the third one. These arrays are called quasinull arrays because they are—according to their array and behaviour—between the traditional and null arrays. It is shown by numerical modelling that, in detecting small‐effect inhomogeneity, these configurations may be more effective than the traditional ones, including the optimized Stummer configuration. Certain γ_11n configurations—especially the γ₁₁₂, γ_113, and γ₁₁₄—produced better results both in horizontal and vertical resolution investigations. Based on the numerical studies, the γ_11n configurations seem to be very promising in problems where the anomalies are similar to the numerically investigated ones, i.e., they can detect and characterize, e.g., tunnels, caves, cables, tubes, abandoned riverbeds, or discontinuity, in a clay layer with greater efficacy than those of the traditional configurations. γ_11n measurements need less data than traditional configurations; therefore, the time demand of electrical resistivity tomography measurements can be shortened by their use.     

PERFORMANCE OF 2000 AND 6000 PSI AIR GUNS: THEORY AND EXPERIMENT*

Air guns have been used in various applications for a number of years. They were first used in coal-mining operations and were operated at up to 16000 psi charge pressures. Later, single air guns, operated at 2000 psi, found application as an oceanographic survey tool. Air gun arrays were first used in offshore seismic exploration in the mid-1960's. These early arrays were several hundred cubic inches in total volume and were operated at 2000 psi; they were either tuned arrays or several large guns of the same size with wave-shape kits. Today's arrays have total volumes greater than 5000 cu in. and are typically operated at 2000 psi. Recently, higher-pressure, lower-volume arrays operated at 4000–5000 psi have been introduced; guns used in these arrays are descendants of the coal-mining gun. On first thought one would equate increased gun pressure linearly with the amplitude of the initial pulse. This is approximately true for the signature radiated by a “free-bubble” (no confining vessel) and recorded broadband. The exact relation depends on the depth at which the gun is operated; from solution of the free-bubble oscillation equation, the relation is If P_c,1= 6014.7 psia, P_c,2= 2014.7 psia and P_{O, 1}=P_{O, 2}= 25.8 psia (corresponding to absolute pressure at 25 ft water depth), then Experiments were conducted offshore California in deep water to determine the performance of several models of air guns at pressures ranging from 2000 to 6000 psi and gun volumes ranging from 5 to 300 cu in. At a given gun pressure, the initial acoustic pulse P_a correlated with gun volume V_c according to the classical relation For 1 ms sampled data the ratio varied between 4.5 and 5.5 dB depending on gun model. Pulse width of the 2000 psi signatures indicated they are compatible with 2 ms sample-rate recording while pulse width of the 6000 psi signatures was greater, indicating they are less compatible with 2 ms sample-rate recording. Conclusions reached were that 2000 psi air guns are more efficient than higher pressure guns and are more compatible with 2 ms sample-rate requirements.     

Hydrochemical variations during flood pulses in the south‐west China peak cluster karst: impacts of CaCO3–H2O–CO2 interactions

High‐resolution measurements of rainfall, water level, pH, conductivity, temperature and carbonate chemistry parameters of groundwater at two adjacent locations within the peak cluster karst of the Guilin Karst Experimental Site in Guangxi Province, China, were made with different types of multiparameter sonde. The data were stored using data loggers recording with 2 min or 15 min resolution. Waters from a large, perennial spring represent the exit for the aquifer's conduit flow, and a nearby well measures water in the conduit‐adjacent, fractured media. During flood pulses, the pH of the conduit flow water rises as the conductivity falls. In contrast, and at the same time, the pH of groundwater in the fractures drops, as conductivity rises. As Ca²⁺ and HCO₃^? were the dominant (>90%) ions, we developed linear relationships (both r² > 0·91) between conductivity and those ions, respectively, and in turn calculated variations in the calcite saturation index (SI_C) and CO₂ partial pressure (P) of water during flood pulses. Results indicate that the P of fracture water during flood periods is higher than that at lower flows, and its SI_C is lower. Simultaneously, P of conduit water during the flood period is lower than that at lower flows, and its SI_C also is lower. From these results we conclude that at least two key processes are controlling hydrochemical variations during flood periods: (i) dilution by precipitation and (ii) water–rock–gas interactions. To explain hydrochemical variations in the fracture water, the water–rock–gas interactions may be more important. For example, during flood periods, soil gas with high CO₂ concentrations dissolves in water and enters the fracture system, the water, which in turn has become more highly undersaturated, dissolves more limestone, and the conductivity increases. Dilution of rainfall is more important in controlling hydrochemical variations of conduit water, because rainfall with higher pH (in this area apparently owing to interaction with limestone dust in the lower atmosphere) and low conductivity travels through the conduit system rapidly. These results illustrate that to understand the hydrochemical variations in karst systems, considering only water–rock interactions is not sufficient, and the variable effects of CO₂ on the system should be evaluated. Consideration of water–rock–gas interactions is thus a must in understanding variations in karst hydrochemistry. Copyright © 2004 John Wiley & Sons, Ltd.     

Application of the redundant-lifting scheme for ground-roll attenuation in near-surface characterization using full-waveform inversion on P-wave seismic data

Seismic modelling of the shallow subsurface (within the first few metres) is often challenging when the data are dominated by ground-roll and devoid of reflection. We showed that, even when transmission is the only available phase for analysis, fine-scale and interpretable P-wave velocity (V_P) and attenuation (Q_P⁻¹) models can still be prepared using full-waveform inversion, with data being preconditioned for ground-roll. To prove this idea, we suppressed the ground-roll in two different ways before full-waveform inversion modelling: first, through a bottom mute; second, through a novel wavelet transform-based method known as the redundant-lifting scheme. The applicability of full-waveform inversion is tested through imaging two buried targets. These include a pair of utility water pipes with known diameters of 0.8 m and burial depths of 1.5 m, respectively. The second target is the poorly documented backfill, which was the former location of the pipe(s). The data for full-waveform inversion are acquired along a 2D profile using a static array of 24, 40 Hz vertical component geophones and a buried point source. The results show that (a) the redundant-lifting scheme better suppresses the ground roll, which in turn provides better images of the targets in full-waveform inversion; and (b) the V_P and Q_P⁻¹ models from full-waveform inversion of redundant-lifting scheme data could detect the two targets adequately.     

Evaporation from three water bodies of different sizes and climates: Measurements and scaling analysis

Evaporation from small reservoirs, wetlands, and lakes continues to be a theoretical and practical problem in surface hydrology and micrometeorology because atmospheric flows above such systems can rarely be approximated as stationary and planar-homogeneous with no mean subsidence (hereafter referred to as idealized flow state). Here, the turbulence statistics of temperature (T) and water vapor (q) most pertinent to lake evaporation measurements over three water bodies differing in climate, thermal inertia and degree of advective conditions are explored. The three systems included Lac Léman in Switzerland (high thermal inertia, near homogeneous conditions with no appreciable advection due to long upwind fetch), Eshkol reservoir in Israel (intermediate thermal inertia, frequent strong advective conditions) and Tilopozo wetland in Chile (low thermal inertia, frequent but moderate advection). The data analysis focused on how similarity constants for the flux-variance approach, C_T/C_q, and relative transport efficiencies R_wT/R_wq, are perturbed from unity with increased advection or the active role of temperature. When advection is small and thermal inertia is large, C_T/C_q < 1 (or R_wT/R_wq > 1) primarily due to the active role of temperature, which is consistent with a large number of studies conducted over bare soil and vegetated surfaces. However, when advection is significantly large, then C_T/C_q > 1 (orR_wT/R_wq < 1). When advection is moderate and thermal inertia is low, then C_T/C_q ∼ 1. This latter equality, while consistent with Monin–Obukhov similarity theory (MOST), is due to the fact that advection tends to increase C_T/C_q above unity while the active role of temperature tends to decrease C_T/C_q below unity. A simplified scaling analysis derived from the scalar variance budget equation, explained qualitatively how advection could perturb MOST scaling (assumed to represent the idealized flow state).