Steady-state equilibrium phase inversion recovery ON-resonant water suppression (IRON) MR angiography in conjunction with superparamagnetic nanoparticles. A robust technique for imaging within a wide range of contrast agent dosages.

PURPOSE: To investigate the ability of inversion recovery ON-resonant water suppression (IRON) in conjunction with P904 (superparamagnetic nanoparticles which consisting of a maghemite core coated with a low-molecular-weight amino-alcohol derivative of glucose) to perform steady-state equilibrium phase MR angiography (MRA) over a wide dose range. MATERIALS AND METHODS: Experiments were approved by the institutional animal care committee. Rabbits (n = 12) were imaged at baseline and serially after the administration of 10 incremental dosages of 0.57-5.7 mgFe/Kg P904. Conventional T1-weighted and IRON MRA were obtained on a clinical 1.5 Tesla (T) scanner to image the thoracic and abdominal aorta, and peripheral vessels. Contrast-to-noise ratios (CNR) and vessel sharpness were quantified. RESULTS: Using IRON MRA, CNR and vessel sharpness progressively increased with incremental dosages of the contrast agent P904, exhibiting constantly higher contrast values than T1 -weighted MRA over a very wide range of contrast agent doses (CNR of 18.8 ± 5.6 for IRON versus 11.1 ± 2.8 for T1 -weighted MRA at 1.71 mgFe/kg, P = 0.02 and 19.8 ± 5.9 for IRON versus -0.8 ± 1.4 for T1-weighted MRA at 3.99 mgFe/kg, P = 0.0002). Similar results were obtained for vessel sharpness in peripheral vessels, (Vessel sharpness of 46.76 ± 6.48% for IRON versus 33.20 ± 3.53% for T1-weighted MRA at 1.71 mgFe/Kg, P = 0.002, and of 48.66 ± 5.50% for IRON versus 19.00 ± 7.41% for T1-weighted MRA at 3.99 mgFe/Kg, P = 0.003). CONCLUSION: Our study suggests that quantitative CNR and vessel sharpness after the injection of P904 are consistently higher for IRON MRA when compared with conventional T1-weighted MRA. These findings apply for a wide range of contrast agent dosages.

Hyperpolarizing gases via dynamic nuclear polarization and sublimation.

A high throughput method was designed to produce hyperpolarized gases by combining low-temperature dynamic nuclear polarization with a sublimation procedure. It is illustrated by applications to 129Xe nuclear magnetic resonance in xenon gas, leading to a signal enhancement of 3 to 4 orders of magnitude compared to the room-temperature thermal equilibrium signal at 7.05 T.

Phase-sensitive dual-inversion recovery for accelerated carotid vessel wall imaging.

OBJECTIVES: Dual-inversion recovery (DIR) is widely used for magnetic resonance vessel wall imaging. However, optimal contrast may be difficult to obtain and is subject to RR variability. Furthermore, DIR imaging is time-inefficient and multislice acquisitions may lead to prolonged scanning times. Therefore, an extension of phase-sensitive (PS) DIR is proposed for carotid vessel wall imaging. METHODS: The statistical distribution of the phase signal after DIR is probed to segment carotid lumens and suppress their residual blood signal. The proposed PS-DIR technique was characterized over a broad range of inversion times. Multislice imaging was then implemented by interleaving the acquisition of 3 slices after DIR. Quantitative evaluation was then performed in healthy adult subjects and compared with conventional DIR imaging. RESULTS: Single-slice PS-DIR provided effective blood-signal suppression over a wide range of inversion times, enhancing wall-lumen contrast and vessel wall conspicuity for carotid arteries. Multislice PS-DIR imaging with effective blood-signal suppression is enabled. CONCLUSIONS: A variant of the PS-DIR method has successfully been implemented and tested for carotid vessel wall imaging. This technique removes timing constraints related to inversion recovery, enhances wall-lumen contrast, and enables a 3-fold increase in volumetric coverage at no extra cost in scanning time.

Implications of the serpentine phase transition on the behaviour of beryllium and lithium-boron of subducted ultramafic rocks

The Totalp-Platta-Malenco ophiolites in the Eastern Central Alps offer a unique opportunity to study the behaviour of Li, Be and B in ultramafic rocks in response to serpentinization and to progressive Alpine metamorphism. These units represent the remnants of a former ocean-continent transition that was intensely serpentinized during exposure on the Jurassic seafloor of the Ligurian Tethys. From north to the south, three isograd reactions (lizardite double right arrow antigorite + brucite; lizardite + talc double right arrow antigorite; lizardite + tremolite double right arrow antigorite + diopside) have been used to quantify the evolution of the light element content of metamorphic minerals. We determined the Li, Be and B concentrations in major silicate minerals from the ultramafic bodies of Totalp, Platta and Malenco by secondary ion mass spectrometry. Mantle minerals have Be concentrations (e.g. <0.001-0.009 mu g/g in olivine) similar to the metamorphic minerals that replace them (e.g. <0.001-0.016 mu g/g in serpentine). The mantle signature of Be is thus neither erased during seafloor alteration nor by progressive metamorphism from prehnite-pumpellyite to epidote-amphibolite facies. In contrast, the Li and B inventories of metamorphic minerals are related to the lizardite-to-antigorite transition. Both elements display higher concentrations in the low-temperature serpentine polymorph lizardite (max. 156 mu/g Li, max. 318 mu g/g B) than in antigorite (max. 0.11 mu g/g Li, max. 12 mu g/g B). Calculated average B/Li ratios for lizardite (similar to 1395) and antigorite (similar to 115) indicate that Li fractionates from B during the lizardite-to-antigorite transition during prograde metamorphism in ultramafic rocks. In subduction zones, this signature is likely to be recorded in the B-rich nature of forearc fluids. Relative to oceanic mantle the Be content of mantle clinopyroxene is much higher, but similar to Be values from mantle xenoliths and subduction-related peridotite massifs. These data support previous hypothesis that the mantle rocks from the Eastern Central Alps have a subcontinental origin. We conclude that Be behaves conservatively during subduction metamorphism of ultramafic rocks, at least at low-temperature, and thus retains the fingerprint of ancient subduction-related igneous events in mantle peridotites. (C) 2010 Elsevier Ltd. All rights reserved.

Phase-sensitive black-blood coronary vessel wall imaging.

Black-blood MR coronary vessel wall imaging may become a powerful tool for the quantitative and noninvasive assessment of atherosclerosis and positive arterial remodeling. Although dual-inversion recovery is currently the gold standard, optimal lumen-to-vessel wall contrast is sometimes difficult to obtain, and the time window available for imaging is limited due to competing requirements between blood signal nulling time and period of minimal myocardial motion. Further, atherosclerosis is a spatially heterogeneous disease, and imaging at multiple anatomic levels of the coronary circulation is mandatory. However, this requirement of enhanced volumetric coverage comes at the expense of scanning time. Phase-sensitive inversion recovery has shown to be very valuable for enhancing tissue-tissue contrast and for making inversion recovery imaging less sensitive to tissue signal nulling time. This work enables multislice black-blood coronary vessel wall imaging in a single breath hold by extending phase-sensitive inversion recovery to phase-sensitive dual-inversion recovery, by combining it with spiral imaging and yet relaxing constraints related to blood signal nulling time and period of minimal myocardial motion.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

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.

Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole georadar data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a time-domain waveform inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little-to-no trade-off between the wavelet estimation and the tomographic imaging procedures.

Juxtaposition of melt impregnation and high-temperature shear zones in the upper mantle; field and petrological constraints from the Lanzo peridotite (Northern Italy)

Results of a field and microstructural study between the northern and the central bodies of the Lanzo plagioclase peridotite massif (NW Italy) indicate that the spatial distribution of deformation is asymmetric across kilometre-scale mantle shear zones. The southwestern part of the shear zone (footwall) shows a gradually increasing degree of deformation from porphyroclastic peridotites to mylonite, whereas the northeastern part (hanging wall) quickly grades into weakly deformed peridotites. Discordant gabbroic and basaltic dykes are asymmetrically distributed and far more abundant in the footwall of the shear zone. The porphyroclastic peridotite displays porphyroclastic zones and domains of igneous crystallization whereas mylonites are characterized by elongated porphyroclasts, embedded between fine-grained, polycrystalline bands of olivine, plagioclase, clinopyroxene, orthopyroxene, spinel, rare titanian pargasite, and domains of recrystallized olivine. Two types of melt impregnation textures have been found: (1) clinopyroxene porphyroclasts incongruently reacted with migrating melt to form orthopyroxene plagioclase; (2) olivine porphyroclasts are partially replaced by interstitial orthopyroxene. The meltrock reaction textures tend to disappear in the mylonites, indicating that deformation in the mylonite continued under subsolidus conditions. The pyroxene chemistry is correlated with grain size. High-Al pyroxene cores indicate high temperatures (11001030C), whereas low-Al neoblasts display lower final equilibration temperatures (860C). The spinel Cr-number [molar Cr/(Cr Al)] and TiO2 concentrations show extreme variability covering almost the entire range known from abyssal peridotites. The spinel compositions of porphyroclastic peridotites from the central body are more variable than spinel from mylonite, mylonite with ultra-mylonite bands, and porphyroclastic rocks of the northern body. The spinel compositions probably indicate disequilibrium and would favour rapid cooling, and a faster exhumation of the central peridotite body, relative to the northern one. Our results indicate that melt migration and high-temperature deformation are juxtaposed both in time and space. Meltrock reaction may have caused grain-size reduction, which in turn led to localization of deformation. It is likely that melt-lubricated, actively deforming peridotites acted as melt focusing zones, with permeabilities higher than the surrounding, less deformed peridotites. Later, under subsolidus conditions, pinning in polycrystalline bands in the mylonites inhibited substantial grain growth and led to permanent weak zones in the upper mantle peridotite, with a permeability that is lower than in the weakly deformed peridotites. Such an inversion in permeability might explain why actively deforming, fine-grained peridotite mylonite acted as a permeability barrier and why ascending mafic melts might terminate and crystallize as gabbros along actively deforming shear zones. Melt-lubricated mantle shear zones provide a mechanism for explaining the discontinuous distribution of gabbros in oceancontinent transition zones, oceanic core complexes and ultraslow-spreading ridges.

Stable isotope geochemistry and phase-equilibria of coesite-bearing whiteschists, Dora Maira massif, Western Alps

Peak metamorphic temperatures for the coesite-pyrope-bearing whiteschists from the Dora Maira Massif, western Alps were determined with oxygen isotope thermometry. The deltaO-18(SMOW) values of the quartz (after coesite) (delta O-18 = 8.1 to 8.6 parts per thousand, n = 6), phengite (6.2 to 6.4 parts per thousand, n = 3), kyanite (6.1 parts per thousand, n = 2), garnet (5.5 to 5.8 parts per thousand, n = 9), ellenbergerite (6.3 parts per thousand, n = 1) and rutile (3.3. to 3.6 parts per thousand, n = 3) reflect isotopic equilibrium. Temperature estimates based on quartz-garnet-rutile fractionation are 700-750-degrees-C. Minimum pressures are 31-32 kb based on the pressure-sensitive reaction pyrope + coesite = kyanite + enstatite. In order to stabilize pyrope and coesite by the temperature-sensitive dehydration reaction talc + kyanite = pyrope + coesite + H2O, the a(H2O) must be reduced to 0.4-0.75 at 700 750-degrees-C. The reduced a(H2O) cannot be due to dilution by CO2, as pyrope is not stable at X (CO2) > 0.02 (T = 750-degrees-C; P = 30 kb). In the absence of a more exotic fluid diluent (e.g. CH4 or N2), a melt phase is required. Granite solidus temperatures are approximately 680-degrees-C/30 kb at a(H2O) = 1.0 and are calculated to be approximately 70-degrees-C higher at a(H2O) = 0.7, consistent with this hypothesis. Kyanite-jadeite-quartz bands may represent a relict melt phase. Peak P-T-f(H2O) estimates for the whiteschist are 34 +/- 2 kb, 700-750-degrees-C and 0.4-0.75. The oxygen isotope fractionation between quartz (deltaO-18 = 11.6%.) and garnet (deltaO-18 = 8.7 parts per thousand) in the surrounding orthognesiss is identical to that in the coesite-bearing unit, suggesting that the two units shared a common, final metamorphic history. Hydrogen isotope measurements were made on primary talc and phengite (deltaD(smow) = -27 to -32 parts per thousand), on secondary talc and chlorite after pyrope (deltaD = - 39 to - 44 parts per thousand) and on the surrounding biotite (deltaD = -64 parts per thousand) and phengite (deltaD = -44 parts per thousand) gneiss. All phases appear to be in near-equilibrium. The very high deltaD values for the primary hydrous phases is consistent with an initial oceanic-derived/connate fluid source. The fluid source for the retrograde talc + chlorite after pyrope may be fluids evolved locally during retrograde melt crystallization. The similar deltaD, but dissimilar deltaO-18 values of the coesite-bearing whiteschists and hosting orthogneiss suggest that the two were in hydrogen isotope equilibrium, but not oxygen isotope equilibrium. The unusual hydrogen and oxygen isotope compositions of the coesite-bearing unit can be explained as the result of metasomatism from slab-derived fluids at depth.

Source-wavelet estimation during full-waveform inversion of ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole ground-penetrating radar (GPR) data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a recently published time-domain inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity of both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little to no trade-off between the wavelet estimation and the tomographic imaging procedures.

Quantification of magmatic and hydrothermal processes in a peralkaline syenite-alkali granite complex based on textures, phase equilibria, and stable and radiogenic isotopes

The Puklen complex of the Mid-Proterozoic Gardar Province, South Greenland, consists of various silica-saturated to quartz-bearing syenites, which are intruded by a peralkaline granite. The primary mafic minerals in the syenites are augite +/- olivine + Fe-Ti oxide + amphibole. Ternary feldspar thermometry and phase equilibria among mafic silicates yield T = 950-750degreesC, a(SiO2) = 0.7-1 and an f(O2) of 1-3 log units below the fayalite-magnetite-quartz (FMQ) buffer at 1 kbar. In the granites, the primary mafic minerals are ilmenite and Li-bearing arfvedsonite, which crystallized at temperatures below 750degreesC and at f(O2) values around the FMQ buffer. In both rock types, a secondary post-magmatic assemblage overprints the primary magmatic phases. In syenites, primary Ca-bearing minerals are replaced by Na-rich minerals such as aegirine-augite and albite, resulting in the release of Ca. Accordingly, secondary minerals include ferro-actinolite, (calcite-siderite)(ss), titanite and andradite in equilibrium with the Na-rich minerals. Phase equilibria indicate that formation of these minerals took place over a long temperature interval from near-magmatic temperatures down to similar to300degreesC. In the course of this cooling, oxygen fugacity rose in most samples. For example, late-stage aegirine in granites formed at the expense of arfvedsonite at temperatures below 300degreesC and at an oxygen fugacity above the haematite-magnetite (HM) buffer. The calculated delta(18)O(melt) value for the syenites (+5.9 to +6.3parts per thousand) implies a mantle origin, whereas the inferred delta(18)O(melt) value of <+5.1parts per thousand for the granitic melts is significantly lower. Thus, the granites require an additional low-delta(18)O contaminant, which was not involved in the genesis of the syenites. Rb/Sr data for minerals of both rock types indicate open-system behaviour for Rb and Sr during post-magmatic metasomatism. Neodymium isotope compositions (epsilonNd(1170 Ma) = -3.8 to -6.4) of primary minerals in syenites are highly variable, and suggest that assimilation of crustal rocks occurred to variable extents. Homogeneous epsilon(Nd) values of -5.9 and -6.0 for magmatic amphibole in the granites lie within the range of the syenites. Because of the very similar neodymium isotopic compositions of magmatic and late- to post-magmatic minerals from the same syenite samples a principally closed-system behaviour during cooling is implied. In contrast, for the granites an externally derived fluid phase is required to explain the extremely low epsilon(Nd) values of about -10 and low delta(18)O between +2.0 and +0.5parts per thousand for late-stage aegirine, indicating an open system in the late-stage history. In this study we show that the combination of phase equilibria constraints with stable and radiogenic isotope data on mineral separates can provide much better constraints on magma evolution during emplacement and crystallization than conventional whole-rock studies.

Epidote in hydrous mafic magmas : near solidus phase relations in the lower crust

Crystallisation of hydrous mafic magmas at high pressure is a subject of numerous petrologic and experimental studies since the last century and is mainly related to the process of continental crust formation and the possible link between mantle derived melts and low pressure granitoids. Albeit the sequence of crystallization is well constrained by experimental studies, the origin of exposed lower crustal rocks exposed on the earth surface is controversial. Ones line of argument is favouring high pressure crystallization of dry or wet mafic magmas, whereas others invoke partial melting of pre-existing crust. Therefore studies involving field, textural and chemical observations of exposed lower crust such as in Kohistan (Pakistan) or Talkeetna (Alaska) are crucial to understand the continental crust formation processes via arc magmatism. Epidote-bearing gabbros are very sparse and always associated with the deep part of continental crust (>30 km) as in the Kohistan Arc Complex (Pakistan) or in the Chelan Complex (western U.S.). Magmatic epidote is restricted to a small temperature interval above the water-saturated solidus of MORB and represent the last crystallizing liquids in lower crustal regions. However, epidote and melt stability at lower crustal pressures are not clearly established.¦The Chelan complex (western U.S.) at the base of the Cascadian Arc is composed mainly by peraluminous tonalité associated with gabbroic and ultramafic rocks and was traditionally interpreted as a migmatitic terrain. However field, chemical and mineralogical observations rather suggest a magmatic origin and point to a protracted crystallization at intermediate to high pressure ~ 1.0 GPa dominated by amphibole fractionation and followed by isobaric cooling down to 650°C. Crystal fractionation modelling using whole rock composition and field constraints is able to generate peraluminous tonalité. The modelled crystallisation sequence and the volume proportions are in agreement with experimental studies performed at these pressures. The Chelan complex was thus not formed during a partial melting event, but represent the sequence of crystallisation occurring at the base of the crust. Massive fractionation of hornblende is able to generate peraluminous tonalité without significant assimilation of crustal rock.¦Similarly to the Chelan complex, the base of the Kohistan arc is composed of cumulates derived by high pressure crystallization of hydrous magma. In garnet gabbros, epidote occurs as magmatic phase, crystallising from hydrous interstitial melt trapped between grain boundaries at lower crustal pressures (Ρ ~ 1.2 GPa) for temperature of (650-700 °C). Trace and REE signature in epidote indicate that epidote was formed through peritectic reaction involving garnet, clinopyroxene and plagioclase. At the beginning of the crystallisation epidote signature is dominated by REE content in the melt, whereas at the end the signature is dominated by reacting phases. Melt in equilibrium with epidote inferred from the partition coefficients available is similar to intrusive tonalité up the section indicating that hydrous melt was extracted from the garnet gabbros. In some gabbros epidote shows single homogeneous compositions, while in others coexisting epidote have different compositions indicating the presence of solvi along the Al-Fe3+ join. The overgowths are only observed in presence of paragonite in the assemblage, suggesting high water content. At high water content, the hydrous solidus is shift to lower temperature and probably intersects the solvi observed along the Al-Fe3+ join. Therefore, several compositions of epidote is stable at high water content.¦-¦La composition chimique de la croûte continentale est considérée comme similaire à celle du magmatisme calco-alcalin de marge continentale active (enrichissement en éléments mobiles dans les fluides, anomalies négatives en Nb, Ta et éléments à haut potentiel électronique, etc...). Cependant la nature andésitique de la croûte continentale (Si02 > 60 wt%), résultant des nombreuses intrusions de granitoïdes dans la croûte supérieure, est sujette à polémique et le lien entre les magmas dérivés du manteau et les roches évoluées de faible profondeur n'est pas clairement établi (fusion partielle de croûte basaltique, cristallisation fractionnée à haute pression, etc...).¦Les affleurements de croûte profonde sont rares mais précieux, car ils permettent d'observer les phénomènes se passant à grande profondeur. Le complexe de Chelan (Washington Cascades) en est un exemple. Formé à environ 30 km de profondeur, il est composé de roches gabbroïques et ultramafiques, ainsi que de tonalités, qui furent souvent interprétés comme le produit de la fusion partielle de la croûte. Cependant, les relations de terrain, la chimie des éléments majeurs et des éléments traces sont cohérentes avec l'évolution d'un complexe magmatique mafique dans la croûte profonde ou moyenne ( 1.0 GPa), dominée par le fractionnement de l'amphibole. Après son emplacement, le complexe a subi un refroidissement isobare jusqu'à des températures de l'ordre de 650 °C, déduit de la composition chimique des minéraux. Un bilan de masse contraint pax les observations de terrain permet de calculer la séquence et les volumes de fractionnement. Les faciès évolués légèrement hyperalumineux observés sur le terrain peuvent être générés par la cristallisation de 3 % de websterite à olivine, 12 % d'hornblendite à pyroxène 33 % d'hornblendite, 19 % de gabbros, 15 % de diorite et 2 % de tonalité. Nous montrons ainsi qu'une série de fractionnement contrôlée par l'amphibole permet de générer des tonalités sans assimilation de matériel crustal et l'exemple de Chelan illustre la viabilité de ce processus dans la formation de croûte continentale.¦Les réactions proches du solidus saturé en H20 dans les systèmes basaltiques à des pressions élevées restent énigmatiques. Diverses expériences tendent à montrer que l'épidote est stable dans ces conditions, mais rarement observée (décrite ?) comme phase primaire dans les systèmes naturels. Les épidotes trouvées dans les gabbros de Jijal (nord-Pakistan) montrent des textures de type .magmatique telles qu'observées dans les roches évoluées. Le contenu en terres rares de ces épidotes est très variable allant de signatures enrichies en terres rares légères impliquant la présence de liquide interstitiel à des signatures complètement déprimées en ces mêmes éléments, évoquant une cristallisation en coexistence avec du grenat. Ces diverses signatures reflètent un chemin de cristallisation en présence de liquide interstitiel et enregistrent des réactions péritectiques impliquant grenat, clinopyroxene et plagioclase à des pressions de ~ 1.2 GPa pour des températures de 650-700 °C. Cependant dans quelques échantillons deux ou trois compositions d'épidotes coexistent démontrant la présence de lacunes d'immiscibilité le long de la solution solide épidote-clinozoïsite. La forte teneur en H20 du liquide magmatique est certainement à l'origine de la coexistence de deux compositions distinctes.

Decreased thermal conductance during the luteal phase of the menstrual cycle in women.

To study the influence of the menstrual cycle on whole body thermal balance and on thermoregulatory mechanisms, metabolic heat production (M) was measured by indirect calorimetry and total heat losses (H) were measured by direct calorimetry in nine women during the follicular (F) and the luteal (L) phases of the menstrual cycle. The subjects were studied while exposed for 90 min to neutral environmental conditions (ambient temperature 28 degrees C, relative humidity 40%) in a direct calorimeter. The values of M and H were not modified by the phase of the menstrual cycle. Furthermore, in both phases the subjects were in thermal equilibrium because M was similar to H (69.7 +/- 1.8 and 72.1 +/- 1.8 W in F and 70.4 +/- 1.9 and 71.4 +/- 1.7 W in L phases, respectively). Tympanic temperature (Tty) was 0.24 +/- 0.07 degrees C higher in the L than in the F phase (P less than 0.05), whereas mean skin temperature (Tsk) was unchanged. Calculated skin thermal conductance (Ksk) was lower in the L (17.9 +/- 0.6 W.m-2.degrees C-1) than in the F phase (20.1 +/- 1.1 W.m-2.degrees C-1; P less than 0.05). Calculated skin blood flow (Fsk) was also lower in the L (0.101 +/- 0.008 l.min-1.m-2) than in the F phase (0.131 +/- 0.015 l.min-1.m-2; P less than 0.05). Differences in Tty, Ksk, and Fsk were not correlated with changes in plasma progesterone concentration. It is concluded that, during the L phase, a decreased thermal conductance in women exposed to a neutral environment allows the maintenance of a higher internal temperature.