Landscape variations in stream water SO4 and δSSO4 in a boreal stream network

Despite reduced anthropogenic deposition during the last decades, deposition sulphate may still play an important role in the biogeochemical cycles of S and many catchments may act as net sources of S that may remain for several decades. The aim of this study is to elucidate the temporal and spatial dynamics of both SO₄²⁻ and δ³⁴S_SO4 in stream water from catchments with varying percentage of wetland and forest coverage and to determine their relative importance for catchment losses of S. Stream water samples were collected from 15 subcatchments ranging in size from 3 to 6780 ha, in a boreal stream network, northern Sweden. In forested catchments (<2% wetland cover) S-SO₄²⁻ concentrations in stream water averaged 1.7 mg L⁻¹ whereas in wetland dominated catchments (>30% wetland cover) the concentrations averaged 0.3 mg L⁻¹. A significant negative relationship was observed between S-SO₄²⁻ and percentage wetland coverage (r² = 0.77, p < 0.001) and the annual export of stream water SO₄²⁻ and wetland coverage (r² = 0.76, p < 0.001). The percentage forest coverage was on the other hand positively related to stream water SO₄²⁻ concentrations and the annual export of stream water SO₄²⁻ (r² = 0.77 and r² = 0.79, respectively). The annual average δ³⁴S_SO4 value in wetland dominated streams was +7.6‰ and in streams of forested catchments +6.7‰. At spring flood the δ³⁴S_SO4 values decreased in all streams by 1‰ to 5‰. The δ³⁴S_SO4 values in all streams were higher than the δ³⁴S_SO4 value of +4.7‰ in precipitation (snow). The export of S ranged from 0.5 kg S ha⁻¹ yr⁻¹ (wetland headwater stream) to 3.8 kg S ha⁻¹ yr⁻¹ (forested headwater stream). With an average S deposition in open field of 1.3 kg S ha⁻¹ yr⁻¹ (2002-2006) the mass balance results in a net export of S from all catchments, except in catchments with >30% wetland. The high temporal and spatial resolution of this study demonstrates that the reducing environments of wetlands play a key role for the biogeochemistry of S in boreal landscapes and are net sinks of S. Forested areas, on the other hand were net sources of S.     

Patterns of Ca/Sr and Sr/Sr variation before and after a whole watershed CaSiO3 addition at the Hubbard Brook Experimental Forest, USA

Forty-one metric tons of the mineral wollastonite (CaSiO₃) was applied to an 11.8 hectare watershed at the Hubbard Brook Experimental Forest (HBEF; White Mountains, New Hampshire, USA) with the goal of restoring the Ca estimated to have been depleted from the soil exchange complex by acid deposition. This experiment provided an opportunity to gain qualitative information on whole watershed hydrologic flow paths by studying the response of stream water chemistry to the addition of Ca. Because the Ca/Sr and ⁸⁷Sr/⁸⁶Sr ratios of wollastonite strongly contrast that of other Ca sources in the watershed, the wollastonite-derived Ca can be identified and its amount estimated in various ecosystem components. Stream water chemistry at the HBEF varies seasonally due to shifts in the proportion of base flow and interflow. Prior to the wollastonite application, seasonal variations in ⁸⁷Sr/⁸⁶Sr ratios indicated that ⁸⁷Sr/⁸⁶Sr was higher during base flow than interflow, due largely to greater amounts of biotite weathering along deeper flow paths. After the application, Ca/Sr and ⁸⁷Sr/⁸⁶Sr changed markedly as the high Ca/Sr and low ⁸⁷Sr/⁸⁶Sr wollastonite dissolved and mixed with stream water. The Ca addition provided information on the response times of various flow paths and ion exchange processes to Ca addition in this small upland watershed. During the first year after the addition, wollastonite applied to the near stream zone dissolved and was partially immobilized by cation exchange sites in the hyporheic zone. In the second and third years after the addition we infer that much of this Ca and Sr was subsequently desorbed from the hyporheic zone and was exported from the watershed in stream flow. In the fourth through ninth years after the addition, Ca and Sr from wollastonite that had dissolved in upland soils was transported to the stream by interflow during wet periods when the ground water table was elevated. Between years three and nine the minimum annual Ca/Sr ratio (in late summer base flow) increased, providing evidence that Ca and Sr had increasingly infiltrated to the deepest flow paths. Strong seasonal variations in Ca/Sr and ⁸⁷Sr/⁸⁶Sr ratios of stream water resulted from the wollastonite addition to upland forest soils, and these ratios have become sensitive to changing flow paths during the annual cycle. Most notably, high flow events now produce large excursions in stream geochemistry toward the high Ca/Sr and low ⁸⁷Sr/⁸⁶Sr ratios of wollastonite. Nine years after the application we estimate that ∼360 kg of Ca from wollastonite has been exported from the watershed in stream flow. The rate of export of Ca from wollastonite dissolution has stabilized at about 11 kg of Ca per year, which accounts for ∼30% of the dissolved Ca in the stream water. Given that 19 metric tons of Ca were applied to the watershed, and assuming this current rate of loss, it should take over 1000 years for this added Ca to be transported from the watershed.     

Formation,fate and leaching of chloroform in coniferous forest soils

Chloroform is a common groundwater pollutant but also a natural compound in forest ecosystems. Leaching of natural chloroform from forest soil to groundwater was followed by regular analysis of soil air and groundwater from multilevel wells at four different sites in Denmark for a period of up to 4 a. Significant seasonal variation in chloroform was observed in soil air 0.5 m below surface ranging at one site from 120 ppb by volume in summer to 20 ppb during winter. With depth, the seasonal variation diminished gradually, ranging from 30 ppb in summer to 20 ppb during winter, near the groundwater table. Chloroform in the shallowest groundwater ranged from 0.5–1.5 μg L⁻¹ at one site to 2–5 μg L⁻¹ at another site showing no clear correlation with season. Comparing changes in chloroform in soil air versus depth with on-site recorded meteorological data indicated that a clear relationship appears between rain events and leaching of chloroform. Chloroform in top soil air co-varied with CO₂ given a delay of 3–4 weeks providing evidence for its biological origin. This was confirmed during laboratory incubation experiments which further located the fermentation layer as the most chloroform producing soil horizon. Sorption of chloroform to soils, examined using ¹⁴C–CHCl₃, correlated with organic matter content, being high in the upper organic rich soils and low in the deeper more minerogenic soils. The marked decrease in chloroform in soil with depth may in part be due to microbial degradation which was shown to occur at all depths by laboratory tests using ¹⁴C–CHCl₃.     

Hydrogeochemistry of Fe and Mn in small boreal streams: The role of seasonality, landscape type and scale

Stream water from a stream network of 15 small boreal catchments (0.03-67 km²) in northern Sweden was analyzed for unfiltered (total) and filtered (<0.4 μm) concentrations of iron (Fe_tot and Fe_<0.4) and manganese (Mn_tot and Mn_<0.4). The purpose was to investigate the temporal and spatial dynamics of Fe, Mn and dissolved organic carbon (DOC) as influenced by snow melt driven spring floods and landscape properties, in particular the proportion of wetland area. During spring flood, concentrations of Fe_tot, Fe_<0.4, Mn_tot, Mn_<0.4 and DOC increased in streams with forested catchments (<2% wetland area). In catchments with high coverage of wetlands (>30% wetland area) the opposite behavior was observed. The hydrogeochemistry of Fe was highly dependent on wetlands as shown by the strong positive correlation of the Fe_tot/Al_tot ratio with wetland coverage (r² = 0.89, p < 0.001). Furthermore, PCA analysis showed that at base flow Fe_tot and Fe_<0.4 were positively associated with wetlands and DOC, whereas they were not associated during peak flow at spring flood. The temporal variation of Fe was likely related to varying hydrological pathways. At peak discharge Fe_tot was associated with variables like silt coverage, which highlights the importance of particulates during high discharge events. For Mn there was no significant correlation with wetlands, instead, PCA analysis showed that during spring flood Mn was apparently more dependent on the supply of minerogenic particulates from silt deposits on the stream banks of some of the streams. The influence of minerogenic particulates on the concentration of, in particular, Mn was greatest in the larger, lower gradient streams, characterized by silt deposits in the near-stream zone. In the small forested streams underlain by till, DOC was of greater importance for the observed concentrations, as indicated by the positive correlation of both Fe_tot and Fe_<0.4 with DOC (r² = 0.77 and r² = 0.76, p < 0.001) at the smallest headwater forest site. In conclusion, wetland area and DOC were important for Fe concentrations in this boreal stream network, whereas silt deposits strongly influenced Mn concentrations. This study highlights the importance of studying stream water chemistry from a landscape perspective in order to address future environmental issues concerning mobility of Fe, Mn and associated trace metals.     

Predictive modelling of the spatial pattern of past and future forest cover changes in India

This study was carried out to simulate the forest cover changes in India using Land Change Modeler. Classified multi-temporal long-term forest cover data was used to generate the forest covers of 1880 and 2025. The spatial data were overlaid with variables such as the proximity to roads, settlements, water bodies, elevation and slope to determine the relationship between forest cover change and explanatory variables. The predicted forest cover in 1880 indicates an area of 10,42,008 km², which represents 31.7% of the geographical area of India. About 40% of the forest cover in India was lost during the time interval of 1880–2013. Ownership of majority of forest lands by non-governmental agencies and large scale shifting cultivation are responsible for higher deforestation rates in the Northeastern states. The six states of the Northeast (Assam, Manipur, Meghalaya, Mizoram, Nagaland, Tripura) and one union territory (Andaman & Nicobar Islands) had shown an annual gross rate of deforestation of >0.3 from 2005 to 2013 and has been considered in the present study for the prediction of future forest cover in 2025. The modelling results predicted widespread deforestation in Northeast India and in Andaman & Nicobar Islands and hence is likely to affect the remaining forests significantly before 2025. The multi-layer perceptron neural network has predicted the forest cover for the period of 1880 and 2025 with a Kappa statistic of >0.70. The model predicted a further decrease of 2305 km² of forest area in the Northeast and Andaman & Nicobar Islands by 2025. The majority of the protected areas are successful in the protection of the forest cover in the Northeast due to management practices, with the exception of Manas, Sonai-Rupai, Nameri and Marat Longri. The predicted forest cover scenario for the year 2025 would provide useful inputs for effective resource management and help in biodiversity conservation and for mitigating climate change.     

Naturally acidic surface and ground waters draining porphyry-related mineralized areas of the Southern Rocky Mountains,Colorado and New Mexico

Acidic, metal-rich waters produced by the oxidative weathering and resulting leaching of major and trace elements from pyritic rocks can adversely affect water quality in receiving streams and riparian ecosystems. Five study areas in the southern Rocky Mountains with naturally acidic waters associated with porphyry mineralization were studied to document variations in water chemistry and processes that control the chemical variations. Study areas include the Upper Animas River watershed, East Alpine Gulch, Mount Emmons, and Handcart Gulch in Colorado and the Red River in New Mexico. Although host-rock lithologies in all these areas range from Precambrian gneisses to Cretaceous sedimentary units to Tertiary volcanic complexes, the mineralization is Tertiary in age and associated with intermediate to felsic composition, porphyritic plutons. Pyrite is ubiquitous, ranging from ∼1 to >5 vol.%. Springs and headwater streams have pH values as low as 2.6, SO₄ up to 3700 mg/L and high dissolved metal concentrations (for example: Fe up to 400 mg/L; Cu up to 3.5 mg/L; and Zn up to 14.4 mg/L). Intensity of hydrothermal alteration and presence of sulfides are the primary controls of water chemistry of these naturally acidic waters. Subbasins underlain by intensely hydrothermally altered lithologies are poorly vegetated and quite susceptible to storm-induced surface runoff. Within the Red River study area, results from a storm runoff study documented downstream changes in river chemistry: pH decreased from 7.80 to 4.83, alkalinity decreased from 49.4 to <1 mg/L, SO₄ increased from 162 to 314 mg/L, dissolved Fe increased from to 0.011 to 0.596 mg/L, and dissolved Zn increased from 0.056 to 0.607 mg/L. Compared to mine drainage in the same study areas, the chemistry of naturally acidic waters tends to overlap but not reach the extreme concentrations of metals and acidity as some mine waters. The chemistry of waters draining these mineralized but unmined areas can be used to estimate premining conditions at sites with similar geologic and hydrologic conditions. For example, the US Geological Survey was asked to estimate premining ground-water chemistry at the Questa Mo mine, and the proximal analog approach was used because a mineralized but unmined area was located adjacent to the mine property. By comparing and contrasting water chemistry from different porphyry mineralized areas, this study not only documents the range in concentrations of constituents of interest but also provides insight into the primary controls of water chemistry.     

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.     

Changes in major element hydrochemistry of the Pecos River in the American Southwest since 1935

The Pecos River, situated in eastern New Mexico and western Texas, receives water from a drainage area of 91 000 km². There are primarily two major water inputs, namely snowmelt from winter storms in the headwater region of the southern Rocky Mountains and runoff from warm-season monsoonal rainfall in the lower valley. The Pecos River suffers from high levels of total dissolved solids (TDS >5000 mg L⁻¹) under normal flow conditions. This not only poses serious problems for agricultural irrigation and safe drinking water supply, but also results in a permanent loss of biodiversity. This study examines changes in stream flow and water chemistry of the Pecos River over the last 70 a to better understand the long-term variability in stream salinity and the role of agricultural practices in salt transfer. A TDS record from the lower Pecos River near Langtry (Texas) back to 1935 was extracted to show a distinct pattern of decadal variability similar to the Pacific Decadal Oscillation (PDO), in which stream salinity is overall above average when the PDO is in positive (warm) phase and below average when the PDO is in negative (cold) phase. This is due to: (1) the dissolved salts contributed to the river are largely from dissolution of NaCl and CaSO₄-bearing minerals (e.g., halite and gypsum) in the upper basin, (2) the amount of the dissolved salts that reach the lower basin is mainly determined by the stream flow yield in the upper basin and (3) the stream flow yield from the upper basin is positively correlated with the PDO index. This further attests that large-scale climatic oscillation is the major source of long-term changes in stream flow and salinity of the Pecos River. On the other hand, there is also a strong indication that the rate of salt export has been affected by reservoir operations and water diversions for agricultural practices.     

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.     

Natural acid rock drainage associated with black shale in the Yukon Territory,Canada

Investigated herein are water and sediment geochemistry, and metal attenuation processes associated with natural acid rock drainage originating from black shale formations in the Macmillan Pass area, Clear Lake prospect and Engineer Creek by the Dempster Highway in the Yukon Territory, Canada. The most metalliferous water having pH 3.0, 150 mg/L Zn, 39 mg/L Ni, 2.8 mg/L Cu and 9.1 mg/L As was found in a tributary stream of Engineer Creek with no known mineral deposits occurring in the vicinity. For all three study areas, the water and sediment geochemistry is significantly affected by the local lithology and prevailing metal attenuation processes. Despite their anomalous acidity and metal contents, the natural acid streams contribute only a small fraction of the contaminant loadings to the major water courses because of their low flows. Dilution, neutralization, sorption and co-precipitation are identified as the major mechanisms attenuating aqueous transport of potentially deleterious metals. However, microbial mediation in metal attenuation is also evident in low-flow systems. The wide variation of water and sediment geochemistry along a flow path renders the establishment of background metal values difficult. In assessing environmental impacts, it may be more practical to consider metal loadings on a watershed scale than to rely on a comparison with operationally defined background concentrations.     

Estimating specific surface area of fine stream bed sediments from geochemistry

Specific surface area (SSA) of headwater stream bed sediments is a fundamental property which determines the nature of sediment surface reactions and influences ecosystem-level, biological processes. Measurements of SSA – commonly undertaken by BET nitrogen adsorption – are relatively costly in terms of instrumentation and operator time. A novel approach is presented for estimating fine (<150 μm) stream bed sediment SSA from their geochemistry – after removal of organic matter – for agricultural headwater catchments across 15,400 km² of central England, UK. From a regional set of 1972 stream bed sediment sites with common characteristics for which geochemical data were available, 60 samples were selected – based on maximising their variation in Al concentrations – and their BET SSA measured by N₂ adsorption. After careful selection of potential regression predictors following a principal component analysis and removal of a subset of samples with the largest Mo concentrations (>2.5 mg kg⁻¹), four elements were identified as significant predictors of SSA (ordered by decreasing predictive power): V > Ca > Al > Rb. The optimum model from these four elements accounted for 73% of the variation in bed sediment SSA (range 6–46 m² g⁻¹) with a root mean squared error of prediction – based on leave-one-out cross-validation – of 6.3 m² g⁻¹. It is believed that V is the most significant predictor because its concentration is strongly correlated both with the quantity of Fe-oxides and clay minerals in the stream bed sediments, which dominate sediment SSA. Sample heterogeneity in SSA – based on triplicate measurements of sub-samples – was a substantial source of variation (standard error = 2.2 m² g⁻¹) which cannot be accounted for in the regression model.     

Effect of precipitation regime on δD values of soil n-alkanes from elevation gradients - Implications for the study of paleo-elevation

We measured δD values of long chain n-alkanes isolated from 30 surface soil samples along two elevation transects on the Tibetan Plateau differing in precipitation regime and water source. The East Asian Monsoon precipitation dominates the wetter regime on the eastern slope (from 1230 to 4300 m) of Gongga Shan on the eastern Tibetan Plateau. Precipitation from the Polar Westerlies dominates the drier region on the slope from 1900 to 5000 m in the West Kunlun Shan on the northwestern Tibetan Plateau. The decrease in δD value with elevation in the wetter region greatly exceeded that in drier region by, −1.9 ± 0.1‰/100 m and −1.4 ± 1.0‰/100 m respectively. The apparent fractionation between leaf wax and precipitation ε_wax-p values in the wetter region (ca. −164‰) were more negative than those in drier region (ca. −125‰ above 3200 m).We also measured δD values in leaves of six common living trees (values from −287‰ to −193‰) from Gongga Shan, ranging from about 2900-4200 m. The abundance-weighted average values of the n-alkanes (δD_wax) show a strong reverse correlation with sample source elevation (R² 0.78 for soils from Gongga Shan; R² 0.85 for soils from West Kunlun Shan above 3200 m), suggesting that n-alkane δD_wax faithfully records the precipitation δD and that the isotopic altitude effect of precipitation controls δD_wax altitudinal gradients in the mountains. The data show a fairly strong monotonic dependency of n-alkane δD values on elevation for the eastern Plateau, but a complex relationship between n-alkane δD values and elevation for the northwestern Plateau. The δD_wax values at sites below 3200 m from the Kunlun Shan area exhibit an unexpected positive correlation with elevation. The study confirms the potential for using sediment δD_wax values to reconstruct paleo-elevation in wetter regions, but suggests caution in applying the approach to dry regions. Our results also show it is essential to consider the intricacy of the pattern of atmospheric circulation and water sources and their influence on the lapse rate of δD values with elevation.     

Water Geochemistry of Three Mountain Streams from Carbonate Watersheds in the Southern French Alps

The chemical composition of water from three streamsflowing through a carbonate watershed wasinvestigated. Although the study area is not spatiallyvery large (a few km²), local inhomogeneitieswithin the lithology appreciably affect the chemicalcomposition of the water and the geochemical gradientas a function of elevation. The main chemical processwhich leads to the observed water chemistry is thedissolution of calcium carbonate by atmospheric andmetabolic CO₂. In the stream La Sigouste, thewater dissolves nearly pure calcite and this reactionproceeds until an equilibrium with respect to calciteis reached. In Le Lauzon stream, local inhomogeneitiesin the lithologic composition prevent theestablishment of an unambiguous weathering budgetwhile, for Le Rif de l'Arc a stoichiometric modelshows that the weathering process is adequatelyrepresented by the dissolution of calcite associatedwith a minor incongruent dissolution of chloritepresent in marls. For the two last streamssupersaturation with respect to calcite is observed.For both streams, supersaturation with respect toatmospheric CO₂ is nearly permanent. The observedaltitudinal gradients of sulfate are interpretedthrough a set of biogeochemical redox reactions.     

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₂.     

Rare earth elements (REE) and yttrium in stream waters, stream sediments, and Fe-Mn oxyhydroxides: Fractionation, speciation, and controls over REE + Y patterns in the surface environment

We have collected ∼500 stream waters and associated bed-load sediments over an ∼400 km² region of Eastern Canada and analyzed these samples for Fe, Mn, and the rare earth elements (REE + Y). In addition to analyzing the stream sediments by total digestion (multi-acid dissolution with metaborate fusion), we also leached the sediments with 0.25 M hydroxylamine hydrochloride (in 0.05 M HCl), to determine the REE + Y associated with amorphous Fe- and Mn-oxyhydroxide phases. We are thus able to partition the REE into “dissolved” (<0.45 μm), labile (hydroxylamine) and detrital sediment fractions to investigate REE fractionation, and in particular, with respect to the development of Ce and Eu anomalies in oxygenated surface environments. Surface waters are typically LREE depleted ([La/Sm]_NASC ranges from 0.16 to 5.84, average = 0.604, n = 410; where the REE are normalized to the North America Shale Composite), have strongly negative Ce anomalies ([Ce/Ce^∗]_NASC ranges from 0.02 to 1.25, average = 0.277, n = 354), and commonly have positive Eu anomalies ([Eu/Eu^∗]_NASC ranges from 0.295 to 1.77, average = 0.764, n = 84). In contrast, the total sediment have flatter REE + Y patterns relative to NASC ([La/Sm]_NASC ranges from 0.352 to 1.12, average = 0.778, n = 451) and are slightly middle REE enriched ([Gd/Yb]_NASC ranges from 0.55 to 3.75, average = 1.42). Most total sediments have negative Ce and Eu anomalies ([Ce/Ce^∗]_NASC ranges from 0.097 to 2.12, average = 0.799 and [Eu/Eu^∗]_NASC ranges from 0.39 to 1.43, average = 0.802). The partial extraction sediments are commonly less LREE depleted than the total sediments ([La/Sm]_NASC ranges from 0.24 to 3.31, average = 0.901, n = 4537), more MREE enriched ([Gd/Yb]_NASC ranges from 0.765 to 6.28, average = 1.97) and Ce and Eu anomalies (negative and positive) are more pronounced.The partial extraction recovered, on average ∼20% of the Fe in the total sediment, ∼80% of the Mn, and 21-29% of the REEs (Ce = 19% and Y = 32%). Comparison between REEs in water, partial extraction and total sediment analyses indicates that REEs + Y in the stream sediments have two primary sources, the host lithologies (i.e., mechanical dispersion) and hydromorphically transported (the labile fraction). Furthermore, Eu appears to be more mobile than the other REE, whereas Ce is preferentially removed from solution and accumulates in the stream sediments in a less labile form than the other REEs + Y. Despite poor statistical correlations between the REEs + Y and Mn in either the total sediment or partial extractions, based on apparent distribution coefficients and the pH of the stream waters, we suggest that either sediment organic matter and/or possibly δ-MnO₂/FeOOH are likely the predominant sinks for Ce, and to a lesser extent the other REE, in the stream sediments.     

Land use land cover dynamics as a function of changing demography and hydrology

This paper describes the spatiotemporal changes pertaining to land use land cover (LULC) and the driving forces behind these changes in Doodhganga watershed of Jhelum Basin. An integrated approach utilizing remote sensing and geographic information system (GIS) was used to extract information pertaining to LULC change. Multi-date LULC maps were generated by analyzing remotely sensed images of three dates which include LandSat TM 1992, LandSat ETM+ 2001 and IRS LISS-III 2005. The LULC information was extracted by adopting on-screen image interpretation technique in a GIS environment at 1:25,000 scale. Based on the analysis, changes were observed in the spatial extent of different LULC types over a period of 13 years. Significant changes were observed in the spatial extent of forest, horticulture, built-up and agriculture. Forest cover in the watershed has decreased by 1.47 %, Agricultural by 0.93 % while as built-up area has increased by 0.92 %. The net decrease in forest cover and agriculture land indicate the anthropogenic interference into surrounding natural ecosystems. From the study it was found that the major driving forces for these changes were population growth and changes in the stream discharge. The changes in the stream discharge were found responsible for the conversion of agricultural land into horticulture, as horticulture has increased by 1.14 % in spatial extent. It has been found that increasing human population together with decreasing stream discharge account for LULC changes in the watershed. Therefore, the existing policy framework needs to focus upon mitigating the impacts of forces responsible for LULC change so as to ensure sustainable development of land resources.     

Sr isotopes as a tracer of weathering processes and dust inputs in a tropical granitoid watershed, Luquillo Mountains, Puerto Rico

Sr isotope data from soils, water, and atmospheric inputs in a small tropical granitoid watershed in the Luquillo Mountains of Puerto Rico constrain soil mineral development, weathering fluxes, and atmospheric deposition. This study provides new information on pedogenic processes and geochemical fluxes that is not apparent in watershed mass balances based on major elements alone. ⁸⁷Sr/⁸⁶Sr data reveal that Saharan mineral aerosol dust contributes significantly to atmospheric inputs. Watershed-scale Sr isotope mass balance calculations indicate that the dust deposition flux for the watershed is 2100 ± 700 mg cm⁻² ka⁻¹. Nd isotope analyses of soil and saprolite samples provide independent evidence for the presence of Saharan dust in the regolith. Watershed-scale Sr isotope mass balance calculations are used to calculate the overall short-term chemical denudation velocity for the watershed, which agrees well with previous denudation rate estimates based on major element chemistry and cosmogenic nuclides. The dissolved streamwater Sr flux is dominated by weathering of plagioclase and hornblende and partial weathering of biotite in the saprock zone. A steep gradient in regolith porewater ⁸⁷Sr/⁸⁶Sr ratio with depth, from 0.70635 to as high as 0.71395, reflects the transition from primary mineral-derived Sr to a combination of residual biotite-derived Sr and atmospherically-derived Sr near the surface, and allows multiple origins of kaolinite to be identified.     

Does reactive surface area depend on grain size? Results from pH 3, 25 °C far-from-equilibrium flow-through dissolution experiments on anorthite and biotite

Laboratory determined mineral weathering rates need to be normalised to allow their extrapolation to natural systems. The principle normalisation terms used in the literature are mass, and geometric- and BET specific surface area (SSA). The purpose of this study was to determine how dissolution rates normalised to these terms vary with grain size. Different size fractions of anorthite and biotite ranging from 180-150 to 20-10 μm were dissolved in pH 3, HCl at 25 °C in flow through reactors under far from equilibrium conditions. Steady state dissolution rates after 5376 h (anorthite) and 4992 h (biotite) were calculated from Si concentrations and were normalised to initial- and final- mass and geometric-, geometric edge- (biotite), and BET SSA. For anorthite, rates normalised to initial- and final-BET SSA ranged from 0.33 to 2.77 × 10⁻¹⁰ mol_feldspar m⁻² s⁻¹, rates normalised to initial- and final-geometric SSA ranged from 5.74 to 8.88 × 10⁻¹⁰ mol_feldspar m⁻² s⁻¹ and rates normalised to initial- and final-mass ranged from 0.11 to 1.65 mol_feldspar g⁻¹ s⁻¹. For biotite, rates normalised to initial- and final-BET SSA ranged from 1.02 to 2.03 × 10⁻¹² mol_biotite m⁻² s⁻¹, rates normalised to initial- and final-geometric SSA ranged from 3.26 to 16.21 × 10⁻¹² mol_biotite m⁻² s⁻¹, rates normalised to initial- and final-geometric edge SSA ranged from 59.46 to 111.32 × 10⁻¹² mol_biotite m⁻² s⁻¹ and rates normalised to initial- and final-mass ranged from 0.81 to 6.93 × 10⁻¹² mol_biotite g⁻¹ s⁻¹. For all normalising terms rates varied significantly (p ? 0.05) with grain size. The normalising terms which gave least variation in dissolution rate between grain sizes for anorthite were initial BET SSA and initial- and final-geometric SSA. This is consistent with: (1) dissolution being dominated by the slower dissolving but area dominant non-etched surfaces of the grains and, (2) the walls of etch pits and other dissolution features being relatively unreactive. These steady state normalised dissolution rates are likely to be constant with time. Normalisation to final BET SSA did not give constant ratios across grain size due to a non-uniform distribution of dissolution features. After dissolution coarser grains had a greater density of dissolution features with BET-measurable but unreactive wall surface area than the finer grains. The normalising term which gave the least variation in dissolution rates between grain sizes for biotite was initial BET SSA. Initial- and final-geometric edge SSA and final BET SSA gave the next least varied rates. The basal surfaces dissolved sufficiently rapidly to influence bulk dissolution rate and prevent geometric edge SSA normalised dissolution rates showing the least variation. Simple modelling indicated that biotite grain edges dissolved 71-132 times faster than basal surfaces. In this experiment, initial BET SSA best integrated the different areas and reactivities of the edge and basal surfaces of biotite. Steady state dissolution rates are likely to vary with time as dissolution alters the ratio of edge to basal surface area. Therefore they would be more properly termed pseudo-steady state rates, only appearing constant because the time period over which they were measured (1512 h) was less than the time period over which they would change significantly.     

Effects of watershed and estuarine characteristics on the abundance of submerged aquatic vegetation in Chesapeake Bay subestuaries

Watershed land use can affect submerged aquatic vegetation (SAV) by elevating nutrient and sediment loading to estuaries. We analyzed the effects of watershed use and estuarine characteristics on the spatial variation of SAV abundance among 101 shallow subestuaries of Chesapeake Bay during 1984–2003. Areas of these subestuaries range from 0.1 to 101 km², and their associated local watershed areas range from 6 to 1664 km². Watershed land cover ranges from 6% to 81% forest, 1% to 64% cropland, 2% to 38% grassland, and 0.3% to 89% developed land. Landscape analyses were applied to develop a number of subestuary metrics (such as subestuary area, mouth width, elongation ratio, fractal dimension of shoreline, and the ratio of local watershed area to subestuary area) and watershed metrics (such as watershed area). Using mapped data from aerial SAV surveys, we calculated SAV coverage for each subestuary in each year during 1984–2003 as a proportion of potential SAV habitat (the area < 2 m deep). The variation in SAV abundance among subestuaries was strongly linked with subestuary and watershed characteristics. A regression tree model indicated that 60% of the variance in SAV abundance could be explained by subestuary fractal dimension, mean tidal range, local watershed dominant land cover, watershed to subestuary area ratio, and mean wave height. Similar explanatory powers were found in wet and dry years, but different independent variables were used. Repeated measures ANOVA with multiple-mean comparison showed that SAV abundance declined with the dominant watershed land cover in the order: forested, mixed-undisturbed, or mixed-developed > mixed-agricultural > agricultural > developed. Change-point analyses indicated strong threshold responses of SAV abundance to point source total nitrogen and phosphorus inputs, the ratio of local watershed area to subestuary area, and septic system density in the local watershed.     

A survey of the inorganic chemistry of bottled mineral waters from the British Isles

The inorganic chemistry of 85 samples of bottled natural mineral waters and spring waters has been investigated from 67 sources across the British Isles (England, Wales, Scotland, Northern Ireland, Republic of Ireland). Sources include boreholes, springs and wells. Waters are from a diverse range of aquifer lithologies and are disproportionately derived from comparatively minor aquifers, the most represented being Lower Palaeozoic (10 sources), Devonian Sandstone (10 sources) and Carboniferous Limestone (9 sources). The waters show correspondingly variable major-ion compositions, ranging from Ca–HCO₃, through mixed-cation–mixed-anion to Na–HCO₃ types. Concentrations of total dissolved solids are mostly low to very low (range 58–800 mg/L). All samples analysed in the study had concentrations of inorganic constituents well within the limits for compliance with European and national standards for bottled waters. Concentrations of NO₃–N reached up to half the limit of 11.3 mg/L, although 62% of samples had concentrations <1 mg/L. Concentrations of Ba were high (up to 1010 μg/L) in two spring water samples. Such concentrations would have been non-compliant had they been classed as natural mineral waters, although no limit exists for Ba in European bottled spring water. In addition, though no European limit exists for U in bottled water, should a limit commensurate with the current WHO provisional guideline value for U in drinking water (15 μg/L) be introduced in the future, a small number of groundwater sources would have concentrations close to this value. Two sources had groundwater U concentrations > 10 μg/L, both being from the Welsh Devonian Sandstone. The highest observed U concentration was 13.6 μg/L.