-
Number concentrations, size distributions, and mixing states of aerosols are essential parameters for accurate estimations of aerosol direct and indirect effects. In this study, we develop an aerosol module, designated the Aerosol Two-dimensional bin module for foRmation and Aging Simulation (ATRAS), that can explicitly represent these parameters by considering new particle formation (NPF), black carbon (BC) aging, and secondary organic aerosol (SOA) processes. A two-dimensional bin representation is used for particles with dry diameters from 40 nm to 10 μm to resolve both aerosol sizes (12 bins) and BC mixing states (10 bins) for a total of 120 bins. The particles with diameters between 1 and 40 nm are resolved using additional eight size bins to calculate NPF. The ATRAS module is implemented in the WRF-Chem model and applied to examine the sensitivity of simulated mass, number, size distributions, and optical and radiative parameters of aerosols to NPF, BC aging, and SOA processes over East Asia during the spring of 2009. The BC absorption enhancement by coating materials is about 50% over East Asia during the spring, and the contribution of SOA processes to the absorption enhancement is estimated to be 10–20% over northern East Asia and 20–35% over southern East Asia. A clear north–south contrast is also found between the impacts of NPF and SOA processes on cloud condensation nuclei (CCN) concentrations: NPF increases CCN concentrations at higher supersaturations (smaller particles) over northern East Asia, whereas SOA increases CCN concentrations at lower supersaturations (larger particles) over southern East Asia. The application of ATRAS in East Asia also shows that the impact of each process on each optical and radiative parameter depends strongly on the process and the parameter in question. The module can be used in the future as a benchmark model to evaluate the accuracy of simpler aerosol models and examine interactions between NPF, BC aging, and SOA processes under different meteorological conditions and emissions.
-
Biogenic NO<sub>x</sub> emissions from soils are a large natural source with substantial uncertainties in global bottom-up estimates (ranging from 4 to 15 Tg N yr<sup>−1</sup>). We reduce this range in emission estimates, and present a top-down soil NO<sub>x</sub> emission inventory for 2005 based on retrieved tropospheric NO<sub>2</sub> columns from the Ozone Monitoring Instrument (OMI). We use a state-of-science soil NO<sub>x</sub> emission inventory (Hudman et al., 2012) as a priori in the GEOS-Chem chemistry transport model to identify 11 regions where tropospheric NO<sub>2</sub> columns are dominated by soil NO<sub>x</sub> emissions. Strong correlations between soil NO<sub>x</sub> emissions and simulated NO<sub>2</sub> columns indicate that spatial patterns in simulated NO<sub>2</sub> columns in these regions indeed reflect the underlying soil NO<sub>x</sub> emissions. Subsequently, we use a mass-balance approach to constrain emissions for these 11 regions on all major continents using OMI observed and GEOS-Chem simulated tropospheric NO<sub>2</sub> columns. We find that responses of simulated NO<sub>2</sub> columns to changing NO<sub>x</sub> emissions are suppressed over low NO<sub>x</sub> regions, and account for these non-linearities in our inversion approach. In general, our approach suggests that emissions need to be increased in most regions. Our OMI top-down soil NO<sub>x</sub> inventory amounts to 10.0 Tg N for 2005 when only constraining the 11 regions, and 12.9 Tg N when extrapolating the constraints globally. Substantial regional differences exist (ranging from −40% to +90%), and globally our top-down inventory is 4–35% higher than the GEOS-Chem a priori (9.6 Tg N yr<sup>−1</sup>). We evaluate NO<sub>2</sub> concentrations simulated with our new OMI top-down inventory against surface NO<sub>2</sub> measurements from monitoring stations in Africa, the USA and Europe. Although this comparison is complicated by several factors, we find an encouraging improved agreement when using the OMI top-down inventory compared to using the a priori inventory. To our knowledge, this study provides, for the first time, specific constraints on soil NO<sub>x</sub> emissions on all major continents using OMI NO<sub>2</sub> columns. Our results rule out the low end of reported soil NO<sub>x</sub> emission estimates, and suggest that global emissions are most likely around 12.9 ± 3.9 Tg N yr<sup>−1</sup>.
-
The chemical composition of organic aerosols in the atmosphere is strongly influenced by human emissions. The effect these have on the environment, human health, and climate change is determined by the molecular nature of these chemical species. The complexity of organic aerosol samples limits the ability to study the chemical composition, and therefore the associated properties and the impacts they have. Many studies have addressed the water-soluble fraction of organic aerosols and have had much success in identifying specific molecular formulas for thousands of compounds present. However, little attention is given to the water-insoluble portion, which can contain most of the fossil material that is emitted through human activity. Here we compare the organic aerosols present in water extracts and organic solvent extracts (pyridine and acetonitrile) of an ambient aerosol sample collected in a rural location that is impacted by natural and anthropogenic emission sources. A semiquantitative method was developed using proton nuclear magnetic resonance spectroscopy to determine that the amount of organic matter extracted by pyridine is comparable to that of water. Electrospray ionization Fourier transform ion cyclotron resonance mass spectra show that pyridine extracts a molecularly unique fraction of organic matter compared to water or acetonitrile, which extract chemically similar organic matter components. The molecular formulas unique to pyridine were less polar, more aliphatic, and reveal formulas containing sulfur to be an important component of insoluble aerosol organic matter.
-
We present decadal trend estimates of surface carbon monoxide (CO) simulated using the atmospheric chemistry general circulation model ECHAM5/MESSy (EMAC; ECHAM5 and MESSy stand for fifth-generation European Centre Hamburg general circulation model and Modular Earth Submodel System, respectively) based on the emission scenarios Representative Concentration Pathways (RCP) 8.5 for anthropogenic activity and Global Fire Emissions Database (GFED) v3.1 for biomass burning from 2001 through 2010. The spatial distribution of the modeled surface CO is evaluated with monthly data from the Measurements Of Pollution In The Troposphere (MOPITT) thermal infrared product. The global means of correlation coefficient and relative bias for the decade 2001–2010 are 0.95 and −4.29%, respectively. We also find a reasonable correlation (<i>R</i> = 0.78) between the trends of EMAC surface CO and full 10-year monthly records from ground-based observation (World Data Centre for Greenhouse Gases, WDCGG). Over western Europe, eastern USA, and northern Australia, the significant decreases in EMAC surface CO are estimated at −35.5 ± 5.8, −59.6 ± 9.1, and −13.7 ± 9.5 ppbv decade<sup>−1</sup>, respectively. In contrast, the surface CO increases by +8.9 ± 4.8 ppbv decade<sup>−1</sup> over southern Asia. A high correlation (<i>R</i> = 0.92) between the changes in EMAC-simulated surface CO and total emission flux shows that the significant regional trends are attributed to the changes in primary and direct emissions from both anthropogenic activity and biomass burning.
-
In order to investigate the impact of different treatments for the contact angle (α) in heterogeneous ice nucleating properties of natural dust and black carbon (BC) particles, we implement the classical-nucleation-theory-based parameterization of heterogeneous ice nucleation (Hoose et al., 2010) in the Community Atmospheric Model version 5 (CAM5) and then improve it by replacing the original single-contact-angle model with the probability-density-function-of-α (α-PDF) model to better represent the ice nucleation behavior of natural dust found in observations. We refit the classical nucleation theory (CNT) to constrain the uncertain parameters (i.e., onset α and activation energy in the single-α model; mean contact angle and standard deviation in the α-PDF model) using recent observation data sets for Saharan natural dust and BC (soot). We investigate the impact of the time dependence of droplet freezing on mixed-phase clouds and climate in CAM5 as well as the roles of natural dust and soot in different nucleation mechanisms. Our results show that, when compared with observations, the potential ice nuclei (IN) calculated by the α-PDF model show better agreement than those calculated by the single-α model at warm temperatures (<i>T</i>; <i>T</i> > −20 °C). More ice crystals can form at low altitudes (with warm temperatures) simulated by the α-PDF model than compared to the single-α model in CAM5. All of these can be attributed to different ice nucleation efficiencies among aerosol particles, with some particles having smaller contact angles (higher efficiencies) in the α-PDF model. In the sensitivity tests with the α-PDF model, we find that the change in mean contact angle has a larger impact on the active fraction at a given temperature than a change in standard deviation, even though the change in standard deviation can lead to a change in freezing behavior. Both the single-α and the α-PDF model indicate that the immersion freezing of natural dust plays a more important role in the heterogeneous nucleation than that of soot in mixed-phase clouds. The new parameterizations implemented in CAM5 induce more significant aerosol indirect effects than the default parameterization.
-
A method to constrain carbon dioxide (CO<sub>2</sub>) emissions from open biomass burning by using satellite observations of co-emitted species and a chemistry-transport model (CTM) is proposed and applied to the case of wildfires in Siberia. CO<sub>2</sub> emissions are assessed by means of an emission model assuming a direct relationship between the biomass burning rate (BBR) and the fire radiative power (FRP) derived from MODIS measurements. The key features of the method are (1) estimating the FRP-to-BBR conversion factors (α) for different vegetative land cover types by assimilating the satellite observations of co-emitted species into the CTM, (2) optimal combination of the estimates of α derived independently from satellite observations of different species (CO and aerosol in this study), and (3) estimation of the diurnal cycle of the fire emissions directly from the FRP measurements. Values of α for forest and grassland fires in Siberia and their uncertainties are estimated using the Infrared Atmospheric Sounding Interferometer (IASI) carbon monoxide (CO) retrievals and MODIS aerosol optical depth (AOD) measurements combined with outputs from the CHIMERE mesoscale chemistry-transport model. The constrained CO emissions are validated through comparison of the respective simulations with independent data of ground-based CO measurements at the ZOTTO site. Using our optimal regional-scale estimates of the conversion factors (which are found to be in agreement with earlier published estimates obtained from local measurements of experimental fires), the total CO<sub>2</sub> emissions from wildfires in Siberia in 2012 are estimated to be in the range from 280 to 550 Tg C, with the optimal (maximum likelihood) value of 392 Tg C. Sensitivity test cases featuring different assumptions regarding the injection height and diurnal variations of emissions indicate that the derived estimates of the total CO<sub>2</sub> emissions in Siberia are robust with respect to the modeling options (the different estimates vary within less than 15% of their magnitude). The CO<sub>2</sub> emission estimates obtained for several years are compared with independent estimates provided by the GFED3.1 and GFASv1.0 global emission inventories. It is found that our "top-down" estimates for the total annual biomass burning CO<sub>2</sub> emissions in the period from 2007 to 2011 in Siberia are by factors of 2.5 and 1.8 larger than the respective bottom-up estimates; these discrepancies cannot be fully explained by uncertainties in our estimates. There are also considerable differences in the spatial distribution of the different emission estimates; some of those differences have a systematic character and require further analysis.
-
Naturally produced very short-lived substances (VSLS) account for almost a quarter of the current stratospheric inorganic bromine, Br<sub>y</sub>. Following VSLS oxidation, bromine radicals (Br and BrO) can catalytically destroy ozone. The extent to which possible increases in surface emissions or transport of these VSLS bromocarbons to the stratosphere could counteract the effect of halogen reductions under the Montreal Protocol is an important policy question. Here, by using a chemistry–climate model, UM-UKCA, we investigate the impact of a hypothetical doubling (an increase of 5 ppt Br<sub>y</sub>) of VSLS bromocarbons on ozone and how the resulting ozone changes depend on the background concentrations of chlorine and bromine. Our model experiments indicate that for the 5 ppt increase in Br<sub>y</sub> from VSLS, the ozone decrease in the lowermost stratosphere of the Southern Hemisphere (SH) may reach up to 10% in the annual mean; the ozone decrease in the Northern Hemisphere (NH) is smaller (4–6%). The largest impact on the ozone column is found in the Antarctic spring. There is a significantly larger ozone decrease following the doubling of the VSLS burden under a high stratospheric chlorine background than under a low chlorine background, indicating the importance of the inter-halogen reactions. For example, the decline in the high-latitude, lower-stratospheric ozone concentration as a function of Br<sub>y</sub> is higher by about 30–40% when stratospheric Cl<sub>y</sub> is ~ 3 ppb (present day), compared with Cl<sub>y</sub> of ~ 0.8 ppb (a pre-industrial or projected future situation). Bromine will play an important role in the future ozone layer. However, even if bromine levels from natural VSLS were to increase significantly later this century, changes in the concentration of ozone will likely be dominated by the decrease in anthropogenic chlorine. Our calculation suggests that for a 5 ppt increase in Br<sub>y</sub> from VSLS, the Antarctic ozone hole recovery date could be delayed by approximately 6–8 years, depending on Cl<sub>y</sub> levels.
-
Global model data from the European Centre for Medium-Range Weather Forecasts (ECMWF) are analyzed for resolved gravity waves (GWs). Based on fitted 3-D wave vectors of individual waves and using the ECMWF global scale background fields, backward ray tracing from 25 km altitude is performed. Different sources such as orography, convection and winter storms are identified. It is found that due to oblique propagation waves spread widely from narrow source regions. Gravity waves which originate from regions of strong convection are frequently excited around the tropopause and have in the ECMWF model low phase and group velocities as well as very long horizontal wavelengths compared to other models and to measurements. While the total amount of momentum flux for convective GWs changes little over season, GWs generated by storms and mountain waves show large day-to-day variability, which has a strong influence also on total hemispheric fluxes; from one day to the next the total hemispheric flux may increase by a factor of 3. Implications of these results for using the ECMWF model in predicting, analyzing and interpreting global GW distributions as well as implications for seamless climate prediction are discussed.
-
Fog-induced visibility reduction is responsible for a variety of hazards in the transport sector. Therefore there is a large demand for an improved understanding of fog formation and thus improved forecasts. Improved fog forecasts require a better understanding of the numerous complex mechanisms during the fog life cycle. During winter 2012/13 a field campaign called ParisFog aiming at fog research took place at SIRTA (Instrumented Site for Atmospheric Remote Sensing Research). SIRTA is located about 20 km southwest of the Paris city center, France, in a semi-urban environment. In situ activation properties of the prevailing fog were investigated by measuring (1) total and interstitial (non-activated) dry particle number size distributions behind two different inlet systems; (2) interstitial hydrated aerosol and fog droplet size distributions at ambient conditions; and (3) cloud condensation nuclei (CCN) number concentration at different supersaturations (SS) with a CCN counter. The aerosol particles were characterized regarding their hygroscopic properties, fog droplet activation behavior and contribution to light scattering for 17 developed fog events. Low particle hygroscopicity with an overall median of the hygroscopicity parameter, κ, of 0.14 was found, likely caused by substantial influence from local traffic and wood burning emissions. Measurements of the aerosol size distribution at ambient RH revealed that the critical wet diameter, above which the hydrated aerosols activate to fog droplets, is rather large (with a median value of 2.6μm) and is highly variable (ranging from 1 to 5μm) between the different fog events. Thus, the number of activated fog droplets was very small and the non-activated hydrated particles were found to contribute significantly to the observed light scattering and thus to the reduction in visibility. Combining all experimental data, the effective peak supersaturation, SS<sub>peak</sub>, a measure of the peak supersaturation during the fog formation, was determined. The median SS<sub>peak</sub> value was estimated to be in the range from 0.031 to 0.046% (upper and lower limit estimations), which is in good agreement with previous experimental and modeling studies of fog.
-
This work reports the first comprehensive analysis of methanol concentrations in rainwater. Methanol concentrations measured in 49 rain events collected between 28 August 2007 and 10 July 2008 in Wilmington, NC, USA, ranged from below the detection limit of 6 nM to 9.3 μM with a volume-weighted average concentration of 1 ± 0.2 μM. Methanol concentrations in rainwater were up to ~200 times greater than concentrations reported previously in marine waters, indicating wet deposition as a potentially significant source of methanol to marine waters. Assuming that these methanol concentrations are an appropriate proxy for global methanol rainwater concentrations, the global methanol wet deposition sink is estimated as 20 Tg yr<sup>−1</sup>, which implies that previous methanol budgets underestimate removal by precipitation. Methanol concentrations in rainwater did not correlate significantly with H<sup>+</sup>, NO<sub>3</sub><sup>−</sup>, and NSS, which suggests that the dominant source of the alcohol to rainwater is not anthropogenic. However, methanol concentrations were strongly correlated with acetaldehyde, which has a primarily biogenic input. The methanol volume-weighted concentration during the summer (2.7 ± 0.9 μM) was ~3 times that of the winter (0.9 ± 0.2 μM), further promoting biogenic emissions as the primary cause of temporal variations of methanol concentrations. Methanol concentrations peaked in rainwater collected during the time period 12 p.m.–6 p.m. Peaking during this period of optimal sunlight implies a possible relationship with photochemical methanol production, but there are also increases in biogenic activity during this time period. Rain events with terrestrial origin had greater concentrations than those of marine origin, demonstrating the significance of the continental source of methanol in rainwater.
