Structural development of the Tso Morari ultra-high pressure nappe of the Ladakh Himalaya

A continental subduction-related and multistage exhumation process for the Tso Morari ultra-high pressure nappe is proposed. The model is constrained by published thermo-barometry and age data, combined with new geological and tectonic maps. Additionally, observations on the structural and metamorphic evolution of the Tso Morari area and the North Himalayan nappes are presented. The northern margin of the Indian continental crust was subducted to a depth of >90 km below Asia after continental collision some 55 Ma ago. The underthrusting was accompanied by the detachment and accretion of Late Proterozoic to Early Eocene sediments, creating the North Himalayan accretionary wedge, in front of the active Asian margin and the 103-50 Ma Ladakh arc batholith. The basic dikes in the Ordovician Tso Morari granite were transformed to eclogites with crystallization of coesite, some 53 Ma ago at a depth of >90 kin (>27 kbar) and temperatures of 500 to 600 degrees C. The detachment and extrusion of the low density Tso Morari nappe, composed of 70% of the Tso Morari granite and 30% of graywackes with some eclogitic dikes, occurred by ductile pure and simple shear deformation. It was pushed by buoyancy forces and by squeezing between the underthrusted Indian lithosphere and the Asian mantle wedge. The extruding Tso Morari nappe reached a depth of 35 km at the base of the North Himalayan accretionary wedge some 48 Ma ago. There the whole nappe stack recrystallized under amphibolite facies conditions of a Barrovian regional metamorphism with a metamorphic field gradient of 20 degrees C/km. An intense schistosity with a W-E oriented stretching lineation L, and top-to-the E shear criteria and crystallization of oriented sillimanite needles after kyanite, testify to the Tso Morari nappe extrusion and pressure drop. The whole nappe stack, comprising from the base to top the Tso Morari, Tetraogal, Karzok and Mata-Nyimaling-Tsarap nappes, was overprinted by new schistosities with a first N-directed and a second NE-directed stretching lineation L-2 and L-3 reaching the base of the North Himalayan accretionary wedge. They are characterized by top-to-the S and SW shear criteria. This structural overprint was related to an early N- and a younger NE-directed underthrusting of the Indian plate below Asia that was accompanied by anticlockwise rotation of India. The warping of the Tso Morari dome started already some 48 Ma ago with the formation of an extruding nappe at depth. The Tso Morari dome reached a depth of 15 km about 40 Ma ago in the eastern Kiagar La region and 30 Ma ago in the western Nuruchan region. The extrusion rate was of about 3 cm/yr between 53 and 48 Ma, followed by an uplift rate of 1.2 mm/yr between 48 and 30 Ma and of only 0.5 mm/yr after 30 Ma. Geomorphology observations show that the Tso Morari dome is still affected by faults, open regional dome, and basin and pull-apart structures, in a zone of active dextral transpression parallel to the Indus Suture zone.

Exploração da Implementação Actual de IPv6 no Kernel do Linux

A presente dissertação visa retratar a exploração do suporte do protocolo internet versão seis (IPv6) no kernel do Linux, conjuntamente com a análise detalhada do estado de implementação dos diferentes aspectos em que se baseia o protocolo. O estudo incide na experimentação do funcionamento em geral do stack, a identificação de inconsistências deste em relação RFC’s respectivos, bem como a simulação laboratorial de cenários que reproduzam casos de utilização de cada uma das facilidades analisadas. O objecto desta dissertação não é explicar o funcionamento do novo protocolo IPv6, mas antes, centrar-se essencialmente na exploração do IPv6 no kernel do Linux. Não é um documento para leigos em IPv6, no entanto, optou-se por desenvolver uma parte inicial onde é abordado o essencial do protocolo: a sua evolução até à aprovação e a sua especificação. Com base no estudo realizado, explora-se o suporte do IPv6 no kernel do Linux, fazendo uma análise detalhada do estudo de implementação dos diferentes aspectos em que se baseia o protocolo. Bem como a realização de testes de conformidade IPv6 em relação aos RFC’s.

Influence of reconstruction parameters during filtered backprojection and ordered-subset expectation maximization in the measurement of the left-ventricular volumes and function during gated SPECT.

A crucial method for investigating patients with coronary artery disease (CAD) is the calculation of the left ventricular ejection fraction (LVEF). It is, consequently, imperative to precisely estimate the value of LVEF--a process that can be done with myocardial perfusion scintigraphy. Therefore, the present study aimed to establish and compare the estimation performance of the quantitative parameters of the reconstruction methods filtered backprojection (FBP) and ordered-subset expectation maximization (OSEM). METHODS: A beating-heart phantom with known values of end-diastolic volume, end-systolic volume, and LVEF was used. Quantitative gated SPECT/quantitative perfusion SPECT software was used to obtain these quantitative parameters in a semiautomatic mode. The Butterworth filter was used in FBP, with the cutoff frequencies between 0.2 and 0.8 cycles per pixel combined with the orders of 5, 10, 15, and 20. Sixty-three reconstructions were performed using 2, 4, 6, 8, 10, 12, and 16 OSEM subsets, combined with several iterations: 2, 4, 6, 8, 10, 12, 16, 32, and 64. RESULTS: With FBP, the values of end-diastolic, end-systolic, and the stroke volumes rise as the cutoff frequency increases, whereas the value of LVEF diminishes. This same pattern is verified with the OSEM reconstruction. However, with OSEM there is a more precise estimation of the quantitative parameters, especially with the combinations 2 iterations × 10 subsets and 2 iterations × 12 subsets. CONCLUSION: The OSEM reconstruction presents better estimations of the quantitative parameters than does FBP. This study recommends the use of 2 iterations with 10 or 12 subsets for OSEM and a cutoff frequency of 0.5 cycles per pixel with the orders 5, 10, or 15 for FBP as the best estimations for the left ventricular volumes and ejection fraction quantification in myocardial perfusion scintigraphy.

Mutation identification by direct comparison of whole-genome sequencing data from mutant and wild-type individuals using k-mers.

Genes underlying mutant phenotypes can be isolated by combining marker discovery, genetic mapping and resequencing, but a more straightforward strategy for mapping mutations would be the direct comparison of mutant and wild-type genomes. Applying such an approach, however, is hampered by the need for reference sequences and by mutational loads that confound the unambiguous identification of causal mutations. Here we introduce NIKS (needle in the k-stack), a reference-free algorithm based on comparing k-mers in whole-genome sequencing data for precise discovery of homozygous mutations. We applied NIKS to eight mutants induced in nonreference rice cultivars and to two mutants of the nonmodel species Arabis alpina. In both species, comparing pooled F2 individuals selected for mutant phenotypes revealed small sets of mutations including the causal changes. Moreover, comparing M3 seedlings of two allelic mutants unambiguously identified the causal gene. Thus, for any species amenable to mutagenesis, NIKS enables forward genetics without requiring segregating populations, genetic maps and reference sequences.

Physical considerations on discrepancies in target volume delineation.

BACKGROUND AND PURPOSE: To compare the delineations and interpretations of target volumes by physicians in different radio-oncology centers. MATERIALS AND METHODS: Eleven Swiss radio-oncology centers delineated volumes according to ICRU 50 recommendations for one prostate and one head and neck case. In order to evaluate the consistency of the volume delineations, the following parameters were determined: 1) the target volumes (GTV, CTV and manually expanded PTV) and their extensions in the three main axes and 2) the correlation of the volume delineated by each pair of centers using the ratio of the intersection to the union (called proximity index). RESULTS: The delineated prostate volume was 105+/-55cm(3) for the CTV and 218+/-44cm(3) for the PTV. The delineated head and neck volume was 46+/-15cm(3) for the GTV, 327+/-154cm(3) for the CTV and 528+/-106cm(3) for the PTV. The mean proximity index for the prostate case was 0.50+/-0.13 for the CTV and 0.57+/-0.11 for the PTV. The proximity index for the head and neck case was 0.45+/-0.09 for the GTV, 0.42+/-0.13 for the CTV and 0.59+/-0.06 for the PTV. CONCLUSIONS: Large discrepancies between all the delineated target volumes were observed. There was an inverse relationship between the CTV volume and the margin between CTV and PTV, leading to less discrepancies in the PTV than is the CTV delineations. There was more spread in the sagittal and frontal planes due to CT pixel anisotropy, which suggests that radiation oncologists should delineate the target volumes not only in the transverse plane, but also in the sagittal and frontal planes to improve the delineation by allowing a consistency check.

Parametric imaging for characterizing focal liver lesions in contrast-enhanced ultrasound.

The differentiation between benign and malignant focal liver lesions plays an important role in diagnosis of liver disease and therapeutic planning of local or general disease. This differentiation, based on characterization, relies on the observation of the dynamic vascular patterns (DVP) of lesions with respect to adjacent parenchyma, and may be assessed during contrast-enhanced ultrasound imaging after a bolus injection. For instance, hemangiomas (i.e., benign lesions) exhibit hyper-enhanced signatures over time, whereas metastases (i.e., malignant lesions) frequently present hyperenhanced foci during the arterial phase and always become hypo-enhanced afterwards. The objective of this work was to develop a new parametric imaging technique, aimed at mapping the DVP signatures into a single image called a DVP parametric image, conceived as a diagnostic aid tool for characterizing lesion types. The methodology consisted in processing a time sequence of images (DICOM video data) using four consecutive steps: (1) pre-processing combining image motion correction and linearization to derive an echo-power signal, in each pixel, proportional to local contrast agent concentration over time; (2) signal modeling, by means of a curve-fitting optimization, to compute a difference signal in each pixel, as the subtraction of adjacent parenchyma kinetic from the echopower signal; (3) classification of difference signals; and (4) parametric image rendering to represent classified pixels as a support for diagnosis. DVP parametric imaging was the object of a clinical assessment on a total of 146 lesions, imaged using different medical ultrasound systems. The resulting sensitivity and specificity were 97% and 91%, respectively, which compare favorably with scores of 81 to 95% and 80 to 95% reported in medical literature for sensitivity and specificity, respectively.

Tertiary Himalayan structures and metamorphism in the Kulu Valley (Mandi-Khoksar transect of the Western Himalaya) - Shikar-Beh-nappe and crystalline nappe

The Crystalline Nappe of the High Himalayan Crystalline has been examined along the Kulu Valley and its vicinity (Mandi-Khoksar transect). This nappe was believed to have undergone deformation related only to its transport towards the SW essentially during the `'Main Central Thrust event''. New data has led to the conclusion that during the Himalayan orogeny, two distinctive phases, related to two opposite transport directions, characterize the evolution of this part of the chain, before the creation of the late NE-vergent backfolding. The first phase corresponds to an early NE-vergent folding and thrusting, creating the Tandi Syncline and the NE-oriented Shikar Beh Nappe stack, with a displacement amplitude of about 50 km. Two schistosities, together with a strong stretching lineation are developed at a deep tectonic level under amphibolite facies conditions (kyanite-staurolite-garnet-two mica schists). At a higher tectonic level and in the southern part of the section (Tandy Syncline and southern Kulu Valley between Kulu and Mandi) one or two schistosities are developed in the greenschist facies grade rocks (garnet-biotite and biotite schists). These structures and the associated Barrovian type metamorphism are all related to the NE-verging Shikar Beh Nappe. The creation of the NE-verging Shikar Beh Nappe may be explained by the reactivation of a SW dipping listric normal fault of the N Indian flexural passive margin, during the early stages of the Himalayan orogeny. In the second phase, the still hot metamorphic rocks of the Shikar Beh Nappe were folded and thrust towards the SW (mainly along the MBT and the MCT with a displacement in excess of 100 km) onto the cold, low-grade metamorphic rocks of the Larji-Kulu-Rampur Window or, near Mandi, on the non-metamorphic sandstones of the Ganges Molasse (Siwaliks). Sense of shear criteria and a strong NE-SW stretching-lineation indicate that the Crystalline Nappe has been overthrusted towards the SW. Thermometry on synkinematically crystallised garnet-biotite and garnet-hornblende pairs reveals the lower amphibolite facies temperature conditions related to the Crystalline Nappe formation. From the muscovite and biotite Rb-Sr cooling ages, the Shikar Beh Nappe emplacement occurred before 32 Ma and the southwestward thrusting of the Crystalline Nappe began before 21 Ma. Our model involving two opposite directions of thrusting goes against the conventional idea of only one main SW-oriented transport direction in the High Himalayan Crystalline Nappes.

Comparison of the polynomial model against explicit measurements of noise components for different mammography systems.

Given the adverse impact of image noise on the perception of important clinical details in digital mammography, routine quality control measurements should include an evaluation of noise. The European Guidelines, for example, employ a second-order polynomial fit of pixel variance as a function of detector air kerma (DAK) to decompose noise into quantum, electronic and fixed pattern (FP) components and assess the DAK range where quantum noise dominates. This work examines the robustness of the polynomial method against an explicit noise decomposition method. The two methods were applied to variance and noise power spectrum (NPS) data from six digital mammography units. Twenty homogeneously exposed images were acquired with PMMA blocks for target DAKs ranging from 6.25 to 1600 µGy. Both methods were explored for the effects of data weighting and squared fit coefficients during the curve fitting, the influence of the additional filter material (2 mm Al versus 40 mm PMMA) and noise de-trending. Finally, spatial stationarity of noise was assessed.Data weighting improved noise model fitting over large DAK ranges, especially at low detector exposures. The polynomial and explicit decompositions generally agreed for quantum and electronic noise but FP noise fraction was consistently underestimated by the polynomial method. Noise decomposition as a function of position in the image showed limited noise stationarity, especially for FP noise; thus the position of the region of interest (ROI) used for noise decomposition may influence fractional noise composition. The ROI area and position used in the Guidelines offer an acceptable estimation of noise components. While there are limitations to the polynomial model, when used with care and with appropriate data weighting, the method offers a simple and robust means of examining the detector noise components as a function of detector exposure.

Globally aligned photomosaic of the Lucky Strike hydrothermal vent field (Mid-Atlantic Ridge, 37°18.5’N): Release of georeferenced data, mosaic construction, and viewing software

We present a georeferenced photomosaic of the Lucky Strike hydrothermal vent field (Mid-Atlantic Ridge, 37°18’N). The photomosaic was generated from digital photographs acquired using the ARGO II seafloor imaging system during the 1996 LUSTRE cruise, which surveyed a ~1 km2 zone and provided a coverage of ~20% of the seafloor. The photomosaic has a pixel resolution of 15 mm and encloses the areas with known active hydrothermal venting. The final mosaic is generated after an optimization that includes the automatic detection of the same benthic features across different images (feature-matching), followed by a global alignment of images based on the vehicle navigation. We also provide software to construct mosaics from large sets of images for which georeferencing information exists (location, attitude, and altitude per image), to visualize them, and to extract data. Georeferencing information can be provided by the raw navigation data (collected during the survey) or result from the optimization obtained from imatge matching. Mosaics based solely on navigation can be readily generated by any user but the optimization and global alignment of the mosaic requires a case-by-case approach for which no universally software is available. The Lucky Strike photomosaics (optimized and navigated-only) are publicly available through the Marine Geoscience Data System (MGDS, http://www.marine-geo.org). The mosaic-generating and viewing software is available through the Computer Vision and Robotics Group Web page at the University of Girona (http://eia.udg.es/_rafa/mosaicviewer.html)

A noninvasive method for measuring the velocity of diffuse hydrothermal flow by tracking moving refractive index anomalies

Diffuse flow velocimetry (DFV) is introduced as a new, noninvasive, optical technique for measuring the velocity of diffuse hydrothermal flow. The technique uses images of a motionless, random medium (e.g.,rocks) obtained through the lens of a moving refraction index anomaly (e.g., a hot upwelling). The method works in two stages. First, the changes in apparent background deformation are calculated using particle image velocimetry (PIV). The deformation vectors are determined by a cross correlation of pixel intensities across consecutive images. Second, the 2-D velocity field is calculated by cross correlating the deformation vectors between consecutive PIV calculations. The accuracy of the method is tested with laboratory and numerical experiments of a laminar, axisymmetric plume in fluids with both constant and temperaturedependent viscosity. Results show that average RMS errors are ∼5%–7% and are most accurate in regions of pervasive apparent background deformation which is commonly encountered in regions of diffuse hydrothermal flow. The method is applied to a 25 s video sequence of diffuse flow from a small fracture captured during the Bathyluck’09 cruise to the Lucky Strike hydrothermal field (September 2009). The velocities of the ∼10°C–15°C effluent reach ∼5.5 cm/s, in strong agreement with previous measurements of diffuse flow. DFV is found to be most accurate for approximately 2‐D flows where background objects have a small spatial scale, such as sand or gravel

Enhanced radar precipitation estimates using a combined clutter and beam blockage correction technique

Weather radar observations are currently the most reliable method for remote sensing of precipitation. However, a number of factors affect the quality of radar observations and may limit seriously automated quantitative applications of radar precipitation estimates such as those required in Numerical Weather Prediction (NWP) data assimilation or in hydrological models. In this paper, a technique to correct two different problems typically present in radar data is presented and evaluated. The aspects dealt with are non-precipitating echoes - caused either by permanent ground clutter or by anomalous propagation of the radar beam (anaprop echoes) - and also topographical beam blockage. The correction technique is based in the computation of realistic beam propagation trajectories based upon recent radiosonde observations instead of assuming standard radio propagation conditions. The correction consists of three different steps: 1) calculation of a Dynamic Elevation Map which provides the minimum clutter-free antenna elevation for each pixel within the radar coverage; 2) correction for residual anaprop, checking the vertical reflectivity gradients within the radar volume; and 3) topographical beam blockage estimation and correction using a geometric optics approach. The technique is evaluated with four case studies in the region of the Po Valley (N Italy) using a C-band Doppler radar and a network of raingauges providing hourly precipitation measurements. The case studies cover different seasons, different radio propagation conditions and also stratiform and convective precipitation type events. After applying the proposed correction, a comparison of the radar precipitation estimates with raingauges indicates a general reduction in both the root mean squared error and the fractional error variance indicating the efficiency and robustness of the procedure. Moreover, the technique presented is not computationally expensive so it seems well suited to be implemented in an operational environment.

A wavelet-based method for the determination of the relative resolution between remotely sensed images

Spatial resolution is a key parameter of all remote sensing satellites and platforms. The nominal spatial resolution of satellites is a well-known characteristic because it is directly related to the area in ground that represents a pixel in the detector. Nevertheless, in practice, the actual resolution of a specific image obtained from a satellite is difficult to know precisely because it depends on many other factors such as atmospheric conditions. However, if one has two or more images of the same region, it is possible to compare their relative resolutions. In this paper, a wavelet-decomposition-based method for the determination of the relative resolution between two remotely sensed images of the same area is proposed. The method can be applied to panchromatic, multispectral, and mixed (one panchromatic and one multispectral) images. As an example, the method was applied to compute the relative resolution between SPOT-3, Landsat-5, and Landsat-7 panchromatic and multispectral images taken under similar as well as under very different conditions. On the other hand, if the true absolute resolution of one of the images of the pair is known, the resolution of the other can be computed. Thus, in the last part of this paper, a spatial calibrator that is designed and constructed to help compute the absolute resolution of a single remotely sensed image is described, and an example of its use is presented.

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.

Correlação linear e espacial entre a resistência do solo ao penetrômetro e a produtividade do feijoeiro irrigado

A resistência do solo ao penetrômetro exerce grande influência sobre o crescimento e desenvolvimento vegetal, uma vez que o crescimento das raízes, assim como o rendimento das culturas, varia de forma inversamente proporcional ao seu valor. Dessa forma, a análise da variabilidade espacial da resistência do solo ao penetrômetro e da produtividade, por meio da geoestatística, pode indicar alternativas de manejo para reduzir os efeitos da variabilidade do solo sobre a produtividade e também melhorar a estimativa de respostas das culturas sob determinadas práticas de manejo. Diante do exposto, o objetivo deste trabalho foi relacionar e caracterizar a variabilidade espacial da resistência do solo ao penetrômetro (RP) e a produtividade do feijoeiro irrigado em sistema de semeadura direta, em duas safras consecutivas. O experimento foi realizado em Latossolo Vermelho distroférrico típico, no campo experimental da Faculdade de Engenharia Agrícola da Unicamp, no município de Campinas-SP, cujas coordenadas geográficas são: 22 ° 48 ' 57 " de latitude sul, 47 ° 03 ' 33 " de longitude oeste e altitude média de 640 m. As avaliações foram realizadas em uma malha regular de amostragem de 3 x 3 m, totalizando 60 pontos amostrais por parcela. A análise da dependência espacial foi avaliada pela geoestatística, e os parâmetros dos semivariogramas utilizados para construir mapas de isolinhas, por meio do interpolador de krigagem do programa Surfer 8.0. A regressão linear simples entre mapas (pixel-a-pixel) mostrou correlação negativa entre os valores de RP e a produtividade; no entanto, a produtividade do feijoeiro irrigado apresentou baixa correlação com a resistência do solo ao penetrômetro em sistema semeadura direta nas duas safras.

Optical-geometrical effects on the photoluminescence spectra of Si nanocrystals embedded in SiO2

We demonstrate that thickness, optical constants, and details of the multilayer stack, together with the detection setting, strongly influence the photoluminescence spectra of Si nanocrystals embedded in SiO2. Due to multiple reflections of the visible light against the opaque silicon substrate, an interference pattern is built inside the oxide layer, which is responsible for the modifications in the measured spectra. This interference effect is complicated by the depth dependence of (i) the intensity of the excitation laser and (ii) the concentration of the emitting nanocrystals. These variations can give rise to apparent features in the recorded spectra, such as peak shifts, satellite shoulders, and even splittings, which can be mistaken as intrinsic material features. Thus, they can give rise to an erroneous attribution of optical bands or estimate of the average particle size, while they are only optical-geometrical artifacts. We have analyzed these effects as a function of material composition (Si excess fraction) and thickness, and also evaluated how the geometry of the detection setup affects the measurements. To correct the experimental photoluminescence spectra and extract the true spectral shape of the emission from Si nanocrystals, we have developed an algorithm based on a modulation function, which depends on both the multilayer sequence and the experimental configuration. This procedure can be easily extended to other heterogeneous systems.