Structures and Concentrations of Surfactants in Gut Fluid of the Marine Polychaete Arenicola Marina

Marine invertebrate deposit feeders secrete surfactants into their gut fluid in concentrations sufficient to induce micelle formation, enhancing solubilization of sedimentary lipids. We isolated and identified 3 related surfactant molecules from the deposit-feeding polychaete lugworm Arenicola marina. Surfactants were isolated and separated by a combination of solvent extraction and thin-layer and gas chromatography. Identification was performed using mass and infrared spectrometry, coupled to various derivatization and hydrolysis reactions. A. marina produces a mixture of related yet distinct anionic surfactants composed of branched, C9, saturated and unsaturated fatty acids that are amide linked to leucine or glycine residues, showing some similarity to crustacean surfactants. The critical micelle concentration of the mixture of these surfactants in gut fluid was about 2 mM, and total concentrations ranged from 5.5 to 19.5 mM. The hydrophilic amide linkage helps to explain previous observations that gut surfactants do not adsorb onto sediment transiting the gut.

Characterization and differentiation of body fluids, putrefaction fluid, and blood using Hounsfield unit in postmortem CT

The purpose of the present study was to evaluate the ranges of Hounsfield unit (HU) found in body fluids, putrefaction fluids, and blood on postmortem CT and how these ranges are affected by postmortem interval, temperatures, and CT beam energy. Body fluids, putrefaction fluids, and blood from a total of 53 corpses were analyzed to determine the ranges of HU values from postmortem CT images that were taken prior to autopsy. The fluids measured in CT images were obtained at autopsy and examined in terms of macroscopic and microscopic appearances. Body fluids and blood were also collected in plastic bottles, which were subjected to CT scans at different beam energies (80-130 kV) and at various fluid temperatures (4 to 40 °C). At a postmortem interval of 1 to 4 days, the ranges of HU values of the serous fluids (13-38 HU) and the nonsedimented blood (40-88 HU) did not overlap. In the sedimented blood, the upper serum layer exhibited HU value ranges that overlapped with those of the serous fluids. The putrefaction fluids exhibited a range of HU values between 80 and -130 HU. Elevated HU values were observed in fluids with accretive cell impurities. HU values decreased slightly with increasing temperature and CT beam energy. We concluded that serous fluids and blood in fresh corpses can be characterized and differentiated from each other based on HU value ranges. In contrast, body fluids in decomposed corpses cannot be differentiated by their HU value ranges. Different beam energies and corpse temperatures had only minor influences on HU value ranges and therefore should not be obstacles to the differentiation and characterization of body fluids and blood.

Modeling episodic fluid-release events in the ductile carbonates of the Glarus thrust

The exposed Glarus thrust displays midcrustal deformation with tens of kilometers of displacement on an ultrathin layer, the principal slip zone (PSZ). Geological observations indicate that this structure resulted from repeated stick-slip events in the presence of highly overpressured fluids. Here we show that the major characteristics of the Glarus thrust movement (localization, periodicity, and evidence of pressurized fluids) can be reconciled by the coupling of two processes, namely, shear heating and fluid release by carbonate decomposition. During this coupling, slow ductile creep deformation raises the temperature through shear heating and ultimately activates the chemical decomposition of carbonates. The subsequent release of highly overpressurized fluids forms and lubricates the PSZ, allowing a ductile fault to move tens of kilometers on millimeter-thick bands in episodic stick-slip events. This model identifies carbonate decomposition as a key process for motion on the Glarus thrust and explains the source of overpressured fluids accessing the PSZ.

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.

Chemical and isotope composition of poor fluid from gas-hydrate samples of ODP Site 146-892

Although the presence of extensive gas hydrate on the Cascadia margin, offshore from the western U.S. and Canada, has been inferred from marine seismic records and pore water chemistry, solid gas hydrate has only been found at one location. At Ocean Drilling Program (ODP) Site 892, offshore from central Oregon, gas hydrate was recovered close to the sediment-water interface at 2-19 m below the seafloor (mbsf) at 670 m water depth. The gas hydrate occurs as elongated platy crystals or crystal aggregates, mostly disseminated irregularly, with higher concentrations occurring in discrete zones, thin layers, and/or veinlets parallel or oblique to the bedding. A 2- to 3-cm thick massive gas hydrate layer, parallel to bedding, was recovered at ~17 mbsf. Gas from a sample of this layer was composed of both CH4 and H2S. This sample is the first mixed-gas hydrate of CH4-H2S documented in ODP; it also contains ethane and minor amounts of CO2. Measured temperatures of the recovered core ranged from 2 to -1.8°C and are 6 to 8 degrees lower than in-situ temperatures. These temperature anomalies were caused by the partial dissociation of the CH4-H2S hydrate during recovery without a pressure core sampler. During this dissociation, toxic levels of H2S (delta34S, +27.4?) were released. The delta13C values of the CH4 in the gas hydrate, -64.5 to -67.5? (PDB), together with deltaD values of -197 to -199? (SMOW) indicate a primarily microbial source for the CH4. The delta18O value of the hydrate H2O is +2.9? (SMOW), comparable with the experimental fractionation factor for sea-ice. The unusual composition (CH4-H2S) and depth distribution (2-19 mbsf) of this gas hydrate indicate mixing between a methane-rich fluid with a pore fluid enriched in sulfide; at this site the former is advecting along an inclined fault into the active sulfate reduction zone. The facts that the CH4-H2S hydrate is primarily confined to the present day active sulfate reduction zone (2-19 mbsf), and that from here down to the BSR depth (19-68 mbsf) the gas hydrate inferred to exist is a >=99% CH4 hydrate, suggest that the mixing of CH4 and H2S is a geologically young process. Because the existence of a mixed CH4-H2S hydrate is indicative of moderate to intense advection of a methane-rich fluid into a near surface active sulfate reduction zone, tectonically active (faulted) margins with organic-rich sediments and moderate to high sedimentation rates are the most likely regions of occurrence. The extension of such a mixed hydrate below the sulfate reduction zone should reflect the time-span of methane advection into the sulfate reduction zone.

Linear global instability of non-orthogonal incompressible swept attachment-line boundary layer flow

Instability of the orthogonal swept attachment line boundary layer has received attention by local1, 2 and global3–5 analysis methods over several decades, owing to the significance of this model to transition to turbulence on the surface of swept wings. However, substantially less attention has been paid to the problem of laminar flow instability in the non-orthogonal swept attachment-line boundary layer; only a local analysis framework has been employed to-date.6 The present contribution addresses this issue from a linear global (BiGlobal) instability analysis point of view in the incompressible regime. Direct numerical simulations have also been performed in order to verify the analysis results and unravel the limits of validity of the Dorrepaal basic flow7 model analyzed. Cross-validated results document the effect of the angle _ on the critical conditions identified by Hall et al.1 and show linear destabilization of the flow with decreasing AoA, up to a limit at which the assumptions of the Dorrepaal model become questionable. Finally, a simple extension of the extended G¨ortler-H¨ammerlin ODE-based polynomial model proposed by Theofilis et al.4 is presented for the non-orthogonal flow. In this model, the symmetries of the three-dimensional disturbances are broken by the non-orthogonal flow conditions. Temporal and spatial one-dimensional linear eigenvalue codes were developed, obtaining consistent results with BiGlobal stability analysis and DNS. Beyond the computational advantages presented by the ODE-based model, it allows us to understand the functional dependence of the three-dimensional disturbances in the non-orthogonal case as well as their connections with the disturbances of the orthogonal stability problem.

Tranversal motion and flow structure of fully nonlinear streaks in a laminar boundary layer

Typical streak computations present in the literature correspond to linear streaks or to small amplitude nonlinear streaks computed using DNS or nonlinear PSE. We use the Reduced Navier-Stokes (RNS) equations to compute the streamwise evolution of fully non-linear streaks with high amplitude in a laminar flat plate boundary layer. The RNS formulation provides Reynolds number independent solutions that are asymptotically exact in the limit $Re \gg 1$, it requires much less computational effort than DNS, and it does not have the consistency and convergence problems of the PSE. We present various streak computations to show that the flow configuration changes substantially when the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-streamwise plane) becomes more important and strongly distorts the streamwise velocity profiles, that end up being quite different from those of the linear case. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks and compare them with available experimental results.

Linear instability of orthogonal compressible leading-edge boundary layer flow

Instability analysis of compressible orthogonal swept leading-edge boundary layer flow was performed in the context of BiGlobal linear theory. 1, 2 An algorithm was developed exploiting the sparsity characteristics of the matrix discretizing the PDE-based eigenvalue problem. This allowed use of the MUMPS sparse linear algebra package 3 to obtain a direct solution of the linear systems associated with the Arnoldi iteration. The developed algorithm was then applied to efficiently analyze the effect of compressibility on the stability of the swept leading-edge boundary layer and obtain neutral curves of this flow as a function of the Mach number in the range 0 ≤ Ma ≤ 1. The present numerical results fully confirmed the asymptotic theory results of Theofilis et al. 4 Up to the maximum Mach number value studied, it was found that an increase of this parameter reduces the critical Reynolds number and the range of the unstable spanwise wavenumbers.

Low-dimensional modeling of streaks in a wedge flow boundary layer

This paper is concerned with the low dimensional structure of optimal streaks in a wedge flow boundary layer, which have been recently shown to consist of a unique (up to a constant factor) three-dimensional streamwise evolving mode, known as the most unstable streaky mode. Optimal streaks exhibit a still unexplored/unexploited approximate self-similarity (not associated with the boundary layer self-similarity), namely the streamwise velocity re-scaled with their maximum remains almost independent of both the spanwise wavenumber and the streamwise coordinate; the remaining two velocity components instead do not satisfy this property. The approximate self-similar behavior is analyzed here and exploited to further simplify the description of optimal streaks. In particular, it is shown that streaks can be approximately described in terms of the streamwise evolution of the scalar amplitudes of just three one-dimensional modes, providing the wall normal profiles of the streamwise velocity and two combinations of the cross flow velocity components; the scalar amplitudes obey a singular system of three ordinary differential equations (involving only two degrees of freedom), which approximates well the streamwise evolution of the general streaks.

Nonlinear description of transversal motion in a laminar boundary layer with streaks

The nonlinear streamwise growth of a spanwise periodic array of steady streaks in a flat plate boundary layer is numerically computed using the well known Reduced Navier- Stokes formulation. It is found that the flow configuration changes substantially when the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-spanwise plane), which is normally not considered, becomes non-negligible in the nonlinear regime, and it strongly distorts the streamwise velocity profiles, which end up being quite different from those predicted by the linear theory. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks, and compare them with available experimental results.

Low-Dimensional Modelling for optimal streaks in the blasius boundary layer

This paper is concerned with the low dimensional structure of optimal streaks in the Blasius boundary layer. Optimal streaks are well known to exhibit an approximate self-similarity, namely the streamwise velocity re-scaled with their maximum remains almost independent of both the spanwise wavenumber and the streamwise coordinate. However, the reason of this self-similar behavior is still unexplained as well as unexploited. After revisiting the structure of the streaks near the leading edge singularity, two additional approximately self-similar relations involving the velocity components and their wall normal derivatives are identified. Based on these properties, we derive a low dimensional model with two degrees of freedom. The comparison with the results obtained from the linearized boundary layer equations shows that this model is consistent and provide good approximations.

Laminar-turbulent transition induced by a discrete roughness element in a supersonic boundary layer

The linear instability and breakdown to turbulence induced by an isolated roughness element in a boundary layer at Mach 2:5, over an isothermal flat plate with laminar adiabatic wall temperature, have been analysed by means of direct numerical simulations, aided by spatial BiGlobal and three-dimensional parabolized (PSE-3D) stability analyses. It is important to understand transition in this flow regime since the process can be slower than in incompressible flow and is crucial to prediction of local heat loads on next-generation flight vehicles. The results show that the roughness element, with a height of the order of the boundary layer displacement thickness, generates a highly unstable wake, which is composed of a low-velocity streak surrounded by a three-dimensional high-shear layer and is able to sustain the rapid growth of a number of instability modes. The most unstable of these modes are associated with varicose or sinuous deformations of the low-velocity streak; they are a consequence of the instability developing in the three-dimensional shear layer as a whole (the varicose mode) or in the lateral shear layers (the sinuous mode). The most unstable wake mode is of the varicose type and grows on average 17% faster tan the most unstable sinuous mode and 30 times faster than the most unstable boundary layer mode occurring in the absence of a roughness element. Due to the high growthrates registered in the presence of the roughness element, an amplification factor of N D 9 is reached within 50 roughness heights from the roughness trailing edge. The independently performed Navier–Stokes, spatial BiGlobal and PSE-3D stability results are in excellent agreement with each other, validating the use of simplified theories for roughness-induced transition involving wake instabilities. Following the linear stages of the laminar–turbulent transition process, the roll-up of the three-dimensional shear layer leads to the formation of a wedge of turbulence, which spreads laterally at a rate similar to that observed in the case of compressible turbulent spots for the same Mach number.

Numerical issues in Lagrangian tracking and topological evolution of fluid particles in wallbounded turbulent flows

The determination of the local Lagrangian evolution of the flow topology in wall-bounded turbulence, and of the Lagrangian evolution associated with entrainment across the turbulent / non-turbulent interface into a turbulent boundary layer, require accurate tracking of a fluid particle and its local velocity gradients. This paper addresses the implementation of fluid-particle tracking in both a turbulent boundary layer direct numerical simulation and in a fully developed channel flow simulation. Determination of the sub-grid particle velocity is performed using both cubic B-spline, four-point Hermite spline and higher-order Hermite spline interpolation. Both wall-bounded flows show similar oscillations in the Lagrangian tracers of both velocity and velocity gradients, corresponding to the movement of particles across the boundaries of computational cells. While these oscillation in the particle velocity are relatively small and have negligible effect on the particle trajectories for time-steps of the order of CFL = 0.1, they appear to be the cause of significant oscillations in the evolution of the invariants of the velocity gradient tensor.

Chemistry and fluid inclusions of the sheeted dike comples of ODP Holes 137-504B and 140-504B

Fluid inclusions in variably altered diabase recovered from Ocean Drilling Program Legs 137 and 140 at Hole 504B, Costa Rica Rift, exhibit fluid salinities up to 3.7 times that of seawater values (11.7 wt% NaCl equivalent) and exhibit uncorrected homogenization temperatures of 125°C to 202°C. The liquid-dominated inclusions commonly are entrapped in zones of secondary plagioclase and may be primary in origin. Fluid salinities are similar to compositions of fluids venting on the seafloor (0.4-7.0 wt% NaCl) and overlap with those measured in metabasalt samples recovered from near the Kane Fracture Zone on the Mid-Atlantic Ridge and from the Troodos ophiolite, Cyprus. The salinity variations may reflect hydration reactions involving formation of secondary mineral assemblages under rock-dominated conditions, which modify the ionic strength of hydrothermal fluids by consuming or liberating water and chloride ion. Rare CO2-CH4-bearing inclusions, subjacent to zones where talc after olivine becomes an important secondary mineral phase (1700 mbsf), may have formed due to local interaction of seawater and olivine at low water to rock ratios. Corrected average fluid inclusion homogenization temperatures exhibit a gradient from 159°C at a depth of 1370 mbsf to 183°C at a depth of 1992 mbsf and are in apparent equilibrium with the present conductive downhole temperatures. These data indicate that fluid inclusions may be used to estimate downhole temperatures if logging data are unavailable. The compositional and thermal evolution of the diabase-hosted fluids may reflect late-stage, off-axis circulation and conductive heating of compositionally modified seawater in the sheeted dike complex at Hole 504B.