Items

Advanced search

Importance of transboundary transport of biomass burning emissions to regional air quality in Southeast Asia during a high fire event

Smoke from biomass and peat burning has a notable impact on ambient air quality and climate in the Southeast Asia (SEA) region. We modeled a large fire-induced haze episode in 2006 stemming mostly from Indonesia using the Weather Research and Forecasting model coupled with chemistry (WRF-Chem). We focused on the evolution of the fire plume composition and its interaction with the urbanized area of the city state of Singapore, and on comparisons of modeled and measured aerosol and carbon monoxide (CO) concentrations. Two simulations were run with WRF-Chem using the complex volatility basis set (VBS) scheme to reproduce primary and secondary aerosol evolution and concentration. The first simulation referred to as WRF-FIRE included anthropogenic, biogenic and biomass burning emissions from the Global Fire Emissions Database (GFED3) while the second simulation referred to as WRF-NOFIRE was run without emissions from biomass burning. To test model performance, we used three independent data sets for comparison including airborne measurements of particulate matter (PM) with a diameter of 10 μm or less (PM10) in Singapore, CO measurements in Sumatra, and aerosol optical depth (AOD) column observations from four satellite-based sensors. We found reasonable agreement between the model runs and both ground-based measurements of CO and PM10. The comparison with AOD was less favorable and indicated the model underestimated AOD, although the degree of mismatch varied between different satellite data sets. During our study period, forest and peat fires in Sumatra were the main cause of enhanced aerosol concentrations from regional transport over Singapore. Analysis of the biomass burning plume showed high concentrations of primary organic aerosols (POA) with values up to 600 μg m−3 over the fire locations. The concentration of POA remained quite stable within the plume between the main burning region and Singapore while the secondary organic aerosol (SOA) concentration slightly increased. However, the absolute concentrations of SOA (up to 20 μg m−3) were much lower than those from POA, indicating a minor role of SOA in these biomass burning plumes. Our results show that about 21% of the total mass loading of ambient PM10 during the July–October study period in Singapore was due to biomass and peat burning in Sumatra, but this contribution increased during high burning periods. In total, our model results indicated that during 35 days aerosol concentrations in Singapore were above the threshold of 50 μg m−3 day−1 indicating poor air quality. During 17 days this was due to fires, based on the difference between the simulations with and without fires. Local pollution in combination with recirculation of air masses was probably the main cause of poor air quality during the other 18 days, although fires from Sumatra and probably also from Kalimantan (Indonesian part of the island of Borneo) added to the enhanced PM10 concentrations. The model versus measurement comparisons highlighted that for our study period and region the GFED3 biomass burning aerosol emissions were more in line with observations than found in other studies. This indicates that care should be taken when using AOD to constrain emissions or estimate ground-level air quality. This study also shows the need for relatively high resolution modeling to accurately reproduce the advection of air masses necessary to quantify the impacts and feedbacks on regional air quality.
Erythemal ultraviolet irradiation trends in the Iberian Peninsula from 1950 to 2011

Erythemal ultraviolet (UVER) irradiation was reconstructed at nine Spanish locations, with series starting around 1950 in five of these places. Each series was checked by applying homogeneity tests in order to discard non-homogeneous series. Available series were used to create an averaged Iberian Peninsula UVER series. Results indicate that annual UVER irradiation in the Iberian Peninsula increased by 155 Jm−2 (6.5%) between 1950 and 2011 due to the decrease observed in atmospheric ozone rather than changes in aerosol and clouds. By contrast, annual UVER irradiation increased by 135 Jm−2 (5.6%) between 1985 and 2011, mainly due to changes in aerosol and clouds.
Integrating laboratory and field data to quantify the immersion freezing ice nucleation activity of mineral dust particles

Data from both laboratory studies and atmospheric measurements are used to develop an empirical parameterization for the immersion freezing activity of natural mineral dust particles. Measurements made with the Colorado State University (CSU) continuous flow diffusion chamber (CFDC) when processing mineral dust aerosols at a nominal 105% relative humidity with respect to water (RHw) are taken as a measure of the immersion freezing nucleation activity of particles. Ice active frozen fractions vs. temperature for dusts representative of Saharan and Asian desert sources were consistent with similar measurements in atmospheric dust plumes for a limited set of comparisons available. The parameterization developed follows the form of one suggested previously for atmospheric particles of non-specific composition in quantifying ice nucleating particle concentrations as functions of temperature and the total number concentration of particles larger than 0.5 μm diameter. Such an approach does not explicitly account for surface area and time dependencies for ice nucleation, but sufficiently encapsulates the activation properties for potential use in regional and global modeling simulations, and possible application in developing remote sensing retrievals for ice nucleating particles. A calibration factor is introduced to account for the apparent underestimate (by approximately 3, on average) of the immersion freezing fraction of mineral dust particles for CSU CFDC data processed at an RHw of 105% vs. maximum fractions active at higher RHw. Instrumental factors that affect activation behavior vs. RHw in CFDC instruments remain to be fully explored in future studies. Nevertheless, the use of this calibration factor is supported by comparison to ice activation data obtained for the same aerosols from Aerosol Interactions and Dynamics of the Atmosphere (AIDA) expansion chamber cloud parcel experiments. Further comparison of the new parameterization, including calibration correction, to predictions of the immersion freezing surface active site density parameterization for mineral dust particles, developed separately from AIDA experimental data alone, shows excellent agreement for data collected in a descent through a Saharan aerosol layer. These studies support the utility of laboratory measurements to obtain atmospherically relevant data on the ice nucleation properties of dust and other particle types, and suggest the suitability of considering all mineral dust as a single type of ice nucleating particle as a useful first-order approximation in numerical modeling investigations.
Exploring the chemical fate of the sulfate radical anion by reaction with sulfur dioxide in the gas phase

The gas phase reaction between SO4−(H2O)n and SO2, n = 0–2, is investigated using ab initio calculations and kinetic modelling. Structures of reactants, transition states and products are reported. Our calculations predict that the SO2SO4−(H2O)n cluster ion, which is formed upon SO2 and SO4−(H2O)n collision, can isomerize to SO3SO3−(H2O)n. The overall reaction is SO2 oxidation by the SO4−(H2O)n anionic cluster. The results show that SO4−(H2O)n is a good SO2 oxidant, especially at low relative humidity, with a reaction rate constant up to 1.5 × 10−10 cm3 molecule−1s−1. At high relative humidity, instead, the re-evaporation of SO2 from the SO2SO</sub4−(H2O)n cluster ion is favoured.
Understanding high wintertime ozone pollution events in an oil- and natural gas-producing region of the western US

Recent increases in oil and natural gas (NG) production throughout the western US have come with scientific and public interest in emission rates, air quality and climate impacts related to this industry. This study uses a regional-scale air quality model (WRF-Chem) to simulate high ozone (O3) episodes during the winter of 2013 over the Uinta Basin (UB) in northeastern Utah, which is densely populated by thousands of oil and NG wells. The high-resolution meteorological simulations are able qualitatively to reproduce the wintertime cold pool conditions that occurred in 2013, allowing the model to reproduce the observed multi-day buildup of atmospheric pollutants and the accompanying rapid photochemical ozone formation in the UB. Two different emission scenarios for the oil and NG sector were employed in this study. The first emission scenario (bottom-up) was based on the US Environmental Protection Agency (EPA) National Emission Inventory (NEI) (2011, version 1) for the oil and NG sector for the UB. The second emission scenario (top-down) was based on estimates of methane (CH4) emissions derived from in situ aircraft measurements and a regression analysis for multiple species relative to CH4 concentration measurements in the UB. Evaluation of the model results shows greater underestimates of CH4 and other volatile organic compounds (VOCs) in the simulation with the NEI-2011 inventory than in the case when the top-down emission scenario was used. Unlike VOCs, the NEI-2011 inventory significantly overestimates the emissions of nitrogen oxides (NOx), while the top-down emission scenario results in a moderate negative bias. The model simulation using the top-down emission case captures the buildup and afternoon peaks observed during high O3 episodes. In contrast, the simulation using the bottom-up inventory is not able to reproduce any of the observed high O3 concentrations in the UB. Simple emission reduction scenarios show that O3 production is VOC sensitive and NOx insensitive within the UB. The model results show a disproportionate contribution of aromatic VOCs to O3 formation relative to all other VOC emissions. The model analysis reveals that the major factors driving high wintertime O3 in the UB are shallow boundary layers with light winds, high emissions of VOCs from oil and NG operations compared to NOx emissions, enhancement of photolysis fluxes and reduction of O3 loss from deposition due to snow cover.
Variability of BVOC emissions from a Mediterranean mixed forest in southern France with a focus on Quercus pubescens

We aimed at quantifying biogenic volatile organic compound (BVOC) emissions in June from three Mediterranean species located at the O3HP site (southern France): Quercus pubescens, Acer monspessulanum and C. coggygria (for isoprene only). As Q. pubescens was shown to be the main BVOC emitter with isoprene representing ≈ 99% of the carbon emitted as BVOC, we mainly focused on this species. C. coggygria was found to be a non-isoprene emitter (no other BVOCs were investigated). To fully understand both the canopy effect on Q. pubescens isoprene emissions and the inter-individual variability (tree to tree and within canopy), diurnal variations of isoprene were investigated from nine branches (seven branches located to the top of canopy at ≈ 4 m above ground level (a.g.l.), and two inside the canopy at ≈ 2 m a.g.l.). The Q. pubescens daily mean isoprene emission rate (ERd) fluctuated between 23 and 98 μgC gDM−1 h−1. Q. pubescens daily mean net assimilation (Pn) ranged between 5.4 and 13.8, and 2.8 and 6.4 μmol CO2 m−2 s−1 for sunlit and shaded branches respectively. Both ERd and isoprene emission factors (Is), assessed according to Guenther et al. (1993) algorithm, varied by a factor of 4.3 among the sunlit branches. While sunlit branches ERd was clearly higher than for shaded branches, there was a non-significant variability of Is (59 to 77 μgC gDM−1 h−1). Diurnal variations of isoprene emission rates (ERs) for sunlit branches were also investigated. ERs were detected at dawn 2 h after Pn became positive and were mostly exponentially dependent on Pn. Diurnal variations of ERs were not equally well described throughout the day by temperature (CT) and light (CL) parameters according to G93 algorithm. Temperature had more impact than photosynthetically active radiation (PAR) on the morning emissions increase, and ER was no longer correlated to CL × CT between solar noon (maximum ER) and mid-afternoon, possibly due to thermal stress of the plant. A comparison between measured and calculated emissions using two isoprene algorithms (G93 and MEGAN – Model of Emissions of Gases and Aerosols from Nature) highlighted the importance of isoprene emission factor Is value used, and some weakness in assessing isoprene emissions under Mediterranean climate conditions (drought) with current isoprene models.
Improved AIOMFAC model parameterisation of the temperature dependence of activity coefficients for aqueous organic mixtures

This study presents a new, improved parameterisation of the temperature dependence of activity coefficients in the AIOMFAC (Aerosol Inorganic–Organic Mixtures Functional groups Activity Coefficients) model applicable for aqueous as well as water-free organic solutions. For electrolyte-free organic and organic–water mixtures the AIOMFAC model uses a group-contribution approach based on UNIFAC (UNIversal quasi-chemical Functional-group Activity Coefficients). This group-contribution approach explicitly accounts for interactions among organic functional groups and between organic functional groups and water. The previous AIOMFAC version uses a simple parameterisation of the temperature dependence of activity coefficients, aimed to be applicable in the temperature range from ~ 275 to ~ 400 K. With the goal to improve the description of a wide variety of organic compounds found in atmospheric aerosols, we extend the AIOMFAC parameterisation for the functional groups carboxyl, hydroxyl, ketone, aldehyde, ether, ester, alkyl, aromatic carbon-alcohol, and aromatic hydrocarbon to atmospherically relevant low temperatures. To this end we introduce a new parameterisation for the temperature dependence. The improved temperature dependence parameterisation is derived from classical thermodynamic theory by describing effects from changes in molar enthalpy and heat capacity of a multi-component system. Thermodynamic equilibrium data of aqueous organic and water-free organic mixtures from the literature are carefully assessed and complemented with new measurements to establish a comprehensive database, covering a wide temperature range (~ 190 to ~ 440 K) for many of the functional group combinations considered. Different experimental data types and their processing for the estimation of AIOMFAC model parameters are discussed. The new AIOMFAC parameterisation for the temperature dependence of activity coefficients from low to high temperatures shows an overall improvement of 28% in comparison to the previous model version, when both versions are compared to our database of experimentally determined activity coefficients and related thermodynamic data. When comparing the previous and new AIOMFAC model parameterisations to the subsets of experimental data with all temperatures below 274 K or all temperatures above 322 K (i.e. outside a 25 K margin of the reference temperature of 298 K), applying the new parameterisation leads to 37% improvement in each of the two temperature ranges considered. The new parameterisation of AIOMFAC agrees well with a large number of experimental data sets. Larger model–measurement discrepancies were found particularly for some of the systems containing multi-functional organic compounds. The affected systems were typically also poorly represented at room temperature and further improvements will be necessary to achieve better performance of AIOMFAC in these cases (assuming the experimental data are reliable). The performance of the AIOMFAC parameterisation is typically better for systems containing relatively small organic compounds and larger deviations may occur in mixtures where molecules of high structural complexity such as highly oxygenated compounds or molecules of high molecular mass (e.g. oligomers) prevail. Nevertheless, the new parameterisation enables the calculation of activity coefficients for a wide variety of different aqueous/water-free organic solutions down to the low temperatures present in the upper troposphere.
Modeling dust as component minerals in the Community Atmosphere Model: development of framework and impact on radiative forcing

The mineralogy of desert dust is important due to its effect on radiation, clouds and biogeochemical cycling of trace nutrients. This study presents the simulation of dust radiative forcing as a function of both mineral composition and size at the global scale, using mineral soil maps for estimating emissions. Externally mixed mineral aerosols in the bulk aerosol module in the Community Atmosphere Model version 4 (CAM4) and internally mixed mineral aerosols in the modal aerosol module in the Community Atmosphere Model version 5.1 (CAM5) embedded in the Community Earth System Model version 1.0.5 (CESM) are speciated into common mineral components in place of total dust. The simulations with mineralogy are compared to available observations of mineral atmospheric distribution and deposition along with observations of clear-sky radiative forcing efficiency. Based on these simulations, we estimate the all-sky direct radiative forcing at the top of the atmosphere as + 0.05 Wm−2 for both CAM4 and CAM5 simulations with mineralogy. We compare this to the radiative forcing from simulations of dust in release versions of CAM4 and CAM5 (+0.08 and +0.17 Wm−2) and of dust with optimized optical properties, wet scavenging and particle size distribution in CAM4 and CAM5, −0.05 and −0.17 Wm−2, respectively. The ability to correctly include the mineralogy of dust in climate models is hindered by its spatial and temporal variability as well as insufficient global in situ observations, incomplete and uncertain source mineralogies and the uncertainties associated with data retrieved from remote sensing methods.
A global survey of cloud overlap based on CALIPSO and CloudSat measurements

Using 2B-CLDCLASS-LIDAR (radar–lidar) cloud classification and 2B-FLXHR-LIDAR radiation products from CloudSat over 4 years, this study evaluates the co-occurrence frequencies of different cloud types, analyzes their along-track horizontal scales and cloud radiative effects (CREs), and utilizes the vertical distributions of cloud types to evaluate cloud-overlap assumptions. The statistical results show that high clouds, altostratus (As), altocumulus (Ac) and cumulus (Cu) tend to coexist with other cloud types. However, stratus (St) (or stratocumulus, Sc), nimbostratus (Ns) and convective clouds are much more likely to exhibit individual features than other cloud types. On average, altostratus-over-stratus/stratocumulus cloud systems have a maximum horizontal scale of 17.4 km, with a standard deviation of 23.5 km. Altocumulus-over-cumulus cloud types have a minimum scale of 2.8 km, with a standard deviation of 3.1 km. By considering the weight of each multilayered cloud type, we find that the global mean instantaneous net CREs of multilayered cloud systems during the daytime are approximately −41.3 and −50.2 W m−2, which account for 40.1 and 42.3% of the global mean total net CREs at the top of the atmosphere (TOA) and at the surface, respectively. The radiative contributions of high-over-altocumulus and high-over-stratus/stratocumulus (or cumulus) in the all multilayered cloud systems are dominant due to their frequency. Considering the overlap of cloud types, the cloud fraction based on the random overlap assumption is underestimated over vast oceans, except in the west-central Pacific Ocean warm pool. Obvious overestimations mainly occur over tropical and subtropical land masses. In view of a lower degree of overlap than that predicted by the random overlap assumption to occur over the vast ocean, particularly poleward of 40° S, the study therefore suggests that a linear combination of minimum and random overlap assumptions may further improve the predictions of actual cloud fractions for multilayered cloud types (e.g., As + St/Sc and Ac + St/Sc) over the Southern Ocean. The establishment of a statistical relationship between multilayered cloud types and the environmental conditions (e.g., atmospheric vertical motion, convective stability and wind shear) would be useful for parameterization design of cloud overlap in numerical models.
Multiday production of condensing organic aerosol mass in urban and forest outflow

Secondary organic aerosol (SOA) production in air masses containing either anthropogenic or biogenic (terpene-dominated) emissions is investigated using the explicit gas-phase chemical mechanism generator GECKO-A. Simulations show several-fold increases in SOA mass continuing for multiple days in the urban outflow, even as the initial air parcel is diluted into the regional atmosphere. The SOA mass increase in the forest outflow is more modest (~50%) and of shorter duration (1–2 days). The multiday production in the urban outflow stems from continuing oxidation of gas-phase precursors which persist in equilibrium with the particle phase, and can be attributed to multigenerational reaction products of both aromatics and alkanes, especially those with relatively low carbon numbers (C4–15). In particular we find large contributions from substituted maleic anhydrides and multi-substituted peroxide-bicyclic alkenes. The results show that the predicted production is a robust feature of our model even under changing atmospheric conditions and different vapor pressure schemes, and contradict the notion that SOA undergoes little mass production beyond a short initial formation period. The results imply that anthropogenic aerosol precursors could influence the chemical and radiative characteristics of the atmosphere over an extremely wide region, and that SOA measurements near precursor sources may routinely underestimate this influence.