Statistical characterization of urban CO2 emission signals observed by commercial airliner measurements.

Cities are responsible for the largest anthropogenic CO2 emissions and are key to effective emission reduction strategies. Urban CO2 emissions estimated from vertical atmospheric measurements can contribute to an independent quantification of the reporting of national emissions and will thus have political implications. We analyzed vertical atmospheric CO2 mole fraction data obtained onboard commercial aircraft in proximity to 36 airports worldwide, as part of the Comprehensive Observation Network for Trace gases by Airliners (CONTRAIL) program. At many airports, we observed significant flight-to-flight variations of CO2 enhancements downwind of neighboring cities, providing advective fingerprints of city CO2 emissions. Observed CO2 variability increased with decreasing altitude, the magnitude of which varied from city to city. We found that the magnitude of CO2 variability near the ground (~1 km altitude) at an airport was correlated with the intensity of CO2 emissions from a nearby city. Our study has demonstrated the usefulness of commercial aircraft data for city-scale anthropogenic CO2 emission studies.

Climate change is considered to be one of the consequences of increased emissions of anthropogenic greenhouse gases during the industrial era, and carbon dioxide (CO 2 ) is the dominant contributor to the enhanced radiative forcing caused by anthropogenic long-lived greenhouse gases. Atmospheric CO 2 mole fraction has increased from ~280 ppm (parts per million) in 1750 1 to > 400 ppm in recent years 2 due to the rapid growth of human activities and population since the beginning of the industrial era. About half of the anthropogenic CO 2 emissions related to fossil fuel combustion and human driven land-use change is taken up by the ocean and the terrestrial biosphere 3,4 . In a top-down approach using atmospheric transport models, the global fossil fuel CO 2 emissions have often been presumed to provide good estimates of the strength of the land and ocean sinks 5 . But recent studies suggest that uncertainties in fossil fuel emission database could lead to significant biases in the optimized estimates of biospheric flux 6,7 . In fact, uncertainties of fossil fuel CO 2 emission estimates are growing because of increasing contributions from developing countries 3,4,8,9 . About 70% of the current anthropogenic CO 2 emissions is considered to come from urban areas that contain over 50% of the world population 10 , and thus accurate quantification of CO 2 emissions from urban areas is of particular importance. The existing and projected rate of urbanization is different for different parts of the world; many cities in Asia and Africa are expected to continue to see rapid growths in population while cities in developed countries have basically already stabilized 11 . For effective mitigation actions against climate change, various independent approaches must be used to reduce the uncertainties associated with citywide greenhouse gas emissions estimates, and one of those approaches can be provided by atmospheric CO 2 observations 12 . For this purpose, atmospheric CO 2 measurements focusing on urban areas have been examined in recent years by means of citywide in-situ ground measurement networks [13][14][15][16] or satellite measurements 17,18 . The former methodology provides dense and accurate data and the latter broad spatial coverage, whilst both also have limitations.
Since 2000s, atmospheric observation instruments onboard commercial airliners have successfully acquired extensive number of trace gas (including CO 2 ) data [19][20][21] . The CONTRAIL (Comprehensive Observation Network for TRace gases by AIrLiners) program, an ongoing project that measures atmospheric CO 2 and other trace gases using aircraft of Japan Airlines 20 , has obtained thousands of vertical profiles of CO 2 over many airports since 2005 [22][23][24][25] . Given that major airports are often located reasonably close to large cities and airlines optimize their flight destinations with priority on connecting between large cities, data collected by commercial airliners might provide useful information that can contribute to better assessments of greenhouse gas emissions from cities, complementing ground and satellite measurements. Here we present some statistical results from analysis of urban CO 2 emission signals contained in vertical profiles of atmospheric CO 2 data over airports. We present detailed analysis of data obtained at Moscow, Tokyo and 34 other airports, to offer a conceptual framework in which urban emission estimates can be obtained consistently with the statistical properties of the measurements. We show that the magnitude of atmospheric CO 2 variability increases with intensity of CO 2 emissions from the neighboring cities, demonstrating the usefulness of commercial aircraft-based measurements for urban emission studies.

Results
Moscow. We analyzed appearance of high CO 2 events in two altitude layers (4.0-4.5 km and 1.0-1.5 km) over Moscow Domodedovo Airport (DME) as a function of wind direction and speed (Fig. 1a,b). These plots indicate wind direction and speed in which high CO 2 events were observed in the measurement area. At altitudes of 4.0-4.5 km, almost no CO 2 enhancements in any wind direction and speed are evident (Fig. 1a), indicating that variability of CO 2 enhancement (defined as excess CO 2 ; see Methods) is small in the free troposphere. In contrast, at 1.0-1.5 km, it is seen that high CO 2 events were associated with winds mainly from the northwest sector (Fig. 1b). At altitudes of < 2 km, our measurement aircraft was primarily in the southeast of the Moscow metropolitan area (pink area in Fig. 1d) in all cases and influenced by CO 2 emissions from the city when advected from northwest. In addition, the fact that high CO 2 events appeared with relatively low wind speed (<15 m s −1 ) suggests that a presence of strong emissions close to the airport is likely. The Moscow metropolis covers an area of ~50 km in the northwest-southeast direction, a distance an air parcel with a speed of ~15 m s −1 traverses in ~1 hour. Longer residence time of an air mass over the metropolitan area would allow city emissions to be more integrated into the air mass, leading to results consistent with the observations of high excess CO 2 at low wind speeds. Histograms of excess CO 2 values at both altitudes (Fig. 1c) shows that the excess CO 2 at the lower altitude undergoes higher www.nature.com/scientificreports www.nature.com/scientificreports/ variability (i.e. larger spread of the histogram), and it is likely that the urban plumes contributed to the increased CO 2 variability with decreasing altitude.
tokyo. Similarly, we analyzed the data over Tokyo Narita Airport (NRT), located ~60 km to the east of the Tokyo agglomeration center (Fig. 2). The flight tracks at <2 km altitudes over NRT are roughly located on the eastern edge of the metropolitan area (Fig. 2d). We found that high CO 2 air masses were associated primarily with air parcels with relatively weak wind speed from the west of the measurement positions (Fig. 2b). The high excess CO 2 values (>10 ppm) were associated with <20 m s −1 westerly and <10 m s −1 easterly. The Greater Tokyo Area is ~100 km wide, the distance a westerly air mass at <20 m s −1 would pass over in > ~1.5 hour. Along with the DME analysis, this result suggests that an adequate residence time of >1 hour over the metropolitan area provides measurable accumulation of CO 2 up to ~1 km altitude in the overlying air mass. At the higher altitudes ( Fig. 2a), transport of high CO 2 was observed in the westerly at relatively high wind speed (>10 m s −1 ). It is again seen that the variability of excess CO 2 is larger in the lower layer (Fig. 2c). cities worldwide. We applied the above analysis to all of the 36 airports worldwide. At ~1 km altitude over many airports, we identified characteristic wind direction sectors where high CO 2 events were observed, and in each case the wind direction corresponded well to the geographical location of the urban agglomeration of the nearby large city (Supplementary Fig. S3-S36). High CO 2 enhancements were usually observed at relatively low wind speeds (roughly < 15 m s −1 ). Vertical distributions of excess CO 2 indicate that, in most cases, the ~4 km altitude is above the boundary layer, and thus the associated CO 2 histograms reflect the free tropospheric conditions. In contrast, the histograms in the 1.0-1.5 km layer show considerable broadening at many airports, indicating the presence of nearby surface fluxes that cause fluctuations in CO 2 in the overlying atmosphere. Therefore, at most airports, the lowermost altitude layers of the CONTRAIL measurements are under substantial influence of CO 2 emissions from the nearby urban area.
The above results point to the possibility of a relationship between the magnitudes of atmospheric CO 2 variability and of associated urban CO 2 emission. We used standard deviation (SD) of excess CO 2 as an indicator of the magnitude of the short-term CO 2 variability observed at all airports (Fig. S2). We observed SD exceeding 5 ppm at heights 1 km and lower over several airports, but at every airport, it decreases to <2 ppm at higher altitudes, i.e. the free troposphere (Fig. S37). Although the SD varies seasonally at heights less than or around 1.0 km over some airports (e.g. NRT), the seasonality in the SD is not common over most airports. For the data analysis below, we therefore calculated the annual SD values, regardless of the data gaps at individual airports.
We found that CO 2 variability at 1.0-1.5 km altitudes is largest over some Asian cities, and moderately large over some cities in Europe and North America (Fig. 3a-c). Large CO 2 variability is associated with urban areas with largest populations in the world; for instance, Tokyo (corresponding airports are NRT and HND) is ranked as the world largest urban agglomeration by population in 2010, and Delhi (DEL), Mexico City (MEX), Shanghai www.nature.com/scientificreports www.nature.com/scientificreports/ (SHA), and Osaka (ITM and KIX), over which we found the large SD values, are all listed in the top 10 population cities 11 . The most airports with small SD and vertical gradient are located in coastal areas, and the associated cities are relatively small in CO 2 emissions (Table S1).
We plotted the airport SD as a function of CO 2 emission from the nearby associated city for 2010 (Fig. 3d). Here we aggregated CO 2 emissions within 50 km to the east, west, north and south from a representative city center such as city hall, main station and local government building based on the ODIAC (Open-source Data Inventory for Anthropogenic Carbon dioxide, version ODIAC2015) dataset 26,27 . It is seen that, in the free troposphere (4.0-4.5 km), the SD is small (<2 ppm) over all the airports, irrespective of the intensity of CO 2 emission from a nearby city (R 2 = 0.16); in contrast, at the low altitude (1.0-1.5 km), the SD has a significant correlation with city CO 2 emission (R 2 = 0.73), suggesting considerable influence of the urban CO 2 emissions on our measurements. A significant correlation (R 2 = 0.68) is likewise observed for the data below 1 km, but with a larger scatter and a less compact correlation plausibly due to higher variability and less number of data. It is notable that there are some exceptional cities whose variabilities show excursions from the general trend; the SD values over Singapore (SIN), Los Angeles (LAX), and Chicago (ORD) are as small as those in the free troposphere; those over Hong Kong (HKG), Osaka (KIX), and Mexico City (MEX) are disproportionally high with large uncertainty. Here we used the ODIAC dataset, as it provides up-to-date gridded estimates of anthropogenic CO 2 emissions. It is noted that the use of different anthropogenic CO 2 emission datasets would not significantly affect the present result, given the fact that the existing CO 2 emission datasets (including ODIAC), with their major differences www.nature.com/scientificreports www.nature.com/scientificreports/ originating in the downscaling of national totals, correlate well with each other at the spatial scale of the current study (100 km × 100 km) or even smaller 28 .
Despite various factors that could contribute to the CO 2 variability over airports, our results strongly suggest that the intensity of the nearby urban CO 2 emissions is the primary component of the magnitude of the observed CO 2 variability (Fig. 3d). Our analyses of the wind pattern (Figs. 1, 2 and Supplementary Figs. S3-S36) indicate that the CO 2 variability at ~1 km altitude is, to a large degree, driven by the advection of high-CO 2 air masses influenced substantially by CO 2 emissions from the nearby city. It is likely that larger CO 2 emissions enhance contrast of atmospheric CO 2 mole fraction between their downwind and other locations, which makes fluctuations observed during multiple flights larger. The observed correlation may suggest that our lowest measurable altitude (1.0-1.5 km) is plausibly the height region where city-scale (~tens of km) emission imprint is distinct, and in which smaller-scale heterogeneity due to spatial pattern of emissions and local meteorology is relatively smoothed out. Previously, ground measurements in Paris showed that more urbanized part of the city showed larger CO 2 variations due to spatiotemporal variations of nearby anthropogenic emissions 15 . Moreover, the present result is consistent with a previous modeling study 23 , which showed, by analyzing the CONTRAIL CO 2 vertical profiles over NRT, that the magnitude of the atmospheric CO 2 variability in the boundary layer is highly sensitive to the magnitude of anthropogenic local fluxes.

Discussion
If we assume constant height of the convective boundary layer, one would expect a linear relationship between the surface CO 2 emission and the enhancement of the CO 2 mole fraction in the downwind boundary layer. In reality however, such emission-mole fraction relationship is complicated by the boundary-layer dynamics which undergoes significant daily and seasonal variations, causing fluctuations in the emission-driven CO 2 accumulation in the boundary layer. Our vertical CO 2 measurements were made under such varying meteorological conditions, likely a contributing factor in the variability in the SD values around the red line in Fig. 3d. In general, every airport has at least two opposite approach/departure directions along the runway, depending on the wind direction; flight routes are determined by airspace over the airport as well as by meteorological condition, thereby differ from flight to flight. We note that our observations might have been biased by specific flight routes under prevailing wind direction. The CO 2 mole fraction would be elevated in downwind of a neighboring city when the aircraft enters "urban CO 2 dome/plume" (Fig. 4a.c). In contrast, when the aircraft is positioned at upwind of the city, CO 2 emissions from the city would be hardly captured, except in cases where the aircraft flies down over the city (Fig. 4b,d). The varying spatial extent of the CO 2 dome/plume is another crucial factor in detecting urban emission signals. It is influenced by wind direction and speed, as well as the development and decay of convective boundary layer at the time of observation 15,24 . In this respect, the local time of the measurements is relevant, as commercial aircraft are scheduled at specific time of a day. Likewise, upwind biospheric CO 2 fluxes that undergo significant diurnal and short-term (shorter than seasonal) variations could contribute to the observed variability. Relative contributions of these various factors would be dependent on the measurement location, which might contribute to the variability around the linear trend in Fig. 3d. Case studies dedicated to data analysis of these contributing factors are required to identify the nature and extent of influence of the local surface fluxes to the observed variations; our earlier work provides an example for Delhi (DEL), India 24 . The variability-emission relationship in CO 2 (Fig. 3d) could not have been found without the commercial airliner in-flight measurements that enable worldwide vertical CO 2 scanning over large cities. We note that the relationship represents atmospheric CO 2 variability obtained from the last 10 years of aircraft data plotted against the 2010 CO 2 emission data. Given the above consideration (Fig. 4), it is reasonable to assume that this relationship would hold for other years of city emission data, and further analyses based on yearly or monthly time scale are required to establish this relationship on a much firmer ground. Currently the vertical profile data over Tokyo (NRT and HND), where the measurement is most frequent, might satisfy such requirement. The value of the results (Fig. 3), as interpreted within our conceptual frame (Fig. 4), becomes more obvious in contributing towards constraining anthropogenic CO 2 emissions for cities in developing countries where uncertainty of emission inventory is large and atmospheric measurements available for top-down emission estimation are sparse.
We also explored the possibility of a relationship of vertical gradient (i.e., CO 2 difference between ~1 km and the free troposphere) against city CO 2 emissions, with a hypothesis that strong ground emissions would elevate the near-surface CO 2 levels while the free troposphere is generally free from surface influence. In other words, vertical CO 2 gradient could represent the aforementioned linear relationship with associated surface emissions. We found large vertical gradients at some airports with strong nearby CO 2 emissions, however, we did not find a significant relationship between the vertical gradient and the city CO 2 emissions for various cities. The calculated vertical gradient was found to be considerably dependent on measurement time, season and location of individual airports. As discussed above, CO 2 variability over cities is influenced by various factors and, among them, the vertical extent of the CO 2 dome/ plume is of particular importance. It depends on the height of the convective boundary layer that varies seasonally and diurnally and influences the meteorological condition during vertical flight measurements. Seasonal coverage of the data and flight schedule are thus also relevant. Furthermore, geographical location of an airport relative to the nearby city varies significantly from city to city; thus, spatial representativeness of the observed vertical gradient differs among airports. Accordingly, for a global characterization of urban CO 2 emissions based on vertical aircraft measurements at world's major airports, as in this study, the vertical gradient would not provide a consistent relationship with corresponding city CO 2 emissions due to uneven measurement opportunities at various airports. We however consider that it could be a good indicator of CO 2 enhancement due to surface emissions for in-depth analyses at individual airports.
In summary, this study has demonstrated that the data taken by commercial airliners in airport proximity contain clear advective fingerprints of urban CO 2 emissions. Better understanding of transport processes of local surface fluxes to the atmosphere for individual measurement locations is also needed. Although the observed variability of CO 2 in the lower troposphere can be caused by many physical processes, the present study showed that the magnitude of the CO 2 variability at an airport reflects the amount of CO 2 emitted into the atmosphere from the corresponding city. The variability-emission relationship deduced from a relatively large number of vertical profile measurements could be helpful for evaluating surface fluxes and vertical propagation processes in the boundary layer simulated by various local to regional scale transport models. Further studies incorporating commercial aircraft data into urban CO 2 emission studies will contribute significantly towards better estimates of CO 2 emissions from major cities worldwide. Although vertical profile simulations by a regional transport model have been previously examined with commercial airliner carbon monoxide data 29 , the high-resolution vertical profile CO 2 data used in this study should be considered for incorporation in urban-scale inversion studies 30,31 or to be used in an independent validation of the emission estimations derived from ground-based data.
Although some recent projects have focused on urban CO 2 monitoring and quantification of CO 2 fluxes from cities, such as Paris 15 , Los Angeles 14 , Indianapolis 13 and Toronto 16 , we still have long ways to go, particularly when it comes to addressing key upcoming megacities that include those in developing countries 11,12 . In this respect, the advantage of commercial airliner is that measurements could provide a great opportunity to study many cities that have not been probed by the current ongoing urban monitoring projects and contribute to the reduction in uncertainties in the estimates of urban CO 2 emissions. The present study with Japan Airlines covers many cities in Asia where CO 2 measurements are sparse and a new set up of secured long-term ground-site measurements is difficult. Although installation of a greenhouse gas measurement system into a European airline is in progress 21 , further implementation, in particular into airlines of the United States, China and United Kingdom whose commercial aviation operations together accounts for about 40% of the world 32 , will significantly extend global monitoring capability of urban areas. Measurements onboard several airlines at same airport will provide vertical profile data at different local-time schedules, which will be also helpful for better quantification of surface fluxes.

experimental.
The CONTRAIL project deploys Continuous CO 2 Measuring Equipment (CME) and its in-flight measurements started on 5 November 2005. We refer to our earlier papers for details 20,22,25,33 . Briefly, the CME unit measures CO 2 mole fractions onboard the aircraft using a non-dispersive infrared gas analyzer (LI-840, LI-COR Biogeosciences). As of October 2019, installation of the CME is certified for eight Boeing 777-200ER and two Boeing 777-300ER aircraft of Japan Airlines (JAL). Once installed, the CME is operated automatically using the aircraft's flight navigation data from the ARINC 429 data bus until it is unloaded from the aircraft about two months later. The measured sample values are compared with two working standard gases (CO 2 in air) installed inside the CME that are traceable to the NIES (National Institute for Environmental Studies)-09 CO 2 scale. The NIES-09 CO 2 scale defines mole fraction of CO 2 in dry synthetic air in μmol mol −1 (reported in ppm in this paper). The results from the Round Robin intercomparison experiment show that the NIES-09 CO 2 scale differs from the WMO-CO2-X2007 scale (http://www.esrl. noaa.gov/gmd/ccgg/wmorr/wmorr_results.php) by <0.1 ppm. The CME data are recorded at 10-s intervals during ascent/descent (~100-m intervals in altitude) and at 1-min intervals during cruise (~15-km intervals horizontally) as well as in-flight aircraft position, static air temperature and wind data received from the ARINC 429 system. To avoid heavy pollution around airports, CME is not operated within 2000 ft (609.6 m) of the ground surface (this altitude was set to 1200 ft in early years). The analytical precision of the CME is estimated to be <0.2 ppm. Data analysis. Here we examined atmospheric CO 2 variations due to emissions from cities and local surface fluxes. It was assumed that the observed CO 2 variation is composed of variations on different timescales e.g. interannual, seasonal and shorter-term variations, where local scale signals of our interest most likely contribute to the last term. It is noted that this term also contains synoptic-scale variations which represent signals of larger scale fluxes in space (i.e. regional). We therefore needed to subtract longer-timescale variations according to our earlier analysis 25 as follows: lat lon alt t  CO lat lon alt t Trend CO at MLO t  ( ,  ,  , Here lat, lon, alt, t are the latitude, longitude, altitude and time of individual CME data points, respectively, and Trend CO 2 at MLO is the long-term trend curve calculated from the flask-based CO 2 data at Mauna Loa (MLO; 19.54°N, 155.58° W, 3397 m.a.s.l.), Hawaii, obtained from NOAA/ESRL/GMD (National Oceanic and Atmospheric Administration/Earth System Research Laboratory/Global Monitoring Division; available at ftp:// aftp.cmdl.noaa.gov/data/) with a digital filtering technique 34 . In general, the long-term CO 2 trend at MLO is representative of the large-scale clean atmosphere and has been used as a reference site; thus, ΔCO 2 in Eq. (1) gives an estimate of the climatological seasonal variations 25,35 (see Supplementary Fig. S1). Next we calculated excess CO 2 as follows: excess CO lat lon alt t CO lat lon alt t median CO airport alt bin t bin ( , , , Here airport refers to the nearby airport where the vertical profile data are taken, alt-bin and t-bin are the altitude and time-of-the-year bins to which individual data points belong. Individual data points are grouped into 500-m altitude and 14-day bins (giving a total of 26 time bins per year). Note that altitude in this study refers to level above ground of the nearby airport (above ground level, a.g.l.). By subtracting the seasonally varying median values according to Eq. (2), excess CO 2 represents detrended and deseasonalized variations and is analyzed below in terms of altitude and geographical locations. We also calculated SD of excess CO 2 in each time bin and found that its seasonality is not common to the data over most airports ( Supplementary Fig. S2). For statistical analysis, we selected 36 airports that contained sufficient number of vertical profiles over the 2005-2016 period; the number of vertical profiles ranged from 50 to >7000 (Supplementary Table S1).

Data availability
The CONTRAIL CME CO 2 data are available on the Global Environmental Database of the Center for Global Environmental Studies of NIES (https://doi.org/10.17595/20180208.001). The data are also available from the ObsPack data product (http://www.esrl.noaa.gov/gmd/ccgg/obspack/) and the World Data Center for Greenhouse Gases (https://gaw.kishou.go.jp/). The ODIAC emission data product is available from the Global Environmental Database operated by the Center for Global Environmental Research, National Institute for Environmental Studies (http://db.cger.nies.go.jp/dataset/ODIAC/).