# Under-reporting of greenhouse gas emissions in U.S. cities

Kevin Robert Gurney, Jianming Liang, Geoffrey Roest, Yang Song, Kimberly Mueller, Thomas Lauvaux

<figure><img src="image_1.png" /></figure>

## To cite this version:

Kevin Robert Gurney, Jianming Liang, Geoffrey Roest, Yang Song, Kimberly Mueller, et al. Under-reporting of greenhouse gas emissions in U.S. cities. Nature Communications, 2021, 12, 10.1038/s41467-020-20871-0. insu-03721920

HAL Id: insu-03721920

https://insu.hal.science/insu-03721920

Submitted on 13 Jul 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-020-20871-0

OPEN

<figure><img src="image_4.png" /></figure>

Check for updates

# Under-reporting of greenhouse gas emissions in U.S. cities

Kevin Robert Gurney ID <fnref n="1" />✉, Jianming Liang², Geoffrey Roest<fnref n="1" />, Yang Song<fnref n="1" />, Kimberly Mueller ID ³ & Thomas Lauvaux⁴

Cities dominate greenhouse gas emissions. Many have generated self-reported emission inventories, but their value to emissions mitigation depends on their accuracy, which remains untested. Here, we compare self-reported inventories from 48 US cities to independent estimates from the Vulcan carbon dioxide emissions data product, which is consistent with atmospheric measurements. We found that cities under-report their own greenhouse gas emissions, on average, by 18.3% (range: -145.5% to +63.5%) - a difference which if extrapolated to all U.S. cities, exceeds California's total emissions by 23.5%. Differences arise because city inventories omit particular fuels and source types and estimate transportation emissions differently. These results raise concerns about self-reported inventories in planning or assessing emissions, and warrant consideration of the new urban greenhouse gas information system recently developed by the scientific community.

<figure><img src="image_5.png" /></figure>

<footnote n="1">¹School of Informatics, Computing, and Cyber Systems, Northern Arizona University, Flagstaff, AZ, USA. ²ESRI, Redlands, CA, USA. ³National Institute for Standards and Technology, Gaithersburg, MD, USA. ⁴Laboratoire des Sciences du Climat et de l'Environnement, Gif-sur-Yvette, France.</footnote>

<footnote>✉email: kevin.gurney@nau.edu</footnote>
Almost three-quarters of fossil fuel carbon dioxide (FFCO₂) emissions, the most important anthropogenic greenhouse gas (GHG), emanate from cities¹. Projections show cities could add over two billion people this century with global urban area tripling by 2030²,³. This outsized role is compounded by the concern that urban development can lock-in resistance to low-carbon pathways through technological, institutional, and behavioral inertia⁴. At the same time, cities are taking a leadership role in climate emissions mitigation activities with cities across the globe pledging to ambitious emissions mitigation goals and playing an increasingly important role in the international climate change negotiating process⁵. A necessary element in urban emissions mitigation is the development of a numerical accounting of GHG emissions following a standardized method⁶. This enables emissions mitigation planning activities and the establishment of emissions baselines with which to establish reduction targets. Cities across the globe have developed estimates of their GHG emissions, which we refer to as “self-reported inventories” (SRIs), mostly following one of a few publicly available protocols⁷⁻⁹. However, there is no systematic, peer-reviewed assessment of SRI quality or accuracy in spite of their importance to establishing baseline emissions, urban mitigation target setting, and mitigation outcome assessment.

Given this assessment gap, we have investigated 48 SRIs published online by cities in the USA and compare these to urban extractions from our recently completed, research-driven Vulcan version 3.0 estimate of US fossil fuel carbon dioxide (FFCO₂) emissions. Vulcan quantifies complete FFCO₂ emissions at points (e.g., factories, powerplants), lines (i.e., roadways), and polygons (i.e., US census block-groups) for the whole US spanning the 2010–2015 time frame using a large number of publicly available input datasets¹⁰. The data inputs are primarily archived from local fuel, activity, and flux monitoring data. The Vulcan v3.0 results have been compared to an emissions estimate based on observations of ¹⁴CO₂ in the USA, an ideal atmospheric tracer for FFCO₂, and is within 1.5% of the atmospheric-based result¹¹. Furthermore, the urban-intensive INFLUX effort in the city of Indianapolis shows Vulcan results within 3% of atmospheric-based results over a 3-year mean timespan, and within 6% annually¹².

The samples of 48 cities included in this comparison were chosen through a combination of previously published work¹³, the presence among the top 100 emitting cities in the USA, and the availability of adequately documented SRI results. The 48 sample cities account for 13.7% of 2015 US urban emissions and 17.7% of 2015 US urban population. Unambiguous city government boundaries were used to extract sector/fuel-specific FFCO₂ emissions from the “native” (points, lines, polygons) Vulcan FFCO₂ emissions. Both the Vulcan and SRI results were inspected to ensure geographic alignment of city boundaries and emission sectors, gases, and scope. For example, were individual emission sectors (e.g. airport, industrial) absent in the SRIs, they were similarly eliminated in Vulcan, to the extent SRI documentation offers clarity on sector representation. The omission of sub-sector elements in the SRI (e.g., individual fuels within a sector, individual factories in the city domain), rarely documented, are considered part of the comparison in this analysis, and hence not adjusted for (see “Methods”).

## Results

The mean relative difference (RD) between the two emission datasets is +18.3% (Vulcan > SRI; calculated as [(Vulcan - SRI)/average(Vulcan, SRI))] with a mean absolute (unsigned) relative difference (MAD) of 29.1% (Fig. 1). This is in the context of a mean RD 95% confidence interval of +5.2%/+31.7% for the Vulcan urban emissions. Of the cities reporting fewer emissions than Vulcan (N = 37), the mean RD value is +30.7%. Furthermore, the summed difference across all 48 cities in units of emitted mass per year is 19,076,760 tC/year, a value nearly equivalent to the 2015 Massachusetts state emissions (Table 1). Were the 18.3% difference extrapolated to all US cities, the total would be 129,219,255 tC/year (Source data), an amount 23.5% larger than the entire 2015 California state emissions. We find a smaller RD value in the Western half of the USA (+11.0%) compared to the East (+25.0%), but a nearly identical MAD value (30.3% and 28.0%, respectively) due to the greater number of negative differences in the Western cities, causing some cancellation of the positive differences. Finally, we find little correlation (Pearson's R = 0.01) between the individual absolute relative difference (AD) values and the magnitude of city emissions.

<figure><img src="image_6.png" /><figcaption>**Fig. 1** Total individual city (N = 48) FFCO₂ emissions and absolute difference (AD = positive RD) between the Vulcan version 3.0 data product and self-reported inventories (SRIs). AD values denoted by color scale and total emissions denoted by circle size. Two regions are expanded (California Bay Area and the Northeast coastal urban corridor) for closer inspection (see “Methods” and Source data for additional details). Background map tiles by Stamen Design (no changes made), under CC BY 3.0, data by © OpenStreetMap contributors, under ODbL.</figcaption></figure>
<table><caption>Table 1 State total FFCO₂⁺ cement production emissions for the year 2015 as generated by Vulcan version 3.0¹⁰.</caption><thead><tr><td>State</td><td>Total emissions (tC)</td></tr></thead><tbody><tr><td>AK</td><td>7698510</td></tr><tr><td>AL</td><td>42704041</td></tr><tr><td>AR</td><td>16561519</td></tr><tr><td>AZ</td><td>29782799</td></tr><tr><td>CA</td><td>104660731</td></tr><tr><td>CO</td><td>31918346</td></tr><tr><td>CT</td><td>10610878</td></tr><tr><td>DC</td><td>1230078</td></tr><tr><td>DE</td><td>4445990</td></tr><tr><td>FL</td><td>78340164</td></tr><tr><td>GA</td><td>42264025</td></tr><tr><td>HI</td><td>4956370</td></tr><tr><td>IA</td><td>21843909</td></tr><tr><td>ID</td><td>6694903</td></tr><tr><td>IL</td><td>67373603</td></tr><tr><td>IN</td><td>55895053</td></tr><tr><td>KS</td><td>19782874</td></tr><tr><td>KY</td><td>35159975</td></tr><tr><td>LA</td><td>47294753</td></tr><tr><td>MA</td><td>20694772</td></tr><tr><td>MD</td><td>19458650</td></tr><tr><td>ME</td><td>5359411</td></tr><tr><td>MI</td><td>49813475</td></tr><tr><td>MN</td><td>33056713</td></tr><tr><td>MO</td><td>37477007</td></tr><tr><td>MS</td><td>18315636</td></tr><tr><td>MT</td><td>8950523</td></tr><tr><td>NC</td><td>42308189</td></tr><tr><td>ND</td><td>13208952</td></tr><tr><td>NE</td><td>14012960</td></tr><tr><td>NH</td><td>4360434</td></tr><tr><td>NJ</td><td>26118536</td></tr><tr><td>NM</td><td>14628211</td></tr><tr><td>NV</td><td>12791132</td></tr><tr><td>NY</td><td>53255701</td></tr><tr><td>OH</td><td>68698200</td></tr><tr><td>OK</td><td>30469338</td></tr><tr><td>OR</td><td>11373707</td></tr><tr><td>PA</td><td>70973620</td></tr><tr><td>RI</td><td>5099102</td></tr><tr><td>SC</td><td>22082573</td></tr><tr><td>SD</td><td>3604415</td></tr><tr><td>TN</td><td>28900311</td></tr><tr><td>TX</td><td>151055606</td></tr><tr><td>UT</td><td>20803673</td></tr><tr><td>VA</td><td>28278324</td></tr><tr><td>VT</td><td>1839216</td></tr><tr><td>WA</td><td>20867652</td></tr><tr><td>WI</td><td>31991556</td></tr><tr><td>WV</td><td>23249400</td></tr><tr><td>WY</td><td>21966262</td></tr><tr><td>US</td><td>1544281777</td></tr></tbody></table>

Figure 2 provides the individual city RD values for the total emissions and the mean across-sector AD (MASAD) in each city. The mean of these MASAD values for all cities is 50.3%. This demonstrates that even for cities with close agreement in total emissions, individual sectors show large discrepancies. For example, the city of Detroit has a total RD of +1.3%, but a MASAD value of 44.2% comprised of a RD of -28.5% for onroad sources and a compensating RD of +38.7% for stationary sources (Fig. 3). The mean of the sector-specific AD values across all cities for the onroad (N = 46), stationary (N = 43), and other transportation (N = 27) sectors is 28.1%, 37.9%, and 82.6%,

respectively (Fig. 3). However, for those cities reporting fewer emissions than Vulcan, the differences are 25.7% (N = 28), 44.9% (N = 23), and 102.2% (N = 11), respectively.

Careful examination of the differences between Vulcan v3.0 and the SRIs suggests that the most common differences are related to SRI omission of petroleum fuel use and point source emissions in the industrial/commercial sector, different accounting perspectives on marine shipping and airborne emissions, and different methods in onroad emission estimation. Many cities collect utility natural gas consumption data on sector-based stationary sources. These data will oftentimes miss petroleum consumption occurring in buildings typically associated with industrial point source processes or smaller fuel supply contracts that are not provided by a centralized utility. SRIs will often isolate the component of marine shipping to local activity in harbor areas or attempt to estimate the resident population utilization of airborne activity. Cities also estimated onroad emissions variously using gas station fuel sales or travel demand modeling to isolate onroad emissions stimulated by activity within the city. Vulcan, by contrast, includes all emission sources within the administrative city boundary, including airborne emissions from taxi/takeoff/landing up to 3000 ft above ground, all marine shipping emissions in waters within the city boundary, and emissions from all roadways within the city regardless of the origin/destination of travel. Further explanation of the individual city/sector differences are difficult to isolate with the SRI information provided. An important caveat to this study is that better documentation of SRI system boundary details could change the numerical difference statistics found here.

Discussion

The development of SRIs was a rational response to the need by cities to tackle GHG mitigation, a necessary and globally relevant effort, given the large role cities play in emissions and their governance purview over many aspects of the urban emission landscape. However, constructing baseline emissions and objective tracking over time using city SRIs raises several challenges, and these are likely responsible for the differences found in this study. The development of an SRI is a costly endeavor, placing a burden on city staff and resources. Data collection, processing, and modeling can present technical challenges for cities, further burdening available resources and often resulting in incomplete estimates or biased outcomes, as demonstrated here. Compounding the problem, independent objective assessment of SRI estimates creates additional technical hurdles and is thus, rarely attempted. These challenges are particularly important when placed in the context of the reduction targets. For example, the city of Indianapolis has indicated that they aim to make a 20% reduction in building GHG emissions between 2025 relative to 2016 values¹⁴. However, with the 26.9% underestimate found here, it will be difficult to know when and if this target is truly achieved or track progress towards it.

A potential solution to this challenge is a more systematic GHG flux information system that combines the type of bottom-up quantification exemplified by the Vulcan data product used here with a tiered (ground, airborne, space-borne) atmospheric observation and model system¹⁵. Rather than being independent of individual city efforts, building and managing a GHG flux information system will require collaboration with city staff, availing of local knowledge and tailored to local needs. An urban-scale prototype of such a system has been piloted in the city of Indianapolis through the INFLUX experiment with similar prototypes being developed at the national scale¹¹,¹². This approach could offer a comprehensive, consistent approach across the entire continental landscape delivering fine spatial (e.g., down to
<figure><img src="image_7.png" /><figcaption>**Fig. 2** Sector-based individual city (N = 48) FFCO₂ emissions relative difference (RD) between the Vulcan version 3.0 urban FFCO₂ emissions and self-reported inventories (SRIs). Black: total emissions RD; error bars: total emissions RD range from Vulcan 95% confidence interval; red: mean across-sector absolute relative difference (MASAD). RD calculated as [(Vulcan - SRI)/mean(Vulcan, SRI)]. Scale capped at 100%. Not all sector categories available in all cities; where individual sectors were missing in the SRIs, the mean only includes those available in both Vulcan and SRI results (see "Methods" and Source data for additional details).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 3** Sector-based individual city (N = 48) FFCO₂ emissions relative difference (RD) between the Vulcan version 3.0 urban FFCO₂ emissions and self-reported inventories (SRIs). Black: total emissions RD; error bars: total emissions RD range from Vulcan confidence interval; red: onroad emissions RD; green: stationary (residential + commercial + industrial) emissions RD; blue: other transportation (airborne + marine shipping + nonroad + railroad) emissions RD. RD calculated as [(Vulcan - SRI)/average (Vulcan, SRI)]. Scale capped at 100% RD. Not all sector categories were available in all cities; where individual sectors were missing in the SRIs, the aggregate category only sums those available in both Vulcan and SRI results (see Source data for additional details).</figcaption></figure>

the neighborhood/street scale) and functional (e.g., sector, sub-sector, combustion process) detail for better targeting of mitigation policy¹⁶. It also offers the ability to track mitigation progress using the most accurate arbiter of impact, atmospheric GHG amounts. Finally, an advantage of these approaches is their ability to conserve mass through spatial scales, allowing for consistency, comparability, and context from the nation to the building.

Accuracy and precision are critical to estimating GHG emissions, whether reported by a city, state, or country. The absence of an accurate emissions assessment (i.e., baseline and ongoing) makes prioritizing mitigation policy options difficult, can lead to misallocation of scarce mitigation resources, and presents challenges to independent assessment and course correction. The results presented here raise serious concerns about the current self-reported approach to quantifying urban GHG emissions in the USA. Similar dynamics may be at play in cities across the globe where SRIs have similarly been reported¹⁷. Fortunately, there is progress on building a systematic emissions quantification system that promises a systematic approach to generating space/time-resolved, atmospherically calibrated emissions information for all cities in collaboration with local authorities. With such a collaborative system, urban GHG mitigation practitioners

can devote time and resource to the activity they have the greatest knowledge and political influence over: the best mitigation strategies for their city.

## Methods

**Cities analyzed.** The cities compared in this study were collected according to a combination of available literature, emissions magnitude, and adequacy of documentation. A recent paper by Nangini et al. provides a list of cities that have completed SRIs, reporting them to the Carbon Disclosure Project (CDP)¹⁸. This list of 64 US cities was combined with the top 50 urban emitters from the 2011 Vulcan estimate. A systematic online search was performed to locate documents that reported SRI results for each of these cities. Many cities did not have published inventories or the inventories were published but the information content was insufficient to isolate those emissions comparable to the Vulcan results. Furthermore, SRI estimates were not considered if they reported for years outside the Vulcan timespan (2010–2015). Four exceptions were made for cities with emissions reported in years very close to the Vulcan timespan: Columbus, OH: 2017; New Orleans, LA: 2007; Longmont, CO: 2016; Salt Lake City, UT: 2009. In these cases, the year closest to the Vulcan timespan (2010 or 2014) was used. The sensitivity of the results to the removal of these four cities was negligible with little impact on the mean statistics (Table 2).
<table><caption>**Table 2 Standard case and sensitivity of statistical results to removal of four cities with SRIs outside the 2010–2015 timespan of the Vulcan output (Columbus, OH: 2017; New Orleans: 2007; Longmont, CO: 2016; Salt Lake City, UT: 2009) and to the division of the Eastern US versus the Western US.**</caption><thead><tr><th>Sensitivity cases</th><th>No. of cities</th><th>Mean RD (%) [stdev]</th><th>MAD (%)</th><th>Mean RD of positive diff (%) [N]</th><th>Mean RD with 95% CI</th></tr></thead><tbody><tr><td>Standard</td><td>48</td><td>+18.3 [38.1]</td><td>+29.1</td><td>+30.7 [37]</td><td>+5.2/+31.7</td></tr><tr><td>Remove four cities</td><td>44</td><td>+18.5 [39.2]</td><td>+29.2</td><td>+30.8 [34]</td><td>+5.4/+31.8</td></tr><tr><td>Eastern US</td><td>25</td><td>+25.0 [31.2]</td><td>+28.0</td><td>+30.1 [22]</td><td>+12.8/+37.8</td></tr><tr><td>Western US</td><td>23</td><td>+11.0 [44.0]</td><td>+30.3</td><td>+31.6 [15]</td><td>-3.0/+25.1</td></tr></tbody></table>

The division of East/West was placed at a longitude of −94.0°. See Source data for additional detail.

<figure><img src="image_9.png" /><figcaption>**Fig. 4 Rank distribution (red) and cumulative emissions (blue) of US city FFCO₂ emissions.** The sample of 48 cities included in this study are denoted with black “X” within the rank distribution. City FFCO₂ emissions are from the Vulcan v3.0 results within US Census Designated Places. The sample of 48 cities used for comparison represents 13.7% of 2015 US urban FFCO₂ emissions and 17.7% of the 2015 US urban population. Of these 48, 28 cities are within the top 100 urban emitters in the USA (see Source data for additional details).</figcaption></figure>

The emissions ranking of the final list of 48 cities in relation to a rank distribution of all cities in the USA, in addition to their share of cumulative urban emissions, is shown in Fig. 4.

**Urban boundaries.** To align the urban boundaries in the annual Vulcan version 3.0 results to the city boundaries in the SRIs, we use the US Census Bureau Designated Places (DPs) reflecting boundary conditions in the year 2010¹⁹. Online SRI documentation rarely defines the precise boundary of the domain considered but it is assumed that the SRI emissions reflect those of the stated city within the governance boundary, unless otherwise described. The US Census DP for the 48 cities considered in this study map exactly to a Google map boundary search for the named cities. The only exception relates to coastal cities in which the DP boundary will extend offshore, while cities may only include emissions from the land area within their administrative boundary. This will create differences in the marine shipping sector.

The Vulcan emissions were regridded at 300 m × 300 m so they can be more finely attributed to the DP entities. The parameters of the 300 m grid are shown in Table 3. The DP shapefile, containing 28,705 polygon features, each representing a DP entity, was rasterized using ArcGIS into the same 300 m × 300 m Vulcan grid. In rasterizing the DP polygons, each grid cell was assigned a DP identity from the polygon that overlaps the center of the cell (https://pro.arcgis.com/en/pro-app/tool-reference/conversion/polygon-to-raster.htm). The 300 m × 300 m gridded Vulcan output was then aggregated by DP entity to generate the emissions values for the 28,705 DP entities. The caveat of this method is that certain small, localized

<table><caption>**Table 3 Parameter values for the 300 m grid used to intersect Vulcan emissions with urban boundaries.**</caption><thead><tr><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>Datum/ellipsoid</td><td>WGS84</td></tr><tr><td>Projection</td><td>Lambert Conformal Conic</td></tr><tr><td>Dimensions</td><td>9551 rows × 15,340 columns</td></tr><tr><td>Resolution</td><td>300 m</td></tr><tr><td>Extent</td><td>xmin = −2,272,425, xmax = 2,329,575, ymin = −1,589,775, ymax = 1,275,525</td></tr><tr><td>Longitude of origin</td><td>97.0°W</td></tr><tr><td>Latitude of origin</td><td>40.0°N</td></tr><tr><td>Reference latitude 1</td><td>33.0°N</td></tr><tr><td>Reference latitude 2</td><td>45.0°N</td></tr></tbody></table>

Both the US Census Bureau DPs and the shapefiles underlying the Vulcan v3.0 emissions were gridded into a 300 m grid. Vulcan emissions were then aggregated by DP entity to generate emissions for each designated place.

sources (e.g., point sources) may fall within a grid cell associated with a certain DP entity, even if the actual point source falls just outside of the urban boundary polygon. Likewise, it is possible that localized sources located just within the urban boundary polygon may be omitted if the DP entity is not assigned to a grid cell that only partly overlaps the polygon. This is minimized by the fact that many city boundaries are located in less dense portions of an urban area.

A population value was also assigned to each DP entity following the same method with the US Census Bureau 2010 population data¹⁹.

**Dataset alignment.** Examination of the 48 SRI results made clear that a direct comparison of the Vulcan output to the SRI output would suffer from a series of categorical mismatches. Hence, we proceeded with a number of adjustments to align the results for comparison. The first of these was to ensure consistency in scope. Three scopes have been identified in the literature: scopes 1, 2, and 3²⁰,²¹. When applied to FFCO₂ emissions, scope 1 reflects emissions that physically emanate from fuel combustion within the geography considered (i.e., city administrative boundary). Scope 2 refers to FFCO₂ emissions resulting from the consumption of electricity occurring within the boundary of the city, regardless of where the emissions physically occur. For example, a given city may have one power production facility within city limits. City residents may consume electricity at such a level as to use all electricity at that in-city facility and electricity from an additional power production facility outside city limits. Scope 1 emissions in the electricity production sector, in this example, would be comprised of that single facility. Scope 2 emissions, by contrast, would include the emissions from both facilities in proportion to the electricity consumed within the city limits. Finally, scope 3 emissions allocate the emissions resulting from the production of goods and services, regardless of where they occur, to the point of consumption. For example, were a factory in our example city producing goods that are exported to consumers outside city limits, those emissions would not be accounted for within the example city, but allocated to the locations where final consumption occurred. Conversely, the quantity of emissions necessary to produce and transport goods to the example city consumers would be counted as within-city emissions.

Vulcan quantifies all CO₂ emissions from the combustion of fossil fuels and cement production within the geographic boundary of the city domain or scope 1 emissions. Vulcan does not include biogenically originated emissions, fugitive emissions, or non-CO₂ GHGs. Because none of the SRIs examined for this study included emissions from cement production (stoichiometric emissions as opposed to fossil fuel combustion used to process cement production), these emissions were removed from the Vulcan results.
The reason that the Vulcan FFCO₂ emissions data product was initially restricted to scope 1 emissions was the use of Vulcan within systems that utilize atmospheric mixing ratios measurements to infer surface fluxes, often referred to as the atmospheric “inversion” approach²²⁻²⁴. Vulcan and other space/time-explicit emissions data products are often used as a boundary condition or “prior” flux to these systems. The inversion systems use models of atmospheric motion to invert atmospheric mixing ratios over a given spatial domain, and hence the surface flux boundary condition or prior flux must reflect fluxes specific to the geography considered—scope 1 emissions only.

Scope 2 and 3 are of understandable interest to urban GHG practitioners. However, to avail of the independent estimate of fluxes that atmospheric inversion systems can provide, scope 1 emissions (which typically comprises the largest emissions among the scope perspectives) are important to isolate within the overall emissions accounting. Independent assessment using atmospheric inversions offers an important element in evaluating inventory estimates and course-correcting mitigation activity.

The city SRIs, in contrast to the Vulcan approach, often will include elements of scopes 1, 2, and 3 combined to comprise the total city inventory. The mixture of scopes is often a reflection of individual city interests and partly a reflection of the protocols currently used by cities to estimate emissions. For example, Vulcan quantifies all FFCO₂ emissions from powerplants residing within the city boundary regardless of where the electricity is sold (scope 1). Most cities, however, will estimate the FFCO₂ emissions driven by consumption of electricity within the boundary of the city, regardless of where those emissions physically occur (scope 2).

In light of these differences, each city SRI was carefully examined in order to isolate the scope 1 sector emissions. Where scope 1 emission sectors were absent in the SRIs, they were similarly eliminated in Vulcan. Typical sectors missing in the SRIs were airport, railroad, and nonroad. Vulcan currently does not report emissions associated with the decay of material in the waste stream, and hence this was removed from the SRIs. Removal of the emissions associated with waste also assists in more closely aligning the emitted gases. Most of the city SRIs report in units of equivalent CO₂, which incorporates non-CO₂ GHGs weighted by their relative influence on climate change. Since the non-CO₂ emissions in the SRIs are dominated by waste processing and Vulcan only quantifies CO₂ emissions from fossil fuel combustion. Nevertheless, there may be a small amount of non-CO₂ emissions within the SRIs that are not quantified in Vulcan and we would, therefore, expect the SRI results to report larger emissions than Vulcan, all else being equal.

The average proportion of the individual city emissions after these alignment steps relative to the original total was 86.4%, indicating that the comparison carried out here captures most of the original SRI emissions. All of the individual adjustments to both the Vulcan and SRI results are described in Source data.

**Comment on Nangini et al.'s study.** Nangini et al. reported on an archive of SRIs that, in the USA, were submitted to the CDP¹³. These SRIs are supplied to the CDP directly with little documentation or information support. Nangini et al. performed some quality control on the submitted data and appended the SRI submissions with a variety of socioeconomic ancillary data (e.g., population, gross domestic product, climate). The paper was primarily aimed at documentation support for the contained data archive that included many cities beyond those in the USA. The USA submissions to the CDP, either as originally deposited or as represented in the Nangini et al. archive, were not used in this study. The individual city estimates were not documented sufficiently making alignment of system boundaries difficult to achieve. Indeed, for those overlapping cities between the current study and those archived in Nangini et al., numerical consistency was difficult to achieve. Hence, this study chose to acquire original SRI estimates with available documentation as the SRI dataset upon which to build the comparison reported here.

A limited urban emissions comparison to an older version of the Vulcan data product (version 2.0) was included in Nangini et al. study. However, Vulcan version 2.0 is resolved at a spatial resolution of 100 km² (compared to 1 km² in version 3.0) and represents emissions in the calendar year 2002. The isolation of urban areas in the Nangini et al. study used circular outlines approximating city geographies. Since none of the SRI estimates used in Nangini et al. were aligned with the year 2002, the combination of calendar year mismatch, low-resolution Vulcan results, lack of sectoral detail, differing system boundaries, and the approximate geographical representation make the comparison performed in that paper inadequate as a comparison analysis.

**Data availability**

The complete Vulcan version 3.0 can be downloaded from the Oak Ridge National Laboratory Distributed Active Archive Center (https://doi.org/10.3334/ORNLDAAC/1741)²⁵. Source data are provided with this paper. All other data and materials are available upon request.

Received: 27 April 2020; Accepted: 15 December 2020;
Published online: 02 February 2021

**References**

1. Intergovernmental Panel on Climate Change. in *Climate Change 2014: Mitigation of Climate Change* (Cambridge Univ. Press, 2015).

2. United Nations Department of Economic and Social Affairs. *World Urbanization Prospects: The 2018 Revision. World Urbanization Prospects: The 2018 Revision* (United Nations Department of Economic and Social Affairs, 2019).

3. Seto, K. C., Guneralp, B. & Hutyra, L. R. Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools. *Proc. Natl Acad. Sci. USA* **109**, 16083–16088 (2012).

4. Seto, K. C. et al. Carbon Lock-in: Types, Causes, And Policy Implications. *Annu. Rev. Environ. Resour.* **41**, 425–452 (2016).

5. ARUP, C40 & University of Leeds. *The Future of Urban Consumption World in a 1.5°C World C40—Headline Report* (ARUP, C40 & University of Leeds, 2019).

6. Hsu, A. et al. A research roadmap for quantifying non-state and subnational climate mitigation action. *Nat. Clim. Chang.* **9**, 11–17 (2019).

7. ICLEI—Local Governments for Sustainability USA. U.S. Community Protocol for Accounting and Reporting of Greenhouse Gas Emissions, Version 1.2. 1–67. www.icleiusa.org/tools/ghg-protocol/community-protocol/us-community-protocol-for-accounting-and-reporting-of-greenhouse-gas-emissions (2019).

8. Schultz, S. et al. *Global Protocol for Community-Scale Greenhouse Gas Emissions*, V2.0. (2014).

9. Global Covenant of Mayors for Climate Energy. *Global Covenant of Mayors Common Reporting Framework Table of Contents* (Global Covenant of Mayors for Climate Energy, 2018).

10. Gurney, K. R. et al. The Vulcan version 3.0 high-resolution fossil fuel CO₂ emissions for the United States. *J. Geophys. Res. Atmos.* **125**, e2020JD032974 (2020).

11. Basu, S. et al. Estimating US fossil fuel CO₂ emissions from measurements of 14C in atmospheric CO₂. *Proc. Natl Acad. Sci. USA* **117**, 13300–13307 (2020).

12. Lauvaux, T. et al. Policy-relevant assessment of urban CO₂ emissions. *Environ. Sci. Technol.* **54**, 10237–10245 (2020).

13. Nangini, C. et al. A global dataset of CO₂ emissions and ancillary data related to emissions for 343 cities. *Sci. Data* **6**, 1–29 (2019).

14. City of Indianapolis Office of Sustainability. *Thrive Indianapolis* (City of Indianapolis Office of Sustainability, 2019).

15. Mueller, K., Robert Gurney, K. & Lauvaux, T. The policy value of urban greenhouse gas information system. *Earth's Fut.* (in preparation).

16. Gurney, K. R. Track urban emissions on a human scale. *Nature* **525**, 179–180 (2015).

17. Global Covenant of Mayors. *Global Covenant of Mayors for Climate and Energy, Data4Cities* (accessed 6 March 2020); www.globalcovenantofmayors.org/our-initiatives/data4cities/.

18. Carbon Disclosure Project. *Carbon Disclosure Project* (accessed 6 March 2020); https://data.cdp.net/browse?category=Emissions.

19. United States Census Bureau. *2010 Census Summary File 1, 2010 Census of Population and Housing. Technical Documentation* (United States Census Bureau, 2012).

20. Chen, G. et al. Review on city-level carbon accounting. *Environ. Sci. Technol.* **53**, 5545–5558 (2019).

21. Lombardi, M., Laiola, E., Tricase, C. & Rana, R. Assessing the urban carbon footprint: an overview. *Environ. Impact Assess. Rev.* **66**, 43–52 (2017).

22. Lauvaux, T. et al. High-resolution atmospheric inversion of urban CO₂ emissions during the dormant season of the Indianapolis flux experiment (INFLUX). *J. Geophys. Res.* **121**, 5213–5236 (2016).

23. Brioude, J. et al. Top-down estimate of surface flux in the Los Angeles Basin using a mesoscale inverse modeling technique: assessing anthropogenic emissions of CO, NOx and CO₂ and their impacts. *Atmos. Chem. Phys.* **13**, 3661–3677 (2013).

24. Vogel, F. R. et al. Regional non-CO₂ greenhouse gas fluxes inferred from atmospheric measurements in Ontario, Canada. *J. Integr. Environ. Sci.* **9**, 41–55 (2012).

25. Gurney, K. R. et al. Vulcan fossil fuel carbon dioxide (FFCO₂) emissions data product, version 3.0, 1 km grid https://doi.org/10.3334/ORNLDAAC/1741 (2019).

**Acknowledgements**

We thank Bruce Hungate, Kate Petersen, and Phil DeCola for editorial input. K.R. Gurney, J. Liang, G. Roest, and Y. Song received support from the National Aeronautics and Space Administration Grant NNX14A20G and the National Institute of Standards and Technology Grant 70NANB16H264N. T. Lauvaux was supported by the French research program Make Our Planet Great Again (project CIUDAD).
**Author contributions**

K.R.G. conceptualized the research, performed analysis and calculations, and wrote the manuscript. J.L. performed analysis and calculations. G.R. and Y.S. contributed to the SRI analysis and manuscript development. K.M. and T.L. contributed to the manuscript development.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41467-020-20871-0.

**Correspondence** and requests for materials should be addressed to K.R.G.

**Peer review information** *Nature Communications* thanks the anonymous reviewer(s) for their contribution to the peer review of this work.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021