-
Extreme haze events have occurred frequently over China in recent years. Although many studies have investigated the formation mechanisms associated with PM<sub>2.5</sub> for heavily polluted regions in China based on observational data, adequately predicting peak PM<sub>2.5</sub> concentrations is still challenging for regional air quality models. In this study, we evaluate the performance of one configuration of the Weather Research and Forecasting model coupled with chemistry (WRF-Chem) and use the model to investigate the sensitivity of heterogeneous reactions on simulated peak sulfate, nitrate, and ammonium concentrations in the vicinity of Beijing during four extreme haze episodes in October 2014 over the North China Plain. The highest observed PM<sub>2.5</sub> concentration of 469 µg m<sup>−3</sup> occurred in Beijing. Comparisons with observations show that the model reproduced the temporal variability in PM<sub>2.5</sub> with the highest PM<sub>2.5</sub> values on polluted days (defined as days in which observed PM<sub>2.5</sub> is greater than 75 µg m<sup>−3</sup>), but predictions of sulfate, nitrate, and ammonium were too low on days with the highest observed concentrations. Observational data indicate that the sulfur/nitric oxidation rates are strongly correlated with relative humidity during periods of peak PM<sub>2.5</sub>; however, the model failed to reproduce the highest PM<sub>2.5</sub> concentrations due to missing heterogeneous/aqueous reactions. As the parameterizations of those heterogeneous reactions are not well established yet, estimates of SO<sub>2</sub>-to-H<sub>2</sub>SO<sub>4</sub> and NO<sub>2</sub>/NO<sub>3</sub>-to-HNO<sub>3</sub> reaction rates that depend on relative humidity were applied, which improved the simulation of sulfate, nitrate, and ammonium enhancement on polluted days in terms of both concentrations and partitioning among those species. Sensitivity simulations showed that the extremely high heterogeneous reaction rates and also higher emission rates than those reported in the emission inventory were likely important factors contributing to those peak PM<sub>2.5</sub> concentrations.
-
Long-term (2009–2015) observations of atmospheric black carbon (BC) mass concentrations were performed using a continuous soot-monitoring system (COSMOS) at Fukue Island, western Japan, to provide information on wet removal rate constraints and the emission strengths of important source regions in East Asia (China and others). The annual average mass concentration was 0.36 µg m<sup>−3</sup>, with distinct seasonality; high concentrations were recorded during autumn, winter, and spring and were caused by Asian continental outflows, which reached Fukue Island in 6–46 h. The observed data were categorized into two classes, i.e., with and without a wet removal effect, using the accumulated precipitation along a backward trajectory (APT) for the last 3 days as an index. Statistical analysis of the observed ΔBC ∕ ΔCO ratios was performed to obtain information on the emission ratios (from data with zero APT only) and wet removal rates (including data with nonzero APTs). The estimated emission ratios (5.2–6.9 ng m<sup>−3</sup> ppb<sup>−1</sup>) varied over the six air mass origin areas; the higher ratios for south-central East China (30–35° N) than for north-central East China (35–40° N) indicated the relative importance of domestic emissions and/or biomass burning sectors. The significantly higher BC ∕ CO emission ratios adopted in the bottom-up Regional Emission inventory in Asia (REAS) version 2 (8.3–23 ng m<sup>−3</sup> ppb<sup>−1</sup>) over central East China and Korea needed to be reduced at least by factors of 1.3 and 2.8 for central East China and Korea, respectively, but the ratio for Japan was reasonable. The wintertime enhancement of the BC emission from China, predicted by REAS2, was verified for air masses from south-central East China but not for those from north-central East China. Wet removal of BC was clearly identified as a decrease in the ΔBC ∕ ΔCO ratio against APT. The transport efficiency (TE), defined as the ratio of the ΔBC ∕ ΔCO ratio with precipitation to that without precipitation, was fitted reasonably well by a stretched exponential decay curve against APT; a single set of fitting parameters was sufficient to represent the results for air masses originating from different areas. An accumulated precipitation of 25.5 ± 6.1 mm reduced the TE to 1∕<i>e</i>. BC-containing particles traveling to Fukue must have already been converted from hydrophobic to hydrophilic particles, because the behavior of TE against APT was similar to that of PM<sub>2.5</sub>, the major components of which are hydrophilic. Wet loss of BC greatly influenced interannual variations in the ΔBC ∕ ΔCO ratios and BC mass concentrations. This long-term data set will provide a benchmark for testing chemical transport/climate model simulations covering East Asia.
-
The 2014–2015 Holuhraun lava-flood eruption of Bárðarbunga volcano (Iceland) emitted prodigious amounts of sulfur dioxide into the atmosphere. This eruption caused a large-scale episode of air pollution throughout Western Europe in September 2014, the first event of this magnitude recorded in the modern era. We gathered chemistry-transport simulations and a wealth of complementary observations from satellite sensors (OMI, IASI), ground-based remote sensing (lidar, sunphotometry, differential optical absorption spectroscopy) and ground-level air quality monitoring networks to characterize both the spatial-temporal distributions of volcanic SO<sub>2</sub> and sulfate aerosols as well as the dynamics of the planetary boundary layer. Time variations of dynamical and microphysical properties of sulfate aerosols in the aged low-tropospheric volcanic cloud, including loading, vertical distribution, size distribution and single scattering albedo, are provided. Retrospective chemistry-transport simulations at low horizontal resolution (25 km × 25 km) capture the correct temporal dynamics of this far-range air pollution event but fail to reproduce the correct magnitude of SO<sub>2</sub> concentration at ground-level. Simulations at higher spatial resolution, relying on two nested domains with finest resolution of 7.3 km × 7.3 km, improve substantially the far-range vertical distribution of the volcanic cloud and subsequently the description of ground-level SO<sub>2</sub> concentrations. However, remaining discrepancies between model and observations are shown to result from an inaccurate representation of the planetary boundary layer (PBL) dynamics. Comparison with lidar observations points out a systematic under-estimation of the PBL height by the model, whichever the PBL parameterization scheme. Such a shortcoming impedes the capture of the overlying Bárðarbunga cloud into the PBL at the right time and in sufficient quantities. This study therefore demonstrates the key role played by the PBL dynamics in accurately modelling large-scale volcanogenic air pollution.
-
The molecular structure of volatile organic compounds determines their oxidation pathway, directly impacting secondary organic aerosol (SOA) formation. This study comprehensively investigates the impact of molecular structure on SOA formation from the photooxidation of 12 different eight- to nine-carbon aromatic hydrocarbons under low-NO<sub><i>x</i></sub> conditions. The effects of the alkyl substitute number, location, carbon chain length and branching structure on the photooxidation of aromatic hydrocarbons are demonstrated by analyzing SOA yield, chemical composition and physical properties. Aromatic hydrocarbons, categorized into five groups, show a yield order of ortho (<i>o</i>-xylene and <i>o</i>-ethyltoluene) > one substitute (ethylbenzene, propylbenzene and isopropylbenzene) > meta (<i>m</i>-xylene and <i>m</i>-ethyltoluene) > three substitute (trimethylbenzenes) > para (<i>p</i>-xylene and <i>p</i>-ethyltoluene). SOA yields of aromatic hydrocarbon photooxidation do not monotonically decrease when increasing alkyl substitute number. The ortho position promotes SOA formation while the para position suppresses aromatic oxidation and SOA formation. Observed SOA chemical composition and volatility confirm that higher yield is associated with further oxidation. SOA chemical composition also suggests that aromatic oxidation increases with increasing alkyl substitute chain length and branching structure. Further, carbon dilution conjecture developed by Li et al. (2016) is extended in this study to serve as a standard method to determine the extent of oxidation of an alkyl-substituted aromatic hydrocarbon.
-
The ability of seven state-of-the-art chemistry–aerosol models to reproduce distributions of tropospheric ozone and its precursors, as well as aerosols over eastern Asia in summer 2008, is evaluated. The study focuses on the performance of models used to assess impacts of pollutants on climate and air quality as part of the EU ECLIPSE project. Models, run using the same ECLIPSE emissions, are compared over different spatial scales to in situ surface, vertical profiles and satellite data. Several rather clear biases are found between model results and observations, including overestimation of ozone at rural locations downwind of the main emission regions in China, as well as downwind over the Pacific. Several models produce too much ozone over polluted regions, which is then transported downwind. Analysis points to different factors related to the ability of models to simulate VOC-limited regimes over polluted regions and NO<sub><i>x</i></sub> limited regimes downwind. This may also be linked to biases compared to satellite NO<sub>2</sub>, indicating overestimation of NO<sub>2</sub> over and to the north of the northern China Plain emission region. On the other hand, model NO<sub>2</sub> is too low to the south and west of this region and over South Korea/Japan. Overestimation of ozone is linked to systematic underestimation of CO particularly at rural sites and downwind of the main Chinese emission regions. This is likely to be due to enhanced destruction of CO by OH. Overestimation of Asian ozone and its transport downwind implies that radiative forcing from this source may be overestimated. Model-observation discrepancies over Beijing do not appear to be due to emission controls linked to the Olympic Games in summer 2008.<br><br>With regard to aerosols, most models reproduce the satellite-derived AOD patterns over eastern China. Our study nevertheless reveals an overestimation of ECLIPSE model mean surface BC and sulphate aerosols in urban China in summer 2008. The effect of the short-term emission mitigation in Beijing is too weak to explain the differences between the models. Our results rather point to an overestimation of SO<sub>2</sub> emissions, in particular, close to the surface in Chinese urban areas. However, we also identify a clear underestimation of aerosol concentrations over northern India, suggesting that the rapid recent growth of emissions in India, as well as their spatial extension, is underestimated in emission inventories. Model deficiencies in the representation of pollution accumulation due to the Indian monsoon may also be playing a role. Comparison with vertical aerosol lidar measurements highlights a general underestimation of scattering aerosols in the boundary layer associated with overestimation in the free troposphere pointing to modelled aerosol lifetimes that are too long. This is likely linked to too strong vertical transport and/or insufficient deposition efficiency during transport or export from the boundary layer, rather than chemical processing (in the case of sulphate aerosols). Underestimation of sulphate in the boundary layer implies potentially large errors in simulated aerosol–cloud interactions, via impacts on boundary-layer clouds.<br><br>This evaluation has important implications for accurate assessment of air pollutants on regional air quality and global climate based on global model calculations. Ideally, models should be run at higher resolution over source regions to better simulate urban–rural pollutant gradients and/or chemical regimes, and also to better resolve pollutant processing and loss by wet deposition as well as vertical transport. Discrepancies in vertical distributions require further quantification and improvement since these are a key factor in the determination of radiative forcing from short-lived pollutants.
-
This paper promotes an understanding of the mineralogical, chemical, and physical interrelationships of re-suspended mineral dusts collected as grab samples from global dust sources. Surface soils were collected from arid regions, including the southwestern USA, Mali, Chad, Morocco, Canary Islands, Cabo Verde, Djibouti, Afghanistan, Iraq, Kuwait, Qatar, UAE, Serbia, China, Namibia, Botswana, Australia, and Chile. The < 38 µm sieved fraction of each sample was re-suspended in a chamber, from which the airborne mineral dust could be extracted, sampled, and analyzed. Instruments integrated into the entrainment facility included two PM<sub>10</sub> and two PM<sub>2.5</sub> filter samplers, a beta attenuation gauge for the continuous measurement of PM<sub>10</sub> and PM<sub>2.5</sub> particulate mass fractions, an aerodynamic particle size analyzer, and a three-wavelength (405, 532, 781 nm) photoacoustic instrument with integrating reciprocal nephelometer for monitoring absorption and scattering coefficients during the dust re-suspension process. Filter sampling media included Teflon<span style="position:relative; bottom:0.5em; " class="text">®</span> membrane and quartz fiber filters for chemical analysis and Nuclepore<span style="position:relative; bottom:0.5em; " class="text">®</span> filters for individual particle analysis by scanning electron microscopy (SEM). The < 38 µm sieved fractions were also analyzed by X-ray diffraction for their mineral content while the > 75, < 125 µm soil fractions were mineralogically assessed by optical microscopy. Presented here are results of the optical measurements, showing the interdependency of single-scattering albedos (SSA) at three different wavelengths and mineralogical content of the entrained dust samples. To explain the elevated concentrations of iron (Fe) and Fe ∕ Al ratios in the soil re-suspensions, we propose that dust particles are to a large extent composed of nano-sized particles of micas, clays, metal oxides, and ions of potassium (K<sup>+</sup>), calcium (Ca<sup>2+</sup>), and sodium (Na<sup>+</sup>) evenly dispersed as a colloid or adsorbed in amorphous clay-like material. Also shown are differences in SSA of the kaolinite/hematite/goethite samples from Mali and those from colloidal soils elsewhere. Results from this study can be integrated into a database of mineral dust properties, for applications in climate modeling, remote sensing, visibility, health (medical geology), ocean fertilization, and impact on equipment.
-
The aerosol size-segregated chemical composition was analyzed at an urban (Bologna) and a rural (San Pietro Capofiume) site in the Po Valley, Italy, during June and July 2012, by ion-chromatography (major water-soluble ions and organic acids) and evolved gas analysis (total and water-soluble carbon), to investigate sources and mechanisms of secondary aerosol formation during the summer. A significant enhancement of secondary organic and inorganic aerosol mass was observed under anticyclonic conditions with recirculation of planetary boundary layer air but with substantial differences between the urban and the rural site. The data analysis, including a principal component analysis (PCA) on the size-resolved dataset of chemical concentrations, indicated that the photochemical oxidation of inorganic and organic gaseous precursors was an important mechanism of secondary aerosol formation at both sites. In addition, at the rural site a second formation process, explaining the largest fraction (22 %) of the total variance, was active at nighttime, especially under stagnant conditions. Nocturnal chemistry in the rural Po Valley was associated with the formation of ammonium nitrate in large accumulation-mode (0.42–1.2 µm) aerosols favored by local thermodynamic conditions (higher relative humidity and lower temperature compared to the urban site). Nocturnal concentrations of fine nitrate were, in fact, on average 5 times higher at the rural site than in Bologna. The water uptake by this highly hygroscopic compound under high RH conditions provided the medium for increased nocturnal aerosol uptake of water-soluble organic gases and possibly also for aqueous chemistry, as revealed by the shifting of peak concentrations of secondary compounds (water-soluble organic carbon (WSOC) and sulfate) toward the large accumulation mode (0.42–1.2 µm). Contrarily, the diurnal production of WSOC (proxy for secondary organic aerosol) by photochemistry was similar at the two sites but mostly affected the small accumulation mode of particles (0.14–0.42 µm) in Bologna, while a shift to larger accumulation mode was observed at the rural site. A significant increment in carbonaceous aerosol concentration (for both WSOC and water-insoluble carbon) at the urban site was recorded mainly in the quasi-ultrafine fraction (size range 0.05–0.14 µm), indicating a direct influence of traffic emissions on the mass concentrations of this range of particles.
-
Downward transport of ozone (O<sub>3</sub>) from the stratosphere can be a significant contributor to tropospheric O<sub>3</sub> background levels. However, this process often is not well represented in current regional models. In this study, we develop a seasonally and spatially varying potential vorticity (PV)-based function to parameterize upper tropospheric and/or lower stratospheric (UTLS) O<sub>3</sub> in a chemistry transport model. This dynamic O<sub>3</sub>–PV function is developed based on 21-year ozonesonde records from World Ozone and Ultraviolet Radiation Data Centre (WOUDC) with corresponding PV values from a 21-year Weather Research and Forecasting (WRF) simulation across the Northern Hemisphere from 1990 to 2010. The result suggests strong spatial and seasonal variations of O<sub>3</sub> ∕ PV ratios which exhibits large values in the upper layers and in high-latitude regions, with highest values in spring and the lowest values in autumn over an annual cycle. The newly developed O<sub>3</sub> ∕ PV function was then applied in the Community Multiscale Air Quality (CMAQ) model for an annual simulation of the year 2006. The simulated UTLS O<sub>3</sub> agrees much better with observations in both magnitude and seasonality after the implementation of the new parameterization. Considerable impacts on surface O<sub>3</sub> model performance were found in the comparison with observations from three observational networks, i.e., EMEP, CASTNET and WDCGG. With the new parameterization, the negative bias in spring is reduced from −20 to −15 % in the reference case to −9 to −1 %, while the positive bias in autumn is increased from 1 to 15 % in the reference case to 5 to 22 %. Therefore, the downward transport of O<sub>3</sub> from upper layers has large impacts on surface concentration and needs to be properly represented in regional models.
-
Extensive undersaturations of carbon tetrachloride (CCl<sub>4</sub>) in Pacific, Atlantic, and Southern Ocean surface waters indicate that atmospheric CCl<sub>4</sub> is consumed in large amounts by the ocean. Observations made on 16 research cruises between 1987 and 2010, ranging in latitude from 60° N to 77° S, show that negative saturations extend over most of the surface ocean. Corrected for physical effects associated with radiative heat flux, mixing, and air injection, these anomalies were commonly on the order of −5 to −10 %, with no clear relationship to temperature, productivity, or other gross surface water characteristics other than being more negative in association with upwelling. The atmospheric flux required to sustain these undersaturations is 12.4 (9.4–15.4) Gg yr<sup>−1</sup>, a loss rate implying a partial atmospheric lifetime with respect to the oceanic loss of 183 (147–241) yr and that ∼ 18 (14–22) % of atmospheric CCl<sub>4</sub> is lost to the ocean. Although CCl<sub>4</sub> hydrolyzes in seawater, published hydrolysis rates for this gas are too slow to support such large undersaturations, given our current understanding of air–sea gas exchange rates. The even larger undersaturations in intermediate depth waters associated with reduced oxygen levels, observed in this study and by other investigators, strongly suggest that CCl<sub>4</sub> is ubiquitously consumed at mid-depth, presumably by microbiota. Although this subsurface sink creates a gradient that drives a downward flux of CCl<sub>4</sub>, the gradient alone is not sufficient to explain the observed surface undersaturations. Since known chemical losses are likewise insufficient to sustain the observed undersaturations, this suggests a possible biological sink for CCl<sub>4</sub> in surface or near-surface waters of the ocean. The total atmospheric lifetime for CCl<sub>4</sub>, based on these results and the most recent studies of soil uptake and loss in the stratosphere is now 32 (26–43) yr.
-
There is an ongoing debate on whether the observed decadal variations in surface solar radiation, known as "dimming and brightening", are a global or just local phenomenon. We investigated this issue using a comprehensive set of long-term sunshine duration records from China, which experienced a rapid growth in urbanization during past decades. 172 pairs of urban and nearby rural stations were analyzed over the period 1960–1989 ("dimming phase") and 1990–2013 ("brightening phase"). There is a large overlap in urban and rural sunshine duration trends for both dimming (≈ 86 %) and brightening (≈ 84 %) phases. This indicates that rather than urban dimming or rural brightening, the global dimming and brightening phenomena are more of national/regional scale in China. In the dimming phase, sunshine duration significantly declined in both urban and rural areas at an average rate of −0.20 h d<sup>−1</sup> decade<sup>−1</sup> and −0.14 h d<sup>−1</sup> decade<sup>−1</sup> respectively, i.e. rural dimming has been around two-thirds of urban dimming. This ratio generally increases from a minimum of 0.39 to a maximum of 0.87 with increasing indices of urbanization, reaching saturation when the urbanization level exceeds 50 %, or the urban population exceeds 20 million persons, or the population density becomes higher than 250 person km<sup>−2</sup>. Urbanization can be treated as a useful indicator for anthropogenic air pollution in studying pollution-driven changes in sunshine duration during the dimming phase when pollution control and monitoring were largely absent. After the transition into the brightening phase, the increasing number of environment-related laws and regulations as well as investments in the abatement of environmental pollution might have helped in counteracting air pollutants generated during the urbanization process. Therefore, in the brightening phase, urbanization no longer simply indicates an increase in air pollution and its effect on sunshine duration becomes insignificant. In conclusion, urbanization can give a general indication of pollution-driven sunshine dimming until pollution regulations become effective.