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.

Results of analyses of fluid inclusions from ODP Site 735B, Leg 118 and 176

Microthermometric and isotopic analyses of fluid inclusions in primitive olivine gabbros, oxide gabbros, and evolved granitic material recovered from Ocean Drilling Program Hole 735B at the Southwest Indian Ridge provide new insights into the evolution of C-O-H-NaCl fluids in the plutonic foundation of the oceanic crust. The variably altered and deformed plutonic rocks span a crustal section of over 1500 m and record a remarkably complex magma-hydrothermal history. Magmatic fluids within this suite followed two chemically distinct paths during cooling through the subsolidus regime: the first path included formation of CO2+CH4+H2O+C fluids with up to 43 mole% CH4; the second path produced hypersaline brines that contain up to 50% NaCl equivalent salinities. Subsequent to devolatilization, respeciation of magmatic CO2, attendant graphite precipitation, and cooling from 800°C to 500°C promoted formation of CH4-enriched fluids. These fluids are characterized by average d13C(CH4) values of -27.1+/-4.3 per mil (N=45) with associated d13C(CO2) compositions ranging from -24.9 per mil to -1.9 per mil (N=39), and average dD values of exsolved vapor of -41+/-12 per mil (N=23). In pods, veins, and lenses of highly fractionated residual material, hypersaline brines formed during condensation and by direct exsolution in the absence of a conjugate vapor phase. Entrapped CO2+CH4+H2O-rich fluids within many oxide-bearing rocks and felsic zones are significantly depleted in 13C (with d13C(CO2) values down to about -25 per mil) and contain CO2 concentrations higher than those predicted by equilibrium devolatilization models. We hypothesize that lower effective pressures in high-temperature shear zones promoted infiltration of highly fractionated melts and compositionally evolved volatiles into focused zones of deformation, significantly weakening the rock strength. In felsic-rich zones, volatile build-up may have driven hydraulic fracturing of gabbroic wall rocks resulting in the formation of magmatic breccias. Comparison of isotopic compositions of fluids in plutonic rocks from 735B, the MARK area of the Mid-Atlantic Ridge, and the Mid-Cayman Rise indicate (1) that the carbon isotope composition of the lower oceanic crust may be far more heterogeneous than previously believed and (2) that carbon-bearing species in the oceanic crust and their distribution at depth are highly variable.

Geochemistry, isotopes, trace elements and fluid inclusion geothermometry at DSDP Hole 83-504B

Stockwork-like metal sulfide mineralizations were found at 910-928 m below seafloor (BSF) in the pillow/dike transition zone of Hole 504B. This is the same interval where most physical properties of the 5.9-m.y.-old crust of the Costa Rica Rift change from those characteristic of Layer 2B to those of Layer 2C. The pillow lavas, breccias, and veins of the stockwork-like zone were studied by transmitted and reflected light microscopy, X-ray diffraction, and electron microprobe analysis. Bulk rock oxygen isotopic analyses as well as isolated mineral oxygen and sulfur isotopic analyses and fluid inclusion measurements were carried out. A complex alteration history was reconstructed that includes three generations of fissures, each followed by precipitation of characteristic hydrothermal mineral parageneses: (1) Minor and local deposition of quartz occurred on fissure walls; adjacent wall rocks were silicified, followed by formation of chlorite and minor pyrite I in the veins, whereas albite, sphene, chlorite and chlorite-expandable clay mixtures, actinolite, and pyrite replaced igneous phases in the host rocks. The hydrothermal fluids responsible for this first stage were probably partially reacted seawater, and their temperatures were at least 200-250° C. (2) Fissures filled during the first stage were reopened and new cracks formed. They were filled with quartz, minor chlorite and chlorite-expandable clay mixtures, traces of epidote, common pyrite, sphalerite, chalcopyrite, and minor galena. During the second stage, hydrothermal fluids were relatively evolved metal- and Si-rich solutions whose temperatures ranged from 230 to 340° C. The fluctuating chemical composition and temperature of the solutions produced a complex depositional sequence of sulfides in the veins: chalcopyrite I, ± Fe-rich sphalerite, chalcopyrite II ("disease"), Fe-poor sphalerite, chalcopyrite III, galena, and pyrite II. (3) During the last stage, zeolites and Mg-poor calcite filled up the remaining spaces and newly formed cracks and replaced the host rock plagioclase. Analcite and stilbite were first to form in veins, possibly at temperatures below 200°C; analcite and earlier quartz were replaced by laumontite at 250°C, whereas calcite formation temperature ranged from 135 to 220°C. The last stage hydrothermal fluids were depleted in Mg and enriched in Ca and 18O compared to seawater and contained a mantle carbon component. This complex alteration history paralleling a complex mineral paragenesis can be interpreted as the result of a relatively long-term evolution of a hydrothermal system with superimposed shorter term fluctuations in solution temperature and composition. Hydrothermal activity probably began close to the axis of the Costa Rica Rift with the overall cooling of the system and multiple fracturing stages due to movement of the crust away from the axis and/or cooling of a magmatic heat source.

Geochemistry of the sheeted dike complex at Pito Deep

Alteration of sheeted dikes exposed along submarine escarpments at the Pito Deep Rift (NE edge of the Easter microplate) provides constraints on the crustal component of axial hydrothermal systems at fast spreading mid-ocean ridges. Samples from vertical transects through the upper crust constrain the temporal and spatial scales of hydrothermal fluid flow and fluid-rock reaction. The dikes are relatively fresh (average extent of alteration is 27%), with the extent of alteration ranging from 0 to >80%. Alteration is heterogeneous on scales of tens to hundreds of meters and displays few systematic spatial trends. Background alteration is amphibole-dominated, with chlorite-rich dikes sporadically distributed throughout the dike complex, indicating that peak temperatures ranged from <300°C to >450°C and did not vary systematically with depth. Dikes locally show substantial metal mobility, with Zn and Cu depletion and Mn enrichment. Amphibole and chlorite fill fractures throughout the dike complex, whereas quartz-filled fractures and faults are only locally present. Regional variability in alteration characteristics is found on a scale of <1-2 km, illustrating the diversity of fluid-rock interaction that can be expected in fast spreading crust. We propose that much of the alteration in sheeted dike complexes develops within broad, hot upwelling zones, as the inferred conditions of alteration cannot be achieved in downwelling zones, particularly in the shallow dikes. Migration of circulating cells along rides axes and local evolution of fluid compositions produce sections of the upper crust with a distinctive character of alteration, on a scale of <1-2 km and <5-20 ka.

Rare earth elements of ROV sample GeoB12338-2 from the Makran accretionary complex

Authigenic carbonates forming at an active methane-seep on the Makran accretionary prism mainly consist of aragonite in the form of microcrystalline, cryptocrystalline, and botryoidal phases. The d13Ccarbonate values are very negative (-49.0 to -44.0 per mill V-PDB), agreeing with microbial methane as dominant carbon source. The d18Ocarbonate values are exclusively positive (+ 3.0 to + 4.5 per mill V-PDB) and indicate precipitation in equilibrium with seawater at bottom water temperatures. The content of rare earth elements and yttrium (REE + Y) determined by laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) and solution ICP-MS varies for each aragonite variety, with early microcrystalline aragonite yielding the highest, cryptocrystalline aragonite intermediate, and later botryoidal aragonite the lowest REE + Y concentrations. Shale-normalised REE + Y patterns of different types of authigenic carbonate reflect distinct pore fluid compositions during precipitation: Microcrystalline aragonite shows high contents of middle rare earth elements (MREE), reflecting REE patterns ascribed to anoxic pore water. Cryptocrystalline aragonite exhibits a seawater-like REE + Y pattern at elevated total REE + Y concentrations, indicating higher concentrations of REEs in pore waters, which were influenced by seawater. Botryoidal aragonite is characterised by seawater-like REE + Y patterns at initial growth stages followed by an increase of light rare earth elements (LREE) with advancing crystal growth, reflecting changing pore fluid composition during precipitation of this cement. Conventional sample preparation involving micro-drilling of carbonate phases and subsequent solution ICP-MS does not allow to recognise such subtle changes in the REE + Y composition of individual carbonate phases. To be able to reconstruct the evolution of pore water composition during early diagenesis, an analytical approach is required that allows to track the changing elemental composition in a paragenetic sequence as well as in individual phases. High-resolution analysis of seep carbonates from the Makran accretionary prism by LA-ICP-MS reveals that pore fluid composition not only evolved in the course of the formation of different phases, but also changed during the precipitation of individual phases.

Chemistry of secondary minerals from the deep sheeted dike complex of ODP Holes 137-504B and 140-504B

Dolerites sampled from the lower sheeted dikes from Hole 504B during Ocean Drilling Program Legs 137 and 140, between 1562.4 and 2000.4 mbsf, were examined to document the mineralogy, petrography, and mineral parageneses associated with secondary alteration, to constrain the thermal history and composition of hydrothermal fluids. The main methods used were mineral chemical analyses by electron microprobe, X-ray diffraction, and cathodoluminescence microscopy. Temperatures of alteration were estimated on the basis of single and/or coexisting mineral chemistry. Permeability is important in controlling the type and extent of alteration in the studied dike section. At the meter-scale, intervals of weakly altered dolerites containing fresh olivine are interpreted as having experienced restricted exposure to hydrothermal fluids. At the centimeter- or millimeter-scale, alteration patches and extensively altered halos adjacent to veins reflect the permeability related to intergranular primary porosity and cracks. Most of the sheeted dike alteration in this case resulted from non-focused, pervasive fluid-rock interaction. This study confirms and extends the previous model for hydrothermal alteration at Hole 504B: hydrothermal alteration at the ridge axis followed by seawater recharge and off-axis alteration. The major new discoveries, all related to higher temperatures of alteration, are: (1) the presence of hydrothermal plagioclase (An80-95), (2) the presence of deuteric and/or hydrothermal diopside, and (3) the general increasing proportion of amphiboles, and particularly magnesio-hornblende with depth. We propose that the dolerites at Hole 504B were altered in five stages. Stage 1 occurred at high temperatures (less than 500° to 700°C) and involved late-magmatic formation of Na- and Ti-rich diopside, the hydrothermal formation of Na, Ti-poor diopside and the hydrothermal formation of an assemblage of An-rich plagioclase + hornblende. Stage 2 occurred at lower temperatures (250°-320°C) and is characterized by the appearance of actinolite, chlorite, chlorite-smectite, and/or talc (in low permeability zones) and albite. During Stage 3, quartz and epidote precipitated from evolved hydrothermal fluids at temperatures between 310° and 320°C. Anhydrite appeared during Stage 4 and likely precipitated directly from heated seawater. Stage 5 occurred off-axis at low temperatures (250°C) with laumontite and prehnite from evolved fluids.

New data on the ophiolitic VMS deposits of Moeche (Cabo Ortegal Complex, NW Spain)

As a result of the variscan collision, several allochtonous complexes were emplaced on the Iberian margin in Devonian times, among them the Cabo Ortegal Complex comprising the Moeche ophiolitic sequence. Copper has been won from several mines (Piquitos I & II, Barqueira, Maruxa) from disseminated ores and thin massive sulphide layers in the Moeche Unit, a strongly deformed meta-volcanic sequence comprising mainly quartz-chlorite schists and mylonites, which defines the top of the ophiolite. The ores were metamorphosed and strongly deformed under brittle conditions (for pyrite), but their textures are often apparently post-deformational, due to very common solution-transfer processes; they are composed mostly of pyrite and chalcopyrite, with minor sphalerite, pyrrhotite, etc., and with traces of native gold and PGE. The geology, mineralogy, and geochemistry of the orebodies relate closely to VMS of the Cu-Zn (Cyprus) type. Fluid inclusion studies allowed an estimation of metamorphic conditions at pressures of 2/2’5 kb and T 325/350ºC. New determinations using the chlorite geothermometer yield temperatures around 320 ºC, corresponding to pressures near 2 kb according to the isochores deduced from the fluid inclusion study, although in the Barqueira mine higher temperatures, up to 350 ºC, are found, corresponding to presssures up to 2’5 kb. Pb isotopic compositions of pyrite point to a double source of Pb, i.e. a main mantle and a subordinate crustal source. The values for 87SR/86Sr in pyrite support this interpretation, but some results suggest later mobilization in an open system, corresponding to solution-transfer. Age determinations of pyrite deduced from the Pb isotope uranogenic graph, ≈ 480 Ma, do not fit with the metamorphic ages published for the Moeche Unit, and might point to the age of Pb extraction from the mantle.

Coupling time-stepping numerical methods and standard aerodynamics codes for instability analysis of flows in complex geometries

The development of a global instability analysis code coupling a time-stepping approach, as applied to the solution of BiGlobal and TriGlobal instability analysis 1, 2 and finite-volume-based spatial discretization, as used in standard aerodynamics codes is presented. The key advantage of the time-stepping method over matrix-formulation approaches is that the former provides a solution to the computer-storage issues associated with the latter methodology. To-date both approaches are successfully in use to analyze instability in complex geometries, although their relative advantages have never been quantified. The ultimate goal of the present work is to address this issue in the context of spatial discretization schemes typically used in industry. The time-stepping approach of Chiba 3 has been implemented in conjunction with two direct numerical simulation algorithms, one based on the typically-used in this context high-order method and another based on low-order methods representative of those in common use in industry. The two codes have been validated with solutions of the BiGlobal EVP and it has been showed that small errors in the base flow do not have affect significantly the results. As a result, a three-dimensional compressible unsteady second-order code for global linear stability has been successfully developed based on finite-volume spatial discretization and time-stepping method with the ability to study complex geometries by means of unstructured and hybrid meshes

Fluid physics under reduced gravity : an overview

Problems related to several fluid physics experiments the case of two experiemnts to be performed under reduced gravity, onboard Space- temperatures,lab (1st Mission), are discussed. Special attention is placed on parallel and preparatory work on erth, wich could throw some light on the opportunity, present interest and limitations of these experiments. The need for strong supporting fundamental research and, in particular, a more precise determination of the paremeters involved, is stressed.

Simulation of wind farms in flat and complex terrain using CFD

Use of computational fluid dynamic (CFD) methods to predict the power production from wind entire wind farms in flat and complex terrain is presented in this paper. Two full 3D Navier–Stokes solvers for incompressible flow are employed that incorporate the k–ε and k–ω turbulence models respectively. The wind turbines (W/Ts) are modelled as momentum absorbers by means of their thrust coefficient using the actuator disk approach. The WT thrust is estimated using the wind speed one diameter upstream of the rotor at hub height. An alternative method that employs an induction-factor based concept is also tested. This method features the advantage of not utilizing the wind speed at a specific distance from the rotor disk, which is a doubtful approximation when a W/T is located in the wake of another and/or the terrain is complex. To account for the underestimation of the near wake deficit, a correction is introduced to the turbulence model. The turbulence time scale is bounded using the general “realizability” constraint for the turbulent velocities. Application is made on two wind farms, a five-machine one located in flat terrain and another 43-machine one located in complex terrain. In the flat terrain case, the combination of the induction factor method along with the turbulence correction provides satisfactory results. In the complex terrain case, there are some significant discrepancies with the measurements, which are discussed. In this case, the induction factor method does not provide satisfactory results.

Modelling wake effects in large wind farms in complex terrain: the problem, the methods and the issues

Computational fluid dynamic (CFD) methods are used in this paper to predict the power production from entire wind farms in complex terrain and to shed some light into the wake flow patterns. Two full three-dimensional Navier–Stokes solvers for incompressible fluid flow, employing k − ϵ and k − ω turbulence closures, are used. The wind turbines are modeled as momentum absorbers by means of their thrust coefficient through the actuator disk approach. Alternative methods for estimating the reference wind speed in the calculation of the thrust are tested. The work presented in this paper is part of the work being undertaken within the UpWind Integrated Project that aims to develop the design tools for next generation of large wind turbines. In this part of UpWind, the performance of wind farm and wake models is being examined in complex terrain environment where there are few pre-existing relevant measurements. The focus of the work being carried out is to evaluate the performance of CFD models in large wind farm applications in complex terrain and to examine the development of the wakes in a complex terrain environment.