Interpretation of fluid inclusions in quartz deformed by weak ductile shearing: reconstruction of differential stress magnitudes and pre-deformation fluid properties

A well developed theoretical framework is available in which paleofluid properties, such as chemical composition and density, can be reconstructed from fluid inclusions in minerals that have undergone no ductile deformation. The present study extends this framework to encompass fluid inclusions hosted by quartz that has undergone weak ductile deformation following fluid entrapment. Recent experiments have shown that such deformation causes inclusions to become dismembered into clusters of irregularly shaped relict inclusions surrounded by planar arrays of tiny, new-formed (neonate) inclusions. Comparison of the experimental samples with a naturally sheared quartz vein from Grimsel Pass, Aar Massif, Central Alps, Switzerland, reveals striking similarities. This strong concordance justifies applying the experimentally derived rules of fluid inclusion behaviour to nature. Thus, planar arrays of dismembered inclusions defining cleavage planes in quartz may be taken as diagnostic of small amounts of intracrystalline strain. Deformed inclusions preserve their pre-deformation concentration ratios of gases to electrolytes, but their H2O contents typically have changed. Morphologically intact inclusions, in contrast, preserve the pre-deformation composition and density of their originally trapped fluid. The orientation of the maximum principal compressive stress (σ1σ1) at the time of shear deformation can be derived from the pole to the cleavage plane within which the dismembered inclusions are aligned. Finally, the density of neonate inclusions is commensurate with the pressure value of σ1σ1 at the temperature and time of deformation. This last rule offers a means to estimate magnitudes of shear stresses from fluid inclusion studies. Application of this new paleopiezometer approach to the Grimsel vein yields a differential stress (σ1–σ3σ1–σ3) of ∼300 MPa∼300 MPa at View the MathML source390±30°C during late Miocene NNW–SSE orogenic shortening and regional uplift of the Aar Massif. This differential stress resulted in strain-hardening of the quartz at very low total strain (<5%<5%) while nearby shear zones were accommodating significant displacements. Further implementation of these experimentally derived rules should provide new insight into processes of fluid–rock interaction in the ductile regime within the Earth's crust.

Experimental Characterization and Quantification of Cement-Bentonite Interaction using Core Infiltration Techniques Coupled with X-ray Tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel canisters and as drift seals. Sand/bentonite (s/b) is foreseen as backfill material of access galleries or as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore-water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predicted significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this thesis was to characterize and quantify the cement/bentonite interactions both spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used to perform X-ray computed tomography (CT) scans without interruption of running experiments. CT scans allowed tracking the evolution of the reaction plume and changes in core volume/diameter/density during the experiments. In total 4 core infiltration experiments were carried out for this study with the compacted and saturated cores consisting of MX-80 bentonite and sand/MX-80 bentonite mixture (s/b; 65/35%). Two different high-pH cementitious pore-fluids were infiltrated: a young (early) ordinary Portland cement pore-fluid (APWOPC; K+–Na+–OH-; pH 13.4; ionic strength 0.28 mol/kg) and a young ‘low-pH’ ESDRED shotcrete pore-fluid (APWESDRED; Ca2+–Na+–K+–formate; pH 11.4; ionic strength 0.11 mol/kg). The experiments lasted between 1 and 2 years. In both bentonite experiments, the hydraulic conductivity was strongly reduced after switching to high-pH fluids, changing eventually from an advective to a diffusion-dominated transport regime. The reduction was mainly induced by mineral precipitation and possibly partly also by high ionic strength pore-fluids. Both bentonite cores showed a volume reduction and a resulting transient flow in which pore-water was squeezed out during high-pH infiltration. The outflow chemistry was characterized by a high ionic strength, while chloride in the initial pore water got replaced as main anionic charge carrier by sulfate, originating from gypsum dissolution. The chemistry of the high-pH fluids got strongly buffered by the bentonite, consuming hydroxide and in case of APWESDRED also formate. Hydroxide got consumed by mineral reactions (saponite and possibly talc and brucite precipitation), while formate being affected by bacterial degradation. Post-mortem analysis showed reaction zones near the inlet of the bentonite core, characterized by calcium and magnesium enrichment, consisting predominately of calcite and saponite, respectively. Silica got enriched in the outflow, indicating dissolution of silicate-minerals, identified as preferentially cristobalite. In s/b, infiltration of APWOPC reduced the hydraulic conductivity strongly, while APWESDRED infiltration had no effect. The reduction was mainly induced by mineral precipitation and probably partly also by high ionic strength pore-fluids. Not clear is why the observed mineral precipitates in the APWESDRED experiment had no effect on the fluid flow. Both s/b cores showed a volume expansion along with decreasing ionic strengths of the outflow, due to mineral reactions or in case of APWESDRED infiltration also mediated by microbiological activity, consuming hydroxide and formate, respectively. The chemistry of the high-pH fluids got strongly buffered by the s/b. In the case of APWESDRED infiltration, formate reached the outflow only for a short time, followed by enrichment in acetate, indicating most likely biological activity. This was in agreement to post-mortem analysis of the core, observing black spots on the inflow surface, while the sample had a rotten-egg smell indicative of some sulfate reduction. Post-mortem analysis showed further in both cores a Ca-enrichment in the first 10 mm of the core due to calcite precipitation. Mg-enrichment was only observed in the APWOPC experiment, originating from newly formed saponite. Silica got enriched in the outflow of both experiments, indicating dissolution of silicate-minerals, identified in the OPC experiment as cristobalite. The experiments attested an effective buffering capacity for bentonite and s/b, a progressing coupled hydraulic-chemical sealing process and also the preservation of the physical integrity of the interface region in this setup with a total pressure boundary condition on the core sample. No complete pore-clogging was observed but the hydraulic conductivity got rather strongly reduced in 3 experiments, explained by clogging of the intergranular porosity (macroporosity). Such a drop in hydraulic conductivity may impact the saturation time of the buffer in a nuclear waste repository, although the processes and geometry will be more complex in repository situation.

Oxygen isotope record of oceanic and high-pressure metasomatism: a P-T-time-fluid path for the Monviso eclogites (Italy)

Fluids are considered a fundamental agent for chemical exchanges between different rock types in the subduction system. Constraints on the sources and pathways of subduction fluids thus provide crucial information to reconstruct subduction processes. The Monviso ophiolitic sequence is composed of mafic, ultramafic and minor sediments that have been subducted to ~80 km depth. In this sequence, both localized fluid flow and channelized fluids along major shear zones have been documented. We investigate the timing and source of the fluids that affected the dominant mafic rocks using microscale U-Pb dating of zircon and oxygen isotope analysis of mineral zones (garnet, zircon and antigorite) in high pressure rocks with variable degree of metasomatic modification. In mafic eclogites, Jurassic zircon cores are the only mineralogical relicts of the protolith gabbros and retain δ18O values of 4.5–6 ‰, typical of mantle melts. Garnet and metamorphic zircon that grew during prograde to peak metamorphism display low δ18O values between 0.2 and 3.8 ‰, which are likely inherited from high-temperature alteration of the protolith on the sea floor. This is corroborated by δ18O values of 3.0 and 3.6 ‰ in antigorite from surrounding serpentinites. In metasomatised eclogites within the Lower Shear Zone, garnet rim formed at the metamorphic peak shows a shift to higher δ18O up to 6‰. The age of zircons in high-pressure veins and metasomatised eclogites constrains the timing of fluid flow at high pressure at around 45–46 Ma. Although the oxygen data do not contradict previous reports of interaction with serpentinite-derived fluids, the shift to isotopically heavier oxygen compositions requires contribution from sediment-derived fluids. The scarcity of metasediments in the Monviso sequence suggests that such fluids were concentrated and fluxed along the Lower Shear Zone in a sufficient amount to modify the oxygen composition of the eclogitic minerals.

Modeled Interaction of Comet 67P/Churyumov-Gerasimenko with the Solar Wind Inside 2 AU

Periodic comets move around the Sun on elliptical orbits. As such comet 67P/Churyumov-Gerasimenko (hereafter 67P) spends a portion of time in the inner solar system where it is exposed to increased solar insolation. Therefore given the change in heliocentric distance, in case of 67P from aphelion at 5.68 AU to perihelion at ~1.24 AU, the comet’s activity—the production of neutral gas and dust—undergoes significant variations. As a consequence, during the inbound portion, the mass loading of the solar wind increases and extends to larger spatial scales. This paper investigates how this interaction changes the character of the plasma environment of the comet by means of multifluid MHD simulations. The multifluid MHD model is capable of separating the dynamics of the solar wind ions and the pick-up ions created through photoionization and electron impact ionization in the coma of the comet. We show how two of the major boundaries, the bow shock and the diamagnetic cavity, form and develop as the comet moves through the inner solar system. Likewise for 67P, although most likely shifted back in time with respect to perihelion passage, this process is reversed on the outbound portion of the orbit. The presented model herein is able to reproduce some of the key features previously only accessible to particle-based models that take full account of the ions’ gyration. The results shown herein are in decent agreement to these hybrid-type kinetic simulations.

Focused fluid transfer through the mantle above subduction zones

Volcanic arcs above subduction zones are enriched in volatiles and fluid-mobile elements with respect to mid-oceanic ridge basalts. There is general consensus that this particular subduction zone signature is generated by fluid-induced extraction of these elements from subducted oceanic crust and its sedimentary cover. However, how these fluids are transferred through the mantle wedge to the locus of partial melting and what modification the fluids will experience is unresolved. Here we investigate the interaction of slab fluids with the mantle wedge through a series of high-pressure experiments. We explore two end-member processes of focused and porous reactive flow of hydrous slab melts through the mantle. Transfer by porous flow leads to the formation of hydrous minerals that sequester fluid-mobile elements and residual fluids characterized by trace element patterns inconsistent with typical arc lavas. In contrast, no hydrous minerals are formed in the reaction zone of experiments mimicking focused flow, and the typical trace element signature acquired during fluid extraction from the slab is preserved, indicating that this is an efficient process for element transfer through the mantle wedge.

Microthermometry and Raman analysis of fluid inclusions in different generations of quartz veins of the Athabasca Basin

The Athabasca Basin (Canada) contains the highest grade unconformity-type uranium deposits in the world. Underlying the Athabasca Group sedimentary rocks of the Dufferin Lake zone are variably graphitic pelitic schists (VGPS), altered to chlorite and hematite (Red/Green Zone: RGZ), and locally bleached near the unconformity during paleoweathering and/or later fluid interaction, leading to a loss of graphite near the unconformity. Fluid inclusions were examined in different generations of quartz veins, using microthermometry and Raman analysis, to characterize and compare the different fluids that interacted with the RGZ and the VGPS. In the VGPS, CH4-, N2- and CO2-rich fluids circulated. CH4- and N2-rich fluids could be the result of the breakdown of graphite to CH4/CO2, whereas N2-rich fluid is interpreted to be the result of breakdown of feldspars/micas to NH4+/N2. In the RGZ, highly saline fluids interpreted to be basinally derived have been recorded. The circulation of the two types of fluids (carbonic and brines) occurred at two different distinct events: 1) during the retrograde metamorphism of the basement rocks before the deposition of the Athabasca Basin for the carbonic fluids, and 2) after the deposition of the Athabasca Basin for the brines. Thus, in addition to possibly be related to graphite depletion in the RGZ, the brines can be linked to uranium mineralization.

Sediment and pore fluid chemistry for IODP Sites 311-U1325 and 311-U1329

Analytical challenges in obtaining high quality measurements of rare earth elements (REEs) from small pore fluid volumes have limited the application of REEs as deep fluid geochemical tracers. Using a recently developed analytical technique, we analyzed REEs from pore fluids collected from Sites U1325 and U1329, drilled on the northern Cascadia margin during the Integrated Ocean Drilling Program (IODP) Expedition 311, to investigate the REE behavior during diagenesis and their utility as tracers of deep fluid migration. These sites were selected because they represent contrasting settings on an accretionary margin: a ponded basin at the toe of the margin, and the landward Tofino Basin near the shelf's edge. REE concentrations of pore fluid in the methanogenic zone at Sites U1325 and U1329 correlate positively with concentrations of dissolved organic carbon (DOC) and alkalinity. Fractionations across the REE series are driven by preferential complexation of the heavy REEs. Simultaneous enrichment of diagenetic indicators (DOC and alkalinity) and of REEs (in particular the heavy elements Ho to Lu), suggests that the heavy REEs are released during particulate organic carbon (POC) degradation and are subsequently chelated by DOC. REE concentrations are greater at Site U1325, a site where shorter residence times of POC in sulfate-bearing redox zones may enhance REE burial efficiency within sulfidic and methanogenic sediment zones where REE release ensues. Cross-plots of La concentrations versus Cl, Li and Sr delineate a distinct field for the deep fluids (z > 75 mbsf) at Site U1329, and indicate the presence of a fluid not observed at the other sites drilled on the Cascadia margin. Changes in REE patterns, the presence of a positive Eu anomaly, and other available geochemical data for this site suggest a complex hydrology and possible interaction with the igneous Crescent Terrane, located east of the drilled transect.

A laboratory-scale experiment and a numerical simulation of unusual spiral plumes in a High-Prandtl-number fluid, links supplementary material

We experimentally and numerically investigated the generation of plumes from a local heat source (LHS) and studied the interaction of these plumes with cellular convective motion (CCM) in a rectangular cavity filled with silicon oil at a Prandtl number (Pr) of approximately two thousand. The LHS is generated using a 0.2-W green laser beam. A roll-type CCM is generated by vertically heating one side of the cavity. The CCM may lead to the formation of an unusual spiral convective plume that resembles a vertical Archimedes spiral. A similar plume is obtained in a direct numerical simulation. We discuss the physical mechanism for the formation of a spiral plume and the application of the results to mantle convection problems. We also estimate the Reynolds (Re) and Rayleigh (Ra) numbers and apply self-similarity theory to convection in the Earth's mantle. Spiral plumes can be used to interpret mantle tomography results over the last decade.

Particle methods parallel implementations by GP-GPU strategies

This paper outlines the problems found in the parallelization of SPH (Smoothed Particle Hydrodynamics) algorithms using Graphics Processing Units. Different results of some parallel GPU implementations in terms of the speed-up and the scalability compared to the CPU sequential codes are shown. The most problematic stage in the GPU-SPH algorithms is the one responsible for locating neighboring particles and building the vectors where this information is stored, since these specific algorithms raise many dificulties for a data-level parallelization. Because of the fact that the neighbor location using linked lists does not show enough data-level parallelism, two new approaches have been pro- posed to minimize bank conflicts in the writing and subsequent reading of the neighbor lists. The first strategy proposes an efficient coordination between CPU-GPU, using GPU algorithms for those stages that allow a straight forward parallelization, and sequential CPU algorithms for those instructions that involve some kind of vector reduction. This coordination provides a relatively orderly reading of the neighbor lists in the interactions stage, achieving a speed-up factor of x47 in this stage. However, since the construction of the neighbor lists is quite expensive, it is achieved an overall speed-up of x41. The second strategy seeks to maximize the use of the GPU in the neighbor's location process by executing a specific vector sorting algorithm that allows some data-level parallelism. Al- though this strategy has succeeded in improving the speed-up on the stage of neighboring location, the global speed-up on the interactions stage falls, due to inefficient reading of the neighbor vectors. Some changes to these strategies are proposed, aimed at maximizing the computational load of the GPU and using the GPU texture-units, in order to reach the maximum speed-up for such codes. Different practical applications have been added to the mentioned GPU codes. First, the classical dam-break problem is studied. Second, the wave impact of the sloshing fluid contained in LNG vessel tanks is also simulated as a practical example of particle methods

Polymer-solvent interaction parameters of SBS rubbers by inverse gas chromatography measurements

The solubility parameters of two SBS commercial rubbers with different structures (lineal and radial), and with slightly different styrene content have been determined by inverse gas chromatography technique. The Flory–Huggins interaction parameters of several polymer–solvent mixtures have also been calculated. The influence of the polymer composition, the solvent molecular weight and the temperature over these parameters have been discussed; besides, these parameters have been compared with previous ones, obtained by intrinsic viscosity measurements. From the Flory–Huggins interaction parameters, the infinite dilution activity coefficients of the solvents have been calculated and fitted to the well-known NRTL model. These NRTL binary interaction parameters have a great importance in modelling the separation steps in the process of obtaining the rubber.

Código híbrido avanzado de motores de plasma de efecto Hall

La propulsión eléctrica constituye hoy una tecnología muy competitiva y de gran proyección de futuro. Dentro de los diversos motores de plasma existentes, el motor de efecto Hall ha adquirido una gran madurez y constituye un medio de propulsión idóneo para un rango amplio de misiones. En la presente Tesis se estudian los motores Hall con geometría convencional y paredes dieléctricas. La compleja interacción entre los múltiples fenómenos físicos presentes hace que sea difícil la simulación del plasma en estos motores. Los modelos híbridos son los que representan un mejor compromiso entre precisión y tiempo de cálculo. Se basan en utilizar un modelo fluido para los electrones y algoritmos de dinámica de partículas PIC (Particle-In- Cell) para los iones y los neutros. Permiten hacer uso de la hipótesis de cuasineutralidad del plasma, a cambio de resolver separadamente las capas límite (o vainas) que se forman en torno a las paredes de la cámara. Partiendo de un código híbrido existente, llamado HPHall-2, el objetivo de la Tesis doctoral ha sido el desarrollo de un código híbrido avanzado que mejorara la simulación de la descarga de plasma en un motor de efecto Hall. Las actualizaciones y mejoras realizadas en las diferentes partes que componen el código comprenden tanto aspectos teóricos como numéricos. Fruto de la extensa revisión de la algoritmia del código HPHall-2 se han conseguido reducir los errores de precisión un orden de magnitud, y se ha incrementado notablemente su consistencia y robustez, permitiendo la simulación del motor en un amplio rango de condiciones. Algunos aspectos relevantes a destacar en el subcódigo de partículas son: la implementación de un nuevo algoritmo de pesado que permite determinar de forma más precisa el flujo de las magnitudes del plasma; la implementación de un nuevo algoritmo de control de población, que permite tener suficiente número de partículas cerca de las paredes de la cámara, donde los gradientes son mayores y las condiciones de cálculo son más críticas; las mejoras en los balances de masa y energía; y un mejor cálculo del campo eléctrico en una malla no uniforme. Merece especial atención el cumplimiento de la condición de Bohm en el borde de vaina, que en los códigos híbridos representa una condición de contorno necesaria para obtener una solución consistente con el modelo de interacción plasma-pared, y que en HPHall-2 aún no se había resuelto satisfactoriamente. En esta Tesis se ha implementado el criterio cinético de Bohm para una población de iones con diferentes cargas eléctricas y una gran dispersión de velocidades. En el código, el cumplimiento de la condición cinética de Bohm se consigue por medio de un algoritmo que introduce una fina capa de aceleración nocolisional adyacente a la vaina y mide adecuadamente el flujo de partículas en el espacio y en el tiempo. Las mejoras realizadas en el subcódigo de electrones incrementan la capacidad de simulación del código, especialmente en la región aguas abajo del motor, donde se simula la neutralización del chorro del plasma por medio de un modelo de cátodo volumétrico. Sin abordar el estudio detallado de la turbulencia del plasma, se implementan modelos sencillos de ajuste de la difusión anómala de Bohm, que permiten reproducir los valores experimentales del potencial y la temperatura del plasma, así como la corriente de descarga del motor. En cuanto a los aspectos teóricos, se hace especial énfasis en la interacción plasma-pared y en la dinámica de los electrones secundarios libres en el interior del plasma, cuestiones que representan hoy en día problemas abiertos en la simulación de los motores Hall. Los nuevos modelos desarrollados buscan una imagen más fiel a la realidad. Así, se implementa el modelo de vaina de termalización parcial, que considera una función de distribución no-Maxwelliana para los electrones primarios y contabiliza unas pérdidas energéticas más cercanas a la realidad. Respecto a los electrones secundarios, se realiza un estudio cinético simplificado para evaluar su grado de confinamiento en el plasma, y mediante un modelo fluido en el límite no-colisional, se determinan las densidades y energías de los electrones secundarios libres, así como su posible efecto en la ionización. El resultado obtenido muestra que los electrones secundarios se pierden en las paredes rápidamente, por lo que su efecto en el plasma es despreciable, no así en las vainas, donde determinan el salto de potencial. Por último, el trabajo teórico y de simulación numérica se complementa con el trabajo experimental realizado en el Pnnceton Plasma Physics Laboratory, en el que se analiza el interesante transitorio inicial que experimenta el motor en el proceso de arranque. Del estudio se extrae que la presencia de gases residuales adheridos a las paredes juegan un papel relevante, y se recomienda, en general, la purga completa del motor antes del modo normal de operación. El resultado final de la investigación muestra que el código híbrido desarrollado representa una buena herramienta de simulación de un motor Hall. Reproduce adecuadamente la física del motor, proporcionando resultados similares a los experimentales, y demuestra ser un buen laboratorio numérico para estudiar el plasma en el interior del motor. Abstract Electric propulsion is today a very competitive technology and has a great projection into the future. Among the various existing plasma thrusters, the Hall effect thruster has acquired a considerable maturity and constitutes an ideal means of propulsion for a wide range of missions. In the present Thesis only Hall thrusters with conventional geometry and dielectric walls are studied. The complex interaction between multiple physical phenomena makes difficult the plasma simulation in these engines. Hybrid models are those representing a better compromise between precision and computational cost. They use a fluid model for electrons and Particle-In-Cell (PIC) algorithms for ions and neutrals. The hypothesis of plasma quasineutrality is invoked, which requires to solve separately the sheaths formed around the chamber walls. On the basis of an existing hybrid code, called HPHall-2, the aim of this doctoral Thesis is to develop an advanced hybrid code that better simulates the plasma discharge in a Hall effect thruster. Updates and improvements of the code include both theoretical and numerical issues. The extensive revision of the algorithms has succeeded in reducing the accuracy errors in one order of magnitude, and the consistency and robustness of the code have been notably increased, allowing the simulation of the thruster in a wide range of conditions. The most relevant achievements related to the particle subcode are: the implementation of a new weighing algorithm that determines more accurately the plasma flux magnitudes; the implementation of a new algorithm to control the particle population, assuring enough number of particles near the chamber walls, where there are strong gradients and the conditions to perform good computations are more critical; improvements in the mass and energy balances; and a new algorithm to compute the electric field in a non-uniform mesh. It deserves special attention the fulfilment of the Bohm condition at the edge of the sheath, which represents a boundary condition necessary to match consistently the hybrid code solution with the plasma-wall interaction, and remained as a question unsatisfactory solved in the HPHall-2 code. In this Thesis, the kinetic Bohm criterion has been implemented for an ion particle population with different electric charges and a large dispersion in their velocities. In the code, the fulfilment of the kinetic Bohm condition is accomplished by an algorithm that introduces a thin non-collisional layer next to the sheaths, producing the ion acceleration, and measures properly the flux of particles in time and space. The improvements made in the electron subcode increase the code simulation capabilities, specially in the region downstream of the thruster, where the neutralization of the plasma jet is simulated using a volumetric cathode model. Without addressing the detailed study of the plasma turbulence, simple models for a parametric adjustment of the anomalous Bohm difussion are implemented in the code. They allow to reproduce the experimental values of the plasma potential and the electron temperature, as well as the discharge current of the thruster. Regarding the theoretical issues, special emphasis has been made in the plasma-wall interaction of the thruster and in the dynamics of free secondary electrons within the plasma, questions that still remain unsolved in the simulation of Hall thrusters. The new developed models look for results closer to reality, such as the partial thermalization sheath model, that assumes a non-Maxwellian distribution functions for primary electrons, and better computes the energy losses at the walls. The evaluation of secondary electrons confinement within the chamber is addressed by a simplified kinetic study; and using a collisionless fluid model, the densities and energies of free secondary electrons are computed, as well as their effect on the plasma ionization. Simulations show that secondary electrons are quickly lost at walls, with a negligible effect in the bulk of the plasma, but they determine the potential fall at sheaths. Finally, numerical simulation and theoretical work is complemented by the experimental work carried out at the Princeton Plasma Physics Laboratory, devoted to analyze the interesting transitional regime experienced by the thruster in the startup process. It is concluded that the gas impurities adhered to the thruster walls play a relevant role in the transitional regime and, as a general recomendation, a complete purge of the thruster before starting its normal mode of operation it is suggested. The final result of the research conducted in this Thesis shows that the developed code represents a good tool for the simulation of Hall thrusters. The code reproduces properly the physics of the thruster, with results similar to the experimental ones, and represents a good numerical laboratory to study the plasma inside the thruster.

Ferromanganese nodules and micro-hardgrounds associated with the Cadiz Contourite Channel (NE Atlantic): Palaeoenvironmental records of fluid venting and bottom currents

Ferromanganese nodule fields and hardgrounds have recently been discovered in the Cadiz Contourite Channel in the Gulf of Cadiz (850–1000 m). This channel is part of a large contourite depositional system generated by the Mediterranean Outflow Water. Ferromanganese deposits linked to contourites are interesting tools for palaeoenviromental studies and show an increasing economic interest as potential mineral resources for base and strategic metals. We present a complete characterisation of these deposits based on submarine photographs and geophysical, petrographic, mineralogical and geochemical data. The genesis and growth of ferromanganese deposits, strongly enriched in Fe vs. Mn (av. 39% vs. 6%) in this contourite depositional system result from the combination of hydrogenetic and diagenetic processes. The interaction of the Mediterranean Outflow Water with the continental margin has led to the formation of Late Pleistocene–Holocene ferromanganese mineral deposits, in parallel to the evolution of the contourite depositional system triggered by climatic and tectonic events. The diagenetic growth was fuelled by the anaerobic oxidation of thermogenic hydrocarbons (δ13CPDB=−20 to −37‰) and organic matter within the channel floor sediments, promoting the formation of Fe–Mn carbonate nodules. High 87Sr/86Sr isotopic values (up to 0.70993±0.00025) observed in the inner parts of nodules are related to the influence of radiogenic fluids fuelled by deep-seated fluid venting across the fault systems in the diapirs below the Cadiz Contourite Channel. Erosive action of the Mediterranean Outflow Water undercurrent could have exhumed the Fe–Mn carbonate nodules, especially in the glacial periods, when the lower core of the undercurrent was more active in the study area. The growth rate determined by 230Thexcess/232Th was 113±11 mm/Ma, supporting the hypothesis that the growth of the nodules records palaeoenvironmental changes during the last 70 ka. Ca-rich layers in the nodules could point to the interaction between the Mediterranean Outflow Water and the North Atlantic Deep Water during the Heinrich events. Siderite–rhodochrosite nodules exposed to the oxidising seabottom waters were replaced by Fe–Mn oxyhydroxides. Slow hydrogenetic growth of goethite from the seawaters is observed in the outermost parts of the exhumed nodules and hardgrounds, which show imprints of the Mediterranean Outflow Water with low 87Sr/86Sr isotopic values (down to 0.70693±0.00081). We propose a new genetic and evolutionary model for ferromanganese oxide nodules derived from ferromanganese carbonate nodules formed on continental margins above the carbonate compensation depth and dominated by hydrocarbon seepage structures and strong erosive action of bottom currents. We also compare and discuss the generation of ferromanganese deposits in the Cadiz Contourite Channel with that in other locations and suggest that our model can be applied to ferromanganiferous deposits in other contouritic systems affected by fluid venting.

Ultra intense laser/plasma interaction at normal incidence: Relativistic mirrors effects, high harmonics generation and absorption.

An analytical study of the relativistic interaction of a linearly-polarized laser-field of w frequency with highly overdense plasma is presented. Very intense high harmonics are generated produced by relativistic mirrors effects due to the relativistic electron plasma oscillation. Also, in agreement with 1D Particle-In-Cell Simulations (PICS), the model self-consistently explains the transition between the sheath inverse bremsstrahlung (SIB) absorption regime and the J×B heating (responsible for the 2w electron bunches), as well as the mean electron energy.

Coronal fluid-dynamics in laser fusion

The fluid-dynamics of the corona ejected by laser-fusion targets in the direct-drive approach (thermal radiation and atomic physics unimportant) is discussed. A two-fluid model involves inverse bremsstrahlung absorption, refraction, different ion and electron temperatures with energy exchange, different ion and electron velocities and magnetic field generation, and their effect on ion-electron friction and heat flux. Four dimensionless parameters determine coronal regimes for one-dimensional flows under uniform irradiation. One additional parameter is involved in two-dimensional problems,including the stability of one-dimensional flows, and the smoothing of nonuniform driving.

E-beam propagation and interaction with atmosphere

This study purports to investigate whether a conductive tether left uninsulated and electrically floating in LEO could serve as an effective e-beam source to produce artificial auroras. An electrically floating tether comes out biased highly negative over most of its length. Ambient ions impacting it with KeV energies liberate secondary electrons, which are locally accelerated through the 2D tether voltage-bias, race down magnetic lines, and result in peak auroral emissions at about 120-160 km altitude. Since no current flows at either tether end, a bare-tether e-beam is fully free of spacecraft charging problems. Beam propagation and beam-atmosphere interactions need be modelled in a simple but quantitative way so as to allow a satisfactory discussion of observational options and their feasibility. The evolution in the energy spectrum of secondary electrons, their pitch distribution, and beam broadening due to collisions with neutrals, which would result in a broader but weaker tether footprint in the E-layer, need be modelled. Relations between particle/energy flux values, and ionization and accompanying emission rates, are considered.