-
Data collected at four ground-level sites are analyzed (1) to determine the surface cloud radiative effect (CRE) induced by cirrus clouds at regional scale for shortwave (CRE<sub>SW</sub>) and longwave (CRE<sub>LW</sub>) fluxes and (2) to derive the sensitivity of surface CRE<sub>SW</sub> to the cloud optical thickness (COT) modulated by the solar zenith angle and the atmospheric turbidity (noted CRE<sub>SW*</sub>) and the sensitivity of surface CRE<sub>LW</sub> to the infrared emissive power of cirrus cloud modulated by the water vapor content (noted CRE<sub>LW*</sub>). The average CRE<sub>SW*</sub> is −120 W m<sup>−2</sup> COT<sup>−1</sup> but it ranges from −80 to −140 m<sup>−2</sup> COT<sup>−1</sup> depending on the solar illumination with a residual variability ranges from +40 and −40 W m<sup>−2</sup> COT<sup>−1</sup> from pristine to turbid conditions, respectively. The CRE<sub>LW*</sub>, that corresponds to the infrared transmissivity of the atmosphere, ranges from 3% to 40% from dry to wet atmospheric conditions, respectively. The subvisible cirrus class (COT<0.03) over mid-latitude sites, that represents 20% of the population, induces a significant increase in surface LW irradiance at the 2–7 W m<sup>−2</sup> level. The semi-transparent cirrus class (0.03<COT<0.3), that represents 45% of the population, will affect the surface SW irradiance by −12 to −25 W m<sup>−2</sup>. Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) and the Atmospheric Infrared Sounder (AIRS) are used here to estimate the surface radiative effect at global scale. Global CRE estimations show very significant zonal and seasonal variability of each component of the CRE<sub>NET</sub>. CRE<sub>NET</sub> is 0.4 W m<sup>−2</sup> during winter/autumn for 15–75° N and 1 W m<sup>−2</sup> for 45–75° S whereas it is near −3 W m<sup>−2</sup> for 15° S–15° N. The summer period shows a cirrus cloud global cooling at all the latitudes except for 75–45° S with a quasi null effect and a peak at −3.6 W m<sup>−2</sup> for 15° S–45° N. The global average cumulative CRE is −2.8, 1.7 and −1.1 W m<sup>−2</sup> for CRE<sub>SW</sub>, CRE<sub>LW</sub>, and CRE<sub>NET</sub>, respectively.
-
The atmospheric chemistry general circulation model EMAC (ECHAM5/MESSy atmospheric chemistry) is used to investigate the effect of height dependent emissions on tropospheric chemistry. In a sensitivity simulation, anthropogenic and biomass burning emissions are released in the lowest model layer. The resulting tracer distributions are compared to those of a former simulation applying height dependent emissions. Although the differences between the two simulations in the free troposphere are small (less than 5%), large differences are present in polluted regions at the surface, in particular for NO<sub>x</sub> (more than 100%), CO (up to 30%) and non-methane hydrocarbons (up to 30%), whereas for OH the differences at the same locations are somewhat lower (15%). Global ozone formation is virtually unaffected by the choice of the vertical distribution of emissions. Nevertheless, local ozone changes can be up to 30%. Model results of both simulations are further compared to observations from field campaigns and to data from measurement stations.
-
The vertical distribution of ice crystals in Arctic boundary-layer mixed-phase (ABM) clouds was investigated by airborne remote-sensing and in situ measurements during the Arctic Study of Tropospheric Aerosol, Clouds and Radiation (ASTAR) campaign in March and April 2007. Information on the spectral absorption of solar radiation by ice and liquid water cloud particles is derived from airborne measurements of solar spectral radiation reflected by these clouds. It is shown by calculation of the vertical weighting function of the measurements that the observed absorption of solar radiation is dominated by the upper cloud layers (50% within 200 m from cloud top). This vertical weighting function is shifted even closer to cloud top for wavelengths where absorption by ice dominates. On this basis an indicator of the vertical distribution of ice crystals in ABM clouds is designed. <br><br> Applying in situ measured microphysical properties, the cloud-top reflectivity was calculated by radiative transfer simulations and compared to the measurements. It is found that ice crystals near cloud top (mixed-phase cloud top layer) are necessary to reproduce the measurements at wavelengths where absorption by ice dominates. The observation of backscatter glories on the cloud top generally indicating liquid water droplets does not contradict the postulated presence of ice crystals. Radiative transfer simulations reproduce the observed glories even if the cloud top layer is of mixed-phase character.
-
In this paper we analyze distributions of water vapour isotopologues in the troposphere using infrared spectra recorded by the Infrared Atmospheric Sounding Interferometer (IASI), which operates onboard the Metop satellite in nadir geometry. The simultaneous uncorrelated retrievals of H<sub>2</sub><sup>16</sup>O and HDO are performed on radiance measurements using a line-by-line radiative transfer model and an inversion procedure based on the Optimal Estimation Method (OEM). The characterizations of the retrieved products in terms of vertical sensitivity and error budgets show that IASI measurements contain up to 6 independent pieces of information on the vertical distribution of H<sub>2</sub><sup>16</sup>O and up to 3.5 for HDO from the surface up to the upper troposphere (0–20 km). Although the purpose of the paper is not validation, a restricted comparison with sonde measurements shows that the retrieved H<sub>2</sub><sup>16</sup>O profiles capture the seasonal/latitudinal variations of the water content, with good accuracy in the lowest layer but with larger uncertainties higher in the free and upper troposphere. Our results then demonstrate the ability of the IASI instrument to monitor atmospheric isotopologic water vapour distributions and to provide information on the partitioning of HDO as compared to H<sub>2</sub><sup>16</sup>O. The derivation of the δ<sub><i>D</i></sub> is challenging and associated with large errors in the uncorrelated retrieval approach chosen here. As a result averaging on the vertical to produce a column-averaged δ<sub><i>D</i></sub> is required to produce meaningful results for geophysical interpretation. As a case study, we analyse concentration distributions and spatio-temporal variations of H<sub>2</sub><sup>16</sup>O and δ<sub><i>D</i></sub> during the October 2007 Krosa super-typhoon over South-East Asia. We show that individual δ<sub><i>D</i></sub> have uncertainties of 37‰ for the vertically averaged values. Using the latter, we suggest that the typhoon produces a so-called amount-effect, where the δ<sub><i>D</i></sub> is negatively correlated to the water amounts as a result of intense depletion of the deuterated species.
-
Number fluxes of particles with diameter larger than 10 nm were measured with the eddy covariance method over the Amazon rain forest during the wet season as part of the LBA (The Large Scale Biosphere Atmosphere Experiment in Amazonia) campaign 2008. The primary goal was to investigate whether sources or sinks dominate the aerosol number flux in the tropical rain forest-atmosphere system. <br><br> During the measurement campaign, from 12 March to 18 May, 60% of the particle fluxes pointed downward, which is a similar fraction to what has been observed over boreal forests. The net deposition flux prevailed even in the absolute cleanest atmospheric conditions during the campaign and therefore cannot be explained only by deposition of anthropogenic particles. The particle transfer velocity <i>v</i><sub><i>t</i></sub> increased with increasing friction velocity and the relation is described by the equation <i>v</i><sub><i>t</i></sub> = 2.4×10<sup>−3</sup>×<i>u</i><sub>*</sub> where <i>u</i><sub>*</sub> is the friction velocity. <br><br> Upward particle fluxes often appeared in the morning hours and seem to a large extent to be an effect of entrainment fluxes into a growing mixed layer rather than primary aerosol emission. In general, the number source of primary aerosol particles within the footprint area of the measurements was small, possibly because the measured particle number fluxes reflect mostly particles less than approximately 200 nm. This is an indication that the contribution of primary biogenic aerosol particles to the aerosol population in the Amazon boundary layer may be low in terms of number concentrations. However, the possibility of horizontal variations in primary aerosol emission over the Amazon rain forest cannot be ruled out.
-
Interactions with water are crucial for the properties, transformation and climate effects of atmospheric aerosols. Here we present a conceptual framework for the interaction of amorphous aerosol particles with water vapor, outlining characteristic features and differences in comparison to crystalline particles. We used a hygroscopicity tandem differential mobility analyzer (H-TDMA) to characterize the hydration and dehydration of crystalline ammonium sulfate, amorphous oxalic acid and amorphous levoglucosan particles (diameter ~100 nm, relative humidity 5–95% at 298 K). The experimental data and accompanying Köhler model calculations provide new insights into particle microstructure, surface adsorption, bulk absorption, phase transitions and hygroscopic growth. The results of these and related investigations lead to the following conclusions: <br><br> (1) Many organic substances, including carboxylic acids, carbohydrates and proteins, tend to form amorphous rather than crystalline phases upon drying of aqueous solution droplets. Depending on viscosity and microstructure, the amorphous phases can be classified as glasses, rubbers, gels or viscous liquids. <br><br> (2) Amorphous organic substances tend to absorb water vapor and undergo gradual deliquescence and hygroscopic growth at lower relative humidity than their crystalline counterparts. <br><br> (3) In the course of hydration and dehydration, certain organic substances can form rubber- or gel-like structures (supramolecular networks) and undergo transitions between swollen and collapsed network structures. <br><br> (4) Organic gels or (semi-)solid amorphous shells (glassy, rubbery, ultra-viscous) with low molecular diffusivity can kinetically limit the uptake and release of water and may influence the hygroscopic growth and activation of aerosol particles as cloud condensation nuclei (CCN) and ice nuclei (IN). Moreover, (semi-)solid amorphous phases may influence the uptake of gaseous photo-oxidants and the chemical transformation and aging of atmospheric aerosols. <br><br> (5) The shape and porosity of amorphous and crystalline particles formed upon dehydration of aqueous solution droplets depend on chemical composition and drying conditions. The apparent volume void fractions of particles with highly porous structures can range up to ~50% or more (xerogels, aerogels). <br><br> (6) For efficient description of water uptake and phase transitions of aerosol particles, we propose not to limit the terms deliquescence and efflorescence to equilibrium phase transitions of crystalline substances. Instead we propose generalized definitions according to which amorphous and crystalline components can undergo gradual or prompt, partial or full deliquescence or efflorescence. <br><br> We suggest that (semi-)solid amorphous phases may be important not only in the upper atmosphere as suggested in recent studies of glass formation at low temperatures. Depending on relative humidity, (semi-)solid phases and moisture-induced glass transitions may also play a role in gas-particle interactions at ambient temperatures in the lower atmosphere.
-
In this paper we investigate some basic properties of the multi-model ensemble systems, which can be deduced from a general characteristic of statistical distributions of the ensemble members with the help of mathematical tools. In particular we show how to find optimal linear combination of model results, which minimizes the mean square error both in the case of uncorrelated and correlated models. By proving basic estimations we try to deduce general properties describing multi-model ensemble systems. We show also how mathematical formalism can be used for investigation of the characteristics of such systems.
-
The Antarctic and sub-Antarctic methanesulphonic acid (MSA) to non sea salt sulphate (nssSO<sub>4</sub>) ratio is simulated with the Laboratoire de Météorologie Dynamique Atmospheric General Circulation Model including an atmospheric sulphur chemistry module. Spatial variations of the MSA/nssSO<sub>4</sub> ratio in different regions have been suggested to be mostly dependent on temperature or sulphur source contributions. Its past variations in ice cores have been interpreted as related to the DMS precursor source location. Our model results are compared with available field measurements in the Antarctic and sub-Antarctic regions. This suggests that the MSA/nssSO<sub>4</sub> ratio in the extra-tropical south hemisphere is mostly dependent on the relative importance of various DMS oxidation pathways. In order to evaluate the effect of a rapid conversion of dimethyl sulphoxide (DMSO) into MSA, not implemented in the model, the MSA+DMSO to nssSO<sub>4</sub> ratio is also discussed. Using this modified ratio, the model mostly captures the seasonal variations of MSA/nssSO<sub>4</sub> at mid and high-southern latitudes. In addition, the model qualitatively reproduces the bell shaped meridional variations of the ratio, which is highly dependent on the adopted relative reaction rates for the DMS+OH addition and abstraction pathways, and on the assumed reaction products of the MSIA+OH reaction. MSA/nssSO<sub>4</sub> ratio in Antarctic snow is fairly well reproduced except at the most inland sites characterized with very low snow accumulation rates. Our results also suggest that atmospheric chemistry plays an important role in the observed decrease of the ratio in snow between coastal regions and central Antarctica. The still insufficient understanding of the DMS oxidation scheme limits our ability to model the MSA/nssSO<sub>4</sub> ratio. Specifically, reaction products of the MSIA+OH reaction should be better quantified, and the impact of a fast DMSO conversion to MSA in spring to fall over Antarctica should be evaluated. A better understanding of BrO source processes is needed in order to include DMS + BrO chemistry in global models. Completing the observations of DMS, BrO and MSA at Halley Bay with DMSO measurements would better constrain the role of BrO in DMS oxidation. Direct measurements of MSA and nssSO<sub>4</sub> dry deposition velocities on Antarctic snow would improve our ability to model MSA and nssSO<sub>4</sub> in ice cores.
-
Air quality in the outflow from urban centers affects millions of people, as well as, natural and managed ecosystems downwind. In locations where there are large sources of biogenic VOCs downwind of urban centers, the outflow is characterized by a high VOC reactivity due to biogenic emissions and low NO<sub>x</sub>. However most field and chamber studies have focused on limiting cases of high NO<sub>x</sub> or of near zero NO<sub>x</sub>. Recent measurements of a wide suite of VOCs, O<sub>3</sub> and meteorological parameters at several locations within the Sacramento urban plume have provided a detailed benchmark for testing our understanding of chemistry in a plume transitioning from high NO<sub>x</sub> to low NO<sub>x</sub> and high VOC reactivity. As an additional simplification, the strong mountain valley circulation in the region makes this urban plume a physical realization of a nearly idealized Lagrangian plume. Here, we describe a model of this plume. We use a Lagrangian model representing chemistry based on the Master Chemical Mechanism (MCM) v3.1 along with mixing and deposition. We discuss the effects of entrainment of background air, the branching ratio for the production of isoprene nitrates and the effects of soil NO<sub>x</sub> emissions on the composition of the evolving plume. The model predicts that after 2–3 h of chemical processing only 45% of the peroxynitrates (ΣPNs) are PAN and that most (69%) RONO<sub>2</sub> are secondary alkyl nitrate products of the reaction of OH with RONO<sub>2</sub>. We find the model is more consistent with the observations if: a) the yield of ΣPNs from large and multi-functional aldehydes is close to zero; and b) the reaction between OH and RONO<sub>2</sub> produces multifunctional nitrates as opposed to either HNO<sub>3</sub> or NO<sub>2</sub> as is typical in most currently adopted reaction mechanisms. Model results also show that adding NO<sub>x</sub> emissions throughout the transect increases the available NO<sub>x</sub> in the downwind regions, but modeled ozone concentrations were little affected by the increased NO<sub>x</sub>.
-
Formaldehyde (HCHO) columns have been retrieved from ground-based Fourier transform infrared (FTIR) campaign measurements in 2004 and 2007 and from UV-Visible MAX-DOAS measurements in 2004–2005 at the NDACC site of Réunion Island (21° S, 55° E). The FTIR and MAX-DOAS daily mean formaldehyde total columns are intercompared in their common measurement period, from August to October 2004. The ground-based data are also compared to correlative SCIAMACHY data. The comparisons account for the vertical sensitivity differences of the data sets, by including their respective averaging kernels. Complete error budgets are also presented. <br><br> The FTIR and MAX-DOAS daily mean total columns agree very well: no significant bias is observed and the standard deviation of the comparisons is only 8%. Both FTIR and MAX-DOAS HCHO total columns are in good agreement with SCIAMACHY values in the 2004–2005 period, with standard deviations of 21% and 31%, respectively. The same seasonal cycle is observed by the different instruments, with a minimum in austral winter and a maximum in February–March. <br><br> The FTIR and MAX-DOAS data are confronted with HCHO columns calculated by a global CTM, the IMAGES model. The model underestimates the HCHO columns by 23–29% in comparison with FTIR, and by 15% in comparison with DOAS. This bias might have multiple causes, including an underestimation of OH concentrations in the model (as indicated by a sensitivity study using prescribed OH fields) and/or an underestimated contribution of large-scale transport of HCHO precursors from Madagascar. The latter hypothesis is comforted by the large observed day-to-day variability of HCHO columns, and by the observation that the peak values of FTIR columns can often be associated with free tropospheric transport patterns from source regions over Madagascar to Réunion Island, according to simulations performed with the Lagrangian particle dispersion model FLEXPART.