ALBERT

Description: In this work, the performance of the Framework for Atmosphere-Canopy Exchange Modelling (FACEM: Pieterse et al., 2007) coupled to a Lagrangian atmospheric transport model is evaluated for carbon dioxide. Before incorporating FACEM into the Lagrangian COMET model (Vermeulen et al., 2006), its performance for the European domain is compared with the Simple Biosphere model (SiB3: Sellers et al., 1996). Overall, FACEM is well correlated to SiB3 (R2≥0.60), but shows less variability for regions with predominantly bare soil. There is no significant overall bias between the models except for the winter conditions and in general for the Iberian peninsula. When coupled to the COMET transport model, both biosphere models yield similar correlations (R2≥0.60) and bias relative to the 1-hourly concentration measurements for the year 2002, performed at three different sites in Europe; Cabauw (Netherlands), Hegyhatsal (Hungary) and Mace Head (Ireland). The overall results indicate that FACEM is comparable to SiB3 in terms of its applicability for atmospheric modelling studies.

Description: The reactivity of aromatic compounds is of great relevance to pure and applied chemical disciplines, yet existing methods for estimating gas-phase rate coefficients for their reactions with free radicals lack accuracy and universality. Here a novel approach is taken, whereby strong relationships between rate coefficients of aromatic hydrocarbons and a Randić-type topological index are investigated, optimized and developed into a method which requires no specialist software or computing power. Measured gas-phase rate coefficients for the reaction of aromatic hydrocarbons with OH radicals were correlated with a calculated Randić-type index, and optimized by including a term for side chain length. Although this method is exclusively for use with hydrocarbons, it is more diverse than any single existing methodology since it incorporates alkenylbenzenes into correlations, and can be extended towards other radical species such as O(3P) (and tentatively NO3, H and Cl). A comparison (with species common to both techniques) is made between the topological approach advocated here and a popular approach based on electrophilic subsituent constants, where it compares favourably. A modelling study was carried out to assess the impact of using estimated rate coefficients as opposed to measured data in an atmospheric model. The difference in model output was negligible for a range of NOx concentrations, which implies that this method has utility in complex chemical models. Strong relationships (e.g.~for OH, R2 = 0.96) between seemingly diverse compounds including benzene, multisubstituted benzenes with saturated, unsaturated, aliphatic and cyclic substitutions and the nonbenzenoid aromatic, azulene suggests that the Randić-type index presented here represents a new and effective way of describing aromatic reactivity, based on a quantitative structure-activity relationship (QSAR).

Description: The Lagrangian transport model COMET has been developed to evaluate emission estimates based on atmospheric concentration observations. This paper describes the model and its application in modelling the methane concentrations at the European stations Cabauw and Macehead. The COMET model captures in most cases both synoptic and diurnal variations of the concentrations as a function of time and in absolute size quite well. The explained variability by COMET of the mixed layer concentration for Cabauw varies from 50% to 84%; for all hourly observations in 2002 the explained variability is 71% with a RMSE of 112 ppb. The explained variability for Macehead is 48%. The most important model parameters were tested for their influence on model performance, but in general the model is not very sensitive to variations within acceptable limits. For a regionally and locally polluted continental site the COMET model shows only a small bias and a moderate random error, and therefore is considered to capture the influence of the sources on the concentration variations quite well. It is therefore concluded that inverse methods and more specifically the COMET model is suitable to be applied in deriving independent estimates of greenhouse gas emissions using Source-Receptor relationships.

Description: Inverse modeling techniques used to quantify surface carbon fluxes commonly assume that the uncertainty of fossil fuel CO2 (FFCO2) emissions is negligible and that intra-annual variations can be neglected. To investigate these assumptions, we analyzed the differences between four fossil fuel emission maps with spatial and temporal differences over Europe and their impact on the model simulated CO2 concentration. Large temporal flux variations characterize the hourly fields (~40% and ~80% for the seasonal and diurnal cycles, peak-to-peak) and annual country totals differ by 10% on average and up to 40% for some countries (i.e., The Netherlands). These emissions have been prescribed to seven different transport models, resulting in 28 different FFCO2 concentrations fields. The modeled FFCO2 concentration time series at surface sites using time-varying emissions show larger seasonal cycles (+2 ppm at the Hungarian tall tower (HUN)) and smaller diurnal cycles in summer (−1 ppm at HUN) than when using constant emissions. The concentration range spanned by all simulations varies between stations, and is generally larger in winter (up to ~10 ppm peak-to-peak at HUN) than in summer (~5 ppm). The contribution of transport model differences to the simulated concentration std-dev is 2–3 times larger than the contribution of emission differences only, at typical European sites used in global inversions. These contributions to the hourly (monthly) std-dev's amount to ~1.2 (0.8) ppm and ~0.4 (0.3) ppm for transport and emissions, respectively. First comparisons of the modeled concentrations with 14C-based fossil fuel CO2 observations show that the large transport differences still hamper a quantitative evaluation/validation of the emission inventories. Changes in the estimated monthly biosphere flux (Fbio) over Europe, using two inverse modeling approaches, are relatively small (less that 5%) while changes in annual Fbio (up to ~0.15 Gt C/yr) are only slightly smaller than the differences in annual emission totals and around 30% of the mean European ecosystem carbon sink. These results point to an urgent need to improve not only the transport models but also the assumed spatial and temporal distribution of fossil fuel emission maps.

Description: The isotopic composition of molecular hydrogen (H2) contains independent information for constraining the global H2 budget. To explore this, we have implemented hydrogen sources and sinks, including their stable isotopic composition and isotope fractionation constants, into the global chemistry transport model TM5. For the first time, a global model now includes a simplified but explicit isotope reaction scheme for the photochemical production of H2. We present a comparison of modelled results for the H2 mixing ratio and isotope composition with available measurements on seasonal to inter annual time scales for the years 2001–2007. The base model results agree well with observations for H2 mixing ratios. For δD[H2], modelled values are slightly lower than measurements. A detailed sensitivity study is performed to identify the most important parameters for modelling the isotopic composition of H2. The results show that on the global scale, the discrepancy between model and measurements can be closed by adjusting the default values of the isotope effects in deposition, photochemistry and the stratosphere-troposphere exchange within the known range of uncertainty. However, the available isotope data do not provide sufficient information to uniquely constrain the global isotope budget. Therefore, additional studies focussing on the isotopic composition near the tropopause and on the isotope effects in the photochemistry and deposition are recommended.

Description: Inverse modeling techniques used to quantify surface carbon fluxes commonly assume that the uncertainty of fossil fuel CO2 (FFCO2) emissions is negligible and that intra-annual variations can be neglected. To investigate these assumptions, we analyzed the differences between four fossil fuel emission inventories with spatial and temporal differences over Europe and their impact on the model simulated CO2 concentration. Large temporal flux variations characterize the hourly fields (~40 % and ~80 % for the seasonal and diurnal cycles, peak-to-peak) and annual country totals differ by 10 % on average and up to 40 % for some countries (i.e., the Netherlands). These emissions have been prescribed to seven different transport models, resulting in 28 different FFCO2 concentrations fields. The modeled FFCO2 concentration time series at surface sites using time-varying emissions show larger seasonal cycles (+2 ppm at the Hungarian tall tower (HUN)) and smaller diurnal cycles in summer (−1 ppm at HUN) than when using constant emissions. The concentration range spanned by all simulations varies between stations, and is generally larger in winter (up to ~10 ppm peak-to-peak at HUN) than in summer (~5 ppm). The contribution of transport model differences to the simulated concentration std-dev is 2–3 times larger than the contribution of emission differences only, at typical European sites used in global inversions. These contributions to the hourly (monthly) std-dev's amount to ~1.2 (0.8) ppm and ~0.4 (0.3) ppm for transport and emissions, respectively. First comparisons of the modeled concentrations with 14C-based fossil fuel CO2 observations show that the large transport differences still hamper a quantitative evaluation/validation of the emission inventories. Changes in the estimated monthly biosphere flux (Fbio) over Europe, using two inverse modeling approaches, are relatively small (less that 5 %) while changes in annual Fbio (up to ~0.15 % GtC yr−1) are only slightly smaller than the differences in annual emission totals and around 30 % of the mean European ecosystem carbon sink. These results point to an urgent need to improve not only the transport models but also the assumed spatial and temporal distribution of fossil fuel emission inventories.

Description: The isotopic composition of molecular hydrogen (H2) produced by photochemical oxidation of methane (CH4) and Volatile Organic Compounds (VOCs) is a key quantity in the global isotope budget of (H2). The many individual reaction steps involved complicate its investigation. Here we present a simplified structure-activity approach to assign isotope effects to the individual elementary reaction steps in the oxidation sequence of CH4 and some other VOCs. The approach builds on and extends the work by Gerst and Quay (2001) and Feilberg et al. (2007b). The description is generalized and allows the application, in principle, also to other compounds. The idea is that the C-H and C-D bonds – seen as reactive sites – have similar relative reaction probabilities in isotopically substituted, but otherwise identical molecules. The limitations of this approach are discussed for the reaction CH4+Cl. The same approach is applied to VOCs, which are important precursors of H2 that need to be included into models. Unfortunately, quantitative information on VOC isotope effects and source isotope signatures is very limited and the isotope scheme at this time is limited to a strongly parameterized statistical approach, which neglects kinetic isotope effects. Using these concepts we implement a full hydrogen isotope scheme in a chemical box model and carry out a sensitivity study to identify those reaction steps and conditions that are most critical for the isotope composition of the final H2 product. The reaction scheme is directly applicable in global chemistry models, which can thus include the isotope pathway of H2 produced from CH4 and VOCs in a consistent way.

Description: Despite the potential of isotope measurements to improve our understanding of the global atmospheric molecular hydrogen (H2) cycle, few H2 isotope data have been published so far. Now, within the EUROpean network for atmospheric HYDRogen Observations and Studies project (EUROHYDROS), weekly to monthly air samples from six locations in a global sampling network have been analysed for H2 mixing ratio (m(H2)) and the stable isotopic composition of the H2 (δ(D,H2), hereafter referred to as δD). The time series thus obtained now cover one to five years for all stations. This is the largest set of ground station observations of δD so far. Annual average δD values are higher at the Southern Hemisphere (SH) than at the Northern Hemisphere (NH) stations; the maximum is observed at Neumayer (Antarctica), and the minimum at the non-arctic NH stations. The maximum seasonal differences in δD range from ≈18 ‰ at Neumayer to ≈45 ‰ at Schauinsland (Southern Germany); in general, seasonal variability is largest at the NH stations. The timing of minima and maxima differs per station as well. In Alert (Arctic Canada), the variations in δD and m(H2) can be approximated as simple harmonic functions with a ≈5-month relative phase shift. This out-of-phase seasonal behaviour of δD and m(H2) can also be detected, but delayed and with a ≈6-month relative phase shift, at Mace Head and Cape Verde. However, no seasonal δD cycle could be observed at Schauinsland, which likely reflects the larger influence of local sources and sinks at this continental station. At the two SH stations, no seasonal cycle could be detected in the δD data. If it is assumed that the sink processes are the main drivers of the observed seasonality in m(H2) and δD on the NH, the relative seasonal variations can be used to estimate the relative sink strength of the two major sinks, deposition to soils and atmospheric oxidation by the hydroxyl (OH) radical. For the NH coastal and marine stations this analysis suggests that the relative contribution of soil uptake to the total annual H2 removal increases with latitude.