Analysis of atmospheric CH4 in Canadian Arctic and estimation of the regional CH4 fluxes

The Canadian Arctic (> 60 N, 60–141W) may undergo drastic changes if the Arctic warming trend continues. For methane (CH4), Arctic reservoirs are large and widespread, and the climate feedbacks from such changes may be potentially substantial. Current bottom-up and topdown estimates of the regional CH4 flux range widely. This study analyzes the recent observations of atmospheric CH4 from five arctic monitoring sites and presents estimates of the regional CH4 fluxes for 2012–2015. The observational data reveal sizeable synoptic summertime enhancements in the atmospheric CH4 that are distinguishable from background variations, which indicate strong regional fluxes (primarily wetland and biomass burning CH4 emissions) around Behchoko and Inuvik in the western Canadian Arctic. Three regional Bayesian inversion modelling systems with two Lagrangian particle dispersion models and three meteorological datasets are applied to estimate fluxes for the Canadian Arctic and show relatively robust results in amplitude and temporal variations across different transport models, prior fluxes, and subregion masking. The estimated mean total CH4 flux for the entire Canadian Arctic is 1.8± 0.6 Tg CH4 yr−1. The flux estimate is partitioned into biomass burning of 0.3± 0.1 Tg CH4 yr−1 and the remaining natural (wetland) flux of 1.5± 0.5 Tg CH4 yr−1. The summer natural CH4 flux estimates clearly show inter-annual variability that is positively correlated with surface temperature anomalies. The results indicate that years with warmer summer conditions result in more wetland CH4 emissions. More data and analysis are required to statistically characterize the dependence of regional CH4 fluxes on the climate in the Arctic. These Arctic measurement sites will aid in quantifying the inter-annual variations and long-term trends in CH4 emissions in the Canadian Arctic. Copyright statement. The works published in this journal are distributed under the Creative Commons Attribution 4.0 License. This license does not affect the Crown copyright work, which is re-usable under the Open Government Licence (OGL). The Creative Commons Attribution 4.0 License and the OGL are interoperable and do not conflict with, reduce or limit each other. © Crown copyright 2019

Abstract. The Canadian Arctic (> 60 • N, 60-141 • W) may undergo drastic changes if the Arctic warming trend continues. For methane (CH 4 ), Arctic reservoirs are large and widespread, and the climate feedbacks from such changes may be potentially substantial. Current bottom-up and topdown estimates of the regional CH 4 flux range widely. This study analyzes the recent observations of atmospheric CH 4 from five arctic monitoring sites and presents estimates of the regional CH 4 fluxes for 2012-2015. The observational data reveal sizeable synoptic summertime enhancements in the atmospheric CH 4 that are distinguishable from background variations, which indicate strong regional fluxes (primarily wetland and biomass burning CH 4 emissions) around Behchoko and Inuvik in the western Canadian Arctic. Three regional Bayesian inversion modelling systems with two Lagrangian particle dispersion models and three meteorological datasets are applied to estimate fluxes for the Canadian Arctic and show relatively robust results in amplitude and temporal variations across different transport models, prior fluxes, and subregion masking. The estimated mean total CH 4 flux for the entire Canadian Arctic is 1.8 ± 0.6 Tg CH 4 yr −1 . The flux estimate is partitioned into biomass burning of 0.3 ± 0.1 Tg CH 4 yr −1 and the remaining natural (wetland) flux of 1.5 ± 0.5 Tg CH 4 yr −1 . The summer natural CH 4 flux estimates clearly show inter-annual variability that is positively correlated with surface temperature anomalies. The results indicate that years with warmer summer conditions result in more wetland CH 4 emissions. More data and analysis are required to statistically characterize the dependence of regional CH 4 fluxes on the climate in the Arctic. These Arctic measurement sites will aid in quantifying the inter-annual varia-

Introduction
Atmospheric methane (CH 4 ) is one of the principal greenhouse gases with a global warming potential (GWP) 34 times stronger than carbon dioxide (CO 2 ) over a 100-year time period and 96 times stronger over a 20-year time period (Gasser et al., 2017). Atmospheric level of CH 4 has doubled since the pre-industrial era, from about 722 to 1803 ppb in 2011 . Natural wetland CH 4 emission in Arctic regions is of much interest to the scientific community because the CH 4 emissions can potentially increase in a warming climate (AMAP, 2015). The Arctic is underlain with continuous permafrost containing large quantities of soil carbon, ∼ 1700 PgC (Tarnocai et al., 2009). Under warming scenarios, this stored carbon may be highly vulnerable to conversion to CH 4 and CO 2 , which can be emitted to the atmosphere. However, there is only low confidence in the exact magnitude of CO 2 and CH 4 emissions caused by the permafrost thawing and whether carbon will decompose aero-bically to release CO 2 or anaerobically to release CH 4 (e.g., McGuire et al., 2009;Schuur et al., 2015;Thornton et al., 2016).
Natural CH 4 flux estimates are highly uncertain in northern high latitudes. There have been many studies on CH 4 emission using both bottom-up and top-down methods. A thorough review of these studies can be found in Saunois et al. (2016). In general, bottom-up flux estimates for the northern high latitudes from biogeochemical CH 4 models have large variations, with the mean estimates being much higher than top-down estimates from inverse modelling (Saunois et al., 2016). For the boreal North America region including Alaska and the Hudson Bay Lowlands (HBL, the second largest boreal wetland in the world), the bottom-up mean estimate is ∼ 32 Tg CH 4 yr −1 , with a wide range from 15 to 60 Tg CH 4 yr −1 . Conversely, the top-down mean estimate is ∼ 12 Tg CH 4 yr −1 with a narrower range from ∼ 7 to 21 Tg CH 4 yr −1 .
Bottom-up estimates from wetland methane models in WETCHIMP show large discrepancies in the spatial distribution of the wetland CH 4 source, as well as its magnitude (Melton et al., 2013). In the higher latitudes, the limited ground-based information has hindered the mapping of wetland area. Recently, remote sensing provided more information, but the high-latitude wetland extent still contains large uncertainties (Olefeldt et al., 2016;Thornton et al., 2016). In addition to uncertainty in wetland extent, other factors affecting high-latitude wetland emissions in different models still remain. In a recent inter-comparison of CH 4 wetland models (Poulter et al., 2017), all models used the same wetland extent, Surface Water Microwave Product Series (SWAMPS) (Schroeder et al., 2015) with Global Lakes and Wetland Database (GLWD) (Lehner and Döll, 2004), and the same meteorological data (CRU-NCEP v4.0 reconstructed climate data) to drive their models. The models showed a large range in estimated CH 4 emission for the North American boreal-Arctic region, compared to other regions in the world. This large range of the CH 4 emissions for the North American boreal-Arctic region clearly shows the uncertainty in our current understanding of the physical and biogeochemical processes that contribute to wetland CH 4 emissions.
Top-down atmospheric inverse models have been developed to infer fluxes with observed atmospheric CH 4 mixing ratios as constraints. Global CH 4 inversion studies estimate global distribution of emissions and sinks from observational sites from around the world (e.g., Bergamaschi et al., 2013;Bruhwiler et al., 2014) but with very limited observational information in northern high latitudes for these studies. Recently ground-based observational coverage in northern high latitudes improved with the expansion of towers and aircraft atmospheric observing platforms in Arctic regions (e. g. Karion et al., 2016;Sasakawa et al., 2010;Chang et at., 2014). These observations have been used for CH 4 flux estimation in specific regions. For North America, previous atmospheric CH 4 studies were mainly focused on Alaska (e.g., Miller et al., 2016;Hartery et al. 2018). Thompson et al. (2017) conducted CH 4 flux estimates for the entire region north of 50 • N, combining recent high-latitude surface observations. Estimated CH 4 emissions for the Canadian Arctic (representing the land region of Canada north of 60 • N) show discrepancies among the inverse studies; the mean annual total CH 4 emission (2006)(2007)(2008)(2009)(2010) is ∼ 1.8 Tg CH 4 by TM5-4DVAR , 0.5 Tg CH 4 by CarbonTracker-CH 4 , and 2.1 Tg CH 4 by FLEXINVERT (Thompson et al., 2017). Differences including model transports, prior fluxes, and observational datasets could affect the inversion results. The previous CH 4 inversion studies used only observations from Alert in the Canadian Arctic, the most northern observational site in the world. Several new sites in the Canadian Arctic (described in the next paragraph) might be helpful in constraining flux estimation.
Canada has a vast Arctic and sub-Arctic region (∼ 39 % of Canada's total land) with wetlands and permafrost. It is important to study the methane cycle and monitor the effects of climate change in this sensitive region as these effects could impact atmospheric CH 4 levels at national, continental, and hemispheric scales. Environment and Climate Change Canada (ECCC) has recently added five greenhouse gas (GHG) measurement sites in northern Canada to monitor the time evolution of Arctic GHGs in view of better constraining flux estimates in the region. In October 2010, ECCC started continuous GHG measurements at Behchoko (BCK, 115.9 • W, 62.8 • N), the first ground-based site of continuous measurements in the Canadian Arctic other than Alert (ALT, 82.5 • N, 62.5 • W) where continuous measurements were implemented in 1988. Following BCK, additional continuous measurement programs were implemented at Churchill (CHL, 58.7 • N, 93.8 • W) in 2011, Inuvik (INU, 68.3  We present the first study to analyze the atmospheric CH 4 mixing ratios from the above new ECCC observational sites located in the Canadian Arctic region. In this study, we address three key questions: (1) what information can these new measurements provide in regards to local and regional sources? (2) What are the estimated CH 4 fluxes in the Canadian Arctic using inverse modelling with these new measurements? (3) Are there any relationships between the Canadian Arctic CH 4 fluxes and climate-environmental variations? This paper is structured as follows: In Sect. 2, the description of the measurement sites as well as the observational data analyses from daily to inter-annual timescales are given. Section 3 describes the inversion model framework. Section 4 presents flux estimates along with flux uncertainties and potential correlations with climate anomalies.

Measurements
ECCC operates six measurement sites in the Canadian Arctic region to monitor the GHG mixing ratios. Alert (ALT) is the most northern GHG monitoring site on the globe. Weekly flask samples for CO 2 measurements began in 1975. Continuous CH 4 measurement started in 1988. The other five Arctic and sub-Arctic sites, Behchoko (BCK), Churchill (CHL), Inuvik (INU), Cambridge Bay (CBY), and Baker Lake (BLK), gradually became operational starting in 2007. BLK is the most recent site in the Canadian Arctic with continuous measurement started in July 2017, which augmented the flask air sampling measurement program (started in 2014). At the four other sites, continuous measurement systems were initiated during the period of 2010-2012. The observations from these four sites were used for the inversion in this study. Further detailed site information of all six measurement sites is in Table 1. A map showing their locations is in Fig. 1. Currently, all ECCC continuous measurements are conducted using an in situ cavity ring-down spectrometer (CRDS, Picarro G1301, G2301, or G2401), and discrete flask air samples are measured using a gas chromatograph equipped with flame ionization detectors (GC-FID, Agilent 6890). Both measurements are traceable to the World Meteorological Organization (WMO) X2004 scale (Dlugokencky et al., 2005). In the following sections, we describe the sites briefly and characterize the observed variations in the CH 4 mixing ratios at the sites.

Site descriptions
Being located thousands of kilometres from major GHG source regions, Alert (ALT, 82.5 • N, 62.5 • W) is often referred to as an Arctic background site. The Alert observatory is located ∼ 6 km away from the main military base camp. The lack of a local source surrounding the site results in no significant diurnal variation in observed atmospheric CH 4 mixing ratios throughout the year. In winter, under weak vertical mixing, well-defined synoptic variations are observed due to intercontinental-scale transport along with mainly anthropogenic CH 4 originating from the Eurasia continent (Worthy et al., 2009). The measurements at Alert typically represent large-scale background conditions and are thus ideal for providing information on the long-term trend and seasonal cycle for the Arctic (Worthy et al., 2009).
Behchoko (BCK,62.8 • N,115.9 • W) is located on the northwest tip of Great Slave Lake. The continuous measurement program started in October 2012. Flask sampling is not conducted at this site. The air sampling intake is installed at the top of a 60 m communication tower. The observational equipment is placed in a small isolated room in a local backup power generation station that is rarely turned on. BCK is located 10 km away from the town of Behchoko that is a community ∼ 80 km northwest of Yellowknife, the cap-ital of the Northwest Territories. Mixed forests, lakes, and ponds surround the BCK site.
Inuvik (INU,68.3 • N,133.5 • W) is located ∼ 120 km south of the coast of the Arctic Ocean. The continuous measurement program started in February 2012. Flask sampling began in May 2012. The measurement system is located in the ECCC upper air weather station building, 5 km southeast of the town of Inuvik. INU is ecologically surrounded by Arctic tundra and geologically located in the east channel of the Mackenzie Delta, where a number of water streams and ponds are formed and vast hydrocarbon deposits are found. Although there are proposed developments of natural gas and pipeline projects, most have been on hold.
Cambridge Bay (CBY, 69.1 • N, 105.1 • W) is on the southeast coast of Victoria Island. CBY is located ∼ 1 km north of the town of Cambridge Bay, the largest port of the Arctic Ocean's Northwest Passage. Both continuous and flask sampling measurements started in December 2012. The measurement system is located in the ECCC upper air weather station building.
Baker Lake (BKL, 64.3 • N, 96.0 • W) is located on the shore of Baker Lake, ∼ 320 km inland of Hudson Bay. Weekly flask air sampling started in June 2014, and the continuous measurement program began in July 2017. The air sampling system is located at the ECCC upper air weather station. BCK is in an Arctic tundra region, surrounded by small lakes.
Churchill (CHL, 58.7 • N, 93.8 • W) is located on the west coast of Hudson Bay. The GHG monitoring program began with flask air sampling in 2007. The continuous observational program started in October 2011. The sampling equipment is placed in the Churchill Northern Studies Research Facility, ∼ 23 km east of the town of Churchill. CHL is situated in an area with Arctic tundra to the north and on the northern perimeter of the Hudson Bay Lowlands, the largest contiguous boreal wetland region in North America. Figure 2 shows a time series of CH 4 mixing ratios, hourly means, afternoon means (using hourly data from 12:00 to 16:00 local time), and values from flask sampling for each site. The fitted curve and long-term trend were generated using the merged data containing both continuous afternoon means and flask data. The curve-fitting method has two harmonics of 1-year and half-year cycles along with two lowand high-pass digital filters with cut-off periods of 4 months and 24 months respectively (Nakazawa et al., 1997).

Temporal variations
Overall, the features of the continuous and flask measurements show similar long-term trends and seasonal cycles. The continuous measurements reveal short timescale variations that are less visible in the flask data records. The diurnal and synoptic variations in atmospheric CH 4 provide information on local-and regional-scale interactions between the atmosphere and the source fluxes (Chan et al., 2004).  All the sites show similar upward trends of atmospheric CH 4 . The growth rates observed at the Canadian Arctic sites are comparable to the global mean growth rates based on the global network of National Oceanic and Atmospheric Administration (NOAA) air sampling (https://www.esrl.noaa. gov/gmd/ccgg/trends_ch4/, last access: 1 February 2019, Fig. S1 in the Supplement). In 2014, the growth rates jumped at all the sites except BCK. In the following year 2015, the growth rates dropped but were still higher than the growth rates prior to 2014. The rapid enhancement in growth rates at the Canadian Arctic sites is consistent with the globally averaged atmospheric CH 4 . The year 2014 growth rate at BCK was also enhanced, but the enhancement was not as high as those at the other Arctic sites. This moderate growth rate for BCK might be an artefact in its long-term component partially due to a 2-month period of missing data (mid-November 2014 to mid-January 2015).

Seasonal and inter-annual variations
Since the long-term trends reflect the global-scale sourcesink changes, the long-term component at ALT is subtracted from all the sites in order to isolate the regional-scale signals in the observed atmospheric CH 4 data records (Fig. 3). The mean seasonal cycles show a maximum in winter and a minimum in summer. All sites show a similar pattern with a maximum in January-February, while the differences are more noticeable in summer. The summer minimum that is typically representative of the large-scale Arctic background can be seen from July to August at ALT. The summer minima at the other Arctic sites are generally enhanced (with higher CH 4 mixing ratios) relative to ALT and very considerably amongst the sites and from year to year. These enhancements and inter-annual variability are due to the superposition of the enhanced atmospheric CH 4 sink and increased wetland emissions during warm seasons. Minima are seen in June at BCK, INU, and CHL, followed by BKL and CBY with ∼ 1-to 1.5month lags. Another interesting feature seen at INU, BCK, and CHL is a secondary maximum in summer (a summer bump), also indicative of the influence of local and regional wetland and biomass burning emissions. As seen in Fig. 3, these secondary maxima summer bumps vary in timing and amplitude from year to year. The bumps were observed at BCK, INU, and CHL in 2012, which were in phase with each Figure 2. Time series of atmospheric CH 4 mixing ratios at Canadian Arctic sites. The observed values are the hourly means (grey dot), the afternoon means (black dot, 12:00-16:00 local time) from the continuous measurements, and the ones from flask sampling (circle in light blue). BCK has only continuous measurements. At BKL, flask air sampling is only available after being initiated in 2014. The red and green curves are fitted curves and long-term trends which are obtained by applying a fitting-curve method to the observed afternoon means.
other, but not at any site in 2013. In 2014, a more substantial summer bump was observed at BCK than 2012 (2014 has strong biomass burning contributions) while the summer bump at CHL was similar to the one in 2012. The cause(s) for the summer bumps at BCK, CHL, and INU might vary year to year, such as local-regional (wetland and forest fires) emission change due to climate anomaly. Another possible cause is inter-annually varying atmospheric transport.

Synoptic and diurnal variability
All measurements of atmospheric CH 4 in the Canadian Arctic show synoptic and daily variations with seasonally changing amplitudes. One quantitative measure of synoptic variability in the observed CH 4 mixing ratios is the monthly standard deviation (SD) of the residual of observed time series relative to their fitted curves. Figure 4 shows the mean seasonality in SD of all 24 hourly data (SD_24) in CH 4 mixing ratios at each site except BKL, as well as the mean seasonality of the afternoon hourly data (SD_PM). Although SD_24 and SD_PM appear similar (some are almost identical) except during summer months, the differences between SD_PM and SD_24 provide a measure of whether the daily variability is reflecting a local-scale change in emission or rather a change in the atmospheric transport processes. The nighttime planetary boundary layer (PBL) is usually shallow, while the daytime boundary layer is usually deeper and well mixed. If there are local CH 4 sources, the emission will be mixed into a shallow PBL at night (yielding a higher mixing ratio) and diluted through a deeper PBL during the day (yielding a lower mixing ratio). The resultant diurnal variations in the CH 4 mixing ratios are evident as larger CH 4 SD_24 compared to SD_PM. In the absence of local sources, SD_24 is comparable to SD_PM.
The most substantial synoptic variations are observed in summer at all sites except ALT (Fig. 4). This indicates that the major regional CH 4 emissions in the continental Canadian Arctic occur in summer. At ALT, the largest synoptic variations are observed in winter. The winter synoptic variations are mainly due to strong long-range transport from other regions, which has been demonstrated for ALT (Worthy et al., 2009).
The diurnal variability of atmospheric CH 4 is mainly caused by a local CH 4 emission signal modulated by daily PBL development or a temporal change in the local source. In summer, the SD_24 values are higher by > 5 ppb than the SD_PM except for ALT. The larger SD_24 in summer supports the existence of local CH 4 sources around the sites, likely wetland CH 4 emissions. In contrast, the fact ALT has identical SD_24 and SD_PM all year round confirms that there is no significant local source near the site.
Similar to the three continental sites (BCK, INU, CHL), CBY also shows the maxima of SD_24 and SD_PM in summer, but they remain lower than BCK, INU, and CHL but higher than ALT. This indicates that there is a weaker local source of CH 4 around CBY than around the three continental sites. In the cold season (September to May), the SD_24 and SD_PM at CBY are almost identical to ALT. It is noticeable that the SD_24 and SD_PM at BCK, INU, and CHL are still higher than ALT until December. These higher SD_24 and SD_PM values in the first half of the cold season might indicate possible CH 4 emissions from the ground. Zona et al. (2016) suggested that there are ongoing CH 4 emissions from the Alaskan Arctic tundra during the "zero curtain" period when the soil temperature is near zero with average air temperature below 0 • C until the surface is completely frozen.
The SD_24 and SD_PM for winter to spring (January to May) at INU remain higher than the other sites. Also, SD_24 at INU becomes higher than SD_PM from April and remains higher over summer. At the other sites, the difference between SD_24 and SD_PM is seen mainly in the summer months (June-August). This higher variability in atmospheric CH 4 at INU in winter and spring, when the surrounding wetland ecosystem is inactive, likely indicates strong local CH 4 sources, such as anthropogenic CH 4 emissions from natural gas well/refinery facilities. During winter, such local CH 4 signals are amplified by the seasonally calm condition (the mean seasonal cycles of wind speed are shown in Fig. S2) and by reduced vertical mixing within the shallow PBL due to the shorter period of daylight in the polar region. Figure S3 shows the deviations (from the fitted curve) of observed hourly and afternoon mean CH 4 at INU and BCK along with wind speed. In April and May, the difference in deviations between hourly CH 4 and afternoon mean CH 4 becomes larger again after the relatively quiet period. This may indicate the signals of local (anthropogenic) emission around INU being amplified as the PBL diurnal variation starts developing due to longer daytime periods. Another possible local source for the large spring SD_24 and SD_PM at INU may be the natural CH 4 emissions originating from lakes and ponds during the spring thaw (Jammet et al., 2015). In contrast, SD_24 and SD_PM at BCK become smaller as the wind speed increases, indicating a lack of local CH 4 source around BCK in spring.
Since ALT is representative of the Arctic background state in synoptic variability, the difference in SD_24 or SD_PM between ALT and each of the other sites gives a measure of the regional source influence to the site. The sizeable regional source influence signals in summer shown in Fig. 4 should be useful in constraining the regional flux estimation modelling in the next sections.

Regional inversion model description
To estimate the regional CH 4 fluxes in the Canadian Arctic, we apply a Bayesian inversion approach, based on the backward simulations by Lagrangian particle dispersion models (LPDMs). In this study, three different transport models and three prior CH 4 flux distributions were used to help estimate the model uncertainties. The following sections describe the various components of our regional inverse modelling.

Transport models and meteorological data
LPDMs simulate an ensemble of air-following particles which are released from the measurement sites. The air particles travel backwards in time for 5 days with the wind field. Previous studies (e.g., Cooper et al., 2010;Gloor et al., 2001;Stohl et al., 2009) have shown 5 days is typically sufficient to capture the surface influence on a measurement site from the surrounding region. The backward trajectories are used to calculate the footprints as the integrated residence times the particles spent inside the PBL at a resolution of 1.0 • × 1.0 • . We use three different regional model setups combining two different LPDMs, FLEXPART and STILT, and three different meteorological data from the European Centre for Mediumrange Weather Forecasts (ECMWF), Japanese Meteorological Agency (JMA), and Weather Research and Forecasting (WRF) model.
LPDMs simulate local contributions for 5 days prior to the measurements at sites. The background condition of atmospheric CH 4 mixing ratios at the endpoints of the particles is provided by a global model, National Institute for Environmental Studies Transport Model (NIES TM) with global CH 4 flux fields. Below are the details of model setups in this study.

LPDM: FLEXPART_EI
The first model setup is FLEXible PARTicle dispersion model (FLEXPART) (Stohl et al., 2005) driven by reanalysis meteorology from the European Centre for Medium-range Weather Forecasts (ECMWF) ERA-Interim (Dee et al., 2011;Uppala et al., 2005). The input meteorological data are at 3hourly time steps and interpolated to 1.0 • × 1.0 • horizontal resolution with 62 vertical layers.

LPDM: FLEXPART_JRA55
The second model setup is also FLEXPART, but driven by the Japanese 55-year Reanalysis (JRA-55) from Japanese Meteorological Agency (JMA; Kobayashi et al., 2015;Harada et al., 2016). JRA-55 is at 6-hourly time steps and TL319 (∼ 0.5625 • , ∼ 55 km) horizontal resolution with 60 vertical layers. For this study, we use the JRA-55 dataset at the lower resolution (∼ 1.25 • ). This model setup was used for a global inverse modelling system by the Global Eulerian-Lagrangian Coupled Atmospheric Model (GELCA) that is a coupled atmospheric model of NIES TM and FLEXPART . The primary meteorological observational data for JRA-55 have been supplied by ECMWF. In addition to the ECMWF data, the observational data obtained by JMA and other sources are also used.

LPDM:WRF-STILT
The third model setup uses the Stochastic, Time-Inverted, Lagrangian Transport Model (STILT) Lin and Gerbig, 2005). The wind fields to drive STILT are from the Weather Research and Forecasting (WRF) model (Skamarock et al., 2008) at 10 km resolution. Detailed descriptions are found elsewhere (Hu et al., 2019;Miller et al., 2014;Henderson et al., 2015). The footprints are aggregated to 1.0 • ×1.0 • horizontal resolution, similar to the other models in this study. The STILT footprint data are provided from CarbonTracker-Lagrange, which is a Lagrangian assimilation framework developed at the NOAA Earth System Research Laboratory (https://www.esrl.noaa.gov/gmd/ccgg/ carbontracker-lagrange/, last access: 1 February 2019).

Global background model: NIES TM
The background or initial condition for the LPDMs is obtained by sampling a global model of CH 4 at the 5-day back endpoint locations of the LPDM particles. The global background field of the CH 4 mixing ratio is simulated by NIES TM version 8.1i (Belikov et al., 2013) with the optimized CH 4 fluxes with the GELCA-CH 4 inversion system Saunois et al., 2016). The GELCA-CH 4 inverse modelling system optimized the monthly CH 4 fluxes for 2000-2015 to assimilate a global network of surface CH 4 measurements available through the GAW World Data Centre for Greenhouse Gases (WDCGG, http://ds.data. jma.go.jp/gmd/wdcgg, last access: 1 February 2019). The prior CH 4 fluxes for the GELCA-CH 4 global inversion are also used for the regional inversion in this study as described in the later section. The NIES TM has 2.5 • × 2.5 • horizontal resolution and 32 vertical layers, driven by JRA-55. For the global simulation, the CH 4 loss in the atmosphere is included; the stratospheric CH 4 loss and OH oxidation schemes are adapted from a model inter-comparison project "TransCom-CH 4 " (Patra et al., 2011).

Prior fluxes
Three cases of prior emissions, labelled as VIS, GEL, and WetC, are used as listed in Table 2. Since the global background atmospheric CH 4 field was calculated with GELCA-CH 4 inversion posterior fluxes, we chose the prior (VIS) and posterior (GEL) fluxes from GELCA as two cases of prior fluxes in our regional inversion. Note that the continuous CH 4 mixing ratio data from the new Canadian Arctic sites were not used in the GELCA-CH 4 inversion. In this study, the mean wetland fluxes for the last 5 years of the GELCA global model were used, and the prior forest fire CH 4 fluxes are detailed in Sect. 3.2.2. The third prior case (WetC) is the same as GEL but with wetland CH 4 fluxes from WetCHARTs (a recent global wetland methane emission model ensemble for use in atmospheric chemical transport models). WetCHARTs provides inter-annually varying monthly wetland CH 4 fluxes for this study period. The details of prior fluxes are described in the following sections.

Wetland CH 4 fluxes
We used the monthly CH 4 wetland fluxes from two different models. The first model is the Vegetation Integrative Simulator for Trace gases (VISIT) (Ito and Inatomi, 2012). VISIT is a process-based model, using GLWD as wetland extent. In addition to wetland CH 4 flux, VISIT calculates soil CH 4 uptake and CH 4 emission through rice cultivation. The wetland fluxes combined with CH 4 fluxes from rice cultivation were optimized through the GELCA-CH 4 global inversion as a natural CH 4 flux. The second model is WetCHARTs 1.0 (Bloom et al., 2017a). WetCHARTs derives wetland CH 4 fluxes as a function of a global scaling factor, wetland extent, carbon heterotrophic respiration, and temperature dependence (Bloom et al., 2017b). We used the ensemble mean fluxes over 18 model sets which are available for 2001-2015, using (1) three global scaling factors, (2) two wetland extents, GLWD and GLOBCOVER, (3) CARDAMOM (the global CARbon Data MOdel fraMework) as terrestrial carbon analysis, and (4) three temperature-dependent CH 4 respiration functions. The WetCHARTs horizontal resolution is 0.5 • × 0.5 • . The modelled CH 4 fluxes are aggregated into 1.0 • × 1.0 • for this study. Figure 5 shows the spatial distribution of three wetland CH 4 fluxes for the summer months (July-August). Overall, they are similar, while WetCHARTs (WetC) has stronger emissions in Northwest Territories than the two wetland fluxes from VIS and GEL, which are based on VISIT.

Forest fire CH 4 fluxes
GFAS (Global Fire Assimilation System) v1.2 (Kaiser et al., 2012) provides biomass burning (BB) emissions by assimilating fire radiative power (FRP) from the Moderate Resolution Imaging Spectrometer (MODIS). The FRP observations are firstly corrected for data gaps and then linked to dry matter combustion rates with CH 4 emission factors. GFAS has a daily temporal resolution and 0.1 • × 0.1 • horizontal resolution. In this study, the daily fire CH 4 emissions are spatially aggregated into 1.0 • × 1.0 • resolution for the regional inversion, though monthly fluxes were used for the GELCA global inversion.

Anthropogenic emission
The anthropogenic CH 4 emissions are provided by EDGAR (Emission Database for Global Atmospheric Research) v4.2FT2010 (http://edgar.jrc.ec.europa.eu, last access: 1 August 2018), except for rice cultivation. EDGARv4.2FT2010 emission, which is originally at 0.1 • × 0.1 • resolution, is Termites 1 GISS GISS (optimized, as natural 1 ) GISS (optimized, as natural 1 ) VIS used the same prior fluxes with those for global GELCA-CH 4 inversion except biomass burning. GELCA-CH 4 inversion optimized CH 4 fluxes for four source types: 1 natural, 2 soil uptake, 3 anthropogenic, and 4 biomass burning, which are also indicated by superscripted numbers. GEL used the posterior fluxes from global GELCA-CH 4 inversion. For VIS and GEL, the 5-year mean of each source type was used. WetC used WetCHARTs extended mean fluxes as wetland CH 4 , while other fluxes were the same as GEL. For all the scenarios, GFAS v1.2 daily fluxes were used as biomass burning.
aggregated into 1.0 • × 1.0 • . Since the EDGARv4.2FT2010 data are available until 2010, the same values for 2010 are used for the years beyond 2010. The CH 4 emission from rice cultivation was replaced with the one from VISIT-CH 4 and then treated as a part of natural fluxes. Since there is no rice field in the Canadian Arctic and the rest of the North American Arctic-boreal region, the influence of CH 4 emission from rice cultivation in the region of interest in this study is negligible. The difference in the optimized anthropogenic emissions in the Canadian Arctic from the prior by the global GELCA inversion is almost negligible (from 0.0247 to 0.0250 Tg CH 4 yr −1 ). Compared to the wetland emissions, the anthropogenic emissions are substantially smaller and localized (see Fig. 5).

Other natural CH 4 fluxes
For other natural CH 4 fluxes, we use a map of climatological termite emissions from Fung et al. (1991) and modelled soil uptake from VISIT-CH 4 . Because of no termite CH 4 emissions in the Canadian Arctic, termite CH 4 emission has no direct impact but it is included in global simulation for the background CH 4 mixing ratio. The prior soil CH 4 uptake is provided by VISIT-CH 4 as oxidative consumption by methanotrophic bacteria in unsaturated lands. Soil CH 4 uptake has large uncertainty regionally and also globally. Kirschke et al. (2013) reported that global soil uptake ranges from 9 to 47 Tg CH 4 yr −1 . In the Canadian Arctic, the VISITmodelled soil uptake is weak (0.094 Tg CH 4 yr −1 ) but spread widely (Fig. 5). In some parts of the eastern Canadian Arctic, soil uptakes exceed other CH 4 emissions, resulting in negative fluxes/net sink of atmospheric CH 4 .

Regional inversion
In this study, we use the Bayesian inversion approach (Tarantola, 1987;Rodgers, 2000;Enting, 2002). The Bayesian inversion used here optimizes the scaling factors of posterior fluxes by minimizing the mismatch between modelled and observed mixing ratios with constraints and given uncertainties using the cost function (J ) minimization method.
where y (N × 1) is the vector of observations (with the background mixing ratio subtracted; see Sect. 3.1.4). N is the number of time points multiplied by number of stations (N is reduced if observations are missing). λ (R × 1) is the vector of the posterior scaling factors to be estimated; R is the number of fluxes to be solved. R is two fluxes per subregions × number of subregions (i.e., two to six in this study). λ prior is the vector of the prior scaling factors, which are all initialized to 1 for all subregions, and K (N × R) is the matrix of contributions on the observations (N) from all the fluxes (R) of subregions. K is a product of two matrices, M(N ×LL) and x (LL×R), M is the modelled transport (or footprints in this study), and x is the spatial distribution of the surface fluxes. LL (lat × long) is the dimension of our domain (1 • × 1 • in latitude by longitude). A linear regularization term has been added which is the second term on the right-hand side of the equation. D and D prior are the error covariance matrices. D is the prior model-observation error/uncertainty matrix (N × N), where the diagonal elements are (σ e ) 2 . D prior is the prior scaling factor uncertainty matrix (R × R), where the diagonal elements are σ prior 2 . We assume that the model-observation mismatch errors are uncor- related with each other and the contributions from the subregions are uncorrelated. All the off-diagonal elements in D and D prior are assumed to be zero. We assigned σ e = 0.33 for the model-observation error (Gerbig et al., 2013;Lin et al., 2004;Zhao et al., 2009) and σ prior = 0.30 for the prior uncertainty (Zhao et al., 2009). We examined the inversion's sensitivity to these uncertainties by doubling their values. The posterior fluxes changed by less than 5 % for all subregions (and the different subregion masks). The results showed the optimized fluxes are not strongly dependent on these prescribed uncertainties. The estimate for λ is calculated according to the expression below.
The posterior error variance and covariance, post , for the estimates of λ is calculated, We optimize the CH 4 flux from biomass burning and separately optimize the remainder flux (consisting of wetland emission, soil uptake, and anthropogenic emission) per subregion at a monthly time resolution.

Domain/subregions
For the Canadian Arctic, which consists of three territories, Northwest Territories (NT), Yukon (YT), and Nunavut (NU), we set up three subregion masks, masks A, B, and C, as shown in Fig. S4. Outside of the Canadian Arctic is treated as one outer region. Regarding the subdivision of the Arctic region, we examined the sensitivity of the flux estimation to the number of subregions. As a starting point, the three territories are treated separately (mask A). Secondly, YT is combined with NT (mask B). There is no existing measurement site in Yukon and no significant CH 4 emissions in prior fluxes. The inversion results in the next section will show that YT could not be reliably constrained as a separate subregion (model uncertainties made the estimated fluxes in YT fluctu-ate from positive to negative). As the third region mask, we solve the fluxes for one region of the entire Canadian Arctic (mask C). Like YT, NU is a weak source region, compared to NT, and weak observational constraint might lead to unrealistic flux estimates. This exercise on the subdivision gives insights into the constraining power of the existing measurements. Table 3 shows all the inversion experiments in this study. We perform 27 experiments in total with three prior emission cases, three different transport models, and three different subregion masks.

Atmospheric measurements
This regional inversion study used the continuous measurements at BCK, INU, CBY, and CHL for 2012-2015 (Fig. 2). First, the afternoon mean values are calculated by averaging the hourly data over 4 h from 12:00 to 16:00 local time so that the observations we use in this study are more regionally representative assuming midday is in a well-mixed planetary boundary layer. Second, the modelled background mixing ratios, which were described in Sect. 3.1.4, are subtracted from the afternoon mean observations. The differences in mixing ratio between observations and background were input into the regional inversion system as local contributions. The observational data examined in Sect. 2 have already been pre-screened for possible contaminations due to mechanical/technical problems during sampling or analyzing processes. Except for the pre-screening, we did not apply any additional data screening or filtering.

Results and discussions
4.1 Comparison of footprints Figure 6 shows the mean footprints (mean emission sensitivities) of all four sites by the three different LPDMs. There are common features, but there are also noticeable seasonal differences and differences between the models. The spatial coverage is similar, but the sensitivity to emissions around sites depends on the models. Among the models, STILT shows the strongest sensitivity near the sites, while FLEX-PART_JRA55 has the weakest sensitivity. All the footprints near the sites for the winter season are stronger than the summer season. The footprint differences among the models are also significant. STILT appears to be more localized around the sites. These differences indicate that choosing multiple implementations for the atmospheric transport will allow us to reflect some of the uncertainties introduced to our inversion estimate by transport models.

Signals in the observations (relative to background)
The regional inversion depends on how well local signals can be detected in the observations. Therefore, we first look at the detectability of local and regional fluxes in the observed atmospheric CH 4 mixing ratios. If the amplitude of local signals is comparable to the background contribution, estimated regional fluxes would be more uncertain because local signals would be difficult to distinguish from the background contributions. In Sect. 2.2.2, we examined the synoptic variability in observed CH 4 . Here we apply the same procedures to the modelled background CH 4 for the sites to see if the local synoptic signal is distinguishable from the background CH 4 . Figure 7 shows the mean monthly SD of modelled background CH 4 to their fitted curves for the case of FLEX-PART_EI (other model setups are analogous), along with those of observed CH 4 (SD_PM in Fig. 4) for the four sites to be used as observational constraints (BCK, INU, CBY and CHL). In summer, all the SD_PM values of the observations are much larger (up to 3 times) than the respective background SDs, indicating strong local influence. However, in winter, both the observation SD_PM and the background SD are comparable. Thus, the observations could provide more constraints on the estimated regional fluxes in summer than in winter.

Comparison of prior and posterior fluxes with different transport models
The inversion experiments outlined in Table 3 were made to estimate the CH 4 fluxes in the Canadian Arctic using atmospheric observations from the four ECCC sites as mentioned in Sect. 3.3.3. We calculated the posterior flux estimates as the mean of the fluxes estimated in the nine experiments in Table 3 (for each set of subregion masks). The variations in the flux results (standard deviation) are used to represent the flux uncertainty due to transport errors (three transport models) and prior flux errors (three prior emission cases). This flux uncertainty is larger than the posterior flux covariance uncertainty estimates, Eq. (3). Figure 8 shows the monthly posterior fluxes with subregion masks A and B. The monthly posterior fluxes with mask C are shown in Fig. S5, along with the aggregated fluxes with masks A and B for the entire Canadian Arctic. As shown in Fig. 8a, the fluxes in NT are dominant, and all the posterior fluxes in NT show a clear seasonal cycle and inter-annual variations that are reflected in the total fluxes for the entire Canadian Arctic (Fig. S5).
In contrast, no clear seasonal pattern is found for NU and YT ( Fig. 8a and b). The inversion model has difficulty optimizing the weak flux regions. As a result, negative mean fluxes, i.e., CH 4 sinks, could appear, especially in YT (Fig. 8a); the negative biomass burning fluxes are "spurious" since the biomass burning CH 4 source cannot be negative. However, a null flux would be consistent within error bars. Next, the differences and similarities in the inversion results from the three transport models are summarized. The differences in the flux estimates by the three different transport models can be seen in Fig. 9. Figure 9 displays the results for YT + NT with mask B by the three different transport models. FLEXPART_JRA55 tends to esti- Table 3. Experiment configurations. Using each of three different masks (A, B, and C in Fig. S3), nine inversion runs were conducted with a combination of three prior flux cases (VIS, GEL, and WetC on Table 2) and three different models (FLEXPART_EI, FLEXPART_JRA55, and WRF-STILT). In total 27 inversion runs were conducted.  mate higher total fluxes than the other models, resulting in higher emissions by ∼ 0.6 Tg CH 4 yr −1 than the average of ∼ 1.8 Tg CH 4 yr −1 . WRF-STILT tends to yield the lowest estimate among the three models, lower by ∼ 0.5 Tg CH 4 yr −1 than the average. The posterior total fluxes by FLEX-PART_EI appear to be moderate. In winter, the FLEXPART-EI fluxes are close to zero, the same with WRF-STILT. These results are consistent with their footprints (mean emission sensitivities) in Fig. 6. Higher footprint sensitivities near the sites tend to yield lower posterior fluxes and vice versa. The inter-model differences in the posterior forest fire fluxes (biomass burning, BB) are quite noticeable in 2014, which is the extreme fire year in NT. Due to the sporadic nature of the fire events, the differences in transport (transport errors) are evident in the modelled prior CH 4 mixing ratios (Fig. S6c) and could lead to substantial differences in Figure 7. The 4-year (2012-2015) mean monthly SD of modelled background CH 4 mixing ratios and SD of observed CH4 mixing ratios (afternoon data only, SD_PM). The background CH 4 mixing ratios are NIES TM modelled mixing ratios weighted by the endpoints of the 5-day back trajectory. the posterior fluxes (Fig. 9). The WRF-STILT estimated BB in 2014 appears to be moderate (0.23 Tg CH 4 yr −1 ), similar to in 2013 (∼ 0.3 Tg CH 4 yr −1 ), while the other two models show the highest BB flux estimates (0.55-0.67 Tg CH 4 yr −1 ) in 2014, comparable to the prior flux GFAS estimates.

Mask
In contrast, the inter-annual variability in total posterior fluxes is very similar among all three transport model results (as shown in Figs. 8 and S5). The inter-annual variability in the transport models (an intra-model result) appears to be consistent, yielding similar posterior flux inter-annual variability. Since all three different transport models capture this inter-annual variability, it appears to be a robust feature of the CH 4 source-sink in the Canadian Arctic.
Another robust feature appears to be the similarity in the results for the total Arctic emission with different numbers of subregions used in the inversion. The subregion with strong signals in the prior fluxes (NT) and strong observational constraints (BCK and INU within NT) yielded posterior flux results with small uncertainties, while subregions with weak signals in the prior fluxes (YT) and weak observational constraint (no observations in YT) yielded large uncertainties in the posterior flux estimates. A weak subregion like YT could be combined with another subregion (NT) without a strong impact on the inversion results. The temporal varia-tions in the inversion results with different numbers of subregions (an intra-model result) seem to also be a robust feature. Given that the strong observational constraints and the strong wetland emissions are both located in the central part of the Canadian Arctic, representing the Canadian Arctic as a single region was able to yield reasonable inversion results.

Comparison with previous estimates
The estimated fluxes for the entire Canadian Arctic in this study are relatively robust in amplitude and temporal variations even with the different prior fluxes and subregion masking. The mean estimated total CH 4 annual flux for the Canadian Arctic is 1.8 ± 0.6 Tg CH 4 yr −1 . Compared with two previous inversion estimates, our estimate is slightly lower than the mean total flux of 2.14 Tg CH 4 yr −1 (average from 2009 to 2013) inferred by FLEXINVERT regional inversion (Thompson et al., 2017) but much higher than the estimate of 0.5 Tg CH 4 yr −1 (average from 2006 to 2010) from the CarbonTracker-CH 4 global inversion   (Fig. 10a).
All the estimated fluxes are seasonally high around July and August (Fig. 10b). The mean summertime maximum of our estimates is quite consistent with the one by Thompson et al. (2017), but our estimated fluxes have a narrow high summer emission period and low wintertime emission compared with the estimates by Thompson et al. (2017). These temporal differences in estimated fluxes might reflect the observational constraints used in the respective inversions. Thompson et al. (2017) employed a similar type of regional inversion but for the entire northern high latitudes (north of 50 • N). Except for the flask measurement data at CHL, none of the Canadian Arctic sites used in this study were included in Thompson et al. (2017). The strong regional CH 4 signals at INU and BCK in this study appear to yield flux estimates with a narrower high summer emission period and lower wintertime wetland emission compared with the estimates by Thompson et al. (2017).
The flux estimate in this study is partitioned into biomass burning (BB), 0.3±0.1 Tg CH 4 yr −1 , and the remaining flux, 1.5 ± 0.5 Tg CH 4 yr −1 . The remaining flux is mainly natural/wetland CH 4 emissions, given that anthropogenic contribution to the total prior fluxes without BB is ∼ 2 % according to the EDGAR prior fluxes. The estimated wetland flux is comparable to the WetCHARTs (wetland) ensemble mean of 1.35 Tg CH 4 yr −1 (Bloom et al. 2017a, b).
The estimated summertime natural CH 4 fluxes show clear inter-annual variability. The higher emissions are estimated for 2012 and 2014 in this study, which is similar to the results from Carbon Arctic Reservoirs Vulnerability Experiment (CARVE) aircraft measurements over Alaska for 2012 to 2014 (Hartery at al., 2018).

Relationship of fluxes with climate anomalies
Inter-annual variations in estimated CH 4 fluxes are examined in relation to climate parameters, specifically with surface air temperature and precipitation from NCEP reanalysis (Kalnay et al., 1996). First, monthly mean values at the subregions as well as the 4-year mean (2012)(2013)(2014)(2015) for each month are calculated; then the monthly anomalies are computed from the monthly mean values and the 4-year mean of the corresponding month. The temperature and precipitation anomalies are aggregated to the respective regions, NT, YT, and NU. On the regional level, climate anomalies in NT and NU are quite similar, though YT is less similar to NT and NU. YT is mainly covered by mountains with little wetland. Furthermore, among the three Arctic territories, NT has the largest wetland extent and most of the forest fire emissions in 2012-2015. Thus, we look into the correlation in monthly anomalies of CH 4 fluxes with summer climate anomalies in NT.
In Fig. 11, the inter-annual variability of wetland CH 4 fluxes exhibits a moderate positive correlation with the surface temperature anomaly (r = 0.55) and only weakly correlated with precipitation anomalies (r = 0.11). This indicates that the hotter summer weather condition stimulates the wetland CH 4 emission, and precipitation appears to have a less immediate or no direct impact on wetland conditions. In prior cases, VIS and GEL, natural CH 4 fluxes (wetland and other fluxes except biomass burning CH 4 flux) are multi-year mean monthly fluxes. Therefore, these prior fluxes have no yearto-year anomalies and no correlation with the meteorological anomalies. Only in WetC, the prior with wetland CH 4 fluxes from WetCHARTs ensemble mean has inter-annual variations, and the correlations with temperature and precipitation anomalies are r = 0.26 and r = 0.90 respectively. The posterior natural fluxes with WetC show slightly higher correlations (r = 0.55 with temperature, r = 0.16 with precipitation) than the mean correlation values. However, overall there is no clear dependency of posterior correlations on the inherent climate anomaly correlations in the prior fluxes. This result indicates that the inter-annual variations in posterior wetland fluxes in this study are mainly determined by the observations, rather than by prior fluxes.
Inter-annual variations in estimated BB CH 4 fluxes show a negative correlation with precipitation (r = −0.41). Also throughout the fire season (June-September), all estimated BB fluxes negatively correlate with precipitation while the prior BB fluxes appear to have no consistent correlations. The inversion results support that dry conditions would enhance the forest fire. The estimated BB fluxes show a weakly negative correlation with surface temperature (r = −0.23) on midsummer average, but the monthly correlations fluctuate from r = −0.40 to r = 0.47 over the fire season. Since the period is limited in this study (2012)(2013)(2014)(2015), these statistical relationships are still not clear. Also, the relationship of CH 4 emissions with climate conditions could be complex and non-linear (with extreme fire events in some years). More data and analysis are required to characterize the dependency of CH 4 fluxes on climate in the Arctic.

Comparison of modelled and observed mixing ratios
The model-observation statistical comparison is shown with the Taylor diagrams of correlation coefficients and normal- ized standard deviation (NSD) by three different transport models for the four Arctic sites (Fig. 12) using the inversion results with mask B and prior flux case WetC. At BCK and INU, the correlation coefficients and NSD for each model are improved by the inversion. At these two sites, the observations contain large synoptic signals from the Canadian Arctic wetland and provide strong constraints to the inversions. At INU, the improvement for STILT is noticeable, especially with NSD. This improvement of STILT is explained further below. At CBY and CHL, no significant changes between the prior and posterior results are seen. This indicates that the regional flux in the Canadian Arctic only weakly influences CBY and CHL. Further investigation has been conducted for INU. Figure S7a shows the time series of modelled prior and posterior mixing ratios by the three transport models and the observed mixing ratio. The Taylor diagrams in Fig. S7b show the comparison results of modelled mixing ratios with the observations by season, summer months (June-September) and winter months (October-May), as well as for the entire period together. The modelled mixing ratios by STILT with the prior fluxes could be much higher than the mixing ratios by the other two models. That results in the higher prior NSD values, especially in the winter season. The inversion Figure 11. CH 4 flux anomalies vs. surface temperature and precipitation anomalies for summer (July and August). The CH 4 fluxes are July and August posterior fluxes for Northwest Territories (NT) from nine inversion experiments with mask A. Regional climate parameter anomalies in NT are monthly deviations from the 4-year (2012-2015) means. was able to improve the results by reducing the fluxes and consequently the posterior NSD.
Another qualitative measure of the goodness of fit of the model to the observations is the reduced chi-square statistics (Drosg, 2009;Hughes and Hase, 2010). For the limiting case of an infinite number of data points, the data are independent and normally distributed, and the value of reduced chi-square should be 1. The overall reduced chi-squares for all our experiments are in a narrow range of 1.23-1.27. Given that the observations and modelled mixing ratios are not normally distributed (more frequent high mixing ratio events than low mixing ratio events) and the limited number of observations, there does not seem to be a strong reason to reject the model results.

Contributions of background CH 4 mixing ratios on the posterior fluxes
We used the same background CH 4 mixing ratios for the different transport models, which are calculated using the particle endpoints from FLEXPART_JRA55. The idea of using the same background CH 4 fields is to focus on the impact of the local and regional transport contribution on regional inversion, separating from the background contribution. One notable feature in the background CH 4 mixing ratios is the relatively large synoptic variability, especially in winter compared to the observed CH 4 , which might contribute to large uncertainties in flux estimation. To examine how sensitive the inversion results are to these temporal variations in the background CH 4 mixing ratios, additional experiments with background CH 4 with smoothing windows of 5 days, 10 days, and 30 days were made (see Fig. S9a).
The examples of the results are shown in Fig. S9b and c. The posterior fluxes are not strongly dependent on the different background CH 4 mixing ratios. Compared with the unsmoothed background case, slightly small values of NSD are found in the model-observation statistical comparison for summer with 30-day smoothing. It seems there are sufficient observations (signal level) above the background CH 4 mixing ratios (noise level) to constrain the inversion results (Fig. 7).

Conclusions
The Canadian Arctic region is one of the potential enhanced CH 4 source regions related to the ongoing global warming (AMAP, 2015). Earth system models differ in their prediction of how the carbon loss in this region will be split up between CO 2 and CH 4 emissions. Even current bottom-up and topdown estimates of CH 4 flux in the region vary widely. This study 1. analyzed the measurements of atmospheric CH 4 mixing ratios from five sites established in the Canadian Arctic by ECCC to characterize the observed variations and examine the detectability of regional fluxes, 2. estimated the regional fluxes for 4 years (2012)(2013)(2014)(2015) with the continuous observational data of atmospheric CH 4 and also examined the relationship of the estimated fluxes with the climate anomalies.
The observational data analysis reveals large synoptic summertime signals in the atmospheric CH 4 , indicating strong regional fluxes (most likely wetland and biomass burning CH 4 emissions) around Behchoko and Inuvik in Northwest Territory, the western Canadian Arctic. The local signals of atmospheric CH 4 also allow inverse models to optimize biomass burning CH 4 flux (emissions due to forest fire), separately from the natural CH 4 fluxes.
The estimated mean total CH 4 annual flux for the Canadian Arctic is 1.8 ± 0.6 Tg CH 4 yr −1 (wetland flux is 1.5 ± 0.5 Tg CH 4 yr −1 , biomass burning flux 0.3 ± 0.1 Tg CH 4 yr −1 ). The mean total flux in this study is comparable to another regional flux inversion result of 2.14 Tg CH 4 yr −1 by Thompson et al. (2017) but much higher than the global inversion result of 0.5 Tg CH 4 yr −1 by CarbonTracker-CH 4 . The strong regional CH 4 signals at INU and BCK appear to yield flux estimates in this study with a narrower high summer emission period and lower wintertime wetland emission compared with the estimates by Thompson et al. (2017).
Clear inter-annual variability is found in all the estimated summertime natural CH 4 fluxes for the Canadian Arctic, mostly due to Northwest Territories. These summertime flux variations are positively correlated with the surface temperature anomaly (r = 0.55). This result indicates that years with warmer summer conditions result in more wetland CH 4 emissions.
With longer data records and more analysis in the Arctic, inversion CH 4 flux estimates could yield more details on CH 4 emission strength and seasonal cycle (onset and termination of wetland emissions) and dependence of wetland fluxes on climate conditions. More knowledge on the flux and climate relationship could help evaluate and improve bottom-up wetland CH 4 flux models.
Next, we will perform a similar study for the CO 2 measurements from these sites to estimate the Canadian Arctic CO 2 fluxes. Estimation of CO 2 and CH 4 fluxes and monitoring how these fluxes change in the future will improve our understanding of the response of the Arctic carbon cycle to climate change and also yield long-term trends in CO 2 and CH 4 emissions in the Canadian Arctic.
Data availability. The data measured at ECCC sites used in this study are available upon request to Doug Worthy (doug.worthy@canada.ca). The Alert data are also periodically updated to WDCGG, http://ds.data.jma.go.jp/gmd/wdcgg (WDCGG, 2019). The model-related data are available upon request to the corresponding author.
Author contributions. MI and DC designed the research and prepared the paper. MI performed data analysis and inversion experiments. DW led the measurement programs and collected the observational data. MI, DC, and EC provided footprint (potential emission sensitivities) information for inversions. All authors contributed to the discussion and interpretation of the results.