Investigating the mesh dependency and upscaling of 3D grain-based models for the simulation of brittle fracture processes in low-porosity crystalline rock

The application of 3D grain-based modelling techniques is investigated in both small and large scale 3DEC models, in order to simulate brittle fracture processes in low-porosity crystalline rock. Mesh dependency in 3D grain-based models (GBMs) is examined through a number of cases to compare Voronoi and tetrahedral grain assemblages. Various methods are used in the generation of tessellations, each with a number of issues and advantages. A number of comparative UCS test simulations capture the distinct failure mechanisms, strength profiles, and progressive damage development using various Voronoi and tetrahedral GBMs. Relative calibration requirements are outlined to generate similar macro-strength and damage profiles for all the models. The results confirmed a number of inherent model behaviors that arise due to mesh dependency. In Voronoi models, inherent tensile failure mechanisms are produced by internal wedging and rotation of Voronoi grains. This results in a combined dependence on frictional and cohesive strength. In tetrahedral models, increased kinematic freedom of grains and an abundance of straight, connected failure pathways causes a preference for shear failure. This results in an inability to develop significant normal stresses causing cohesional strength dependence. In general, Voronoi models require high relative contact tensile strength values, with lower contact stiffness and contact cohesional strength compared to tetrahedral tessellations. Upscaling of 3D GBMs is investigated for both Voronoi and tetrahedral tessellations using a case study from the AECL’s Mine-by-Experiment at the Underground Research Laboratory. An upscaled tetrahedral model was able to reasonably simulate damage development in the roof forming a notch geometry by adjusting the cohesive strength. An upscaled Voronoi model underestimated the damage development in the roof and floor, and overestimated the damage in the side-walls. This was attributed to the discretization resolution limitations.

Analysis of optimal control problems for the incompressible MHD equations and implementation in a finite element multiphysics code

This thesis deals with the study of optimal control problems for the incompressible Magnetohydrodynamics (MHD) equations. Particular attention to these problems arises from several applications in science and engineering, such as fission nuclear reactors with liquid metal coolant and aluminum casting in metallurgy. In such applications it is of great interest to achieve the control on the fluid state variables through the action of the magnetic Lorentz force. In this thesis we investigate a class of boundary optimal control problems, in which the flow is controlled through the boundary conditions of the magnetic field. Due to their complexity, these problems present various challenges in the definition of an adequate solution approach, both from a theoretical and from a computational point of view. In this thesis we propose a new boundary control approach, based on lifting functions of the boundary conditions, which yields both theoretical and numerical advantages. With the introduction of lifting functions, boundary control problems can be formulated as extended distributed problems. We consider a systematic mathematical formulation of these problems in terms of the minimization of a cost functional constrained by the MHD equations. The existence of a solution to the flow equations and to the optimal control problem are shown. The Lagrange multiplier technique is used to derive an optimality system from which candidate solutions for the control problem can be obtained. In order to achieve the numerical solution of this system, a finite element approximation is considered for the discretization together with an appropriate gradient-type algorithm. A finite element object-oriented library has been developed to obtain a parallel and multigrid computational implementation of the optimality system based on a multiphysics approach. Numerical results of two- and three-dimensional computations show that a possible minimum for the control problem can be computed in a robust and accurate manner.

Keratinocyte transcriptional regulation of the human c-Myc promoter occurs via a novel Lef/Tcf binding element distinct from neoplastic cells

The proto-oncogene c-Myc is involved in early neoplastic transformations. Two consensus Lef/Tcf binding elements (TBE) were found to be prerequisite for transcriptional transactivation by the armadillo proteins beta-catenin and plakoglobin (PG) together with Tcf4 in human neoplastic cells. In epidermal keratinocytes, c-Myc was reported to be repressed by Lef-1 and PG. Using reporter gene assays, here we demonstrate that deletion of the two consensus TBE fails to abrogate transcriptional regulation by Lef-1/PG in wildtype and beta-catenin-/- keratinocytes, while it reduces transcription in pre-neoplastic PG-/- keratinocytes. We identified a TBE sequence variant downstream of the major transcriptional initiation site that binds Lef-1 in vitro and in vivo, and its mutation compromised transcriptional regulation by Lef-1/PG. Collectively, this study demonstrates that the two consensus TBE's reported in neoplastic cells are dispensable for c-Myc regulation in normal keratinocytes, which instead use a novel TBE sequence variant. This unprecedented finding may have important implications for armadillo target genes involved in carcinogenesis.

Element concentrations of plagioclase and clinopyroxene in poikilitic olivine gabbros from ODP Hole 153-923A (Table 4)

Ocean Drilling Program Hole 923A, located on the western flank of the Mid-Atlantic Ridge south of the Kane Fracture Zone, recovered primitive gabbros that have mineral trace element compositions inconsistent with growth from a single parental melt. Plagioclase crystals commonly show embayed anorthitic cores overgrown by more albitic rims. Ion probe analyses of plagioclase cores and rims show consistent differences in trace element ratios, indicating variation in the trace element characteristics of their respective parental melts. This requires the existence of at least two distinct melt compositions within the crust during the generation of these gabbros. Melt compositions calculated to be parental to plagioclase cores are depleted in light rare earth elements, but enriched in yttrium, compared to basalts from this region of the Mid-Atlantic Ridge, which are normal mid-ocean ridge basalt (N-MORB). Clinopyroxene trace element compositions are similar to those predicted to be in equilibrium with N-MORB. However, primitive clinopyroxene crystals are much more magnesian than those produced in one-atmosphere experiments on N-MORB, suggesting that the major element composition of the melt was unlike N-MORB. These data require that the diverse array of melt compositions generated within the mantle beneath mid-ocean ridges are not always fully homogenised during melt extraction from the mantle and that the final stage of mixing can occur efficiently within crustal magma chambers. This has implications for the process of melt extraction from the mantle and the liquid line of descent of MORB

Major and trace element analyses of ODP Leg 120 volcanic rocks (Table 1)

Widespread Lower Cretaceous magmatism occurred along the Indian-Australian/Antarctic margins, and in the juvenile Indian Ocean, during the rifting of eastern Gondwana. The formation of this magmatic province probably began around 120-130 Ma with the eruption of basalts on the Naturaliste Plateau and at Bunbury, western Australia. On the northeast margin of India, activity began around 117 Ma with the Rajmahal continental basalts and associated lamprophyre intrusions. The formation of the Kerguelen Plateau in the Indian Ocean began no later than 114 Ma. Ultramafic lamprophyres (alnoites) were emplaced in the Prince Charles Mountains near the Antarctic continental margin at ~ 110 Ma. These events are considered to be related to a major mantle plume, the remnant of which is situated beneath the region of Kerguelen and Heard islands at the present day. Geochemical data are presented for each of these volcanic suites and are indicative of complex interactions between asthenosphere-derived magmas and the continental lithosphere. Kerguelen Plateau basalts have Sr and Nd isotopic compositions lying outside the field for Indian Ocean mid-ocean ridge basalts (MORB) but, with the exception of Site 738 at the southern end of the plateau, within the range of more recent hotspot basalts from Kerguelen and Heard Islands. However, a number of the plateau tholeiites are characterized by lower 206Pb/204Pb ratios than are basalts from Kerguelen Island, and many also have anomalously high La/Nb ratios. These features suggest that the source of the Kerguelen Plateau basalts suffered contamination by components derived from the Gondwana continental lithosphere. An extreme expression of this lithospheric signature is shown by a tholeiite from Site 738, suggesting that the southernmost part of the Kerguelen Plateau may be underlain by continental crust. The Rajmahal tholeiites mostly fall into two distinct geochemical groups. Some Group I tholeiites have Sr and Nd isotopic compositions and incompatible element abundances, similar to Kerguelen Plateau tholeiites from Sites 749 and 750, indicating that the Kerguelen-Heard mantle plume may have directly furnished Rajmahal volcanism. However, their elevated 207Pb/204Pb ratios indicate that these magmas did not totally escape contamination by continental lithosphere. In contrast to the Group I tholeiites, significant contamination is suggested for Group II Rajmahal tholeiites, on the basis of incompatible element abundances and isotopic compositions. The Naturaliste Plateau and the Bunbury Basalt samples show varying degrees of enrichment in incompatible elements over normal MORB. The Naturaliste Plateau samples (and Bunbury Basalt) have high La/Nb ratios, a feature not inconsistent with the notion that the plateau may consist of stretched continental lithosphere, near the ocean-continent divide.

Major and trace element compositions of the Buendnerschiefer samples

The Bündnerschiefer of the Swiss-Italian Alps is a large sedimentary complex deposited on the Piemonte-Liguria and Valais oceans and associated continental margins from the upper Jurassic to Eocene. It is made of a large variety of sequences associated or not with an ophiolitic basement. The Bündnerschiefer makes an accretionary prism that developed syn-tectonically from the onset of alpine subduction, and it records orogenic metamorphism following episodes of HP metamorphism. The Bündnerschiefer shares important similarities with the Otago schists of New Zealand and with the Wepawaug schists of Connecticut, both of which form accretionary prisms and have an orogenic metamorphic imprint. With the aim of testing the hypothesis of mobility of chemical components as a function of metamorphic grade, in this work I present fifty-five bulk chemical analyses of various lithological facies of the Bündnerschiefer collected along the well-studied field gradient of the Lepontine dome of Central Switzerland, in the Prättigau half window of East Switzerland, and in the Tsaté Nappe of Valle d'Aosta (Italy). The dataset includes the concentration of major components, large ion lithophile elements (Rb, Sr, Ba, Cs), high field strength elements (Zr, Ti, Nb, Th, U, Ta, Hf), fluid-mobile light elements (B, Li), volatiles (CO2, S), REEs, and Y, V, Cr, Co, Sn, Pb, Cu, Zn, Tl, Sb, Be, and Au. These data are compared against the compositions of the global marine sediment reservoir, typical crustal reservoirs, and against the previously measured compositions of Otago and Wepawaug schists. Results reveal that, irrespective of their metamorphic evolution, the bulk chemical compositions of orogenic metasediments are characterized by mostly constant compositional ratios (e.g., K2O/Al2O3, Ba/Al2O3, Sr/CaO, etc.), whose values in most cases are undistinguishable from those of actual marine sediments and other crustal reservoirs. For these rocks, only volatile concentrations decrease dramatically as a function of metamorphic temperature, and significant deviations from the reservoir signatures are evident for SiO2, B, and Li. These results are interpreted as an indication of residual enrichment in the sediments, a process taking place during syn-metamorphic dehydration from the onset of metamorphism in a regime of chemical immobility. Residual enrichment increased the absolute concentrations of the chemical components of these rocks, but did not modify significantly their fundamental ratios. This poor compositional modification of the sediments indicates that orogenic metamorphism in general does not promote significant mass transfer from accretionary prisms. In contrast, mass transfer calculations carried out in a shear zone crosscutting the Bündnerschiefer shows that significant mass transfer occurs within these narrow zones, resulting in gains of H2O, SiO2, Al2O3, K2O, Ba, Y, Rb, Cu, V, Tl, Mo, and Ce during deformation and loss of Na2O, CO2, S, Ni, B, U, and Pb from the rock. These components were presumably transported by an aquo-carbonic fluid along the shear zone. These distinct attitudes to mobilize chemical elements from orogenic sediments may have implications for a potentially large number of geochemical processes in active continental margins, from the recycling of chemical components at plate margins to the genesis of hydrothermal ore deposits.

Element oxides in basaltic glass fragments from ODP Hole 203-1243B

Sand-sized basaltic glass fragments were recovered in the liner of Core 203-1243B-19R, the deepest recovery from Hole 1243B. Microprobe analysis of 582 glassy cuttings cluster into five compositionally distinct groups, most of which are unlike the lithologic units described on board ship. Drilling operations intended to sweep cuttings from the caving hole and differences between the cuttings and geochemically distinct lithologic units of the upper part of the basement indicate that the cuttings came mainly, if not entirely, from the lower part of the hole. They give information about the part of Hole 1243B that had poor core recovery. Enriched mid-ocean-ridge basalt (MORB) from the upper part of the hole and transitional MORB from two groups of cuttings from sources low in the hole may be a trace of the Galápagos plume on the Pacific plate or may be a normal consequence of eruptions from two distinct magmas on fast-spreading crust.

Mineral composition and major element oxides of gabbros from ODP Hole 179-1105A, southwest Indian Ridge

Hole 1105A penetrated 158 m of gabbros at a site offset 1.3 km east-northeast from Hole 735B on the Atlantis Bank near the Atlantis II Fracture Zone. A total of 118 m of dominantly medium- to coarse-grained intercalated Fe-Ti oxide gabbro and olivine gabbro was recovered from Hole 1105A that shows many petrographic features similar to those recovered from the upper part of Hole 735B. The main rock types are distinguished based on the constituent cumulus phases, with the most primitive gabbros consisting of olivine, plagioclase, and clinopyroxene. The inferred crystallization order is subsequently Fe-Ti oxides (ilmenite and titanomagnetite), followed by orthopyroxene, then apatite, and finally biotite. Orthopyroxene appears to replace olivine in a narrow middle interval. The magmatic evolution is likewise reflected in the mineral compositions. Plagioclase varies from An66 to An28. Olivine varies from Fo78 to Fo35. The gap in olivine crystallization occurs between Fo46 and Fo40 and coincides approximately with the appearance of orthopyroxene (~En50). The clinopyroxenes show large compositional variation in Mg/(Mg + Fe total) from 0.84 to 0.51. The nonquadrilateral cations of clinopyroxene similarly show large variations with Ti increasing for the olivine gabbros and decreasing for the Fe-Ti oxide gabbros with the decrease in Mg/(Mg + Fe total). The apatites are mainly flourapatites. The compositional variation in the gabbros is interpreted as a comagmatic suite resulting from fractional crystallization. Pyroxene geothermometry suggests equilibration temperatures from 1100°C and below. The coexisting Fe-Ti oxide minerals indicate subsolidus equilibration temperatures from 900°C for olivine gabbros to 700°C for the most evolved apatite-bearing gabbros. The cryptic variation in the olivine gabbros defines two or three lenses, 40 to 60 m thick, each characterized by a distinct convex zoning with a lower segment indicating upward reverse fractionation, a central maximum, and an upper segment showing normal fractionation. The Fe-Ti oxide gabbros show cryptic variations independent of the host olivine gabbros and reveal a systematic upward normal fractionation trend transgressing host olivine gabbro boundaries. Forward fractional crystallization modeling, using a likely parental magma composition from the Atlantis II Fracture Zone (MgO = 7.2 wt%; Mg/[Mg + Fe2+] = 0.62), closely matches the compositions of coexisting olivine, plagioclase, and clinopyroxene. This modeling suggests cosaturation of olivine, plagioclase, and clinopyroxene from 1155°C and the addition of Fe-Ti oxides from 1100°C. The liquid line of descent initially shows increasing FeO with moderately increasing SiO2. After saturation of Fe-Ti oxides, the liquid strongly decreases in FeO and TiO2 and increases in SiO2, reaching dacitic compositions at ~10% liquid remaining. The calculations indicate that formation of olivine gabbros can be accounted for by <65% fractionation and that only the residual 35% liquid was saturated in Fe-Ti oxides. The modeling of the solid fractionation products shows that both the olivine gabbro and the Fe-Ti oxide gabbros contain very small amounts of trapped liquid (<5%). The implications are that the gabbros represent crystal mush that originated in a recharging and tapping subaxial chamber. Compaction and upward melt migration in the crystal mush appear to have been terminated with relatively large amounts of interstitial liquid remaining in the upper parts of the cumulate mush. This termination may have been caused by tectonic disturbances, uplift, and associated withdrawal of magma into the subaxial dike and sill system. Prolonged compaction and cooling of the trapped melt in the mush formed small differentiated bodies and lenses by pressure release migration and crystallization along syntectonic channels. This resulted in differentiation products along lateral and vertical channelways in the host gabbro that vary from olivine gabbro, to Fe-Ti oxide gabbro, gabbronorite, and apatite gabbros and show large compositional variations independent of the host olivine gabbros.

Element concentrations and isotopic composition of tholeiite and whole rocks samples from the east Mariana Basin

Studies of seafloor magnetic anomaly patterns suggest the presence of Jurassic oceanic crust in a large area in the western Pacific that includes the East Mariana, Nauru and Pigafetta Basins. Sampling of the igneous crust in this area by the Deep Sea Drilling Program (DSDP) and the Ocean Drilling Program (ODP) allows direct evaluation of the age and petrogenesis of this crust. ODP Leg 129 drilled a 51 m sequence of basalt pillows and massive flows in the central East Mariana Basin. 40Ar/39Ar ages determined in this study for two Leg 129 basalts average 114.6 +/- 3.2 Ma. This age is in agreement with the Albian-late Aptian paleontologic age of the overlying sediments, but is distinctively younger than the Jurassic age predicted by magnetic anomaly patterns in the basin. Compositionally, the East Mariana Basin basalts are uniformly low-K tholeiites that are depleted in highly incompatible elements compared to moderately incompatible ones, which is typical of mid-ocean ridge basalts (MORB) erupted near hotspots. The Sr, Nd and Pb isotopic compositions of the tholeiites (87Sr/86Sr init = 0.70360-0.70374; 143Nd/144Nd init = 0.512769-0.512790; 206Pb/204Pb meas = 18.355-18.386) also overlap with some Indian Ocean Ridge MORB, although they are distinct from the isotopic compositions of Jurassic basalts drilled in the Pigafetta Basin, the oldest Pacific MORB. The isotopic compositions of the East Mariana Basin tholeiites are also similar to those of intraplate basalts, and in particular, to the isotopic signature of basalts from the nearby Ontong Java and Manihiki Plateaus. The East Mariana Basin tholeiites also share many petrologic and isotopic characteristics with the oceanic basement drilled in the Nauru Basin at DSDP Site 462. In addition, the new 110.8 +/- 1.0 Ma 40Ar/39Ar age for two flows from the bottom of Site 462 in the Nauru Basin is indistinguishable from the age of the East Mariana Basin flows. Thus, while magnetic anomaly patterns predict that the igneous basement in the Nauru and East Mariana Basins is Jurassic in age, the geochemical and chronological results discussed here suggest that the basement formed during a Cretaceous rifting event within the Jurassic crust. This magmatic and tectonic event was created by the widespread volcanism responsible for the genesis of the large oceanic plateaus of the western Pacific.

Mineralogy and major element oxide composition of sediments from ODP Hole 172-1063D

Dansgaard-Oeschger (D-O) cycles in sediment at Site 1063 are characterized by distinct fluctuations in physical properties. Stadials are marked by low bulk density and interstadials by high bulk density. Compressional (P-)wave velocity is in phase with bulk density over some but not all depth intervals. Four of the D-O cycles straddling the oxygen isotope Stage 4/5 boundary have been studied in detail to understand the origin of the physical properties changes. Sediment on the Bermuda Rise is comprised of three main components: calcite, aluminosilicate minerals, and biogenic silica. Calcite concentrations vary from 1% to 43% of bulk sediment and are highest during interstadials. Aluminosilicate concentrations vary from 52% to 92% of bulk sediment and are highest during stadials. The major element ratios Al2O3/TiO2 and K2O/Al2O3 show increases across bulk density cycles, suggesting a change in the composition of aluminosilicates. This interpretation is supported by mineralogical analyses, which show a subtle change in clay composition. Biogenic silica concentrations vary from 0% to 23% of bulk sediment and are also highest during stadials. However, the abundance of silica varies significantly from one D-O cycle to another. Silt and fine sand abundance also increase during the first of the four stadials. This coarsening of sediment coincides with the increase in biogenic silica. The low grain density and high porosity associated with biogenic silica result in intervals of low bulk-sediment density. The abundance of biogenic silica closely matches P-wave velocity, suggesting that silica imparts a greater rigidity to the sediment.

Major element oxides in ODP Site 164-997 sediments

Since being first discovered in the Blake-Bahama region of the west Atlantic in the 1970s (Hollister, Ewing, et al., 1972, doi:10.2973/dsdp.proc.11.1972), submarine gas hydrates have been identified in the continental margin worldwide. Ocean Drilling Program (ODP) Leg 164 was the first drilling designated to study the occurrence and distribution of natural gas hydrates in Blake Ridge where a well developed, distinct BSR (Bottom Simulating Reflector) has been identified (Paull, Matsumoto, Wallace, et al., 1996, doi:10.2973/odp.proc.ir.164.1996). It has been reported there is a prominent discrepancy between the BSR and the base of gas hydrate stability (Paull, Matsumoto, Wallace, et al., 1996, doi:10.2973/odp.proc.ir.164.1996; Ruppel, 1997, doi:10.1130/0091-7613(1997)025<0699:ACTOAT>2.3.CO;2), though theoretically they should be at the same depth. Natural gas hydrate in marine sediments coexists with sediment particles, so detailed delineation of sediment geochemistry will be of benefit to solve this apparent discrepancy. The main objectives of this study are to supply background data of the major chemical compositions of sediments from a hydrated sediment section.

Trace-element geochemistry of tholeiitic basalts from DSDP Hole 55-433C

Hole 433C, a multiple re-entry hole drilled in 1862 meters of water on Suiko Seamount in the central Emperor Seamounts, penetrated 387.5 meters of lava flows overlain by 163.0 meters of sediments. The recovered volcanic rocks consist of three flow units (1-3) of alkalic basalt underlain by more than 105 flows or flow lobes (Flow Units 4-67) of tholeiitic basalt. This study reports trace-element, including rare-earth element (REE), data for 25 samples from 24 of the least altered tholeiitic flows. These data are used to evaluate the origin and evolution of tholeiitic basalts from Suiko Seamount and to evaluate changes in the mantle source between the time when Suiko Seamount formed, 64.7 ± 1.1 m.y. ago (see Dalrymple et al., 1980), and the present day. Stearns (1946), Macdonald and Katsura (1964) and Macdonald (1968) have established that chemically distinct lavas erupt during four eruptive stages of development of a Hawaiian volcano. These stages, from initial to final, are shield-building, caldera-filling, post-caldera, and post-erosional. The lavas of the shield-building stage are tholeiitic basalts, which erupt rapidly and in great volume. The shield-building stage is quickly followed by caldera collapse and by the caldera-filling stage, during which the caldera is filled by tholeiitic and alkalic lavas. During the post-caldera stage, a relatively thin veneer of alkalic basalts and associated differentiated lavas are erupted, sometimes accompanied by minor eruptions of tholeiitic lava. After a period of volcanic quiescence and erosion, lavas of the nephelinitic suite, which include both alkalic basalts and strongly SiO2-undersaturated nephelinitic basalts, may erupt from satellite vents during the post-erosional stage. Many Hawaiian volcanoes develop through all four stages; but individual volcanoes have become extinct before the cycle is complete. We interpret the tholeiitic lavas drilled on Suiko Seamount to have erupted during either the shield-building or the caldera-filling stage, and the overlying alkalic flows to have erupted during either the caldera-filling or the post-caldera stage (see Kirkpatrick et al., 1980).

Ash layer characteristics and major element oxide content of glasses from ODP Sites 162-982 and 162-985

The Cenozoic volcanic activity on Iceland has been recorded in North Atlantic sediments drilled during several Ocean Drilling Program (ODP)/Deep Sea Drilling Project legs (Legs 104, 151, 152, 162, and 163). Leg 162 (North Atlantic-Arctic Gateways II) recovered ash layers at Sites 982, 985, and 907 (Jansen, Raymo, Blum, et al., 1996, doi:10.2973/odp.proc.ir.162.1996). The revisited Site 907 was first drilled during Leg 151, and the ash from this site has been described in detail by Lacasse et al. (1996, doi:10.2973/odp.proc.sr.151.122.1996) and Werner et al. (1996, doi:10.2973/odp.proc.sr.151.123.1996). Site 982 is located within the Hatton-Rockall Basin on the Rockall Plateau, which is situated west of the British Isles. Site 985 is located northeast of Iceland at the foot of the eastern slope of the Iceland Plateau, adjacent to the Norwegian Basin. Here we report chemical analyses of Neogene tephra layers from Holes 982A, 983B, 982C, 985A, and 985B. The sedimentary sequence at Site 982 spans the lower Miocene-Holocene; Site 985 recovered sediments spanning the upper Oligocene-Holocene. Twenty-two distinct ash layers and ash-bearing sediments were sampled in Holes 982A-982C (Cores 162-982A-16H through 24H, 162-982B-14H through 56X, and 162-982C-15H through 27H), and 59 ash layers were sampled in Holes 985A and 985B (Cores 162-985A-11H through 59X, and 162-985B-11H through 14H). Almost 50% of the sampled ash is strongly altered (predominantly from Site 985). A cluster of altered thin layers in the lower Pliocene of Site 985 (top of Unit III) is remarkable.

Rhenium and platinum-group element concentrations and Os isotope data for Philippine Sea Plate basalts from DSDP Legs 31, 58, and 59

Platinum-group elements (PGE), rhenium and osmium isotope data are reported for basalts from Deep Sea Drilling Project cores in the Philippine Sea Plate (PSP). Lithophile trace element and isotopic characteristics indicate a range of source components including DMM, EMII and subduction-enriched mantle. MORB-like basalts possess smooth, inclined chondrite-normalised PGE patterns with high palladium-PGE/iridium-PGE ratios, consistent with previously published data for MORB, and with the inferred compatibility of PGE. In contrast, while basalts with EMII-type lithophile element chemistry possess high Pt/Ir ratios, many have much lower Pd/Ir and unusually high Ru/Ir of >10. Similarly, back-arc samples from the Shikoku and Parece-Vela basins have very high Ru/Ir ratios (>30) and Pd/Ir as low as 1.1. Such extreme Pd/Ir and Ru/Ir ratios have not been previously reported in mafic volcanic suites and cannot be easily explained by variable degrees of melting, fractional crystallisation or by a shallow-level process such as alteration or degassing. The data appear most consistent with sampling of at least two mantle components with distinct PGE compositions. Peridotites with the required PGE characteristics (i.e. low Pd, but relatively high Ru and Re) have not been documented in oceanic mantle, but have been found in sub-continental mantle lithosphere and are the result of considerable melt depletion and selective metasomatic enrichment (mainly Re). The long-term presence of subduction zones surrounding the Philippine Sea Plate makes this a prime location for metasomatic enrichment of mantle, either through fluid enrichment or infiltration by small melt fractions. The Re-Os isotope data are difficult to interpret with confidence due to low Os concentrations in most samples and the uncertainty in sample age. Data for Site 444A (Shikoku Basin) give an age of 17.7+/-1.3 Ma (MSWD = 14), consistent with the proposed age of basement at the site and thus provides the first robust radiometric age for these samples. The initial 187Os/188Os of 0.1298+/-0.0069 is consistent with global MORB, and precludes significant metasomatic enrichment of Os by radiogenic slab fluids. Re-Os data for Sites 446A (two suites, Daito Basin) and 450 (Parece-Vela Basin) indicate ages of 73, 68 and 43 Ma, which are respectively, 30, 17 and >12 Ma older than previously proposed ages. The alkalic and tholeiitic suites from Site 446A define regression lines with different 187Os/188Osinitial (0.170+/-0.033 and 0.112+/-0.024, respectively) which could perhaps be explained by preferential sampling of interstitial, metasomatic sulphides (with higher time-integrated Re/Os ratios) by smaller percentage alkalic melts. One sample, with lithophile elements indistinguishable from MORB, is Os-rich (146 pg/g) and has an initial 187Os/188Os of 0.1594, which is at the upper limit of the accepted OIB range. Given the Os-rich nature of this sample and the lack of evidence for subduction or recycled crust inputs, this osmium isotope ratio likely reflects heterogeneity in the DMM. The dataset as a whole is a striking indication of the possible PGE and Os isotope variability within a region of mantle that has experienced a complex tectonic history.

Trace element contents in zircons from oceanic crust

We present newly acquired trace element compositions for more than 300 zircon grains in 36 gabbros formed at the slow-spreading Mid-Atlantic and Southwest Indian Ridges. Rare earth element patterns for zircon from modern oceanic crust completely overlap with those for zircon crystallized in continental granitoids. However, plots of U versus Yb and U/Yb versus Hf or Y discriminate zircons crystallized in oceanic crust from continental zircon, and provide a relatively robust method for distinguishing zircons from these environments. Approximately 80% of the modern ocean crust zircons are distinct from the field defined by more than 1700 continental zircons from Archean and Phanerozoic samples. These discrimination diagrams provide a new tool for fingerprinting ocean crust zircons derived from reservoirs like that of modern mid-ocean ridge basalt (MORB) in both modern and ancient detrital zircon populations. Hadean detrital zircons previously reported from the Acasta Gneiss, Canada, and the Narryer Gneiss terrane, Western Australia, plot in the continental granitoid field, supporting hypotheses that at least some Hadean detrital zircons crystallized in continental crust forming magmas and not from a reservoir like modern MORB.