-
We report the first measurements of particulate lead (Pb) from Aerodyne Aerosol Mass Spectrometers, which were deployed in and around Mexico City during the Megacity Initiative: Local and Global Research Observations (MILAGRO)/Mexico City Metropolitan Area 2006 (MCMA-2006) field campaigns. The high resolution mass spectrometer of one of the AMS instruments (HR-AMS) and the measured isotopic ratios unequivocally prove the detection of Pb in ambient particles. A substantial fraction of the lead evaporated slowly from the vaporizer of the instruments, which is indicative of species with low volatility at 600 °C. A model was developed in order to estimate the ambient particulate Pb entering the AMS from the signals in the "open" and the "closed" (or "background") mass spectrum modes of the AMS. The model suggests the presence of at least two lead fractions with ~25% of the Pb signal exhibiting rapid evaporation (1/e decay constant, τ<0.1 s) and ~75% exhibiting slow evaporation (τ~2.4 min) at the T0 urban supersite and a different fraction (70% prompt and 30% slow evaporation) at a site northwest from the metropolitan area (PEMEX site). From laboratory experiments with pure Pb(NO<sub>3</sub>)<sub>2</sub> particles, we estimated that the Pb ionization efficiency relative to nitrate (RIE<sub>Pb</sub>) is 0.5. Comparison of time series of AMS Pb with other measurements carried out at the T0 supersite during MILAGRO (using Proton Induced X-ray Emission (PIXE), Inductively-Coupled Plasma Mass Spectrometry (ICP-MS) and single-particle counts from an Aerosol Time-of-Fight Mass Spectrometer (ATOFMS)) shows similar levels (for PIXE and ICP-MS) and substantial correlation. During part of the campaign, sampling at T0 was alternated every 10 min with an Aerosol Concentrator, which enabled the detection of signals for PbCl<sup>+</sup> and PbS<sup>+</sup> ions. PbS<sup>+</sup> displays the signature of a slowly evaporating species, while PbCl<sup>+</sup> appears to arise only from fast evaporation, which is likely due to the higher vapor pressure of the compounds generating PbCl<sup>+</sup>. This is consistent with the evaporation model results. Levels of particulate Pb measured at T0 were similar to previous studies in Mexico City. Pb shows a diurnal cycle with a maximum in the early morning, which is typical of primary urban pollutants. Pb shows correlation with Zn, consistent with previous studies, while the sources of Pb appear to be at least partially disjoint from those of particulate chloride. Back trajectory analysis of the T0 Pb data suggests the presence of sources inside the urban area SSW and N of T0, with different chemical forms of Pb being associated with different source locations. High signals due to particulate lead were also detected in the PEMEX site; again, no correlation between Pb and chloride plumes was observed, suggesting mostly different sources for both species.
-
Measurements of the sensitivity of organic aerosol (OA, and its components) mass to changes in temperature were recently reported by Huffman et al.~(2009) using a tandem thermodenuder-aerosol mass spectrometer (TD-AMS) system in Mexico City and the Los Angeles area. Here, we use these measurements to derive quantitative estimates of aerosol volatility within the framework of absorptive partitioning theory using a kinetic model of aerosol evaporation in the TD. OA volatility distributions (or "basis-sets") are determined using several assumptions as to the enthalpy of vaporization (Δ<i>H</i><sub>vap</sub>). We present two definitions of "non-volatile OA," one being a global and one a local definition. Based on these definitions, our analysis indicates that a substantial fraction of the organic aerosol is comprised of non-volatile components that will not evaporate under any atmospheric conditions; on the order of 50–80% when the most realistic Δ<i>H</i><sub>vap</sub> assumptions are considered. The sensitivity of the total OA mass to dilution and ambient changes in temperature has been assessed for the various Δ<i>H</i><sub>vap</sub> assumptions. The temperature sensitivity is relatively independent of the particular Δ<i>H</i><sub>vap</sub> assumptions whereas dilution sensitivity is found to be greatest for the low (Δ<i>H</i><sub>vap</sub> = 50 kJ/mol) and lowest for the high (Δ<i>H</i><sub>vap</sub> = 150 kJ/mol) assumptions. This difference arises from the high Δ<i>H</i><sub>vap</sub> assumptions yielding volatility distributions with a greater fraction of non-volatile material than the low Δ<i>H</i><sub>vap</sub> assumptions. If the observations are fit using a 1 or 2-component model the sensitivity of the OA to dilution is unrealistically high. An empirical method introduced by Faulhaber et al. (2009) has also been used to independently estimate a volatility distribution for the ambient OA and is found to give results consistent with the high and variable Δ<i>H</i><sub>vap</sub> assumptions. Our results also show that the amount of semivolatile gas-phase organics in equilibrium with the OA could range from ~20% to 400% of the OA mass, with smaller values generally corresponding to the higher Δ<i>H</i><sub>vap</sub> assumptions. The volatility of various OA components determined from factor analysis of AMS spectra has also been assessed. In general, it is found that the fraction of non-volatile material follows the pattern: biomass burning OA < hydrocarbon-like OA < semivolatile oxygenated OA < low-volatility oxygenated OA. Correspondingly, the sensitivity to dilution and the estimated amount of semivolatile gas-phase material for the OA factors follows the reverse order. Primary OA has a substantial semivolatile fraction, in agreement with previous results, while the non-volatile fraction appears to be dominated by oxygenated OA produced by atmospheric aging. The overall OA volatility is thus controlled by the relative contribution of each aerosol type to the total OA burden. Finally, the model/measurement comparison appears to require OA having an evaporation coefficient (γ<sub><i>e</i></sub>) substantially greater than 10<sup>−2</sup>; at this point it is not possible to place firmer constraints on γ<sub><i>e</i></sub> based on the observations.
-
Variations in the chemical composition of atmospheric aerosols alter their hygroscopicity and can lead to changes in the cloud-active fraction of the aerosols, or cloud condensation nuclei (CCN) number concentration. To investigate the importance of this effect under different atmospheric conditions, cloud droplet formation was simulated with a Lagrangian parcel model. Initial values of updraft speed and temperature were systematically varied along with aerosol number concentration, size and hygroscopicity (represented by the hygroscopicity parameter, κ). A previous study classifies the sensitivity of CCN activity to compositional changes based on the supersaturation reached in the parcel model. We found that these classifications could not be generalized to a range of aerosol size distribution median radii. Instead, variations in sensitivity with size depend on the location of the dry critical radius for droplet activation relative to the size distribution median radius. The parcel model output was used to construct droplet activation lookup tables based on κ that were implemented in the Regional Atmospheric Modeling System (RAMS) microphysical scheme. As a first application of this system, aerosol hygroscopicity and size were varied in a series of RAMS mesoscale simulations designed to investigate the sensitivity of a mixed-phase orographic cloud case to the parameter variations. Observations from a recent field campaign in northwestern Colorado provided the basis for the aerosol field initializations. Model results show moderate sensitivity in the distribution of total case precipitation to extreme changes in κ, and minimal sensitivity to observed changes in estimated κ. The impact of varying aerosol hygroscopicity diminished with increasing median radius, as expected from the parcel model results. The conclusions drawn from these simulations could simplify similar research in other cloud regimes by defining the need, or lack of need, for detailed knowledge of aerosol composition.
-
Atmospheric deposition was measured at three sites in the Pensacola Bay watershed, Florida, between November 2004 and December 2007. Mercury deposition in the Pensacola Bay watershed was similar to that from nearby Mercury Deposition Network sites along the Northern Gulf of Mexico coast. Mercury deposition during the summer months is higher than other months due to higher concentrations in rainfall throughout the region. Deposition of constituents like H<sup>+</sup>, sulfate, nitrate, ammonium, chloride and sodium, were much higher in Pensacola Bay that at National Atmospheric Deposition Program (NADP) sites. Chloride and sodium deposition are higher because Pensacola Bay sites are closer to the Gulf of Mexico which is a source of sea salt aerosols. Acid rain constituents, H<sup>+</sup>, sulfate, nitrate and ammonium are most likely higher at Pensacola Bay sites because these sites are much closer to emission sources of these constituents than NADP sites, particularly two Florida NADP sites, FL14 and FL23, which are located in rural counties far from major industrial activities.
-
It has been established that observed local and regional levels of secondary organic aerosols (SOA) in polluted areas cannot be explained by the oxidation and partitioning of anthropogenic and biogenic VOC precursors, at least using current mechanisms and parameterizations. In this study, the 3-D regional air quality model CHIMERE is applied to estimate the potential contribution to SOA formation of recently identified semi-volatile and intermediate volatility organic precursors (S/IVOC) in and around Mexico City for the MILAGRO field experiment during March 2006. The model has been updated to include explicitly the volatility distribution of primary organic aerosols (POA), their gas-particle partitioning and the gas-phase oxidation of the vapors. Two recently proposed parameterizations, those of Robinson et al. (2007) ("ROB") and Grieshop et al. (2009) ("GRI") are compared and evaluated against surface and aircraft measurements. The 3-D model results are assessed by comparing with the concentrations of OA components from Positive Matrix Factorization of Aerosol Mass Spectrometer (AMS) data, and for the first time also with oxygen-to-carbon ratios derived from high-resolution AMS measurements. The results show a substantial enhancement in predicted SOA concentrations (2–4 times) with respect to the previously published base case without S/IVOCs (Hodzic et al., 2009), both within and downwind of the city leading to much reduced discrepancies with the total OA measurements. Model improvements in OA predictions are associated with the better-captured SOA magnitude and diurnal variability. The predicted production from anthropogenic and biomass burning S/IVOC represents 40–60% of the total measured SOA at the surface during the day and is somewhat larger than that from commonly measured aromatic VOCs, especially at the T1 site at the edge of the city. The SOA production from the continued multi-generation S/IVOC oxidation products continues actively downwind. Similar to aircraft observations, the predicted OA/ΔCO ratio for the ROB case increases from 20–30 μg sm<sup>−3</sup> ppm<sup>−1</sup> up to 60–70 μg sm<sup>−3</sup> ppm<sup>−1</sup> between a fresh and 1-day aged air mass, while the GRI case produces a 30% higher OA growth than observed. The predicted average O/C ratio of total OA for the ROB case is 0.16 at T0, substantially below observed value of 0.5. A much better agreement for O/C ratios and temporal variability (<i>R</i><sup>2</sup>=0.63) is achieved with the updated GRI treatment. Both treatments show a deficiency in regard to POA ageing with a tendency to over-evaporate POA upon dilution of the urban plume suggesting that atmospheric HOA may be less volatile than assumed in these parameterizations. This study highlights the important potential role of S/IVOC chemistry in the SOA budget in this region, and highlights the need for further improvements in available parameterizations. The agreement observed in this study is not sufficient evidence to conclude that S/IVOC are the major missing SOA source in megacity environments. The model is still very underconstrained, and other possible pathways such as formation from very volatile species like glyoxal may explain some of the mass and especially increase the O/C ratio.
-
Secondary organic aerosol (SOA) formation in the atmosphere is currently often modeled using a multiple lumped "two-product" (<I>N</I>·2p) approach. The <I>N</I>·2p approach neglects: 1) variation of activity coefficient (ζ<sub><I>i</I></sub>) values and mean molecular weight <span style="border-top: 1px solid #000; color: #000;">MW</span> in the particulate matter (PM) phase; 2) water uptake into the PM; and 3) the possibility of phase separation in the PM. This study considers these effects by adopting an (<I>N</I>·2p)<sup>ζ<sub>p</sub><span style="border-top: 1px solid #000; color: #000;">MW</span>,ζ</sup> approach (θ is a phase index). Specific chemical structures are assigned to 25 lumped SOA compounds and to 15 representative primary organic aerosol (POA) compounds to allow calculation of ζ<sub><I>i</I></sub> and <span style="border-top: 1px solid #000; color: #000;">MW</span> values. The SOA structure assignments are based on chamber-derived 2p gas/particle partition coefficient values coupled with known effects of structure on vapor pressure <I>p</I><sub>L,<I>i</I></sub><sup>o</sup> (atm). To facilitate adoption of the (<I>N</I>·2p)<sup>ζ<sub><I>p</I></sub><span style="border-top: 1px solid #000; color: #000;">MW</span>,θ</sup> approach in large-scale models, this study also develops CP-Wilson.1 (Chang-Pankow-Wilson.1), a group-contribution ζ<sub><I>i</I></sub>-prediction method that is more computationally economical than the UNIFAC model of Fredenslund et al. (1975). Group parameter values required by CP-Wilson.1 are obtained by fitting ζ<sub><I>i</I></sub> values to predictions from UNIFAC. The (<I>N</I>·2p)<sup>ζ<sub><I>p</I></sub><span style="border-top: 1px solid #000; color: #000;">MW</span>,θ</sup> approach is applied (using CP-Wilson.1) to several real α-pinene/O<sub>3</sub> chamber cases for high reacted hydrocarbon levels (ΔHC≈400 to 1000 μg m<sup>−3</sup>) when relative humidity (RH) ≈50%. Good agreement between the chamber and predicted results is obtained using both the (<I>N</I>·2p)<sup>ζ<sub><I>p</I></sub><span style="border-top: 1px solid #000; color: #000;">MW</span>,θ</sup> and <I>N</I>·2p approaches, indicating relatively small water effects under these conditions. However, for a hypothetical α-pinene/O<sub>3</sub> case at ΔHC=30 μg m<sup>−3</sup> and RH=50%, the (<I>N</I>·2p)<sup>ζ<sub><I>p</I></sub><span style="border-top: 1px solid #000; color: #000;">MW</span>,θ</sup> approach predicts that water uptake will lead to an organic PM level that is more double that predicted by the </I>N</I>·2p approach. Adoption of the (<I>N</I>·2p)<sup>ζ<sub><I>p</I></sub><span style="border-top: 1px solid #000; color: #000;">MW</span>,θ</sup> approach using reasonable lumped structures for SOA and POA compounds is recommended for ambient PM modeling.
-
A statistical cirrus cloud scheme that accounts for mesoscale temperature perturbations is implemented in a coupled aerosol and atmospheric circulation model to better represent both subgrid-scale supersaturation and cloud formation. This new scheme treats the effects of aerosol on cloud formation and ice freezing in an improved manner, and both homogeneous freezing and heterogeneous freezing are included. The scheme is able to better simulate the observed probability distribution of relative humidity compared to the scheme that was implemented in an older version of the model. Heterogeneous ice nuclei (IN) are shown to decrease the frequency of occurrence of supersaturation, and improve the comparison with observations at 192 hPa. Homogeneous freezing alone can not reproduce observed ice crystal number concentrations at low temperatures (<205 K), but the addition of heterogeneous IN improves the comparison somewhat. Increases in heterogeneous IN affect both high level cirrus clouds and low level liquid clouds. Increases in cirrus clouds lead to a more cloudy and moist lower troposphere with less precipitation, effects which we associate with the decreased convective activity. The change in the net cloud forcing is not very sensitive to the change in ice crystal concentrations, but the change in the net radiative flux at the top of the atmosphere is still large because of changes in water vapor. Changes in the magnitude of the assumed mesoscale temperature perturbations by 25% alter the ice crystal number concentrations and the net radiative fluxes by an amount that is comparable to that from a factor of 10 change in the heterogeneous IN number concentrations. Further improvements on the representation of mesoscale temperature perturbations, heterogeneous IN and the competition between homogeneous freezing and heterogeneous freezing are needed.
-
The heterogeneous oxidation of laboratory Secondary Organic Aerosol (SOA) particles by OH radicals was investigated. SOA particles, produced by reaction of α-pinene and O<sub>3</sub>, were exposed to OH radicals in a flow tube, and particle chemical composition, size, and hygroscopicity were measured to assess modifications due to oxidative aging. Aerosol Mass Spectrometer (AMS) mass spectra indicated that the degree of oxidation of 200 nm diameter SOA particles was significantly enhanced due to OH-initiated oxidation, as evidenced by the increase in the fraction of <i>m/z</i> 44 fragment of total organic mass concentration (F44). F44 values of the SOA particles, initially in the range F44=0.04–0.07, increased by up to ΔF44~0.01 under equivalent atmospheric aging timescales of 2 weeks, assuming a 24-h average OH concentration of 10<sup>6</sup> cm<sup>−3</sup>. Particle O/C ratios calculated from F44 values, initially in the range O/C=0.25–0.35, rose by a maximum of ΔO/C~0.04 units for 2 weeks of aging. Particle densities also increased with heterogeneous oxidation, consistent with the observed increase in the degree of oxidation. Minor reductions in particle size, with volume losses of up to 10%, were observed due to volatilization of oxidation products. The SOA particles activated more readily to form cloud droplets with an increase in the κ hygroscopicity parameter of up to a factor of two for the equivalent of 2 weeks of OH atmospheric exposure. These results indicate that OH heterogeneous oxidation of typical SOA needs to be considered as an atmospheric organic aerosol aging mechanism, most likely of higher relative importance away from VOC source regions, where other aging mechanisms are less dominant.
-
Temperature inversions are one of the dominant features of the Arctic atmosphere and play a crucial role in various processes by controlling the transfer of mass and moisture fluxes through the lower troposphere. It is therefore essential that they are accurately quantified, monitored and simulated as realistically as possible over the Arctic regions. In the present study, the characteristics of inversions in terms of frequency and strength are quantified for the entire Arctic Ocean for summer and winter seasons of 2003 to 2008 using the AIRS data for the clear-sky conditions. The probability density functions (PDFs) of the inversion strength are also presented for every summer and winter month. <br><br> Our analysis shows that although the inversion frequency along the coastal regions of Arctic decreases from June to August, inversions are still seen in almost each profile retrieved over the inner Arctic region. In winter, inversions are ubiquitous and are also present in every profile analysed over the inner Arctic region. When averaged over the entire study area (70° N–90° N), the inversion frequency in summer ranges from 69 to 86% for the ascending passes and 72–86% for the descending passes. For winter, the frequency values are 88–91% for the ascending passes and 89–92% for the descending passes of AIRS/AQUA. The PDFs of inversion strength for the summer months are narrow and right-skewed (or positively skewed), while in winter, they are much broader. In summer months, the mean values of inversion strength for the entire study area range from 2.5 to 3.9 K, while in winter, they range from 7.8 to 8.9 K. The standard deviation of the inversion strength is double in winter compared to summer. The inversions in the summer months of 2007 were very strong compared to other years. The warming in the troposphere of about 1.5–3.0 K vertically extending up to 400 hPa was observed in the summer months of 2007.
-
Processes occurring in the tropical upper troposphere (UT), the Tropical Transition Layer (TTL), and the lower stratosphere (LS) are of importance for the global climate, for stratospheric dynamics and air chemistry, and for their influence on the global distribution of water vapour, trace gases and aerosols. In this contribution we present aerosol and trace gas (in-situ) measurements from the tropical UT/LS over Southern Brazil, Northern Australia, and West Africa. The instruments were operated on board of the Russian high altitude research aircraft M-55 "Geophysica" and the DLR Falcon-20 during the campaigns TROCCINOX (Araçatuba, Brazil, February 2005), SCOUT-O3 (Darwin, Australia, December 2005), and SCOUT-AMMA (Ouagadougou, Burkina Faso, August 2006). The data cover submicron particle number densities and volatility from the COndensation PArticle counting System (COPAS), as well as relevant trace gases like N<sub>2</sub>O, ozone, and CO. We use these trace gas measurements to place the aerosol data into a broader atmospheric context. Also a juxtaposition of the submicron particle data with previous measurements over Costa Rica and other tropical locations between 1999 and 2007 (NASA DC-8 and NASA WB-57F) is provided. The submicron particle number densities, as a function of altitude, were found to be remarkably constant in the tropical UT/LS altitude band for the two decades after 1987. Thus, a parameterisation suitable for models can be extracted from these measurements. Compared to the average levels in the period between 1987 and 2007 a slight increase of particle abundances was found for 2005/2006 at altitudes with potential temperatures, Θ, above 430 K. The origins of this increase are unknown except for increases measured during SCOUT-AMMA. Here the eruption of the Soufrière Hills volcano in the Caribbean caused elevated particle mixing ratios. The vertical profiles from Northern hemispheric mid-latitudes between 1999 and 2006 also are compact enough to derive a parameterisation. The tropical profiles all show a broad maximum of particle mixing ratios (between Θ≈340 K and 390 K) which extends from below the TTL to above the thermal tropopause. Thus these particles are a "reservoir" for vertical transport into the stratosphere. The ratio of non-volatile particle number density to total particle number density was also measured by COPAS. The vertical profiles of this ratio have a maximum of 50% above 370 K over Australia and West Africa and a pronounced minimum directly below. Without detailed chemical composition measurements a reason for the increase of non-volatile particle fractions cannot yet be given. However, half of the particles from the tropical "reservoir" contain compounds other than sulphuric acid and water. Correlations of the measured aerosol mixing ratios with N<sub>2</sub>O and ozone exhibit compact relationships for the tropical data from SCOUT-AMMA, TROCCINOX, and SCOUT-O3. Correlations with CO are more scattered probably because of the connection to different pollution source regions. We provide additional data from the long distance transfer flights to the campaign sites in Brazil, Australia, and West-Africa. These were executed during a time window of 17 months within a period of relative volcanic quiescence. Thus the data represent a "snapshot picture" documenting the status of a significant part of the global UT/LS fine aerosol at low concentration levels 15 years after the last major (i.e., the 1991 Mount Pinatubo) eruption. The corresponding latitudinal distributions of the measured particle number densities are presented in this paper to provide data of the UT/LS background aerosol for modelling purposes.