Are the Fenno-Scandinavian Arctic wetlands a significant regional source of formic acid? Atmosphere

: This study presents the ﬁrst gaseous formic acid (HC(O)OH) concentration measurements collected over the Fenno-Scandinavian wetlands (67.9–68.0 ◦ N, 22.1–27.8 ◦ E) as part of the MAMM (Methane and other greenhouse gases in the Arctic-Measurements, process studies and Modelling) aircraft campaigns conducted in August and September 2013. A boundary layer box model approach has been used to calculate a regionally representative (~240 km 2 ) surface ﬂux for HC(O)OH of 0.0098 ( ± 0.0057) mg [HCOOH] · m − 2 · h − 1 . A surface-type classiﬁcation map was used to estimate proportional source contributions to the observed HC(O)OH ﬂux over the measurement region. The removal of expected source contributions (using available literature parameterisations) from the calculated surface ﬂux identiﬁed that 75% remained unaccounted for. This may suggest that HC(O)OH emission from wetland within the Fenno-Scandinavian region could contribute up to 29 times higher per unit area than previous theoretical HC(O)OH globally-averaged wetland estimates, highlighting a need for further constrained wetland studies of HC(O)OH emission to better understand its potentially signiﬁcant impact on the Arctic HC(O)OH budget and consequent impacts on oxidative capacity.


Introduction
Formic acid (HC(O)OH) is considered to be the most abundant carboxylic acid, ubiquitous in the atmosphere. It has been detected in the gaseous [1], liquid [2] and the particulate phases [3] and seen across urban [4], rural [5] and marine environments [5]. It is considered to be a significant component of rain water acidity in remote regions, thereby influencing pH dependent aqueous reactions [3]. Aerosol containing HC(O)OH can act as efficient cloud condensation nuclei (CCN), contributing to the total aerosol loading and thereby influencing the total indirect radiative forcing on the planet [5]. Gaseous HC(O)OH can also influence the oxidative capacity of the atmosphere by the removal of hydroxyl radicals via reaction Equation (1) [6].
Current global models have been found to substantially under-predict HC(O)OH levels in the high latitudes [7,8], and although numerous sources have been identified, their relative contribution remains highly uncertain. Source types include direct emissions from anthropogenic and biogenic processes, and secondary production via alkene ozonolysis [9,10]. The dominant sources in the high latitudes are expected to be tree [11] and soil emission [12], vehicle emission [4,13,14], long range transport of biomass burning emission [15], and ocean emission [5,16]. Previous studies have suggested oceanic algae may also contribute [5,17], however the lack of understanding of the mechanism, reaction pathways, and the large range of different algae species means that the assessment of algae emission is particularly challenging. The principal sinks for HC(O)OH are wet and dry deposition, thus the atmospheric lifetime for HC(O)OH is estimated to be 1.6 days in the boundary layer and approximately 3.2 days in the free troposphere [5].

Arctic Wetland Source of HC(O)OH
HC(O)OH concentrations in the high latitudes have shown elevated levels in the summer which would suggest the missing source contributions to modelled HC(O)OH predictions are likely to be of biogenic origin [7,18]. One such source yet to be quantitatively evaluated is potential emission from wetland soils. Soil bacteria are a well-documented source of HC(O)OH and an active component in wetlands [19]. To our knowledge, no emission flux studies have previously been undertaken, measuring direct HC(O)OH emission from wetlands. Paulot et al. [5] suggested that HC(O)OH emission from soil is directly proportional to the emitted NO concentration from the same source. Using the ratio 0.95 nmol (HC(O)OH) ·nmol −1 (NO) obtained from data collected over dry savannah soil [12], and globally-modelled NO data collected over known regions of wetlands, Paulot et al. [5] inferred a globally-averaged HC(O)OH emission factor of 2.5 × 10 −4 mg [HCOOH] ·m −2 ·h −1 at 303 K for wetland soils. The wet conditions associated with wetland soils are expected to reduce HC(O)OH emission, whereas increased soil acidity and soil temperature are likely to promote emission, although field studies are clearly required to assess these theories [5]. The mean surface temperatures in the high northern latitudes have been shown to increase at twice the rate of the global average over the last 20 years, and these higher temperatures are expected to continue to rise in the future [20,21]. Although HC(O)OH emission from wetlands is poorly understood at present, the increasing temperatures in the Arctic are likely to have significant consequences for temperature dependent sources such as Arctic wetland soils. Wetland is considered the most biologically diverse of all ecosystems and represents a broad spectrum of different types of environments, which include swamps, marshes, bogs and fens, therefore it is apparent that wetlands must be assessed based on their specific environment. Wetlands are defined based on three components; soil, vegetation and hydrology [22]. Of these components, vegetation and hydrology are sensitive to daily conditions, as well as seasonal climatic fluctuations. It is these components that are likely to influence the processes governing HC(O)OH emission. The high affinity of HC(O)OH to water would suggest high levels of water content are likely to reduce or prohibit HC(O)OH emission, therefore hydrology is likely to have a significant role in the process. The Fenno-Scandinavian region spanning Finland and Sweden represents one of the largest wetland areas in Europe and thereby a potentially significant source of HC(O)OH impacting the Arctic atmosphere. As temperatures rise in the high latitudes and the ground thaw reaches greater soil depths, organic chemical species within the wetland soil are likely to become more available for release into the atmosphere [23]. Organic matter tends to accumulate in wetland soils to a greater extent than upland soils because of its high rate of production relative to its rate of decomposition [24]. Furthermore, SoilGrids1KM [25] estimates the Fenno-Scandinavian wetlands to be organically rich (100-300 g·kg −1 organic carbon content) and estimates a mean soil pH ranging between 4.6 and 5.6. Based on previous wetland estimates, these conditions may be considered favourable for the release of HC(O)OH [5]. HC(O)OH emissions in the Fenno-Scandinavian wetlands however, have not been evaluated and the magnitude of this source to the regional budget remains unknown. This study represents the first flux derived from gaseous HC(O)OH measured concentrations over a wetland region in the Fenno-Scandinavian Arctic, with a view to assessing the impact of wetlands on the regional HC(O)OH budget.

Campaign Outline
As part of the MAMM (Methane and other greenhouse gases in the Arctic-Measurements, process studies and Modelling) aircraft campaign, measurements were collected within the Fenno-Scandinavian Arctic, a region proportionally occupied by wetlands. The data presented in this study incorporate two aircraft campaigns based out of Kiruna, Sweden during the periods of 15-18 August and 19-23 September 2013. These were used to derive HC(O)OH mass budget fluxes, with the aim of understanding and reducing the uncertainty of regional emission sources with a particular focus on the Fenno-Scandinavian wetlands.

CIMS Measurement Technique
A Chemical Ionisation Mass Spectrometer (CIMS) (Georgia Institute of Technology, Atlanta, GA, USA) was deployed for in-situ measurement of HC(O)OH on board the FAAM BAe-146 research aircraft at a frequency of 1 Hz. The CIMS instrument was constructed by Georgia Institute of Technology as described by Nowak et al. [26] and set up on the FAAM BAe-146 aircraft as described by Le Breton et al. [9]. This technique utilised iodide anions for the detection of HC(O)OH as described by Slusher et al. [27]. This method followed the same procedure as described by Le Breton et al. [9] and Jones et al. [17] whereby a flow of methyl iodide (CH 3 I) gas mixture was passed over a radioactive 210 Po source, initiating a dissociative electron attachment reaction with CH 3 I to form the reagent anion I − ·(H 2 O) as shown by Equation (2). This anion is passed into the ion molecule region where it combines with HC(O)OH and HCN species via the adduct reactions Equations (3)

Formic Acid and HCN Calibration
Following the same approach as Jones et al. [17], water-specific background calibrations were performed for HC(O)OH and HCN before each flight. Set flows of dry N 2 (ranging 0-7380 SCCM dry N 2 ) were combined with scrubbed ambient air, which altered the I − ·H 2 O ion counts, as well as the I − ·HC(O)OH and I − ·HCN background signal. These measurements produce a linear relationship between I − ·H 2 O ion counts and the detected ion signal for the species; I − ·HC(O)OH and I − ·HCN. This linear relationship is then applied to the measured data to produce a water-corrected instrument background signal of each species for each data point. Water-specific HC(O)OH sensitivity calibrations were also performed at the start, mid-point and end of the campaigns, to determine the HC(O)OH sensitivity corresponding to the I − ·H 2 O signal. Firstly, a water-specific background calibration was conducted to establish the water-corrected instrument background for HC(O)OH. Secondly a HC(O)OH gas standard was added (10 SCCM) to a set flow of scrubbed ambient air (7380 SCCM), where the change in I − ·HC(O)OH counts from its background is used to calculate a sensitivity. I − ·H 2 O counts are then altered by the addition of set flows of dry N 2 where the HC(O)OH sensitivity change is monitored in relation to the change in I − ·H 2 O counts. This calibration shows a linear relationship between I − ·H 2 O counts and HC(O)OH sensitivity. HC(O)OH gas mixtures were prepared using a custom-made manifold and cylinder concentrations were validated after the campaigns using GC-FID analysis. This methodology is described in greater detail by Jones et al. [17].
HCN was calibrated relative to that of HC(O)OH, where the sensitivity of HCN was determined from laboratory calibrations, where a flow of 5 SCCM (BW Technologies HCN calibration cylinder, 10 ppm) was diluted in a flow of ambient air at a stable count of~670,000 for I − ·H 2 O signal, m/z = 145. A second laboratory calibration was performed for HC(O)OH, equivalent to the same concentration as HCN in the same conditions at approximately 670,000 I − ·H 2 O signal. This equated to a sensitivity ratio of 24:1 at I − ·H 2 O equal to 670,000 counts. HCN sensitivity was found to be independent of I − ·H 2 O counts. Therefore, HCN ion count data were normalised for the campaigns based on HC(O)OH sensitivity at 670,000 I − ·H 2 O counts.

CH 4 Measurements (FGGA and WAS Sampling)
In-situ CH 4 measurements (1 Hz) were obtained using a commercially available cavity-enhanced absorption spectroscopy (CEAS, Model RMT-200, Los Gatos Research Inc., San Jose, CA, USA) that was modified for aircraft use O'Shea et al. (2013) provides a full description of the configuration of this instrument on board the FAAM Bae-146 research aircraft, together with the associated data processing and quality control. CH 4 uncertainty is calculated at ±1.3 ppb; 1 Hz precision is ±2.5 ppb [28]. Measurements of CH 4 were also made by analysing whole air samples using cavity-ring down spectroscopy (Model G1301, Picarro Inc., Santa Clara, CA, USA). The mean bias of the whole air samples (400 samples) relative to the in-situ measurements was −0.5 (1 standard deviation = 4.6) ppb for CH 4 . Whole Air Samples (WAS) were collected in stainless steel flasks as described by Lewis et al. [29] and analysed post-flight in the laboratory for the δ 13 C isotopic ratio of CH 4 . The δ 13 C (CH 4 ) isotopic signature analysis was performed by continuous-flow gas chromatography/isotopic-ratio mass spectrometry (CF-GC-IRMS), with a precision of 0.1 , to derive a Keeling plot [30,31], and used to interpret the dominant source type. BAe-146 Core Instruments) CO concentration measurements are collected on board the FAAM BAe-146 aircraft using an AL 5002 VUV Fast Fluorescence CO Analyser (Aero-Laser GmbH, Garmisch-Partenkirchen, Germany) with an uncertainty of ±5% [32]. O 3 data were collected by a TE49C UV Photometric Ozone Analyser (Thermo-Fischer Scientific, Waltham, MA, USA) as part of the core instrumentation on board the aircraft, with an uncertainty of ±3 ppb [33]. Relative humidity is calculated on the aircraft by the FAAM CORE General Eastern GE 1011B Chilled Mirror Hygrometer. Additionally, temperature was recorded using the FAAM CORE Rosemount/Goodrich type 102 Total Air Temperature sensors [34].

Flux Calculation
Previously, aircraft mass budget calculations have been performed successfully to derive regionally representative fluxes of trace species [34][35][36][37][38][39]. For this approach to be applicable, the composition of the air inflow-outflow to a particular region needs to be characterised to determine the enhancement of a species associated with a given region. We employ the same approach as O'Shea et al. [34], whereby measurements are collected along the prevailing wind vector. If it can be assumed that within the PBL (with depth Z 1 ) no vertical gradients exist in the species of interest, a net flux of the species S can be determined as shown by Equation (5).
U (m·s −1 ) represents the mean wind speed, and n (molecules·m −3 ) is the atmospheric number density, which is integrated from the surface to the top of the PBL (dz). The dS (molecules molecules −1 ) term represents the species enhancement for "S" over the distance cosθ. dx parallel to the prevailing wind, where θ is the difference between the mean wind vector and the direction taken by the aircraft along transect x. This enhancement is represented by a concentration gradient determined by averaging the data to 30 s intervals along transect x. The measured variation in the wind speed, HC(O)OH concentration gradient in relation to distance x, wind direction and boundary layer height, defined as one times the standard deviation (1σ) of each term, and are used to calculate an overall flux uncertainty for each flux estimate. The east-west transect was carried out at an angle of 92.6 (±2.6) • which was roughly aligned against the prevailing wind bearing of 315.6 (±26.5) • . The data obtained were reversed in terms of direction to align with the wind direction to ascertain the angle alignment difference of 41.1 • , used in the flux calculation. The mean wind speed during transect x was calculated as 2.9 (±0.9) m·s −1 , as shown in Table 1. This calculation assumes the near-surface emissions are immediately mixed throughout the PBL. Therefore, to validate this assumption, the time taken for air to advect across transect x must be significantly greater than the PBL turnover time.

Case Study B796b
Following the same approach as Karion et al. [40], the mean PBL turnover time (t t ) for species to mix from the surface to the top of the boundary layer, can be estimated by t t = Z 1 /w*, where Z 1 is the PBL height and w* is the convective velocity scale. The convective velocity scale is calculated using Equation (6).
In relation to Equation (6), g signifies the standard acceleration due to gravity at the Earth's surface (9.8 m·s −2 ), w θ V is the surface buoyancy flux (0.0191 K·m·s −1 ), and θ v (284 K) represents the mean virtual potential temperature obtained from the profile ascent (in this case between 13:33-13:42 UTC). The surface buoyancy flux data were collected from the Sodankylä ground research site located at coordinates 67.369 • N, 26.654 • E, at time 13:30 UTC. Based on this calculation the PBL turnover time was estimated to be~22 min. This is significantly shorter than the time taken for air to advect across transect x (~22 h and 53 min). Therefore, the assumption that emissions are immediately mixed throughout the PBL over the length scale of transect x can be considered reasonable.
The PBL height during transect x was derived from potential temperature data recorded by thermodynamic instrumentation on flight B796b, using the ascent and descent vertical profiles at 13:33-13:42 UTC and 17:23-17:32 UTC, respectively. In examining the ascent profile shown in Figure 1

Case Study B797
Following the approach described in Section 2.2, two transect-average flux concentration measurements of HC(O)OH were obtained from flight B797 on the 17 August 2013 over the Fenno-Scandinavian wetlands and used to calculate an overall flux. The first flux measurement was recorded between the coordinates 67.9° N, 26.6° E and 68.0° N, 24.5° E following an east-to-west trajectory over a distance of ~78 km at 08:31:30-08:43 UTC. The average wind direction of 331 (±7)° was aligned against the aircraft bearing of 87 (±5)°. Following the same approach to flight B796b, 180° was added to the aircraft bearing in order to ascertain the alignment difference (63°) between the prevailing wind and the aircraft bearing. The average wind speed during this transect was 6.2 (±0.5) m·s −1 at an average temperature of 283.7 (±0.5) K comparable with flight B796b conditions. The second set of data was obtained on flight B797 over the same region, between coordinates 68.0° N, 24.5° E and 68.0° N, 28.0° E at 08:45-09:07:30 UTC, travelling a distance of ~147 km in a west-to-east direction. The average wind speed during this period was 5.8 (±0.6) m·s −1 , with an average temperature of 284.1 (±0.4) K. An average wind direction of 348.6 (±11.9)° was aligned within a 90° angle to the aircraft bearing of 269.3 (±15.0)°, resulting in an angle difference of 79.3°. The two data sets were collected over the same latitude and longitude, with a minimal time difference between runs. Therefore, to provide a more representative flux, the data sets were combined to calculate an overall flux. Potential temperature data obtained from a descent profile flight on B797 were used to estimate a PBL top of 1310 m with an uncertainty of ±100 m, as shown in Figure 2.

Case Study B797
Following the approach described in Section 2.2, two transect-average flux concentration measurements of HC(O)OH were obtained from flight B797 on the 17 August 2013 over the Fenno-Scandinavian wetlands and used to calculate an overall flux. The first flux measurement was recorded between the coordinates 67.9 • N, 26.6 • E and 68.0 • N, 24.5 • E following an east-to-west trajectory over a distance of~78 km at 08:31:30-08:43 UTC. The average wind direction of 331 (±7) • was aligned against the aircraft bearing of 87 (±5) • . Following the same approach to flight B796b, 180 • was added to the aircraft bearing in order to ascertain the alignment difference (63 • ) between the prevailing wind and the aircraft bearing. The average wind speed during this transect was 6.2 (±0.5) m·s −1 at an average temperature of 283.7 (±0.5) K comparable with flight B796b conditions. The second set of data was obtained on flight B797 over the same region, between coordinates 68.0 • N, 24.5 • E and 68.0 • N, 28.0 • E at 08:45-09:07:30 UTC, travelling a distance of~147 km in a west-to-east direction. The average wind speed during this period was 5.8 (±0.6) m·s −1 , with an average temperature of 284.1 (±0.4) K. An average wind direction of 348.6 (±11.9) • was aligned within a 90 • angle to the aircraft bearing of 269.3 (±15.0) • , resulting in an angle difference of 79.3 • . The two data sets were collected over the same latitude and longitude, with a minimal time difference between runs. Therefore, to provide a more representative flux, the data sets were combined to calculate an overall flux. Potential temperature data obtained from a descent profile flight on B797 were used to estimate a PBL top of 1310 m with an uncertainty of ±100 m, as shown in Figure 2.

Case Study B804
The

Case Study B804
The

Case Study B804
The

Flight B796b: HC(O)OH Flux
This study recorded a HC(O)OH concentration change of approximately 143 ppt ranging from a low of~120 ppt to a high of~262 ppt with a R 2 of 0.5, as shown by Figure 4. Using Equation (5), the calculated HC(O)OH emission flux over the Fenno-Scandinavian wetlands region was estimated to be 0.0098 (±0.0057) mg [HCOOH] ·m −2 ·h −1 on 16 August 2013. The uncertainty in the total flux was determined by propagating the uncertainties associated with the individual terms. The individual uncertainties and total uncertainties associated with each case study are illustrated in Table 1. These included the uncertainties in the observed (fitted) HC(O)OH concentration gradient, the known variability in the wind speed, wind direction and the boundary layer mixing height. HCN concentrations for transect x were comparatively low at 64 (±17) ppt. Following the same approach by Le Breton et al. [41], the 6-sigma method for identifying airmasses of biomass burning origin was applied to the HCN measurements from flight B796b. This method identified an airmass with a HCN concentrations exceeding 109 ppt would be indicative of a biomass burning plume. HCN concentrations during this flight were no greater than 102 ppt suggesting the observed airmass during the flux measurement contained no biomass burning influence. CO and O 3 concentrations remained relatively constant at background levels of 83.9 (±0.4) ppb and 29.1 (±0.9) ppb respectively as shown by Figure 5. No rainfall occurred during these measurement periods thereby removing the effect of wet deposition, which would suggest suitable conditions for capturing a surface emission flux.
The mean relative humidity during this transect was 82.5 (±8.4)%. Interestingly, relative humidity revealed a partial correlation with HC(O)OH concentrations (R 2 = 0.45) along transect x, with an observed sharp decrease of 26.5% from approximately 24.2 • E to 26.1 • E longitude, as shown by Figure 6. This sharp decrease in relative humidity may represent a potential loss mechanism for gaseous HC(O)OH via uptake onto water vapour droplets. However, a recent study by Veres et al. [42] revealed gaseous HC(O)OH concentrations were relatively unchanged under a range of controlled relative humidity conditions (ranging between 0% and 100% relative humidity), therefore for the purpose of this flux calculation this potential loss pathway was not considered.

Flight B796b: HC(O)OH Flux
This study recorded a HC(O)OH concentration change of approximately 143 ppt ranging from a low of ~120 ppt to a high of ~262 ppt with a R 2 of 0.5, as shown by Figure 4. Using Equation (5), the calculated HC(O)OH emission flux over the Fenno-Scandinavian wetlands region was estimated to be 0.0098 (±0.0057) mg[HCOOH]·m −2 ·h −1 on 16 August 2013. The uncertainty in the total flux was determined by propagating the uncertainties associated with the individual terms. The individual uncertainties and total uncertainties associated with each case study are illustrated in Table 1. These included the uncertainties in the observed (fitted) HC(O)OH concentration gradient, the known variability in the wind speed, wind direction and the boundary layer mixing height. HCN concentrations for transect x were comparatively low at 64 (±17) ppt. Following the same approach by Le Breton et al. [41], the 6-sigma method for identifying airmasses of biomass burning origin was applied to the HCN measurements from flight B796b. This method identified an airmass with a HCN concentrations exceeding 109 ppt would be indicative of a biomass burning plume. HCN concentrations during this flight were no greater than 102 ppt suggesting the observed airmass during the flux measurement contained no biomass burning influence. CO and O3 concentrations remained relatively constant at background levels of 83.9 (±0.4) ppb and 29.1 (±0.9) ppb respectively as shown by Figure 5. No rainfall occurred during these measurement periods thereby removing the effect of wet deposition, which would suggest suitable conditions for capturing a surface emission flux.
The mean relative humidity during this transect was 82.5 (±8.4)%. Interestingly, relative humidity revealed a partial correlation with HC(O)OH concentrations (R 2 = 0.45) along transect x, with an observed sharp decrease of 26.5% from approximately 24.2° E to 26.1° E longitude, as shown by Figure 6. This sharp decrease in relative humidity may represent a potential loss mechanism for gaseous HC(O)OH via uptake onto water vapour droplets. However, a recent study by Veres et al. [42] revealed gaseous HC(O)OH concentrations were relatively unchanged under a range of controlled relative humidity conditions (ranging between 0% and 100% relative humidity), therefore for the purpose of this flux calculation this potential loss pathway was not considered.    Methane (CH4) is a well-documented emission from wetlands to the atmosphere [34,43] and represents a useful tracer for wetland emissions. Regional methane sources from wetlands have successfully been identified through Keeling plot analysis of δ 13 C in CH4 from WAS samples [43][44][45]. Methane measurements and WAS samples were also collected during the same observation period. The resulting Keeling analysis revealed a source signature of −72.8 ± 1.9‰, as displayed by Figure 7. Based on previously measured source analyses over the Fenno-Scandinavian wetlands during the summer, this δ 13 C (in CH4) source signature is within the expected range (−75 to −69‰) for a high northern latitude wetland CH4 emission source, thereby adding further credence to the suggestion of a wetland source contributing to the observed HC(O)OH flux. Figures 5 and 8 illustrate HC(O)OH concentrations partially anti-correlate with CH4 concentrations along transect x with an R 2 of 0.37. Influences from sources unique to the individual species may explain these observations. Additionally, several strains of methanogens have been found to use HC(O)OH as a substrate for producing CH4 [46] which may contribute to the detected anti-correlation. CH4 emissions from wetlands are expected to be greater in saturated soil conditions [47], whereas gaseous HC(O)OH emissions are likely to be reduced in saturated conditions due to the hydrophilic properties of HC(O)OH. Therefore, the effects of water may also play a crucial role in governing the emission processes of these species. Constrained laboratory studies monitoring the effect of soil moisture content on HC(O)OH emission would be required to test this hypothesis.  Methane (CH4) is a well-documented emission from wetlands to the atmosphere [34,43] and represents a useful tracer for wetland emissions. Regional methane sources from wetlands have successfully been identified through Keeling plot analysis of δ 13 C in CH4 from WAS samples [43][44][45]. Methane measurements and WAS samples were also collected during the same observation period. The resulting Keeling analysis revealed a source signature of −72.8 ± 1.9‰, as displayed by Figure 7. Based on previously measured source analyses over the Fenno-Scandinavian wetlands during the summer, this δ 13 C (in CH4) source signature is within the expected range (−75 to −69‰) for a high northern latitude wetland CH4 emission source, thereby adding further credence to the suggestion of a wetland source contributing to the observed HC(O)OH flux. Figures 5 and 8 illustrate HC(O)OH concentrations partially anti-correlate with CH4 concentrations along transect x with an R 2 of 0.37. Influences from sources unique to the individual species may explain these observations. Additionally, several strains of methanogens have been found to use HC(O)OH as a substrate for producing CH4 [46] which may contribute to the detected anti-correlation. CH4 emissions from wetlands are expected to be greater in saturated soil conditions [47], whereas gaseous HC(O)OH emissions are likely to be reduced in saturated conditions due to the hydrophilic properties of HC(O)OH. Therefore, the effects of water may also play a crucial role in governing the emission processes of these species. Constrained laboratory studies monitoring the effect of soil moisture content on HC(O)OH emission would be required to test this hypothesis. Methane (CH 4 ) is a well-documented emission from wetlands to the atmosphere [34,43] and represents a useful tracer for wetland emissions. Regional methane sources from wetlands have successfully been identified through Keeling plot analysis of δ 13 C in CH 4 from WAS samples [43][44][45]. Methane measurements and WAS samples were also collected during the same observation period. The resulting Keeling analysis revealed a source signature of −72.8 ± 1.9 , as displayed by Figure 7. Based on previously measured source analyses over the Fenno-Scandinavian wetlands during the summer, this δ 13 C (in CH 4 ) source signature is within the expected range (−75 to −69 ) for a high northern latitude wetland CH 4 emission source, thereby adding further credence to the suggestion of a wetland source contributing to the observed HC(O)OH flux. Figures 5 and 8 illustrate HC(O)OH concentrations partially anti-correlate with CH 4 concentrations along transect x with an R 2 of 0.37. Influences from sources unique to the individual species may explain these observations. Additionally, several strains of methanogens have been found to use HC(O)OH as a substrate for producing CH 4 [46] which may contribute to the detected anti-correlation. CH 4 emissions from wetlands are expected to be greater in saturated soil conditions [47], whereas gaseous HC(O)OH emissions are likely to be reduced in saturated conditions due to the hydrophilic properties of HC(O)OH. Therefore, the effects of water may also play a crucial role in governing the emission processes of these species. Constrained laboratory studies monitoring the effect of soil moisture content on HC(O)OH emission would be required to test this hypothesis.

Source Attribution
To address the influence from other known sources of HC(O)OH, an underlying surface-type classification map [48], was used to identify the vegetation and land use types within the flight transect area, as detailed in Table 2 and illustrated in Figure 9. This map characterised the sampling region for flight B796b to consist of coniferous forests (26.9%), peat bogs (24.3%), mixed forests (22.2%), woodland shrub (12.2%), water bodies (7.1%), broad-leaf forest (5.1%), moors and heathland (1.8%), water courses (0.3%) and pastures and agricultural land (0.1%). In order to hypothesise the contribution of other sources in the measurement region, emission factors of known sources were deducted from the flux estimate.

Source Attribution
To address the influence from other known sources of HC(O)OH, an underlying surface-type classification map [48], was used to identify the vegetation and land use types within the flight transect area, as detailed in Table 2 and illustrated in Figure 9. This map characterised the sampling region for flight B796b to consist of coniferous forests (26.9%), peat bogs (24.3%), mixed forests (22.2%), woodland shrub (12.2%), water bodies (7.1%), broad-leaf forest (5.1%), moors and heathland (1.8%), water courses (0.3%) and pastures and agricultural land (0.1%). In order to hypothesise the contribution of other sources in the measurement region, emission factors of known sources were deducted from the flux estimate.

Source Attribution
To address the influence from other known sources of HC(O)OH, an underlying surface-type classification map [48], was used to identify the vegetation and land use types within the flight transect area, as detailed in Table 2 and illustrated in Figure 9. This map characterised the sampling region for flight B796b to consist of coniferous forests (26.9%), peat bogs (24.3%), mixed forests (22.2%), woodland shrub (12.2%), water bodies (7.1%), broad-leaf forest (5.1%), moors and heathland (1.8%), water courses (0.3%) and pastures and agricultural land (0.1%). In order to hypothesise the contribution of other sources in the measurement region, emission factors of known sources were deducted from the flux estimate.
Although not a conclusive measure, this method offers insight into possible contributions of un-quantified sources such as wetlands or of sources not yet known. Established emissions sources in the region include tree emission and soil emission. The dominant tree species in this region of Fenno-Scandinavia is the Norwegian Spruce tree, a well-established source of HC(O)OH with a standard emission factor of 0.0040 mg [HCOOH] ·m −2 ·h −1 [11]. Assuming coniferous and mixed tree areas (49.1%) are entirely occupied by Norwegian Spruce trees and the emission factor is independent of temperature, we can assume an emission factor contribution (49.1%) from tree species of 0.0020 mg [HCOOH] ·m −2 ·h −1 . Furthermore, soil is another well-known source of HC(O)OH emission and based on Figure 9 and Table 2, there is likely a strong contribution from soil emission in the flux region. Paulot et al. [5] devised a temperature dependent calculation for estimating a HC(O)OH emission factor from soil as shown below, where T is expressed in Celsius for this calculation.
Based on the average temperature recorded during concentration measurements along transect x (283 K, 10 • C), this study hypothesises a soil emission factor of 6.4 × 10 −4 mg (HCOOH] ·m −2 ·h −1 in this region. In this discussion, we assume a soil contribution of 68.3%, resulting in an overall estimated contribution of 4.4 × 10 −4 mg [HCOOH] ·m −2 ·h −1 . This assumes all land represents soil emission of HC(O)OH except for areas identified with a water surface ( Table 2; peat bogs: 24.3%, water bodies: 7.1% and water courses; 0.3%). The estimated emission factors used in this investigation are likely to represent an overestimation of their contribution to the observed flux. However, by deducting these emission factors in proportion to their percentage contribution, an estimated 0.0074 mg [HCOOH] ·m −2 ·h −1 remains unaccounted for, and the hypothesis is that wetland emission may contribute much more significantly than previously thought. Paulot et al. [5] estimates a wetland contribution of 2.5 × 10 −4 mg [HCOOH] ·m −2 ·h −1 , less than 4% of the residual emission rate observed in these measurements, thus suggesting wetlands or other sources, not yet known, may be contributing significantly to the observed flux. Although not a conclusive measure, this method offers insight into possible contributions of un-quantified sources such as wetlands or of sources not yet known. Established emissions sources in the region include tree emission and soil emission. The dominant tree species in this region of Fenno-Scandinavia is the Norwegian Spruce tree, a well-established source of HC(O)OH with a standard emission factor of 0.0040 mg[HCOOH]·m −2 ·h −1 [11]. Assuming coniferous and mixed tree areas (49.1%) are entirely occupied by Norwegian Spruce trees and the emission factor is independent of temperature, we can assume an emission factor contribution (49.1%) from tree species of 0.0020 mg[HCOOH]·m −2 ·h −1 . Furthermore, soil is another well-known source of HC(O)OH emission and based on Figure 9 and Table 2, there is likely a strong contribution from soil emission in the flux region. Paulot et al. [5] devised a temperature dependent calculation for estimating a HC(O)OH emission factor from soil as shown below, where T is expressed in Celsius for this calculation.
Based on the average temperature recorded during concentration measurements along transect x (283 K, 10 °C), this study hypothesises a soil emission factor of 6.4 × 10 −4 mg(HCOOH]·m −2 ·h −1 in this region. In this discussion, we assume a soil contribution of 68.3%, resulting in an overall estimated contribution of 4.4 × 10 −4 mg[HCOOH]·m −2 ·h −1 . This assumes all land represents soil emission of HC(O)OH except for areas identified with a water surface ( Table 2; peat bogs: 24.3%, water bodies: 7.1% and water courses; 0.3%). The estimated emission factors used in this investigation are likely to represent an overestimation of their contribution to the observed flux. However, by deducting these emission factors in proportion to their percentage contribution, an estimated 0.0074 mg[HCOOH]·m −2 ·h −1 remains unaccounted for, and the hypothesis is that wetland emission may contribute much more significantly than previously thought. Paulot et al. [5] estimates a wetland contribution of 2.5 × 10 −4 mg[HCOOH]·m −2 ·h −1 , less than 4% of the residual emission rate observed in these measurements, thus suggesting wetlands or other sources, not yet known, may be contributing significantly to the observed flux.    Table 2. The key for the numbers and colours associated with underlying surface-type classification map (Figure 9), which is used to identify prevalent vegetation and land use types within transect x for flight B696b. * Represents surface types within the flight track for transect x that were too small to be observed on the classification map illustrated in Figure 9.

Flight B797: HC(O)OH Flux
The combined data for flight B797 were used to calculate an overall flux emission rate of 0.0096 (±0.0079) mg [HCOOH] ·m −2 ·h −1 , using the same approach as the flux emission rate calculated for case B796b. Given the difficulties associated with changing meteorological conditions for calculating a flux over the same region between consecutive days, the agreement between flux concentration measurements is excellent. HCN concentrations remained low throughout this transect with an average concentration of 44 (±17) ppt. Furthermore, HCN concentrations remained below the 6-sigma threshold of 88 ppt, suggesting the airmass was not influenced from biomass burning [41]. The HC(O)OH concentration gradient in this flux ranged from a low of 186 ppt in the west to 247 ppt in the east, with an average concentration of 215 (±12) ppt. Concentrations are within a similar range and follow the same trend as the previous flight, with larger concentrations in the east, downwind of the Fenno-Scandinavian wetlands. The concentration change in this flight, however, was much more variable resulting in an R 2 of only 0.03, significantly less than for the previous case (flight B796b). O 3 concentrations followed a similar trend to HC(O)OH, increasing across the west-to-east transect, with a mean of 26.4 (±1.8) ppb. A HC(O)OH enhancement was also observed at longitude 26.5 • E correlating with a 6 ppb rise in O 3 , which could be indicative HC(O)OH production via ozonolysis, however without additional NO x measurements this process cannot be investigated further. Alternatively, based on the overlaying land map of the region, this area at longitude 26.5 • E is classified as the largest area occupied by wetlands, and instead may represent the largest contribution from wetlands emission. The average CO concentration was relatively low at 89.2 (±1.2) ppb during the flight, suggesting local anthropogenic pollution was minimal. The angle difference between the average wind direction and the average aircraft bearing during flight B797 concentration measurements was 74 • . Although this is considered in the flux calculation and uncertainty, the magnitude would suggest that the HC(O)OH concentrations were not directly downwind of transect x, and thus likely influenced by additional sources from outside of the associated flux region. Therefore, considering these uncertainties, flight B797's flux has been included qualitatively in this study to investigate repeatability within the flux concentrations between consecutive days and to better interpret local sources within transect x. Despite this, HC(O)OH concentrations followed a similar concentration range and trend observed in flight B796b, increasing in the eastern area of the measurement region. δ 13 C (CH 4 ) analysis of 10 WAS bottles collected during flight B797, including one sample taken along transect x, produced a Keeling plot isotopic signature of −71.1 ± 0.6 . This is in excellent agreement with isotopic source signatures of boreal wetland emission from this region. As shown in Figure 10a,b, CH 4 measurements did not correlate with HC(O)OH levels, particularly after longitude 26 • E, similar to observations from flight B796b. The mean relative humidity during this combined flux measurement was 87.7 (±2.6)%, and remained relatively unchanged throughout this transect with no correlation with HC(O)OH measurements, suggesting stable conditions. No rainfall occurred during measurements, thus indicating HC(O)OH loss via wet deposition was insignificant. Although there remains a high degree of uncertainty in this second flux estimate, the similarity in HC(O)OH levels between this flux measurement and flight B796b (within a day of each other) adds further credence to the suggestion that these observations are the result of a relatively local source. δ 13 C (CH4) analysis of 10 WAS bottles collected during flight B797, including one sample taken along transect x, produced a Keeling plot isotopic signature of −71.1 ± 0.6‰. This is in excellent agreement with isotopic source signatures of boreal wetland emission from this region. As shown in Figure 10a,b, CH4 measurements did not correlate with HC(O)OH levels, particularly after longitude 26° E, similar to observations from flight B796b. The mean relative humidity during this combined flux measurement was 87.7 (±2.6)%, and remained relatively unchanged throughout this transect with no correlation with HC(O)OH measurements, suggesting stable conditions. No rainfall occurred during measurements, thus indicating HC(O)OH loss via wet deposition was insignificant. Although there remains a high degree of uncertainty in this second flux estimate, the similarity in HC(O)OH levels between this flux measurement and flight B796b (within a day of each other) adds further credence to the suggestion that these observations are the result of a relatively local source.

Flight B804: HC(O)OH Flux
Aircraft measurements were collected between longitudes 23.8° E and 26.7° E, however the data taken between 23.8° E and 25.0° E longitudes revealed sharp enhancements in O3 (~6 ppb), CH4 (~12 ppb) and CO (~6 ppb), as shown by Figure 11. HC(O)OH concentrations were variable and elevated (~60 ppt) between these longitudinal points, with a mean concentration of 98 (±17) ppt. The 6 ppb enhancements of O3 and CO may suggest an anthropogenic influence on the observed measurements. A four-day HYSPLIT backward trajectory at the peak O3 and CO levels (at 16:54 UTC) indicated the air mass originated from a southeast direction, passing through central continental Russia. The enhancements in these trace species may suggest influence from an additional direct source (likely of anthropogenic origin) other than the known sources in this region. CH4 measurements were also found to be initially high, and decreased at longitude 25.0° E. These initial elevations may have interfered with expected background HC(O)OH concentrations. As illustrated in Figure 12, relative humidity measurements were found to shift to a higher average at approximately 25.1° E longitude along transect x, suggesting a possible change of air mass. The average O3 and CO concentrations during the remaining leg of transect x (25.1° E and 26.7° E longitudes) were relatively stable with background concentrations of 25.4 (±1.4) ppb and 107.3 (±0.5) ppb respectively. Relative humidity measurements also remained relatively stable for the remaining leg of this transect with no correlation with HC(O)OH measurements and a mean relative humidity of 91.6 (±2.1)%. Therefore, to provide a more representative surface flux over the region, the aircraft measurements were restricted to data collected between the coordinates 68.

Flight B804: HC(O)OH Flux
Aircraft measurements were collected between longitudes 23.8 • E and 26.7 • E, however the data taken between 23.8 • E and 25.0 • E longitudes revealed sharp enhancements in O 3 (~6 ppb), CH 4 (~12 ppb) and CO (~6 ppb), as shown by Figure 11. HC(O)OH concentrations were variable and elevated (~60 ppt) between these longitudinal points, with a mean concentration of 98 (±17) ppt. The 6 ppb enhancements of O 3 and CO may suggest an anthropogenic influence on the observed measurements. A four-day HYSPLIT backward trajectory at the peak O 3 and CO levels (at 16:54 UTC) indicated the air mass originated from a southeast direction, passing through central continental Russia. The enhancements in these trace species may suggest influence from an additional direct source (likely of anthropogenic origin) other than the known sources in this region. CH 4 measurements were also found to be initially high, and decreased at longitude 25.0 • E. These initial elevations may have interfered with expected background HC(O)OH concentrations. As illustrated in Figure 12 Figure 12. No rainfall occurred during the measurement period suggesting wet deposition was minimal. HCN measurements remained low throughout transect x with an average concentration of 55 (±19) ppt between longitudes 23.8 • E and 26.7 • E, indicative of no biomass burning influence to the observations. Furthermore, the HCN 6-sigma threshold limit for a biomass burning plume was calculated to be 132 ppt, above the observed HCN measurements for the transect, suggesting that biomass burning had no significant influence on the observed measurements [41]. The HC(O)OH concentration measurements between the defined longitudes resulted in a near surface flux of −0.0146 (±0.0029) mg [HCOOH] ·m −2 ·h −1 , indicating a reversal of this region from a source to a sink. The concentration gradient obtained during these measurements generated an R 2 value of 0.08, significantly lower than the R 2 value obtained in flight B796b flux calculation. Nevertheless, the processes governing HC(O)OH emission within wetlands are highly uncertain due to a lack of wetland-specific emission measurements. Furthermore, the poorly defined classification of the wetland soils in this region creates further difficulties in understanding the important variables impacting HC(O)OH emission from such sources. It is widely accepted that HC(O)OH emission from savannah soils is temperature dependent, and therefore may represent an important factor governing wetland emission. The temperatures recorded between these flux measurements were relatively unchanged suggesting the change of this wetlands region from a source to a sink is not the result of a temperature influence, although this may represent an additional factor in the emission process yet to be investigated [5,12]. δ 13 C (CH 4 ) data for WAS samples collected below 1000 m altitude during case study B804 were used to calculate an isotopic source signature of −66.6 (±7.2) using a Keeling plot. The large uncertainty associated with this isotopic source signature is a consequence of the relatively low CH 4 enhancements observed during flight B804. Although this value is higher than the expected isotopic signature range for Fenno-Scandinavian wetlands, there is no isotopic evidence from other flights in this campaign to suggest an alternative δ 13 C (CH 4 ) emission source. Relative humidity measurements were found to be 9.1% higher on average than case study B796b. Although this is not expected to influence gaseous HC(O)OH levels, the differences observed between these flux concentration measurements might suggest it could be an influencing factor within the wetlands emission process. Surface based measurements over a prolonged period of time would be needed in order to assess such an influence on HC(O)OH emission from wetlands. The lower HC(O)OH concentrations observed during flight B804 in comparison to HC(O)OH concentrations in the previous case studies would suggest the observed flux was not a result of an additional source contribution in the background location but rather a reduction in the dominant emission source. As discussed, HC(O)OH's favourability to partition into the liquid phase in the presence of water would suggest soil moisture content is likely to be an important factor influencing HC(O)OH emission from wetlands. Therefore we hypothesise that a change in soil moisture content is likely responsible for the apparent change of this region from a source to a sink, although this theory cannot be confirmed based on the results collected here.
Ultimately, the transition from a source to a sink for this region, as highlighted between case studies B796b and B804, further complicates our understanding of the HC(O)OH emission process from wetlands soils, however, it would appear that wetlands may represent a potentially significant source contributing to the Arctic HC(O)OH budget. burning influence to the observations. Furthermore, the HCN 6-sigma threshold limit for a biomass burning plume was calculated to be 132 ppt, above the observed HCN measurements for the transect, suggesting that biomass burning had no significant influence on the observed measurements [41]. The HC(O)OH concentration measurements between the defined longitudes resulted in a near surface flux of −0.0146 (±0.0029) mg[HCOOH]·m −2 ·h −1 , indicating a reversal of this region from a source to a sink. The concentration gradient obtained during these measurements generated an R 2 value of 0.08, significantly lower than the R 2 value obtained in flight B796b flux calculation. Nevertheless, the processes governing HC(O)OH emission within wetlands are highly uncertain due to a lack of wetland-specific emission measurements. Furthermore, the poorly defined classification of the wetland soils in this region creates further difficulties in understanding the important variables impacting HC(O)OH emission from such sources. It is widely accepted that HC(O)OH emission from savannah soils is temperature dependent, and therefore may represent an important factor governing wetland emission. The temperatures recorded between these flux measurements were relatively unchanged suggesting the change of this wetlands region from a source to a sink is not the result of a temperature influence, although this may represent an additional factor in the emission process yet to be investigated [5,12]. δ 13 C (CH4) data for WAS samples collected below 1000 m altitude during case study B804 were used to calculate an isotopic source signature of −66.6 (±7.2)‰ using a Keeling plot. The large uncertainty associated with this isotopic source signature is a consequence of the relatively low CH4 enhancements observed during flight B804. Although this value is higher than the expected isotopic signature range for Fenno-Scandinavian wetlands, there is no isotopic evidence from other flights in this campaign to suggest an alternative δ 13 C (CH4) emission source. Relative humidity measurements were found to be 9.1% higher on average than case study B796b. Although this is not expected to influence gaseous HC(O)OH levels, the differences observed between these flux concentration measurements might suggest it could be an influencing factor within the wetlands emission process. Surface based measurements over a prolonged period of time would be needed in order to assess such an influence on HC(O)OH emission from wetlands. The lower HC(O)OH concentrations observed during flight B804 in comparison to HC(O)OH concentrations in the previous case studies would suggest the observed flux was not a result of an additional source contribution in the background location but rather a reduction in the dominant emission source. As discussed, HC(O)OH's favourability to partition into the liquid phase in the presence of water would suggest soil moisture content is likely to be an important factor influencing HC(O)OH emission from wetlands. Therefore we hypothesise that a change in soil moisture content is likely responsible for the apparent change of this region from a source to a sink, although this theory cannot be confirmed based on the results collected here.
Ultimately, the transition from a source to a sink for this region, as highlighted between case studies B796b and B804, further complicates our understanding of the HC(O)OH emission process from wetlands soils, however, it would appear that wetlands may represent a potentially significant source contributing to the Arctic HC(O)OH budget.

Conclusions
This work presents the first flux derived from gaseous HC(O)OH measured concentrations over a wetland region in the Fenno-Scandinavian area. A best-case scenario (flight B796b) with the most favourable meteorological conditions led to an estimated emission of 0.0098 (±0.0057) mg[HCOOH]·m −2 ·h −1 . The removal of (upper limit estimation) emission factors of known sources in the region (tree and soil emission) from the overall flux calculation revealed 0.0074 mg[HCOOH]·m −2 ·h −1 still unaccounted for. Due to the lack of knowledge of wetland HC(O)OH production processes and release pathways, it is difficult to assign a HC(O)OH emission rate solely to wetland emission, however it would appear that wetland emission is likely to be an important contributing factor. A recent estimation of HC(O)OH wetland emission approximates an emission factor of 2.5 × 10 −4 mg[HCOOH]·m −2 ·h −1 [5], whereas work from this study estimates an emission factor over 29 times greater than this. This highlights the need for additional research on HC(O)OH release from wetlands, of which constrained ground-based flux studies are best suited to complement this work. The additional flux measurement in September (case study B804) identified a change in this region from a surface source to a surface sink. The average temperature and relative humidity between the August and September flux concentration measurements were very similar suggesting these components are unlikely to be the driving factors behind the observed change in this region from a source to a sink.
Delineation of wetland soils in this region may also provide further insight into the emission process. The lack of surface measurements within the wetland region, however, make it difficult to hypothesise a single cause for this change, however, based on the available information, we hypothesise that soil moisture content is likely an important factor governing this process change. From this work, it is clear that further studies are required to fully assess wetlands as a source of HC(O)OH and to investigate its impact on the regional atmospheric budget, particularly in the Arctic region. This work would suggest seasonal, ground-based flux chamber studies in this wetland region are required to better constrain and interpret gaseous HC(O)OH emission from this source type. However, based on the observed HC(O)OH flux from case study B796b, wetlands may represent a much greater emission source of HC(O)OH in the Arctic summer, than previously thought.

Conclusions
This work presents the first flux derived from gaseous HC(O)OH measured concentrations over a wetland region in the Fenno-Scandinavian area. A best-case scenario (flight B796b) with the most favourable meteorological conditions led to an estimated emission of 0.0098 (±0.0057) mg [HCOOH] ·m −2 ·h −1 . The removal of (upper limit estimation) emission factors of known sources in the region (tree and soil emission) from the overall flux calculation revealed 0.0074 mg [HCOOH] ·m −2 ·h −1 still unaccounted for. Due to the lack of knowledge of wetland HC(O)OH production processes and release pathways, it is difficult to assign a HC(O)OH emission rate solely to wetland emission, however it would appear that wetland emission is likely to be an important contributing factor. A recent estimation of HC(O)OH wetland emission approximates an emission factor of 2.5 × 10 −4 mg [HCOOH] ·m −2 ·h −1 [5], whereas work from this study estimates an emission factor over 29 times greater than this. This highlights the need for additional research on HC(O)OH release from wetlands, of which constrained ground-based flux studies are best suited to complement this work. The additional flux measurement in September (case study B804) identified a change in this region from a surface source to a surface sink. The average temperature and relative humidity between the August and September flux concentration measurements were very similar suggesting these components are unlikely to be the driving factors behind the observed change in this region from a source to a sink.
Delineation of wetland soils in this region may also provide further insight into the emission process. The lack of surface measurements within the wetland region, however, make it difficult to hypothesise a single cause for this change, however, based on the available information, we hypothesise that soil moisture content is likely an important factor governing this process change. From this work, it is clear that further studies are required to fully assess wetlands as a source of HC(O)OH and to investigate its impact on the regional atmospheric budget, particularly in the Arctic region. This work would suggest seasonal, ground-based flux chamber studies in this wetland region are required to better constrain and interpret gaseous HC(O)OH emission from this source type. However, based on the observed HC(O)OH flux from case study B796b, wetlands may represent a much greater emission source of HC(O)OH in the Arctic summer, than previously thought. Rebecca Fisher and Dave Lowry. Thomas Bannan's contribution concerns the development and operation of the CIMS instrument, which provided the formic acid data used in the study. All authors participated in the data analysis and interpretation and contributed to the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.