Hydro-climate and biogeochemical processes control watershed organic carbon inflows: Development of an in-stream organic carbon module coupled with a process-based hydrologic model

Organic carbon and mercury exports from pan-Arctic rivers in a thawing permafrost context - A review

Climate change affects more than elsewhere the northern circumpolar permafrost region. This zone comprises large rivers flowing mainly to the Arctic Ocean, delivering about 10 % of the global riverine water flux. These pan-Arctic Rivers drive the dynamics of northern organic carbon (OC) and mercury (Hg) cycling. Permafrost degradation may release substantial amounts of OC and Hg, with potential regional and global impacts. In this review, we summarise the main findings in the last three decades about the role of the pan-Arctic Rivers in OC and Hg cycling and the effect of climate change on these dynamics. Total DOC and POC fluxes delivered by the pan-Arctic rivers presently reach 34.4 ± 1.2 TgC·yr-1 and 7.9 ± 0.5 TgC·yr-1, while the export of Hg reaches 38.9 ± 1.7 Mg·yr-1. This review highlights future challenges for the scientific community in evaluating spatial and temporal dynamics of the processes involved in OC and Hg cycling in permafrost-affected areas. Permafrost thawing could lead to greater fluxes of OC and Hg with ill-known resulting risks for food chains. Within this context, efforts should be made to study OC effects on Hg methylation. Moreover, assessing the spatial variability of OC and Hg mobilisation and transport within the pan-Arctic watersheds may help understand the future OC and Hg cycling dynamics in the northern circumpolar permafrost region.

Identification of surface water - groundwater nitrate governing factors in Jianghuai hilly area based on coupled SWAT-MODFLOW-RT3D modeling approach

A comprehensive understanding of the key controlling factors on NO3-N spatiotemporal distribution in surface and groundwater is of great significance to nitrogen pollution control and water resources management in watershed. Hence, the coupled SWAT-MODFLOW-RT3D model was employed to simulate nitrate (NO3-) fate and transport in Huashan watershed system. The model was calibrated using a combination of stream discharge, groundwater levels, NO3-N in-stream loading and groundwater NO3-N concentrations. The simulation revealed the significant spatiotemporal variations in surface water-groundwater nitrate interactions. The annual average percolation of NO3- from rivers to groundwater was 171.5 kg/km2 and the annual average discharge NO3- content from groundwater into rivers was 451.9 kg/km2 over the simulation period. The highest percolation of NO3- from rivers to groundwater occurred in April and the highest discharge NO3- content from groundwater into rivers occurred in July. Grassland and agriculture land contributed more nitrate contents in river water and groundwater compared to bare land and forest in the study area and the water exchange was the primary driving force for nitrate interactions in the surface water-groundwater system. Sensitivity analysis indicated that river runoff and groundwater levels were most influenced by the SCS runoff curve number f (CN2) and aquifer hydraulic conductivity (K), which, in turn, significantly affected nitrate transport. Regarding water quality parameters, the denitrification exponential rate coefficient (CDN) had the most pronounced impact on NO3-N in-stream loading and groundwater NO3-N concentrations. This study underscores the central role of surface-groundwater (SW-GW) interactions in watershed-scale nitrate research and suggests that parameters with higher sensitivity should be prioritized in analogous watershed modeling.

Assessing fluvial organic carbon flux and its response to short climate variability and damming on a large-scale tropical Asian river basin

Fluvial organic carbon (OC) transfer is an essential resource for downstream ecosystems. Multiple factors affect its transfer process, e.g., climate or anthropogenic activities. Quantifying OC fluxes with fine spatiotemporal resolution is challenging in anthropised catchments. This study aims to quantify daily OC dynamics and to assess the impacts of short climate variability and damming on OC spatiotemporal transfer processes in a large tropical Asian river basin (the Red River) for an extended period (2003-2013) by combining empirical equations with modelling outputs. Firstly, empirical equations for calculating dissolved (DOC) and particulate OC (POC) concentrations were calibrated based on in-situ sampling data. Then, simulated daily discharge (Q) and suspended sediment concentrations were used to quantify daily OC fluxes. Results show that the parameters of the DOC and POC equations well represent the subbasins characteristics, underlining the effects of soil OC content, mean annual Q and Chlorophyll a. DOC and POC exports reached 222 and 406 kt yr-1 at the basin outlet, accounting for 0.38 % of the total OC (TOC) exported by Asian rivers to the ocean. However, the specific yields of DOC (1.62 t km-2 yr-1) and POC (2.96 t km-2 yr-1) of the Red River basin were ~ 1.5 times those of other Asian basins. By comparing a reference scenario (without dams) to current conditions, we estimated 12 % and 88 % decreases in DOC and POC fluxes between 2008-2013 and 2003-2007, mainly due to damming. This study shows that climate variability may not impact OC dynamics in rivers as it explained <2 % of the variations. However, dam management, especially recent ones operating since 2008, deeply influences OC variations as the POC/TOC ratio decreased from 86 % to 47 %. Damming significantly decreased POC exports due to sediment retention, altering the equilibrium of OC cycling downstream, which may impact the food chain.

Roles of dissolved organic matter (DOM) in shaping the distribution pattern of heavy metal in the Yangtze River

Dissolved organic matter (DOM) strongly influences the solid-liquid partitioning and migration characteristics of heavy metals, yet little is known about the metal distribution and risk with the participation of DOM in large riverine systems. This study investigated the spatiotemporal distribution of 14 heavy metals and DOM along the entire Yangtze River (over 6000 km), and highlighted the critical roles of DOM in regulating the environmental behaviors of heavy metals. Significant spatial variations of metal contents were observed, with the river source and lower reach remarkably different from the upper-middle reaches. Heavy metals in the Yangtze River were mainly from the natural sources with minor anthropogenic disturbance. We found DOM could promote the conversion of metals from solid to liquid phase and DOM with higher aromaticity showed higher metal affinities. Although low ecological risks were observed in the Yangtze River, potential risks of metal leaching warrant attention, especially for As, Cd and Sb in the middle-lower reaches with higher DOM content and aromaticity. This study established a source-to-sea investigative approach to evaluate the influences of DOM features on metal partitioning, which is crucial for the risk control and sustainable management of large rivers.

Investigating hydrological transport pathways of dissolved organic carbon in cold region watershed based on a watershed biogeochemical model

Dissolved organic carbon (DOC) is a significant component of regional and global carbon cycles and an important surface water quality indicator. DOC affects the processes of solubility, bioavailability and transport for a number of contaminants, such as heavy metals. Therefore, it is crucial to understand DOC fate and transport in the watershed and the transport pathways of DOC load. We modified a previously developed watershed-scale organic carbon model by incorporating the DOC load from glacier melt runoff and used the modified model to simulate periodic daily DOC load in the upper Athabasca River Basin (ARB) in the cold region of western Canada. The calibrated model achieved an overall acceptable performance for simulating daily DOC load with model uncertainties mainly from the underestimation of peak loads. Parameter sensitivity analysis indicates that the fate and transport of DOC load in upper ARB are mainly controlled by DOC production in the soil layers, DOC transport at the soil surface, and reactions in the stream system. The modeling results indicated that the DOC load is mainly from the terrestrial sources and the stream system was a negligible sink in the upper ARB. It also indicated that rainfall-induced surface runoff was the major transport pathway of DOC load in the upper ARB. However, the DOC loads transported by glacier melt runoff were negligible and only accounted for 0.02% of the total DOC loads. In addition, snowmelt-induced surface runoff and lateral flow contributed 18.7% of total DOC load, which is comparable to the contribution from the groundwater flow. Our study investigated the DOC dynamics and sources in the cold region watershed in western Canada and quantified the contribution of different hydrological pathways to DOC load, which could provide a useful reference and insight for understanding watershed-scale carbon cycle processes.

Modelling Subarctic watershed dissolved organic carbon response to hydroclimatic regime

Shifts in hydroclimatic regimes associated with global climate change may impact freshwater availability and quality. In high latitudes of the northern hemisphere, where vast quantities of carbon are stored terrestrially, explaining landscape-scale carbon (C) budgets and associated pollutant transfer is necessary for understanding the impact of changing hydroclimatic regimes. We used a dynamic modelling approach to simulate streamflow, DOC concentration, and DOC export in a northern Canadian catchment that has undergone notable climate warming, and will continue to for the remainder of this century. The Integrated Catchment model for Carbon (INCA-C) was successfully calibrated to a multi-year period (2012-2016) that represents a range in hydrologic conditions. The model was subsequently run over 30-year periods representing baseline and two future climate scenarios. Average discharge is predicted to decrease under an elevated temperature scenario (22-27 % of baseline) but increase (116-175 % of baseline) under an elevated temperature and precipitation scenario. In the latter scenario the nival hydroclimatic regime is expected to shift to a combined nival and pluvial regime. Average DOC flux over 30 years is predicted to decrease (24-27 % of baseline) under the elevated temperature scenario, as higher DOC concentrations are offset by lower runoff. Under the elevated temperature and precipitation scenario, results suggest an increase in carbon export of 64-81 % above baseline. These increases are attributed to greater connectivity of the catchment. The largest increase in DOC export is expected to occur in early winter. These predicted changes in DOC export, particularly under a climate that is warmer and wetter could be part of larger ecosystem change and warrant additional monitoring efforts in the region.

Integrated water quality modeling in a river-reservoir system to support watershed management

Water resources planning and management are dependable on an adequate integration of physical, chemical, biological, and socio-economic realities of multiple water users. The dynamic of water quantity and quality in rivers are affected by several conditions, such as land use, soil characteristics, and meteorological and hydrological processes. Among these, the presence of hydraulic structures, such as dams and reservoirs, are often responsible for hydrodynamic and geomorphological alterations, leading to impacts on water quality and ecological behavior. In this context, this study presents the combination of a solution of the one-dimensional flow and transport and fate of contaminants in rivers (SihQual model) with a continuously stirred tank reactor (CSTR) to represent the reservoir. The first model solves the hydrodynamic and water quality equations under unsteady state, while the second approach considers the reservoir as a complete mixed system with inputs that vary over time. The main goal is to provide an integrated analysis for planning and management in a watershed where water has multiple purposes of use. The case study is the Iguaçu watershed, where the main river is affected by several dams for hydroelectric power generation. The simulation covers 542 km of the Iguaçu River and the Foz do Areia reservoir (flooded area of 139.5 km²), using data from 12 monitoring stations. The region also has issues of water quality impairment and water scarcity events due to deforestation, and urban and agricultural activities, exemplifying challenges throughout the world. Results show that the integrated models can reproduce the expected variability in different systems, although calibration challenges arise in multiscale modeling. The data indicate that, overall, the lentic environment is able to deplete organic matter and phosphorus, in comparison to levels in the fluvial flow. Nonetheless, experiments show that the river-reservoir system may be highly sensible to external and internal changes, such as water availability throughout time and pollution from the main tributary, as well as outlet discharges and transformation processes in the water column, leading to a possibility of critical events. Therefore, the study highlights how planning and managing actions in the watershed can benefit from an integrated river-reservoir analysis.

Anthropogenically driven climate and landscape change effects on inland water carbon dynamics: What have we learned and where are we going?

Inland waters serve as important hydrological connections between the terrestrial landscape and oceans but are often overlooked in global carbon (C) budgets and Earth System Models. Terrestrially derived C entering inland waters from the watershed can be transported to oceans but over 83% is either buried in sediments or emitted to the atmosphere before reaching oceans. Anthropogenic pressures such as climate and landscape changes are altering the magnitude of these C fluxes in inland waters. Here, we synthesize the most recent estimates of C fluxes and the differential contributions across inland waterbody types (rivers, streams, lakes, reservoirs, and ponds), including recent measurements that incorporate improved sampling methods, small waterbodies, and dried areas. Across all inland waters, we report a global C emission estimate of 4.40 Pg C/year (95% confidence interval: 3.95-4.85 Pg C/year), representing a 13% increase from the most recent estimate. We also review the mechanisms by which the most globally widespread anthropogenically driven climate and landscape changes influence inland water C fluxes. The majority of these drivers are expected to influence terrestrial C inputs to inland waters due to alterations in terrestrial C quality and quantity, hydrological pathways, and biogeochemical processing. We recommend four research priorities for the future study of anthropogenic alterations to inland water C fluxes: (1) before-and-after measurements of C fluxes associated with climate change events and landscape changes, (2) better quantification of C input from land, (3) improved assessment of spatial coverage and contributions of small inland waterbodies to C fluxes, and (4) integration of dried and drawdown areas to global C flux estimates. Improved measurements of inland water C fluxes and quantification of uncertainty in these estimates will be vital to understanding both terrestrial C losses and the "moving target" of inland water C emissions in response to rapid and complex anthropogenic pressures.

Non-stationary response of rain-fed spring wheat yield to future climate change in northern latitudes

The non-stationary response of crop growth to changes in hydro-climatic variables makes yield projection uncertain and the design and implementation of adaptation strategies debatable. This study simulated the time-varying behavior of the underlying cause-and-effect mechanisms affecting spring wheat yield (SWY) under various climate change and nitrogen (N) application scenarios in the Red Deer River basin in agricultural lands of the western Canadian Prairies. A calibrated and validated Soil and Water Assessment Tool and Analysis of Variance decomposition methods were utilized to assess the contribution of crop growth parameters, Global Climate Models, Representative Concentration Pathways, and downscaling techniques to the total SWY variance for the 2040-2064 period. The results showed that the cause-and-effect mechanisms, driving crop yield, shifted from water stress (W-stress) dominated (27 days of W-stress days) during the historical period to nitrogen stress (N-stress) dominated (27 to 35 N-stress days) in the future period. It was shown that while higher precipitation, warmer weather, and early snowmelts, along with elevated CO₂ may favor SWY in cold regions in the future (up to 50% more yields in some sub-basins), the yield potentials may be limited by N-stress (only up to 0.7% yield increase in some sub-basins). The N-stress might be partially related to the N deficiency in the soil, which can be compensated by N fertilizer application. However, inadequate N uptake due to limited evapotranspiration under elevated atmospheric CO₂ might pose restrictions to SWY potentials even in the least N deficient regions. This study uncovers important information on the understanding of spatiotemporal variability of hydrogeochemical processes driving crop yields and the non-stationary response of yields to changing climate. The results also underscore spatiotemporal variability of N-stress due to N deficiency in the soil or N uptake by crops, both of which may restrain SWY by changes in atmospheric CO₂ concentrations in the future.

Bacterial Communities in Stream Biofilms in a Degrading Grassland Watershed on the Qinghai-Tibet Plateau

Grassland is among the largest terrestrial biomes and is experiencing serious degradation, especially on the Qinghai-Tibet Plateau (QTP). However, the influences of grassland degradation on microbial communities in stream biofilms are largely unknown. Using 16S rRNA gene sequencing, we investigated the bacterial communities in stream biofilms in sub-basins with different grassland status in the Qinghai Lake watershed. Grassland status in the sub-basins was quantified using the normalized difference vegetation index (NDVI). Proteobacteria, Bacteroidetes, Cyanobacteria, and Verrucomicrobia were the dominant bacterial phyla. OTUs, 7,050, were detected in total, within which 19 were abundant taxa, and 6,922 were rare taxa. Chao 1, the number of observed OTUs, and phylogenetic diversity had positive correlations with carbon (C), nitrogen (N), and/or phosphorus (P) in biofilms per se. The variation of bacterial communities in stream biofilms was closely associated with the rate of change in NDVI, pH, conductivity, as well as C, N, P, contents and C:N ratio of the biofilms. Abundant subcommunities were more influenced by environmental variables relative to the whole community and to rare subcommunities. These results suggest that the history of grassland degradation (indicated as the rate of change in NDVI) influences bacterial communities in stream biofilms. Moreover, the bacterial community network showed high modularity with five major modules (>50 nodes) that responded differently to environmental variables. According to the module structure, only one module connector and 12 module hubs were identified, suggesting high fragmentation of the network and considerable independence of the modules. Most of the keystone taxa were rare taxa, consistent with fragmentation of the network and with adverse consequences for bacterial community integrity and function in the biofilms. By documenting the properties of bacterial communities in stream biofilms in a degrading grassland watershed, our study adds to our knowledge of the potential influences of grassland degradation on aquatic ecosystems.

Impacts of Hydrological Processes on Stream Temperature in a Cold Region Watershed Based on the SWAT Equilibrium Temperature Model

Variance in stream temperature from historical norms, which reflects the impacts from both hydrological and meteorological factors, is a significant indicator of the stream ecosystem health. Therefore, it is imperative to study the hydrological processes controlling stream temperature in the watershed. The impacts of hydrological processes on stream temperature in the cold region of Western Canada were investigated based on the previously developed Soil and Water Assessment Tool (SWAT) equilibrium temperature model. The model was calibrated and validated for streamflow and stream temperature based on the observations and a global parameter sensitivity analysis conducted to identify the most important hydrological process governing the stream temperature dynamics. The precipitation and air temperature lapse rates were found to be the most sensitive parameters controlling the stream temperature, followed by the parameters regulating the processes of soil water dynamics, surface runoff, and channel routing. Our analysis showed an inverse relationship between streamflow volume and stream temperature, and different runoff components have different impacts on temporal regimes of stream temperatures. This study elaborates on the response of the stream temperature to changes in hydrological processes at the watershed scale and indicates that hydrological processes should be taken into account for prediction of stream temperatures.