Articles | Volume 14, issue 2
Atmos. Meas. Tech., 14, 945–959, 2021
Atmos. Meas. Tech., 14, 945–959, 2021

Research article 08 Feb 2021

Research article | 08 Feb 2021

Quantifying fugitive gas emissions from an oil sands tailings pond with open-path Fourier transform infrared measurements

Quantifying fugitive gas emissions from an oil sands tailings pond with open-path Fourier transform infrared measurements
Yuan You1,a, Samar G. Moussa1, Lucas Zhang2, Long Fu2, James Beck3, and Ralf M. Staebler1 Yuan You et al.
  • 1Air Quality Research Division, Environment and Climate Change Canada (ECCC), Toronto, M3H 5T4, Canada
  • 2Alberta Environment and Parks, Edmonton, T5J, 5C6, Canada
  • 3Suncor Energy Inc., Calgary, T2P 3Y7, Canada
  • anow at: Department of Physics, University of Toronto, Toronto, M5S 1A7, Canada

Correspondence: Ralf M. Staebler (


Fugitive emissions from tailings ponds contribute significantly to facility emissions in the Alberta oil sands, but details on chemical emission profiles and the temporal and spatial variability of emissions to the atmosphere are sparse, since flux measurement techniques applied for compliance monitoring have their limitations. In this study, open-path Fourier transform infrared spectroscopy was evaluated as a potential alternative method for quantifying spatially representative fluxes for various pollutants (methane, ammonia, and alkanes) from a particular pond, using vertical-flux-gradient and inverse-dispersion methods. Gradient fluxes of methane averaged 4.3 g m−2 d−1 but were 44 % lower than nearby eddy covariance measurements, while inverse-dispersion fluxes agreed to within 30 %. With the gradient fluxes method, significant NH3 emission fluxes were observed (0.05 g m−2 d−1, 42 t yr−1), and total alkane fluxes were estimated to be 1.05 g m−2 d−1 (881 t yr−1), representing 9.6 % of the facility emissions.

1 Introduction

Tailings from the oil sands industrial processes in Alberta's Athabasca oil sands consist of a mixture of water, sand, non-recovered bitumen, and additives from the bitumen extraction processes (Small et al., 2015). These tailings are deposited into large engineered tailings ponds on site. Separation of processed water from remaining tailings occurs continuously in the tailings pond, and the processed water is recycled (Canada's oil sands tailings ponds:, last access: 29 September 2019). The total liquid surface area covered by tailings ponds in the Athabasca oil sands was 103 km2 in 2016 and continues to grow (Alberta Environment and Parks, 2016).

Emissions to the atmosphere from tailings ponds include methane (CH4), carbon dioxide (CO2), reduced sulfur compounds, volatile organic compounds (VOCs), and polycyclic aromatic hydrocarbons (PAHs) (Siddique et al., 2007; Simpson et al., 2010; Yeh et al., 2010; Siddique et al., 2011, 2012; Galarneau et al., 2014; Small et al., 2015; Bari and Kindzierski, 2018; Zhang et al., 2019; Foght et al., 2017; Penner and Foght, 2010). Emissions from tailings ponds vary with pond conditions, such as pond age and solvent additives in the ponds, and can contribute significantly to total facility emissions (Small et al., 2015).

Very few studies focusing on emissions of air pollutants from tailings ponds have been published (Galarneau et al., 2014; Small et al., 2015; Zhang et al., 2019). Compounds of particular interest include alkanes and ammonia (NH3). Alkanes are part of the solvents used in the extraction process (Small et al., 2015) and can dominate VOCs emissions from oil sands facilities (Li et al., 2017). Previously reported VOCs emissions by facilities had large uncertainties, especially from fugitive sources, due to limitations of the methods used to estimate emissions for compliance monitoring purposes (Li et al., 2017). VOCs in the atmosphere are important because of their effects on ambient ozone and secondary aerosol formation (Field et al., 2015; Kroll and Seinfeld, 2008). Emissions of NH3 from tailings ponds to the atmosphere have not been published, although NH3 has been observed in the oil sands region (Bytnerowicz et al., 2010; Whaley et al., 2018). NH3 emissions have important environmental implications, such as forming atmospheric aerosols with sulfuric acid (Kürten, et al., 2016) and affecting nitrogen deposition in the ecosystem (Makar, et al., 2018). This information is important for model simulations of critical loads of acidifying deposition in the ecosystem (Makar, et al., 2018). This field measurement project provided a great opportunity to continuously measure and to quantify tailings pond emissions over more than a month, especially for NH3 and total alkanes.

Open-path Fourier transform infrared (OP-FTIR) spectroscopy has been considered a good candidate for an alternative method to monitor fugitive emissions from industrial or hazardous waste area sources, since the method is non-intrusive, integrates over long path lengths, and has the ability to quantify several different gases of interest simultaneously and continuously (Marshall et al., 1994), without sample line issues. It has previously been used to quantify mole fractions of various air pollutants from different sources such as forest fires (Griffith et al., 1991; Yokelson et al., 1996, 1997; Goode et al., 1999; Yokelson, 1999; Yokelson et al., 2007; Burling et al., 2010; Johnson et al., 2010; Akagi et al., 2013; Yokelson et al., 2013; Akagi et al., 2014; Paton-Walsh et al., 2014; Smith et al., 2014), volcanoes (Horrocks et al., 1999; Oppenheimer and Kyle, 2008), industrial sites (Wu et al., 1995), harbors (Wiacek et al., 2018), and road vehicles (Bradley et al., 2000; Grutter et al., 2003; You et al., 2017). OP-FTIR measurements with vertically separated paths have previously been conducted to derive the emission rate of air pollutants. Schäfer et al. (2012) deployed two OP-FTIR spectrometers with parallel paths 2.2 m vertically apart at a grassland in Fuhrberg, Germany, to measure nitrous oxide (N2O) emissions with a flux-gradient method and showed the calculated flux is comparable to the chamber measurements at the same grassland. Flesch et al. (2016) deployed OP-FTIR measurement with one spectrometer and two paths vertically separated by about 1 m on average (slant path configuration) at a cattle field in Alberta, Canada. They derived emission rates of N2O and NH3 by flux-gradient and inverse-dispersion methods, demonstrating the capability of OP-FTIR systems to measure emission rates of N2O and NH3. Following the flux-gradient method in Flesch et al. (2016), Bai et al. (2018) measured the flux of N2O, NH3, CH4, and CO2 from a vegetable farm in Australia using an OP-FTIR system with two paths vertically separated by 0.5 m on average. At the same vegetable farm, Bai et al. (2019) measured emission rates of N2O using flux chambers and OP-FTIR slant path configuration with flux-gradient methods and showed a large variation of the ratio of N2O fluxes with these two measurements. Inverse-dispersion models have also been applied to OP-FTIR measurements to quantify emission rates in previous studies (Flesch et al., 2004, 2005; Bai et al., 2014; Hu et al., 2016; Shonkwiler and Ham, 2018).

Longer continuous coverage with a greater height difference between paths is one distinguishing feature of this study compared to previous research. The motivation of this work is to quantify emission rates of pollutants from one specific tailings pond by combining OP-FTIR measurements with micrometeorological methods. Emissions of CH4, NH3, and total alkanes as well as a comparison of gradient and inverse-dispersion methods are presented in this study.

2 Open-path FTIR field measurements and methods for deriving fluxes

2.1 Site and measurement setup

The main site of this study was on the south shore of Suncor Pond 2/3 (Fig. 1; 56590.90′′ N, 1113030.30′′ W, 305 m a.s.l.). Sensible heat and momentum fluxes were measured on a mobile tower with sonic anemometers (model CSAT3, Campbell Scientific, USA) at 8, 18, and 32 m above ground. Vertical gradients of gaseous pollutant mole fractions were measured by drawing air from 8, 18, and 32 m on the tower to instrumentation housed in a trailer on the ground. A fourth sample inlet at 4 m was on the roof of the main trailer beside the mobile tower. CH4 and H2O mole fractions at these four levels were measured sequentially by cavity ring-down spectroscopy (CRDS) (model G2204, Picarro, USA). CH4 mole fractions at 4 m were used to calibrate the mole fraction from OP-FTIR retrievals. At 18 m, CH4 mole fractions were also measured by another CRDS (model G2311-f, Picarro, USA) at 10 Hz, to be combined with sonic anemometer measurements to calculate the eddy covariance (EC) flux. Meteorological parameters including temperature and relative humidity (RH) were measured at the same three levels on the tower and 1 m above ground. A propeller anemometer (model 05103-10, Campbell Scientific, USA) on the roof of the main trailer at 4 m above ground provided an additional measurement of wind speed and direction. Measurements were conducted from 28 July to 5 September 2017. The FTIR spectrometer was located 10 m to the east of the flux tower, and the paths were along the south shore of the pond. This paper focuses on derived fluxes from the measurement of OP-FTIR. Other experimental details of the project can be found in You et al. (2020).

Figure 1Map of the site: open-path FTIR on the south shore marked in red, open-path FTIR of Alberta Environment and Parks (AEP) on the north shore marked in yellow, and the outfall on the west side. The colored rose plot shows 50 % and 80 % contribution distances for eddy covariance fluxes at 18 m using the Flux Footprint Prediction (FFP) model (Kljun et al., 2015). The unit of contribution distances is in meters. The white circle labels the 2.3 km distance from the main site, equivalent to 100 times the height of the highest point of the top FTIR path. This figure is adapted from You et al. (2020), Fig. 1.

2.2 Open-path Fourier transform infrared spectrometer (OP-FTIR) system

The FTIR measurements were taken with a commercial Open Path FTIR Spectrometer (Open Path Air Monitoring System (OPS), Bruker, Germany), which was set up at 1.7 m above the ground in a trailer. The infrared source is an air-cooled Globar. The emitted radiation is directed through the interferometer where it is modulated, travels along the measurement path (200 m horizontal distance) to a retroreflector array that reflects the radiation, travels back to the spectrometer, and enters a Stirling-cooled mercury cadmium telluride (MCT) detector (monostatic configuration). Three retroreflectors were employed in this study: one near ground level (1 m) on a tripod and two at higher elevations on basket lifts, resulting in heights of reflectors of approximately 1, 11, and 23 m above ground. Three paths with these three retroreflectors are referred to as bottom, middle, and top paths. The bottom retroreflector was approximately twice the size of the upper two (59 vs. 30 reflector cubes). All retroreflectors were cleaned with an alcohol solution once during the study, and the bottom mirror was rinsed with de-ionized water three times. Return signal strength decreased by around 65 % during the 5-week study due to reflector deterioration, presumably mostly due to impaction by particulate matter. This reflector deterioration also decreased the signal-to-noise ratio by around 67 %, based on spectral retrievals for CH4, but did not affect the mean mole fractions measured.

In this study, spectra were measured at a resolution of 0.5 cm−1 with 250 scans co-added to increase signal-to-noise ratio, resulting in roughly a 1 min temporal resolution. Stray light spectra were recorded regularly by pointing the spectrometer away from the retroreflectors. This stray light spectrum accounts for radiation back to the detector from internal reflections inside the spectrometer, i.e., not from the retroreflector array, and was subtracted from all the measurement spectra before performing further analysis.

Table 1Spectral windows of OP-FTIR spectra for retrieving mole fractions of pollutants in this study.

a Threshold correlation coefficient is a input for OPUS_RS when performing fitting analysis of FTIR spectra. When the correlation coefficient between measured spectrum and reference spectrum with the defined spectral window is below this threshold, that pollutant is not identified and the mole fraction is reported as zero in OPUS_RS (You et al., 2017). b Butane and octane mixing ratio are quantified as two surrogates to quantify a total alkane mixing ratio = butane + octane (Thoma et al., 2010). c Detection limit is calculated by converting 3σ of the noise of the measurements with a retroreflector distance of 225 m by Bruker to 3σ of the noise with 200 m in this study.

Download Print Version | Download XLSX

Spectral fitting was performed with OPUS_RS (Bruker), which uses a non-linear curve-fitting algorithm (You et al., 2017). Spectral windows and interference gases for each gas (Table 1) were determined by optimizing capture of the absorption features while minimizing interferences. To further improve fittings, baselines were optimized through either linear or Gaussian fits under given spectral windows and interfering gases. For CH4 and CO2, temperature-dependent reference files were used for fitting and retrieving mole fractions. For other pollutants, reference spectra at 296K were used and retrieved mole fractions were corrected for air density using measured ambient temperature and pressure. For the bottom path, the retrieved mole fractions were corrected for the temperature effect given the temperature difference between the temperatures at 8 m (used as the input temperature) and 1 m. For the top path, the retrieved mole fractions were corrected for the temperature effect given the temperature difference between the input temperature at 8 m and the temperature at 12 m (linearly estimated by using temperature at 8 and 18 m). The corrected mole fractions from the three paths were then converted to dry mole fractions using the RH, temperature, and pressure measured at 1, 8, and 18 m. The dry CH4 mole fractions from FTIR were then calibrated against CH4 mole fractions from point cavity ring-down spectrometer (CRDS) measurements (Picarro G2204) at 4 m (Sect. S1.1 and Fig. S1 in the Supplement). These calibrated CH4 mole fractions from the FTIR were then used in flux calculations.

We also attempted to retrieve several other pollutants from measured FTIR spectra but encountered insufficient signal-to-noise ratios, given the existing mole fractions at this location, variability in ambient H2O vapor, etc. These pollutants include toluene, benzene, xylenes, sulfur dioxide, dimethyl sulfide, carbonyl sulfide, formic acid, and hydrogen cyanide. Given the complex mixture of interfering gas signatures at this site, the detection limits of this open-path system were insufficient for flux calculations.

2.3 Method of deriving gradient flux

2.3.1 Method of deriving eddy diffusivity of gas (Kc)

Gradient flux estimates are derived from the vertical gradient of mole fractions and the associated turbulence, given by

(1) F c = - K c c z ,

where Fc is the flux for a pollutant gas c, and cz is the vertical gradient of mole fractions of gas c. Kc is the eddy diffusivity of the gas c, a transfer coefficient characterizing turbulent transport (Monin and Obukhov, 1954) and relating the vertical gradient of gas c to its flux. To calculate the gradient flux of pollutants measured by the OP-FTIR system, Kc has to be calculated first. In this study, Kc is calculated with a variation on the modified Bowen ratio method (Meyers et al., 1996; Bolinius et al., 2016), which was described in detail in You et al. (2020). The key steps are repeated here. With the EC flux of CH4 measured at 18 m and mole fraction gradient measured at 8 and 32 m by CRDS, Kcwas calculated using Eq. (1). The limitation of using this Kc is that when the vertical gradient cz is close to zero, both Fc and Kc are poorly defined. A much better behaved eddy diffusivity is that for momentum, calculated by combining the wind profile with the momentum flux (Eq. 2). To construct a continuous time series of Kc for CH4, we used Km and the relationship between Km and Kc in Eq. (3).


Kc can thus be calculated from Km if the Schmidt number Sc is known. You et al. (2020) showed the details of two approaches to calculate Sc. In this paper, we took the result of stability-corrected variable Sc from the second approach in You et al. (2020).

2.3.2 Gradient flux based on OP-FTIR measurements

Kc is a function of height above the surface (Monin and Obukhov, 1954). The Kc calculated above applies to a height of 18 m and needs to be adjusted to the OP-FTIR path heights. In this study, the vertical profiles of the CH4 mole fractions varied over time and mostly showed linear vertical profiles when the wind was from the pond (Sect. S1.2). In the following calculation, the vertical profiles of CH4 and other gases are considered linear over the entire project. Therefore, the representative average height of the FTIR top path is taken as the height of the middle point (at 12 m). KcFTIR for gradient flux calculated from the top and bottom FTIR paths has been adjusted linearly based on the Kc(8 m,32 m) calculated from point measurements at 8 and 32 m on the tower:

(4) K c FTIR K c (8 m,32 m) = K m FTIR S c K m (8 m,32 m) S c = K m FTIR K m (8 m,32 m) .

Km is a function φm, which in turn is a function of stability (z/L), where z is the height above ground and L the Obukhov length:


Here u is the friction velocity. From Eqs. (4)–(6), KmFTIR/Km(8 m,32 m) is calculated for each half-hour period. The FTIR gradient flux is then calculated by multiplying this KcFTIR by the mole fraction gradient between top and bottom path of the FTIR in Eq. (1). The difference between assuming linear and logarithmic vertical profiles of the mole fractions is discussed in the Sect. S1.2. The logarithmic vertical profile assumption resulted in an average flux which is 19 % greater than the average gradient flux calculated with linear vertical profiles.

In addition to calculating gradient fluxes by using the CH4 mole fraction gradient between top and bottom paths, gradient fluxes of CH4 were also calculated by using the mole fraction gradient between middle and bottom paths. Results show that the average gradient fluxes with middle-bottom paths gradient are 29 % less than that with top-bottom paths (Sect. S1.2, Table S1). These results suggest the gradient fluxes in this study may include uncertainties from the vertical profile of Kc, which has not been studied much in the past. Lastly, the approach developed by Flesch et al. (2016) was applied and compared to the approach shown above. Fluxes were within 70 %. Kc values for NH3 and total alkanes are assumed to be the same as Kc for CH4.

2.3.3 Inverse-dispersion fluxes

Inverse-dispersion models (IDMs) are a useful method for deriving emission rate estimates based on line-integrated or point mole fraction measurements downwind of a defined source. IDMs require inputs including the mole fraction measurements, the surface turbulence statistics between the measurements and the source, and the background mole fraction (mole fraction upwind of the source) of the pollutant. In this study, we used WindTrax 2.0 (Thunder Beach Scientific,, last access: 26 January 2021; Flesch et al., 1995, 2004), which is based on a backward Lagrangian stochastic (bLS) model. Details on IDM calculations and resulting CH4 fluxes were presented in You et al. (2020). Key information is repeated here. The emission rate is derived through

(7) Q = C - C b C Q sim ,

where C is the mole fraction of the pollutant measured at a given location, Cb is the background mole fraction and (C/Q)sim is the simulated ratio (calculated by the bLS model) of the pollutant mole fraction at the site to the emission rate at the specified source. In this study, the inputs of surface turbulence statistics are u and L calculated from the ultrasonic anemometer measurements at 8 m, ambient temperature measured at 8 m, and the mean wind direction measured by the propeller anemometer at 4 m. All surface turbulence statistics inputs were half-hour averages. Half-hour periods when u< 0.15 m s−1 were excluded (Flesch et al., 2004). IDM fluxes of NH3 and total alkane are shown in this work for the first time for comparison with gradient fluxes. Meteorological inputs for calculating inverse-dispersion fluxes of NH3 and total alkanes are the same as for CH4 fluxes. CH4, NH3, and total alkanes mole fraction inputs were mole fractions from the bottom path of OP-FTIR.

3 Results and discussion

3.1 Meteorological conditions

The measurement site including the OP-FTIR was at the south shore of the pond (Fig. 1); therefore, the northern wind (wind direction (WD)  286 or WD  76 ) was defined as the wind coming from the pond (You et al., 2020). The wind came from the north for about 22 % of the entire measurement period (Fig. S2 in You et al., 2020). There was no significant diurnal variation in wind direction during the study period. Detailed ambient temperature, water surface temperature, wind speed, and other meteorological parameters can be found in You et al. (2020). As discussed in You et al. (2020), the warm pond surface resulted in continuing convective turbulence at night, resulting in continuing transport of pollutants from the pond into the atmosphere without significant diurnal variation.

Gradient and IDM fluxes for NH3 and total alkanes are averages for half-hour periods when the wind came from the pond. The half-hour fluxes were binned into 16 wind direction sectors, and the area-weighted averages of fluxes from the pond were calculated as described in You et al. (2020).

3.2 Methane

Path-integrated mole fractions and associated gradient fluxes of CH4 from OP-FTIR are presented here to test if the gradient fluxes derived from the mole fractions with OP-FTIR are comparable to CH4 fluxes from EC and IDM methods (You et al., 2020). The area-weighted flux statistics from different methods are summarized in Table 2. The path-integrated measurement from the FTIR bottom path clearly indicates that the CH4 mole fraction was elevated when the wind was from the pond direction, while it was steady near 2 ppm when the wind was from other directions (Figs. S2 and S3). In addition, a clear vertical gradient (Fig. S3), with mole fractions along the bottom path on the order of 0.5 to 1 ppm higher than mole fractions from the top path, identified the pond as the CH4 source. The fact that the CH4 mole fraction increased when the wind was from the pond direction, and decreased with height, clearly points to the pond as the dominant local source.

Table 2Summary of fluxes from OP-FTIR measurements. Results are area-weight-averaged fluxes from the pond.

 Errors with the mean fluxes are calculated with a top-down approach: the average of standard deviations of fluxes from five periods when the fluxes displayed high stationarity.

Download Print Version | Download XLSX

For comparison, vertical profiles of the CH4 mole fraction by point measurements on the nearby tower are given in Sect. S1.2. A linear vertical extrapolation of the profiles to the point where the mole fraction reaches 2.0 ppm (background levels) indicated a median plume height of 64 m (Figs. S4 and S5).

Figure 2Gradient flux of CH4 from FTIR binned by wind direction in 22.5 bins. Lower and upper bounds of the box plot are the 25th and 75th percentile; the line in the box marks the median and the black square labels the mean; the whiskers label the 10th and 90th percentile.


The gradient flux derived from the OP-FTIR shows that the flux was minimal when the wind was from other directions, except for the sector centered at 270 (Fig. 2), which represented a mix of pond and shoreline influences. The average and interquartile ranges of fluxes in wind direction sectors centered at 315, 337.5, and 0 are comparable. This gradient flux result is consistent with the EC fluxes measured on the adjacent flux tower (You et al., 2020), and these results also suggest that the pond is the main source of measured CH4 fluxes. However, the sector centered at 292.5 shows average flux 22 % and 73 % greater than sectors centered at 315 and 337.5. This is different from the EC fluxes, which showed closer agreement between the 292.5, 315, 337.5, and 0 sectors (Fig. 7 in You et al., 2020). The footprint of the EC fluxes measured on the adjacent flux tower at 18 m was calculated by using the Flux Footprint Prediction (FFP) model in Kljun et al. (2015), and results showed the 80 % contribution distance was typically within 1 km, which is closer to the main site than the north edge of pond liquid surface (Fig. 1; You et al., 2020). The outfall was about 1.4 km from the main site. The discrepancy suggests that the footprint of the gradient method incorporated emissions from the outfall more clearly than the smaller footprint of the EC method.

FTIR CH4 gradient fluxes and EC fluxes showed a linear correlation, but on average the gradient fluxes were lower than the EC fluxes by 44 % (r2= 0.63) (Fig. S7). In agreement with EC, the gradient flux showed no significant diurnal variations when the wind was from the pond (Fig. S8, with a relative standard deviation of 38 %). To investigate the difference between CH4 gradient fluxes derived from FTIR and EC fluxes, the latter were examined in relation to meteorological conditions, similar to the analysis presented in You et al. (2020). The analysis in You et al. (2020) showed weak correlation between EC flux and friction velocity (u) or wind speed, consistent with the observation of the weak correlation between gradient flux and wind speed (Fig. S5 in You et al., 2020). As described in Sect. S1.2 (Figs. S4, S5, and S6), CH4 vertical profiles were closer to linear when the wind speed was less than 6 m s−1 and were more logarithmic with wind speed greater than 6 m s−1. The sector centered at 292.5 was often associated with wind speeds greater than 6 m s−1 (21 % of the time). The approximation of linear vertical profile could have underestimated the CH4 flux with periods of high winds by 19 %. This weak dependence of gradient flux on wind speed may reflect the uncertainty in the vertical profile of Kc in the calculation of gradient flux.

Model calculations by Horst (1999) showed that the estimated footprint of a gradient flux measurement at the geometric mean height of the gradient is similar to the footprint of EC flux at that same height for homogeneous upwind area sources. However, mole fraction footprints are significantly larger than perturbation (flux) footprints (Schmid, 1994), and some CH4 sources on the far shore (e.g., the outfall) may have contributed to the upper-path CH4 mole fraction. This decreased the vertical gradient difference and thus the derived flux relative to the EC flux with its smaller footprint, since the latter is more likely to represent water surface emissions only. As an approximate estimation, the footprint of the path-integrated mole fraction of the top path is about 2.3 km (23 m × 100; Flesch et al., 2016), and this covers the whole pond including the north edge (Fig. 1).

Background mole fractions, upwind of the source under investigation, must be provided for the bLS calculations of CH4 fluxes. We quantified these using two methods. First, the background mole fraction was determined with the FTIR measurements at the south of the pond, as follows: for most of the days, it was taken as the minimum CH4 mole fraction from the FTIR bottom path on each day while the wind direction was between 180 and 240. On 7 and 30 August, there was no half-hour period when the wind was from this sector, and the background mole fraction was chosen as the minimum mole fraction for the day. For 1 August, there was also no half-hour period for this sector, and the minimum of the day was 2.40 ppm, significantly greater than the minimum mole fraction of other days. Therefore, the background mole fraction of the previous day, 1.92 ppm, was used for 1 August.

Alberta Environment and Parks (AEP) conducted OP-FTIR measurements (RAM2000 G2; KASSAY FSI, ITT Corp., Mohrsville, PA, USA) at the north side of the pond (Fig. 1), quantifying CH4 to be used as the second estimate of background mole fractions. For most of the days, the half-hour averaged mole fractions were directly used as the background mole fractions. From 3 August 22:00 to 4 August 13:30, 6 August 08:00 to 17:00, and 23 August 01:30 to 02:00 UTC there were no data from AEP, so background mole fractions for these periods were picked as the interpolation of mole fractions before and after this period. In this approach, the bLS flux results can be negative when the AEP mole fraction is greater than the mole fraction from the measurements at the south shore, possibly due to influences by other emission sources in the surrounding area, gas diffusion under low wind speeds, plume inertia when wind direction changes suddenly, or instrument mismatch differences.

CH4 IDM fluxes with background determined from the first approach (using measurements at the south of the pond) agreed with IDM fluxes with background determined from the second approach (using measurements at the north of the pond), with a linear regression r2 of 0.92 and a slope of 0.90; there was a 20 % difference between average fluxes from the two approaches (Fig. S10; Table S1). These results confirm that the CH4 flux estimate from this inverse-dispersion approach is consistent and that the first approach to determining backgrounds is appropriate. In the following results and discussion, IDM fluxes with background mole fractions from the first method are used.

IDM and EC flux showed moderate comparison (slope = 0.69, r= 0.62; Fig. S7 in You et al., 2020), although the mean IDM fluxes are 30 % smaller than the EC flux.

3.3 NH3

The mole fraction of NH3 was elevated when the wind was from the pond but was mainly below 5 ppb when the wind was from the south (Fig. 3a). NH3 gradient fluxes were significant when the wind came from the pond direction (Fig. 4a).

Figure 3Normalized rose plot of NH3 (a) and total alkane (b) mole fractions from FTIR bottom path. Colors represent the mole fraction in part per billion. The length of each colored segment presents the time fractions of that mole fraction range in each direction bin. The radius of the black open sectors indicates the frequency of wind in each direction bin; the angle represents wind direction – straight up is north and straight left is west.


Figure 4Gradient flux of NH3 (a) and total alkane (b) from FTIR top-bottom path binned by wind direction in 22.5 bins. Lower and upper bounds of the box plot are the 25th and 75th percentile; the line in the box marks the median and the black square labels the mean; the whiskers label the 10th and 90th percentile.


Figure 5(a) NH3 gradient flux compared to CH4 gradient flux; (b) NH3 IDM flux compared to CH4 IDM flux.


The time series of mole fraction vertical gradient of NH3 and CH4 were similar (Figs. S3 and S12). The NH3 gradient flux and CH4 gradient flux showed good correlation (r2= 0.7, Fig. 5a). The diurnal variation in the NH3 gradient flux (relative standard deviation = 72 %) was stronger than for the CH4 gradient flux (relative standard deviation = 36 %), with greater fluxes from 13:00 to 18:00 MDT (mountain daylight savings time) (Fig. S11a). Previous studies showed tailings pond waters contained elevated NH3 concentrations, which makes them potential sources of NH3 to the atmosphere (Allen, 2008; Risacher et al., 2018). NH3 in the pond is mainly produced through nitrate and/or nitrite reduction during microbial activities (Barton and Fauque, 2009; Collins et al., 2016). In addition, some of these nitrate and/or nitrite reduction microbes may also produce reduced sulfur (Barton and Fauque, 2009). The water sample collected from Pond 2/3 on August 2017 during this study was alkaline (pH = 8.0 ± 0.5), which also supports the emission of NH3.

The average flux in the sector centered at 292.5 was 33 % and 57 % more than the average flux in the sectors centered at 315 and 337.5, but the median fluxes in these three sectors are within 14 %. These suggest the high average flux in 281–304 is skewed by big spikes which are associated with the outfall (with wind directions in 281–304), but the majority of NH3 fluxes in the 281–304 wind sector correlated well with CH4 fluxes which were less affected by the outfall. Although hydrotreating processes in upgraders remove most of the sulfur and nitrogen from the bitumen residue, a small amount of NH3 might still be carried with the processed water and tailings (Bytnerowicz et al., 2010) and transported with the outfall liquid into the pond. The negative fluxes observed for the 67.5 sector may be due to elevated NH3 plumes originating from the upgrader facility 3 km upwind in this direction, resulting in a negative gradient and thus deposition to the pond under some circumstances.

IDM fluxes of NH3 were calculated the same way as CH4 and show a weak correlation with CH4 IDM fluxes (Fig. 5b). The NH3 background mole fraction was based on the mean daily minimum, approximately 1 ppb (Fig. 3). Vertical profiles of NH3 mole fraction (Fig. S12) with northerly wind also show roughly linear profiles similar to CH4. Profiles of sectors centered at 292.5 and 315 are linear. Therefore, the outfall on average did not significantly contribute to the NH3 profile; i.e., the pond surface was the main source of NH3. NH3 fluxes from the gradient method were significantly less than from the IDM method. This difference is mostly due to the input background NH3 mole fraction. The background NH3 mole fraction was not measured and could be greater than 1 ppb if there was any source to the north of the pond. Assuming the NH3 gradient flux as a reference, different backgrounds were tested in IDM to match the mean gradient flux. A background of 7 ppb of NH3 was required to close the gap between gradient and IDM fluxes. This seems large but cannot be verified, since there was no ground level measurement of NH3 near the north of the pond. This illustrates the advantage of using either EC or gradient flux measurements, which are based on mole fraction fluctuations or gradients at a single location and are therefore independent of upwind background mole fraction.

3.4 Total alkane

Total alkane derived from FTIR spectra used butane and octane as two surrogates in this study, following the method in EPA OTM 10 (Thoma et al., 2010). Results only reflect alkanes which have similar absorption features to butane and octane and cannot accurately represent other VOCs. Total alkane fluxes from the pond were evident (Fig. 4b). A comparison to CH4 fluxes showed only a weak correlation (r2= 0.3, Fig. S13), unlike the correlation between NH3 and CH4 (Fig. 5a). This difference can be explained by sources of alkane and CH4 at this site. Figure 4b shows the average flux from the sector centered at 292.5 is 2.4 g m−2 d−1, which is 2.09 and 2.29 times of the average fluxes from the sectors centered at 315 and 337.5. Figure 2 shows that the average fluxes of CH4 from the sector centered at 292.5 is 1.22 and 1.72 times of the average fluxes from sectors centered at 315 and 337.5. These results indicate that there was an enhanced contribution (27 %) from around 281 to 304 to total alkane flux measured at the site but not to the observed CH4 flux. This enhancement of alkane flux is likely due to the outfall, which was at the edge of the pond, 1.9 km from the site at 295. The liquid mixture flowing into the pond contained naphthenic solvent, which includes a mixture of alkanes, alkenes, and aromatic hydrocarbons, and since the outfall is at a temperature of approximately 33 C, enhanced evaporation of volatile components can be expected from this area. The outfall also introduces some mechanical mixing in the upper layers of the pond water, which may contribute to elevated emission rates. The diurnal variation of total alkane gradient flux when the wind came from the pond direction was also weak (the standard deviation of average fluxes at each hour is comparable to the interquartile ranges, Fig. S14). The vertical profiles of total alkane mole fraction with northeastern winds were vertically invariant (Fig. S15). With northwestern winds, profiles showed a decreased mole fraction from the bottom to middle path and a minimal decrease or even increase from the middle to the top path. These total alkane vertical profiles with northern wind are different from CH4 or NH3 profiles with the same periods, suggesting there were additional sources other than the pond surface to the measured total alkane flux, such as the outfall and industrial activities upwind. IDM flux of total alkanes agreed moderately with gradient flux, with a difference of 21 % (26 %) in the pond average (median) flux.

3.5 Methanol (CH3OH)

Unlike pollutants studied above, the CH3OH mole fraction did not show significant enhancement when wind was from the pond (Fig. S16), suggesting the pond was not significantly contributing to CH3OH compared to potential sources surrounding the pond. The calculation of the gradient flux of CH3OH was attempted in the same way as for the other pollutants, and the flux was on the order of 1 mg m−2 d−1 but with an uncertainty that made it not statistically different from zero. The lifetime of CH3OH is around 10 d (Simpson et al., 2011; Shephard et al., 2015), and the main source is vegetation (Millet et al., 2008). The mole fractions observed at the site were consistent with satellite measurements representative of the general oil sands region (Shephard et al., 2015) and with an airborne study of VOCs (Simpson et al., 2010). In addition to emissions from vegetation, CH3OH observed in the oil sands region could be due to transport from biomass burning (Simpson et al., 2011; Bari and Kindzierski, 2018) and local traffic (Rogers et al., 2006; You et al., 2017).

3.6 Comparison of calculated fluxes to reported emissions and approaches

As a test of the robustness of these results, fluxes of CH4, NH3, and total alkane are also calculated from the slant path method described by Flesch et al. (2016). Calculation inputs and results are summarized in the Sect. S5. Compared to results from our flux-gradient approach, CH4, NH3, and total alkane fluxes with the slant path flux-gradient method were reasonably comparable (24 %, 25 %, and 30 % lower, respectively; Table S1). The difference in fluxes from the two approaches could be due to differences in assumptions regarding the vertical profiles. In our gradient flux approach, we used linear vertical mole fraction profiles of pollutants to calculate the difference of mean heights between the two paths. In Flesch et al. (2016) the gradient flux for pairs of points along the two paths was integrated along the entire path length assuming the flux was uniform horizontally. In that approach, the dependence of Kc on height is incorporated explicitly, assuming a logarithmic wind profile including a stability correction. In our gradient flux approach, Kc is derived from a measured and stability-corrected Km and therefore did not require a wind profile shape assumption.

To facilitate a transparent comparison of the emission results from this study to reported facility-wide emissions of Suncor, we present emission rates based on simple extrapolation of the measured August emissions to the whole year. Other possible seasonal emission profiles have been reported (Small et al., 2015); using these to convert from August emissions to annual average values would for example result in a scale factor of 0.92 (mass transfer model), 0.64 (mass transfer model adjusted for ice cover), or 0.42 (thawing degree-day model) (Cumulative Environmental Management Association, 2011). The seasonally invariant total emission estimates for NH3 and total alkanes from Pond 2/3 to the air were 42 and 881 t yr−1 in 2017, converted from gradient fluxes results in g m−2 d−1. However, NH3 emissions from Pond 2/3 have not been reported in the past, because NH3 was not being measured as part of compliance flux chamber monitoring. Therefore, the facility-wide NH3 emissions reported to the Government of Canada National Pollutant Release Inventory (NPRI) (0.82 t yr−1) in 2017 did not include fugitive emissions from tailings ponds. The solvents entering this tailings pond are naphtha additives with octane, nonane, and heptane as the biggest contributors. Li et al. (2017) quantified alkane emissions (including n-alkanes, branched alkanes, and cycloalkanes) as 36.2 t d−1 from the whole facility using airborne measurements in 2013, which is 73.9 % of 49 t d−1 total VOC emission. Alkanes account for 54 % of total VOC emitted from Pond 2/3 in 2017 (internal communication with Samar Moussa). If we use reported facility-wide VOC annual emissions in 2017 (17 242 t yr−1, NPRI) to estimate facility-wide total alkane annual emission, we obtain 17 242 × 53 % = 9138 t yr−1, and total alkane emissions from Pond 2/3 contribute 9.6 % to facility-wide emissions for 2017. The fugitive NH3 emissions from Pond 2/3 in this study were 42 t yr−1, a number that is 51 times the process-related emission number reported for the facility to NPRI. Negligible volume of ammonia may be carried over to the pond through the naphtha recovery unit and froth treatment tailings (FTT) line, and it is believed that the ammonia is mainly generated from the biogenic activities in the mature fine tailings (MFT) layer of the pond. The majority of H2S, NH3, and CH4 emissions are related to microbiological activities as evidenced in this study.

4 Conclusions and implications

We have shown that OP-FTIR is an effective method to quantify mole fractions and vertical gradients of CH4, NH3, and total alkanes continuously and simultaneously for an area source such as a tailings pond. Benefits are the integration of mole fractions over long path lengths, thus providing a spatially representative average, and the avoidance of sample line issues that can be serious problems for sticky gases such as NH3. Results from the gradient flux method and IDM calculations suggest OP-FTIR is a useful tool for deriving emission rates of CH4, NH3, and total alkane from this type of fugitive area source. For the two approaches for determining background mole fractions of CH4 with the IDM method, i.e., upwind background measurement vs. local background estimation, the area-weight-averaged fluxes of CH4 were within 20 %. FTIR CH4 gradient fluxes and EC fluxes showed a linear correlation, but on average the gradient fluxes were lower than the EC fluxes by 44 %. CH4 IDM and EC flux showed moderate comparison, and the mean IDM fluxes are 30 % smaller than the EC mean flux. NH3 gradient flux and IDM flux showed a difference of more than 50 %, which suggests that there may have been sources of NH3 upwind (north) of the pond that were not captured by assuming that southern and northern background NH3 were similar, thus illustrating a limitation of the IDM method. CH4, NH3, and total alkane fluxes were also calculated using the slant path flux-gradient method (Flesch et al., 2016), to compare to the our flux-gradient approach, and results from the slant path method were reasonably comparable (24 %, 25 %, and 30 % lower, respectively).

The NH3 emissions results in this study are the first to quantify NH3 fugitive fluxes from a tailings pond and clearly showed that Pond 2/3 is a significant source of NH3, most likely through microbial activities in the pond. This suggests that at least some tailings ponds in the oil sands could be significant sources of NH3, compared to process-related facility emissions. Further measurements of NH3 emissions from tailings ponds are recommended to elucidate our understanding of the mechanisms behind NH3 emissions and to improve the total facility emission estimates reported to NPRI.

Total alkane gradient fluxes from OP-FTIR measurements clearly showed that the pond is a significant source of total alkane, and annual alkane emissions extrapolated from gradient flux method represented 9.6 % of facility emissions. The outfall area contributed significantly (27 %) to pond alkane emissions, showing spatial variability of alkane emissions from the pond. Observed CH3OH mole fractions show that the pond was not likely a significant source of CH3OH. This study demonstrated the applicability of OP-FTIR combined with gradient flux or inverse-dispersion methods for determining emission fluxes of multiple gases simultaneously, with high temporal resolution and comprehensive spatial coverage.

Data availability

The supplement related to this article is available online at:

Author contributions

YY and RMS wrote the manuscript; SGM, LZ, LF, and JB contributed data and comments.

Competing interests

James Beck is an employee of Suncor Energy. The other authors declare that they have no conflict of interest.


The authors thank the technical team of Andrew Sheppard, Roman Tiuliugenev, Raymon Atienza, and Raj Santhaneswaran for their invaluable contributions throughout; Julie Narayan for spatial analysis; Stewart Cober for management; and Stoyka Netcheva for home base logistical support. We thank Suncor and its project team (Dan Burt et al.), AECOM (April Kliachik, Peter Tkalec) and SGS (Nathan Grey, Ardan Ross) for site logistics support. This work was partially funded under the Oil Sands Monitoring Program and is a contribution to the program but does not necessarily reflect the position of the program. We also acknowledge funding from the Program for Energy Research and Development (Natural Resources Canada) and from the Climate Change and Air Pollution Program (ECCC).

Financial support

This research has been supported by the Oil Sands Monitoring Program, the Program for Energy Research and Development (Natural Resources Canada), and the Climate Change and Air Pollution Program (ECCC).

Review statement

This paper was edited by Daniela Famulari and reviewed by two anonymous referees.


Akagi, S. K., Yokelson, R. J., Burling, I. R., Meinardi, S., Simpson, I., Blake, D. R., McMeeking, G. R., Sullivan, A., Lee, T., Kreidenweis, S., Urbanski, S., Reardon, J., Griffith, D. W. T., Johnson, T. J., and Weise, D. R.: Measurements of reactive trace gases and variable O3 formation rates in some South Carolina biomass burning plumes, Atmos. Chem. Phys., 13, 1141–1165,, 2013. 

Akagi, S. K., Burling, I. R., Mendoza, A., Johnson, T. J., Cameron, M., Griffith, D. W. T., Paton-Walsh, C., Weise, D. R., Reardon, J., and Yokelson, R. J.: Field measurements of trace gases emitted by prescribed fires in southeastern US pine forests using an open-path FTIR system, Atmos. Chem. Phys., 14, 199–215,, 2014. 

Alberta Environment and Parks: Total Area of the Oil Sands Tailings Ponds over Time, available at: (last access: 22 September 2019), 2016. 

Allen, E. W.: Process water treatment in Canada's oil sands industry: I. Target pollutants and treatment objectives, J. Environ. Eng. Sci., 7, 123–138,, 2008. 

Bai, M., Suter, H., Lam, S. K., Sun, J., and Chen, D.: Use of open-path FTIR and inverse dispersion technique to quantify gaseous nitrogen loss from an intensive vegetable production site, Atmos. Environ., 94, 687–691,, 2014. 

Bai, M., Suter, H., Lam, S. K., Davies, R., Flesch, T. K., and Chen, D.: Gaseous emissions from an intensive vegetable farm measured with slant-path FTIR technique, Agr. Forest Meterol., 258, 50–55,, 2018. 

Bai, M., Suter, H., Lam, S. K., Flesch, T. K., and Chen, D.: Comparison of slant open-path flux gradient and static closed chamber techniques to measure soil N2O emissions, Atmos. Meas. Tech., 12, 1095–1102,, 2019. 

Bari, M. A. and Kindzierski, W. B.: Ambient volatile organic compounds (VOCs) in communities of the Athabasca oil sands region: Sources and screening health risk assessment, Environ. Pollut., 235, 602–614,, 2018. 

Barton, L. L. and Fauque, G. D.: Chapter 2: Biochemistry, physiology and biotechnology of sulfate-reducing bacteria, Adv. Appl. Microbiol., 68, 41–98,, 2009. 

Bolinius, D. J., Jahnke, A., and MacLeod, M.: Comparison of eddy covariance and modified Bowen ratio methods for measuring gas fluxes and implications for measuring fluxes of persistent organic pollutants, Atmos. Chem. Phys., 16, 5315–5322,, 2016. 

Bradley, K. S., Brooks, K. B., Hubbard, L. K., Popp, P. J., and Stedman, D. H.: Motor vehicle fleet emissions by OP-FTIR, Environ. Sci. Tech., 34, 897–899,, 2000. 

Burling, I. R., Yokelson, R. J., Griffith, D. W. T., Johnson, T. J., Veres, P., Roberts, J. M., Warneke, C., Urbanski, S. P., Reardon, J., Weise, D. R., Hao, W. M., and de Gouw, J.: Laboratory measurements of trace gas emissions from biomass burning of fuel types from the southeastern and southwestern United States, Atmos. Chem. Phys., 10, 11115–11130,, 2010. 

Bytnerowicz, A., Fraczek, W., Schilling, S., and Alexander, D.: Spatial and temporal distribution of ambient nitric acid and ammonia in the Athabasca Oil Sands Region, Alberta, J. Limnol., 69, 11–21,, 2010. 

Collins, C. E. V., Foght, J. M., and Siddique, T.: Co-occurrence of methanogenesis and N2 fixation in oil sands tailings, Sci. Total Environ., 565, 306–312,, 2016. 

Cumulative Environmental Management Association: Protocol for Updating and Preparing a Modelling Emission Inventory, available at: (last access: 26 February 2020), 2011. 

Environment and Climate Change Canada: Emissions from tailings ponds to the atmosphere, oil sands region, available at: oil-sands-region/emissions-from-tailings-ponds-to-the-atmosphere-oil-sands-region/, last access: 26 January 2021. 

Field, R. A., Soltis, J., McCarthy, M. C., Murphy, S., and Montague, D. C.: Influence of oil and gas field operations on spatial and temporal distributions of atmospheric non-methane hydrocarbons and their effect on ozone formation in winter, Atmos. Chem. Phys., 15, 3527–3542,, 2015. 

Flesch, T. K., Wilson, J. D., and Yee, E.: Backward-time Lagrangian stochastic dispersion models and their application to estimate gaseous emissions, J. Appl. Meteorol., 34, 1320–1332,<1320:BTLSDM>2.0.CO;2, 1995. 

Flesch, T. K., Wilson, J. D., Harper, L. A., Crenna, B. P., and Sharpe, R. R.: Deducing ground-to-air emissions from observed trace gas concentrations: A field trial, J. Appl. Meteorol., 43, 487–502,<0487:DGEFOT>2.0.CO;2, 2004. 

Flesch, T. K., Wilson, J., Harper, L., and Crenna, B.: Estimating gas emissions from a farm with an inverse-dispersion technique, Atmos. Environ., 39, 4863–4874,, 2005. 

Flesch, T. K., Baron, V. S., Wilson, J. D., Griffith, D. W. T., Basarab, J. A., and Carlson, P. J.: Agricultural gas emissions during the spring thaw: Applying a new measurement technique, Agr. Forest Meterol., 221, 111–121,, 2016. 

Foght, J. M., Gieg, L. M., and Siddique, T.: The microbiology of oil sands tailings: Past, present, future, FEMS Microbiol. Ecol., 93, fix034,, 2017. 

Galarneau, E., Hollebone, B. P., Yang, Z., and Schuster, J.: Preliminary measurement-based estimates of PAH emissions from oil sands tailings ponds, Atmos. Environ., 97, 332–335,, 2014. 

Goode, J. G., Yokelson, R. J., Susott, R. A., and Ward, D. E.: Trace gas emissions from laboratory biomass fires measured by open-path Fourier transform infrared spectroscopy: Fires in grass and surface fuels, J. Geophys. Res., 104, 21237–21245,, 1999. 

Government of Canada National Pollutant Release Inventory:, last access: 7 January 2020. 

Griffith, D. W. T., Mankin, W. G., Coffey, M. T., Ward, D. E., and Riebau, A.: “FTIR remote sensing of biomass burning emissions of CO2, CO, CH4, CH2O, NO, NO2, NH3, and N2O.” Global biomass burning: atmospheric, alimate, and biospheric implications, MIT Press, Cambridge, MA, USA, 1991. 

Grutter, M., Flores, E., Basaldud, R., and Ruiz-Suarez, L. G.: Open-path FTIR spectroscopic studies of the trace gases over Mexico City, Atmos. Ocean. Opt., 16, 232–236, 2003. 

Horrocks, L., Burton, M., Francis, P., and Oppenheimer, C.: Stable gas plume composition measured by OP-FTIR spectroscopy at Masaya Volcano, Nicaragua, 1998–1999, Geophys. Res. Lett,, 26, 3497–3500,, 1999. 

Horst, T. W.: The footprint for estimation of atmosphere-surface exchange fluxes by profile techniques, Bound.-Lay. Meteorol., 90, 171–188,, 1999. 

Hu, N., Flesch, T. K., Wilson, J. D., Baron, V. S., and Basarab, J. A.: Refining an inverse dispersion method to quantify gas sources on rolling terrain, Agr. Forest Meterol., 225, 1–7,, 2016. 

Johnson, T. J., Profeta, L. T. M., Sams, R. L., Griffith, D. W. T., and Yokelson, R. L.: An infrared spectral database for detection of gases emitted by biomass burning, Vib. Spectrosc., 53, 97–102,, 2010. 

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

Kroll, J. and Seinfeld, J. H.: Chemistry of secondary organic aerosol: Formation and evolution of low-volatility organics in the atmosphere, Atmos. Environ., 42, 3593–3624,, 2008. 

Kürten, A., Bianchi, F., Almeida, J., Kupiainen-Määttä, O., Dunne, E. M., Duplissy, J., Williamson, C., Barmet, P., Breitenlechner, M., Dommen, J., Donahue, N. M., Flagan, R. C., Franchin, A., Gordon, H., Hakala, J., Hansel, A., Heinritzi, M., Ickes, L., Jokinen, T., Kangasluoma, J., Kim, J., Kirkby, J., Kupc, A., Lehtipalo, K., Leiminger, M., Makhmutov, V., Onnela, A., Ortega, I. K., Petäjä, T., Praplan, A. P., Riccobono, F., Rissanen, M. P., Rondo, L., Schnitzhofer, R., Schobesberger, S., Smith, J. N., Steiner, G., Stozhkov, Y., Tomé, A., Tröstl, J., Tsagkogeorgas, G., Wagner, P. E., Wimmer, D., Ye, P., Baltensperger, U., Carslaw, K., Kulmala, M., and Curtius, J.: Experimental particle formation rates spanning tropospheric sulfuric acid and ammonia abundances, ion production rates, and temperatures, J. Geophys. Res., 121, 12377–12400,, 2016. 

Li, S. M., Leithead, A., Moussa, S. G., Liggio, J., Moran, M. D., Wang, D., Hayden, K., Darlington, A., Gordon, M., Staebler, R., Makar, P. A., Stroud, C. A., McLaren, R., Liu, P. S. K., O'Brien, J., Mittermeier, R. L., Zhang, J., Marson, G., Cober, S. G., Wolde, M., and Wentzell, J. J. B.: Differences between measured and reported volatile organic compound emissions from oil sands facilities in Alberta, Canada, P. Natl. Acad. Sci. USA., 114, E3756–E3765,, 2017. 

Makar, P. A., Akingunola, A., Aherne, J., Cole, A. S., Aklilu, Y.-A., Zhang, J., Wong, I., Hayden, K., Li, S.-M., Kirk, J., Scott, K., Moran, M. D., Robichaud, A., Cathcart, H., Baratzedah, P., Pabla, B., Cheung, P., Zheng, Q., and Jeffries, D. S.: Estimates of exceedances of critical loads for acidifying deposition in Alberta and Saskatchewan, Atmos. Chem. Phys., 18, 9897–9927,, 2018. 

Marshall, T. L., Chaffin, C. T., Hammaker, R. M., and Fateley, W. G.: An introduction to open-path FT-IR atmospheric monitoring, Environ. Sci. Technol., 28, 224–232,, 1994. 

Meyers, T. P., Hall, M. E., Lindberg, S. E., and Kim, K.: Use of the modified bowen-ratio technique to measure fluxes of trace gases, Atmos. Environ., 30, 3321–3329,, 1996. 

Millet, D. B., Jacob, D. J., Custer, T. G., de Gouw, J. A., Goldstein, A. H., Karl, T., Singh, H. B., Sive, B. C., Talbot, R. W., Warneke, C., and Williams, J.: New constraints on terrestrial and oceanic sources of atmospheric methanol, Atmos. Chem. Phys., 8, 6887–6905,, 2008. 

Monin, A. S. and Obukhov, A. M.: Basic laws of turbulent mixing in the surface layer of the atmosphere, Contrib. Geophys. Inst. Acad. Sci. USSR, 24, 163–187, 1954. 

Oppenheimer, C. and Kyle, P. R.: Probing the magma plumbing of Erebus volcano, Antarctica, by open-path FTIR spectroscopy of gas emissions, J. Volcanol. Geotherm. Res., 177, 743–754,, 2008. 

Paton-Walsh, C., Smith, T. E. L., Young, E. L., Griffith, D. W. T., and Guérette, É.-A.: New emission factors for Australian vegetation fires measured using open-path Fourier transform infrared spectroscopy – Part 1: Methods and Australian temperate forest fires, Atmos. Chem. Phys., 14, 11313–11333,, 2014. 

Penner, T. J. and Foght, J. M.: Mature fine tailings from oil sands processing harbour diverse methanogenic communities, Can. J. Microbiol., 56, 459–470,, 2010. 

Risacher, F. F., Morris, P. K., Arriaga, D., Goad, C., Nelson, T. C., Slater, G. F., and Warren, L. A.: The interplay of methane and ammonia as key oxygen consuming constituents in early stage development of Base Mine Lake, the first demonstration oil sands pit lake, Appl. Geochem., 93, 49–59,, 2018. 

Rogers, T. M., Grimsrud, E. P., Herndon, S. C., Jayne, J. T., Kolb, C. E., Allwine, E., Westberg, H., Lamb, B. K., Zavala, M., Molina, L. T., Molina, M. J., and Knighton, W. B.: On-road measurements of volatile organic compounds in the Mexico City metropolitan area using proton transfer reaction mass spectrometry, Int. J. Mass Soectrom., 252, 26–37,, 2006. 

Schäfer, K., Grant, R. H., Emeis, S., Raabe, A., von der Heide, C., and Schmid, H. P.: Areal-averaged trace gas emission rates from long-range open-path measurements in stable boundary layer conditions, Atmos. Meas. Tech., 5, 1571–1583,, 2012. 

Schmid, H. P.: Source areas for scalars and scalar fluxes. Bound.-Lay. Meteorol., 67, 293–318,, 1994. 

Shephard, M. W., McLinden, C. A., Cady-Pereira, K. E., Luo, M., Moussa, S. G., Leithead, A., Liggio, J., Staebler, R. M., Akingunola, A., Makar, P., Lehr, P., Zhang, J., Henze, D. K., Millet, D. B., Bash, J. O., Zhu, L., Wells, K. C., Capps, S. L., Chaliyakunnel, S., Gordon, M., Hayden, K., Brook, J. R., Wolde, M., and Li, S.-M.: Tropospheric Emission Spectrometer (TES) satellite observations of ammonia, methanol, formic acid, and carbon monoxide over the Canadian oil sands: validation and model evaluation, Atmos. Meas. Tech., 8, 5189–5211,, 2015. 

Shonkwiler, K. B. and Ham, J. M.: Ammonia emissions from a beef feedlot: Comparison of inverse modeling techniques using long-path and point measurements of fenceline NH3, Agr. Forest Meterol., 258, 29–42,, 2018. 

Siddique, T., Fedorak, P. M., MacKinnon, M. D., and Foght, J. M.: Metabolism of BTEX and Naphtha Compounds to Methane in Oil Sands Tailings, Environ. Sci. Tech., 41, 2350–2356,, 2007. 

Siddique, T., Penner, T., Semple, K., and Foght, J. M.: Anaerobic biodegradation of longer-chain n-alkanes coupled to methane production in oil sands tailings, Environ. Sci. Tech., 45, 5892–5899,, 2011. 

Siddique, T., Penner, T., Klassen, J., Nesbø, C., and Foght, J. M.: Microbial communities involved in methane production from hydrocarbons in oil sands tailings, Environ. Sci. Tech., 46, 9802–9810,, 2012. 

Simpson, I. J., Blake, N. J., Barletta, B., Diskin, G. S., Fuelberg, H. E., Gorham, K., Huey, L. G., Meinardi, S., Rowland, F. S., Vay, S. A., Weinheimer, A. J., Yang, M., and Blake, D. R.: Characterization of trace gases measured over Alberta oil sands mining operations: 76 speciated C2–C10 volatile organic compounds (VOCs), C2–C10 volatile organic compounds (VOCs), CO2, CH4, CO, NO, NO2, NOy, O3 and SO2, , Atmos. Chem. Phys., 10, 11931–11954,, 2010. 

Simpson, I. J., Akagi, S. K., Barletta, B., Blake, N. J., Choi, Y., Diskin, G. S., Fried, A., Fuelberg, H. E., Meinardi, S., Rowland, F. S., Vay, S. A., Weinheimer, A. J., Wennberg, P. O., Wiebring, P., Wisthaler, A., Yang, M., Yokelson, R. J., and Blake, D. R.: Boreal forest fire emissions in fresh Canadian smoke plumes: C1-C10 volatile organic compounds (VOCs), CO2, CO, NO2, NO, HCN and CH3CN, Atmos. Chem. Phys., 11, 6445–6463,, 2011. 

Small, C. C., Cho, S., Hashisho, Z., and Ulrich, A. C.: Emissions from oil sands tailings ponds: Review of tailings pond parameters and emission estimates, J. Petrol. Sci. Eng., 127, 490–501,, 2015. 

Smith, T. E. L., Paton-Walsh, C., Meyer, C. P., Cook, G. D., Maier, S. W., Russell-Smith, J., Wooster, M. J., and Yates, C. P.: New emission factors for Australian vegetation fires measured using open-path Fourier transform infrared spectroscopy – Part 2: Australian tropical savanna fires, Atmos. Chem. Phys., 14, 11335–11352,, 2014. 

Thoma, E. D., Green, R. B., Hater, G. R., Goldsmith, C. D., Swan, N. D., Chase, M. J., and Hashmonay, R. A.: Development of EPA OTM 10 for landfill applications, J. Environ. Eng., 136, 769–776,, 2010. 

Whaley, C. H., Makar, P. A., Shephard, M. W., Zhang, L., Zhang, J., Zheng, Q., Akingunola, A., Wentworth, G. R., Murphy, J. G., Kharol, S. K., and Cady-Pereira, K. E.: Contributions of natural and anthropogenic sources to ambient ammonia in the Athabasca Oil Sands and north-western Canada, Atmos. Chem. Phys., 18, 2011–2034,, 2018. 

Wiacek, A., Li, L., Tobin, K., and Mitchell, M.: Characterization of trace gas emissions at an intermediate port, Atmos. Chem. Phys., 18, 13787–13812,, 2018. 

Wu, R. T., Chang, S.-Y., Chung, Y. W., Tzou, H. C., and Tso, T.-L.: FTIR remote sensor measurements of air pollutants in the petrochemical industrial park, Proc. SPIE 2552, SPIE's 1995 International Symposium on Optical Science, Engineering, and Instrumentation, 9–14 July 1995, San Diego, USA, 719–727,, 1995. 

Yeh, S., Jordaan, S. M., Brandt, A. R., Turetsky, M. R., Spatari, S., and Keith, D. W.: Land use greenhouse gas emissions from conventional oil production and oil sands, Environ. Sci. Tech., 44, 8766–8772,, 2010. 

Yokelson, R. J.: Emissions of formaldehyde, acetic acid, methanol, and other trace gases from biomass fires in North Carolina measured by airborne Fourier Transform Infrared spectroscopy, J. Geophys. Res.-Atmos., 104, 30109–30125,, 1999. 

Yokelson, R. J., Griffith, D. W. T., and Ward, D. E.: Open-path Fourier Transform Infrared studies of large-scale laboratory biomass fires, J. Geophys. Res.-Atmos., 101, 21067–21080,, 1996. 

Yokelson, R. J., Susott, R., Ward, D. E., Reardon, J., and Griffith, D. W. T.: Emissions from smoldering combustion of biomass measured by open-path Fourier Transform Infrared spectroscopy, J. Geophys. Res.-Atmos., 102, 18865–18877,, 1997. 

Yokelson, R. J., Karl, T., Artaxo, P., Blake, D. R., Christian, T. J., Griffith, D. W. T., Guenther, A., and Hao, W. M.: The Tropical Forest and Fire Emissions Experiment: overview and airborne fire emission factor measurements, Atmos. Chem. Phys., 7, 5175–5196,, 2007. 

Yokelson, R. J., Burling, I. R., Gilman, J. B., Warneke, C., Stockwell, C. E., de Gouw, J., Akagi, S. K., Urbanski, S. P., Veres, P., Roberts, J. M., Kuster, W. C., Reardon, J., Griffith, D. W. T., Johnson, T. J., Hosseini, S., Miller, J. W., Cocker III, D. R., Jung, H., and Weise, D. R.: Coupling field and laboratory measurements to estimate the emission factors of identified and unidentified trace gases for prescribed fires, Atmos. Chem. Phys., 13, 89–116,, 2013.  

You, Y., Staebler, R. M., Moussa, S. G., Su, Y., Munoz, T., Stroud, C., Zhang, J., and Moran, M. D.: Long-path measurements of pollutants and micrometeorology over Highway 401 in Toronto, Atmos. Chem. Phys., 17, 14119–14143,, 2017. 

You, Y., Staebler, R. M., Moussa, S. G., Beck, J., and Mittermeier, R. L.: Methane emissions from an oil sands tailings pond: A quantitative comparison of fluxes derived by different methods, Atmos. Meas. Tech. Discuss. [preprint],, in review, 2020. 

Zhang, L., Cho, S., Hashisho, Z., and Brown, C.: Quantification of fugitive emissions from an oil sands tailings pond by eddy covariance, Fuel, 237, 457–464,, 2019. 

Short summary
Tailings ponds in the Alberta oil sands represent an insufficiently characterized source of fugitive emissions of pollutants to the atmosphere. In this study, a novel approach of using a Fourier transform infrared spectrometer along with measurements of atmospheric turbulence is shown to present a practical, non-intrusive method of quantifying emission rates for ammonia, alkanes, and methane. Results from a 1-month field study are presented and discussed.