A preliminary assessment of water partitioning and ecohydrological coupling in northern headwaters using stable isotopes and conceptual runoff models

Funded by European Research Council ERC. Grant Number: GA 335910 VEWA Swedish Science Foundation (SITES) Future Forest Formas (ForWater) SKB the Kempe foundation Environment Canada the Garfield Weston Foundation the Natural Sciences and Engineering Research Council of Canada (NSERC) the Northwest Territories Cumulative Impacts Monitoring ProgramPeer reviewedPublisher PD

Current forest carbon fixation fuels stream CO2 emissions

Stream CO2 emissions contribute significantly to atmospheric climate forcing. While there are strong indications that groundwater inputs sustain these emissions, the specific biogeochemical pathways and timescales involved in this lateral CO2 export are still obscure. Here, via an extensive radiocarbon (14C) characterisation of CO2 and DOC in stream water and its groundwater sources in an old-growth boreal forest, we demonstrate that the 14C-CO2 is consistently in tune with the current atmospheric 14C-CO2 level and shows little association with the 14C-DOC in the same waters. Our findings thus indicate that stream CO2 emissions act as a shortcut that returns CO2 recently fixed by the forest vegetation to the atmosphere. Our results expose a positive feedback mechanism within the C budget of forested catchments, where stream CO2 emissions will be highly sensitive to changes in forest C allocation patterns associated with climate and land-use changes

Reading tea leaves worldwide: decoupled drivers of initial litter decomposition mass‐loss rate and stabilization

The breakdown of plant material fuels soil functioning and biodiversity. Currently, process understanding of global decomposition patterns and the drivers of such patterns are hampered by the lack of coherent large‐scale datasets. We buried 36,000 individual litterbags (tea bags) worldwide and found an overall negative correlation between initial mass‐loss rates and stabilization factors of plant‐derived carbon, using the Tea Bag Index (TBI). The stabilization factor quantifies the degree to which easy‐to‐degrade components accumulate during early‐stage decomposition (e.g. by environmental limitations). However, agriculture and an interaction between moisture and temperature led to a decoupling between initial mass‐loss rates and stabilization, notably in colder locations. Using TBI improved mass‐loss estimates of natural litter compared to models that ignored stabilization. Ignoring the transformation of dead plant material to more recalcitrant substances during early‐stage decomposition, and the environmental control of this transformation, could overestimate carbon losses during early decomposition in carbon cycle models

Biomass offsets little or none of permafrost carbon release from soils, streams, and wildfire : an expert assessment

As the permafrost region warms, its large organic carbon pool will be increasingly vulnerable to decomposition, combustion, and hydrologic export. Models predict that some portion of this release will be offset by increased production of Arctic and boreal biomass; however, the lack of robust estimates of net carbon balance increases the risk of further overshooting international emissions targets. Precise empirical or model-based assessments of the critical factors driving carbon balance are unlikely in the near future, so to address this gap, we present estimates from 98 permafrost-region experts of the response of biomass, wildfire, and hydrologic carbon flux to climate change. Results suggest that contrary to model projections, total permafrost-region biomass could decrease due to water stress and disturbance, factors that are not adequately incorporated in current models. Assessments indicate that end-of-the-century organic carbon release from Arctic rivers and collapsing coastlines could increase by 75% while carbon loss via burning could increase four-fold. Experts identified water balance, shifts in vegetation community, and permafrost degradation as the key sources of uncertainty in predicting future system response. In combination with previous findings, results suggest the permafrost region will become a carbon source to the atmosphere by 2100 regardless of warming scenario but that 65%-85% of permafrost carbon release can still be avoided if human emissions are actively reduced.Peer reviewe

Spatial heterogeneity of the spring flood acid pulse in a boreal stream network

Spatial and temporal patterns in streamwater acidity are ecologically important, but difficult to measure in parallel. Here we present the spatial distribution of streamwater chemistry relevant to acidity from 60 stream sites distributed throughout a 67 km2 boreal catchment, sampled during a period of winter baseflow (high pH) and during a spring flood episode (low pH). Sites were grouped based on pH level and pH change from winter baseflow to spring flood. The site attributes of each pH group were then assessed in terms of both stream chemistry and subcatchment landscape characteristics. Winter baseflow pH was high throughout most of the stream network (median pH 6.4), but during the spring flood episode stream sites experienced declines in pH ranging from 0–1.6 pH units, resulting in pH ranging from 4.3–6.3. Spring flood pH was highest in larger, lower altitude catchments underlain by fine sorted sediments, and lowest in small, higher altitude catchments with a mixture of peat wetlands and forested till. Wetland-dominated headwater catchments had low but stable pH, while the spring flood pH drop was largest in a group of catchments of intermediate size which contained well-developed coniferous forest and a moderate proportion of peat wetlands. There was a trend with distance downstream of higher pH, acid neutralizing capacity (ANC) and base cation concentrations together with lower dissolved organic carbon (DOC, strongly negatively correlated with pH). This apparent scale- dependence of stream chemistry could be explained by a number of environmental factors which vary predictably with altitude, catchment area and distance downstream— most notably, a shift in surficial sediment type from unsorted till and peat wetlands to fine sorted sediments at lower altitudes in this catchment. As a result of the combination of spatial heterogeneity in landscape characteristics and scale-related processes, boreal catchments like this one can be expected to experience high spatial variability both in terms of chemistry at any given point in time, and in the change experienced during high discharge episodes. Although chemistry patterns showed associations with landscape characteristics, considerable additional variability remained, suggesting that the modeling of dynamic stream chemistry from map parameters will continue to present a challenge

The role of landscape properties, storage and evapotranspiration on variability in streamflow recessions in a boreal catchment

Streamflow recession analysis provides valuable insights into catchment functioning that can be related to runoff generation, storage retention and baseflow dynamics. As an integrated characteristic, recession analysis is particularly useful in catchment comparison studies to help explain drivers of spatial and temporal variability in hydrological behavior. Here, five years of hourly streamflow data from 14, partly nested, catchments within a 68 km2 boreal forest landscape in Northern Sweden were used to explore spatiotemporal variation in hydrological processes through recession analysis. The aim of this study was to better understand spatial variation in runoff generation and storage-discharge dynamics across the landscape, as well as the relation to landscape properties. Due to high collinearity between variables, partial least square regression was used to quantify the associations between recession characteristics and catchment properties, as well as to identify key variables controlling recession behavior. We analyzed recession characteristics using both an aggregated approach including all recession data and individual recession events. The analyses based on individual recession events, indicated that catchment topography, quantified by indices such as mean slope or elevation above the stream network, is a primary control on the recession behavior during relatively high flows, whereas catchment area gains importance when flows are relatively low. The proportion of sediment and deep soils also controlled recession behavior. Furthermore, we found that recession characteristics are influenced by both evapotranspiration (ET) and proxies of antecedent catchment storage, but that the patterns were different depending on catchment properties. ET was less influential in catchments with deeper soils and larger catchment area. Shifts in recession rates were primarily related to variation in storage, with faster streamflow recessions occurring during periods with low storage. The results demonstrate the influence of catchment properties on recession behavior, and we found great value in analyzing individual recession events for an increased understanding of spatial and temporal recession characteristics. When recession properties were lumped together, the relationships to catchment characteristics were obscured. This indicates the value of more detailed analyses, at least under the strongly seasonal hydroclimatic conditions of this site

Specific discharge variability in a boreal landscape

Specific discharge variations within a mesoscale catchment were studied on the basis of three synoptic sampling campaigns. These were conducted during stable flow conditions within the Krycklan catchment study area in northern Sweden. During each campaign, about 80 individual locations were measured for discharge draining from catchment areas ranging between 0.12 and 67 km2. These discharge samplings allowed for the comparison between years within a given season (September 2005 versus September 2008) and between seasons within a given year (May 2008 versus September 2008) of specific discharge across this boreal landscape. There was considerable variability in specific discharge across this landscape. The ratio of the interquartile range (IQR) defined as the difference between the 75th and 25th percentiles of the specific discharges to the median of the specific discharges ranged from 37% to 43%. Factor analysis was used to explore potential relations between landscape characteristics and the specific discharge observed for 55 of the individual locations that were measured in all three synoptic sampling campaigns. Percentage wet area (i.e., wetlands, mires, and lakes) and elevation were found to be directly related to the specific discharge during the drier September 2008 sampling while potential annual evaporation was found to be inversely related. There was less of a relationship determined during the wetter post spring flood May 2008 sampling and the late summer rewetted September 2005 sampling. These results indicate the ability of forests to “dry out” parts of the catchment over the summer months while wetlands “keep wet” other parts. To demonstrate the biogeochemical implications of such spatiotemporal variations in specific discharge, we estimate dissolved organic carbon (DOC) exports with available data for the May 2008 and September 2008 samplings using both the spatially variable observed specific discharges and the spatially constant catchment average values. The average absolute difference in DOC export for the various subcatchments between using a variable and using a constant specific discharge was 28% for the May 2008 sampling and 20% for the September 2008 sampling