Articles | Volume 17, issue 13
https://doi.org/10.5194/bg-17-3409-2020
https://doi.org/10.5194/bg-17-3409-2020
Research article
 | 
06 Jul 2020
Research article |  | 06 Jul 2020

Carbon dioxide and methane fluxes from different surface types in a created urban wetland

Xuefei Li, Outi Wahlroos, Sami Haapanala, Jukka Pumpanen, Harri Vasander, Anne Ojala, Timo Vesala, and Ivan Mammarella
Abstract

Many wetlands have been drained due to urbanization, agriculture, forestry or other purposes, which has resulted in a loss of their ecosystem services. To protect receiving waters and to achieve services such as flood control and storm water quality mitigation, new wetlands are created in urbanized areas. However, our knowledge of greenhouse gas exchange in newly created wetlands in urban areas is currently limited. In this paper we present measurements carried out at a created urban wetland in Southern Finland in the boreal climate.

We conducted measurements of ecosystem CO2 flux and CH4 flux (FCH4) at the created storm water wetland Gateway in Nummela, Vihti, Southern Finland, using the eddy covariance (EC) technique. The measurements were commenced the fourth year after construction and lasted for 1 full year and two subsequent growing seasons. Besides ecosystem-scale fluxes measured by the EC tower, the diffusive CO2 and CH4 fluxes from the open-water areas (FwCO2 and FwCH4, respectively) were modelled based on measurements of CO2 and CH4 concentration in the water. Fluxes from the vegetated areas were estimated by applying a simple mixing model using the above-mentioned fluxes and the footprint-weighted fractional area. The half-hourly footprint-weighted contribution of diffusive fluxes from open water ranged from 0 % to 25.5 % in 2013.

The annual net ecosystem exchange (NEE) of the studied wetland was 8.0 g C-CO2 m−2 yr−1, with the 95 % confidence interval between −18.9 and 34.9 g C-CO2 m−2 yr−1, and FCH4 was 3.9 g C-CH4 m−2 yr−1, with the 95 % confidence interval between 3.75 and 4.07 g C-CH4 m−2 yr−1. The ecosystem sequestered CO2 during summer months (June–August), while the rest of the year it was a CO2 source. CH4 displayed strong seasonal dynamics, higher in summer and lower in winter, with a sporadic emission episode in the end of May 2013. Both CH4 and CO2 fluxes, especially those obtained from vegetated areas, exhibited strong diurnal cycles during summer with synchronized peaks around noon. The annual FwCO2 was 297.5 g C-CO2 m−2 yr−1 and FwCH4 was 1.73 g C-CH4 m−2 yr−1. The peak diffusive CH4 flux was 137.6 nmol C-CH4 m−2 s−1, which was synchronized with the FCH4.

Overall, during the monitored time period, the established storm water wetland had a climate-warming effect with 0.263 kg CO2-eq m−2 yr−1 of which 89 % was contributed by CH4. The radiative forcing of the open-water areas exceeded that of the vegetation areas (1.194 and 0.111 kg CO2-eq m−2 yr−1, respectively), which implies that, when considering solely the climate impact of a created wetland over a 100-year horizon, it would be more beneficial to design and establish wetlands with large patches of emergent vegetation and to limit the areas of open water to the minimum necessitated by other desired ecosystem services.

1 Introduction

Wetlands provide many beneficial ecosystem services such as flood control and water quality mitigation, natural habitat for flora and fauna, and recreational opportunities (Mitsch and Gosselink, 2015). Many wetlands have been drained globally for agriculture, forestry and other purposes, including urbanization at the cost of losing wetland ecosystem services (Vasander et al., 2003). Migration from rural areas to cities will increase in even greater numbers in the near future, and the United Nations report (United Nations, 2016) has predicted that 75 % of the world population will be living in cities by 2030. There is an urgent need for more sustainable urbanism, and one effective measure is to create functional and connected wetland networks in cities (Lucas et al., 2015; Mungasavalli and Viraraghavan, 2006).

Wetlands can take up carbon dioxide (CO2) through emergent and submerged vegetation, but they are also important sources of methane (CH4), a greenhouse gas more potent than CO2 when considered over a 100-year horizon (IPCC, 2013). The exchange of greenhouse gases (GHG) such as CO2 and CH4 between atmosphere and ecosystem have direct influence on the atmospheric concentration of these gases; thus besides the ecosystem services that wetlands provide, the GHG budget of constructed wetlands should be accounted for according to international agreements such as the Paris Agreement.

Reports on boreal wetlands, such as peatlands, have shown that large carbon storage remains in the soil due to anaerobic conditions limiting microbial decomposition and thus offering a global cooling effect (Frolking et al., 2006). However, in newly constructed urban wetlands on mineral soil the gas exchange may be very different from natural wetlands: (1) the cooling effect of a wetland may be reduced, or it becomes a source of carbon due to the early successional stage of the wetland. When an urban wetland is newly created by wetting a landscape, it takes time for vegetation to establish itself in the new environment. The low coverage of vegetation at the initial phase of wetland establishment can lead to low CO2 sequestration on an ecosystem scale. (2) Wetlands in close proximity to urban centres receive a significant amount of nutrients and dissolved organic carbon from runoff (Lu et al., 2009; Vohla et al., 2007; Valkama et al., 2017). At the areas with emergent vegetation, CO2 is absorbed by photosynthetic activity during daytime and growing season and is released through respiration. At open-water surfaces, the net production of CO2 is a result of photosynthesis by algae, cyanobacteria as well as submerged aquatic plants, and respiration of organic carbon. When the CO2 concentration in the water exceeds atmospheric equilibrium, the surface becomes a source of CO2. CH4 can be produced through anaerobic metabolism in wetland soil and can be transported to the atmosphere by plant-mediated pathways through aerenchyma, sediment ebullition and diffusive fluxes at the water–atmosphere interface. In open water, the transport is dominated by diffusion, whereas in the vegetated areas the plant-mediated transport is most prominent.

Urban wetlands have received extensive attention globally, and their societal and economical importance have been evaluated (Salminen et al., 2013), whereas their climate impact is still largely overlooked except for only a few studies (e.g. Morin et al., 2014a, b). The thus far review of GHG emission in constructed wetlands for wastewater treatment reported that the average CO2 emission was 92.3 mg CO2-C m−2 h−1 and that the CH4 emission ranged from 1.6 to 27 mg CH4-C m−2 h−1 from free water surface (Mander et al., 2014). All of the studies were based on static chamber measurements during a short period so that the annual carbon balance of the ecosystem could not be assessed. In contrast to static chamber measurements, the eddy covariance (EC) method provides continuous measurements of GHG exchange at ecosystem scale, presenting the net result of fluxes as exchange in different source areas contributing simultaneously within the footprint extent (Baldocchi, 2003). It is worth noting that one of the assumptions of the EC method is surface homogeneity, yet in many study sites the landscape setting is far from ideal. The change of source area due to changes in wind provides difficulties in estimating GHG emissions in spatially heterogeneous sites, especially in short-term flux measurements (Baldocchi et al., 2012). Therefore, for heterogeneous sites such as urban wetlands, accurate footprint modelling and a surface area map at high spatial resolution are important in identifying the source area, and a land-surface-specific analysis is vital to reveal the diel pattern, sink/source strength of the wetland.

The objective of this study is to investigate how CO2 and CH4 surface–atmosphere exchange vary with seasonality and spatial heterogeneity and what the annual radiative forcing of these gases is in a created urban wetland at the Nummela suburb, municipality of Vihti, Southern Finland. The studied Gateway wetland was designed and implemented to serve the purposes of storm water quality treatment, creating an urban park, as well as supporting biodiversity. Besides taking advantage of ecosystem-scale EC measurements, we also parse the variability of gas exchange induced by surface heterogeneity (the open-water and the vegetated areas) using diffusional flux modelling and footprint modelling overlapped on a high-resolution surface map. To illustrate how the urban wetland functions as a source or a sink of GHG equivalents, we calculate separately the sustained global warming potential (SGWP) of CO2 and CH4 over a 100-year horizon in each surface type.

2 Materials and methods

2.1 Site description

Our study site is a created storm water wetland Gateway, located by an eutrophicated lake Enäjärvi in the district of Nummela, municipality of Vihti, Southern Finland (60.3272 N, 24.3369 E). Southern Finland experiences a climate with a 30-year mean air temperature of 4.6 C and an annual precipitation rate of 627 mm in the period of year 1981–2010 (Pirinen et al., 2012).

The wetland was constructed in 2010 at the mouth of a 550 ha largely urbanized (35 % impervious) watershed of Kilsoi stream. It was excavated over 6 weeks in early winter 2010 on an abandoned agricultural field growing meadow vegetation. All of the old drainage ditches were blocked as amphibian habitats, which also ensured only one inlet route receiving water from Kilsoi stream and one outlet route discharging water to the nearby lake Enäjärvi. Lake Enäjärvi is an eutrophicated lake. The internal phosphorus load from human activities and the runoff from its catchments have resulted in regular cyanobacterial blooms and fish kills in the lake (Varis et al., 1989; Salonen and Varjo, 2000).

The wetland park has a total area of 7 ha within which – during mean water flow conditions – a 0.5 ha inundated wetland is located. This storm water treatment wetland consists of an inlet stilling pond, a meandering shallow water area with three habitat islands and an outlet pond. The average water depth in the ponds is 1.5 m; within emergent vegetation patches water depth ranges between 0.3 and 0.5 m. There are also submerged macrophytes in the open water as the water is shallow; thus in the paper we refer to the area with emergent vegetation as the “vegetated area” and to the area covered by water in the absence of emergent macrophytes as the “open-water area”. The outlet bottom dam sets a low water level (WL) to 50.04 m above the Baltic Sea level (N60+ coordinate system). Herbaceous vegetation has been allowed to fully self-establish after the construction of the wetland. Annual monitoring of vegetation carried out in summers 2010, 2011 and 2012 indicated rapid self-establishment of vegetation which was rich in taxa and dominated by native species (Wahlroos et al., 2015). At the frequently inundated area (elevation levels of 50–50.35 m), vegetation was arranged in dense patches with different dominating wetland plant species: Typha latifolia L., Iris pseudacorus L., Carex spp. or Juncus effuses L. At the major less frequently inundated area (elevation levels of 50.35–50.45 m), the wet meadow species Filipendula ulmaria L. (Maxim.), Lysimachia vulgaris L. and Lythrum salicaria L. with the three species coexisting at a 1 : 1 : 1 ratio formed the plant community. Drier areas (elevation levels of 50.45–50.60 m) were mostly colonized by dry meadow species such as Poa spp. and Calamagrostis spp., including patches dominated by Cirsium species (Fig. 1). Note that the area with water level lower than 49.5 m is defined as the open-water area while the rest is defined as the vegetated area in this study.

https://www.biogeosciences.net/17/3409/2020/bg-17-3409-2020-f01

Figure 1The landscape classification of the Nummela wetland. Wetland subareas specified according to mean water level are shown with different colours. The arrows indicate the direction of the water flow. The black dots indicate the inflow and outflow measuring stations and the location of the EC tower.

Download

2.2 Water and micrometeorological measurements

Water monitoring stations were set up at the inlet (60.3283 N, 24.3356 E) and at the outlet (60.3281 N, 24.3377 E) of the wetland. During the 2012–2013 and 2013–2014 monitoring periods, water temperature as well as water turbidity, oxygen concentration, conductivity and pH were measured at the inlet and outlet monitoring station with the YSI-6600 series multiparameter sonde (YSI Inc., Yellow Springs, OH, USA). Measurements were conducted continuously with a 10 min interval. Water level at the outlet was measured continuously with a pressure gauge (STS Sensor Technik Sirnach AG, Switzerland). At the outlet monitoring station, the concentrations of dissolved carbon dioxide ([CO2]) and dissolved methane ([CH4]) were measured with Contros HydroC™ CO2 and HydroC™ CH4 sensors (CONTROS Systems and Solutions GmbH, Germany). In 2014, the same sensors were also installed at the inlet monitoring station to measure [CO2] and [CH4]. Dissolved CO2 and CH4 molecules diffuse from the water column into the detection chamber through a thin-film composite membrane where the concentration of CO2 and CH4 is determined by means of IR absorption spectrometry and Tunable Diode Laser Absorption Spectroscopy, respectively. NO3-N and total phosphorus (TP) data have been previously published in Wahlroos et al. (2015). Briefly, NO3-N was measured with Scan sensors (Scan GmbH, Austria) and TP was estimated based on turbidity data measured at 10 min intervals.

Local weather conditions were recorded with a Vaisala WXT weather transmitter (WXT520, Vaisala Oyj, Finland) at the inlet monitoring station. Rainfall, wind speed and direction, temperature, and relative humidity were recorded continuously at a 10 min interval. Photosynthetic photon flux density (PPFD) was measured with a PQS1 PAR quantum sensor (Kipp & Zonen, the Netherlands). Due to instrument failure we obtained PPFD data only from 26 January to 7 April and from 22 July to 29 December 2013. The gaps were filled with PPFD data from another meteorological station nearby (6038 N, 2358 E) in Lettosuo, Finland. The prevailing wind directions were southwest and northeast, and the average of half-hourly average wind speed was 1.13 m s−1 from January to December 2013 with higher wind speed in winter than in summer. The average daily air temperature was 5.9, with the minimum and maximum daily temperatures of −24.4 and 23.3 in 2013. During the winter 2012–2013, there was ice coverage from the beginning of December 2012 to the end of March 2013. In contrast, winter was mild and warm in 2014, and there was practically no snow cover during a winter period (December 2013–March 2014).

2.3 Greenhouse gas measurements by the EC tower technique and gap-filling procedure

To understand the whole-ecosystem exchange of CO2 and CH4 in the wetland, a 2.9 m eddy covariance tower was established in the autumn of 2012 on the southern side of the wetland. The operational period of the EC tower was the entire calendar year of 2013 (from 1 January to 31 December 2013) and the peak growing season in 2014 (from 1 June to 31 August 2014). The EC set-up included a 3D-sonic anemometer (uSonic-3, Metek, Elmshorn, Germany) to measure the three wind speed components and sonic temperature, a gas analyser (LI-7200, LI-COR Inc., Lincoln, Nebraska, USA) which measures CO2 and H2O mixing ratio, and a TDL gas analyser (TGA100A, Campbell Scientific Inc., USA) to measure CH4 mixing ratio. Data from the analysers were collected on a computer at the frequency of 10 Hz. The post-processing of the EC flux data has been done with EddyUH post-processing software (Mammarella et al., 2016). The fluxes were calculated as 30 min covariances between the vertical wind velocity and the gas mixing ratio using block averaging. The raw data were despiked according to standard methods (Vickers and Mahrt, 1997). Coordinate rotations were conducted by performing a two-step rotation to make the x axis along the mean wind direction and the mean vertical wind velocity zero within each 30 min block. The time lag between the anemometer and gas analyser signals, resulting from the transport through the inlet tube, was determined for each 30 min interval by maximizing the cross-correlation function between vertical wind speed and the scalar (CO2 and CH4). The fluxes were corrected for high-frequency loss due to the limited frequency response of the EC system and low-frequency loss due to the limited averaging time period used for calculating the fluxes. Theoretically and experimentally determined co-spectral transfer functions at low and high frequency were used in the correction (Mammarella et al., 2009).

After calculating the fluxes, data collected from periods when the sonic anemometer showed signs of freezing (mean temperature <0.5C and standard deviation of temperature >1.5) were discarded. The data collected during weak turbulence with friction velocity below 0.1 m s−1 have been removed. The measurement points with flux stationarity greater than 1 were omitted to ensure the quality of the co-variances. Fluxes were further filtered according to the wind direction. Since the patchy forest to the southeast of the EC tower (from 100 to 200) and the highway to the west (from 200 to 280) could potentially lead to flow distortion and an additional source of CO2 and CH4, only fluxes from 280 to 100 were accepted for further analysis. The percentage of 30 min fluxes excluded from this analysis was 72 % for CO2 and 73 % for CH4 in 2013, whereas in 2014 the percentage for data exclusion was 54 % for CO2 and 68 % for CH4.

We used an artificial neural network (ANN) technique to gap-fill half-hourly flux data using meteorological variables (Moffat et al., 2007; Papale et al., 2006). Those variables included radiation, air temperature, water temperature, water level, wind speed, relative humidity, time of the day, season, and dissolved CO2 and CH4 concentration in the water. We tested the model performance with different ANN architectures, starting from the architecture with the most complexity, and then reduced the variables to find the simplest ANN architecture with good performance (more than 5 % loss in model accuracy with additional variable reduction). For CO2, water level and wind speed were found to have trivial contributions to the ANN model, and thus they were removed from the model input, while for CH4, only wind speed was removed for the same reason. We found that dissolved gas concentrations greatly improved the model prediction as they captured the variation of diffusive fluxes from the water (Fig. S1 in the Supplement). Ancillary meteorological variables in general had good data coverage, and short gaps (up to several hours) were gap-filled by linear interpolation. The only exception was dissolved gas concentration, which had a long measurement breakage in 2013 (day of year 214–254). Fluxes were therefore gap-filled with two separate ANNs, one with dissolved gas concentration and one without. During the above-mentioned period with long gaps, the ANNs modelled without dissolved gas concentration were used to gap-fill the half-hourly flux data.

The Levenberg–Marquardt algorithm was used in the learning process of ANN. The optimized number of neurons in the hidden layer was determined by training the network 100 times with varying number of neurons (from 3 to 15), and 10 neurons were considered to be sufficient after evaluating the performance of the network using root-mean-square error (RMSE) (data not shown). The entire dataset was divided into three parts: two-thirds of the data were used to train the networks, one-sixth were used for testing the networks and the remaining one-sixth were used for validating the networks. Since the training of the networks can be biased towards periods with greater data coverage (e.g. daytime conditions), the environmental variables were first divided into five natural clusters using a k-mean clustering algorithm in MATLAB (MathWorks, Inc., Natick, MA, USA), and then the data used for training, testing and validation were proportionally extracted from each cluster. After each data extraction, the network was reinitialized for 10 times to avoid local minima, the initialization with the lowest RMSE was selected and the resulting network was saved. We repeated the whole process of data extraction and initialization for 20 times, and we used the median of these 20 predictions to gap-fill the missing flux values. The uncertainty of the ANN gap-filling procedure was presented using a 95 % confidence interval of the 20 ANN predictions.

In order to be confident in our gap-filling results, we also applied alternative gap-filling methods to EC fluxes using parameterization based on biological principles (see Supplement). The results on annual cumulative fluxes were not significantly different from the ones gap-filled using ANN; thus we only report the results from ANN in the following text.

2.4 Diffusive gas exchange

We calculated diffusive gas exchange F from open water according to the boundary layer model:

(1) F = k c aq - c eq ,

where k is the gas transfer velocity (cm h−1), caq is the gas concentration in surface water (mol m−3) and ceq is the gas concentration that surface water would have when it reaches equilibrium with the air (mol m−3). caq and ceq can be obtained according to the solubility of the gas:

(2)caq=10-3kHpχwater,(3)ceq=10-3kHpχair,

where kH is Henry's law constant for the respective gas (mol L−1 atm−1), p is air pressure (atm), χwater is the gas mixing ratio in surface water (ppm) and χair is the gas mixing ratio in the air (ppm). In this study, χwater was obtained from the outlet monitoring station as it was located most of the time in the flux footprint area and it had longer data coverage than from the inlet monitoring station. The gas transfer velocity k can be calculated as the formula below (Cole and Caraco, 1998):

(4) k = 2.07 + 0.215 U 10 1.7 Sc 600 - 0.5 ,

where U10 is the horizontal wind speed extrapolated to 10 m using the theoretical log wind profile equation (m s−1, approximately U10=1.15U, where U is the measured wind speed at 2.9 m height in the study site) and Sc is the temperature-dependent Schmidt number of the respective gas. When a gas concentration measurement was not available, linear interpolation was applied to obtain monthly and annual diffusive GHG fluxes from the open water.

Although the above-mentioned Cole–Caraco (CC) method is the most simple and most often used model for gas transfer velocity, the limitation of the CC method is that it considers wind as the sole factor to cause the water turbulence and to drive the gas exchange. More complicated models were suggested to include the effect of buoyancy-flux-driven turbulence (Heiskanen et al., 2014; Tedford et al., 2014), which could not be applied in the current study as we did not have the measurements of some of the required model inputs, e.g. the net shortwave and longwave radiation from the water body. However, we believe that the CC method is suitable to be applied here. First, the water body in our study is located in an open area where the contribution of wind shear to the turbulence in the surface mixed layer is relatively high. During the study period in 2013, the average wind speed was 1.57 m s−1 with a maximum of 7.1 m s−1. Secondly, the k600 estimated in our study was on average 0.66 m d−1, well situated within the range of the k600 which was directly measured by a floating chamber or gas tracer for small lakes and ponds (Holgerson et al., 2017). Thirdly, the estimated air–water fluxes of CH4 and CO2 based on the CC model were also well within the range of the diffusive gas fluxes over small lakes from other studies (Erkkilä et al., 2018; Mammarella et al., 2015). Finally, the parameterization of Cole and Caraco has been similarly applied to connected small open-water pools, and good agreement between the model estimation and the measurements was found (Cole et al., 2010; McNicol et al., 2017). Therefore, we decided to use the CC model to estimate diffusive fluxes from the water, bearing in mind that the calculated fluxes could potentially be underestimated.

2.5 Estimating zone fluxes and radiative forcing

By combining the EC tower and diffusive flux from the open water, the following model can be derived:

(5) F EC = F water × f water + F veg × f veg ,

where FEC is the flux measured by the EC tower, and Fwater and Fveg stand for the fluxes from open-water and vegetated areas, respectively. fwater and fveg are the footprint-weighted spatial fraction of open-water and vegetated areas. In this study, ebullition was neither measured nor calculated, so the flux from water was only represented by the diffusive flux.

Specifically, we first modelled the half-hourly flux footprint with a parameterization of a three-dimensional backward Lagrangian footprint model (Kljun et al., 2015) in MATLAB. Periods in which the wind came from the patchy forest to the southeast of the EC tower (between 100 and 200) and the highway to the west (between 200 and 280) were eliminated in the footprint analysis. Secondly, a land cover classification map of vegetated and open-water zones was delineated manually using a high-resolution aerial image (data from the National Land Survey of Finland Topographic Database 06/2013) with an image manipulation software (Gimp 2.10.6, https://www.gimp.org, last access: 22 November 2018). Thirdly, the flux footprints were aligned and combined with the land cover classification map to calculate half-hourly fveg and fwater within 90 % footprint contour lines. Specifically, we assigned each footprint pixel within the 90 % footprint area to either open-water area or vegetated area on the land cover classification map while the footprint of the pixels outside the 90 % footprint area was regarded as zero. fwater was calculated as the sum of the footprint within the open-water area to the total footprints, while fveg was calculated as the sum of footprint within the vegetated area to the total footprints. None of the 90 % footprint contour lines exceeded the map area, and the sum of fveg and fwater equaled to 1. In order to obtain the long-term aggregated footprint of carbon fluxes, we calculated also the monthly and annual aggregated footprint climatology during the study period.

The uncertainty of the vegetation and water fraction comes from two sources. Firstly, the delineation of the distinct surface types was conducted based on a land surface map of the growing season in 2013, which neglected the change in the spatial extent of the vegetation throughout the presented GHG monitoring time. Secondly, although the footprint model used here is proven to be robust and general, there are uncertainties in the model predictions. To be more confident in the footprint estimation, it would be good to compare our results with footprint estimates based on large eddy simulations; however, that was out of the scope of the current study. With only one EC tower we could not cross-check the results as was done in another study (Matthes et al., 2014). However, we chose to follow a simple approach dividing the landscape into vegetation and open water because we did not observe significant vegetation expansion to cover more of the open water during the growing season in 2013, and thus the area of open water was relatively constant during the monitored time periods. Vegetation establishment at the Gateway wetland was very rapid in the first growing season, and by the summer 2012 emergent vegetation had densely established at the intended shallow wetland areas. Furthermore, the clear effect of the footprint-weighted fraction of open water on the synchronization between EC CH4 measurements and diffusive CH4 flux from water was nicely demonstrated in our analysis (Fig. S3).

To better understand the influence of greenhouse gas fluxes in this urban wetland, we calculated the sustained global warming potential (SGWP) for CO2 and CH4 over a 100-year horizon in each surface type. The difference between SGWP and global warming potential (GWP) is that SGWP accounts for the effect of GHG remaining in the atmosphere during the period. Since CH4 is a more potent greenhouse gas, we multiply the emission of CH4 by a factor of 45 to convert it to kg CO2-eq m−2 yr−1 (Neubauer and Megonigal, 2015). However, for an easy comparison between our results and those from other studies using the conventional method, we calculated also CH4 fluxes as CO2 equivalents using a GWP of 34 following the Fifth Assessment Report of IPCC (2013).

2.6 Statistical analysis

Differences in the fluxes and environmental variables between the two peak growing seasons (summer 2013 and 2014) were evaluated using the t test. Cumulative annual GHG fluxes measured by the EC tower are reported as the median of the 20 ANN predictions, and uncertainty is presented as the 95 % confidence interval of the 20 ANN predictions. As diffusive GHG fluxes were calculated from gas concentration meteorological parameters, no standard error is reported for the cumulative annual fluxes from the open water. All statistical analyses were performed in MATLAB.

We also conducted wavelet coherence analysis to explore the temporal correlations between fluxes and environmental variables on the multi-temporal scales (Grinsted et al., 2004; Torrence and Webster, 1998). Since the fluxes are gap-filled using some of the environmental variables, simply applying the wavelet coherence analysis to all the variables can overstate the correlations. Therefore, we only conducted wavelet coherence analysis between gap-filled ecosystem flux time series and those independent environmental variables which were not used in the gap-filling procedure (concentration of NO3-N and TP), while Pearson correlations (r) were determined between non-gap-filled fluxes and the other environmental variables.

3 Results

3.1 Ecosystem seasonality and environmental variables

Daily average PPFD ranged from 0.9 to 691.5 µmol m−2 s−1 in 2013, with the highest value appearing in July. June had the highest monthly average PPFD with 486.1 µmol m−2 s−1, followed by July and August with 470.2 and 430.6 µmol m−2 s−1, respectively. The PPFD during the peak growing season in 2014 was on average 361.8 µmol m−2 s−1, which is lower than that during the same period in 2013 (Fig. 2a).

https://www.biogeosciences.net/17/3409/2020/bg-17-3409-2020-f02

Figure 2The daily average of (a) photosynthetic photon flux density (PPFD), (b) air and water temperature (Tair and Twater), (c) water level, (d) rainfall, (e) CO2 concentration ([CO2]), (f) CH4 concentration ([CH4]), (g) concentration of total phosphorus (TP) and (h) concentration of NO3-N of the Nummela wetland from January 2013 to August 2014.

Download

Mean daily water temperature (Twater) ranged from 0 C in March to 23.7 C in June with an annual average of 7.9 C in 2013 and from 0 C in February to 21.4 C in July in 2014. Mean daily air temperature (Tair) had more fluctuation and ranged from −15.6C in January to 23.3 C in June 2013 and from −19.0C in January to 23.4 C in July 2014 (Fig. 2b). The open-water area experienced an ice-covered period between 1 January and 31 March 2013, while the winter 2013–2014 was so mild and warm that there was practically no snow cover during December 2013–March 2014. Comparing the temperature between the two peak growing seasons, both Twater and Tair were higher in June 2013 while Tair was lower in July 2013 than in 2014. In August, there was no significant temperature difference between the 2 years. Four seasons were classified for the ecosystem based on the trend in Tair and Twater. In spring (April and May), the daily temperature started to increase, the vegetation showed a sign of early growing season and the warm temperature unfroze the lake. In summer, the peak growing season (June–August), vegetation exhibited the maximum growth, and the temperatures reached the annual maxima. In autumn (September and October), daily temperatures began to drop and the vegetation showed signs of early senescence. In winter (January to March and November and December), temperatures reached the annual minimum and vegetation was inactive in carbon sequestration. Precipitation was higher in August 2014 than in the preceding August, almost twice as high as that of 2013.

WL was higher in the winter and lower in the summer in 2013. The daily average of WL varied between 50.06 m in July 2013 and 50.4 m in April 2013. There was a spring peak in 2013 when the highest WL was observed due to snowmelt, while in 2014 no such event appeared due to the mild winter 2013–2014 without an ice-covered period (Fig. 2c). The average daily WL from January to August was similar (50.13 and 50.15 cm for 2013 and 2014, respectively). However, during the peak growing season, it was on average 5.7 cm higher in 2014 than in 2013.

The annual rainfall in 2013 (snowfall not included) was 363.6 mm, which happened mostly during summer and autumn (Fig. 2d). The maximum daily-averaged rainfall was in August (26.7 mm d−1), while monthly-averaged rainfall was highest in November with 73.8 mm per month followed by August with 68.3 mm per month. In 2014, an exceptionally high amount of rainfall was observed in August (125.7 mm per month), while the amount of rainfall in the other months was similar to 2013.

The daily-averaged CO2 concentration in the water ([CO2]) in 2013 had large variation with the maximum (461 µmol L−1) and the minimum (21.6 µmol L−1) (Fig. 2e). [CO2] was higher in 2014 with an average of 262.6 µmol L−1 from January to August than in 2013 with an average of 211.5 µmol L−1. It also exhibited seasonal variation with high concentration in summer (360.3 µmol L−1) and low concentration in winter (223.4 µmol L−1). The [CO2] measured in the inflow was generally lower than that in the outflow, and they were well correlated (r=0.84). [CH4] in the outflow was on average 5 times higher in 2014 than in 2013. The average annual concentration was 0.81 µmol L−1 in 2013 and 2.25 µmol L−1 in 2014. There were peak [CH4] episodes in the outflow in May 2013 with a maximum of 5.43 µmol L−1. During the summer months in 2014 there were even higher outflow [CH4] peaks, with a maximum of 16.83 µmol L−1. The [CH4] had a mean of 0.42 µmol L−1 in the inflow which was lower than that in the outflow, and there were no prominent [CH4] peaks observed in the inflow. [CH4] in the inflow and outflow were weakly correlated (r=0.2) (Fig. 2f).

The median concentration of total phosphorus (TP) measured at the outflow monitoring station was 56 µg L−1, and the median NO3-N concentration was 0.69 mg L−1 in 2013 (Fig. 2g, h). In the annual perspective, TP and NO3-N concentration consisted of several runoff peaks occurring after rain or snow melting events. This wetland serves as a nutrient removal measure as it improved water quality by retaining P and N from runoff before the release to the receiving lake, where the annual TP reduction was 10 % and NO3-N reduction was 3 % from the original concentration in 2013 (Wahlroos, 2019).

3.2 Flux footprint mapping

A footprint distribution was modelled for each half hour when an eddy flux measurement was collected at the EC tower. The open-water area accounted for 10 % to 16 % of the total wetland area within the footprint, while the rest was comprised of wetland vegetation. When weighted with footprint distribution, fwater ranged from 0 % to 25.5 % and fveg from 74.5 % to 100 %. The first quartile, median and third quartile of fwater were 0.09 %, 14.1 % and 17.9 %, respectively, and those of fveg were 82.1 %, 85.9 % and 91.3 %, respectively.

The monthly cumulative footprint was slightly different for CO2 and for CH4 due to the different missing flux values. However, the difference on average was so small (7 %) and the footprint of CO2 was used in further analysis. The flux footprints were shown to be northeast to the EC mast due to the wind direction filtering, meaning only half-hourly data with wind directions from the wetland area were considered in the analysis (Fig. S2). The monthly average of the 90 % footprint area covered a minimum of 0.69 ha to a maximum of 2.28 ha with a mean of 1.3 ha. The mean extent of the 90 % flux footprints was 128 m. After applying the flux footprint function, the monthly average of the footprint-weighted spatial fraction of open water showed lower value in summer and higher value in winter, ranging from 11.3 % to 21.4 % with a mean of 13.3 % in 2013. In 2014, during the peak growing season, on average 13.8 % of the wetland area was open water and the mean fwater was 10 %.

3.3CO2 and CH4 fluxes

3.3.1 Ecosystem CO2 and CH4 fluxes

Ecosystem CO2 and CH4 fluxes measured by the EC tower showed the ecosystem was nearly CO2 neutral and it was a small CH4 source in 2013.

Daily average of NEE was near zero during winter time (January to March, on average 0.37 µmol C-CO2 m−2 s−1), it was slightly positive in spring and it became negative from the end of May until the end of August, indicating the ecosystem was a CO2 sink during this period, with a maximum negative value of −5.14µmol C-CO2 m−2 s−1 in June. Daily average NEE was highest in September with a maximum of 3.29 µmol C-CO2 m−2 s−1. In October, November and December, NEE remained low but still positive (on average 0.77 g µmol C-CO2 m−2 s−1), demonstrating the milder winter between 2013 and 2014 (Fig. 3). NEE exhibited strong seasonality in 2013, which was negative during June, July and August, meaning the ecosystem was a CO2 sink while the rest of year it was a CO2 source. NEE was lowest in June and highest in September. The cumulative NEE in 2013 was 8 g C-CO2 m−2 yr−1, with the 95 % confidence interval between −18.9 and 34.9 g C-CO2 m−2 yr−1 (Fig. 3).

https://www.biogeosciences.net/17/3409/2020/bg-17-3409-2020-f03

Figure 3Daily average of net ecosystem exchange of CO2 (NEE, µmol m−2 s−1) in 2013. Filled dots indicate measurement (when available half-hourly measurement data ≥10) and circles indicate gap-filled data (when available half-hourly measurement data <10). The insert shows cumulative NEE (g C m−2) in the ecosystem and the red line indicates the zero reference line. Error bounds (marked in grey) on cumulative NEE reflect the 95 % confidence interval for the gap-filling procedure.

Download

Daily-averaged CH4 was low but not negligible from January to April (on average 5.1 nmol C-CH4 m−2 s−1), with a sudden rise at the end of May reaching a maximum of 48.9 nmol C-CH4 m−2 s−1. During summer months the ecosystem exhibited relatively high CH4 emission (on average 15.4 nmol C-CH4 m−2 s−1), which is not comparable with the emission episode in May but higher than winter months. In autumn (September and October) the daily-average CH4 was 8.8 nmol C-CH4 m−2 s−1, and after that it gradually decreased throughout the rest of the year with an average of 5.5 nmol C-CH4 m−2 s−1. The cumulative CH4 for 2013 was 3.9 g C-CH4 m−2 yr−1, with the 95 % confidence interval between 3.75 and 4.07 g C-CH4 m−2 yr−1 (Fig. 4).

https://www.biogeosciences.net/17/3409/2020/bg-17-3409-2020-f04

Figure 4Daily average of ecosystem CH4 flux measured by the EC tower and cumulative CH4 emission in 2013. Filled dots indicate measurement (when available half-hourly measurement data ≥10) and circles indicate gap-filled data (when available half-hourly measurement data <10). The insert shows cumulative CH4 emission, with the error bounds in grey reflecting the 95 % confidence interval for the gap-filling procedure.

Download

Comparing the peak growing season between 2013 and 2014, the 30 min NEE ranged from −20.0µmol C-CO2 m−2 s−1 in June to 18.5 µmol C-CO2 m−2 s−1 in September 2013. During the peak growing season 2014, NEE had the lowest value of −22.6µmol C-CO2 m−2 s−1 in June. The monthly NEEs of the peak growing season were −84.1, −76.1 and −22.2 g C-CO2 m−2 per month in June, July and August 2013, and they were −97.6, −47.5 and −19.6 g C-CO2 m−2 per month in 2014. The average CH4 emissions in June, July and August were 24.4, 10.8 and 11 nmol m−2 s−1 in 2013, and they were 15.5, 21.3 and 21.3 nmol m−2 s−1 in 2014, respectively.

3.3.2 Diffusive CO2 and CH4 fluxes from the open-water area

Diffusive CO2 and CH4 fluxes from the open water were estimated based on wind speed, [CO2] and [CH4] (see Sect. 2.4). The variation of diffusive fluxes demonstrated a pattern driven by both wind speed in the short term and gas concentration dynamics in the water in the long term. Diffusive CO2 fluxes ranged from −0.07 to 4.09 µmol CO2 m−2 s−1 with a mean of 1.04 µmol CO2 m−2 s−1 in 2013, indicating CO2 oversaturation in the water. From June to September the averaged flux (1.27 µmol CO2 m−2 s−1) was higher than that of the other months (Fig. 5a), corresponding to the higher [CO2] in the water during summer months (Fig. 2e). The monthly-averaged diffusive CO2 flux during peak growing season in 2014 was 2.34, 2.71 and 1.99 µmol CO2 m−2 s−1 for June, July and August, which is significantly higher than during the same period in 2013 due to the high [CO2] in the open water (Fig. 2e).

https://www.biogeosciences.net/17/3409/2020/bg-17-3409-2020-f05

Figure 5Monthly average of (a) diffusive CO2 and (b) CH4 flux from the open water in 2013. Error bar indicates the standard error of the mean. From January to March there was an ice-covered period.

Download

The average diffusive CH4 emissions in 2013 was 4.9 nmol C-CH4 m−2 s−1, where a peak emission appeared in late May with the highest flux of 137.6 nmol C-CH4 m−2 s−1. Monthly-averaged CH4 diffusive fluxes showed an increasing trend towards the end of the year with large variation in May due to the peak concentration episode. This phenomenon was mainly driven by the increasing dissolved CH4 concentration in the outflow in 2013. The monthly-averaged diffusive CH4 flux during peak growing season in 2014 was 20.9, 18.9 and 13.5 nmol CH4 m−2 s−1 for June, July and August, respectively, and they were significantly higher than the same period in 2013 due to the high [CH4] in the open water (Fig. 2f).

3.3.3 Diel patterns in CO2 and CH4 fluxes

Only non-gap-filled data were used for determination of diel patterns in both gas fluxes. CO2 and CH4 fluxes from the vegetated areas (Fveg) were calculated for each 30 min interval according to Eq. (5). As expected, CO2 flux showed a strong diel pattern in summer with CO2 uptake during daytime and release in the night, which was controlled by photosynthetic activity (Fig. 6a). The summer peak CO2 uptake reached 11.5 µmol m−2 s−1 for the whole constructed wetland ecosystem and 15.2 µmol m−2 s−1 for the vegetated areas. The CO2 flux from the vegetated areas had higher maximum uptake than the EC measurements carried out over the whole constructed wetland. In the winter, the CO2 fluxes from both tower and vegetation were similar, being on average 0.46 and 0.55 µmol m−2 s−1 respectively (Fig. 6b).

https://www.biogeosciences.net/17/3409/2020/bg-17-3409-2020-f06

Figure 6Mean diel pattern of the half-hourly net CO2 and CH4 fluxes in summer (a, c) and in winter (b, d). The dashed lines represent the standard deviation. Red lines indicate measurement from the EC tower and the blue lines show the fluxes modelled for the vegetated area.

Download

CH4 flux also showed diel patterns in the summer with much larger variability than those from CO2 flux. CH4 emission in general was higher in daytime than in nighttime. In the daytime in summer, CH4 flux from the vegetated area was higher than the flux measured from the tower while there was no difference during the nighttime (Fig. 6c, d). The summer peak daytime flux from the tower (18.9 nmol m−2 s−1) and the vegetated area (24.7 nmol m−2 s−1) was 2.4 times and 3.3 times higher than the nighttime flux (7.5 nmol m−2 s−1), respectively. This can be understood as daytime CH4 flux being linked with photosynthesis while nighttime CH4 flux is controlled by other processes like diffusion, ebullition and convection between the soil, water and atmosphere. In winter there was a small (on average 4.6 nmol m−2 s−1) but constantly positive CH4 flux without an obvious diel pattern.

3.4 Environmental variables with fluxes

Only non-gap-filled flux data were used in the Pearson correlation analysis between environmental variables and flux pairs. Radiation, Tair and Twater all had a high negative correlation coefficient (r) with NEE and a high positive r with CH4 flux in 2013, corresponding to the results of ANN model parameter selection. Radiation was best correlated with NEE and Twater was best correlated with CH4 (Table 1). The correlations were rather weak (small r or even the opposite sign of r) during 2014 due to the short measuring period and narrow ranges of the variables. Water level was positively correlated with NEE and negatively correlated with CH4, which was counter-intuitive, possibly because it was masked by temperature variation as the water level was in general higher in winter and lower in summer. [CO2] and [CH4] were not correlated with either NEE or CH4, although they were shown to be important parameters in ANN model selection.

NO3-N did not show consistent correlation with any of the fluxes (Fig. 7a, b). The variation of TP was negatively leading the change in NEE at 1 d scale (more TP leads to more CO2 uptake; Fig. 7c) where the time lag varies between 1 and 5 h (data not shown). TP had a positive correlation with CH4 flux (more TP leads to more CH4 emission) at 1 d scale (Fig. 7d) and TP led the CH4 flux by ∼2 h (data not shown).

Table 1Pearson correlation coefficient (r) between the daily averages of environmental variables and fluxes in 2013 and 2014. NEE – net ecosystem exchange; Tair – air temperature; Twater – water temperature; PPFD – photosynthetic photon flux density; WL – water level; [CO2] and [CH4] – CO2 and CH4 concentration measured in the outlet; * indicates only peak growing seasons (June, July and August) are included in the analysis.

Download Print Version | Download XLSX

https://www.biogeosciences.net/17/3409/2020/bg-17-3409-2020-f07

Figure 7Wavelet coherence analysis and the phase difference between ecosystem fluxes, the net ecosystem exchange (NEE) and the CH4 flux (FCH4), and nutrient concentration in the water, NO3-N and total phosphorus from January to December 2013. The colour represents the power of the coherence from 0 to 1. The phase difference is indicated by black arrows which only show up where the coherence is greater than or equal to 0.5. Rightwards arrow () indicates in-phase (two time series in synchrony) and arrows in other direction indicate out of phase (representing lags between time series); i.e. leftwards arrow () indicates anti-phase, downwards arrow () indicates the first series (fluxes) leads by quarter cycle and upwards arrow () indicates the second series (NO3-N and total phosphorus) leads by quarter cycle. White dash contour lines indicate the cone of influence.

Download

3.5 Estimating radiative forcing from different zones

To obtain the climate forcings from each land surface type, we calculated the annual cumulative fluxes from the vegetated area based on Eq. (5) using footprint-weighted spatial fraction, ecosystem fluxes and diffusive fluxes from the open water (see Sect. 2.5). The annual median value of footprint-weighted spatial extent was used to calculate the annual fluxes, which showed the open-water area was a CO2 source (297.5 g C-CO2 m−2 yr−1) and the vegetated area was a CO2 sink (−39.5 C-CO2 m−2 yr−1). Both the open-water area and the vegetated area were CH4 sources, but the CH4 emission from the vegetated area was higher than the open-water area, being 4.26 and 1.73 g C-CO2 m−2 yr−1, respectively (Table 2).

Table 2Annual CO2 and CH4 exchange, their sustained global warming potential (SGWP) and global warming potential (GWP) from different surface zones in the Nummela wetland in 2013. “Ecosystem”, “water” and “vegetation” represent flux, SGWP and GWP measured or calculated from the ecosystem by the EC tower, from the open-water and vegetated areas, respectively. The numbers in the square bracket represent the 95 % confidence interval of the average. No error bounds are reported for flux, SGWP and GWP from open water as they are modelled using gas concentration in the water and meteorological measurements.

Download Print Version | Download XLSX

Open water has contributed a large amount of CO2 emission into the atmosphere through diffusion (1.09 kg CO2-eq m−2 yr−1), whereas the CH4 emission was relatively small (0.104 kg CO2-eq m−2 yr−1). The vegetated area was a small sink of CO2 but the cooling effect of vegetation by CO2 uptake was relatively small (−0.145 kg CO2-eq m−2 yr−1) compared to its CH4 emission (0.256 kg CO2-eq m−2 yr−1). Overall, the ecosystem had a small warming effect of 0.263 kg CO2-eq m−2 yr−1, of which 89 % was contributed by CH4 (Table 2).

4 Discussion

4.1 The GHG fluxes from an urban storm water wetland ecosystem

The studied urban wetland ecosystem was a small carbon source over the full-year studied period in 2013. Due to the scarcity of studies on urban wetlands using the EC method, we compare our results to restored wetlands which can be considered to be proxy ecosystems to urban wetlands, with both including wetting practice in an ecosystem which has been drained or dry previously. The annual CO2 balance of 8 g C-CO2 m−2 yr−1 from the ecosystem, or −39.5 g C-CO2 m−2 yr−1 from the vegetated area (Table 2), was small compared to a restored wetland in western Denmark where the annual CO2 balance ranged from −286 to −53 g C-CO2 m−2 yr−1 (Herbst et al., 2013), and the annual CH4 balance of 3.9 g C-CH4 m−2 yr−1 was less than half of the annual CH4 emission (between 9 and 13 g C-CH4 m−2 yr−1) in that study. Over a network of restored freshwater wetlands in California, the CO2 sequestration can be up to nearly 700 g C m−2 yr−1 and CH4 emission up to 63 g C m−2 yr−1 (Hemes et al., 2018). It is not surprising that the studied ecosystem appeared to be CO2 neutral as it was recently constructed. The herbaceous vegetation has been allowed to fully self-establish without human intervention, and at the early successional stage, plant diversity and biomass were still increasing each year (Wahlroos, 2019). With the vegetation being more developed, a greater CO2 uptake from the vegetated area can be expected in the following years. The low CH4 emission observed in this study may be due to the depletion of organic matter in the bottom soil from agricultural uses; thus it provided little substrate for anaerobic microbial activity to produce CH4. With the accumulation of organic matter in the anoxic wetland sediment, CH4 production may increase in the future. Certain chemical compounds like Fe in mineral soils can also inhibit CH4 production, leading to much-lower ecosystem-scale CH4 flux (Chamberlain et al., 2018). In the meantime, methane-oxidizing bacteria (methanotroph) regulate CH4 consumption at the soil–water interface. With the ecosystem being used previously as cropland, the physical disturbance of soil may have greatly reduced the methanotroph communities so that the CH4 oxidation may also be low in the soil (Smith et al., 2000; Saggar et al., 2008). Furthermore, after the initial establishing phase, the ecosystem productivity can also be reduced at a water treatment wetland due to the standing litter that inhibits the generation of new vegetation growth. It was shown that in a restored freshwater wetland the ecosystem was a net CO2 sink (-804±131 g C-CO2 m−2 yr−1) in 2002–2003, 6 years after the restoration, but nearly CO2 neutral in 2010–2011 due to the reduced photosynthetic plants (Anderson et al., 2016). Thus, given the urban wetland is sustained for a sufficiently long period, it is still unclear whether the CO2 uptake from the vegetated zone would compensate for its CH4 emission, not considering the large GHG emission from the open-water zone. Thus, similar studies as the present one should be conducted at a later stage after the construction of the wetland to fully reveal the GHG balance of the ecosystem over time.

Overall, the ecosystem CO2 and CH4 fluxes measured by the EC tower ranged from −5.33 to 3.4 g C-CO2 m−2 d−1 and from 1.0 to 55.2 mg C-CH4 m−2 d−1, respectively. They are consistent with the flux ranges provided by other studies on GHG fluxes in restored wetlands (Anderson et al., 2016; Knox et al., 2015; Matthes et al., 2014; Morin et al., 2014b; Herbst et al., 2013), although for both gases they tend to be on the lower end. NEE exhibited seasonal variation so that the ecosystem was a CO2 sink between June and August. The highest NEE appeared in September possibly because photosynthesis greatly reduced due to plant senescence while ecosystem respiration remained relatively high because of the warm temperature. Previous studies have found good agreement between CH4 emission and photosynthesis as plants provide substrates for methanogenesis (Rinne et al., 2018), which was not observed in the daily average of gas fluxes in this study (Figs. 2 and 3) as the peak CH4 flux appeared in May and peak gross primary productivity appeared in July (data now shown). Nonetheless, both CH4 and CO2 fluxes, especially those obtained from the vegetated areas, exhibited strong diurnal cycles during summer with synchronized peaks around noon (Fig. 6a, c). This finding reflects that short-term CH4 emission from vegetation is linked with photosynthesis by providing labile carbon from root exudate and by gas transport through aerenchyma and open stomata, while long-term CH4 emission may be determined by complex processes related to environmental variables, e.g. temperature and redox potential (Linden et al., 2014).

4.2 Parsing GHG fluxes from heterogeneous land surfaces

We found that the open-water area was constantly a source of CO2 and CH4 to the atmosphere during the studied period as the [CO2] and [CH4] in the water generally exceeded the atmosphere equilibrium except during the ice-covered period (Fig. 5). The annual average of [CO2] in the surface water in 2013 was 0.3 % in our study, comparable to 0.4 % in another temperate restored wetland (McNicol et al., 2017), while the seasonal pattern (higher in summer and fall) was the opposite of what McNicol et al. (2017) found. We also found that both [CO2] and [CH4] were higher in 2014 than 2013 (Figs. 2d, 1e). The O2 concentration ([O2]) and O2 balance ([O2]outlet [O2]inlet) measured by another study on the same wetland (Wahlroos, 2019) could partially explain the observed phenomenon. The relatively high water temperature and oxic conditions in the water in fall 2013 have allowed high decomposition of detritus leading to high [CO2] (Wahlroos, 2019). The long period of hypoxia during summer 2014 could explain the 3-fold increase in [CH4] as the condition was more favourable for CH4 production. The negative O2 balance in summer 2014 indicated strong O2 consumption by microbial decomposition producing CO2 in the water. As the long-term diffusive fluxes (daily and monthly) were mainly driven by gas concentration in the water, it was straightforward to understand high diffusive CO2 and CH4 fluxes in 2014 compared to 2013. Interestingly, the ecosystem CH4 emission in 2013 was well synchronized with the diffusive CH4 flux by capturing sporadic emission episodes from the water (Fig. S3a, c), while they were not synchronized in summer 2014 although several stronger diffusive peaks happened (Fig. S3b, d). When footprint-weighted contribution was accounted for, it clearly revealed that the synchronization of CH4 emission from ecosystem and water was closely related to the flux footprint distribution. When there was high flux contribution from the open water (20 %–25 %), high diffusive CH4 was also reflected in ecosystem flux measured by EC. This has further proved the application of footprint analysis is essential in explaining gas exchange from heterogeneous surfaces using EC data.

It is worth noting that in our study we only classified the surface landscapes into “open water” and “vegetation” but neglected the difference in sink/source strength from different plant types within the vegetation zone (Fig. 1). We did not account for the dissimilarity between vegetation types because the characteristics in gas exchange are much more distinct between open water and vegetation, which was the focus of this study. For the same reason, ebullition was not considered in this study either, as ebullition was shown to have only minor significance in a restored wetland, accounting for less than 0.1 % of ecosystem CO2 flux and 4.1 % of ecosystem CH4 flux (McNicol et al., 2017). However, for a proper downscaling analysis of EC data, the subareas of different plant types and ebullition should also be taken into account.

4.3 Climate impact of urban wetland and implications for management

In the present study, the urban boreal wetland had an overall SGWP of 0.263 kg CO2-eq m−2 yr−1 which was comparable to or higher than other restored wetlands in the boreal region (Herbst et al., 2013) and within the range of inter-annual variation or lower than restored wetlands in the temperate zone (McNicol et al., 2017; Anderson et al., 2016). Different from other studies, the urban wetland was CO2 neutral and a CH4 source. It is worth noting that the paramount contribution of CH4 in ecosystem SGWP was mainly driven by the large footprint-weighted spatial area of the vegetation (see Sect. 3.2). In fact, The SGWP of GHG emission from open water (1.194 kg CO2-eq m−2 yr−1) was 10 times as large as that from vegetation (0.111 kg CO2-eq m−2 yr−1) (Table 2). The implication of this result is that, during wetland restoration, it would be more beneficial to have large patches of emergent vegetation at least from the GHG emission point of view. Similar results have been obtained by other studies as well, indicating that open water has more climate-warming impact than emergent vegetation due to the large diffusive fluxes from open water (Stefanik and Mitsch, 2014; McNicol et al., 2017). The climate impact of natural wetland depends on the net balance between the cooling effect of CO2 uptake by vegetation and the warming effect of other GHG emissions, mainly CH4 (Bridgham et al., 2013). In wetlands constructed in urban areas, the extent of open water, which is a significant emitter of CO2, should also be taken into consideration when evaluating the role of urban wetland in global climate change.

Our results also showed that total phosphorus enhanced both CO2 uptake and CH4 emission, which have contradictory climate impacts to the ecosystem (Fig. 7b, d). Although it is out of the scope of our study, it would be very interesting to understand the mechanisms and quantify the magnitude and the duration of these enhancements induced by nutrient input. Previous studies have found that nutrient inputs can influence the identity of the key primary producer (submerged plants versus phytoplankton) in the water, which is crucial in shaping the CH4 emission from shallow water (West et al., 2016; Davidson et al., 2018). Submerged plants may decrease CH4 production in the lake by transporting oxygen to the sediment and providing good habitat for CH4-oxidizing bacteria (Heilman and Carlton, 2001), while phytoplankton was shown to significantly increase CH4 ebullition by changing the quality of the dissolved organic carbon which promotes methanogenesis (West et al., 2016) and/or by altering the sediment texture and redox conditions favouring the release of bubbles. As a result, we suggest controlling the nutrient input to the water of the newly established wetland to limit the abundance of phytoplankton as well as to support the existence of submerged plants.

5 Conclusions

Urban wetlands have received global attention as a nature-based urban runoff management solution for sustainable cities, as they provide cost-efficient flood control and water quality mitigation as well as many ecological and cultural services. In the meantime, the climate impact of urban wetlands should also be considered. Wetting a landscape may enhance the CO2 sequestration in the ecosystem, whereas CH4 can be emitted due to the anaerobic conditions in the soil after wetting. Furthermore, heterogeneity induced in newly created urban wetlands may contribute differently to the overall climate impact.

In the present study, for the first time a full annual carbon balance of an urban storm water wetland in the boreal region was evaluated, and the radiative forcing from heterogeneous landscapes was presented. We found that, during the monitored period at the study wetland, both the open-water area and the vegetated area within the created wetland were carbon sources, and thus the urban wetland had a net climate-warming effect in the fourth year after the wetland establishment. However, if the same carbon content from the contributing watershed would have reached the receiving lake without treatment at the studied in-stream created wetland, conversion to CH4 would likely have exceeded emissions observed at the wetland. The radiative forcing effect of the open-water area exceeded the vegetated area, which indicated that limiting open-water surfaces and setting a design preference for areas of emergent vegetation in the establishment of urban wetlands can be a beneficial practice when considering only the climate impact of a created urban wetland. In the meantime, we also emphasize that the value of urban wetlands should not be determined solely by GHG radiative forcing. The values of urban wetlands in other areas e.g. flood control, pollutant removal, biodiversity, recreation and education, are also of paramount importance to human society.

Data availability

Eddy covariance, gas concentration and meteorological data are available from the DRYAD database at https://doi.org/10.5061/dryad.9s4mw6mb6 (Li et al., 2019).

Supplement

The supplement related to this article is available online at: https://doi.org/10.5194/bg-17-3409-2020-supplement.

Author contributions

IM, OW, HV, AO and TV designed the field study. SH, IM and JP carried out eddy covariance measurements, automatic gas concentration measurements in the open-water and manual field measurements. XL and IM participated in eddy covariance data processing and analysis. XL analysed the results and prepared the manuscript with contributions from all co-authors.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

We thank Mikko Yli-Rosti and Kiril Aspila for assistance for the maintenance of the field measurements. This research was supported by the EU Life+11 ENV/FI/911 Urban Oases project grant, Academy of Finland, Academy Professor projects (312571 and 282842), ICOS Finland (by Academy of Finland 281255 and University of Helsinki), the Maa-ja vesitekniikan tuki ry, the Ministry of the Environment of Finland and the municipality of Vihti. In memoriam: The greenhouse gas exchange measurements at the Gateway wetland were made possible due to the creative and caring support by the late Eero Nikinmaa to the Urban Oases project.

Financial support

This research has been supported by the Academy of Finland (grant nos. 312571 and 282842), ICOS Finland (grant no. 281255) and the LIFE+ PROJECT (grant no. ENV/FI/911).

Open access funding provided by Helsinki University Library.

Review statement

This paper was edited by Steven Bouillon and reviewed by Dennis Baldocchi and one anonymous referee.

References

Anderson, F. E., Bergamaschi, B., Sturtevant, C., Knox, S., Hastings, L., Windham-Myers, L., Detto, M., Hestir, E. L., Drexler, J., Miller, R. L., Matthes, J. H., Verfaillie, J., Baldocchi, D., Snyder, R. L., and Fujii, R.: Variation of energy and carbon fluxes from a restored temperate freshwater wetland and implications for carbon market verification protocols, J. Geophys. Res.-Biogeo., 121, 777–795, https://doi.org/10.1002/2015jg003083, 2016. 

Baldocchi, D., Detto, M., Sonnentag, O., Verfaillie, J., Teh, Y. A., Silver, W., and Kelly, N. M.: The challenges of measuring methane fluxes and concentrations over a peatland pasture, Agr. Forest Meteorol., 153, 177–187, https://doi.org/10.1016/j.agrformet.2011.04.013, 2012. 

Baldocchi, D. D.: Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: past, present and future, Glob. Change Biol., 9, 479–492, https://doi.org/10.1046/j.1365-2486.2003.00629.x, 2003. 

Bridgham, S. D., Cadillo-Quiroz, H., Keller, J. K., and Zhuang, Q. L.: Methane emissions from wetlands: biogeochemical, microbial, and modeling perspectives from local to global scales, Glob. Change Biol., 19, 1325–1346, https://doi.org/10.1111/gcb.12131, 2013. 

Chamberlain, S. D., Anthony, T. L., Silver, W. L., Eichelmann, E., Hemes, K. S., Oikawa, P. Y., Sturtevant, C., Szutu, D. J., Verfaillie, J. G., and Baldocchi, D. D.: Soil properties and sediment accretion modulate methane fluxes from restored wetlands, Glob. Change Biol., 24, 4107–4121, https://doi.org/10.1111/gcb.14124, 2018. 

Cole, J. J. and Caraco, N. F.: Atmospheric exchange of carbon dioxide in a low-wind oligotrophic lake measured by the addition of SF6, Limnol. Oceanogr., 43, 647–656, https://doi.org/10.4319/lo.1998.43.4.0647, 1998. 

Cole, J. J., Bade, D. L., Bastviken, D., Pace, M. L., and Van de Bogert, M.: Multiple approaches to estimating air-water gas exchange in small lakes, Limnol. Oceanogr.-Meth., 8, 285–293, https://doi.org/10.4319/lom.2010.8.285, 2010. 

Davidson, T. A., Audet, J., Jeppesen, E., Landkildehus, F., Lauridsen, T. L., Sondergaard, M., and Syvaranta, J.: Synergy between nutrients and warming enhances methane ebullition from experimental lakes, Nat. Clim. Change, 8, 156–160, 2018. 

Erkkilä, K.-M., Ojala, A., Bastviken, D., Biermann, T., Heiskanen, J. J., Lindroth, A., Peltola, O., Rantakari, M., Vesala, T., and Mammarella, I.: Methane and carbon dioxide fluxes over a lake: comparison between eddy covariance, floating chambers and boundary layer method, Biogeosciences, 15, 429–445, https://doi.org/10.5194/bg-15-429-2018, 2018. 

Frolking, S., Roulet, N., and Fuglestvedt, J.: How northern peatlands influence the Earth's radiative budget: Sustained methane emission versus sustained carbon sequestration, J. Geophys. Res.-Biogeo., 111, G01008, https://doi.org/10.1029/2005jg000091, 2006. 

Grinsted, A., Moore, J. C., and Jevrejeva, S.: Application of the cross wavelet transform and wavelet coherence to geophysical time series, Nonlin. Processes Geophys., 11, 561–566, https://doi.org/10.5194/npg-11-561-2004, 2004. 

Heilman, M. and Carlton, R.: Methane oxidation associated with submersed vascular macrophytes and its impact on plant diffusive methane flux, Biogeochemistry, 52, 207–224, 2001. 

Heiskanen, J. J., Mammarella, I., Haapanala, S., Pumpanen, J., Vesala, T., Macintyre, S., and Ojala, A.: Effects of cooling and internal wave motions on gas transfer coefficients in a boreal lake, Tellus B, 66, 1, https://doi.org/10.3402/tellusb.v66.22827, 2014. 

Hemes, K. S., Chamberlain, S. D., Eichelmann, E., Knox, S. H., and Baldocchi, D. D.: A Biogeochemical Compromise: The High Methane Cost of Sequestering Carbon in Restored Wetlands, Geophys. Res. Lett., 45, 6081–6091, https://doi.org/10.1029/2018gl077747, 2018. 

Herbst, M., Friborg, T., Schelde, K., Jensen, R., Ringgaard, R., Vasquez, V., Thomsen, A. G., and Soegaard, H.: Climate and site management as driving factors for the atmospheric greenhouse gas exchange of a restored wetland, Biogeosciences, 10, 39–52, https://doi.org/10.5194/bg-10-39-2013, 2013. 

Holgerson, M. A., Farr, E. R., and Raymond, P. A.: Gas transfer velocities in small forested ponds, J. Geophys. Res.-Biogeo., 122, 5, 1011–1021, 2017. 

IPCC: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1–1535, 2013. 

Kljun, N., Calanca, P., Rotach, M. W., and Schmid, H. P.: A simple two-dimensional parameterisation for Flux Footprint Prediction (FFP), Geosci. Model Dev., 8, 3695–3713, https://doi.org/10.5194/gmd-8-3695-2015, 2015. 

Knox, S. H., Sturtevant, C., Matthes, J. H., Koteen, L., Verfaillie, J., and Baldocchi, D.: Agricultural peatland restoration: effects of land-use change on greenhouse gas (CO2 and CH4) fluxes in the Sacramento-San Joaquin Delta, Glob. Change Biol., 21, 750–765, https://doi.org/10.1111/gcb.12745, 2015. 

Li, X., Wahlroos, O., Haapanala, S., Pumpanen, J., Vasander, H., Ojala, A., Vesala, T., and Mammarella, I.: Data from: Carbon dioxide and methane fluxes from different surface types in a created urban wetland, Dryad, Dataset, https://doi.org/10.5061/dryad.9s4mw6mb6, 2019. 

Linden, A., Heinonsalo, J., Buchmann, N., Oinonen, M., Sonninen, E., Hilasvuori, E., and Pumpanen, J.: Contrasting effects of increased carbon input on boreal SOM decomposition with and without presence of living root system of Pinus sylvestris L, Plant Soil, 377, 145–158, https://doi.org/10.1007/s11104-013-1987-3, 2014. 

Lu, S. Y., Wu, F. C., Lu, Y., Xiang, C. S., Zhang, P. Y., and Jin, C. X.: Phosphorus removal from agricultural runoff by constructed wetland, Ecol. Eng., 35, 402–409, 2009. 

Lucas, R., Earl, E. R., Babatunde, A. O., and Bockelmann-Evans, B. N.: Constructed wetlands for stormwater management in the UK: a concise review, Civ. Eng. Environ. Syst., 32, 251–268, https://doi.org/10.1080/10286608.2014.958472, 2015. 

Mammarella, I., Launiainen, S., Gronholm, T., Keronen, P., Pumpanen, J., Rannik, U., and Vesala, T.: Relative Humidity Effect on the High-Frequency Attenuation of Water Vapor Flux Measured by a Closed-Path Eddy Covariance System, J. Atmos. Ocean. Tech., 26, 1856–1866, https://doi.org/10.1175/2009jtecha1179.1, 2009. 

Mammarella, I., Nordbo, A., Rannik, U., Haapanala, S., Levula, J., Laakso, H., Ojala, A., Peltola, O., Heiskanen, J., Pumpanen, J., and Vesala, T.: Carbon dioxide and energy fluxes over a small boreal lake in Southern Finland, J. Geophys. Res.-Biogeo., 120, 1296–1314, 2015. 

Mammarella, I., Peltola, O., Nordbo, A., Järvi, L., and Rannik, Ü.: Quantifying the uncertainty of eddy covariance fluxes due to the use of different software packages and combinations of processing steps in two contrasting ecosystems, Atmos. Meas. Tech., 9, 4915–4933, https://doi.org/10.5194/amt-9-4915-2016, 2016. 

Mander, U., Dotro, G., Ebie, Y., Towprayoon, S., Chiemchaisri, C., Nogueira, S. F., Jamsranjav, B., Kasak, K., Truu, J., Tournebize, J., and Mitsch, W. J.: Greenhouse gas emission in constructed wetlands for wastewater treatment: A review, Ecol. Eng., 66, 19–35, https://doi.org/10.1016/j.ecoleng.2013.12.006, 2014. 

Matthes, J. H., Sturtevant, C., Verfaillie, J., Knox, S., and Baldocchi, D.: Parsing the variability in CH4 flux at a spatially heterogeneous wetland: Integrating multiple eddy covariance towers with high-resolution flux footprint analysis, J. Geophys. Res.-Biogeo., 119, 1322–1339, https://doi.org/10.1002/2014jg002642, 2014. 

McNicol, G., Sturtevant, C. S., Knox, S. H., Dronova, I., Baldocchi, D. D., and Silver, W. L.: Effects of seasonality, transport pathway, and spatial structure on greenhouse gas fluxes in a restored wetland, Glob. Change Biol., 23, 2768–2782, https://doi.org/10.1111/gcb.13580, 2017. 

Mitsch, W. J. and Gosselink, J. G.: Wetlands, 5th Edn., John Wiley & Sons Inc., Hoboken, NJ, 2015. 

Moffat, A. M., Papale, D., Reichstein, M., Hollinger, D. Y., Richardson, A. D., Barr, A. G., Beckstein, C., Braswell, B. H., Churkina, G., Desai, A. R., Falge, E., Gove, J. H., Heimann, M., Hui, D. F., Jarvis, A. J., Kattge, J., Noormets, A., and Stauch, V. J.: Comprehensive comparison of gap-filling techniques for eddy covariance net carbon fluxes, Agr. Forest Meteorol., 147, 209–232, https://doi.org/10.1016/j.agrformet.2007.08.011, 2007. 

Morin, T. H., Bohrer, G., Frasson, R., Naor-Azreli, L., Mesi, S., Stefanik, K. C., and Schafer, K. V. R.: Environmental drivers of methane fluxes from an urban temperate wetland park, J. Geophys. Res.-Biogeo., 119, 2188–2208, https://doi.org/10.1002/2014jg002750, 2014a. 

Morin, T. H., Bohrer, G., Naor-Azrieli, L., Mesi, S., Kenny, W. T., Mitsch, W. J., and Schafer, K. V. R.: The seasonal and diurnal dynamics of methane flux at a created urban wetland, Ecol. Eng., 72, 74–83, https://doi.org/10.1016/j.ecoleng.2014.02.002, 2014b. 

Mungasavalli, D. P. and Viraraghavan, T.: Constructed wetlands for stormwater management: A review, Fresen. Environ. Bull., 15, 1363–1372, 2006. 

Neubauer, S. C. and Megonigal, J. P.: Moving Beyond Global Warming Potentials to Quantify the Climatic Role of Ecosystems, Ecosystems, 18, 1000–1013, https://doi.org/10.1007/s10021-015-9879-4, 2015. 

Papale, D., Reichstein, M., Aubinet, M., Canfora, E., Bernhofer, C., Kutsch, W., Longdoz, B., Rambal, S., Valentini, R., Vesala, T., and Yakir, D.: Towards a standardized processing of Net Ecosystem Exchange measured with eddy covariance technique: algorithms and uncertainty estimation, Biogeosciences, 3, 571–583, https://doi.org/10.5194/bg-3-571-2006, 2006. 

Pirinen, P., Simola, H., Aalto, J., Kaukoranta, J.-P., Karlsson, P., and Ruuhela, R.: Tilastoja Suomen ilmastosta 1981–2010, Finnish Meteorological Institute, Helsinki, Finland, 2012. 

Rinne, J., Tuittila, E. S., Peltola, O., Li, X. F., Raivonen, M., Alekseychik, P., Haapanala, S., Pihlatie, M., Aurela, M., Mammarella, I., and Vesala, T.: Temporal Variation of Ecosystem Scale Methane Emission From a Boreal Fen in Relation to Temperature, Water Table Position, and Carbon Dioxide Fluxes, Global Biogeochem. Cy., 32, 1087–1106, https://doi.org/10.1029/2017gb005747, 2018. 

Saggar, S., Tate, K. R., Giltrap, D. L., and Singh, J.: Soil-atmosphere exchange of nitrous oxide and methane in New Zealand terrestrial ecosystems and their mitigation options: a review, Plant Soil, 309, 25–42, https://doi.org/10.1007/s11104-007-9421-3, 2008. 

Salminen, O., Ahponen, H., Valkama, P., Vessman, T., Rantakokko, K., Vaahtera, E., Taylor, A., Vasander, H., and Eero N.: TEEB Nordic case: Benefits of green infrastructure – socio-economic importance of constructed wetlands (Nummela, Finland), in: Socio-economic importance of ecosystem services in the Nordic Countries – Scoping assessment in the context of The Economics of Ecosystems and Biodiversity (TEEB), edited by: Kettunen, M., Vihervaara, P., Kinnunen, S., D'Amato, D., Badura, T., Argimon, M., and Ten Brink, P., Nordic Council of Ministers, Copenhagen, 2013. 

Salonen, V.-P. and Varjo, E.: Vihdin Enäjärven kunnostuksen vaikutus pohjasedimentin ominaisuuksiin [The effects of restoration actions at the Lake Enäjärvi in Vihti, Finland on bottom sediment characteristics], Geologi, 52, 159–163, 2000 (Finnish). 

Smith, K. A., Dobbie, K. E., Ball, B. C., Bakken, L. R., Sitaula, B. K., Hansen, S., Brumme, R., Borken, W., Christensen, S., Prieme, A., Fowler, D., Macdonald, J. A., Skiba, U., Klemedtsson, L., Kasimir-Klemedtsson, A., Degorska, A., and Orlanski, P.: Oxidation of atmospheric methane in Northern European soils, comparison with other ecosystems, and uncertainties in the global terrestrial sink, Glob. Change Biol., 6, 791–803, https://doi.org/10.1046/j.1365-2486.2000.00356.x, 2000. 

Stefanik, K. C. and Mitsch, W. J.: Metabolism and methane flux of dominant macrophyte communities in created riverine wetlands using open system flow through chambers, Ecol. Eng., 72, 67–73, https://doi.org/10.1016/j.ecoleng.2013.10.036, 2014. 

Tedford, E. W., MacIntyre, S., Miller, S. D., and Czikowsky, M. J.: Similarity scaling of turbulence in a temperate lake during fall cooling, J. Geophys. Res.-Oceans, 119, 4689–4713, https://doi.org/10.1002/2014jc010135, 2014. 

Torrence, C. and Compo, G. P.: A practical guide to wavelet analysis, B. Am. Meteorol. Soc., 79, 61–78, 1998. 

United Nations: Department of Economic and Social Affairs: Global Sustainable Development Report 2016, New York, 2016. 

Valkama, P., Makinen, E., Ojala, A., Vahtera, H., Lahti, K., Rantakokko, K., Vasander, H., Nikinmaa, E., and Wahlroos, O.: Seasonal variation in nutrient removal efficiency of a boreal wetland detected by high-frequency on-line monitoring, Ecol. Eng., 98, 307–317, 2017. 

Varis, O., Sirvio, H., and Kettunen, J.: Multivariate analysis of lake phytoplankton and environmental factors, Arch. Hydrobiol., 117, 163–175, 1989. 

Vasander, H., Tuittila, E. S., Lode, E., Lundin, L., Ilomets, M., Sallantaus, T., Heikkila, R., Pitkanen, M. L., and Laine, J.: Status and restoration of peatlands in northern Europe, Wetl. Ecol. Manag., 11, 51–63, https://doi.org/10.1023/a:1022061622602, 2003.  

Vickers, D. and Mahrt, L.: Quality control and flux sampling problems for tower and aircraft data, J. Atmos. Ocean. Tech., 14, 512–526, 1997. 

Vohla, C., Alas, R., Nurk, K., Baatz, S., and Mander, U.: Dynamics of phosphorus, nitrogen and carbon removal in a horizontal subsurface flow constructed wetland, Sci. Total Environ., 380, 66–74, 2007. 

Wahlroos, O.: Life+ Urban Oases final project report, available at: http://www.helsinki.fi/urbanoases, last access: 30 June 2019. 

Wahlroos, O., Valkama, P., Mäkinen, E., Ojala, A., Vasander, H., Väänänen, V.-M., Halonen, A., Lindén, L., Nummi, P., Ahponen, H., Lahti, K., Vessman, T., Rantakokko, K., and Nikinmaa, E.: Urban wetland parks in Finland: improving water quality and creating endangered habitats, International Journal of Biodiversity Science, Ecosystem Services & Management, 11, 46–60, https://doi.org/10.1080/21513732.2015.1006681, 2015. 

West, W. E., Creamer, K. P., and Jones, S. E.: Productivity and depth regulate lake contributions to atmospheric methane (Special issue: Methane Emissions from Oceans, Wetlands, and Freshwater Habitats: New Perspectives and Feedbacks on Climate), Limnol. Oceanogr., 61, S51–S61, 2016. 

Download
Short summary
We measured CO2 and CH4 fluxes and quantified the global warming potential of different surface areas in a recently created urban wetland in Southern Finland. The ecosystem has a small net climate warming effect which was mainly contributed by the open-water areas. Our results suggest that limiting open-water areas and setting a design preference for areas of emergent vegetation in the establishment of urban wetlands can be a beneficial practice when considering solely the climate impact.
Altmetrics
Final-revised paper
Preprint