To simulate global mercury (Hg) dynamics in chemical transport models (CTMs), surface-atmosphere exchange of gaseous elemental mercury, Hg 0 , is often parameterized based on resistance-based dry deposition schemes coupled with a re-emission function, mainly from soils. Despite extensive use of this approach, direct evaluations of this implementation against field observations of net Hg 0 exchange are lacking. In this study, we evaluate an existing net exchange parameterization (referred to here as the base model) by comparing modeled fluxes of Hg 0 to fluxes measured in the field using micrometeorological techniques. Comparisons were performed in two terrestrial ecosystems: a grassland site in Switzerland and an Arctic tundra site in Alaska, U.S., each including summer and winter seasons. The base model included the dry deposition and soil re-emission parameterizations from Zhang et al. (2003) and the global CTM GEOS-Chem, respectively. Comparisons of modeled and measured Hg 0 fluxes showed large discrepancies, particularly in the summer months when the base model overestimated daytime net deposition by approximately 9 and 2 ng m −2 h −1 at the grassland and tundra sites, respectively. In addition, the base model was unable to capture a measured nighttime net Hg 0 deposition and wintertime deposition. We conducted a series of sensitivity analyses and recommend that Hg simulations using CTMs: (i) reduce stomatal uptake of Hg 0 over grassland and tundra in models by a factor 5–7; (ii) increase nighttime net Hg 0 deposition, e.g. , by increasing ground and cuticular uptake by reducing the respective resistance terms by factors of 3–4 and 2–4, respectively; and (iii) implement a new soil re-emission parameterization to produce larger daytime emissions and lower nighttime emissions. We also compared leaf Hg 0 uptake over the growing season estimated by the dry deposition model against foliar Hg measurements, which revealed good agreement with the measured leaf Hg concentrations after adjusting the base model as suggested above. We conclude that the use of resistance-based models combined with the new soil re-emission flux parameterization is able to reproduce observed diel and seasonal patterns of Hg 0 exchange in these ecosystems. This approach can be used to improve model parameterizations for other ecosystems if flux measurements become available.
more »
« less
A Vertically Resolved Canopy Improves Chemical Transport Model Predictions of Ozone Deposition to North Temperate Forests
Abstract Dry deposition is the second largest tropospheric ozone (O3) sink and occurs through stomatal and nonstomatal pathways. Current O3uptake predictions are limited by the simplistic big‐leaf schemes commonly used in chemical transport models (CTMs) to parameterize deposition. Such schemes fail to reproduce observed O3fluxes over terrestrial ecosystems, highlighting the need for more realistic treatment of surface‐atmosphere exchange in CTMs. We address this need by linking a resolved canopy model (1D Multi‐Layer Canopy CHemistry and Exchange Model, MLC‐CHEM) to the GEOS‐Chem CTM and use this new framework to simulate O3fluxes over three north temperate forests. We compare results with in situ measurements from four field studies and with standalone, observationally constrained MLC‐CHEM runs to test current knowledge of O3deposition and its drivers. We show that GEOS‐Chem overpredicts observed O3fluxes across all four studies by up to 2×, whereas the resolved‐canopy models capture observed diel profiles of O3deposition and in‐canopy concentrations to within 10%. Relative humidity and solar irradiance are strong O3flux drivers over these forests, and uncertainties in those fields provide the largest remaining source of model deposition biases. Flux partitioning analysis shows that: (a) nonstomatal loss accounts for 60% of O3deposition on average; (b) in‐canopy chemistry makes only a small contribution to total O3fluxes; and (c) the CTM big‐leaf treatment overestimates O3‐driven stomatal loss and plant phytotoxicity in these temperate forests by up to 7×. Results motivate the application of fully online vertically explicit canopy schemes in CTMs for improved O3predictions.
more »
« less
- PAR ID:
- 10618739
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Atmospheres
- Volume:
- 129
- Issue:
- 24
- ISSN:
- 2169-897X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract. During the Program for Research on Oxidants: PHotochemistry, Emissions, and Transport (PROPHET) campaign from 21 July to 3 August 2016,field experiments on leaf-level trace gas exchange of nitric oxide (NO), nitrogen dioxide (NO2), and ozone (O3) were conducted for thefirst time on the native American tree species Pinus strobus (eastern white pine), Acer rubrum (redmaple), Populus grandidentata (bigtooth aspen), and Quercus rubra (red oak) in a temperate hardwood forest inMichigan, USA. We measured the leaf-level trace gas exchange rates andinvestigated the existence of an NO2 compensation point, hypothesizedbased on a comparison of a previously observed average diurnal cycle ofNOx (NO2+NO) concentrations with that simulated using amulti-layer canopy exchange model. Known amounts of trace gases wereintroduced into a tree branch enclosure and a paired blank referenceenclosure. The trace gas concentrations before and after the enclosures weremeasured, as well as the enclosed leaf area (single-sided) and gas flow rate to obtain the trace gas fluxes with respect to leaf surface. There was nodetectable NO uptake for all tree types. The foliar NO2 and O3uptake largely followed a diurnal cycle, correlating with that of the leafstomatal conductance. NO2 and O3 fluxes were driven by theirconcentration gradient from ambient to leaf internal space. The NO2 loss rate at the leaf surface, equivalently the foliar NO2 deposition velocity toward the leaf surface, ranged from 0 to 3.6 mm s−1 for bigtooth aspen and from 0 to 0.76 mm s−1 for red oak, both of which are∼90 % of the expected values based on the stomatalconductance of water. The deposition velocities for red maple and white pineranged from 0.3 to 1.6 and from 0.01 to 1.1 mm s−1, respectively, and were lower than predicted from the stomatal conductance, implying amesophyll resistance to the uptake. Additionally, for white pine, theextrapolated velocity at zero stomatal conductance was 0.4±0.08 mm s−1, indicating a non-stomatal uptake pathway. The NO2compensation point was ≤60 ppt for all four tree species andindistinguishable from zero at the 95 % confidence level. This agrees withrecent reports for several European and California tree species butcontradicts some earlier experimental results where the compensation pointswere found to be on the order of 1 ppb or higher. Given that the sampledtree types represent 80 %–90 % of the total leaf area at this site, theseresults negate the previously hypothesized important role of a leaf-scaleNO2 compensation point. Consequently, to reconcile these findings,further detailed comparisons between the observed and simulated in- and above-canopy NOx concentrations and the leaf- and canopy-scaleNOx fluxes, using the multi-layer canopy exchange model withconsideration of the leaf-scale NOx deposition velocities as well asstomatal conductances reported here, are recommended.more » « less
-
null (Ed.)Abstract. At the leaf level, stomata control the exchange of water and carbon across the air–leaf interface. Stomatal conductance is typically modeledempirically, based on environmental conditions at the leaf surface. Recently developed stomatal optimization models show great skills at predictingcarbon and water fluxes at both the leaf and tree levels. However, how well the optimization models perform atlarger scales has not been extensively evaluated. Furthermore, stomatal models are often used with simple single-leaf representations of canopy radiative transfer (RT), such asbig-leaf models. Nevertheless, the single-leaf canopy RT schemes do not have the capability to model optical properties of the leaves nor the entirecanopy. As a result, they are unable to directly link canopy optical properties with light distribution within the canopy to remote sensing dataobserved from afar. Here, we incorporated one optimization-based and two empirical stomatal models with a comprehensive RT model in the landcomponent of a new Earth system model within CliMA, the Climate Modelling Alliance. The model allowed us to simultaneously simulate carbon and waterfluxes as well as leaf and canopy reflectance and fluorescence spectra. We tested our model by comparing our modeled carbon and water fluxes andsolar-induced chlorophyll fluorescence (SIF) to two flux tower observations (a gymnosperm forest and an angiosperm forest) and satellite SIFretrievals, respectively. All three stomatal models quantitatively predicted the carbon and water fluxes for both forests. The optimization model,in particular, showed increased skill in predicting the water flux given the lower error (ca. 14.2 % and 21.8 % improvement for thegymnosperm and angiosperm forests, respectively) and better 1:1 comparison (slope increases from ca. 0.34 to 0.91 for the gymnosperm forest andfrom ca. 0.38 to 0.62 for the angiosperm forest). Our model also predicted the SIF yield, quantitatively reproducing seasonal cycles for bothforests. We found that using stomatal optimization with a comprehensive RT model showed high accuracy in simulating land surface processes. Theever-increasing number of regional and global datasets of terrestrial plants, such as leaf area index and chlorophyll contents, will helpparameterize the land model and improve future Earth system modeling in general.more » « less
-
Abstract. We implement the GEOS-Chem chemistry module as a chemical mechanism in version 2 of the Community Earth System Model (CESM). Our implementation allowsthe state-of-the-science GEOS-Chem chemistry module to be used with identical emissions, meteorology, and climate feedbacks as the CAM-chemchemistry module within CESM. We use coupling interfaces to allow GEOS-Chem to operate almost unchanged within CESM. Aerosols are converted at eachtime step between the GEOS-Chem bulk representation and the size-resolved representation of CESM's Modal Aerosol Model (MAM4). Land-type informationneeded for dry-deposition calculations in GEOS-Chem is communicated through a coupler, allowing online land–atmosphere interactions. Wet scavengingin GEOS-Chem is replaced with the Neu and Prather scheme, and a common emissions approach is developed for both CAM-chem and GEOS-Chem in CESM. We compare how GEOS-Chem embedded in CESM (C-GC) compares to the existing CAM-chem chemistry option (C-CC) when used to simulate atmosphericchemistry in 2016, with identical meteorology and emissions. We compare the atmospheric composition and deposition tendencies between the twosimulations and evaluate the residual differences between C-GC and its use as a stand-alone chemistry transport model in the GEOS-Chem HighPerformance configuration (S-GC). We find that stratospheric ozone agrees well between the three models, with differences of less than 10 % inthe core of the ozone layer, but that ozone in the troposphere is generally lower in C-GC than in either C-CC or S-GC. This is likely due to greatertropospheric concentrations of bromine, although other factors such as water vapor may contribute to lesser or greater extents depending on theregion. This difference in tropospheric ozone is not uniform, with tropospheric ozone in C-GC being 30 % lower in the Southern Hemisphere whencompared with S-GC but within 10 % in the Northern Hemisphere. This suggests differences in the effects of anthropogenic emissions. Aerosolconcentrations in C-GC agree with those in S-GC at low altitudes in the tropics but are over 100 % greater in the upper troposphere due todifferences in the representation of convective scavenging. We also find that water vapor concentrations vary substantially between the stand-aloneand CESM-implemented version of GEOS-Chem, as the simulated hydrological cycle in CESM diverges from that represented in the source NASA Modern-Era Retrospective analysis for Research and Applications (Version 2; MERRA-2)reanalysis meteorology which is used directly in the GEOS-Chem chemistrytransport model (CTM). Our implementation of GEOS-Chem as a chemistry option in CESM (including full chemistry–climate feedback) is publicly available and is beingconsidered for inclusion in the CESM main code repository. This work is a significant step in the MUlti-Scale Infrastructure for Chemistry andAerosols (MUSICA) project, enabling two communities of atmospheric researchers (CESM and GEOS-Chem) to share expertise through a common modelingframework, thereby accelerating progress in atmospheric science.more » « less
-
Abstract. The impact of biomass burning (BB) on the atmospheric burden of volatile organic compounds (VOCs) is highly uncertain. Here we apply the GEOS-Chemchemical transport model (CTM) to constrain BB emissions in the western USA at ∼ 25 km resolution. Across three BB emission inventorieswidely used in CTMs, the inventory–inventory comparison suggests that the totals of 14 modeled BB VOC emissions in the western USA agree with eachother within 30 %–40 %. However, emissions for individual VOCs can differ by a factor of 1–5, driven by the regionally averaged emissionratios (ERs, reflecting both assigned ERs for specific biome and vegetation classifications) across the three inventories. We further evaluate GEOS-Chemsimulations with aircraft observations made during WE-CAN (Western Wildfire Experiment for Cloud Chemistry, Aerosol Absorption and Nitrogen) andFIREX-AQ (Fire Influence on Regional to Global Environments and Air Quality) field campaigns. Despite being driven by different global BBinventories or applying various injection height assumptions, the model–observation comparison suggests that GEOS-Chem simulations underpredictobserved vertical profiles by a factor of 3–7. The model shows small to no bias for most species in low-/no-smoke conditions. We thus attribute thenegative model biases mostly to underestimated BB emissions in these inventories. Tripling BB emissions in the model reproduces observed verticalprofiles for primary compounds, i.e., CO, propane, benzene, and toluene. However, it shows no to less significant improvements for oxygenatedVOCs, particularly for formaldehyde, formic acid, acetic acid, and lumped ≥ C3 aldehydes, suggesting the model is missing secondarysources of these compounds in BB-impacted environments. The underestimation of primary BB emissions in inventories is likely attributable tounderpredicted amounts of effective dry matter burned, rather than errors in fire detection, injection height, or ERs, as constrained by aircraftand ground measurements. We cannot rule out potential sub-grid uncertainties (i.e., not being able to fully resolve fire plumes) in the nestedGEOS-Chem which could explain the negative model bias partially, though back-of-the-envelope calculation and evaluation using longer-term groundmeasurements help support the argument of the dry matter burned underestimation. The total ERs of the 14 BB VOCs implemented in GEOS-Chem onlyaccount for half of the total 161 measured VOCs (∼ 75 versus 150 ppb ppm−1). This reveals a significant amount of missing reactiveorganic carbon in widely used BB emission inventories. Considering both uncertainties in effective dry matter burned (× 3) and unmodeledVOCs (× 2), we infer that BB contributed to 10 % in 2019 and 45 % in 2018 (240 and 2040 Gg C) of the total VOC primaryemission flux in the western USA during these two fire seasons, compared to only 1 %–10 % in the standard GEOS-Chem.more » « less
An official website of the United States government

