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.

Effects of anomalous propagation conditions on weather radar observations

In this chapter our objective is to provide an overview of the effects of anomalous propagation conditions on weather radar observations, based mostly on studies performed by the authors during the last decade, summarizing results from recent publications, presentations, or unpublished material. We believe this chapter may be useful as an introductory text for graduate students, or researchers and practitioners dealing with this topic. Throughout the text a spherical symmetric atmosphere is assumed and the focus is on the occurrence of ground and sea clutter and subsequent problems for weather radar applications. Other related topics such as long-path, over-the-horizon propagation and detection of radar targets (either clutter or weather systems) at long ranges is not considered here; however readers should be aware of the potential problems these phenomena may have as range aliasing may cause these echoes appear nearer than they are ¿ for more details see the discussion about second trip echoes by Zrnic, this volume.

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.

Real-time full field laser Doppler imaging.

We present a full field laser Doppler imaging instrument, which enables real-time in vivo assessment of blood flow in dermal tissue and skin. This instrument monitors the blood perfusion in an area of about 50 cm(2) with 480 × 480 pixels per frame at a rate of 12-14 frames per second. Smaller frames can be monitored at much higher frame rates. We recorded the microcirculation in healthy skin before, during and after arterial occlusion. In initial clinical case studies, we imaged the microcirculation in burned skin and monitored the recovery of blood flow in a skin flap during reconstructive surgery indicating the high potential of LDI for clinical applications. Small animal imaging in mouse ears clearly revealed the network of blood vessels and the corresponding blood perfusion.

Bayesian Markov-chain-Monte-Carlo inversion of time-lapse cross hole ground-penetrating radar data to characterize the vadose zone at the Arrenaes field site, Denmark

The ground-penetrating radar (GPR) geophysical method has the potential to provide valuable information on the hydraulic properties of the vadose zone because of its strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR traveltime data can allow for a significant reduction in uncertainty regarding subsurface van Genuchten-Mualem (VGM) parameters. Much of the previous work on the stochastic estimation of VGM parameters from crosshole GPR data has considered the case of steady-state infiltration conditions, which represent only a small fraction of practically relevant scenarios. We explored in detail the dynamic infiltration case, specifically examining to what extent time-lapse crosshole GPR traveltimes, measured during a forced infiltration experiment at the Arreneas field site in Denmark, could help to quantify VGM parameters and their uncertainties in a layered medium, as well as the corresponding soil hydraulic properties. We used a Bayesian Markov-chain-Monte-Carlo inversion approach. We first explored the advantages and limitations of this approach with regard to a realistic synthetic example before applying it to field measurements. In our analysis, we also considered different degrees of prior information. Our findings indicate that the stochastic inversion of the time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions compared with the corresponding priors, which in turn significantly improves knowledge of soil hydraulic properties. Overall, the results obtained clearly demonstrate the value of the information contained in time-lapse GPR data for characterizing vadose zone dynamics.

Polarimetric radar interferometry for improved mine detection and surface clutter rejection

A recently developed technique, polarimetric radar interferometry, is applied to tackle the problem of the detection of buried objects embedded in surface clutter. An experiment with a fully polarimetric radar in an anechoic chamber has been carried out using different frequency bands and baselines. The processed results show the ability of this technique to detect buried plastic mines and to measure their depth. This technique enables the detection of plastic mines even if their backscatter response is much lower than that of the surface clutter.

Study of fish swimming activity using Acoustical Doppler Velocimetry (ADV) techniques

Peer reviewed

Reconstruction of the join time-delay Doppler-scale reflectivity density in the wide-band regime. A frame theory based approach

The inverse scattering problem concerning the determination of the joint time-delayDoppler-scale reflectivity density characterizing continuous target environments is addressed by recourse to the generalized frame theory. A reconstruction formula,involving the echoes of a frame of outgoing signals and its corresponding reciprocalframe, is developed. A ‘‘realistic’’ situation with respect to the transmission ofa finite number of signals is further considered. In such a case, our reconstruction formula is shown to yield the orthogonal projection of the reflectivity density onto a subspace generated by the transmitted signals.

Transcranial Doppler for predicting delayed cerebral ischemia after subarachnoid hemorrhage.

OBJECTIVE: Transcranial Doppler (TCD) is widely used to monitor the temporal course of vasospasm after subarachnoid hemorrhage (SAH), but its ability to predict clinical deterioration or infarction from delayed cerebral ischemia (DCI) remains controversial. We sought to determine the prognostic utility of serial TCD examination after SAH. METHODS: We analyzed 1877 TCD examinations in 441 aneurysmal SAH patients within 14 days of onset. The highest mean blood flow velocity (mBFV) value in any vessel before DCI onset was recorded. DCI was defined as clinical deterioration or computed tomographic evidence of infarction caused by vasospasm, with adjudication by consensus of the study team. Logistic regression was used to calculate adjusted odds ratios for DCI risk after controlling for other risk factors. RESULTS: DCI occurred in 21% of patients (n = 92). Multivariate predictors of DCI included modified Fisher computed tomographic score (P = 0.001), poor clinical grade (P = 0.04), and female sex (P = 0.008). After controlling for these variables, all TCD mBFV thresholds between 120 and 180 cm/s added a modest degree of incremental predictive value for DCI at nearly all time points, with maximal sensitivity by SAH day 8. However, the sensitivity of any mBFV more than 120 cm/s for subsequent DCI was only 63%, with a positive predictive value of 22% among patients with Hunt and Hess grades I to III and 36% in patients with Hunt and Hess grades IV and V. Positive predictive value was only slightly higher if mBFV exceeded 180 cm/s. CONCLUSION: Increased TCD flow velocities imply only a mild incremental risk of DCI after SAH, with maximal sensitivity by day 8. Nearly 40% of patients with DCI never attained an mBFV more than 120 cm/s during the course of monitoring. Given the poor overall sensitivity of TCD, improved methods for identifying patients at high risk for DCI after SAH are needed.

Aplicación de imágenes radar de satélite a la cartografía de la vegetación en zonas boreales

El presente estudio se enmarca en el proyecto europeo SIBERIA. Trata de explorar el uso de imágenes radar de satélite (ERS y JERS) para la actualización de la cartografía de vegetación de zonas boreales. Se dispone de 8 imágenes de amplitud y coherencia tomadas en 1998, así como de un inventario de vegetación georreferenciado de dos pequeñas zonas. Se proponen tres tipos de clasificaciones supervisadas por el método de máxima verosimilitud. La primera con las imágenes de satélite, la segunda añadiendo algunas imágenes texturales, y la tercera utilizando sólo las imágenes de los componentes principales más significativos. Se siguen los criterios establecidos en el proyecto SIBERIA para la obtención de áreas de entrenamiento. Se propone una doble validación, por una parte vía matrices de confusión a partir de áreas de verdad-terreno obtenidas por el mismo método que las áreas de entrenamiento, y por otra parte contrastando y correlacionando las clasificaciones con los parámetros de inventario disponibles para dos pequeñas áreas de verdad-terreno. Los resultados indican una sensible mejora en la clasificación con la incorporación de imágenes texturales (la precisión aumenta de un 66% a un 75%), y señalan el parámetro biomasa como el mejor correlacionado con las clasificaciones derivadas (coeficiente de correlación r de hasta 0,49). Diferentes fuentes de error permiten augurar un margen de mejora para posteriores estudios.