Thermal evolution and physics of melt extraction on the ureilite parent body

We develop a physical model of the thermal history of the ureilite parent body (UPB) that numerically tracks the history of its heating, hydration, dehydration, partial melting and smelting as a function of its formation time and the initial values of its composition, formation temperature and water ice content. Petrologic and chemical data from the main group (non-polymict) ureilite meteorites, which sample the interior of the UPB between depths corresponding to pressures in the range 3-10 MPa, are used to constrain the model. We find that to achieve the ∼30% melting inferred for ureilites from all sampled depths, the UPB must have had a radius between ∼80 and ∼130 km and must have accreted about 0.55 Ma after CAI formation. Melting began in the body at ∼1 Ma after CAI, and the time at which 30% melting was reached varied with depth in the asteroid but was always between ∼4.5 and ∼5.8 Ma after CAI. The total rate at which melt was produced in the UPB varied from more than 100 m³ s⁻¹ in the very early stages of melting at ∼1 Ma after CAI to ∼5 m³ s⁻¹ between 2 and 3 Ma after CAI, decreasing to extremely small values as the end of melting was approached beyond ∼5 Ma. Although the initial period of high melt production occupied only a short time around 1 Ma after CAI, it corresponded to ∼half (16%) of total silicate melting, and all strictly basaltic (i.e. plagioclase-saturated) melts must have been produced during this period.A very efficient melt transport network, consisting of a hierarchy of veins and larger pathways (dikes), developed quickly at the start of melting, ensuring rapid (timescales of months) transport of any single parcel of melt to shallow levels, thus ensuring that chemical interaction between melts and the rocks through which they subsequently passed was negligible. Volatile (mainly carbon monoxide) production due to smelting began at the start of silicate melting in the shallowest parts of the UPB and at later times at greater depths. Except at the very start and very end of melting, the volatile content of the melts produced was always high - generally between 15 and 35 mass % - and most of the melt produced was erupted at the surface of the UPB with speeds well in excess of the escape velocity and was lost into space. However, we show that 30% melting at the 3 MPa pressure level was only possible if ∼15% of the total melt produced in the asteroid was retained as a small number (∼5) of very extensive, sill-like intrusions centered at a depth of ∼7 km below the surface, near the base of the ∼8 km thick outer crust of the asteroid that was maintained at temperatures below the basalt solidus by conductive heat loss to the surface. The horizontal extents of these sills occupied about 75% of the surface area of the UPB, and the sills acted as buffers between the steady supply of melt from depth and the intermittent explosive eruption of the melt into space. We infer that samples from these intrusions are preserved as the rare feldspathic (loosely basaltic) clasts in polymict ureilites, and show that the cooling histories of the sills are consistent with these clasts reaching isotopic closure at ∼5 Ma after CAI, as given by ²⁶Al-²⁶Mg, ⁵³Mn-⁵³Cr and Pb-Pb age dates.     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.     

Flexural isostatic response of the Alps to increased Quaternary erosion recorded by foreland basin remnants,SE France

We test the hypothesis that flexural isostatic compensation of the mass removed by enhanced Quaternary erosion is responsible for uplift of the Western European Alps and their forelands. We use two well‐preserved and well‐dated (1.8 Ma) abandonment surfaces of foreland basin remnants in SE France (the Chambaran and Valensole plateaux) as passive benchmarks for tilting of the foreland. Estimating their initial slope from morphometric scaling relationships, we determine bulk post‐depositional tilting of 0.5–0.8% for these surfaces. The calculated isostatic response of the Alpine lithosphere to erosional unloading, using the method recently proposed by Champagnac et al. [Geology 35 (2007) 195–198] , yields a predicted tilting of 0.3–0.4% in the considered areas, explaining approximately half of the determined post‐depositional tilting. Such long‐term deformation being insensitive to cyclic loading/unloading because of glaciations, we suspect the other half to be related to as yet undetermined long‐wavelength and long‐lived tectonic process(es).     

The Dynamics of Covering Behaviour in Dominant Echinoid Populations from American and European West Coasts

Abstract.  The so-called 'covering behaviour' of echinoids is known to vary with habitat according to the types of covering items available, but the full extent of the role played by habitat characteristics in this behaviour is not known. Moreover, whether or not different species inhabiting similar environmental conditions and habitats also show similar patterns of this particular behaviour has yet to be investigated. In this study, two prominent west coast echinoid species, Paracentrotus lividus and Strongylocentrotus purpuratus , were examined with respect to covering behaviour and spatial organisation. Populations of both species are found in exposed intertidal zones, occupying boreholes within tidepools. Important similarities were found in the spatial organisation of both species' populations, with lowest densities in upper shore pools. However, the size distributions of the two populations differed significantly; P. lividus within upper shore pools were significantly larger than those in mid or lower shore pools, whereas S. purpuratus on the lower shore were significantly larger than those at other shore heights. We found significant between-species variability in covering item use, although the number of covering items available at both sites was not significantly different. For example, greater densities of covering items were used by P. lividus than by S. purpuratus . We also found that the percent of echinoid surface area covered varied significantly in space for both species. These data emphasise the complex and multifaceted nature of covering behaviour, but most importantly, they also strongly suggest that habitat characteristics can only partly explain the extent of this behavioural variability amongst echinoid species.     

Mercury speciation in surface and deep waters of the Mediterranean Sea

Mercury speciation and its distribution in surface and deep waters of the Mediterranean Sea were studied during two oceanographic cruises on board the Italian research vessel URANIA in summer 2003 and spring 2004 as part of the Med Oceaneor and MERCYMS projects. The study included deep water profiles of dissolved gaseous Hg (DGM), reactive Hg (RHg), total Hg (THg), monomethyl Hg (MeHg) and dimethyl Hg (DMeHg) in open ocean waters. Average concentrations of measured Hg species were characterized by seasonal and spatial variations. Overall average THg concentrations ranged between 0.41 and 2.65 pM (1.32 ± 0.48 pM) and were comparable to those obtained in previous studies of the Mediterranean Sea. A significant fraction of Hg was present as “reactive” Hg (average 0.33 ± 0.32 pM). Dissolved gaseous Hg (DGM), which consists mainly of Hg⁰, represents a considerable proportion of THg (average 20%, 0.23 ± 0.11 pM). The portion of DGM typically increased towards the bottom, especially in areas with strong tectonic activity (Alboran Sea, Strait of Sicily, Tyrrhenian Sea), indicating its geotectonic origin. No dimethyl Hg was found in surface waters down to the depth of 40 m. Below this depth, its average concentration was 2.67 ± 2.9 fM. Dissolved fractions of total Hg and MeHg were measured in filtered water samples and were 0.68 ± 0.43 pM and 0.29 ± 0.17 pM for THg and MeHg respectively. The fraction of Hg as MeHg was in average 43%, which is relatively high compared to other ocean environments. The concentrations reported in this study are among the lowest found in marine environments and the quality of analytical methods are of key importance. Speciation of Hg in sea water is of crucial importance as THg concentrations alone do not give adequate data for understanding Hg sources and cycling in marine environments. For example, photoinduced transformations are important for the presence of reactive and elemental mercury in the surface layers, biologically mediated reactions are important for the production/degradation of MeHg and DGM in the photic zones of the water column, and the data for DGM in deep sea indicate the natural sources of Hg in geotectonicaly active areas of the Mediterranean Sea.     

Vertical habitat partitioning by large calanoid copepods in the oceanic subarctic Pacific during Spring

The copepods Neocalanus plumchrus, N. flemingeri, N. cristatus, and Eucalanus bungii dominate the net zooplankton throughout the subarctic Pacific Ocean. All four species have an extensive seasonal ontogenetic vertical migration, completing most or all of their feeding and somatic growth in spring and early summer. We used stratified tows with MOCNESS and BIONESS instrumented net systems to resolve their upper ocean vertical distributions in May and June of 1984, 1987 and 1988. In each year the feeding copepodite stages of all four species were concentrated above the permanent halocline (roughly from 0 to 150m). However, the four species showed strong vertical species zonation and segregation within this layer. We consistently found a near-surface pair (N. plumchrus and N. flemingeri) and a subsurface pair (N. cristatus and E. bungii). The boundary between these groups shifts vertically, but was sharply defined and was very often coincident with a weak and transient thermocline marking the base of the layer actively mixed by surface wind and wave energy. Diel vertical migration was very limited during our sampling periods.The data suggest that the vertical distribution patterns of the copepods could be set by responses to the local intensity of turbulent mixing in the watercolumn. N. plumchrus and N. flemingeri occupied a stratum characterized by strong turbulence. N. cristatus and E. bungii occupied a stratum that was a local minimum in turbulence profiles. The depth of the boundary between the species pairs was deeper when winds and surface energy inputs were strong. The vertical partition pattern may also be determined by a difference in feeding strategy between the species pairs. N. plumchrus and N. flemingeri may feed on the enhanced protozoan population of the mixed layer, while N. cristatus and E. bungii feed on particle aggregates settling from above.