Convective patterns under the Indo-Atlantic << box >>

Using fluid mechanics, we reinterpret the mantle images obtained from global and regional tomography together with geochemical, geological and paleomagnetic observations, and attempt to unravel the pattern of convection in the Indo-Atlantic "box" and its temporal evolution over the last 260 Myr. The << box >> presently contains a) a broad slow seismic anomaly at the CMB which has a shape similar to Pangea 250 Myr ago, and which divides into several branches higher in the lower mantle, b) a "superswell, centered on the western edge of South Africa, c) at least 6 "primary hotspots" with long tracks related to traps, and d) numerous smaller hotspots. In the last 260 Myr, this mantle box has undergone 10 trap events, 7 of them related to continental breakup. Several of these past events are spatially correlated with present-day seismic anomalies and/or upwellings. Laboratory experiments show that superswells, long-lived hotspot tracks and traps may represent three evolutionary stages of the same phenomenon, i.e. episodic destabilization of a hot, chemically heterogeneous thermal boundary layer, close to the bottom of the mantle. When scaled to the Earth's mantle, its recurrence time is on the order of 100-200 Myr. At any given time, the Indo-Atlantic box should contain 3 to 9 of these instabilities at different stages of their development, in agreement with observations. The return flow of the downwelling slabs, although confined to two main << boxes >> (Indo-Atlantic and Pacific) by subduction zone geometry, may therefore not be passive, but rather take the form of active thermochemical instabilities. (c) 2005 Elsevier B.V. All rights reserved.

Constrains on the structure of Maio Island (Cape Verde) by a three-dimensional gravity model: Imaging partially exhumed magma chambers

We propose a 3-D gravity model for the volcanic structure of the island of Maio (Cape Verde archipelago) with the objective of solving some open questions concerning the geometry and depth of the intrusive Central Igneous Complex. A gravity survey was made covering almost the entire surface of the island. The gravity data was inverted through a non-linear 3-D approach which provided a model constructed in a random growth process. The residual Bouguer gravity field shows a single positive anomaly presenting an elliptic shape with a NWSE trending long axis. This Bouguer gravity anomaly is slightly off-centred with the island but its outline is concordant with the surface exposure of the Central Igneous Complex. The gravimetric modelling shows a high-density volume whose centre of mass is about 4500 m deep. With increasing depth, and despite the restricted gravimetric resolution, the horizontal sections of the model suggest the presence of two distinct bodies, whose relative position accounts for the elongated shape of the high positive Bouguer gravity anomaly. These bodies are interpreted as magma chambers whose coeval volcanic counterparts are no longer preserved. The orientation defined by the two bodies is similar to that of other structures known in the southern group of the Cape Verde islands, thus suggesting a possible structural control constraining the location of the plutonic intrusions.

Deep crustal structure of the Parana Basin from receiver functions and Rayleigh-wave dispersion: Evidence for a fragmented cratonic root

The deep crustal structure of the Parana Basin of southern Brazil is investigated by analyzing P- and PP-wave receiver functions at 17 Brazilian Lithosphere Seismic Project stations within the basin. The study area can be described as a typical Paleozoic intracratonic basin that hosts one of the largest Large Igneous Province of the world and makes a unique setting for investigating models of basin subsidence and their interaction with mantle plumes. Our study consists of (1) an analysis of the Moho interaction phases in the receiver functions to obtain the thickness and bulk Vp/Vs ratio of the basin`s underlying crust and (2) a joint inversion with Rayleigh-wave dispersion velocities from an independent tomographic study to delineate the detailed S-wave velocity variation with depth. The results of our analysis reveal that Moho depths and bulk Vp/Vs ratios (including sediments) vary between 41 and 48 km and between 1.70 and 1.76, respectively, with the largest values roughly coinciding with the basin`s axis, and that S-wave velocities in the lower crust are generally below 3.8 km/s. Select sites within the basin, however, show lower crustal S-wave velocities slightly above 3.9 km/s suggestive of underplated mafic material. We show that these observations are consistent with a fragmented cratonic root under the Parana basin that defined a zone of weakness for the initial Paleozoic subsidence of the basin and which allowed localized mafic underplating of the crust along the suture zones by Cenozoic magmatism.

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.

Geochemistry and oxygen and iron isotope ratios of mafic rocks

Recycling of oceanic crust into the deep mantle via subduction is a widely accepted mechanism for creating compositional heterogeneity in the upper mantle and for explaining the distinct geochemistry of mantle plumes. The oxygen isotope ratios (d18O) of some ocean island basalts (OIB) span values both above and below that of unmetasomatised upper mantle (5.5 ± 0.4 per mil) and provide support for this hypothesis, as it is widely assumed that most variations in d18O are produced by near-surface low-temperature processes. Here we show a significant linear relationship between d18O and stable iron isotope ratios (d57Fe) in a suite of pristine eclogite xenoliths. The d18O values of both bulk samples and garnets range from values within error of normal mantle to significantly lighter values. The observed range and correlation between d18O and d57Fe is unlikely to be inherited from oceanic crust, as d57Fe values determined for samples of hydrothermally altered oceanic crust do not differ significantly from the mantle value and show no correlation with d18O. It is proposed that the correlated d57Fe and d18O variations in this particular eclogite suite are predominantly related to isotopic fractionation by disequilibrium partial melting although modification by melt percolation processes cannot be ruled out. Fractionation of Fe and O isotopes by removal of partial melt enriched in isotopically heavy Fe and O is supported by negative correlations between bulk sample d57Fe and Cr content and bulk sample and garnet d18O and Sc contents, as Cr and Sc are elements that become enriched in garnet- and pyroxene-bearing melt residues. Melt extraction could take place either during subduction, where the eclogites represent the residues of melted oceanic lithosphere, or could take place during long-term residence within the lithospheric mantle, in which case the protoliths of the eclogites could be of either crustal or mantle origin. This modification of both d57Fe and d18O by melting processes and specifically the production of low-d18O signatures in mafic rocks implies that some of the isotopically light d18O values observed in OIB and eclogite xenoliths may not necessarily reflect near-surface processes or components.

Isotope ratios and trace element concentrations of basalts from DSDP Leg 15 sites (Appendix 1)

Formation of the Cretaceous Caribbean plateau, including the komatiites of Gorgona, has been linked to the currently active Galápagos hotspot. We use Hf-Nd isotopes and trace element data to characterise both the Caribbean plateau and the Galápagos hotspot, and to investigate the relationship between them. Four geochemical components are identified in the Galápagos mantle plume: two 'enriched' components with epsilon-Hf and epsilon-Nd similar to enriched components observed in other mantle plumes, one moderately enriched component with high Nb/Y, and a fourth component which most likely represents depleted MORB source mantle. The Caribbean plateau basalt data form a linear array in Hf-Nd isotope space, consistent with mixing between two mantle components. Combined Hf-Nd-Pb-Sr-He isotope and trace element data from this study and the literature suggest that the more enriched Caribbean end member corresponds to one or both of the enriched components identified on Galápagos. Likewise, the depleted end member of the array is geochemically indistinguishable from MORB and corresponds to the depleted component of the Galápagos system. Enriched basalts from Gorgona partially overlap with the Caribbean plateau array in epsilon-Hf vs. epsilon-Nd, whereas depleted basalts, picrites and komatiites from Gorgona have a high epsilon-Hf for a given epsilon-Nd, defining a high-epsilon-Hf depleted end member that is not observed elsewhere within the Caribbean plateau sequences. This component is similar, however, in terms of Hf-Nd-Pb-He isotopes and trace elements to the depleted plume component recognised in basalts from Iceland and along the Reykjanes Ridge. We suggest that the Caribbean plateau represents the initial outpourings of the ancestral Galápagos plume. Absence of a moderately enriched, high Nb/Y component in the older Caribbean plateau (but found today on the island of Floreana) is either due to changing source compositions of the plume over its 90 Ma history, or is an artifact of limited sampling. The high-epsilon-Hf depleted component sampled by the Gorgona komatiites and depleted basalts is unique to Gorgona and is not found in the Caribbean plateau. This may be an indication of the scale of heterogeneity of the Caribbean plateau system; alternatively Gorgona may represent a separate oceanic plateau derived from a completely different Pacific plume, such as the Sala y Gomez.

(Table 1) Bulk geochemistry of some lavas at DSDP Hole 55-433C

According to Wilson's (1963a, b) hypothesis, the volcanoes of the Hawaiian-Emperor Chain are formed as the Pacific lithospheric plate moves over a source of magma in the mantle. Morgan (1971, 1972) proposed that these "hot spots" resulted from "mantle plumes" that rise vertically from the core/mantle boundary and that are fixed about the deep mantle and rotating globe poles. The age of volcanoes increases with distance away from the recent "hot spot" beneath Kilauea volcano. The Hawaiian-Emperor bend indicates that the direction of motion of the Pacific plate changed about 40 m.y. ago.

Chemical composition of basalts and basaltic glasses from the Sierra Leone Fracture Zone region

Major and trace element (including REE) geochemistry of basalts and chilled basaltic glasses from the MAR axial zone in the vicinity of the Sierra Leone FZ (5-7°10'N) has been studied. Associations of basalts of various compositions with particular ocean-floor geological structural features have been analyzed as well. Three basaltic varieties have been discriminated. Almost ubiquitous are high-Mg basalts (Variety 1) that are derivatives of N-MORB tholeiitic melts and that are produced in the axial zone of spreading. Variety 2 is alkaline basalts widespread on the southwestern flank of the MAR crest zone in the Sierra Leone region, likely generated through deep mantle melting under plume impact. Variety 3 is basalts derivative from T- and P-MORB-like tholeiitic melts and originating through addition of deeper mantle material to depleted upper mantle melts. Magma generation parameters, as calculated from chilled glass compositions, are different for depleted tholeiites (44-55 km, 1320-1370°C) and enriched tholeiites (45-78 km, 1330-1450°C). Mantle plume impact is shown to affect not only tholeiitic basalt compositions but also magma generation conditions in the axial spreading zone, resulting in higher Ti and Na concentrations in melts parental to rift-related basalts occurring near the plume. T- and P-MORBs are also developed near areas where mantle plumes are localized. High-Mg basalts are shown to come in several types with distinctive Ti and Na contents. Nearly every single MAR segment (bounded by sinistral strike slips and the Bogdanov Fracture Zone) is featured by its own basalt type suggesting that it has formed above an asthenospheric diapir with its unique magma generation conditions. These conditions are time variable. Likely causes of temporal and spatial instability of the mantle upwelling beneath this portion of the MAR are singular tectonic processes and plume activity. In sulfide-bearing rift morphostructures (so-called "Ore area'' and the Markov Basin), basalts make up highly evolved suites generated through olivine and plagioclase fractionation, which is suggestive of relatively long-lived magma chambers beneath the sulfide-bearing rift morphostructures. Functioning of these chambers is a combined effect of singular geodynamic regime and plume activity. In these chambers melts undergo deep differentiation leading to progressively increasing concentration of sulfide phase, eventually to be supplied to the hydrothermal plumbing system.

Major and trace element content of the 17 lithostratigraphic units recovered during DSDP Leg 49 (Table 1)

IPOD Leg 49 recovered basalts from 9 holes at 7 sites along 3 transects across the Mid-Atlantic Ridge: 63°N (Reykjanes), 45°N and 36°N (FAMOUS area). This has provided further information on the nature of mantle heterogeneity in the North Atlantic by enabling studies to be made of the variation of basalt composition with depth and with time near critical areas (Iceland and the Azores) where deep mantle plumes are thought to exist. Over 150 samples have been analysed for up to 40 major and trace elements and the results used to place constraints on the petrogenesis of the erupted basalts and hence on the geochemical nature of their source regions. It is apparent that few of the recovered basalts have the geochemical characteristics of typical "depleted" midocean ridge basalts (MORB). An unusually wide range of basalt compositions may be erupted at a single site: the range of rare earth patterns within the short section cored at Site 413, for instance, encompasses the total variation of REE patterns previously reported from the FAMOUS area. Nevertheless it is possible to account for most of the compositional variation at a single site by partial melting processes (including dynamic melting) and fractional crystallization. Partial melting mechanisms seem to be the dominant processes relating basalt compositions, particularly at 36°N and 45°N, suggesting that long-lived sub-axial magma chambers may not be a consistent feature of the slow-spreading Mid-Atlantic Ridge. Comparisons of basalts erupted at the same ridge segment for periods of the order of 35 m.y. (now lying along the same mantle flow line) do show some significant inter-site differences in Rb/Sr, Ce/Yb, 87Sr/86Sr, etc., which cannot be accounted for by fractionation mechanisms and which must reflect heterogeneities in the mantle source. However when hygromagmatophile (HYG) trace element levels and ratios are considered, it is the constancy or consistency of these HYG ratios which is the more remarkable, implying that the mantle source feeding a particular ridge segment was uniform with respect to these elements for periods of the order of 35 m.y. and probably since the opening of the Atlantic. Yet these HYG element ratios at 63°N are very different from those at 45°N and 36°N and significantly different from the values at 22°N and in "MORB". The observed variations are difficult to reconcile with current concepts of mantle plumes and binary mixing models. The mantle is certainly heterogeneous, but there is not simply an "enriched" and a "depleted" source, but rather a range of sources heterogeneous on different scales for different elements - to an extent and volume depending on previous depletion/enrichment events. HYG element ratios offer the best method of defining compositionally different mantle segments since they are little modified by the fractionation processes associated with basalt generation.

Element analyses and isotope composition of basalts from the West Valley segment, Juan de Fuca Ridge, northeast Pacific

The 50 km-long West Valley segment of the northern Juan de Fuca Ridge is a young, extension-dominated spreading centre, with volcanic activity concentrated in its southern half. A suite of basalts dredged from the West Valley floor, the adjacent Heck Seamount chain, and a small near-axis cone here named Southwest Seamount, includes a spectrum of geochemical compositions ranging from highly depleted normal (N-) MORB to enriched (E-) MORB. Heck Seamount lavas have chondrite-normalized La/Sm en -0.3, 87Sr/86Sr = 0.70235 - 0.70242, and 206Pb/204Pb = 18.22 - 18.44, requiring a source which is highly depleted in trace elements both at the time of melt generation and over geologic time. The E-MORB from Southwest Seamount have La/Sm en -1.8, 87Sr/86Sr = 0.70245 - 0.70260, and 206Pb/204Pb = 18.73 - 19.15, indicating a more enriched source. Basalts from the West Valley floor have chemical compositions intermediate between these two end-members. As a group, West Valley basalts from a two-component mixing array in element-element and element-isotope plots which is best explained by magma mixing. Evidence for crustal-level magma mixing in some basalts includes mineral-melt chemical and isotopic disequilibrium, but mixing of melts at depth (within the mantle) may also occur. The mantle beneath the northern Juan de Fuca Ridge is modelled as a plum-pudding, with "plums" of enriched, amphibole-bearing peridotite floating in a depleted matrix (DM). Low degrees of melting preferentially melt the "plums", initially removing only the amphibole component and producing alkaline to transitional E-MORB. Higher degrees of melting tap both the "plums" and the depleted matrix to yield N-MORB. The subtly different isotopic compositions of the E-MORBs compared to the N-MORBs require that any enriched component in the upper mantle was derived from a depleted source. If the enriched component crystallized from fluids with a DM source, the "plums" could evolve to their more evolved isotopic composition after a period of 1.5-2.0 Ga. Alternatively, the enriched component could have formed recently from fluids with a lessdepleted source than DM, such as subducted oceanic crust. A third possibility is that enriched material might be dispersed as "plums" throughout the upper mantle, transported from depth by mantle plumes.

Chemistry of peridotites from ODP 109-670A

During Leg 109 of the Ocean Drilling Program, about 100 m of serpentinized peridotites were drilled on the western wall of the M.A.R. axial rift valley, 45 km south of the Kane Fracture Zone. The present study reports petrological and mineralogical data obtained from 29 small pieces of these ultramafic rocks, including about 60% serpentinized harzburgites, 26% serpentinized lherzolites, 14% serpentinized dunites, and one sample of olivine websterite. Modal analyses show that all these rocks are plagioclase-free four-phase peridotites equilibrated in the spinel lherzolite facies. The estimated average modal composition of the sample set is about 80% olivine, 14% opx, 5% cpx, and 1% spinel, that is, a cpx-poor lherzolite. The well developed porphyroclastic structures and mineralogical characteristics of these rocks indicate their affinity with the group of residual mantle tectonites, among the abyssal peridotites. Features typical of magmatic cumulates are lacking. The high contents in Al2O3 of the cpx (average 5.4%) and of the opx (average 4.3%) porphyroclasts, the low Cr# of the spinels (average 22.9%), and the rather high content in modal cpx (about 5%), indicate a moderate percentage of melting, of the order of 10%-15%. Site 670 peridotites plot close to the least depleted mantle rocks collected in the oceans in most diagrams used to define the average trend of the ocean-floor peridotites. Microprobe traverses across the cores of the exsolved opx and cpx porphyroclasts permitted the recalculation of the magmatic compositions of these pyroxenes: the 'primitive' opx were equilibrated at about 1300°C, probably at the end of the main melting episodes, whereas the 'primitive' cpx show lower equilibration temperatures, at about 1200°C, reflecting a more complex thermal history. The subsolidus evolution is well recorded, from 1200°C to about 950CC, by the exsolved pyroxenes and the olivine and spinel phases. Unusually high blocking temperatures, close to 1000°C, indicate that the peridotite body was cooled very rapidly between 1000°C and the beginning of serpentinization. Oxygen fugacities, calculated for 10 kb and at the blocking temperatures indicated by the olivine/spinel geothermometer, are close to the usual fugacities calculated in oceanic peridotites and basalts (of the order of 10**-10 to 10**-11, on the QFM buffer). Site 670 peridotites have compositions close to those of the peridotites collected in the Kane Fracture Zone area, and obviously belong to the moderately depleted mantle peridotites which characterize abyssal peridotites collected away from mantle plumes and oceanic islands. In particular, they differ from the highly residual harzburgites collected along the M.A.R. over the Azores bulge.

Rare-earth elements and isotope ratios at DSDP Leg 82 Holes

We report 48 analyses of rare-earth elements (REE) and 15 143Nd/144Nd and 87Sr/86Sr analyses for basalts from the eight holes drilled during Leg 82. Discrete and distinct REE patterns and 143Nd/144Nd ratios characterize the eight holes, with little variation observed downhole except in Holes 561 and 558, thus suggesting dominantly long-term temporal and large-scale spatial variations in the mantle source of these basalts beneath the Mid-Atlantic Ridge over the last 35 Ma of its spreading activity. There is a good inverse correlation between 143Nd/144Nd and (La/Sm)EF with one exception in Hole 558 (approximately 35 Ma), the latter suggesting a recent (35 Ma) light REE depletion event, perhaps caused by dynamic or fractional melting. Short-term temporal and small-scale spatial mantle source variability is also evident in Hole 561 (approximately 18 Ma), which has rapid fluctuations in REE patterns and 143Nd/144Nd ratios (suggesting rapid transfer of magma from the time of melting) and is evidence contrary to the presence of a well-mixed magma chamber at this particular site and time. The mantle source variations noted can be interpreted within two extreme models. The first model invokes a convecting mantle depleted in large ion lithophile elements (LILE) and containing lumps (or veins) of LILE-enriched material of various shapes and sizes, passively and randomly distributed throughout. A second more restrictive model considers the interaction of fixed mantle plumes and the LILE-depleted asthenosphere flowing towards a migrating Mid- Atlantic Ridge (MAR) axis. With the exception of Hole 558 and the uncertainties of reconstructions of absolute plate movements in the region, the observed variations can be explained by two hot spots; the nearly ridge-centered Azores hot spot (plume) and another hot spot located beneath the African plate that may be affecting the source of basalts currently erupting at the MAR axis at 35°N and which, in the past, would have produced the New England chain of seamounts on the North American plate and (later) the Atlantis-Great Meteor chain on the African plate. Basalts erupted south of the Hayes Fracture Zone have not been affected by either of these two hot spots over the last 35 Ma and appear to have been continuously derived from the LILE-depleted source. Subaxial flow downridge from the Azores plume appears to have started 9 Ma, on the basis of the southward converging V-shaped time-transgressive ridges branching from the Pico and Corves Island, or not earlier than 16 Ma, on the basis of the geochemical results. Variations within Hole 558 remains unexplained by the latter model, unless we hypothesize a third hot spot.

Mineralogy and geochemistry of a basalt from ODP Hole 119-738C

The Lower Cretaceous tholeiitic basalt cored at Site 738, on the southernmost part of the Kerguelen Plateau, shows anomalous Sr, Nd, and Pb isotopic compositions compared to other lavas from Kerguelen Island and the Kerguelen Plateau. The strongly negative value of eNd (- 8.5) and high 207Pb/204Pb ratio (15.71) reflect a long-term evolution in a source high in Nd/Sm and µ. These geochemical properties, not observed in the Indian Ocean mantle plumes (St. Paul, Kerguelen Islands), have been reported for alkali lavas erupted in East Antarctica, thus suggesting that they originate from the Gondwana subcontinental lithosphere.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.