Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.

Identification and characterization of novel molecular mechanisms involved in initiation and progression of myelination

Schwann cells synthesize a large amount of membrane that form a specialized structure called myelin that surrounds axons and facilitate the transmission of electrical signal along neurons in peripheral nervous system (PNS). Previous studies demonstrated that both Schwann cell differentiation and de-differentiation (in the situation of a nerve injury or demyelinating disease) are regulated by cell-intrinsic regulators including several transcription factors. In particular, the de-differentiation of mature Schwann cells is driven by the activation of multiple negative regulators of myelination including Sox2, c-Jun, Notch and Pax3, all usually expressed in immature Schwann cells and suppressed at the onset of myelination. In order to identify new regulators of myelination involved in the development of the PNS, we analyzed the gene-expression profiling data from developing PNS and from three models of demyelinating neuropathies. This analysis led to the identification of Sox4, a member of the Sox family of transcription factors, as a potential candidate. To characterize the molecular function of Sox4 in PNS, we generated two transgenic lines of mice, which overexpress Sox4 specifically in Schwann cells. Detailed analysis of these mice showed that the overexpression of Sox4 in Schwann cells causes a delay in progression of myelination between post-natal day 2 (P2) and P5. Our in vitro analysis suggested that Sox4 cDNA can be overexpressed while the protein translation is tightly regulated. Interestingly, we observed that Sox4 protein is stabilized in nerves of the CMT4C mouse, a model of the human neuropathy. We therefore crossed Sox4 transgenic mice with CMT4C mice and we observed that Sox4 overexpression exacerbated the neuropathy phenotype in these mice. While recognized as being crucial for the normal function of both neurons and myelinating glial cells, the processes that regulate the beginning of myelination and the nature of the neuro-glial cross-talk remains mostly unknown. In order to gain insight into the molecular pathways involved in the interactions between neurons and associated glial cells, we developed a neuron-glia co-culture system based on microfluidic chambers and successfully induced myelination in this system by ascorbic acid. Importantly, we observed that in addition to acting on Schwann cells, ascorbic acid also modulate neuronal/axonal NRG1/ErbB2-B3 signalling. The experimental setting used in our study thus allowed us to discover a novel phenomena of propagation for myelination in vitro. The further characterization of this event brought us to identify other compounds able to induce myelination: ADAMs secretases inhibitor GM6001 and cyclic-AMP. The results generated during my thesis project are therefore not only important for the advancement of our understanding of how the PNS works, but may also potentially help to develop new therapies aiming at improvement of PNS myelination under disease conditions. - Les cellules de Schwann synthétisent une grande quantité de membrane formant une structure spécialisée appelée myéline qui entoure les axones et facilite la transmission du signal électrique le long des neurones du système nerveux périphérique (SNP). Des études antérieures ont démontré que la différenciation et la dédifférenciation des cellules de Schwann (dans la situation d'une lésion nerveuse ou d'une maladie démyélinisante) sont régulées par des régulateurs cellulaires intrinsèques, incluant plusieurs facteurs de transcription. En particulier, la dédifférenciation des cellules de Schwann matures est contrôlée par l'activation de plusieurs régulateurs négatifs de la myélinisation dont Sox2, c-Jun, Notch et Pax3, tous habituellement exprimés dans des cellules de Schwann immatures et supprimés au début de la myélinisation. Afin d'identifier de nouveaux régulateurs de myélinisation impliqués dans le développement du SNP, nous avons analysé le profil d'expression génique durant le développement du SNP ainsi que dans trois modèles de neuropathies démyélinisantes. Cette analyse a mené à l'identification de Sox4, un membre de la famille des facteurs de transcription Sox, comme étant un candidat potentiel. Dans le but de caractériser la fonction moléculaire de Sox4 dans le SNP, nous avons généré deux lignées transgéniques de souris qui surexpriment Sox4 spécifiquement dans les cellules de Schwann. L'analyse détaillée de ces souris a montré que la surexpression de Sox4 dans les cellules de Schwann provoque un retard dans la progression de la myélinisation entre le jour postnatal 2 (P2) et P5. Notre analyse in vitro a suggéré que l'ADNc de Sox4 peut être surexprimé alors que la traduction des protéines est quand à elle étroitement régulée. De façon intéressante, nous avons observé que la protéine Sox4 est stabilisée dans les nerfs des souris CMT4C, un modèle de neuropathie humaine. Nous avons donc croisé les souris transgéniques Sox4 avec des souris CMT4C et avons observé que la surexpression de Sox4 exacerbe le phénotype de neuropathie chez ces souris. Bien que reconnus comme étant cruciaux pour le fonctionnement normal des neurones et des cellules gliales myélinisantes, les processus qui régulent le début de la myélinisation ainsi que la nature des interactions neurone-glie restent largement méconnus. Afin de mieux comprendre les mécanismes moléculaires impliqués dans les interactions entre les neurones et les cellules gliales leur étant associés, nous avons développé un système de co-culture neurone-glie basé sur des chambres microfluidiques et y avons induit avec succès la myélinisation avec de l'acide ascorbique. Étonnamment, nous avons remarqué que, en plus d'agir sur les cellules de Schwann, l'acide ascorbique module également la voie de signalisation neuronale/axonale NRG1/ErbB2-B3. Le protocole expérimental utilisé dans notre étude a ainsi permis de découvrir un nouveau phénomène de propagation de la myélinisation in vitro. La caractérisation plus poussée de ce phénomène nous a menés à identifier d'autres composés capables d'induire la myélinisation: L'inhibiteur de sécrétases ADAMs GM6001 et l'AMP cyclique. Les résultats obtenus au cours de mon projet de thèse ne sont donc pas seulement importants pour l'avancement de notre compréhension sur la façon dont le SNP fonctionne, mais peuvent aussi potentiellement aider à développer de nouvelles thérapies visant à l'amélioration de la myélinisation du SNP dans des conditions pathologiques.

Quantifying network properties in multi-electrode recordings: spatiotemporal characterization and inter-trial variation of evoked gamma oscillations in mouse somatosensory cortex in vitro.

Linking the structural connectivity of brain circuits to their cooperative dynamics and emergent functions is a central aim of neuroscience research. Graph theory has recently been applied to study the structure-function relationship of networks, where dynamical similarity of different nodes has been turned into a "static" functional connection. However, the capability of the brain to adapt, learn and process external stimuli requires a constant dynamical functional rewiring between circuitries and cell assemblies. Hence, we must capture the changes of network functional connectivity over time. Multi-electrode array data present a unique challenge within this framework. We study the dynamics of gamma oscillations in acute slices of the somatosensory cortex from juvenile mice recorded by planar multi-electrode arrays. Bursts of gamma oscillatory activity lasting a few hundred milliseconds could be initiated only by brief trains of electrical stimulations applied at the deepest cortical layers and simultaneously delivered at multiple locations. Local field potentials were used to study the spatio-temporal properties and the instantaneous synchronization profile of the gamma oscillatory activity, combined with current source density (CSD) analysis. Pair-wise differences in the oscillation phase were used to determine the presence of instantaneous synchronization between the different sites of the circuitry during the oscillatory period. Despite variation in the duration of the oscillatory response over successive trials, they showed a constant average power, suggesting that the rate of expenditure of energy during the gamma bursts is consistent across repeated stimulations. Within each gamma burst, the functional connectivity map reflected the columnar organization of the neocortex. Over successive trials, an apparently random rearrangement of the functional connectivity was observed, with a more stable columnar than horizontal organization. This work reveals new features of evoked gamma oscillations in developing cortex.

Geophysical Investigations for Hydrogeological Characterization at the Kappelen Test Site (Switzerland)

Vertical electric soundings, 2D resistivity imaging and several logging measurements were performed at Kappelen test site to identify the various geolelectric facies that allowed determining the tabular and horizontal structure of the aquifer. The surface-based geoelectric methods allowed for a reliable characterization of the overall structure and the geometry of the aquifer, while geophysical logging methods allowed for inferring detailed hydrogeophysical characteristics, such as the electrical resistivity, total porosity, global and matrix density and hydraulic conductivity. The synoptic interpretation and integration of this broad and diverse database allows for constraining the key hydrological characteristics and hence forms the basis for the detailed hydraulic modelling of flow and transport process.

Localization and characterization of an active fault in an urbanized area in central Guatemala by means of geoelectrical imaging

The Polochic and Motagua faults define the active plate boundary between the North American and Caribbean plates in central Guatemala. A splay of the Polochic Fault traverses the rapidly growing city of San Miguel Uspantan that is periodically affected by destructive earthquakes. This fault splay was located using a 2D electrical resistivity tomography (ERT) survey that also characterized the fault damage zone and evaluated the thickness and nature of recent deposits upon which most of the city is built. ERT images show the fault as a similar to 50 m wide, near-vertical low-resistivity anomaly, bounded within a few meters by high resistivity anomalies. Forward modeling reproduces the key aspects of the observed electrical resistivity data with remarkable fidelity thus defining the overall location, geometry, and internal structure of the fault zone as well as the affected lithologies. Our results indicate that the city is constructed on a similar to 20 m thick surficial layer consisting of poorly consolidated, highly porous, water-logged pumice. This soft layer is likely to amplify seismic waves and to liquefy upon moderate to strong ground shaking. The electrical conductivity as well as the major element chemistry of the groundwater provides evidence to suggest that the local aquifer might, at least in part, be fed by water rising along the fault. Therefore, the potential threat posed by this fault splay may not be limited to its seismic activity per se, but could be compounded its potential propensity to enhance seismic site effects by injecting water into the soft surficial sediments. The results of this study provide the basis for a rigorous analysis of seismic hazard and sustainable development of San Miguel Uspantan and illustrate the potential of ERT surveying for paleoseismic studies.

Detection and Characterization of Preferential Flow Paths in the Downstream Area of a Hazardous Landfill

We have used surface-based electrical resistivity tomography to detect and characterize preferential hydraulic pathways in the immediate downstream area of an abandoned, hazardous landfill. The landfill occupies the void left by a former gravel pit and its base is close to the groundwater table and lacking an engineered barrier. As such, this site is remarkably typical of many small- to medium-sized waste deposits throughout the densely populated and heavily industrialized foreland on both sides of the Alpine arc. Outflows of pollutants lastingly contaminated local drinking water supplies and necessitated a partial remediation in the form of a synthetic cover barrier, which is meant to prevent meteoric water from percolating through the waste before reaching the groundwater table. Any future additional isolation of the landfill in the form of lateral barriers thus requires adequate knowledge of potential preferential hydraulic pathways for outflowing contaminants. Our results, inferred from a suite of tomographically inverted surfaced-based electrical resistivity profiles oriented roughly perpendicular to the local hydraulic gradient, indicate that potential contaminant outflows would predominantly occur along an unexploited lateral extension of the original gravel deposit. This finds its expression as a distinct and laterally continuous high-resistivity anomaly in the resistivity tomograms. This interpretation is ground-truthed through a litholog from a nearby well. Since the probed glacio-fluvial deposits are largely devoid of mineralogical clay, the geometry of hydraulic and electrical pathways across the pore space of a given lithological unit can be assumed to be identical, which allows for an order-of-magnitude estimation of the overall permeability structure. These estimates indicate that the permeability of the imaged extension of the gravel body is at least two to three orders-of-magnitude higher than that of its finer-grained embedding matrix. This corroborates the preeminent role of the high-resistivity anomaly as a potential preferential flow path.

Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data

Predictive groundwater modeling requires accurate information about aquifer characteristics. Geophysical imaging is a powerful tool for delineating aquifer properties at an appropriate scale and resolution, but it suffers from problems of ambiguity. One way to overcome such limitations is to adopt a simultaneous multitechnique inversion strategy. We have developed a methodology for aquifer characterization based on structural joint inversion of multiple geophysical data sets followed by clustering to form zones and subsequent inversion for zonal parameters. Joint inversions based on cross-gradient structural constraints require less restrictive assumptions than, say, applying predefined petro-physical relationships and generally yield superior results. This approach has, for the first time, been applied to three geophysical data types in three dimensions. A classification scheme using maximum likelihood estimation is used to determine the parameters of a Gaussian mixture model that defines zonal geometries from joint-inversion tomograms. The resulting zones are used to estimate representative geophysical parameters of each zone, which are then used for field-scale petrophysical analysis. A synthetic study demonstrated how joint inversion of seismic and radar traveltimes and electrical resistance tomography (ERT) data greatly reduces misclassification of zones (down from 21.3% to 3.7%) and improves the accuracy of retrieved zonal parameters (from 1.8% to 0.3%) compared to individual inversions. We applied our scheme to a data set collected in northeastern Switzerland to delineate lithologic subunits within a gravel aquifer. The inversion models resolve three principal subhorizontal units along with some important 3D heterogeneity. Petro-physical analysis of the zonal parameters indicated approximately 30% variation in porosity within the gravel aquifer and an increasing fraction of finer sediments with depth.

Detection and characterization of hydraulically active fractures in a carbonate aquifer: Results from self-potential, temperature and fluid conductivity logging along a 260-m-deep borehole in the Combioula hydrothermal system in the southwestern Swiss Alps

A geophysical and geochemical study has been conducted in a fractured carbonate aquifer located at Combioula in the southwestern Swiss Alps with the objective to detect and characterize hydraulically active fractures along a 260-m-deep borehole. Hydrochemical analyses, borehole diameter, temperature and fluid electrical conductivity logging data were integrated in order to relate electrokinetic self-potential signals to groundwater flow inside the fracture network. The results show a generally good, albeit locally variable correlation of variations of the self-potential signals with variations in temperature, fluid electrical conductivity and borehole diameter. Together with the hydrochemical evidence, which was found to be critical for the interpretation of the self-potential data, these measurements not only made it possible to detect the hydraulically active fractures but also to characterize them as zones of fluid gain or fluid loss. The results complement the available information from the corresponding litholog and illustrate the potential of electrokinetic self-potential signals in conjunction with temperature, fluid electrical conductivity and hydrochemical analyses for the characterization of fractured aquifers, and thus may offer a perspective for an effective quantitative characterization of this increasingly important class of aquifers and geothermal reservoirs.

3D crosshole ERT for aquifer characterization and monitoring of infiltrating river water

The hydrogeological properties and responses of a productive aquifer in northeastern Switzerland are investigated. For this purpose, 3D crosshole electrical resistivity tomography (ERT) is used to define the main lithological structures within the aquifer (through static inversion) and to monitor the water infiltration from an adjacent river. During precipitation events and subsequent river flooding, the river water resistivity increases. As a consequence, the electrical characteristics of the infiltrating water can be used as a natural tracer to delineate preferential flow paths and flow velocities. The focus is primarily on the experiment installation, data collection strategy, and the structural characterization of the site and a brief overview of the ERT monitoring results. The monitoring system comprises 18 boreholes each equipped with 10 electrodes straddling the entire thickness of the gravel aquifer. A multi-channel resistivity system programmed to cycle through various four-point electrode configurations of the 180 electrodes in a rolling sequence allows for the measurement of approximately 15,500 apparent resistivity values every 7 h on a continuous basis. The 3D static ERT inversion of data acquired under stable hydrological conditions provides a base model for future time-lapse inversion studies and the means to investigate the resolving capability of our acquisition scheme. In particular, it enables definition of the main lithological structures within the aquifer. The final ERT static model delineates a relatively high-resistivity, low-porosity, intermediate-depth layer throughout the investigated aquifer volume that is consistent with results from well logging and seismic and radar tomography models. The next step will be to define and implement an appropriate time-lapse ERT inversion scheme using the river water as a natural tracer. The main challenge will be to separate the superposed time-varying effects of water table height, temperature, and salinity variations associated with the infiltrating water.

Probabilistic electrical resistivity tomography of a CO2 sequestration analog

Electrical resistivity tomography (ERT) is a well-established method for geophysical characterization and has shown potential for monitoring geologic CO2 sequestration, due to its sensitivity to electrical resistivity contrasts generated by liquid/gas saturation variability. In contrast to deterministic inversion approaches, probabilistic inversion provides the full posterior probability density function of the saturation field and accounts for the uncertainties inherent in the petrophysical parameters relating the resistivity to saturation. In this study, the data are from benchtop ERT experiments conducted during gas injection into a quasi-2D brine-saturated sand chamber with a packing that mimics a simple anticlinal geological reservoir. The saturation fields are estimated by Markov chain Monte Carlo inversion of the measured data and compared to independent saturation measurements from light transmission through the chamber. Different model parameterizations are evaluated in terms of the recovered saturation and petrophysical parameter values. The saturation field is parameterized (1) in Cartesian coordinates, (2) by means of its discrete cosine transform coefficients, and (3) by fixed saturation values in structural elements whose shape and location is assumed known or represented by an arbitrary Gaussian Bell structure. Results show that the estimated saturation fields are in overall agreement with saturations measured by light transmission, but differ strongly in terms of parameter estimates, parameter uncertainties and computational intensity. Discretization in the frequency domain (as in the discrete cosine transform parameterization) provides more accurate models at a lower computational cost compared to spatially discretized (Cartesian) models. A priori knowledge about the expected geologic structures allows for non-discretized model descriptions with markedly reduced degrees of freedom. Constraining the solutions to the known injected gas volume improved estimates of saturation and parameter values of the petrophysical relationship. (C) 2014 Elsevier B.V. All rights reserved.