Derivation of soil-specific streaming potential electrical parameters from hydrodynamic characteristics of partially saturated soils

Water movement in unsaturated soils gives rise to measurable electrical potential differences that are related to the flow direction and volumetric fluxes, as well as to the soil properties themselves. Laboratory and field data suggest that these so-called streaming potentials may be several orders of magnitudes larger than theoretical predictions that only consider the influence of the relative permeability and electrical conductivity on the self potential (SP) data. Recent work has improved predictions somewhat by considering how the volumetric excess charge in the pore space scales with the inverse of water saturation. We present a new theoretical approach that uses the flux-averaged excess charge, not the volumetric excess charge, to predict streaming potentials. We present relationships for how this effective excess charge varies with water saturation for typical soil properties using either the water retention or the relative permeability function. We find large differences between soil types and the predictions based on the relative permeability function display the best agreement with field data. The new relationships better explain laboratory data than previous work and allow us to predict the recorded magnitudes of the streaming potentials following a rainfall event in sandy loam, whereas previous models predict values that are three orders of magnitude too small. We suggest that the strong signals in unsaturated media can be used to gain information about fluxes (including very small ones related to film flow), but also to constrain the relative permeability function, the water retention curve, and the relative electrical conductivity function.

Optimal 3-T MRI for depiction of the finger A2 pulley: comparison between T1-weighted, fat-saturated T2-weighted and gadolinium-enhanced fat-saturated T1-weighted sequences

OBJECTIVE: To compare three spin-echo sequences, transverse T1-weighted (T1WI), transverse fat-saturated (FS) T2-weighted (T2WI), and transverse gadolinium-enhanced (Gd) FS T1WI, for the visualisation of normal and abnormal finger A2 pulley with magnetic resonance (MR) imaging at 3 tesla (T). MATERIALS AND METHODS: Sixty-three fingers from 21 patients were consecutively investigated. Two musculoskeletal radiologists retrospectively compared all sequences to assess the visibility of normal and abnormal A2 pulleys and the presence of motion or ghost artefacts. RESULTS: Normal and abnormal A2 pulleys were visible in 94% (59/63) and 95% (60/63) on T1WI sequences, in 63% (40/63) and 60% (38/63) on FS T2WI sequences, and in 87% (55/63) and 73% (46/63) on Gd FS T1WI sequences when read by the first and second observer, respectively. Motion and ghost artefacts were higher on FS T2WI sequences. Seven among eight abnormal A2 pulleys were detected, and were best depicted with Gd FS T1WI sequences in 71% (5/7) and 86% (6/7) by the first and the second observer, respectively. CONCLUSION: In 3-T MRI, the comparison between transverse T1WI, FS T2WI, and Gd FS T1WI sequences shows that transverse T1WI allows excellent depiction of the A2 pulley, that FS T2WI suffers from a higher rate of motion and ghost artefacts, and transverse Gd FS T1WI is the best sequence for the depiction of abnormal A2 pulley.

Athletic injuries of the extensor carpi ulnaris subsheath: MRI findings and utility of gadolinium-enhanced fat-saturated T1-weighted sequences with wrist pronation and supination

Purpose: To report the magnetic resonance imaging (MRI) findings in athletic injuries of the extensor carpi ulnaris (ECU) subsheath, assessing the utility of gadolinium-enhanced (Gd) fat-saturated (FS) T1-weighted sequences with wrist pronation and supination. Methods and Materials: Sixteen patients (13 males, 3 females; mean age 30.3 years) with athletic injuries of the ECU subsheath sustained between January 2003 and June 2009 were included in this retrospective study. Initial and follow‑up 1.5-T wrist MRIs were performed with transverse T1-weighted and STIR sequences in pronation, and Gd FS T1-weighted sequences with wrist pronation and supination. Two radiologists assessed the type of injury (A to C), ECU tendon stability, associated lesions and rated pulse sequences using a three-point scale: 1 = poor, 2 = good and 3 = excellent. Results: Gd-enhanced FS T1-weighted transverse sequences in supination (2.63) and pronation (2.56) were most valuable, compared with STIR (2.19) and T1 weighted (1.94). Nine type A, one type B and six type C injuries were found. There were trends towards diminution in size, signal intensity and enhancement of associated pouches on follow‑up MRI and tendon stabilisation within the ulnar groove. Conclusion: Gd-enhanced FS T1-weighted sequences with wrist pronation and supination are most valuable in assessing and follow‑up athletic injuries of the ECU subsheath on 1.5-T MRI.

Seismic attenuation and velocity dispersion in heterogeneous partially saturated porous rocks

Using a numerical approach, we explore wave-induced fluid flow effects in partially saturated porous rocks in which the gas-water saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighbouring regions can exhibit different levels of saturation. To determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. The multiscale nature of the fluid saturation is taken into account by locally computing the physical properties of an effective fluid, which are then used for the larger-scale simulations. We consider two sets of numerical experiments to analyse such effects in heterogeneous partially saturated porous media, where the saturation field is determined by variations in porosity and clay content, respectively. In both cases we also evaluate the seismic responses of corresponding binary, patchy-type saturation patterns. Our results indicate that significant attenuation and modest velocity dispersion effects take place in this kind of media for both binary patchy-type and spatially continuous gas saturation patterns and in particular in the presence of relatively small amounts of gas. The numerical experiments also show that the nature of the gas distribution patterns is a critical parameter controlling the seismic responses of these environments, since attenuation and velocity dispersion effects are much more significant and occur over a broader saturation range for binary patchy-type gas-water distributions. This analysis therefore suggests that the physical mechanisms governing partial saturation should be accounted for when analysing seismic data in a poroelastic framework. In this context, heterogeneities associated with the dry frame properties, which do not play important roles in wave-induced fluid flow processes per se, should be taken into account since they may determine the kind of gas distribution pattern taking place in the porous rock.

Seismic attenuation in heterogeneous partially saturated rocks

We study wave-induced fluid flow effects in porous rocks partially saturated with gas and water, where the saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighboring regions can exhibit different levels of saturation. In order to determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. We consider numerical experiments to analyze such effects in heterogeneous partially saturated porous media, where the saturation field is determined by realistic variations in porosity. Our results indicate that the spatially continuous nature of gas saturation inherent to this study is a critical parameter controlling the seismic response of these environments, which in turn suggests that the physical mechanisms governing partial saturation should be accounted for when analyzing seismic data in a poro-elastic context.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

Athletic injuries of the extensor carpi ulnaris subsheath: MRI findings and utility of gadolinium-enhanced fat-saturated T1-weighted sequences with wrist pronation and supination.

OBJECTIVE: To report the magnetic resonance imaging (MRI) findings in athletic injuries of the extensor carpi ulnaris (ECU) subsheath, assessing the utility of gadolinium-enhanced (Gd) fat-saturated (FS) T1-weighted sequences with wrist pronation and supination. METHODS: Sixteen patients (13 male, three female; mean age 30.3 years) with athletic injuries of the ECU subsheath sustained between January 2003 and June 2009 were included in this retrospective study. Initial and follow-up 1.5-T wrist MRIs were performed with transverse T1-weighted and STIR sequences in pronation, and Gd FS T1-weighted sequences with wrist pronation and supination. Two radiologists assessed the type of injury (A to C), ECU tendon stability, associated lesions and rated pulse sequences using a three-point scale: 1=poor, 2=good and 3=excellent. RESULTS: Gd-enhanced FS T1-weighted transverse sequences in supination (2.63) and pronation (2.56) were most valuable, compared with STIR (2.19) and T1-weighted (1.94). Nine type A, one type B and six type C injuries were found. There were trends towards diminution in size, signal intensity and enhancement of associated pouches on follow-up MRI and tendon stabilisation within the ulnar groove. CONCLUSION: Gd-enhanced FS T1-weighted sequences with wrist pronation and supination are most valuable in assessing and follow-up athletic injuries of the ECU subsheath on 1.5-T MRI.

Self-potentials in partially saturated media: the importance of explicit modeling of electrode effects

Self-potential (SP) data are of interest to vadose zone hydrology because of their direct sensitivity to water flow and ionic transport. There is unfortunately little consensus in the literature about how to best model SP data under partially saturated conditions, and different approaches (often supported by one laboratory data set alone) have been proposed. We argue that this lack of agreement can largely be traced to electrode effects that have not been properly taken into account. A series of drainage and imbibition experiments were considered in which we found that previously proposed approaches to remove electrode effects were unlikely to provide adequate corrections. Instead, we explicitly modeled the electrode effects together with classical SP contributions using a flow and transport model. The simulated data agreed overall with the observed SP signals and allowed decomposing the different signal contributions to analyze them separately. After reviewing other published experimental data, we suggest that most of them include electrode effects that have not been properly taken into account. Our results suggest that previously presented SP theory works well when considering the modeling uncertainties presently associated with electrode effects. Additional work is warranted to not only develop suitable electrodes for laboratory experiments but also to assure that associated electrode effects that appear inevitable in longer term experiments are predictable, so that they can be incorporated into the modeling framework.

Impact of controlled changes in grain size and pore space characteristics on the hydraulic conductivity and spectral induced polarization response of "proxies" of saturated alluvial sediments

Understanding the influence of pore space characteristics on the hydraulic conductivity and spectral induced polarization (SIP) response is critical for establishing relationships between the electrical and hydrological properties of surficial unconsolidated sedimentary deposits, which host the bulk of the world's readily accessible groundwater resources. Here, we present the results of laboratory SIP measurements on industrial-grade, saturated quartz samples with granulometric characteristics ranging from fine sand to fine gravel, which can be regarded as proxies for widespread alluvial deposits. We altered the pore space characteristics by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples. In general, the results indicate a clear connection between the SIP response and the granulometric as well as pore space characteristics. In particular, we observe a systematic correlation between the hydraulic conductivity and the relaxation time of the Cole-Cole model describing the observed SIP effect for the entire range of considered grain sizes. The results do, however, also indicate that the detailed nature of these relations depends strongly on variations in the pore space characteristics, such as, for example, the degree of compaction. The results of this study underline the complexity of the origin of the SIP signal as well as the difficulty to relate it to a single structural factor of a studied sample, and hence raise some fundamental questions with regard to the practical use of SIP measurements as site- and/or sample-independent predictors of the hydraulic conductivity.

Erosion and channel change as factors of landslides and valley formation in Champlain sea clays: The Chacoura river, Quebec, Canada

The Champlain Sea clays of Eastern Canada are incised by numerous rivers. Their slopes have been modified by landslides: on the Chacoura River near Trois-Rivières (Quebec), several large landslide scars, more or less recent, are visible. The role of erosion (channel incision, lateral channel migration and erosion of slopes due to agricultural drainage) as a trigger of these landslides is important. The aim of this study is to understand how erosion and landslides are related to valley development. From a detailed analysis of aerial photographs and DEMs, a map of the phenomena has been drawn by identifying various elements such as landslides, limits of the slope, position of the channel, and the area covered by forest. It is shown that channel change and erosion are strongly linked to landslides by the fact that they change the bank morphology in an unstable way. A slide in itself is a natural way for the slope to achieve stability. But when it occurs in a stream, it creates a disturbance to the stream flow enhancing local erosion which may change the river path and generate more erosion downstream or upstream resulting in more slides. Cross-valley sections and a longitudinal profile show that landslides are a major factor of valley formation. It appears that the upper part of the Chacoura River valley is still unaffected by landslides and has V-shaped sections. The lower part has been subject to intense erosion and many landslide scars can be seen. This shows that the valley morphology is transient, and that future activity is more likely to occur in the upper part of the river. Therefore the identification of areas prone to erosion will help determine the possible location of future large landslides just like the ones that occurred in the lower part.

Application des méthodes géophysiques à l'étude de deux sites (Kappelen et Grenchen) de l'aquifère poreux complexe du Seeland

L'aquifère du Seeland représente une richesse en ressources hydriques qu'il est impératif de préserver contre tout risque de détérioration. Cet aquifère prolifique est constitué principalement de sédiments alluviaux post-glaciaires (graviers, sables et limons). Il est soumis aux contraintes environnementales des pratiques d'agriculture intensive, du réseau routier, des villes et de leurs activités industrielles. La connaissance optimale de ces ressources est donc primordiale pour leur protection. Dans cette optique, deux sites Kappelen et Grenchen représentatifs de l'aquifère du Seeland ont été étudiés. L'objectif de ce travail est de caractériser d'un point de vue hydrogéophysique l'aquifère au niveau de ces deux sites, c'est-à-dire, comprendre la dynamique des écoulements souterrains par l'application des méthodes électriques de surface associées aux diagraphies en intégrant des méthodes hydrogéologiques. Pour le site de Kappelen, les méthodes électriques de surface ont permis d'identifier les différents faciès géoélectriques en présence et de mettre en évidence leur disposition en une structure tabulaire et horizontale. Il s'agit d'un aquifère libre constitué d'une série de graviers allant jusqu'à 15 m de profondeur reposant sur de la moraine argileuse. Les diagraphies électriques, nucléaires et du fluide ont servis à la détermination des caractéristiques pétrophysiques et hydrauliques de l'aquifère qui contrôlent son comportement hydrodynamique. Les graviers aquifères de Kappelen présentent deux minéraux dominants: quartz et calcite. Les analyses minéralogiques indiquent que ces deux éléments constituent 65 à 75% de la matrice. La porosité totale obtenue par les diagraphies nucléaires varie de 20 à 30 %, et de 22 à 29 % par diagraphies électrique. Avec les faibles valeurs de Gamma Ray ces résultats indiquent que l'aquifère des graviers de Kappelen est dépourvu d'argile minéralogique. La perméabilité obtenue par diagraphies du fluide varie de 3.10-4 à 5.10-2 m/s, et par essais de pompage de 10-4 à 10-2 m/s. Les résultats des analyses granulométriques indiquent une hétérogénéité granulométrique au niveau des graviers aquifères. La fraction de sables, sables très fins, silts et limons constitue de 10 à 40 %. Ces éléments jouent un rôle important dans le comportement hydraulique de l'aquifère. La porosité efficace de 11 à 25% estimée à partir des résultats des analyses granulométriques suppose que les zones les plus perméables correspondent aux zones les plus graveleuses du site. Etablie sur le site de Kappelen, cette méthodologie a été utilisée sur le site de Grenchen. Les méthodes électriques de surface indiquent que l'aquifère captif de Grenchen est constitué des sables silteux comprenant des passages sableux, encadrés par des silts argileux imperméables. L'aquifère de Grenchen est disposé dans une structure relativement tabulaire et horizontale. Son épaisseur totale peut atteindre les 25 m vers le sud et le sud ouest ou les passages sableux sont les plus importants. La détermination des caractéristiques pétrophysiques et hydrauliques s'est faite à l'aide des diagraphies. Les intensités Gamma Ray varient de 30 à 100 cps, les plus fortes valeurs n'indiquent qu'une présence d'éléments argileux mais pas de bancs d'argile. Les porosités totales de 15 à 25% et les densités globales de 2.25 à 2.45 g/cm3 indiquent que la phase minérale (matrice) est composée essentiellement de quartz et de calcaire. Les densités de matrice varient entre 2.65 et 2.75 g/cm3. La perméabilité varie de 2 10-6 à 5 10-4 m/s. La surestimation des porosités totales à partir des diagraphies électriques de 25 à 42% est due à la présence d'argiles. -- The vast alluvial Seeland aquifer system in northwestern Switzerland is subjected to environmental challenges due to intensive agriculture, roads, cities and industrial activities. Optimal knowledge of the hydrological resources of this aquifer system is therefore important for their protection. Two representative sites, Kappelen and Grenchen, of the Seeland aquifer were investigated using surface-based geoelectric methods and geophysical borehole logging methods. By integrating of hydrogeological and hydrogeophysical methods, a reliable characterization of the aquifer system at these two sites can be performed in order to better understand the governing flow and transport process. At the Kappelen site, surface-based geoelectric methods allowed to identify various geoelectric facies and highlighted their tabular and horizontal structure. It is an unconfined aquifer made up of 15 m thick gravels with an important sandy fraction and bounded by a shaly glacial aquitard. Electrical and nuclear logging measurements allow for constraining the petrophysical and hydrological parameters of saturated gravels. Results indicate that in agreement with mineralogical analyses, matrix of the probed formations is dominated by quartz and calcite with densities of 2.65 and 2.71 g/cc, respectively. These two minerals constitute approximately 65 to 75 % of the mineral matrix. Matrix density values vary from 2.60 to 2.75 g/cc. Total porosity values obtained from nuclear logs range from 20 to 30 % and are consistent with those obtained from electrical logs ranging from 22 to 29 %. Together with the inherently low natural gamma radiation and the matrix density values obtained from other nuclear logging measurements, this indicates that at Kappelen site the aquifer is essentially devoid of clay. Hydraulic conductivity values obtained by the Dilution Technique vary between 3.10-4 and 5.10-2 m/s, while pumping tests give values ranging from 10-4 to 10-2 m/s. Grain size analysis of gravel samples collected from boreholes cores reveal significant granulometric heterogeneity of these deposits. Calculations based on these granulometric data have shown that the sand-, silt- and shale-sized fractions constitute between 10 and 40 % of the sample mass. The presence of these fine elements in general and their spatial distribution in particular are important as they largely control the distribution of the total and effective porosity as well as the hydraulic conductivity. Effective porosity values ranging from 11 to 25% estimated from grain size analyses indicate that the zones of higher hydraulic conductivity values correspond to the zones dominated by gravels. The methodology established at the Kappelen site was then applied to the Grenchen site. Results from surface-based geoelectric measurements indicate that it is a confined aquifer made up predominantly of shaly sands with intercalated sand lenses confined impermeable shally clay. The Grenchen confined aquifer has a relatively tabular and horizontal structure with a maximum thickness of 25 m in the south and the southwest with important sand passages. Petrophysical and hydrological characteristics were performed using electrical and nuclear logging. Natural gamma radiation values ranging from 30 to 100 cps indicate presence of a clay fraction but not of pure clay layers. Total porosity values obtained from electrical logs vary form 25 to 42%, whereas those obtained from nuclear logs values vary from 15 to 25%. This over-estimation confirms presences of clays. Density values obtained from nuclear logs varying from 2.25 to 2.45 g/cc in conjunction with the total porosity values indicate that the dominating matrix minerals are quartz and calcite. Matrix density values vary between 2.65 and 2.75 g/cc. Hydraulic conductivity values obtained by the Dilution Technique vary from 2 10-6 to 5 10-4 m/s.