Search

There are growing concerns about the climate impacts of absorbing organic carbon (also known as brown carbon, BrC) in the environment, yet its chemical composition and association with the light absorption capabilities remain poorly understood. This study characterized water-soluble and water-insoluble organic carbon (WSOC and WIOC) from residential solid fuel combustion at the molecular level and evaluated their quantitative relationship with mass absorption efficiency (MAE). The MAE values at λ = 365 nm from biomass burning were significantly higher than those from coal combustion (p < 0.05). Thousands of peaks were identified in the m/z range of 150–800, with the most intense ion peaks occurring between m/z 200–500 for WSOC and m/z 600–800 for WIOC, respectively. The CHO group predominated in the WSOC extract from biomass burning emissions, while sulfur-containing compounds (SOCs) including CHOS and CHONS were more intense in the WIOC extract, particularly from coal emissions. Emissions of the CHON group were positively correlated with the fuel nitrogen content (r = 0.936; p < 0.05), explaining their higher abundance in coal emissions compared to biomass. The SOC emissions were more predominant during flaming phases, as indicated by a positive correlation with modified combustion efficiency (MCE) (r = 0.750; p < 0.05). The unique formulas of coal combustion aerosols were in the lower H/C and O/C regions, with higher unsaturated compounds in the van Krevelen (VK) diagram. In the WIOC extract, coal combustion emissions contained significantly higher fractions of condensed aromatics (32 %–59 %) compared to only 4.3 %–9.7 % in biomass burning emissions. In contrast, the CHOS group in biomass burning emissions was characterized by larger condensed aromatic compound fractions than those in coal combustion. Moreover, the CHOS aromatic compound fractions were positively correlated with MAE values in both WSOC (r = 0.714; p < 0.05) and WIOC extracts (r = 0.929; p < 0.001), suggesting that these compounds significantly contributed to MAE variabilities across different fuels.

Associations between tropospheric ozone (O3) and climate variations have been extensively investigated worldwide. However, given the lack of historical O3 monitoring data, the knowledge gaps regarding the influences of climate variations on long-term O3 trends in China remain. The present study used a tropospheric O3 dataset from the summers of 1999 to 2017 simulated by an atmospheric chemistry model to explore the linkage between summer O3 and a dominant atmospheric circulation system – the Western Pacific Subtropical High (WPSH) pressure – on an interannual basis in China. During this period, both WPSH strength and O3 concentrations in eastern and central China illustrated a growing trend. An EOF analysis was conducted to examine significant summer O3 characteristics and patterns and their potential connections with the WPSH. We find that the correlation between the first principal component of summer ozone concentration in the EOF analysis and the WPSH reached 0.56 (P≤0.01) in China from 1999 to 2017. We show that the WPSH determines interannual fluctuations of summer O3, whereas O3 precursor emissions contribute primarily to the O3 long-term trend. Our results reveal that the WPSH plays a vital role in O3 perturbation in the eastern seaboard regions and inland China. Precursor emissions made more significant contributions of up to 60 % to increasing O3 trends in the inland urban agglomerations than coastal regions in eastern and southern China. The strongest contribution of meteorological conditions associated with the WPSH to summer O3 occurred in the Yangtze River Delta (YRD), accounting for over 9 % to ozone perturbations from 1999 to 2017. We find that the effect of the WPSH on regional O3 depends on the spatial proximity to the WPSH. We attributed the effects of the WPSH on O3 interannual variations to the changes in air temperature, precipitation, and winds associated with the WPSH's intensity and positions.

Reactive organic gases (ROGs), as important precursors of secondary pollutants, are not well resolved as their chemical complexity has challenged their quantification in many studies. Here, a near-complete speciation of ROG emissions from residential combustion was developed by the combination of proton transfer reaction time-of-flight mass spectrometry (PTR-ToF-MS) with a gas chromatography system equipped with a mass spectrometer and a flame ionization detector (GC-MS/FID), including 1049 species in all. Among them, 125 identified species, ∼ 90 % of the total ROG masses, were applied to evaluate their emission characteristics through real combustion samplings in rural households of China. The study revealed that with 55 species, mainly oxygenated species, higher hydrocarbons with ≥8 carbon atoms, and nitrogen-containing species, previously un- and under-characterized, ROG emissions from residential coal and biomass combustion were underestimated by 44.3 % ± 11.8 % and 22.7 % ± 3.9 %, respectively, which further amplified the underestimation of secondary organic aerosol formation potential (SOAP) as high as 70.3 % ± 1.6 % and 89.2 % ± 1.0 %, respectively. The hydroxyl radical reactivity (OHR) of ROG emissions was also undervalued significantly. The study provided a feasible method for the near-complete speciation of ROGs in the atmosphere and highlighted the importance of acquiring completely speciated measurement of ROGs from residential emissions, as well as other processes.

Polychlorinated biphenyls (PCBs) are nowadays mainly unintentionally produced, since their use and manufacture has been banned. PCBs, especially those dioxin-like compounds, have been proven harmful to ecosystems and human health. Detailed information about the spatiotemporally distribution of unintentionally produced PCBs (UP-PCBs) is crucial for understanding the environmental fate and associated health risks. However, research estimating UP-PCBs emission in China is limited. In this study, a source-specific and year-varying emission factor (EF) dataset was used to compile provincial emissions of 12 dioxin-like UP-PCBs congeners for 66 sources in the Chinese mainland from 1960 to 2019. In addition to that, emissions for 2019 were allocated into 0.1∘×0.1∘ grids. The results indicated that historical national UP-PCBs emission represented an increasing trend until around 1995 and then showed an overall decreasing trend from 1995 to 2019. Cement production was the largest UP-PCBs emission source in the Chinese mainland. Geographically, East and North China contributed the largest UP-PCBs emission across the Chinese mainland. High emission densities occurred mainly in the densely populated and well-developed regions such as Yangtze River Delta Urban Agglomeration and Pearl River Delta. Furthermore, highly positive correlations of emission densities with population and GDP densities were identified.

Aerosol acidity plays a key role in regulating the chemistry and toxicity of atmospheric aerosol particles. The trend of aerosol pH and its drivers is crucial in understanding the multiphase formation pathways of aerosols. Here, we reported the first trend analysis of aerosol pH from 2011 to 2019 in eastern China, calculated with the ISORROPIA model based on observed gas and aerosol compositions. The implementation of the Air Pollution Prevention and Control Action Plan led to −35.8 %, −37.6 %, −9.6 %, −81.0 % and 1.2 % changes of PM2.5, SO42-, NHx, non-volatile cations (NVCs) and NO3- in the Yangtze River Delta (YRD) region during this period. Different from the drastic changes of aerosol compositions due to the implementation of the Air Pollution Prevention and Control Action Plan, aerosol pH showed a minor change of −0.24 over the 9 years. Besides the multiphase buffer effect, the opposite effects from the changes of SO42- and non-volatile cations played key roles in determining this minor pH trend, contributing to a change of +0.38 and −0.35, respectively. Seasonal variations in aerosol pH were mainly driven by the temperature, while the diurnal variations were driven by both temperature and relative humidity. In the future, SO2, NOx and NH3 emissions are expected to be further reduced by 86.9 %, 74.9 % and 41.7 % in 2050 according to the best health effect pollution control scenario (SSP1-26-BHE). The corresponding aerosol pH in eastern China is estimated to increase by ∼0.19, resulting in 0.04 less NO3- and 0.12 less NH4+ partitioning ratios, which suggests that NH3 and NOx emission controls are effective in mitigating haze pollution in eastern China.

Sulfur dioxide (SO2) is a major atmospheric pollutant and precursor of sulfate aerosols, which influences air quality, cloud microphysics, and climate. Therefore, better understanding the conversion of SO2 to sulfate is essential to simulate and predict sulfur compounds more accurately. This study evaluates the effects of in-cloud aqueous-phase chemistry on SO2 oxidation in the Community Earth System Model version 2 (CESM2). We replaced the default parameterized SO2 aqueous-phase reactions with detailed HOx, Fe, N, and carbonate chemistry in cloud droplets and performed a global simulation for 2014–2015. Compared with the observations, the results incorporating detailed cloud aqueous-phase chemistry greatly reduced SO2 overestimation. This overestimation was reduced by 0.1–10 ppbv (parts per billion by volume) in most of Europe, North America, and Asia and more than 10 ppbv in parts of China. The biases in annual simulated SO2 mixing ratios decreased by 46 %, 41 %, and 22 % in Europe, the USA, and China, respectively. Fe chemistry and HOx chemistry contributed more to SO2 oxidation than N chemistry. Higher concentrations of soluble Fe and higher pH values could further enhance the oxidation capacity. This study emphasizes the importance of detailed in-cloud aqueous-phase chemistry for the oxidation of SO2. These mechanisms can improve SO2 simulation in CESM2 and deepen understanding of SO2 oxidation and sulfate formation.

China's fossil-fuel CO2 (FFCO2) emissions accounted for approximately 28 % of the global total FFCO2 in 2016. An accurate estimate of China's FFCO2 emissions is a prerequisite for global and regional carbon budget analyses and the monitoring of carbon emission reduction efforts. However, significant uncertainties and discrepancies exist in estimations of China's FFCO2 emissions due to a lack of detailed traceable emission factors (EFs) and multiple statistical data sources. Here, we evaluated China's FFCO2 emissions from nine published global and regional emission datasets. These datasets show that the total emissions increased from 3.4 (3.0–3.7) in 2000 to 9.8 (9.2–10.4) Gt CO2 yr−1 in 2016. The variations in these estimates were largely due to the different EF (0.491–0.746 t C per t of coal) and activity data. The large-scale patterns of gridded emissions showed a reasonable agreement, with high emissions being concentrated in major city clusters, and the standard deviation mostly ranged from 10 % to 40 % at the provincial level. However, patterns beyond the provincial scale varied significantly, with the top 5 % of the grid level accounting for 50 %–90 % of total emissions in these datasets. Our findings highlight the significance of using locally measured EF for Chinese coal. To reduce uncertainty, we recommend using physical CO2 measurements and use these values for dataset validation, key input data sharing (e.g., point sources), and finer-resolution validations at various levels.

The hygroscopic behavior of black carbon (BC)-containing particles (BCPs) has a significant impact on global and regional climate change. However, the mechanism and factors controlling the hygroscopicity of BCPs from different carbon sources are not well understood. Here, we systematically measured the equilibrium and kinetics of water uptake by 15 different BCPs (10 herb-derived BCPs, 2 wood-derived BCPs, and 3 soot-type BCPs) using a gravimetric water vapor sorption method combined with in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS). In the gravimetric analysis, the sorption–desorption equilibrium isotherms were measured under continuous-stepwise water vapor pressure conditions, while the kinetics was measured at a variety of humidity levels obtained by different saturated aqueous salt solutions. The equilibrium water uptake of the tested group of BCPs at high relative humidity (>80 %) positively correlated to the dissolved mineral content (0.01–13.0 wt %) (R2=0.86, P=0.0001), the content of the thermogravimetrically analyzed organic carbon (OCTGA, 4.48–15.25 wt %) (R2=0.52, P=0.002), and the content of the alkali-extracted organic carbon (OCAE, 0.14–8.39 wt %) (R2=0.80, P=0.0001). In contrast, no positive correlation was obtained with the content of total organic carbon or elemental carbon. Among the major soluble ionic constituents, chloride and ammonium were each correlated with the equilibrium water uptake at high relative humidity. Compared with the herbal BCPs and soot, the woody BCPs had much lower equilibrium water uptake, especially at high relative humidity, likely due to the very low dissolved mineral content and OC content. The DRIFTS analysis provided generally consistent results at low relative humidity. The kinetics of water uptake (measured by pseudo-second-order rate constant) correlated to the content of OCTGA and OCAE as well as the content of chloride and ammonium at low relative humidity (33 %) but to the porosity of BCPs at high relative humidity (94 %). This was the first study to show that BCPs of different types and sources had greatly varying hygroscopic properties.

The Atmospheric Pollution and Human Health in a Chinese Megacity (APHH-Beijing) programme is an international collaborative project focusing on understanding the sources, processes and health effects of air pollution in the Beijing megacity. APHH-Beijing brings together leading China and UK research groups, state-of-the-art infrastructure and air quality models to work on four research themes: (1) sources and emissions of air pollutants; (2) atmospheric processes affecting urban air pollution; (3) air pollution exposure and health impacts; and (4) interventions and solutions. Themes 1 and 2 are closely integrated and support Theme 3, while Themes 1–3 provide scientific data for Theme 4 to develop cost-effective air pollution mitigation solutions. This paper provides an introduction to (i) the rationale of the APHH-Beijing programme and (ii) the measurement and modelling activities performed as part of it. In addition, this paper introduces the meteorology and air quality conditions during two joint intensive field campaigns – a core integration activity in APHH-Beijing. The coordinated campaigns provided observations of the atmospheric chemistry and physics at two sites: (i) the Institute of Atmospheric Physics in central Beijing and (ii) Pinggu in rural Beijing during 10 November–10 December 2016 (winter) and 21 May–22 June 2017 (summer). The campaigns were complemented by numerical modelling and automatic air quality and low-cost sensor observations in the Beijing megacity. In summary, the paper provides background information on the APHH-Beijing programme and sets the scene for more focused papers addressing specific aspects, processes and effects of air pollution in Beijing.

Parameterizations that impact wet removal of black carbon (BC) remain uncertain in global climate models. In this study, we enhance the default wet deposition scheme for BC in the Community Earth System Model (CESM) to (a) add relevant physical processes that were not resolved in the default model and (b) facilitate understanding of the relative importance of various cloud processes on BC distributions. We find that the enhanced scheme greatly improves model performance against HIPPO observations relative to the default scheme. We find that convection scavenging, aerosol activation, ice nucleation, evaporation of rain or snow, and below-cloud scavenging dominate wet deposition of BC. BC conversion rates for processes related to in-cloud water–ice conversion (i.e., riming, the Bergeron process, and evaporation of cloud water sedimentation) are relatively smaller, but have large seasonal variations. We also conduct sensitivity simulations that turn off each cloud process one at a time to quantify the influence of cloud processes on BC distributions and radiative forcing. Convective scavenging is found to have the largest impact on BC concentrations at mid-altitudes over the tropics and even globally. In addition, BC is sensitive to all cloud processes over the Northern Hemisphere at high latitudes. As for BC vertical distributions, convective scavenging greatly influences BC fractions at different altitudes. Suppressing BC droplet activation in clouds mainly decreases the fraction of column BC below 5 km, whereas suppressing BC ice nucleation increases that above 10 km. During wintertime, the Bergeron process also significantly increases BC concentrations at lower altitudes over the Arctic. Our simulation yields a global BC burden of 85 Gg; corresponding direct radiative forcing (DRF) of BC estimated using the Parallel Offline Radiative Transfer (PORT) is 0.13 W m−2, much lower than previous studies. The range of DRF derived from sensitivity simulations is large, 0.09–0.33 W m−2, corresponding to BC burdens varying from 73 to 151 Gg. Due to differences in BC vertical distributions among each sensitivity simulation, fractional changes in DRF (relative to the baseline simulation) are always higher than fractional changes in BC burdens; this occurs because relocating BC in the vertical influences the radiative forcing per BC mass. Our results highlight the influences of cloud microphysical processes on BC concentrations and radiative forcing.

Combining measurements of atmospheric CO2 and its radiocarbon (14CO2) fraction and transport modeling in atmospheric inversions offers a way to derive improved estimates of CO2 emitted from fossil fuel (FFCO2). In this study, we solve for the monthly FFCO2 emission budgets at regional scale (i.e., the size of a medium-sized country in Europe) and investigate the performance of different observation networks and sampling strategies across Europe. The inversion system is built on the LMDZv4 global transport model at 3.75° × 2.5° resolution. We conduct Observing System Simulation Experiments (OSSEs) and use two types of diagnostics to assess the potential of the observation and inverse modeling frameworks. The first one relies on the theoretical computation of the uncertainty in the estimate of emissions from the inversion, known as posterior uncertainty, and on the uncertainty reduction compared to the uncertainty in the inventories of these emissions, which are used as a prior knowledge by the inversion (called prior uncertainty). The second one is based on comparisons of prior and posterior estimates of the emission to synthetic true emissions when these true emissions are used beforehand to generate the synthetic fossil fuel CO2 mixing ratio measurements that are assimilated in the inversion. With 17 stations currently measuring 14CO2 across Europe using 2-week integrated sampling, the uncertainty reduction for monthly FFCO2 emissions in a country where the network is rather dense like Germany, is larger than 30 %. With the 43 14CO2 measurement stations planned in Europe, the uncertainty reduction for monthly FFCO2 emissions is increased for the UK, France, Italy, eastern Europe and the Balkans, depending on the configuration of prior uncertainty. Further increasing the number of stations or the sampling frequency improves the uncertainty reduction (up to 40 to 70 %) in high emitting regions, but the performance of the inversion remains limited over low-emitting regions, even assuming a dense observation network covering the whole of Europe. This study also shows that both the theoretical uncertainty reduction (and resulting posterior uncertainty) from the inversion and the posterior estimate of emissions itself, for a given prior and true estimate of the emissions, are highly sensitive to the choice between two configurations of the prior uncertainty derived from the general estimate by inventory compilers or computations on existing inventories. In particular, when the configuration of the prior uncertainty statistics in the inversion system does not match the difference between these prior and true estimates, the posterior estimate of emissions deviates significantly from the truth. This highlights the difficulty of filtering the targeted signal in the model–data misfit for this specific inversion framework, the need to strongly rely on the prior uncertainty characterization for this and, consequently, the need for improved estimates of the uncertainties in current emission inventories for real applications with actual data. We apply the posterior uncertainty in annual emissions to the problem of detecting a trend of FFCO2, showing that increasing the monitoring period (e.g., more than 20 years) is more efficient than reducing uncertainty in annual emissions by adding stations. The coarse spatial resolution of the atmospheric transport model used in this OSSE (typical of models used for global inversions of natural CO2 fluxes) leads to large representation errors (related to the inability of the transport model to capture the spatial variability of the actual fluxes and mixing ratios at subgrid scales), which is a key limitation of our OSSE setup to improve the accuracy of the monitoring of FFCO2 emissions in European regions. Using a high-resolution transport model should improve the potential to retrieve FFCO2 emissions, and this needs to be investigated.

The response of surface ozone (O3) concentrations to basin-scale warming and cooling of Northern Hemisphere oceans is investigated using the Community Earth System Model (CESM). Idealized, spatially uniform sea surface temperature (SST) anomalies of ±1 °C are individually superimposed onto the North Pacific, North Atlantic, and North Indian oceans. Our simulations suggest large seasonal and regional variability in surface O3 in response to SST anomalies, especially in the boreal summer. The responses of surface O3 associated with basin-scale SST warming and cooling have similar magnitude but are opposite in sign. Increasing the SST by 1 °C in one of the oceans generally decreases the surface O3 concentrations from 1 to 5 ppbv. With fixed emissions, SST increases in a specific ocean basin in the Northern Hemisphere tend to increase the summertime surface O3 concentrations over upwind regions, accompanied by a widespread reduction over downwind continents. We implement the integrated process rate (IPR) analysis in CESM and find that meteorological O3 transport in response to SST changes is the key process causing surface O3 perturbations in most cases. During the boreal summer, basin-scale SST warming facilitates the vertical transport of O3 to the surface over upwind regions while significantly reducing the vertical transport over downwind continents. This process, as confirmed by tagged CO-like tracers, indicates a considerable suppression of intercontinental O3 transport due to increased tropospheric stability at lower midlatitudes induced by SST changes. Conversely, the responses of chemical O3 production to regional SST warming can exert positive effects on surface O3 levels over highly polluted continents, except South Asia, where intensified cloud loading in response to North Indian SST warming depresses both the surface air temperature and solar radiation, and thus photochemical O3 production. Our findings indicate a robust linkage between basin-scale SST variability and continental surface O3 pollution, which should be considered in regional air quality management.

Methane (CH4) has a 28-fold greater global warming potential than CO2 over 100 years. Atmospheric CH4 concentration has tripled since 1750. Anthropogenic CH4 emissions from China have been growing rapidly in the past decades and contribute more than 10 % of global anthropogenic CH4 emissions with large uncertainties in existing global inventories, generally limited to country-scale statistics. To date, a long-term CH4 emission inventory including the major sources sectors and based on province-level emission factors is still lacking. In this study, we produced a detailed annual bottom-up inventory of anthropogenic CH4 emissions from the eight major source sectors in China for the period 1980–2010. In the past 3 decades, the total CH4 emissions increased from 24.4 [18.6–30.5] Tg CH4 yr−1 in 1980 (mean [minimum–maximum of 95 % confidence interval]) to 44.9 [36.6–56.4] Tg CH4 yr−1 in 2010. Most of this increase took place in the 2000s decade with averaged yearly emissions of 38.5 [30.6–48.3] Tg CH4 yr−1. This fast increase of the total CH4 emissions after 2000 is mainly driven by CH4 emissions from coal exploitation. The largest contribution to total CH4 emissions also shifted from rice cultivation in 1980 to coal exploitation in 2010. The total emissions inferred in this work compare well with the EPA inventory but appear to be 36 and 18 % lower than the EDGAR4.2 inventory and the estimates using the same method but IPCC default emission factors, respectively. The uncertainty of our inventory is investigated using emission factors collected from state-of-the-art published literatures. We also distributed province-scale emissions into 0.1° × 0.1° maps using socioeconomic activity data. This new inventory could help understanding CH4 budgets at regional scale and guiding CH4 mitigation policies in China.

We develop and examine a microphysics-based black carbon (BC) aerosol aging scheme that accounts for condensation, coagulation, and heterogeneous chemical oxidation processes in a global 3-D chemical transport model (GEOS-Chem) by interpreting the BC measurements from the HIAPER Pole-to-Pole Observations (HIPPO, 2009–2011) using the model. We convert aerosol mass in the model to number concentration by assuming lognormal aerosol size distributions and compute the microphysical BC aging rate (excluding chemical oxidation aging) explicitly from the condensation of soluble materials onto hydrophobic BC and the coagulation between hydrophobic BC and preexisting soluble particles. The chemical oxidation aging is tested in the sensitivity simulation. The microphysical aging rate is ∼ 4 times higher in the lower troposphere over source regions than that from a fixed aging scheme with an e-folding time of 1.2 days. The higher aging rate reflects the large emissions of sulfate–nitrate and secondary organic aerosol precursors hence faster BC aging through condensation and coagulation. In contrast, the microphysical aging is more than 5-fold slower than the fixed aging in remote regions, where condensation and coagulation are weak. Globally, BC microphysical aging is dominated by condensation, while coagulation contribution is largest over eastern China, India, and central Africa. The fixed aging scheme results in an overestimate of HIPPO BC throughout the troposphere by a factor of 6 on average. The microphysical scheme reduces this discrepancy by a factor of ∼ 3, particularly in the middle and upper troposphere. It also leads to a 3-fold reduction in model bias in the latitudinal BC column burden averaged along the HIPPO flight tracks, with largest improvements in the tropics. The resulting global annual mean BC lifetime is 4.2 days and BC burden is 0.25 mg m−2, with 7.3 % of the burden at high altitudes (above 5 km). Wet scavenging accounts for 80.3 % of global BC deposition. We find that, in source regions, the microphysical aging rate is insensitive to aerosol size distribution, condensation threshold, and chemical oxidation aging, while it is the opposite in remote regions, where the aging rate is orders of magnitude smaller. As a result, global BC burden and lifetime show little sensitivity (< 5 % change) to these three factors.

Improving the ability of global models to predict concentrations of black carbon (BC) over the Pacific Ocean is essential to evaluate the impact of BC on marine climate. In this study, we tag BC tracers from 13 source regions around the globe in a global chemical transport model, Model for Ozone and Related Chemical Tracers, version 4 (MOZART-4). Numerous sensitivity simulations are carried out varying the aging timescale of BC emitted from each source region. The aging timescale for each source region is optimized by minimizing errors in vertical profiles of BC mass mixing ratios between simulations and HIAPER Pole-to-Pole Observations (HIPPO). For most HIPPO deployments, in the Northern Hemisphere, optimized aging timescales are less than half a day for BC emitted from tropical and midlatitude source regions and about 1 week for BC emitted from high-latitude regions in all seasons except summer. We find that East Asian emissions contribute most to the BC loading over the North Pacific, while South American, African and Australian emissions dominate BC loadings over the South Pacific. Dominant source regions contributing to BC loadings in other parts of the globe are also assessed. The lifetime of BC originating from East Asia (i.e., the world's largest BC emitter) is found to be only 2.2 days, much shorter than the global average lifetime of 4.9 days, making the contribution from East Asia to the global BC burden only 36 % of that from the second largest emitter, Africa. Thus, evaluating only relative emission rates without accounting for differences in aging timescales and deposition rates is not predictive of the contribution of a given source region to climate impacts. Our simulations indicate that the lifetime of BC increases nearly linearly with aging timescale for all source regions. When the aging rate is fast, the lifetime of BC is largely determined by factors that control local deposition rates (e.g., precipitation). The sensitivity of lifetime to aging timescale depends strongly on the initial hygroscopicity of freshly emitted BC. Our findings suggest that the aging timescale of BC varies significantly by region and season and can strongly influence the contribution of source regions to BC burdens around the globe. Therefore, improving parameterizations of the aging process for BC is important for enhancing the predictive skill of global models. Future observations that investigate the evolution of the hygroscopicity of BC as it ages from different source regions to the remote atmosphere are urgently needed.

Rapid urbanization throughout eastern China is imposing an irreversible effect on local climate and air quality. In this paper, we examine the response of a range of meteorological and air quality indicators to urbanization. Our study uses the Weather Research and Forecasting model coupled with chemistry (WRF/Chem) to simulate the climate and air quality impacts of four hypothetical urbanization scenarios with fixed surface pollutant emissions during the month of July from 2008 to 2012. An improved integrated process rate (IPR) analysis scheme is implemented in WRF/Chem to investigate the mechanisms behind the forcing–response relationship at the process level. For all years, as urban land area expands, concentrations of CO, elemental carbon (EC), and particulate matter with aerodynamic diameter less than 2.5 microns (PM2.5) tend to decrease near the surface (below ~ 500 m), but increase at higher altitudes (1–3 km), resulting in a reduced vertical concentration gradient. On the other hand, the O3 burden, averaged over all newly urbanized grid cells, consistently increases from the surface to a height of about 4 km. Sensitivity tests show that the responses of pollutant concentrations to the spatial extent of urbanization are nearly linear near the surface, but nonlinear at higher altitudes. Over eastern China, each 10 % increase in nearby urban land coverage on average leads to a decrease of approximately 2 % in surface concentrations for CO, EC, and PM2.5, while for O3 an increase of about 1 % is simulated. At 800 hPa, pollutants' concentrations tend to increase even more rapidly with an increase in nearby urban land coverage. This indicates that as large tracts of new urban land emerge, the influence of urban expansion on meteorology and air pollution would be significantly amplified. IPR analysis reveals the contribution of individual atmospheric processes to pollutants' concentration changes. It indicates that, for primary pollutants, the enhanced sink (source) caused by turbulent mixing and vertical advection in the lower (upper) atmosphere could be a key factor in changes to simulated vertical profiles. The evolution of secondary pollutants is further influenced by the upward relocation of precursors that impact gas-phase chemistry for O3 and aerosol processes for PM2.5. Our study indicates that dense urbanization has a moderate dilution effect on surface primary airborne contaminants, but may intensify severe haze and ozone pollution if local emissions are not well controlled.

Atmospheric deposition of iron (Fe) plays an important role in controlling oceanic primary productivity. However, the sources of Fe in the atmosphere are not well understood. In particular, the combustion sources of Fe and the subsequent deposition to the oceans have been accounted for in only few ocean biogeochemical models of the carbon cycle. Here we used a mass-balance method to estimate the emissions of Fe from the combustion of fossil fuels and biomass by accounting for the Fe contents in fuel and the partitioning of Fe during combustion. The emissions of Fe attached to aerosols from combustion sources were estimated by particle size, and their uncertainties were quantified by a Monte Carlo simulation. The emissions of Fe from mineral sources were estimated using the latest soil mineralogical database to date. As a result, the total Fe emissions from combustion averaged for 1960–2007 were estimated to be 5.3 Tg yr−1 (90% confidence of 2.3 to 12.1). Of these emissions, 1, 27 and 72% were emitted in particles < 1 μm (PM1), 1–10 μm (PM1-10), and > 10 μm (PM> 10), respectively, compared to a total Fe emission from mineral dust of 41.0 Tg yr−1 in a log-normal distribution with a mass median diameter of 2.5 μm and a geometric standard deviation of 2. For combustion sources, different temporal trends were found in fine and medium-to-coarse particles, with a notable increase in Fe emissions in PM1 since 2000 due to an increase in Fe emission from motor vehicles (from 0.008 to 0.0103 Tg yr−1 in 2000 and 2007, respectively). These emissions have been introduced in a global 3-D transport model run at a spatial resolution of 0.94° latitude by 1.28° longitude to evaluate our estimation of Fe emissions. The modelled Fe concentrations as monthly means were compared with the monthly (57 sites) or daily (768 sites) measured concentrations at a total of 825 sampling stations. The deviation between modelled and observed Fe concentrations attached to aerosols at the surface was within a factor of 2 at most sampling stations, and the deviation was within a factor of 1.5 at sampling stations dominated by combustion sources. We analysed the relative contribution of combustion sources to total Fe concentrations over different regions of the world. The new mineralogical database led to a modest improvement in the simulation relative to station data even in dust-dominated regions, but could provide useful information on the chemical forms of Fe in dust for coupling with ocean biota models. We estimated a total Fe deposition sink of 8.4 Tg yr−1 over global oceans, 7% of which originated from the combustion sources. Our central estimates of Fe emissions from fossil fuel combustion (mainly from coal) are generally higher than those in previous studies, although they are within the uncertainty range of our estimates. In particular, the higher than previously estimated Fe emission from coal combustion implies a larger atmospheric anthropogenic input of soluble Fe to the northern Atlantic and northern Pacific Oceans, which is expected to enhance the biological carbon pump in those regions.

Long-range transport of black carbon (BC) is a growing concern as a result of the efficiency of BC in warming the climate and its adverse impact on human health. We study transpacific transport of BC during HIPPO-3 using a combination of inverse modeling and sensitivity analysis. We use the GEOS-Chem chemical transport model and its adjoint to constrain Asian BC emissions and estimate the source of BC over the North Pacific. We find that different sources of BC dominate the transport to the North Pacific during the southbound (29 March 2010) and northbound (13 April 2010) measurements in HIPPO-3. While biomass burning in Southeast Asia (SE) contributes about 60% of BC in March, more than 90% of BC comes from fossil fuel and biofuel combustion in East Asia (EA) during the April mission. GEOS-Chem simulations generally resolve the spatial and temporal variation of BC concentrations over the North Pacific, but are unable to reproduce the low and high tails of the observed BC distribution. We find that the optimized BC emissions derived from inverse modeling fail to improve model simulations significantly. This failure indicates that uncertainties in BC removal as well as transport, rather than in emissions, account for the major biases in GEOS-Chem simulations of BC over the North Pacific. The aging process, transforming BC from hydrophobic into hydrophilic form, is one of the key factors controlling wet scavenging and remote concentrations of BC. Sensitivity tests on BC aging (ignoring uncertainties of other factors controlling BC long range transport) suggest that in order to fit HIPPO-3 observations, the aging timescale of anthropogenic BC from EA may be several hours (faster than assumed in most global models), while the aging process of biomass burning BC from SE may occur much slower, with a timescale of a few days. To evaluate the effects of BC aging and wet deposition on transpacific transport of BC, we develop an idealized model of BC transport. We find that the mid-latitude air masses sampled during HIPPO-3 may have experienced a series of precipitation events, particularly near the EA and SE source region. Transpacific transport of BC is sensitive to BC aging when the aging rate is fast; this sensitivity peaks when the aging timescale is in the range of 1–1.5 d. Our findings indicate that BC aging close to the source must be simulated accurately at a process level in order to simulate better the global abundance and climate forcing of BC.

Multiannual satellite measurements of tropospheric NO2 columns are used for evaluation of CO2 emission changes in China in the period from 1996 to 2008. Indirect top-down annual estimates of CO2 emissions are derived from the satellite NO2 column measurements by means of a simple inverse modeling procedure involving simulations performed with the CHIMERE mesoscale chemistry–transport model and the CO2-to-NOx emission ratios from the Emission Database for Global Atmospheric Research (EDGAR) global anthropogenic emission inventory and Regional Emission Inventory in Asia (REAS). Exponential trends in the normalized time series of annual emissions are evaluated separately for the periods from 1996 to 2001 and from 2001 to 2008. The results indicate that the both periods manifest strong positive trends in the CO2 emissions, and that the trend in the second period was significantly larger than the trend in the first period. Specifically, the trends in the first and second periods are best estimated to be in the range from 3.7 to 8.3 and from 11.0 to 13.2% per year, respectively, taking into account statistical uncertainties and differences between the CO2-to-NOx emission ratios from the EDGAR and REAS inventories. Comparison of our indirect top-down estimates of the CO2 emission changes with the corresponding bottom-up estimates provided by the EDGAR (version 4.2) and Global Carbon Project (GCP) glomal emission inventories reveals that while acceleration of the CO2 emission growth in the considered period is a common feature of both kinds of estimates, nonlinearity in the CO2 emission changes may be strongly exaggerated in the global emission inventories. Specifically, the atmospheric NO2 observations do not confirm the existence of a sharp bend in the emission inventory data time series in the period from 2000 to 2002. A significant quantitative difference is revealed between the bottom-up and indirect top-down estimates of the CO2 emission trend in the period from 1996 to 2001 (specifically, the trend was not positive according to the global emission inventories, but is strongly positive in our estimates). These results confirm the findings of earlier studies that indicated probable large uncertainties in the energy production and other activity data for China from international energy statistics used as the input information in the global emission inventories. For the period from 2001 to 2008, some quantitative differences between the different kinds of estimates are found to be in the range of possible systematic uncertainties associated with our estimation method. In general, satellite measurements of tropospheric NO2 are shown to be a useful source of information on CO2 sources collocated with sources of nitrogen oxides; the corresponding potential of these measurements should be exploited further in future studies.

High-resolution mapping of fuel combustion and CO2 emission provides valuable information for modeling pollutant transport, developing mitigation policy, and for inverse modeling of CO2 fluxes. Previous global emission maps included only few fuel types, and emissions were estimated on a grid by distributing national fuel data on an equal per capita basis, using population density maps. This process distorts the geographical distribution of emissions within countries. In this study, a sub-national disaggregation method (SDM) of fuel data is applied to establish a global 0.1° × 0.1° geo-referenced inventory of fuel combustion (PKU-FUEL) and corresponding CO2 emissions (PKU-CO2) based upon 64 fuel sub-types for the year 2007. Uncertainties of the emission maps are evaluated using a Monte Carlo method. It is estimated that CO2 emission from combustion sources including fossil fuel, biomass, and solid wastes in 2007 was 11.2 Pg C yr−1 (9.1 Pg C yr−1 and 13.3 Pg C yr−1 as 5th and 95th percentiles). Of this, emission from fossil fuel combustion is 7.83 Pg C yr−1, which is very close to the estimate of the International Energy Agency (7.87 Pg C yr−1). By replacing national data disaggregation with sub-national data in this study, the average 95th minus 5th percentile ranges of CO2 emission for all grid points can be reduced from 417 to 68.2 Mg km−2 yr−1. The spread is reduced because the uneven distribution of per capita fuel consumptions within countries is better taken into account by using sub-national fuel consumption data directly. Significant difference in per capita CO2 emissions between urban and rural areas was found in developing countries (2.08 vs. 0.598 Mg C/(cap. × yr)), but not in developed countries (3.55 vs. 3.41 Mg C/(cap. × yr)). This implies that rapid urbanization of developing countries is very likely to drive up their emissions in the future.

Secondary organic aerosols (SOA) exert a significant influence on ambient air quality and regional climate. Recent field, laboratorial and modeling studies have confirmed that in-cloud processes contribute to a large fraction of SOA production with large space-time heterogeneity. This study evaluates the key factors that govern the production of cloud-process SOA (SOAcld) on a global scale based on the GFDL coupled chemistry-climate model AM3 in which full cloud chemistry is employed. The association between SOAcld production rate and six factors (i.e., liquid water content (LWC), total carbon chemical loss rate (TCloss), temperature, VOC/NOx, OH, and O3) is examined. We find that LWC alone determines the spatial pattern of SOAcld production, particularly over the tropical, subtropical and temperate forest regions, and is strongly correlated with SOAcld production. TCloss ranks the second and mainly represents the seasonal variability of vegetation growth. Other individual factors are essentially uncorrelated spatiotemporally to SOAcld production. We find that the rate of SOAcld production is simultaneously determined by both LWC and TCloss, but responds linearly to LWC and nonlinearly (or concavely) to TCloss. A parameterization based on LWC and TCloss can capture well the spatial and temporal variability of the process-based SOAcld formation (R2 = 0.5) and can be easily applied to global three dimensional models to represent the SOA production from cloud processes.

A global-scale three dimensional atmospheric transport and chemistry model was applied to simulate transpacific transport of Benzo[a]pyrene (BaP) emitted from Asia. The model results were compared with observations at six monitoring sites. The annual mean and seasonal variation of transport patterns and the contributions of different Asian source regions to transpacific transport flux were investigated. The episodic nature of transpacific transport was also systematically explored. Interannual variability of transpacific transport of BaP was also assessed during the period of 1948–2007. Results showed that strong enhancements of modeled BaP occurred in an area bounded by 70–80° E and 100–120° E. Air containing these elevated BaP concentrations was then delivered eastward by westerly winds. When approaching the West Coast of North America, the descending atmospheric motion carried BaP-laden air into the lower atmosphere. The transpacific transport flux was 1.6 times higher in the winter than in the summer. East Asian emission dominates the transpacific transport flux with a contribution of about 97%. Near ground concentration of BaP induced by Asian sources in North America varied between 1–20 pg m−3. A case study for observation at Cheeka Peak Observatory during March 2002–May 2002 reveals the importance of warm conveyor belt for transpacific transport. The number of days with transpacific transport flux with a factor of 0.5, 1.0, 1.5, and 2.0 larger than the running mean were 9.4%, 0.72%, 0.06% and 0.01%, respectively, implying a mild contribution of episodic transport to the long-term mean transport flux. Significant interannual fluctuation of transpacific transport of BaP was found, including a general decreasing trend during 1948–2007, and especially after the 1970s. The transpacific transport was found to be positively correlated with the Southern Oscillation Index.