Interactive comment on “ The dispersion characteristics of air pollution from world ’ s megacities ”

This paper presents results of forward simulations of megacity plumes with FLEXPART. The simulations are performed for a passive tracer as well as for Black Carbon with wet and dry deposition processes. The model results are used for 4 different sub-projects: to calculate transport distance, speed and height from the different megacities, to calculate deposition in the Arctic and in the Antarctic, and to calculate population exposure from individual megacity plumes both inside and outside of the source region.


Introduction
More than 50 % of the human population is living in urban areas (United Nations Population Fund report unit http://www.unfpa.org/swp/2007/english/introduction.html) and this fraction is still increasing, although urban areas currently cover only 0.5 % of the Earth's land area (Schneider et al., 2009).Extreme examples of urbanization are the so-called megacities with a population of several millions of people; ten million are often used as a threshold to define a megacity (Gurjar and Lelieveld, 2005).This large concentration of human activity causes unprecedented pressure on the environment, including the quality of air both inside the city (e.g., Gurjar et Published by Copernicus Publications on behalf of the European Geosciences Union.al., 2008) and on a regional and global scale.Some case studies (e.g., Guttikunda et al., 2003;Stohl et al., 2003;Real et al., 2008;Mena-Carrasco et al., 2009;Voss et al., 2010) have investigated regional and long-range transport of megacity plumes and the associated chemical processes, but more must be done to quantify the relative importance of their local versus large-scale impacts.
In a relatively recent paper Lawrence et al. (2007), called L07 in the remainder of the paper, investigated the differences and common features in the regional dispersion characteristics of plumes from a collection of the world's largest cities.In that work three artificial tracers, with a characteristic decaying time of 1, 10 and 100 days, were released and transported using a 3-D global Eulerian chemistry transport model.This study gave significant information on the common regional features of dispersion from these megacities and identified cases where the dispersion was dominated by local conditions.It is well known that for studying dispersion characteristics of plumes originating from a relatively localized source, a modeling approach based on 3-D Lagrangian particle models is a viable alternative to Eulerian models.Lagrangian models minimize artificial numerical diffusivity effects and can account for various subgrid-scale processes in an accurate way.On the other end, the Lagrangian nature of the model does not make it well suited to investigate the influence of the background on local megacity conditions, and this aspect is not considered here, although it can be an important aspect of the local pollution within the megacities, as recently discussed in Butler et al. (2012) and Crippa et al. (2013).
Based on these premises we decided to use the Lagrangian particle model FLEXPART (Stohl et al., 2005) for simulating the dispersion from a group of the world's megacities during the years 2003-2005, as established within the framework of the European Union project MEGAPOLI (Megacities: Emissions, urban, regional and Global Atmospheric POLlution and climate effects, and Integrated tools for assessment and mitigation, http://megapoli.info/).Our investigation has a similar scope as the study of L07 but uses a Lagrangian particle model with a 4-D sampling grid.The extension of the grid into the time domain allowed the use of additional metrics to describe the dispersion process of the plumes.Moreover, we specifically investigate long-range transport into the sensitive polar regions and into the upper layer of the troposphere and the stratosphere.Finally, we answer some of the questions pointed out by L07 by investigating the behavior of a passive tracer versus an aerosol-like tracer subject to removal processes by dry and wet deposition, and inferring information about population exposure by combining modeled near-ground concentration and population maps.
We mention that for this study, which compares several metrics for many megacities, only a modeling approach seems feasible.Consistent measurements collected for all megacities are unavailable to date.Moreover in this study the plumes are followed across many scales, and such plumes cannot be isolated in a straightforward way, in measurement data, from emissions from larger regions.Therefore, a detailed comparison with measurements for every megacity seems unfeasible and is not part of this study.Apart from the already mentioned L07, a pure modeling approach to study various impacts of emissions from an ensemble of the world's megacities has been used in several recent studies: Butler et al. (2009) studied the impact of megacities on global atmospheric chemistry, Butler et al. (2012) the megacities' ozone air quality under different future scenarios, Folberth et al. (2012) the influence of megacities on global radiative forcing, and Kunkel et al. (2012) investigated the global remote deposition of aerosol originating from megacities with an analysis similar to that of L07.Kunkel et al. (2012) in particular used five monomodal aerosols with diameters ranging from 0.1 µm to 10 µm and obtained similar deposition patterns for all the aerosols smaller than 2.5 µm.
In the next section a brief description of the FLEXPART model, the emission inventory and the meteorological input data will be given.In Sect. 3 the mean transport characteristics of the megacity plumes will be discussed.In Sect. 4 some metrics originally proposed by L07 will be calculated and compared with their original results.The contribution of megacities to pollution in the polar regions will be investigated in Sect. 5.In Sect.6, a metric of the large-scale versus local impacts of a megacity on human population will be defined.

Model description
For this study, the Lagrangian particle dispersion model FLEXPART (Stohl et al., 2005) was run for the period November 2002 to December 2005.FLEXPART is a widely used transport model and has been extensively validated with measurements.Some examples of validation relevant to the present study are against (1) large-scale tracer experiments in Stohl et al. (1998), (2) aircraft measurements of a large-scale intercontinental megacity pollution plume in Stohl et al. (2003), (3) aircraft measurements of a pollution plume transported from eastern Asia to Europe in Stohl et al. (2007), (4) an in situ observation of Asian pollution transported into the Arctic lowermost stratosphere in Roiger et al. (2011) and (5) a measurement of the transport of American megacity pollution into the North Atlantic in Neuman et al. (2006).The code and the documentation are freely available (see http://transport.nilu.no/flexpart/).It is worth recalling here that (i) FLEXPART uses the well-mixed criteria, introduced by Thomson (1987), to define stochastic differential equations with a Gaussian velocity distribution to simulate the effect of planetary boundary layer turbulence on particle dispersion and accounts for the vertical air density gradient (Stohl and Thomson, 1999), (ii) it includes a horizontal mesoscale meandering parameterization through an additional Langevin equation for the mesoscale velocity, (iii) it simulates convective exchanges using the scheme of Emanuel and Zivkovic-Rothman (1999), and (iv) the model includes dry/wet deposition processes.Thanks to the Lagrangian particle representation the model does not suffer from significant numerical diffusivity and the grid is used only to extract information from the particle field, which in general preserves very well any flow structures generated during the advection and dispersion process.This last characteristic, and the superior accuracy in the representation of subgrid dispersion processes, differentiates this class of models from Eulerian (fully grid-based) chemistry transport models, and makes FLEXPART particularly suitable to study dispersion from localized sources in the planetary boundary layer.
The model is driven here by operational meteorological data from the European Centre for Medium-Range Weather Forecasts (ECMWF) retrieved with a horizontal resolution of 1 • × 1 • and with 60 vertical levels in a hybrid coordinate system and a time resolution of 3 h.The model output grid has a horizontal resolution of 1 • × 1 • and 25 vertical levels ranging from 100 m in the layer closest to the ground to 3 km at an elevation above ground of 27-30 km.The grid extends also in the time space, thus partitioning the particles among 10 time (age) classes based on the time elapsed since the particle release, i.e., the time of emission.The age class durations used are expressed in days as 0 ≤ t < 1, 1 ≤ t < 3, 3 ≤ t < 6, 6 ≤ t < 9, 9 ≤ t < 14, 14 ≤ t < 20, 20 ≤ t < 27, 27 ≤ t < 34, 34 ≤ t < 41 and 41 ≤ t < 48.Note that, because of their different duration, the mass content of a perfectly conserved tracer differs for the different age classes.When a particle reaches 48 days of elapsed time, it is removed from the simulation; therefore even the conserved scalar (BCtr, see below) has a maximum time of residence in the atmosphere of 48 days.The model output is stored every 5 days as a 5-day average concentration c (x, y, z, t) or, written in a discretized form, c ij k,it , where it is used as index for the age and i, j, k are indices for the longitude, latitude and elevation above terrain level (ATL).
Aerosols are removed in FLEXPART by dry and wet deposition.The dry deposition is simulated according to a resistance approach (Wesely and Hicks, 1977) considering the aerodynamic resistance, the quasilaminar sublayer resistance, and gravitational settling (Slinn, 1982).For the aerosols considered in our simulations (see below), we assumed a logarithmic diameter distribution with a geometric mean diameter of 2.5 × 10 −7 m, a geometric standard deviation of 1.25, and a density of 1.4 × 10 3 kg m −3 .Wet deposition distinguishes between in-cloud and below-cloud scavenging (McMahon and Denison, 1979;Hertel et al., 1995).Scavenging coefficients are determined primarily as a function of cloud liquid water content and precipitation rate, respectively.Note that FLEXPART has no aerosol processing scheme and that details such as the conversion of some aerosols from an initially hydrophobic to a hydrophilic state (Vignati et al., 2010) are ignored.

Megacity emissions and deposition timescales
The emission inventory is extracted from the original Representative Concentration Pathway (RCP) 8.5 emission scenario for the year 2005 developed by the MESSAGE (Model for Energy Supply Strategy Alternatives and their General Environmental impact) modeling team and IIASA (Integrated Assessment Framework at the International Institute for Applied Systems Analysis), Austria (see http://www.iiasa.ac.at and Riahi et al., 2007).The global emission data used here has a 1 • × 1 • resolution and was obtained from the original data set at 0.5 • × 0.5 • .The mask used to identify the megacities' emissions was developed by Butler et al. (2009) and subsequently extended as reported in Folberth et al. (2012) to include extended agglomerations such as the Po Valley and other megacities, as agreed in the framework of the MEGAPOLI project.Any simulated megacity occupies one or more grid cells of this emission data set as listed in Table 1 and illustrated in Fig. 1 and, it is possible to see that both single-grid-cell as well as multiple-grid-cell megacities are included.From the RCP 8.5 data set, the emissions of several species can be extracted.It should be noted in Fig. 1 and in Table 1 that the definition of the megacities is not perfectly consistent since extended areas, like the agglomerate of New York, Boston and Washington (here just called the New York megacity for simplicity), the Dhaka megacity and the Po Valley are compared to more localized areas, like Chicago or Paris.However, the emission and population are consistently extracted using the same mask.
A first estimate on the relative importance of megacity emissions on a global scale can be deduced immediately from Fig. 2, where the cumulative emission contribution for several chemical species (black carbon, BC; carbon monoxide, CO; nitrogen oxides, NO x ; sulfur dioxide, SO 2 ) for all the megacities up to a given degree of latitude is scaled with the total cumulative species emission (including nonmegacity emission) to the same degree of latitude.For example the total cumulated emission starting from the South Pole and up to a specific latitude, i lat , could be written as , where E (i, j ) is the emission in grid cell i, j .Similarly, the cumulative contribution for all the cities starting from the South Pole to latitude i lat could be written as where N city is the number of cities considered, N cells the number of cells representing the city and δ is a function with a value of one within the city limits and zero outside.
The data have been separated in the Northern and Southern Hemisphere (NH and SH), and the cumulative contribution Table 1.Several quantities for all the megacities being part of this study are reported.T li and T lo are two decaying timescales obtained for the deposition process (see Fig. 3 and the text for details).PWC is the population-weighted ground level concentration and measures the concentration to which the population inside the city is exposed on average.GLC × Pop is the product of population and ground level concentration.Other quantities are self explanatory.S latitude megacities in the SH contribute more to total emissions than their counterparts in the NH.In the NH, there is a strong seasonal variation at higher latitudes especially for CO and BC, which is due to intense open biomass burning in spring and summer, so that megacity emissions contribute less in summer.

Properties
In the autumn and winter months, a single city (St.Petersburg) accounts for almost 10 % of the total emissions north of about 60 • N for the BC considered in this study.It can be of interest to note that, adding to the list of cities an emitting area extending from a lower left corner of (59 • N, 10 • E) to an upper right corner (61 • N, 26 • E), the emission fraction goes up to 36 % of the total emission over 60 • N for BC in autumn and winter.This area contains five large urban areas in the Baltic and Scandinavian countries: Tallinn (Estonia), Helsinki (Finland), Stockholm and Uppsala (Sweden), and Oslo (Norway).Since on the scale of this study these cities are relatively close to the investigated megacity of St. Petersburg, the results obtained for this latter city could be of guidance for the behavior of emissions that in the autumn and winter contribute 36 % of the total world emissions north of 60 • N.These cities' contribution is lower in the summer, but considering that during these months the concentrations of BC in the Arctic lower atmosphere are lower due to stronger deposition processes (e.g., Hirdman et al., 2010), the winter peak of these urban emissions could be significant.However, it is important to mention here that emissions from individual cities can differ, often by a factor of two, between different emission inventories (Butler et al., 2008) and this uncertainty is implicitly contained in our analysis.
The further discussion in this paper will be based only on the emissions of one species, BC, which we consider as representative for all combustion-related emissions.This species will be used for two different artificial scalar tracers: one considered as a conserved tracer (BCtr) and not subject to any removal process and another one (BCdp) with dry and wet deposition processes included in the simulations.
For studying the impact of pollution on human population we used the GPWv3 data set of the Center for International Earth Science Information Network (CIESIN), Columbia University -available at http://sedac.ciesin.columbia.edu/gpw (January 2011).The population of a city was extracted from the GPWv3 data set, using the same megacity mask used for the emission; the resulting population for the cities is listed in Table 1.
Before proceeding with a detailed analysis of the plume characteristics it is useful to explore the effects of the modeled wet and dry deposition processes on the suspended mass in the atmosphere as a function of time after emission ("age").The deposition processes reduce the mass content significantly with age but differently for the different cities considered.In most figures in this paper, megacities have been divided into NH and SH megacities, and the results of these two groups will be in general plotted separately to improve the clarity and readability of the figures.

Dispersion characteristics of transported plume 6
In this section we start the discussion about the megacity plumes with an analysis of their 7 physical dispersion characteristics both horizontally and vertically.To characterize the horizontal 8 dispersion of the plumes, the model output was first averaged over the time period of interest 9 (e.g. three winter seasons or all three years) 10 a logarithmic scale, while Fig. 3c, d shows the same ratio on a linear scale.This indicates how dry and wet deposition processes deplete the atmospheric mass of BCdp with time.From Fig. 3c, d it is apparent that, initially, deposition processes are more active for European and tropical cities.However, after a few days the plumes originating from megacities located at midlatitudes in Asia are those characterized by the higher deposition rates, related to strong scavenging in the North Pacific storm track (Fig. 3a, b).For the last few age classes, European and North American plumes have again higher deposition rates.This complex behavior reflects the different three-dimensional distributions of the plumes originating from different geographical regions.For very large ages, all the plumes should converge towards the same deposition rates typical of the Northern and Southern Hemisphere, while initially the deposition rate is strongly influenced by the local processes close to the source.This can be better appreciated defining two timescales obtained by fitting an exponential decay to the data on a logarithmic scale, T lo , (based on Fig. 3a, b) and to the data on a linear scale, T li (based on Fig. 3c, d).These timescales are reported in Table 1.T lo is an estimate of the average e-folding deposition timescale for the 48-day time window and the source location.T li being more related to the initially higher values of BCdp/BCtr, gives an idea about the time needed for the mass emitted from a specific source to be depleted to a small fraction of its initial value, and therefore it should strongly depend on the local processes in the boundary layer close to the source.T li varies between one and ten days, while T lo varies between about six and nine days, but with most of the cases between about six and eight days, indicating the approximate averaged global e-folding lifetime of the BCdp aerosols in FLEXPART.
An important feature of including age classes in the analysis is that many of the results reported here involving a mass assigned to the age classes can be extended to consider the effect of a superimposed decay process (e.g., radioactive decay, or chemical destruction) just by applying, a posteriori, the decay -e.g., exponential characterized by an assigned timescale.This is true for both BCtr and BCdp.In the first case the superimposed decay is applied to a substance not undergoing deposition, and in the second case to a substance undergoing wet and dry deposition.

Dispersion characteristics of transported plume
In this section we start the discussion about the megacity plumes with an analysis of their physical dispersion characteristics both horizontally and vertically.To characterize the horizontal dispersion of the plumes, the model output was first averaged over the time period of interest (e.g., three winter seasons or all three years): temporary asymptotic distance.In fact Figure 4(c) shows that for BCtr all the plumes but the 7 Lima and Jakarta plumes are transported more rapidly compared to the NH plumes.This 8 behavior is likely related to the stronger westerlies in the SH.Jakarta shows a typical behavior of 9 a tropical megacity as previously discussed for the NH, while Lima shows the slowest horizontal transport speed of all the megacity plumes investigated.(c, d).The smaller panels on the left are a zoom for 0-3 days using a logarithmic scale.
Here N indicates the number of 5-day output intervals in the period of interest, x, y, z are the three space dimensions, and t indicates the elapsed time from the release (i.e., the age).
From this time-averaged concentration a vertically integrated concentration field was introduced: Here z k indicates the vertical grid level extension and N z (= 25) is the number of vertical levels in the FLEXPART output grid.From this field an advection measure was obtained as the average distance of the mass from the source location, d (t), in the following called transport distance, where d ij is used to indicate the distance, along the Earth's surface, between the center of the grid point with indices (i, j ) and the source position.Since some sources are composed of more than one single grid cell element, the source position was defined as the grid element containing the maximum ground level concentration averaged over the extension of the first age class (i.e., 0 < t < 1 day).ij is the area of a grid element with indices (i, j ), and ( cz ) ij is the concentration associated with the grid element.It should be noted again that, since the grid extends also in the time dimension, any of these measures is defined as a function of age for the intervals defined above.In a similar manner, the averaged vertical position, h (t), of the horizontally integrated concentration ( cxy ) k as a function of time can be introduced, where k is used for the vertical grid cell index.
The results for the distance from the source location, d (t), for BCtr are shown in Fig. 4a and c, respectively, for the NH and the SH, and those for BCdp are shown in Fig. 4b and  d.The smaller left panels are a zoom into the age classes 0-1 and 1-3 days and use a logarithmic scale for the plot.For BCtr in the NH, it can be seen that initially the plumes subject to the fastest advection are those emitted at higher latitude and associated therefore with European and North American megacities, while the plumes emitted by tropical or subtropical megacities have the slowest initial transport.However, the situation changes, and for elapsed time greater than 10 days all the higher latitude cities show comparatively short transport distances.Conversely, the tropical and subtropical megacity plumes then become associated with the   (2) cities at lower latitudes have a larger possible asymptotic distance for a uniform scalar distribution in the hemisphere; and (3) the tropics are characterized by lower near-surface wind speeds but a high frequency of deep convection, so that the dispersion is at first slower than for higher-latitude cities but after a few days significant amounts of tracer in the tropics are transported by deep convection to the upper troposphere, where winds are faster.This is illustrated in Fig. 5a and c, which shows that tropical cities and particularly southeast Asian cities (e.g., Manila) are associated with a very fast increase of the plume average elevation.
The zoom figures for the horizontal transport behavior between 1 and 3 days (Fig. 4) show that in the NH the cities with the initial fastest plume transport are London, New York, Osaka and Tokyo.These are all cities close to the ocean which is therefore an important factor in determining a fast early phase of plume transport.Among the European megacities the Po Valley shows a unique behavior since the tracer is trapped there by mountains on three sides, which strongly limits the early phase of transport.One consequence of this is the high pollution levels usually registered in this area (e.g., Bigi et al., 2012).However, overall in the NH the megacity with the slowest horizontal transport up to 10 days is the tropical city of Lagos, which also has severe air pollution problems (Baumbach et al., 1995).
Removal processes lead to a reduction of BCdp mass in the atmosphere.The lifetime T li of the BCdp tracer for the various megacities ranges from 1.0 days for Jakarta to 10.5 days for Cairo (Table 1), reflecting the very different efficiencies of wet scavenging in tropical and desert regions.This large difference in the scavenging has also been reported and quantified, based on climate classes, in Kunkel et al. (2012): for small aerosol particles, a large remote deposition fraction over 1000 km was found for cities located in arid regions (average of about 75 %), while a much lower value (36 %) was found for cities in tropical regions.
Removal processes also change the transport characteristics of the BCdp plumes, and these changes are sometimes difficult to understand (Fig. 4b, d).For example, the preferential removal of mass during upward transport associated with precipitation means that initially the vertical transport of BCdp is slower compared to BCtr (Fig. 5).As a consequence, initially the horizontal transport is dominated by the slower low-level winds, leading also to shorter horizontal transport distances of BCdp compared to BCtr (see left panels in Fig. 4).After a few days, however, due to preferential wet scavenging in the lower troposphere versus lack of precipitation in the upper troposphere and lower stratosphere, for the remaining BCdp the fraction located in the he averaged (over all cities) horizontal distance of centers of mass, D(t), and center of mass, H(t), for BCdp over BCtr.upper troposphere and stratosphere increases.Thus, the centroid altitude of the BCdp tracer increases throughout all remaining age classes (Fig. 5b, d), and the horizontal transport distance also catches up with that of BCtr (Fig. 4b).This is summarized in Fig. 6a and b, which shows the ratio, for BCdp versus BCtr, of the horizontal transport distances (D (t) = Notice that these results for H (t) do not imply that a larger fraction of the total emissions can reach the upper troposphere for BCdp than for BCtr.However, after some time a larger relative fraction of the remaining BCdp than of the remaining BCtr resides in the upper troposphere, as demonstrated in Fig. 7, which shows the fraction of the total BCtr and BCdp mass suspended in the atmosphere which is contained between 9 km and the top of our sampling domain (30 km).For BCdp, this fraction increases throughout the 48 days considered in our study, and for some cities it approaches the value of 1 (Fig. 7b and d).
The SH megacity plumes have a similar behavior compared to their NH counterparts.Midlatitude megacity plumes show initially a faster horizontal transport, while ultimately the tropical and subtropical latitude megacity plume have a larger transport distance.However, the transport for midlatitude cities is here faster than in the NH, and the plumes reach faster their temporary asymptotic distance.In fact Fig. 4c shows that for BCtr all the plumes but the Lima and Jakarta plumes are transported more rapidly compared to the NH plumes.This behavior is likely related to the stronger west-erlies in the SH.Jakarta shows a typical behavior of a tropical megacity as previously discussed for the NH, while Lima shows the slowest horizontal transport speed of all the megacity plumes investigated.
4 Statistical measures of regional pollution export L07 defined several metrics to characterize the dispersion of the megacity plumes and their potential impact on a regional scale.Here we will calculate within our modeling framework two of their metrics of pollutant export over long range (EL-Rcol and ELR1km) and compare our results with theirs.EL-Rcol is defined as the fraction of the total BC mass in the atmosphere which is more than 1000 km away from the source at any altitude.ELR1km is the same as ELRcol but for the mass contained in the first 1 km ATL.Contrary to L07 we considered only whole grid cells and did not split grid cells according to their fractions within or beyond the 1000 km radius from the source; however, our grid has finer resolution (1 • versus 1.9 • ), which compensates for this less accurate treatment.Also, the simulations of L07 were for the year 1995, while our simulations are for the years 2003-2005.Also, as before, we will investigate the effect of the deposition, thus extending the results of L07.
Our age class treatment generalizes these metrics to any decay time within the 48-day range considered.The ELR measure has been calculated as cumulated mass exported up to a specific age normalized with the total mass of the substance suspended in the atmosphere cumulated up to the same age.This means that the value of the last age class corresponds to the total exported mass fraction.
L07 report in detail the ELR1km values for a scalar with a 10-day decaying timescale.Figure 8a shows that the value of ELR1km for BCtr becomes, after about one week, almost constant for most cities.This means that all tracers with a decay timescale larger than about one week will behave similarly to a purely passive tracer without decay.This simply reflects the partitioning of the tracer between the lowest 1 km and the rest of the atmosphere, as basically all of the tracer will be exported beyond 1000 km after a week.However, using the discretized time dependence information present in the age class partitioning we can explicitly apply the exact exponential decay process with a 10-day timescale to allow a more precise comparison of our results to those of L07.Table 2 compares the ELR1km results for the 48-day age class for a conserved tracer with those for a tracer with a 10-day lifetime (BC10).The results of L07 for BC10 are also reported.It is possible to see that the difference between BC10 and the conserved tracer, BCtr, is minor in most of the cases.However, in general and as expected, the results obtained for the decaying scalar agree better with the results reported by L07.Considering that the results of L07 were for a different year and used different meteorological input data (NCEP reanalysis) and a different transport model (MATCH), the L07 defined several metrics to characterize the dispersion of the megacity plumes and their 7 potential impact on a regional scale.Here we will calculate within our modeling framework two 8 of their metrics of pollutant Export over Long-Range (ELRcol and ELR1km), and compare our 9 results with theirs.ELRcol is defined as the fraction of the total BC mass in the atmosphere 10 which is more than 1000 km away from the source at any altitude.ELR1km is the same as 11 ELRcol but for the mass contained in the first 1 km ATL.Contrary to L07 we considered only 12 whole grid cells and did not split grid cells according to their fractions within or beyond the 13 1000 km radius from the source, however, our grid has finer resolution (1° versus 1.9°), which 14 compensates for this less accurate treatment.Also, the simulations of L07 were for the year 1995 15 agreement between the two studies is very good.The values of ELR1km obtained in our study range from 6.3 % to 30.1 % compared with 3.2 % to 34 % in L07.
A scatter plot (Fig. 9) shows that the ELR1km values for European, North American and midlatitude Asian cities agree reasonably well in the two studies.Southern Hemisphere megacities have generally higher values in the present calculation based on FLEXPART.The worst agreement is achieved for the Indian subcontinent, Asian tropical and South American megacities, although the differences are not systematic.Jakarta stands out with ELR1km values of 7 % and 3.2 % in FLEXPART and MATCH, respectively.This difference must be related to the different convection schemes used in FLEXPART (Emanuel and Zivcovic-Rothman scheme, implemented by Forster et al., 2007) and MATCH (plume ensemble scheme, proposed by Lawrence and Rasch, 2005).The low fraction in MATCH means that in this model the BC10 tracer must have been vented out of the boundary layer very efficiently, with only small compensating downward fluxes.This result is obtained, in somewhat less extreme form, also for other tropical cities such as Manila, Bangkok, Sao Paulo and Rio de Janeiro.On the other hand, FLEXPART has lower ELR1km values than MATCH for the cities on the Indian subcontinent.This is probably due to different representation of the winter monsoon in the two models, where FLEXPART shows a relatively rapid lowlevel southward transport to the intertropical convergence zone, followed by rapid lifting (Stohl et al., 2002).
Figure 8b shows the metric ELR1km computed for BCdp.The export efficiency for the BCdp tracer depends not only on the combined action of horizontal wind and vertical transport, but also on the removal by wet and dry deposition.In general the effect of deposition is to reduce the value of ELR1km, and the reduction is stronger for higher local deposition rates, i.e., low value of T li .However, the megacities in the Indian subcontinent plus the city of Cairo have a noticeable increase of ELR1km for BCdp with respect to that of BCtr.A possible mechanism for this happening is if precipitations act mainly on rising air masses outside the 1000 km  radius from the source.This seems consistent with the winter monsoon behavior for the Indian subcontinent, as described above.The city of Cairo is instead characterized by an arid local climate (it is the city with the highest T li ), and therefore, since there are very low levels of precipitations close to the city, it is plausible that precipitations act mainly on rising air masses outside the 1000 km radius.This is in agreement with the analysis of the deposition for small aerosol particles reported in Kunkel et al. (2012) which, for Cairo, found that deposition happens very far away from the source after the plume has travelled above the Sahara and has reached the inner tropical convergence zone.In general this factor adds variability to the possible range of ELR1km for BCdp, with values spanning about a factor of ten, compared to only a factor of four for BCtr.The two lowest total values of ELR1km for BCtr and BCdp are found for Jakarta and Lagos (see Table 2), and it is interesting to examine why.Lagos is the city with the slowest horizontal transport of the plume during the first 3 days after tracer release (Fig. 4a, b, left).Jakarta on the other hand is the second slowest in the range 3 < t < 14 days (Fig. 4a, b, right).This difference explains also why Jakarta has an initial peak in the value of ELR1km for BCtr (Fig. 8a).When deposition is accounted for, Jakarta is the city with the shortest T li , and therefore the deposition generates a stronger decrease of ELR1km compared to that of any other city (including Lagos).As a consequence the difference between the two lowest values of ELR1km becomes more marked for BCdp.It is also interesting to discuss why Mexico City has initially the lowest value of ELR1km and why it has a quite low overall ELR1km, despite having a quite fast overall horizontal transport (Fig. 4a, b, right).This city is characterized by slow initial transport (t < 1 day, see Fig. 4a, b) and has the highest initial plume elevation among all cites (Fig. 5); as a consequence, the plume is initially transported slowly and subsequently, when the transport becomes faster, it is mostly transported at an elevation above 1 km ATL.The two top low-level exporters irrespective of deposition are the cities of Moscow and Saint Petersburg, and these are cities with a quite fast initial transport (Fig. 4) but also and especially very low average elevation (Fig. 5).
As previously reported by L07, results for the total mass exported to a distance of more than 1000 km (ELRcol) are considerably less variable.In particular for the 48-day age class we find a very narrow range of values from 85 % to 95 % (Fig. 10).This is expected since in this case the values should converge approximately towards the same asymptotic limit of 98.8 % for a scalar well mixed in the hemisphere, where an area bounded by the 1000 km radius covers 1.2 % of the hemisphere.The values of ELRcol for BCtr in the twoweek age class goes from about 60 % to 85 %.This is indicative of the value of a scalar subject to a superimposed decay of a similar duration, and indeed this range is quite similar to the cumulated value of ELRcol for BC10, which ranges from about 59 % to 82 %.The range of values of BC10 is again in good agreement with the results of L07 for the 10-day decaying scalar for which they reported a range from 62 % to 84 %.
It is interesting to note that Lima, which has clearly the lowest value of ELRcol for any age, was not among the cities having a very low ELR1km.This city has an extremely slow horizontal transport (Fig. 4), which is perfectly reflected in the low ELRcol value.However, contrary to Jakarta and Lagos, which have a very low ELR1km, Lima is not characterized by a very fast uplifting of the plume (Fig. 5), and this increases the value of ELR1km.The deposition modifies again the picture by widening the range of values.However, in this case the interpretation is simpler with respect to ELR1km since deposition systematically decreases the value of ELRcol.The cities with the highest values of ELRcol are those characterized by a high value of T li , and Cairo (which is the city with the highest T li ) becomes the top exporter.The fact that Cairo was not the top low-level exporter for BCdp underlines again that for ELR1km it is very important that the plume is transported at low elevation, which is not the case for Cairo.

Air pollution and BC deposition in the polar regions
Analysis of BC concentration in the polar regions is important due to the positive radiative forcing of this species both in the atmosphere and when deposited on snow or ice and also due to the high efficacy of this forcing for climate warming (e.g., Hansen and Nazarenko, 2004;Flanner et al., 2009).We note that a recent paper by Folberth et al. (2012) gives some evidence that the megacities' global radiative forcing is modest, also because of compensating effects of the emitted pollutants.However, this, as discussed by same authors, does not mean that there is not a regional impact on climate (e.g., of the northern European megacities in the Arctic).Our simulations allow us to quantify the efficiency with which BCdp and BCtr originating from megacities in different parts of the world can reach the polar atmosphere and deposit on the ice/snow cover.This study is somewhat complementary to the works of Stohl (2006), Stohl and Sodemann (2010) and Hirdman et al. (2010), who studied the possible sources of BC in the Arctic and Antarctic atmosphere from a receptororiented perspective.Here, we take a source-oriented view by following the emissions from the world's megacities.Recently Kunkel et al. (2012) carried out an extensive study about aerosol deposition from large urban areas, but they kept a more general perspective than our present study, without any particular focus on the polar regions.

Air pollution
Figure 11a and b show the total BCtr and BCdp mass in the lowest 1 km ATL of the NH polar atmosphere (north of 66 • N), scaled for the mass emitted in the respective age class interval, while Fig. 11c, d provides the same information for the SH (south of 66 • S).As discussed extensively in Stohl (2006), the contribution to mass in this layer in the NH depends mainly on the latitude of the source, with largest contributions from northern Eurasia.This is confirmed here for the BCtr results.Specifically, we find that releases from Saint Petersburg give the largest fractional contributions of all megacities and for all age classes.As much as 10 % of the emissions can be located in the Arctic already as early as about 2 days after the emission.In fact, the contribution decreases after about 2 weeks, as the Saint Petersburg plume is initially preferably transported to the north and subsequently disperses into the entire NH.Similar behavior could be expected for pollution export from other Baltic and Scandinavian cities that are relatively close to Saint Petersburg.This group of cities is consequently likely to be the largest "localized" contributor of anthropogenic pollution in the lower layer of the Arctic atmosphere for substances with lifetimes on the order of days to a few weeks.Furthermore, while Moscow makes a smaller BCtr fractional contribution to the pollution in the Arctic atmosphere during the first 3 days after emission, the contributions at later times are almost identical to those of Saint Petersburg.The age class sampling allows appreciating that, for any substance with a decaying timescale of less than about 10 days, the contributions from non-European cities are very small.Only after this elapsed time do contributions from North American megacities start to be more relevant, and soon after also from east Asian midlatitude cities. Contributions from tropical locations are, even for the last age class considered (48 days), an order of magnitude smaller than contributions from Saint Petersburg.Table 3 (column 1 and 3) reports the cumulative amount of mass irrespective of the age classes and expressed both as mass and as a percentage of total emitted mass.The table confirms that cities in Europe have significantly larger fractional contributions, but in terms of absolute numbers (i.e., taking into account also the different emission strengths) the extended New York megacity is the largest contributor and Tokyo, Seoul and Beijing have numbers comparable to those of most European cities.However, we remind the reader that, as shown by the age class plots, these cities have smaller contributions relative to European cities when the lifetime is significantly shorter than 48 days.When taking into account deposition processes (BCdp), the fractional contributions from all cities become much smaller (Fig. 11b).However, the reductions are particularly large for non-European cities, and, thus, significant contributions to the BCdp loading in the lowest 1 km of the Arctic atmosphere are confined to northern European sources with small contributions from North American cities.Saint Pe-tersburg and Moscow give similarly high contributions in the age classes going from 2 to 20 days, although both London and the Ruhr Valley are also relevant.The cumulative numbers reported in Table 3 (column 2 and 4) show that, with our definitions of the megacities' extension, for BCdp all the megacities in Asia contribute less than the sole city of Saint Petersburg and that the extended New York megacity has a similar contribution to Saint Petersburg despite having almost 20 times the population.
Figure 12 reports the fraction of mass of BCdp and BCtr in the Arctic free troposphere and stratosphere, 2-9 km (a, b) and 9-30 km (c, d) above sea level, respectively.For Bctr, the contributions in the free troposphere are similar to those observed in the first layer (0-1 km) but with higher values from midlatitude Asian and North American cities, which are as large as those of European cities as early as two weeks after the release.Also, with the increase of the plume age, the contributions from sources in the tropics become noticeable.This again suggests that, if a superimposed decay process is investigated, a decaying timescale of about a week would give clearly the largest contributions from European and North American cities, while longer decaying timescales would give substantial contributions from Asian megacities.If deposition is accounted for (Fig. 12b), the contribution from European megacities becomes more significant even for longer superimposed decaying timescales.However, only in the case of substances with a superimposed decay of less than 3 days would the contribution from midlatitude Asian cities be negligible.This is in contrast to the 0-1 km ATL layer, where European sources are totally dominant.
Moving to the stratosphere, Fig. 12c shows that the most efficient transport of BCtr mass into the Arctic stratosphere is from the midlatitude Asian and North American megacities.The differences between these and other sources are better appreciated on a linear scale (insert in Fig. 12c).Looking in more detail it is possible to see that initially (age class of up to 3 days) the most efficient in contributing to the BCtr mass in this upper layer are Northern European cities but soon thereafter midlatitude Asian and North American cities become noticeably more efficient.This suggests a quite complex scenario in case a superimposed decay is considered, with dominant contribution from different regions of the NH depending on the lifetime of the chemical species.Similar lifetime dependence could be observed in the case of BCdp, although in this case the fractions of mass reaching these layers of the atmosphere are much lower.Interestingly the Po Valley megacity becomes in this case the most efficient source in transporting mass into the Arctic stratosphere.In Table 3 (last two columns) the total mass transported in the polar stratosphere for BCtr and BCdp is reported, and it is possible to see how the Po Valley stands out as the largest contributor to BCdp in the Arctic stratosphere.
The current results -reporting a contribution mainly from European megacities to the mass of BC in the lowest kilometer of the Arctic atmosphere and a progressively increasing  reports (1) the cumulative value of mass of BCdp deposited (Dep66) in the polar regions (|Lat|>66 • ); (2) the mass of BCtr and Bcdp suspended (Sus66) in the polar atmosphere up to an elevation h<1 km above terrain; (3) The mass of BCdp deposited over land during winter in the northern hemisphere for Lat>40 • (Dep40l) (4) the mass in the polar stratosphere (Strat66).The mass is reported in tons, kg and as a fraction (%, ‰) of the total averaged emitted mass within the 48-day age class limit.contribution of North American and east Asian megacities when moving up to the free troposphere and the stratosphere -are in agreement with the findings of Hirdman et al. (2010), who reported a substantial contribution to BC measurements at lower altitude stations only from northern Eurasian sources, while significant contributions from more remote sources (e.g., east Asian) were observed for measurements on top of the Greenland ice sheet.

Deposition
For quantifying BC deposition in the polar regions we consider the accumulated mass over an area of interest during a period (i.e., year, season).When normalized quantities are shown, this is done considering the total mass available for deposition during the period considered, given by the initial mass suspended in the atmosphere at the beginning of the period plus the mass emitted during the period.Figure 13 shows the normalized (left panels) and unnormalized (right panels) deposition north/south of 66 • N/66 • S within 48 days from the release.The values cumulated over all age classes are reported in Table 3. Saint Petersburg has the largest fraction of its BCdp emissions deposited in the Arctic (9.4 %), followed by Moscow (5.6 %) and the Ruhr Valley and London (both about 4 %).However, the absolute deposition in the Arctic is comparable for these cities, and also the Po Valley and Paris cause significant deposition values.The North American megacities of New York and Los Angeles give normalized deposition values that are more than an order of magnitude smaller, but because of the larger BC emissions of these American cities, the actual deposition values are only about 2 to 3 times lower than the contributions of European megacities.Absolute and relative contributions from megacities on other continents are much lower.This means, as shown in Table 3, that the sole city of Saint Petersburg (population 2.6 million based on the GPWv3 database and with our definitions of his extension) produces a BCdp deposition in the Arctic that is marginally larger than the BCdp deposition from all the megacities in Asia (population roughly 500 million).This seems to confirm again that most of the Arctic BC deposition from megacities must be attributed to northern Eurasian sources.
Looking in detail to the deposition fractions for Asian megacities in Table 3, it is interesting to note that Teheran, with a value above 1 %, stands out with respect to all other cities of comparable or more northern latitude.This is related to the arid conditions that allow for a considerable fraction of the mass to be exported away from the source (T li = 8.2 days).This result agrees well with the analysis of Kunkel et al. (2012) that found for cities in arid climate the largest remote aerosol deposition potential.
The BCdp deposition over land during the winter season for latitudes north of 40 • is also reported in Table 3 (columns  7 and 8).During winter, a significant fraction of land above this latitude may be covered by snow, and the BC deposited may have an influence in accelerating the melting of the snow in spring (e.g., Quinn et al., 2008;Flanner et al., 2009).For this region, the fractional contributions from European and North American megacities are comparable, and the Asian megacities of Beijing, Tianjin, and Teheran give significant fractional contributions followed by Cairo in northern Africa.Similar considerations are valid for the absolute mass deposition, and in this case New York City emerges as the largest contributor to deposited BCdp.

The Southern Hemisphere
Although there are just eight megacities in the SH, according to the emission inventories used, they represent about 12 % of the total emission (see Fig. 2) in the SH from the South Pole to a latitude of 25 • S. The export of BCtr and BCdp from SH megacities to the Antarctic region is much smaller than from NH megacities to the Arctic, and even export efficiencies relative to the emission strength are lower.In fact, due to the generally longer distances of the SH megacities to the Antarctic region, the fraction of BCtr mass transported to the lower layer of the Antarctic atmosphere (0-1 km) reaches maximum values of 2 × 10 −2 (Fig. 11c) in the last age class.Values for ages of less than about 20 days are lower and there is no peak in the early age class, again due to the considerable distance of the sources from the polar region.Moreover, these values are reached only for Buenos Aires and Sydney, while the other megacities, and especially Lima and Jakarta, give substantially lower contributions.For the depositing tracer BCdp, no more than a fraction of 10 −4 of the emitted mass in any age class reaches the Antarctic atmosphere (Fig. 11d).Summary numbers, cumulating the contribution from all age classes, are reported in Table 3 for both BCtr and BCdp.They show that for BCtr the total megacity emissions in the SH generate a suspended mass in the Antarctic atmosphere comparable to that produced by the single city of St. Petersburg in the Arctic.If deposition is furthermore taken into account, the suspended mass is entirely negligible with less than 100 kg of suspended mass overall.
The fraction of BCdp mass depositing south of 66 • S is also very low.Our simulations estimate a total annual deposition of BCdp in the Antarctic by all the megacities of about 45 tons (see Table 3), and most of it comes from just four megacities, Sydney, Johannesburg, Santiago and Buenos Aires, with a marginal contribution from Sao Paulo and Rio de Janeiro.The comparison with the total deposition of 2330 tons of BCdp from NH megacities in the Arctic shows that this deposition is very small, also considering that the total emission from NH megacities is just about 5.5 times the SH megacity emissions.However, this 45 tons of BCdp deposition is comparable to the local emissions, below 60 • S, of BC from shipping (100 tons) as reported by Stohl and Sodemann (2010) based on the emission inventory of Bond et al. (2007), with tourist ships contributing about 20 tons (Stohl and Sodemann, 2010).Also, if we consider the 12 % megacity emissions representative of all SH emissions south of 25 • S, we can estimate that less than 400 tons of BCdp are deposited in the South Pole region from all continental sources.

Exposure of human population
In this section, population (P ) and annual mean ground level concentration (GLC) maps are combined to derive simple but quantitative measures of the exposure of the human population to air pollution from megacities.Particular focus will be on the relative importance of population exposure inside and outside the megacity.We use the GLC multiplied by population (CmP) as a measure of total population exposure to air pollution and define this metric as a function of time (i.e., for any age class), (5) Here i, j are grid indices, indicating that the metric is integrated over the entire globe, and it is the index of the age class considered.Therefore, CmP is a function of time since emission, and the concentration used is given by the summation of the concentration associated with the current and the previous age classes.Figure 14 shows the values for all the age classes, and it can be seen that when deposition is considered the CmP exposure metric for BCdp (Fig. 14, right) is quite flat and determined mainly by the exposure on the very first day.This is due to a combination of three factors: (i) the highest population density at the location of the release (i.e., within the megacity), (ii) all BCdp being initially located near the surface, and (iii) the quick loss of BCdp mass with time.This suggests that most of the population exposure to BCdp occurs within or very close to the megacity.Also for the conserved tracer BCtr, CmP values do not grow substantially after the first day (Fig. 14, left), suggesting that even for longlived pollutants most of the population exposure occurs in the city or its vicinity.This is also the reason why the CmP values for BCtr and BCdp are quite similar for most cities, suggesting that population exposure is not very sensitive to a substance's lifetime if it is longer than a day.
In Table 1 we also report the population-weighted concentration for the city (PWC city ) for BCtr and BCdp: i.e., the same as the CmP metric but computed only for the grid cells defining each megacity, with GLC integrated over all the age classes (N a = 10) and normalized with the megacity's total population (P city ).This metric measures the concentration to which the population inside the city is exposed on average.Some interesting features arise when comparing the megacities.A further parameter of interest is CmP normalized by CmP obtained when integrating only over the megacity grid cells, called here CmP c .This metric (CmP / CmP c ) measures the total population exposure relative to the exposure inside the city.These values are plotted as a function of the age in Fig. 15.First, we note that for the nondepositing tracer BCtr, several cities have an impact on surrounding regions which is at least half as large as on the population of the city itself (i.e., CmP / CmP c > 1.5), but few cities have an impact outside the city's boundaries larger than the impact inside the city boundary (i.e., CmP / CmP c > 2) and only three cities reach this thresholds if deposition is taken into account (Fig. 15, right panel).For BCtr, Saint Petersburg shows the largest CmP / CmP c values but only for ages greater than 10 days.This is due both to the fact that it is surrounded, but not closely surrounded, by relatively densely populated areas and that it has the smallest population among the cities considered (thus reducing the impact inside the city boundaries).Southern Hemisphere megacities show generally somewhat lower CmP / CmP c values, with minima obtained for Buenos Aires and Sao Paulo, since these are densely populated areas surrounded by relatively sparsely populated areas and are close to the ocean.In the NH, Tokyo stands out with the lowest value of the ratio, since most of its pollution is exported to the ocean.If deposition is accounted for (BCdp), Sydney also has a very low CmP / CmP c ratio, while the cities showing highest values of CmP / CmP c are the densely populated Chinese megacities of Tianjin and Beijing, which are actually surrounded by regions with extremely high population density outside the city boundaries, and the city of Karachi.This last city shows high impact only after about 5 days; thus this should be due to the Karachi's plume impacting densely populated Indian urban areas.It is interesting to note that the cities with a high CmP / CmP c ratio also have a quite high index of low-level BCdp export (ELR1km, Fig. 8b).
In general and especially for the BCdp tracer, it seems that the impact of megacities on population outside their boundaries relative to their impact inside city boundaries is low, with only three of the megacities (Tianjin, Beijing and Karachi) showing an impact comparable to or greater than the impact inside the cities, while most of the cities have CmP / CmP c < 1.4.It is also of some interest to note that based on the time related behavior of this index one could speculate whether the impact on population is transnational.For example when deposition is involved, the Chinese megacities of Tianjin and Beijing have an index CmP / CmP c > 1.6 already in the first age class (less than a day) and this value does not grow substantially with age.Therefore, most of their influence is on national population.Conversely, Karachi reaches a similar index only in the second age class (3 days) so it is probably related to transboundary pollution transport.

Summary and conclusions 15
A study of the characteristics of air pollution plumes dispersing from a group of thirty-six of 16 world's megacities has been performed using the Lagrangian particle model FLEXPART.The 17 study focuses on plumes generated by BC emissions.BC was chosen since: i) this is a species 18 Based on the age class spectrum it is also easily possible to extrapolate the behavior for decaying scalars, but in the case of BCdp any superimposed decay with characteristic time greater than about ten days would basically leave the results unchanged.

Summary and conclusions
A study of the characteristics of air pollution plumes dispersing from a group of 36 of the world's megacities has been performed using the Lagrangian particle model FLEX-PART.The study focuses on plumes generated by BC emissions.BC was chosen since (i) this is a species well suited for a Lagrangian investigation because it can be represented in a linear modeling framework, (ii) it is representative of combustion-related emissions and (iii) it is an important air pollutant and short-lived climate forcer.Based on BC emissions two artificial tracers were modeled, a purely conserved tracer with a cut-off lifetime at 48 days -BCtr -and one resembling more closely the properties of true BC and subject to wet and dry deposition -BCdp.Four aspects of megacities' plumes have been investigated: (1) physical dispersion characteristics, e.g., transport distance and average altitude; (2) pollution export in the lower troposphere over more than one thousand kilometers (ELR1km); (3) concentrations and deposition of BC in the polar regions; and (4) exposure of human population with a focus on comparing the inside-thecity versus outside-the-city exposure values.
Although local conditions can be very important, the main factor influencing the physical dispersion characteristics was found to be the latitude of the emitting cities; initially, plumes from higher-latitude cities travel faster, but after about three days the higher elevation reached by the plumes at tropical/subtropical latitude generated a complementary faster horizontal advection.The proximity to the ocean was found to be an important factor in determining a fast early phase of plume transport.Among the European megacities the Po Valley shows a unique behavior since the tracer is trapped there by mountains on three sides, which strongly limits the early phase of transport, but, overall, in the NH the megacity with the slowest horizontal transport up to 10 days is the tropical city of Lagos.However, Lima is the city showing the slowest horizontal transport.
It was found that the cities showing a high level of regional pollution export in the lower troposphere (ELR1km) are generally high-latitude cities since these are characterized by an initially fast transport and by a low plume elevation.The influence of deposition was mainly to reduce the value of ELR1km, and this reduction was found to be stronger for cities with low values of the local characteristic deposition timescale T li .However, a few exceptions (Indian subcontinent megacities and Cairo) were observed where the value of ELR1km increased due to the effect of deposition.This behavior could be explained with the possibility of lowlevel transport at a distance larger than 1000 km followed by a preferential deposition in uplifting air masses.A good agreement between our simulations and previous results by Lawrence et al. (2007) for ELR1km was observed for a nondepositing scalar with a 10-day decaying timescale.
The contribution to lower-tropospheric concentrations and deposition in the Arctic must be mainly attributed to northern European megacities (St.Petersburg, Moscow, London and Ruhr Valley).Based on the emission inventory used, the sole city of St. Petersburg contributes up to 10 % of the BC emission north of 60 • N during winter, and when deposition is considered (BCdp) its contribution to suspended mass of BCdp in the lower 1km of the Arctic atmosphere and deposition on the Arctic surface is larger than that of the whole group of Asian megacities (population about 500 million).When considering the upper layer of the Arctic troposphere (2-9 km) and stratosphere (9-30 km), the contribution of remote sources becomes progressively more important.In the stratosphere the contribution from Asian megacities is the largest, followed by those from North American megacities and European megacities, with the Po Valley standing out among this latter group.However, all the megacities can be considered significant for pollution in this layer of Arctic atmosphere.
The export of BCtr and BCdp from Southern Hemisphere megacities to the Antarctic region is much smaller than from Northern Hemisphere megacities to the Arctic, and even the export efficiency relative to the emission strength is lower.Our simulations estimate a total annual deposition of BCdp in the Antarctic by all the Southern Hemisphere megacities of about 45 tons, compared to the total deposition of 2330 tons of BCdp from Northern Hemisphere megacities in the Arctic.Note that the contribution by all megacities to the BCdp deposition in the Antarctic is, therefore, estimated to be smaller than the BC emissions coming from the local shipping activities.
We remind the reader, however, that emissions from individual cities can be very different between different emission inventories, and this uncertainty is implicit in our analysis.
Regarding the exposure of the human population, we found that for the passive tracer BCtr ten megacities have an impact outside their boundaries greater than their impact inside the cities' boundaries.If BCdp is considered, only three megacities have larger impacts outside than inside their boundaries.These three cities are Beijing, Tianjin and Karachi, and they were shown to have quite high values of low-level pollution export -ELR1km.Therefore the combined conditions of high ELR1km and high population densities outside cities' boundaries were necessary to have a higher impact on population outside the city than that on population inside the city.In general, however, it was found that most of the impact on population is related to pollution inside the city boundaries.

Figure 3 .
Figure 3.Total mass of of BCdp normalized with that of BCtr (BCdp/BCtr) as function of the 3 age for the northern hemisphere (NH, a) and the southern hemisphere (SH, b).Lower panels (c, 4 d) as above but using a linear scale.5

Fig. 3 .
Fig. 3. Total mass of BCdp normalized with that of BCtr (BCdp/BCtr) as function of the age for the Northern Hemisphere (NH, a) and the Southern Hemisphere (SH, b).Lower panels (c, d) as above but using a linear scale.

Figure 4 .Fig. 4 .
Figure 4. Transport distance ( ഥ ሺሻሻ from the source megacity location for the vertically integrated concentration of BCtr and BCdp for the NH) (a, b) and SH (c, d).The smaller panels on the left are a zoom for 0-3 days using a logarithmic scale.

1Figure 5 . 4 Figure 6 .
Figure 5. Altitude (above sea level) of the center of mass of horizontally integrated concentration 2 ( ഥ ሺሻሻ for BCtr (a, c) and BCdp (b, d) in the NH (a, b) and in the SH (c, d). 3

Fig. 5 .
Fig. 5. Altitude (above sea level) of the center of mass of horizontally integrated concentration ( h(t)) for BCtr (a, c) and BCdp (b, d) in the NH (a, b) and in the SH (c, d).
above sea level) of the center of mass of horizontally integrated concentration c) and BCdp (b, d) in the NH (a, b) and in the SH (c, d).

Fig. 6 .
Fig. 6.Ratio of the averaged (over all cities) horizontal distance of centers of mass, D(t), and vertical position of center of mass, H (t), for BCdp over BCtr.
N cities i=1 di (t)/N cities ) and of the average megacities plume altitude (H (t) = N cities i=1 hi (t)/N cities ) for all cities.

1Figure 7 .
Figure 7. Mass in the upper troposphere and stratosphere (9-30 km) normalized by the total mass 2 suspended in the atmosphere as a function of tracer age for BCtr and BCdp, for the NH (a, b) 3 and SH (c, d). 4 5 4. Statistical measures of regional pollution export 6

Fig. 7 . 19 Figure 8 .Fig. 8 .Fig. 9 .
Fig. 7. Mass in the upper troposphere and stratosphere (9-30 km) normalized by the total mass suspended in the atmosphere as a function of tracer age for BCtr and BCdp, for the NH (a, b) and SH (c, d).

4 5Figure 10 .Fig. 10 .
Figure 10.Fraction of mass transported over 1000km away from the source over the whole air 6 column (ELRcol) for BCtr and for BCdp.The normalization is the same used before in Figure 8. 7

Fig. 11 .
Figure 11.BCtr (left panels) and BCdp (right panels) mass content in a layer 0-1000 m ATL in 2 the Arctic (a, b) and Antarctic (c, d) regions, scaled for the total mass emitted in any the 3 respective age class.4

1Figure 12 .Fig. 12 .
Figure 12.BCtr (left panels) and BCdp (right panels) mass content in the Arctic region scaled for 2 the total mass emitted in the respective age class for layers of the atmosphere extending from 2 3 to 9 km (a, b) and from 9 to 30 km (c, d) above sea level.The insert in panel (c) shows the same 4 information but on a linear scale.5 6 7 8

7 Figure 13 .Fig. 13 .
Figure 13.Deposition of BC at latitudes north of 66 o N (a, b) and south of 66 o S (c, d), relative to 8 the mass available for deposition (left) and absolute value (right).9 10 lower and there is no peak in the early age class, again due to the considerable distance of the 11 sources from the polar region.Moreover, these values are reached only for Buenos Aires and 12

Fig. 14 .
Fig. 14.The CmP metric for BCtr (left) and BCdp (right) as a function of tracer age.

Fig. 15 .
Fig. 15.CmP / CmP c as a function of tracer age, for BCtr (left) and BCdp (right).In the left panel, CmP / CmP c values for the city of St. Petersburg exceed the scale and reach about 4.5 at 48 days.
Location and extension of the megacities investigated.is calculated starting from the North and the South Pole, respectively.No data are shown south of 35 • S and north of 60 • N since there are no megacities.Further, the data are separated into autumn/winter emissions and spring/summer emissions.It can be seen that the 35

Table 2 .
Comparison of ELR1km for a conserved tracer with no decay (BCtr) and one subject to deposition (BCdp), for the BCtr tracer with a superimposed exponential decay process with a timescale of 10 days (BC10, see text for clarification) and a tracer with a 10-day decay timescale computed by L07 using the model MATCH for the year 1995.

Table 3 .
The Table The values computed for the 48-day age class are reported in Table1for both BCtr and BCdp.Notice that since the lifetime of BCdp is substantially shorter than 48 days for all cities, the metric measures the total impact of BCdp.The table shows that the largest values of CmP are reached for the megacities in Asia, with Dhaka setting the records of 149 and 104 (10 10 ng m −3 people), respectively, for BCtr and BCdp.The first non-Asian megacity is the extended New York area which ranks seventh, with values of 40 and 22 for BCtr and BCdp, while the first-ranking European city is Moscow with CmP values of 19 and 14. very low CmP/CmPc ratio while the cities showing highest values of CmP/CmPc are the densely 5 populated Chinese megacities of Tianjin and Beijing, which are actually surrounded by regions 6 with extremely high population density outside the city boundaries, and the city of Karachi.This 7 last city shows high impact only after about 5 days, thus this should be due to the Karachi's 8 plume impacting on densely populated Indian urban areas.It is interesting to note that the cities 9 with high CmP/CmPc ratio also have a quite high index of low level BCdp export (ELR1km, The CmP metric for BCtr (left) and BCdp (right) as a function of tracer age.13 14In general and especially for the BCdp tracer, it seems that the impact of megacities on 15 population outside their boundaries relative to their impact inside city boundaries is low, with 16 only three of the megacities (Tianjin, Beijing and Karachi) showing an impact comparable or 17 greater than the impact inside the cities, while most of the cities have CmP/CmPc < 1.4.It is 18 The highest value of PWC can be found in Asian megacities, with Seoul having values of 1621 and 1535 ng m −3 , respectively, for BCtr and BCdp.The lowest values are registered for Istanbul, respectively 64 and 61 ng m −3 for BCtr and BCdp.It should be noted that for cities covering only one cell this value is the same as the maximum concentration in the domain, while it can be significantly lower for larger megacities.
Therefore, most of their influence is on national population.On the contrary, Karachi reaches a 5 similar index only in the second age class (3 days) so it is probably related to transboundary 6 pollution transport.7 Based on the age class spectrum it is also easily possible to extrapolate the behavior for decaying 8 scalars but in case of BCdp any superimposed decay with characteristic time greater than about 9 ten days would basically leave the results unchanged.10 11 Figure 15.CmP/CmPc as a function of tracer age, for BCtr (left) and BCdp (right).In the left 12 panel, CmP/CmPc values for the city of St. Petersburg exceed the scale and reach about 4.5 at 48 13 days.14