Metamorphic decarbonation, silicate weathering and the long-term carbon cycle

The seawater ⁸⁷Sr/⁸⁶Sr curve implies a 50–100 Myr episodicity in weathering rate which requires a corresponding variation in CO₂ degassing from the solid earth to the atmosphere. It is proposed that this is caused by orogenesis, which both produces CO₂ as a result of metamorphic decarbonation reactions, and consumes extra CO₂ as a consequence of erosion-enhanced weathering. Global climate on the geological time-scale is therefore contTolled by the difference between the relatively large and variable orogenic-moderated degassing and weathering CO₂ fluxes.     

Fluxes of Sr into the headwaters of the Ganges

Himalayan weathering is recognized as an important agent in modifying sea water chemistry, but there are significant uncertainties in our understanding of Himalayan riverine fluxes. This paper examines causes of the variability, including that of the seasons, by analysis of downstream variations in Sr, ⁸⁷Sr, and major ions in the mainstream, in relation to the composition of tributary streams from subcatchments with differing geologic substrates.Water samples were collected over four periods spanning the premonsoon, monsoon, and postmonsoon seasons. Uncertainties in the relative fluxes have been estimated, using Monte Carlo techniques, from the short-term variability of mainstream chemistry and the scatter of tributary compositions. The results show marked seasonal variations in the relative inputs related to high monsoon rainfall in the High and Lesser Himalaya, contrasting with the major contribution from glacial melt waters from the Tibetan Sedimentary Series (TSS) at times of low rainfall. Much of the spread in previously published estimates of the sources of Sr in Himalayan rivers may result from these seasonal variations in Sr fluxes.The annual fluxes of Sr into the headwaters of the Ganges are derived from the three main tectonic units in the proportions 35 ± 1% from the TSS, 27 ± 3% from the High Himalayan Crystalline Series (HHCS), and 38 ± 8% from the Lesser Himalaya. The particularly elevated ⁸⁷Sr/⁸⁶Sr ratios characteristic of the HHCS and the Lesser Himalaya enhance their influence on seawater Sr-isotope composition. The TSS contributes 13 ± 1%, the HHCS 30 ± 3%, and the Lesser Himalaya 57 ± 11% of the ⁸⁷Sr flux in excess of the seawater ⁸⁷Sr/⁸⁶Sr ratio of 0.709.     

Dynamics and role of the Durban cyclonic eddy in the KwaZulu-Natal Bight ecosystem

The semi-permanent Durban Eddy is a mesoscale, lee-trapped, cold-core cyclonic circulation that occurs off the east coast of South Africa between Durban in the north and Sezela, some 70 km to the south. When present, strong north-eastward countercurrents reaching 100 cm s^–1 are found inshore. It is hypothesised that the cyclone is driven by the strong south-westward flowing Agulhas Current offshore of the regressing shelf edge near Durban. Analysis of ADCP data and satellite imagery shows the eddy to be present off Durban approximately 55% of the time, with an average lifespan of 8.6 days, and inter-eddy periods of 4 to 8 days. After spin-up the eddy breaks loose from its lee position and propagates downstream on the inshore boundary of the Agulhas Current. The eddy is highly variable in occurrence, strength and downstream propagation speeds. There is no detectable seasonal cycle in eddy occurrence, with the Natal Pulse causing more variability than any seasonal signal. A thermistor array deployed in the eddy centre, together with ship CTD data, indicates upward doming of the thermal structure in the eddy core associated with cooler water and nutrients being moved higher in the water column, stimulating primary production. Together with the use of satellite imagery, our findings indicate a second mechanism of upwelling, viz. divergent upwelling in the northern limb of the eddy. Satellite-tracked surface drifters released in the eddy demonstrated the potential for nutrient-rich eddy water to be transported northwards along the inshore regions of the KwaZulu-Natal (KZN) Bight, thus contributing to the functioning of the bight ecosystem, as well as southwards along the KZN and Transkei coasts – both by the eddy migrating downstream and by eddy water being recirculated into the inshore boundary of the Agulhas Current itself.     

Sharks caught in the KwaZulu-Natal bather protection programme,South Africa. 14. The smooth hammerhead shark Sphyrna zygaena (Linnaeus)

The current study provides long-term catch-rate, biological and feeding data for smooth hammerhead sharks, Sphyrna zygaena, caught in South Africa’s KwaZulu-Natal bather protection programme. In total, 2 512 S. zygaena were caught in net installations between 1978 and 2014, and 72 S. zygaena were caught on drumlines between 2007 and 2014. There was no significant log-linear year trend in the net catch rate over time (slope = 0.0054, t = 1.808, p = 0.07). However, there was a significant temporal increase in mean size of the captured sharks (slope = 0.0012, t = 3.502, p < 0.001). A quasi-Poisson generalised additive mixed model showed that increasing latitude, winter months, colder sea temperatures and the deployment of drumlines all had a significant positive effect on the catch rate of sharks in nets. The size frequency of the catch was unimodal, with significantly more females caught in the nets and more males on the drumlines. The majority (93.1%) of all sharks caught were immature and measured between 80 and 120 cm precaudal length. Teleosts and cephalopods dominated the sharks’ diet in terms of all dietary indices. The prey species consumed indicate that immature S. zygaena are feeding primarily within the pelagic zone of shallow coastal habitats.     

A strontium,neodymium and oxygen isotope study of hydrothermal metamorphism and crustal anatexis in the Trois Seigneurs Massif,Pyrenees, France

Nd, Sr, and O isotope analyses have been made on metamorphic and igneous rocks and minerals from a 310–340 Ma Hercynian-age metamorphic terrane in the Pyrenees, France. Lower Paleozoic shales and phyllites have ⁸⁷Sr/⁸⁶Sr values of 0.707–0.717 at 310 Ma, but model values at 310 Ma of 0.709–0.736 (based on assumed depositional age of 450 Ma and an initial ⁸⁷Sr/⁸⁶Sr=0.707). On a regional scale, ⁸⁷Sr/⁸⁶Sr was homogenized to about 0.713 to 0.717 in the higher-grade pelitic schists during metamorphism. Much of this ⁸⁷Sr/⁸⁶Sr exchange occurred at very low grades (below the biotite isograd), but significant changes also accompanied the δ ¹⁸O lowering of the phyllites (+13 to +16) during their transformation to andalusite- and sillimanite-grade schists (δ ¹⁸O=+11 to +12); all of these effects are attributed to pervasive interactions with hydrothermal fluids (Wickham and Taylor 1985). The data also show that a syn-metamorphic plutonic complex, dominated by a biotite granite body, was derived by mixing of a relatively mafic magmatic end-member (⁸⁷Sr/⁸⁶Sr～ 0.7025–0.7050 and δ ¹⁸O～ +7.5 to +8.0) with two metasedimentary sources, both having ⁸⁷Sr/⁸⁶Sr～0.715 and δ ¹⁸O～ +10.0 to +12.0, but with one being more homogeneous than the other. The more homogeneous component and the (mantle-derived?) magmatic end-member dominate at low structural levels within the complex. The less homogeneous end-member that dominates at high levels is clearly derived from the local Paleozoic pelitic schists. A Rb-Sr age of 330±20 Ma was obtained on hornblende from a deep level within the complex, which fixes this age for the regional metamorphism, as well. Although a post-metamorphic granodiorite magma body at Trois Seigneurs also displays heterogeneities in δ ¹⁸O and ⁸⁷Sr/⁸⁶Sr (and thus does not give a clear-cut Rb-Sr isochron), the data are consistent with an emplacement age between 260 and 310 Ma, similar to ages of other late granodiorites in the Pyrenees. ¹⁴³Nd/ ¹⁴⁴Nd is very uniform within the Hercynian crust, both at Trois Seigneurs (? _Nd=?3 to ?7) and elsewhere in the Pyrenees; almost all igneous lithologies have depleted-mantle, mid-Proterozoic model ages, consistent with efficient recycling of crustal material following original crustal accretion in this area at about 1600 Ma or earlier. Rb-Sr mineral ages exhibit a complex cooling history reflecting late Hercynian and Mesozoic thermal events. Our results show that profound homogenization of the ⁸⁷Sr/⁸⁶Sr and ¹⁸O/¹⁶O ratios of large volumes of the crust can occur during regional metamorphism and crustal anatexis, particularly in regions undergoing extensional tectonics. Such processes can significantly modify the isotopic compositions of the protoliths of granitic magmas; this may explain why many peraluminous Hercynian granitoids of Western Europe have anomalously low (⁸⁷Sr/⁸⁶Sr) initial values compared to their probable sedimentary parent rocks.     

Constraining the potential temperature of the Archaean mantle: A review of the evidence from komatiites

The maximum potential temperature of the Archaean mantle is poorly known, and is best constrained by the MgO contents of komatiitic liquids, which are directly related to eruptive temperatures. However, most Archaean komatiites are significantly altered and it is difficult to assess the composition of the erupted liquid. Relatively fresh lavas from the SASKMAR suite, Belingwe Greenstone Belt, Zimbabwe (2.7 Ga) include chills of 25.6 wt.% MgO, and olivines ranging to Fo_93.6, implying eruption at around 1520°C. A chill sample from Alexo Township, Ontario (also 2.7 Ga) is 28 wt.% MgO, and associated olivines range to Fo_94.1, implying eruption at 1560°C. However, inferences of erupted liquids containing 32–33 wt.% MgO, from lavas in the Barberton Greenstone Belt, South Africa (3.45 Ga) and from the Perseverance Complex, Western Australia (2.7 Ga) may be challenged on the grounds that they contain excess (cumulate) olivine, or were enriched in Mg during alteration or metamorphism. Re-interpretation of olivine compositions from these rocks shows that they most likely contained a maximum of 29 wt.% MgO corresponding to an eruption temperature of 1580°C. This composition is the highest liquid MgO content of an erupted lava that can be supported with any confidence. The hottest modern magma, on Gorgona Island (0.155 Ga) contained 18–20% MgO and erupted at circa 1400°C.

If 1580°C is taken as the temperature of the most magnesian known eruption, then the source mantle from which the liquids rose would have been at up to 2200°C at pressures of 18 GPa corresponding to a mantle potential temperature of 1900°C. These temperatures are in excess of the mantle temperatures predicted by secular cooling models, and thus komatiites can only be formed in hot rising convective jets in the mantle. This result requires that Archaean mantle jets may have been 300°C hotter than the Archaean ambient mantle temperature. This temperature difference is similar to the 200–300°C temperature difference between present day jets and ambient mantle temperatures. An important subsidiary result of this study is the confirmation that spinifex rocks may be cumulates and do not necessarily represent liquid compositions.     

On the kinetics of textural equilibration in forsterite marbles

In the aureole of the Beinn an Dubhaich granite, Skye, the minimum observed median forsterite-calcite-calcite dihedral angle varies from 110° at the olivinein isograd to about 165° (the equilibrium value) at the granite-limestone contact. Laboratory experiments were performed to investigate the kinetics of this textural change. It was found that the rate of change of the forsterite-calcite-calcite dihedral angle followed approximately first-order kinetics with an activation energy of 48±4 kJ mol^-1 for fluid-present conditions, and 90 ±4 kJ mol^-1 for fluid-absent conditions. Scanning ion imaging demonstrated that, at least in the early stages of textural change, solution-reprecipitation of the calcite was the rate determining step in the fluid-present runs. The latter result and the value of the activation energy are both consistent with the activation energy found by previous authors for (albeit zeroth order) silicate-aqueous solution solution-reprecipitation reactions. The value of activation energy for the dry data does not correspond to those for either grain boundary or volume diffusion of oxygen or magnesium in forsterite. The mechanism for textural equilibration in the fluid-absent runs is uncertain. Application of the experimentally-derived rate equation to the Beinn an Dubhaich marbles gave activation energies higher than those obtained experimentally. It is concluded from consideration of grain growth effects that activation energies derived from the Beinn an Dubhaich marbles probably reflect textural equilibration under predominantly fluid-absent conditions.     

The transport of heat and matter by fluids during metamorphism

Non-dimensional solutions to the equations for the combined advective and diffusive one-dimensional transport of heat and solute in a layer are derived for fixed temperature/concentration on the boundaries and initial conditions of a linear gradient across the layer or a step function at the lower boundary. The solutions allow distinction of regimes in which advective or diffusive transport of either heat or solute predominate as a function of fluid flux, time and a length scale. The much lower diffusive coefficients for solute than heat results in a significant range of length scales and fluid flux rates characterised by advection of matter and diffusion of heat. The advective velocity of a component is a function of its fluid:rock partition coefficient. The most rapidly transported tracers which partition largely into the fluid phase, such as He, will travel orders of magnitude faster than heat or compatible solutes such as oxygen. Geochemical profiles in boundary layer regions where both advective and diffusive transport are significant are shown to be particularly informative as to properties of the rocks related to fluid flow such as porosity, permeability, time scales and fluid flux rates. The importance of advection can be directly estimated from the asymmetry of the geochemical profiles across individual layers.     

Strontium and oxygen isotope profiles across marble-silicate contacts,Lizzies Basin,East Humboldt Range,Nevada: constraints on metamorphic permeability contrasts and fluid flow

 Sr isotope profiles across marble-silicate rock contacts are used in conjunction with previously published oxygen isotope profiles to constrain fluid movement, porosity and permeability contrasts in migmatitic rocks from Lizzies Basin in the East Humboldt Range, Nevada. The ¹⁸O/¹⁶O systematics in the high-grade sequence have been interpreted to reflect infiltration of ∼2×10² m³/m² of a relatively low ¹⁸O hydrous fluid through the sequence, but with preservation of δ¹⁸O anomalies in thin marble bands due to a 30-fold lower porosity in the marble compared with silicate rocks (Wickham and Peters 1992). The Sr isotope profiles confirm that tracer exchange between marble and silicate rock was primarily by diffusion, and in one case, indicate that porosities differed by less than a factor of four in the ∼10 cm boundary layer which exhibits diffusive modification of ⁸⁷Sr/⁸⁶Sr ratios. This contrasts with modelling of the oxygen isotope profiles which imply porosity contrasts >10 for one marble band and >50 for a second marble band. Either strontium and oxygen isotope diffusion reflect different events (possible if fluid Sr contents varied with time) or porosity varied substantially with the silicate rocks. Oxygen isotope profiles in the deeper part of the metamorphic section in which δ¹⁸O values of silicate rocks have been homogenised and lowered, indicate similar diffusion distances (and thus porosity-time evolution) to oxygen isotopic profiles higher in the section. Comparison of strontium and oxygen isotope diffusion distances constrains fluid Sr contents to between ∼50 and ∼500 ppm deep in the section, but less than ∼10 ppm higher in the section. The difference is related to release of relatively saline, Sr-rich fluids, by the abundant leucogranites and associated skarns deep in the section (cf. Peters and Wickham 1995). Received: 9 December 1994/Accepted: 13 April 1995