A Small Baseline DIFSAR Approach for Investigating Deformations on Full Resolution SAR Interferograms

This paper presents a differential synthetic apertureradar (SAR) interferometry (DIFSAR) approach for investigatingdeformation phenomena on full-resolution DIFSAR interferograms.In particular, our algorithm extends the capabilityof the small-baseline subset (SBAS) technique that relies onsmall-baseline DIFSAR interferograms only and is mainly focusedon investigating large-scale deformations with spatial resolutionsof about 100 100 m. The proposed technique is implemented byusing two different sets of data generated at low (multilook data)and full (single-look data) spatial resolution, respectively. Theformer is used to identify and estimate, via the conventional SBAStechnique, large spatial scale deformation patterns, topographicerrors in the available digital elevation model, and possibleatmospheric phase artifacts; the latter allows us to detect, onthe full-resolution residual phase components, structures highlycoherent over time (buildings, rocks, lava, structures, etc.), as wellas their height and displacements. In particular, the estimation ofthe temporal evolution of these local deformations is easily implementedby applying the singular value decomposition technique.The proposed algorithm has been tested with data acquired by theEuropean Remote Sensing satellites relative to the Campania area(Italy) and validated by using geodetic measurements.

Efecto de los usos del suelo en la emisión de dióxido de carbono del suelo a la atmósfera en un agroecosistema semiárido del Valle del Ebro

Los cambios en los usos del suelo han contribuido de manera importante al incremento de gases de efecto invernadero en la atmósfera, especialmente de dióxido de carbono, aumentando sus emisiones desde 1970 en un 80%. Estos cambios causan la alteración de los suelos provocando un impacto sobre el ciclo del carbono, aumentando las tasas de descomposición de la fracción orgánica creando así un flujo de CO2 a la atmosfera. Entre las recomendaciones del Panel Intergubernamental de expertos sobre el Cambio Climático (IPCC, en inglés), y contemplado en el Protocolo de Kyoto, se encuentra el proceso de secuestro de carbono en suelos, que implica la eliminación del CO2 atmosférico por parte de las plantas y su almacenamiento como materia orgánica del suelo. Para poder favorecer dicho proceso, en un determinado tipo de ecosistema, es fundamental conocer cuáles son los factores que gobiernan la respiración del suelo y el impacto que tienen los diferentes usos en la emisión de CO2. En el presente trabajo se han estudiado 4 usos del suelo representativos del secano aragonés: un monocultivo de cebada en siembra directa (NT), un suelo abandonado labrado (AC), un suelo abandonado no alterado (AU) y un suelo forestal (FR) con el objetivo de conocer sus tasas de respiración, la influencia de diferentes parámetros edáficos en ellas, y proponer cambios en el uso del suelo que ayuden a mitigar estas emisiones. Además, se ha dedicado un apartado para conocer cómo influyen diferentes técnicas de fertilización nitrogenada (mineral y orgánica) en la respiración de un monocultivo de cebada en siembra directa. En cuanto a los usos, los resultados obtenidos tanto in situ como en laboratorio muestran una mayor respiración en AC, siendo los valores más bajos los de NT y FR. Una de las principales conclusiones es que la supresión del laboreo y del periodo de barbecho largo, así como la conversión de tierras abandonadas y marginales a cultivos y zonas forestales se presentan en este tipo de ecosistemas como prácticas de secuestro de carbono. En el estudio de aplicación de fertilizantes, no se observó ningún cambio en la respiración del suelo después de la aplicación de nitrógeno mineral. En cambio, el suelo fertilizado con purín sí que mostró picos de emisión durante las siguientes horas a la incorporación de éste, debido fundamentalmente a su alto contenido de carbono lábil.

Declining labor force attachment and downward trends in unemployment and participation

The US labor market witnessed two apparently unrelated secular movements in thelast 30 years: a decline in unemployment between the early 1980s and the early 2000s,and a decline in participation since the early 2000s. Using CPS micro data and a stock-flow accounting framework, we show that a substantial, and hitherto unnoticed, factorbehind both trends is a decline in the share of nonparticipants who are at the margin ofparticipation. A lower share of marginal nonparticipants implies a lower unemploymentrate, because marginal nonparticipants enter the labor force mostly through unemployment,while other nonparticipants enter the labor force mostly through employment.

Transformations in occluded light fraction organic matter in a clayey oxisol: evidence from 13C-CPMAS-NMR and delta13C Signature

We hypothesised that, during occlusion inside granular aggregates of oxide-rich soils, the light fraction organic matter would undergo a strong process of decomposition, either due to the slow process of aggregate formation and stabilisation or due to digestion in the macro- and meso-fauna guts. This process would favour the accumulation of recalcitrant materials inside aggregates. The aim of this study was to compare the dynamics and the chemical composition of free and occluded light fraction organic matter in a natural cerrado vegetation (woodland savannah) and a nearby pasture (Brachiaria spp.) to elucidate the transformations during occlusion of light fraction in aggregates of a clayey Oxisol. Nuclear Magnetic Resonance of the 13C, with Cross Polarisation and Magic Angle Spinning (13C-CPMAS-NMR), and 13C/12C isotopic ratio were combined to study organic matter composition and changes in carbon dynamics, respectively. The occluded light fraction had a slower turnover than the free light fraction and the heavy fraction. Organic matter in the occluded fraction also showed a higher degree of decomposition. The results confirm that processes of soil organic matter occlusion in the typical "very fine strong granular" structure of the studied oxide-rich soil led to an intense transformation, selectively preserving stable organic matter. The small amount of organic material stored as occluded light faction, as well as its stability, suggests that this is not an important or manageable sink for sequestration of atmospheric CO2.

Bayesian image reconstruction with space-variant noise suppression

In this paper we present a Bayesian image reconstruction algorithm with entropy prior (FMAPE) that uses a space-variant hyperparameter. The spatial variation of the hyperparameter allows different degrees of resolution in areas of different statistical characteristics, thus avoiding the large residuals resulting from algorithms that use a constant hyperparameter. In the first implementation of the algorithm, we begin by segmenting a Maximum Likelihood Estimator (MLE) reconstruction. The segmentation method is based on using a wavelet decomposition and a self-organizing neural network. The result is a predetermined number of extended regions plus a small region for each star or bright object. To assign a different value of the hyperparameter to each extended region and star, we use either feasibility tests or cross-validation methods. Once the set of hyperparameters is obtained, we carried out the final Bayesian reconstruction, leading to a reconstruction with decreased bias and excellent visual characteristics. The method has been applied to data from the non-refurbished Hubble Space Telescope. The method can be also applied to ground-based images.

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.

A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

Nitrogen application timing and soil inorganic nitrogen dynamics under no-till oat/maize sequential cropping

The timing of N application to maize is a key factor to be considered in no-till oat/maize sequential cropping. This study aimed to evaluate the influence of pre-planting, planting and sidedress N application on oat residue decomposition, on soil N immobilisation and remineralisation and on N uptake by maize plants in no-till oat/maize sequential cropping. Undisturbed soil cores of 10 and 20 cm diameter were collected from the 0-15 cm layer of a no-till Red Latossol, when the oat cover crop was in the milk-grain stage. Two greenhouse experiments were conducted simultaneously. Experiment A, established in the 10 cm diameter cores and without plant cultivation, was used to asses N dynamics in soil and oat residues. Experiment B, established in the 20 cm diameter cores and with maize cultivation, was used to assess plant growth and N uptake. An amount of 6.0 Mg ha-1 dry matter of oat residues was spread on the surface of the cores. A rate of 90 kg N ha-1 applied as ammonium sulphate in both experiments was split in pre-planting, planting and sidedress applications as follows: (a) 00-00-00 (control), (b) 90-00-00 (pre-planting application, 20 days before planting), (c) 00-90-00 (planting application), (d) 00-30-60 (split in a planting and a sidedress application 31 days after emergence), (e) 00-00-00* (control, without oat residue) and (f) 90-00-00* (pre-planting application, without oat residue). The N concentration and N content in oat residues were not affected during decomposition by N fertilisation. Most of the fertiliser NH4+-N was converted into NO3--N within 20 days after application. A significant decrease in NO3--N contents in the 0-4 cm layer was observed in all treatments between 40 and 60 days after the oat residue placement on the soil surface, suggesting the occurrence of N immobilisation in this period. Considering that most of the inorganic N was converted into NO3- and that no immobilisation of the pre planting fertiliser N occurred at the time of its application, it was possible to conclude that pre-planting applied N was prone to losses by leaching. On the other hand, with split N applications, maize plants showed N deficiency symptoms before sidedress application. Two indications for fertiliser-N management in no-till oat/maize sequential cropping could be suggested: (a) in case of split application, the sidedress should be earlier than 30 days after emergence, and (b) if integral application is preferred to save field operations, this should be done at planting.

Biomass and nutrient cycling in pure and mixed stands of native tree species in southeastern Bahia, Brazil

The objective of this paper is to study selected components of the nutrient cycle of pure and mixed stands of native forest species of Atlantic Forest in southeastern Brazil. Tree diameter, height, above-ground biomass, and nutrient content were determined in 22-year-old stands. Litterfall, litter decomposition, and nutrient concentration were evaluated from August 1994 to July 1995. The following species were studied: Peltogyne angustiflora, Centrolobium robustum, Arapatiella psilophylla, Sclerolobium chrysophyllum, Cordia trichotoma, Macrolobium latifolium. The litter of a natural forest and a 40-year-old naturally regenerated second-growth forest was sampled as well. The mixed-species outmatched pure stands in height, stem volume and total biomass (29.4 % more). The greatest amount of forest litter was observed in the natural forest (9.3 Mg ha-1), followed by the mixed-species stand (7.6 Mg ha-1) and secondary forest (7.3 Mg ha-1), and least litterfall was measured in the pure C. robustum stand (5.5 Mg ha-1). Litterfall seasonality varied among species in pure stands (CV from 44.7 to 91.4 %), unlike litterfall in the mixed-tree stand, where the variation was lower (CV 31.2 %). In the natural and second-growth forest, litterfall varied by 57.8 and 34.0 %, respectively. The annual rate of nutrient return via litterfall varied widely among forest ecosystems. Differences were detected between forest ecosystems in both the litter accumulation and quantity of litterlayer nutrients. The highest mean nutrient accumulation in above-ground biomass was observed in mixed-species stands. The total nutrient accumulation (N + P + K+ Ca + Mg) ranged from 0.97 to 1.93 kg tree-1 in pure stands, and from 1.21 to 2.63 kg tree-1 in mixed-species stands. Soil fertility under mixed-species stands (0-10 cm) was intermediate between the primary forest and pure-stand systems. The litterfall rate of native forest species in a mixed-species system is more constant, resulting in a more continuous decomposition rate. Consequently, both nutrient availability and quantity of organic matter in the soil are higher and the production system ecologically more sustainable.

Feminist Perspectives on Motherhood and Reproduction

Motherhood and reproduction have been at the core of the feminist discourse about women's rights ever since its onset. For the first and second feminist movements, the right to abortion and the public recognition of motherhood have been main issues in the discourse on reproduction. Since the last two dec- ades of the 20th century, the potentials of assisted reproductive technologies (ART) have opened up new venues of feminist discourse.In this paper we sketch the main feminist lines of argumentation regarding motherhood and reproduction since the 1970s, and we identify specific shifts in their recurrent issues. We argue that an essential contribution of feminism to the understanding of motherhood as a structuring category has been its insis- tence on the distinction between biological and social motherhood. Feminist discourse shows how ART has further decomposed biological motherhood and has altered the meaning of motherhood and reproduction. Feminist analysis maintains that despite the rhetoric of choice surrounding ART, these technolo- gies have not increased women's reproductive freedom. The decomposition of biological motherhood, the medical, legal, and commercial development of re- production, and the change in the social perception of motherhood have rather established new forms of control over female reproduction.

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.

On some numerical aspects of an active strain model in cardiac mechanics

We are interested in the development, implementation and testing of an orthotropic model for cardiac contraction based on an active strain decomposition. Our model addresses the coupling of a transversely isotropic mechanical description at the cell level, with an orthotropic constitutive law for incompressible tissue at the macroscopic level. The main differences with the active stress model are addressed in detail, and a finite element discretization using Taylor-Hood and MINI elements is proposed and illustrated with numerical examples.

Characterization of magmatic sulfur in the Aegean island arc by means of the δ34S values of fumarolic H2S, elemental S, and hydrothermal gypsum from Nisyros and Milos islands

A delta(34)S value of +6.3 +/- 1.5% was estimated for the rhyodacitic degassing magma present underneath the hydrothermal system of Nisyros, based on the S isotope ratios of H2S in fumarolic vapors. This value was estimated by modeling the irreversible water-rock mass transfers occurring during the generation of the hydrothermal liquid which separates these fumarolic vapors. The S isotope ratio of the rhyodacitic degassing magma of Nisyros is consistent with fractional crystallization of a parent basaltic magma with an initial delta(34)S value of +4% (+/-at least 1.5%). This positive value could be explained by mantle contamination due to by either transference of fluids derived from subducted materials or involvement of altered oceanic crust, whereas contribution of biogenic sulfides from sediments seems to be negligible or nil. This conclusion agrees with the lack of N-2 and CO2 from thermal decomposition of organic matter contained in subducted sediments, which is a characteristic of the whole Aegean arc system. Since hydrothermal S at Milos and Santorini has isotope ratios similar to those determined at Nisyros, it seems likely that common controlling processes are active throughout the Aegean island arc. (C) 2002 Elsevier, Science B.V. All rights reserved.

Mineralogical, K-Ar, stable and Sr isotope systematics of K-white micas during very low-grade metamorphism of limestones (Helvetic nappes, Western Switzerland)

Mineralogical, K-Ar, Rb-Sr and stable isotope analyses have been carried out on K-white micas from Helvetic Malm limestones in order to examine their evolution during very low- to low-grade Alpine metamorphism, associated with intense ductile deformation. Metamorphic temperatures were estimated al approximately 300-degrees-C from stable isotopes (quartz-calcite thermometry), occurrence of chloritoid, and `'epizonal'' illite crystallinity index. K-white micas consist of variable mixtures of 2M, phengite and muscovite, as revealed by detailed X-ray diffraction analyses using peak decomposition of the (060, 331) spectra. K-Ar apparent ages display a strong grain-size dependence in which mainly fine-grained size fractions (< 2 mum) record Alpine ages (37-15 Ma). However, these ages provide a relative rather than an absolute chronology of the diachronous Alpine metamorphic evolution of the Helvetic nappes. The resetting of the K-Ar isotopic system of K-white micas to Alpine metamorphic conditions reflects an apparent combination of crystallization/recrystallization and radiogenic Ar-40 diffusion loss. The oxygen isotope compositions of micas (+ 15 to + 22 parts per thousand) are intermediate between detrital and O-18-enriched values expected for micas neoformed within an abundant marine carbonate matrix. No isotopic equilibrium has been reached between calcite and micas. The variable depletion of hydrogen isotope compositions (- 126 to - 82 parts per thousand) is influenced by the interaction with organic matter under closed-system conditions. Organic matter, if not removed, may also represent a serious source of error in K-Ar age determination, by introducing radiogenic Ar-40 contamination. Sr-87/Sr-86 isotope ratios of micas range from 0.70879 to 0.70902 with one outlier at 0.71794. The low values reflect Sr exchange with calcite occurring during crystallization/recrystallization of micas under closed-system conditions.

The effect of substrate on high-temperature annealing of GaN epilayers: Si versus sapphire

We have studied the effects of rapid thermal annealing at 1300¿°C on GaN epilayers grown on AlN buffered Si(111) and on sapphire substrates. After annealing, the epilayers grown on Si display visible alterations with craterlike morphology scattered over the surface. The annealed GaN/Si layers were characterized by a range of experimental techniques: scanning electron microscopy, optical confocal imaging, energy dispersive x-ray microanalysis, Raman scattering, and cathodoluminescence. A substantial Si migration to the GaN epilayer was observed in the crater regions, where decomposition of GaN and formation of Si3N4 crystallites as well as metallic Ga droplets and Si nanocrystals have occurred. The average diameter of the Si nanocrystals was estimated from Raman scattering to be around 3¿nm. Such annealing effects, which are not observed in GaN grown on sapphire, are a significant issue for applications of GaN grown on Si(111) substrates when subsequent high-temperature processing is required.