Seismic wave velocity and anisotropy of serpentinized peridotite in the Oman ophiolite

Shallow seismic measurements in harzburgite from the Oman ophiolite performed in a zone where the maximum horizontal anisotropy is expected (vertical foliation and horizontal lineation) point to a dominant dependence of seismic properties on fracturing.

Optical microscopy studies show that microcracks are guided by the serpentine (lizardite) penetrative network oriented subparallel to the harzburgite foliation and subperpendicular to the mineral lineation, and that serpentine (lizardite) vein filling has a maximum concentration of (001) planes parallel to the veins walls. The calculated elastic properties of the oriented alteration veins filled with serpentine in an anisotropic matrix formed by oriented crystals of olivine and orthopyroxene are compared with seismic velocities measured on hand specimens.

Laboratory ultrasonic data indicate that open microcracks are closed at 100 MPa pressure, e.g. (J. Geophys. Res. 65, (1960) 1083) and (Proc. ODP Sci. Results Leg 118, (1990) 227). Above this pressure, laboratory measurements and modeling show that P-compressional and S-shear wave velocities are mainly controlled by the mineral preferred orientation. Veins sealed with serpentine are effective in slightly lowering P and S velocities and increasing anisotropy. The penetrative lizardite network does not affect directly the geometry of seismic anisotropy, but contributes indirectly in the fact that this network controls the microcrack orientations.

Comparison between seismic measurements of peridotite and gabbro in the same conditions suggest that P- and S-waves anisotropies are a possible discriminating factor between the two lithologies in the suboceanic lithosphere.     

Phylogeographical patterns of the gastropod Siphonaria pectinata (Linnaeus, 1758)

Siphonaria pectinata is a marine gastropod distributed across the shores of the Western and Eastern Atlantic. Although previous studies have indicated that the species could be a recent introduction to Central and North America from Western Africa and Europe, other authors have questioned this statement. In fact, a recent study has shown that individuals from both sides of the Atlantic are genetically disconnected, rejecting the idea of a recent invasion. However, genetic information on this species is very scarce and based only on a single study using a limited number of specimens and sites. The data reported in the present study represent one step forward, collecting information from 119 specimens from 22 sites that cover the entire distributional range of the species. Populations from the Western Atlantic separated from eastern ones around 1.4?1.9 Mya, whereas the centre of dispersion is located on the eastern side of the Atlantic. Siphonaria pectinata should be considered native to the Western Atlantic and cannot be considered as a species with an amphi‐atlantic distribution owing to the lack of genetic connection between both sides of the Atlantic, as previous authors have pointed out. The presence of a star‐shaped phylogeny is discussed in the context of Plio‐Pleistocene glaciations.     

goGPS: open-source MATLAB software

GPS Solutions - goGPS is a positioning software application designed to process single-frequency code and phase observations for absolute or relative positioning. Published under a free and...     

Bioaccumulation and depuration of Zn and Cd in mangrove oysters (Crassostrea rhizophorae, Guilding, 1828) transplanted to and from a contaminated tropical coastal lagoon

In order to study Zn and Cd accumulation and depuration, a set of oysters, Crassostrea rhizophorae, were transplanted to a metal contaminated coastal lagoon and another one was harvested there and transplanted to a non-polluted site. C. rhizophorae oysters and Perna perna mussels native from both sites were collected in order to monitor variability of metal concentration in resident populations. After three months exposure, oysters transplanted to the polluted site accumulated fourfold Zn (307-1319 microgg(-1)) without reaching the concentration level of resident oysters (9770 microgg(-1)). Cadmium concentrations had a slight but significant decrease during the same period (1.25-0.54 microgg(-1)). Oysters transplanted to the non-polluted site, showed threefold Zn depuration (6727-2404 microgg(-1)), while Cd had no significant variation (0.90-1.45 microgg(-1)). Results showed that transplanted oysters do not reach heavy metal concentrations in indigenous populations suggesting transplanted organisms would be better used to evaluate bioavailability instead of environmental concentrations.     

Shallow water gravity waves: A note on the particle orbits

When surface gravity waves of small amplitude progress in shallow water of constant mean depth, the fluid particle orbits are observed to be oval, where the longer axis of the oval is parallel to the flat bottom, and at the bottom the orbits are straight lines. Potential flow, upon which the standard wave theory is based, predicts that the oval orbits are ellipses, but by a rather lengthy mathematical procedure that is founded on the questionable assumption of irrotationality. Using a more elementary and physical method, that does not employ the irrotational assumption, the elliptical orbits can be understood much more easily. The elementary method features a balance of two oppositely directed forces on each fluid particle: the outward centrifugal force and the inward pressure force.     

Shoaling Surface Gravity Waves Cause a Force and a Torque on the Bottom

Freely propagating surface gravity waves are observed to slow down and to stop at a beach when the bottom has a relatively gentle upward slope toward the shore and the frequency range of the waves covers the most energetic wind waves (sea and swell). Essentially no wave reflection can be seen and the measured reflected energy is very small compared to that transmitted shoreward. One consequence of this is that the flux of the wave’s linear momentum decreases in the direction of wave propagation, which is equivalent to a time rate of change of the momentum. It takes a force to cause the time rate of change of the momentum. Therefore, the bottom exerts a force on the waves in order to decrease the momentum flux. By Newton’s third law (action equals reaction) the waves then impart an equal but opposite force to the bottom. In shallow (but finite) water depths the wave force per unit bottom area is calculated, for normal angle of incidence to the beach, to be directly proportional to the square of the wave amplitude and to the bottom slope and inversely proportional to the mean depth; it is independent of the wave frequency. Constants of proportionality are: 1/4, the fluid density and the acceleration of gravity. Swell attenuation near coasts and some characteristics of sand movement in the near-shore region are not inconsistent with the algebraic structure of the wave force formula. Since the force has a depth variation which is significantly faster than that of the dimensions of the particle orbits in the vertical direction, the bottom induces a torque on the fluid particles that decreases the angular momentum flux of the waves. By an extension of Newton’s third law, the waves also exert an equal but opposite torque on the bottom. And because the bottom force on the waves exists over a horizontal distance, it does work on the waves and decreases their energy flux. Thus, theoretically, the fluxes of energy, angular and linear momentum are not conserved for shoaling surface gravity waves. Mass flux, associated with the Stokes drift, is assumed to be conserved, and the wave frequency is constant for a steady medium.     

Comet Hale-Bopp in outburst: Imaging the dynamics of icy particles with HST/NICMOS

Comet Hale-Bopp was imaged at wavelengths from 1.87 to 2.22 μm by HST/NICMOS in post-perihelion observations starting on UT 1997 August 27.95. Diffraction-limited (∼02) images were obtained at high signal-to-noise (∼1500) to probe the composition and dynamics of the inner coma and also the size and activity of the nucleus. The velocities of several unusual morphological features over a 1.7 h period, indicate that a significant outburst occurred 7.4 h prior to these images while the comet was at a heliocentric distance of 2.49 AU. Similar features are also apparent after re-analysis of pre-perihelion ground-based images. The inner coma (radius ?2500 km) is dominated by an “arc” feature, which expanded and became more diffuse with time. This feature can be modeled as the bright central portion of a “jet of outburst” from a near-equatorial region of the nucleus. Less prominent, time-variable linear and circular morphologies are also apparent. The expansion rates of both the arc feature and the circular morphologies imply a common origin and also suggest a grain size distribution with two broad maxima. In addition, several static linear features extend to the edge of the field of view (21,100 km). Radial brightness profiles are highly asymmetric and only approach a ρ⁻¹ decline at distances ?15,000 km. Images in a narrow-band filter at 2.04 μm exhibit a ∼4% absorption feature relative to nearly simultaneous images at wavelengths of 2.22, 1.90, and 1.87 μm. This absorption is attributed to H₂O ice in the coma grains. The spatial distribution and expansion velocity of the absorption at 2.04 μm indicate that these grains are associated with the outburst. The constancy of the absorption feature indicates no appreciable sublimation over 1.7 h. The unresolved nucleus has a flux density consistent with a 40±10 km diameter assuming a 4% geometric albedo.     

Monitoring GOCE gradiometer calibration parameters using accelerometer and star sensor data: methodology and first results

The Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite, launched on 17 March 2009, is designed to measure the Earth’s mean gravity field with unprecedented accuracy at spatial resolutions down to 100?km. The accurate calibration of the gravity gradiometer on-board GOCE is of utmost importance for achieving the mission goals. ESA’s baseline method for the calibration uses star sensor and accelerometer data of a dedicated calibration procedure, which is executed every 2?months. In this paper, we describe a method for monitoring the evolution of calibration parameter during that time. The method works with star sensor and accelerometer data and does not require gravity field models, which distinguishes it from other existing methods. We present time series of calibration parameters estimated from GOCE data from 1 November 2009 to 17 May 2010. The time series confirm drifts in the calibration parameters that are present in the results of other methods, including ESA’s baseline method. Although these drifts are very small, they degrade the gravity gradients, leading to the conclusion that the calibration parameters of the ESA’s baseline method need to be linearly interpolated. Further, we find a correction of ?36 × 10^?6 for one calibration parameter (in-line differential scale factor of the cross-track gradiometer arm), which improves the gravity gradient performance. The results are validated by investigating the trace of the calibrated gravity gradients and comparing calibrated gravity gradients with reference gradients computed along the GOCE orbit using the ITG-Grace-2010s gravity field model.     

A global method for coupling transport with chemistry in heterogeneous porous media

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection–diffusion PDEs coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper, a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton–Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that one be able to solve chemical equilibrium problems (and compute derivatives) without having to know the solution method. An additional advantage of the Newton–Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.