The Banhado Alkaline Complex, Southeastern Brazil: source and evolution of potassic SiO2-undersaturated high-Ca and low-Ca magmatic series

The Cretaceous Banhado alkaline complex in southeastern Brazil presents two potassic SiO2-undersaturated series. The high-Ca magmatic series consist of initially fractionated olivine (Fo(92-91)) + diopside (Wo(48-43)En(49-35)Ae(0-7)), as evidenced by the presence of xenocrysts and xenoliths. In that sequence, diopside (Wo(47-38)En(46-37)Ae(0-8)) + phlogopite + apatite + perovskite (Prv(> 92)) crystallized to form the phlogopite melteigite and led to the Ca enrichment of the magma. Diopside (Wo(47-41)En(32-24) Ae(3-14)) continued to crystallize as an early mafic mineral, followed by nepheline (Ne(74.8-70.1)Ks(26.3-21.2)Qz(7.6-0.9)) and leucite (Lc(65-56)) and subsequently by melanite and potassic feldspar (Or(85-99)Ab(1-7)) to form melanite ijolites, wollastonite-melanite urtites and melanite-nepheline syenites. Melanite-pseudoleucite-nepheline syenites are interpreted to be a leucite accumulation. Melanite nephelinite dykes are believed to represent some of the magmatic differentiation steps. The low-Ca magmatic series is representative of a typical fractionation of aegirine-augite (Wo(36-29)En(25-4)Ae(39-18)) + alkali feldspar (Or(57-96)Ab(3-43)) + nepheline (Ne(76.5-69.0)Ks(19.9-14.4)Qz(15.1-7.7)) + titanite from phonolite magma. The evolution of this series from potassic nepheline syenites to sodic sodalite syenites and sodalitolites is attributed to an extensive fractionation of potassic feldspar, which led to an increase of the NaCl activity in the melt during the final stages forming sodalite-rich rocks. Phonolite dykes followed a similar evolutionary process and also registered some crustal assimilation. The mesocratic nepheline syenites showed interactions with phlogopite melteigites, such as compatible trace element enrichments and the presence of diopside xenocrysts, which were interpreted to be due to a mixing/mingling process of phonolite and nephelinite magmas. The geochemical data show higher TiO2 and P2O5 contents and lower SiO2 contents for the high-Ca series and different LILE evolution trends and REE chondrite-normalized patterns as compared to the low-Ca series. The Sr-87/Sr-86, Nd-143/Nd-144, Pb-206/Pb-204 and Pb-208/Pb-204 initial ratios for the high-Ca series (0.70407-0.70526, 0.51242-0.51251, 17.782-19.266 and 38.051-39.521, respectively) were slightly different from those of the low-Ca series (0.70542-0.70583, 0.51232-0.51240, 17.758-17.772 and 38.021-38.061, respectively). For both series, a CO2-rich potassic metasomatized lithospheric mantle enriched the source with rutile-bearing phlogopite clinopyroxenite veins. Kamafugite-like parental magma is attributed to the high-Ca series with major contributions from the melting of the veins. Potassic nephelinite-like parental magma is assigned to the low-Ca series, where the metasomatized wall-rock played a more significant role in the melting process.

Crustal heterogeneities and complexities of fault processes

In this work we study the relation between crustal heterogeneities and complexities in fault processes. The first kind of heterogeneity considered involves the concept of asperity. The presence of an asperity in the hypocentral region of the M = 6.5 earthquake of June 17-th, 2000 in the South Iceland Seismic Zone was invoked to explain the change of seismicity pattern before and after the mainshock: in particular, the spatial distribution of foreshock epicentres trends NW while the strike of the main fault is N 7◦ E and aftershocks trend accordingly; the foreshock depths were typically deeper than average aftershock depths. A model is devised which simulates the presence of an asperity in terms of a spherical inclusion, within a softer elastic medium in a transform domain with a deviatoric stress field imposed at remote distances (compressive NE − SW, tensile NW − SE). An isotropic compressive stress component is induced outside the asperity, in the direction of the compressive stress axis, and a tensile component in the direction of the tensile axis; as a consequence, fluid flow is inhibited in the compressive quadrants while it is favoured in tensile quadrants. Within the asperity the isotropic stress vanishes but the deviatoric stress increases substantially, without any significant change in the principal stress directions. Hydrofracture processes in the tensile quadrants and viscoelastic relaxation at depth may contribute to lower the effective rigidity of the medium surrounding the asperity. According to the present model, foreshocks may be interpreted as induced, close to the brittle-ductile transition, by high pressure fluids migrating upwards within the tensile quadrants; this process increases the deviatoric stress within the asperity which eventually fails, becoming the hypocenter of the mainshock, on the optimally oriented fault plane. In the second part of our work we study the complexities induced in fault processes by the layered structure of the crust. In the first model proposed we study the case in which fault bending takes place in a shallow layer. The problem can be addressed in terms of a deep vertical planar crack, interacting with a shallower inclined planar crack. An asymptotic study of the singular behaviour of the dislocation density at the interface reveals that the density distribution has an algebraic singularity at the interface of degree ω between -1 and 0, depending on the dip angle of the upper crack section and on the rigidity contrast between the two media. From the welded boundary condition at the interface between medium 1 and 2, a stress drop discontinuity condition is obtained which can be fulfilled if the stress drop in the upper medium is lower than required for a planar trough-going surface: as a corollary, a vertically dipping strike-slip fault at depth may cross the interface with a sedimentary layer, provided that the shallower section is suitably inclined (fault "refraction"); this results has important implications for our understanding of the complexity of the fault system in the SISZ; in particular, we may understand the observed offset of secondary surface fractures with respect to the strike direction of the seismic fault. The results of this model also suggest that further fractures can develop in the opposite quadrant and so a second model describing fault branching in the upper layer is proposed. As the previous model, this model can be applied only when the stress drop in the shallow layer is lower than the value prescribed for a vertical planar crack surface. Alternative solutions must be considered if the stress drop in the upper layer is higher than in the other layer, which may be the case when anelastic processes relax deviatoric stress in layer 2. In such a case one through-going crack cannot fulfil the welded boundary conditions and unwelding of the interface may take place. We have solved this problem within the theory of fracture mechanics, employing the boundary element method. The fault terminates against the interface in a T-shaped configuration, whose segments interact among each other: the lateral extent of the unwelded surface can be computed in terms of the main fault parameters and the stress field resulting in the shallower layer can be modelled. A wide stripe of high and nearly uniform shear stress develops above the unwelded surface, whose width is controlled by the lateral extension of unwelding. Secondary shear fractures may then open within this stripe, according to the Coulomb failure criterion, and the depth of open fractures opening in mixed mode may be computed and compared with the well studied fault complexities observed in the field. In absence of the T-shaped decollement structure, stress concentration above the seismic fault would be difficult to reconcile with observations, being much higher and narrower.

Active tectonics in the NW-German Basin: Evidence from correlations between the modern landscape and deep geological structures (Lower Saxony, river Hunte)

KurzfassungIm Einzugsgebiet der Hunte (NW-Deutsches Becken, Niedersachsen) wurde untersucht, ob die Landschaftsgenese durch tektonische Bewegungen der Oberkruste beeinflußt ist. Krustenbewegungen führten im Bereich einer Hauptschollengrenze zu einer Hebung der weichselzeitlichen Niederterrasse (durchschnittliche Hebungssrate von ~0,5 mm/a über die letzten 12000 Jahre). Tektonischer Einfluß auf die heutige Landoberfläche ist über einem permischen Salzkissen zu verzeichnen, wo sich das Gefälle der holozänen Aue umkehrt. Krustenbewegungen haben mit großer Wahrscheinlichkeit Vorzugsrichtungen verursacht, die an der Tertiärbasis und in der heutigen Landschaft nachweisbar sind (0-5° und 90-95°). Das Abfließen der Hunte nach Norden scheint durch eine aktive, nordwärts gerichtete Kippung des NW-Deutschen Beckens verursacht zu sein. Hohe lineare Korrelationskoeffizienten zwischen Tiefenlage der Tertiärbasis und Höhenlage der heutigen Landoberfläche weisen auf eine aktive Kippung des Beckens hin. Beckensubsidenz hat möglicherweise die Akkumulation der weichselzeitlichen Niederterrasse gesteuert, da eine Übereinstimmung zwischen rezenter Beckensubsidenz und durchschnittlicher Sedimentationsrate des Niederterrassenkörpers besteht. Untersuchungen an einer geschlossenen Hohlform deuten auf eine aktive Sackungsstruktur hin, da sich Anomalien des geologischen Untergrundes mit der topographischen Lage der Struktur decken.

Petrogenesis of dunite, wehrlite and websterite xenoliths from Kimberley, South Africa: Origin as mantle peridotites of cumulates

Dunite, wehrlite and websterite xenoliths occur amongst a large abundance of mantle xenoliths in kimberlites of the Kimberley cluster in South Africa. Up to know they have mostly been neglected. On the basis of texture, major and trace elements, oxygen isotopes as well as Re-Os isotope characteristics, they can be subdivided into two groups. A coarse-grained mantle peridotite group, comprising dunite, wehrlite and websterite xenoliths, that are similar to fertile peridotites and represent upper mantle assemblages that are differently influenced by mantle metasomatism. And a cumulate group, containing fine-grained Fe-rich dunite xenoliths that represent cumulates of flood basalt magmatism related to ~183 Ma Karoo and ~2.7 Ga Ventersdorp events in southern Africa. Dunite, wehrlite and websterite xenoliths have preserved a complex history of melt depletion and metasomatic re-enrichment events, which gives information about the different re-enrichment stages of the subcratonic lithospheric mantle and the spatial differences within the Kaapvaal craton upper mantle. Websterite xenoliths comprise orthopyroxene (40-85 Vol. %), clinopyroxene (5-42 Vol. %), garnet (4-10 Vol. %) and subordinately olivine, while dunite and wehrlite xenoliths contain predominantly olivine (65-100 Vol %) and subordinately orthopyroxene, clinopyroxene and garnet. High melt depletion and a dunitic to harzburgitic protolith composition are reflected by high forsterite (Fo90-92) and high olivine NiO contents (2800-5000 ppm) and high orthopyroxene Mg# (Mg/(Mg+Fe)) of 0.91-0.93. Re-depletion ages of predominantly 2.9 Ga reflect a minimum age of melt depletion. Melt depletion ceased in conjunction with collision of the Kimberley block with the Witwatersrand block ~2.9 Ga ago. Subduction related re-fertilisation of the previously depleted mantle xenoliths is documented by i) amoeboid textured orthopyroxene, clinopyroxene and garnet, which crystallized in schlieren along olivine grain boundaries, ii) high whole-rock SiO2, Al2O3, CaO, TiO2, FeO contents, iii) low oxygen isotope ratios in clinopyroxene and garnet of 4.8-5.4 ‰ and 4.7-5.3 ‰, respectively and iv) trace element compositions of wehrlitic clinopyroxene and garnet in equilibrium with high-pressure partial melts of eclogite. Trace element disequilibrium of orthopyroxene with clinopyroxene and garnet indicates a separate origin for orthopyroxene, on one side as primary mantle orthopyroxene in dunite and wehrlite xenoliths and on the other side as reaction product with Si-rich melts produced by partial melting of eclogite. This reaction triggered replacement of olivine by orthopyroxene in the surrounding mantle and produced the typical Si-rich composition of Kaapvaal mantle peridotites. Partial melting of eclogite at higher temperatures produced a second metasomatic melt with lower SiO2, but higher Al2O3, CaO, FeO, Ti, Zr, Hf and a low oxygen isotope ratio. This melt triggered clinopyroxene and locally garnet and rutile crystallization in percolation veins, replacing olivine and orthopyroxene in the Kaapvaal upper mantle. Additionally, websterite xenoliths have experienced late stage cryptic metasomatism by the host kimberlite melt, changing the trace element composition of clinopyroxene, orthopyroxene and garnet to different extent. Hence websterite and most fertile lherzolite xenoliths have experienced three metasomatic events: i) reaction with high-Si melt, ii) percolation of subduction related silica melt with lower SiO2 content and iii) cryptic metasomatism by kimberlite. In contrast, dunite and wehrlite xenoliths have only experienced the second metasomatic event. They represent mantle lithologies further away from metasomatising agents. The Fe-rich dunites comprise olivine neoblasts with subordinate olivine porphyroclasts and parallel-orientated needles of ilmenite, which may enclose spinel. The lower forsterite and NiO contents of olivine in Fe-rich dunites compared to mantle peridotite xenoliths (Fo87-89 vs. Fo93-95 and 1300-2800ppm vs. 2200-3900 ppm, respectively), rules out a restitic origin. Cr-rich spinels are remnants of the original cumulate mineralogy that survived a late stage metasomatic overprint related to the production of the host kimberlite, producing ilmenite and phlogopite in some samples. Olivine porphyroclasts and neoblasts have different trace element compositions, the latter having high Ti, V, Cr and Ni and low Zn, Zr and Nb contents, indicating contrasting origins for neoblasts and porphyroclasts. The dunites have high 187Os/188Os ratios (0.11-0.15) indicating young (Phanerozoic) model ages for most samples, whereas three samples show isotopic mixtures between Phanerozoic neoblasts and ancient porphyroclastic material. Most Fe-rich dunite xenoliths can be interpreted as cumulates of fractional crystallization of Karoo magmatism, whereas the porphyroclasts are interpreted to be remnants from the much earlier Archaean Ventersdorp magmatism.

Magma-Driven Multiple Dike Propagation and Fracture Toughness of Crustal Rocks

Dike swarms consisting of tens to thousands of subparallel dikes are commonly observed at Earth's surface, raising the possibility of simultaneous propagation of two or more dikes at various stages of a swarm's development. The behavior of multiple propagating dikes differs from that of a single dike owing to the interacting stress fields associated with each dike. We analyze an array of parallel, periodically spaced dikes that grow simultaneously from an overpressured source into a semi-infinite, linear elastic host rock. To simplify the analysis, we assume steady state (constant velocity) magma flow and dike propagation. We use a perturbation method to analyze the coupled, nonlinear problem of multiple dike propagation and magma transport. The stress intensity factor at the dike tips and the opening displacements of the dike surfaces are calculated. The numerical results show that dike spacing has a profound effect on the behavior of dike propagation. The stress intensity factors at the tips of parallel dikes decrease with a decrease in dike spacing and are significantly smaller than that for a single dike with the same length. The reduced stress intensity factor indicates that, compared to a single dike, propagation of parallel dikes is more likely to be arrested under otherwise the same conditions. It also implies that fracture toughness of the host rock in a high confining pressure environment may not be as high as inferred from the propagation of a single dike. Our numerical results suggest fracture toughness values on the order of 100 MPa root m. The opening displacements for parallel dikes are smaller than that for a single dike, which results in higher magma pressure gradients in parallel dikes and lower flux of magma transport.

Influence of Martian crustal magnetic anomalies on the emission of energetic neutral hydrogen atoms

We analyze the data on hydrogen energetic neutral atoms (ENAs) emissions from the dayside of Mars, recorded by a Neutral Particle Detector of the Analyzer of Space Plasmas and Energetic Atoms aboard Mars Express from 14 March to 9 July 2004. We first identify and analyze events of the ENA flux enhancement coinciding with the presence of the crustal magnetic anomalies on the dayside of Mars. We then backtrace the ENA emissions to the lower altitudes (source region) and build up an average map of the flux intensities in the geographic coordinates with all the available data. The map shows a peak-to-valley ENA flux enhancement of 40%–90% close to the crustal magnetic anomaly regions. These results suggest the influence of the magnetic anomalies on the ENA emission from the dayside of Mars. The enhancement may result from the deviation of the highly directional plasma flow above anomalies toward the detectors such that more charge exchange ENAs would be recorded. Alternatively, higher exospheric densities above the anomalies would also result in an increase of the charge exchange ENA flux.

Geochemistry of the lower sheeted dike complex of ODP Holes 137-504B and 140-504B

Sixty-three samples representing 379 m of sheeted dikes from Deep Sea Drilling Project/Ocean Drilling Program Site 504B have been analyzed for major and selected trace elements by X-ray fluorescence. The samples range from microcrystalline aphyric basalts to moderately phyric (2%-10% phenocrysts) diabase that are typically multiply saturated with plagioclase, olivine, and clinopyroxene, in order of relative abundance. All analyzed samples are classified as Group D compositions with moderate to slightly elevated compatible elements (MgÆ-value = 0.65% ± 0.03%; Al2O3 = 15.5% ± 0.8%; CaO = 13.0% ± 0.3%; Ni = 114 ± 29 ppm), and unusually depleted levels of moderate to highly incompatible elements (Nb < 1 ppm; Zr = 44 ± 7 ppm; Rb < 0.5 ppm; Ba ~ 1 ppm; P2O5 = 0.07% ± 0.02%). These compositions are consistent with a multistage melting of a normal ocean ridge basaltic mantle source followed by extensive fractionation of olivine, plagioclase, and clinopyroxene. Leg 140 aphyric to sparsely phyric (0%-2% phenocrysts) basalts and diabases are compositionally indistinguishable from similarly phyric samples at higher levels in the hole. An examination of the entire crustal section, from the overlying volcanics through the sheeted dikes observed in Leg 140, reveals no significant trends indicating the enrichment or depletion of Costa Rica Rift Zone source magmas over time. Similarly, significant trends toward increased or decreased differentiation cannot be identified, although compositional patterns reflecting variable amounts of phenocryst addition are apparent at various depths. Below ? 1700 mbsf to the bottom of the Leg 140 section, there is a broadly systematic pattern of Zn depletion with depth, the result of high-temperature hydrothermal leaching. This zone of depletion is thought to be a significant source of Zn for the hydrothermal fluids depositing metal sulfides at ridge-crest hydrothermal vents and the sulfide-mineralization zone, located in the transition between pillow lavas and sheeted dikes. Localized zones of intense alteration (60%-95% recrystallization) are present on a centimeter to meter scale in many lithologic units. Within these zones, normally immobile elements Ti, Zr, Y, and rare-earth elements are strongly depleted compared with "fresher" samples centimeters away. The extent of compositional variability of these elements tends to obscure primary igneous trends if the highly altered samples are not identified or removed. At levels up to 40% (or possibly 60%) recrystallization, Ti, Zr, and Y retain their primary signatures. Although the mechanisms are unclear, it is possible that these intense alteration zones are a source of Y and rare-earth elements for the typically rare-earth-element-enriched hydrothermal vent fluids of mid-ocean ridges.

The petrology of the lower series volcanics of ODP Hole 104-642E

Between 1086.6 and 1229.4 m below seafloor at Site 642 on the Outer Vøring Plateau, a series of intermediate volcanic extrusive flow units and volcaniclastic sediments was sampled. A mixed sequence of dacitic subaerial flows, andesitic basalts, intermediate volcaniclastics, subordinate mid-ocean ridge basalt, (MORB) lithologies, and intrusives was recovered, in sharp contrast to the more uniform tholeiitic T-type MORB units of the overlying upper series. This lower series of volcanics is composed of three chemically distinct groups, (B, A2, A1), rather than the two previously identified. Flows of the dacitic group (B) have trace-element and initial Sr isotope signatures which indicate that their source magma derived from the partial melting of a component of continental material in a magma chamber at a relatively high level in the crust. The relative proportions of crustal components in this complex melt are not known precisely. The most basic group (A2) probably represents a mixture of this material with MORB-type tholeiitic melt. A third group (A1), of which there was only one representative flow recovered, is chemically intermediate between the two groups above, and may suggest a repetition of, or a transition phase in, the mixing processes.

Geochemistry of Qeqertaa xenoliths

This study presents electron microprobe data for dunite xenoliths from a lamprophyre dyke located on the island of Qeqertaa, West Greenland. The minimum age of this dyke is Palaeoproterozoic and it experienced amphibolite facies metamorphism and deformation during that era. The samples consist of nearly 200 xenoliths with a size range of 0.5-8 cm. These dunite xenoliths have olivine Mg#, that range from 80.3 to 94.6 (n = 579) with a mean of 92.6. Orthopyroxene is found in three xenoliths and garnet in five others. The latter suggests the depth of the Qeqertaa xenolith suite to be near the diamond stability-field, which is substantiated by the finding of diamonds in bulk samples of the Qeqertaa dyke. This further indicates the presence of a lithospheric mantle domain dominated by high-Mg# dunite to this depth in Palaeoproterozoic time. Cr-rich spinel, in the 0.1-0.2 mm size range, is found within and between olivine grains in individual xenoliths. These Cr-spinels yield Fe-Mg exchange temperatures of 400-600°C. However, the presence of intermediate spinel compositions spanning the lower temperature solvus suggests that equilibration temperatures were >550°C. Fe3+#, expressed as 100xFe3+/(Fe3++Al+Cr)), is shown to be a useful parameter in order to screen for altered spinel (Fe3+#>10) with disturbed Mg# and Cr#. The screened spinel data (Fe3+#<10) show a distinctly different trend in terms of spinel Cr# versus Mg#, compared to unmetamorphosed xenoliths in Tertiary lavas and dikes from Ubekendt Ejland and Wiedemann Fjord, respectively, also located within the North Atlantic craton. This difference likely reflects amphibolite facies metamorphic resetting of the Qeqertaa xenolith suite by Fe-Mg exchange. Given the similarity of the Qeqertaa xenolith suite with the Ubekendt and Wiedemann suites, in terms of their olivine Mg# and spinel Cr# distribution, high-Mg# dunite is likely to be an important component of the subcontinental lithospheric mantle beneath the North Atlantic craton and appears to have spanned a vertical distance of at least 150 km in this region, even during the Palaeoproterozoic.

(Figures S2a-c) Phase diagrams of the upper crust, the lower crust and the upper mantle based on averaged compositions of the oceanic crust and mantle

The age of the subducting Nazca Plate off Chile increases northwards from 0 Ma at the Chile Triple Junction (46°S) to 37 Ma at the latitude of Valparaíso (32°S). Age-related variations in the thermal state of the subducting plate impact on (a) the water influx to the subduction zone, as well as on (b) the volumes of water that are released under the continental forearc or, alternatively, carried beyond the arc. Southern Central Chile is an ideal setting to study this effect, because other factors for the subduction zone water budget appear constant. We determine the water influx by calculating the crustal water uptake and by modeling the upper mantle serpentinization at the outer rise of the Chile Trench. The water release under forearc and arc is determined by coupling FEM thermal models of the subducting plate with stability fields of water-releasing mineral reactions for upper and lower crust and hydrated mantle. Results show that both the influx of water stored in, and the outflux of water released from upper crust, lower crust and mantle vary drastically over segment boundaries. In particular, the oldest and coldest segments carry roughly twice as much water into the subduction zone as the youngest and hottest segments, but their release flux to the forearc is only about one fourth of the latter. This high variability over a subduction zone of < 1500 km length shows that it is insufficient to consider subduction zones as uniform entities in global estimates of subduction zone fluxes. This article is protected by copyright. All rights reserved.

(Table T1) Geochemical composition of xenoliths and host rocks from ODP Hole 193-1188A and dredge haul MD-138, PACMANUS field

A number of intensely altered, dark xenoliths with palimpsest quench textures were recorded within altered dacitic host rocks at Site 1189 (Roman Ruins, PACMANUS) during Ocean Drilling Program (ODP) Leg 193. Several of these displayed puzzling marginal fringes, apparently of altered plagioclase with variolitic texture, protruding into adjacent host rocks. Despite their alteration, the xenoliths were interpreted as fragments of rapidly chilled, possibly olivine-bearing basalts incorporated into the dacitic magmas either within the crustal plumbing system or during eruption at the seafloor (figures F15, F16, F17, F42, and F43 in Shipboard Scientific Party, 2002, doi:10.2973/odp.proc.ir.193.104.2002). An additional example of formerly spinifex-textured xenolith, the first from Site 1188 (Snowcap) and the first from the upper cristobalite-bearing zone of alteration, has been revealed by postcruise studies. Furthermore, a pristine sample of the parent lithology has been found within a dredge haul (MD-138, Binatang-2000 Cruise of Franklin; 3°43.60'S, 151°40.35'E, 1688 meters below sea level) from the Satanic Mills hydrothermal field at PACMANUS, near ODP Site 1191. The purpose of this report is to document these discoveries and thereby to confirm and build on shipboard interpretations. To my knowledge, similar xenoliths have never before been found in modern island arc or backarc volcanic sequences. Spinifex textures are most common in Archean komatiites, some of which are bimodally associated with calc-alkaline felsic volcanic rocks.