Inferring transport characteristics in a fractured rock aquifer by combining single-hole ground-penetrating radar reflection monitoring and tracer test data

Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Integration of ground-penetrating radar and microgravimetric methods to map shallow caves

Ground-penetrating radar (GPR) and microgravimetric surveys have been conducted in the southern Jura mountains of western Switzerland in order to map subsurface karstic features. The study site, La Grande Rolaz cave, is an extensive system in which many portions have been mapped. By using small station spacing and careful processing for the geophysical data, and by modeling these data with topographic information from within the cave, accurate interpretations have been achieved. The constraints on the interpreted geologic models are better when combining the geophysical methods than when using only one of the methods, despite the general limitations of two-dimensional (2D) profiling. For example, microgravimetry can complement GPR methods for accurately delineating a shallow cave section approximately 10 X 10 mt in size. Conversely, GPR methods can be complementary in determining cavity depths and in verifying the presence of off-line features and numerous areas of small cavities and fractures, which may be difficult to resolve in microgravimetric data.

Strokes and TIAs during and after Carotid Artery Doppler: Cause or Coincidence?

The aim of the present study was to explore the prevalence of acute cerebrovascular symptoms temporally related to carotid Doppler examination (DEx), in order to increase the awareness and recording of such events and to discuss possible mechanisms. All adult patients who complained of acute onset neurologic symptoms during or shortly after a carotid DEx, between 01/2003 and 12/2011 in the University Hospital of Lausanne were prospectively collected. We identified four consecutive patients with acute onset neurologic symptoms during or shortly after a carotid DEx among approximately 13,500 patients who underwent carotid DEx in our facility during the nine-year period (0.015% of all examined carotids). Clinical data, imaging reports and CTA (CT angiography) or/and ultrasound images are presented for each patient. Ischemic cerebrovascular events during or immediately after cervical Doppler could be due to chance or to several physical factors. They should be promptly recognized by Doppler personnel and properly treated, but do not put into question the overwhelming benefits of Doppler in cerebrovascular patients.

Reduced atrial emptying after orthotopic heart transplantation masquerading as restrictive transmitral Doppler flow pattern?

BACKGROUND: An elevated early (E) to late (A) diastolic filling velocities ratio, typically seen in advanced diastolic dysfunction, has also been observed after cardioversion of atrial fibrillation as a consequence of the depressed left atrial (LA) contractility. We hypothesized that the impaired LA contractile function demonstrated after orthotopic cardiac transplantation (OCT) could also lead to this "pseudorestrictive" pattern. METHOD: E/A ratio related to the tissue Doppler early mitral annular velocity (Ea) as preload-independent index of LV relaxation was evaluated in all consecutive OCT patients between 2005 and 2007. RESULTS: The study population comprised 48 patients 97 ± 77 months after OCT. Thirty-two patients (67%) had an E/A ratio > 2. LV systolic function and myocardial relaxation assessed by the Ea velocity were similar compared to patients with normal ratio (61 ± 6% vs. 60 ± 12%, P = 0.854 and 15 ± 4 cm/s vs. 14 ± 3 cm/s, r = 0.15, P = 0.323, respectively). On the other hand, the proportion of the recipient and donor LA cuffs as estimated by the recipient/global LA area ratio and the LA emptying fraction significantly correlated with the E/A ratio (r = 0.40, P = 0.005 and r =-0.33, P = 0.022, respectively). CONCLUSION: Our study shows that there is a high prevalence of elevated E/A ratio after standard OCT which seems mainly related to reduced LA contractility. Recognition of this "pseudorestrictive" pattern may avoid misdiagnosis of diastolic dysfunction.

FluxEXPLORER: a new high-speed laser Doppler imaging system for the assessment of burn injuries.

BACKGROUND: Deep burn assessment made by clinical evaluation has an accuracy varying between 60% and 80% and will determine if a burn injury will need tangential excision and skin grafting or if it will be able to heal spontaneously. Laser Doppler Imaging (LDI) techniques allow an improved burn depth assessment but their use is limited by the time-consuming image acquisition which may take up to 6 min per image. METHODS: To evaluate the effectiveness and reliability of a newly developed full-field LDI technology, 15 consecutive patients presenting with intermediate depth burns were assessed both clinically and by FluxExplorer LDI technology. Comparison between the two methods of assessment was carried out. RESULTS: Image acquisition was done within 6 s. FluxEXPLORER LDI technology achieved a significantly improved accuracy of burn depth assessment compared to the clinical judgement performed by board certified plastic and reconstructive surgeons (P < 0.05, 93% of correctly assessed burns injuries vs. 80% for clinical assessment). CONCLUSION: Technological improvements of LDI technology leading to a decreased image acquisition time and reliable burn depth assessment allow the routine use of such devices in the acute setting of burn care without interfering with the patient's treatment. Rapid and reliable LDI technology may assist clinicians in burn depth assessment and may limit the morbidity of burn patients through a minimization of the area of surgical debridement. Future technological improvements allowing the miniaturization of the device will further ease its clinical application.

Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole georadar data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a time-domain waveform inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little-to-no trade-off between the wavelet estimation and the tomographic imaging procedures.

Advances in Near-surface Seismology and Ground-penetrating Radar

Advances in Near-surface Seismology and Ground-penetrating Radar (SEG Geophysical Developments Series No. 15) is a collection of original papers by renowned and respected authors from around the world. Technologies used in the application of near-surface seismology and ground-penetrating radar have seen significant advances in the last several years. Both methods have benefited from new processing tools, increased computer speeds, and an expanded variety of applications. This book, divided into four sections ? ?Reviews,? ?Methodology,? ?Integrative Approaches,? and ?Case Studies? ? captures the most significant cutting-edge issues in active areas of research, unveiling truly pertinent studies that address fundamental applied problems. This collection of manuscripts grew from a core group of papers presented at a postconvention workshop, ?Advances in Near-surface Seismology and Ground-penetrating Radar,? held during the 2009 SEG Annual Meeting in Houston, Texas. This is the first cooperative publication effort between the near-surface communities of SEG, AGU, and EEGS. It will appeal to a large and diverse audience that includes researchers and practitioners inside and outside the near-surface geophysics community.

Source-wavelet estimation during full-waveform inversion of ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole ground-penetrating radar (GPR) data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a recently published time-domain inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity of both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little to no trade-off between the wavelet estimation and the tomographic imaging procedures.

Estimating vadose zone hydraulic properties using ground penetrating radar: The impact of prior information

A number of geophysical methods, such as ground-penetrating radar (GPR), have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, the stochastic inversion of such data within a coupled geophysical-hydrological framework may allow for the effective estimation of vadose zone hydraulic parameters and their corresponding uncertainties. A critical issue in stochastic inversion is choosing prior parameter probability distributions from which potential model configurations are drawn and tested against observed data. A well chosen prior should reflect as honestly as possible the initial state of knowledge regarding the parameters and be neither overly specific nor too conservative. In a Bayesian context, combining the prior with available data yields a posterior state of knowledge about the parameters, which can then be used statistically for predictions and risk assessment. Here we investigate the influence of prior information regarding the van Genuchten-Mualem (VGM) parameters, which describe vadose zone hydraulic properties, on the stochastic inversion of crosshole GPR data collected under steady state, natural-loading conditions. We do this using a Bayesian Markov chain Monte Carlo (MCMC) inversion approach, considering first noninformative uniform prior distributions and then more informative priors derived from soil property databases. For the informative priors, we further explore the effect of including information regarding parameter correlation. Analysis of both synthetic and field data indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when we combine these data with a realistic, informative prior.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

Transcranial Doppler Pulsatility Index: What it is and What it Isn't.

BACKGROUND: Transcranial Doppler (TCD) pulsatility index (PI) has traditionally been interpreted as a descriptor of distal cerebrovascular resistance (CVR). We sought to evaluate the relationship between PI and CVR in situations, where CVR increases (mild hypocapnia) and decreases (plateau waves of intracranial pressure-ICP). METHODS: Recordings from patients with head-injury undergoing monitoring of arterial blood pressure (ABP), ICP, cerebral perfusion pressure (CPP), and TCD assessed cerebral blood flow velocities (FV) were analyzed. The Gosling pulsatility index (PI) was compared between baseline and ICP plateau waves (n = 20 patients) or short term (30-60 min) hypocapnia (n = 31). In addition, a modeling study was conducted with the "spectral" PI (calculated using fundamental harmonic of FV) resulting in a theoretical formula expressing the dependence of PI on balance of cerebrovascular impedances. RESULTS: PI increased significantly (p < 0.001) while CVR decreased (p < 0.001) during plateau waves. During hypocapnia PI and CVR increased (p < 0.001). The modeling formula explained more than 65% of the variability of Gosling PI and 90% of the variability of the "spectral" PI (R = 0.81 and R = 0.95, respectively). CONCLUSION: TCD pulsatility index can be easily and quickly assessed but is usually misinterpreted as a descriptor of CVR. The mathematical model presents a complex relationship between PI and multiple haemodynamic variables.

Transcranial Doppler ultrasound in the acute phase of aneurysmal subarachnoid hemorrhage.

BACKGROUND: Angiographic studies suggest that acute vasospasm within 48 h of aneurysmal subarachnoid hemorrhage (SAH) predicts symptomatic vasospasm. However, the value of transcranial Doppler within 48 h of SAH is unknown. METHODS: We analyzed 199 patients who had at least 1 middle cerebral artery (MCA) transcranial Doppler examination within 48 h of SAH onset. Abnormal MCA mean blood flow velocity (mBFV) was defined as >90 cm/s. Delayed cerebral ischemia (DCI) was defined as clinical deterioration or radiological evidence of infarction due to vasospasm. RESULTS: Seventy-six patients (38%) had an elevation of MCA mBFV >90 cm/s within 48 h of SAH onset. The predictors of elevated mBFV included younger age (OR = 0.97 per year of age, p = 0.002), admission angiographic vasospasm (OR = 5.4, p = 0.009) and elevated white blood cell count (OR = 1.1 per 1,000 white blood cells, p = 0.003). Patients with elevated mBFV were more likely to experience a 10 cm/s fall in velocity at the first follow-up than those with normal baseline velocities (24 vs. 10%, p < 0.01), suggestive of resolving spasm. DCI developed in 19% of the patients. An elevated admission mBFV >90 cm/s during the first 48 h (adjusted OR = 2.7, p = 0.007) and a poor clinical grade (Hunt-Hess score 4 or 5, OR = 3.2, p = 0.002) were associated with a significant increase in the risk of DCI. CONCLUSION: Early elevations of mBFV correlate with acute angiographic vasospasm and are associated with a significantly increased risk of DCI. Transcranial Doppler ultrasound may be an early useful tool to identify patients at higher risk to develop DCI after SAH.