Items

Advanced search

HFC-152a and HFC-134a emission estimates and characterization of CFCs, CFC replacements, and other halogenated solvents measured during the 2008 ARCTAS campaign (CARB phase) over the South Coast Air Basin of California

This work presents results from the NASA Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) study. Whole air samples were obtained on board research flights that flew over California during June 2008 and analyzed for selected volatile organic compounds, including several halogenated species. Samples collected over the South Coast Air Basin of California (SoCAB), which includes much of Los Angeles (LA) County, were compared with samples from inflow air masses over the Pacific Ocean. The levels of many halocarbon species were enhanced significantly over the SoCAB, including compounds regulated by the Montreal Protocol and subsequent amendments. Emissions estimates of HFC-152a (1,1-difluoroethane, CH3CHF2; 0.82 ± 0.11 Gg) and HFC-134a (1,1,1,2-tetrafluoroethane, CH2FCF3; 1.16 ± 0.22 Gg) in LA County for 2008 were obtained using the observed HFC:carbon monoxide (CO) enhancement ratio. Emission rates also were calculated for the SoCAB (1.60 ± 0.22 Gg yr−1 for HFC-152a and 2.12 ± 0.28 Gg yr−1 for HFC-134a) and then extrapolated to the United States (32 ± 4 Gg yr−1 for HFC-152a and 43 ± 6 Gg yr−1 for HFC-134a) using population data. In addition, emission rates of the two HFCs in LA County and SoCAB were calculated by a second method that utilizes air quality modeling. Emissions estimates obtained using both methods differ by less than 25% for the LA County and less than 45% for the SoCAB.
A statistical subgrid-scale algorithm for precipitation formation in stratiform clouds in the ECHAM5 single column model

Cloud properties are usually assumed to be homogeneous within the cloudy part of the grid-box, i.e. subgrid-scale inhomogeneities in cloud cover and/or microphysical properties are often neglected. However, precipitation formation is initiated by large particles. Thus mean values are not representative and could lead to a delayed onset of precipitation. For a more physical description of the subgrid-scale structure of clouds we introduce a new statistical sub-column algorithm to study the impact of cloud inhomogeneities on stratiform precipitation. Each model column is divided into N independent sub-columns with sub-boxes in each layer, which are completely clear or cloudy. The cloud cover is distributed over the sub-columns depending on the diagnosed cloud fraction. Mass and number concentrations of cloud droplets and ice crystals are distributed randomly over the cloudy sub-columns according to prescribed probability distributions. Shapes and standard deviations of the distributions are obtained from aircraft observations. We have implemented this sub-column algorithm into the ECHAM5 global climate model to take subgrid variability of cloud cover and microphysical properties into account. Simulations with the Single Column Model version of ECHAM5 were carried out for one period of the Mixed-Phase Polar Arctic Cloud Experiment (MPACE) campaign as well as for the Eastern Pacific Investigation of climate Processes (EPIC) campaign. Results with the new algorithm show an earlier onset of precipitation for the EPIC campaign and a higher conversion of liquid to ice for the MPACE campaign, which reduces the liquid water path in better agreement with the observations than the original version of the ECHAM5 model.
Measurement and simulation of the 16/17 April 2010 Eyjafjallajökull volcanic ash layer dispersion in the northern Alpine region

The spatial structure and the progression speed of the first ash layer from the Icelandic Eyjafjallajökull volcano which reached Germany on 16/17 April is investigated from remote sensing data and numerical simulations. The ceilometer network of the German Meteorological Service was able to follow the progression of the ash layer over the whole of Germany. This first ash layer turned out to be a rather shallow layer of only several hundreds of metres thickness which was oriented slantwise in the middle troposphere and which was brought downward by large-scale sinking motion over Southern Germany and the Alps. Special Raman lidar measurements, trajectory analyses and in-situ observations from mountain observatories helped to confirm the volcanic origin of the detected aerosol layer. Ultralight aircraft measurements permitted the detection of the arrival of a second major flush of volcanic material in Southern Germany. Numerical simulations with the Eulerian meso-scale model MCCM were able to reproduce the temporal and spatial structure of the ash layer. Comparisons of the model results with the ceilometer network data on 17 April and with the ultralight aircraft data on 19 April were satisfying. This is the first example of a model validation study from this ceilometer network data.
Evaluating the effects of microphysical complexity in idealised simulations of trade wind cumulus using the Factorial Method

The effect of microphysical and environmental factors on the development of precipitation in warm idealised cloud is explored using a kinematic modelling framework. A simple one-dimensional column model is used to drive a suite of microphysics schemes including a flexible multi-moment bulk scheme (including both single and dual moment cloud liquid water) and a state-of-the-art bin-resolved scheme with explicit treatments of liquid and aerosol. The Factorial Method is employed to quantify and compare the sensitivities of each scheme under a set of controlled conditions, in order to isolate the effect of additional microphysical complexity in terms of the impact on surface precipitation. At relatively low updraught speeds, the sensitivity of the bulk schemes was found to depend on the assumptions made with regards the treatment of droplet activation. It was possible to achieve a much closer agreement between the single and dual moment bulk schemes by tuning the specified droplet number concentration in the single moment scheme, suggesting that a diagnostic representation of droplet number may be an acceptable alternative to the more expensive prognostic option. However the effect of changes in CCN concentration were found to produce a relatively stronger effect on precipitation in the bulk schemes compared to the bin scheme; this is believed to be a consequence of differences in the treatment of drop growth by collision and coalescence. Collectively, these results demonstrate the usefulness of the Factorial Method as a model development tool for quantitatively comparing and contrasting the behaviour of microphysics schemes of differing levels of complexity within a specified parameter space.
Dry deposition of reactive nitrogen to European ecosystems: a comparison of inferential models across the NitroEurope network

Inferential models have long been used to determine pollutant dry deposition to ecosystems from measurements of air concentrations and as part of national and regional atmospheric chemistry and transport models, and yet models still suffer very large uncertainties. An inferential network of 55 sites throughout Europe for atmospheric reactive nitrogen (Nr) was established in 2007, providing ambient concentrations of gaseous NH3, NO2, HNO3 and HONO and aerosol NH4+ and NO3− as part of the NitroEurope Integrated Project. Network results providing modelled inorganic Nr dry deposition to the 55 monitoring sites are presented, using four existing dry deposition routines, revealing inter-model differences and providing ensemble average deposition estimates. Dry deposition is generally largest over forests in regions with large ambient NH3 concentrations, exceeding 30–40 kg N ha−1 yr−1 over parts of the Netherlands and Belgium, while some remote forests in Scandinavia receive less than 2 kg N ha−1 yr−1. Turbulent Nr deposition to short vegetation ecosystems is generally smaller than to forests due to reduced turbulent exchange, but also because NH3 inputs to fertilised, agricultural systems are limited by the presence of a substantial NH3 source in the vegetation, leading to periods of emission as well as deposition. Differences between models reach a factor 2–3 and are often greater than differences between monitoring sites. For soluble Nr gases such as NH3 and HNO3, the non-stomatal pathways are responsible for most of the annual uptake over many surfaces, especially the non-agricultural land uses, but parameterisations of the sink strength vary considerably among models. For aerosol NH4+ and NO3− discrepancies between theoretical models and field flux measurements lead to much uncertainty in dry deposition rates for fine particles (0.1–0.5 μm). The validation of inferential models at the ecosystem scale is best achieved by comparison with direct long-term micrometeorological Nr flux measurements, but too few such datasets are available, especially for HNO3 and aerosol NH4+ and NO3−.
Sources of carbonaceous aerosol in the Amazon basin

The quantification of sources of carbonaceous aerosol is important to understand their atmospheric concentrations and regulating processes and to study possible effects on climate and air quality, in addition to develop mitigation strategies. In the framework of the European Integrated Project on Aerosol Cloud Climate Interactions (EUCAARI) fine (Dp < 2.5 μm) and coarse (2.5 μm < Dp <10 μm) aerosol particles were sampled from February to June (wet season) and from August to September (dry season) 2008 in the central Amazon basin. The mass of fine particles averaged 2.4 μg m−3 during the wet season and 4.2 μg m−3 during the dry season. The average coarse aerosol mass concentration during wet and dry periods was 7.9 and 7.6 μg m−3, respectively. The overall chemical composition of fine and coarse mass did not show any seasonality with the largest fraction of fine and coarse aerosol mass explained by organic carbon (OC); the average OC to mass ratio was 0.4 and 0.6 in fine and coarse aerosol modes, respectively. The mass absorbing cross section of soot was determined by comparison of elemental carbon and light absorption coefficient measurements and it was equal to 4.7 m2 g−1 at 637 nm. Carbon aerosol sources were identified by Positive Matrix Factorization (PMF) analysis of thermograms: 44% of fine total carbon mass was assigned to biomass burning, 43% to secondary organic aerosol (SOA), and 13% to volatile species that are difficult to apportion. In the coarse mode, primary biogenic aerosol particles (PBAP) dominated the carbonaceous aerosol mass. The results confirmed the importance of PBAP in forested areas. The source apportionment results were employed to evaluate the ability of global chemistry transport models to simulate carbonaceous aerosol sources in a regional tropical background site. The comparison showed an overestimation of elemental carbon (EC) by the TM5 model during the dry season and OC both during the dry and wet periods. The overestimation was likely due to the overestimation of biomass burning emission inventories and SOA production over tropical areas.
Cloud condensation nuclei in polluted air and biomass burning smoke near the mega-city Guangzhou, China – Part 2: Size-resolved aerosol chemical composition, diurnal cycles, and externally mixed weakly CCN-active soot particles

Size-resolved chemical composition, mixing state, and cloud condensation nucleus (CCN) activity of aerosol particles in polluted mega-city air and biomass burning smoke were measured during the PRIDE-PRD2006 campaign near Guangzhou, China, using an aerosol mass spectrometer (AMS), a volatility tandem differential mobility analyzer (VTDMA), and a continuous-flow CCN counter (DMT-CCNC). The size-dependence and temporal variations of the effective average hygroscopicity parameter for CCN-active particles (κa) could be parameterized as a function of organic and inorganic mass fractions (forg, finorg) determined by the AMS: κa,p=κorg·forg + κinorg·finorg. The characteristic κ values of organic and inorganic components were similar to those observed in other continental regions of the world: κorg≈0.1 and κinorg≈0.6. The campaign average κa values increased with particle size from ~0.25 at ~50 nm to ~0.4 at ~200 nm, while forg decreased with particle size. At ~50 nm, forg was on average 60% and increased to almost 100% during a biomass burning event. The VTDMA results and complementary aerosol optical data suggest that the large fractions of CCN-inactive particles observed at low supersaturations (up to 60% at S≤0.27%) were externally mixed weakly CCN-active soot particles with low volatility (diameter reduction <5% at 300 °C) and effective hygroscopicity parameters around κLV≈0.01. A proxy for the effective average hygroscopicity of the total ensemble of CCN-active particles including weakly CCN-active particles (κt) could be parameterized as a function of κa,p and the number fraction of low volatility particles determined by VTDMA (φLV): κt,p=κa,p−φLV·(κa,p−κLV). Based on κ values derived from AMS and VTDMA data, the observed CCN number concentrations (NCCN,S≈102–104 cm−3 at S = 0.068–0.47%) could be efficiently predicted from the measured particle number size distribution. The mean relative deviations between observed and predicted CCN concentrations were ~10% when using κt,p, and they increased to ~20% when using only κa,p. The mean relative deviations were not higher (~20%) when using an approximate continental average value of κ≈0.3, although the constant κ value cannot account for the observed temporal variations in particle composition and mixing state (diurnal cycles and biomass burning events). Overall, the results confirm that on a global and climate modeling scale an average value of κ≈0.3 can be used for approximate predictions of CCN number concentrations in continental boundary layer air when aerosol size distribution data are available without information about chemical composition. Bulk or size-resolved data on aerosol chemical composition enable improved CCN predictions resolving regional and temporal variations, but the composition data need to be highly accurate and complemented by information about particle mixing state to achieve high precision (relative deviations <20%).
Chemical composition and mixing-state of ice residuals sampled within mixed phase clouds

During an intensive campaign at the high alpine research station Jungfraujoch, Switzerland, in February/March 2006 ice particle residuals within mixed-phase clouds were sampled using the Ice-counterflow virtual impactor (Ice-CVI). Size, morphology, chemical composition, mineralogy and mixing state of the ice residual and the interstitial (i.e., non-activated) aerosol particles were analyzed by scanning and transmission electron microscopy. Ice nuclei (IN) were identified from the significant enrichment of particle groups in the ice residual (IR) samples relative to the interstitial aerosol. In terms of number lead-bearing particles are enriched by a factor of approximately 25, complex internal mixtures with silicates or metal oxides as major components by a factor of 11, and mixtures of secondary aerosol and carbonaceous material (C-O-S particles) by a factor of 2. Other particle groups (sulfates, sea salt, Ca-rich particles, external silicates) observed in the ice-residual samples cannot be assigned unambiguously as IN. Between 9 and 24% of all IR are Pb-bearing particles. Pb was found as major component in around 10% of these particles (PbO, PbCl2). In the other particles, Pb was found as some 100 nm sized agglomerates consisting of 3–8 nm sized primary particles (PbS, elemental Pb). C-O-S particles are present in the IR at an abundance of 17–27%. The soot component within these particles is strongly aged. Complex internal mixtures occur in the IR at an abundance of 9–15%. Most IN identified at the Jungfraujoch station are internal mixtures containing anthropogenic components (either as main or minor constituent), and it is concluded that admixture of the anthropogenic component is responsible for the increased IN efficiency within mixed phase clouds. The mixing state appears to be a key parameter for the ice nucleation behaviour that cannot be predicted from the sole knowledge of the main component of an individual particle.
Origin and radiative forcing of black carbon transported to the Himalayas and Tibetan Plateau

The remote and high elevation regions of central Asia are influenced by black carbon (BC) emissions from a variety of locations. BC deposition contributes to melting of glaciers and questions exist, of both scientific and policy interest, as to the origin of the BC reaching the glaciers. We use the adjoint of the GEOS-Chem model to identify the location from which BC arriving at a variety of locations in the Himalayas and Tibetan Plateau originates. We then calculate its direct and snow-albedo radiative forcing. We analyze the seasonal variation in the origin of BC using an adjoint sensitivity analysis, which provides a detailed map of the location of emissions that directly contribute to black carbon concentrations at receptor locations. We find that emissions from northern India and central China contribute the majority of BC to the Himalayas, although the precise location varies with season. The Tibetan Plateau receives most BC from western and central China, as well as from India, Nepal, the Middle East, Pakistan and other countries. The magnitude of contribution from each region varies with season and receptor location. We find that sources as varied as African biomass burning and Middle Eastern fossil fuel combustion can significantly contribute to the BC reaching the Himalayas and Tibetan Plateau. We compute radiative forcing in the snow-covered regions and find the forcing due to the BC induced snow-albedo effect to vary from 5–15 W m−2 within the region, an order of magnitude larger than radiative forcing due to the direct effect, and with significant seasonal variation in the northern Tibetan Plateau. Radiative forcing from reduced snow albedo likely accelerates glacier melting. Our analysis may help inform mitigation efforts to slow the rate of glacial melt by identifying regions that make the largest contributions to BC deposition in the Himalayas and Tibetan Plateau.
Evaluating a 3-D transport model of atmospheric CO2 using ground-based, aircraft, and space-borne data

We evaluate the GEOS-Chem atmospheric transport model (v8-02-01) of CO2 over 2003–2006, driven by GEOS-4 and GEOS-5 meteorology from the NASA Goddard Global Modeling and Assimilation Office, using surface, aircraft and space-borne concentration measurements of CO2. We use an established ensemble Kalman Filter to estimate a posteriori biospheric+biomass burning (BS + BB) and oceanic (OC) CO2 fluxes from 22 geographical regions, following the TransCom-3 protocol, using boundary layer CO2 data from a subset of GLOBALVIEW surface sites. Global annual net BS + BB + OC CO2 fluxes over 2004–2006 for GEOS-4 (GEOS-5) meteorology are −4.4 ± 0.9 (−4.2 ± 0.9), −3.9 ± 0.9 (−4.5 ± 0.9), and −5.2 ± 0.9 (−4.9 ± 0.9) PgC yr−1, respectively. After taking into account anthropogenic fossil fuel and bio-fuel emissions, the global annual net CO2 emissions for 2004–2006 are estimated to be 4.0 ± 0.9 (4.2 ± 0.9), 4.8 ± 0.9 (4.2 ± 0.9), and 3.8 ± 0.9 (4.1 ± 0.9) PgC yr−1, respectively. The estimated 3-yr total net emission for GEOS-4 (GEOS-5) meteorology is equal to 12.5 (12.4) PgC, agreeing with other recent top-down estimates (12–13 PgC). The regional a posteriori fluxes are broadly consistent in the sign and magnitude of the TransCom-3 study for 1992–1996, but we find larger net sinks over northern and southern continents. We find large departures from our a priori over Europe during summer 2003, over temperate Eurasia during 2004, and over North America during 2005, reflecting an incomplete description of terrestrial carbon dynamics. We find GEOS-4 (GEOS-5) a posteriori CO2 concentrations reproduce the observed surface trend of 1.91–2.43 ppm yr−1 (parts per million per year), depending on latitude, within 0.15 ppm yr−1 (0.2 ppm yr−1) and the seasonal cycle within 0.2 ppm (0.2 ppm) at all latitudes. We find the a posteriori model reproduces the aircraft vertical profile measurements of CO2 over North America and Siberia generally within 1.5 ppm in the free and upper troposphere but can be biased by up to 4–5 ppm in the boundary layer at the start and end of the growing season. The model has a small negative bias in the free troposphere CO2 trend (1.95–2.19 ppm yr−1) compared to AIRS data which has a trend of 2.21–2.63 ppm yr−1 during 2004–2006, consistent with surface data. Model CO2 concentrations in the upper troposphere, evaluated using CONTRAIL (Comprehensive Observation Network for TRace gases by AIrLiner) aircraft measurements, reproduce the magnitude and phase of the seasonal cycle of CO2 in both hemispheres. We generally find that the GEOS meteorology reproduces much of the observed tropospheric CO2 variability, suggesting that these meteorological fields will help make significant progress in understanding carbon fluxes as more data become available.