Tree growth potential and its relationship with soil moisture conditions across a heterogeneous boreal forest landscape.

Forest growth varies across landscapes due to the intricate relationships between various environmental drivers and forest management. In this study, we analysed the variation of tree growth potential across a landscape scale and its relation to soil moisture. We hypothesised that soil moisture conditions drive landscape-level variation in site quality and that intermediate soil moisture conditions demonstrate the highest potential forest production. We used an age-independent difference model to estimate site quality in terms of maximum achievable tree height by measuring the relative change in Lorey's mean height for a five year period across 337 plots within a 68 km2 boreal landscape. We achieved wall-to-wall estimates of site quality by extrapolating the modelled relationship using repeated airborne laser scanning data collected in connection to the field surveys. We found a clear decrease in site quality under the highest soil moisture conditions. However, intermediate soil moisture conditions did not demonstrate clear site quality differences; this is most likely a result of the nature of the modelled soil moisture conditions and limitations connected to the site quality estimation. There was considerable unexplained variation in the modelled site quality both on the plot and landscape levels. We successfully demonstrated that there is a significant relationship between soil moisture conditions and site quality despite limitations associated with a short study period in a low productive region and the precision of airborne laser scanning measurements of mean height.

A drained nutrient-poor peatland forest in boreal Sweden constitutes a net carbon sink after integrating terrestrial and aquatic fluxes.

Northern peatlands provide a globally important carbon (C) store. Since the beginning of the 20th century, however, large areas of natural peatlands have been drained for biomass production across Fennoscandia. Today, drained peatland forests constitute a common feature of the managed boreal landscape, yet their ecosystem C balance and associated climate impact are not well understood, particularly within the nutrient-poor boreal region. In this study, we estimated the net ecosystem carbon balance (NECB) from a nutrient-poor drained peatland forest and an adjacent natural mire in northern Sweden by integrating terrestrial carbon dioxide (CO2 ) and methane (CH4 ) fluxes with aquatic losses of dissolved organic C (DOC) and inorganic C based on eddy covariance and stream discharge measurements, respectively, over two hydrological years. Since the forest included a dense spruce-birch area and a sparse pine area, we were able to further evaluate the effect of contrasting forest structure on the NECB and component fluxes. We found that the drained peatland forest was a net C sink with a 2-year mean NECB of -115 ± 5 g C m-2 year-1 while the adjacent mire was close to C neutral with 14.6 ± 1.7 g C m-2 year-1 . The NECB of the drained peatland forest was dominated by the net CO2 exchange (net ecosystem exchange [NEE]), whereas NEE and DOC export fluxes contributed equally to the mire NECB. We further found that the C sink strength in the sparse pine forest area (-153 ± 8 g C m-2 year-1 ) was about 1.5 times as high as in the dense spruce-birch forest area (-95 ± 8 g C m-2 year-1 ) due to enhanced C uptake by ground vegetation and lower DOC export. Our study suggests that historically drained peatland forests in nutrient-poor boreal regions may provide a significant net ecosystem C sink and associated climate benefits.

Consequences of rewetting and ditch cleaning on hydrology, water quality and greenhouse gas balance in a drained northern landscape.

Drainage for forestry has created ~ 1 million km of artificial waterways in Sweden, making it one of the largest human-induced environmental disturbances in the country. These extensive modifications of both peatland and mineral soil dominated landscapes still carry largely unknown, but potentially enormous environmental legacy effects. However, the consequences of contemporary ditch management strategies, such as hydrological restoration via ditch blocking or enhancing forest drainage to promote biomass production via ditch cleaning, on water resources and greenhouse gas (GHG) fluxes are unclear. To close the gap between science and management, we have developed a unique field research platform to experimentally evaluate key environmental strategies for drained northern landscapes with the aim to avoid further environmental degeneration. The Trollberget Experimental Area (TEA) includes replicated and controlled treatments applied at the catchment scale based on a BACI approach (before-after and control-impact). The treatments represent the dominant ecosystem types impacted by ditching in Sweden and the boreal zone: (1) rewetting of a drained peatland, (2) ditch cleaning in productive upland forests and (3) leaving these ditches unmanaged. Here we describe the TEA platform, report initial results, suggest ways forward for how to best manage this historical large-scale alteration of the boreal landscape, as well as warn against applying these treatments broadly before more long-term results are reported.

Soil moisture controls the partitioning of carbon stocks across a managed boreal forest landscape.

Boreal forests sequester and store vast carbon (C) pools that may be subject to significant feedback effects induced by climatic warming. The boreal landscape consists of a mosaic of forests and peatlands with wide variation in total C stocks, making it important to understand the factors controlling C pool sizes in different ecosystems. We therefore quantified the total C stocks in the organic layer, mineral soil, and tree biomass in 430 plots across a 68 km2 boreal catchment. The organic layer held the largest C pool, accounting for 39% of the total C storage; tree and mineral C pools accounted for 38% and 23%, respectively. The size of the soil C pool was positively related to modelled soil moisture conditions, especially in the organic soil layer (R2 = 0.50). Conversely, the tree C pool exhibited a unimodal relationship: storage was highest under intermediate wetness conditions. The magnitude and variation in the total soil C stocks observed in this work were comparable to those found at the national level in Sweden, suggesting that C accumulation in boreal landscapes is more sensitive to local variation resulting primarily from differences in soil moisture conditions than to regional differences in climate, nitrogen deposition, and parent material.

Leveraging research infrastructure co-location to evaluate constraints on terrestrial carbon cycling in northern European forests.

Integrated long-term, in-situ observations are needed to document ongoing environmental change, to "ground-truth" remote sensing and model outputs and to predict future Earth system behaviour. The scientific and societal value of in-situ observations increases with site representativeness, temporal duration, number of parameters measured and comparability within and across sites. Research Infrastructures (RIs) can support harmonised, cross-site data collection, curation and publication. Integrating RI networks through site co-location and standardised observation methods can help answers three questions about the terrestrial carbon sink: (i) What are present and future carbon sequestration rates in northern European forests? (ii) How are these rates controlled? (iii) Why do the observed patterns exist? Here, we present a conceptual model for RI co-location and highlight potential insights into the terrestrial carbon sink achievable when long-term in-situ Earth observation sites participate in multiple RI networks (e.g., ICOS and eLTER). Finally, we offer recommendations to promote RI co-location.

On the uncertainty in estimates of the carbon balance recovery time after forest clear-cutting.

An essential metric for describing carbon dynamics in managed forest landscapes is the recovery time of the carbon balance after clear-cutting. Here, we demonstrate how the age-dependent mathematical trajectory is affected by both the selected model and data availability, leading to considerable uncertainty in the modelling of the net ecosystem production (NEP) over stand age. We further show that the initial carbon loss estimates associated with the timing of the source-sink transition (SST) are significant, but may have a limited effect on the total carbon sequestration at the end of the standard (RP, 120 years) or optimal (OCS) rotation periods.

Seasonal trends and retention of polycyclic aromatic compounds (PACs) in a remote sub-Arctic catchment.

Atmospheric deposition is the most dominant source of polycyclic aromatic hydrocarbons (PAHs) in remote and pristine areas. Despite low bioaccumulation potential, PAHs and their persistent transformation products (PAH-derivatives) are chemicals of concern as they can harm human and animal health through chronic low dose exposure. In this study, atmospheric deposition fluxes of polycyclic aromatic compounds (PACs) were measured on a seasonal basis (3-month periods) from 2012 to 2016 in a remote subarctic forest catchment in northern Europe. The target PACs included 19 PAHs and 15 PAH-derivatives (oxygenated, nitrogenated, and methylated PAHs). The deposition fluxes of Æ©PAHs and Æ©PAH-derivatives were in the same range and averaged 530 and 500 ng m2 day-1, respectively. The fluxes were found to be higher with a factor of 2.5 for Æ©PAHs and a factor of 3 for Æ©PAH-derivatives during cold (<0 °C) in comparison to warm (>10 °C) periods. PAHs and PAH-derivatives showed similar seasonal patterns, which suggests that these two compound classes have similar sources and deposition mechanisms, and that the source strength of the PAH-derivatives in air follows that of the PAHs. The terrestrial export of PACs via the outlet of the catchment stream was estimated to be 1.1% for Æ©PAHs and 1.7% for Æ©PAH-derivatives in relation to the annual amounts deposited to the catchment, which suggests that boreal forests are sinks for PACs derived from the atmosphere. Some individual PACs showed higher export than others (i.e. chrysene, cyclopenta(c,d)pyrene, carbazole, quinoline, and benzo(f)quinoline), with high export PACs mostly characterized by low molecular weight and low hydrophobicity (2-3 benzene rings; log Kow<6.0).

Landscape-variability of the carbon balance across managed boreal forests.

Boreal forests are important global carbon (C) sinks and, therefore, considered as a key element in climate change mitigation policies. However, their actual C sink strength is uncertain and under debate, particularly for the actively managed forests in the boreal regions of Fennoscandia. In this study, we use an extensive set of biometric- and chamber-based C flux data collected in 50 forest stands (ranging from 5 to 211 years) over 3 years (2016-2018) with the aim to explore the variations of the annual net ecosystem production (NEP; i.e., the ecosystem C balance) across a 68 km2 managed boreal forest landscape in northern Sweden. Our results demonstrate that net primary production rather than heterotrophic respiration regulated the spatio-temporal variations of NEP across the heterogeneous mosaic of the managed boreal forest landscape. We further find divergent successional patterns of NEP in our managed forests relative to naturally regenerating boreal forests, including (i) a fast recovery of the C sink function within the first decade after harvest due to the rapid establishment of a productive understory layer and (ii) a sustained C sink in old stands (131-211 years). We estimate that the rotation period for optimum C sequestration extends to 138 years, which over multiple rotations results in a long-term C sequestration rate of 86.5 t C ha-1 per rotation. Our study highlights the potential of forest management to maximize C sequestration of boreal forest landscapes and associate climate change mitigation effects by developing strategies that optimize tree biomass production rather than heterotrophic soil C emissions.

The emerging role of drought as a regulator of dissolved organic carbon in boreal landscapes.

One likely consequence of global climate change is an increased frequency and intensity of droughts at high latitudes. Here we use a 17-year record from 13 nested boreal streams to examine direct and lagged effects of summer drought on the quantity and quality of dissolved organic carbon (DOC) inputs from catchment soils. Protracted periods of drought reduced DOC concentrations in all catchments but also led to large stream DOC pulses upon rewetting. Concurrent changes in DOC optical properties and chemical character suggest that seasonal drying and rewetting trigger soil processes that alter the forms of carbon supplied to streams. Contrary to expectations, clearest drought effects were observed in larger watersheds, whereas responses were most muted in smaller, peatland-dominated catchments. Collectively, our results indicate that summer drought causes a fundamental shift in the seasonal distribution of DOC concentrations and character, which together operate as primary controls over the ecological and biogeochemical functioning of northern aquatic ecosystems.

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.

Current water quality guidelines across North America and Europe do not protect lakes from salinization.

Human-induced salinization caused by the use of road deicing salts, agricultural practices, mining operations, and climate change is a major threat to the biodiversity and functioning of freshwater ecosystems. Yet, it is unclear if freshwater ecosystems are protected from salinization by current water quality guidelines. Leveraging an experimental network of land-based and in-lake mesocosms across North America and Europe, we tested how salinization-indicated as elevated chloride (Cl-) concentration-will affect lake food webs and if two of the lowest Cl- thresholds found globally are sufficient to protect these food webs. Our results indicated that salinization will cause substantial zooplankton mortality at the lowest Cl- thresholds established in Canada (120 mg Cl-/L) and the United States (230 mg Cl-/L) and throughout Europe where Cl- thresholds are generally higher. For instance, at 73% of our study sites, Cl- concentrations that caused a ≥50% reduction in cladoceran abundance were at or below Cl- thresholds in Canada, in the United States, and throughout Europe. Similar trends occurred for copepod and rotifer zooplankton. The loss of zooplankton triggered a cascading effect causing an increase in phytoplankton biomass at 47% of study sites. Such changes in lake food webs could alter nutrient cycling and water clarity and trigger declines in fish production. Current Cl- thresholds across North America and Europe clearly do not adequately protect lake food webs. Water quality guidelines should be developed where they do not exist, and there is an urgent need to reassess existing guidelines to protect lake ecosystems from human-induced salinization.

The role of terrestrial productivity and hydrology in regulating aquatic dissolved organic carbon concentrations in boreal catchments.

The past decades have witnessed an increase in dissolved organic carbon (DOC) concentrations in the catchments of the Northern Hemisphere. Increasing terrestrial productivity and changing hydrology may be reasons for the increases in DOC concentration. The aim of this study is to investigate the impacts of increased terrestrial productivity and changed hydrology following climate change on DOC concentrations. We tested and quantified the effects of gross primary production (GPP), ecosystem respiration (RE) and discharge on DOC concentrations in boreal catchments over 3 years. As catchment characteristics can regulate the extent of rising DOC concentrations caused by the regional or global environmental changes, we selected four catchments with different sizes (small, medium and large) and landscapes (forest, mire and forest-mire mixed). We applied multiple models: Wavelet coherence analysis detected the delay-effects of terrestrial productivity and discharge on aquatic DOC variations of boreal catchments; thereafter, the distributed-lag linear models quantified the contributions of each factor on DOC variations. Our results showed that the combined impacts of terrestrial productivity and discharge explained 62% of aquatic DOC variations on average across all sites, whereas discharge, gross primary production (GPP) and RE accounted for 26%, 22% and 3%, respectively. The impact of GPP and discharge on DOC changes was directly related to catchment size: GPP dominated DOC fluctuations in small catchments (<1 km2 ), whereas discharge controlled DOC variations in big catchments (>1 km2 ). The direction of the relation between GPP and discharge on DOC varied. Increasing RE always made a positive contribution to DOC concentration. This study reveals that climate change-induced terrestrial greening and shifting hydrology change the DOC export from terrestrial to aquatic ecosystems. The work improves our mechanistic understanding of surface water DOC regulation in boreal catchments and confirms the importance of DOC fluxes in regulating ecosystem C budgets.

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.

How tree species, tree size, and topographical location influenced tree transpiration in northern boreal forests during the historic 2018 drought.

Trees in northern latitude ecosystems are projected to experience increasing drought stress as a result of rising air temperatures and changes in precipitation patterns in northern latitude ecosystems. However, most drought-related studies on high-latitude boreal forests (>50°N) have been conducted in North America, with few studies quantifying the response in European and Eurasian boreal forests. Here, we tested how daily whole-tree transpiration (Q, Liters day-1 ) and Q normalized for mean daytime vapor pressure deficit (QDZ , Liters day-1 kPa-1 ) were affected by the historic 2018 drought in Europe. More specifically, we examined how tree species, size, and topographic position affected drought response in high-latitude mature boreal forest trees. We monitored 30 Pinus sylvestris (pine) and 30 Picea abies (spruce) trees distributed across a topographic gradient in northern Sweden. In general, pine showed a greater QDZ control compared to spruce during periods of severe drought (standardized precipitation-evapotranspiration index: SPEI < -1.5), suggesting that the latter are more sensitive to drought. Overall, QDZ reductions (using non-drought QDZ as reference) were less pronounced in larger trees during severe drought, but there was a species-specific pattern: QDZ reductions were greater in pine trees at high elevations and greater in spruce trees at lower elevations. Despite lower QDZ during severe drought, drought spells were interspersed with small precipitation events and overcast conditions, and QDZ returned to pre-drought conditions relatively quickly. This study highlights unique species-specific responses to drought, which are additionally driven by a codependent interaction among tree size, relative topographic position, and unique regional climate conditions.

The undetected loss of aged carbon from boreal mineral soils.

The boreal forest is among the largest terrestrial biomes on earth, storing more carbon (C) than the atmosphere. Due to rapid climatic warming and enhanced human development, the boreal region may have begun transitioning from a net C sink to a net source. This raises serious concern that old biogenic soil C can be re-introduced into the modern C cycle in near future. Combining bio-decay experiments, mixing models and the Keeling plot method, we discovered a distinct old pre-bomb organic carbon fraction with high biodegradation rate. In total, 34 ± 12% of water-extractable organic carbon (WEOC) in podzols, one of the dominating boreal soil types, consisted of aged (~ 1000 year) labile C. The omission of this aged (i.e., Δ14C depleted) WEOC fraction in earlier studies is due to the co-occurrence with Δ14C enriched modern C formed following 1950s nuclear bomb testing masking its existence. High lability of aged soil WEOC and masking effects of modern Δ14C enriched C suggests that the risk for mobilization and re-introduction of this ancient C pool into the modern C cycle has gone undetected. Our findings have important implications for earth systems models in terms of climate-carbon feedbacks and the future C balance of the boreal forest.

Assessment of a portable UV-Vis spectrophotometer's performance in remote areas: Stream water DOC, Fe content and spectral data.

This paper presents data for the assessment of a portable UV-Vis spectrophotometer's performance on predicting stream water DOC and Fe content. The dataset contains DOC and Fe concentrations by laboratory methods, in-situ and ex-situ spectral absorbances, monitoring environmental indexes such as water depth, temperature, turbidity and voltage. The records in Yli-Nuortti river (Cold station, Finland) took place during the hydrological year 2018-2019 and in Krycklan (C4 and C5, Sweden) during the hydrological years 2016-2019. The data analyses were conducted with 'pls' and 'caret' package in R. The correlation coefficient (R), root-mean-square deviation (RMSD), standard deviation (STD) and bias were used to check the performance of the models. This dataset can be combined with datasets from other regions around the world to build more universal models. For discussion and more information of the dataset creation, please refer to the full-length article "Assessment of a portable UV-Vis spectrophotometer's performance for stream water DOC and Fe content monitoring in remote areas" [1].

Peatland drainage - a missing link behind increasing TOC concentrations in waters from high latitude forest catchments?

Total Organic Carbon (TOC) concentrations in stream waters from peat-covered catchments have increased over the last 15-25 years, resulting in large-scale brownification of lakes and rivers in high latitudes. While this increase has primarily been attributed to decreased acid deposition and climatic warming in most regions, we studied whether peatland drainage in forested catchments has contributed to the increasing TOC concentrations. We analysed the spatial variability of average TOC concentrations from a total of 133 peatland dominated catchments in Sweden and Finland, of which 62 were pristine and 71 were drained during the last century. In addition, we performed a trend analysis on 37 catchments for which long-term data were available. We found about 14 mg l-1 higher TOC concentrations in streams discharging from drained than undrained sites in southern latitudes, and about 8 mg l-1 higher concentrations from drained sites in northern latitudes. Trend analysis did not indicate significant differences in TOC concentration trends between drained and undrained catchments but indicated that tree stand volume correlated with increasing trends. This supports earlier findings in that the general increase in forest cover and biomass that has occurred in high latitudes during the last decades is another factor that has contributed to brownification.

Functionally reversible impacts of disturbances on lake food webs linked to spatial and seasonal dependencies.

Increasing human impact on the environment is causing drastic changes in disturbance regimes and how they prevail over time. Of increasing relevance is to further our understanding on biological responses to pulse disturbances (short duration) and how they interact with other ongoing press disturbances (constantly present). Because the temporal and spatial contexts of single experiments often limit our ability to generalize results across space and time, we conducted a modularized mesocosm experiment replicated in space (five lakes along a latitudinal gradient in Scandinavia) and time (two seasons, spring and summer) to generate general predictions on how the functioning and composition of multitrophic plankton communities (zoo-, phyto- and bacterioplankton) respond to pulse disturbances acting either in isolation or combined with press disturbances. As pulse disturbance, we used short-term changes in fish presence, and as press disturbance, we addressed the ongoing reduction in light availability caused by increased cloudiness and lake browning in many boreal and subarctic lakes. First, our results show that the top-down pulse disturbance had the strongest effects on both functioning and composition of the three trophic levels across sites and seasons, with signs for interactive impacts with the bottom-up press disturbance on phytoplankton communities. Second, community composition responses to disturbances were highly divergent between lakes and seasons: temporal accumulated community turnover of the same trophic level either increased (destabilization) or decreased (stabilization) in response to the disturbances compared to control conditions. Third, we found functional recovery from the pulse disturbances to be frequent at the end of most experiments. In a broader context, these results demonstrate that top-down, pulse disturbances, either alone or with additional constant stress upon primary producers caused by bottom-up disturbances, can induce profound but often functionally reversible changes across multiple trophic levels, which are strongly linked to spatial and temporal context dependencies. Furthermore, the identified dichotomy of disturbance effects on the turnover in community composition demonstrates the potential of disturbances to either stabilize or destabilize biodiversity patterns over time across a wide range of environmental conditions.

Assessment of a portable UV-Vis spectrophotometer's performance for stream water DOC and Fe content monitoring in remote areas.

Quantification of dissolved organic carbon (DOC) and iron (Fe) in surface waters is critical for understanding the water quality dynamics, brownification and carbon balance in the northern hemisphere. Especially in the remote areas, sampling and laboratory analysis of DOC and Fe content at a sufficient temporal frequency is difficult. Ultraviolet-visible (UV-Vis) spectrophotometry is a promising tool for water quality monitoring to increase the sampling frequency and applications in remote regions. The aim of this study was (1) to investigate the performance of an in-situ UV-Vis spectrophotometer for detecting spectral absorbances in comparison with a laboratory benchtop instrument; (2) to analyse the stability of DOC and Fe estimates from UV-Vis spectrophotometers among different rivers using multivariate methods; (3) to compare site-specific calibration of models to pooled models and investigate the extrapolation of DOC and Fe predictions from one catchment to another. This study indicates that absorbances that were measured by UV-Vis sensor explained 96% of the absorbance data from the laboratory benchtop instrument. Among the three tested multivariate methods, multiple stepwise regression (MSR) was the best model for both DOC and Fe predictions. Accurate and unbiased models for multiple watersheds for DOC were built successfully, and these models could be extrapolated from one watershed to another even without site-specific calibration for DOC. However, for Fe the combination of different datasets was not possible.