Mantle heterogeneity beneath the southern Mid-Atlantic Ridge: trace element evidence for contamination of ambient asthenospheric mantle

We report new trace element data for an extensive suite of quench basalt glasses dredged from the southern Mid-Atlantic Ridge (MAR) between 40°S and 52.5°S. Ratios between highly incompatible trace elements are strongly correlated and indicate a systematic distribution of incompatible element enriched mid-ocean ridge basalt (MORB) (E-type: Zr/Nb=5.9-19, Y/Nb=0.9-8.4, (La/Sm)_n=1.0-2.9) and incompatible element depleted MORB (N-type: Zr/Nb=30-69, Y/Nb=11-29, (La/Sm)_n=0.48-0.79) along this section of the southern MAR. A notable feature of N-type MORB from the region is the higher than usual Ba/Nb (4-9), La/Nb (1.2-2.4) and primitive mantle normalised K/Nb ratios (>1). Ba/Nb ratios in E-type MORB samples from 47.5 to 49°S are especially elevated (>10). The occurrence and geographic distribution of E-type MORB along this section of the southern MAR can be correlated with the ridge-centred Shona and off-axis Discovery mantle plumes. In conjunction with published isotope data for a subset of the same sample suite [Douglass et al., J. Geophys. Res. 104 (1999) 2941], a model is developed whereby prior to the breakup of Gondwana and the opening of the South Atlantic Ocean, the underlying asthenospheric mantle was locally contaminated by fluids/melts rising from the major Mesozoic subduction zone along the south-southwest boundary of Gondwana, leaving a subduction zone geochemical imprint (elevated (K/Nb)_n and ⁸⁷Sr/⁸⁶Sr ratios, decreased ¹⁴³Nd/¹⁴⁴Nd ratios). Subsequent impingement of three major mantle plume heads (Tristan/Gough, Discovery, Shona) resulted in heating and thermal erosion of the lowermost subcontinental lithosphere and dispersal into the convecting asthenospheric mantle. With the opening of the ocean basin, continued plume upwelling led to plume-ridge interactions and mixing between geochemically enriched mantle derived from the Shona and Discovery mantle plumes, material derived from delamination of the subcontinental lithosphere, and mildly subduction zone contaminated depleted asthenospheric mantle.     

Oxygen isotope evidence for the origin of enriched mantle beneath the mid-Atlantic ridge

Geochemical variations in mid-ocean ridge basalts have been attributed to differing proportions of compositionally distinct mantle components in their sources, some of which may be recycled crust. Oxygen isotopes are strongly fractionated by near-surface interactions of rocks with the hydrosphere, and thus provide a tracer of near-surface materials that have been recycled into the mantle. We present here oxygen isotope analyses of basaltic glasses from the mid-Atlantic ridge south of and across the Azores platform. Variations in δ¹⁸O in these samples are subtle (range of 0.47‰) and may partly reflect shallow fractional crystallization; we present a method to correct for these effects. Relatively high fractionation-corrected δ¹⁸O in these samples is associated with geochemical indices of enrichment, including high La/Sm, Ce/Pb, and ⁸⁷Sr/⁸⁶Sr and low ¹⁴³Nd/¹⁴⁴Nd. Our results suggest two first-order conclusions about these enriched materials: (1) they are derived (directly or indirectly) from recycled upper oceanic crustal rocks and/or sediments; and (2) these materials are present in the north Atlantic MORB sources in abundances of less than 10% (average 2–5%). Modeling of variations of δ¹⁸O with other geochemical variables further indicates that the enriched component is not derived from incorporation of sediment or bulk altered oceanic crust, from metasomatism of the mantle by hydrous or carbonate-rich fluids, or from partial melting of subducted sediment. Instead, the data appear to require a model in which the enriched component is depleted mantle that has been metasomatized by small-degree partial melts of subducted, dehydrated, altered oceanic crust. The age of this partial melting is broadly constrained to 250 Ma. Reconstructed plate motions suggest that the enriched component in the north Atlantic mantle may have originated by subduction along the western margin of Pangea.     

Trapping of aeolian sediments and build‐up of the ice cover of a dry‐based Antarctic lake

Among the perennially frozen lakes of the Dry Valleys of South Victoria Land (Antarctica), some are dry‐based, i.e. frozen to the bottom. One of these is studied here by a multiparametric investigation (isotopic composition in δD and δ¹⁸O, ions, gas and ice texture analyses). A sediment layer about 10 cm thick appearing at a depth of 3·5 m is also studied by grain size, X‐ray diffraction and scanning electron microscope analyses. The information retrieved indicates that this ice‐block lake results from a build‐up in two steps and explains how aeolian sediments were included as a layer into the ice. Copyright © 2002 John Wiley & Sons, Ltd.     

Simulating the effects of geologic heterogeneity and transient boundary conditions on streambed temperatures — Implications for temperature-based water flux calculations

Analytical solutions to the one-dimensional heat transport equation for steady-state conditions can provide simple means to quantify groundwater surface water exchange. The errors in exchange flux calculations that are introduced when the underlying assumptions of homogeneous sediments and constant temperature boundary conditions are violated were systematically evaluated in a simulation study. Temperatures in heterogeneous sediments were simulated using a numerical model. Heterogeneity in the sediments was represented by discrete, binary geologic units. High contrasts between the hydraulic conductivities (K) of the geologic units were found to lead to large errors, while the influence of the structural arrangement of the units was smaller. The effects of transient temperature boundary conditions were investigated using an analytical equation. Errors introduced by transient boundary conditions were small for Darcy-velocities > 0.1 m d^− 1 in the period near maximum and minimum annual surface water temperatures. For smaller fluxes, however, errors can be large. Assuming steady-state conditions and vertical flow in homogeneous sediments is acceptable at certain times of the year and for medium to high exchange fluxes, but pronounced geologic heterogeneity can lead to large errors.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Anisotropy and depth‐related heterogeneity of hydraulic conductivity in a bog peat. II: modelling the effects on groundwater flow

Anisotropy and heterogeneity of hydraulic conductivity (K) are seldom considered in models of mire hydrology. We investigated the effect of anisotropy and heterogeneity on groundwater flow in bog peat using a steady‐state groundwater model. In five model simulations, four sets of K data were used. The first set comprised measured K values from an anisotropic and heterogeneous bog peat. These data were aggregated to produce the following simplified data sets: an isotropic and heterogeneous distribution of K; an isotropic and homogeneous distribution; and an anisotropic and homogeneous distribution. We demonstrate that, where anisotropy and heterogeneity exist, groundwater flow in bog peat is complex. Fine‐scale variations in K have the potential to influence patterns and rates of groundwater flow. However, for our data at least, it is heterogeneity and not anisotropy that has the greater influence on producing complex patterns of groundwater flow. We also demonstrate that patterns and rates of groundwater flow are simplified and reduced when measured K values are aggregated to create a more uniform distribution of K. For example, when measured K values are aggregated to produce isotropy and homogeneity, the rate of modelled seepage is reduced by 28%. We also show that when measured K values are used, the presence of a drainage ditch can increase seepage through a modelled cross‐section. Our work has implications for the accurate interpretation of hydraulic head data obtained from peat soils, and also the understanding of the effect of drainage ditches on patterns and rates of groundwater flow. Copyright © 2002 John Wiley & Sons, Ltd.     

Geochemistry and Sr–Nd–Pb isotopic constraints on the petrogenesis of Cenozoic lavas from the Pali Aike and Morro Chico area (52°S), southern Patagonia,South America

Geochemical and isotopic analyses (Sr–Nd–Pb) of late Miocene to Quaternary plateau lavas from the Pali Aike and Morro Chico areas (52°S) were undertaken to constrain the melting processes and mantle sources that contributed to magma generation and the geodynamic evolution of southernmost Patagonia, South America. The Pali Aike and Morro Chico lavas are alkaline (Pali Aike, 45–49 wt.% SiO₂; 4.3–5.9 wt.% Na₂O+K₂O) and subalkaline (Morro Chico, 50.5–50.8 wt.% SiO₂; 4.0–4.4 wt.% Na₂O+K₂O), relatively primitive (Pali Aike, 9.5–13.7 wt.% MgO; Morro Chico, 7.6–8.8 wt.% MgO) mafic volcanic rocks that have typical intraplate ocean island basalt‐like signatures. Incompatible trace element ratios and isotopic ratios of the Pali Aike and Morro Chico lavas differ from those of the majority of Neogene southern Patagonian slab window lavas in showing more enriched characteristics and are similar to high‐μ (HIMU)‐like basalts. The rare earth element (REE) modeling to constrain mantle melting percentages suggests that these lavas were produced by low degrees of partial melting (1.0–2.0% for Pali Aike lavas and about 2.6–2.7% for Morro Chico lavas) of a garnet lherzolite mantle source. The major systematic variations of Sr–Nd–Pb isotopes in southern Patagonian lavas are related to geographic location. The Pali Aike and Morro Chico lavas from the southernmost part of Patagonia have lower ⁸⁷Sr/⁸⁶Sr and higher ¹⁴³Nd/¹⁴⁴Nd and ²⁰⁶Pb/²⁰⁴Pb ratios, relative to most of the southern Patagonian lavas erupted north of 49.5°S, pointing to a HIMU‐like signature. An isotopically depleted and HIMU‐like asthenospheric domain may have been the main source of magmas in the southernmost part of Patagonia (e.g. Pali Aike, Morro Chico, and Camusu Aike volcanic field), suggesting the presence of a major discontinuity in the isotopic composition of the asthenosphere in southern Patagonia. On the basis of geochemical and isotope data and the available geological and geotectonic reconstructions, a link between the HIMU asthenospheric mantle domain beneath southernmost Patagonia and the HIMU mega‐province of the southwestern Pacific Ocean is proposed.