Isotopic temporal and spatial variations of tropical rivers in Thailand reflect monsoon precipitation signals

Stable isotopic compositions (δ¹⁸O and d-excess) from 25 rivers in Thailand were analysed monthly during 2013–2015. Results indicated that monsoon precipitation fundamentally influences the river isotopes. The overland flow supplied from monsoon precipitation and human-altered flow regimes produces considerable isotopic variability. Spatial and temporal variations were observed among four principal geographical regions. The seasonality of monsoon precipitation in mountainous Thailand produced large variations in isotopic compositions because most rainfall occurred during the southwest monsoon, and dry conditions prevailed during the northeast monsoon. The northern and northeastern regions are mountainous, highland areas. Low δ¹⁸O values were found in these regions, likely because of altitude effects on precipitation. Conversely, monsoonal precipitation continually supplies rivers in southern Thailand all year round, producing higher and more consistent δ¹⁸O values than in the other regions. The Chao Phraya plain in the central region experienced enrichment of δ¹⁸O river runoff related to evaporation in irrigation systems. Larger catchment areas and longer residence times resulted in more pronounced evaporation effects, producing lower values of d-excess and local river water line slopes compared with precipitation. The isotopic differences between river waters and precipitation were utilized to determine river recharge elevations and water transit time. The methods presented here can be used to explore hydrological interactions in other tropical river basins.     

Controls on stable isotopic characteristics of water vapour over Thailand

This study examined the weekly water vapour isotopic composition (δ¹⁸O_v) in Thailand. The water vapour was cryogenically collected from eight sites across the country. Two observational samples were collected over one 24-h period each week (a daytime and a night-time sample), from September 2013 to September 2014. The primary aim was to investigate the environmental factors influencing water vapour isotopes. The results revealed differences in water vapour isotopic values between day and night samples. Three periods of depleted δ¹⁸O_v were associated with large-scale convective systems in September, December, and May. The statistical relationship between the climate variables and water vapour isotopes indicated that the amount of precipitation and relative humidity were the primary controls on both diurnal and seasonal isotopic variability. The temperature did not affect the δ¹⁸O_v, mainly because the atmospheric processes are a function of vertical convection rather than temperature in tropical regions. The water vapour deuterium excess (d-excess) showed greater variability in 2013 than in 2014. The d-excess variation reflected the differences in convection occurring in the day and night. In addition, the vapour phase data were combined with the local meteoric water line to identify the local water vapour line and the interaction between the isotopic composition of water vapour and liquid water. The water vapour isotopic patterns paralleled the precipitation isotopes on rainy days because of equilibrium isotopic exchange. Water vapour and precipitation were isotopically similar under low humidity but showed greater differences from each other under wetter conditions. The study results provide insight into water vapour isotopic characteristics in tropical regions and constrain the role of large-scale atmospheric processes relative to isotopic variability of water vapour in Thailand and nearby countries.     

Vapour source and spatiotemporal variation of precipitation isotopes in Southwest Spain

The δ²H and δ¹⁸O composition of 77 precipitation samples collected between January 2014 and April 2019 from two sites across the Guadalquivir Basin, SW Spain, were analysed. The first site is located in an urban area of Seville at 100 km distance to the Atlantic coast and the second site is located in a dune area of the Doñana National Park a few kilometres from the coast. Sampling was performed within intervals of at least 14 days if rain occurred but frequently intervals were longer according to the rainfall incidence. Samples from both sites are available for the period February 2016 to June 2018 with six samples containing identical rain events at both locations. Precipitation weighted averages and local meteoric water lines produced by weighting and non-weighting regression methods are presented for its use in hydrological applications. Results show a remarkably high variability in δ²H and δ¹⁸O values and precipitation weighted average d-excess values of 11.8‰ and 13‰ at the sites Plaza de España and Doñana, respectively. Temperature and amount effects were found to be weak. A significant influence of secondary evaporation for single rainfall events during summer was identified by enriched isotopic signatures with reduced d-excess values plotting close or below the global meteoric water line. Backward trajectory analysis of 115 days with daily rainfall above 3 mm yield a predominant Atlantic Ocean vapour source with negligible Mediterranean influence and therefore, d-excess variability is attributed to the different ocean surface conditions of relative humidity and sea surface temperature. Parallel sampling indicate very similar isotopic signatures at both sites and point to the existence of thermal effects of the Plaza de España site in Seville city during the summer season.     

Spatial and temporal distributions of stable isotopes in precipitation over Thailand

Spatial and temporal variations of the isotopic composition of precipitation over Thailand were investigated. The local meteoric water line for Thailand deviates slightly from the global meteoric water line, with lower slopes (7.62 ± 0.07, 7.59 ± 0.08) and intercepts (6.42 ± 0.39, 6.22 ± 0.42) using ordinary and precipitation weighted methods. Differences in spatial and temporal δ¹⁸O distributions between the tropical monsoon and tropical savanna climate zones were found due to differing moisture source contributions and seasonal precipitation patterns. The temporal data reveals that the northeast monsoon rains originate from isotopically-enriched local moisture with isotope values of −9.36 to −0.09‰ (mean − 3.73 ± 0.42‰), whereas the southwest monsoon clouds had a more significant rainout effect from Rayleigh distillation, with isotope values of −9.56 to −1.78‰ (mean − 5.40 ± 0.38‰). The precipitation amount at each site was negatively correlated with δ¹⁸O (−0.24 to −3.20‰ per 100 mm, R² = 0.1–0.9). Furthermore, δ¹⁸O was negatively correlated with geography (latitude, altitude) for the southwest monsoon periods, as expected based on other observed correlations. However, an inverse correlation was seen in the northeast monsoon due to differing moisture transportation as part of the continental effect. The correlation coefficient (R) was higher in the southwest monsoon (−0.84 for latitude effect, −0.64 for altitude effect) than the northeast monsoon (0.67 for latitude effect, 0.35 for altitude effect). The spatial pattern of isotopic composition reflects the southwest monsoon more clearly than the northeast monsoon, but the two monsoons also have a cancelling impact on orographic patterns. An agreement of the δ¹⁸O and deuterium excess (d-excess) was a negative correlation and found to reflect precipitation sources and re-evaporation processes. The d-excess was slightly higher for the northeast monsoon, bringing moisture from the Pacific Ocean and travelling across the continent before reaching the observed stations. By contrast, the d-excess was relatively lower for the Indian Ocean's moisture in the southwest monsoon.     

Mining lakes as groundwater‐dominated hydrological systems: assessment of the water balance of Mining Lake Plessa 117 (Lusatia,Germany) using stable isotopes

In the present study, the stable isotopes δ¹⁸O and δ²H were used for assessment of the water balance in a heterogeneously structured catchment area in the Lusatian Lignite Mining District, in particular, for estimation of the annual groundwater inflow and outflow (I_GW and O_GW) of Mining Lake Plessa 117. The application of stable isotopes was possible since the water exchange in the catchment area had reached steady‐state conditions after the abandonment of mining activities in 1968 and the filling of the voids and aquifers by re‐rising groundwater in the years thereafter. Diverging slopes of the Evaporation Line and the Global Meteoric Water Line manifested as evaporation from the lake catchment area. The calculated isotope water balance was compared with the commonly used surface water balance, which is unable to differentiate between I_GW and O_GW, and with a local groundwater model. The groundwater model calculated an I_GW of about 811 000 m³ yr^?1 and an O_GW close to zero, whereas the isotope water balance showed fluxes of about 914 000 and 140 000 m³ yr^?1, respectively. Considering the contribution of the groundwater inflow to the total annual input into the lake (ΔI_T) and the mean residence time (τ), where the groundwater model and the isotope water balance calculated 42 and 47% for ΔI_T and 4·3 and 3·9 years for τ, respectively, it was shown that both water balance calculation methods led to comparable results despite the differences in I_GW and O_GW. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydrologic characteristics of lake‐ and stream‐side riparian wetted margins in the McMurdo Dry Valleys,Antarctica

Water is a limiting factor for life in the McMurdo Dry Valleys (MDV), Antarctica. The active layer (seasonally thawed soil overlying permafrost) accommodates dynamic hydrological and biological processes for 10–16 weeks per year. Wetted margins (visually wetted areas with high moisture content) adjacent to lakes and streams are potential locations of great importance in the MDV because of the regular presence of liquid water, compared with the rest of the landscape where liquid water is rare. At 11 plots (four adjacent to lakes, seven adjacent to streams), soil particle size distribution, soil electrical conductivity, soil water content and isotopic signature, width of the wetted margin, and active layer thaw depth were characterised to determine how these gradients influence physicochemical properties that determine microbial habitat and biogeochemical cycling. Sediments were generally coarse‐grained in wetted margins adjacent to both lakes and streams. Wetted margins ranged from 1·04 to 11·01 m in average length and were found to be longer at lakeside sites than streamside. Average thaw depths ranged from 0·12 to 0·85 m, and were found to be deepest under lake margins. Lake margins also had much higher soil electrical conductivity, steeper topographic gradients, but more gradual soil moisture gradients than stream margins. Patterns of soil water δ¹⁸O and δD distribution indicate capillary action and evaporation from wetted margins; margin pore waters generally demonstrated isotopic enrichment with distance from the shore, indicating evaporation of soil water. Lake margin pore waters were significantly more negative in D_XS (D_XS = δD‐8δ¹⁸O) than streamside pore waters, indicating a longer history of evaporation there. Differences between lake and stream margins can be explained by the more consistent availability of water to lake margins than stream margins. Differences in margin characteristics between lakes and streams have important consequences for the microbial habitat of these margins and their functional role in biogeochemical cycling at these terrestrial–aquatic interfaces. Copyright © 2009 John Wiley & Sons, Ltd.     

Hydrological regulation of chemical weathering and dissolved inorganic carbon biogeochemical processes in a monsoonal river

To better understand the mechanisms relating to hydrological regulations of chemical weathering processes and dissolved inorganic carbon (DIC) behaviours, high-frequency sampling campaigns and associated analyses were conducted in the Yu River, South China. Hydrological variability modifies the biogeochemical processes of dissolved solutes, so major ions display different behaviours in response to discharge change. Most ions become diluted with increasing discharge because of the shortened reactive time between rock and water under high-flow conditions. Carbonate weathering is the main source of major ions, which shows strong chemostatic behaviour in response to changes in discharge. Ions from silicate weathering exhibit a significant dilution effect relative to the carbonate-sourced ions. Under high temperatures, the increased soil CO₂ influx from the mineralisation of organic material shifts the negative carbon isotope ratios of DIC (δ¹³C_DIC) during the high-flow season. The δ¹³C_DIC values show a higher sensitivity than DIC contents in response to various hydrological conditions. Results from a modified isotope-mixing model (IsoSource) demonstrate that biological carbon is a dominant source of DIC and plays an important role in temporal carbon dynamics. Furthermore, this study provides insights into chemical weathering processes and carbon dynamics, highlighting the significant influence of hydrological variability to aid understanding of the global carbon cycle.     

Changes in carbon and nitrogen cycling in a floodplain lake over recent decades linked to littoral expansion,declining riverine influx,and eutrophication

This study aimed to understand changes in the biogeochemical processing of organic matter (OM) in response to multiple stressors (e.g., littoral area expansion, wastewater input, and hydrological regulation) in East Dongting Lake (Central China) over the past 60 years, using analyses of total organic carbon (TOC), total nitrogen (TN), C/N ratios, δ¹³C, δ¹⁵N, and diatoms from 2 sediment cores collected from the littoral and central parts of the lake. OM mainly originated from phytoplankton and C₃ plant‐derived soil OM based on the ranges of C/N ratios (from 7 to 11) and δ¹³C (between ?27‰ and ?23‰). Littoral area expansion due to siltation caused an increasing influx of terrestrial soil OM in the 1980s and the 1990s, subsequently lowering δ¹³C values and rising C/N ratios in both sediment cores. Meanwhile, higher δ¹⁵N was linked to a high influx of isotopically heavy nitrate from urban and agricultural wastewaters. After 2000, slight decreases in TOC and TN in the littoral area were attributable to reducing inputs of external OM, likely linked to declining sediment influx from the upper reaches resulting from the Three Gorges Dam impoundment. Contrasting increases in TOC, TN, and C/N ratios in the central part indicated a high influx of terrestrial soil OM due to the declining distance from the shoreline with littoral area expansion. Declining δ¹⁵N values after 2000 indicated an increase in N₂‐fixing cyanobacteria with eutrophication. Changes in diatom assemblages in both the littoral and central zones reflected nutrient enrichment and hydrological alterations. These results indicate that littoral expansion, declining riverine influx, and anthropogenic nutrient inputs are potential driving forces for the biogeochemical processing of OM in floodplain lakes. This study provides sedimentary biogeochemical clues for tracking past limnological conditions of floodplain lakes that are subjected to increasing disturbances from hydrological regulation and eutrophication.     

On the biogeochemical response of a glacierized High Arctic watershed to climate change: revealing patterns,processes and heterogeneity among micro‐catchments

Our novel study examines landscape biogeochemical evolution following deglaciation and permafrost change in Svalbard by looking at the productivity of various micro‐catchments existing within one watershed. It also sheds light on how moraine, talus and soil environments contribute to solute export from the entire watershed into the downstream marine ecosystem. We find that solute dynamics in different micro‐catchments are sensitive to abiotic factors such as runoff volume, water temperature, geology, geomorphological controls upon hydrological flowpaths and landscape evolution following sea level and glacial changes. Biotic factors influence the anionic composition of runoff because of the importance of microbial SO₄^2? and NO₃^? production. The legacy of glaciation and its impact upon sea level changes is shown to influence local hydrochemistry, allowing Cl^? to be used as a tracer of thawing permafrost that has marine origins. However, we show that a ‘glacial signal’ dominates solute export from the watershed. Therefore, although climatically driven change in the proglacial area has an influence on local ecosystems, the biogeochemical response of the entire watershed is dominated by glacially derived products of rapid chemical weathering. Consequently, only the study of micro‐catchments existing within watersheds can uncover the landscape response to contemporary climate change. Copyright © 2014 John Wiley & Sons, Ltd.     

Identifying surface water evaporation loss of inland river basin based on evaporation enrichment model

Accurately quantifying the evaporation loss of surface water is essential for regional water resources management, especially in arid and semi-arid areas where water resources are already scarce. The long-term monitoring of stable isotopes (δ¹⁸O and δ²H) in water can provide a sensitive indicator of water loss by evaporation. In this study, we obtained surface water samples of Shiyang River Basin from April to October between 2017 and 2019. The spatial and temporal characteristics of stable isotopes in surface water show the trend of enrichment in summer, depletion in spring, enrichment in deserts and depletion in mountains. The Local Evaporation Line (LEL) obtained by the regression of δ²H and δ¹⁸O in surface water has been defined by the lines: δ²H = 7.61δ¹⁸O + 14.58 for mountainous area, δ²H = 4.19δ¹⁸O − 17.85 for oasis area, δ²H = 4.08δ¹⁸O − 18.92 for desert area. The slope of LEL shows a gradual decrease from mountain to desert, indicating that the evaporation of surface water is gradually increasing. The evaporation loss of stable isotopes in surface water is 24.82% for mountainous area, 32.19% for oasis area, and 70.98% for desert area, respectively. Temperature and air humidity are the main meteorological factors affecting the evaporation loss, and the construction of reservoirs and farmland irrigation are the main man-made factors affecting the evaporation loss.     

Improved high-resolution global and regionalized isoscapes of δ18O, δ2H and d-excess in precipitation

Patterns of δ¹⁸O and δ²H in Earth's precipitation provide essential scientific data for use in hydrological, climatological, ecological and forensic research. Insufficient global spatial data coverage promulgated the use of gridded datasets employing geostatistical techniques (isoscapes) for spatiotemporally coherent isotope predictions. Cluster-based isoscape regionalization combines the advantages of local or regional prediction calibrations into a global framework. Here we present a revision of a Regionalized Cluster-Based Water Isotope Prediction model (RCWIP2) incorporating new isotope data having extensive spatial coverage and a wider array of predictor variables combined with high-resolution gridded climatic data. We introduced coupling of δ¹⁸O and δ²H (e.g., d-excess constrained) in the model predictions to prevent runaway isoscapes when each isotope is modelled separately and cross-checked observed versus modelled d-excess values. We improved model error quantification by adopting full uncertainty propagation in all calculations. RCWIP2 improved the RMSE over previous isoscape models by ca. 0.3 ‰ for δ¹⁸O and 2.5 ‰ for δ²H with an uncertainty <1.0 ‰ for δ¹⁸O and < 8 ‰ for δ²H for most regions of the world. The determination of the relative importance of each predictor variable in each ecoclimatic zone is a new approach to identify previously unrecognized climatic drivers on mean annual precipitation δ¹⁸O and δ²H. The improved RCWIP2 isoscape grids and maps (season, monthly, annual, regional) are available for download at https://isotopehydrologynetwork.iaea.org .     

Partitioning resistance to overland flow on rough mobile beds

For overland flows transporting predominantly bed load over rough mobile beds without rainfall, resistance to flow f may be divided into four components: surface resistance f_s, form resistance f_f, wave resistance f_w, and bed‐mobility resistance f_m. In this study it is assumed that f = f_s + f_f + f_w + f_m, and an equation is developed for each component. The equations for f_s and f_f are borrowed from the literature, while those for f_w and f_m are developed from two series of flume experiments in which the beds are covered with various concentrations of large‐scale roughness elements. The first series consists of 65 experiments on fixed beds, while the second series contains 194 experiments on mobile beds. All experiments were performed on the same slope (S = 0·114) and with the same size of sediment (D = 0·00074 m). The equations for f_w and f_m are derived by a combination of dimensional analysis and regression analysis. The analyses reveal that the major controls of f_w and f_m are the Froude number F and the concentration of the roughness elements C_r. When the equations for f_w and f_m are summed, the C_r terms cancel out, leaving f_w+m = 0·63F^?2. An equation is developed that predicts total f, and the contributions of f_s, f_f, f_w and f_m to f are computed from the series 1 and 2 experiments. An analysis of the first series reveals that in clear‐water flows over fixed beds, f_w accounts for 52 per cent of f. A similar analysis of the second series indicates that in sediment‐laden flows over mobile beds f_w comprises 37 per cent and f_m 32 per cent of f, so that together f_w and f_m account for almost 70 per cent of f. Finally, regression analyses indicate that where F > 0·5, f_w and f_m each vary with F^?2 and f_w/f_m = 1·18. The equation developed here for predicting total f applies only to the range of hydraulic, sediment, and bed roughness conditions represented by the experimental data. With additional data from a broader range of conditions the same methodology as employed here could be used to develop a more general equation. Copyright © 2006 John Wiley & Sons, Ltd.     

Continuous real‐time analysis of the isotopic composition of precipitation during tropical rain events: Insights into tropical convection

To investigate stable isotopic variability of precipitation in Singapore, we continuously analysed the δ‐value of individual rain events from November 2014 to August 2017 using an online system composed of a diffusion sampler coupled to Cavity Ring‐Down Spectrometer. Over this period, the average value (δ¹⁸O_Avg), the lowest value (δ¹⁸O_Low), and the initial value (δ¹⁸O_Init) varied significantly, ranging from ?0.45 to ?15.54‰, ?0.9 to ?17.65‰, and 0 to ?13.13‰, respectively. All 3 values share similar variability, and events with low δ¹⁸O_Low and δ¹⁸O_Avg values have low δ¹⁸O_Init value. Individual events have limited intraevent variability in δ‐value (Δδ) with the majority having a Δδ below 4‰. Correlation of δ¹⁸O_Low and δ¹⁸O_Avg with δ¹⁸O_Init is much higher than that with Δδ, suggesting that convective activities prior to events have more control over δ‐value than on‐site convective activities. The d‐excess of events also varies considerably in response to the seasonal variation in moisture sources. A 2‐month running mean analysis of δ¹⁸O reveals clear seasonal and interannual variability. Seasonal variability is associated with the meridional movement of the Intertropical Convergence Zone and evolution of the Asian monsoon. El Niño–Southern Oscillation is a likely driver of interannual variability. During 2015–2016, the strongest El Niño year in recorded history, the majority of events have a δ¹⁸O value higher than the weighted average δ¹⁸O of daily precipitation. δ¹⁸O shows a positive correlation with outgoing longwave radiation in the western Pacific and the Asian monsoon region, and also with Oceanic Niño Index. During El Niño, the convection centre shifts eastward to the central/eastern Pacific, weakening convective activities in Southeast Asia. Our study shows that precipitation δ‐value contains information about El Niño–Southern Oscillation and the Intertropical Convergence Zone, which has a significant implication for the interpretation of water isotope data and understanding of hydrological processes in tropical regions.     

Use of stable water isotopes to identify hydrological processes of meteoric water in montane catchments

This study analyzes the stable isotopic compositions of hydrogen and oxygen (δ²H, δ¹⁸O) in montane meteoric waters including precipitation and stream water of central Taiwan to identify hydrological processes in montane catchments. Results of precipitation demonstrate that monsoon and altitude effects are two principal processes affecting δ and deuterium excess (d^E) values of inland precipitation in central Taiwan. Furthermore, slope and intercept values of summer and winter local meteoric water line are modified by secondary evaporation effects such as moisture recycling and raindrop evaporation. Additionally, stream water's results indicate that differences in δ values among stream waters reflect isotopic altitude effect whereby lower values are more evident in stream water originating from high‐elevation catchments than low‐elevation catchments. Comparison of the isotopic results between precipitation and stream water indicates that summer precipitation containing recycled moisture is the most important water source for the studied stream waters and indicates that catchment effect and base flow contribution are the two major hydrological processes affecting mountain stream hydrology. The hydrological processes identified by the isotopic study re‐stress the important role of forests in mountain hydrology. Copyright © 2015 John Wiley & Sons, Ltd.     

Influence of below-cloud evaporation on stable isotopes of precipitation in the Yellow River source region

The study of below-cloud evaporation effects under clouds in the Yellow River source region is of great significance for regional water resource generation as well as for water resource security in the arid and semi-arid regions of northern China. In this study, we quantitatively assessed the evapotranspiration effect in the Yellow River source region from March to November based on the improved Stewart model. The study concluded that: (1) below-cloud evaporation was slightly higher in summer than in other seasons (residual fractions of raindrop evaporation were 80.57% in summer, 81.12% in spring, and 84.2% in autumn, respectively); and (2) sub-cloud evaporation diminishes with increasing altitude (residual fractions of raindrop evaporation were 83.09% in the western part of the area, 81.82% in the central part of the area, and 81.36% in the eastern part of the area, respectively). (3) The total linear index between study areas f and ∆d is 2.24, where f > 95%, it is 1.19; that is, the evaporation of raindrops increases by 1% and the reduction in the excess of mercury by about 2‰. (4) Local meteorological factors (temperature, precipitation, and relative humidity) and raindrop diameter have a cross-influence on below-cloud evaporation, with relative humidity having the most significant effect, with the highest correlation coefficient of 3.03 when relative humidity is less than 70%. The results of the study can provide a parameter basis for hydrological and climatic models in the Yellow River Basin.     

Effects of permafrost degradation on thermokarst lake hydrochemistry in the Qinghai-Tibet Plateau,China

Thermokarst lakes play a key role in the hydrological and biogeochemical cycles of permafrost regions. Current knowledge regarding the changes caused by permafrost degradation to the hydrochemistry of lakes in the Qinghai-Tibet Plateau (QTP) is limited. To address this gap, a systematic investigation of thermokarst lake water, suprapermafrost water, ground ice, and precipitation was conducted in the hinterland of the QTP. The thermokarst lake water in the QTP was identified to be of the Na-HCO₃-Cl type. The mean concentrations of HCO₃⁻ and Na⁺ were 281.8 mg L⁻¹ (146.0–546.2 mg L⁻¹) and 73.3 mg L⁻¹ (9.2–345.8 mg L⁻¹), respectively. The concentrations of Li⁺, NH₄⁺, K⁺, F⁻, NO₂⁻, and NO₃⁻ were relatively low. Freeze-out fractionation concentrated the dissolved solids within the lake water during winter, which was deeply deepened on lake depth and lake ice thickness. Owing to solute enrichment, the ground ice was characterized by high salinity. Conversely, repeated replenishment via precipitation led to lower solute concentrations in the ground ice near the permafrost table compared to that within the permafrost. Although lower solute concentration existed in precipitation, the soil leaching and saline ground ice melting processes enhanced the solute load in suprapermafrost water, which is considered an important water and solute resource in thermokarst lakes. The influencing mechanism of permafrost degradation on thermokarst lake hydrochemistry is presumably linked to: (1) the liberation of soluble materials sequestered in ground ice; (2) the increase of solutes in suprapermafrost water and soil pore water; and (3) the changes in lake morphometry. These results have major implications on the understanding of the effects of ground ice melting on ecosystem functions, biogeochemical processes, and energy balance in a rapidly changing climate.     

Evaluating sulfur dynamics during storm events for three watersheds in the northeastern USA: a combined hydrological,chemical and isotopic approach

Concerns related to climate change have resulted in an increasing interest in the importance of hydrological events such as droughts in affecting biogeochemical responses of watersheds. The effects of an unusually dry summer in 2002 had a marked impact on the biogeochemistry of three watersheds in the north‐eastern USA. Chemical, isotopic and hydrological responses with particular emphasis on S dynamics were evaluated for Archer Creek (New York), Sleepers River (Vermont) and Cone Pond (New Hampshire) watersheds. From 1 August to 14 September 2002, all three watersheds had very low precipitation (48 to 69 mm) resulting in either very low or no discharge (mean 0·015, 0·15 and 0·000 mm day^?1 for Archer Creek, Sleepers River and Cone Pond, respectively). From 15 September to 31 October 2002, there was a substantial increase in precipitation totals (212, 246 and 198 mm, respectively) with increased discharge. Archer Creek was characterized by a large range of SO₄^2? concentrations (152 to 389 µeq L^?1, mean = 273 µeq L^?1) and also exhibited the greatest range in δ³⁴S values of SO₄^2? (?1·4 to 8·8 ‰ ). Sleepers River's SO₄^2? concentrations ranged from 136 to 243 µeq L^?1 (mean = 167 µeq L^?1) and δ³⁴S values of SO₄^2? ranged from 4·0 to 9·0 ‰ . Cone Pond's SO₄^2? concentrations (126–187 µeq L^?1, mean = 154 µeq L^?1) and δ³⁴S values (2·4 to 4·3 ‰ ) had the smallest ranges of the three watersheds. The range and mean of δ¹⁸O‐SO₄^2? values for Archer Creek and Cone Pond were similar (3·0 to 8·9 ‰ , mean = 4·5 ‰ ; 3·9 to 6·3 ‰ , mean = 4·9 ‰ ; respectively) while δ¹⁸O‐SO₄^2? values for Sleepers River covered a larger range with a lower mean (1·2 to 10·0 ‰ , mean = 2·5). The difference in Sleepers River chemical and isotopic responses was attributed to weathering reactions contributing SO₄^2?. For Archer Creek wetland areas containing previously reduced S compounds that were reoxidized to SO₄^2? probably provided a substantial source of S. Cone Pond had limited internal S sources and less chemical or isotopic response to storms. Differences among the three watersheds in S biogeochemical responses during these storm events were attributed to differences in S mineral weathering contributions, hydrological pathways and landscape features. Further evaluations of differences and similarities in biogeochemical and hydrological responses among watersheds are needed to predict the impacts of climate change. Copyright © 2008 John Wiley & Sons, Ltd.     

Variable nitrate concentration–discharge relationships in a forested watershed

The relationship between solute concentrations and discharge can inform an integrated understanding of hydrological and biogeochemical processes at watershed scales. Recent work from multiple catchments has shown that there is typically little variation in concentration relative to large variations in discharge. This pattern has been described as chemostatic behavior. Pond Branch, a forested headwater catchment in Maryland, has been monitored for stream nitrate (NO₃^?) concentrations at weekly intervals for 14 years. In the growing season and autumn of 2011 a high‐frequency optical NO₃^? sensor was used to supplement the long‐term weekly data. In this watershed, long‐term weekly data show that NO₃^? concentrations decrease with increasing discharge whereas 6 months of 15‐minute sensor observed concentrations reveal a more chemostatic behavior. High‐frequency NO₃^? concentrations from the sensor collected during different storm events reveal variable concentration–discharge patterns highlighting the importance of high resolution data and ecohydrological drivers in controlling solute export for biologically reactive solutes such as NO₃^?.     

Spatial and temporal variability of stable isotopes (δ18O and δ2H) in surface waters of arid,mountainous Central Asia

Characterization of spatial and temporal variability of stable isotopes (δ¹⁸O and δ²H) of surface waters is essential to interpret hydrological processes and establish modern isotope–elevation gradients across mountainous terrains. Here, we present stable isotope data for river waters across Kyrgyzstan. River water isotopes exhibit substantial spatial heterogeneity among different watersheds in Kyrgyzstan. Higher river water isotope values were found mainly in the Issyk‐Kul Lake watershed, whereas waters in the Son‐Kul Lake watershed display lower values. Results show a close δ¹⁸O–δ²H relation between river water and the local meteoric water line, implying that river water experiences little evaporative enrichment. River water from the high‐elevation regions (e.g., Naryn and Son‐Kul Lake watershed) had the most negative isotope values, implying that river water is dominated by snowmelt. Higher deuterium excess (average d = 13.9‰) in river water probably represents the isotopic signature of combined contributions from direct precipitation and glacier melt in stream discharge across Kyrgyzstan. A significant relationship between river water δ¹⁸O and elevation was observed with a vertical lapse rate of 0.13‰/100 m. These findings provide crucial information about hydrological processes across Kyrgyzstan and contribute to a better understanding of the paleoclimate/elevation reconstruction of this region.     

Identification of groundwater recharge sources and processes in a heterogeneous alluvial aquifer: results from multi‐level monitoring of hydrochemistry and environmental isotopes in a riverside agricultural area in Korea

We evaluated sources and pathways of groundwater recharge for a heterogeneous alluvial aquifer beneath an agricultural field, based on multi‐level monitoring of hydrochemistry and environmental isotopes of a riverside groundwater system at Buyeo, Korea. Two distinct groundwater zones were identified with depth: (1) a shallow oxic groundwater zone, characterized by elevated concentrations of NO₃^? and (2) a deeper (>10–14 m from the ground surface) sub‐oxic groundwater zone with high concentrations of dissolved Fe, silica, and HCO₃^?, but little nitrate. The change of redox zones occurred at a depth where the aquifer sediments change from an upper sandy stratum to a silty stratum with mud caps. The δ¹⁸O and δ²H values of groundwater were also different between the two zones. Hydrochemical and δ¹⁸O? δ²H data of oxic groundwater are similar to those of soil water. This illustrates that recharge of oxic groundwater mainly occurs through direct infiltration of rain and irrigation water in the sandy soil area where vegetable cropping with abundant fertilizer use is predominant. Oxic groundwater is therefore severely contaminated by agrochemical pollutants such as nitrate. In contrast, deeper sub‐oxic groundwater contains only small amounts of dissolved oxygen (DO) and NO₃^?. The ³H contents and elevated silica concentrations in sub‐oxic groundwater indicate a somewhat longer mean residence time of groundwater within this part of the aquifer. Sub‐oxic groundwater was also characterized by higher δ¹⁸O and δ²H values and lower d‐excess values, indicating significant evaporation during recharge. We suggest that recharge of sub‐oxic groundwater occurs in the areas of paddy rice fields where standing irrigation and rain water are affected by strong evaporation, and that reducing conditions develop during subsequent sub‐surface infiltration. This study illustrates the existence of two groundwater bodies with different recharge processes within an alluvial aquifer. Copyright © 2009 John Wiley & Sons, Ltd.