Domain Editor-in-Chief: Joel D. Blum; Department of Earth & Environmental Sciences, University of Michigan, Ann Arbor, Michigan, United States
Guest Editor: David Long; Michigan State University, East Lansing, United States

## 1 Introduction

While it is clear that fossil fuel combustion and land use change are driving observed increases in global atmospheric CO2 concentrations (e.g., Conway et al., 1994; Ciais et al., 1995; Keeling et al., 1995; Miller et al., 2012; Francey et al., 2013; Tans, 2014), only 40%–50% of emitted CO2 remains in the atmosphere due to uptake by oceanic and terrestrial sinks (Canadell et al., 2007; Knorr, 2009). Quantitatively modeling uptake of fossil fuel CO2 by the oceanic and terrestrial sinks is challenging because seasonal changes in the biosphere (Conway et al., 1994; Miller et al., 2012), the El Niño-Southern Oscillation cycle (Keeling et al., 1995), and other phenomena drive diurnal and seasonal variations in atmospheric CO2 concentrations. Quantifying the balance between various CO2 sources and sinks has been done most successfully at the global and hemispheric scales, usually across annual time periods (Ciais et al., 1995; Francey et al., 2013).

Even greater challenges exist to quantify carbon sources and sinks at regional and local scales, primarily because the processes occurring at such scales are less well-characterized. At the local scale, fewer studies have focused on urban centers compared to natural and agricultural systems (Peters and McFadden, 2012), in part due to the inherent complexity of urban regions. This knowledge gap is notable given that urban regions contribute ∼70% of anthropogenic CO2 emissions (Churkina, 2008) and contain the majority of the human population (53.6%, UN, 2014). Because the percentage of the global population living in urban centers is projected to increase to 66.4% by 2050 (UN, 2014), presumably with a commensurate increase in anthropogenic CO2 emissions, studies focused on urban carbon cycling are essential for more accurately predicting future climate change forcings and feedbacks. To parameterize urban carbon cycle models, it is necessary to analyze the spatial and temporal variation of local carbon sources (fossil fuel emissions and ecosystem respiration) and sinks (photosynthetic uptake) and to understand how local phenomena relate to larger regional and global trends. Establishing a comprehensive urban carbon cycle model will require data from cities and metropolitan regions representing diverse climates, population densities, and fossil fuel usage patterns (Pataki et al., 2006a; Pataki et al., 2007; Nordbo et al., 2012).

Measuring atmospheric CO2 concentrations ([CO2]) and/or CO2 fluxes at ≤1 minute intervals provides one method for quantifying carbon cycling at daily and seasonal time scales. This approach has been used for decades to study local carbon cycling in forested ecosystems (e.g., Keeling, 1958, 1961; Flanagan et al., 1996; Bowling et al., 2003; Pataki et al., 2003b; Bowling et al., 2005, 2009; Andrews et al., 2014). When collected at multiple, widely separated locations, high-frequency measurements can be employed in conjunction with model simulations to examine carbon cycling at regional to continental scales (e.g., Ciais et al., 1995; Ballantyne et al., 2011; Miller et al., 2012; Francey et al., 2013). Recently, an increasing number of studies have used high-frequency measurements to characterize urban carbon cycling (Reid and Steyn, 1997; Grimmond et al., 2002; Nemitz et al., 2002; Takahashi et al., 2002; Pataki et al., 2003a, 2006b; Vogt et al., 2006; Coutts et al., 2007; George et al., 2007; Miyaoka et al., 2007; Bergeron and Strachan, 2011; Christen et al., 2011; Helfter et al., 2011; Rice and Bostrom, 2011; Büns and Kuttler, 2012; Contini et al., 2012; Lietzke and Vogt, 2013).

At or near the surface in urban areas, average [CO2] are typically highest in densely populated centers, lower in neighboring suburbs, and near background levels in surrounding rural areas (Idso et al., 2001; George et al., 2007; Pataki et al., 2007; Büns and Kuttler, 2012). Additionally, high [CO2] can extend 100s of meters above the surface and downwind of urban areas (Mays et al., 2009; Cambaliza et al., 2014; Turnbull et al., 2015). Both anthropogenic activity and natural processes influence seasonal and daily carbon cycles in urban ecosystems, with resulting variations in urban [CO2] and CO2 fluxes (e.g., Reid and Steyn, 1997; Grimmond et al., 2002; Pataki et al., 2007; Helfter et al., 2011). Land use also affects average [CO2] in urban areas (Henninger and Kuttler, 2010).

However, measuring [CO2] alone provides limited information about carbon sources and sinks. Net CO2 fluxes can be calculated, but contributions from fossil fuels and ecosystem processes cannot be distinguished (Turnbull et al., 2006; Miller et al., 2012). Complementary analyses of other gases, such as CO or NOx, or application of the 14C tracer, can be used to estimate CO2 contributions from fossil fuel combustion versus ecosystem respiration (Turnbull et al., 2006; Miller et al., 2012; Lopez et al., 2013). Stable carbon isotope ratios (δ13C) of CO2 can also identify contributions from photosynthesis, soil respiration, petroleum combustion, and natural gas combustion in urban areas (Clark-Thorne and Yapp, 2003; Pataki et al., 2003a; Widory and Javoy, 2003; Pataki et al., 2005, 2006b; Bush et al., 2007; Pataki et al., 2007; Newman et al., 2008; Rice and Bostrom, 2011; Wada et al., 2011; Górka and Lewicka-Szczebak, 2013; Newman et al., 2013). In most previous studies, however, the number of δ13C measurements has been constrained to several data points per day because air must be manually collected in flasks and subsequently processed through vacuum extraction or gas chromatograph systems before analysis using a gas source isotope ratio mass spectrometer (Clark-Thorne and Yapp, 2003; Pataki et al., 2003a, 2005, 2007; Newman et al., 2008; Górka and Lewicka-Szczebak, 2013). Thus, while broad patterns for urban carbon cycles have begun to emerge, long-term, continuous records for [CO2] and δ13C in urban areas are relatively scarce. In North America, continuous records of urban [CO2] that are six months or longer have only been published for Baltimore, MD, USA (George et al., 2007), Phoenix, AZ, USA (Idso et al., 2002), and Salt Lake City, UT, USA (Pataki et al., 2003a, 2007; McKain et al., 2012). Worldwide, δ13C values collected for at least six months only exist for Dallas, TX, USA (Clark-Thorne and Yapp, 2003), Los Angeles, CA, USA (Newman et al., 2008), Salt Lake City (Pataki et al., 2003a, 2007), and Wroclaw, Poland (Górka and Lewicka-Szczebak, 2013). The only published continuous record of urban [CO2] for the midcontinent of North America is from Chicago, IL, USA, for June through August, 1995 (Grimmond et al., 2002). Many factors influence urban [CO2], e.g., seasonal climate variations, ecosystem type, fossil fuel usage patterns, regional topography, and population density. Several of these factors are likely distinct for the midcontinent compared to other previously studied cities.

The advent of wavelength scanned-cavity ring down spectrometers (WS-CRDS) now permits the simultaneous analysis of [CO2] and δ13C values at rates of >10,000 data points per day. Application of WS-CRDS is rapidly growing (Wahl et al., 2006; Rice and Bostrom, 2011; Lauvaux et al., 2012; Richardson et al., 2012; Newman et al., 2013; Sturm et al., 2013; Vogel et al., 2013). In the present study, we used a Picarro G1101-i WS-CRDS to generate a nearly continuous seven-month record for [CO2] and δ13C values for the Chicago region, which is the third largest urban center in the United States (U.S. Census Bureau, 2010). The record comprises ∼2x106 paired measurements, consolidated into hourly averages, and also includes complementary meteorological data obtained from the Illinois Environmental Protection Agency and the National Weather Service. The only similar study focused on [CO2] in Chicago was conducted in 1995 (Grimmond et al., 2002). Our dataset represents one of the longest records of [CO2] and δ13C published for any city worldwide. By incorporating information for various background locations and other cities, we use the dataset to examine how natural and anthropogenic sources of CO2 in a midcontinent urban setting change on a seasonal to daily basis as a function of key variables, such as wind speed, wind direction, and fossil fuel usage patterns.

## 2 Methods

### 2.1 Instrumentation and sample site location

In August 2011, a Picarro G1101-i WS-CRDS was installed on the roof of Scott Hall on the Northwestern University campus in Evanston, IL (Fig. 1, N42° 03′ 06.36″, W87° 40′ 39.16″). The instrument was operated nearly continuously through February 2012. Identifying the optimal location for measuring meteorological conditions and sampling gases in urban areas involves tradeoffs (World Meteorological Organization, 2008), such as maximizing the sampling height while maintaining accessibility for equipment maintenance. Urban areas can exhibit large diurnal variations in [CO2] (Reid and Steyn, 1997; Grimmond et al., 2002; Pataki et al., 2006b; Helfter et al., 2011; Rice and Bostrom, 2011; Büns and Kuttler, 2012). However, compared to the large range in diurnal variations, only relatively small concentration differences exist as a function of measurement height within the roughness sublayer when urban air samples are collected more than a few meters above the ground surface (Büns and Kuttler, 2012; Lietzke and Vogt, 2013). In the present study, the air intake was located on the rooftop of Scott Hall, ∼18 m above street level. Scott Hall is located ∼0.6 km west of Lake Michigan and sits at the intersection of two major roads (Fig. 1C). Evanston is a suburb of Chicago with 74,486 residents and a population density of 3687 persons km−2, ∼11% higher than the population-weighted density for the Chicago metropolitan area (U.S. Census Bureau, 2010). The city of Chicago, immediately south of Evanston, has ∼2.7 million residents (U.S. Census Bureau, 2010). The Chicago metropolitan area contains ∼9.5 million residents (U.S. Census Bureau, 2010). The Scott Hall sampling site is 3.6 km north of the Evanston–Chicago boundary and ∼20 km north of downtown Chicago (Fig. 1B).

doi: 10.12952/journal.elementa.000052.f001.
Figure 1.

#### Location maps.

Location maps: (A) Chicago within the USA, (B) the sampling site in relationship to downtown Chicago, and (C) the Northwestern University CO2 (NUCO2) sampling site, the IL EPA site where wind direction and wind speed data were collected, and the Northwestern University steam plant, all in Evanston, IL.

The Chicago, IL region has a continental climate with warm–hot summers and cold winters. From 1981 to February 2014, the mean annual temperature at O’Hare International Airport, ∼20 km southwest of the sampling site, was 9.9°C, with average temperatures of 22.1, 11.4, and 3.1°C in the summer (June–August), fall (September–November), and winter (December–February) (National Weather Service, 2014). Average annual precipitation during the same time period was 937 mm with 306, 242, and 147 mm in the summer, fall, and winter (National Weather Service, 2014). During the study period, the average temperatures in summer, fall, and winter were 23.3, 12.2, and 0.44°C, and total precipitation was 484, 225, and 156 mm (National Weather Service, 2014).

### 2.2 Sampling, calibration, quality control, and data treatment

The air intake and gas transfer system was modified from a NCAR design (Stephens and Watt, 2006). A vacuum pump (KNF 811 KTP) containing a Teflon-coated diaphragm drew ambient air through a 30 µm filter attached at the inlet. Air was pumped ∼10 m to a manifold using Synflex metal/plastic composite tubing with aluminum lining. Ambient air, as well as air from quality control check standards, flowed through two Nafion dryers (Perma Pure, MD–110–96F–4) before entering the instrument and then flowed out of the system through three molecular sieve dryers (Scott Specialty Gases [now Air Liquide America Specialty Gases], Model 44R, 13X molecular sieve with 4X indicator). Consequently, sample air had a low and relatively constant water content, <0.13% on average, with August–September values typically around 0.10%–0.15% and October–February values typically around 0.05%–0.10%, and the standard air was slightly humidified to the same level as sample air. Consistent water content increases the precision of WS-CRDS carbon isotope measurements (Rella, 2011). CO2 concentrations and δ13C values were measured ∼8 times min−1.

To calibrate the instrument, two NOAA gas standards with certified CO2 concentrations tied to the WMO scale and δ13C values (relative to VPDB) determined by INSTAAR at University of Colorado were analyzed ∼60 times over several months (Table S1). Two check standards were repeatedly analyzed in June and July 2011 to further confirm the calibration and assess instrumental stability. The check standards were tanks of CO2 in synthetic air (N2 + O2) from commercial suppliers (Airgas or Scott Specialty Gases [now Air Liquide America Specialty Gases]). During the data collection period, the NOAA gas standards were periodically analyzed to confirm the quality of the calibration. Within instrumental uncertainty, the measured [CO2] and δ13C values agreed with the NOAA and INSTAAR values. Also during the data collection period, two check standards similar to the ones described above were analyzed once or twice daily, with three minutes for flushing and five minutes for measurement (Table S1). Measured [CO2] for the check standards agreed with the supplier values and showed little evidence of drift (Table S1, Fig. S1). Flask samples were collected from each check standard for analysis of δ13C values by IRMS in the SIRFER lab at the University of Utah (Table S1). The δ13C values measured with the G1101-i showed little evidence for drift beyond normal instrumental uncertainty but were 0.54‰ ± 0.15‰ (1σ) lower than those determined by IRMS (Table S1). This difference is similar to the 0.43‰ offset observed in a previous study and can be attributed to the lack of Ar in the synthetic balance air (Friedrichs et al., 2010). The NOAA standards contain Ar. Because isotopic analysis of the synthetic-air check standards was used only to monitor instrumental drift, the offset has no bearing on the measured atmospheric values.

Data reduction and analysis were performed using R (R Core Team, 2013), implemented with the standard functionality, as well as the Lubridate (Grolemund and Wickham, 2011) and Openair libraries (Carslaw and Ropkins, 2012; Carslaw, 2013). Openair was employed for atmospheric analysis and plotting, including calculation and plotting of hourly averages for diurnal plots and polar plots. Hourly averages are reported on a 24–hour period, on the basis of local time – Central Daylight Time from August to Sunday, November 6, 2011 and Central Standard Time from Sunday, November 6 to the end of the sampling period in February 2012.

### 2.3 Ambient meteorological data

Wind speed and direction were collected by the Illinois Environmental Protection Agency (IL EPA) at a meteorological station located ∼1.1 km from the study site (N 42°03′42.61″, W87°40′28.84″, Fig. 1C). Temperature and other weather data were recorded at O’Hare International Airport by the National Weather Service (NWS). Wind speed and direction data from O’Hare were used to fill the following intervals when the IL EPA site did not collect data: August 1, 2011 until August 9, 2011 at 14:00, and November 4 until December 8 (with the exception of November 28, 15:00 to November 29, 23:59). Wind speed and wind direction data from the NWS were reported in bins of ∼0.4–0.6 m s−1 and 10°, respectively. The IL EPA wind direction data were partitioned into the same bins as the NWS data. The IL EPA and the O’Hare wind directions and speeds agreed well, with O’Hare reporting a 13° bias to the west and higher wind speeds by 0.3 m s−1 (see Fig. S2).

### 2.4 Keeling plots and mixing calculations

Measured [CO2] and δ13C values were used to construct Keeling plots (Keeling, 1958, 1961; Pataki et al., 2003b). As described by Pataki et al. (2003b), the measured atmospheric CO2 concentration ([CO2]m) represents the sum of the background CO2concentration ([CO2]b) and the concentration of CO2 emitted locally ([CO2]L):

${\left[{\mathrm{CO}}_{2}\right]}_{m}={\left[{\mathrm{CO}}_{2}\right]}_{b}+{\left[{\mathrm{CO}}_{2}\right]}_{L}$
(1)

Following the conservation of mass and the relationship between concentrations and isotope ratios (Pataki et al., 2003b), the δ13C value of local CO2 sources is calculated with the equation:

${\delta }^{\mathit{13}}{C}_{m}={\left[{\mathrm{CO}}_{2}\right]}_{b}\left({\delta }^{\mathit{13}}{C}_{b}-{\delta }^{\mathit{13}}{C}_{L}\right)\left(1/{\left[{\mathrm{CO}}_{2}\right]}_{m}\right)+{\delta }^{\mathit{13}}{C}_{L},$
(2)

where δ13Cm, δ 13CL, and δ13Cb represent the carbon isotope composition of measured, local, and background CO2, respectively. Because δ13CL is the intercept in a linear equation, the value can be calculated from regression of δ13Cm versus 1/[CO 2]m (Keeling, 1958, 1961; Pataki et al., 2003b). Ordinary least squares regression was used because the method introduces smaller uncertainties for calculated δ13CL values than other regression methods (Zobitz et al., 2006).

The δ13CL values were used to estimate relative contributions from various local CO2 sources according to equations presented in Newman et al. (2008):

${f}_{r}+{f}_{p}+{f}_{g}=1$
(3)
${\delta }^{\mathit{13}}{C}_{L}={f}_{r}{\delta }^{\mathit{13}}{C}_{r}+{f}_{p}{\delta }^{\mathit{13}}{C}_{p}+{f}_{g}{\delta }^{\mathit{13}}{C}_{g},$
(4)

where fr, fp, and fg represent the fraction of [CO2]L from ecosystem respiration, petroleum combustion (almost exclusively gasoline in Evanston), and natural gas combustion, respectively, and δ13Cr, δ13Cp, and δ13Cg represent the corresponding carbon isotope compositions.

Data from the NOAA Ameriflux tall tower site at Park Falls, WI, USA, located ∼500 km northwest of Evanston, were used to estimate background end members (Andrews et al., 2014). Average diurnal, monthly, and seasonal [CO2] (Tables 1 and 2) were calculated from hourly [CO2] data for the Park Falls 30 and 396 m inlets (Andrews et al., 2014). Monthly average δ13C values from Park Falls 396 m were used for Keeling plots and mixing calculations (White and Vaughn, 2011). Carbon isotope data are not available for Park Falls 30 m. Monthly reported [CO2] and δ13C values for Mauna Loa, HI, USA (Tans, 2014) were also used to compare results from this study with previous work (Table 1). Table 2 contains literature data for average [CO2] and δ13C values for other mixing end-members, namely petroleum and natural gas combustion, as well as ecosystem respiration, which includes plant and soil respiration (Bakwin et al., 1998; Widory and Javoy, 2003).

10.12952/journal.elementa.000052.t001

Table 1.

#### Monthly averages for background CO2 and δ13Ca

Park Falls, WI - 30 m Park Falls, WI - 396 m Mauna Loa, HI
Month CO2 (ppm) δ13C (‰) CO2 (ppm) δ13C (‰) CO2 (ppm) δ13C (‰)
August 2011 388.8 nmb 376.3 -7.75 390.0 -8.24
September 2011 391.7 nm 386.1 -8.04 388.6 -8.17
October 2011 398.4 nm 394.3 -8.40 389.2 -8.19
November 2011 400.8 nm 399.4 -8.64 390.5 -8.26
December 2011 401.9 nm 400.6 -8.70 392.2 -8.32
January 2012 401.7 nm 400.6 -8.73 393.0 -8.35
February 2012 401.5 nm 400.7 -8.77 393.3 -8.38
Whole period average 397.9 nm 394.0 -8.43 390.9 -8.27
10.12952/journal.elementa.000052.t002

Table 2.

#### End members for CO2 and δ13C in Keeling plots

CO2 (ppm) 10,000 *1/CO2 (ppm−1) δ13C (‰) Sourceb
Background - Park Falls, WI (396 m)
Aug-Sep 2011 381.1 26.2 -7.89 A
Oct–Nov 2011 396.8 25.2 -8.52 A
Dec 2011–Feb 2012 400.6 25.0 -8.73 A
CO2 sources
Natural gas 7757.1 1.3 -39.1 B
Petroleum 32340 0.3 -28.7 B
Respirationa 7757.1 1.3 -24.7 C

## 3 Results and discussion

Evanston consistently displayed higher [CO2] and lower δ13C values than the background sites. Below, measurements from Evanston and the background sites are compared in progressively shorter time periods, from the entire sample period to seasonal, and finally, diurnal. The role of wind speed and direction is also discussed. Finally, seasonal variations in local CO2 sources are quantified using Keeling plots and mass-balance mixing equations.

### 3.1 Average [CO2] and δ13C compared to background

The average [CO2] in Evanston from August 2011–February 2012 was 16, 20, and 23 ppm higher than the Park Falls 30 m, Park Falls 396 m, and Mauna Loa background sites, respectively (Tables 1, 3). The average [CO2] difference between Evanston and Mauna Loa determined here is similar to the 22 ppm difference reported in a previous study conducted in Chicago during June–August, 1995 (Grimmond et al., 2002). The differences between Evanston and the background sites are in the middle of the range reported in other studies that compared urban and background [CO2]. At the low end, [CO2] in Portland, Oregon, USA was only ∼5 ppm higher than a nearby rural area (Rice and Bostrom, 2011). Cities with [CO2] ranging from 15–35 ppm higher than background, similar to Evanston, include Essen, Germany (Büns and Kuttler, 2012); surburban or foothill sites in Salt Lake City (Pataki et al., 2005, 2007); Tokyo, Japan (Moriwaki et al., 2006); and Vancouver, Canada (Reid and Steyn, 1997). Concentration differences of 50 ppm or more between urban and background sites were observed in densely populated sections of Los Angeles, downtown Salt Lake City, and Phoenix, all cities surrounded by mountains (Idso et al., 2001; Pataki et al., 2007; Newman et al., 2008) and in densely populated Baltimore, where samples were collected only 2 m above ground level (George et al., 2007). These data suggest that the difference between urban and background [CO2] observed in Evanston is fairly typical, although some cities have smaller or larger differences depending on the extent of local emissions, local meteorological conditions, or geological/geographic characteristics.

10.12952/journal.elementa.000052.t003

Table 3.

#### Monthly averages for CO2, δ13C, and meteorological data

Month CO2 (ppm) 1 σ δ13C (‰) 1 σ n T (°C) Wind speed (m/s) Wind direction (degrees) Wind direction
August 2011 397.3 18.6 -8.29 1.0 680 23.1 3.01 314 NW
September 2011 402.9 20.3 -8.65 0.9 653 16.8 3.80 337 NNW
October 2011 417.8 23.0 -9.72 1.2 682 12.8 3.53 256 WSW
November 2011 414.1 13.4 -9.64 0.7 660 7.2 5.40 245 WSW
December 2011 419.0 12.2 -9.97 0.7 618 2.0 3.57 258 WSW
January 2012 426.9 22.2 -10.02 1.2 503 -0.7 3.65 269 W
February 2012 419.6 14.2 -9.54 0.8 667 0.8 3.64 275 W
Whole period average 413.9 10.3 -9.40 0.7 4463 10.5 3.80 314 NW

The average δ13C value in Evanston was 0.97‰ and 1.13‰ lower than Park Falls 396 m and Mauna Loa, respectively (Tables 1, 3). Several studies have reported that δ13C values in urban air are lower than background values (Clark-Thorne and Yapp, 2003; Pataki et al., 2003a; Widory and Javoy, 2003; Pataki et al., 2005, 2006a, 2006b, 2007; Rice and Bostrom, 2011; Górka and Lewicka-Szczebak, 2013; Lopez et al., 2013). The only previous study to report and compare average urban and background δ13C values over a several-month period found that urban δ13C values were 1.7‰ lower than background in the early 2000s (Newman et al., 2008). The Newman et al. (2008) study was conducted in Pasadena, California, USA, which receives marine air blowing directly off the Pacific Ocean, and thus, different results are expected for a midcontinent setting where seasonal changes in ecosystem productivity yield seasonal differences in δ13C values.

### 3.2 Seasonal [CO2] and δ13C variations

Mean [CO2] was lowest during the late summer and early fall (August–September) and highest during the winter (December–February) (Fig. 2, Table 3). The ∼22 ppm difference between the [CO2] minimum in August–September and the maximum in December–February was larger for Evanston (Table 3) than the differences for the background sites: 19, 13, and 4 ppm for Park Falls 396 m, Park Falls 30 m, and Mauna Loa, respectively (Table 1). The overall seasonal pattern for Evanston is consistent with the Park Falls background site (Table 1), as well as broader northern hemisphere [CO2] trends, which are driven by increased photosynthesis during summer and decreased photosynthesis during winter (Lauvaux et al., 2012; Andrews et al., 2014; Tans, 2014). The seasonal pattern is also consistent with measurements made in other northern hemisphere urban areas (Pataki et al., 2003a, 2007; Kordowski and Kuttler, 2010; Helfter et al., 2011).

doi: 10.12952/journal.elementa.000052.f002.
Figure 2.

#### Hourly CO2 concentrations and δ13C values across sampling period.

Hourly (A) CO2 concentrations ([CO2]) and (B) δ13C values from August 2011 to February 2012. Gaps represent downtime for instrument maintenance. The thick line represents a locally-weighted smoothed trend. Diurnal differences occur as high-frequency, large amplitude changes. Seasonal differences are highlighted by the smoothed line, with lower [CO2] in the summer and higher [CO2] in the winter. δ13C values show the opposite trend. Monthly (C) [CO2] and (D) δ13C values where the dark line represents the median, the box represents the interquartile range (IQR, 25%–75%), the whiskers represent all values within 1.5 IQR, and the points represent values outside 1.5 IQR.

During the summer, [CO2] in Evanston was ∼10 ppm higher than Park Falls 30 m, which is roughly the same height as the Scott Hall station. During the winter, [CO2] was ∼20 ppm higher in Evanston. The smaller difference during summer suggests that photosynthesis and respiration strongly influence [CO2] in both rural and urban midcontinent settings, while the larger difference during winter points to the importance of urban fossil fuel emissions (see section 3.5). A similar seasonal pattern was observed in Salt Lake City, one of the few urban locations with high resolution winter data (Pataki et al., 2007).

The δ13C values for Evanston and the background sites inversely correlate with [CO2], with a summer maximum and a winter minimum (Fig. 2, Tables 1, 3). As with [CO2], the difference between summer and winter δ13C values was larger in Evanston and Park Falls 396 m compared to Mauna Loa. The greater amplitude of the seasonal δ13C variation in Evanston (1.6‰) than Park Falls 396 m (1.0‰) is also consistent with CO2 contributions from urban fossil fuel emissions, including differing emission sources between summer and winter.

### 3.3 Diurnal [CO2] and δ13C variations

Across the entire sampling period, average [CO2] and δ13C values varied diurnally in Evanston, with low [CO2] and high δ13C values in the afternoon, increasing [CO2] and decreasing δ13C values during the evening, and high [CO2] and low δ13C values in the early morning (Fig. 3). The average diurnal variation was ∼20 ppm for [CO2] and 1‰–1.1‰ for δ13C values. Evanston appears similar to other urban areas, where average [CO2] are lowest during the daytime hours, when a relatively turbulent atmosphere mixes surface air containing local CO2 emissions with background air containing lower [CO2] (Reid and Steyn, 1997; Grimmond et al., 2002; Pataki et al., 2006b; Helfter et al., 2011; Rice and Bostrom, 2011; Büns and Kuttler, 2012). The same pattern also occurs in rural areas (Flanagan et al., 1996; Bowling et al., 2003, 2005). During evening through early morning hours, the [CO2] rises because the atmosphere is generally more stable (Reid and Steyn, 1997; Grimmond et al., 2002; Pataki et al., 2006b; Helfter et al., 2011; Rice and Bostrom, 2011; Büns and Kuttler, 2012). On an hourly basis, urban CO2 fluxes show less variation than [CO2] (Grimmond et al., 2002) and sometimes exhibit complex patterns. For example, in some urban areas, the highest hourly CO2 fluxes occur in the afternoon, when [CO2] are lowest (Vogt et al., 2006; Helfter et al., 2011).

doi: 10.12952/journal.elementa.000052.f003.
Figure 3.

#### Diurnal variations in CO2 concentrations and δ13C values for each season.

Hourly mean values for CO2 concentrations ([CO2]) and δ13C values during (A, D) summer, (B, E) fall, and (C, F) winter. [CO2] averages shown in grey are background measurements from 30 m height at the background site in Park Falls, WI, USA for (A) summer, (B) fall, and (C) winter. The bold line shows hourly means, and the shaded regions represent 95% confidence intervals. Evanston [CO2] increases, and δ13C values decrease, from summer to winter. The amplitude of the average diurnal variation in Evanston decreases from summer to winter. [CO2] at the background site also increases from summer to winter, although to a smaller degree. Average diurnal variations for Evanston and Park Falls are similar in the summer (25–30 ppm) but much different in the winter: ∼15 ppm for Evanston and <1.25 ppm for Park Falls.

While Evanston and other urban areas exhibit diurnal variations throughout the year, the amplitude and timing of the variations change seasonally. Understanding combinations of anthropogenic and natural factors that drive seasonally-changing diurnal variations is important for elucidating why urban carbon cycles vary throughout the year. These factors are discussed in more detail in the following sections.

#### 3.3.1 Summer diurnal variations

Summer diurnal variations for [CO2] and δ13C values approached 30 ppm and 1.5‰ (Fig. 3A, 3D). The highest [CO2] and lowest δ13C values occurred between 6:00 and 7:00, while the lowest [CO2] and highest δ13C values occurred from 12:00–18:00 (Fig. 3A, 3D). The Park Falls 30 m background site displayed a 25 ppm diurnal range, with [CO2] peaking in the morning and rising in the evening (Fig. 3A). CO2 concentrations for the Park Falls 396 m site exhibited a smaller diurnal range of 5 ppm because samples were collected above the well-mixed planetary boundary layer and thus represent a larger scale, regional average. The daily average and diurnal amplitudes for [CO2] were similar for summer weekdays and weekends (Fig. 4A). A weekday morning [CO2] peak that occurred an hour later than the weekend peak (Fig. 4A) and a somewhat larger diurnal [CO2] variation in Evanston compared to Park Falls 30 m are consistent with fossil fuel emissions, likely from weekday rush hour traffic.

doi: 10.12952/journal.elementa.000052.f004.
Figure 4.

#### Weekday versus weekend diurnal variations in CO2 concentrations for each season.

Hourly mean CO2 concentrations ([CO2]) for weekday (red) versus weekends (grey) for (A) summer, (B) fall, and (C) winter. The bold line shows hourly means, and the shaded regions represent 95% confidence intervals. In the summer, weekday and weekend hourly averages are indistinguishable. In the fall and winter, weekday hourly averages are higher than weekend averages.

Summer diurnal ranges observed in Evanston are similar to those reported for other northern hemisphere, mid-latitude cities. For example, the 30 ppm summer diurnal range for [CO2] in Evanston is similar to the 35 ppm diurnal range observed during the summer of 1995 in Chicago (Grimmond et al., 2002) and to the 20–30 ppm diurnal range measured in Basel, Switzerland (Vogt et al., 2006), London, England (Sparks and Toumi, 2010), Portland (Rice and Bostrom, 2011), Vancouver (Reid and Steyn, 1997), and a suburban site in Salt Lake City (Pataki et al., 2007). A morning [CO2] peak and an afternoon [CO2] low have been observed in other urban areas (e.g., Reid and Steyn, 1997; Grimmond et al., 2002; Vogt et al., 2006; Coutts et al., 2007; Pataki et al., 2007; Büns and Kuttler, 2012; Contini et al., 2012; Lietzke and Vogt, 2013). In Melbourne, Australia, the only urban study from the southern hemisphere, a smaller summer diurnal range of 8–12 ppm was observed (Coutts et al., 2007).

#### 3.3.2 Fall diurnal variations

Fall diurnal variations were ∼20 ppm for [CO2] and 1.2‰ for δ13C values (Fig. 3B, 3E). Compared to summer, the highest [CO2] and lowest δ13C values shifted later to 7:00–9:00. The shift partly reflects the change from daylight savings time to standard time, which occurred approximately halfway through the fall season. The afternoon period of low [CO2] was shorter relative to summer, and the rise toward higher evening [CO2] was steeper and began earlier, around 17:00 (Fig. 3B, 3E). Similar to summer, the fall morning [CO2] peak was synchronous with the peak recorded at the Park Falls 30 m site, but the late afternoon [CO2] rise in Evanston was earlier and more rapid (Fig. 3B).

While fall weekdays and weekends displayed the same average daily [CO2] within uncertainty, average hourly [CO2] for weekdays was higher than weekends from 11:00 to 19:00, and the late afternoon [CO2] rise was earlier and more abrupt on weekdays (Fig. 4B). Higher hourly [CO2] during the midday and late afternoon suggest contributions from local vehicular traffic. The late afternoon [CO2] increase likely reflects emissions during the evening rush hour (Fig. 4B). Based on 2010 measurements, average traffic volumes from three nearby measurement locations – the Evanston intersection where sampling occurred and two adjacent intersections located north and south of the sampling site–were 13,545 vehicles day−1, with 24% (3260 vehicles) of the daily traffic occurring from 16:00–19:00 (Illinois Department of Transportation, 2010). These vehicle counts are likely minimum estimates because the data were collected during the summer when Northwestern University was out of session, resulting in less traffic than during the academic year. The late afternoon [CO2] increase in Evanston was not evident during the summer (Fig. 3A, 4A), which could reflect several factors, including longer periods of daylight causing photosynthetic activity to continue through the early evening, a longer duration of atmospheric mixing, and less activity at the university compared to the academic year. The suggestion that photosynthetic drawdown of CO2 during the summer can be large enough to mask the evening rush hour signal highlights the significance of local plant growth in urban carbon cycling. Some studies attributed afternoon [CO2] increases to the end of the day rush hour (Lietzke and Vogt, 2013; Newman et al., 2013), whereas others did not observe an afternoon increase (Reid and Steyn, 1997; McKain et al., 2012). Because afternoon [CO2] increases are relatively small, their detection partly depends on measurement height. For example, in Basel, an afternoon increase was observed for measurements taken at 3 and 19 m above ground level but not at 39 m above ground level, where air is more well-mixed (Lietzke and Vogt, 2013).

#### 3.3.3 Winter diurnal variations

Winter diurnal variations were ∼15 ppm for [CO2] and 1‰ for δ13C values (Fig. 3C, 3F). Average daily [CO2] for weekdays was 2.7 ppm higher than weekends. Average hourly [CO2] was higher during weekdays than weekends during the time period 5:00–13:00 (Fig. 4C). Similar to the fall, the late afternoon [CO2] rise in winter was more abrupt during weekdays (Fig. 4C). About 42% of the total daily traffic passed during morning and afternoon rush hours: 2453 vehicles (18%) from 7:00–10:00 and 3260 vehicles (24%) from 16:00–19:00 (Illinois Department of Transportation, 2010). In contrast to Evanston, winter diurnal variations observed at both Park Fall sites were <1.25 ppm (Fig. 3C). The small diurnal variation at Park Falls demonstrates minimal CO2 contributions from ecosystem respiration, as expected for the winter season. The larger diurnal variations in Evanston compared to Park Falls therefore point to fossil fuel emissions as the major source of CO2. The increase in [CO2] from 0:00–6:00 almost certainly reflects natural gas emissions from building heating because vehicular traffic was low. Only ∼360 vehicles, or 2.6% of the daily total traffic, passed during this time (Illinois Department of Transportation, 2010). The timing and size of the weekday versus weekend morning peaks in Evanston (Fig. 4A, 4C), as well as the similarity of the summer and fall morning peaks in Evanston and Park Falls but near absence of a winter morning peak at Park Falls (Figs. 3AC), collectively suggest that local ecosystem respiration controls much of the size and timing of morning [CO2] in summer and that fossil fuel emissions control the morning [CO2] peak in winter.

The winter diurnal range in Evanston is similar to the 15–20 ppm diurnal range observed at several other northern hemisphere urban sites, such as Basel (Lietzke and Vogt, 2013); Essen (Büns and Kuttler, 2012); Lecce, Italy (Contini et al., 2012); and London (Sparks and Toumi, 2010). Similar to fall and winter patterns in Evanston, weekday peaks in Basel occurred later and were larger than weekend peaks (Lietzke and Vogt, 2013). Suburban and downtown sites in Salt Lake City exhibited average winter diurnal variations of 25–30 ppm and 40–50 ppm, respectively (Pataki et al., 2007). The larger diurnal variations recorded for Salt Lake City likely reflect the influence of mountains, which surround the city and reduce atmospheric mixing (Pataki et al., 2005).

Both urban and background sites exhibit diurnal variations, with the amplitude and timing changing seasonally. However, diurnal ranges are consistently larger in Evanston compared to the background sites, especially during the winter. The seasonal pattern of diurnal variations in Evanston is similar to the patterns documented for many other northern hemisphere cities. While local factors, like topography for Salt Lake City, can impact diurnal variations, the overall similarity between Evanston and other northern hemisphere cities suggests common controlling variables, such as meteorological conditions and seasonal changes in natural versus fossil fuel [CO2] sources. These factors are discussed in more detail in the following two sections.

### 3.4 [CO2] as a function of wind speed and direction

Although considerable variation exists, the highest [CO2] generally occurred at the lowest wind speeds regardless of season (Fig. 5A). As wind speeds increased, the concentration and extent of variation decreased. Daily variations in wind speed (horizontal advection of air) and the thickness of the well-mixed planetary boundary layer (vertical mixing) strongly influenced [CO2] and δ13C values (Fig. 3, 5A). For example, from 12:00 to 16:00, mean [CO2] was lowest (Fig. 3) while wind speed was highest (4.30 m s−1 on average). In contrast, from 03:00–07:00, [CO2] was highest (Fig. 3) while mean wind speed was lowest (3.44 m s−1 on average). The thickness of the well-mixed layer also affects [CO2] by changing the fraction of the atmosphere that receives direct inputs of CO2 emissions from the surface. In urban areas, the thickness of the well-mixed layer increases from sunrise until early to mid-afternoon and then decreases to a minimum at night (Helfter et al., 2011; Strong et al., 2011; Newman et al., 2013). In four cities located in, or around, the midwestern US, the thickness of the well-mixed layer was 300–500 m in the morning and 1700–1900 m in the afternoon (Holzworth, 1967). As the thickness of the well-mixed layer increases, low [CO2] air from higher in the atmosphere can be entrained, thereby lowering locally measured [CO2] (Holzworth, 1967; Strong et al., 2011). This type of diurnal relationship between [CO2] and atmospheric mixing is common in both urban and rural areas (e.g., Flanagan et al., 1996; Reid and Steyn, 1997; Grimmond et al., 2002; Bowling et al., 2003, 2005; Vogt et al., 2006; Coutts et al., 2007; Büns and Kuttler, 2012).

doi: 10.12952/journal.elementa.000052.f005.
Figure 5.

#### CO2 concentrations as a function of wind speed, direction, and time of day.

(A) CO2 concentrations ([CO2]) versus wind speed for summer, fall, and winter. [CO2] concentrations at zero wind speed are not shown. While variability exists, [CO2] are highest at low wind speeds and are lower and less variable at higher wind speeds. (B) Polar annulus plot of [CO2] versus wind direction and hour of the day. The 0 inside the annulus represents midnight, and 23 outside the annulus represents 23:00. Higher concentrations occur when winds are from the northwest, west, and south, with the highest concentrations occurring during the evening through the early morning hours. A secondary high occurs for southeast winds during the midday period.

Wind direction is also a critical factor. The correspondence between hourly [CO2] and wind direction demonstrates that CO2 emissions from nearby populated areas contributed to high [CO2]. In Evanston, densely populated areas are located to the northwest, west, and south, whereas Lake Michigan is located to the east and north (Fig. 1). Higher [CO2] occurred when winds blew from the west or south (Fig. 5B). Lower concentrations dominated when winds blew from the east or north over Lake Michigan (Fig. 5B). When winds blew from the west-northwest to southwest, [CO2] was highest at night and early morning (21:00 to 08:00) and lowest during the day (Fig. 5B). When winds blew from the east to southeast, the lowest concentrations occurred during the evening or midday (Fig. 5B). One exception occurred during the midday, when winds blew from the southeast over the nearby intersection of two major roads and the population center of Chicago.

Fundamental patterns emerge when [CO2] is examined as a function of both wind speed and direction. Figure 6 provides a polar plot, where mean [CO2] is plotted for each vector of wind speed and wind direction (Carslaw and Ropkins, 2012; Carslaw, 2013). The overall pattern with wind speed and direction was consistent across each of the three seasons (Fig. S3). Local [CO2] was at, or above, the average concentration for the entire measurement period when winds blew from the south to southwest at speeds of 10 m s−1 or less, or from the west at speeds of 7 m s−1 or less (Fig. 6). Additional small regions yielding above average concentrations were seen to the north-northwest and northeast at higher wind speeds (8–12 m s−1, Fig. 6). The higher concentrations to the north-northwest likely reflect winds blowing over nearby neighborhoods and suburbs. Higher concentrations were recorded when winds blew from the northeast at velocities of 9–14 m s−1 (Fig. 6), which could reflect CO2 emissions from the Northwestern University steam plant, ∼0.45 km northeast of the sampling site (Fig. 1). In 2010, the plant had a capacity of 180 million BTU hour−1 in the most commonly used boiler and 125 million BTU hour−1 in secondary boilers that supplemented output during high demand periods, especially during winter (personal communication, Northwestern University Facilities and Maintenance). Based on estimates of the total CO2 emissions within Evanston (City of Evanston, 2008) and the US Energy Information Administration coefficient of 53.2 kg CO2 emitted per million BTU for natural gas, the steam plant emits 17%–26% of the CO2 within the Evanston borders.

doi: 10.12952/journal.elementa.000052.f006.
Figure 6.

#### Mean CO2 concentration as a function of wind direction and speed.

Polar plot of [CO2] versus wind speed and direction for August 2011–February 2012. Each pixel represents the mean [CO2] for each given wind speed and direction. The wind speed increases with distance from the center. The white areas represent wind speed and directions with little or no data. The highest [CO2] are observed at wind speeds <3–5 m s−1 and when winds are from the south or west. Above-average [CO2] are also observed at wind speeds of 10–15 m s−1 when the wind blows from the northeast, which is the direction of the natural gas-fired university central utility plant. Below-average [CO2] are generally observed at higher wind speeds and also when the wind blows from the east off Lake Michigan.

### 3.5 CO2 sources and sinks

Keeling plots showing δ13C values versus 1/[CO2] were used to estimate the fraction of CO2 contributed from fossil fuel and plant respiration sources. Summer data plot between the mixing lines for ecosystem respiration and petroleum combustion (Fig. 7A). Local CO2 drawdown and attendant carbon isotope fractionation by photosynthesis clearly occurred during summer, as many δ13C values are more positive than the background value (Fig. 7A). Fall and winter data plot between the mixing lines for petroleum and natural gas combustion, with fall data plotting closer to the petroleum line and winter data plotting closer to the natural gas line (Figs. 7B and 7C). Several studies have identified local fossil fuel CO2 emissions as the likely source of low δ13C values measured for urban CO2 (Clark-Thorne and Yapp, 2003; Pataki et al., 2003a; Widory and Javoy, 2003; Pataki et al., 2005, 2006a, 2006b, 2007; Rice and Bostrom, 2011; Górka and Lewicka-Szczebak, 2013; Lopez et al., 2013).

doi: 10.12952/journal.elementa.000052.f007.
Figure 7.

#### Keeling plots of δ13C versus 1/[CO2] for determining CO2 sources.

Keeling plots of δ13C versus 1/[CO2]x 10,000 (ppm−1) for (A) summer, (B) fall, and (C) winter. Note that higher [CO2] plot to the left. Dashed lines show the linear regression fit to the data for each season (see Table 5). The black dots represent average background values for each period measured at the Ameriflux tall tower site (396 m elevation) in Park Falls, WI (Andrews et al., 2014). Black solid lines represent mixing between the background value and either petroleum (upper line) or natural gas (lower line) combustion. The dotted line represents mixing between the background value and ecosystem respiration. Each tick mark along the petroleum and natural gas mixing lines represents a 10% contribution from the end-member source. For example, in the summer Keeling plot, the tick mark at approximately -10‰ and 23.75 indicates that 10% of the CO2 originates from a petroleum source and 90% originates from the background atmosphere. In summer, hours with CO2 in excess of background plot between the ecosystem respiration and petroleum mixing lines. In fall, hours with excess CO2 plots between the petroleum and natural gas mixing lines, and in winter, the hourly data are closer to the natural gas mixing line.

The 1/[CO2] and δ13C values highly correlate for multiple time periods, including each season, afternoon hours (12:00–6:00), and night hours (00:00–06:00 local time) (Table 4). Afternoon values for the slope and intercept are the same within statistical uncertainty across the sampling period. From summer to winter, the slope of the best–fit line became steeper, and the y-intercept, δ13CL, became more negative by 7‰ (Table 4, Fig. 7). The nighttime data reflect periods when the well-mixed layer was at its lowest elevation and advective transport was minimal. Several studies have used the δ13CL from nighttime data to determine the δ13C value of ecosystem respiration (Pataki et al., 2003b; Zobitz et al., 2006; Bowling et al., 2009) and to determine CO2 sources in urban areas (Pataki et al., 2003a, 2006b, 2007). Similarly, we used nighttime data to determine the δ13CL value for local CO2 sources representing a mixture of ecosystem respiration and fossil fuel emissions. Given that the Evanston dataset represents the longest continuous δ13C record yet published, with hundreds of nighttime data points for each season, the estimated δ13CL values are likely robust.

10.12952/journal.elementa.000052.t004

Table 4.

#### Linear fits for Keeling plots

Month Slope 1 σ δ13CL (intercept) 1 σ R2 n
All hours
Aug–Sep 2011 0.79 0.01 -28.3 0.2 0.88 1332
Oct–Nov 2011 0.91 0.01 -31.6 0.2 0.90 1342
Dec 2011–Feb 2012 1.05 0.01 -34.7 0.2 0.90 1788
Entire period 0.94 0.005 -32.2 0.1 0.90 4462
Afternoons (12:00–16:00)
Aug-Sep 2011 0.79 0.03 -28.1 0.8 0.75 196
Oct–Nov 2011 0.76 0.03 -27.9 0.6 0.82 208
Dec 2011–Feb 2012 0.82 0.03 -29.2 0.7 0.72 313
Entire period 0.88 0.01 -30.6 0.3 0.86 717
Nights (00:00–06:00)
Aug–Sep 2011 0.78 0.02 -28.1 0.4 0.89 312
Oct–Nov 2011 0.95 0.02 -32.8 0.4 0.90 330
Dec 2011–Feb 2012 1.08 0.01 -35.7 0.3 0.93 468
Dec 2011–Feb 2012, Easterly windsa 0.94 0.06 -32.2 1.5 0.76 76
Entire period 0.96 0.01 -32.7 0.2 0.89 1110
10.12952/journal.elementa.000052.t005

Table 5.

#### Scenarios for contributions of ecosystem respiration & fossil fuel CO2 sources

Month Measured nighttime δ13C value (‰) Fossil fuel δ13C valuea (‰) Petroleum contribution (%) Natural gas contribution (%) Total fossil fuel contribution (%) Ecosystem respiration contribution (%)
Summer
Scenario 1 – 100% natural gas -28.1 -39.1 0.0 100.0 23.7 76.3
Scenario 2 – 75% natural gas -28.1 -36.5 25.0 75.0 28.8 71.2
Scenario 3 – 50% natural gas -28.1 -33.9 50.0 50.0 37.0 63.0
Scenario 4 – 25% natural gas -28.1 -31.3 75.0 25.0 51.5 48.5
Scenario 5 – 0% natural gas -28.1 -28.7 100.0 0.0 85.0 15.0
Winter
Scenario 1 – 0% respiration -35.7 -35.7 32.7 67.3 100.0 0.0
Scenario 2 – 10% respiration -35.7 -36.9 20.9 79.1 90.0 10.0
Scenario 3 – 20% respiration -35.7 -38.5 6.2 93.8 80.0 20.0
Scenario 4 – 30% respiration -34.8b -39.1 0.0 100.0 70.0 30.0

Simple mixing model scenarios were used to constrain contributions from ecosystem respiration and different fossil fuel sources in the summer and winter seasons (Table 5). Both ecosystem respiration and building heating likely contributed CO2 during the fall. Therefore, we were unable to make simplifying assumptions necessary for computing mixing model scenarios for this season. The summer and winter scenarios were used to bracket the combinations of fossil fuel emissions and respiration inputs that reasonably reproduce the measured δ13CL values. Note that the statistical uncertainty of the linear regression (0.2‰–0.4‰) used to determine δ 13CL values yields a 3%–5% uncertainty in the contributions of different CO2 sources (Table 5).

For the summer, scenarios 1–5 show the respiration contribution across a range of fossil fuel contributions from 100% natural gas to 100% petroleum (Table 5). In scenarios 1–4, ecosystem respiration contributes at least 48.5% to δ13CL. However, scenario 1 and 2 are unlikely because the estimates for natural gas usage are too high for summer, yet scenario 5 also seems unlikely because even during summer, some natural gas combustion is necessary for heating water, particularly on campus and in dense residential areas surrounding the sampling location. Given that natural gas demand is relatively low during warm months, scenario 4, and to a lesser extent, scenario 3, are the most plausible, with plant respiration contributing 48.5%–63% to δ13CL (Table 5). An uncertainty of <15% seems reasonable given that uncertainties can exceed 30% in cases where a single source does not contribute at least 70% of the local CO2 (Pataki et al., 2007).

The summer contribution from ecosystem respiration in Evanston lies at the middle of the range estimated for other cities using δ13C, as well as CO and 14C. On the low end, respiration contributed 20%–30% of the excess CO2 during late spring in Nagoya, Japan (Wada et al., 2011) and 35%–50% in late summer in Salt Lake City (Pataki et al., 2007). Respiration in the Los Angeles basin has been estimated to be as low as 50% during the late spring (Newman et al., 2013) and as high as 75% for a mid-spring day at the end of the local rainy season (Djuricin et al., 2010). Evanston was similar to Dallas and Salt Lake City, in that petroleum combustion was the major fossil fuel contributor in the summer (Table 5, Clark-Thorne and Yapp, 2003; Pataki et al., 2007).

In winter, the ecosystem respiration component was almost certainly ≤20% and likely 0%–10% (Table 5). When ecosystem respiration is modeled as 30% in winter (scenario 4), then the modeled nighttime δ13CL value was 0.9‰ more positive than the measured value (Table 5). Low respiration is consistent with low winter temperatures, which slow metabolic processes. Diurnal [CO2] variations for Park Falls 30 m were only ∼1 ppm for the same period, which provides additional evidence for minimal respiration (Fig. 3C). Thus, the Evanston winter CO2 end member was dominated by contributions from local petroleum and natural gas combustion rather than respiration (Table 5). The signal from fossil fuel combustion in densely populated areas to the west and south can be seen by comparing the overall average winter nighttime δ13CL value of –35.7‰ with the δ13CL value of –32.2‰ when winds blew from the east over Lake Michigan (0–180 degrees, Table 4). Air coming from the direction of Lake Michigan contains less CO2 from fossil fuel emissions, resulting in a more positive δ13CL value. In winter, natural gas combustion contributed more than 67% of the excess CO2 in scenario 1 with no respiration and 76% in scenario 2 where respiration constituted 10% of the high [CO2] end member (Table 5). Scenario 3 is unlikely because the contribution of petroleum combustion to local CO2 almost certainly exceeds 6%, given the estimated results for summer. Large contributions from natural gas combustion occurred during winter, even though the temperature was slighter warmer than average, 0.44°C for the study period versus –3.1°C from 1981–2014 (National Weather Service, 2014). Scenarios suggesting large natural gas contributions during the winter and modest contributions during the summer are generally consistent with an emissions inventory for the broader Chicago region, where annual greenhouse gas emissions (primarily CO2) from natural gas usage were ∼50% higher than emissions from the transportation sector (McGraw et al., 2010). Similarly, a greenhouse gas inventory for Evanston found that emissions from natural gas combustion were 2.7 times higher than emissions from gasoline or diesel combustion (City of Evanston, 2008). The scenarios are also consistent with CO2 emissions estimates based on fuel consumption for the approximate latitude and longitude of Chicago, which indicate that emissions from natural gas combustion vary by 70%–80% between summer and winter, while emissions from petroleum combustion vary by <10% (Gregg et al., 2009). The net result is that overall emissions are highest during winter months and that natural gas combustion is the largest source (Gregg et al., 2009).

Our estimates for winter ecosystem respiration in Evanston are similar to other northern cities where winter temperatures are cold enough to require building heating: 11% in Nagoya, Japan (Wada et al., 2011), 5%–30% in Salt Lake City (Pataki et al., 2007), and 20%–30% in Paris, France (Lopez et al., 2013). Like Evanston, natural gas was the major fossil fuel contributor in the winter in Salt Lake City (Pataki et al., 2006b, 2007) and Paris (Lopez et al., 2013). In Wroclaw, carbon isotope data revealed that coal was the major fossil fuel contributor in the winter (Górka and Lewicka-Szczebak, 2013). Evanston respiration values are lower than the 52%–66% estimated for Los Angeles (Djuricin et al., 2010), which is expected given that warmer temperatures and seasonally higher rainfall in Los Angeles support plant growth in the winter. In Los Angeles and Portland, the winter contribution from petroleum combustion was 25%–50% higher than the contribution from natural gas combustion (Djuricin et al., 2010; Rice and Bostrom, 2011). In Los Angeles, winter natural gas usage is low by comparison to Evanston, and the Portland study did not consider ecosystem respiration. Clearly, climate, extent of vegetation cover, and local energy infrastructure affect the relative contributions of ecosystem respiration and fossil fuels to urban carbon cycles, with seasonal variations expected for cities located in temperate climates.

## 4 Conclusions

This study employed a Picarro G1101-i WS-CRDS to generate nearly continuous records for the concentration and carbon isotope composition (δ13C) of atmospheric CO2 in Evanston, IL for the time period August 2011 through February 2012. Evanston is adjacent to Chicago, IL, which is the third largest city in the US. Evanston displayed high [CO2] and low δ13C values compared to background sites in Park Falls, WI, and Mauna Loa, HI. Similar to previous studies conducted in other urban areas (Clark-Thorne and Yapp, 2003; Pataki et al., 2006b, 2007; Lietzke and Vogt, 2013), the patterns for Evanston reflect the combined influence of natural and anthropogenic factors, especially ecosystem respiration and local fossil fuel CO2 emissions.

While overall seasonal variations in [CO2] and δ13C largely followed broader midcontinental trends, relatively large amplitude shifts in local δ13C values point to seasonally changing CO2 sources in Evanston. Similarities in the timing and size of average diurnal fluctuations in [CO2] for Evanston and Park Falls during the summer suggest that ecosystem respiration influenced [CO2] in Evanston. Larger average diurnal [CO2] fluctuations in Evanston compared to Park Falls during fall and winter point to the importance of local fossil fuel CO2 emissions. In particular, relatively high [CO2] and low δ13C values measured during the fall and winter suggest elevated contributions from natural gas combustion.

On any given day, a combination of local natural CO2 cycling, anthropogenic activity, and atmospheric conditions controlled [CO2] and δ13C values. For hourly and daily periods, wind speed and direction strongly influenced atmospheric [CO2] and δ13C values. The highest [CO2] and lowest δ13C values were observed at wind speeds <3 m s−1 and when winds blew from the northwest, west, and south over densely populated areas. Also, a large local CO2 source, the Northwestern University central steam plant, was detectable when wind blew from the northeast, forcing emissions from the plant toward the measurement site.

Keeling plots were used to constrain seasonally changing CO2 sources. In summer, nighttime measurements suggest that local ecosystem respiration contributed 50%–60% of the CO2 in excess of background and that fossil fuel combustion, primarily petroleum combustion, contributed the remainder. In winter, ecosystem respiration only contributed 0%–10% of the excess CO2 and fossil fuel emissions, primarily natural gas combustion, contributed the remaining 90%–100%. This study confirms the utility of δ13C for estimating ecosystem and fossil fuel contributions to urban CO2. Overall, the seasonal pattern in Evanston/Chicago is similar to that observed for other northern hemisphere, mid-latitude cities where local vegetation provides a CO2 sink and source during summer and fossil fuel emissions, from natural gas in the case of Evanston, dominate CO2 contributions during winter.

## Supplemental material

Table S1.
CO2 concentration and δ13C value calibration and check standard data
File Type: PDF
File Size: 0.03 MB

Figure S1.
Time series of CO2 concentrations and δ13C values for check standards
File Type: TIF
File Size: 1.72 MB

Figure S2.
Frequency plot of differences in wind direction and speed between weather stations
File Type: TIF
File Size: 0.39 MB

Figure S3.
Seasonal polar plots with mean CO2 concentrations as a function of wind direction and speed
File Type: TIF
File Size: 3.52 MB

## Data accessibility statement

The NOAA CO2 concentration and δ13C values used in this study are publicly available from the NOAA or Ameriflux websites. The climate data are publically available from the National Weather Service and from the Illinois Environmental Protection Agency. The CO2 concentrations and δ13C values generated in this study along with the compiled meteorological data are available from the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.h8d0n.

© 2015 Moore and Jacobson. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

## Contributions

Contributed to conception and design: JM, ADJ

Contributed to acquisition of data: JM

Contributed to analysis and interpretation of data: JM, ADJ

Drafted and/or revised the article: JM, ADJ

Approved the submitted version for publication: JM, ADJ

## Competing interests

The authors do not have any competing interests.

## Funding information

The project was mostly supported by a David and Lucile Packard Foundation Fellowship awarded to ADJ and the Institute for Sustainability and Energy at Northwestern (ISEN). We are grateful for additional funding provided by several sources at Northwestern University, including the Weinberg College of Arts and Sciences, the McCormick School of Engineering, and the Vice President’s Office for Research. JM thanks the Jess and Mildred Fisher Endowed Chair and the Faculty Development and Research Committee at Towson University for financial support.