A Random Evolution Inclusion of Subdifferential Type in Hilbert Spaces

In this paper we study a nonlinear evolution inclusion of subdifferential type in Hilbert spaces. The perturbation term is Hausdorff continuous in the state variable and has closed but not necessarily convex values. Our result is a stochastic generalization of an existence theorem proved by Kravvaritis and Papageorgiou in [6].

Generation and evolution of magma beneath the East Pacific Rise: Constraints from U-series disequilibrium and plagioclase-hosted melt inclusions

Mid-ocean ridge basalt (MORB) samples from the East Pacific Rise (EPR 12 degrees 50'N) were analyzed for U-series isotopes and compositions of plagioclase-hosted melt inclusions. The Ra-226 and Th-230 excesses are negatively correlated; the Ra-226 excess is positively correlated with Mg# and Sm/Nd, and is negatively correlated with La/Sm and Fe-8; the Th-230 excess is positively correlated with Fe-8 and La/Sm and is negatively correlated with Mg# and Sm/Nd. Interpretation of these correlations is critical for understanding the magmatic process. There are two models (the dynamic model and the "two-porosity" model) for interpreting these correlations, however, some crucial parameters used in these models are not ascertained. We propose instead a model to explain the U-series isotopic compositions based on the control of melt density variation. For melting either peridotite or the "marble-cake" mantle, the FeOt content, Th-230 excess and La/Sm ratio increases and Sm/Nd decreases with increasing pressure. A deep melt will evolve to a higher density and lower Mg# than a shallow melt, the former corresponds to a long residence time, which lowers the Ra-226 excess significantly. This model is supported by the existence of low Ra-226 excesses and high Th-230 excesses in MORBs having a high Fe-8 content and high density. The positive correlation of Ra-226 excess and magma liquidus temperature implies that the shallow melt is cooled less than the deep melt due to its low density and short residence time. The correlations among Fe-8, Ti-8 and Ca-8/Al-8 in plagioclase-hosted melt inclusions further prove that MORBs are formed from melts having a negative correlation in melting depths and degrees. The negative correlation of Ra-226 excess vs. chemical diversity index (standard deviation of Fe-8, Ti-8 and Ca-8/Al-8) of the melt inclusions is in accordance with the influence of a density-controlled magma residence time. We conclude that the magma density variation exerts significant control on residence time and U-series isotopic compositions. (c) 2010 Elsevier B.V. All rights reserved.

P-T-Fluid evolution and graphite deposition during retrograde metamorphism in Ribeira Fold Belt, SE Brazil: Oxygen fugacity, fluid inclusions and C-O-H isotopic evidence

Combined fluid inclusion (FI) microthermometry, Raman spectroscopy, X-ray diffraction, C-O-H isotopes and oxygen fugacities of granulites from central Ribeira Fold Belt, SE Brazil, provided the following results: i) Magnetite-Hematite fO(2) estimates range from 10(-11.5) bar (QFM + 1) to 10(-18.3) bar (QFM - 1) for the temperature range of 896 degrees C-656 degrees C, implying fO(2) decrease from metamorphic peak temperatures to retrograde conditions; ii) 5 main types of fluid inclusions were observed: a) CO(2) and CO(2)-N(2) (0-11 mol%) high to medium density (1.01-0.59 g/cm(3)) FI; b) CO(2) and CO(2)-N(2) (0-36 mol%) low density (0.19-0.29 g/cm(3)) FI; c) CO(2) (94-95 mol%)-N(2) (3 mol%)-CH(4) (2-3 mol%)-H(2)O (water phi(v) (25 degrees C) = 0.1) FI; d) low-salinity H(2)O-CO(2) FI; and e) late low-salinity H(2)O FI; iii) Raman analyses evidence two graphite types in khondalites: an early highly ordered graphite (T similar to 450 degrees C) overgrown by a disordered kind (T similar to 330 degrees C); iv) delta(18)O quartz results of 10.3-10.7 parts per thousand, imply high-temperature CO(2) delta(18)O values of 14.4-14.8 parts per thousand, suggesting the involvement of a metamorphic fluid, whereas lower temperature biotite delta(18)O and delta D results of 7.5-8.5 parts per thousand and -54 to -67 parts per thousand respectively imply H(2)O delta(18)O values of 10-11 parts per thousand and delta D(H2O) of -23 to -36 parts per thousand suggesting delta(18)O depletion and increasing fluid/rock ratio from metamorphic peak to retrograde conditions. Isotopic results are compatible with low-temperature H(2)O influx and fO(2) decrease that promoted graphite deposition in retrograde granulites, simultaneous with low density CO(2), CO(2)-N(2) and CO(2)-N(2)-CH(4)-H(2)O fluid inclusions at T = 450-330 degrees C. Graphite delta(13)C results of -10.9 to -11.4 parts per thousand imply CO(2) delta(13)C values of -0.8 to -1.3 parts per thousand suggesting decarbonation of Cambrian marine carbonates with small admixture of lighter biogenic or mantle derived fluids. Based on these results, it is suggested that metamorphic fluids from the central segment of Ribeira Fold Belt evolved to CO(2)-N(2) fluids during granulitic metamorphism at high fO(2), followed by rapid pressure drop at T similar to 400-450 degrees C during late exhumation that caused fO(2) reduction induced by temperature decrease and water influx, turning carbonic fluids into CO(2)-H(2)O (depleting biotite delta(18)O and delta D values), and progressively into H(2)O. When fO(2) decreased substantially by mixture of carbonic and aqueous fluids, graphite deposited forming khondalites. (C) 2010 Elsevier Ltd. All rights reserved.

Magmatic evolution of topaz-bearing granite stocks within the Wiborg rapakivi granite batholith

This study brings new insights into the magmatic evolution of natural F-enriched peraluminous granitic systems. The Artjärvi, Sääskjärvi and Kymi granite stocks within the 1.64 Ga Wiborg rapakivi granite batholith have been investigated by petrographic, geochemical, experimental and melt inclusion methods. These stocks represent late-stage leucocratic and weakly peraluminous intrusive phases typical of rapakivi granites worldwide. The Artjärvi and Sääskjärvi stocks are multiphase intrusions in which the most evolved phase is topaz granite. The Kymi stock contains topaz throughout and has a well-developed zoned structure, from the rim to the center: stockscheider pegmatite equigranular topaz granite porphyritic topaz granite. Geochemically the topaz granites are enriched in F, Li, Be, Ga, Rb, Sn and Nb and depleted in Mg, Fe, Ti, Ba, Sr, Zr and Eu. The anomalous geochemistry and mineralogy of the topaz granites are essentially magmatic in origin; postmagmatic reactions have only slightly modified the compositions. The Kymi equigranular topaz granite shows the most evolved character, and the topaz granites at Artjärvi and Sääskjärvi resemble the less evolved porphyritic topaz granite of the Kymi stock. Stockscheiders are found at the roof contacts of the Artjärvi and Kymi stocks. The stockscheider at Artjärvi is composed of biotite-rich schlieren and pegmatite layers parallel to the contact. The schlieren layering is considered to have formed by velocity-gradient sorting mechanism parallel to the flow, which led to the accumulation of mafic minerals along the upper contact of the topaz granite. Cooling and contraction of the topaz granite formed fractures parallel to the roof contact and residual pegmatite magmas were injected along the fractures and formed the pegmatite layers. The zoned structure of the Kymi stock is the result of intrusion of highly evolved residual melt from deeper parts of the magma chamber along the fractured contact between the porphyritic granite crystal mush and country rock. The equigranular topaz granite and marginal pegmatite (stockscheider) crystallized from this evolved melt. Phase relations of the Kymi equigranular topaz granite have been investigated utilizing crystallization experiments at 100 to 500 MPa as a function of water activity and F content. Fluorite and topaz can crystallize as liquidus phases in F-rich peraluminous systems, but the F content of the melt should exceed 2.5 - 3.0 wt % to facilitate crystallization of topaz. In peraluminous F-bearing melts containing more than 1 wt % F, topaz and muscovite are expected to be the first F-bearing phases to crystallize at high pressure, whereas fluorite and topaz should crystallize first at low pressure. Overall, the saturation of fluorite and topaz follows the reaction: CaAl2Si2O8 (plagioclase) + 2[AlF3]melt = CaF2 (fluorite) + 2Al2SiO4F2 (topaz). The obtained partition coefficient for F between biotite and glass D(F)Bt/glass is 1.89 to 0.80 (average 1.29) and can be used as an empirical fluormeter to determine the F content of coexisting melts. In order to study the magmatic evolution of the Kymi stock, crystallized melt inclusions in quartz and topaz grains in the porphyritic and the equigranular topaz granites and the marginal pegmatite were rehomogenized and analyzed. The homogenization conditions for the melt inclusions from the granites were 700 °C, 300 MPa, and 24 h, and for melt inclusions from the pegmatite, 700 °C, 100 MPa, and 24/96 h. The majority of the melt inclusions is chemically similar to the bulk rocks (excluding H2O content), but a few melt inclusions in the equigranular granite show clearly higher F and low K2O contents (on average 11.6 wt % F, 0.65 wt % K2O). The melt inclusion compositions indicate coexistence of two melt fractions, a prevailing peraluminous and a very volatile-rich, possibly peralkaline. Combined petrological, experimental and melt inclusion studies of the Kymi equigranular topaz granite indicate that plagioclase was the liquidus phase at nearly water-saturated (fluid-saturated) conditions and that the F content of the melt was at least 2 wt %. The early crystallization of biotite and the presence of muscovite in crystallization experiments at 200 MPa contrasts with the late-stage crystallization of biotite and the absence of muscovite in the equigranular granite, indicating that crystallization pressure may have been lower than 200 MPa for the granite.

Ductile deformation of single inclusions in simple shear with a finite-strain hyperelastoviscoplastic rheology

We examine the 2D plane-­strain deformation of initially round, matrix-­bonded, deformable single inclusions in isothermal simple shear using a recently introduced hyperelastoviscoplastic rheology. The broad parameter space spanned by the wide range of effective viscosities, yield stresses, relaxation times, and strain rates encountered in the ductile lithosphere is explored systematically for weak and strong inclusions, the effective viscosity of which varies with respect to the matrix. Most inclusion studies to date focused on elastic or purely viscous rheologies. Comparing our results with linear-­viscous inclusions in a linear-­viscous matrix, we observe significantly different shape evolution of weak and strong inclusions over most of the relevant parameter space. The evolution of inclusion inclination relative to the shear plane is more strongly affected by elastic and plastic contributions to rheology in the case of strong inclusions. In addition, we found that strong inclusions deform in the transient viscoelastic stress regime at high Weissenberg numbers (≥0.01) up to bulk shear strains larger than 3. Studies using the shapes of deformed objects for finite-­strain analysis or viscosity-­ratio estimation should establish carefully which rheology and loading conditions reflect material and deformation properties. We suggest that relatively strong, deformable clasts in shear zones retain stored energy up to fairly high shear strains. Hence, purely viscous models of clast deformation may overlook an important contribution to the energy budget, which may drive dissipation processes within and around natural inclusions.

Solubility and Diffusivity of Water in Lunar Basalt, and Chemical Zonation in Olivine-Hosted Melt Inclusions

I report the solubility and diffusivity of water in lunar basalt and an iron-free basaltic analogue at 1 atm and 1350 °C. Such parameters are critical for understanding the degassing histories of lunar pyroclastic glasses. Solubility experiments have been conducted over a range of fO₂ conditions from three log units below to five log units above the iron-wüstite buffer (IW) and over a range of pH₂/pH₂O from 0.03 to 24. Quenched experimental glasses were analyzed by Fourier transform infrared spectroscopy (FTIR) and secondary ionization mass spectrometry (SIMS) and were found to contain up to ~420 ppm water. Results demonstrate that, under the conditions of our experiments: (1) hydroxyl is the only H-bearing species detected by FTIR; (2) the solubility of water is proportional to the square root of pH₂O in the furnace atmosphere and is independent of fO₂ and pH₂/pH₂O; (3) the solubility of water is very similar in both melt compositions; (4) the concentration of H₂ in our iron-free experiments is <3 ppm, even at oxygen fugacities as low as IW-2.3 and pH₂/pH₂O as high as 24; and (5) SIMS analyses of water in iron-rich glasses equilibrated under variable fO₂ conditions can be strongly influenced by matrix effects, even when the concentrations of water in the glasses are low. Our results can be used to constrain the entrapment pressure of the lunar melt inclusions of Hauri et al. (2011).

Diffusion experiments were conducted over a range of fO₂ conditions from IW-2.2 to IW+6.7 and over a range of pH₂/pH₂O from nominally zero to ~10. The water concentrations measured in our quenched experimental glasses by SIMS and FTIR vary from a few ppm to ~430 ppm. Water concentration gradients are well described by models in which the diffusivity of water (D^*_water) is assumed to be constant. The relationship between D^*_water and water concentration is well described by a modified speciation model (Ni et al. 2012) in which both molecular water and hydroxyl are allowed to diffuse. The success of this modified speciation model for describing our results suggests that we have resolved the diffusivity of hydroxyl in basaltic melt for the first time. Best-fit values of D^*_water for our experiments on lunar basalt vary within a factor of ~2 over a range of pH₂/pH₂O from 0.007 to 9.7, a range of fO₂ from IW-2.2 to IW+4.9, and a water concentration range from ~80 ppm to ~280 ppm. The relative insensitivity of our best-fit values of D^*_water to variations in pH₂ suggests that H₂ diffusion was not significant during degassing of the lunar glasses of Saal et al. (2008). D^*_water during dehydration and hydration in H₂/CO₂ gas mixtures are approximately the same, which supports an equilibrium boundary condition for these experiments. However, dehydration experiments into CO₂ and CO/CO₂ gas mixtures leave some scope for the importance of kinetics during dehydration into H-free environments. The value of D^*_water chosen by Saal et al. (2008) for modeling the diffusive degassing of the lunar volcanic glasses is within a factor of three of our measured value in our lunar basaltic melt at 1350 °C.

In Chapter 4 of this thesis, I document significant zonation in major, minor, trace, and volatile elements in naturally glassy olivine-hosted melt inclusions from the Siqueiros Fracture Zone and the Galapagos Islands. Components with a higher concentration in the host olivine than in the melt (MgO, FeO, Cr₂O₃, and MnO) are depleted at the edges of the zoned melt inclusions relative to their centers, whereas except for CaO, H₂O, and F, components with a lower concentration in the host olivine than in the melt (Al₂O₃, SiO₂, Na₂O, K₂O, TiO₂, S, and Cl) are enriched near the melt inclusion edges. This zonation is due to formation of an olivine-depleted boundary layer in the adjacent melt in response to cooling and crystallization of olivine on the walls of the melt inclusions concurrent with diffusive propagation of the boundary layer toward the inclusion center.

Concentration profiles of some components in the melt inclusions exhibit multicomponent diffusion effects such as uphill diffusion (CaO, FeO) or slowing of the diffusion of typically rapidly diffusing components (Na₂O, K₂O) by coupling to slow diffusing components such as SiO₂ and Al₂O₃. Concentrations of H2O and F decrease towards the edges of some of the Siqueiros melt inclusions, suggesting either that these components have been lost from the inclusions into the host olivine late in their cooling histories and/or that these components are exhibiting multicomponent diffusion effects.

A model has been developed of the time-dependent evolution of MgO concentration profiles in melt inclusions due to simultaneous depletion of MgO at the inclusion walls due to olivine growth and diffusion of MgO in the melt inclusions in response to this depletion. Observed concentration profiles were fit to this model to constrain their thermal histories. Cooling rates determined by a single-stage linear cooling model are 150–13,000 °C hr^-1 from the liquidus down to ~1000 °C, consistent with previously determined cooling rates for basaltic glasses; compositional trends with melt inclusion size observed in the Siqueiros melt inclusions are described well by this simple single-stage linear cooling model. Despite the overall success of the modeling of MgO concentration profiles using a single-stage cooling history, MgO concentration profiles in some melt inclusions are better fit by a two-stage cooling history with a slower-cooling first stage followed by a faster-cooling second stage; the inferred total duration of cooling from the liquidus down to ~1000 °C is 40 s to just over one hour.

Based on our observations and models, compositions of zoned melt inclusions (even if measured at the centers of the inclusions) will typically have been diffusively fractionated relative to the initially trapped melt; for such inclusions, the initial composition cannot be simply reconstructed based on olivine-addition calculations, so caution should be exercised in application of such reconstructions to correct for post-entrapment crystallization of olivine on inclusion walls. Off-center analyses of a melt inclusion can also give results significantly fractionated relative to simple olivine crystallization.

All melt inclusions from the Siqueiros and Galapagos sample suites exhibit zoning profiles, and this feature may be nearly universal in glassy, olivine-hosted inclusions. If so, zoning profiles in melt inclusions could be widely useful to constrain late-stage syneruptive processes and as natural diffusion experiments.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Petrogenesis of Basaltic Andesites at Tofua Volcano

Tofua volcano is situated midway along the Tonga oceanic arc and has undergone two phases of ignimbrite-forming activity. The eruptive products are almost entirely basaltic andesites (52 center dot 5-57 wt % SiO2) with the exception of a volumetrically minor pre-caldera dacite. The suite displays a strong tholeiitic trend with K2O <1 wt %. Phenocryst assemblages typically comprise plagioclase + clinopyroxene +/- orthopyroxene with microlites of Ti-magnetite. Olivine (Fo(83-88)) is rare and believed to be dominantly antecrystic. An increase in the extent and frequency of reverse zoning in phenocrysts, sieve-textured plagioclase and the occurrence of antecrystic phases in post-caldera lavas record a shift to dynamic conditions, allowing the interaction of magma batches that were previously distinct. Pyroxene thermobarometry suggests crystallization at 950-1200 degrees C and 0 center dot 8-1 center dot 8 kbar. Volatile measurements of glassy melt inclusions indicate a maximum H2O content of 4 center dot 16 wt % H2O, and CO2-H2O saturation curves indicate that crystallization occurred at two levels, at depths of 4-5 center dot 5 km and 1 center dot 5-2 center dot 5 km. Major and trace element models suggest that the compositions of the majority of the samples represent a differentiation trend whereby the dacite was produced by 65% fractional crystallization of the most primitive basaltic andesite. Trace element models suggest that the sub-arc mantle source is the residuum of depleted Indian mid-ocean ridge basalt mantle (IDMM-1% melt), whereas radiogenic isotope data imply addition of 0 center dot 2% average Tongan sediment melt and a fluid component derived from the subducted altered Pacific oceanic crust. A horizontal array on the U-Th equiline diagram and Ra excesses of up to 500% suggest fluid addition to the mantle wedge within the last few thousand years. Time-integrated (Ra-226/Th-230) vs Sr/Th and Ba/Th fractionation models imply differentiation timescales of up to 4500 years for the dacitic magma compositions at Tofua.

Possible Host Adaptation as an Evolution Factor of Cowpea aphid-borne mosaic virus Deduced by Coat Protein Gene Analysis

Cowpea aphid-borne mosaic virus (CABMV) causes major diseases in cowpea and passion flower plants in Brazil and also in other countries. CABMV has also been isolated from leguminous species including, Cassia hoffmannseggii, Canavalia rosea, Crotalaria juncea and Arachis hypogaea in Brazil. The virus seems to be adapted to two distinct families, the Passifloraceae and Fabaceae. Aiming to identify CABMV and elucidate a possible host adaptation of this virus species, isolates from cowpea, passion flower and C.hoffmannseggii collected in the states of Pernambuco and Rio Grande do Norte were analysed by sequencing the complete coat protein genes. A phylogenetic tree was constructed based on the obtained sequences and those available in public databases. Major Brazilian isolates from passion flower, independently of the geographical distances among them, were grouped in three different clusters. The possible host adaptation was also observed in fabaceous-infecting CABMV Brazilian isolates. These host adaptations possibly occurred independently within Brazil, so all these clusters belong to a bigger Brazilian cluster. Nevertheless, African passion flower or cowpea-infecting isolates formed totally different clusters. These results showed that host adaptation could be one factor for CABMV evolution, although geographical isolation is a stronger factor.

Mineral chemistry and structural relationships of inclusions in diamond samples

Diamant ist das härteste Mineral – und dazu ein Edelstein -, das unter höchstem Druck und hohen Temperaturen in tiefen kontinentalen Regionen der Erde kristallisiert. Die Mineraleinschlüsse in Diamanten werden durch die physikalische Stabilität und chemische Beständigkeit der umgebenden – eigentlich metastabilen -Diamant-Phase geschützt. Aufgrund der koexistierenden Phasenkombination ermöglichen sie, die Mineral-Entwicklung zu studieren, während deren der Einschlüssen und die Diamanten kristallisierten. rnDie Phasenkombinationen von Diamant und Chrom-Pyrop, Chrom-Diopsid, Chromit, Olivin, Graphit und Enstatit nebeneinander (teilweise in Berührungsexistenz) mit Chrom-Pyrop Einschlüssen wurden von neunundzwanzig Diamant-Proben von sechs Standorten in Südafrika (Premier, Koffiefontein, De Beers Pool, Finsch, Venetia und Koingnaas Minen) und Udachnaya (Sibirien/Russland) identifiziert und charakterisiert. Die Mineraleinschlüsse weisen z.T. kubo-oktaedrische Form auf, die unabhängig von ihren eigenen Kristallsystemen ausgebildet werden können. Das bedeutet, dass sie syngenetische Einschlüsse sind, die durch die sehr hohe Formenergie des umgebenden Diamanten morphologisch unter Zwang stehen. Aus zweidiemnsionalen Messungen der ersten Ordnung von charakteristischen Raman-Banden lassen sich relative Restdrucke in Diamanten zwischen Diamant und Einschlussmineral gewinnen; sie haben charakteristische Werte von ca. 0,4 bis 0,9 GPa um Chrom-Pyrop-Einschlüsse, 0,6 bis 2,0 GPa um Chrom-Diopsid-Einschlüsse, 0,3 bis 1,2 GPa um Olivin-Einschlüsse, 0,2 bis 1,0 GPa um Chromit-Einschlüsse, beziehungsweise 0,5 GPa um Graphit Einschlüsse.rnDie kristallstrukturellen Beziehung von Diamanten und ihren monomineralischen Einschlüssen wurden mit Hilfe der Quantifizierung der Winkelkorrelationen zwischen der [111] Richtung von Diamanten und spezifisch ausgewählten Richtungen ihrer mineralischen Einschlüsse untersucht. Die Winkelkorrelationen zwischen Diamant [111] und Chrom-Pyrop [111] oder Chromit [111] zeigen die kleinsten Verzerrungen von 2,2 bis zu 3,4. Die Chrom-Diopsid- und Olivin-Einschlüsse zeigen die Missorientierungswerte mit Diamant [111] bis zu 10,2 und 12,9 von Chrom-Diopsid [010] beziehungsweise Olivin [100].rnDie chemische Zusammensetzung von neun herausgearbeiteten (orientiertes Anschleifen) Einschlüssen (drei Chrom-Pyrop-Einschlüsse von Koffiefontein-, Finsch- und Venetia-Mine (zwei von drei koexistieren nebeneinander mit Enstatit), ein Chromit von Udachnaya (Sibirien/Russland), drei Chrom-Diopside von Koffiefontein, Koingnaas und Udachnaya (Sibirien/Russland) und zwei Olivin Einschlüsse von De Beers Pool und Koingnaas) wurden mit Hilfe EPMA und LA-ICP-MS analysiert. Auf der Grundlage der chemischen Zusammensetzung können die Mineraleinschlüsse in Diamanten in dieser Arbeit der peridotitischen Suite zugeordnet werden.rnDie Geothermobarometrie-Untersuchungen waren aufgrund der berührenden Koexistenz von Chrom-Pyrop- und Enstatit in einzelnen Diamanten möglich. Durchschnittliche Temperaturen und Drücke der Bildung sind mit ca. 1087 (± 15) C, 5,2 (± 0,1) GPa für Diamant DHK6.2 von der Koffiefontein Mine beziehungsweise ca. 1041 (± 5) C, 5,0 (± 0,1) GPa für Diamant DHF10.2 von der Finsch Mine zu interpretieren.rn