Imaging and quantifying salt-tracer transport in a riparian groundwater system by means of 3D ERT monitoring

Determining groundwater flow paths of infiltrated river water is necessary for studying biochemical processes in the riparian zone, but their characterization is complicated by strong temporal and spatial heterogeneity. We investigated to what extent repeat 3D surface electrical resistance tomography (ERT) can be used to monitor transport of a salt-tracer plume under close to natural gradient conditions. The aim is to estimate groundwater flow velocities and pathways at a site located within a riparian groundwater system adjacent to the perialpine Thur River in northeastern Switzerland. Our ERT time-lapse images provide constraints on the plume's shape, flow direction, and velocity. These images allow the movement of the plume to be followed for 35 m. Although the hydraulic gradient is only 1.43 parts per thousand, the ERT time-lapse images demonstrate that the plume's center of mass and its front propagate with velocities of 2x10(-4) m/s and 5x10(-4) m/s, respectively. These velocities are compatible with groundwater resistivity monitoring data in two observation wells 5 m from the injection well. Five additional sensors in the 5-30 m distance range did not detect the plume. Comparison of the ERT time-lapse images with a groundwater transport model and time-lapse inversions of synthetic ERT data indicate that the movement of the plume can be described for the first 6 h after injection by a uniform transport model. Subsurface heterogeneity causes a change of the plume's direction and velocity at later times. Our results demonstrate the effectiveness of using time-lapse 3D surface ERT to monitor flow pathways in a challenging perialpine environment over larger scales than is practically possible with crosshole 3D ERT.

Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data

Time-lapse geophysical monitoring and inversion are valuable tools in hydrogeology for monitoring changes in the subsurface due to natural and forced (tracer) dynamics. However, the resulting models may suffer from insufficient resolution, which leads to underestimated variability and poor mass recovery. Structural joint inversion using cross-gradient constraints can provide higher-resolution models compared with individual inversions and we present the first application to time-lapse data. The results from a synthetic and field vadose zone water tracer injection experiment show that joint 3-D time-lapse inversion of crosshole electrical resistance tomography (ERT) and ground penetrating radar (GPR) traveltime data significantly improve the imaged characteristics of the point injected plume, such as lateral spreading and center of mass, as well as the overall consistency between models. The joint inversion method appears to work well for cases when one hydrological state variable (in this case moisture content) controls the time-lapse response of both geophysical methods. Citation: Doetsch, J., N. Linde, and A. Binley (2010), Structural joint inversion of time-lapse crosshole ERT and GPR traveltime data, Geophys. Res. Lett., 37, L24404, doi: 10.1029/2010GL045482.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.

Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

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.

3D geoelectrical monitoring to follow hydric soil behaviour during rainfall

RESUME Durant les dernières années, les méthodes électriques ont souvent été utilisées pour l'investigation des structures de subsurface. L'imagerie électrique (Electrical Resistivity Tomography, ERT) est une technique de prospection non-invasive et spatialement intégrée. La méthode ERT a subi des améliorations significatives avec le développement de nouveaux algorithmes d'inversion et le perfectionnement des techniques d'acquisition. La technologie multicanale et les ordinateurs de dernière génération permettent la collecte et le traitement de données en quelques heures. Les domaines d'application sont nombreux et divers: géologie et hydrogéologie, génie civil et géotechnique, archéologie et études environnementales. En particulier, les méthodes électriques sont souvent employées dans l'étude hydrologique de la zone vadose. Le but de ce travail est le développement d'un système de monitorage 3D automatique, non- invasif, fiable, peu coûteux, basé sur une technique multicanale et approprié pour suivre les variations de résistivité électrique dans le sous-sol lors d'événements pluvieux. En raison des limitations techniques et afin d'éviter toute perturbation physique dans la subsurface, ce dispositif de mesure emploie une installation non-conventionnelle, où toutes les électrodes de courant sont placées au bord de la zone d'étude. Le dispositif le plus approprié pour suivre les variations verticales et latérales de la résistivité électrique à partir d'une installation permanente a été choisi à l'aide de modélisations numériques. Les résultats démontrent que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle et plus apte à détecter les variations latérales et verticales de la résistivité électrique, et cela malgré la configuration non-conventionnelle des électrodes. Pour tester l'efficacité du système proposé, des données de terrain ont été collectées sur un site d'étude expérimental. La technique de monitorage utilisée permet de suivre le processus d'infiltration 3D pendant des événements pluvieux. Une bonne corrélation est observée entre les résultats de modélisation numérique et les données de terrain, confirmant par ailleurs que le dispositif pôle-dipôle offre une meilleure résolution que le dispositif pôle-pôle. La nouvelle technique de monitorage 3D de résistivité électrique permet de caractériser les zones d'écoulement préférentiel et de caractériser le rôle de la lithologie et de la pédologie de manière quantitative dans les processus hydrologiques responsables d'écoulement de crue. ABSTRACT During the last years, electrical methods were often used for the investigation of subsurface structures. Electrical resistivity tomography (ERT) has been reported to be a useful non-invasive and spatially integrative prospecting technique. The ERT method provides significant improvements, with the developments of new inversion algorithms, and the increasing efficiency of data collection techniques. Multichannel technology and powerful computers allow collecting and processing resistivity data within few hours. Application domains are numerous and varied: geology and hydrogeology, civil engineering and geotechnics, archaeology and environmental studies. In particular, electrical methods are commonly used in hydrological studies of the vadose zone. The aim of this study was to develop a multichannel, automatic, non-invasive, reliable and inexpensive 3D monitoring system designed to follow electrical resistivity variations in soil during rainfall. Because of technical limitations and in order to not disturb the subsurface, the proposed measurement device uses a non-conventional electrode set-up, where all the current electrodes are located near the edges of the survey grid. Using numerical modelling, the most appropriate arrays were selected to detect vertical and lateral variations of the electrical resistivity in the framework of a permanent surveying installation system. The results show that a pole-dipole array has a better resolution than a pole-pole array and can successfully follow vertical and lateral resistivity variations despite the non-conventional electrode configuration used. Field data are then collected at a test site to assess the efficiency of the proposed monitoring technique. The system allows following the 3D infiltration processes during a rainfall event. A good correlation between the results of numerical modelling and field data results can be observed since the field pole-dipole data give a better resolution image than the pole-pole data. The new device and technique makes it possible to better characterize the zones of preferential flow and to quantify the role of lithology and pedology in flood- generating hydrological processes.

The role of ceramide trihexoside (globotriaosylceramide) in the diagnosis and follow-up of the efficacy of treatment of Fabry disease: a review of the literature.

Fabry disease is caused by a deficiency of a-galactosidase A which leads to the progressive intra-lysosomal accumulation of ceramide trihexoside (CTH), also known as globotriaosylceramide (Gb3), in different cell types and body fluids. The clinical manifestations are multisystemic and predominantly affect the heart, kidney and central nervous system. The role of CTH in the pathophysiological process of Fabry disease is not established, and the link between the degree of accumulation and disease manifestations is not systematic. The use of CTH as a diagnostic tool has been proposed for several decades. The recent introduction of a specific treatment for Fabry disease in the form of enzyme replacement therapy (ERT) has led to the need for a biological marker, in place of a clinical sign, for evaluating the efficacy of treatment and also as a tool for following the long term effects of treatment. The ideal biomarker must adhere to strict criteria, and there should be a correlation between the degree of clinical efficacy of treatment and a change in its concentration. This review of the literature assesses the utility of CTH as a diagnostic tool and as a marker of the efficacy of ERT in patients with Fabry disease. Several techniques have been developed for measuring CTH; the principles and the sensitivity thresholds of these methods and the units used to express the results should be taken into consideration when interpreting data. The use of CTH measurement in Fabry disease should be re-evaluated in light of recent published data.

Quantitative integration of geophysical and hydraulic data : from the local towards the regional scale range

Des progrès significatifs ont été réalisés dans le domaine de l'intégration quantitative des données géophysique et hydrologique l'échelle locale. Cependant, l'extension à de plus grandes échelles des approches correspondantes constitue encore un défi majeur. Il est néanmoins extrêmement important de relever ce défi pour développer des modèles fiables de flux des eaux souterraines et de transport de contaminant. Pour résoudre ce problème, j'ai développé une technique d'intégration des données hydrogéophysiques basée sur une procédure bayésienne de simulation séquentielle en deux étapes. Cette procédure vise des problèmes à plus grande échelle. L'objectif est de simuler la distribution d'un paramètre hydraulique cible à partir, d'une part, de mesures d'un paramètre géophysique pertinent qui couvrent l'espace de manière exhaustive, mais avec une faible résolution (spatiale) et, d'autre part, de mesures locales de très haute résolution des mêmes paramètres géophysique et hydraulique. Pour cela, mon algorithme lie dans un premier temps les données géophysiques de faible et de haute résolution à travers une procédure de réduction déchelle. Les données géophysiques régionales réduites sont ensuite reliées au champ du paramètre hydraulique à haute résolution. J'illustre d'abord l'application de cette nouvelle approche dintégration des données à une base de données synthétiques réaliste. Celle-ci est constituée de mesures de conductivité hydraulique et électrique de haute résolution réalisées dans les mêmes forages ainsi que destimations des conductivités électriques obtenues à partir de mesures de tomographic de résistivité électrique (ERT) sur l'ensemble de l'espace. Ces dernières mesures ont une faible résolution spatiale. La viabilité globale de cette méthode est testée en effectuant les simulations de flux et de transport au travers du modèle original du champ de conductivité hydraulique ainsi que du modèle simulé. Les simulations sont alors comparées. Les résultats obtenus indiquent que la procédure dintégration des données proposée permet d'obtenir des estimations de la conductivité en adéquation avec la structure à grande échelle ainsi que des predictions fiables des caractéristiques de transports sur des distances de moyenne à grande échelle. Les résultats correspondant au scénario de terrain indiquent que l'approche d'intégration des données nouvellement mise au point est capable d'appréhender correctement les hétérogénéitées à petite échelle aussi bien que les tendances à gande échelle du champ hydraulique prévalent. Les résultats montrent également une flexibilté remarquable et une robustesse de cette nouvelle approche dintégration des données. De ce fait, elle est susceptible d'être appliquée à un large éventail de données géophysiques et hydrologiques, à toutes les gammes déchelles. Dans la deuxième partie de ma thèse, j'évalue en détail la viabilité du réechantillonnage geostatique séquentiel comme mécanisme de proposition pour les méthodes Markov Chain Monte Carlo (MCMC) appliquées à des probmes inverses géophysiques et hydrologiques de grande dimension . L'objectif est de permettre une quantification plus précise et plus réaliste des incertitudes associées aux modèles obtenus. En considérant une série dexemples de tomographic radar puits à puits, j'étudie deux classes de stratégies de rééchantillonnage spatial en considérant leur habilité à générer efficacement et précisément des réalisations de la distribution postérieure bayésienne. Les résultats obtenus montrent que, malgré sa popularité, le réechantillonnage séquentiel est plutôt inefficace à générer des échantillons postérieurs indépendants pour des études de cas synthétiques réalistes, notamment pour le cas assez communs et importants où il existe de fortes corrélations spatiales entre le modèle et les paramètres. Pour résoudre ce problème, j'ai développé un nouvelle approche de perturbation basée sur une déformation progressive. Cette approche est flexible en ce qui concerne le nombre de paramètres du modèle et lintensité de la perturbation. Par rapport au rééchantillonage séquentiel, cette nouvelle approche s'avère être très efficace pour diminuer le nombre requis d'itérations pour générer des échantillons indépendants à partir de la distribution postérieure bayésienne. - Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending corresponding approaches beyond the local scale still represents a major challenge, yet is critically important for the development of reliable groundwater flow and contaminant transport models. To address this issue, I have developed a hydrogeophysical data integration technique based on a two-step Bayesian sequential simulation procedure that is specifically targeted towards larger-scale problems. The objective is to simulate the distribution of a target hydraulic parameter based on spatially exhaustive, but poorly resolved, measurements of a pertinent geophysical parameter and locally highly resolved, but spatially sparse, measurements of the considered geophysical and hydraulic parameters. To this end, my algorithm links the low- and high-resolution geophysical data via a downscaling procedure before relating the downscaled regional-scale geophysical data to the high-resolution hydraulic parameter field. I first illustrate the application of this novel data integration approach to a realistic synthetic database consisting of collocated high-resolution borehole measurements of the hydraulic and electrical conductivities and spatially exhaustive, low-resolution electrical conductivity estimates obtained from electrical resistivity tomography (ERT). The overall viability of this method is tested and verified by performing and comparing flow and transport simulations through the original and simulated hydraulic conductivity fields. The corresponding results indicate that the proposed data integration procedure does indeed allow for obtaining faithful estimates of the larger-scale hydraulic conductivity structure and reliable predictions of the transport characteristics over medium- to regional-scale distances. The approach is then applied to a corresponding field scenario consisting of collocated high- resolution measurements of the electrical conductivity, as measured using a cone penetrometer testing (CPT) system, and the hydraulic conductivity, as estimated from electromagnetic flowmeter and slug test measurements, in combination with spatially exhaustive low-resolution electrical conductivity estimates obtained from surface-based electrical resistivity tomography (ERT). The corresponding results indicate that the newly developed data integration approach is indeed capable of adequately capturing both the small-scale heterogeneity as well as the larger-scale trend of the prevailing hydraulic conductivity field. The results also indicate that this novel data integration approach is remarkably flexible and robust and hence can be expected to be applicable to a wide range of geophysical and hydrological data at all scale ranges. In the second part of my thesis, I evaluate in detail the viability of sequential geostatistical resampling as a proposal mechanism for Markov Chain Monte Carlo (MCMC) methods applied to high-dimensional geophysical and hydrological inverse problems in order to allow for a more accurate and realistic quantification of the uncertainty associated with the thus inferred models. Focusing on a series of pertinent crosshole georadar tomographic examples, I investigated two classes of geostatistical resampling strategies with regard to their ability to efficiently and accurately generate independent realizations from the Bayesian posterior distribution. The corresponding results indicate that, despite its popularity, sequential resampling is rather inefficient at drawing independent posterior samples for realistic synthetic case studies, notably for the practically common and important scenario of pronounced spatial correlation between model parameters. To address this issue, I have developed a new gradual-deformation-based perturbation approach, which is flexible with regard to the number of model parameters as well as the perturbation strength. Compared to sequential resampling, this newly proposed approach was proven to be highly effective in decreasing the number of iterations required for drawing independent samples from the Bayesian posterior distribution.

Multiple-scale hydrological and geophysical data integration through non-linear Bayesian sequential simulation

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the scale of a field site represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed downscaling procedure based on a non-linear Bayesian sequential simulation approach. The main objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity logged at collocated wells and surface resistivity measurements, which are available throughout the studied site. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariatekernel density function. Then a stochastic integration of low-resolution, large-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities is applied. The overall viability of this downscaling approach is tested and validated by comparing flow and transport simulation through the original and the upscaled hydraulic conductivity fields. Our results indicate that the proposed procedure allows obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Kidney transplantation in patients with Fabry disease

Introduction and Aims: Fabry disease is an X-linked lysosomal storage disorder caused by absence or deficient activity of the lysosomal enzyme alpha-galactosidase A. Renal manifestations occur early in life in a significant proportion of children, in many women and in almost all men with Fabry disease. These manifestations ultimately progress to end-stage renal disease in nearly all males and in some female patients. Data on kidney transplantation in patients with Fabry disease who are receiving enzyme replacement therapy (ERT), however, are scarce. Methods: We examined the clinical characteristics of kidney transplant recipients (KTRs) in the Fabry Outcome Survey (FOS) - a European database of patients with Fabry disease that was established to monitor the safety and outcome of ERT. Results: Of the 752 patients enrolled in FOS up to October 2005, 34 (4.5%) were reported to be KTRs. The mean age of these 32 male and 2 female patients was 45 ± 9 years, the median time since the transplant was 9 years, the median estimated glomerular filtration rate (eGFR) was 46 mL/min/1.73 m2 and the median level of proteinuria was 180 mg/24 hours. ERT was well tolerated, with mild infusion-related reactions reported in only one patient. Amongst these patients, 53% were reported to have hypertension, 71% left ventricular hypertrophy, 27% cardiac valve disease and 27% arrhythmia. A total of 23 (68%) of the patients (1 female, 22 males) were receiving ERT with agalsidase alfa (Replagal; Shire Human Genetic Therapies, UK), with a median duration of treatment of 2.5 years. There were no differences in age or time since transplantation between treated and untreated patients. The median eGFRs were 46 and 49 mL/min/1.73 m2 and the median levels of proteinuria were 200 and 160 mg/24 hours, respectively. Conclusions: KTRs represent a significant minority of individuals enrolled in a large international registry of patients with Fabry disease (FOS). Approximately two-thirds of KTRs with Fabry disease enrolled in FOS receive ERT with agalsidase alfa, which is well tolerated. Comparison of treated and untreated patients has the potential to examine effects of ERT on the progression of renal and cardiovascular disease.

SOURCE, FATE AND FUNCTION OF EARLY GLIAL CELLS EXPRESSING NKX2.1 IN MOUSE EMBRYONIC BRAINS

The brain tissue is made of neuronal and glial cells generated in the germinal layer bordering the ventricles. These cells divide, differentiate and migrate following specific pathways. The specification of GABAergic interneurons and glutamatergic neurons has been broadly studied but little is known about the origin, the fate and the function of early glial cells in the embryonic telencephalon. It has been commonly accepted since long that the glial cells and more particularly the astrocytes were generated after neurogenesis from the dorsal telencephalon. However, our work shows that, unlike what was previously thought, numerous glial cells (astroglia and polydendrocytes) are generated during neurogenesis in the early embryonic stages from E14.5 to E16.5, and originate from the ventral Nkx2.1-expressing precursors instead. NK2 homeobox 1 (Nkx2.1) is a member of the NK2 family of homeodomaincontaining transcription factors. The specification of the MGE precursors requires the expression of the Nkx2.1 homeobox gene. Moreover, Nkx2.1 is previously known to regulate the specification of GABAergic interneurons and early oligodendrocytes in the ventral telencephalon. Here, in my thesis work, I have discovered that, in addition, Nkx2.1 also regulates astroglia and polydendrocytes differentiation. The use of Nkx2.1 antibody and Nkx2.1 riboprobe have revealed the presence of numerous Nkx2.1-positive cells that express astroglial markers (like GLAST and GFAP) in the entire embryonic brain. Thus, to selectively fate map MGE-derived GABAergic interneurons and glia, we crossed Nkx2.1-Cre mice, Glast-Cre ERT+/- inducible mice and NG2-Cre mice with the Cre reporter Rosa26-lox-STOP-lox-YFP (Rosa26-YFP) mice. The precise origin of Nkx2.1-positive astroglia has been directly ascertained by combining glial immunostaining and focal electroporation of the pCAG-GS-EGFP plasmids into the subpallial domains of organotypic slices, as well as, by using in vitro neurosphere experiments and in utero electroporation of the pCAG-GS-tomato plasmid into the ventral pallium of E14.5 Nkx2.1-Cre+/Rosa-YFP+/- embryos. We have, thus, confirmed that the three germinal regions of the ventral telencephalon i.e. the MGE, the AEP/POA and the triangular septal nucleus are able to generate early astroglial cells. Moreover, immunohistochemistry for several astroglial cells and polydendrocyte markers, both in the Nkx2.1-/- and control embryos and in the neurospheres, has revealed a severe loss of both glial cell types in the Nkx2.1 mutants. We found that the loss of glia corresponded to a decrease of Nkx2.1-derived precursor division capacity and glial differentiation. There was a drastic decrease of BrdU+ dividing cells labeled for Nkx2.1 in the MGE*, the POA* and the septal nucleus* of Nkx2.1 mutants. In addition, we noticed that while some remaining Nkx2.1+ precursors still succeeded to give rise to post-mitotic neurons in vitro and in vivo in the Nkx2.1-/-, they completely lost the capacity to differentiate in astrocytes. Altogether, these observations indicate for the first time that the transcription factor Nkx2.1 regulates the proliferation and differentiation of precursors in three subpallial domains that generate early embryonic astroglia and polydendrocytes. Furthermore, in order to investigate the potential function of these early Nkx2.1- derived glia, we have performed multiple immunohistochemical stainings on Nkx2.1-/- and wild-type animals, and Nkx2.1-Cre mice that were crossed to Rosa-DTA+/- mice in which the highly toxic diphtheria toxin aided to selectively deplete a majority of the Nkx2.1-derived cells. Interestingly, in these two mutants, we observed a drastic and significant loss of GFAP+, GLAST+, NG2+ and S100ß+ astroglial cells at the telencephalic midline and in the medial cortical areas. This cells loss could be directly correlated with severe axonal guidance defects observed in the corpus callosum (CC), the hippocampal commissure (HIC), the fornix (F) and the anterior commissure (AC). Axonal guidance is a key step allowing neurons to form specific connections and to become organized in a functional network. The contribution of guidepost cells inside the CC and the AC in mediating the growth of commissural axons have until now been attributed to specialized midline guidepost astroglia. Previous published results in our group have unravelled that, during embryonic development, the CC is populated in addition to astroglia by numerous glutamatergic and GABAergic guidepost neurons that are essential for the correct midline crossing of callosal axons. Therefore, the relative contribution of individual neuronal or glial populations towards the guidance of commissural axons remains largely to be investigated to understand guidance mechanisms further. Thus, we crossed Nkx2.1-Cre mice with NSE-DTA+/- mice that express the diphtheria toxin only in neurons and allowed us to selectively deplete Nkx2.1-derived GABAergic neurons. Interestingly, in the Nkx2.1-/- mice, the CC midline was totally disorganized and the callosal axons partly lost their orientation, whereas in the Nkx2.1Cre+/Rosa-DTA+/- and the Nkx2.1Cre+/NSE-DTA+/- mice, the axonal organization of the CC was not affected. In the three types of mice, hippocampal axons of the fornix were not properly fasciculated and formed disoriented bundles through the septum. Additionally, the AC formation was completely absent in Nkx2.1-/- mice and the AC was divided into two/three separate paths in the Nkx2.1Cre+/Rosa-DTA+/- mice that project in wrong territories. On the other hand, the AC didn't form or was reduced to a relatively narrower tract in the Nkx2.1Cre+/NSE-DTA+/- mice as compared to wild-type AC. These results clearly indicate that midline Nkx2.1-derived cells play a major role in commissural axons pathfinding and that both Nkx2.1-derived guidepost neurons and glia are necessary elements for the correct development of these commissures. Furthermore, during our investigations on Nkx2.1-/- and Nkx2.1Cre+/Rosa-DTA+/- mice, we noticed similar and severe defects in the erythrocytes distribution and the blood vessels network morphology in the embryonic brain of both mutants. As the Cre-mediated recombination was never observed to occur in the blood vessels of Nkx2.1-Cre mice, we inferred that the vessels defects observed were due to the loss of Nkx2.1-derived cells and not to the cells autonomous effects of Nkx2.1 in regulating endothelial cell precursors. Thereafter, the respective contribution of individual Nkx2.1-regulated neuronal or glial populations in the blood vessels network building were studied with the use of transgenic mice strains. Indeed, the use of Nkx2.1Cre+/NSE-DTA+/- mice indicated that the Nkx2.1-derived neurons were not implicated in this process. Finally, to discriminate between the two Nkx2.1-derived glial cell populations, the GLAST+ astroglia and the NG2+ polydendrocytes, an NG2-Cre mouse strain crossed to the Rosa-DTA+/- mice was used. In that mutant, the blood vessel network and the erythrocytes distribution were similarly affected as observed in Nkx2.1Cre+/Rosa-DTA+/- animals. Therefore, this result indicates that most probably, the NG2+ polydendrocytes are involved in helping to build the vessels network in the brain. Taken altogether, these observations show that during brain development, Nkx2.1- derived embryonic glial cells act as guidepost cells on the guidance of axons as well as forming vessels. Both Nkx2.1-regulated guidepost GABAergic neurons and glia collaborate to guide growing commissural axons, while polydendrocytes are implicated in regulating brain angiogenesis. - Le tissu cérébral est composé de cellules neuronales et gliales générées dans les couches germinales qui bordent les ventricules. Ces cellules se divisent, se différencient et migrent selon des voies particulières. La spécification des interneurones GABAergiques et des neurones glutamatergiques a été largement étudiée, par contre, l'origine, le destin et la fonction des cellules gliales précoces du télencéphale embryonnaire restent peu élucidées. Depuis longtemps, il était communément accepté que les cellules gliales, et plus particulièrement les astrocytes, sont générés après la neurogénèse à partir du télencéphale dorsal. Toutefois, notre travail montre que de nombreuses cellules gliales sont générées à partir de précurseurs ventraux qui expriment le gène Nkx2.1, entre E14.5 et E16.5, c'est-à dire,à des stades embryonnaires très précoces. Le gène NK2 homéobox 1 (Nkx2.1) appartient à une famille de facteurs de transcription appelée NK2. Il s'agit de protéines qui contiennent un homéo-domaine. La spécification des précurseurs de la MGE requiert l'expression du gène homéobox Nkx2.1. De plus, la fonction du gène Nkx2.1 dans la régulation de la spécification des interneurones GABAergiques et des oligodendrocytes dans le télencéphale ventral était déjà connue. Au cours de mon travail de thèse, j'ai également mis en évidence que, Nkx2.1 régule aussi les étapes de prolifération et de différenciation de divers sous-types de cellules gliales soit de type astrocytes ou bien polydendrocytes. L'utilisation d'un anticorps contre la protéine Nkx2.1 ainsi qu'une sonde à ribonucléotides contre l'ARN messager du gène Nkx2.1 ont révélé la présence de nombreuses cellules positives pour Nkx2.1 qui exprimaient des marqueurs astrocytaires (comme GLAST et GFAP) dans le télencéphale embryonnaire. Afin de déterminer de manière sélective le sort des interneurones GABAergiques, des polydendrocytes et des astrocytes dérivés de la MGE, nous avons croisé soit des souris Nkx2.1-Cre, des souris Glast-Cre ERT+/- inductibles ou bien des souris NG2-Cre avec des souris Rosa26-lox-STOP-lox-YFP (Rosa26-YFP) Cre rapportrices. L'origine précise des astroglies positives pour Nkx2.1 a été directement établie en combinant une coloration immunologique pour les glies et une électroporation focale d'un plasmide pCAG-GS-EGFP dans les domaines subpalliaux de tranches organotypiques, puis également, par des cultures de neurosphères in vitro et des expériences d'électroporation in utero d'un plasmide pCAG-GS-tomato dans le pallium ventral d'embryons Nkx2.1-Cre+/Rosa- YFP+/- au stade E14.5. Nous avons donc confirmé que les trois régions germinales du télencéphale ventral, c'est-à-dire, la MGE, l'AEP/POA et le noyau triangulaire septal sont capables de générer des cellules astrogliales. D'autre part, l'immunohistochimie pour plusieurs marqueurs d'astrocytes ou de polydendrocytes, dans les embryons Nkx2.1-/- et contrôles ainsi que dans les neurosphères, a révélé une sévère perte de ces deux types gliaux chez les mutants. Nous avons trouvé que la perte de glies correspondait à une diminution de la capacité de division des précurseurs dérivés de Nkx2.1, ainsi que l'incapacité de ces précurseurs de se différencier en cellules gliales. Nous avons en effet observé une diminution importante des cellules BrdU+ en division exprimant Nkx2.1dans la MGE*, la POA* et le noyau septal* des mutants pour Nkx2.1. D'autre part, nous avons pu mettre en évidence aussi bien in vitro, qu'in vivo, que certains précurseurs Nkx2.1+ chez le mutant gardent la capacité à se différencier en neurones tandis qu'ils perdent celle de se différencier en cellules gliales. Prises dans leur ensemble, ces observations indiquent pour la première fois que le facteur de transcription Nkx2.1 régule les étapes de prolifération et de différentiation des précurseurs des trois domaines subpalliaux qui génèrent les astroglies et polydendrocytes embryonnaires précoces. Par la suite, dans le but de comprendre la fonction potentielle de ces glies précoces, nous avons procédé à de multiples colorations immunohistochimiques sur des animaux Nkx2.1-/- et sauvages, ainsi que sur des souris Nkx2.1-Cre croisées à des souris Rosa-DTA+/- dans lesquelles la toxine diphthérique hautement toxique a permis de supprimer sélectivement la majorité des cellules dérivées de Nkx2.1. De manière intéressante, nous avons observé dans ces deux mutants, une perte drastique et significative de cellules astrogliales GFAP+, GLAST+ et polydendrocytaires NG2+ et S100ß+ dans le télencéphale, à la midline et dans les aires corticales médianes. Ces pertes ont pu être directement corrélées avec des défauts de guidage axonal observés dans le corps calleux (CC), la commissure hippocampique (HIC), le fornix (F) et la commissure antérieure (AC). Le guidage axonal est une étape clé permettant aux neurones de former des connections spécifiques et de s'organiser dans un réseau fonctionnel. La contribution des cellules « guidepost » dans le CC et dans la AC comme médiateurs de la croissance des axones commissuraux à jusqu'à aujourd'hui été attribuée spécifiquement à des astroglies « guidepost » de la midline. Des résultats publiés précédemment dans notre groupe, ont permis de montrer que, pendant le développement embryonnaire, le CC est peuplé en plus de la glie par de nombreux neurones « guidepost » glutamatergiques et GABAergiques qui sont essentiels pour le croisement correct des axones callosaux à la midline. Ainsi, la contribution relative des populations individuelles neuronales ou gliales pour le guidage des axones commissuraux demande à être approfondie afin de mieux comprendre les mécanismes de guidage. A ces fins, nous avons croisé des souris Nkx2.1-Cre avec des souris NSE-DTA+/- qui expriment la toxine diphthérique uniquement dans les neurones et ainsi, nous avons pu sélectivement supprimer les neurones dérivés de domaines Nkx2.1+. Dans les souris Nkx2.1-/-,nous avons découvert que le CC était désorganisé avec des axones callosaux perdant partiellement leur orientation, alors que dans les souris Nkx2.1Cre+/Rosa-DTA+/- et Nkx2.1Cre+/NSE-DTA+/-, l'organisation axonale n'était pas affectée. De plus, les faisceaux hippocampiques du fornix étaient défasciculés dans les trois types de mutants. Par ailleurs, la formation de la commissure antérieure (AC) était complètement absente dans les souris Nkx2.1-/- d'une part, et d'autre part, celle-ci était divisée en deux à trois voies séparées dans les souris Nkx2.1Cre+/Rosa-DTA+/-. Finalement, la AC était soit absente, soit réduite de manière ne former plus qu'un faisceau relativement plus étroit dans les souris Nkx2.1Cre+/NSE-DTA+/- en comparaison avec la AC sauvage. Ces derniers résultats indiquent clairement que les cellules dérivées de Nkx2.1 à la midline, jouent un rôle majeur dans le guidage des axones commissuraux et que, autant les neurones, que les astrocytes « guidepost » dérivés de Nkx2.1, sont des éléments nécessaires au développement correct de ces commissures. En outre, lors de nos investigations sur les souris Nkx2.1-/- et Nkx2.1Cre+/Rosa-DTA+/-, nous avons remarqués des défauts sévères et similaires dans la distribution des erythrocytes et dans la morphologie du réseau de vaisseaux sanguins dans le cerveau embryonnaire des deux mutants précités. Puisque nous n'avons jamais observé de recombinaison de la Cre recombinase dans les vaisseaux sanguins des souris Nkx2.1Cre, nous en avons déduit que les défauts de vaisseaux observés étaient dus à la perte de cellules dérivées de Nkx2.1. Il existerait donc en plus de la fonction cellulaire autonome de Nkx2.1 reconnue pour régulée directement la spécification des cellules endothéliales, une fonction indirecte de Nkx2.1. Afin de déterminer la contribution respective des populations individuelles neuronales ou gliales régulées par Nkx2.1 dans la construction du réseau de vaisseaux sanguins, nous avons utilisé diverses lignées de souris transgéniques. L'utilisation de souris Nkx2.1Cre+/NSE-DTA+/- a indiqué que les neurones dérivés de Nkx2.1 n'étaient pas impliqués dans ce processus. Finalement, afin de discriminer entre les deux populations de cellules gliales dérivées de Nkx2.1, les astroglies et les polydendrocytes, nous avons croisé une lignée de souris NG2-Cre avec des souris Rosa-DTA+/-. Dans ce dernier mutant, le réseau de vaisseaux sanguins du cortex ainsi que la distribution des erythrocytes étaient affectés de la même manière que dans le cortex des souris Nkx2.1Cre+/Rosa-DTA+/-. Par conséquent, ce résultat indique que très probablement, les polydendrocytes NG2+ sont impliqués dans la mise en place du réseau de vaisseaux dans le cerveau. Prises dans leur ensemble, ces observations montrent que durant le développement embryonnaire du cerveau, des sous-populations de glies régulées par Nkx2.1 jouent un rôle de cellules « guidepost » dans le guidage des axones, ainsi que des vaisseaux. Les polydendrocytes sont impliquées dans la régulation de l'angiogenèse tandis que, autant les neurones GABAergiques que les astrocytes collaborent dans le guidage des axones commissuraux en croissance.

Clinical results of enzyme replacement therapy in Fabry disease: a comprehensive review of literature.

Enzyme replacement therapy (ERT) has been used to treat Fabry disease - a progressive lysosomal storage disorder - since 2001. Two preparations of the enzyme alpha-galactosidase A are available in Europe: agalsidase alpha, produced in a human cell line, and agalsidase beta, produced in Chinese hamster ovary cells. To review critically the published evidence for the clinical efficacy of these two enzyme preparations. A systematic literature search was undertaken to identify open or randomised controlled trials published on Fabry disease since 2001. Eleven trials fulfilled the criteria for inclusion in this review, of a total of 586 references on Fabry disease. To date, no direct comparisons exists between the two available enzyme preparations. Significant clinical benefits compared with placebo, however, have been demonstrated with ERT, with positive effects on the heart, kidneys, nervous system and quality of life. The quality of most of these publications was less than optimal. Further prospective studies are required to confirm the long-term clinical benefits of ERT. More studies are also needed on the effects of ERT in women and on the use of ERT early in the course of Fabry disease, to prevent organ damage. Large national and international outcomes databases will also be invaluable in evaluating treatment effects and safety.

Cystatin C as a marker of early changes of renal function in Fabry nephropathy.

BACKGROUND: A sensitive, feasible and reproducible marker for renal function is necessary to evaluate the clinical efficacy of enzyme replacement therapy (ERT) in Fabry nephropathy. Serum creatinine has some limitations and cystatin C has been proposed, in other nephropathies, as a useful marker of renal function. The use of cystatin C as a marker of glomerular filtration rate (GFR) was investigated in Fabry patients receiving ERT. METHODS: Renal function was evaluated with serum creatinine, serum cystatin C and estimated GFR (through Modification of Diet in Renal Disease [MDRD], Cockcroft-Gault [C&G] and Hoek formulae) in 21 Fabry patients receiving ERT with agalsidase alfa for 3 years and in 13 Fabry patients receiving agalsidase alfa for 4 years. RESULTS: During years of ERT while serum creatinine remained stable, cystatin C values showed a significant, increasing trend right from the first year of ERT. CONCLUSIONS: In Fabry disease, cystatin C is a sensitive and reliable marker of renal function, and it should be taken into account when evaluating GFR trends during ERT.

Rate and processes of river network rearrangement during incipient faulting: The case of the Cahabon river, Guatemala

Deeply incised river networks are generally regarded as robust features that are not easily modified by erosion or tectonics. Although the reorganization of deeply incised drainage systems has been documented, the corresponding importance with regard to the overall landscape evolution of mountain ranges and the factors that permit such reorganizations are poorly understood. To address this problem, we have explored the rapid drainage reorganization that affected the Cahabon River in Guatemala during the Quaternary. Sediment-provenance analysis, field mapping, and electrical resistivity tomography (ERT) imaging are used to reconstruct the geometry of the valley before the river was captured. Dating of the abandoned valley sediments by the Be-10-Al-26 burial method and geomagnetic polarity analysis allow us to determine the age of the capture events and then to quantify several processes, such as the rate of tectonic deformation of the paleovalley, the rate of propagation of post-capture drainage reversal, and the rate at which canyons that formed at the capture sites have propagated along the paleovalley. Transtensional faulting started 1 to 3 million years ago, produced ground tilting and ground faulting along the Cahabon River, and thus generated differential uplift rate of 0.3 +/- 0.1 up to 0.7 +/- 0.4 mm . y(-1) along the river's course. The river responded to faulting by incising the areas of relative uplift and depositing a few tens of meters of sediment above the areas of relative subsidence. Then, the river experienced two captures and one avulsion between 700 ky and 100 ky. The captures breached high-standing ridges that separate the Cahabon River from its captors. Captures occurred at specific points where ridges are made permeable by fault damage zones and/or soluble rocks. Groundwater flow from the Cahabon River down to its captors likely increased the erosive power of the captors thus promoting focused erosion of the ridges. Valley-fill formation and capture occurred in close temporal succession, suggesting a genetic link between the two. We suggest that the aquifers accumulated within the valley-fills, increased the head along the subterraneous system connecting the Cahabon River to its captors, and promoted their development. Upon capture, the breached valley experienced widespread drainage reversal toward the capture sites. We attribute the generalized reversal to combined effects of groundwater sapping in the valley-fill, axial drainage obstruction by lateral fans, and tectonic tilting. Drainage reversal increased the size of the captured areas by a factor of 4 to 6. At the capture sites, 500 m deep canyons have been incised into the bedrock and are propagating upstream at a rate of 3 to 11 mm . y(-1) deepening at a rate of 0.7 to 1 5 mm . y(-1). At this rate, 1 to 2 million years will be necessary for headward erosion to completely erase the topographic expression of the paleovalley. It is concluded that the rapid reorganization of this drainage system was made possible by the way the river adjusted to the new tectonic strain field, which involved transient sedimentation along the river's course. If the river had escaped its early reorganization and had been given the time necessary to reach a new dynamic equilibrium, then the transient conditions that promoted capture would have vanished and its vulnerability to capture would have been strongly reduced.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.