Jukuri, open repository of the Natural Resources Institute Finland (Luke) All material supplied via Jukuri is protected by copyright and other intellectual property rights. Duplication or sale, in electronic or print form, of any part of the repository collections is prohibited. Making electronic or print copies of the material is permitted only for your own personal use or for educational purposes. For other purposes, this article may be used in accordance with the publisher’s terms. There may be differences between this version and the publisher’s version. You are advised to cite the publisher’s version. This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): Lauri Heiskanen, Juha-Pekka Tuovinen, Henriikka Vekuri, Aleksi Räsänen, Tarmo Virtanen, Sari Juutinen, Annalea Lohila, Juha Mikola & Mika Aurela Title: Meteorological responses of carbon dioxide and methane fluxes in the terrestrial and aquatic ecosystems of a subarctic landscape Year: 2023 Version: Published version Copyright: The Author(s) 2023 Rights: CC BY 4.0 Rights url: http://creativecommons.org/licenses/by/4.0/ Please cite the original version: Heiskanen, L., Tuovinen, J.-P., Vekuri, H., Räsänen, A., Virtanen, T., Juutinen, S., Lohila, A., Mikola, J., and Aurela, M.: Meteorological responses of carbon dioxide and methane fluxes in the terrestrial and aquatic ecosystems of a subarctic landscape, Biogeosciences, 20, 545–572, https://doi.org/10.5194/bg-20-545-2023, 2023. Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 © Author(s) 2023. This work is distributed under the Creative Commons Attribution 4.0 License. R esearch article Meteorological responses of carbon dioxide and methane fluxes in the terrestrial and aquatic ecosystems of a subarctic landscape Lauri Heiskanen1, Juha-Pekka Tuovinen1, Henriikka Vekuri1, Aleksi Räsänen2,3, Tarmo Virtanen2, Sari Juutinen2,4, Annalea Lohila1,5, Juha Mikola2,6, and Mika Aurela1 1Climate System Research Unit, Finnish Meteorological Institute, Helsinki, Finland 2Ecosystems and Environment Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland 3Bioeconomy and Environment Unit, Natural Resources Institute Finland (LUKE), Oulu, Finland 4Department of Geographical and Historical Studies, University of Eastern Finland, Joensuu, Finland 5Institute for Atmospheric and Earth System Research (INAR)/Physics, Faculty of Science, University of Helsinki, Helsinki, Finland 6Bioeconomy and Environment Unit, Natural Resources Institute Finland (LUKE), Helsinki, Finland Correspondence: Lauri Heiskanen (lauri.heiskanen@helsinki.fi) Received: 11 March 2022 – Discussion started: 6 April 2022 Revised: 14 December 2022 – Accepted: 11 January 2023 – Published: 3 February 2023 Abstract. The subarctic landscape consists of a mosaic of forest, peatland, and aquatic ecosystems and their ecotones. The carbon (C) exchange between ecosystems and the at- mosphere through carbon dioxide (CO2) and methane (CH4) fluxes varies spatially and temporally among these ecosys- tems. Our study area in Kaamanen in northern Finland cov- ered 7 km2 of boreal subarctic landscape with upland forest, open peatland, pine bogs, and lakes. We measured the CO2 and CH4 fluxes with eddy covariance and chambers between June 2017 and June 2019 and studied the C flux responses to varying meteorological conditions. The landscape area was an annual CO2 sink of −45± 22 and −33± 23 gCm−2 and a CH4 source of 3.0± 0.2 and 2.7± 0.2 gCm−2 dur- ing the first and second study years, respectively. The pine forest had the largest contribution to the landscape-level CO2 sink, −126± 21 and −101± 19 gCm−2, and the fen to the CH4 emissions, 7.8± 0.2 and 6.3± 0.3 gCm−2, dur- ing the first and second study years, respectively. The lakes within the area acted as CO2 and CH4 sources to the at- mosphere throughout the measurement period, and a lake located downstream from the fen with organic sediment showed 4-fold fluxes compared to a mineral sediment lake. The annual C balances were affected most by the rainy peak growing season in 2017, the warm summer in 2018, and a heatwave and drought event in July 2018. The rainy pe- riod increased ecosystem respiration (ER) in the pine for- est due to continuously high soil moisture content, and ER was on a level similar to the following, notably warmer, sum- mer. A corresponding ER response to abundant precipitation was not observed for the fen ecosystem, which is adapted to high water table levels, and thus a higher ER sum was observed during the warm summer 2018. During the heat- wave and drought period, similar responses were observed for all terrestrial ecosystems, with decreased gross primary productivity and net CO2 uptake, caused by the unfavourable growing conditions and plant stress due to the soil moisture and vapour pressure deficits. Additionally, the CH4 emis- sions from the fen decreased during and after the drought. However, the timing and duration of drought effects varied between the fen and forest ecosystems, as C fluxes were af- fected sooner and had a shorter post-drought recovery time in the fen than forest. The differing CO2 flux response to weather variations showed that terrestrial ecosystems can have a contrasting impact on the landscape-level C balance in a changing climate, even if they function similarly most of the time. Published by Copernicus Publications on behalf of the European Geosciences Union. 546 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 1 Introduction A typical boreal subarctic landscape consists of a mosaic of land cover types including forests, open and forested peat- lands, and waterbodies, with each ecosystem acting as a source or sink of atmospheric carbon (C), depending on its characteristics and weather conditions. To grasp a full pic- ture of the C balance at the landscape scale, the carbon diox- ide (CO2) and methane (CH4) exchange between the main ecosystems and the atmosphere needs to be assessed. This is also because part of the C fixed in terrestrial systems can be emitted back to the atmosphere via aquatic ecosystems or transported out from the catchment. The net atmosphere– ecosystem CO2 exchange is a result of gross primary produc- tion (GPP) of plants and ecosystem respiration (ER), which is composed of autotrophic and heterotrophic respiration (Chapin et al., 2011). The CH4 emissions from ecosystems are mainly produced by methanogens in the waterlogged, anoxic zone in soils, while the CH4 uptake is caused by ox- idation by methanotrophs in aerobic soil conditions (Le Mer and Roger, 2001). As a result of the ongoing climate change, the subarctic regions warm rapidly, 2 to 3 times as fast as the rest of the world (Masson-Delmotte et al., 2018). The higher tempera- tures affect the C cycle of subarctic ecosystems by length- ening the growing season and shortening the snow and ice cover periods. Higher temperatures and a longer growing season enable greater CO2 uptake through photosynthesis (Aurela et al., 2004; Silfver et al., 2020), but they also al- low greater CO2 emissions due to the increased heterotrophic and autotrophic respiration (Kätterer et al., 1998). However, the variation in ecosystem C dynamics is also highly de- pendent on water balance, as too low a moisture content in soils decreases the GPP and decomposition of organic matter (Jones, 2013; Meyer et al., 2018), and methanogenesis re- quires anoxic soil conditions (Chapin et al., 2011). In aquatic ecosystems, the longer ice-free periods, with increased C in- put from surrounding terrestrial ecosystems, lead to higher annual C emissions to the atmosphere (Cole et al., 2007; Guo et al., 2020). In a longer than 1-year time span, the effects of climate change are observed as vegetation composition shifts, such as Arctic shrubification that enhances CO2 up- take and alters ecosystem respiration and nutrient availability (Mekonnen et al., 2021). Concurrently, the potential of severe weather events, such as heatwaves and droughts, to affect the C cycle increases within the region (Masson-Delmotte et al., 2018). Boreal upland forests have been shown to act as either annual net sinks or sources of atmospheric CO2, depending on species composition, climatic conditions, and the occur- rence of extreme weather events, fires, and herbivores and pathogens (Bond-Lamberty et al., 2007; Kljun et al., 2007; Lindroth et al., 2008; Aurela et al., 2015; Hadden and Grelle, 2017), and these ecosystems usually act as small CH4 sinks (Dinsmore et al., 2017). Undrained boreal peatlands have been, in the long term, net CO2 sinks via the accumulation of organic C, but, similar to forests, their annual CO2 balance is sensitive to environmental conditions (Bubier et al., 1998; Aurela et al., 2004; Lindroth et al., 2007). Additionally, these peatlands release CH4 into the atmosphere (Chapin et al., 2011). Peatlands have more C stored in the soil but less C in living plant matter compared to upland forests, and thus, the same disturbances on the ecosystems can lead to differing effects on the C balance. Globally, inland waters act as net CO2 and CH4 sources to the atmosphere (Tranvik et al., 2009). Boreal lakes vary greatly in their area and depth, input of organic and inor- ganic C from the terrestrial ecosystems, ice-free period dura- tion, and the amount of C in sediments and aquatic vegeta- tion, which all affect the magnitude of the C exchange with the atmosphere (Bastviken et al., 2004; Wik et al., 2018; Den- feld et al., 2020). Additionally, waterbodies provide paths for the lateral flow of water, C, and other matter between ecosys- tems (Cole et al., 2007). The landscape-level C balance comprised of diverse ecosystems is evidently dependent on the ecosystem compo- sition, and the regional effect of an environmental change or disturbance depends on the integrated response of individual ecosystems. As an example of such a disturbance, a heatwave and drought event encompassed northwestern Europe during the summer of 2018 (Lehtonen and Pirinen, 2019a, b). The drought conditions had different effects on the CO2 exchange of boreal forests and peatlands in the area (Rinne et al., 2020; Lindroth et al., 2020; Matkala et al., 2021). There was a large variation in how different forests reacted to the heatwave and drought, with some showing a substantial decrease and some no change or even an increase in CO2 uptake (Lindroth et al., 2020; Matkala et al., 2021). Most of the studied peatlands turned momentarily into CO2 sources during the drought, which lowered their annual CO2 uptake significantly and, simultaneously, the CH4 emissions decreased (Rinne et al., 2020). Thus, the landscape-level C balance depends on mul- tiple responses, whose net effects are not easily predictable. In this study, we assess the ecosystem–atmosphere ex- change of CO2 and CH4 within a subarctic landscape based on eddy covariance and chamber measurements. The stud- ied area in Kaamanen in northern Finland includes an upland pine forest, pine bogs, a mesotrophic flark fen, and shallow lakes of glacial origin. We focus on the temporal variation in C exchange during 2 full years (June 2017–June 2019). The meteorological conditions during the first year were similar to the long-term average, with the exception of the higher- than-average precipitation sum in summer. The second year included the above-mentioned heatwave and drought period. Motivated by the plant community-level study of Heiskanen et al. (2021), who found that the diversity of plant commu- nities constrained C loss from the Kaamanen fen during the 2018 drought, we study if a similar pattern can be observed at the ecosystem level. Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 547 Figure 1. The study area in Kaamanen with land cover type classi- fication and flux measurement site locations. Here, we address questions on how sensitive the C fluxes of different ecosystems are to changes in environmental con- ditions and how this is reflected on the landscape-level C ex- change. The more specific scientific questions addressed are as follows: (1) what are the contributions of each ecosys- tem to the landscape-level CO2 and CH4 fluxes? (2) Do the CO2 and CH4 fluxes of different ecosystems show similar responses to varying meteorological conditions? 2 Materials and methods 2.1 Study site The studied landscape is located in Kaamanen in northern Finland (69◦8′ N, 27◦16′ E; 155 ma.s.l. – above sea level ) within the subarctic climate zone and the northern boreal vegetation zone. The annual mean temperature in the re- gion was −0.4 ◦C and the mean annual precipitation sum 472 mm in 1981–2010 (Pirinen et al., 2012). Even though the study area is located within the sporadic permafrost zone, no permafrost has been found there anymore in recent decades (Fronzek et al., 2010). The growing season is short and lasts for 150–180 d (Aurela et al., 2002). The study area cov- ers 7 km2 around the flux measurement sites (Fig. 1) and consists of five main ecosystem types, namely an upland pine forest, a patterned mesotrophic flark fen, a treed pine bog, a sparsely treed pine bog, and lakes and a connecting stream (Fig. 2). The ecosystems within the study area are homoge- neous in their species composition and provide a representa- tive sample of the landscape. The areal coverage of different ecosystem types in the landscape was estimated with a land cover classification Table 1. Relative areal coverage of the land cover types within the study area. Land cover type Areal coverage (%) Upland pine forest 29.2 Treed pine bog 9.1 Sparsely treed pine bog 13.5 Fen 26.1 Waterbody, mineral sediment 12.5 Waterbody, organic sediment 7.8 Non-vegetated 1.8 utilising remote sensing data and field observations, fol- lowing the methodology described by Räsänen and Virta- nen (2019). An aerial orthophoto was segmented, and 55 fea- tures for each segment were calculated from multi-source re- mote sensing data, including orthophoto, aerial laser scan- ning, and satellite imagery, and a supervised random forest classification (Breiman, 2001) was carried out using field data for model training and validation (Sect. S1 in the Sup- plement). The pine forest, pine bog, fen, and lake ecosystems each encompass roughly an equal area (Table 1). The ecosys- tems in the area are pristine with the exception of some forest logging, thinning, and selective removal of birches in part of the pine-dominated forests. The mesotrophic patterned flark fen is an open peatland ecosystem characterised by a mosaic of string and flark mi- croforms. The string formations are 0.5–1 m high and can re- main frozen until the late summer (Aurela et al., 2001). The fen vegetation has also a clear patterning and can be parti- tioned to five distinct plant community types (PCTs) that dif- fer in their vegetation composition, water table level, and car- bon fluxes such as flark, Trichophorum tussock, string mar- gin, string top, and tall sedge fen (Maanavilja et al., 2011; Räsänen et al., 2019; Heiskanen et al., 2021). In the flarks, the plant communities are dominated by sedges, including Trichophorum tussocks, and brown mosses. The top of the strings act as ombrotrophic bog-like surfaces within the fen and are covered mainly by dwarf shrubs, herbs, mosses, and lichens, while string margins are populated by dwarf birch (Betula nana) and other dwarf shrubs, sedges, and Sphag- num mosses. The fifth PCT, a tall sedge fen that is covered by tall sedges, deciduous shrubs, and forbs, can be found in the riparian areas of lakes and small streams. The species com- position of these PCTs was described in detail by Maanav- ilja et al. (2011) and Räsänen et al. (2019) and in the pale- orecords by Piilo et al. (2020). The bedrock under the peat- land slopes towards the south. The average peat thickness increases from 1 m in the northern part of the fen towards the south, where it is up to 4 m (Piilo et al., 2020). Most of the aboveground biomass of the peatland resides in shrubs and mosses, and forbs and graminoids contribute increasingly to https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 548 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape Figure 2. The main ecosystem types within the study area are a (a) Scots pine forest, (b) treed pine bog, (c) sparsely treed pine bog, (d) fen, and (e) lake. Photo credits: (a) Lauri Heiskanen and (b–e) Tarmo Virtanen. the total leaf area in the mid- and late growing season (Ta- ble A1 in the Appendix). The dominant tree species in the upland mineral soil for- est is Scots pine (Pinus sylvestris), with a few downy birch (Betula pubescens) trees also present (Table A2). The forest around the flux measurement site was logged about 50 years ago, but most of the pine forests within the study area are pine-dominated old-growth forest with an uneven age distri- bution. The average heights of all trees and the main canopy were approximately 8 and 11 m, respectively (Table A3). The field layer is dominated by evergreen dwarf shrubs (e.g. Vac- cinium vitis-idaea and Calluna vulgaris), and the ground layer is covered by mosses and lichens. The soil type in the forest is sandy podzol, which has higher bulk density, lower nitrogen (N) content, and higher C : N ratio than the soil in the peatlands (Table A4; Sect. S2). Mosses and lichens contribute about 10 % to the total aboveground biomass (Ta- bles A1 and A5). The tree height and density decrease from the upland pine forest to the pine bog ecosystems. The average tree height in the pine bog was 5 m, and in the treed pine bog and sparsely treed pine bog, the aboveground tree biomass was about 50 % and 16 %, respectively, of the biomass in the pine forest (Ta- bles A3 and A5). The field layer vegetation of the pine bog ecosystems is similar to that of the string tops (Table A1), with evergreen (Ledum palustre) and deciduous (Vaccinium uliginosum and Betula nana) shrubs, herbs (Rubus chamae- morus), and graminoids (mostly Carex spp.). The pine bog ground layer is formed by mosses and a few lichens. The soil carbon content is higher in the sparsely treed bog patches within the fen than in bog vegetation with thinner organic layer at the forest edge. Soil N content and C : N ratio are similar to those found in pine forest and string tops (Ta- ble A4). The waterbodies within the study area vary in size from small streams and ponds to lakes larger than 25 ha (Lake Ulkujärvi in the northwest and Lake Jänkäjärvi in the south; Fig. 1). The lakes in the area are shallow, with the depth rang- ing from less than 1 m to a few metres, and have a sandy or organic sediment bottom. Most waterbodies in the study area belong to the same catchment, with water flowing from north to south. The aquatic vegetation is sparse, consisting of macrophytes, mainly horsetail (Equisetum fluviatile), sedges (Carex spp.), Menyanthes trifoliata, and benthic algae, near the shore. A few-metres-wide stream flows through the peat- land and connects Lake Ulkujärvi to Lake Jänkäjärvi. Addi- tionally, water flows on the surface across the fen and through the peat layers and eskers. Each year, during the spring thaw, the fen and part of the pine bogs become flooded. 2.2 Flux measurement methods 2.2.1 Pine forest CO2 flux measurements The CO2 flux measurements in the Scots pine forest were set up at the beginning of the growing season in 2017, and the data acquisition was started on 8 June 2017. The eddy co- variance (EC) measurements were conducted on a 14 m tall tower that was located on the southern edge of the for- est so that, in the wind sector 250–65◦, the forest coverage was higher than 80 %. The EC system consisted of a three- axis sonic anemometer (USA-1; METEK Meteorologische Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 549 Messtechnik GmbH, Germany) and a closed-path infrared gas analyser for CO2 and H2O mixing ratios (LI-7200; LI- COR Biosciences, USA). The inlet tube for the gas analyser was ∼ 18 m long and had an inner tube diameter of 3.1 mm. The flow rate was 5–6 Lmin−1. The sampling frequency for the EC flux data was 10 Hz. Standard methods were used for calculating the half-hourly turbulent fluxes (Aubinet et al., 2012), with block averaging and a double rotation of the coordinate system (McMillen, 1988). The high-frequency signal attenuation flux losses were taken into account using an experimental transfer func- tion with a half-power frequency of 1.9 Hz (Laurila et al., 2005). The half-hourly averaged data were accepted based on the following screening criteria: relative stationarity< 100 % (Foken and Wichura, 1996), amount of recorded data per 30 min> 17 400, number of signal spikes per 30 min< 360, and mean CO2 mixing ratio within 340–550 ppm (parts per million). A friction velocity (u∗) limit of 0.24 ms−1 was used for screening the periods of insufficient turbulence. Data from the wind direction sector 250–65◦ were used to calcu- late the CO2 balances of the forest. As the net ecosystem exchange of CO2 in forests (NEEforest) can be significantly affected by the flux due to storage change below the measurement height, this was es- timated according to Montagnani et al. (2018) and added to the measured eddy flux (FCO2 ) as follows: NEEforest = FCO2 + ρd 1CO2 1t h, (1) where ρd is the mean dry air density, 1CO2 1t is the change in CO2 dry molar fraction at the EC measurement height dur- ing the 30 min averaging period, and h is the measurement height. The CO2 concentration profile below the measure- ment height was assumed to be constant. 2.2.2 Fen CO2 and CH4 flux measurements The CO2 and CH4 fluxes of the Kaamanen fen were mea- sured with both the EC and flux chamber methods. The EC measurements were used for the fen ecosystem flux analysis here, and the PCT-specific chamber-based data were used for estimating the CO2 fluxes of the pine bog ecosystems. Even though the fen is comprised of a mosaic of different PCTs, the contribution of each type is similar in all wind directions within a 100–150 m radius, and thus, the footprint variation is not expected to bias the ecosystem balances derived from EC measurements. The measurements were conducted on a 5 m tall tower with a three-axis sonic anemometer (USA-1; METEK Me- teorologische Messtechnik GmbH, Germany), a closed-path infrared gas analyser for CO2 and H2O mixing ratios (LI- 7000; LI-COR Biosciences, USA), and a laser-based gas analyser for the CH4 mixing ratio (RMT-200; Los Gatos Re- search, USA). The heated inlet tubes (inner diameter 3.1 Table 2. Area coverage of the plant community types of the peat- land inside a 200 m radius around the eddy covariance tower. Plant community type Area coverage (%) Flark 37 Trichophorum tussock 10 Tall sedge fen 17 String margin 14 String top 16 Pine bog 6 Table 3. Coefficients for scaling the chamber-based ER and GPP fluxes. Time period ER ratio GPP ratio June 2017–October 2017 1.170 1.514 November 2017–April 2018 1.012 1.337 May 2018–June 2019 0.854 1.160 and 8 mm, for LI-7000 and RMT-200, respectively) were 6 m long and had a flow rate of 6 and 15 Lmin−1 for the LI- 7000 and RMT-200, respectively. The sampling frequency was 10 Hz for both analysers. The same flux calculation and data processing methods were used as with the pine forest EC data, except for the discarded wind direction sector (260– 315◦), the u∗ limit (0.1 ms−1), and the relative stationarity limit (30 %; Foken and Wichura, 1996); in addition, a mean CH4 mixing ratio of 1.8–2.8 ppm was required. The manual chamber measurements were conducted bi- weekly in 12 June–11 October 2017 and 31 May–4 Septem- ber 2018 on four main plant communities that grow on the microtopography gradient (flark, Trichophorum tussock, string margin, and string top). A total of 17 chamber plots were used, with four or five replicates on each PCT. For de- termining the ER flux and the light response of CO2 flux, a transparent chamber was used with one to three shading ele- ments over the chamber. The chamber was closed for 2 min during each measurement, and the CO2 and CH4 fluxes were calculated from the mixing ratio change measured inside the chamber. The EC and chamber measurements at the fen and the PCT-specific fluxes are presented in more detail by Heiskanen et al. (2021). For modelling the string top CO2 fluxes for the time span from June 2017 to June 2019, to be used for pine bog (Sect. 2.4), the chamber-based fen ecosystem fluxes were scaled to match the EC-based fen fluxes. This was done sim- ilarly to Piilo et al. (2020), by first upscaling the chamber fluxes of the main PCTs of the fen to the ecosystem scale according to their relative areas within 200 m from the EC tower (Table 2). The PCT-specific fluxes were presented in more detail by Heiskanen et al. (2021). The mean ratios between the EC-based and upscaled chamber-based ER and GPP fluxes were used for scaling the https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 550 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape chamber flux data (Table 3). For ER, the mean ratio was cal- culated from all data, while the mean GPP ratio was calcu- lated from the daytime data (from 06:00 to 18:00 LST, local standard time). 2.2.3 Eddy covariance data gap-filling and partitioning The EC flux time series of both the pine forest and fen ecosystems had substantial gaps due to equipment failures, wind sector exclusions, and quality control filtering applied during the post-processing of the data. In total, 86 % and 91 % of the forest CO2 flux data were gap-filled in the time series of 11 June 2017–10 June 2018 and 11 June 2018– 10 June 2019, respectively. During the growing seasons, i.e. in 11 June–31 October 2017 and 11 June–31 October 2018, the gap-filling percentage was 77 % and 84 %, respec- tively. In total, 64 % and 63 % of the fen CO2 flux data were gap-filled in the time series of the first and second years, respectively. The contribution of each data filtering step is shown in Table S1 in the Supplement. For gap-filling the pine forest CO2 fluxes during the win- ter, when the data coverage was lowest but the fluxes rela- tively small and stable, we averaged the measurement data pooled into two soil temperature categories, i.e. over and un- der −2 ◦C at the 10 cm depth. The mean fluxes were 0.0173 and 0.0103 mgCO2 m−2 s−1 in the warmer and colder cate- gory, respectively. The winter period was determined based on the timing of frost at 10 cm depth, which occurred in 31 September 2017–24 April 2018 and 30 September 2018– 20 April 2019. The growing season pine forest CO2 flux data and the entire fen CO2 flux data time series were gap-filled using a machine learning method called extreme gradient boost- ing (XGBoost; Chen and Guestrin, 2016). XGBoost is a decision-tree-based ensemble method in which trees are built in a sequential manner so that each tree corrects for the er- rors in the previous trees. XGBoost has been shown to per- form well, even with long gaps in the EC data (Zhu et al., 2022; Irvin et al., 2021). The environmental variables that were used to predict NEE were photosynthetic photon flux density (PPFD), air temperature, relative humidity, vapour pressure deficit (VPD), soil temperature at−10 cm, soil tem- perature at −5 cm, and soil moisture at −10 cm. First, we optimised the hyperparameters of the model using a grid search to control for overfitting and to select the best model architecture. The hyperparameters that were optimised control for the maximum depth of a tree (“max_depth”), the subsample ratio of columns when constructing each tree (“colsample_bytree”), learn- ing rate (“learning_rate”), and the minimum number of samples required to create a new node in a tree (“min_child_weight”). The selected hyperparameters were 0.8 for colsample_bytree, 0.05 for learning_rate, 20 for max_depth, and 9 for min_child_weight. The squared er- ror was used as the loss function. Model performance for the more sparse pine forest data was evaluated using 10- fold cross-validation, resulting in the coefficient of deter- mination (R2) of 0.88± 0.02 and the mean squared error of 0.003± 0.0006 mgCO2 m−2 s−1 (Fig. S3 in the Supple- ment). Further details of the model performance and justi- fication for gap-filling the EC flux time series with a high proportion of missing data are presented in Sect. S3. The ER and GPP fluxes were partitioned based on environ- mental response functions that were fitted to the gap-filled NEE data. A rectangular hyperbola was used to model the dependency of GPP on PPFD (e.g. Whiting, 1994), while the exponential model of Lloyd and Taylor (1994) was used for ER. For the pine forest flux parameterisation, air tempera- ture was used during the growing season, while soil temper- ature (Ts) at the 10 cm depth was used for the other periods; for the fen, Ts at the 10 cm depth was used for the whole year. For fitting the parameters of the ER function, night- time data were used (PPFD< 30 µmolm−2 s−1). The activa- tion energy parameter (E0) of the respiration function was fitted first with a 91 d moving window, after which the base respiration rate at 10 ◦C (R10) was fitted to the data within a moving window of 7 d. The fitting of the GPP function was performed in a moving window of 3 d, using the day- time flux data (PPFD> 30 µmolm−2 s−1). For gap-filling the CH4 fluxes measured at the fen, a simple moving average interpolation of the half-hour fluxes was used. A moving av- erage window of± 1, 2, 4, 8, 16, or 32 d was used, depending on the length of the gaps in data. In total, 69 % and 70 % of the CH4 flux data were gap-filled in the time series of the first and second years, respectively (Table S1 in the Supplement). 2.2.4 Lake CO2 and CH4 flux measurements The diffusive fluxes of CO2 and CH4 were measured on the intermediate-sized (1.5 ha) Lake Jänkälampi, with min- eral sediment (MS), and on the larger Lake Jänkäjärvi, with organic sediment (OS; Fig. 1). On Jänkäjärvi, CH4 ebulli- tion was also measured. The lakes are hereafter referred to as the MS and OS lakes, respectively. The MS lake fluxes were measured on 5 d during June–October in 2017 and an- other 5 d during June–September 2018. The OS lake fluxes were measured bi-weekly during June–August 2017. The flux measurements were conducted with floating flux chambers, while the CH4 ebullition was determined from both chamber and floating bubble collector data. The diffu- sive fluxes of the MS lake were measured with an opaque aluminium chamber (60 cm× 60 cm× 30 cm) at 20 m from the north shore. The chamber air was mixed with a battery- driven fan. The closure time was 7 min, after which the chamber was ventilated for 3 min. The changes in CO2, CH4, and H2O mixing ratios inside the chamber were mea- sured using a closed-path laser-based gas analyser (Picarro G2401; Picarro, Inc., USA) connected to the chamber via a 50 m long inlet tube (Teflon; inside diameter 3.1 mm). On the OS lake, floating chambers with a volume of 8 L and an area Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 551 of 0.05 m2 were applied, using a 30–60 min closure. Four samples (30 mL) of the gas space were drawn using polyethy- lene syringes, and the samples were stored in 12 mL glass vials flushed with sample air. Samples were analysed within a month, using a gas chromatograph equipped with electron capture, thermal conductivity, and flame ionisation detectors (Agilent 7890B, with a Gilson GX-271 autosampler). We tested CO2 fluxes against those determined using an online CO2 sensor (K33 ELG CO2 module; Senseair) in three of the chambers and obtained similar flux values. In all, 5 to 10 chambers were deployed at a time. The chambers were lined with a rope, anchored at both ends, so that they cap- tured the open-water section of the lake from near-shore to- wards the lake centre. The chambers captured some ebulli- tion events, and these were included as bubble flux. In addi- tion, funnel bubble collectors with an area of 0.03 m2 were floated. Auxiliary flux measurements from three lakes out- side the study area were conducted with the same method as for the OS lake. These measurements were utilised to assess the representativeness of the flux data collected on the MS and OS lakes. The diffusive fluxes of CO2 and CH4 were calculated from the change in the CO2 and CH4 mixing ratios inside the closed chamber. The rate of change was calculated with lin- ear regression based on ordinary least squares, similar to the fen ecosystem chamber measurements by Heiskanen et al. (2021). The measurement data of diffusive fluxes were screened to reject cases with a non-linear concentration change and disturbances due to chamber leakage and ebullition events. The number of accepted/total data of the MS lake fluxes was 54/56 in 2017 and 62/75 in 2018. The number of ac- cepted/total OS lake data was 42/50 for diffusive CO2 fluxes and 49/50 for diffusive CH4 fluxes. There were ebullition events in 7 % of the individual flux measurements, and this ebullition frequency was used to estimate the total ebullition- induced flux over the ice-free period. 2.3 Abiotic and biotic environmental measurements In addition to the C flux measurements, meteorological and environmental variables were measured continuously at and close to the flux measurement locations in the pine forest, fen, and lake ecosystems. Air temperature, precipitation sum, and snow depth were measured at the weather station of the Finnish Meteorological Institute, located between the pine forest and fen EC measurement sites. Air temperature and humidity (Vaisala; HMP230), global and reflected radiation (Kipp & Zonen; CM7), and downward and upward PPFD (Kipp & Zonen; PQS1) were measured at the forest and fen sites at heights of 14 and 3 m, respectively. The water vapour pressure deficit was calculated from air temperature and rel- ative humidity according to Jones (2013). Forest soil temper- ature was measured at 5 and 10 cm depths (Onset; HOBO) and soil moisture at 10 cm depth (Onset; HOBO). Ts profiles at the fen were measured from a string (at 10, 30, 50, 75, and 105 cm depth) and flark (at 10, 30, and 50 cm depth; IKES Pt100 sensors). The water table level (WTL) was measured manually from the flux chamber positions at the fen and OS lake, while the WTL of the MS lake was measured with a pressure sensor (Onset; HOBO U20) at the floating flux chamber position. Lake water temperature (Onset; HOBO Pendant) was measured at the same position at 10 cm from the lake bottom. The ice-free period of the flux measurement lakes was de- termined from air temperature data and repeat digital pho- tographs at the fen site (Linkosalmi et al., 2022). In October 2018, the lake freezing was recorded with a temperature and pressure logger (Onset; HOBO). The freezing occurred af- ter the air temperature remained continuously below 0 ◦C for 3 d, with no subsequent exceedances. The timing of thaw was determined from the snowmelt and ice thaw at the fen. We defined meteorological drought based on the atmo- spheric VPD as the period during which the daily maximum VPD (VPDmax) exceeded 20 hPa (Lindroth et al., 2007; Au- rela et al., 2007). The Standardised Precipitation Evapotranspiration Index (SPEI), which takes into account both precipitation and po- tential evapotranspiration in determining drought conditions (Vicente-Serrano et al., 2010), was used as a climatological reference for the study period. Monthly SPEI data covering the years 1950–2018 for the 0.5◦× 0.5◦ grid cell including Kaamanen were extracted from the global SPEI database (SPEIbase v2.6; https://spei.csic.es/database.html, last ac- cess: 19 October 2021). 2.4 Estimation of pine bog fluxes In the subarctic aapa mires, pine bogs typically form nar- row zones bordering forests and peatlands. Conducting direct EC flux measurements on them is challenging, as the fetch is too limited for the EC method. Thus, the fluxes of the two pine bog ecosystems within the study area were not directly measured but, to enable regional upscaling, were modelled based on the fluxes measured in the forest and fen ecosys- tems. The use of forest and fen fluxes was considered appro- priate because, first, the tree species composition in pine for- est and pine bog is similar, and, second, the peat soil ground layer vegetation of pine bog is similar to that of the string top PCT at the fen (Table A1). For estimating the pine bog fluxes, the total forest flux was assumed to consist of the flux of trees and the flux of the other forest ecosystem elements, Fforest = Fforest, trees+Fforest, other, (2) and, similarly, for the pine bog ecosystem, Fpine bog = Fpine bog, trees+Fpine bog, other. (3) https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 https://spei.csic.es/database.html 552 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape The “trees” fluxes were assumed to be proportional to tree biomass, i.e. Fpine bog, trees = a×Fforest, trees, (4) where a is the tree biomass ratio between pine bog and pine forest (a= 0.52 and 0.16 for treed and sparsely treed pine bogs, respectively; Table A5; Sect. S4). Assuming further that the ground layer (“other”) fluxes in both the forest and pine bog equal the flux of the fen’s driest PCT, i.e. string top, Fforest, other = Fpine bog, other = Fstring top. (5) Equations (2)–(5) show that the GPP and ER fluxes for the pine bog can be estimated as follows: Fpine bog = a×Fforest+ (1− a)×Fstring top, (6) where Fforest and Fstring top are obtained from measurements. 2.5 Upscaling fluxes to landscape level The CO2 and CH4 fluxes of different ecosystems were up- scaled to the landscape-level by taking into account the areal contributions inside the 7 km2 study area (Tables 1 and A7) and summing the half-hourly area-weighted NEE, GPP, ER, and CH4 flux estimates of each ecosystem. The pine forest and pine bog daily average growing season CH4 flux esti- mates were taken from the literature (Bubier et al., 2005; Dinsmore et al., 2017), which were further calculated for 180 d long growing seasons at Kaamanen. The annual CO2 and CH4 balances for lakes were derived from the results of Kou et al. (2022), who modelled these fluxes for the Kaamanen catchment with the Arctic Lake Bio- chemistry Model (ALBM). The ALBM is a one-dimensional process-based model operating on a daily time step (Tan et al., 2015, 2017). For Kaamanen, it was calibrated with the flux measurements described in Sect. 2.2.4 and driven with local meteorological data (Kou et al., 2022). 2.6 Estimating flux uncertainty The uncertainty in the EC-based CO2 and CH4 fluxes of the pine forest and fen ecosystems were estimated by tak- ing into account the most significant error sources. First, the uncertainty related to low-turbulence screening was esti- mated by filtering the data with 100 bootstrapped u∗ thresh- olds and gap-filling each of the resulting time series. The re- lated uncertainty was defined as the standard deviation of the 100 gap-filled 30 min NEE values. Second, to estimate the random measurement uncertainty, we binned the measured data into 0.2 (pine forest) or 0.1 mgCO2 m−2 s−1 (fen) wide bins and calculated the standard deviation of model residu- als for each bin. The linear relationship between this stan- dard deviation and the magnitude of the flux was estimated separately for negative and positive fluxes and used to esti- mate the uncertainty in the measured data (Richardson et al., 2008). Third, the uncertainty due to gap-filling with XG- Boost was estimated following the procedures detailed by Irvin et al. (2021). This uncertainty component was estimated from an ensemble of 10 models that were fitted using boot- strapped data, and the uncertainty intervals were calibrated using Platt scaling as an additional post-processing step. To estimate the uncertainty in the annual balances, the uncertainty related to u∗ was determined as the standard deviation of the 100 bootstrapped balances. Also, the gap- filling uncertainties were determined using a bootstrapping approach. For winter, we resampled the wintertime data 100 times with replacement and calculated alternative bal- ances. For the growing season, we used the ensemble of 10 XGBoost models to calculate alternative balances. In both cases, the gap-filling uncertainty was determined as the stan- dard deviation of the balances. For the flux-chamber-based monthly CO2 and CH4 flux sums of the lake ecosystem, the measurement uncertainty was estimated as the standard deviation of the flux data from the repeated chamber closures. The uncertainty in the annual CO2, CH4, and C flux sums were accumulated as the root sum square of individual uncertainties. For the modelled pine bog fluxes, both for treed and sparsely treed bogs, the uncertainty was combined from the pine forest string top flux uncertainties proportionally to their contribution to the derived flux. The string top flux uncer- tainty was estimated by combining the uncertainty due to the estimated parameters of environmental response functions and the flux variation among the chamber plots (Heiskanen et al., 2021). The landscape-level flux uncertainty was calculated by combining the total uncertainties in each ecosystem, exclud- ing lakes due to differing data sets, proportionally to their areal coverages (Table 1). All uncertainty estimates reported for C balances refer to the 95 % confidence intervals. 2.7 Statistical analysis The estimates of ER, GPP, NEE, and CH4 fluxes of pine forest and fen ecosystems were compared between the two growing seasons using the Z test, for which we produced full flux time series with uncertainty estimates. Welch’s t test was used for comparing the monthly average lake fluxes between the 2 years, as the compared data comprised of a sample of measured fluxes. For identifying the main environmental drivers of the pine forest CO2 flux, a linear regression model was estimated to explain the 5 d averaged GPP and ER fluxes derived from the measured EC data. The FR data were normalised to 10 ◦C (denoted as FR10), while the FGPP data were normalised to a near-optimal PPFD level of 1200 µmolm−2 s−1 (denoted as FGPP1200; Laurila et al., 2001). Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 553 Figure 3. Monthly precipitation sum and mean air temperature at Kaamanen, and the corresponding 30-year averages (Pirinen et al., 2012) measured at the Ivalo weather station (68◦36′ N, 27◦25′ E; 59 km south of Kaamanen). The following fixed explanatory 5 d variables were tested in the models: precipitation sum, average VPDmax, soil mois- ture and temperature at a depth of 10 cm, air temperature, and daily maximum PPFD. Ts was included in the final analysis, while air temperature was excluded because these variables were strongly correlated (Pearson’s R= 0.96 in the FR data andR= 0.93 in the FGPP1200 data), and Ts showed higher ab- solute correlation with the response variables. The MODIS normalised difference vegetation index was also tested as an explanatory variable, but it was excluded from the final anal- ysis, as it strongly correlated with Ts (Pearson’s R= 0.87 in the FR data and R= 0.89 in the FGPP1200 data). For the final models, the explanatory variables were chosen with a step- wise procedure to both directions by minimising Akaike’s information criterion value. To evaluate the relative impact of each explanatory variable, the standardised regression coef- ficients were calculated. Data analyses were conducted in R (R Core Team, 2021) with the MASS (Venables and Ripley, 2002) and MuMIN (Barton, 2020) packages. 3 Results and discussion 3.1 Environmental conditions We studied CO2 and CH4 fluxes across different ecosys- tem types in Kaamanen from June 2017 to June 2019. The growing seasons (May–October) of 2017 and 2018 differed from each other in terms of their meteorological conditions, which affected the CO2 and CH4 exchange between the at- mosphere and ecosystems. The early growing season in 2017 was colder than the corresponding period in 1981–2010 on average (Pirinen et al., 2012), with May being 1.9 ◦C and June 1.4 ◦C colder (Fig. 3). The summer months of June– August were rainy in 2017, which was reflected in cloudi- ness that decreased the incoming solar radiation, with July in particular being cloudy (Fig. 4). In contrast, the spring of 2018 was warm, with May being 3.8 ◦C warmer than the 30-year May average. In July 2018, a widespread drought and heatwave event in northwestern Europe reached Kaa- manen and caused the monthly mean air temperature to rise to 18.6 ◦C, which was 5.3 ◦C higher than the 30-year July av- erage. Additionally, the monthly precipitation sum was less than half of the 30-year average (34 and 72 mm, respec- tively). After the drought, the following August and Septem- ber were also warmer than the corresponding 30-year aver- ages, by 1.5 and 1.9 ◦C, respectively, but the precipitation sum was then double the 30-year average. Spring 2019 was characterised by a dry and warm April (3.7 ◦C above the long-term average) and rainy weather in May and June. Lake water levels were visibly lower in 2018, exposing the shore- line sediments of shallow lakes, for instance, in the MS lake of this study. The forest soil freeze and thaw occurred each year at the end of October and in mid-May, respectively (Fig. 5a). In late autumn 2018, the forest Ts decreased from 8 to 1 ◦C in just 20 d between 23 September and 13 October, which was due to the lack of snow cover. In the fen, the surface peat of the strings froze concurrently with the forest soil (Fig. 5b). How- ever, the thaw occurred slightly earlier in the strings than in the forest, where trees shaded the surface. The continuous water saturation in the flark peat led to weaker peat tempera- ture responses to changing air temperature than in the dryer strings. The lake water temperature, measured only during the ice-free periods of 2017 and 2018 (Fig. 5c), was affected, in addition to air temperature, by the incoming flow to the lake and by the water level. Following the air temperatures, the forest soil, fen peat, and lake water temperatures were higher in July 2018 than in 2017. The ice-free period of the measured lakes lasted for 167 and 176 d (from May to October) during the first and second study years, respectively (Table 4). The meteorological drought lasted from 2 July to 1 August 2018 (Fig. 6a). The drought limit of daily VPDmax= 20 hPa was exceeded on 13 d in total in 2018, while no exceedances took place in 2017. The average daily maximum VPD in 2 July–1 August 2018 was 16.9 hPa, while, during the same https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 554 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape Figure 4. Weekly average photosynthetic photon flux density (PPFD) measured at Kaamanen. Figure 5. Daily average (a) forest soil temperature measured at 10 cm depth and snow depth, (b) fen peat temperature measured from strings and flarks at 10 cm depth, and (c) MS lake water temperature at the floating chamber position 10 cm from the bottom. Table 4. Estimated ice-free period start and end dates for the studied lakes. Period start Period end 25 May 2017 26 Oct 2017 10 May 2018 30 Oct 2018 5 May 2019 19 Oct 2019 period in 2017, it was 9.4 hPa. The corresponding daily air temperatures were 18.5 and 14.2 ◦C, respectively. The high- est daily maximum VPD of 31 hPa was observed on 18 July 2018. VPD responded rapidly to the reduced amount of water in the environment, while the water table level in the fen, the water depth of lakes, and the forest soil moisture decreased gradually as the drought developed (Fig. 6). Within the fen, the WTL varied strongly among the plant community types (Fig. 6b). In the flark and Trichophorum tussock PCTs, WTL was close to or at the peat surface most of the time, while in elevated strings it was deeper, away from the peat surface, fluctuating between −10 and −20 cm in string margins and between −30 and −70 cm in string tops. The drought decreased the WTL in all plant communities si- multaneously, and by mid-August, the WTL had recovered to a normal level, except for the string tops. The water depth of the MS lake was within 30–40 cm (Fig. 6c), and it dropped in June–July 2018 to less than 20 cm at the centre of the lake. This drawdown was associated with a shoreline retreat of 6 to 8 m. However, after the first rain- Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 555 Figure 6. (a) Daily maximum vapour pressure deficit. (b) Mean water table level at the main plant communities of the fen. Error bars represent the standard deviation. (c) Water table level at the floating chamber position in Lake Jänkälampi. (d) Forest soil moisture at 10 cm depth. Figure 7. The 1-month (July) and 5-month (May–September) Standardised Precipitation Evapotranspiration Index (SPEI) at Kaamanen in 1950–2018. Positive SPEI values denote moist conditions and negative values drought conditions. The years 2017 and 2018 are marked with shading. https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 556 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape fall in 3–5 August, the water depth quickly reverted to the pre-drought level. Our field observations and visual interpre- tation of satellite and drone imagery suggest that the water level drop magnitude was uniform in the waterbodies within the study area. The annual cycle of the soil moisture measured in the pine forest showed a low moisture content (< 0.1 m3 m−3) during winter and fluctuations after the soil thaw (Fig. 6d). During both springs, the water availability in soil was at its maxi- mum in 10–14 May due to melting snow. The drought during summer 2018 was clearly indicated by the soil moisture data that showed a drastic drop from 0.2 to 0.1 m3 m−3 between 23 June and 22 July. After the rainfall in the beginning of August, the soil moisture rose in a few weeks to the average growing season level. Thus, even though the meteorological drought was already over at the beginning of August, the wa- ter availability to vegetation was not yet fully recovered. In terms of the drought index (SPEI), the drought event in July 2018 was the eighth most severe drought in Kaamanen between January 1950 and December 2018 (Fig. 7). How- ever, the drought event was relatively short, and the 5-month SPEI between May and September was close to neutral. 3.2 Ecosystem CO2 and CH4 fluxes 3.2.1 Pine forest fluxes The pine forest acted as a net CO2 sink during both study years. The annual CO2 balances (with the 95 % confidence intervals) during the first and second study years were −126± 21 and −101± 19 gCm−2, respectively (Fig. 8b; Table A6). For other evergreen needleleaf forests in north- ern Fennoscandia, the observed balances have been smaller in magnitude; the Scots pine forest at Värriö in northern Fin- land was a CO2 sink in 2012–2014 (−48 to −7 gCm−2) and a small source in 2015 (14 gCm−2; Kulmala et al., 2019), while the Norway spruce forest at Kenttärova in northern Fin- land had a close-to-neutral balance (on average −2 gCm−2; Aurela et al., 2015). There were two periods when the CO2 fluxes behaved dif- ferently between the 2 study years, due to differing mete- orological conditions. First, even though the summer was warmer in 2018 than in 2017, and thus one could expect en- hanced respiration, the ER fluxes of the pine forest ecosystem were on average similar during the growing seasons. This was due to the rainy period in June–August 2017 (Fig. 3), as a result of which the forest soil remained saturated or nearly saturated in water during the whole growing season (Fig. 6d). As the forest soil moisture content was continuously close to the maximum water-holding capacity, the effect of abun- dant precipitation emerged as increasing lake water table levels (Fig. 6c). We suspect that the stronger ER tempera- ture response observed in June–August 2017 was caused by enhanced heterotrophic soil respiration (Fig. A1 in the Ap- pendix), which is known to increase with soil moisture until near-saturation with water (Orchard and Cook, 1983; Moy- ano et al., 2012; McElligott et al., 2017; Du et al., 2020). The second period of dissimilar behaviour in CO2 fluxes was observed when the drought and heatwave discussed above limited GPP fluxes; compared to the previous year, the GPP sum in 22 July–17 August was 35 gCm−2 lower in 2018 (Z test; p= 0.003; Fig. 8d). The daily maximum VPD surpassed the 20 hPa limit, indicating meteorological drought, for the first time on 2 July and the last time on 1 August 2018 (Fig. 6a), during which period the average air temperature was 5 ◦C higher than in the previous year (Fig. 3). In July 2018, the forest soil moisture at 10 cm depth also dropped by 50 % from the normal growing season level, decreasing to 0.1 m3 m−3 on 22 July (Fig. 6d). Soil moisture recovered to a normal level 3 weeks later on 13 August. There is an evident discrepancy between the timing of the meteorological drought and the reduced ecosystem CO2 fluxes. The response to drought in plants occurs both in roots, where water availability through soil is a key fac- tor, and in the stomatal openings in leaves, where the gas exchange is affected by stomatal control due to VPD. The tighter stomatal control is due to the high VPD limits C as- similation (Martín-Gómez et al., 2017), and in Scots pine trees, this has been estimated to occur at a VPD of 8 hPa, while at approximately 28 hPa the stomata are fully closed (Büker et al., 2012). The drought seemed to have an effect on the CO2 exchange until 18 August, even though soil mois- ture and VPD hardly had any direct effect anymore, as only after this date were the magnitude of ER and GPP fluxes not continuously lower compared to the previous year. This lagged effect could be caused by embolism, defoliation, or root degradation (Aguadé et al., 2015) or by weakened my- corrhizal symbiosis in the roots (Muilu-Mäkelä et al., 2015). Gao et al. (2017) also found that the GPP of a southern boreal Scots pine forest was suppressed during a severe soil mois- ture drought in the summer of 2006, as the plants regulated their stomata due to high VPD. According to the linear regression model, the main en- vironmental factors affecting the radiation-normalised GPP flux (FGPP1200) of the pine forest were Ts, VPDmax and soil moisture (Table 5). The magnitude of FGPP1200 increased with increasing Ts and soil moisture and decreased with in- creasing VPDmax. The positive correlation with soil moisture originated from the anomalously low moisture levels during the 2018 drought (Fig. 6d). The temperature-normalised res- piration (FR10) was affected by Ts, VPDmax, soil moisture, and precipitation sum (Table 5). FR10 increased with increas- ing Ts, soil moisture, and precipitation sum and decreased with increasing VPDmax. 3.2.2 Fen fluxes The fen ecosystem was a small net CO2 sink of −14± 17 gCm−2 and a small net CO2 source of 9± 11 gCm−2 during the first and second years, re- Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 557 Figure 8. Pine forest gap-filled fluxes. Weekly averaged (a) and cumulative (b) CO2 flux. Weekly averaged ecosystem respiration (c) and gross primary productivity (d). Figure 9. Fen gap-filled fluxes. Weekly averaged CO2 (a), ER (c), GPP (d), and CH4 (e) flux. Cumulative CO2 (b) and CH4 (f) flux. https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 558 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape Table 5. Standardised regression coefficients (± standard error) for the explanatory variables of the linear regression model for EC- based FGPP1200 and FR10. The adjusted coefficient of determina- tion (R2 a ), normalised root mean square error (nRMSE; RMSE di- vided by the range of response variable values), and degrees of free- dom (df) are also shown. Regression df coefficient FGPP1200 (R2 a = 0.78; nRMSE= 0.12) Ts (−10 cm) −1.02± 0.11 34 VPDmax 0.28± 0.11 34 Soil moisture (−10 cm) −0.17± 0.08 34 FR10 (R2 a = 0.68; nRMSE= 0.15) Ts (−10 cm) 0.87± 0.18 31 VPDmax −0.78± 0.18 31 Soil moisture (−10 cm) 0.24± 0.12 31 Precipitation sum 0.21± 0.11 31 spectively (Fig. 9b; Table A6). There were three periods during which either CO2 or CH4 fluxes diverged between the years, namely the start of the growing season, the warmer- than-average growing season of 2018, and the drought and heatwave event in 2018. The warm growing season increased the ER flux sum by 32 gCm−2 (Z test; p= 0.124) during 11 June–23 Septem- ber 2018 (mean air temperature 13.2 ◦C) compared to 2017 (10.5 ◦C). Half of this difference, 16 gCm−2 (Z test; p= 0.043), accumulated in just 26 d, 17 July–12 August, when the temperature difference between the years was largest (13.4 and 17.9 ◦C in 2017 and 2018, respectively). The earlier start for the growing season in 2018 resulted in a higher CO2 uptake in 11–30 June compared to the pre- vious year (Fig. 9; Table A6); the balances of this period were −5± 13 and −20± 17 gCm−2 in 2017 and 2018, re- spectively. This was due to the nearly doubled GPP (Z test; p< 0.001), which was 27± 3 and 49± 4 gCm−2 in 2017 and 2018, respectively. The increase in the net CO2 uptake in northern mires due to earlier snowmelt and warm spring tem- peratures has been reported previously by, for example, Au- rela et al. (2004) and Sagerfors et al. (2008). The differences in CO2 exchange between the microforms of the fen were studied in more detail by Heiskanen et al. (2021), who found that the magnitude of ER and GPP increased gradually from the wettest PCT, i.e. flark, to the driest one, i.e. string top. However, the variation in the net CO2 uptake among the mi- croforms was weaker than in the two flux components. The increased GPP due to warm spring weather was observed in all main microforms. However, the simultaneous increase in ER led to a significantly (p< 0.05) higher net CO2 uptake only in string tops, when comparing the early growing sea- sons of 2017 and 2018 (Fig. A2). The higher CO2 uptake during the early growing season of 2018 was offset by the decreased uptake due to the drought and heatwave event in 8 July–4 August 2018. The GPP sum was 17 gCm−2 smaller (Z test; p= 0.229) during this pe- riod in 2018 than in 2017 (Fig. 9d). As discussed above, the drought was observed as higher-than-average temperatures (Fig. 3), an elevated VPD (Fig. 6a), and water level draw- down by 5–20 cm at the fen microforms (Fig. 6b). These anomalies likely caused drought stress in the mire plants (Alm et al., 1999). The drought stress reduced GPP, and the high temperatures increased ER compared to the previous weeks and the same time period during the previous year (Fig. 9c and d). The CO2 uptake decreased during this period, and at the end of the period the fen ecosystem even turned into a CO2 source. Unlike the increased CO2 uptake during the early growing season, which could be allocated largely to the string plant communities, the drought affected plant communities in all microforms (Heiskanen et al., 2021). The CO2 exchange in flarks, including Trichophorum tussocks, was immediately affected by the lowering WTL in July and August 2018 (Fig. 6b), while the drier string communities were affected to a lesser extent. While the fen was an annual net CO2 sink, it also acted as a CH4 source to the atmosphere during both study years. The annual CH4 balance was 7.8± 0.2 and 6.3± 0.3 gCm−2 during the first and second study years, respectively (Fig. 9f). The flark, Trichophorum tussock, and string margin PCTs contributed 98 % of the emissions, with string margins ac- counting for 44 % of the emissions during the growing sea- son of 2017 and all three having similar emissions in 2018 (Heiskanen et al., 2021). The drought decreased the annual CH4 emissions, mostly during 21 July–28 August 2018, when the emissions were 0.8 gCm−2 lower (Z test; p< 0.001) than during the same period in the previous year (Fig. 9e). Notably, the decrease in CH4 emissions occurred a few weeks after the meteoro- logical drought begun and continued well after the WTL had reverted to the pre-drought level. The decrease in CH4 emis- sions is likely due to both the reduced release of carbon com- pounds by plant roots and increased oxic soil zone, which re- duced CH4 production and increased CH4 oxidation (Strack and Waddington, 2007; White et al., 2008; Deppe et al., 2010). The lagged recovery of CH4 flux coinciding with the GPP recovery indicates the link between methanogen- esis and plant root exudates. Unfortunately, the uncertain- ties in the monthly CH4 balances derived from the manual chamber measurements were large, and the lower emissions could not be allocated to any specific PCT. However, the CH4 emissions from flarks were larger in 2018 than 2017, which suggests that the drought-induced decrease in emis- sions occurred for the Trichophorum tussock and string mar- gin PCTs that had a lower WTL (Fig. 6b; Heiskanen et al., 2021). The drought, which covered large parts of northwest- ern Europe, did not affect the CO2 and CH4 exchange at Kaa- manen as much as it did at more southern mires, where the Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 559 Table 6. Modelled annual CO2 and CH4 balances (gCm−2) of the MS and OS lakes during the 2 study years (Kou et al., 2022). 11 Jun 2017– 11 Jun 2018– 10 Jun 2018 10 Jun 2019 Mineral sediment lake CO2 flux 24.4 32.0 CH4 diffusive flux 1.0 1.1 CH4 ebullition flux 0.3 0.3 Total C flux 25.6 33.4 Organic sediment lake CO2 flux 99.4 94.2 CH4 diffusive flux 2.9 3.3 CH4 ebullition flux 0.8 0.9 Total C flux 103.1 98.4 drought duration was longer, although the water level draw- down was similar in magnitude (Rinne et al., 2020). The total carbon balance, i.e. the sum of the CO2 and CH4 fluxes, showed that the fen was an annual carbon sink, −7± 17 gC m−2, and a carbon source to the atmo- sphere, 15± 11 gC m−2, in the first and second study years, respectively (Table A6). The lower net CO2 uptake dur- ing the drought period contributed most of the difference between the years. Previous studies at the Kaamanen fen show that the ecosystem has been on average a larger an- nual CO2 sink, −22 gCm−2 (from −4 to −53 gCm−2) dur- ing 1997–2002 (Aurela et al., 2004), and a similar average annual CH4 source of 6 gCm−2 in 1995, 1997–1998, and 2011–2016 (Hargreaves et al., 2001; Piilo et al., 2020). 3.2.3 Lake fluxes The measured CO2 fluxes of the OS and MS lake were on average 0.45 and 0.06 gCm−2 d−1 during the ice-free period, and the mean CH4 fluxes were 0.014 and 0.003 gCm−2 d−1, respectively (Fig. 10). The large difference in the flux rates between the lakes was most likely caused by the fact that the OS lake is situated immediately downstream of the fen (Fig. 1). Due to the location, the OS lake receives transported organic carbon (OC) in the streamflow and has a higher lake water total organic C content than the MS lake (9 vs. 5 mgL−1) and a higher C content in its sediment (30.3 % C of dry mass) than in the MS lake sediment (5.5–22.6 % C, de- pending on location; Table A8). Based on the ALBM model calculations (Kou et al., 2022), calibrated with the flux ob- servations at the MS and OS lakes (Sect. 2.5), the annual CO2 emissions were on average 28 gCm−2 from the MS lake and 97 gCm−2 from the OS lake. The modelled CH4 emis- sions were 1.4 gC m−2 from the MS lake and 4.0 gCm−2 from the OS lake (Table 6). The CO2 emissions of the OS lake (0.2–0.8 gCm−2 d−1; Fig. 10a) were similar in magnitude, for instance, to those (0.3–1.4 gCm−2 d−1) observed for Lake Stortjärn in Swe- den, which is a northern boreal lake located adjacent to a mire (Denfeld et al., 2020). In lakes with less organic car- bon in the water, the net CO2 emissions are lower, yet we did not observe CO2 uptake by either of the studied lakes. In a similar shallow subarctic lake next to a fen at Abisko in Sweden, Jammet et al. (2017) found high CO2 emissions of 33.3 gCm−2 during the spring period of 41 d, but the photosynthetic CO2 uptake during summer reduced the net emission of the ice-free period to 8.9 gCm−2. Lohila et al. (2015) estimated an annual CO2 balance of 33 gCm−2 for the shallow parts of the large subarctic Lake Pallasjärvi in Finland and recorded a small CO2 uptake during midday in the summer months. Also, small annual CO2 emissions, 11.5 gCm−2, were estimated for the small (9.6 ha) Lake Kipojärvi located near Kaamanen, surrounded by both an es- ker and a peatland (Juutinen et al., 2013). This shows that the CO2 balance of lakes varies substantially from lake to lake in the subarctic region. The annual CH4 emissions of Lake Kipojärvi, 3.4 gCm−2, were similar to the modelled OS lake emissions. The av- erage daily diffusive CH4 fluxes varied from 0.001 to 0.020 gCm−2 d−1 in this study and were similar to those reported in previous studies of boreal lakes. Denfeld et al. (2020) observed a range of 0.001–0.008 gCm−2 d−1 for Lake Stortjärn during the ice-free period, while Rasilo et al. (2015) found spatially highly variable CH4 diffusive fluxes of 0.008± 0.020 gCm−2 d−1 across 224 boreal lakes in Canada during summer. The length of the ice-free period can affect the annual emission of different boreal lakes considerably, as the ice cover prevents gas exchange with the atmosphere (Guo et al., 2020). Our estimate of ebullition was 22 % of the total CH4 emis- sion. This is smaller than observed in other studies of bo- real lakes, especially considering that our study lakes are mostly shallow (< 2 m deep) and thus likely to show more frequent CH4 ebullition than deep lakes (Bastviken et al., 2004). Contributions of 40 %–80 % of the total CH4 emis- sions have been estimated for subarctic lakes (Bastviken et al., 2004; Wik et al., 2013; Jansen et al., 2020), and the an- nual CH4 emissions of 1.0–8.3 gCm−2 have been observed (Bastviken et al., 2004; Thornton et al., 2015; Jammet et al., 2017). In general, there is high temporal and spatial varia- tion in ebullition, as there are both seep and non-seep areas where CH4 bubbles emerge at the lake bottom, which com- plicates the ebullition estimation (Walter et al., 2006). Wik et al. (2018) found that the spatial ebullition potential was af- fected by the coarse detritus, buried aquatic vegetation, and redeposited peat rather than the amount of total organic car- bon or CH4 in the sediment. Both lakes showed the largest CO2 emissions in early June (Fig. 10a). However, a similar development of fluxes was not observed with CH4 (Fig. 10b). The dissimilarity between the https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 560 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape Figure 10. Weekly average lake CO2 (a) and CH4 (b) diffusive fluxes. Error bars represent the standard deviation of individual chamber measurements in each week. Figure 11. Treed pine bog gap-filled fluxes. Weekly averaged CO2 flux (a) and cumulative CO2 flux (b). Weekly averaged ER (c) and GPP (d) flux. gases is likely due to the differing production and accumu- lation processes. CO2 emissions are driven by the amount of incoming dissolved and particulate carbon and the rate of decomposition of that organic carbon (Kortelainen et al., 2006), while CH4 production is controlled by the substrates available for the methanogens to reduce it to CH4 in anoxic conditions in sediments (Chapin et al., 2011). CH4 concen- tration has been observed to be roughly constant with depth during the early summer, while CO2 is strongly stratified, with a larger concentration towards the bottom of the lake (Denfeld et al., 2020). The breakdown of thermal stratifica- tion and turnover mixing occurs around the start and end of the ice-free period (e.g. López Bellido et al., 2009). At the MS lake, the flux measurements also covered the au- tumn turnover mixing in September, when CO2 emissions increased compared to the previous month. The daily aver- age CO2 flux on the MS lake was lower in July 2018 than in 2017 (Welch’s t test; p= 0.01; Fig. 10), but no significant changes were observed in the measured CH4 fluxes between the years. Both lake types emitted 96 % of the total C efflux as CO2 (Table 6). For comparison, in an extensive study of Alaskan Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 561 subarctic lakes, the non-yedoma lakes emitted about 85 % as CO2 and 15 % as CH4 (Sepulveda-Jauregui et al., 2015), and in the above-mentioned 224-lake study, the CH4 diffu- sive fluxes contributed 8± 23 % to the total lake C emissions (Rasilo et al., 2015). 3.2.4 Pine bog fluxes The treed pine bog ecosystem between the forest and fen ecosystems covers 9.1 % of the studied landscape, and the sparsely treed pine bog covers 13.5 % and is located within the fen ecosystem in the driest parts of peatland. The an- nual CO2 balance of the treed pine bog ecosystem was similar in the first and second study years, −92± 102 and −92± 110 gCm−2, respectively (Fig. 11b; Table A6). The corresponding estimates for the sparsely treed pine bog ecosystem were lower and more uncertain at −48± 134 and −68± 145 gCm−2 (Fig. 12b; Table A6). The larger uncer- tainties in these balances compared to the forest and fen ecosystems were due to the contribution of the modelled string top fluxes, which were based on a limited number of chamber measurements (Sect. 2.4). For both pine bog ecosystems, the GPP rates were higher during the latter year from mid-June to mid-July (Figs. 11d and 12d). This led to larger net CO2 uptake during the first part of growing season in 2018 than in 2017 (Figs. 11a, b and 12a, b). However, this difference in fluxes was more prominent in the sparsely treed pine bog. The drought event decreased the magnitude of the GPP rate similarly in both ecosystems between mid-July and mid-August 2018 com- pared to the previous 2 summer months and the same period in 2017 (Figs. 11d and 12d). The pine bog ecosystem fluxes were derived from the pine forest and fen ecosystem fluxes. String tops, which were used as a proxy for a treeless pine bog, are relatively dry microsites, with a typical water table depth of 40–60 cm (Fig. 6b). As part of the bog ecosystems are wetter than the typical string tops, and as the flark PCTs at Kaamanen had lower CO2 fluxes than the string tops (Heiskanen et al., 2021), the proxy approach may bias the pine bog flux esti- mates. However, the use of string margin fluxes instead of the string top fluxes would not make any significant differ- ence. The fluxes of the pine bog ground layer have been studied previously within a nearby Kipojärvi catchment by Juutinen et al. (2013), who found that in 2006 the ground layer was an annual CO2 sink of −130± 91 gCm−2, with an annual GPP sum of −456± 77 gCm−2 and ER sum of 326± 48 gCm−2. The ER sum was similar to that measured at the Kaamanen fen on the string top and margin PCTs (Heiskanen et al., 2021), but the magnitude of the GPP sum and thus the CO2 sink were greater for the Kipojärvi catch- ment (Fig. A2b). This suggests that our pine bog flux esti- mates might underestimate the true CO2 sink of this ecosys- tem. 3.3 Upscaled landscape-level fluxes By upscaling the ecosystem balances to our study area of 7 km2, we obtained an annual landscape-wide CO2 bal- ance of −45± 22 and −33± 23 gCm−2 for the 2 study years. The corresponding CH4 balances were 3.0± 0.2 and 2.7± 0.2 gCm−2, and the sum of CO2 and CH4 bal- ances were −42± 22 and −31± 23 gCm−2, respectively (Fig. 13b, f; Table A6). The pine forest ecosystem con- tributed to the total landscape C balance with a large CO2 sink and a minimal CH4 sink. In the evergreen forest ecosystem, the CO2 uptake period was longer than that of the fen ecosystem by 12 d at the beginning of the growing season and by 30 d at the end of the growing season (Figs. 8d and 9d), and during the first half of the growing season, the magnitude of CO2 uptake was larger in the forest (Figs. 8a and 9a). The lake ecosystem was the complete opposite, with substantial CO2 and CH4 emissions, while the fen ecosys- tem balances were between forest and lake, with annual net CO2 uptake and CH4 emissions, although the fen had larger CH4 emissions per unit area than the lakes. The pine bog ecosystems, which were modelled using pine forest and fen fluxes, most likely acted as CO2 sinks and CH4 sources. The ecosystem-specific CH4 emissions were much lower than the magnitude of CO2 exchange, with the exception of the fen ecosystem. The fen CH4 fluxes largely determined the landscape-level CH4 emissions (Fig. 14). The average an- nual terrestrial C uptake of the landscape was 338± 156 t C, of which 24 % was released back to the atmosphere by the lakes, resulting in a net C uptake of 256± 156 t C within the landscape (Fig. 14). The studied landscape consisted of roughly an equal area of pine forest, pine bog, open peatland, and lakes (Table 1). Kou et al. (2022) studied the effect of the land-cover-type misclassification to the ecosystem and landscape C balances at Kaamanen using biogeochemical and atmospheric mod- els. Their study area encompassed the entire 32.8 km2 catch- ment at Kaamanen. In a related manner, the 7 km2 area of this study partially coincided with the most southern part of the catchment. The catchment consisted of proportionally more pine forest (53 %) and less open peatland (16 %), fewer lakes (13 %), and fewer pine bogs (9 %). In contrast to this study, birch-dominated ecosystems (1 %) and mixed forests (6 %) were also present in the area. Thus, the 7 km2 study area represented almost all ecosystems found in the greater landscape area around Kaamanen. The average annual land- scape CO2 balance was −104 gCm−2 and the CH4 bal- ance 2.4 gCm−2. The larger annual landscape CO2 uptake on catchment-scale highlights the effect of the proportion- ally larger forested area that was effective in CO2 sequester- ing. Aurela et al. (2015) estimated an annual landscape-level CO2 balance of −5 gCm−2 for a 1963 km2 area in Pallas in northwestern Finland that was comprised of 71 % pine for- est, 12 % open wetland, 6 % water surfaces, and 2 % tree- less fell tops. Within this area, the 105 km2 Lake Pallasjärvi https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 562 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape Figure 12. Sparsely treed pine bog gap-filled fluxes. Weekly averaged CO2 flux (a) and cumulative CO2 flux (b). Weekly averaged ER (c) and GPP (d) flux. Figure 13. Landscape-level gap-filled fluxes. Weekly averaged CO2 (a), ER (c), GPP (d), and CH4 (e) flux. Cumulative CO2 (b) and CH4 (f) flux. Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 563 Figure 14. Annual CO2, CH4, and C flux sums for different ecosystems and the landscape (Fig. 1) upscaled with the corresponding area. The CH4 flux estimate for forest is from Dinsmore et al. (2017) and for pine bog from Bubier et al. (2005). The error bars denote the 95 % confidence interval. Uncertainty estimates were not available for the modelled lake balances (Kou et al., 2022). catchment included fewer forests and peatlands, but more lake and fell area, and the annual landscape-level CO2 bal- ance was estimated to be 15 gCm−2. The annual CO2 up- take in Pallas was lower than in Kaamanen in spite of the much larger proportion of forests. In a study temporally over- lapping our study, Chi et al. (2020) reported landscape-level C balances for a typical northern boreal forest landscape in Svartberget, Sweden, where forests covered 87 % of the area; the annual landscape-level CO2 balance was −37 gCm−2 in October 2016–September 2017 and −108 gCm−2 in Octo- ber 2017–September 2018. The latter year showed a drought response similar to Kaamanen; a higher annual CO2 up- take was observed in the forest, while the mire ecosystem (9 % of the area) turned from an annual CO2 sink to a source. The annual CO2 balances of Svartberget were comparable to Kaamanen, whose pine forest had an annual CO2 bal- ance of −78 and −119 gCm−2 during the first and sec- ond study years (Fig. 8d). In a synthesis by Virkkala et al. (2021), the mean annual NEE of 41 boreal biomes was on average −46 gCm−2 (uplands −47 gCm−2 and wetlands −38 gCm−2), i.e. the CO2 sink was lower in uplands and larger in wetlands than the mean landscape-scale sink in Kaa- manen. These results not only show that there can be large spatial variation in C exchange among boreal ecosystems, e.g. forests, but also part of the variation between landscape C balances derive from the ecosystem composition. The temporal variation in the landscape-level CO2 flux (Fig. 13) obviously depended on the flux variation in in- dividual ecosystems, and different ecosystems showed dif- ferent environmental responses. During the 2 study years, there were three periods when the ecosystem-specific fluxes clearly deviated from each other, and these differences were reflected in the landscape-level C balance. First, the warmer- than-average early growing season in 2018 increased the CO2 uptake of the fen ecosystem (Fig. 9a) and probably also the pine bog fluxes (Figs. 11a and 12a), or conversely, the colder-than-average early growing season in 2017 led to lower CO2 uptake by these ecosystems. However, equally large variation in the early growing season CO2 uptake was not observed in the pine forest fluxes, as the evergreen pine forest phenology differs from that of the largely deciduous mire. Second, the rainy peak growing season in 2017 in- creased the ER rates of the pine forest ecosystem compared to the same period in the next year, which was much warmer, and thus, the ER fluxes were similar (Figs. 8 and A1). Con- versely, at the fen, the warm summer of 2018 increased the ER sum compared to the cool, rainy summer in 2017. This was due to the inherently different water balance between uplands and peatlands. Third, the drought period in summer 2018 decreased the CO2 uptake in both the pine forest and fen ecosystems (Figs. 8b and 9b), and the pine bogs likely showed a similar response. https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 564 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape The lakes in the landscape acted as CO2 sources to the at- mosphere throughout the growing seasons, with peak emis- sions linked with the spring thaw. However, the temporal variation in the lake emissions had only a minor effect on the temporal variation in the landscape-scale CO2 fluxes, as the flux magnitude did not vary as much on the lakes as on the terrestrial ecosystems during the growing seasons (Figs. 8, 9, and 10). Integrated over the study area, the total lake emis- sions had a considerable effect on the landscape-scale annual balances (Fig. 14). As the pine forest and pine bog CH4 fluxes used for up- scaling were adopted from the literature (Dinsmore et al., 2017; Bubier et al., 2005), only the fen and lake ecosys- tem CH4 fluxes affect the annual variation in the upscaled landscape-level fluxes (Figs. 13 and 14). The CH4 emis- sions from the fen decreased due to the drought in 2018 (Fig. 9e), thus decreasing the annual emissions compared to the first study year (Table A6). However, as the modelled lake CH4 emissions were slightly higher during the second year (Table 6), there was only a minor difference in the annual landscape-level CH4 balance (Fig. 13f; Table A6). During the 2 study years, the meteorological conditions were not optimal for C sequestration, as the fen on average acted as a weaker CO2 sink than in some previous years (Au- rela et al., 2007) and even a source in the second year. Thus, the landscape-level CO2 uptake would probably be higher in more optimal conditions. However, the conditions favour- ing C sequestration differ between ecosystems, as all ter- restrial ecosystems are likely to sequester more CO2 during longer growing seasons, but mires can simultaneously emit more CH4, and lakes emit more of both CO2 and CH4. For estimating the relative radiative impact of the CO2 and CH4 exchanges in Kaamanen, the CH4 flux was translated to CO2-equivalent flux by multiplying it by the sustained global warming potential (SGWP) coefficient (Neubauer and Megonigal, 2015; Neubauer, 2021). If the CO2 and CH4 flux observed during the 2 study years continued for the next 500 years (SGWP= 14 gCO2 eq.gCH4 −1), then the net CO2-equivalent flux of the fen would be positive. However, in longer timescales, the fen will eventually have a negative radiative balance, even with a low annual CO2 uptake. In ad- dition to the time-dependent contribution of fens, the annual net CO2-equivalent flux of the landscape was affected by up- take of both CO2 and CH4 by forests, i.e. a systematically negative CO2-equivalent flux, and emissions of both CO2 and CH4 from lakes, i.e. a systematically positive CO2- equivalent flux. Assuming the present-day fluxes for each ecosystem, upscaling suggests that the net CO2-equivalent flux of the landscape would be initially positive but turn neg- ative soon after 100 years (SGWP= 45 gCO2 eq.gCH4 −1 for the 100-year time horizon). 4 Conclusions We estimated the ecosystem–atmosphere exchange of CO2 and CH4 for the main ecosystems in a subarctic landscape during 2 full years. The 7 km2 study area consisted of pine forest, patterned flark fen, two pine bog ecosystems, and two lake types. For the terrestrial ecosystems, C exchange was most sensitive to changes in temperature and moisture con- ditions, while on the lake it depended on the amount of avail- able carbon in sediment and the length of the ice-free period. The lakes in the study area released 24 % of the C that was sequestered by the landscape during the 2-year study period back to the atmosphere, and there was a 4-fold difference in the CO2 emissions between organic- and mineral sedi- ment lakes. Thus, more measurements are needed to accu- rately define the role of lake emission variability. Similarly, the CH4 fluxes were much greater from the organic sediment lake, but the overall impact of lake fluxes on the landscape scale was smaller than for CO2. Obviously, the great differ- ence in the observed C fluxes between the lake types should be considered when estimating regional-scale fluxes in a het- erogeneous environment such as a northern boreal landscape. There were three periods when the C fluxes of the terres- trial ecosystems were clearly different between the 2 study years due to meteorological conditions. In the pine forest, the CO2 fluxes were affected by the rainy weather in sum- mer 2017 increasing the ER rate to the level observed dur- ing the warmer growing season in the following year. At the fen, however, the warm growing season in 2018 resulted in a higher ER sum than in the previous growing season, high- lighting the differing water balance of the fen and upland forest ecosystems. The warmer-than-average early growing season in 2018 advanced the plant growth at the fen, thus increasing the CO2 uptake of this ecosystem. All terres- trial ecosystems were affected by a short but severe drought event in July 2018, which decreased the GPP rates and thus CO2 uptake. However, both the onset of and recovery from drought effects occurred more rapidly at the fen than in the pine forest. Additionally, during the drought, the CH4 emis- sions from the fen decreased due to water level drawdown and possibly also due to decreased plant root carbon input. The CO2 flux responses to changing environmental condi- tions were reflected in the landscape-level CO2 fluxes, even though only the drought affected the CO2 fluxes of all ter- restrial ecosystems similarly. For this reason, none of the ecosystems alone controlled the changes in the landscape- level CO2 exchange. In contrast to the CO2 fluxes, it appears that the landscape-level CH4 flux and its variation can almost entirely be estimated based on the fen data. Both study years had periods of nonoptimal C sequestra- tion conditions, but still the landscape remained a CO2 sink. This indicates that the multitude of ecosystems contributes positively to the landscape resilience to C loss. Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 565 Appendix A: Ancillary environmental information, flux, and statistical tests Table A1. Ground layer biomass and leaf area index (mean± standard deviation) and the number of field measurement points. Ground layer Biomass (g of dry matter per metre squared) Leaf area index (m2 m−2) n Land cover class Shrubs Forbs and Mosses and Total Shrubs Forbs and Total graminoids lichens graminoids Pine forest 144.1± 73.0 3.2± 7.8 398.1± 200.7 545.4± 189.5 0.39± 0.15 0.02± 0.04 0.41± 0.17 75 Pine bog 284.9± 225.6 38.2± 30.2 398.4± 172.9 721.6± 324.5 0.51± 0.22 0.22± 0.16 0.73± 0.25 36 Fen 149.4± 90.9 69.6± 36.9 146.7± 169.1 365.6± 231.2 0.21± 0.10 0.37± 0.19 0.59± 0.21 223 String top∗ 247.5± 147.4 39.5± 29.9 261.8± 263.7 548.8± 318.8 0.50± 0.20 0.21± 0.14 0.71± 0.19 49 ∗ String top CO2 fluxes were used for estimating the ground layer flux of pine bogs. Table A2. Pine forest basal area, stand density, and diameter at breast height (DBH; mean± standard deviation) derived from field measure- ments. Tree species Basal area No. of trees per ha DBH (cm) Scots pine 13.2± 6.4 545± 827 15± 6 Downy birch 0.3± 0.7 103± 243 7± 6 Table A3. Average tree height (mean± standard deviation) and the number of field measurement points. Ecosystem Average tree height (m) n Pine forest 7.7± 2.0 43 Pine bog 4.9± 2.0 17 Table A4. Soil and peat properties (mean± standard deviation) for pine forest, pine bog, and fen ecosystems. Except for the pine forest organic layer, pH was always measured at 30 cm depth. In pine bog and fen ecosystems, bulk density, C and N content, and C : N ratio are the mean of 0–5 and 15–20 cm peat layers. Ecosystem pH Bulk density Soil C content Soil N content C : N ratio No. of sample (gcm−3) (mgcm−3) (mgcm−3) plots Pine forest 6.0± 0.4 11 Organic layer 4.6± 0.2 0.136± 0.101 52.9± 44.9 1.2± 0.9 43.5± 8.4 Eluvial layer 1.037± 0.265 20.6± 11.9 0.5± 0.3 45.7± 11.2 Illuvial layer top 1.219± 0.110 16.5± 10.4 0.5± 0.4 36.4± 9.3 Illuvial layer bottom 1.337± 0.145 5.0± 4.1 0.2± 0.2 43.8± 40.1 50 cm below organic layer 1.511± 0.124 2.1± 3.3 0.1± 0.2 42.4± 30.6 100 cm below organic layer 1.588± 0.100 1.1± 0.9 0.01± 0.02 256± 506 Pine bog 5.0± 0.3 0.146± 0.246 43.8± 29.6 1.2± 1.1 44.0± 11.8 8 Fen String top 4.6± 0.4 0.095± 0.030 50.9± 15.9 1.5± 0.8 38.9± 10.5 20 String margin 5.7± 0.5 0.094± 0.046 44.7± 21.5 2.5± 1.9 24.3± 13.0 16 Trichophorum tussock 5.9± 0.2 0.133± 0.042 57.2± 11.3 3.4± 0.8 17.1± 1.8 18 Flark 5.8± 0.2 0.098± 0.017 44.2± 8.6 2.7± 0.6 16.5± 2.2 18 https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 566 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape Table A5. Tree biomass and leaf area index derived from field and remote sensing measurements. Ecosystem Tree biomass Tree leaf area index (g of dry matter per metre squared) (m2 m−2) Pine forest 3165 1.97 Treed pine bog 1649 1.02 Sparsely treed pine bog 514 0.32 Table A6. Annual CO2, CH4, and C balances of the five ecosystems and landscape. Uncertainty estimates were not available for the lake balances (Kou et al., 2022), and thus, they are not taken into account in the landscape-scale uncertainty estimates. 11 Jun 2017–10 Jun 2018 11 Jun 2018–10 Jun 2019 CO2 balance (gCm−2) Pine forest −126± 10 −101± 9 Treed pine bog −92± 102 −92± 110 Sparsely treed pine bog −48± 134 −68± 145 Fen −14± 7 9± 5 Lakes 53 56 Landscape −45± 59 −33± 63 CH4 balance (gCm−2) Pine forest −0.2±< 0.1 a −0.2±< 0.1 a Pine bog 1.6± 1.4 b 1.6± 1.4 b Fen 7.8± 0.2 6.3± 0.3 Lakes 3.3 3.7 Landscape 3.0± 0.7 2.7± 0.7 C balance (gCm−2) Pine forest −126± 10 −101± 9 Treed pine bog −90± 102 −90± 110 Sparsely treed pine bog −46± 134 −66± 145 Fen −7± 7 15± 5 Lakes 56 60 Landscape −42± 59 −31± 63 a Dinsmore et al. (2017). b Bubier et al. (2005). Table A7. Data sources for different ecosystems. Ecosystem Eddy covariance Flux chamber Pine bog flux model Arctic Lake Biogeochemistry Estimate from measurements measurements (Sect. 2.4) Model (Kou et al., 2022) the literature Pine forest CO2 CH4 a Fen CO2, CH4 CO2, CH4 Lakes CO2, CH4 CO2, CH4 Pine bog CO2 CH4 b a Dinsmore et al. (2017). b Bubier et al. (2005). Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 567 Figure A1. Temperature response of ecosystem respiration in the pine forest based on the nighttime EC data for 15 June–15 August 2017 and 2018. The lines show the fitted temperature response ER= R10× e E0( 1 56.02− 1 T−227.13 ), with E0= 400 K−1 (Lloyd and Taylor, 1994). The shaded areas indicate the one-sigma confidence intervals. Figure A2. String top gap-filled fluxes. Weekly averaged (a) and cumulative CO2 flux (b) and weekly averaged ER (c) and GPP (d) flux, which is used for estimating pine bog ecosystem fluxes. https://doi.org/10.5194/bg-20-545-2023 Biogeosciences, 20, 545–572, 2023 568 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape Table A8. Average bulk density, C and N content, and C : N ratio (± standard deviation) in the top 20 cm of lake sediments. Lake Bulk density Soil C content of Soil N content of C : N ratio No. of core (gcm−3) dry mass (%) dry mass (%) samples Jänkäjärvi (organic sediment lake) 0.13± 0.03 30.3± 1.6 1.6± 0.1 18.4± 1.0 1 Jänkälampi (mineral sediment lake) Northern basin 0.91± 0.47 5.5± 5.4 0.3± 0.3 19.9± 6.2 4 Southern basin 0.13± 0.10 22.6± 3.1 1.4± 0.4 15.6± 2.7 1 Data availability. The measured flux and ancillary meteo- rological and environmental data are available on Zenodo (https://doi.org/10.5281/zenodo.7437333; Heiskanen et al., 2022). Supplement. The supplement related to this article is available on- line at: https://doi.org/10.5194/bg-20-545-2023-supplement. Author contributions. MA, JPT, and LH designed the study. Ter- restrial ecosystem field flux measurements and maintenance were carried out by LH and MA. Aquatic ecosystem field flux measure- ments were carried out by LH, MA, and SJ. AR, TV, and SJ mea- sured and analysed the vegetation and land cover data. JM measured and analysed the soil data. SJ and LH collected the lake sediment data, and SJ and JM took part in the sediment data analysis. HV was responsible for the eddy covariance data gap-filling and uncertainty analysis. AR and TV conducted the linear regression model analy- sis. The rest of the data analysis was carried out by LH, JPT, and MA. LH wrote the paper, with contributions from all co-authors. Competing interests. The contact author has declared that none of the authors has any competing interests. Disclaimer. Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Acknowledgements. We thank Jani Antila, Holtti Hakonen, and Tuuli Lehtosalo, for field assistance, and Jan Weckström and Maija Heikkilä, for granting ancillary lake sediment chemistry data. We thank Dan Kou, for providing the lake flux modelling results. We thank Tuomas Laurila and Timo Penttilä, for expert support on planning the measurement site set-up. Financial support. This research has been supported by the Academy of Finland (grant nos. 296888 and 296423) as part of CAPTURE-project and the Maj and Tor Nessling foundation (grant no. 202000391). Review statement. This paper was edited by Jamie Shutler and re- viewed by two anonymous referees. References Aguadé, D., Poyatos, R., Gómez, M., Oliva, J., and Martínez- Vilalta, J.: The role of defoliation and root rot pathogen infec- tion in driving the mode of drought-related physiological decline in Scots pine (Pinus sylvestris L.), Tree Physiol., 35, 229–242, https://doi.org/10.1093/treephys/tpv005, 2015. Alm, J., Schulman, L., Walden, J., Nykänen, H., Mar- tikainen, P. J., and Silvola, J.: Carbon Balance of a Bo- real Bog during a Year with an Exceptionally Dry Sum- mer, Ecology, 80, 161–174, https://doi.org/10.1890/0012- 9658(1999)080[0161:CBOABB]2.0.CO;2, 1999. Aubinet, M., Vesala, T., and Papale, D.: Eddy Covariance: A Prac- tical Guide to Measurement and Data Analysis, Springer Nether- lands, Dordrecht, https://doi.org/10.1007/978-94-007-2351-1, 2012. Aurela, M., Laurila, T., and Tuovinen, J.-P.: Seasonal CO2 balances of a subarctic mire, J. Geophys. Res.-Atmos., 106, 1623–1637, https://doi.org/10.1029/2000JD900481, 2001. Aurela, M., Laurila, T., and Tuovinen, J.-P.: Annual CO2 bal- ance of a subarctic fen in northern Europe: Importance of the wintertime efflux, J. Geophys. Res.-Atmos., 107, 4607, https://doi.org/10.1029/2002JD002055, 2002. Aurela, M., Laurila, T., and Tuovinen, J.-P.: The tim- ing of snow melt controls the annual CO2 balance in a subarctic fen, Geophys. Res. Lett., 31, L16119, https://doi.org/10.1029/2004GL020315, 2004. Aurela, M., Riutta, T., Laurila, T., Tuovinen, J.-P., Vesala, T., Tuit- tila, E., Rinne, J., Haapanala, S., and Laine, J.: CO2 exchange of a sedge fen in southern Finland-the impact of a drought period, Tellus B, 59, 826–837, https://doi.org/10.1111/j.1600- 0889.2007.00309.x, 2007. Aurela, M., Lohila, A., Tuovinen, J.-P., Hatakka, J., Penttilä, T., and Laurila, T.: Carbon dioxide and energy flux measurements in four northern-boreal ecosystems at Pallas, Boreal Environ. Res., 20, 455–473, 2015. Barton, K.: MuMIn: Multi-Model Inference, R package version 1.43.17 [code], https://CRAN.R-project.org/package=MuMIn (last access 1 October 2021), 2020. Bastviken, D., Cole, J., Pace, M., and Tranvik, L.: Methane emis- sions from lakes: Dependence of lake characteristics, two re- gional assessments, and a global estimate, Global Biogeochem. Cy., 18, GB4009, https://doi.org/10.1029/2004GB002238, 2004. Biogeosciences, 20, 545–572, 2023 https://doi.org/10.5194/bg-20-545-2023 https://doi.org/10.5281/zenodo.7437333 https://doi.org/10.5194/bg-20-545-2023-supplement https://doi.org/10.1093/treephys/tpv005 https://doi.org/10.1890/0012-9658(1999)080[0161:CBOABB]2.0.CO;2 https://doi.org/10.1890/0012-9658(1999)080[0161:CBOABB]2.0.CO;2 https://doi.org/10.1007/978-94-007-2351-1 https://doi.org/10.1029/2000JD900481 https://doi.org/10.1029/2002JD002055 https://doi.org/10.1029/2004GL020315 https://doi.org/10.1111/j.1600-0889.2007.00309.x https://doi.org/10.1111/j.1600-0889.2007.00309.x https://CRAN.R-project.org/package=MuMIn https://doi.org/10.1029/2004GB002238 L. Heiskanen et al.: Meteorological responses of CO2 and CH4 fluxes in ecosystems of a subarctic landscape 569 Bond-Lamberty, B., Peckham, S. D., Ahl, D. E., and Gower, S. T.: Fire as the dominant driver of central Cana- dian boreal forest carbon balance, Nature, 450, 89–92, https://doi.org/10.1038/nature06272, 2007. Breiman, L.: Random forests, Mach. Learn., 45, 5–32, 2001. Bubier, J. L., Crill, P. M., Moore, T. R., Savage, K., and Varner, R. K.: Seasonal patterns and controls on net ecosystem CO2 ex- change in a boreal peatland complex, Global Biogeochem. Cy., 12, 703–714, https://doi.org/10.1029/98GB02426, 1998. Bubier, J., Moore, T., Savage, K., and Crill, P.: A compari- son of methane flux in a boreal landscape between a dry and a wet year, Global Biogeochem. Cy., 19, GB1023, https://doi.org/10.1029/2004GB002351, 2005. Büker, P., Morrissey, T., Briolat, A., Falk, R., Simpson, D., Tuovi- nen, J.-P., Alonso, R., Barth, S., Baumgarten, M., Grulke, N., Karlsson, P. E., King, J., Lagergren, F., Matyssek, R., Nunn, A., Ogaya, R., Peñuelas, J., Rhea, L., Schaub, M., Uddling, J., Werner, W., and Emberson, L. D.: DO3SE modelling of soil moisture to determine ozone flux to forest trees, Atmos. Chem. Phys., 12, 5537–5562, https://doi.org/10.5194/acp-12- 5537-2012, 2012. Chapin, F. S., Matson, P. A., and Vitousek, P.: Princi- ples of Terrestrial Ecosystem Ecology, 2nd Edn., Springer, https://doi.org/10.1007/978-1-4419-9504-9, 2011. Chen, T. and Guestrin, C.: Xgboost: A scalable tree boosting sys- tem, Proceedings of the 22nd ACM SIGKDD international con- ference on knowledge discovery and data mining, 13–17 August 2016, San Francisco, California, USA, 785–794, 2016. Chi, J., Nilsson, M. B., Laudon, H., Lindroth, A., Wallerman, J., Fransson, J. E. S., Kljun, N., Lundmark, T., Ottosson Löfve- nius, M., and Peichl, M.: The Net Landscape Carbon Bal- ance—Integrating terrestrial and aquatic carbon fluxes in a man- aged boreal forest landscape in Sweden, Glob. Change Biol., 26, 2353–2367, https://doi.org/10.1111/gcb.14983, 2020. Cole, J. J., Prairie, Y. T., Caraco, N. F., McDowell, W. H., Tranvik, L. J., Striegl, R. G., Duarte, C. M., Kortelainen, P., Downing, J. A., Middelburg, J. J., and Melack, J.: Plumbing the global carbon cycle: Integrating inland waters into the terrestrial carbon bud- get, Ecosystems, 10, 171–184, https://doi.org/10.1007/s10021- 006-9013-8, 2007. Denfeld, B. A., Lupon, A., Sponseller, R. A., Laudon, H., and Karlsson, J.: Heterogeneous CO2 and CH4 patterns across space and time in a small boreal lake, Inland Waters, 10, 348–359, https://doi.org/10.1080/20442041.2020.1787765, 2020. Deppe, M., Knorr, K.-H., McKnight, D. M., and Blodau, C.: Effects of short-term drying and irrigation on CO2 and CH4 production and emission from mesocosms of a north- ern bog and an alpine fen, Biogeochemistry, 100, 89–103, https://doi.org/10.1007/s10533-010-9406-9, 2010. Dinsmore, K.