Examining the competing effects of contemporary land management vs. land cover changes on global air quality

Our work explores the impact of two important dimensions of land system changes, land use and land cover change (LULCC) as well as direct agricultural reactive nitrogen (Nr) emissions from soils, on ozone (O3) and fine particulate matter (PM2.5) in terms of air quality over contemporary (1992 to 2014) timescales. We account for LULCC and agricultural Nr emissions changes with consistent remote sensing products and new global emission inventories respectively estimating their impacts on global surface O3 and PM2.5 concentrations as well as Nr deposition using the GEOS-Chem global chemical transport model. Over this time period, our model results show that agricultural Nr emission changes cause a reduction of annual mean PM2.5 levels over Europe and northern 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 emission changes only lead to minor changes (up to ±0.6 ppbv) in annual mean surface O3 levels, mainly over China, India and Myanmar. Meanwhile, our model result suggests a stronger impact of LULCC on surface O3 over the time period across South America; the combination of changes in dry deposition and isoprene emissions results in −0.8 to +1.2 ppbv surface ozone 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. In certain 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×106 km2) 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 importance of land system changes for air quality over multidecadal timescales.

Abstract. Our work explores the impact of two important dimensions of land system changes, land use and land cover change (LULCC) as well as direct agricultural reactive nitrogen (N r ) emissions from soils, on ozone (O 3 ) and fine particulate matter (PM 2.5 ) in terms of air quality over contemporary (1992 to 2014) timescales. We account for LULCC and agricultural N r emissions changes with consistent remote sensing products and new global emission inventories respectively estimating their impacts on global surface O 3 and PM 2.5 concentrations as well as N r deposition using the GEOS-Chem global chemical transport model. Over this time period, our model results show that agricultural N r emission changes cause a reduction of annual mean PM 2.5 levels over Europe and northern Asia (up to −2.1 µg m −3 ) while increasing PM 2.5 levels in India, China and the eastern US (up to +3.5 µg m −3 ). Land cover changes induce small reductions in PM 2.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., NO 2 , SO 2 ). Agricultural N r emission changes only lead to minor changes (up to ±0.6 ppbv) in annual mean surface O 3 levels, mainly over China, India and Myanmar. Meanwhile, our model result suggests a stronger impact of LULCC on surface O 3 over the time period across South America; the combination of changes in dry deposition and isoprene emissions results in −0.8 to +1.2 ppbv surface ozone 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. In certain regions, the combined effects of LULCC and agricultural N r emis-sion changes on O 3 and PM 2.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 N r emissions leads to a net increase in global land area (+3.67 × 10 6 km 2 ) that potentially faces exceedance of the critical N r load (>5 kg N ha −1 yr −1 ). Our result demonstrates the impacts of contemporary LULCC and agricultural N r emission changes on PM 2.5 and O 3 in terms of air quality, as well as the importance of land system changes for air quality over multidecadal timescales.

Introduction
The broad term "land use and land cover change (LULCC)" encapsulates both the anthropogenic (e.g., agricultural expansion) and natural (e.g., ecological succession) dimensions of terrestrial biome changes (Reick et al., 2013), which alter the physical and ecophysiological properties of the land surface. These perturbations alter the transfer and uptake of air pollutants by ecosystems and can also have large impacts on the emission of biogenic volatile organic compounds (BVOCs), which play vital roles in tropospheric ozone (O 3 ) and secondary organic aerosol (SOA) formation (Fu and Tai, 2015;Ganzeveld et al., 2010;Heald and Geddes, 2016;Heald and Spracklen, 2015;Squire et al., 2014;Wu et al., 2012a).
Agricultural activities, in addition to being a large driver of LULCC (e.g., Ellis, 2015;Ellis et al., 2013;Goldewijk et al., 2017;Kaplan et al., 2011), also introduce an enormous amount of reactive nitrogen into the soil (Galloway et al., 2008), which can be emitted into the atmosphere either as oxidized or reduced nitrogen. The reactive nitrogen Published by Copernicus Publications on behalf of the European Geosciences Union. 16480 A. Y. H. Wong and J. A. Geddes: Effects of contemporary land management vs. land cover changes oxides emitted from soil, NO x (≡ NO + NO 2 ), enhance O 3 production when volatile organic compounds (VOCs) are relatively abundant (i.e., NO x -limited regimes) but suppresses O 3 production when the concentration of VOCs is relatively low (i.e., VOC-limited regimes) (Sillman et al., 1990). Reactive nitrogen also contributes to aerosol formation. Ammonia (NH 3 ) can combine with nitrate and sulfate ions to form secondary inorganic aerosol, while the emissions of NO x can oxidize further and contribute to particulate nitrate formation (Ansari and Pandis, 1998). Indeed, agricultural emissions are the dominant global anthropogenic source of NH 3 (Hoesly et al., 2018) and have been identified as a major contributor to global premature mortality due to particulate matter (PM) pollution (Lelieveld et al., 2015). Trends in atmospheric reactive nitrogen also affect nitrogen deposition (e.g., Geddes and Martin, 2017), with potentially negative impacts on biodiversity (e.g., Bobbink et al., 2010;Payne et al., 2017;Wal-lisDeVries and Bobbink, 2017) and eutrophication of aquatic ecosystems (e.g., Fenn et al., 2003). These ecosystem impacts may contribute to economic loss comparable to the benefits of extra crop output from LULCC and agricultural emissions (Paulot and Jacob, 2014;Sobota et al., 2015).
Even while land cover at a particular location may not change, modifications in human management of the land (e.g., intensification of agriculture, irrigation practices, fertilizer application, selective harvesting) may still be associated with changes in pollutant emission and uptake. An obvious example would be a region where direct agricultural emissions may have changed without any concomitant changes in land cover. Reducing NH 3 emission, particularly from the agricultural sector, has been explored as a potent strategy of controlling PM pollution (Giannadaki et al., 2018;Pinder et al., 2007;Pozzer et al., 2017). Bauer et al. (2016) suggest that agricultural emissions are the main source of presentday PM 2.5 (fine particulate with an aerodynamic diameter less than 2.5 µm) over the eastern US, Europe and northern China. However, as anthropogenic NO x and SO 2 emissions are expected to be lower in the future, some aerosol formation chemistry is expected to become less sensitive to NH 3 emissions.
The potential impacts of LULCC and agricultural emission changes on air quality have been explored previously. To date, this work has focused on future projections in land use (Bauer et al., 2016;Ganzeveld et al., 2010;Hardacre et al., 2013;Heald et al., 2008;Squire et al., 2014;Tai et al., 2013;Wu et al., 2012b), contrasted preindustrial estimates of land cover and agricultural emissions with present-day conditions Hollaway et al., 2017), or has been regional in focus (Fu and Tai, 2015;Geddes et al., 2015;Silva et al., 2016). For example, Wu et al. (2012) propose that LULCC induced by climate, CO 2 abundance and agriculture could significantly affect surface O 3 in the future, mainly through modulating dry deposition and isoprene emissions. Over more contemporary timescales (e.g., across the last 20-30 years), Fu et al. (2016), Fu and Tai (2015), and Silva et al. (2016) find that LULCC could have impacts on O 3 and PM in terms of air quality over China and Southeast Asia.
Given the large spatial scale of LULCC (e.g., Hansen et al., 2013;Li et al., 2018) and agricultural emission changes (e.g., Crippa et al., 2018;Hoesly et al., 2018;Xu et al., 2019) over recent decades, these two land system changes could be substantially contributing to global trends in O 3 and PM pollution. While changes in land cover and agricultural emissions actually occur simultaneously across the globe, they are rarely considered together in simulations of air quality from chemical transport models. The importance of studying these combined processes at the same time was highlighted by Ganzeveld et al. (2010) in their analysis of air quality impacts from future land use and land cover changes. In this study, for example, opposing effects on O 3 were simulated, with decreases in tropical forest soil NO x emissions being compensated for by increases in soil NO x emissions associated with agriculture. Still, this work did not explore the concomitant changes in ammonia emissions that would be expected with the changes in agricultural activity. It remains unclear to what extent LULCC may either amplify or offset the impacts of some of the associated agricultural emission changes, how this may vary regionally, and to what extent these land system impacts may compare to concomitant changes resulting from other direct anthropogenic emissions (e.g., emissions from industrial and transport sectors).
Consistent long-term land records of land cover derived from satellite remote sensing observations and global anthropogenic emission inventories have become readily available. This opens an opportunity for a more holistic and observationally constrained assessment of the impacts on global O 3 and PM in terms of air quality from contemporary changes in LULCC and agricultural emissions simultaneously, which has been advocated by Ganzeveld et al. (2010), and a comparison of these to the effects of direct anthropogenic emissions. In this study, we model the effects of contemporary LULCC and agriculture emissions changes on global surface O 3 and PM 2.5 levels and gauge their importance relative to changes in other direct anthropogenic emissions over the same period of time. We also highlight the effect of agricultural emissions changes on nitrogen deposition on land ecosystems. Through our chemical transport model predictions, we aim to identify potential global hotspots of contemporary land changes that may be substantially altering trends in air quality and nitrogen deposition.

Methods
To simulate global changes in surface O 3 and PM 2.5 concentrations due to LULCC, agricultural emissions, and direct anthropogenic emissions over 1992 to 2014, we use the GEOS-Chem chemical transport model (version 12.7.0, available at https://doi.org/10.5281/zenodo.3634864). We choose  199219921992Land cover 19921992Anthropogenic emissions 19921992 our timeframe due to the availability of consistent highresolution remote sensing products (PFT and LAI maps) and concurrent global emission inventories. We define "direct anthropogenic" and "agricultural" emissions separately in more detail below.
We perform five sets of simulations summarized in Table 1: (1) a "baseline" scenario in which land cover, agricultural emissions and direct anthropogenic emissions are all set to 1992 levels; (2) an "anthropogenic emission" scenario in which direct anthropogenic emissions are updated to 2014 levels; (3) an "anthropogenic emissions and land cover change" scenario in which anthropogenic emissions remain updated to 2014, with land cover inputs now prescribed based on updated 2014 data ; and (4) an "anthropogenic emissions, land cover and agricultural emission change" scenario in which direct anthropogenic emissions and land cover inputs remain updated to 2014, with agricultural emissions also updated to 2014 levels. To test the chemical sensitivity of our results, (5) is performed with anthropogenic emissions held at 1992 levels, but land cover change and agricultural emissions are updated to 2014 levels.
The role of direct anthropogenic emission changes can be evaluated by comparing simulations (1) and (2); the additional role played by land cover changes over this time period is evaluated by comparing simulations (2) and (3); and finally the additional impact of agricultural emission changes is evaluated by comparing simulations (3) and (4). The latter two effects will be the focus of this paper, but we compare these to the role of direct anthropogenic emission changes for context. Since changes in surface ozone and PM 2.5 should be sensitive to the NO x -to-VOC ratio and availability of NO − 3 and SO 2− 4 ions, the sensitivity of the effects from land cover change and agricultural emission changes to anthropogenic emission changes can be quantified by evaluating simulation (5).
We use assimilated meteorological fields from Modern-Era Retrospective analysis for Research and Applications Version 2 (MERRA-2) (Gelaro et al., 2017) to drive GEOS-Chem. All simulations are carried out at 2 • latitude by 2.5 • longitude resolution over the globe, using identical meteorological fields from 2011 to 2014 in order to exclude meteorological variability from the analysis. The output from 2011 is discarded as spin-up. The GEOS-Chem model simulates O 3 chemistry with a comprehensive HO x −NO x −VOC−O 3 −BrO x chemical mechanism Mao et al., 2013). Gaseous dry deposition follows Wang et al. (1998) and Wesely (1989), while particle deposition follows Zhang et al. (2001). In GEOS-Chem, the surface exchange modules are unidirectional (which implies that the effects of bidirectional exchanges of trace gases are not explicitly modeled). In certain regions for which the Community Emission Data System (CEDS) inventory scales the calculated emissions to a regional inventory, the extent of accounting for bidirectional exchange may depend on the underlying assumptions in the regional inventory modeling. For example, agricultural ammonia emissions from NEI for the United States include considering bidirectional ammonia exchange modeling from the Community Multiscale Air Quality Modeling System (CMAQ) (U.S. EPA, 2018). However, we cannot comment with certainty on how this is treated elsewhere across the globe, so we assume that neglecting bidirectional exchange of ammonia (and other species for which an atmospheric compensation point may exist) introduces some uncertainty in our simulation (which we discuss in a subsequent section). Wet deposition is described by Liu et al. (2001) with updates from Amos et al. (2012) and Wang et al. (2011Wang et al. ( , 2014. The recent update from Luo et al. (2019) on wet deposition parameterization is also included to improve model-observation agreement for sulfatenitrate-ammonium (SNA) aerosol. The thermodynamics and gas-aerosol partitioning of the NH 3 −H 2 SO 4 −HNO 3 system are simulated by ISORROPIA II module (Fountoukis and Nenes, 2007). A simple yield-based secondary organic aerosol (SOA) estimate is also included (Kim et al., 2015). Other types of aerosol represented in the model include sea salt, dust, primary black carbon (BC) and organic carbon (OC). The total PM 2.5 mass is then calculated at 35 % relative humidity for consistency with the measurement standard in the US.
We use anthropogenic and agricultural emissions based on the Community Emission Data System (CEDS) inventory (Hoesly et al., 2018), which contains the estimates of anthropogenic NO x , non-methane volatile organic compounds (NMVOCs), CO, BC, OC, SO 2 and NH 3 emissions harmonized from a wide range of global and regional inventories. In this inventory, emissions are from six major sectors: energy production, industry, transportation, RCO (residential, commercial, other), agriculture and waste. For this study, "agricultural emissions" specifically refer to NO x and NH 3 emitted from fertilizer application and manure management, which correspond directly to agricultural nitrogen input. We do not consider the changes in agriculture for other trace species (e.g., CH 4 , SO 2 , CO). For simplicity, we assume that agricultural emissions from fertilizer application in CEDS represent "above-canopy" emissions to the atmosphere (instead of making assumptions about the implicit treatment of canopy reduction over each region). We note that the fertilizer emissions represent only a fraction of the total agricultural NH 3 emissions we are considering here (e.g., which also include livestock operation) so that uncertainty in canopy reduction will only affect a fraction of the total. Likewise, fertilizer NO x emissions are small compared to the total soil NO x emissions (for which canopy reduction is accounted for online in the Hudman et al., 2012, parameterization).
Biogenic volatile organic compound emissions are calculated by the Model of Emissions of Gases and Aerosols from Nature (MEGAN) v2.1 (Guenther et al., 2012). Soil NO x emission follows Hudman et al. (2012), with fertilizer emissions zeroed out to avoid double counting with the agricultural NO x emission in the CEDS inventory. Fire (Global Fire Emissions Database v4.1; Van Der Werf et al., 2017) and lightning (Murray et al., 2012) emissions are held constant at the 2014 level.
We use the European Space Agency Climate Change Initiative (ESA CCI) land cover map  to characterize LULCC and drive the biosphere-atmosphere emission fluxes in our simulations. The ESA CCI land cover map is a consistent global annual land cover time series derived from satellite observations from the AVHRR, MERIS, SOPTVGT and PROBA-V instruments. It has a native spatial resolution of 300 m following the United Nations Land Cover Classification System. Time-consistent land surface characterization also requires leaf area index (LAI) data. We use the Global Land Surface Satellite (GLASS) product  (retrieved from http://globalchange.bnu.edu.cn/, last access: 15 May 2020), which is a global LAI time series combining AVHRR and MODIS observations. The 3-year average (1991-1993 average LAI for 1992 land cover, 2013-2015 average LAI for 2014 land cover) is used as input for LAI to GEOS-Chem to reduce the possible effect of interannual variability.
This satellite-derived land surface characterization on its own is not directly compatible with the input to the vegetation-related modules in GEOS-Chem; it thus requires further harmonization (dry deposition, BVOC emissions, soil NO x emissions), which is a common problem for simulations involving land change (e.g., Geddes et al., 2016). We first aggregate and process the ESACCI land cover map with the tool and crosswalk table provided with the land cover product to derive the percentage coverage of plant functional type (PFT) at 0.05 • resolution, which is the native resolution of GLASS LAI. The dominant surface type can be readily mapped to the 11 deposition surface types in the Wesely (1989) dry deposition model. We adopt the approach of Geddes et al. (2016) to replace roughness length (z 0 ) from assimilated meteorology with that prescribed for each deposition surface type. We ignore changes in displacement height as they are expected to be much less important than the changes in z o (Text S1 in the Supplement). To derive the MODIS-Köppen type land map (Steinkamp and Lawrence, 2011) required for the soil NO x module, we first translate the PFT map according to the International Geosphere-Biosphere Programme (IGBP) land cover classification system (available at: http://www.eomf.ou.edu/static/IGBP.pdf, last access: 15 May 2020). We use global monthly temper-ature climatology (Matsuura and Willmott, 2012) to further differentiate the land types by climate with criteria outlined by Kottek et al. (2006). Finally, the ESA CCI PFT map is converted to a Community Land Model (CLM) PFT map, which is required for the MEGAN BVOC emissions module, by the temperature criteria specified by Bonan et al. (2002). As the method of deriving the C 3 -to-C 4 grass ratio was subsequently updated (Lawrence and Chase, 2007), this ratio is directly taken from the CLM land surface dataset.
In the Supplement, we provide an evaluation of the annual mean simulated SNA aerosol mass concentration and surface O 3 mixing ratios from simulation (4) (representative of 2014 conditions) with globally available observations from the same time period. In general, the model captures the spatial distributions of individual SNA species reasonably well (Fig. S1 in the Supplement). The model is able to capture regional annual means of individual SNA species (Table S1 in the Supplement) over Europe. Over the US and China, where annual means of all SNA species are underestimated by 21 %-55 %, and in regions covered by the Acid Deposition Monitoring Network in East Asia (Japan, Korea and southeast Asia) where SO 2− 4 is underestimated by 36 %, we expect the model may underestimate the sensitivity of SNA concentration to NH 3 emission perturbations. This may imply that results from our study should be interpreted as conservative. Figure S2 shows the reasonable agreement on annual mean surface O 3 between our model output and the gridded observation dataset from Sofen et al. (2016) (mean bias = +1.81 ppbv and mean absolute error = 3.97 ppbv). Our model therefore captures the present-day annual means of surface SNA and O 3 concentrations, providing a basis for our subsequent analyses. We also provide definitions for geographical regions, which largely follow Integrated Modeling of Global Environmental Change (IMAGE) 2.4 classifications, in Table S2.
3 Changes in land cover, biospheric fluxes and agricultural emissions Table 2 shows the changes in the global coverage of the major land cover types from 1992 to 2014 derived by the ESA CCI land cover product. The coverage of managed grass (including cropland and pasture) and built-up area, both of which are unmistakably related to human activities, has increased mainly at the expense of forest coverage. This is consistent with a global trend in deforestation over this period. Figure 1 shows the spatial distribution of changes in fractional coverage of the major land cover types. Expansion of agricultural land at the expense of broadleaf forest coverage is most notable in South America and Southeast Asia, which is well-documented in other studies based on remote sensing (Hansen et al., 2013) and national surveys (Keenan et al., 2015). The expansion of agricultural land over this time period is also observed in central Asia, southern China and Africa, but usually at the expense of land types other than broadleaf forests (mainly primary grassland and needleleaf forests). Meanwhile, transitions from agricultural land to forests and built-up areas are observed in northern China and eastern Europe, consistent with the findings of Potapov et al. (2015) and Lai et al. (2016). Figure 2 shows the global changes in the 3-year (2012-2014 minus 1991-1993) annual mean LAI calculated from the GLASS LAI dataset. Over southern China and South America, the area with regionally consistent deforestation experience general increases in LAI, while the opposite effect is observed in the Sahel and the former Soviet Union. In Europe, LAI increases in most parts despite a fairly consistent retraction of agricultural land being observed. The agricultural expansion and deforestation over Southeast Asia are mostly concurrent with the LAI decreases. LAI increases notably in northern China where agricultural land decreases. The fact that LAI change can be driven by factors other than changes in land cover type (e.g., temperature, precipitation, atmospheric CO 2 level) (e.g., Zhu et al., 2016) may explain the regionally divergent trend response of LAI to agricultural land use change. For example, the general increase in LAI in China is not only driven by changes in biome types, but also the greening within cropland (mainly attributable to agricultural intensification) and forests (mainly attributable to ambitious tree planting programs) . Similarly, some deforested land in South America might have been cultivated intensively, resulting in an increase rather decrease in LAI. We also note that since the relationship between satellite-derived surface reflectance and retrieved LAI depends on land cover, the use of a static land cover map in long-term LAI retrievals (Claverie et al., 2016;Zhu et al., 2013) may not fully capture the effect of LULCC on LAI (Fang et al., 2013). In particular, Fang et al. (2013) show that LAI could be substantially overestimated when grasses and crops are misclassified as forest. We may therefore overestimate dry deposition velocity over regions with significant deforestation. Such an impact on biogenic emissions is secondary as biogenic emissions are expected to be much more sensitive to land cover type than LAI (e.g., Guenther et al., 2012).
These changes in land cover produce changes in the biogenic fluxes of reactive trace gases between the Earth's surface and atmosphere derived by GEOS-Chem. Figure 3a shows the calculated changes in annual mean isoprene emission due to land cover change over 1992 to 2014 and suggests that global isoprene emission could have decreased by 5.12 Tg yr −1 (−1.5 %). The largest local reductions in isoprene emissions (up to 30 %) are observed in parts of South America, where deforestation from highly isoprene-emitting broadleaf forests is most strongly observed. We note that the decrease in isoprene emission simulated in Southeast Asia does not agree with the result from Silva et al. (2016), since our remote sensing data do not have a separate land cover class for oil palm plantations, which have expanded dramati-cally in the region. Our model therefore may not capture the full effects of LULCC on isoprene emission and its effect on PM 2.5 and O 3 over the region. Elsewhere in the world, the signals of land cover change in isoprene emissions are mostly small and follow the local patterns of changes in LAI. Changes in monoterpene (<5 ng m −2 s −1 ) and sesquiterpene (<1 ng m −2 s −1 ) emissions are relatively small. Figure 3b shows the changes in annual mean soil NO emission due to LULCC, which represent the change in soil emission driven purely by LAI (which affects canopy uptake) and land cover changes (which affects both the biome-based emission factor and canopy uptake) (i.e., without considering the changes in nitrogen input). LULCC leads to a small signal of +0.04 Tg yr −1 (+0.6 %) in global soil NO emissions. The magnitude of changes in soil NO emissions induced by LULCC is comparable to that in agricultural NO emission inventories (see below) over certain regions (e.g., South America, Australia, Africa). Relatively large increases in soil NO are simulated over western Africa due to both cropland expansion and LAI reduction, which leads to a smaller canopy reduction factor and larger emission factor. Figure 3c shows the changes in annual mean O 3 dry deposition velocity (v d ), which also closely follow the pattern of LAI changes. Slight increases in v d are observed in China, India, the southeastern US, Central America, South America, Europe and southern Africa. In Southeast Asia v d decreases concurrently with deforestation and reduction in LAI. In central Brazil, the increase in LAI is offset by the deforestation of tropical evergreen broadleaf forests that have higher v d than other land types (Song-Miao Fan et al., 1990;Wang et al., 1998), leading to small overall change in v d . Likewise, despite deforestation observed further south, these losses are offset by strong increases in LAI so that v d increases by up to 0.1 cm s −1 . Significant changes in the v d of O 3 due to LAI also imply that v d of other relevant trace gases (e.g., NO 2 , SO 2 ) would also be perturbed by land cover change in our model, which will be discussed briefly in the subsequent section. Figure 4 shows the changes in agricultural NH 3 and NO x emissions between 1992 and 2014, which consists mostly of emissions from fertilizer application and manure management (Hoesly et al., 2018). According to the CEDS inventory, global direct agricultural NH 3 emissions have increased by 7.6 Tg N yr −1 since 1992, equivalent to a 19 % increase in total anthropogenic NH 3 emissions. Direct agricultural soil NO x emissions have increased by 0.37 Tg N yr −1 since 1992, and while this is a substantial increase in agricultural soil NO x emissions (26 %), it represents only a 1 % increase in total anthropogenic NO x emissions.
The increases in agricultural emissions are most substantial over South Asia, followed by China, parts of the Middle East, Southeast Asia and South America, and to a lesser degree Central America, North America and the Sahel. The sharpest decline of agricultural emissions is observed in Europe and the former Soviet Union, followed by milder de-   (Skjøth and Hertel, 2013). According to the CEDS inventory, changes in agricultural emissions dominate the trend of total NH 3 emissions in all major regions except Africa, where a large part of the NH 3 emissions trend is attributable to the waste management and RCO (residential, commercial, other) sectors (Hoesly et al., 2018) (Fig. S3).
In contrast, the increase in agricultural emissions of NO x does not contribute significantly to the global increase in total NO x emissions over our period of concern. We note that the hotspots of change in managed land cover and of change in agricultural emissions do not always overlap. For example, agricultural emissions increase significantly over northern China and northern India, while the cropland coverage over those regions does not increase correspondingly over this same period. Such agricultural intensification in turn significantly contributes to the positive LAI trend over the above regions . Similarly, agricultural emissions have declined over Kazakhstan, while the area of managed land has not decreased significantly. This highlights a degree of independence between land management and LULCC, with both being components of land change but having potentially distinct spatial patterns and impacts on air quality. This also highlights the importance of treating both in our chemical transport model simulations as they occur contemporaneously around the globe and may have different impacts on air quality. 4 Impact of LULCC and agricultural emission changes on surface PM 2.5 Figure 5 shows the modeled impacts of LULCC, changes in agricultural emissions and the combined effects of both on annual mean surface PM 2.5 (under 2014 anthropogenic emissions). We have calculated the impacts of LULCC on PM 2.5 ( PM 2.5,LULCC ) as the difference in PM 2.5 predicted by simulation (3) and simulation (2), the impacts of agricultural emission changes on PM 2.5 ( PM 2.5,agr_emis ) as the difference in PM 2.5 predicted by simulation (4) and simulation (3), and the impacts of these combined ( PM 2.5,LULCC+agr_emis ) as the difference in PM 2.5 predicted in simulation (4) and simulation (2) (see Table 1). The effect of LULCC on PM 2.5 (Fig. 5a) is mainly through perturbing BVOC emissions as they are a precursor to SOA. Over parts of South America and Southeast Asia, where isoprene emissions drop significantly due to deforestation, PM 2.5 is reduced by up to 0.7 µg m −3 . Land cover changes also lead to changes in the dry deposition velocity of some SNA precursor gases for which stomatal uptake is an important deposition pathway (e.g., NO 2 and SO 2 , Fig. S4). Indeed, over India and China, where our model suggests high levels of SNA aerosol precursors, contemporary LULCC enhances dry deposition of these constituents, which reduces PM 2.5 overall by up to 0.3 µg m −3 , similar to the finding of Fu et al. (2016).
We find that agricultural emissions generally have a larger impact on annual mean surface PM 2.5 level (Fig. 5b) than LULCC. The largest increases in annual mean surface PM 2.5 due to changes in agricultural emissions over 1992 to 2014 occur across China (+0.7 µg m −3 ) and India (+1.6 µg m −3 ). Over some hotspots in the two countries (e.g., northwestern India and the North China Plain), the local changes in   PM 2.5 exceed 3.5 µg m −3 , supporting the previously emphasized importance of controlling NH 3 emissions for PM in terms of air quality in China (Fu et al., 2017), but potentially India as well. Some moderate increases (<2 µg m −3 in most locations) in annual mean PM 2.5 concentrations are also observed in the Middle East, North America, Central America and South America.
The largest decreases (up to 2.1 µg m −3 ) in annual mean PM 2.5 due to changes in agricultural emissions are simulated in central and eastern Europe and the former Soviet Union. Despite comparable reductions in agricultural NH 3 emissions, decreases in PM 2.5 over western Europe are smaller because of weaker sensitivity of SNA aerosol to NH 3 emissions, which is consistent with the finding of Lee et al. (2015) and Pozzer et al. (2017). In general, reductions in annual mean PM 2.5 due to agricultural emission changes simulated over western Europe are weaker than over central and eastern Europe and the former Soviet Union. Figure 5c shows the combined effect of agricultural emissions and LULCC on annual mean surface PM 2.5 , which we have already shown is mostly dominated by the effect of agricultural emissions. Nevertheless, we find that the effects of LULCC are able to partially offset the increase in PM 2.5 due to agricultural emissions changes over China and India. These offsets are occurring in densely populated areas so that the effects on population-weighted average (method described in Text S2 in the Supplement) PM 2.5 concentrations (see below), and therefore potentially exposure, may be noteworthy. This is discussed in further detail below.
We note that the difference between Fig. S5a and 5c illustrates how PM 2.5,LULCC+agr_emis is sensitive to the anthropogenic emissions background. We find that surface PM 2.5 over the US, Europe and the former Soviet Union is less sensitive to NH 3 emissions under the 2014 anthropogenic emissions background, since both SO 2 and NO x emissions in these regions have decreased significantly (>42 % for NO x and >58 % for SO 2 ) over 1992 to 2014. The opposite is simulated over China and India, where SO 2 and NO x emissions have increased by >50 %. Table 3 summarizes the simulated effects of LULCC and agricultural emission changes on PM 2.5 in terms of air quality and compares their magnitudes with the concomitant effects from direct anthropogenic emission changes ( PM 2.5,anth ) over the same time period. We additionally compare area-averaged and population-weighted global and regional metrics. While the resolution of our simulations does not capture urban-scale gradients and nonlinearities in urban chemistry, the use of population weighting allows us to explore whether signals of change in land cover or land management are concentrated over areas of high population or whether they are primarily observed over less populated areas.
Globally, our model results estimate that the global population-weighted change in PM 2.5 resulting from LULCC and agricultural emission changes (+0.70 µg m −3 ) is on the order of ∼ 10 % of the change in PM 2.5 resulting from direct anthropogenic emissions (+7.99 µg m −3 ) over 1992 to 2014. Regionally, the largest impact of land change ( PM 2.5,LULCC+agr_emis ) on population-weighted annual mean surface PM 2.5 is simulated over central and eastern Europe (−1.01 µg m −3 ), the former Soviet Union (−1.00 µg m −3 ), South Asia (+1.71 µg m −3 ), and China (+1.45 µg m −3 ). In most regions, the difference between population-weighted PM 2.5,agr_emis and PM 2.5,LULCC+agr_emis is very small (< ∼ 0.05 µg m −3 ) except in China (0.12 µg m −3 ). Generally, the impacts of land change on population-weighted PM 2.5 have the same sign as the impacts of direct anthropogenic emissions. The only exception to this occurs over North America where anthropogenic NO x and SO 2 emissions have declined, but agricultural emissions have increased. This suggests that the increase in agricultural emissions over North America has partially canceled out the effects of other emission controls on PM 2.5 , though this effect is small so far (∼ 5 %). In other regions, population-weighted PM 2.5,LULCC+agr_emis is generally on the order of 5 % to 12 % of changes due to direct anthropogenic emissions (e.g., in central, eastern and western Europe Europe). Notably, over the former Soviet Union, the Middle East and Central America, PM 2.5,LULCC+agr_emis is much more comparable to the effect of anthropogenic emission changes (24 %, 42 % and 208 %, respectively).
Our result shows that the impact of LULCC and land management changes on PM 2.5 is mainly from agricultural emission changes, while LULCC can result in additional impacts in regions with high SNA precursor emissions (e.g., India, China) through modulating dry deposition. The magnitude of population-weighted PM 2.5,LULCC+agr_emis suggests that land change may contribute significantly to regional and global changes in human PM 2.5 exposure and that the effects of these changes are not isolated to low-population regions. Particularly, over the regions experiencing rapid change in land use intensity (e.g., the former Soviet Union) or slow change in anthropogenic emissions (e.g., Central America, the Middle East), the effects of land changes on particulate air pollution could be comparable (24 % to 208 %) to the effects of direct anthropogenic emission changes.
5 Impact on surface O 3 Figure 6 shows the modeled impacts of LULCC, changes in agricultural emissions and the combined effects of both on annual mean surface O 3 (under 2014 anthropogenic emissions). These changes are calculated identically as for PM 2.5 above: the impact of LULCC on O 3 ( O 3,LULCC ) is the difference in O 3 predicted by simulation (3) and simulation (2), the impact of agricultural emission changes on O 3 ( O 3,agr_emis ) is the difference in PM 2.5 predicted by simulation (4) and simulation (3), and the impact of these combined ( O 3,LULCC +agr_emis) is the difference in PM 2.5 predicted in simulation (4) and simulation (2) (see Table 1). We also use predictions of surface HNO 3 /H 2 O 2 ratios (Fig. S6) as a proxy for VOC-vs. NO x -sensitive chemical O 3 production (Peng et al., 2006;Sillman, 1995) in our discussion of the results.
The modeled response of surface O 3 to LULCC ( O 3,LULCC ) (Fig. 6a) involves several distinct processes (dry deposition, soil NO x and BVOC emissions). Over parts of North America and Central America, the increase in dry deposition velocity (v d ) reduces annual mean surface ozone by up to 0.5 ppbv overall. In central Brazil, deforestation of tropical rainforests leads to a significant reduction in isoprene emissions, reducing surface ozone by up to 0.8 ppbv in this NO x -limited environment (Fig. S6). In contrast, modeled surface ozone decreases by up to 1.2 ppbv further south, where strong increases in LAI lead to increases in v d (up to 0.06 cm s −1 ). The modeled reduction of surface ozone (up to 1 ppbv) over central African rainforests is also likely attributable to increased v d as neither soil NO x nor isoprene emissions change much in the region. However, in other parts of Africa, up to 0.6 ppbv of surface ozone increases are simulated, mainly because of the relatively large increase in soil NO emission. In southern China, up to 0.5 ppbv reduction in surface ozone is simulated, which is likely attributable to the increase in v d and slightly offset by the small increase in isoprene emission under this NO x -saturated environment (Fig. S6). Small surface O 3 changes, mainly due to transport, are also simulated over the Atlantic Ocean.
Overall, the role of agricultural emission changes in fertilizer-associated NO x plays a minor role in surface O 3 changes (Fig. 6b). An exception to this is observed in the large increase in agricultural NO x emissions, which reduces surface O 3 by up to 0.6 ppbv over NO x -saturated India and China but increases surface O 3 in NO x -limited parts of Southeast Asia by a similar magnitude. Slight increases in surface O 3 levels due to increased agricultural NO x emissions are also simulated over parts of eastern Africa and South America. Whether the effect of agricultural emissions strengthens (e.g., China and the Sahel) or offsets (e.g., over Table 3. Changes in area-averaged and population-weighted (in parentheses) annual mean surface PM 2.5 concentrations (in µg m −3 ) due to anthropogenic emissions alone ( PM 2.5(anth) ), LULCC ( PM 2.5(LULCC) ), agricultural emissions ( PM 2.5(agr_emis) ), and the combined effects of LULCC and agricultural emissions ( PM 2.5(LULCC+agr_emis) ) together. Results only from regions with PM 2.5(LULCC) , PM 2.5(agr_emis) or PM 2.5(LULCC+agr_emis) >0.2 µg m −3 are shown.  Table S2. southern Brazil and India) the effect of LULCC is largely region-dependent. As shown in Fig. 6c, LULCC tends to dominate the impacts on surface O 3 over most regions in the world (unlike PM 2.5 for which the effects of agricultural emission changes dominate). Similar to PM 2.5 , we find that the changes in the anthropogenic emission background over 1992 to 2014 are strong enough to alter the sensitivity of O 3 to land change. As indicated by Fig. S6, Asia was less NO x -saturated, while western Europe and the coastal United States were more NO xsaturated in 1992 than in 2014. For example, the increase in soil NO emission over India is more likely to increase rather than decrease the surface ozone concentration (Fig. S7a), leading to different modeled effects on surface ozone. Table 4 shows the change in area-and populationweighted annual mean afternoon surface O 3 due to the effects of anthropogenic emissions ( O 3,anth , Fig. S7b), O 3,LULCC , O 3,agr_emis and O 3,LULCC+agr_emis . In most regions, O 3,LULCC+agr_emis is positive. However, this is offset by the negative population-weighted average O 3,LULCC+agr_emis over the most populous regions (South Asia and China), resulting in very small globally averaged population-weighted O 3,LULCC+agr_emis .

Region
The magnitudes of population-weighted O 3 (within ±0.5 ppbv) display less regional variability than that of PM 2.5 . Over eastern Africa, western Africa and southern Africa, area-averaged O 3,LULCC+agr_emis generally has magnitudes similar to population-weighted O 3,LULCC+agr_emis . In other regions, the differences between area-and population-weighted O 3,LULCC+agr_emis are more substantial. The largest discrepancies between areaand population-weighted O 3,LULCC+agr_emis are found over China, where increases in surface O 3 are predicted over less populated western China, while reductions in surface O 3 are simulated over more densely populated eastern China. In South America, there are large sub-regional signals of O 3,LULCC+agr_emis , but these positive and negative signals largely offset each other, resulting in both small areaweighted and population-weighted O 3,LULCC+agr_emis .
Over China, western Africa, eastern Africa, southern Africa, the former Soviet Union and the Middle East, the magnitudes of population-weighted O 3,LULCC+agr_emis are more than 20 % of that of O 3,anth , implying that contemporary land system changes could be a regionally important component in contemporary trends of surface O 3 . The effects of agricultural emission changes and LULCC can either noticeably enhance (e.g., over the Middle East, Japan and Korea, China) or offset (e.g., over South Asia) each other because of the dependence of O 3,agr_emis+land_cover on regional NO x -VOC chemistry and details of LULCC, indicating the complexity of diagnosing the effect of land change on surface O 3 at regional and global scale.
Our result suggests that contemporary agricultural emission changes and LULCC each have distinct effects on surface O 3 , with LULCC generally stronger in magnitude. Both of the effects are dependent on local NO x -VOC chemistry, as agricultural emission changes perturb NO x emissions, while LULCC tends to affect BVOC emissions. In addition, LULCC is also able to affect surface O 3 (and other precursors) directly through dry deposition and LAI changes over our period of concern. These effects are found to affect O 3 pollution over densely populated regions (e.g., China) and could be comparable to the magnitudes of O 3 changes due to anthropogenic emissions over specific regions (e.g., the former Soviet Union, eastern Africa, western Africa), indicating the importance of land change in studying long-term changes in surface O 3 .

Impact on nitrogen deposition
Finally, we estimate the effect of these land changes on nitrogen deposition estimates. Figure 7 shows the global Table 4. Changes in area-averaged and population-weighted (in parentheses) annual mean surface O 3 concentrations (in ppbv) due to anthropogenic emissions alone ( PM 2.5(anth) ), LULCC ( PM 2.5(LULCC) ), agricultural emissions ( PM 2.5(agr_emis) ), and the combined effects of LULCC and agricultural emissions ( PM 2.5(LULCC+agr_emis) ) together. Results only from regions with population-weighted average O 3(LULCC) , O 3(agr_emis) or O 3(LULCC+agr_emis) >0.2 ppb are shown.   impact of LULCC and agricultural emission changes on total nitrogen deposition ( N dep ), and Table 5 summarizes the regional and global results. The largest increase and decrease in nitrogen deposition (N dep ) are simulated over South Asia (+1.91 Tg N yr −1 ) and the former Soviet Union (−1.28 Tg N yr −1 ), respectively. Notable increases in N dep are also simulated over China (+1.55 Tg N yr −1 ), South America (+1.24 Tg N yr −1 ), North America (+0.66 Tg N yr −1 ), western Africa (+0.39 Tg N yr −1 ) and eastern Africa (+0.41 Tg N yr −1 ). Figure 7 also illustrates the simulated changes over 1992 to 2014 in areas with nitrogen deposition (N dep ) exceeding 5 kg N ha −1 yr −1 , which is a proxy for possible exceedance of critical N dep loads for terrestrial and fresh water (Moriarty, 1988). Globally, there is a net increase in land area with N dep >5 kg N ha −1 yr −1 of 3.67 × 10 6 km 2 . The increase is mostly simulated over the Americas, Africa, the Middle East and China, which is partially offset by the large de- Table 5. Changes in total nitrogen deposition ( N dep ) and land area that has nitrogen deposition >5 kg N ha −1 yr −1 ( Area crit ), which is a proxy for potential risk of critical nitrogen deposition load exceedance. Only regions with significant N dep (> 0.25 Tg N yr −1 ) or Area crit (> 10 5 km 2 ) are shown. +7.20 +3673 * The definitions and abbreviations of all regions can be found in Table S2. crease over the former Soviet Union. Meanwhile, despite agricultural changes that lead to notable N dep , over most of Europe, the eastern US, China, South Asia and Southeast Asia, nitrogen input from other sources is large enough that this signal alone does not lead to substantial changes in N dep exceedances of 5 kg N ha −1 yr −1 . However, over parts of North America, South America, Africa and China, agricultural changes are simulated to increase N dep from below to above 5 kg N ha −1 yr −1 . This implies that these natural ecosystems at the edge of these areas are at risk of nitrogen exceedances due to agricultural changes. In contrast, the substantial reduction of N dep in parts of the former Soviet Union may have significantly reduced the risk of nitrogen exceedance in natural ecosystems from agricultural sources.

Discussion and conclusions
In this work, we have explored how changes in the global land system, through LULCC and agricultural emission changes, may have impacted contemporary global air quality over 1992 to 2014. We model the effects of contemporary LULCC and agricultural emission changes, individually then in combination, on surface O 3 and PM 2.5 using the GEOS-Chem chemical transport model (CTM). With a uniquely consistent framework, we are able to integrate direct information from global emission inventories (CEDS) with updated land surface remote sensing products (ESA CCI land cover and GLASS LAI). This allows us to avoid invoking extra assumptions on land management practices (e.g., constant N r input, emissions or emission factors over time) and biophysical properties of PFTs (e.g., constant PFT-specific LAI over time).
We find that changes in agricultural emissions are simulated to increase the annual mean surface PM 2.5 concentrations in China and India by up to 3.5 µg m −3 and to decrease in Europe by up to 3.5 µg m −3 . Our simulation suggests that though PM 2.5 is mainly attributable to changes in agricultural emissions at the global scale, LULCC over India and China can lead to enhanced dry deposition of certain PM 2.5 precursor gases (SO 2 and NO 2 ), thus partially offsetting (∼ 10 %) the increase in PM 2.5 from agricultural regions. This implies a potentially important role of LULCC in determining the SNA aerosol level over certain heavily polluted regions. Also, LULCC reduces BVOC emissions over Amazonia, which leads to reductions in PM 2.5 by up to 0.7 µg m −3 . In a future with decreasing anthropogenic NO x and SO 2 emissions, which could diminish the importance of agricultural emissions for PM 2.5 formation (Bauer et al., 2016), LULCC may become increasingly important in the overall effect of land change on PM 2.5 . Noticeable changes (>1 µg m −3 ) in populationweighted PM 2.5,LULCC+agr_emis are simulated over China (+1.45 µg m −3 ), South Asia (+1.71 µg m −3 ), central and eastern Europe (−1.00 µg m −3 ), and the former Soviet Union (−1.01 µg m −3 ), indicating the potential impact of land change on long-term public health through modulating the PM 2.5 level at regional scale. Our results suggest that contemporary (1996-2014) LULCC and agricultural emission changes contribute to changes in PM 2.5 at regional and global scales that range from the order of 5 % to 10 % of changes in PM 2.5 resulting from direct anthropogenic emissions over the same time period and up to ∼ 25 % or more in the former Soviet Union and the Middle East specifically.
In contrast, the effect of LULCC is generally stronger than that of agricultural emission change in simulations of surface O 3 . We find that the role of LULCC over 1992 to 2014 is regionally significant enough to induce changes in BVOC emissions and dry deposition, which affect surface O 3 , but that the overall effects largely offset each other on the global scale, leading to very small population-weighted O 3,LULCC+agr_emis . This finding is consistent with that of Ganzeveld et al. (2010), even though the timeframe of study  is different. The effects of both agricultural emission changes and LULCC, through NO x and BVOC emissions, are sensitive to the regional ozone production regime. The increase in agricultural emissions reduces O 3 over NO x -saturated parts of China and South Asia by up to 0.6 ppbv, while the reduction in BVOC emissions increases surface O 3 over VOC-limited Amazonia by up to 1.2 ppbv; enhancements of dry deposition reduce O 3 over parts of China, North America and South America by up to 1.2 ppbv. Overall, the largest populationweighted O 3,LULCC+agr_emis is simulated over western Africa (+0.42 ppbv) and eastern Africa (+0.47 ppbv). We find that the ratio between O 3,LULCC+agr_emis and O 3,anth varies widely depending on region, with some having O 3,LULCC+agr_emis comparable (>20 %) to O 3,anth . These results show the complexity and importance of land change in mediating long-term changes in surface O 3 . We also find that both the modeled O 3,LULCC+agr_emis and PM 2.5,LULCC+agr_emis are sensitive to the changes in anthropogenic emissions suggested by the CEDS inventory over 1992 to 2014, as the changes in NO x , SO 2 and VOC emissions are large enough to considerably perturb atmospheric HNO 3 and H 2 SO 4 production, as well as the ozone production regime in many regions (e.g., Asia and western Europe). This highlights the necessity of accurate and relevant emission inventories when evaluating the impacts of land change on air quality (e.g., Bauer et al., 2016).
The increased atmospheric reactive nitrogen (+7.20 Tg yr −1 ) due to agricultural emissions is mostly found to be deposited near source regions as the atmospheric lifetime of NH 3 is generally short, which implies the potential risk of excessive nitrogen input over natural ecosystems near regions with increases in agricultural emissions.
Our work suggests that, at contemporary timescales (on the order of ∼ 20 years), the effect of land change on air quality can sometimes be important relative to the air quality changes induced by trends in direct anthropogenic emissions. We also find that agricultural emission changes have stronger effects on PM 2.5 , while LULCC has stronger effects on O 3 . This finding is comparable to that from Heald and Geddes (2016), which suggests much more comparable changes in biogenic SOA (mostly induced by LULCC) and particulate nitrate (mostly induced by agricultural emission changes), as well as stronger surface ozone changes induced by land change over 1850-2000. This shows that both the magnitudes and relative contributions from different components of land change effects on air quality vary significantly with the timescale of study, as well as the potential importance at longer timescales (e.g., multidecadal, centennial), despite the relatively small signal that we obtain here.
We find the effects of agricultural emissions and LULCC to be largely linearly additive over contemporary timescales, which may be attributable to two factors: (1) LULCC mainly impacts O 3 precursors, while agricultural emissions mainly impact SNA precursors, and these are often spatially segregated; (2) LULCC and agriculture-related changes in surface fluxes of O 3 and SNA precursors are not large enough to change their respective chemical production regime. At longer timescales when land change signals are stronger, the effects of LULCC and agricultural emissions may be nonlinear.
We note several important limitations and opportunities for development. We were only able to evaluate our simulation extensively over Europe, North American and East Asia. In most other regions where such evaluation of SNA speciation is not feasible, the sensitivity of SNA formation to NH 3 emissions can be a major source of uncertainty. Given that the changes in agricultural emissions have occurred at a global scale, effort to monitor SNA speciation outside North America and Europe (e.g., Weagle et al., 2018) is neces-sary for understanding the sensitivity of PM 2.5 to agricultural emissions to a global extent. A better understanding of both the sources and sinks of HNO 3 (e.g., Heald et al., 2012;Holmes et al., 2019;Luo et al., 2019;Petetin et al., 2016) as well as nitrate partitioning (e.g., Vasilakos et al., 2018) is important for modeling SNA aerosol and its sensitivity to NH 3 emissions. Agricultural NO x and NH 3 emissions estimates also carry large uncertainty due their biological nature and resulting dependence on environmental conditions, which are not explicitly considered in the construction of bottom-up anthropogenic emission inventories (Crippa et al., 2018;Hoesly et al., 2018). Bidirectional exchanges of NO 2 (Breuninger et al., 2013;Chaparro-Suarez et al., 2011;Lerdau et al., 2000) and NH 3 (Bash et al., 2013;Massad et al., 2010;Wichink Kruit et al., 2012;Zhang et al., 2010) are not explicitly modeled (although in some regions they may be implicitly accounted for in the regional scaling performed by CEDS), which introduces some uncertainty in the accuracy of surface flux modeling. Zhu et al. (2015) implemented a bidirectional NH 3 exchange model in GEOS-Chem and found no substantial improvement with observations in the modeled NH 3 concentration, NH + 4 wet deposition and nitrate aerosol concentration compared to the default GEOS-Chem unidirectional exchange framework. This indicates that the unidirectional framework may still be sufficiently accurate in simulating global air quality compared to the bidirectional framework, which requires more observations to properly parameterize at global scale. In the case of NO 2 , we make the assumption that in most regions we are interested in (Fig. S9), the ambient concentrations of NO 2 exceed an ecosystem compensation point (0.05-0.6 ppb) (e.g., Breuninger et al., 2013) so that we can assume deposition would dominate. The simplistic representation of dry deposition in general, particularly the lack of dependence of stomatal conductance on atmospheric and soil water content, may not adequately capture the effects of LULCC, as biomes can have differential responses to meteorological and hydrological conditions. The inherent inconsistency of long-term LAI time series derived from reflectance measured by different instruments (Jiang et al., 2017) and the use of static land cover maps also introduce uncertainty in the LAI retrieval (Fang et al., 2013) and the subsequently computed LAI changes and trends, and these have been shown to be important to changes in simulated O 3 in this study and elsewhere (Wong et al., 2019). Though the use of PFT-based emission factors in regional and global modeling is generally justifiable (Guenther et al., 2012), we cannot rule out the possibility of intra-PFT variabilities of BVOC emission factors affecting the accuracies our results, which is exemplified by the inability of our model to capture the palm-driven isoprene emission increase over Southeast Asia  as discussed in Sect. 3. Finally, the meteorological feedbacks (e.g., changes in sensible heat, latent heat, air temperature, boundary layer height) and the subsequent effects on atmospheric chemistry and transport from LULCC and agricul-tural emissions are not considered in our study, which could potentially be important (e.g., Wang et al., 2020).
Our study helps demonstrate the possible magnitudes and regional patterns of the impacts of contemporary LULCC and agricultural emission changes on PM 2.5 and O 3 , and it suggests that the combination of these factors should not be neglected in the study of regional and global air quality changes over multidecadal timescales. Our results confirm the potential importance of controlling agricultural emissions for improving air quality in terms of PM 2.5 , which could be practical as there are numerous feasible options for reducing agricultural emissions through optimizing livestock and crop production systems (e.g., Ti et al., 2019). Incentivizing these and other practices that improve agricultural nitrogen use efficiency (e.g., including livestock production with cropping, synchronizing nitrogen supply with crop demand) (e.g., Fageria and Baligar, 2005;Langholtz et al., 2021) can be one of the keys to mitigate the air quality impacts of reactive nitrogen input without compromising agricultural productivity (e.g., Guo et al., 2020). Furthermore, as increasing reactive nitrogen input and land use change are the two of the main strategies to meet the global demand for biomassbased products in the future (Foley et al., 2011), the distinct yet significant impacts of agricultural emissions and land use change on O 3 , PM 2.5 and nitrogen deposition should be investigated as part of the overall environmental impacts of land system changes, especially when there is a trade-off between increasing land input and cropland expansion (e.g., Lotze-Campen et al., 2010;Mauser et al., 2015). This could benefit agricultural policy activities by appropriately considering the externalities and socioeconomic costs of different options and scenarios for agricultural expansion.
Code availability. The source code of the GEOS-Chem model is publicly available (https://doi.org/10.5281/zenodo.3634864, the International GEOS-Chem User Community, 2020). The GEOS-Chem model output and other source code used in the project can be obtained by contacting the corresponding author (jgeddes@bu.edu).
Author contributions. AYHW and JAG developed the ideas for this study, formulated the methods and designed the model experiments together. AYHW performed the chemical transport model simulations and data analysis, with input and feedback from JAG. Paper preparation was performed by AYHW, and the paper was reviewed, edited and approved by JAG.
Competing interests. The authors declare that they have no conflict of interest.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.