Articles | Volume 19, issue 7
Research article
03 Apr 2019
Research article |  | 03 Apr 2019

Estimating vehicle carbon dioxide emissions from Boulder, Colorado, using horizontal path-integrated column measurements

Eleanor M. Waxman, Kevin C. Cossel, Fabrizio Giorgetta, Gar-Wing Truong, William C. Swann, Ian Coddington, and Nathan R. Newbury

We performed 7.5 weeks of path-integrated concentration measurements of CO2, CH4, H2O, and HDO over the city of Boulder, Colorado. An open-path dual-comb spectrometer simultaneously measured time-resolved data across a reference path, located near the mountains to the west of the city, and across an over-city path that intersected two-thirds of the city, including two major commuter arteries. By comparing the measured concentrations over the two paths when the wind is primarily out of the west, we observe daytime CO2 enhancements over the city. Given the warm weather and the measurement footprint, the dominant contribution to the CO2 enhancement is from city vehicle traffic. We use a Gaussian plume model combined with reported city traffic patterns to estimate city emissions of on-road CO2 as (6.2±2.2)×105 metric tons (t) CO2 yr−1 after correcting for non-traffic sources. Within the uncertainty, this value agrees with the city's bottom-up greenhouse gas inventory for the on-road vehicle sector of 4.5×105 t CO2 yr−1. Finally, we discuss experimental modifications that could lead to improved estimates from our path-integrated measurements.

1 Introduction

Measurements of greenhouse gases, especially CO2 and CH4, are critical for monitoring, verification, and reporting as countries and cities work towards decreasing their carbon emissions. Measurements on the city-scale are critical because cities contribute a large fraction of global emissions (Marcotullio et al., 2013; Seto et al., 2014). However, quantification of city greenhouse gas emissions is challenging, especially for CO2 since it has a high background and numerous point and diffuse sources including traffic, power plants, and animal and plant respiration. Emissions of pollutants are typically determined using two methods: (1) a top-down approach using atmospheric measurements over a specific site or area to adjust a prior model and (2) bottom-up inventories that calculate emissions based on sector activity and sector emissions factors. Here we demonstrate a technique for top-down measurements that uses an open-path sensor rather than a point sensor and apply it to measure city-scale CO2 emissions.

Quantification of CO2 fluxes from cities has been determined from eddy covariance flux measurements with a point sensor located on a tower in or near a city (Nemitz et al., 2002; Velasco et al., 2005, 2014; Coutts et al., 2007; Bergeron and Strachan, 2011). However, for a single sensor, the relatively small footprint of the eddy covariance flux measurements limits the utility of this technique for large cities, as do violations of the horizontal homogeneity assumption (Järvi et al., 2018). To overcome this limitation, tower networks of point sensors can measure CO2 at multiple sites within a city and at background sites outside the city (McKain et al., 2012; Lauvaux et al., 2013, 2016; Bréon et al., 2015; Shusterman et al., 2016; Staufer et al., 2016; Mueller et al., 2017; Verhulst et al., 2017; Mitchell et al., 2018; Sargent et al., 2018). To distinguish the small enhancements compared to the large background, these networks often use expensive, high-precision cavity ring-down (CRDS) instruments, resulting in a high cost. The BEACO2N network (Shusterman et al., 2016), on the other hand, has a much lower cost per sensor. It requires calibration for quantitative results, but the high density of the point sensors can provide lower sensitivity to systematics (Turner et al., 2016). All of these methods use an inversion to determine the total emissions and thus rely on well-known priors and high-resolution mesoscale atmospheric models.

More recently, several other approaches have also been applied to city-scale emissions. Aircraft mass balance measurements (White et al., 1976; Ryerson et al., 2001) have been used to determine city emissions (Mays et al., 2009; Heimburger et al., 2017). However, the use of an aircraft is costly and labor intensive and, therefore, not suited to long-term continuous measurements. Column measurements from the Total Carbon Column Observation Network (TCCON) were used to calculate total South Coast Air Basin (SoCAB) CO and CH4 emissions but not CO2 (Wunch et al., 2009). In addition, data from the Orbiting Carbon Observatory satellite (OCO-2) were recently combined with TCCON data to estimate CO2 emissions from the Los Angeles (LA) basin (Hedelius et al., 2018).

As an alternative to these approaches, horizontal, kilometer-scale, open-path instruments could in principle be used to determine CO2 emissions from cities. Such instruments are capable of continuous measurements over a large area with a single instrument (e.g., Wong et al., 2016; Dobler et al., 2017; Coburn et al., 2018). These sensors also have the advantage of being insensitive to small changes in local meteorology and are not subject to the same representation errors as point sensors (Ciais et al., 2010). Several such systems have been deployed. A laser absorption spectrometer system (GreenLITE) has mapped CO2 concentrations over Paris but has not yet quantified emissions (Dobler et al., 2017). The California Laboratory of Atmospheric Remote Sensing Fourier Transform Spectrometer (CLARS-FTS) is a downward-looking slant column Fourier transform spectrometer (FTS) that scans across 28 measurement targets in the LA Basin to measure CO2, CH4, and O2 (Wong et al., 2015). Based on the measured CH4:CO2 ratio and the bottom-up CO2 inventory from the California Air Resources Board, researchers have calculated the LA Basin CH4 emissions (Wong et al., 2016) but not yet the CO2 emissions.

Here we present the quantification of city CO2 emissions using open-path measurements made with a dual frequency comb spectrometer. While dual-comb spectroscopy is a relatively new technique, it has a unique set of attributes that make it attractive for open path measurements (Rieker et al., 2014; Coddington et al., 2016; Waxman et al., 2017; Coburn et al., 2018). Dual-comb spectroscopy (DCS) is a high-resolution, broadband technique spanning hundreds of wave numbers but with a resolution that exceeds even high-end Fourier transform infrared spectrometers (FTIRs) leading to a negligible instrument lineshape (Coddington et al., 2016). This allows for simultaneous measurements of multiple species and path-integrated temperature with low systematic uncertainty and without the need for instrument calibration. Additionally, the eye-safe, high-brightness, single-transverse-mode output of a frequency comb allows for beam paths exceeding 10 km, while the speed and parallelism of the measurement suppress any spectral distortion from the inevitable turbulence-induced power fluctuations over such a path (Rieker et al., 2014; Waxman et al., 2017).

Figure 1 shows the measurement layout for an initial campaign to quantify CO2 emissions from Boulder, Colorado. Here we take the light from a dual comb spectrometer near the edge of the city and simultaneously measure two paths: a reference path that points west-southwest towards the mountains and an over-city path that crosses the city to the northeast, covering the main traffic arteries of the city with sensitivity to traffic emissions. We acquire time-resolved data at 5 min resolution of CO2, CH4, H2O and isotopologues over 7.5 weeks. The dry mole fraction of CO2 shows a diurnal cycle consistent with a morning build-up from traffic followed by a midday decline due to the rising boundary layer. In addition, there is a distinct difference between the weekday and weekend cycles for CO2, consistent with traffic patterns. In order to demonstrate the utility of this method for emissions quantification, we perform a preliminary estimate of the CO2 emissions from traffic. To do this, we filter the data for days when the wind is out of the west and not too strong so that there is a measurable daytime enhancement in CO2 between the reference path and over-city path. Given the weather, beam path location, and observation times, the dominant contribution will be from traffic rather than residential or industrial emissions. We apply a Gaussian plume model to calculate the city emissions based on the expected distributed source (due to traffic) and the path-averaged concentrations. After adjusting for small expected contributions from residential sources and a local utility plant, the measured emission value is scaled to annual city-wide emissions based on city traffic count data. We estimate (6.2±2.2)×105 metric tons (t) CO2 yr−1, compared to the bottom-up City of Boulder inventory estimate of 4.46×105 t CO2 yr−1. Finally, we discuss improvements to this estimate, which could be realized by more advantageous beam paths that sample a larger spatial and temporal fraction of the full city emissions and by a more detailed inventory model.

Figure 1Measurement layout. The two measurement paths are shown by red (reference) and black (over-city) lines. The two weather stations that provided wind speed and direction data are given by the green diamonds. The colored circles are turning movement count (TMC) locations, which are used as a proxy for the traffic source locations. Both color and size represent the number of traffic counts at each location. Dominant wind directions for the campaign overall (aqua) and the test case days (purple for 22 October and blue for 25 October) are given by colored arrows.


2 Experimental data

2.1 DCS measurements

The DCS system was located on the top floor of the National Institute of Standards and Technology (NIST) building in Boulder, Colorado. This instrument has been described previously (Truong et al., 2016; Waxman et al., 2017). The light from the combs is split to generate two combined dual-comb outputs, one of which is transmitted over the reference path and one of which is transmitted over the city path (see Fig. 1). Here, we transmit 2–10 mW of light spanning 1.561 to 1.656 µm (6410 to 6040 cm−1), which includes absorption lines from CO2, CH4, H2O and HDO. The returning light from each path is detected and digitized to yield the transmitted optical spectrum at a point spacing of 0.0067 cm−1 (1.5 pm) and with effectively perfect (10 ppb) frequency accuracy and narrow instrument lineshape (4×10-6 cm−1). A typical spectrum from the reference path is shown in Fig. 2. A fit of this transmitted spectrum yields the path-averaged gas concentrations. The absolute frequency accuracy and high frequency resolution of the dual-comb spectrometers translates to high precision and accuracy in the retrieved concentrations. Further, DCS spectra are undistorted by turbulence due to the simultaneous acquisition of all spectral channels and the fast sample rate of the instrument (1.6 ms spectrum−1, averaged up to 5 min here) (Rieker et al., 2014).

Figure 2Typical 32 s spectrum measured over the 2 km reference path. CO2 bands are observed in the 6350 and 6225 cm−1 regions, while CH4 and H2O are measured between 6150 and 6050 cm−1. The larger, slowly varying structure is from the comb intensity profile. The atmospheric absorption appears as the small and narrow dips.


In previous work (Waxman et al., 2017), we confirmed the high precision and accuracy possible with open-path DCS. Two DCS instruments, constructed by different teams, measured atmospheric air over adjacent paths over a 2-week period. The retrieved path-averaged gas concentrations agreed to better than 0.6 ppm (0.14 %) for CO2 and 7 ppb (0.35 %) for CH4 across the full 2-week period, where the analysis of the two DCS instruments used a common spectral database (HITRAN 2008, Rothman et al., 2009) to retrieve the concentrations from the absorption spectrum. In the work here, a single DCS instrument probes the concentrations across two different open paths simultaneously, which should further suppress any systematic offsets to below 0.45 ppm (Waxman et al., 2017). In addition, Waxman et al. (2017) compared the two DCS instruments to a stationary cavity ring-down (CRDS) point sensor whose inlet was approximately at the midpoint of the open path. This comparison actually took place over the reference path during the first 2 weeks of the present work. During that time, we found a roughly constant difference of 3.4 ppm CO2 and 17 ppb CH4 between the DCS and CRDS systems. At present, we attribute this offset to differences in the calibration scheme as the DCS is tied to the HITRAN database while the CRDS is tied to the manometric (or gravimetric depending on the gas) World Meteorological Organization (WMO) scale. Similar level offsets have been observed in comparison with the TCCON open-path FTS instrument and point-sensor-based vertical columns resulting in the TCCON CO2 scaling factor of 0.9898 (4.08 ppm for a mixing ratio of 400 ppm) (Wunch et al., 2017). This offset does not affect the results here as it is common to both the reference and over-city paths.

The reference and over-city paths had different path lengths and, therefore, used slightly different telescopes and launch powers. For the reference path, 2 mW of dual-comb light is launched from a 5 cm (2 in.) homebuilt off-axis telescope (Cossel et al., 2017; Waxman et al., 2017). The light travels to a 6.35 cm (2.5 in.) retroreflector located on a hilltop 1 km to the southwest of NIST and then is reflected back to a detector that is co-located with the launch telescope for a 1950.17±0.15 m round-trip path. Return powers vary constantly with air turbulence but we collect about 200 µW for a typical 10 dB link loss. For the city path, 10 mW of dual-comb light is launched from a modified 25.4 cm (10 in.) diameter astronomical telescope to a 12.7 cm (5 in.) retroreflector located on a building roof 3.35 km to the northeast for a 6730.66±0.15 m round-trip path. We collect about 100 µW for a typical 20 dB link loss. Round-trip path distances were measured with a laser range finder. Telescope tracking of the retroreflector is implemented to compensate for thermal drifts via a co-aligned 850 nm light-emitting diode (LED) and silicon charge-coupled device camera (Cossel et al., 2017; Waxman et al., 2017).

The measured spectra are analyzed as described in Rieker et al. (2014) and Waxman et al. (2017) at 32 s intervals. Briefly, we fit a seventh-order polynomial and HITRAN data to the measured spectrum in 100 GHz (0.333 cm−1) sections to remove the underlying structure from the comb themselves (as opposed to the atmospheric absorption). We fit the resulting absorption spectrum twice: once in the region from 6171 to 6271 cm−1 (1.595 to 1.620 µm) to obtain the path-averaged temperature from the 1.6 µm CO2 band, and once over the entire spectrum to obtain 12CO2, 13CO2, CH4, H2O, and HDO concentrations using the retrieved temperature. We then use the retrieved H2O concentration to correct the wet CO2 and CH4 mole fractions to dry mole fractions, hereafter referred to as XCO2 and XCH4, given in units of ppm and ppb (micromole of CO2 per mole of dry air, and nanomole of CH4 per mole of dry air). The correction equations are XCO2=CO2/(1-H2O) and XCH4=CH4/(1-H2O).

The variations in the retrieved concentrations are due to statistical uncertainty, systematic uncertainty (discussed above), and the true variations in the gas concentrations. Figure 8 of Waxman et al. (2017) quantified the statistical uncertainty in terms of the Allan deviation over the 2 km reference path for both XCH4 and XCO2. Figure 3 provides an Allan deviation for just XCO2 over both the ∼6.7 km city and ∼2 km reference paths, as calculated from a relatively “flat” 1000 s period of this measurement campaign on the night of 3 October 2016. As expected, the statistical uncertainty over both paths improves in relation to the square root of integration time until reaching a floor, which we attribute to real variations in the atmospheric gas concentrations. At 30 s, the statistical uncertainty of XCO2is 0.76 ppm for the reference path and 0.64 ppm for the over-city path, finally dropping to 0.21 and 0.15 ppm, respectively, at about 15 min. In most subsequent figures, we show results at a 5 min averaging time for which the statistical uncertainty is well under 0.3 ppm of XCO2 for both paths and, therefore, well below the typical atmospheric variations. Note that the uncertainty also improves with path length, as expected due to the stronger absorption. The lower uncertainty over the city path reflects the expected improvement from the 3.4× longer path length lessened by the reduction in return signal power for the longer path length.

Figure 3Statistical uncertainty as quantified by the Allan deviations for XCO2 over both the reference path (red triangles) and city path (black squares) from a well-mixed 3 h time period on the night of 3 October 2016.


2.2 Meteorological measurements

Meteorological data including pressure, wind direction, and wind speed are obtained from meteorological stations located at NCAR-Mesa and NCAR-Foothills (, last access: 28 March 2019), which are approximately the endpoints of our measurement paths (see Fig. 1), as well as from a 3-D sonic anemometer located at NIST. The path-averaged air temperature was retrieved from the CO2 spectra as described above. Finally, we obtain solar insolation from the ATOC weather station located in central Boulder (, last access: 29 March 2019).

2.3 Traffic data

We measure a subset of Boulder traffic, so we use the city traffic data to determine the fraction covered by our footprint (see Fig. 1). Traffic data from the City of Boulder are freely available at: (last access: 28 March 2019). The city provides two types of traffic data that are useful in this work: the arterial count program (ART) and the turning movement count (TMC) data.

ART measures traffic at 18 major intersections in Boulder for 5 days (1 work week, Monday through Friday) every year in 1 h bins to create a diurnal cycle. The traffic counts for 2016 are shown in Fig. 4. We use these data to scale our selected measurement time periods to a full day as discussed in Sect. 3.3.4. Note that there is only a 10 %–20 % “peak” in traffic counts at the standard commuter times with generally high traffic levels from 07:00 to  19:00, which agrees with the traffic emissions reported by the Hestia inventory model for the similar city of Salt Lake City, UT (Mitchell et al., 2018).

Figure 4City-wide traffic counts from the Boulder arterial count program (ART), normalized to a peak of unity.


TMC measures the number of vehicles at 140 intersections in Boulder for 1 work day per year during the hours of 07:45–08:45, 12:00–13:00, and 16:45–17:45. One-third of each of these sites is measured every year. We have scaled the 2014 and 2015 data to 2016 traffic levels by using total vehicle mile values available from the City of Boulder. We approximate city vehicle emissions by using the TMC locations as our source locations with a source strength scaled based on the location's fractional traffic count.

3 Results and discussion

3.1 DCS measurements

All 7.5 weeks of DCS measurements of CO2, CH4, H2O, and HDO are shown in Fig. 5. HDO is not used here but is shown for completeness (note that the HDO concentration is scaled by the isotopic abundance in HITRAN). We have insufficient precision to measure time-resolved 13CO2 concentrations over the 2 km path. However, there are very clear enhancements in the over-city path relative to the reference path for the other trace gases, especially for CO2. These enhancements are observed primarily at night when the boundary layer is lower. For example, on 13 October the CO2 enhancement reaches 129 ppm and the CH4 enhancement reaches 265 ppb. Daytime enhancements occur when the wind speed is very low and intermittent (typically below 5 m s−1), which allows emitted gases to build up over the city. When the wind increases to steady, moderate speeds, the concentrations drop quickly as the emissions are flushed out of the city. The H2O retrieval is important as accurate knowledge of the time-dependent water concentration is needed to calculate the dry CO2 and CH4 mole fractions (see Sect. 2.1). Also, the correlation of the water concentration between the two paths indicates the two paths sense the same air mass, which is further substantiated in Fig. 7 and is central to attributing their different CO2 concentration to local urban sources.

Figure 5A total of 7.5 weeks of dual-comb spectroscopy data for the reference path (red) and the over-city path (black) smoothed to 5 min time intervals. Enhancements in the over-city path relative to the reference path are observed in CO2 and CH4 but not in H2O or HDO. Note that the HDO concentration includes the HITRAN isotopic scaling.


3.2 Diurnal cycles

The diurnal cycle of XCO2 and XCH4 for both the reference path and the over-city path are shown in Fig. 6 for weekdays (midnight to midnight Monday through Friday) and weekends (midnight to midnight Saturday and Sunday). We choose to include Monday as a weekday and Saturday as a weekend because the influence of emissions from the previous day is expected to be low. The diurnal cycles of the wind direction and the wind speed measured at NCAR-Foothills are also shown in the top panel of Fig. 6. All diurnal cycles are the median values over the full 7.5 weeks of measurements and the bars reflect the 25 %/75 % quartile values.

Figure 6Diurnal cycle analysis. Data are the median of the full 7.5 weeks. (a) The mean direction in which the wind is blowing (black trace, left axis) and the wind speed (gray trace, right axis) are both from the NCAR-Foothills measurement station; shaded regions reflect the 25th to 75th quartiles; (b) the weekend and (c) weekday median XCO2 values for the over-city path (blue triangles) and reference path (red squares). Uncertainty bars represent the 25%–75% range of values encountered. Panels (d) and (e) represent the same data for XCH4. The vertical dashed black line marks 09:00 local time and the yellow shaded region highlights the region from sunrise to sunset on 22 October 2016.


The diurnal cycle of the reference path CO2 is nearly flat and nearly identical for both weekends and weekdays. It has a slight maximum between 09:00 and 10:00, with average values of 410 to 420 ppm. The diurnal cycle of the city path CO2 shows a different trend with a stronger diurnal variation. Overnight from about 18:00 to 09:00, there is an enhancement in the CO2 relative to the reference path as the CO2 from the city sources builds up due to the low winds out of the west and a presumed collapsing nighttime boundary layer. During the weekdays, this enhancement increases in the morning consistent with the rise in traffic. After the morning, the combination of the presumed rising boundary layer, increased wind speed, and shift in average wind direction out of the west (270) to the southeast (135) result in a drop in the city path CO2. Moreover, this shift in wind direction means that the reference path no longer samples the clean air from the direction of the mountains but rather sees a very similar CO2 enhancement to the city path. Fortunately, as discussed below, there are days when the wind does not shift direction, so there is a measured enhancement of the city path compared to the reference path. In the early evening, as the wind speed drops and the wind direction shifts back to out of the west, the enhancement of the city path over the reference path reappears and continues overnight as the boundary layer presumably drops. In general, the CO2 mixing ratios tend to be higher on the weekdays, sometimes exceeding 500 ppm, while weekend mixing ratios are entirely below 490 ppm. This difference is reflected in the median values as well, which reach about 440 ppm during the weekdays but only 430 ppm during the weekend.

The diurnal cycle of the reference path CH4 is relatively flat for both weekends and weekdays at just over 1.9 ppm, with a slight peak between 09:00 and 10:00. The diurnal cycle of the city path CH4 shows an enhancement, relative to the reference path, between midnight and about 09:00. We attribute this enhancement to sources of CH4 within the city, combined again with low nighttime winds and collapsing boundary layer. These sources may be leaking natural gas infrastructure such as observed in Boston (Phillips et al., 2013; McKain et al., 2015; Hendrick et al., 2016), Washington, D.C. (Jackson et al., 2014), and Indianapolis (Lamb et al., 2016). Unlike for CO2, the CH4 diurnal cycle appears unrelated to traffic (nor would we expect it to be for clean-burning vehicles) as it does not increase during high traffic times.

3.3 Estimate for CO2 emissions due to traffic

3.3.1 Measurement day selections

To select test case days to estimate the city emissions, we filter the XCO2 time series for time periods with daytime enhancement and a moderate wind strength predominantly out of the west (270). Given that the prevailing daytime winds are from the southeast (135) and often strong, this limits the test case days significantly. However, as is clear from Fig. 1, for these wind conditions, the city path samples a significant fraction of the traffic emissions and the reference path samples no traffic emissions. We only consider daytime enhancements because the nighttime boundary layer behavior is significantly more complicated than a well-mixed daytime stable boundary layer. We find 2 days that meet these criteria: Saturday 22 October 2016 from 11:00 to 16:00 and Tuesday 25 October 2016 from 07:00 to 16:00. Both days have moderate wind speeds (on average, 5 m s−1) as measured at both meteorological sites. There are additional days with daytime enhancement in XCO2, but the wind direction is variable. Additionally, there are many days with no daytime enhancement in XCO2 because the high wind speeds (6 m s−1 or higher) prevented buildup of CO2. We use 22 October as a proxy for all weekend days and 25 October as a proxy for all weekdays. The XCO2 and XCH4 mixing ratios, as well as wind speed and wind direction, are shown in Fig. 7 for these 2 case study days.

Figure 7Footprint calculations and time series data for the 2 case study days. Left column – Saturday, 22 October 2016 data; right column – Tuesday, 25 October 2016 data. Upper panels (a, d): footprints for the reference path. Middle panels (b, e): footprints for the over-city path. The footprints are averaged over the respective time windows and open paths. Lower panels (c, f): wind and CO2 data at 5 min time intervals. Reference and over-city measurement paths are shown in red and black, respectively. Data plots show XCO2 over the reference path (red) and city path (black), wind speed and wind direction measurements taken at NCAR-Mesa (blue) and NCAR-Foothills (orange), and the calculated Q(t). On 25 October, Q(t) data near 14:00 have been removed since the reference path wind direction is out of the southeast to the east, resulting in city contamination along the reference path. All data are smoothed to 5 min time intervals.


In order to confirm that the reference path measured clean background air and the over-city path measured city emissions, we calculated footprints for the two test case time periods using the Stochastic Time-Inverted Lagrangian Transport (STILT-R) model (Fasoli et al., 2018). The input meteorology file consisted of a uniform wind field with wind data from the NCAR-Foothills lab, boundary layer height from the North American Regional Reanalysis (NARR), and uniform turbulent velocity variance calculated from the Pasquill stability class (determined from wind speed and solar insolation) from the ground up to the boundary layer. We also used hyper near-field scaling described in Fasoli et al. (2018). Average footprints for the two time periods are shown in Fig. 7. The footprint for the reference path covers undeveloped areas extending from the near foothills into the mountains. The footprint for the over-city path also has contributions from the same general mountain region. In addition, this path has sensitivity to an extended area within the city and, therefore, to a large fraction of the traffic emissions. Note that the open-path geometry leads to a much larger extended footprint for this path than would be the case for a single point sensor located at the same height within the city.

The variability in the reference CO2 on both days is a real atmospheric effect. In processing, any data are removed if the signal power is low, which is indicative of poor telescope alignment or strong weather-related attenuation over the beam path, so the variability is not due to variable signal strength. We attribute this variability to the smaller footprint of the reference path relative to the over-city path, as seen in Fig. 7. If the CO2 in the air is not fully mixed, then the temporal and spatial variability will be more evident in the path with the smaller footprint.

To convert from the measured enhancement to an emissions rate, we require a model that connects the source strength to the plume concentration. Since we do not have a high-resolution, spatially resolved inventory for Boulder similar to the Hestia model for Salt Lake City (Mitchell et al., 2018), we use the existing Boulder traffic inventory (see Sect. 2.3) in conjunction with a Gaussian plume model.

3.3.2 Gaussian plume calculations

The standard Gaussian plume model that includes total reflection at the Earth's surface is as follows (Seinfeld and Pandis, 2006):


where (x,y,z) is the location in space for which the plume concentration is being calculated; (x0,y0,H) is the emissions location; c(x,y,z,t) is the concentration at location (x,y,z) and time t; q is the emissions strength (usually in kg s−1); σy and σz are the plume variances in the y and z direction as a function of travel distance and Pasquill stability class (Seinfeld and Pandis, 2006); and u is the wind speed in m s−1. The wind is assumed to be in the x direction. The plume variances are calculated as follows:

(2) σ y = exp I y + J y ( ln Δ x ) + K y ( ln Δ x ) 2 ,


(3) σ z = exp I z + J z ( ln Δ x ) + K z ( ln Δ x ) 2 ,

where Iy, Jy, Ky, Iz, Jz, and Kz are from a look-up table based on the Pasquill stability class, which depends on the wind speed and solar insolation (Seinfeld and Pandis, 2006) and Δx is the x distance relative to the plume origin. This plume model does not include any reflection at the boundary layer height; however, due to the small spatial scales, this effect is negligible here.

We modify this equation in several ways: (1) since we measure the column-integrated concentration over a finite beam path at an angle to the wind direction, we integrate the plume concentration along this beam path and then normalize to the length of the beam path; (2) we sum over the emissions locations in the city that contribute emissions to our measurements. Thus, our overall measurement equation is as follows:


where (c-c0) is our path-integrated concentration enhancement measurement (in t m−3 and t is metric tons; 1 t =1000 kg) along our path s, which goes from 0 to L; Q is the total city emissions in t h−1; L is our path length in m; (xj,yj) are the source emissions locations; fj is the fraction of traffic at source location (xj,yj) relative to traffic over all locations in the city from the TMC database; u is the wind speed in m s−1; θ is the angle of the beam path with respect to the wind direction; and σy and σz are the plume dispersions in m in the y and z directions, which depend on the sources distance from the beam path. In writing Eq. (3), we assume the wind is in the +x^ direction (this assumption is relaxed below). We assume that all plume emission locations are vehicle tailpipes at 1 m above the ground, and the beam path runs 15 m above ground so all measurement heights are at 15 m above ground.

Grid rotation for variable wind directions

To calculate Eq. (3), we grid the emissions locations using UTM (Universal Transverse Mercator) coordinates obtained from Google Earth, where we then define north as +y^ and east as +x^. We translate the coordinate system such that the DCS path begins at the origin (0,0) and travels a distance L at angle θ with respect to the x axis. Eq. (3) is then valid provided the wind is directly in the +x^ direction. More generally, the wind is at a time-varying small angle ϕ(t) with respect to +x^. Therefore, we apply a rotation about the origin (Prussin et al., 2015):


to generate new traffic coordinates (xj,yj) and a new parameterized DCS beam path of (s cos(θ), s sin(θ)), where θ=θ-ϕ(t). In this new coordinate system, the wind is along the +x^ direction and Eq. (3) holds with the substitutions θθ and yjyj and where the σy and σz are calculated based on the distance Δx=|xj-(yj/tanθ)|.

Time-dependent estimate of Q(t)

The rotated Eq. (3) can be solved for Q in terms of the measured or estimated values of c(t)−c0(t), u(t), Δϕ(t), σy(t), σz(t), θ, L, and fi, where the first five quantities are time dependent. The resulting, time-dependent Q(t) for each test case day is shown in the bottom panels of Fig. 7 and has a mean value and standard deviation of Q22Oct=31±17 t CO2 h−1 for 22 October and Q25Oct=165±45 t CO2 h−1 for 25 October for the 5 min averaged data as shown.

Uncertainty in Q(t)

Seven measured parameters factor into the emissions calculation of Q(t) for the 2 days. These are given in Table 1 along with the instrumental measurement precision and the observed variability. Note that solar insolation is used solely in the determination of the Pasquill stability class (Seinfeld and Pandis, 2006). The stability class is relatively insensitive to the variations in solar insolation observed on the 2 test case days. As can be seen in the table, the uncertainty is dominated by the natural variability in parameters like wind speed, wind direction, and CO2 concentration rather than the DCS spectrometer precision. The observed variability over the 5–9 h period is typically at least a factor of 2 larger than the instrument precision. The variability in these parameters leads to the observed variability in Q(t). We use the mean of Q(t) as our emissions value and the standard deviation (at 5 min time-averaging) as its uncertainty. In using this standard deviation as a measure of the uncertainty, we attempt to capture the uncertainty associated with the discrepancies between, for example, the weather station measurements of wind direction and speed relative to the true wind direction (which results in greater or fewer number of plumes from the given traffic locations intercepting the measurement path). This variability appears in Q(t) as the nominal measured wind direction varies. Future systems with redundant, distributed DCS beam paths would provide a superior estimate of all of these uncertainties.

Table 1Parameters used to calculate the emission rate from Eq. (3). The measurement precision refers to the instrument uncertainty in the measurement quantity. The variability refers to the observed environmental variability over the measurement period. The variability from the enhancement, the wind direction, and the wind speed drive the observed variability in the estimated Q(t). The distance from a given source location to the DCS measurement path, Δxj , varies with location and has a 5 m uncertainty.

Download Print Version | Download XLSX

In addition, there are assumptions and possible uncertainties inherent to the Gaussian plume model. First, the model does not include the effects of buildings, trees, or other objects that could break up the plume between the emissions location and the beam path. Second, we assume that all CO2 emissions come from the discrete locations shown in Fig. 1, while in reality the emissions are likely substantially more diffuse. The assumption of discrete emissions simplifies modeling and is feasible due to the city traffic data but may result in a bias due to the coarse distribution of traffic measurements. Third, we approximate the measurement height at 15 m above ground although the beam height differs over the path since Boulder is not perfectly flat. Finally, we use standard Iy, Jy, Ky, Iz, Jz, and Kz values that were derived for rural areas (Turner, 1970), which may be different to urban or suburban areas. However, the greatest differences between rural and urban conditions are expected to be at night (Turner, 1970).

Further, we ran plume calculations in STILT-R using both wind fields derived from the local meteorological stations shown in Fig. 1 and using the North American Mesoscale Forecast System (NAM,, last access: 28 March 2019). The High-Resolution Rapid Refresh (HRRR,, last access: 28 March 2019) and North American Regional Reanalysis (NARR,, last access: 28 March 2019) wind projections did not match the measured winds at the meteorological stations. These calculations produced emissions values ranging between 55 and 770 t h−1, depending on the wind fields and vertical dispersion parameterization used. This brackets our emission calculations by approximately a factor of 3 in each direction and shows how sensitive these kilometer-scale measurements are to vertical dispersion.

3.3.3 Corrections for non-traffic sources of CO2

There are a number of non-traffic sources of CO2 that could contribute to our measured XCO2 enhancement, including local power plants, residential emissions, and biological activity. These non-traffic sources should have a relatively minor contribution for several reasons. First, the footprint of the over-city path does not overlap the large power plant to the east of the Boulder city limits. Second, the temperature during the 2 test case days was 24 and 20 C (68 and 75 F) on 22 and 25 October, leading to minimal residential and commercial heating. Third, the measurements occurred in October after leaf senescence so there should be negligible biological activity. Nevertheless, as discussed below, we do adjust our measurements to account for the relatively minor contribution from non-traffic sources before scaling up to an estimate of the annual traffic emissions.

We first consider power plants. There are two power generation facilities on the Department of Commerce (DOC) campus located near the NIST building that houses the dual-comb spectrometer: the site's Central Utilities Plant (CUP), and the National Oceanic and Atmospheric Administration (NOAA) building's boilers. To calculate their average CO2 emissions, we used available fuel consumption data (October 2016 monthly average for the CUP and mid-November to mid-December 2016 average for the NOAA boilers; October data were unavailable) and the EPA emissions factor (EPA, 1995). We then modeled the CUP and boiler plume emissions using WindTrax (Flesch et al., 1995, 2004) with wind speed and direction data from the NCAR-Mesa site. We find that due to the moderate wind speeds (∼5 m s−1) during our case study days and the height mismatch between the emission stacks and our measurement path over the DOC campus, there is negligible enhancement over the reference path. Given the location of the emission sources and the wind direction during our measurement periods, the emissions also do not cross the over-city beam path. Therefore, we apply no correction for these two power plant emissions.

The University of Colorado also has a power plant that falls within the main footprint associated with the over-city beam path, shown in Fig. 7, and whose emissions are expected to intersect our over-city beam path. The EPA Greenhouse Gas Reporting Program (GHGRP,, last access: 28 March 2019) lists the 2017 emission from the power plant as 2.7×104 t CO2 or an average of 3.1 t h−1 (no breakdown by season or hour is provided). We apply this correction to our previous daily values and add a conservative uncertainty equal to this correction in quadrature with the previous uncertainty. The new adjusted values are then 28±17 t CO2 h−1 for 22 October and 162±45 t CO2 h−1 for 25 October.

The large Valmont power station lies just outside the city limits to the east of Boulder; however, given its location and the dominant selected westerly wind, emissions from this source do not reach our beam paths. There are no other power generation facilities within the city that report to the GHGRP, so we make no further corrections based on power plants.

In addition, there are also likely diffuse emissions from residential and commercial furnaces and water heaters that use natural gas. The City of Boulder Community Greenhouse Gas Emissions Inventory reports 20 % of the city emissions, or 3.18×105 t CO2, were from natural gas in 2016 (, last access: 28 March 2019). The natural gas usage varies strongly by month with building heating requirements. Although our measurements occurred in October, the measurement days were quite warm (20–24 C) so that residential and commercial building heating was unlikely and the use of an annual average would overestimate any contribution. Instead, we scale the natural gas usage according to the monthly breakdown provided by the United States Energy Information Administration database for Colorado (, last access: 28 March 2019). The mean daytime (approximately sunrise to sunset, 07:00 to 18:00) temperature in October was 18.2 C while the mean temperature (including day and night) for October was 15.7 C. Our daytime-only measurements, therefore, had a mean temperature that was much closer to the mean temperature (day and night) of September, which was 19.2 C. Therefore, we scale the Boulder annual natural gas consumption by the September 2016 natural gas usage, which was 2.4 % of the Colorado annual total (, last access: 28 March 2019). The estimated total emissions from residential and commercial natural gas usage in Boulder over our measurement days is then 10.2 t CO2 h−1. We apply this correction to our measured values and include a (conservative) uncertainty equal to this correction. The new adjusted values are then Q22Oct,adj=18±20 t CO2 h−1 for 22 October and Q25Oct,adj=152±46 t CO2 h−1 for 25 October.

Once leaf senescence has completed, neither plants nor soil respiration contribute to CO2 signal (Matyssek et al., 2013). The National Phenology Network (USA National Phenology Network, 2018) data show that for the site nearest to Boulder (64 km north of Boulder), the leaf fall dates were 15 September 2016 for box elder trees and 6 October 2016 for the eastern cottonwood. Thus, by our measurement dates leaf senescence should be fully complete and plants will not contribute to the city CO2 enhancement. We note that a wide range of biogenic contributions to CO2 have been noted in the literature (Gurney et al., 2017; Mitchell et al., 2018; Sargent et al., 2018).

3.3.4 Scaling to annual emissions

In order to compare with the city inventory, we scale our results to an annual total. To do this, we use the hourly traffic data of Fig. 4 to scale Q22Oct,adj and Q25Oct,adj to a daily emission. Based on Fig. 4, 34% of the total traffic counts occur during the 5 h measurement period on 22 October and 52 % of the total traffic counts occur during the 8 h measurement period on 25 October (excluding the 13:00 to 14:00 period). The daily emissions are then Q22Oct,day=Q22Oct,adj×(5h)/(0.34) and Q25Oct,day=Q25Oct,adj×(8h)/(0.52). The traffic data in Fig. 4 are based on weekday measurement, and we assume that the hourly distribution is the same for weekends; this may lead to a slight overestimate in the weekend data where a larger fraction of emissions occurs between 11:00 and 16:00 than on weekdays. We then scale to annual emissions by assuming that the emissions on 22 October are representative of all 112 weekend and holiday days and the emissions on 25 October are representative of all 253 workdays. Including their uncertainty, this calculation yields (6.2±1.8)×105 t CO2 yr−1.

The scaling relies heavily on the traffic count data supplied by the city of Boulder, which do not have an associated uncertainty value. A comparison of these data over several years shows a typical 7 % statistical variation at a given TMC location after removing a linear trend. We assume this reflects day-to-day fluctuations in traffic. In addition, there will be seasonal variations, which are not captured in the extrapolation from our 2 test case days to the annual emissions. Due to the lack of seasonal data for Boulder traffic, we use the detailed Hestia traffic inventory for Salt Lake City, UT, given in Fig. 2 of Mitchell et al. (2018). These data show a variation of ±18 % in traffic emissions between “summer” and “winter” months. Combined in quadrature with the 7 % statistical uncertainty in the TMC traffic count data, this leads to an additional ∼20 % uncertainty to the scaled annual estimate. As noted earlier, we have not applied any additional uncertainty on the reliance on the TMC data as a proxy for emissions locations.

Including the additional uncertainty on the scaling to annual emissions, we estimate an annual emission rate of (6.2±2.2)×105 t CO2 yr−1 for traffic carbon emissions for Boulder, CO.

4 Comparison with city estimates

The city vehicle emissions estimate comes from total vehicle miles traveled based on data from the transportation department, miles per gallon inputs from the EPA state inventory tool, and vehicle type distribution from the Colorado Department of Public Health and the Environment (Kimberlee Rankin, City of Boulder, personal communication). The City of Boulder estimates total vehicle emissions of 4.50×105 t CO2 in 2016 (, last access: 28 March 2019). On-road emissions account for greater than 99 % of the transportation emissions, so we have scaled this value down by 1 % for an on-road emissions value of 4.46×105 t CO2. We assume that all traffic emissions are CO2 rather than a mix of CO2 and CH4. There is no uncertainty provided by the city on this value.

In comparison, we estimate (6.2±2.2)×105 t CO2 yr−1, which is 139 % of the city estimate but agrees within the given uncertainty. Interestingly, other studies have also found that emissions measurements were higher than the reported inventory values. Brioude et al. (2013) found top-down aircraft estimates of Los Angeles county and the South Coast Air Basin CO2 were 1.45 times larger than the Vulcan 2005 inventory (Gurney et al., 2009). An earlier aircraft campaign over Sacramento, CA, found an average CO2 emission, with 100 % uncertainty, that was 15 %–20 % higher than the Vulcan estimate (Turnbull et al., 2011). Lauvaux et al. (2016) compared Indianapolis city CO2 emissions measured by a network of CRDS instruments to the HESTIA inventory (Gurney et al., 2012) during INFLUX (Davis et al., 2017). They found that despite the building-scale resolution in the HESTIA inventory, it still underestimated the annual CO2 flux by 20 %. An updated version of HESTIA predicted very similar emissions estimates for on-road, residential, and commercial sectors, so the discrepancy was attributed to missing sources of CO2, including animal (primarily human and companion animal) respiration, biofuel combustion, and biosphere respiration (Gurney et al., 2017).

4.1 Improvements in future measurements

Future improvements should include additional and different beam paths, selected based on prevailing wind directions. Our initial assumption that the mountain path would generally act as a reference path was incorrect since the prevailing daytime winds during this time of year are not out of the west but rather the southeast. An east–west running beam north of the city and one south of the city would allow us to utilize a larger fraction of the data as the predominant midday wind direction during the fall is out of the north to the northeast (see Fig. 1). Even longer beam paths would also interrogate a larger fraction of the city and measure a correspondingly larger fraction of the vehicle emissions. Vertically resolved data from, e.g., a series of stacked retroreflectors would better test the assumption of vertically dispersing Gaussian plumes.

Additionally, more extensive modeling to cover variable wind directions and speeds would allow the incorporation of a much larger fraction of the data than the 2 days selected here. An inversion-based model similar to Lauvaux et al. (2013) could potentially be applied to a small city like Boulder; however, this would depend heavily on the quality of the bottom-up emissions inventory used to generate the priors. Indeed, one of the major future improvements would be to generate a detailed Hestia inventory of Boulder, CO, similar to that generated for Salt Lake City, UT (Mitchell et al., 2018).

5 Conclusions

We demonstrate the use of an open-path dual frequency comb spectroscopy system for quantifying city emissions of carbon dioxide. We send light over two paths: a reference path that samples the concentration of gases entering the city from the west, and an over-city path that measures the concentrations of gases after the air mass has crossed approximately two-thirds of the city, including two major commuter arteries. The measured diurnal cycle shows a significant traffic-related enhancement in the carbon dioxide signal during weekdays in the over-city path compared to the reference path. We select 2 case study days with appropriate wind conditions and apply Gaussian plume modeling to estimate the total vehicular carbon emission. We then scale these results up to annual city-wide emissions using traffic data from the City of Boulder. We find overall traffic-related carbon emissions that are approximately 1.4 times greater than the city's bottom-up traffic emissions inventory but with an uncertainty that encompasses the city inventory estimate. Further improvements to this method should include improved design of reference and over-city paths and a more detailed inventory model for Boulder CO, which together should further reduce the overall uncertainty in the estimate.

Data availability

As per NIST regulations, all data are archived at NIST and available upon request.

Appendix A: Modification of the Gaussian plume equation

Equation (1) is the standard Gaussian plume equation as discussed in Sect. 3.3.2 (Seinfeld and Pandis, 2006). It is reproduced here:


where the standard variables are as defined in Sect. 3.3.2.

A1 Path-integrated substitutions

The DCS returns the average concentration along a line path. We denote distance along this path by the variable s, where s runs from 0 to L. This path is assumed to lie in the xy plane at an angle θ with respect to the x axis (which is assumed to be the wind direction in the standard Gaussian plume equation). With these definitions, the contribution to the DCS signal from the plume is as follows:




A2 Accounting for multiple point sources

Rather than a single source at (x0, y0), we have multiple sources at locations (xj, yj), each with a source strength fjq, where fj is the fractional source strength out of the total value q. We now sum over all sources to find the total enhancement. We also change the units of q from kg s−1 to t yr−1 and thus change the emissions variable to Q to indicate the unit change. This gives the following equation:


A3 Height substitutions

We assume that the point source emissions locations are 1 m above ground (z=1), and city topographic data indicate that our beam path is approximately 15 m above ground (H=15). These substitutions finally lead to Eq. (3) in the main text.

Author contributions

EMW, KCC, IC, and NRN designed the experiment. WCS helped build the hardware for the open-path measurements. EMW, KCC, and GWT ran the experiment. FG wrote the processing code for the data analysis. EMW processed the data and did the Gaussian plume modeling. KCC did the STILT-R modeling. EMW, KCC, IC, and NRN cowrote the manuscript.

Competing interests

The authors declare that they have no conflict of interest.


We thank Kimberlee Rankin, Randall Rutsch, Bill Cowern, and Chris Hagelin from the City of Boulder for city inventory and traffic information; Anna Karion for assistance with STILT-R modeling; and Dave Plusquellic and Caroline Alden for assistance with the manuscript. This work was funded by Defense Advanced Research Program Agency DSO SCOUT program and James Whetstone and the NIST Special Programs Office. Eleanor M. Waxman and Kevin C. Cossel are partially supported by National Research Council postdoctoral fellowships.

Review statement

This paper was edited by Ronald Cohen and reviewed by three anonymous referees.


Bergeron, O. and Strachan, I. B.: CO2 sources and sinks in urban and suburban areas of a northern mid-latitude city, Atmos. Environ., 45, 1564–1573,, 2011. 

Bréon, F. M., Broquet, G., Puygrenier, V., Chevallier, F., Xueref-Remy, I., Ramonet, M., Dieudonné, E., Lopez, M., Schmidt, M., Perrussel, O., and Ciais, P.: An attempt at estimating Paris area CO2 emissions from atmospheric concentration measurements, Atmos. Chem. Phys., 15, 1707–1724,, 2015. 

Brioude, J., Angevine, W. M., Ahmadov, R., Kim, S.-W., Evan, S., McKeen, S. A., Hsie, E.-Y., Frost, G. J., Neuman, J. A., Pollack, I. B., Peischl, J., Ryerson, T. B., Holloway, J., Brown, S. S., Nowak, J. B., Roberts, J. M., Wofsy, S. C., Santoni, G. W., Oda, T., and Trainer, M.: Top-down estimate of surface flux in the Los Angeles Basin using a mesoscale inverse modeling technique: assessing anthropogenic emissions of CO, NOx and CO2 and their impacts, Atmos. Chem. Phys., 13, 3661–3677,, 2013. 

Ciais, P., Rayner, P., Chevallier, F., Bousquet, P., Logan, M., Peylin, P., and Ramonet, M.: Atmospheric inversions for estimating CO2 fluxes: methods and perspectives, Clim. Change, 103, 69–92,, 2010. 

Coburn, S., Alden, C. B., Wright, R., Cossel, K., Baumann, E., Truong, G.-W., Giorgetta, F., Sweeney, C., Newbury, N. R., Prasad, K., Coddington, I., and Rieker, G. B.: Regional trace-gas source attribution using a field-deployed dual frequency comb spectrometer, Optica, 5, 320–327,, 2018. 

Coddington, I., Newbury, N., and Swann, W.: Dual-comb spectroscopy, Optica, 3, 414–426,, 2016. 

Cossel, K. C., Waxman, E. M., Giorgetta, F. R., Cermak, M., Coddington, I. R., Hesselius, D., Ruben, S., Swann, W. C., Truong, G.-W., Rieker, G. B., and Newbury, N. R.: Open-path dual-comb spectroscopy to an airborne retroreflector, Optica, 4, 724–728,, 2017. 

Coutts, A. M., Beringer, J., and Tapper, N. J.: Characteristics influencing the variability of urban CO2 fluxes in Melbourne, Australia, Atmos. Environ., 41, 51–62,, 2007. 

Davis, K. J., Deng, A., Lauvaux, T., Miles, N. L., Richardson, S. J., Sarmiento, D. P., Gurney, K. R., Hardesty, R. M., Bonin, T. A., Brewer, W. A., Lamb, B. K., Shepson, P. B., Harvey, R. M., Cambaliza, M. O., Sweeney, C., Turnbull, J. C., Whetstone, J., and Karion, A.: The Indianapolis Flux Experiment (INFLUX): A test-bed for developing urban greenhouse gas emission measurements, Elem. Sci. Anth., 5, 21,, 2017. 

Dobler, J. T., Zaccheo, T. S., Pernini, T. G., Blume, N., Broquet, G., Vogel, F., Ramonet, M., Braun, M., Staufer, J., Ciais, P., and Botos, C.: Demonstration of spatial greenhouse gas mapping using laser absorption spectrometers on local scales, J. Appl. Remote Sens., 11, 014002,, 2017. 

EPA: AP 42, Fifth Edition Compilation of Air Pollutant Emissions Factors, Volume 1: Stationary Point and Area Sources, available at: (last access: 28 March 2019), 1995. 

Fasoli, B., Lin, J. C., Bowling, D. R., Mitchell, L., and Mendoza, D.: Simulating atmospheric tracer concentrations for spatially distributed receptors: updates to the Stochastic Time-Inverted Lagrangian Transport model's R interface (STILT-R version 2), Geosci. Model Dev., 11, 2813–2824,, 2018. 

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. 

Gurney, K. R., Mendoza, D. L., Zhou, Y., Fischer, M. L., Miller, C. C., Geethakumar, S., and Du Can, S. D. L. R.: High Resolution Fossil Fuel Combustion CO2 Emission Fluxes for the United States, Environ. Sci. Technol., 43, 5535–5541,, 2009. 

Gurney, K. R., Razlivanov, I., Song, Y., Zhou, Y., Benes, B., and Abdul-Massih, M.: Quantification of Fossil Fuel CO2 Emissions on the Building/Street Scale for a Large U.S. City, Environ. Sci. Technol., 46, 12194–12202,, 2012. 

Gurney, K. R., Liang, J., Patarasuk, R., O'Keeffe, D., Huang, J., Hutchins, M., Lauvaux, T., Turnbull, J. C., and Shepson, P. B.: Reconciling the differences between a bottom-up and inverse-estimated FFCO2 emissions estimate in a large US urban area, Elem. Sci. Anth., 5, 44,, 2017. 

Hedelius, J. K., Liu, J., Oda, T., Maksyutov, S., Roehl, C. M., Iraci, L. T., Podolske, J. R., Hillyard, P. W., Liang, J., Gurney, K. R., Wunch, D., and Wennberg, P. O.: Southern California megacity CO2, CH4, and CO flux estimates using ground- and space-based remote sensing and a Lagrangian model, Atmos. Chem. Phys., 18, 16271–16291,, 2018. 

Heimburger, A. M. F., Harvey, R. M., Shepson, P. B., Stirm, B. H., Gore, C., Turnbull, J., Cambaliza, M. O. L., Salmon, O. E., Kerlo, A.-E. M., Lavoie, T. N., Davis, K. J., Lauvaux, T., Karion, A., Sweeney, C., Brewer, W. A., Hardesty, R. M., and Gurney, K. R.: Assessing the optimized precision of the aircraft mass balance method for measurement of urban greenhouse gas emission rates through averaging, Elem. Sci. Anth., 5, 26,, 2017. 

Hendrick, M. F., Ackley, R., Sanaie-Movahed, B., Tang, X., and Phillips, N. G.: Fugitive methane emissions from leak-prone natural gas distribution infrastructure in urban environments, Environ. Pollut., 213, 710–716,, 2016. 

Jackson, R. B., Down, A., Phillips, N. G., Ackley, R. C., Cook, C. W., Plata, D. L., and Zhao, K.: Natural Gas Pipeline Leaks Across Washington, DC, Environ. Sci. Technol., 48, 2051–2058,, 2014. 

Järvi, L., Rannik, Ü., Kokkonen, T. V., Kurppa, M., Karppinen, A., Kouznetsov, R. D., Rantala, P., Vesala, T., and Wood, C. R.: Uncertainty of eddy covariance flux measurements over an urban area based on two towers, Atmos. Meas. Tech., 11, 5421–5438,, 2018. 

Lamb, B. K., Cambaliza, M. O. L., Davis, K. J., Edburg, S. L., Ferrara, T. W., Floerchinger, C., Heimburger, A. M. F., Herndon, S., Lauvaux, T., Lavoie, T., Lyon, D. R., Miles, N., Prasad, K. R., Richardson, S., Roscioli, J. R., Salmon, O. E., Shepson, P. B., Stirm, B. H., and Whetstone, J.: Direct and Indirect Measurements and Modeling of Methane Emissions in Indianapolis, Indiana, Environ. Sci. Technol., 50, 8910–8917,, 2016. 

Lauvaux, T., Miles, N. L., Richardson, S. J., Deng, A., Stauffer, D. R., Davis, K. J., Jacobson, G., Rella, C., Calonder, G.-P., and DeCola, P. L.: Urban Emissions of CO2 from Davos, Switzerland: The First Real-Time Monitoring System Using an Atmospheric Inversion Technique, J. Appl. Meteorol. Clim., 52, 2654–2668,, 2013. 

Lauvaux, T., Miles, N. L., Deng, A., Richardson, S. J., Cambaliza, M. O., Davis, K. J., Gaudet, B., Gurney, K. R., Huang, J., O'Keefe, D., Song, Y., Karion, A., Oda, T., Patarasuk, R., Razlivanov, I., Sarmiento, D., Shepson, P., Sweeney, C., Turnbull, J., and Wu, K.: High-resolution atmospheric inversion of urban CO2 emissions during the dormant season of the Indianapolis Flux Experiment (INFLUX), J. Geophys. Res.-Atmos., 121, 2015JD024473,, 2016. 

Marcotullio, P. J., Sarzynski, A., Albrecht, J., Schulz, N., and Garcia, J.: The geography of global urban greenhouse gas emissions: an exploratory analysis, Clim. Change, 121, 621–634,, 2013. 

Matyssek, R., Clarke, N., Cudlin, P., Mikkelsen, T. N., Tuovinen, J.-P., Wieser, G., and Paoletti, E.: Climate Change, Air Pollution and Global Challenges: Understanding and Perspectives from Forest Research, Elsevier, London, UK, available at: (last access: 19 December 2018), 2013. 

Mays, K. L., Shepson, P. B., Stirm, B. H., Karion, A., Sweeney, C., and Gurney, K. R.: Aircraft-Based Measurements of the Carbon Footprint of Indianapolis, Environ. Sci. Technol., 43, 7816–7823,, 2009. 

McKain, K., Wofsy, S. C., Nehrkorn, T., Eluszkiewicz, J., Ehleringer, J. R., and Stephens, B. B.: Assessment of ground-based atmospheric observations for verification of greenhouse gas emissions from an urban region, P. Natl. Acad. Sci. USA, 109, 8423–8428,, 2012. 

McKain, K., Down, A., Raciti, S. M., Budney, J., Hutyra, L. R., Floerchinger, C., Herndon, S. C., Nehrkorn, T., Zahniser, M. S., Jackson, R. B., Phillips, N., and Wofsy, S. C.: Methane emissions from natural gas infrastructure and use in the urban region of Boston, Massachusetts, P. Natl. Acad. Sci. USA, 112, 1941–1946,, 2015. 

Mitchell, L. E., Lin, J. C., Bowling, D. R., Pataki, D. E., Strong, C., Schauer, A. J., Bares, R., Bush, S. E., Stephens, B. B., Mendoza, D., Mallia, D., Holland, L., Gurney, K. R., and Ehleringer, J. R.: Long-term urban carbon dioxide observations reveal spatial and temporal dynamics related to urban characteristics and growth, P. Natl. Acad. Sci. USA, 115, 2912–2917,, 2018. 

Mueller, K., Yadav, V., Lopez-Coto, I., Karion, A., Gourdji, S., Martin, C., and Whetstone, J.: Siting background towers to characterize incoming air for urban greenhouse gas estimation: a case study in the Washington DC/Baltimore Area, J. Geophys. Res.-Atmos., 123, 2910–2926,, 2017. 

Nemitz, E., Hargreaves, K. J., McDonald, A. G., Dorsey, J. R., and Fowler, D.: Micrometeorological Measurements of the Urban Heat Budget and CO2 Emissions on a City Scale, Environ. Sci. Technol., 36, 3139–3146,, 2002. 

Phillips, N. G., Ackley, R., Crosson, E. R., Down, A., Hutyra, L. R., Brondfield, M., Karr, J. D., Zhao, K., and Jackson, R. B.: Mapping urban pipeline leaks: Methane leaks across Boston, Environ. Pollut., 173, 1–4,, 2013. 

Prussin, A. J., Marr, L. C., Schmale, D. G., Stoll, R., and Ross, S. D.: Experimental validation of a long-distance transport model for plant pathogens: Application to Fusarium graminearum, Agric. Forest Meteorol., 203, 118–130,, 2015. 

Rieker, G. B., Giorgetta, F. R., Swann, W. C., Kofler, J., Zolot, A. M., Sinclair, L. C., Baumann, E., Cromer, C., Petron, G., Sweeney, C., Tans, P. P., Coddington, I., and Newbury, N. R.: Frequency-comb-based remote sensing of greenhouse gases over kilometer air paths, Optica, 1, 290–298,, 2014. 

Rothman, L. S., Gordon, I. E., Barbe, A., Benner, D. C., Bernath, P. E., Birk, M., Boudon, V., Brown, L. R., Campargue, A., Champion, J. P., Chance, K., Coudert, L. H., Dana, V., Devi, V. M., Fally, S., Flaud, J. M., Gamache, R. R., Goldman, A., Jacquemart, D., Kleiner, I., Lacome, N., Lafferty, W. J., Mandin, J. Y., Massie, S. T., Mikhailenko, S. N., Miller, C. E., Moazzen-Ahmadi, N., Naumenko, O. V., Nikitin, A. V., Orphal, J., Perevalov, V. I., Perrin, A., Predoi-Cross, A., Rinsland, C. P., Rotger, M., Simeckova, M., Smith, M. A. H., Sung, K., Tashkun, S. A., Tennyson, J., Toth, R. A., Vandaele, A. C., and Vander Auwera, J.: The HITRAN 2008 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 110, 533–572,, 2009. 

Ryerson, T. B., Trainer, M., Holloway, J. S., Parrish, D. D., Huey, L. G., Sueper, D. T., Frost, G. J., Donnelly, S. G., Schauffler, S., Atlas, E. L., Kuster, W. C., Goldan, P. D., Hübler, G., Meagher, J. F., and Fehsenfeld, F. C.: Observations of Ozone Formation in Power Plant Plumes and Implications for Ozone Control Strategies, Science, 292, 719–723,, 2001. 

Sargent, M., Barrera, Y., Nehrkorn, T., Hutyra, L. R., Gately, C. K., Jones, T., McKain, K., Sweeney, C., Hegarty, J., Hardiman, B., Wang, J. A., and Wofsy, S. C.: Anthropogenic and biogenic CO2 fluxes in the Boston urban region, P. Natl. Acad. Sci. USA, 115, 7491–7496,, 2018. 

Seinfeld, J. H. and Pandis, S. N.: Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, Wiley, Hoboken, New Jersey, USA, 2006. 

Seto, K. C., Bigio, A., Bento, A., Cervero, R., and Christensen, P.: Human Settlements, Infrastructure, and Spatial Planning, in: Climate Change 2014: Mitigation of Coimate Change, Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, 78 pp., 2014. 

Shusterman, A. A., Teige, V. E., Turner, A. J., Newman, C., Kim, J., and Cohen, R. C.: The BErkeley Atmospheric CO2 Observation Network: initial evaluation, Atmos. Chem. Phys., 16, 13449–13463,, 2016. 

Staufer, J., Broquet, G., Bréon, F.-M., Puygrenier, V., Chevallier, F., Xueref-Rémy, I., Dieudonné, E., Lopez, M., Schmidt, M., Ramonet, M., Perrussel, O., Lac, C., Wu, L., and Ciais, P.: The first 1-year-long estimate of the Paris region fossil fuel CO2 emissions based on atmospheric inversion, Atmos. Chem. Phys., 16, 14703–14726,, 2016. 

Truong, G.-W., Waxman, E. M., Cossel, K. C., Baumann, E., Klose, A., Giorgetta, F. R., Swann, W. C., Newbury, N. R., and Coddington, I.: Accurate frequency referencing for fieldable dual-comb spectroscopy, Opt. Express, 24, 30495–30504,, 2016. 

Turnbull, J. C., Karion, A., Fischer, M. L., Faloona, I., Guilderson, T., Lehman, S. J., Miller, B. R., Miller, J. B., Montzka, S., Sherwood, T., Saripalli, S., Sweeney, C., and Tans, P. P.: Assessment of fossil fuel carbon dioxide and other anthropogenic trace gas emissions from airborne measurements over Sacramento, California in spring 2009, Atmos. Chem. Phys., 11, 705–721,, 2011. 

Turner, A. J., Shusterman, A. A., McDonald, B. C., Teige, V., Harley, R. A., and Cohen, R. C.: Network design for quantifying urban CO2 emissions: assessing trade-offs between precision and network density, Atmos. Chem. Phys., 16, 13465–13475,, 2016. 

Turner, D. B.: Workbook of Atmospheric Dispersion Estimates, available at: (last access: 5 June 2017), 1970. 

USA National Phenology Network: Plant and Animal Phenology Data, Data type: Site Pheometrics, 09/01/2016-11/31/2016 for Region: Colorado, USA-NPN, Tucson, Arizona, USA,, 2018. 

Velasco, E., Pressley, S., Allwine, E., Westberg, H., and Lamb, B.: Measurements of CO2 fluxes from the Mexico City urban landscape, Atmos. Environ., 39, 7433–7446,, 2005. 

Velasco, E., Perrusquia, R., Jiménez, E., Hernández, F., Camacho, P., Rodríguez, S., Retama, A., and Molina, L. T.: Sources and sinks of carbon dioxide in a neighborhood of Mexico City, Atmos. Environ., 97, 226–238,, 2014. 

Verhulst, K. R., Karion, A., Kim, J., Salameh, P. K., Keeling, R. F., Newman, S., Miller, J., Sloop, C., Pongetti, T., Rao, P., Wong, C., Hopkins, F. M., Yadav, V., Weiss, R. F., Duren, R. M., and Miller, C. E.: Carbon dioxide and methane measurements from the Los Angeles Megacity Carbon Project – Part 1: calibration, urban enhancements, and uncertainty estimates, Atmos. Chem. Phys., 17, 8313–8341,, 2017. 

Waxman, E. M., Cossel, K. C., Truong, G.-W., Giorgetta, F. R., Swann, W. C., Coburn, S., Wright, R. J., Rieker, G. B., Coddington, I., and Newbury, N. R.: Intercomparison of open-path trace gas measurements with two dual-frequency-comb spectrometers, Atmos. Meas. Tech., 10, 3295–3311,, 2017. 

White, W. H., Anderson, J. A., Blumenthal, D. L., Husar, R. B., Gillani, N. V., Husar, J. D., and Wilson, W. E.: Formation and transport of secondary air pollutants: ozone and aerosols in the St. Louis urban plume, Science, 194, 187–189,, 1976. 

Wong, C. K., Pongetti, T. J., Oda, T., Rao, P., Gurney, K. R., Newman, S., Duren, R. M., Miller, C. E., Yung, Y. L., and Sander, S. P.: Monthly trends of methane emissions in Los Angeles from 2011 to 2015 inferred by CLARS-FTS observations, Atmos. Chem. Phys., 16, 13121–13130,, 2016. 

Wong, K. W., Fu, D., Pongetti, T. J., Newman, S., Kort, E. A., Duren, R., Hsu, Y.-K., Miller, C. E., Yung, Y. L., and Sander, S. P.: Mapping CH4:CO2 ratios in Los Angeles with CLARS-FTS from Mount Wilson, California, Atmos. Chem. Phys., 15, 241–252,, 2015. 

Wunch, D., Wennberg, P. O., Toon, G. C., Keppel-Aleks, G., and Yavin, Y. G.: Emissions of greenhouse gases from a North American megacity, Geophys. Res. Lett., 36, L15810,, 2009. 

Wunch, D. Toon, G. C., Sherlock, V., Deutscher, N. M., Liu, C., Feist, D. G., and Wennberg, P. O.: Documentation for the 2014 TCCON Data Release, CaltechDATA,, 2017. 

Short summary
In this work we measure carbon dioxide for the city of Boulder, Colorado, using a novel laser-based instrument. We then use a model to determine the strength of the emissions from the city. Based on our measurement location, we attribute the majority of these emissions to vehicles. We compare our emissions to the city vehicle emissions inventory with good agreement.
Final-revised paper