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.     

Regolith mass balance inferred from combined mineralogical, geochemical and geophysical studies: Mule Hole gneissic watershed, South India

The aim of this study is to propose a method to assess the long-term chemical weathering mass balance for a regolith developed on a heterogeneous silicate substratum at the small experimental watershed scale by adopting a combined approach of geophysics, geochemistry and mineralogy. We initiated in 2003 a study of the steep climatic gradient and associated geomorphologic features of the edge of the rifted continental passive margin of the Karnataka Plateau, Peninsular India. In the transition sub-humid zone of this climatic gradient we have studied the pristine forested small watershed of Mule Hole (4.3 km²) mainly developed on gneissic substratum. Mineralogical, geochemical and geophysical investigations were carried out (i) in characteristic red soil profiles and (ii) in boreholes up to 60 m deep in order to take into account the effect of the weathering mantle roots. In addition, 12 Electrical Resistivity Tomography profiles (ERT), with an investigation depth of 30 m, were generated at the watershed scale to spatially characterize the information gathered in boreholes and soil profiles. The location of the ERT profiles is based on a previous electromagnetic survey, with an investigation depth of about 6 m. The soil cover thickness was inferred from the electromagnetic survey combined with a geological/pedological survey.Taking into account the parent rock heterogeneity, the degree of weathering of each of the regolith samples has been defined using both the mineralogical composition and the geochemical indices (Loss on Ignition, Weathering Index of Parker, Chemical Index of Alteration). Comparing these indices with electrical resistivity logs, it has been found that a value of 400 Ohm m delineates clearly the parent rocks and the weathered materials. Then the 12 inverted ERT profiles were constrained with this value after verifying the uncertainty due to the inversion procedure. Synthetic models based on the field data were used for this purpose. The estimated average regolith thickness at the watershed scale is 17.2 m, including 15.2 m of saprolite and 2 m of soil cover.Finally, using these estimations of the thicknesses, the long-term mass balance is calculated for the average gneiss-derived saprolite and red soil. In the saprolite, the open-system mass-transport function τ indicates that all the major elements except Ca are depleted. The chlorite and biotite crystals, the chief sources for Mg (95%), Fe (84%), Mn (86%) and K (57%, biotite only), are the first to undergo weathering and the oligoclase crystals are relatively intact within the saprolite with a loss of only 18%. The Ca accumulation can be attributed to the precipitation of CaCO₃ from the percolating solution due to the current and/or the paleoclimatic conditions. Overall, the most important losses occur for Si, Mg and Na with −286 × 10⁶ mol/ha (62% of the total mass loss), −67 × 10⁶ mol/ha (15% of the total mass loss) and −39 × 10⁶ mol/ha (9% of the total mass loss), respectively. Al, Fe and K account for 7%, 4% and 3% of the total mass loss, respectively. In the red soil profiles, the open-system mass-transport functions point out that all major elements except Mn are depleted. Most of the oligoclase crystals have broken down with a loss of 90%. The most important losses occur for Si, Na and Mg with −55 × 10⁶ mol/ha (47% of the total mass loss), −22 × 10⁶ mol/ha (19% of the total mass loss) and −16 × 10⁶ mol/ha (14% of the total mass loss), respectively. Ca, Al, K and Fe account for 8%, 6%, 4% and 2% of the total mass loss, respectively.Overall these findings confirm the immaturity of the saprolite at the watershed scale. The soil profiles are more evolved than saprolite but still contain primary minerals that can further undergo weathering and hence consume atmospheric CO₂.     

Mg isotope constraints on soil pore-fluid chemistry: Evidence from Santa Cruz, California

Mg isotope ratios (²⁶Mg/²⁴Mg) are reported in soil pore-fluids, rain and seawater, grass and smectite from a 90 kyr old soil, developed on an uplifted marine terrace from Santa Cruz, California. Rain water has an invariant ²⁶Mg/²⁴Mg ratio (expressed as δ²⁶Mg) at −0.79 ± 0.05‰, identical to seawater δ²⁶Mg. Detrital smectite (from the base of the soil profile, and therefore unweathered) has a δ²⁶Mg value of 0.11‰, potentially enriched in ²⁶Mg by up to 0.3‰ compared to the bulk silicate Earth Mg isotope composition (although within the range of all terrestrial silicates). The soil pore-waters show a continuous profile with depth for δ²⁶Mg, ranging from −0.99‰ near the surface to −0.43‰ at the base of the profile. Shallow pore-waters (<1 m) have δ²⁶Mg values that are similar to, or slightly lower than the rain waters. This implies that the degree of biological cycling of Mg in the pore-waters is relatively small and is quantified as <32%, calculated using the average Mg isotope enrichment factor between grass and rain (δ²⁶Mg_grass-δ²⁶Mg_rain) of 0.21‰. The deep pore-waters (1-15 m deep) have δ²⁶Mg values that are intermediate between the smectite and rain, ranging from −0.76‰ to −0.43‰, and show a similar trend with depth compared to Sr isotope ratios. The similarity between Sr and Mg isotope ratios confirms that the Mg in the pore-waters can be explained by a mixture between rain and smectite derived Mg, despite the fact that Mg and Sr concentrations may be buffered by the exchangeable reservoir. However, whilst Sr isotope ratios in the pore-waters span almost the complete range between mineral and rain inputs, Mg isotopes compositions are much closer to the rain inputs. If Mg and Sr isotope ratios are controlled uniquely by a mixture, the data can be used to estimate the mineral weathering inputs to the pore-waters, by correcting for the rain inputs. This isotopic correction is compared to the commonly used chloride correction for precipitation inputs. A consistent interpretation is only possible if Mg isotope ratios are fractionated either by the precipitation of a secondary Mg bearing phase, not detected by conventional methods, or selective leaching of ²⁴Mg from smectite. There is therefore dual control on the Mg isotopic composition of the pore-waters, mixing of two inputs with distinct isotopic compositions, modified by fractionation. The data provide (1) further evidence for Mg isotope fractionation at the surface of the Earth and (2) the first field evidence of Mg isotope fractionation during uptake by natural plants. The coherent behaviour of Mg isotope ratios in soil environments is encouraging for the development of Mg isotope ratios as a quantitative tracer of both weathering inputs of Mg to waters, and the physicochemical processes that cycle Mg, a major cation linked to the carbon cycle, during continental weathering.     

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

     

Silicate weathering of soil-mantled slopes in an active Alpine landscape

Despite being located on high, steep, actively uplifting, and formerly glaciated slopes of the Swiss Central Alps, soils in the upper Rhone Valley are depleted by up to 50% in cations relative to their parent bedrock. This depletion was determined by a mass loss balance based on Zr as a refractory element. Both Holocene weathering rates and physical erosion rates of these slopes are unexpectedly low, as measured by cosmogenic ¹⁰Be-derived denudation rates. Chemical depletion fractions, CDF, range from 0.12 to 0.48, while the average soil chemical weathering rate is 33 ± 15 t km⁻² yr⁻¹. Both the cosmogenic nuclide-derived denudation rates and model calculations suggest that these soils have reached a weathering steady-state since deglaciation 15 ky ago. The weathering signal varies with elevation and hillslope morphology. In addition, the chemical weathering rates decrease with elevation indicating that temperature may be a dominant controlling factor on weathering in these high Alpine basins. Model calculations suggest that chemical weathering rates are limited by reaction kinetics and not the supply rate of fresh material. We compare hillslope and catchment-wide weathering fluxes with modern stream cation flux, and show that high relief, bare-rock slopes exhibit much lower chemical weathering rates despite higher physical erosion rates. The low weathering fluxes from rocky, rapidly eroding slopes allow for the broader implication that mountain building, while elevating overall denudation rates, may not cause increased chemical weathering rates on hillslopes. In order for this sediment to be weathered, intermediate storage, for instance in floodplains, is required.     

Hafnium and neodymium isotopes in surface waters of the eastern Atlantic Ocean: Implications for sources and inputs of trace metals to the ocean

We present hafnium (Hf) and neodymium (Nd) isotopic compositions and concentrations in surface waters of the eastern Atlantic Ocean between the coast of Spain and South-Africa. These data are complemented by Hf and Nd isotopic and concentration data, as well as rare earth element (REE) concentrations, in Saharan dust.Hafnium concentrations range between a maximum of 0.52 pmol/kg in the area of the Canary Islands and a minimum value of 0.08 pmol/kg in the southern Angola Basin. Neodymium concentrations also show a local maximum in the area of the Canary Islands (26 pmol/kg) but are even higher between ∼20°N and ∼4°N reaching maximum concentrations of 35 pmol/kg. These elevated concentrations provide evidence of inputs from weathering of the Canary Islands and from the partial dissolution of dust from the Sahara/Sahel region. The inputs from ocean island weathering are also reflected in radiogenic Hf and Nd isotopes.The Hf isotopic compositions of dust samples themselves are highly variable, ranging between ε_Hf = −20 and −0.6. The combined Hf and Nd isotopic compositions of dust plot close to the “terrestrial array” during periods of appreciable dust load in the atmosphere. During low atmospheric dust loading combined Hf and Nd isotopic compositions similar to seawater are observed. Most of the variability can be explained in terms of variable degrees of zircon loss from the dust samples, which in turn is linked to sorting during atmospheric transport to the eastern Atlantic Ocean and possibly presorting by sedimentary redistribution on the continent. In addition, increasing relative proportions of radiogenic clay minerals with decreasing grain size may contribute to the radiogenic Hf isotopic compositions observed.While the Nd isotopic composition in the surface ocean reflects the Nd isotopic composition of the Saharan dust adjacent to the Sahara/Sahel region, the release of Hf from that dust appears to be incongruent and results in surface ocean Hf isotopic compositions which are ∼10 ε_Hf more radiogenic than the bulk dust. Radiogenic Hf appears to be released from clays and possibly from trace apatite. Rare earth element patterns of dust samples indicate the presence of apatite but provide no evidence for ferromanganese grain coatings, suggesting that such coatings are insignificant in the release of Hf and Nd from Saharan dust to the surface ocean.The Nd isotopic composition of the surface waters becomes less radiogenic south of the equator, most likely reflecting the release of Nd from Congo river sediments. The release of Hf from Saharan dust and the Congo river sediments, however, does not produce distinct Hf isotopic signatures in the surface ocean, implying that the mobile fraction of Hf integrated over large continental areas is isotopically uniform. The Hf isotopic uniformity in the surface ocean means that the limited variability in deep water isotopic compositions is consistent with a short deep water residence time and reflects homogenous continental inputs rather than efficient deep water homogenization.     

Mineral weathering and elemental transport during hillslope evolution at the Susquehanna/Shale Hills Critical Zone Observatory

Located in the uplands of the Valley and Ridge physiographic province of Pennsylvania, the Susquehanna/Shale Hills Critical Zone Observatory (SSHO) is a tectonically quiescent, first-order catchment developed on shales of the Silurian Rose Hill Formation. We used soil cores augered at the highest point of the watershed and along a subsurface water flowline on a planar hillslope to investigate mineral transformations and physical/chemical weathering fluxes. About 25 m of bedrock was also drilled to estimate parent composition. Depletion of carbonate at tens of meters of depth in bedrock may delineate a deep carbonate-weathering front. Overlying this, extending from ∼6 m below the bedrock-soil interface up into the soil, is the feldspar dissolution front. In the soils, depletion profiles for K, Mg, Si, Fe, and Al relative to the bedrock define the illite and chlorite reaction fronts. When combined with a cosmogenic nuclide-derived erosion rate on watershed sediments, these depletion profiles are consistent with dissolution rates that are several orders of magnitudes slower for chlorite (1-5 × 10⁻¹⁷ mol m⁻² s⁻¹) and illite (2-9 × 10⁻¹⁷ mol m⁻² s⁻¹) than observed in the laboratory. Mineral reactions result in formation of vermiculite, hydroxy-interlayered vermiculite, and minor kaolinite. During weathering, exchangeable divalent cations are replaced by Al as soil pH decreases.The losses of Mg and K in the soils occur largely as solute fluxes; in contrast, losses of Al and Fe are mostly as downslope transport of fine particles. Physical erosion of bulk soils also occurs: results from a steady-state model demonstrate that physical erosion accounts for about half of the total denudation at the ridgetop and midslope positions. Chemical weathering losses of Mg, Na, and K are higher in the upslope positions likely because of the higher degree of chemical undersaturation in porewaters. Chemical weathering slows down in the valley floor and Al and Si even show net accumulation. The simplest model for the hillslope that is consistent with all observations is a steady-state, clay weathering-limited system where soil production rates decrease with increasing soil thickness.     

Stream geochemistry, chemical weathering and CO2 consumption potential of andesitic terrains, Dominica, Lesser Antilles

Recent studies of chemical weathering of andesitic-dacitic material on high-standing islands (HSIs) have shown these terrains have some of the highest observed rates of chemical weathering and associated CO₂ consumption yet reported. However, the paucity of stream gauge data in many of these terrains has limited determination of chemical weathering product fluxes. In July 2006 and March 2008, stream water samples were collected and manual stream gauging was performed in watersheds throughout the volcanic island of Dominica in the Lesser Antilles. Distinct wet and dry season solute concentrations reveal the importance of seasonal variations on the weathering signal. A cluster analysis of the stream geochemical data shows the importance of parent material age on the overall delivery of solutes. Observed Ca:Na, HCO₃:Na and Mg:Na ratios suggest crystallinity of the parent material may also play an important role in determining weathering fluxes. From total dissolved solids concentrations and mean annual discharge calculations we calculate chemical weathering yields of (6-106 t km⁻² a⁻¹), which are similar to those previously determined for basalt terrains. Silicate fluxes (3.1-55.4 t km⁻² a⁻¹) and associated CO₂ consumption (190-1575 × 10³ mol km⁻² a⁻¹) determined from our study are among the highest determined to date. The calculated chemical fluxes from our study confirm the weathering potential of andesitic-dacitic terrains and that additional studies of these terrains are warranted.     

Effects of temperature on silicate weathering: Solute fluxes and chemical weathering in a temperate rain forest watershed,Jamieson Creek,British Columbia

Chemical weathering of silicate minerals has long been known as a sink for atmospheric CO₂, and feedbacks between weathering and climate are believed to affect global climate. While warmer temperatures are believed to increase rates of weathering, weathering in cool climates can be accelerated by increased mineral exposure due to mechanical weathering by ice. In this study, chemical weathering of silicate minerals is investigated in a small temperate watershed. The Jamieson Creek watershed is covered by mature coniferous forest and receives high annual precipitation (4000 mm), mostly in the form of rainfall, and is underlain by quartz diorite bedrock and glacial till. Analysis of pore water concentration gradients indicates that weathering in hydraulically unsaturated ablation till is dominated by dissolution of plagioclase and hornblende. However, a watershed scale solute mass balance indicates high relative fluxes of K and Ca, indicating preferential leaching of these solutes possibly from the relatively unweathered lodgement till. Weathering rates for plagioclase and hornblende calculated from a watershed scale solute mass balance are similar in magnitude to rates determined using pore water concentration gradients.When compared to the Rio Icacos basin in Puerto Rico, a pristine tropical watershed with similar annual precipitation and bedrock, but with dissimilar regolith properties, fluxes of weathering products in stream discharge from the warmer site are 1.8 to 16.2-fold higher, respectively, and regolith profile-averaged plagioclase weathering rates are 3.8 to 9.0-fold higher. This suggests that the Arrhenius effect, which predicts a 3.5- to 9-fold increase in the dissolution rate of plagioclase as temperature is increased from 3.4° to 22 °C, may explain the greater weathering fluxes and rates at the Rio Icacos site. However, more modest differences in K and Ca fluxes between the two sites are attributed to accelerated leaching of those solutes from glacial till at Jamieson Creek. Our findings suggest that under conditions of high rainfall and favorable topography, weathering rates of silicate minerals in warm tropical systems will tend to be higher than in cool temperate systems, even if the temperate system is has been perturbed by an episode of glaciation that deposits regolith high in fresh mineral surface area.     

Natural and anthropogenic budgets of a small watershed in the massif central (France): Chemical and strontium isotopic characterization of water and sediments

A small watershed (160 km²) located in the Massif Central (France) has been chemically, isotopically and hydrologically studied through its dissolved load, bed sediments and soils. This watershed is underlain by basaltic bedrock and associated soils in which the vegetation is dominated mainly by meadows.Dissolved concentrations of major ions (Cl, SO₄, NO₃, HCO₃, Ca, Na, Mg, K, Al and Si), trace elements (Rb and Sr) and strontium isotopes have been determined for two different hydrologic periods on the main stream of the Allanche river and its tributaries.The major objectives of this study were to characterize the chemical and isotopic signatures of each reservoir occurring in the watershed. Changes in chemical and isotopic signatures are interpreted in terms of fluctuations of the different components inputs: rainwater, weathering products, anthropogenic addition.Water quality may be influenced by natural inputs (rainwater, weathering processes) and anthropogenic additions (fertilizers, road salts, etc.). Precipitation serves as a major vehicle for dissolved chemical species in addition to the hydrosystem and, in order to constrain rain inputs, a systematic study of rainwaters is carried out over a one year period using an automatic collector. Corrections of rainwater addition using chloride as an atmospheric input reference were computed for selected elements and the Sr/Sr ratio. After such corrections, the geochemical budget of the watershed was determined and the role of anthropogenic additions evaluated through the relationship between strontium isotopes and major and trace element ratios. Thus, 10% of Ca and Na originate in rainwater input, 40 to 80% in fertilizer additions and 15 to 50% in rock weatheringThe cationic denudation rates for this watershed are around 0.3 g s^–1 km² during low water discharge and 0.6 g s^–1 km² in high water stage. This led to a chemical denudation rate of 5.3 mm/1000 years.For solid matter, the normalization of chemical species relative to parent rocks shows the depletion or enrichment in soils and sediments. The use of K and Ca as mobile reference illustrates the weathering state of soils and sediments relative to parent rocks. This weathering state for bed sediments range from 15 to 45% for the K normalization and from 2 to 50% for the Ca normalization. For the soils, the weathering state ranges from 15 to 57% for the K normalization and from 17 to 90% for the Ca normalization.     

Whole-watershed mercury balance at Sagehen Creek, Sierra Nevada, CA

Few data are available on mercury (Hg) dynamics at high-elevation mountain sites. In this project, a whole-watershed approach was used to quantify major fluxes and pools of Hg in Sagehen basin, a closed basin in the Sierra Nevada mountains in California. Over a period spanning 9 months (January-September 2009), we estimated wet deposition inputs to the watershed at 3.8 μg m⁻². Dry deposition added additional Hg in the range of 0.30-2.45 μg m⁻² during this time period, and was the dominant deposition process during summer time. Seasonal snowpack accounted for only half of the Hg deposited by wet deposition. We suggest that photo-induced reduction of Hg(II) in snow and subsequent volatilization was responsible for this loss. Thus, snowpacks in the Sierra Nevada mountains likely reduce the effective atmospheric mercury flux via wet deposition due to significant emission fluxes prior to snowmelt. As such, wet Hg deposition could be of lesser importance as a Hg source in snow-dominated systems. Finally, stream runoff collected at the outlet of the watershed could account for only 4% of total Hg wet deposition suggesting that a large fraction of mercury deposition was sequestered in the ecosystem, specifically in the soils.     

Carbonate and silicate weathering in two presently glaciated, crystalline catchments in the Swiss Alps

We present a weathering mass balance of the presently glaciated Rhône and Oberaar catchments, located within the crystalline Aar massif (central Switzerland). Annual chemical and physical weathering fluxes are calculated from the monthly weighted means of meltwater samples taken from July, 1999 to May, 2001 and are corrected for precipitation inputs. The meltwater composition issuing from the Oberaar and Rhône catchments is dominated by calcium, which represents 81% and 55% of the total cation flux respectively (i.e. 555 and 82-96 keq km⁻² yr⁻¹). The six to seven times higher Ca²⁺ denudation flux from the Oberaar catchment is attributed to the presence of a strongly foliated gneissic zone. The gneissic zone has an elevated calcite content (as reflected by the 4.6 times higher calcite content of the suspended sediments from Oberaar compared to Rhône) and a higher mechanical erosion rate (resulting in a higher flux of suspended sediment). The mean flux of suspended calcite of the Oberaar meltwaters during the ablation period is 7 times greater than that of the Rhône meltwaters. Taking the suspended calcite as a proxy for the total (including sub-glacial sediments) weathering calcite surface area, it appears that the available surface area is an important factor in controlling weathering rates. However, we also observe an increased supply of protons for carbonate dissolution in the Oberaar catchment, where the sulphate denudation flux is six times greater. Carbonic acid is the second important source of protons, and we calculate that three times as much atmospheric CO₂ is drawn down (short term) in the Oberaar catchment. Silica fluxes from the two catchments are comparable with each other, but are 100 kmol km² yr⁻¹ lower than fluxes from physically comparable, non-glaciated basins.     

Sr and Nd isotope compositions of atmospheric mineral dust at the summit of Mt. Sefuri, north Kyushu, southwest Japan: A marker of the dust provenance and seasonal variability

We analyzed Rb-Sr-Nd isotope ratios of mineral dust in total aerosol load collected with rainwater continuously from 1998 to 2006 at the summit of Mt. Sefuri, northern Kyushu, southwestern Japan. During this period, the total mass of the dust generally increased in late winter, peaked in early spring, and then decreased.⁸⁷Sr/⁸⁶Sr in atmospheric mineral dust varied from 0.7096 to 0.7180, and εNd(0)_CHUR from −19.9 to −3.5. During heavy deposition periods, the dust had high ⁸⁷Sr/⁸⁶Sr isotope ratios and low to middle εNd(0)_CHUR values, respectively. These compositions are comparable to those of sand and loess in arid areas of Northeast China, Takla Makan and Western Beijing. Such particles were transported by westerlies from those areas to northern Kyushu in winter and spring. In summer and autumn, the isotopic compositions of the dust varied greatly; however, during light deposition periods, the Sr isotope composition was low. In these seasons, the contributions to the dust from Japanese soils and volcanic ash, transported by southern winds, were relatively larger than in winter and spring because of decreased mineral dust particle transport from the continent. Nevertheless, fine sandy desert particles and loess in general accounted for most mineral dust deposition in northern Kyushu year-round, even in summer. Local soils derived from weathered granite and volcanic ash were minor components only.The net mass of water-insoluble inorganic matter in the collected mineral dust fluctuated from year-to-year; deposition on Mt. Sefuri was relatively large before 2001, decreased from 2002 to 2005, and increased greatly in spring 2006. These year-to-year differences probably reflected changes in the strength of the westerlies and in climate conditions in the arid source areas.     

The isotopic evolution of atmospheric Pb in central Ontario since AD 1800, and its impacts on the soils, waters, and sediments of a forested watershed, Kawagama Lake

A peat core from an ombrotrophic bog documents the isotopic evolution of atmospheric Pb in central Ontario since AD 1804 ± 53 (²¹⁰Pb dating). Despite the introduction of unleaded gasoline in the mid-1970’s, the ratio ²⁰⁶Pb/²⁰⁷Pb in atmospheric deposition has not increased as expected, but rather continues to decline. In fact, snowpack sampling (2005 and 2009) and rainwater samples (2008) show that the isotopic composition of atmospheric Pb today is often far less radiogenic than the gasoline lead that had been used in Canada in the past. The peat, snow, and rainwater data presented here are consistent with the Pb isotope data for aerosols collected in Dorset in 1984 and 1986 which were traced by Sturges and Barrie (1989) to emissions from the Noranda smelter in northern Quèbec, Canada’s largest single source of atmospheric Pb. Understanding atmospheric Pb deposition in central Ontario, therefore, requires not only consideration of natural sources and past contributions from leaded gasoline, but also emissions from metal smelting and refining.Lead in the streams which enter Kawagama Lake today (²⁰⁶Pb/²⁰⁷Pb = 1.16 − 1.19) represents a mixture between the natural values (1.191 − 1.201 estimated using pre-industrial lake sediments) and the values found in the humus layer of the surrounding forest soils (²⁰⁶Pb/²⁰⁷Pb = 1.15 − 1.19). In the lake itself, however, Pb is much less radiogenic (²⁰⁶Pb/²⁰⁷Pb as low as 1.09) than in the streams, with the dissolved fraction less radiogenic than particulate material. The evolution of Pb isotope ratios within the watershed apparently reflects preferential removal by sedimentation of comparatively dense, radiogenic, terrestrial particles (derived from the mineral fraction of soils) from the humus particles with lower ratios of ²⁰⁶Pb/²⁰⁷Pb (because of atmospheric Pb contamination). Despite the contemporary enrichments of Pb in rain and snow, concentrations of dissolved Pb in the lake are extremely low (sometimes below 10 ng/l), with Pb concentrations and Pb/Sc ratios approaching “natural” values because of efficient binding to particles, and their subsequent removal in the watershed.     

A threshold in soil formation at Earth’s arid-hyperarid transition

The soils of the Atacama Desert in northern Chile have long been known to contain large quantities of unusual salts, yet the processes that form these soils are not yet fully understood. We examined the morphology and geochemistry of soils on post-Miocene fans and stream terraces along a south-to-north (27° to 24° S) rainfall transect that spans the arid to hyperarid transition (21 to ∼2 mm rain y⁻¹). Landform ages are ? 2 My based on cosmogenic radionuclide concentrations in surface boulders, and Ar isotopes in interbedded volcanic ash deposits near the driest site indicate a maximum age of 2.1 My. A chemical mass balance analysis that explicitly accounts for atmospheric additions was used to quantify net changes in mass and volume as a function of rainfall. In the arid (21 mm rain y⁻¹) soil, total mass loss to weathering of silicate alluvium and dust (−1030 kg m⁻²) is offset by net addition of salts (+170 kg m⁻²). The most hyperarid soil has accumulated 830 kg m⁻² of atmospheric salts (including 260 kg sulfate m⁻² and 90 kg chloride m⁻²), resulting in unusually high volumetric expansion (120%) for a soil of this age. The composition of both airborne particles and atmospheric deposition in passive traps indicates that the geochemistry of the driest soil reflects accumulated atmospheric influxes coupled with limited in-soil chemical transformation and loss. Long-term rates of atmospheric solute addition were derived from the ion inventories in the driest soil, divided by the landform age, and compared to measured contemporary rates. With decreasing rainfall, the soil salt inventories increase, and the retained salts are both more soluble and present at shallower depths. All soils generally exhibit vertical variation in their chemistry, suggesting slow and stochastic downward water movement, and greater climate variability over the past 2 My than is reflected in recent (∼100 y) rainfall averages. The geochemistry of these soils shows that the transition from arid to hyperarid rainfall levels marks a fundamental geochemical threshold: in wetter soils, the rate and character of chemical weathering results in net mass loss and associated volumetric collapse after 10⁵ to 10⁶ years, while continuous accumulation of atmospheric solutes in hyperarid soils over similar timescales results in dramatic volumetric expansion. The specific geochemistry of hyperarid soils is a function of atmospheric sources, and is expected to vary accordingly at other hyperarid sites. This work identifies key processes in hyperarid soil formation that are likely to be independent of location, and suggests that analogous processes may occur on Mars.     

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.     

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.     

The role of trace minerals in chemical weathering in a high-elevation granitic watershed (Estibère, France): chemical and mineralogical evidence

We investigated chemical weathering in a high elevation granitic environment in three selected watersheds located in the Pyrenees (France). The sites were located on glacial deposits derived from similar Hercynian (∼300 Ma) granites characterized by the occurrence of zoned plagioclases and trace calcic phases (epidote, prehnite, sphene, apatite). The surface waters at those sites show high Ca/Na molar ratios (>1) which could not be explained by the dissolution of the major plagioclase (oligoclase) present in the rocks. The coupled approach of investigating stream water chemistry and the mineralogy and chemistry of rocks and soils allowed us to explore the role of the weathering of trace calcic minerals in calcium export at the watershed scale. The weathering of the trace calcic minerals which represent ∼ 1% of the total rock volume are responsible for more than 90% of the calcium export at the sites. Annual cationic fluxes (∼ 23.10⁴ eq/km²/yr) calculated for the Estibère watershed are among the highest reported for high elevation systems draining granitic rocks and ∼ 80% of this annual cationic flux can be attributed to the weathering of trace calcic phases. Calculations based on isotopic values (⁸⁷Sr/⁸⁶Sr) go in the same direction. Except apatite, the trace calcic phases appear to be mainly silicates, thus the type of chemical weathering observed in the Estibère watershed may have an influence on atmospheric CO₂ consumption by granite weathering. However, comparison with other watersheds draining granitic environments worldwide, and with the two other sites in the Pyrenees, indicate that the role of trace calcic phases is important in most young environments exposed to chemical weathering (e.g., high elevation catchments on glacial deposits). Other factors such as the date of glacial retreat, the physical denudation rate, the hydrological functioning of the watershed and the nature and structure of the soil cover are also important.     

Lithium isotope systematics in a forested granitic catchment (Strengbach, Vosges Mountains, France)

Over the last decade it has become apparent that Li isotopes may be a good proxy to trace silicate weathering. However, the exact mechanisms which drive the behaviour of Li isotopes in surface environments are not totally understood and there is a need to better calibrate and characterize this proxy. In this study, we analysed the Li concentrations and isotopic compositions in the various surface reservoirs (soils, rocks, waters and plants) of a small forested granitic catchment located in the Vosges Mountains (Strengbach catchment, France, OHGE http://ohge.u-strasbg.fr). Li fluxes were calculated in both soil profiles and at the basin scale and it was found that even in this forested basin, atmospheric inputs and litter fall represented a minor flux compared to input derived from the weathering of rocks and soil minerals (which together represent a minimum of 70% of dissolved Li). Li isotope ratios in soil pore waters show large depth dependent variations. Average dissolved δ⁷Li decreases from −1.1‰ to −14.4‰ between 0 and −30 cm, but is +30.7‰ at −60 cm. This range of Li isotopic compositions is very large and it encompasses almost the entire range of terrestrial Li isotope compositions that have been previously reported. We interpret these variations to result from both the dissolution and precipitation of secondary phases. Large isotopic variations were also measured in the springs and stream waters, with δ⁷Li varying from +5.3‰ to +19.6‰. δ⁷Li increases from the top to the bottom of the basin and also covaries with discharge at the outlet. These variations are interpreted to reflect isotopic fractionations occurring during secondary phase precipitation along the water pathway through the rocks. We suggest that the dissolved δ⁷Li increases with increasing residence time of waters through the rocks, and so with increasing time of interaction between waters and solids. A dissolution precipitation model was used to fit the dissolved Li isotopic compositions. It was found that the isotopic compositions of springs and stream waters are explicable by an isotopic fractionation of −5‰ to −14‰ (best fit −10.8‰), in agreement with Li incorporation into clay. In soil solutions, it was found that isotopic fractionation during secondary precipitation is larger (at least −23‰), suggesting a major role for different secondary phases, such as iron oxides that maybe incorporate Li with a higher isotope fractionation.     

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.