Nd-, O-, and H-isotopic evidence for complex, closed-system fluid evolution of the peralkaline Ilimaussaq intrusion, South Greenland

Relatively homogeneous oxygen isotope compositions of amphibole, clinopyroxene, and olivine separates (+5.2 to +5.7parts per thousand relative to VSMOW) and neodymium isotope compositions (epsilon(Nd(T)) = -0.9 to -1.8 for primary magmatic minerals and epsilon(Nd(T)) = -0.1 and -0.5 for mineral separates from late-stage pegmatites and hydrothermal veins) from the alkaline to agpaitic llimaussaq intrusion, South Greenland, indicate a closed system evolution of this igneous complex and support a mantle derivation of the magma. In contrast to the homogeneous oxygen and neodymium isotopic data, deltaD values for hand-picked amphibole separates vary between -92 and -232parts per thousand and are among the most deuterium-depleted values known from igneous amphiboles. The calculated fluid phase coexisting with these amphiboles has a homogeneous oxygen isotopic composition within the normal range of magmatic waters, but extremely heterogeneous and low D/H ratios, implying a decoupling of the oxygen- and hydrogen isotope systems. Of the several possibilities that can account for such unusually low deltaD values in amphiboles (e.g., late-stage hydrothermal exchange with meteoric water, extensive magmatic degassing, contamination with organic matter, and/or effects of Fe-content and pressure on amphibole-water fractionation) the most likely explanation for the range in deltaD values is that the amphiboles have been influenced by secondary interaction and reequilibration with D-depleted fluids obtained through late-magmatic oxidation of internally generated CH(4) and/or H(2). This interpretation is consistent with the known occurrence of abundant magmatic CH(4) in the Ilimaussaq rocks and with previous studies on the isotopic compositions of the rocks and fluids. Copyright (C) 2004 Elsevier Ltd.

Nd and Sr isotope compositions in modern and fossil bones - Proxies for vertebrate provenance and taphonomy

Rare earth elements (REE), while not essential for the physiologic functions of animals, are ingested and incorporated in ppb concentrations in bones and teeth. Nd isotope compositions of modern bones of animals from isotopically distinct habitats demonstrate that the (143)Nd/(144)Nd of the apatite can be used as a fingerprint for bedrock geology or ambient water mass. This potentially allows the provenance and migration of extant vertebrates to be traced, similar to the use of Sr isotopes. Although REE may be enriched by up to 5 orders of magnitude during diagenesis and recrystallization of bone apatite, in vivo (143)Nd/(144)Nd may be preserved in the inner cortex of fossil bones or enamel. However, tracking the provenance of ancient or extinct vertebrates is possible only for well-preserved archeological and paleontological skeletal remains with in vivo-like Nd contents at the ppb-level. Intra-bone and -tooth REE analysis can be used to screen for appropriate areas. Large intra-bone Nd concentration gradients of 10(1)-10(3) are often measured. Nd concentrations in the inner bone cortex increase over timescales of millions of years, while bone rims may be enriched over millenial timescales. Nevertheless, epsilon(Nd) values are often similar within one epsilon(Nd) unit within a single bone. Larger intra-bone differences in specimens may either reflect a partial preservation of in vivo values or changing epsilon(Nd) values of the diagenetic fluid during fossilization. However, most fossil specimens and the outer rims of bones will record taphonomic (143)Nd/(144)Nd incorporated post mortem during diagenesis. Unlike REE patterns, (143)Nd/(144)Nd are not biased by fractionation processes during REE-uptake into the apatite crystal lattice, hence the epsilon(Nd) value is an important tracer for taphonomy and reworking. Bones and teeth from autochthonous fossil assemblages have small variations of +/- 1 epsilon(Nd) unit only. In contrast, fossil bones and teeth from over 20 different marine and terrestrial fossil sites have a total range of epsilon(Nd) values from -13.0 to 4.9 (n = 80), often matching the composition of the embedding sediment. This implies that the surrounding sediment is the source of Nd in the fossil bones and that the specimens of this study seem not to have been reworked. Differences in epsilon(Nd) values between skeletal remains and embedding sediment may either indicate reworking of fossils and/or a REE-uptake from a diagenetic fluid with non-sediment derived epsilon(Nd) values. The latter often applies to fossil shark teeth, which may preserve paleo-seawater values. Complementary to epsilon(Nd) values, (87)Sr/(86)Sr can help to further constrain the fossil provenance and reworking. (C) 2011 Elsevier Ltd. All rights reserved.

Constraints on Miocene oceanography and climate in the Western and Central Paratethys: O-, Sr-, and Nd-isotope compositions of marine fish and mammal remains

The Paratethys evolved as a marginal sea during the Alpine-Himalayan orogeny in the Oligo-Miocene. Sediments from the northern Alpine Molasse Basin, the Vienna, and the Pannonian Basins located in the western and central part of the Paratethys thus provide unique information on regional changes in climate and oceanography during a period of active Alpine uplift Oxygen isotope compositions of well-preserved phosphatic fossils recovered from the sediments support deposition under sub-tropical to warm-temperate climate with water temperatures of 14 to 28 degrees C for the Miocene. delta(18)O values of fossil shark teeth are similar to those reported for other Miocene marine sections and, using the best available estimates of their biostratigraphic age, show a variation until the end of the Badenian similar to that reported for composite global record. The (87)Sr/(86)Sr isotope ratios of the fossils follow the global Miocene seawater trend, albeit with a much larger scatter. The deviations of (87)Sr/(86)Sr in the samples from the well-constrained seawater curve are interpreted as due to local input of terrestrially-derived Sr. Contribution of local sources is also reflected in the epsilon(Nd) values, consistent with input from ancient crystalline rocks (e.g., Bohemian Massif and/or Mesozoic sediments with epsilon(Nd) < -9. On the other hand, there is evidence for input from areas with Neogene volcanism as suggested by samples with elevated epsilon(Nd) values >-7. Excluding samples showing local influence on the water column, an average epsilon(Nd) value of -7.9 +/- 0.5 may be inferred for the Miocene Paratethys. This value is indistinguishable from the epsilon(Nd) value of the contemporaneous Indian Ocean, supporting a dominant role of this ocean in the Western and Central Paratethys. (C) 2008 Elsevier B.V. All rights reserved.

Sr-Nd isotopic evidence for crustal contamination in the Niquelandia complex, Goias, Central Brazil

The Niquelandia complex is a Neoproterozoic mafic-ultramafic intrusion resulting from fractional crystallization of primary picritic basalt intrusions. It consists of two layered sequences: a lower and larger one (LS), where four stratigraphic units exhibit an upward decrease of ultramafic layers and increase of gabbroic layers; an upper, smaller sequence (US), separated from LS by a high-temperature shear zone and consisting of two stratigraphic units (gabbros + anorthosites and amphibolites). Nd and Sr isotopic analyses and rare earth element (REE) profiles provide evidence that the complex suffered important crustal contamination. The LS isotopic array trends from a DM region with positive epsilon Nd and moderately positive epsilon Sr towards a field occupied by crustal xenoliths, especially abundant in the upper LS (negative epsilon Nd and large, positive E:Sr). Each LS stratigraphic unit is distinct from the next underlying unit, showing lower epsilon Nd and higher epsilon Sr, suggesting inputs of fresh magma and mixing with the contaminated, residual magma. The US is characterised by a relatively high variation of epsilon Nd and constant epsilon Sr. REE patterns vary within each unit from LREE depleted to LREE enriched in the samples having lower epsilon Nd and higher epsilon Sr. The contamination process has been modelled by using the EC-AFC algorithms from [Spera, F.J., Bohrson, W.A., 2001. Energy-constrained open-system magmatic processes 1: general model and energy-constrained assimilation and fractional crystallization (EC-AFC) formulation. J. Petrology 42, 999-1018]. The differences between the LS and US isotopic arrays are consistent with contamination by the same crustal component, provided that its melting degree was higher in LS than in US. The different degrees of anatexis are explained by the heat budget released from the magma, higher in LS (because of its larger mass) than in US. Comparison of the correlations between isotopes and incompatible trace element ratios of the models and of the gabbros shows some differences, which are demonstrably related with the variable amount of cumulus phases and trapped melt in the gabbros. (c) 2007 Elsevier Ltd. All rights reserved.

Trace element geochemistry and Sr-Nd characteristics of Mesoproterozoic mafic intrusive rocks from Rondonia, Brazil, SW Amazonian Craton: Petrogenetic and tectonic inferences

The Amazonian Craton comprises an Archean domain surrounded by four successively younger Proterozoic tectonic provinces. Within the Rio-Negro-Juruena province the Serra da Providencia Intrusive Suite (1.60 and 1.53 Ga) consists of A-type rapakivi granites, charnockites and mangerites genetically associated with diabase dikes, gabbros and amphibolites lites. The original mafic melts were derived from a depleted mantle source (epsilon(Nd(T)) + 2.5 to +2.8; epsilon(Sr(T)) - 12.1). Underplated mafic magma induced melting of a short-lived fielsic crust, thus originating coeval felsic-inafic magmatism in a continental intraplate setting. The Colorado Complex, assigned to the Rondonian-San Ignacio province, comprises 1.35-1.36 Ga intrusive bimodal magmatism represented by monzonite gneisses associated with amphibolite, gabbro and metadiabase dikes intercalated with metasediments with detrital zircon that yield U-Pb ages of 1.35 to 1.42 Ga. Mafic samples display juvenile signatures (epsilon(Nd(T)) 0.0 to +5.2; epsilon(Sr(T)) -5.0 to -30.7) and are less contaminated than the Serra da Previdencia and Nova Brasiladndia ones. The generation of the basaltic magma is related to the subduction of an oceanic slab below the peridotite wedge (intraoceanic arc setting). Fluids and/or small melts from the slab impregnated the mantle. The Nova Brasilandia Sequence (Sunsas-Aguapei province) comprises a metasedimentary sequence intruded by 1.10-1.02 Ga metadiabases, gabbros, meta-gabbros, and amphibolites associated with granitic plutons (bimodal magmatism). The original tholeiitic magmas, derived from a depleted source (epsilon(Nd(T)) = +3.1 to +5.0), in a proto-oceanic setting, underwent subsequent contamination by the host rocks, as indicated by the isotopic and trace element data.

Nd and Pb isotope signatures on the Southeastern South American upper margin: Implications for sediment transport and source rocks

Neodymium and lead isotope values in sediment samples were used to interpret sediment transport and source rocks on the Southeastern South American upper margin. The sediments of the Argentinian margin exhibit an average epsilon(Nd) value of -1.9, indicating the influence of the Andean rocks as sediment sources. Sediments from the Rio de La Plata estuary show an average epsilon(Nd) value of -9.6 which is similar to that of the Southern Brazilian Upper Margin. Finally, sediments of Southeastern Brazil, which are associated with the transport of the Brazil Current exhibit an average epsilon(Nd) of -13.0. The Pb isotope signatures also confirm the differentiation of source rocks in the sedimentation of the study area. In addition, Pb isotopes helped to establish the extent of the influence of the Rio de La Plata on the sedimentation of the Southern Brazilian margin. In terms of Pb isotopes the sediments from the Rio de La Plata estuary and Southern Brazil are more radiogenic than those of Southeastem Brazil and the Argentinian margin. (c) 2007 Elsevier B.V. All rights reserved.

Comparação estratigráfica e isotópica (Sm-Nd) das regiões do Cerro Morado e La Torre, bacia de Ischigualasto-Villa Unión, noroeste argentino

A bacia de Ischigualasto-Villa Unión (Stipanicic e Bonaparte, 1972), que consiste de uma bacia alongada no sentido NW-SE e preenchida por rochas sedimentares Triássicas, foi gerada durante o Eo-Triássico pela reativação de descontinuidades crustais pré-existentes, associadas ao Gondwanides, na região central andina. Esta limitada a oeste, pelo vale do Rio Bermejo, no qual se encontra a zona de sutura do Vale Fértil, que se caracteriza pela borda ativa da bacia. Aqui, apresentar-se-á uma revisão da estratigrafia da bacia de Ischigualasto-Villa Unión, detalhando-se duas áreas chaves do registro estratigráfico inferior da bacia, na região do Cerro Morado e La Torre. Nestas áreas utilizou-se de técnicas estratigráficas e isotópicas (Sm-Nd em rocha total), com o objetivo de estabelecer uma correlação entre seus registros estratigráficos e suas áreas-fonte, além de caracterizar o comportamento tectônico. Ambos os perfis apresentam na porção inferior o preenchimento por um sistema aluvial-fluvial (Formação Tarjados sensu Romer & Jensen, 1966), subdividido em 4 associações de fácies: leques aluviais, canais fluviais, barras longitudinais e planície de inundação; e, na porção superior, o preenchimento por um sistema lacustre-deltáico (Formação Los Chañares sensu Frenguelli, 1944, no Cerro Morado e, Formação Lomas Blancas sensu Bossi e Herbst, 1968, em La Torre), subdividido em 4 associações de fácies: lacustre, pró delta, frente deltáica e planicie deltáica. Sobre as últimas unidades ocorre um lacólito de diabásio (Formação Los Baldecitos) A formação Tarjados apresenta as maiores idades TDM na região do Cerro Morado, variando entre 1,33 e 1,60 Ga, enquanto que em La Torre os valores oscilam entre 1,03 e 1,30 Ga. Isso é interpretado aqui como uma variação da área fonte entre as duas regiões na época de sua sedimentação, cujas regiões poderiam estar compartimentadas em termos tectônicos a época de sua deposição. A formação Los Chañares possui idades TDM variando entre 0,90 e 1,36 Ga, e a formação Lomas Blancas possui idades TDM entre 0,99 e 1,35 Ga, com dominância de valores próximos a 1 Ga nos dois blocos. Os diabásios da formação Los Baldecitos possuem idades TDM entre 0,51 e 0,33 Ga, cujo dado é interpretado com uma extração mantélica próxima a 130 Ma de um manto modificado pela orogenia do Ciclo Gondwânico, ou até mesmo pelo processo de slab break-off da crosta oceânica consumida pela colisão do micro-continente Cuyania com a proto-margem do Gondwana (Ramos, 1999). Todas as rochas dessa formação obtiveram valores de Epsilon Nd positivos. A integração dos dados estratigráficos e isotópicos nos permite interpretar as duas áreas como meio-grábens distintos na etapa inicial de preenchimento, em função da discrepância entre as idades TDM entre as áreas. A partir da transgressão lacustre as áreas-fonte são uniformizadas e o alto estrutural pré-existente é transgredido.

The effect of hydrothermal alteration on the Sr and Nd isotopic signatures of the Barra do Itapirapua Carbonatite, Southern Brazil

The Cretaceous Barra do Itapirapua carbonatite in southern Brazil experienced extensive postmagmatic hydrothermal alteration. In this article, Sr and Nd isotope ratios of coexisting samples of hydrothermally overprinted and of preserved, nonoverprinted carbonatite are presented. Hydrothermal alteration caused strong REE enrichment, leading to the formation of minerals of the bastnaesite group. In the overprinted carbonatite, Nd contents reach 4000 ppm, two orders of magnitude higher than in the fresh carbonatite, but epsilon(Nd) varies only within a range of 3.4 units. In contrast, Sr was leached from the carbonatite during the postmagmatic alteration; hence values of around 10,000 ppm in the fresh carbonatite drop to about 1000 ppm in the overprinted samples. Leaching is accompanied by a variation of Sr isotopic composition toward more radiogenic values, resulting in an increase of 15 units in epsilon(Sr). Variation of Sr isotopic composition is related to postmagmatic alteration and is decoupled from the variation of Nd isotopic composition, ruling out heterogeneities in the mantle source as the main cause of isotopic variability in the data set. Furthermore, this cannot be explained by bulk crustal contamination. A two-step model is proposed in which (1) a REE-rich, carbonatite-derived hydrothermal fluid overprinted the pristine carbonatite, causing REE-enrichment with a relative small change of isotopic composition; and (2) crust-derived hydrothermal fluids percolated the cooling carbonatite, leaching the original Sr from the carbonatite and introducing a more radiogenic Sr isotopic signature. The amounts of carbonatite-derived Nd with primitive, carbonatite-like Nd isotope ratios introduced during the first stage of hydrothermal alteration are high enough to buffer the effect of crust-derived Nd on the Nd isotopic composition of the overprinted carbonatite.

High-pressure granulites from Carire, Borborema Province, NE Brazil: Tectonic setting, metamorphic conditions and U-Pb, Lu-Hf and Sm-Nd geochronology

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Geological evolution of the basement rocks in the east-central part of the Rondonia Tin Province, SW Amazonian craton, Brazil: U-Pb and Sm-Nd isotopic constraints

On the basis of geologic, petrologic, and U-Pb geochronologic data the basement rocks in the east-central part of the Rondonia Tin Province (RTP, southwestern Amazonian craton) are grouped into five lithologic associations: (1) tonalitic gneiss (1.75 Ga); (2) enderbitic granulite (1.73 Ga); (3) paragneiss; (4) granitic and charnockitic augen gneisses (1.57-1.53 Ga); and (5) fine-grained granitic gneiss and charnockitic granulite (1.43-1.42 Ga). The first three are related to development of the Paleoproterozoic Rio Negro-Juruena Province and represent the oldest crust in the region. The tonalitic gneisses and enderbitic granulites show calc-alkaline affinities and Nd isotopic compositions (initial epsilon(Nd) = +0-1 to -1.5; T-DM of 2.2-2.1 Ga) that suggest a continental arc margin setting for the original magmas. The paragneisses yield T-DM values of 2.2-2.1 Ga suggesting that source material was primarily derived from the Ventuari-Tapajos and Rio Negro-Juruena crusts, but detrital zircon ages and an intrusive granitoid bracket deposition between 1.67 and 1.57 Ga. The granitic and charnockitic augen gneisses show predominantly A-type and within-plate granite affinities, but also some volcanic arc granite characteristics. The initial epsilon(Nd) values (+0.6 to +2.0) indicate mixing of magmas derived from depleted mantle and older crustal sources. These rocks are correlated to the 1.60-1.53 Ga Serra da Providencia intrusive suite that reflects inboard magmatism coeval with the Cachoeirinha orogen located to the southeast. The fine-grained granitic gneiss and charnockitic granulites represent the first record of widespread magmatism at 1.43-1.42 Ga in northern Rondonia. Their geochemical signatures and the slightly positive initial epsilon(Nd) values (+0.7 to +1.2) are very similar to those of the most evolved granites of the calc-alkaline Santa Helena batholith farther southeast. U-Pb monazite and Sm-Nd whole-rock-garnet ages demonstrate that a high-grade tectonometa-morphic episode occurred in this region at 1.33-1.30 Ga. This episode attained upper-amphibolite conditions and is interpreted as the peak of the Rondonian-San Ignacio orogeny. The U-Pb and Sm-Nd data presented here and data published on rapakivi granites elsewhere indicate that the east-central part of the RTP is a poly-orogenic region characterized by successive episodes of magmatism, metamorphism, and deformation between 1.75 and 0.97 Ga. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Contrasted crustal sources as defined by whole-rock and Sr-Nd-Pb isotope geochemistry of neoproterozoic early post-collisional granitic magmatism within the Southern Brazilian Shear Belt, Camboriu, Brazil

The early phase of post-collisional granitic magmatism in the Camboriu region, south Brazil, is represented by the porphyritic biotite +/- hornblende Rio Pequeno Granite (RPG; 630-620 Ma) and the younger (similar to 610 Ma), equigranular, biotite +/- muscovite Serra dos Macacos Granite (SMG). The two granite types share some geochemical characteristics, but the more felsic SMG constitutes a distinctive group not related to RPG by simple fractionation processes, as indicated by its lower FeOt, TiO2, K2O/Na2O and higher Zr Al2O3, Na2O, Ba and Sr when compared to RPG of similar SiO2 range. Sr-Nd-Pb isotopes require different sources. The SMG derives from old crustal sources, possibly related to the Paleoproterozoic protoliths of the Camboriu Complex, as indicated by strongly negative epsilon Nd-t (-23 to -24) and unradiogenic Pb (e.g., Pb-206/Pb-204 = 16.0-16.3; Pb-207/Pb-204 = 15.3-15.4) and confirmed by previous LA-MC-ICPMS data showing dominant zircon inheritance of Archean to Paleoproterozoic age. In contrast, the RPG shows less negative epsilon Nd-t (-12 to -15) and a distinctive zircon inheritance pattern with no traces of post-1.6 Ga sources. This is indicative of younger sources whose significance in the regional context is still unclear; some contribution of mantle-derived magmas is indicated by coeval mafic dykes and may account for some of the geochemical and isotopic characteristics of the least differentiated varieties of the RPG. The transcurrent tectonics seems to have played an essential role in the generation of mantle-derived magmas despite their emplacement within a low-strain zone. It may have facilitated their interaction with crustal melts which seem to be to a large extent the products of reworking of Paleoproterozoic orthogneisses from the Camboriu Complex. (C) 2012 Elsevier Ltd. All rights reserved.

Petrology and Sr-Nd characteristics of the Nova Lacerda dike swarm, SW Amazonian Craton: new insights regarding its subcontinental mantle source and Mesoproterozoic geodynamics

The NNW-trending Nova Lacerda tholeiitic dike swarm in Mato Grosso State, Central Brazil, intrudes the Nova Lacerda granite (1.46 Ga) and the Jauru granite-greenstone terrain (ca. 1.79-1.77 Ga). The swarm comprises diabases I and II and amphibolites emplaced at ca. 1.38 Ga. Geochemical data indicate that these are evolved tholeiites characterized by high LILE/HSFE and LREE/HSFE ratios. Isotopic modelling yields positive epsilon(Nd)(T) values (+0.86 to +2.65), whereas values for epsilon(Sr)(T) range from positive to negative (+1.96 to -5.56). Crustal contamination did not play a significant petrogenetic role, as indicated by a comparison of isotopic data (Sr-Nd) from both dikes and country rocks, and by the relationship between isotopic and geochemical parameters (SiO2, K2O, Rb/Sr, and La/Yb) of the dikes. We attribute the origin of these tholeiites to fractional crystallization of evolved melts derived from a heterogeneous mantle source. Comparison of the geochemical and isotopic data of the studied swarm and other tholeiitic Mesoproterozoic mafic intrusions of the SWAmazonian Craton the Serra da Providencia, Colorado, and Nova Brasilandia bimodal suites - indicates that parental melts of the Nova Lacerda swarm were derived from the most enriched mantle source. This enrichment was probably caused by the stronger influence of the EMI component on the DMM end-member. These data, coupled with trace element bulk-rock geochemistry of the country rocks, and comparisons with the Colorado Complex of similar age, suggest a continental-margin arc setting for the emplacement of the Nova Lacerda dikes.

(Table 1) Nd and Sr isotopic analyses of ODP Leg 127 basalts

Sr and Nd isotopic compositions are reported for basaltic rocks collected during ODP Leg 127 from the Yamato Basin, a rifted backarc basin in the Japan Sea. The basalts are classified into two groups in terms of Nd isotopic composition: the upper sills at Site 797 are characterized by higher 143Nd/144Nd ratios (0.513083-0.513158, epsilon-Nd = 8.68-10.14) and the basalts from Site 794 and the lower sills at Site 797 have lower 143Nd/144Nd ratios (0.512684-0.512862, epsilon-Nd = 0.90-4.37). All of the basalts show higher Sr isotopic compositions than those of the mantle array, which is attributed to seawater alteration. The basalts with lower Nd isotopic values ranging in age from 20.6 to 17.3 Ma have tapped an enriched subcontinental upper mantle (SCUM) with the minor involvement of a depleted asthenospheric mantle (AM). Subsequent change in composition through the physical replacement of SCUM by AM yielded the basalts of the upper sills of higher Nd isotopic compositions. This event within the upper mantle was associated with the breakup of the overlying lithosphere during the rifting of the Japan Sea backarc basin.

Nd and Sr isotopic composition of Cretaceous MORB

Chemical and isotopic (Nd and Sr) compositions have been determined for 12 Cretaceous basaltic samples (108 Ma old) from Holes 417D and 418A of Legs 51,52 and 53. We have found that: (1) The chemical compositions are typical of MORB. They do not vary systematically with the stratigraphic positions of the analyzed samples; thus, the chemical evolution is independent of the eruption sequence that occurred at this Cretaceous ridge. (2) REE patterns for all rocks are characterized by a strong LREE depletion with (La/Sm)N = 0.38-0.50; no significant Eu anomalies are found; HREE are nearly flat or slightly depleted towards Yb-Lu and have 12-18 * chondritic abundances. Combining the results of previous studies, it suggests that no significant temporal and spatial variation in magma chemistry (especially for LIL elements) has occurred in the 'normal' ridge segments over the last 150 Ma. (3) lsotopically, 143Nd/144Nd ratios vary from 0.513026 to 0.513154, corresponding to epsilon-Nd(0) = +7.5 to +10, and they fall in the typical range of MORB. However, these rocks have unexpectedly high 87Sr/86Sr ratios (0.70355-0.70470) which are attributed to the result of seawater-rock interaction. (4) The Nd model ages (Tin), ranging from 1.53 to 2.47 (average 2.06) AE, suggest that the upper mantle source(s) underwent a large scale chemical differentiation leading to LREE and other LIL element depletion about 2 AE ago, assuming a simple two-stage model. More realistically, the variation in Tm(Nd) or epsilon-Nd could be derived from mixing of heterogeneous mantle sources that were a consequence of continuous mantle differentiation and continental formation. (5) Because of the low mg values (0.52-0.63), the analyzed basaltic rocks do not represent primary liquids of mantle melting. The variation in La/Sm ratios and TiO2 are not compatible with a model in which all rocks are genetically related by a simple fractional crystallization. Rather, it is proposed that the basaltic rocks might have been derived from some heterogeneous upper mantle source with or without later magmatic mixing, and followed by some shallow-level fraetionations.