Items

Advanced search

A new model of the global biogeochemical cycle of carbonyl sulfide – Part 2: Use of carbonyl sulfide to constrain gross primary productivity in current vegetation models

Clear analogies between carbonyl sulfide (OCS) and carbon dioxide (CO2) diffusion pathways through leaves have been revealed by experimental studies, with plant uptake playing an important role for the atmospheric budget of both species. Here we use atmospheric OCS to evaluate the gross primary production (GPP) of three dynamic global vegetation models (Lund–Potsdam–Jena, LPJ; National Center for Atmospheric Research – Community Land Model 4, NCAR-CLM4; and Organising Carbon and Hydrology In Dynamic Ecosystems, ORCHIDEE). Vegetation uptake of OCS is modeled as a linear function of GPP and leaf relative uptake (LRU), the ratio of OCS to CO2 deposition velocities of plants. New parameterizations for the non-photosynthetic sinks (oxic soils, atmospheric oxidation) and biogenic sources (oceans and anoxic soils) of OCS are also provided. Despite new large oceanic emissions, global OCS budgets created with each vegetation model show exceeding sinks by several hundred Gg S yr−1. An inversion of the surface fluxes (optimization of a global scalar which accounts for flux uncertainties) led to balanced OCS global budgets, as atmospheric measurements suggest, mainly by drastic reduction (up to −50 %) in soil and vegetation uptakes. The amplitude of variations in atmospheric OCS mixing ratios is mainly dictated by the vegetation sink over the Northern Hemisphere. This allows for bias recognition in the GPP representations of the three selected models. The main bias patterns are (i) the terrestrial GPP of ORCHIDEE at high northern latitudes is currently overestimated, (ii) the seasonal variations of the GPP are out of phase in the NCAR-CLM4 model, showing a maximum carbon uptake too early in spring in the northernmost ecosystems, (iii) the overall amplitude of the seasonal variations of GPP in NCAR-CLM4 is too small, and (iv) for the LPJ model, the GPP is slightly out of phase for the northernmost ecosystems and the respiration fluxes might be too large in summer in the Northern Hemisphere. These results rely on the robustness of the OCS modeling framework and, in particular, the choice of the LRU values (assumed constant in time) and the parameterization of soil OCS uptake with small seasonal variations. Refined optimization with regional-scale and seasonally varying coefficients might help to test some of these hypothesis.
A comparison of measured HONO uptake and release with calculated source strengths in a heterogeneous forest environment

Vertical mixing ratio profiles of nitrous acid (HONO) were measured in a clearing and on the forest floor in a rural forest environment. For the forest floor, HONO was found to predominantly deposit, whereas for the clearing, net deposition dominated only during nighttime and net emissions were observed during daytime. For selected days, net fluxes of HONO were calculated from the measured profiles using the aerodynamic gradient method. The emission fluxes were in the range of 0.02 to 0.07 nmol m−2 s−1 and thus were in the lower range of previous observations. These fluxes were compared to the strengths of postulated HONO sources. Laboratory measurements of different soil samples from both sites revealed an upper limit for soil biogenic HONO emission fluxes of 0.025 nmol m−2 s−1. HONO formation by light-induced NO2 conversion was calculated to be below 0.03 nmol m−2 s−1 for the investigated days, which is comparable to the potential soil fluxes. Due to light saturation at low irradiance, this reaction pathway was largely found to be independent of light intensity, i.e. it was only dependent on ambient NO2. We used three different approaches based on measured leaf nitrate loadings for calculating HONO formation from HNO3 photolysis. While the first two approaches based on empirical HONO formation rates yielded values in the same order of magnitude as the estimated fluxes, the third approach based on available kinetic data of the postulated pathway failed to produce noticeable amounts of HONO. Estimates based on reported cross sections of adsorbed HNO3 indicate that the lifetime of adsorbed HNO3 was only about 15 min, which would imply a substantial renoxification. Although the photolysis of HNO3 was significantly enhanced at the surface, the subsequent light-induced conversion of the photolysis product NO2 did not produce considerable amounts of HONO. Consequently, this reaction might occur via an alternative mechanism. By explicitly calculating HONO formation based on available kinetic data and simple parameterizations, we showed that (a) for low NOx the light-induced conversion of NO2 on humic acids is already light saturated by the early morning, (b) HONO formation from photolysis of adsorbed HNO3 appears to proceed via an alternative mechanism and (c) estimates of HONO emissions from soil are very sensitive to mass transfer and acidic soils do not necessarily favour HONO emissions.
Microphysical processing of aerosol particles in orographic clouds

An explicit and detailed treatment of cloud-borne particles allowing for the consideration of aerosol cycling in clouds has been implemented into COSMO-Model, the regional weather forecast and climate model of the Consortium for Small-scale Modeling (COSMO). The effects of aerosol scavenging, cloud microphysical processing and regeneration upon cloud evaporation on the aerosol population and on subsequent cloud formation are investigated. For this, two-dimensional idealized simulations of moist flow over two bell-shaped mountains were carried out varying the treatment of aerosol scavenging and regeneration processes for a warm-phase and a mixed-phase orographic cloud. The results allowed us to identify different aerosol cycling mechanisms. In the simulated non-precipitating warm-phase cloud, aerosol mass is incorporated into cloud droplets by activation scavenging and released back to the atmosphere upon cloud droplet evaporation. In the mixed-phase cloud, a first cycle comprises cloud droplet activation and evaporation via the Wegener–Bergeron–Findeisen (WBF) process. A second cycle includes below-cloud scavenging by precipitating snow particles and snow sublimation and is connected to the first cycle via the riming process which transfers aerosol mass from cloud droplets to snowflakes. In the simulated mixed-phase cloud, only a negligible part of the total aerosol mass is incorporated into ice crystals. Sedimenting snowflakes reaching the surface remove aerosol mass from the atmosphere. The results show that aerosol processing and regeneration lead to a vertical redistribution of aerosol mass and number. Thereby, the processes impact the total aerosol number and mass and additionally alter the shape of the aerosol size distributions by enhancing the internally mixed/soluble Aitken and accumulation mode and generating coarse-mode particles. Concerning subsequent cloud formation at the second mountain, accounting for aerosol processing and regeneration increases the cloud droplet number concentration with possible implications for the ice crystal number concentration.
The influences of mass loading and rapid dilution of secondary organic aerosol on particle volatility

The thermally induced evaporation of secondary organic aerosol (SOA) has been characterized for SOA formed from the dark ozonolysis of α-pinene at initial mass concentrations ranging from 1 to 800 μg m−3. Temperature-dependent particle size distributions were measured using a thermodenuder and the resulting mass thermograms were compared between the SOA formed at the various SOA mass concentrations. Negligible differences were observed between the mass thermograms for SOA concentrations < 300 μg m−3. At higher SOA concentrations, the observed mass thermograms indicated the SOA was actually slightly less volatile than the SOA at lower concentrations; this is likely an artifact due to either saturation of the gas phase or to recondensation during cooling. The thermograms observed when the SOA was formed at high concentrations (> 380 μg m−3) and then rapidly isothermally diluted to low concentrations (1–20 μg m−3) were identical to those for the SOA that was initially formed at low concentrations. The experimental results were compared to a kinetic model that simulates particle evaporation upon heating in a thermodenuder for a given input volatility distribution and particle composition. Three cases were considered: (1) the SOA was composed of semi-volatile monomer species with a volatility distribution based on that derived previously from consideration of SOA growth experiments; (2) the initial SOA was composed almost entirely of non-volatile dimers that decompose upon heating into their semi-volatile monomer units, which can then evaporate; and (3) where a volatility distribution was derived by fitting the model to the observed mass thermograms. It was found that good agreement is obtained between model predictions and the observations when the particle composition is dominated by either compounds of low volatility or by dimers. These same models were used to simulate isothermal evaporation of the SOA and were found to be broadly consistent with literature observations that indicate that SOA evaporation occurs with multiple timescales. The use of the semi-volatile monomer volatility distribution fails to reproduce the observed evaporation. The presence of dimers and larger oligomers in secondary organic aerosol formed from products of the reaction of α-pinene and O3 has been well established in laboratory studies. However, the timescale and relative importance of the formation of oligomers or low-volatility compounds in the growth and evaporation of SOA has been debated. This study provides further support that low-volatility compounds and oligomers are formed in α-pinene + O3 in high abundances and suggests that their formation occurs rapidly upon particle formation.
Particulate organic nitrates observed in an oil and natural gas production region during wintertime

Organic nitrates in both gas and condensed (aerosol) phases were measured during the Uintah Basin Winter Ozone Study from January to February in 2012. A high degree of correlation between total aerosol volume at diameters less than 500 nm and the particulate organic nitrate concentration indicates that organic nitrates are a consistent, if not dominant, fraction of fine aerosol mass. In contrast, a similar correlation with sub-2.5 μm aerosol volume is weaker. The C : N atomic ratio inferred from field measurements of PM2.5 and particulate organic nitrate is 34 : 1. Calculations constrained by the observations indicate that both condensation of gas-phase nitrates and heterogeneous reactions of NO3 / N2O5 are responsible for introducing organic nitrate functionality into the aerosol and that the source molecules are alkanes. Extrapolating the results to urban aerosol suggests organic nitrate production from alkanes may be a major secondary organic aerosol source.
Ash iron mobilization through physicochemical processing in volcanic eruption plumes: a numerical modeling approach

It has been shown that volcanic ash fertilizes the Fe-limited areas of the surface ocean through releasing soluble iron. As ash iron is mostly insoluble upon the eruption, it is hypothesized that heterogeneous in-plume and in-cloud processing of the ash promote the iron solubilization. Direct evidences concerning such processes are, however, lacking. In this study, a 1-D numerical model is developed to simulate the physicochemical interactions of the gas–ash–aerosol in volcanic eruption plumes focusing on the iron mobilization processes at temperatures between 600 and 0 °C. Results show that sulfuric acid and water vapor condense at ~ 150 and ~ 50 °C on the ash surface, respectively. This liquid phase then efficiently scavenges the surrounding gases (> 95 % of HCl, 3–20 % of SO2 and 12–62 % of HF) forming an extremely acidic coating at the ash surface. The low pH conditions of the aqueous film promote acid-mediated dissolution of the Fe-bearing phases present in the ash material. We estimate that 0.1–33 % of the total iron available at the ash surface is dissolved in the aqueous phase before the freezing point is reached. The efficiency of dissolution is controlled by the halogen content of the erupted gas as well as the mineralogy of the iron at ash surface: elevated halogen concentrations and presence of Fe2+-carrying phases lead to the highest dissolution efficiency. Findings of this study are in agreement with the data obtained through leaching experiments.
Regional-scale transport of air pollutants: impacts of Southern California emissions on Phoenix ground-level ozone concentrations

In this study, WRF-Chem is utilized at high resolution (1.333 km grid spacing for the innermost domain) to investigate impacts of southern California anthropogenic emissions (SoCal) on Phoenix ground-level ozone concentrations ([O3]) for a pair of recent exceedance episodes. First, WRF-Chem control simulations, based on the US Environmental Protection Agency (EPA) 2005 National Emissions Inventories (NEI05), are conducted to evaluate model performance. Compared with surface observations of hourly ozone, CO, NOX, and wind fields, the control simulations reproduce observed variability well. Simulated [O3] are comparable with the previous studies in this region. Next, the relative contribution of SoCal and Arizona local anthropogenic emissions (AZ) to ozone exceedances within the Phoenix metropolitan area is investigated via a trio of sensitivity simulations: (1) SoCal emissions are excluded, with all other emissions as in Control; (2) AZ emissions are excluded with all other emissions as in Control; and (3) SoCal and AZ emissions are excluded (i.e., all anthropogenic emissions are eliminated) to account only for Biogenic emissions and lateral boundary inflow (BILB). Based on the USEPA NEI05, results for the selected events indicate the impacts of AZ emissions are dominant on daily maximum 8 h average (DMA8) [O3] in Phoenix. SoCal contributions to DMA8 [O3] for the Phoenix metropolitan area range from a few ppbv to over 30 ppbv (10–30 % relative to Control experiments). [O3] from SoCal and AZ emissions exhibit the expected diurnal characteristics that are determined by physical and photochemical processes, while BILB contributions to DMA8 [O3] in Phoenix also play a key role. Finally, ozone transport processes and pathways within the lower troposphere are investigated. During daytime, pollutants (mainly ozone) near the Southern California coasts are pumped into the planetary boundary-layer over the Southern California desert through the mountain chimney and pass channel effects, aiding eastward transport along the desert air basins in southern California and finally, northeastward along the lower Gila River basin in Arizona, thereby affecting Phoenix air quality during subsequent days. This study indicates that local emission controls in Phoenix need to be augmented with regional emission reductions to attain the federal ozone standard, especially if a more stringent standard is adopted in the future.
Impacts of an unknown daytime HONO source on the mixing ratio and budget of HONO, and hydroxyl, hydroperoxyl, and organic peroxy radicals, in the coastal regions of China

Many field experiments have found high nitrous acid (HONO) mixing ratios in both urban and rural areas during daytime, but these high daytime HONO mixing ratios cannot be explained well by gas-phase production, HONO emissions, and nighttime hydrolysis conversion of nitrogen dioxide (NO2) on aerosols, suggesting that an unknown daytime HONO source (Punknown) could exist. The formula Punknown ≈ 19.60[NO2] · J(NO2) was obtained using observed data from 13 field experiments across the globe. The three additional HONO sources (i.e., the Punknown, nighttime hydrolysis conversion of NO2 on aerosols, and HONO emissions) were coupled into the WRF-Chem model (Weather Research and Forecasting model coupled with Chemistry) to assess the Punknown impacts on the concentrations and budgets of HONO and peroxy (hydroxyl, hydroperoxyl, and organic peroxy) radicals (ROx) (= OH + HO2 + RO2) in the coastal regions of China. Results indicated that the additional HONO sources produced a significant improvement in HONO and OH simulations, particularly in the daytime. High daytime average Punknown values were found in the coastal regions of China, with a maximum of 2.5 ppb h−1 in the Beijing–Tianjin–Hebei region. The Punknown produced a 60–250 % increase of OH, HO2, and RO2 near the ground in the major cities of the coastal regions of China, and a 5–48 % increase of OH, HO2, and RO2 in the daytime meridional-mean mixing ratios within 1000 m above the ground. When the three additional HONO sources were included, the photolysis of HONO was the second most important source in the OH production rate in Beijing, Shanghai, and Guangzhou before 10:00 LST with a maximum of 3.72 ppb h−1 and a corresponding Punknown contribution of 3.06 ppb h−1 in Beijing, whereas the reaction of HO2 + NO (nitric oxide) was dominant after 10:00 LST with a maximum of 9.38 ppb h−1 and a corresponding Punknown contribution of 7.23 ppb h−1 in Beijing. The whole ROx cycle was accelerated by the three additional HONO sources, especially the Punknown. The daytime average OH production rate was enhanced by 0.67 due to the three additional HONO sources; [0.64], due to the Punknown, to 4.32 [3.86] ppb h−1, via the reaction of HO2 + NO, and by 0.49 [0.47] to 1.86 [1.86] ppb h−1, via the photolysis of HONO. The OH daytime average loss rate was enhanced by 0.58 [0.55] to 2.03 [1.92] ppb h−1, via the reaction of OH + NO2, and by 0.31 [0.28] to 1.78 [1.64] ppb h−1, via the reaction of OH + CO (carbon monoxide) in Beijing, Shanghai, and Guangzhou. Similarly, the three additional HONO sources produced an increase of 0.31 [0.28] (with a corresponding Punknown contribution) to 1.78 [1.64] ppb h−1, via the reaction of OH + CO, and 0.10 [0.09] to 0.63 [0.59] ppb h−1, via the reaction of CH3O2 (methylperoxy radical) + NO in the daytime average HO2 production rate, and 0.67 [0.61] to 4.32 [4.27] ppb h−1, via the reaction of HO2 + NO in the daytime average HO2 loss rate in Beijing, Shanghai, and Guangzhou. The above results suggest that the Punknown considerably enhanced the ROx concentrations and accelerated ROx cycles in the coastal regions of China, and could produce significant increases in concentrations of inorganic aerosols and secondary organic aerosols and further aggravate haze events in these regions.
Investigation of post-depositional processing of nitrate in East Antarctic snow: isotopic constraints on photolytic loss, re-oxidation, and source inputs

Snowpits along a traverse from coastal East Antarctica to the summit of the ice sheet (Dome Argus) are used to investigate the post-depositional processing of nitrate (NO3−) in snow. Seven snowpits from sites with accumulation rates between 24 and 172 kg m−2 a−1 were sampled to depths of 150 to 300 cm. At sites from the continental interior (low accumulation, < 55 kg m−2 a−1), nitrate mass fraction is generally > 200 ng g−1 in surface snow and decreases quickly with depth to < 50 ng g−1. Considerably increasing values of δ15N of nitrate are also observed (16–461 &permil; vs. air N2), particularly in the top 20 cm, which is consistent with predicted fractionation constants for the photolysis of nitrate. The δ18O of nitrate (17–84 &permil; vs. VSMOW (Vienna Standard Mean Ocean Water)), on the other hand, decreases with increasing δ15N, suggestive of secondary formation of nitrate in situ (following photolysis) with a low δ18O source. Previous studies have suggested that δ15N and δ18O of nitrate at deeper snow depths should be predictable based upon an exponential change derived near the surface. At deeper depths sampled in this study, however, the relationship between nitrate mass fraction and δ18O changes, with increasing δ18O of nitrate observed between 100 and 200 cm. Predicting the impact of post-depositional loss, and therefore changes in the isotopes with depth, is highly sensitive to the depth interval over which an exponential change is assumed. In the snowpits collected closer to the coast (accumulation > 91 kg m−2 a−1), there are no obvious trends detected with depth and instead seasonality in nitrate mass fraction and isotopic composition is found. In comparison to the interior sites, the coastal pits are lower in δ15N (−15–71 &permil; vs. air N2) and higher in δ18O of nitrate (53–111 &permil; vs. VSMOW). The relationships found amongst mass fraction, δ15N, δ18O and Δ17O (Δ17O = δ17O–0.52 × δ18O) of nitrate cannot be explained by local post-depositional processes alone, and are instead interpreted in the context of a primary atmospheric signal. Consistent with other Antarctic observational and modeling studies, the isotopic results are suggestive of an important influence of stratospheric ozone chemistry on nitrate formation during the cold season and a mix of tropospheric sources and chemistry during the warm season. Overall, the findings in this study speak to the sensitivity of nitrate isotopic composition to post-depositional processing and highlight the strength of combined use of the nitrogen and oxygen isotopes for a mechanistic understanding of this processing.
Current model capabilities for simulating black carbon and sulfate concentrations in the Arctic atmosphere: a multi-model evaluation using a comprehensive measurement data set

The concentrations of sulfate, black carbon (BC) and other aerosols in the Arctic are characterized by high values in late winter and spring (so-called Arctic Haze) and low values in summer. Models have long been struggling to capture this seasonality and especially the high concentrations associated with Arctic Haze. In this study, we evaluate sulfate and BC concentrations from eleven different models driven with the same emission inventory against a comprehensive pan-Arctic measurement data set over a time period of 2 years (2008–2009). The set of models consisted of one Lagrangian particle dispersion model, four chemistry transport models (CTMs), one atmospheric chemistry-weather forecast model and five chemistry climate models (CCMs), of which two were nudged to meteorological analyses and three were running freely. The measurement data set consisted of surface measurements of equivalent BC (eBC) from five stations (Alert, Barrow, Pallas, Tiksi and Zeppelin), elemental carbon (EC) from Station Nord and Alert and aircraft measurements of refractory BC (rBC) from six different campaigns. We find that the models generally captured the measured eBC or rBC and sulfate concentrations quite well, compared to previous comparisons. However, the aerosol seasonality at the surface is still too weak in most models. Concentrations of eBC and sulfate averaged over three surface sites are underestimated in winter/spring in all but one model (model means for January–March underestimated by 59 and 37 % for BC and sulfate, respectively), whereas concentrations in summer are overestimated in the model mean (by 88 and 44 % for July–September), but with overestimates as well as underestimates present in individual models. The most pronounced eBC underestimates, not included in the above multi-site average, are found for the station Tiksi in Siberia where the measured annual mean eBC concentration is 3 times higher than the average annual mean for all other stations. This suggests an underestimate of BC sources in Russia in the emission inventory used. Based on the campaign data, biomass burning was identified as another cause of the modeling problems. For sulfate, very large differences were found in the model ensemble, with an apparent anti-correlation between modeled surface concentrations and total atmospheric columns. There is a strong correlation between observed sulfate and eBC concentrations with consistent sulfate/eBC slopes found for all Arctic stations, indicating that the sources contributing to sulfate and BC are similar throughout the Arctic and that the aerosols are internally mixed and undergo similar removal. However, only three models reproduced this finding, whereas sulfate and BC are weakly correlated in the other models. Overall, no class of models (e.g., CTMs, CCMs) performed better than the others and differences are independent of model resolution.