Stirring in 3-d spherical models of convection in the Earth's mantle

On a global scale basalts from mid-ocean ridges are strikingly more homogeneous than basalts from intraplate volcanism. The observed geochemical heterogeneity argues strongly for the existence of distinct reservoirs in the Earth's mantle. It is an unresolved problem of Geodynamics as to how these findings can be reconciled with large-scale convection. We review observational constraints, and investigate stirring properties of numerical models of mantle convection. Conditions in the early Earth may have supported layered convection with rapid stirring in the upper layers. Material that has been altered near the surface is transported downwards by small-scale convection. Thereby a layer of homogeneous depleted material develops above pristine mantle. As the mantle cools over Earth history, the effects leading to layering become reduced and models show the large-scale convection favoured for the Earth today. Laterally averaged, the upper mantle below the lithosphere is least affected by material that has experienced near-surface differentiation. The geochemical signature obtained during the previous episode of small-scale convection may be preserved there for the longest time. Additionally, stirring is less effective in the high viscosity layer of the central lower mantle [1, 2], supporting the survival of medium-scale heterogeneities there. These models are the first, using 3-d spherical geometry and mostly Earth-like parameters, to address the suggested change of convective style. Although the models are still far from reproducing our planet, we find that proposal might be helpful towards reconciling geochemical and geophysical constraints.

A refined solution to Earth's hidden niobium: implications for evolution of continental crust and mode of core formation

New high-precision niobium (Nb) and tantalum (Ta) concentration data are presented for early Archaean metabasalts, metabasaltic komatiites and their erosion products (mafic metapelites) from SW Greenland and the Acasta gneiss complex, Canada. Individual datasets consistently show sub-chondritic Nb/Ta ratios averaging 15.1+/-11.6. This finding is discussed with regard to two competing models for the solution of the Nb-deficit that characterises the accessible Earth. Firstly, we test whether Nb could have sequestered into the core due to its slightly siderophile (or chalcophile) character under very reducing conditions, as recently proposed from experimental evidence. We demonstrate that troilite inclusions of the Canyon Diablo iron meteorite have Nb and V concentrations in excess of typical chondrites but that the metal phase of the Grant, Toluca and Canyon Diablo iron meteorites do not have significant concentrations of these lithophile elements. We find that if the entire accessible Earth Nb-deficit were explained by Nb in the core, only ca. 17% of the mantle could be depleted and that by 3.7 Ga, continental crust would have already achieved ca. 50% of its present mass. Nb/Ta systematics of late Archaean metabasalts compiled from the literature would further require that by 2.5 Ga, 90% of the present mass of continental crust was already in existence. As an alternative to this explanation, we propose that the average Nb/Ta ratio (15.1+/-11.6) of Earth's oldest mafic rocks is a valid approximation for bulk silicate Earth. This would require that ca. 13% of the terrestrial Nb resided in the Ta-free core. Since the partitioning of Nb between silicate and metal melts depends largely on oxygen fugacity and pressure, this finding could mean that metal/silicate segregation did not occur at the base of a deep magma ocean or that the early mantle was slightly less reducing than generally assumed. A bulk silicate Earth Nb/Ta ratio of 15.1 allows for depletion of up to 40% of the total mantle. This could indicate that in addition to the upper mantle, a portion of the lower mantle is depleted also, or if only the upper mantle were depleted, an additional hidden high Nb/Ta reservoir must exist. Comparison of Nb/Ta systematics between early and late Archaean metabasalts supports the latter idea and indicates deeply subducted high Nb/Ta eclogite slabs could reside in the mantle transition zone or the lower mantle. Accumulation of such slabs appears to have commenced between 2.5 and 2.0 Ga. Regardless of these complexities of terrestrial Nb/Ta systematics, it is shown that the depleted mantle Nb/Th ratio is a very robust proxy for the amount of extracted continental crust, because the temporal evolution of this ratio is dominated by Th-loss to the continents and not Nb-retention in the mantle. We present a new parameterisation of the continental crust volume versus age curve that specifically explores the possibility of lithophile element loss to the core and storage of eclogite slabs in the transition zone. (C) 2003 Elsevier Science B.V. All rights reserved.

The source of the Great Dyke, Zimbabwe, and its tectonic significance: Evidence from Re-Os isotopes

Re-Os data for chromite separates from 10 massive chromitite seams sampled along the 550-km length of the 2.58-Ga Great Dyke layered igneous complex, Zimbabwe, record initial 187Os/188Os ratios in the relatively narrow range between 0.1106 and 0.1126. This range of initial 187Os/188Os values is only slightly higher than the value for the coeval primitive upper mantle (0.1107) as modeled from the Re-Os evolution of chondrites and data of modern mantle melts and mantle derived xenoliths. Analyses of Archean granitoid and gneiss samples from the Zimbabwe Craton show extremely low Os concentrations (3-9 ppt) with surprisingly unradiogenic present-day 187Os/188Os signatures between 0.167 and 0.297. Only one sample yields an elevated 187Os/188Os ratio of 1.008. Using these data, the range of crustal contamination of the Great Dyke magma would be minimally 0%-33% if the magma source was the primitive upper mantle, whereas the range estimated from Nd and Pb isotope systematics is 5%-25%. If it is assumed that the primary Great Dyke magma derived from an enriched deep mantle reservoir (via a plume), a better agreement can be obtained. A significant contribution from a long-lived subcontinental lithospheric mantle (SCLM) reservoir with subchondritic Re/Os to the Great Dyke melts cannot be reconciled with the Os isotope results at all. However, Os isotope data on pre-Great Dyke ultramafic complexes of the Zimbabwe Craton and thermal modeling show that such an SCLM existed below the Zimbabwe Craton at the time of the Great Dyke intrusion. It is therefore concluded that large melt volumes such as that giving rise to the Great Dyke were able to pass lithospheric mantle keels without significant contamination in the late Archean. Because the ultramafic-mafic melts forming the Great Dyke must have originated below the SCLM (which extends to at least a 200-km depth ), the absence of an SCLM signature precludes a subduction-related magma-generation process.

Emergent anisotropy and flow alignment in viscous rock

A novel class of nonlinear, visco-elastic rheologies has recently been developed by MUHLHAUS et al. (2002a, b). The theory was originally developed for the simulation of large deformation processes including folding and kinking in multi-layered visco-elastic rock. The orientation of the layer surfaces or slip planes in the context of crystallographic slip is determined by the normal vector the so-called director of these surfaces. Here the model (MUHLHAUS et al., 2002a, b) is generalized to include thermal effects; it is shown that in 2-D steady states the director is given by the gradient of the flow potential. The model is applied to anisotropic simple shear where the directors are initially parallel to the shear direction. The relative effects of textural hardening and thermal softening are demonstrated. We then turn to natural convection and compare the time evolution and approximately steady states of isotropic and anisotropic convection for a Rayleigh number Ra=5.64x10(5) for aspect ratios of the experimental domain of 1 and 2, respectively. The isotropic case has a simple steady-state solution, whereas in the orthotropic convection model patterns evolve continuously in the core of the convection cell, which makes only a near-steady condition possible. This near-steady state condition shows well aligned boundary layers, and the number of convection cells which develop appears to be reduced in the orthotropic case. At the moderate Rayleigh numbers explored here we found only minor influences in the change from aspect ratio one to two in the model domain.

Histological structure of the cancellous bone layer in Bothriolepis canadensis (Antiarchi, Placodermi)

The Placodermi are extinct basal gnathostomes which had extensive dermal and perichondral bone, but which lacked the endochondral bone which characterizes the more derived bony fishes. Thin sections of bone from a specimen of the antiarch placoderm Bothriolepis canadensis, from the Escuminac Formation ( Frasnian, Upper Devonian), Quebec, Canada, reveal that part of the cancellous layer in its dermal and endoskeletal bone formed from perichondral bone trabeculae growing around cartilage spheres. The resultant structure mimics that of osteichthyan endochondral bone. The layout and dimensions of this polygonal mosaic patterning of the bone trabeculae and flattened cartilage spheres resemble those of the prismatic layers of calcified cartilage in chondrichthyans. If the lack of endoskeletal bone in chondrichthyans is a derived character, then the structure identified in B. canadensis could represent a 'template' for the formation of prismatic calcified cartilage in the absence of bone.

New dinoflagellate cyst species from the upper Callovian-lower Oxfordian Rigaudella aemula zone, Timor sea, northwestern Australia

Four new species of dinoflagellate cysts are described from Callovian to lower Oxfordian (Jurassic) sediments of the Timor Sea, northwestern Australia. These comprise Evansia? lacryma, Egmontodinium elongatum, Leptodinium? ancoralium, and Nannoceratopsis reticulata. They are rare to common constituents of the Rigaudella aemula dinoflagellate cyst Interval Zone, and may prove useful for regional biostratigraphic correlation. (c) 2005 Elsevier B.V. All rights reserved.

Optimal sampling strategy for estimation of spatial genetic structure in tree populations

Fine-scale spatial genetic structure (SGS) in natural tree populations is largely a result of restricted pollen and seed dispersal. Understanding the link between limitations to dispersal in gene vectors and SGS is of key interest to biologists and the availability of highly variable molecular markers has facilitated fine-scale analysis of populations. However, estimation of SGS may depend strongly on the type of genetic marker and sampling strategy (of both loci and individuals). To explore sampling limits, we created a model population with simulated distributions of dominant and codominant alleles, resulting from natural regeneration with restricted gene flow. SGS estimates from subsamples (simulating collection and analysis with amplified fragment length polymorphism (AFLP) and microsatellite markers) were correlated with the 'real' estimate (from the full model population). For both marker types, sampling ranges were evident, with lower limits below which estimation was poorly correlated and upper limits above which sampling became inefficient. Lower limits (correlation of 0.9) were 100 individuals, 10 loci for microsatellites and 150 individuals, 100 loci for AFLPs. Upper limits were 200 individuals, five loci for microsatellites and 200 individuals, 100 loci for AFLPs. The limits indicated by simulation were compared with data sets from real species. Instances where sampling effort had been either insufficient or inefficient were identified. The model results should form practical boundaries for studies aiming to detect SGS. However, greater sample sizes will be required in cases where SGS is weaker than for our simulated population, for example, in species with effective pollen/seed dispersal mechanisms.

A new estimate for the composition of weathered young upper continental crust from alluvial sediments, Queensland, Australia

We present new major element, trace element and Nd-isotope data for 30 alluvial sediments collected from 25 rivers in Queensland, E Australia. Samples were chosen to represent drainage from the region's most important lithologies, including Tertiary intraplate volcanic rocks, a Cretaceous igneous province (and sedimentary rocks derived thereof) as well as Proterozoic blocks. In most chemical and isotopic aspects, the alluvial sediments represent binary or ternary mixing relationships, with absolute abundances implied to reflect the proportion of lithologies in the catchments. When averaged, the studied sediments differ from other proxies of upper continental crust (UCC) mainly in their relative middle rare earth element enrichment (including an elevated Sm/Nd ratio), higher relative Eu abundance and higher Nb/Ta ratio. These features are inherited from eroded Tertiary intraplate basalts, which commonly form topographic highs in the studied region. Despite the high degree of weathering strong to excellent coherence between similarly incompatible elements is found for all samples. From this coherence, we suggest revisions of the following upper crustal element ratios: Y/Ho = 26.2, Yb/Tm = 6.37, Th/W = 7.14, Th/Tl = 24 and Zr/Hf = 36.9. Lithium, Rb, Cs and Be contents do not seem depleted relative to UCC, which may reflect paucity of K-feldspar in the eroded catchments. Nickel, Cr, Pb, Cu and Zn concentrations are elevated in polluted rivers surrounding the state capital. River sediments in the Proterozoic Georgetown Inlier are elevated in Pb, Cu and Zn but this could be a natural phenomenon reflecting abundant sulphide mineralisation in the area. Except for relative Sr concentrations, which broadly anticorrelate with mean annual rainfall in catchments, there is no obvious relationship between the extent of weathering and climate types, which range from and to tropical. The most likely explanation for this observation is that the weathering profiles in many catchments are several Myr old, established during the much wetter Miocene period. The studied sediment compositions (excluding those from the Proterozoic catchments) are used to propose a new trace element normalisation termed MUQ (MUd from Queensland), which serves as an alternative to UCC proxies derived from sedimentary rocks. Copyright (C) 2005 Elsevier Ltd

Theoretical and numerical analysis of large-scale heat transfer problems with temperature-dependent pore-fluid densities

Purpose - In many scientific and engineering fields, large-scale heat transfer problems with temperature-dependent pore-fluid densities are commonly encountered. For example, heat transfer from the mantle into the upper crust of the Earth is a typical problem of them. The main purpose of this paper is to develop and present a new combined methodology to solve large-scale heat transfer problems with temperature-dependent pore-fluid densities in the lithosphere and crust scales. Design/methodology/approach - The theoretical approach is used to determine the thickness and the related thermal boundary conditions of the continental crust on the lithospheric scale, so that some important information can be provided accurately for establishing a numerical model of the crustal scale. The numerical approach is then used to simulate the detailed structures and complicated geometries of the continental crust on the crustal scale. The main advantage in using the proposed combination method of the theoretical and numerical approaches is that if the thermal distribution in the crust is of the primary interest, the use of a reasonable numerical model on the crustal scale can result in a significant reduction in computer efforts. Findings - From the ore body formation and mineralization points of view, the present analytical and numerical solutions have demonstrated that the conductive-and-advective lithosphere with variable pore-fluid density is the most favorite lithosphere because it may result in the thinnest lithosphere so that the temperature at the near surface of the crust can be hot enough to generate the shallow ore deposits there. The upward throughflow (i.e. mantle mass flux) can have a significant effect on the thermal structure within the lithosphere. In addition, the emplacement of hot materials from the mantle may further reduce the thickness of the lithosphere. Originality/value - The present analytical solutions can be used to: validate numerical methods for solving large-scale heat transfer problems; provide correct thermal boundary conditions for numerically solving ore body formation and mineralization problems on the crustal scale; and investigate the fundamental issues related to thermal distributions within the lithosphere. The proposed finite element analysis can be effectively used to consider the geometrical and material complexities of large-scale heat transfer problems with temperature-dependent fluid densities.

Upper Devonian biostratigraphy of the Iberian Pyrite Belt, southwest Spain Part One: Miospores

The Iberian Pyrite Belt (IPB), which forms part of the Variscan orogenic massif, is renowned for the magnitude and extent of its massive sulfide mineralization. The stratigraphic record of the IPB consists of Upper Palaeozoic sedimentary and igneous rocks. In ascending order, these comprise the thick Phyllite-Quartzite Group attributed to the Middle and Upper Devonian and characterized by shales and quartzites with conglomeratic and carbonate intercalations towards the top; the appreciably thinner Volcano-Sedimentary Complex, a heterogeneous uppermost Devonian-Mississippian unit embodying diverse volcanic, subvolcanic, and sedimentary rocks that host the massive sulfide deposits; and the shaly and sandy, turbiditic Culm Group (Carboniferous). This entire succession was folded and faulted during the Asturian phase of the Variscan Orogeny that gave rise to a thin-skinned type structure. The present study constitutes a detailed blostratigraphic investigation of palynologically productive samples representative of the Phyllite-Quartzite Group and the basal (anoxic) portion of the Volcano-Sedimentary Complex. These were collected from surface and mine exposures variously located in the Spanish part of the IPB; out of 282 samples processed, 117 proved to be productive palynologically. The aim of this project is to provide comprehensive palynostratigraphic data applicable to precise dating and correlation of the IPB's stratigraphic succession (i.e., of the two sampled lithostratigraphic units), which has hitherto been investigated biostratigraphically on a relatively localized basis. The results are incorporated in two successive parts. The first of these, i. e., the present paper, focuses on the systematic analysis of the terrestrial (miospore) component of the palynological assemblages. The second part, devoted to the marine, organic-walled microphytoplankton (acritarchs and prasinophytes), will evaluate the stratigraphic significance of the IPB palynofloras and their application to elucidating the geological history of the region. In the systematic-descriptive section, which occupies the bulk of this paper, 55 species of trilete miospores are described and are allocated among 34 genera, two of which (Cristicavatispora and Epigruspora) are newly instituted herein. The majority of the species are either positively identifiable or closely affiliable with previously named species. The nine newly established species are as follows: Camptozonotriletes confertus, Indotriradites diversispinosus, Cristicavatispora dispersa (type species), Epigruspora regularis (type species), Ancyrospora? implicata, Endosporites tuberosus, Rugospora explicata, Spelaeotriletes plicatus, and Teichertospora iberica.

Microsensors technology used to measure productivity and other biochemical processes in the upper regions of aquatic sediments

Carbonate sediments are dynamic three-dimensional environments where the surface layers are constantly moving and mixing due to the energy of the water column. It is also an environment of dynamic biological, chemical and physical interaction and modification. The biological community can actively influence changes to sediment characteristics and associated biochemistry. Bioturbation resulting from macrofaunal activity disrupts sediment structure and biochemical arrangements and reduces the critical shear forces required to move sediment particles, adding to the dynamic and complex physical and biogeochemical nature of the sediment. Laboratory studies using both planner optodes and glass needle microsensors were used to measure abiotic sediment characteristics such as the depth distribution and concentrations of PAR. The biochemical nature of coral reef sediment were also investigated, specifically the quantification and the distribution of dissolved oxygen within coarse and fine-grained sediments under regimes of light and darkness. Results highlighted the significant contribution microalgal productivity and bioturbation has on distribution of dissolved oxygen in the upper sediment layers. On the reef flat a shallow water lander system was employed to measure concentrations of O2, pH, S, Ca and temperature over periods of 24 to 48 hours in coarse and fine-grained sediments. Similarities between laboratory and in situ results where evident, however the in situ environment was more dynamic and the distribution and concentrations of dissolved oxygen were more complex and correlated to irradiance, temperature and biological activity. Microsensor technology provides us with the opportunity to study, at very high resolutions, the upper irradiated; photosynthetically active regions of aquatic sediments along with anoxic processes deeper in sub-euphotic regions of the sediments.