-
We carried out column model simulations to study particle fluxes and deposition and to evaluate different particle formation mechanisms at a boreal forest site in Finland. We show that kinetic nucleation of sulphuric acid cannot be responsible for new particle formation alone as the simulated vertical profile of particle number concentration does not correspond to observations. Instead organic induced nucleation leads to good agreement confirming the relevance of the aerosol formation mechanism including organic compounds emitted by the biosphere. <br><br> The simulation of aerosol concentration within the atmospheric boundary layer during nucleation event days shows a highly dynamical picture, where particle formation is coupled with chemistry and turbulent transport. We have demonstrated the suitability of our turbulent mixing scheme in reproducing the most important characteristics of particle dynamics within the boundary layer. Deposition and particle flux simulations show that deposition affects noticeably only the smallest particles in the lowest part of the atmospheric boundary layer.
-
An analysis of cloud overlap based on high temporal and vertical resolution retrievals of cloud condensate from a suite of ground instruments is performed at a mid-latitude atmospheric observation facility. Two facets of overlap are investigated: cloud fraction overlap, expressed in terms of a parameter "<i>α</i>" indicating the relative contributions of maximum and random overlap, and overlap of horizontal distributions of condensate, expressed in terms of the correlation coefficient of condensate ranks. The degree of proximity to the random and maximum overlap assumptions is also expressed in terms of a decorrelation length, a convenient scalar parameter for overlap parameters assumed to decay exponentially with separation distance. Both cloud fraction overlap and condensate overlap show significant seasonal variations with a clear tendency for more maximum overlap in the summer months. More maximum overlap is also generally observed when the domain size used to define cloud fractions increases. These tendencies also exist for rank correlations, but are significantly weaker. Hitherto unexplored overlap parameter dependencies are investigated by analyzing mean parameter differences at fixed separation distance within different layers of the atmospheric column, and by searching for possible systematic relationships between alpha and rank correlation. We find that for the same separation distance the overlap parameters are significantly distinct in different atmospheric layers, and that random cloud fraction overlap is usually associated with more randomly overlapped condensate ranks.
-
Ice formation in the atmosphere by homogeneous and heterogeneous nucleation is one of the least understood processes in cloud microphysics and climate. Here we describe our investigation of the marine environment as a potential source of atmospheric IN by experimentally observing homogeneous ice nucleation from aqueous NaCl droplets and comparing against heterogeneous ice nucleation from aqueous NaCl droplets containing intact and fragmented diatoms. Homogeneous and heterogeneous ice nucleation are studied as a function of temperature and water activity, <i>a</i><sub>w</sub>. Additional analyses are presented on the dependence of diatom surface area and aqueous volume on heterogeneous freezing temperatures, ice nucleation rates, ω<sub>het</sub>, ice nucleation rate coefficients, <i>J</i><sub>het</sub>, and differential and cumulative ice nuclei spectra, <i>k(T)</i> and <i>K(T)</i>, respectively. Homogeneous freezing temperatures and corresponding nucleation rate coefficients are in agreement with the water activity based homogeneous ice nucleation theory within experimental and predictive uncertainties. Our results confirm, as predicted by classical nucleation theory, that a stochastic interpretation can be used to describe the homogeneous ice nucleation process. Heterogeneous ice nucleation initiated by intact and fragmented diatoms can be adequately represented by a modified water activity based ice nucleation theory. A horizontal shift in water activity, Δ<i>a</i><sub>w, het</sub> = 0.2303, of the ice melting curve can describe median heterogeneous freezing temperatures. Individual freezing temperatures showed no dependence on available diatom surface area and aqueous volume. Determined at median diatom freezing temperatures for <i>a</i><sub>w</sub> from 0.8 to 0.99, ω<sub>het</sub><u>~</u>0.11<sup>+0.06</sup><sub>−0.05</sub> s<sup>−1</sup>, <i>J</i><sub>het</sub><u>~</u>1.0<sup>+1.16</sup><sub>−0.61</sub>×10<sup>4</sup> cm<sup>−2</sup> s<sup>−1</sup>, and <i>K</i><u>~</u>6.2<sup>+3.5</sup><sub>−4.1</sub> ×10<sup>4</sup> cm<sup>−2</sup>. The experimentally derived ice nucleation rates and nuclei spectra allow us to estimate ice particle production which we subsequently use for a comparison with observed ice crystal concentrations typically found in cirrus and polar marine mixed-phase clouds. Differences in application of time-dependent and time-independent analyses to predict ice particle production are discussed.
-
We report on the first multi-year springtime measurements of PAN in the free troposphere over the US Pacific Northwest. The measurements were made at the summit of Mount Bachelor (43.979° N, 121.687° W; 2.7 km a.s.l.) by gas chromatography with electron capture detector during spring 2008, 2009 and 2010. This dataset provides an observational estimate of the month-to-month and springtime interannual variability of PAN mixing ratios in this region. Springtime seasonal mean (1 April–20 May) PAN mixing ratios at Mount Bachelor varied from 100 pptv to 152 pptv. The standard deviation of the three seasonal means was 28 pptv, 21 % of the springtime mean. We summarize the interannual variability in three factors expected to drive PAN variability: biomass burning, transport efficiency over the central and eastern Pacific, and transport temperature. <br></br> Zhang et al. (2008) used the GEOS-Chem global chemical transport model to show that rising Asian NO<sub>x</sub> emissions from 2000 to 2006 resulted in a relatively larger positive trend in PAN than O<sub>3</sub> over western North America. However the model results only considered monotonic changes in Asian emissions, whereas other factors, such as biomass burning, isoprene emissions or climate change can induce greater variability in the atmospheric concentrations and thus extend the time needed for trend detection. We combined the observed variability in PAN and O<sub>3</sub> at Mount Bachelor with a range of possible future trends in these species to determine the observational requirements to detect such trends. Though the relative increase in PAN is expected to be larger than that of O<sub>3</sub>, PAN is more variable. If PAN mixing ratios are currently increasing at a rate of 4 % per year due to rising Asian emissions, we would detect a trend with 13 years of measurements at a site like Mount Bachelor. If the corresponding trend in O<sub>3</sub> is 1 % per year, the trends in O<sub>3</sub> and PAN would be detected on approximately the same timescale.
-
Our ability to predict future climate change relies on our understanding of current and future CO<sub>2</sub> fluxes, particularly on a regional scale (100–1000 km). CO<sub>2</sub> regional sources and sinks are still poorly understood. Inverse transport modeling, a method often used to quantify these fluxes, relies on atmospheric CO<sub>2</sub> measurements. One of the main challenges for the transport models used in the inversions is to properly reproduce CO<sub>2</sub> vertical gradients between the boundary layer and the free troposphere, as these gradients impact on the partitioning of the calculated fluxes between the different model regions. Vertical CO<sub>2</sub> profiles are very well suited to assess the performances of the models. In this paper, we conduct a comparison between observed and modeled CO<sub>2</sub> profiles recorded during two CAATER campaigns that occurred in May 2001 and October 2002 over Western Europe, as described in a companion paper. We test different combinations between a global transport model (LMDZt), a mesoscale transport model (CHIMERE), and different sets of biospheric fluxes, all chosen with a diurnal cycle (CASA, SiB2 and ORCHIDEE). The vertical profile comparison shows that: 1) in most cases the influence of the biospheric flux is small but sometimes not negligible, ORCHIDEE giving the best results in the present study; 2) LMDZt is most of the time too diffuse, as it simulates a too high boundary layer height; 3) CHIMERE better reproduces the observed gradients between the boundary layer and the free troposphere, but is sometimes too variable and gives rise to incoherent structures. We conclude there is a need for more vertical profiles to conduct further studies to improve the parameterization of vertical transport in the models used for CO<sub>2</sub> flux inversions. <br><br> Furthermore, we use a modeling method to quantify CO<sub>2</sub> fluxes at the regional scale from a chosen observing point, coupling influence functions from the transport model LMDZt (that works quite well at the synoptic scale) with information on the space-time distribution of fluxes. This modeling method is compared to a dual tracer method (the so-called Radon method) for a case study on 25 May 2001 during which simultaneous well-correlated in situ CO<sub>2</sub> and Radon 222 measurements have been collected. Both methods give a similar result: a flux within the Radon 222 method uncertainty (35%), that is an atmospheric CO<sub>2</sub> sink of −4.2 to −4.4 gC m<sup>−2</sup> day<sup>−1</sup>. We have estimated the uncertainty of the modeling method to be at least 33% on average, and even more for specific individual events. This method allows the determination of the area that contributed to the CO<sub>2</sub> observed concentration. In our case, the observation point located at 1700 m a.s.l. in the north of France, is influenced by an area of 1500×700 km<sup>2</sup> that covers the Benelux region, part of Germany and western Poland. Furthermore, this method allows deconvolution between the different contributing fluxes. In this case study, the biospheric sink contributes 73% of the total flux, fossil fuel emissions for 27%, the oceanic flux being negligible. However, the uncertainties of the influence function method need to be better assessed. This could be possible by applying it to other cases where the calculated fluxes can be checked independently, for example at tall towers where simultaneous CO<sub>2</sub> and Radon 222 measurements can be conducted. The use of optimized fluxes (from atmospheric inversions) and of mesoscale models for atmospheric transport may also significantly reduce the uncertainties.
-
In this work we present an analysis of the occurrence of nucleation events during more than three years of measurements at two different rural altitude sites, the puy de Dôme research station (1465 m a.s.l.) and the Opme station (660 m a.s.l.), central France. The collected database is a unique combination of scanning mobility particle sizer (10–400 nm), air ion spectrometers (from 0.8 to 42 nm for NTP-conditions), and, neutral clusters and air ion spectrometers (from 0.8 to 42 nm for NTP-conditions) measurements at these two different altitudes nearly located research stations, from February 2007 to June 2010. The seasonality of the frequency of nucleation events was studied at the puy de Dôme station and maximum of events frequency was found during early spring and early autumn. During the measurement period, neither the particle formation rates (<span style="border-top: 1px solid #000; color: #000;"><i>J</i></span><sub>2</sub>= 1.382 ± 0.195 s<sup>−1</sup>) nor the growth rates (<span style="border-top: 1px solid #000; color: #000;">GR</span><sub>1.3−20 nm</sub> = 6.20 ± 0.12 nm h<sup>−1</sup>) differ from one site to the other on average. Hovewer, we found that, on 437 sampling days in common to the two sites, the nucleation frequency was higher at the puy de Dôme station (35.9 %, 157 days) than at the low elevation station of Opme (20.8 %, 91 days). LIDAR measurements and the evolution of the potential equivalent temperature revealed that the nucleation could be triggered either (i) within the whole low tropospheric column at the same time from the planetary boundary layer to the top of the interface layer (29.2 %, 47 events), (ii) above the planetary boundary layer upper limit (43.5 %, 70 events), and (iii) at low altitude and then transported, conserving dynamic and properties, at high altitude (24.8 %, 40 events). This is the first time that the vertical extent of nucleation can be studied over a long observational period, allowing for a rigorous statistical analysis of the occurrence of nucleation over the whole lower troposphere. This work highlights the fact that nucleation can occur over a large vertical extent, at least the whole low tropospheric column, and also the fact that it occurs twice as frequently as actually detected in the planetary boundary layer. The role of sulfuric acid and ions in the nucleation process was investigated at the altitude station and no correlation was found between nucleation events and the estimated sulfuric acid concentrations. However, the contribution of ion-induced nucleation was found to be relatively high (12.49 ± 2.03 % of the total nucleation rate).
-
The source contributions to carbonaceous PM<sub>2.5</sub> aerosol were investigated at a European background site at the edge of the Po Valley, in Northern Italy, during the period January–December 2007. Carbonaceous aerosol was described as the sum of 8 source components: primary (1) and secondary (2) biomass burning organic carbon, biomass burning elemental carbon (3), primary (4) and secondary (5) fossil organic carbon, fossil fuel burning elemental carbon (6), primary (7) and secondary (8) biogenic organic carbon. The mass concentration of each component was quantified using a set of macro tracers (organic carbon OC, elemental carbon EC, and levoglucosan), micro tracers (arabitol and mannitol), and <sup>14</sup>C measurements. This was the first time that <sup>14</sup>C measurements covered a full annual cycle with daily resolution. This set of 6 tracers, together with assumed uncertainty ranges of the ratios of OC-to-EC, and the reference fraction of modern carbon in the 8 source categories, provides strong constraints to the source contributions to carbonaceous aerosol. The uncertainty of contributions was assessed with a Quasi-Monte Carlo (QMC) method accounting for the variability of OC and EC emission factors, the uncertainty of reference fractions of modern carbon, and the measurement uncertainty. <br></br> During winter, biomass burning composed 64 % (±15 %) of the total carbon (TC) concentration, while in summer secondary biogenic OC accounted for 50 % (±16 %) of TC. The contribution of primary biogenic aerosol particles was negligible during the entire year. Moreover, aerosol associated with fossil sources represented 27 % (±16 %) and 41 % (±26 %) of TC in winter and summer, respectively. The contribution of secondary organic aerosol (SOA) to the organic mass (OM) was significant during the entire year. SOA accounted for 30 % (±16 %) and 85 % (±12 %) of OM during winter and summer, respectively. While the summer SOA was dominated by biogenic sources, winter SOA was mainly due to biomass burning and fossil sources. This indicates that the oxidation of semi-volatile and intermediate volatility organic compounds co-emitted with primary organics is a significant source of SOA, as suggested by recent model results and Aerosol Mass Spectrometer measurements. Comparison with previous global model simulations, indicates a strong underestimate of wintertime primary aerosol emissions in this region. The comparison of source apportionment results in different urban and rural areas showed that the sampling site was mainly affected by local aerosol sources during winter and regional air masses from the nearby Po Valley in summer. This observation was further confirmed by back-trajectory analysis applying the Potential Source Contribution Function method to identify potential source regions.
-
We report simultaneous measurements of volatile organic compound (VOC) mixing ratios including C<sub>6</sub> to C<sub>8</sub> aromatics, isoprene, monoterpenes, acetone and organic aerosol mass loadings at a rural location in southwestern Ontario, Canada by Proton-Transfer-Reaction Mass Spectrometry (PTR-MS) and Aerosol Mass Spectrometry (AMS), respectively. During the three-week-long Border Air Quality and Meteorology Study in June–July 2007, air was sampled from a range of sources, including aged air from the polluted US Midwest, direct outflow from Detroit 50 km away, and clean air with higher biogenic input. After normalization to the diurnal profile of CO, a long-lived tracer, diurnal analyses show clear photochemical loss of reactive aromatics and production of oxygenated VOCs and secondary organic aerosol (SOA) during the daytime. Biogenic VOC mixing ratios increase during the daytime in accord with their light- and temperature-dependent sources. Long-lived species, such as hydrocarbon-like organic aerosol and benzene show little to no photochemical reactivity on this timescale. From the normalized diurnal profiles of VOCs, an estimate of OH concentrations during the daytime, measured O<sub>3</sub> concentrations, and laboratory SOA yields, we calculate integrated local organic aerosol production amounts associated with each measured SOA precursor. Under the assumption that biogenic precursors are uniformly distributed across the southwestern Ontario location, we conclude that such precursors contribute significantly to the total amount of SOA formation, even during the period of Detroit outflow. The importance of aromatic precursors is more difficult to assess given that their sources are likely to be localized and thus of variable impact at the sampling location.
-
Noctilucent clouds (NLC) have been measured by the Rayleigh/Mie/Raman-lidar at the ALOMAR research facility in Northern Norway (69° N, 16° E). From 1997 to 2010 NLC were detected during more than 1850 h on 440 different days. Colocated MF-radar measurements and calculations with the Leibniz-Institute Middle Atmosphere (LIMA-) model are used to characterize the background atmosphere. Temperatures as well as horizontal winds at 83 km altitude show distinct differences during NLC observations compared to when NLC are absent. The seasonally averaged temperature is lower and the winds are stronger westward when NLC are detected. The wind separation is a robust feature as it shows up in measurements as well as in model results and it is consistent with the current understanding that lower temperatures support the existence of ice particles. For the whole 14-year data set there is no statistically significant relation between NLC occurrence and solar Lyman-α radiation. On the other hand NLC occurrence and temperatures at 83 km show a significant anti-correlation, which suggests that the thermal state plays a major role for the existence of ice particles and dominates the pure Lyman-α influence on water vapor during certain years. We find the seasonal mean NLC altitudes to be correlated to both Lyman-α radiation and temperature. NLC above ALOMAR are strongly influenced by atmospheric tides. The cloud water content varies by a factor of 2.8 over the diurnal cycle. Diurnal and semidiurnal amplitudes and phases show some pronounced year-to-year variations. In general, amplitudes as well as phases vary in a different manner. Amplitudes change by a factor of more than 3 and phases vary by up to 7 h. Such variability could impact long-term NLC observations which do not cover the full diurnal cycle.
-
We use observations of fire radiative power (FRP) from the Moderate Resolution Imaging Spectroradiometer~(MODIS) and tropospheric NO<sub>2</sub> column measurements from the Ozone Monitoring Instrument (OMI) to derive NO<sub>2</sub> wildfire emission coefficients (g MJ<sup>−1</sup>) for three land types over California and Nevada. Retrieved emission coefficients were 0.279±0.077, 0.342±0.053, and 0.696±0.088 g MJ<sup>−1</sup> NO<sub>2</sub> for forest, grass and shrub fuels, respectively. These emission coefficients reproduce ratios of emissions with fuel type reported previously using independent methods. However, the magnitude of these coefficients is lower than prior estimates. While it is possible that a negative bias in the OMI NO<sub>2</sub> retrieval over regions of active fire emissions is partly responsible, comparison with several other studies of fire emissions using satellite platforms indicates that current emission factors may overestimate the contributions of flaming combustion and underestimate the contributions of smoldering combustion to total fire emissions. Our results indicate that satellite data can provide an extensive characterization of the variability in fire NO<sub>x</sub> emissions; 67 % of the variability in emissions in this region can be accounted for using an FRP-based parameterization.