Dissolved rhenium in the Yamuna river system and the Ganga in the Himalaya: role of black shale weathering on the budgets of Re, Os, and U in rivers and CO2 in the atmosphere

Extensive measurements of dissolved Re and major ion abundances in the Yamuna River System (YRS), a major tributary of the Ganga, have been performed along its entire stretch in the Himalaya, from its source near the Yamunotri Glacier to its outflow at the foothills of the Himalaya at Saharanpur. In addition, Re analysis has been made in granites and Precambrian carbonates, some of the major lithologies of the drainage basin. These data, coupled with those available for black shales in the Lesser Himalaya, allow an assessment of these lithologies’ contributions to the Re budget of the YRS.The Re concentrations in the YRS range from 0.5 to 35.7 pM with a mean of 9.4 pM, a factor of ∼4 higher than that reported for its global average concentration in rivers. Dissolved Re and ΣCations∗ (= Na∗+K+Ca+Mg) are strongly correlated in the YRS, indicating that they are released to these waters in roughly the same proportion throughout their course. The Re/ΣCations∗ in most of these rivers are one to two orders of magnitude higher than the (Re/Na+K+Mg+Ca) measured in granites of the Yamuna basin. This leads to the conclusion that, on average, granites/crystallines make only minor contributions to the dissolved Re budget of the YRS on a basin-wide scale, though they may be important for rivers with low dissolved Re. Similarly, Precambrian carbonates of the Lesser Himalaya do not seem to be a major contributor to dissolved Re in these rivers, as their Re/(Ca+Mg) is much less than those in the rivers. The observation that Re concentrations in rivers flowing through black shales and in groundwaters percolating through phosphorite-black shale-carbonate layers in phosphorite mines are high, and that Re and SO₄ are significantly correlated in YRS, seems to suggest that the bulk of the dissolved Re is derived from black shale/carbonaceous sediments. Material balance considerations, based on average Re of 30 ng g⁻¹ in black shales from the Lesser Himalaya, require that its abundance in the drainage basin of the YRS needs to be a few percent to yield average Re of 9.4 pM. Furthermore, the positive correlation between Re and ΣCations∗ would require that these Re-rich sediments (e.g., black shales) and Re-poor lithologies (e.g., crystallines, Precambrian carbonates) contribute Re and cations in roughly the same proportion throughout the drainage basin. The available data on the abundance and distribution of black shales in the basin are not adequate to test if these requirements can be met.The annual fluxes of dissolved Re at the base of the Himalaya from the Yamuna are ∼150 mol at Batamandi and ∼100 mol at Saharanpur, compared to ∼120 mol from the Ganga at Rishikesh. The total flux from the Yamuna and the Ganga account for ∼0.4% of the global riverine Re flux, much higher than their contribution to global water discharge. This is also borne out from the mobilization rate of Re: ∼1 to 3 g km⁻² y⁻¹ in the Ganga and Yamuna basins in the Himalaya, compared to the global average of ∼0.1 g km⁻² y⁻¹.Black shale weathering can also significantly influence the budgets of Os and U in rivers and CO₂ in rivers and the atmosphere. Using dissolved Re in rivers as a proxy, it is estimated that ∼(6-9) × 10⁸ kg y⁻¹ of black shales are being weathered in the Ganga and Yamuna basins in the Himalaya. Weathering of such amounts of black shales can account for the reported concentrations of Os and U in these rivers. Furthermore, if the weathering results in the conversion of organic carbon in the black shales to CO₂, it would release ∼2 × 10⁵ mol of CO₂ km⁻² y⁻¹ in the Yamuna and Ganga basins in the Himalaya, comparable to the CO₂ consumption from silicate weathering.     

Sr and Sr/Sr in the Yamuna River System in the Himalaya: sources, fluxes, and controls on sr isotope composition

Sr and ⁸⁷Sr/⁸⁶Sr have been measured in the Yamuna river headwaters and many of its tributaries (YRS) in the Himalaya. These results, with those available for major ions in YRS rivers and in various lithologies of their basin, have been used to determine their contributions to riverine Sr and its isotopic budget. Sr in the YRS ranges from 120 to 13,400 nM, and ⁸⁷Sr/⁸⁶Sr from 0.7142 to 0.7932. Streams in the upper reaches, draining predominantly silicates, have low Sr and high ⁸⁷Sr/⁸⁶Sr whereas those draining the lower reaches exhibit the opposite resulting from differences in drainage lithology. ⁸⁷Sr/⁸⁶Sr shows significant co-variation with SiO₂/TDS and (Na^* + K)/TZ⁺ (indices of silicate weathering) in YRS waters, suggesting the dominant role of silicate weathering in contributing to high radiogenic Sr. This is also consistent with the observation that streams draining largely silicate terrains have the highest ⁸⁷Sr/⁸⁶Sr, analogous to that reported for the Ganga headwaters. Evaluation of the significance of other sources such as calc-silicates and trace calcites in regulating Sr budget of these rivers and their high ⁸⁷Sr/⁸⁶Sr needs detailed work on their Sr and ⁸⁷Sr/⁸⁶Sr. Preliminary calculations, however, indicate that they can be a significant source to some of the rivers.It is estimated that on an average, ∼25% of Sr in the YRS is derived from silicate weathering. In the lower reaches, the streams receive ∼15% of their Sr from carbonate weathering whereas in the upper reaches, calc-silicates can contribute significantly (∼50%) to the Sr budget of rivers. These calculations reveal the need for additional sources for rivers in the lower reaches to balance their Sr budget. Evaporites and phosphorites are potential candidates as judged from their occurrence in the drainage basin. In general, Precambrian carbonates, evaporites, and phosphorites “dilute” the high ⁸⁷Sr/⁸⁶Sr supplied by silicates, thus making Sr isotope distribution in YRS an overall two end member mixing. Major constraints in quantifying contributions of Sr and ⁸⁷Sr/⁸⁶Sr from different sources to YRS rivers are the wide range in Sr and ⁸⁷Sr/⁸⁶Sr of major lithologies, limited data on Sr and ⁸⁷Sr/⁸⁶Sr in minor phases and on the behavior of Sr, Na, and Ca during weathering and transport.The Ganga and the Yamuna together transport ∼0.1% of the global Sr flux at the foothills of the Himalaya which is in the same proportion as their contribution to global water discharge. Dissolved Sr flux from the Yamuna and its mobilization rate in the YRS basin is higher than those in the Ganga basin in the Himalaya, a result consistent with higher physical and chemical erosion rates in the YRS.     

Effect of groundwater pH and ionic strength on strontium sorption in aquifer sediments: Implications for Sr mobility at contaminated nuclear sites

Strontium-90 is a beta emitting radionuclide produced during nuclear fission, and is a problem contaminant at many nuclear facilities. Transport of ⁹⁰Sr in groundwaters is primarily controlled by sorption reactions with aquifer sediments. The extent of sorption is controlled by the geochemistry of the groundwater and sediment mineralogy. Here, batch sorption experiments were used to examine the sorption behaviour of ⁹⁰Sr in sediment–water systems representative of the UK Sellafield nuclear site based on groundwater and contaminant fluid compositions. In experiments with low ionic strength groundwaters (<0.01 mol L⁻¹), pH variation is the main control on sorption. The sorption edge for ⁹⁰Sr was observed between pH 4 and 6 with maximum sorption occurring (K_d ∼ 10³ L kg⁻¹) at pH 6–8. At ionic strengths above 10 mmol L⁻¹, and at pH values between 6 and 8, cation exchange processes reduced ⁹⁰Sr uptake to the sediment. This exchange process explains the lower ⁹⁰Sr sorption (K_d ∼ 40 L kg⁻¹) in the presence of artificial Magnox tank liquor (IS = 29 mmol L⁻¹). Strontium K-edge EXAFS spectra collected from sediments incubated with Sr²⁺ in either HCO₃-buffered groundwater or artificial Magnox tank liquor, revealed a coordination environment of ∼9 O atoms at 2.58–2.61 Å after 10 days. This is equivalent to the Sr²⁺ hydration sphere for the aqueous ion and indicates that Sr occurs primarily in outer sphere sorption complexes. No change was observed in the Sr sorption environment with EXAFS analysis after 365 days incubation. Sequential extractions performed on sediments after 365 days also found that ∼80% of solid associated ⁹⁰Sr was exchangeable with 1 M MgCl₂ in all experiments. These results suggest that over long periods, ⁹⁰Sr in contaminated sediments will remain primarily in weakly bound surface complexes. Therefore, if groundwater ionic strength increases (e.g. by saline intrusion related to sea level rise or by design during site remediation) then substantial remobilisation of ⁹⁰Sr is to be expected.     

Major ion chemistry of the Yarlung Tsangpo-Brahmaputra river: Chemical weathering, erosion, and CO2 consumption in the southern Tibetan plateau and eastern syntaxis of the Himalaya

The Yarlung Tsangpo-Brahmaputra river drains a large portion of the Himalaya and southern Tibetan plateau, including the eastern Himalayan syntaxis, one of the most tectonically active regions on the globe. We measured the solute chemistry of 161 streams and major tributaries of the Tsangpo-Brahmaputra to examine the effect of tectonic, climatic, and geologic factors on chemical weathering rates. Specifically, we quantify chemical weathering fluxes and CO₂ consumption by silicate weathering in southern Tibet and the eastern syntaxis of the Himalaya, examine the major chemical weathering reactions in the tributaries of the Tsangpo-Brahmaputra, and determine the total weathering flux from carbonate and silicate weathering processes in this region. We show that high precipitation, rapid tectonic uplift, steep channel slopes, and high stream power generate high rates of chemical weathering in the eastern syntaxis. The total dissolved solids (TDS) flux from the this area is greater than 520 tons km⁻² yr⁻¹ and the silicate cation flux more than 34 tons km⁻² yr⁻¹. In total, chemical weathering in this area consumes 15.2 × 10⁵ mol CO₂ km⁻² yr⁻¹, which is twice the Brahmaputra average. These data show that 15-20% of the total CO₂ consumption by silicate weathering in the Brahmaputra catchment is derived from only 4% of the total land area of the basin. Hot springs and evaporite weathering provide significant contributions to dissolved Na⁺ and Cl⁻ fluxes throughout southern Tibet, comprising more than 50% of all Na⁺ in some stream systems. Carbonate weathering generates 80-90% of all dissolved Ca²⁺ and Mg²⁺ cations in much of the Yarlung Tsangpo catchment.     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

Constraining the marine strontium budget with natural strontium isotope fractionations (Sr/Sr∗, δSr) of carbonates, hydrothermal solutions and river waters

We present strontium (Sr) isotope ratios that, unlike traditional ⁸⁷Sr/⁸⁶Sr data, are not normalized to a fixed ⁸⁸Sr/⁸⁶Sr ratio of 8.375209 (defined as δ^88/86Sr = 0 relative to NIST SRM 987). Instead, we correct for isotope fractionation during mass spectrometry with a ⁸⁷Sr-⁸⁴Sr double spike. This technique yields two independent ratios for ⁸⁷Sr/⁸⁶Sr and ⁸⁸Sr/⁸⁶Sr that are reported as (⁸⁷Sr/⁸⁶Sr∗) and (δ^88/86Sr), respectively. The difference between the traditional radiogenic (⁸⁷Sr/⁸⁶Sr normalized to ⁸⁸Sr/⁸⁶Sr = 8.375209) and the new ⁸⁷Sr/⁸⁶Sr∗ values reflect natural mass-dependent isotope fractionation. In order to constrain glacial/interglacial changes in the marine Sr budget we compare the isotope composition of modern seawater ((⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_Seawater) and modern marine biogenic carbonates ((⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_Carbonates) with the corresponding values of river waters ((⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_River) and hydrothermal solutions ((⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_HydEnd) in a triple isotope plot. The measured (⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_River values of selected rivers that together account for ∼18% of the global Sr discharge yield a Sr flux-weighted mean of (0.7114(8), 0.315(8)‰). The average (⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_HydEnd values for hydrothermal solutions from the Atlantic Ocean are (0.7045(5), 0.27(3)‰). In contrast, the (⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_Carbonates values representing the marine Sr output are (0.70926(2), 0.21(2)‰). We estimate the modern Sr isotope composition of the sources at (0.7106(8), 0.310(8)‰). The difference between the estimated (⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_input and (⁸⁷Sr/⁸⁶Sr∗, δ^88/86Sr)_output values reflects isotope disequilibrium with respect to Sr inputs and outputs. In contrast to the modern ocean, isotope equilibrium between inputs and outputs during the last glacial maximum (10-30 ka before present) can be explained by invoking three times higher Sr inputs from a uniquely “glacial” source: weathering of shelf carbonates exposed at low sea levels. Our data are also consistent with the “weathering peak” hypothesis that invokes enhanced Sr inputs resulting from weathering of post-glacial exposure of abundant fine-grained material.     

Silicon isotopes in meteorites and planetary core formation

The silicon (Si) isotope compositions of 42 meteorite and terrestrial samples have been determined using MC-ICPMS with the aim of resolving the current debate over their compositions and the implications for core formation. No systematic δ³⁰Si differences are resolved between chondrites (δ³⁰Si = −0.49 ± 0.15‰, 2σ_SD) and achondrites (δ³⁰Si = −0.47 ± 0.11‰, 2σ_SD), although enstatite chondrites are consistently lighter (δ³⁰Si = −0.63 ± 0.07‰, 2σ_SD) in comparison to other meteorite groups. The data reported here for meteorites and terrestrial samples display an average difference Δ³⁰Si_{BSE−meteorite∗} = 0.15 ± 0.10‰, which is consistent within uncertainty with previous studies. No effect from sample heterogeneity, preparation, chemistry or mass spectrometry can be identified as responsible for the reported differences between current datasets. The heavier composition of the bulk silicate Earth is consistent with previous conclusions that Si partitioned into the metal phase during metal-silicate equilibration at the time of core formation. Fixing the temperature of core formation to the peridotite liquidus and using an appropriate metal silicate fractionation factor (ε ∼0.89), the Δ³⁰Si_{BSE−meteorite∗} value from this study indicates that the Earth core contains at least 2.5 and possibly up to 16.8 wt% Si.     

Chemical erosion in the eastern Himalaya: Major ion composition of the Brahmaputra and δC of dissolved inorganic carbon

Major ion composition of waters, δ¹³C of its DIC (dissolved inorganic carbon), and the clay mineral composition of bank sediments in the Brahmaputra River System (draining India and Bangladesh) have been measured to understand chemical weathering and erosion and the factors controlling these processes in the eastern Himalaya. The time-series samples, collected biweekly at Guwahati, from the Brahmaputra mainstream, were also analyzed for the major ion composition. Clay mineralogy and chemical index of alteration (CIA) of sediments suggest that weathering intensity is relatively poor in comparison to that in the Ganga basin. This is attributed to higher runoff and associated physical erosion occurring in the Brahmaputra basin. The results of this study show, for the first time, spatial and temporal variations in chemical and silicate erosion rates in the Brahmaputra basin. The subbasins of the Brahmaputra watershed exhibit chemical erosion rates varying by about an order of magnitude. The Eastern Syntaxis basin dominates the erosion with a rate of ∼300 t km⁻² y⁻¹, one of the highest among the world river basins and comparable to those reported for some of the basaltic terrains. In contrast, the flat, cold, and relatively more arid Tibetan basin undergoes much slower chemical erosion (∼40 t km⁻² y⁻¹). The abundance of total dissolved solids (TDS, 102-203 mg/L) in the time-series samples collected over a period of one year shows variations in accordance with the annual discharge, except one of them, cause for which is attributable to flash floods. Na* (Na corrected for cyclic component) shows a strong positive correlation with Si, indicating their common source: silicate weathering. Estimates of silicate cations (Na_sil+K_sil+Ca_sil+Mg_sil) suggest that about half of the dissolved cations in the Brahmaputra are derived from silicates, a proportion higher than that for the Ganga system. The CO₂ consumption rate due to silicate weathering in the Brahmaputra watershed is ∼6 × 10⁵ moles km⁻² y⁻¹; whereas that in the Eastern Syntaxis subbasin is ∼19 × 10⁵ moles km⁻² y⁻¹, similar to the estimates for some of the basaltic terrains. This study suggests that the Eastern Syntaxis basin of the Brahmaputra is one of most intensely chemically eroding regions of the globe; and that runoff and physical erosion are the controlling factors of chemical erosion in the eastern Himalaya.     

Geochemical fluxes and weathering of volcanic terrains on high standing islands: Taranaki and Manawatu-Wanganui regions of New Zealand

Sediment fluxes from high standing oceanic islands (HSIs) such as New Zealand are some of the highest known [Milliman J. D. and Syvitski J. P. M. (1992) Geomorphic/tectonic control of sediment discharge to the ocean: the importance of small mountainous rivers. J. Geol.100, 525-544]. Recent geochemical work has suggested that along with their extremely high physical weathering yields, many New Zealand watersheds also have very high chemical weathering yields. In New Zealand, the magnitude of both the physical and chemical weathering yields is related to the lithology of the watershed. Most of the previous work on this topic has been undertaken in Southern Alps watersheds of schist and greywacke and in East Cape watersheds of semi-consolidated marine sediments and greywacke. We recently sampled North Island watersheds in the Taranaki and Manawatu-Wanganui regions which have been subjected to volcanism since the Miocene. We sampled watersheds that contain both volcanic and sedimentary rocks. A series of water and sediment samples was collected and analyzed for major, minor and trace elements. This was done to quantify the weathering intensities in the watersheds and to establish the relationship between physical and chemical weathering yields in volcanic lithologies. Our results reveal distinct chemical signatures for the different regions. Waters draining the Taranaki region volcanics are significantly enriched in K⁺, and depleted in Ca²⁺ and Sr²⁺ compared to waters draining the Manawatu-Wanganui region volcanics, which also traverse expanses of sedimentary siltstones and mudstones. The Ca²⁺ and Sr²⁺ depletions may reflect the relative absence of CaCO₃ in the Taranaki region watersheds. In addition, sediment samples from the Taranaki region show significant enrichment in Ti, Al, Ca, Fe, Mn, Mg, Ca, and P and depletion in Si and Rb compared to those of the Manawatu-Wanganui region. From total dissolved solids concentrations and mean annual water discharge, we calculate chemical weathering yields of 60-240 tons km⁻² a⁻¹. These weathering yields fall within the middle to upper range of those previously documented for the Southern Alps (93-480 tons km⁻² a⁻¹) and East Cape (62-400 tons km⁻² a⁻¹). Calculated silicate weathering yields of 12-33.6 tons km⁻² a⁻¹ and CO₂ consumption of 852-2390 × 10³ mol km⁻² a⁻¹ for the rivers draining the Taranaki volcanic region are higher than those previously reported for watersheds hosted in sedimentary and metamorphosed rock terrains on HSIs. CO₂ consumption is found to be within the range previously measured for the basaltic terrains of the Deccan Traps (580-2450 × 10³ mol km⁻² a⁻¹) and Réunion Island(1300-4400 × 10³ mol km⁻² a⁻¹). Our calculated chemical weathering yields demonstrate the importance of HSIs, particularly those with volcanic terrains, when considering global geochemical fluxes.     

Sr fluxes and Sr/Sr characterization of river waters from a karstic versus granitic watershed in the Yangtze River

The watershed in the southern Jiangxi Province (Jiangxi Province is called simply Gan) (SGW) and the watershed in the central Guizhou Province (Guizhou Province is called simply Qian) (CQW) are two subtropical watersheds of the Yangtze River in China. Both watersheds have similar latitudes and climate, but distinct differences in basin lithology. These similarities and differences provide a good natural laboratory in which to investigate weathering processes and Sr end-members in river waters. This work aims to identify and contrast the sources, fluxes and controls on Sr isotopic composition in the river waters of these two areas. Results showed that the ⁸⁷Sr/⁸⁶Sr in the SGW waters ranged from 0.716501 to 0.724931, with dissolved Sr averaging 27 μg l^− 1. Rhyolites and granites are two major sources for the dissolved Sr. The SGW waters receive 42% of their Sr from silicates weathering, 32% from carbonates and 3.2% from evaporites. ⁸⁷Sr/⁸⁶Sr in the CQW waters has a lesser variation from 0.707694 to 0.710039, but higher Sr contents (average of 208 μg l^− 1). Dolomite, limestone and dolomitic limestone are major sources of Sr in the waters. The CQW waters receive 69% of their Sr from carbonates, 1.7% from silicates and 0.9% from evaporites. The chemical erosion rate and Sr flux in the CQW are 122 t km^− 2 a^− 1 and 0.079 t km^− 2 a^− 1, respectively, which are higher than those of the SGW (56 t km^− 2 a^− 1 and 0.021 t km^− 2 a^− 1, respectively). These data suggest that the intensive carbonates weathering occurred in the karstic area in the upper-reach of the Yangtze River exert great influence on the high Sr concentration and low Sr isotopic ratios in the River.     

Magnesium isotope systematics of the lithologically varied Moselle river basin, France

Magnesium and strontium isotope signatures were determined during different seasons for the main rivers of the Moselle basin, northeastern France. This small basin is remarkable for its well-constrained and varied lithology on a small distance scale, and this is reflected in river water Sr isotope compositions. Upstream, where the Moselle River drains silicate rocks of the Vosges mountains, waters are characterized by relatively high ⁸⁷Sr/⁸⁶Sr ratios (0.7128-0.7174). In contrast, downstream of the city of Epinal where the Moselle River flows through carbonates and evaporites of the Lorraine plateau, ⁸⁷Sr/⁸⁶Sr ratios are lower, down to 0.70824.Magnesium in river waters draining silicates is systematically depleted in heavy isotopes (δ²⁶Mg values range from −1.2 to −0.7‰) relative to the value presently estimated for the continental crust and a local diorite (−0.5‰). In comparison, δ²⁶Mg values measured in soil samples are higher (∼0.0‰). This suggests that Mg isotope fractionation occurs during mineral leaching and/or formation of secondary clay minerals. On the Lorraine plateau, tributaries draining marls, carbonates and evaporites are characterized by low Ca/Mg (1.5-3.2) and low Ca/Sr (80-400) when compared to local carbonate rocks (Ca/Mg = 29-59; Ca/Sr = 370-2200), similar to other rivers draining carbonates. The most likely cause of the Mg and Sr excesses in these rivers is early thermodynamic saturation of groundwater with calcite relative to magnesite and strontianite as groundwater chemistry progressively evolves in the aquifer. δ²⁶Mg of the dissolved phases of tributaries draining mainly carbonates and evaporites are relatively low and constant throughout the year (from −1.4‰ to −1.6‰ and from −1.2‰ to −1.4‰, respectively), within the range defined for the underlying rocks. Downstream of Epinal, the compositions of the Moselle River samples in a δ²⁶Mg vs. ⁸⁷Sr/⁸⁶Sr diagram can be explained by mixing curves between silicate, carbonate and evaporite waters, with a significant contribution from the Vosgian silicate lithologies (>70%). Temporal co-variation between δ²⁶Mg and ⁸⁷Sr/⁸⁶Sr for the Moselle River throughout year is also observed, and is consistent with a higher contribution from the Vosges mountains in winter, in terms of runoff and dissolved element flux. Overall, this study shows that Mg isotopes measured in waters, rocks and soils, coupled with other tracers such as Sr isotopes, could be used to better constrain riverine Mg sources, particularly if analytical uncertainties in Mg isotope measurements can be improved in order to perform more precise quantifications.     

Sr and Nd isotopes as tracers of clastic sources in Lake Le Bourget sediment (NW Alps, France) during the Little Ice Age: Palaeohydrology implications

Geochemical methods (major elements and Sr, Nd isotopes) have been used to (1) characterize Lake Le Bourget sediments in the French Alps, (2) identify the current sources of the clastic sediments and estimate the source variability over the last 600 years. Major element results indicate that Lake Le Bourget sediments consist of 45% clastic component and 55% endogenic calcite. In addition, several individual flood levels have been identified during the Little Ice Age (LIA) on the basis of their higher clastic content (> 70%).Potential sources of Lake Le Bourget clastic sediments have been investigated from Sr and Nd isotope compositions. The sediments from the Sierroz River and Leysse River which are mainly derived from the Mesozoic Calcareous Massifs are characterised by lower ⁸⁷Sr/⁸⁶Sr ratios and slightly lower ?_Nd(0) ratios than the Arve River sediments which are derived from the Palaeozoic Mont-Blanc External Crystalline Massifs. The Rhône River appears to have been the main source of clastic sediments into the lake for the last 600 years, as evidenced by a similar Sr and Nd isotopic compositions analyzed in core B16 sediments (⁸⁷Sr/⁸⁶Sr = 0.719, ?_Nd(0) = − 10) and in the sediments of the Rhône River (⁸⁷Sr/⁸⁶Sr = 0.719, ?_Nd(0) = − 9.6).The isotopic signatures of flood events and background samples from core B16 in Lake Le Bourget are also similar. This indicates that prior to ∼ 1800, the inputs into the lake have remained relatively homogeneous with the proportion of clastic component mainly being a function of the palaeohydrology of the Rhone River. Early human modification (deforestation and agriculture) of the lake catchment before the 1800s appears to have had little influence on the source of clastic sediments.     

Oxygen- and magnesium-isotope compositions of calcium-aluminum-rich inclusions from Rumuruti (R) chondrites

We report oxygen- and magnesium-isotope compositions of Ca,Al-rich inclusions (CAIs) from several Rumuruti (R) chondrites measured in situ using a Cameca ims-1280 ion microprobe. On a three-isotope oxygen diagram, δ¹⁷O vs. δ¹⁸O, compositions of individual minerals in most R CAIs analyzed fall along a slope-1 line. Based on the variations of Δ¹⁷O values (Δ¹⁷O = δ¹⁷O − 0.52 × δ¹⁸O) within individual inclusions, the R CAIs are divided into (i) ¹⁶O-rich (Δ¹⁷O ∼ −23-26‰), (ii) uniformly ¹⁶O-depleted (Δ¹⁷O ∼ −2‰), and (iii) isotopically heterogeneous (Δ¹⁷O ranges from −25‰ to +5‰). One of the hibonite-rich CAIs, H030/L, has an intermediate Δ¹⁷O value of −12‰ and a highly fractionated composition (δ¹⁸O ∼ +47‰). We infer that like most CAIs in other chondrite groups, the R CAIs formed in an ¹⁶O-rich gaseous reservoir. The uniformly ¹⁶O-depleted and isotopically heterogeneous CAIs subsequently experienced oxygen-isotope exchange during remelting in an ¹⁶O-depleted nebular gas, possibly during R chondrite chondrule formation, and/or during fluid-assisted thermal metamorphism on the R chondrite parent asteroid.Three hibonite-bearing CAIs and one spinel-plagioclase-rich inclusion were analyzed for magnesium-isotope compositions. The CAI with the highly fractionated oxygen isotopes, H030/L, shows a resolvable excess of ²⁶Mg (²⁶Mg^∗) corresponding to an initial ²⁶Al/²⁷Al ratio of ∼7 × 10⁻⁷. Three other CAIs show no resolvable excess of ²⁶Mg (²⁶Mg^∗). The absence of ²⁶Mg^∗ in the spinel-plagioclase-rich CAI from a metamorphosed R chondrite NWA 753 (R3.9) could have resulted from metamorphic resetting. Two other hibonite-bearing CAIs occur in the R chondrites (NWA 1476 and NWA 2446), which appear to have experienced only minor degrees of thermal metamorphism. These inclusions could have formed from precursors with lower than canonical ²⁶Al/²⁷Al ratio.     

Chemical weathering of a marine terrace chronosequence, Santa Cruz, California. Part II: Solute profiles, gradients and the comparisons of contemporary and long-term weathering rates

The spatial and temporal changes in hydrology and pore water elemental and ⁸⁷Sr/⁸⁶Sr compositions are used to determine contemporary weathering rates in a 65- to 226-kyr-old soil chronosequence formed from granitic sediments deposited on marine terraces along coastal California. Soil moisture, tension and saturation exhibit large seasonal variations in shallow soils in response to a Mediterranean climate. These climate effects are dampened in underlying argillic horizons that progressively developed in older soils, and reached steady-state conditions in unsaturated horizons extending to depths in excess of 15 m. Hydraulic fluxes (q_h), based on Cl mass balances, vary from 0.06 to 0.22 m yr⁻¹, resulting in fluid residence times in the terraces of 10-24 yrs.As expected for a coastal environment, the order of cation abundances in soil pore waters is comparable to sea water, i.e., Na > Mg > Ca > K > Sr, while the anion sequence Cl > NO₃ > HCO₃ > SO₄ reflects modifying effects of nutrient cycling in the grassland vegetation. Net Cl-corrected solute Na, K and Si increase with depth, denoting inputs from feldspar weathering. Solute ⁸⁷Sr/⁸⁶Sr ratios exhibit progressive mixing of sea water-dominated precipitation with inputs from less radiogenic plagioclase. While net Sr and Ca concentrations are anomalously high in shallow soils due to biological cycling, they decline with depth to low and/or negative net concentrations. Ca/Mg, Sr/Mg and ⁸⁷Sr/⁸⁶Sr solute and exchange ratios are similar in all the terraces, denoting active exchange equilibration with selectivities close to unity for both detrital smectite and secondary kaolinite. Large differences in the magnitudes of the pore waters and exchange reservoirs result in short-term buffering of the solute Ca, Sr, and Mg. Such buffering over geologic time scales can not be sustained due to declining inputs from residual plagioclase and smectite, implying periodic resetting of the exchange reservoir such as by past vegetational changes and/or climate.Pore waters approach thermodynamic saturation with respect to albite at depth in the younger terraces, indicating that weathering rates ultimately become transport-limited and dependent on hydrologic flux. Contemporary rates R_solute are estimated from linear Na and Si pore weathering gradients b_solute such that

     

Competitive adsorption of strontium and fulvic acid at the muscovite-solution interface observed with resonant anomalous X-ray reflectivity

Molecular-scale distributions of Sr²⁺ and fulvic acid (FA) adsorbed on the muscovite (0 0 1) surface were investigated using in situ specular X-ray reflectivity (XR) and resonant anomalous X-ray reflectivity (RAXR). The total amount of Sr²⁺ adsorbed from a 1 × 10⁻² mol/kg SrCl₂ and 100 mg/kg Elliott Soil Fulvic Acid II (ESFA II) solution at pH 5.5 compensated 81 ± 5% of the muscovite surface charge, less than previously measured (118 ± 5%) in an ESFA II-free solution with the same Sr concentration and pH. Inner-sphere (IS) and outer-sphere (OS) Sr²⁺ constituted 87% of the total adsorbed species in IS:OS proportions of 19:81 compared to 42:58 in the solution without FA, suggesting that adsorbed FA competes with the IS Sr²⁺ for surface sites. The coverage of both IS and OS Sr²⁺ decreased even more in a pH 3.5 solution containing the same concentration of FA and 0.5 × 10⁻² mol/kg Sr(NO₃)₂, whereas a significant amount of Sr²⁺ accumulated farther from the surface in the FA layer. The amount of Sr²⁺ incorporated in the ∼10 Å thick FA layer decreased by 79% with decreasing FA concentration (100 → 1 mg/kg) and increasing Sr²⁺ concentration (0.5 × 10⁻² → 1 × 10⁻² mol/kg) and pH (3.5 → 3.6). These results indicate not only that adsorbed FA molecules (and perhaps also H₃O⁺) displace Sr²⁺ near the muscovite surface, but also that the sorbed FA film provides binding sites for additional Sr²⁺ away from the surface. When a muscovite crystal pre-coated with FA after reaction in a 100 mg/kg ESFA II solution for 50 h was subsequently reacted with a 0.5 × 10⁻² mol/kg Sr(NO₃)₂ and 100 mg/kg ESFA II solution at pH 3.7, a significant fraction of Sr²⁺ was distributed in the outer part of the FA film similar to that observed on fresh muscovite reacted at pH 3.5 with a pre-mixed Sr-FA solution at the same concentrations. However, this Sr²⁺ sorbed in the pre-adsorbed organic film was more widely distributed and had a lower coverage, suggesting that pre-sorbed FA may undergo fractionation and/or conformational changes that diminish its capacity, and that of the muscovite (0 0 1) surface, for adsorbing the aqueous Sr cation.     

Seasonal variations of physical and chemical erosion: A three-year survey of the Rhone River (France)

Numerous studies of weathering fluxes have been carried out on major world rivers during the last decade, to estimate CO₂ consumption rates, landscape evolution and global erosion rates. For obvious logistical reasons, most of these studies were based on large scale investigations carried out on short timescales. By comparison, much less effort has been devoted to long term monitoring, as a means to verify the temporal variability of the average characteristics, their trends, and the representativeness of short-term investigations. Here we report the results of a three-year survey (November 2000 to December 2003) of the major and trace element composition of dissolved and suspended matter in the lower Rhone River (France), the largest river of the Mediterranean area. Subsurface water samples were collected in Arles, about 48 km upstream of the estuary, twice a month routinely, and at higher frequency during flood events.During each flood event, the suspended particulate matter (SPM) show the usual trend of clockwise hysteresis with higher SPM concentrations on the rising limb of the flood than at the same discharge on the falling limb. We show that the annual average SPM flux of the Rhone River to the Mediterranean Sea (7.3 ± 0.6 × 10⁶ tons yr⁻¹) was largely controlled by the flood events (83% of the solid discharge occurred in less than 12% of the time), and that the precision on the total output flux depends strongly on the precise monitoring of SPM variations during the floods.The chemical composition of water and SPM are characterized by the predominance of Ca²⁺ due to the abundance of carbonate rocks in the Rhone watershed. Chemical budgets have been calculated to derive the contributions of atmospheric deposition, carbonate, silicate and evaporite weathering, and anthropogenic inputs. The chemical weathering rate of carbonates is estimated to be 89 ± 5 t km⁻² yr⁻¹ compared to 14.4 ± 3 t km⁻² yr⁻¹ from silicates. By contrast, the physical erosion rate of silicates is about 51 t km⁻² yr⁻¹ against 19 t km⁻² yr⁻¹ for carbonates.The steady-state model of Gaillardet et al. (1995) has been applied to the chemical composition of dissolved and solid products. The results show that the Rhone River currently exports much less material than produced at steady-state by weathering in its watershed. The sediment flux inferred from the steady-state calculation (21-56 × 10⁶ t yr⁻¹) is on the same order as that estimated in literature for the 19th and the beginning of the 20th centuries. This imbalance may suggest that the Rhone is under a transient erosion regime following climate change (i.e. significant decrease of the flooding frequency since the beginning of the 19th century). On the other hand, the imbalance may also be due to the trapping of alluvion by the numerous dams on the river and its tributaries.Our data corroborate with previous studies that suggest a strong coupling between chemical and physical erosion fluxes, during the hydrological seasonal cycle of the Rhone River. The correlation between physical and chemical transport rates is, however, clearly different from that reported for global annual averages in large world rivers.