Search

Soil moisture plays a key role in soil nutrient and carbon cycling; plant productivity; and energy, water, and greenhouse gas exchanges between the land and the atmosphere. The knowledge on drivers of spatiotemporal soil moisture dynamics in subarctic landscapes is limited. In this study, we used the Spatial Forest Hydrology (SpaFHy) model, in situ soil moisture data, and Sentinel-1 synthetic aperture radar (SAR)-based soil moisture estimates to explore spatiotemporal controls of soil moisture in a subarctic headwater catchment in northwestern Finland. The role of groundwater dynamics and lateral flow in soil moisture was studied through three groundwater model conceptualizations: (i) omission of groundwater storage and lateral flow, (ii) conceptual TOPMODEL approach based on topographic wetness index, and (iii) explicit 2D lateral groundwater flow. The model simulations were compared against continuous point soil moisture measurements, distributed manual measurements, and novel SAR-based soil moisture estimates available at high spatial and temporal resolutions. Based on model scenarios and model–data comparisons, we assessed when and where the lateral groundwater flow shapes shallow soil moisture and under which conditions soil moisture variability is driven more by local ecohydrology, i.e., the balance of infiltration, drainage, and evapotranspiration. The choice of groundwater flow model was shown to have a strong impact on modeled soil moisture dynamics within the catchment. All model conceptualizations captured the observed soil moisture dynamics in the upland forests, but accounting for the lateral groundwater flow was necessary to reproduce the saturated conditions common in the peatlands and occasionally in lowland forest grid cells. We further highlight the potential of integrating multi-scale observations with land surface and hydrological models. The results have implications for ecohydrological and biogeochemical processes, as well as for modeling hydrology and Earth system feedbacks in subarctic and boreal environments.

As soil microbial respiration is the major component of land CO2 emissions, differences in the functional dependence of respiration on soil moisture among Earth system models (ESMs) contributes significantly to the uncertainties in their projections. Using soil organic C (SOC) stocks and CO2 data from a boreal forest–mire ecotone in Finland and Bayesian data assimilation, we revised the original precipitation-based monotonic saturation dependency of the Yasso07 soil carbon model using the non-monotonic Ricker function based on soil volumetric water content. We fit the revised functional dependency of moisture to the observed microbial respiration and SOC stocks and compared its performance against the original Yasso07 model and the version used in the JSBACH land surface model with a reduction constant for decomposition rates in wetlands. The Yasso07 soil C model coupled with the calibrated unimodal Ricker moisture function with an optimum in well-drained soils accurately reconstructed observed SOC stocks and soil CO2 emissions and clearly outperformed previous model versions on paludified organo-mineral soils in forested peatlands and water-saturated organic soils in mires. The best estimate of the posterior moisture response of decomposition used both measurements of SOC stocks and CO2 data from the full range of moisture conditions (from dry and xeric to wet and water-saturated soils). We observed unbiased residuals of SOC and CO2 data modelled with the moisture optimum in well-drained soils, suggesting that this modified function accounts more precisely for the long-term SOC change dependency according to ecosystem properties as well as the contribution of short-term CO2 responses including extreme events. The optimum moisture for decomposition in boreal forests was found in well-drained soils instead of the mid-range between dry and water-saturated conditions as is commonly assumed among soil C and ESMs. Although the unimodal moisture modifier with an optimum in well-drained soils implicitly incorporates robust biogeochemical mechanisms of SOC accumulation and CO2 emissions, it needs further evaluation with large-scale data to determine if its use in land surface models will decrease the uncertainty in projections.

The snowpack has a major influence on the land surface energy budget. Accurate simulation of the snowpack energy and radiation budget is challenging due to, e.g., effects of vegetation and topography, as well as limitations in the theoretical understanding of turbulent transfer in the stable boundary layer. Studies that evaluate snow, hydrology and land surface models against detailed observations of all surface energy balance components at high latitudes are scarce. In this study, we compared different configurations of the SURFEX land surface model against surface energy flux, snow depth and soil temperature observations from four eddy-covariance stations in Finland. The sites cover two different climate and snow conditions, representing the southern and northern subarctic zones, as well as the contrasting forest and peatland ecosystems typical for the boreal landscape. We tested different turbulent flux parameterizations implemented in the Crocus snowpack model. In addition, we examined common alternative approaches to conceptualize soil and vegetation, and we assessed their performance in simulating surface energy fluxes, snow conditions and soil thermal regime. Our results show that a stability correction function that increases the turbulent exchange under stable atmospheric conditions is imperative to simulate sensible heat fluxes over the peatland snowpacks and that realistic peat soil texture (soil organic content) parameterization greatly improves the soil temperature simulations. For accurate simulations of surface energy fluxes, snow and soil conditions in forests, an explicit vegetation representation is necessary. Moreover, we demonstrate the high sensitivity of surface fluxes to a poorly documented parameter involved in snow cover fraction computation. Although we focused on models within the SURFEX platform, the results have broader implications for choosing suitable turbulent flux parameterization and model structures depending on the potential use cases for high-latitude land surface modeling.

Pristine boreal mires are known as substantial sinks of carbon dioxide (CO2) and net emitters of methane (CH4). Bogs constitute a major fraction of pristine boreal mires. However, the bog CO2 and CH4 balances are poorly known, having been largely estimated based on discrete and short-term measurements by manual chambers and seldom using the eddy-covariance (EC) technique. Eddy-covariance (EC) measurements of CO2 and CH4 exchange were conducted in the Siikaneva mire complex in southern Finland in 2011–2016. The site is a patterned bog having a moss–sedge–shrub vegetation typical of southern Eurasian taiga, with several ponds near the EC tower. The study presents a complete series of CO2 and CH4 EC flux (FCH4) measurements and identifies the environmental factors controlling the ecosystem–atmosphere CO2 and CH4 exchange. A 6-year average growing season (May–September) cumulative CO2 exchange of −61 ± 24 g C m−2 was observed, which partitions into mean total respiration (Re) of 167 ± 33 (interannual range 146–197) g C m−2 and mean gross primary production (GPP) of 228 ± 46 (interannual range 193–257) g C m−2, while the corresponding FCH4 amounts to 7.1 ± 0.7 (interannual range 6.4–8.4) g C m−2. The contribution of October–December CO2 and CH4 fluxes to the cumulative sums was not negligible based on the measurements during one winter. GPP, Re and FCH4 increased with temperature. GPP and FCH4 did not show any significant decline even after a substantial water table drawdown in 2011. Instead, GPP, Re and FCH4 were limited in the cool, cloudy and wet growing season of 2012. May–September cumulative net ecosystem exchange (NEE) of 2013–2016 averaged at about −73 g C m−2, in contrast with the hot and dry year 2011 and the wet and cool year 2012. Suboptimal weather likely reduced the net sink by about 25 g C m−2 in 2011 due to elevated Re, and by about 40 g C m−2 in 2012 due to limited GPP. The cumulative growing season sums of GPP and CH4 emission showed a strong positive relationship. The EC source area was found to be comprised of eight distinct surface types. However, footprint analyses revealed that contributions of different surface types varied only within 10 %–20 % with respect to wind direction and stability conditions. Consequently, no clear link between CO2 and CH4 fluxes and the EC footprint composition was found despite the apparent variation in fluxes with wind direction.

High-resolution thermal infrared (TIR) imaging is opening up new vistas in biosphere–atmosphere heat exchange studies. The rapidly developing unmanned aerial systems (UASs) and specially designed cameras offer opportunities for TIR survey with increasingly high resolution, reduced geometric and radiometric noise, and prolonged flight times. A state-of-the-art science platform is assembled using a Matrice 210 V2 drone equipped with a Zenmuse XT2 thermal camera and deployed over a pristine boreal peatland with the aim of testing its performance in a heterogeneous sedge-fen ecosystem. The study utilizes the capability of the UAS platform to hover for prolonged times (about 20 min) at a height of 500 m a.g.l. while recording high frame rate (30 Hz) TIR videos of an area of ca. 430 × 340 m. A methodology is developed to derive thermal signatures of near-ground coherent turbulent structures impinging on the land surface, surface temperature spectra, and heat fluxes from the retrieved videos. The size, orientation, and movement of the coherent structures are computed from the surface temperature maps, and their dependency on atmospheric conditions is examined. A range of spectral and wavelet-based approaches are used to infer the properties of the dominant turbulent scene structures. A ground-based eddy-covariance system and an in situ meteorological setup are used for reference.

Although it is well known that evapotranspiration (ET) represents an important water flux at local to global scales, few studies have quantified the magnitude and relative importance of ET and its individual flux components in high-latitude forests. In this study, we combined empirical sapflux, throughfall, and eddy-covariance measurements with estimates from a process-based model to partition the water balance in a northern boreal forested catchment. This study was conducted within the Krycklan catchment, which has a rich history of hydrological measurements, thereby providing us with the unique opportunity to compare the absolute and relative magnitudes of ET and its flux components to other water balance components. During the growing season, ET represented ca. 85 % of the incoming precipitation. Both empirical results and model estimates suggested that tree transpiration (T) and evaporation of intercepted water from the tree canopy (IC) represented 43 % and 31 % of ET, respectively, and together were equal to ca. 70 % of incoming precipitation during the growing season. Understory evapotranspiration (ETu) was less important than T and IC during most of the study period, except for late autumn, when ETu was the largest ET flux component. Overall, our study highlights the importance of trees in regulating the water cycle of boreal catchments, implying that forest management impacts on stand structure as well as climate change effects on tree growth are likely to have large cascading effects on the way water moves through these forested landscapes.

Vegetation is known to have strong influence on evapotranspiration (ET), a major component of terrestrial water balance. Yet hydrological models often describe ET by methods unable to include the variability of vegetation characteristics in their predictions. To take advantage of the increasing availability of high-resolution open GIS data on land use, vegetation and soil characteristics in the boreal zone, a modular, spatially distributed model for predicting ET and other hydrological processes from grid cell to catchment level is presented and validated. An improved approach to upscale stomatal conductance to canopy scale using information on plant type (conifer/deciduous) and stand leaf-area index (LAI) is proposed by coupling a common leaf-scale stomatal conductance model with a simple canopy radiation transfer scheme. Further, a generic parametrization for vegetation-related hydrological processes for Nordic boreal forests is derived based on literature and data from a boreal FluxNet site. With the generic parametrization, the model was shown to reproduce daily ET measured using an eddy-covariance technique well at 10 conifer-dominated Nordic forests whose LAI ranged from 0.2 to 6.8 m2 m−2. Topography, soil and vegetation properties at 21 small boreal headwater catchments in Finland were derived from open GIS data at 16 m × 16 m grid size to upscale water balance from stand to catchment level. The predictions of annual ET and specific discharge were successful in all catchments, located from 60 to 68∘ N, and daily discharge was also reasonably well predicted by calibrating only one parameter against discharge measurements. The role of vegetation heterogeneity in soil moisture and partitioning of ET was demonstrated. The proposed framework can support, for example, forest trafficability forecasting and predicting impacts of climate change and forest management on stand and catchment water balance. With appropriate parametrization it can be generalized outside the boreal coniferous forests.

Alkylamines are important precursors in secondary aerosol formation in the boreal forest atmosphere. To better understand the behavior and sources of two alkylamines, dimethylamine (DMA) and diethylamine (DEA), we estimated the magnitudes of soil–atmosphere fluxes of DMA and DEA using a gradient-diffusion approximation based on measured concentrations in soil solution and in the canopy air space. The ambient air concentration of DMA used in this study was a sum of DMA and ethylamine. To compute the amine fluxes, we first estimated the soil air space concentration from the measured soil solution amine concentration using soil physical (temperature, soil water content) and chemical (pH) state variables. Then, we used the resistance analogy to account for gas transport mechanisms in the soil, soil boundary layer, and canopy air space. The resulting flux estimates revealed that the boreal forest soil with a typical long-term mean pH 5.3 is a possible source of DMA (170 ± 51 nmol m−2 day−1) and a sink of DEA (−1.2 ± 1.2 nmol m−2 day−1). We also investigated the potential role of fungi as a reservoir for alkylamines in boreal forest soil. We found high DMA and DEA concentrations both in fungal hyphae collected from field humus samples and in fungal pure cultures. The highest DMA and DEA concentrations were found in fungal strains belonging to decay and ectomycorrhizal fungal groups, indicating that boreal forest soil and, in particular, fungal biomass may be important reservoirs for these alkylamines.

Canopy structure is one of the most important vegetation characteristics for land–atmosphere interactions, as it determines the energy and scalar exchanges between the land surface and the overlying air mass. In this study we evaluated the performance of a newly developed multi-layer energy budget in the ORCHIDEE-CAN v1.0 land surface model (Organising Carbon and Hydrology In Dynamic Ecosystems – CANopy), which simulates canopy structure and can be coupled to an atmospheric model using an implicit coupling procedure. We aim to provide a set of acceptable parameter values for a range of forest types. Top-canopy and sub-canopy flux observations from eight sites were collected in order to conduct this evaluation. The sites crossed climate zones from temperate to boreal and the vegetation types included deciduous, evergreen broad-leaved and evergreen needle-leaved forest with a maximum leaf area index (LAI; all-sided) ranging from 3.5 to 7.0. The parametrization approach proposed in this study was based on three selected physical processes – namely the diffusion, advection, and turbulent mixing within the canopy. Short-term sub-canopy observations and long-term surface fluxes were used to calibrate the parameters in the sub-canopy radiation, turbulence, and resistance modules with an automatic tuning process. The multi-layer model was found to capture the dynamics of sub-canopy turbulence, temperature, and energy fluxes. The performance of the new multi-layer model was further compared against the existing single-layer model. Although the multi-layer model simulation results showed few or no improvements to both the nighttime energy balance and energy partitioning during winter compared with a single-layer model simulation, the increased model complexity does provide a more detailed description of the canopy micrometeorology of various forest types. The multi-layer model links to potential future environmental and ecological studies such as the assessment of in-canopy species vulnerability to climate change, the climate effects of disturbance intensities and frequencies, and the consequences of biogenic volatile organic compound (BVOC) emissions from the terrestrial ecosystem.

Drainage intensity and forest management in peatlands affect carbon dioxide (CO 2 ) emissions to the atmosphere and export of dissolved organic carbon (DOC) to water courses. The peatland carbon (C) balance results from a complex network of ecosystem processes from where lateral C fluxes have typically been ignored. Here, we present a new version of the SUSI Peatland simulator, the first advanced process-based ecosystem model that compiles a full C balance in drained forested peatland including DOC formation, transport and biodegradation. SUSI considers site, stand and terrain characteristics as well as the interactions and feedbacks between ecosystem processes and offers novel ways to evaluate and mitigate adverse environmental impacts with thorough management planning. Here, we extended SUSI by designing and parameterizing a mass-balance based decomposition module (ESOM) based on literature findings and tested the ESOM performance against an independent dataset measured in the laboratory using peat columns collected from Finland, Estonia, Sweden and Ireland. ESOM predicted the CO 2 emissions and changes in DOC concentrations with a reasonable accuracy for the peat columns. We applied the new SUSI for drained peatland sites and found that reducing the depth to which ditches are cleaned by 0.3 m decreased the annual DOC export by 34 (17 %), 29 (19 %) and 7 (5 %) kg ha -1 in Finland, Estonia and Sweden, respectively, using typical ditch spacing for these countries. Correspondingly, site annual C sink increased by 305, 409 and 32 kg ha -1 in Finland, Estonia and Sweden, respectively. Our results also indicated that terrain slope can markedly alter the water residence time and consequently DOC biodegradation and export to ditches. We conclude that DOC export can be decreased and site C sink increased by reducing the depth to which ditches are cleaned or by increasing the ditch spacing., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Encroachment of vascular plants (VP) in temperate raised bogs, as a consequence of altered hydrological conditions and nutrient input, is widely observed. Effects of such vegetation shift on water and carbon cycles are, however, largely unknown and identification of responsible plant physiological traits is challenging. Process-based modelling offers the opportunity of gaining insights into ecosystem functioning beyond observations, and to infer decisive trait shifts of plant functional groups. We adapted the Soil-Vegetation-Atmosphere Transfer model pyAPES to a temperate raised bog site by calibration against measured peat temperature, water table and surface CO 2 fluxes. We identified the most important traits determining CO 2 fluxes by conducting Morris sensitivity analysis (MSA) under changing conditions throughout the year and simulated VP encroachment. We further investigated transferability of results to other sites by extending MSA to parameter ranges derived from literature review. We found highly variable intra-annual plant traits importance determining ecosystem CO 2 fluxes, but only a partial shift of importance of photosynthetic processes from moss to VP during encroachment. Ecosystem respiration was dominated by peat respiration. Overall, carboxylation rate, base respiration rate and temperature sensitivity (Q 10 ) were most important for determining bog CO 2 balance and parameter ranking was robust even under the extended MSA., (© 2024. The Author(s).)

Continuous cover forestry (CCF) has gained interest as an alternative to even-aged management particularly on drained peatland forests. However, relatively little is known about the physiological response of suppressed trees when larger trees are removed as a part of CCF practices. Consequently, studies concentrating on process-level modeling of the response of trees to selection harvesting are also rare. Here, we compared, modeled and measured harvest response of previously suppressed Norway spruce (Picea abies) trees to a selection harvest. We quantified the harvest response by collecting Norway spruce tree-ring samples in a drained peatland forest site and measuring the change in stable carbon and oxygen isotopic ratios of wood formed during 2010-20, including five post-harvest years. The measured isotopic ratios were compared with ecosystem-level process model predictions for ${\kern0em }^{13}$C discrimination and ${\kern0em }^{18}$O leaf water enrichment. We found that the model predicted similar but lower harvest response than the measurements. Furthermore, accounting for mesophyll conductance was important for capturing the variation in ${\kern0em }^{13}$C discrimination. In addition, we performed sensitivity analysis on the model, which suggests that the modeled ${\kern0em }^{13}$C discrimination is sensitive to parameters related to CO2 transport through stomata to the mesophyll., (© The Author(s) 2023. Published by Oxford University Press.)

The transport processes of methane (CH 4 ) in tree stems remain largely unknown, although they are critical in assessing the whole-forest CH4 dynamics. We used a physically based dynamic model to study the spatial and diurnal dynamics of stem CH 4 transport and fluxes. We parameterised the model using data from laboratory experiments with Pinus sylvestris and Betula pendula and compared the model to experimental data from a field study. Stem CH 4 flux in laboratory and field conditions were explained by the axial advective CH 4 transport from soil with xylem sap flow and the radial CH 4 diffusion through the stem conditions. Diffusion resistance caused by the bark permeability did not significantly affect gas transport or stem CH 4 flux in the laboratory experiments. The role of axial diffusion of CH 4 in trees was unresolved and requires further studies. Due to the transit time of CH 4 in the stem, the diurnal dynamics of stem CH 4 fluxes can deviate markedly from the diurnal dynamics of sap flow., (© 2023 The Authors. Plant, Cell & Environment published by John Wiley & Sons Ltd.)

Land-based mitigation measures are needed to achieve climate targets. One option is the mitigation of currently high greenhouse gas (GHG) emissions of nutrient-rich drained peatland forest soils. Continuous cover forestry (CCF) has been proposed as a measure to manage this GHG emission source; however, its emission reduction potential and impact on timber production at regional and national scales have not been quantified. To quantify the potential emission reduction, we simulated four management scenarios for Finnish forests: (i) The replacement of clear-cutting by selection harvesting on nutrient-rich drained peatlands (CCF) and (ii) the current forest management regime (BAU), and both at two harvest levels, namely (i) the mean annual harvesting (2016-2018) and (ii) the maximum sustainable yield. The simulations were conducted at the stand scale with a forest simulator (MELA) coupled with a hydrological model (SpaFHy), soil C model (Yasso07) and empirical GHG exchange models. Simulations showed that the management scenario that avoided clear-cutting on nutrient-rich drained peatlands (i.e. CCF) produced approximately 1 Tg CO 2 eq. higher carbon sinks annually compared with BAU at equal harvest level for Finland. This emission reduction can be attributed to the maintenance of a higher biomass sink and to the mitigation of soil emissions from nutrient-rich drained peatland sites., (© 2023. Springer Nature Limited.)

Nonvascular photoautotrophs (NVP), including bryophytes, lichens, terrestrial algae, and cyanobacteria, are increasingly recognized as being essential to ecosystem functioning in many regions of the world. Current research suggests that climate change may pose a substantial threat to NVP, but the extent to which this will affect the associated ecosystem functions and services is highly uncertain. Here, we propose a research agenda to address this urgent question, focusing on physiological and ecological processes that link NVP to ecosystem functions while also taking into account the substantial taxonomic diversity across multiple ecosystem types. Accordingly, we developed a new categorization scheme, based on microclimatic gradients, which simplifies the high physiological and morphological diversity of NVP and world-wide distribution with respect to several broad habitat types. We found that habitat-specific ecosystem functions of NVP will likely be substantially affected by climate change, and more quantitative process understanding is required on: (1) potential for acclimation; (2) response to elevated CO 2 ; (3) role of the microbiome; and (4) feedback to (micro)climate. We suggest an integrative approach of innovative, multimethod laboratory and field experiments and ecophysiological modelling, for which sustained scientific collaboration on NVP research will be essential., (© 2022 The Authors New Phytologist © 2022 New Phytologist Foundation.)

We explore needle sugar isotopic compositions (δ 18 O and δ 13 C) in boreal Scots pine (Pinus sylvestris) over two growing seasons. A leaf-level dynamic model driven by environmental conditions and based on current understanding of isotope fractionation processes was built to predict δ 18 O and δ 13 C of two hierarchical needle carbohydrate pools, accounting for the needle sugar pool size and the presence of an invariant pinitol pool. Model results agreed well with observed needle water δ 18 O, δ 18 O and δ 13 C of needle water-soluble carbohydrates (sugars + pinitol), and needle sugar δ 13 C (R 2  = 0.95, 0.84, 0.60, 0.73, respectively). Relative humidity (RH) and intercellular to ambient CO 2 concentration ratio (C i /C a ) were the dominant drivers of δ 18 O and δ 13 C variability, respectively. However, the variability of needle sugar δ 18 O and δ 13 C was reduced on diel and intra-seasonal timescales, compared to predictions based on instantaneous RH and C i /C a , due to the large needle sugar pool, which caused the signal formation period to vary seasonally from 2 d to more than 5 d. Furthermore, accounting for a temperature-sensitive biochemical 18 O-fractionation factor and mesophyll resistance in 13 C-discrimination were critical. Interpreting leaf-level isotopic signals requires understanding on time integration caused by mixing in the needle sugar pool., (© 2022 The Authors. New Phytologist © 2022 New Phytologist Foundation.)

The terrestrial net ecosystem productivity (NEP) has increased during the past three decades, but the mechanisms responsible are still unclear. We analyzed 17 years (2001-2017) of eddy-covariance measurements of NEP, evapotranspiration (ET) and light and water use efficiency from a boreal coniferous forest in Southern Finland for trends and inter-annual variability (IAV). The forest was a mean annual carbon sink (252 [ ± 42] gC m - 2 a - 1 ), and NEP increased at rate +6.4-7.0 gC m - 2 a - 1 (or ca. +2.5% a - 1 ) during the period. This was attributed to the increasing gross-primary productivity GPP and occurred without detectable change in ET. The start of annual carbon uptake period was advanced by 0.7 d a - 1 , and increase in GPP and NEP outside the main growing season contributed ca. one-third and one-fourth of the annual trend, respectively. Meteorological factors were responsible for the IAV of fluxes but did not explain the long-term trends. The growing season GPP trend was strongest in ample light during the peak growing season. Using a multi-layer ecosystem model, we showed that direct CO 2 fertilization effect diminishes when moving from leaf to ecosystem, and only 30-40% of the observed ecosystem GPP increase could be attributed to CO 2 . The increasing trend in leaf-area index (LAI), stimulated by forest thinning in 2002, was the main driver of the enhanced GPP and NEP of the mid-rotation managed forest. It also compensated for the decrease of mean leaf stomatal conductance with increasing CO 2 and LAI, explaining the apparent proportionality between observed GPP and CO 2 trends. The results emphasize that attributing trends to their physical and physiological drivers is challenged by strong IAV, and uncertainty of LAI and species composition changes due to the dynamic flux footprint. The results enlighten the underlying mechanisms responsible for the increasing terrestrial carbon uptake in the boreal zone., (© 2022 The Authors. Global Change Biology published by John Wiley & Sons Ltd.)

This study explores hydrological response of urban catchment in Southern Finland to climate change and urbanization. Process-based urban hydrological modelling and statistical analysis are applied to various urbanization and climate scenarios. Future changes in precipitation and temperature under Representative Concentration Pathways 4.5 and 8.5 (RCP4.5 and RCP8.5, respectively) clearly influence urban streamflow all year-round. We found snowpack shrinks during 2061 to 2090, snowmelt becomes earlier and the amount of melted runoff is reduced under both climate scenarios. The most significant runoff increase occurs in winter with the growth rates of 79% and 127%, respectively. It is also found that the dominant portion of urban streamflow shifts from summer to autumn in the future under both RCP4.5 and RCP8.5. Results indicate that urbanization has direct impact on hydrological response due to the change of imperviousness, but climate change will have more significant impact on seasonal distribution of urban streamflow. Additionally, urbanization impacts shrink monthly streamflow differences along with climate change., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Maintaining and improving surface water quality requires knowledge of nutrient and sediment loads due to past and future land-use practices, but historical data on land cover and its changes are often lacking. In this study, we tested whether land-use-specific export coefficients can be used together with satellite images (Landsat) and/or regional land-use statistics to estimate riverine nutrient loads and concentrations of total nitrogen (TN), total phosphorus (TP), and suspended solids (SS). The study area, Simojoki (3160 km 2 ) in northern Finland, has been intensively drained for peatland forestry since the 1960s. We used different approaches at multiple sub-catchment scales to simulate TN, TP, and SS export in the Simojoki catchment. The uncertainty in estimates based on specific export coefficients was quantified based on historical land-use changes (derived from Landsat data), and an uncertainty boundary was established for each land-use. The uncertainty boundary captured at least 60% of measured values of TN, TP, and SS loads or concentrations. However, the uncertainty in estimates compared with measured values ranged from 7% to 20% for TN, 0% to 18% for TP, and 13% to 43% for SS for different catchments. Some discrepancy between predicted and measured loads and concentrations was expected, as the method did not account for inter-annual variability in hydrological conditions or river processes. However, combining historical land-use change estimates with simple export coefficients can be a practical approach for evaluating the influence on water quality of historical land-use changes such as peatland drainage for forest establishment., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

The FLUXNET2015 dataset provides ecosystem-scale data on CO 2 , water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

Earth observing systems are now routinely used to infer leaf area index (LAI) given its significance in spatial aggregation of land surface fluxes. Whether LAI is an appropriate scaling parameter for daytime growing season energy budget, surface conductance (G s ), water- and light-use efficiency and surface-atmosphere coupling of European boreal coniferous forests was explored using eddy-covariance (EC) energy and CO 2 fluxes. The observed scaling relations were then explained using a biophysical multilayer soil-vegetation-atmosphere transfer model as well as by a bulk G s representation. The LAI variations significantly alter radiation regime, within-canopy microclimate, sink/source distributions of CO 2 , H 2 O and heat, and forest floor fluxes. The contribution of forest floor to ecosystem-scale energy exchange is shown to decrease asymptotically with increased LAI, as expected. Compared with other energy budget components, dry-canopy evapotranspiration (ET) was reasonably 'conservative' over the studied LAI range 0.5-7.0 m 2 m -2 . Both ET and G s experienced a minimum in the LAI range 1-2 m 2 m -2 caused by opposing nonproportional response of stomatally controlled transpiration and 'free' forest floor evaporation to changes in canopy density. The young forests had strongest coupling with the atmosphere while stomatal control of energy partitioning was strongest in relatively sparse (LAI ~2 m 2 m -2 ) pine stands growing on mineral soils. The data analysis and model results suggest that LAI may be an effective scaling parameter for net radiation and its partitioning but only in sparse stands (LAI <3 m 2 m -2 ). This finding emphasizes the significance of stand-replacing disturbances on the controls of surface energy exchange. In denser forests, any LAI dependency varies with physiological traits such as light-saturated water-use efficiency. The results suggest that incorporating species traits and site conditions are necessary when LAI is used in upscaling energy exchanges of boreal coniferous forests., (© 2016 John Wiley & Sons Ltd.)

The majority of C, N and P in boreal lakes are in organic form. Organically bound nutrients are released through biodegradation or photodegradation which affects the water quality, eutrophication and greenhouse gas emissions of lakes. We tested whether open land-use data combined with land-use-specific export coefficients can be used to predict total organic carbon (TOC), dissolved organic nitrogen (DON) and dissolved organic phosphorus (DOP) loading and lake water concentrations. Using data from 12 lake catchments in eastern Finland, we found that land use and management of the catchment explained a substantial proportion of the variations in TOC (r 2  = 0.78), DON (r 2  = 0.55) and DOP (r 2  = 0.80) concentrations between lakes. The computation does not account for in-lake processes, which are reflected as mismatch between the predicted and observed concentrations. However, this simple practical approach is useful in ranking lakes according to their water quality. The results indicated that natural sources dominate TOC, DON and DOP exports; the background leachings accounted for 57-99 %, 48-96 % and 55-99 % of TOC, DON and DOP export, respectively. The proposed method has promise as a practical decision support tool for assessing the impacts of land use on water quality. The results showed that possibilities to control TOC, DON and DOP loading to surface waters are limited to catchments where the peatland proportion is low and anthropogenic sources significant.