skip to main content


Title: Development of an ecophysiology module in the GEOS-Chem chemical transport model version 12.2.0 to represent biosphere–atmosphere fluxes relevant for ozone air quality
Abstract. Ground-level ozone (O3) is a major air pollutant that adversely affects human health and ecosystem productivity. Removal of troposphericO3 by plant stomatal uptake can in turn cause damage to plant tissues with ramifications for ecosystem and crop health. In manyatmospheric and land surface models, the functionality of stomata opening is represented by a bulk stomatal conductance, which is oftensemi-empirically parameterized and highly fitted to historical observations. A lack of mechanistic linkage to ecophysiological processes such asphotosynthesis may render models inadequate to represent plant-mediated responses of atmospheric chemistry to long-term changes in CO2,climate, and short-lived air pollutant concentrations. A new ecophysiology module was thus developed to mechanistically simulate land−atmosphereexchange of important gas species in GEOS-Chem, a chemical transport model widely used in atmospheric chemistry studies. The implementation not onlyallows for dry deposition to be coupled with plant ecophysiology but also enables plant and crop productivity and functions to respond dynamically toatmospheric chemical changes. We conduct simulations to evaluate the effects of the ecophysiology module on simulated dry deposition velocity andconcentration of surface O3 against an observation-derived dataset known as SynFlux. Our estimated stomatal conductance and dry depositionvelocity of O3 are close to SynFlux with root-mean-squared errors (RMSEs) below 0.3 cm s−1 across different plant functionaltypes (PFTs), despite an overall positive bias in surface O3 concentration (by up to 16 ppbv). Representing ecophysiology wasfound to reduce the simulated biases in deposition fluxes from the prior model but worsen the positive biases in simulated O3concentrations. The increase in positive concentration biases is mostly attributable to the ecophysiology-based stomatal conductance being generallysmaller (and closer to SynFlux values) than that estimated by the prior semi-empirical formulation, calling for further improvements in non-stomataldepositional and non-depositional processes relevant for O3 simulations. The estimated global O3 deposition flux is864 Tg O3 yr−1 with GEOS-Chem, and the new module decreases this estimate by 92 Tg O3 yr−1. Estimated global grossprimary production (GPP) without O3 damage is 119 Pg C yr−1. O3-induced reduction in GPP is 4.2 Pg C yr−1(3.5 %). An elevated CO2 scenario (580 ppm) yields higher global GPP (+16.8 %) and lower global O3depositional sink (−3.3 %). Global isoprene emission simulated with a photosynthesis-based scheme is 317.9 Tg C yr−1, which is31.2 Tg C yr−1 (−8.9 %) less than that calculated using the MEGAN(Model of Emissions of Gases and Aerosols from Nature) emission algorithm. This new model development dynamicallyrepresents the two-way interactions between vegetation and air pollutants and thus provides a unique capability in evaluating vegetation-mediatedprocesses and feedbacks that can shape atmospheric chemistry and air quality, as well as pollutant impacts on vegetation health, especially for anytimescales shorter than the multidecadal timescale.  more » « less
Award ID(s):
1848372
NSF-PAR ID:
10432275
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Geoscientific Model Development
Volume:
16
Issue:
9
ISSN:
1991-9603
Page Range / eLocation ID:
2323 to 2342
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract. Our work explores the impact of two important dimensions of landsystem changes, land use and land cover change (LULCC) as well as directagricultural reactive nitrogen (Nr) emissions from soils, on ozone(O3) and fine particulate matter (PM2.5) in terms of air quality overcontemporary (1992 to 2014) timescales. We account for LULCC andagricultural Nr emissions changes with consistent remote sensingproducts and new global emission inventories respectively estimating theirimpacts on global surface O3 and PM2.5 concentrations as well as Nrdeposition using the GEOS-Chem global chemical transport model. Over thistime period, our model results show that agricultural Nr emissionchanges cause a reduction of annual mean PM2.5 levels over Europe andnorthern Asia (up to −2.1 µg m−3) while increasing PM2.5 levels in India, China and the eastern US (up to +3.5 µg m−3). Land cover changes induce small reductions in PM2.5 (up to −0.7 µg m−3) over Amazonia, China and India due to reduced biogenic volatile organic compound (BVOC) emissions and enhanced deposition of aerosol precursor gases (e.g., NO2, SO2). Agricultural Nr emissionchanges only lead to minor changes (up to ±0.6 ppbv) in annual meansurface O3 levels, mainly over China, India and Myanmar. Meanwhile, ourmodel result suggests a stronger impact of LULCC on surface O3 over the time period across South America; the combination of changes in drydeposition and isoprene emissions results in −0.8 to +1.2 ppbv surfaceozone changes. The enhancement of dry deposition reduces the surface ozone level (up to −1 ppbv) over southern China, the eastern US and central Africa. The enhancement of soil NO emission due to crop expansion also contributes to surface ozone changes (up to +0.6 ppbv) over sub-Saharan Africa. Incertain regions, the combined effects of LULCC and agricultural Nr emission changes on O3 and PM2.5 air quality can be comparable (>20 %) to anthropogenic emission changes over the same time period. Finally, we calculate that the increase in global agricultural Nr emissions leads to a net increase in global land area (+3.67×106km2) that potentially faces exceedance of the critical Nr load (>5 kg N ha−1 yr−1). Our result demonstrates the impacts of contemporary LULCC and agricultural Nr emission changes on PM2.5 and O3 in terms of air quality, as well as the importanceof land system changes for air quality over multidecadal timescales. 
    more » « less
  4. Abstract. We use the GEOS-Chem chemical transport model to examine theinfluence of bromine release from blowing-snow sea salt aerosol (SSA) onspringtime bromine activation and O3 depletion events (ODEs) in theArctic lower troposphere. We evaluate our simulation against observations oftropospheric BrO vertical column densities (VCDtropo) from the GOME-2 (second Global Ozone Monitoring Experiment)and Ozone Monitoring Instrument (OMI) spaceborne instruments for 3 years (2007–2009), as well asagainst surface observations of O3. We conduct a simulation withblowing-snow SSA emissions from first-year sea ice (FYI; with a surface snowsalinity of 0.1 psu) and multi-year sea ice (MYI; with a surface snowsalinity of 0.05 psu), assuming a factor of 5 bromide enrichment of surfacesnow relative to seawater. This simulation captures the magnitude ofobserved March–April GOME-2 and OMI VCDtropo to within 17 %, as wellas their spatiotemporal variability (r=0.76–0.85). Many of the large-scalebromine explosions are successfully reproduced, with the exception of eventsin May, which are absent or systematically underpredicted in the model. Ifwe assume a lower salinity on MYI (0.01 psu), some of the bromine explosionsevents observed over MYI are not captured, suggesting that blowing snow overMYI is an important source of bromine activation. We find that the modeledatmospheric deposition onto snow-covered sea ice becomes highly enriched inbromide, increasing from enrichment factors of ∼5 inSeptember–February to 10–60 in May, consistent with composition observations of freshly fallen snow. We propose that this progressive enrichment indeposition could enable blowing-snow-induced halogen activation to propagateinto May and might explain our late-spring underestimate in VCDtropo.We estimate that the atmospheric deposition of SSA could increase snow salinityby up to 0.04 psu between February and April, which could be an importantsource of salinity for surface snow on MYI as well as FYI covered by deepsnowpack. Inclusion of halogen release from blowing-snow SSA in oursimulations decreases monthly mean Arctic surface O3 by 4–8 ppbv(15 %–30 %) in March and 8–14 ppbv (30 %–40 %) in April. We reproduce atransport event of depleted O3 Arctic air down to 40∘ Nobserved at many sub-Arctic surface sites in early April 2007. While oursimulation captures 25 %–40 % of the ODEs observed at coastal Arctic surfacesites, it underestimates the magnitude of many of these events and entirelymisses 60 %–75 % of ODEs. This difficulty in reproducing observed surfaceODEs could be related to the coarse horizontal resolution of the model, theknown biases in simulating Arctic boundary layer exchange processes, thelack of detailed chlorine chemistry, and/or the fact that we did not includedirect halogen activation by snowpack chemistry. 
    more » « less
  5. Atmospheric chemistry models generally assume organic aerosol (OA) to be photochemically inert. Recent mechanisms for the oxidation of biogenic isoprene, a major source of secondary organic aerosol (iSOA), produce excessive OA in the absence of subsequent OA reactivity. At the same time, models underestimate atmospheric concentrations of formic and acetic acids for which OA degradation could provide a source. Here we show that the aqueous photooxidation of an isoprene-derived organosulfate (2-methyltriolsulfate or MTS), an important iSOA component, produces formic and acetic acids in high yields and at timescales competitive with deposition. Experimental data are well fit by a kinetic model in which three sequential oxidation reactions of the isoprene organosulfate produce two molar equivalents of formic acid and one of acetic acid. We incorporate this chemistry and that of 2-methyltetrol, another ubiquitous iSOA component, into the GEOS-Chem global atmospheric chemistry model. Simulations show that photooxidation and subsequent revolatilization of this iSOA may account for up to half of total iSOA loss globally, producing 4 Tg a−1 each of formic and acetic acids. This reduces model biases in gas-phase formic acid and total organic aerosol over the Southeast United States in summer by ∼30% and 60% respectively. While our study shows the importance of adding iSOA photochemical sinks into atmospheric models, uncertainties remain that warrant further study. In particular, improved understanding of reaction dependencies on particle characteristics and concentrations of particle-phase OH and other oxidants are needed to better simulate the effects of this chemistry on the atmospheric budgets of organic acids and iSOA. 
    more » « less