Groundwater travel times predict DOC in streams and riparian soils across a heterogeneous boreal landscape.

Dissolved organic carbon (DOC) in surface waters is an important component of the boreal landscape carbon budget and a critical variable in water quality. A dominant terrestrial DOC source in the boreal landscape is the riparian zone. These near stream areas play a key role in regulating DOC transport between land and aquatic ecosystems. The groundwater dynamics at this interface have been considered a major controlling variable for DOC export to streams. This study focuses on the regulating role of groundwater levels and mean travel times (MTT) on riparian DOC concentrations and, subsequently, stream DOC. This is done by comparing them as explanatory variables to capture the spatial and intra-annual variability of the stream and riparian groundwater DOC. We used a physically based 3D hydrological model, Mike SHE, to simulate DOC concentrations of the riparian zones for 14 sub-catchments within the Krycklan catchment (Sweden). The model concept assumes that DOC concentrations will be higher in groundwater moving through shallow flow paths. In the model, this can be linked to the position of the groundwater table at a point of observation or the travel time, which will generally be shorter for water that has travelled through shallow and more conductive soil layers. We compared the results with both observed stream and groundwater concentrations. The analysis revealed that the correlation between modelled and observed annual averages of stream DOC increased from r = 0.08 to r = 0.87 by using MTT instead of groundwater level. MTT also better captured the observed spatial variability in riparian DOC concentrations and more successfully represented seasonal variability of stream DOC. We, therefore, suggest that MTT is a better predictor than groundwater level for riparian DOC concentration because it can capture a greater variety of catchment heterogeneities, such as variation in soil properties, catchment size, and input from deep groundwater sources.

Scaling relations reveal global and regional differences in morphometry of reservoirs and natural lakes.

Water bodies provide essential ecosystem services linked to morphometric features that might differ between natural lakes and reservoirs. We use the HydroLAKES global dataset to quantitatively compare large (area > 1 km2) reservoirs and natural lakes in terms of scaling exponents between morphometric measures (volume, area, shore length). These exponents are further compared to those expected from geometrical assumptions and constraints. Lakes cover a larger range of volumes for the same range of surface areas than reservoirs, and have a larger volume-area scaling exponent. The volume-area scaling exponent for reservoirs (but not natural lakes) and the area-shore length exponent for all water bodies follow the predictions for self-affine surfaces. Land cover and terrain influence the scaling relations more for lakes than for reservoirs. These morphometric differences may be used to model the impact of reservoirs and lakes on hydrological processes and associated ecosystem services at regional to global scales.

Hydrological control of water quality - Modelling base cation weathering and dynamics across heterogeneous boreal catchments.

Linking biogeochemical processes to water flow paths and solute travel times is important for understanding internal catchment functioning and control of water quality. Base cation weathering is a process closely linked to key factors affecting catchment functioning, including water pathways, soil contact time, and catchment characteristics, particularly in silicate-dominated areas. However, common process-based weathering models are often calibrated and applied for individual soil profiles, which can cause problems when trying to extrapolate results to catchment scale and assess consequences for stream water and groundwater quality. Therefore, in this work, base cation export was instead modelled using a fully calibrated 3D hydrological model (Mike SHE) of a boreal catchment, which was expanded by adding a relatively simple but still reasonably flexible and versatile weathering module including the base cations Na, K, Mg, and Ca. The results were evaluated using a comprehensive dataset of water chemistry from groundwater and stream water in 14 nested sub-catchments, representing different catchment sizes and catchment characteristics. The strongest correlations with annual and seasonal observations were found for Ca (r = 0.89-0.93, p < 0.05), Mg (r = 0.90-0.95, p < 0.05), and Na (r = 0.80-0.89, p < 0.05). These strong correlations suggest that catchment hydrology and landscape properties primarily control weathering rates and stream dynamics of these solutes. Furthermore, catchment export of Mg, Ca, and K was strongly connected to travel times of discharging stream water (r = 0.78-0.83). Conversely, increasing Na export was linked to a reduced areal proportion of mires (r = -0.79). The results suggest that a significant part (~45%) of the catchment stream export came from deep-soil weathering sources (>2.5 m). These results have implications for terrestrial and aquatic water quality assessments. If deep soils are present, focusing mainly on the shallow soil could lead to misrepresentation of base cation availability and the acidification sensitivity of groundwater and water recipients such as streams and lakes.