Rising importance of agricultural nitrogen oxide emissions in China’s future PM2.5 pollution mitigation

Rising importance of agricultural nitrogen oxide emissions in China’s future PM2.5 pollution mitigation

Introduction

Fine particulate matter (PM2.5), which describes inhalable particles with an aerodynamic diameter equal to or less than 2.5 μm, is a major pollutant that has been proven to increase susceptibility to respiratory diseases1. In China, severe PM2.5 pollution contributed to 1.1 million excess deaths in 20152. The anthropogenic emission control strategies implemented by the Chinese government, such as the ‘Three-year Action Plan Fighting for a Bule Sky’ issued in 20183, have significantly reduced PM2.5 pollution4,5. However, most of the Chinese population is still exposed to PM2.5 levels much higher than that recommended by the World Health Organization’s Air Quality Guidelines (AQG, annual mean PM2.5 of 5 μg m−3)6.

As an important PM2.5 precursor, the nitrogen oxides (NOx = NO + NO2) emissions from combustion sources have been effectively reduced with the implementation of clear air actions in China7,8. Besides combustion sources, agricultural soil is also recognized as an important NOx source and can have a significant influence on ozone air quality9,10,11, even in regions that already have a high level of NOx emissions from fuel combustion, such as North China Plain (NCP)12,13,14. Although the NOx emissions from agricultural soils are mainly induced by agricultural fertilizer application14,15, they are usually treated as natural emissions. Recent studies suggest that reducing agricultural NH3 emissions through agricultural nitrogen management is needed to further control PM2.5 pollution in China15,16,17,18,19. Agricultural nitrogen management can also reduce NOx emissions20. However, the influence of controlling agricultural NOx emissions on mitigating PM2.5 pollution is rarely recognized in air quality strategies.

The agricultural emissions of NH3 and NOx contribute to the formation of secondary inorganic aerosols (sulfate, nitrate, and ammonium aerosols, SIAs) through nonlinear chemistry. SO2 and NOx emitted to the atmosphere can be oxidized to sulfuric acid (H2SO4) and nitric acid (HNO3), respectively. NH3 first forms ammonium sulfate aerosol by reacting with H2SO4, and then the excess NH3 will form ammonium nitrate aerosol by reacting with HNO3. When NH3 is excessive, which is commonly found in agricultural-intensive areas, controlling NH3 emissions will first reduce the NH3 concentrations rather than SIAs21. As the SIA formations are more sensitive to HNO3 than NH3 over the NH3-excessive regions, we propose here the effectiveness of controlling agricultural emissions on PM2.5 reduction might be underestimated if agricultural NOx emissions are not considered.

In this study, we investigate the role of agricultural NOx emissions in mitigating PM2.5 pollution during the crop-growing season over eastern China at present and under future scenarios by conducting a set of atmospheric model simulations. We show that over the North China Plain, including agricultural NOx emissions will shift the optimal reactive nitrogen (Nr) control pathway toward controlling more agricultural emissions and will achieve more effective PM2.5 mitigation. Given that most agricultural nitrogen management practices aimed at controlling NH₃ emissions also regulate agricultural NOₓ emissions, changes in agricultural NOx emissions associated with these agricultural management strategies should also be considered when developing emission control strategies for intensive agricultural areas like North China Plain.

Results

Anthropogenic NOx emissions from agricultural and combustion sectors

Figure 1a and Supplementary Fig. 1 show the spatial distribution and seasonal variation of agricultural NOx emissions (ENOx_agri) over China in 2019, as estimated by the Berkeley–Dalhousie Soil NOx Parameterization (BDSNP)22 with fertilizer application provided by Chen et al.23 (Supplementary Fig. 2; Methods). The annual total Chinese ENOx_agri estimated by BDSNP is 0.3 Tg N yr−1, consistent with the ENOx_agri given by the sixth Coupled Model Intercomparison Project (CMIP6) scenarios (also 0.3 Tg N yr−1, Methods). The highest agricultural NOx emissions are found over the North China Plain (NCP), reaching 9 kg N ha−1 yr−1, and contributing 35% of the total Chinese agricultural NOx emission in 2019. In other regions of China, agricultural NOx emissions are typically less than 1 kg N ha−1 yr−1. Significant seasonal variations exist, with the ENOx_agri mainly occurring over May-August, which accounts for 65% of the annual total ENOx_agri (Supplementary Fig. 1). Higher ENOx_agri during May-August can be attributed to high temperature and intensive fertilizer application during these months10,24.

Fig. 1: The spatial distribution and future changes in agricultural NOx emissions in China.
Rising importance of agricultural nitrogen oxide emissions in China’s future PM2.5 pollution mitigation

a The spatial distribution of agricultural NOx emissions in China was estimated by BDSNP in 2019. The total agricultural NOx emission over China is shown inset. b, c The bar plots show the anthropogenic NOx emissions from each emission sector in 2019 and 2030–2060 under SSP2-45 (b) and SSP1-26 (c) scenarios (left axis). The lines show the percentage contribution of agricultural NOx emissions to annual (gray) and May-August (black) total anthropogenic NOx emissions (right axis).

Full size image

We evaluate the contribution of ENOx_agri to total anthropogenic NOx emissions in China in 2019, and the future by comparing with NOx emissions from combustion sectors (ENOx_comb, including industrial, power plant, transportation, and residential burning) in Fig. 1. The ENOx_comb in 2019 is derived from the Multi-resolution Emission Inventory for China (MEIC)7,25,26, and in the future (Fig. 1) derived by applying the anthropogenic emission changes (including the agricultural sector) given in the CMIP6 scenarios (Methods). In 2019, the ENOx_agri (0.3 Tg N yr−1) contributes to 5% of annual total anthropogenic NOx emissions (6.6 Tg N yr−1). In the future, the emission mitigation scenarios in CMIP6 (SSP1-26 and SSP2-45) predict substantial decreases in ENOx_comb but only slight changes in ENOx_agri. The SSP1-26 scenario (strict mitigation pathway scenario), which predicts a sharp reduction of NOx emissions from industry and transportation, shows ENOx_agri accounting for 11% of the annual total anthropogenic NOx emissions in 2040 and 27% in 2060. The SSP2-45 scenario (middle road scenario), predicts lower decreases in industry and transportation NOx emissions compared to the SSP1-26 scenario, resulting in ENOx_agri contributions of 7% in 2040 and 10% in 2060.

The contribution of agricultural NOx emissions is more significant during the crop-growing season (May-August) (Fig.1 and Supplementary Fig. 3). In 2019, the contribution from ENOx_agri (53 Gg N month−1) reached 10% of total anthropogenic NOx emissions during May-August, which is significantly higher than the residential sector. The ENOx_agri would contribute 43% and 19% of the total anthropogenic NOx emissions under SSP1-26 and SSP2-45 scenarios in 2060, respectively. Remarkably, under the SSP1-26 scenario, ENOx_agri (50 Gg N month−1) during May-August in 2060 would be even higher than the transportation (19 Gg N month−1) and industry (31 Gg N month−1) sectors, making ENOx_agri a primary contributor to anthropogenic NOx emissions in China (Supplementary Fig. 3).

The contribution of agricultural NOx emissions to PM2.5 pollution

We use the atmospheric chemistry transport model GEOS-Chem to simulate the PM2.5 components in 2019 and the future under SSP2-45 and SSP1-26 scenarios during crop-growing season over East Asia (Methods). The base simulation can generally capture the spatial distribution of the PM2.5 concentrations in 2019 observed by the surface network (Supplementary Fig. 4). Our analysis mainly focuses on the SIAs, which are directly impacted by NOx emissions. The contributions of ENOx_agri and ENOx_comb to SIAs are calculated as the difference between model simulations with all emissions turned on (base simulations in Supplementary Table 1) and simulations with corresponding emissions turned off (sensitivity simulations Group I and II in Supplementary Table 1).

Figure 2 shows the spatial distribution of simulated surface SIA concentrations during the crop-growing season (May-August) of 2019 and under future emission scenarios. In 2019, the model simulated SIA concentrations are more than 20 μg m−3 over the majority of eastern China except for northeast China and southern coastal regions (5–15 μg m−3). With the reduction of NOx (Fig. 1) and SO2 emissions (Supplementary Fig. 5) from combustion sectors, the model simulated significant decreases in SIA concentrations in the future. However, despite these reductions, the formation of SIAs still hinders eastern China from meeting the AQG level (annual mean PM2.5 of 5 μg m−3)6. Under the SSP2-45 scenario, the SIA concentrations are more than 20 μg m−3 in 2040 and 10 μg m−3 in 2060 over the NCP, implying that SIAs alone (excluding primary aerosols and secondary organic aerosols) could result in PM2.5 levels exceeding the AQG level. Under the SSP1-26 scenario, the SIA concentrations in 2060 would decrease to 2–5 μg m−3 over most regions of eastern China. Given that SIAs account for 30-50% of PM2.5 in eastern China27,28,29, further reducing SIAs is essential to ensure PM2.5 concentrations below the AQG level.

Fig. 2: The SIA concentrations over eastern China and contribution of agricultural and combustion NOx emissions.
figure 2

ae Spatial distribution of SIA concentrations over eastern China simulated by GEOS-Chem averaged over the crop-growing season (May-August) in 2019 (a) and under SSP2-45 and SSP1-26 emission scenarios for 2040 and 2060 (be). fj The corresponding contribution of agricultural NOx emissions (ENOx_agri) and combustion NOx emissions (ENOx_comb) to SIA concentrations over East China (EC, red) and North China Plain (NCP, yellow). The bar plots show the median values, and the error bars show the ranges.

Full size image

We compare the contribution of ENOx_agri and ENOx_comb to SIA concentrations over NCP and the whole of eastern China (EC). The spatial distributions of the contribution of ENOx_agri and ENOx_comb to SIA concentrations are shown in Supplementary Fig. 6, and the median values and ranges over NCP and EC are summarized in Fig. 2. The contribution of ENOx_agri to SIA concentrations in 2019 during the crop-growing season is usually below 10% across most regions in eastern China. Over the NCP, where the highest ENOx_agri levels occur, the contribution to surface SIA concentrations reaches about 10-20%. In comparison, the contribution of ENOx_comb can reach more than 50% over most regions in eastern China, much larger than the contribution from ENOx_agri. However, model simulations for future emission conditions indicate a significant increase in the contribution of ENOx_agri to SIAs, particularly over the NCP. In 2060, the ENOx_agri contributions will further increase to 15% over most of eastern China. Over the NCP, the contribution of ENOx_agri will exceed 20% under the SSP2-45 scenario and 40% under the SSP1-26 scenario. Under the SSP1-26 scenario, ENOx_agri contributions will be even larger than those from total ENOx_comb. Additionally, we compare the contributions of ENOx_agri with NOx emissions from the industrial sector (ENOx_ind), which is the largest combustion source in China (Fig. 1). As shown in Supplementary Fig. 7, under the SSP1-26 scenario, the contribution of ENOx_agri will become roughly equal to that of ENOx_ind from 2040 onwards across eastern China. Over the NCP, the contribution of ENOx_agri will be twice that of ENOx_ind in 2060. Under the SSP2-45 scenario, where industrial emissions are less reduced, the contribution of ENOx_agri will exceed half of ENOx_ind in 2060 over the NCP. Thus, as more stringent measures are expected to be implemented to reduce the emissions from fuel combustion in the future, there will be a growing significance of ENOx_agri in the formation of SIAs. ENOx_agri will serve as a pivotal factor in SIA formations during crop-growing season over eastern China.

The role of controlling agricultural NOx emissions on PM2.5 pollution mitigation

Despite the future reduction of emissions from fossil fuel combustion, eastern China still experiences a substantial level of SIA formations. Since the agricultural emissions show slight changes under the future emission scenarios (Fig. 1 and Supplementary Fig. 8), controlling agricultural emissions could be crucial for further mitigating PM2.5 pollution and achieving more stringent air quality targets such as the AQG level6. We further evaluate the efficacy of controlling ENOx_agri on PM2.5 pollution by comparing the reduction of SIA concentrations resulting from agricultural NH3 and NOx emission controls.

Agricultural nitrogen management practices, such as using enhanced-efficiency fertilizers, reducing fertilizer application rate, machine application, and improving manure management, have the potential to reduce NH3 emissions by ~40%20,30,31. A recent study32 also reported that NH3 emissions could be mitigated by 49% for wheat and 39% for maize under a county-wide NH3 reduction campaign in the NCP. In addition to NH3, agricultural NOx emissions can also be reduced by about 40% through improved nitrogen management practices20. We posit that this potential for a 40% reduction in agricultural NH3 and NOx will persist in the future, given the relatively small changes in agricultural emissions predicted by future scenarios (Fig. 1 and Supplementary Fig. 8). Here, we focus our analyses on May as a representative month of the crop-growing season and the results for June-August are shown in the Supplementary Materials.

Figure 3 shows the decline in SIA concentrations resulting from a 40% reduction of agricultural emissions over NCP. The results for the entire eastern China are shown in Supplementary Fig. 9. In May 2019, reducing 40% of agricultural NH3 emissions would lead to an up to 4.6 μg m−3 decrease in SIA concentrations over eastern China (Supplementary Fig. 9). Over the NCP, where the agricultural emissions are most intensive, the SIA concentrations only decrease by 0.5–3.3 μg m−3 (Fig. 3). We estimate the G ratio33 (Supplementary Text 1) to describe the chemical regime of SIA formations, of which the value large than 1 indicates excessive NH3 for SIA formations (Supplementary Fig. 10). The limited response of SIA concentrations to NH3 emission reductions over NCP can be attributed to highly saturated NH3 simulated by GEOS-Chem sensitivity simulations, where the G values are larger than 1 even with a 40% reduction of ENH3_agri. In the future, the CMIP6 scenarios predict minor changes in agricultural NH3 emissions (Supplementary Fig. 8). NH3 will be more excessive in eastern China with reductions of NOx and SO2 from fuel combustion under future scenarios (Supplementary Fig. 10), leading to a significant decline in the SIA concentration decreases associated with a 40% reduction in ENH3_agri (Fig. 3). Both the SSP2-45 and SSP1-26 scenarios predict that the corresponding mean SIA decreases in 2060 would be less than 1.0 μg m−3 over NCP (Fig. 3).

Fig. 3: The effectiveness of controlling agricultural NOx and NH3 emission in SIA mitigation.
figure 3

a The decrease in SIA concentrations in response to a 40% reduction in agricultural NH3 emissions (gray) and a 40% reduction in both agricultural NH3 and NOx emissions (blue) over the North China Plain in May 2019 and under SSP2-45 and SSP1-26 emission scenarios. The solid dots are the median values, and the error bars are the ranges over the North China Plain. The pie plots show the corresponding contribution of population-weighted SIA concentration reductions led by a 40% reduction of agricultural NOx (yellow) and NH3 emissions (gray). b The efficiency of controlling agricultural NH3 emissions (({beta }_{{{mathrm{NH}}}3}), gray) and agricultural NOx emissions (({beta }_{{NOx}}), yellow) in reducing SIA concentrations in May 2019 and under SSP2-45 (left) and SSP1-26 (right) emission scenarios for 2040 and 2060. The error bars show the ranges estimated by assuming a factor of 2 (50–200%) uncertainties in agricultural NOx emissions. c The same as (b) but for absolute efficiencies (({beta }_{{{mathrm{NH}}}3}^{{prime} }) and ({beta }_{{{mathrm{NOx}}}}^{{prime} })).

Full size image

When both ENH3_agri and ENOx_agri are concurrently reduced by 40%, the response of SIA concentrations to agricultural emission reductions will be more significant than reducing ENH3_agri alone, especially over the NCP (Fig. 3 and Supplementary Fig. 9). This reflects a substantial contribution from ENOx_agri reduction on mitigating SIA formations. In May 2019, ENOx_agri reduction contributed 15% of the overall decrease in SIAs resulting from the total reduction in agricultural Nr (both NH3 and NOx) emissions averaged over the NCP (Fig. 3). Under future scenarios, there will be a noticeable increase in the contribution of ENOx_agri (Fig. 3 and Supplementary Fig. 9). By 2060, the contribution of ENOx_agri reduction to the decrease in SIA concentrations led by total agricultural emission reduction is 34% under the SSP2-45 scenario and 46% under the SSP1-26 scenario averaged over the NCP, implying that reducing ENOx_agri by 40% can lead to similar SIA mitigation as ENH3_agri. Over most of the other regions in eastern China, the contribution of ENOx_agri also increases from around 5–15% in 2019 to 20–30% and 30–40% in 2060, respectively, under SSP2-45 and SSP1-26 scenarios (Supplementary Fig. 9).

The rising importance of ENOx_agri can be further quantified using the efficiency of controlling ENH3_agri (({beta }_{{{mathrm{NH}}}3})) and ENOx_agri (({beta }_{{{mathrm{NOx}}}})) to mitigate SIA formations over NCP21,34. Efficiencies ({beta }_{{{mathrm{NH}}}3}) and ({beta }_{{{mathrm{NOx}}}}) are calculated as the ratio of relative changes in population-weighted regional mean SIA concentrations ((delta {C}_{{{mathrm{SIA}}}})) over the North China Plain to the relative changes in ENH3_agri (left({beta }_{{{mathrm{NH}}}3}=left(frac{delta {C}_{{{mathrm{SIA}}}}}{{C}_{{{mathrm{SIA}}}}}right)/left(frac{delta {E}_{{{mathrm{NH}}}3{_agri}}}{{E}_{{{mathrm{NH}}}3{_agri}}}right)right)) and ENOx_agri (left({beta }_{{{mathrm{NOx}}}}=left(frac{delta {C}_{{{mathrm{SIA}}}}}{{C}_{{{mathrm{SIA}}}}}right)/left(frac{delta {E}_{{{mathrm{NOx}}_agri}}}{{E}_{{{mathrm{NOx}}_agri}}}right)right)), respectively. As shown in Fig. 3, in May 2019, ({beta }_{{{mathrm{NH}}}3}) is 0.35%/%, which is about five times of ({beta }_{{{mathrm{NOx}}}}) (0.07%/%). Projecting to future emission scenarios, ({beta }_{{{mathrm{NH}}}3}) shows only slight variations, remaining around 0.35%/%. Different from ({beta }_{{{mathrm{NH}}}3}), ({beta }_{{{mathrm{NOx}}}}) increases with the reduction of NOx and SO2 emissions from combustion sources. In 2060, the ({beta }_{{{mathrm{NOx}}}}) 0.18%/% (0.10–0.29%/% considering a factor of two uncertainty in agricultural NOx emissions) under SSP2-45 and 0.36%/% (0.21–0.55%/%) under SSP1-26, become comparable to ({beta }_{{{mathrm{NH}}}3}).

We also calculate the absolute efficiency ({beta }_{{NH}3}^{{prime} }) and ({beta }_{{NOx}}^{{prime} }) (({beta }_{{{mathrm{NH}}}3}^{{prime} }=delta {C}_{{{mathrm{SIA}}}}/delta {E}_{{{mathrm{NH}}}3{_agri}}), ({beta }_{{{mathrm{NOx}}}}^{{prime} }=delta {C}_{{{mathrm{SIA}}}}/delta {E}_{{{mathrm{NOx}}_agri}})) over NCP in Fig. 3, which describe the SIA changes in response to unit mass changes in agricultural NH3 and NOx emissions. ({beta }_{{{mathrm{NH}}}3}^{{prime} }) is much smaller than ({beta }_{{{mathrm{NOx}}}}^{{prime} }), reflecting the lower sensitivity of SIAs to NH3 emissions due to the high saturation of NH3 over the North China Plain. However, the relative efficiency ({beta }_{{{mathrm{NH}}}3}) is higher than ({beta }_{{NOx}}), as agricultural NH3 emissions are much larger in magnitude compared to NOx. The decrease of ({beta }_{{{mathrm{NH}}}3}^{{prime} }) reflect that as combustion emissions are reduced, NH3 becomes more excessive for SIA formation. In 2060, ({beta }_{{{mathrm{NOx}}}}^{{prime} }) is 0.089 μg m−3/Gg N (0.085–0.090 μg m−3/Gg N) and 0.062 μg m−3/Gg N (0.059–0.065 μg m−3/Gg N) for the month of May under SSP2-45 and SSP1-26 scenarios, respectively. In comparison, ({beta }_{{{mathrm{NH}}}3}^{{prime} }) is only 0.008 μg m−3/Gg N under SSP2-45 and 0.004 μg m−3/Gg N under SSP1-26. The model simulations for other months of the crop-growing season under the SSP1-26 scenario in 2060 also show comparable ({beta }_{{{mathrm{NH}}}3}) and ({beta }_{{{mathrm{NOx}}}}) and much higher ({beta }_{{{mathrm{NOx}}}}^{{prime} }) than ({beta }_{{{mathrm{NH}}}3}^{{prime} }) (Supplementary Fig. 11). The results indicate that during the crop-growing season, as NH3 is saturated for future SIA formations over most of eastern China, the SIA concentrations are projected to be much more sensitive to per unit reduction of ENOx_agri than ENH3_agri.

Implications for future emission control strategies

Choosing an optimal pathway for controlling combustion and agricultural Nr emissions can help achieve the most efficient PM2.5 mitigation34. As aforementioned, the ENOx_agri are usually treated as natural sources and are not considered in emission control strategies. However, our model sensitivity simulations show that in the future, controlling ENOx_agri can lead to similar SIA reductions to that controlling ENH3_agri. The results imply that ENOx_agri should not be neglected when designing Nr control strategies. Here we further examine the influence of ENOx_agri on the design of efficient emission control pathways for combustion and agricultural reactive nitrogen emissions.

To estimate the optimal Nr emission control pathways, we conducted sensitivity simulations Group V-VII, which gradually reduced Nr emissions from agricultural and combustion sources by 20, 40, 60, and 80% (Supplementary Text 2 and Supplementary Table 1). We focus the analysis on NCP during the crop-growing season. To evaluate the influence from ENOx_agri, we estimate the population-weighted regional mean SIA reductions by only reducing ENH3_agri as well as reducing both ENH3_agri and ENOx_agri in the agricultural Nr experiments. In the combustion Nr experiments, we consider the combined reductions in NOx and NH3 emissions from combustion sources. These experiments are designed to represent the most simplified case, in which the NH3 and NOx are reduced synchronously for each source. The optimal pathway is estimated for 2019 and for the case that the anthropogenic emissions are changed to 2040 levels predicted by the SSP1-26 scenario. The influence from ENOx_agri is evaluated by estimating the optimal pathways with/without including ENOx_agri in agricultural Nr emission control.

We first compare the decline in SIAs due to reducing Nr from combustion sources and agricultural emissions with/without considering ENOx_agri. As shown by Fig. 4a, when ignoring ENOx_agri, the decline in population-weighted regional mean SIA concentrations over NCP due to reducing agricultural Nr emissions is smaller than reducing combustion Nr emissions. We define R as the ratio of SIAs decline due to controlling agricultural and combustion Nr emissions, which represents the relative importance of controlling agricultural Nr emissions in SIA mitigations. We estimate that without considering ENOx_agri, R gradually increased from around 0.5 to 0.8, with the emission reduction increasing from 20 to 80% for both the 2019 and 2040 cases. Such an increase in R reflects that with a deeper reduction in ENH3_agri, the SIA changes from NH3 saturated to NH3 limited situation (Supplementary Fig. 12), and the SIA formations will be more sensitive to NH3 emissions21. However, controlling reactive nitrogen emissions from combustion sources is more effective than agricultural emissions (R < 1), even with 80% emission reductions. Including ENOx_agri in agricultural Nr emissions reduction can obviously raise the relative importance of controlling agricultural emissions. Particularly, when anthropogenic emissions are reduced to the 2040 level, R for 20% and 80% Nr emission reduction becomes 0.8 and 1.0, making the agricultural Nr emission control able to achieve similar SIA mitigations as combustion Nr for a wide range (20–80%) of Nr emission control process.

Fig. 4: Influence of agricultural NOx emissions on the design of optimal reactive nitrogen emission control pathways over the North China Plain.
figure 4

a, d The bar plots show the decline in population-weighted SIA concentrations ((Delta {PWC})) due to 20–80% reductions in emissions from agricultural NH3 (gray), agricultural Nr (NOx and NH3) (yellow), and combustion Nr emissions (red) (left axis). The blue (green) lines show the ratio of SIAs decrease due to reductions of agricultural Nr (NH3) and combustion Nr emissions (right axis). b, c, e, f Diagrams illustrating the effectiveness of controlling Nr emissions from agricultural and combustion sources in mitigating population SIA concentrations. The background color indicates the population-weighted SIA concentrations over NCP, and the red lines show the optimal emission reduction pathways. Panels (b, e) are based on estimates that include only agricultural NH3 emission reductions, while panels (c, f) include both NH3 and NOx emissions in agricultural emission reductions. The arrows indicate the relative magnitude of SIA reductions contributed by reductions in agricultural and combustion reactive nitrogen emissions. The top panels (ac) use 2019 emissions, and the bottom panels (df) use emissions reduced to the 2040 level predicted by the SSP1-26 scenario.

Full size image

We examine the influence of ENOx_agri on designing the optimal Nr emission control pathway following Liu et al.34 (Supplementary Text 2) in Fig. 4. The optimal pathway is estimated following the direction of the gradients of the SIA concentrations mitigation. If ENOx_agri is not considered, under both 2019 and the 2040 anthropogenic emission conditions predicted by the SSP1-26 scenario, the optimal emission control pathway tends to reduce much more combustion Nr emissions than agricultural NH3 emissions. The gradients suggest that to achieve the most efficient SIA mitigation, a 20–80% reduction in combustion Nr emissions should be along with 10–26% and 8–26% agricultural NH3 emissions, respectively, under 2019 and 2040 anthropogenic emissions. Such reduction of Nr emissions will respectively achieve 3.6–13.7 μg m−3 and 1.6–5.3 μg m−3 population-weighted mean SIA mitigations. When including ENOx_agri reduction, we find that the optimal pathway significantly shifts toward controlling agricultural emissions. Particularly, for the 2040 anthropogenic emission case, we estimate that the agricultural Nr emissions should almost be reduced simultaneously with combustion Nr emission reductions to achieve the most efficient SIA reductions, where the 20–80% reduction in combustion Nr emissions is suggested along with 14–59% agricultural Nr reduction. The associated population-weighted mean SIA mitigations are increased to 2.0–6.3 μg m−3, which are 16–38% higher than the case that ENOx_agri is neglected.

Discussion

In summary, our study reveals the crucial role of controlling NOx emissions originating from agricultural activities in mitigating future PM2.5 pollution over eastern China during crop-growing season, particularly in regions characterized by intensive agricultural practices, such as NCP. Despite the substantially lower magnitude of agricultural NOx (0.3 Tg N yr-1) emissions compared to NH3(8.5 Tg N yr-1), we find that the efficiency of reducing agricultural NOx emissions in mitigating SIAs will become comparable to reducing NH3 emissions over the North China Plain when fossil fuel emissions decrease to 2060 levels predicted by future emission scenarios, reflecting the nonlinear chemistry mechanisms of SIA formations. As a result, when designing optimal emission control pathways, including agricultural NOx emissions will shift the emission control strategies more towards reducing agricultural sources, and will achieve more efficient SIA mitigations.

In this study, we consider a simplified scenario, in which agricultural NOx and NH3 are synchronously controlled. This assumption is reasonable for a quick evaluation of the implications of agricultural NOx emissions on emission control strategies, since most of the nitrogen management practices, such as reducing the overuse of chemical fertilizer, can collectively reduce NH3 and NOx emissions. By synthesizing a wealth of field measurements, previous studies have summarized the impact of 11 different reactive nitrogen management practices on nitrogen loss20. These practices vary in terms of technical thresholds, mitigation efficacy, and implementation costs. Most of the practices affect both NH₃ and NOₓ emissions, but the reduction in NH3 and NOx emissions in response to a single practice can differ significantly (Supplementary Table 2). For example, using enhanced-efficiency fertilizers can reduce NH₃ emissions by 69% and NOₓ emissions by 46%, while irrigation management (e.g., reduced irrigation or drip irrigation) can reduce NH₃ emissions by 36% and NOₓ emissions by 93%20. The increased contribution of agricultural NOx emissions does not only occur in Eastern China. As predicted by SSP scenarios, over the anthropogenic emissions hotspots, such as India, Southeast Asia, and Western Europe, the contribution of agricultural sources to total anthropogenic sources will reach ~20% by 2060, similar to the contribution over Eastern China (Supplementary Fig. 13). Thus, to meet more stringent air quality standards and achieve continuous air quality improvement over agricultural-intensive areas in Eastern China and in other regions, we call for a detailed consideration of how nitrogen management practices can influence agricultural NOx emissions for designing future agricultural emission control strategies.

Methods

GEOS-Chem model description

We use the global three-dimensional chemistry transport model GEOS-Chem (version 12.9.3) to simulate the PM2.5 concentrations and their response to emission changes over China during the crop-growing season in 2019. We use the nested version of GEOS-Chem with a horizontal resolution of 0.625°(longitude) × 0.5°(latitude) over East Asia and 2.5°(longitude) × 2°(latitude) over the rest of the world. The model incorporates a comprehensive state-of-the-art O3-NOx-hydrocarbon-aerosol-halogen chemistry in the troposphere35,36. To estimate the partition of gas phase H2SO4-HNO3-NH3 with SIAs, the model employed the ISORROPIA-II thermodynamical model37. The chemistry, transportation, and deposition of chemical species as well as the biogenic emissions including soil NOx are driven by MERRA-2 meteorological data provided by the Global Modeling and Assimilation Office (GMAO) at the National Aeronautics and Space Administration (NASA). The MERRA-2 data has a temporal resolution of 1 h for surface variables and boundary layer height and 3 h for other three-dimensional variables.

The anthropogenic emissions other than agricultural NH3 over China are from the Multi-resolution Emission Inventory for China (MEIC; http://meicmodel.org/)7,38. The agricultural NH3 emissions are from the bottom-up emission inventory given by Zhang et al.7 and Chen et al.23. The natural emissions include NOx emissions from soil estimated by the Berkeley–Dalhousie Soil NOx emission Parameterization (BDSNP)22 and from lightning parameterized based on the cloud-top height39, biogenic volatile organic compound (VOC) emissions estimated by the Model of Emissions of Gases and Aerosols from Nature (MEGAN v2.1) algorithm40, and biomass burning given by the Global Fire Emissions Database version 4 (GFED4)41. By employing the aforementioned emission inventories, the model can generally capture the spatial distribution of PM2.5 concentration over China during the crop-growing season (May-August) when compared with the surface observations collected from the Ministry of Ecology and Environment of China (Supplementary Fig. 4).

Agricultural NOx emission parameterization

The agricultural NOx emissions are estimated as a part of soil NOx emissions by the BDSNP that was implemented in GEOS-Chem22. This parameterization characterizes the soil NOx emissions by considering the available soil nitrogen as well as meteorological conditions, including temperature and water-filled pore space (defined as the ratio of volumetric soil water content to total soil porosity). The dependence of agricultural NOx emissions on soil organic matter42 is not included in BDSNP parameterization. The BDSNP also considers the pulse of soil NOx emissions when the water-stressed bacteria is reactivated by irrigation or rainfalls and the canopy uptake of NOx. Lu et al.12 have demonstrated that the BDSNP parameterization can generally capture the spatial distribution of soil NOx emission in China when compared with available field measurements and satellite observation of tropospheric NO2 column density.

The available soil nitrogen includes the natural nitrogen pool as well as the external nitrogen input from fertilizer application and atmospheric deposition. Thus, agricultural NOx emissions can be calculated using the BDSNP parameterization by distinguishing the soil NOx emissions resulting from fertilizer nitrogen input. In Hudman et al.22, the spatial distribution of chemical fertilizer and manure is from Potter et al.43, with 37% of the manure added to soil nitrogen content. Since Potter et al.43 estimated the chemical fertilizer application and manure based on data reported from 1994 to 2001, we updated the dataset based on Chen et al.23, which estimated the monthly chemical fertilizer and manure for China between 2005 and 2016. Unlike the assumption made by Hudman et al.22 that 75% of the fertilizer is applied around the green-up day, Chen et al.23 considered detailed fertilizer application practices in China. In addition, the manure applied to the soil is estimated by the process-based mass-flow method. Annually, the total fertilizer application (chemical fertilizer and manure applied to soil) from Chen et al.23 (35 Tg yr-1) is higher than that given by Potter et al.43 (Supplementary Fig. 2). Besides direct fertilizer application, agricultural activities also indirectly contribute to soil nitrogen content through enhanced atmospheric deposition of reactive nitrogen. Through a sensitivity simulation that turns off all agricultural reactive nitrogen emissions, we estimate that soil NOx emissions driven by agricultural-induced nitrogen deposition account for only 3% of those from direct fertilizer application. Using the BDSNP parameterization, we calculate total soil NOx emissions in China to be 0.68 Tg N yr⁻¹, of which 0.33 Tg N is induced by agricultural fertilizer application, identified as agricultural NOx emissions in this study. The agricultural NOx emissions in China calculated with the updated fertilizer input data are 43% higher than the estimate obtained using the original BDSNP parameterization (Supplementary Fig. 14). The largest differences are observed over the North China Plain, where our updated estimates can exceed the original BDSNP by as much as 4 kg N ha-1 yr-1.

Future emission scenarios

We investigate the contribution of agricultural NOx emissions on PM2.5 pollution in the future based on the anthropogenic emission changes given by the sixth Coupled Model Intercomparison Project (CMIP6) scenarios44. The CMIP6 scenarios encompass projections of climate forcing in 2100 and future social evolutions by linking the Representative Concentration Pathways (RCPs) with the Shared Socioeconomic Pathways (SSP1-5). The Tier 1 scenarios of CMIP6 provide a comprehensive range of future climate-forcing targets, which include two baseline scenarios (SSP3-70 and SSP5-85) and two mitigation scenarios (SSP1-26 and SSP2-45). In this study, we chose the emission mitigation scenarios taking into account emission control policies in China. The SSP2-45 and SSP1-26 present a middle road and a strict mitigation pathway scenario, respectively. We only apply the annual scale factor for each emission sector estimated for the integrated assessment model (IAM) regions44 to GEOS-Chem emissions to represent future anthropogenic emission changes under various climate targets and emission control strategies.

Model experiments

We investigate the role of agricultural NOx emissions on PM2.5 pollution over eastern China through a series of GEOS-Chem model experiments as listed in Supplementary Table 1. The base simulations use the anthropogenic emission inventories described above, and future PM2.5 changes use the anthropogenic emissions for 2040 and 2060 predicted by SSP2-45 and SSP1-26 scenarios. To assess the contribution of agricultural and combustion NOx emissions to PM2.5 pollution, we conduct sensitivity simulations Group I and II, which are the same as base simulations but turn off agricultural NOx emissions and combustion NOx emissions, respectively. We further examine the role of agricultural NOx emissions in controlling present and future PM2.5 pollution by conducting simulations sensitivity simulations in Groups III and IV. Compared with the standard simulations, simulations in Group III reduce agricultural NH3 emissions by 40%, and simulations in Group IV reduce both agricultural NH3 and NOx emissions by 40%. To evaluate how the uncertainties in ENOx_agri will influence the efficiency of controlling agricultural NOx on mitigating SIA concentrations in the future, we also conduct simulations that assume a factor of 2 uncertainty (50–200%) in ENOx_agri following Lu et al.12 for sensitivity simulations Group III and IV for 2060. To estimate the optimal reactive Nr emission control pathways, we conduct sensitivity simulations in Groups V, VI, VII, which gradually reduce the Nr emissions from agricultural/combustion sources by each 20% step. We calculate the population-weighted regional mean SIA concentrations (PWC) based on the gridded SIA concentrations (({C}_{{{mathrm{SIA}}},i,j})) and the corresponding population densities (({P}_{i,j})) from LandScan Global 201945 .

$${{mathrm{PWC}}}=frac{sum _{i,j}{C}_{{{mathrm{SIA}}},i,j}times {P}_{i,j}}{sum _{i,j}{P}_{i,j}}$$

We should note here since we only focus on the changes in anthropogenic emissions, all above simulations are conducted with the meteorological conditions fixed to 2019. The uncertainties from emissions and meteorological conditions are discussed in the Supplementary Materials.

Related Articles

Nitrogen dioxide (NO2) Meteorology and predictability for air quality management using TROPOMI

Nitrogen dioxide (NO₂) is a critical air pollutant and key indicator for air quality. Due to limited monitoring, we leveraged TROPOMI NO₂ and NASA POWER meteorological datasets to evaluate the meteorological drivers on NO₂ tropospheric column concentrations and to develop predictive models for NO₂ levels over Ghana. Employing an 8:2 ratio for model training and testing, NO₂ and meteorology relationships were assessed by seasonality indices and correlation analyses. Results indicate marked seasonal variability in NO₂ columns, prominent during the dry season. Wind speed, relative humidity, and precipitation significantly reduce NO₂, whereas temperature correlated positively in the southern forested zone. Predictive models demonstrate varying efficacy across climatic zones, with mean percentage differences ranging 9.87 to 37.76% and agreement index up to 0.96. The Random Forest and XGBoost models showed outstanding performance, with correlation reaching 0.92. This results presents a scalable methodology for NO₂ monitoring providing insights for air quality management.

Tomato growth stage modulates bacterial communities across different soil aggregate sizes and disease levels

Soil aggregates contain distinct physio-chemical properties across different size classes. These differences in micro-habitats support varied microbial communities and modulate the effect of plant on microbiome, which affect soil functions such as disease suppression. However, little is known about how the residents of different soil aggregate size classes are impacted by plants throughout their growth stages. Here, we examined how tomato plants impact soil aggregation and bacterial communities within different soil aggregate size classes. Moreover, we investigated whether aggregate size impacts the distribution of soil pathogen and their potential inhibitors. We collected samples from different tomato growth stages: before-planting, seedling, flowering, and fruiting stage. We measured bacterial density, community composition, and pathogen abundance using qPCR and 16 S rRNA gene sequencing. We found the development of tomato growth stages negatively impacted root-adhering soil aggregation, with a gradual decrease of large macro-aggregates (1–2 mm) and an increase of micro-aggregates (<0.25 mm). Additionally, changes in bacterial density and community composition varied across soil aggregate size classes. Furthermore, the pathogen exhibited a preference to micro-aggregates, while macro-aggregates hold a higher abundance of potential pathogen-inhibiting taxa and predicted antibiotic-associated genes. Our results indicate that the impacts of tomatoes on soil differ for different soil aggregate size classes throughout different plant growth stages, and plant pathogens and their potential inhibitors have different habitats within soil aggregate size classes. These findings highlight the importance of fine-scale heterogeneity of soil aggregate size classes in research on microbial ecology and agricultural sustainability, further research focuses on soil aggregates level could help identify candidate tax involved in suppressing pathogens in the virtual micro-habitats.

Bank lending and environmental quality in Gulf Cooperation Council countries

To achieve economies with net-zero carbon emissions, it is essential to develop a robust green financial intermediary channel. This study seeks empirical evidence on how domestic bank lending to sovereign and private sectors in Gulf Cooperation Council (GCC) countries impacts carbon dioxide and greenhouse gas emissions. We employ PMG-ARDL model to panel data comprising six countries in GCC over twenty years for carbon dioxide emissions and nineteen years for greenhouse gas emissions. Our findings reveal a long-term positive impact of both bank lending variables on carbon dioxide and greenhouse gas emissions. In addition, lending to the government shows a negative short-term effect on greenhouse gas emissions. The cross-country results demonstrate the presence of a long-run effect of explanatory variables on both types of emissions, except for greenhouse gas in Saudi Arabia. The sort-term impact of the explanatory variables on carbon dioxide and greenhouse gas emissions is quite diverse. Not only do these effects differ across countries, but some variables have opposing effects on the two types of emissions within a single country. The findings of this study present a new perspective for GCC economies: neglecting total greenhouse gas emissions and concentrating solely on carbon dioxide emissions means missing critical information for devising effective strategies to combat threats of environmental degradation and achieve net-zero goals.

Higher temperature sensitivity of forest soil methane oxidation in colder climates

Forest soils, serving as an important sink for atmospheric methane (CH4), modulate the global CH4 budget. However, the direction and magnitude of the forest soil CH4 sink under warming remain uncertain, partly because the temperature response of microbial CH4 oxidation varies substantially across geographical scales. Here, we reveal the spatial variation in the response of forest soil microbial CH4 oxidation to warming, along with the driving factors, across 84 sites spanning a broad latitudinal gradient in eastern China. Our results show that the temperature sensitivity of soil microbial CH4 oxidation significantly declines with increasing site mean annual temperature, with a range of 0.03 to 0.77 μg CH4 g–1 soil d–1 °C–1. Moreover, soil resources and type II methanotrophs play crucial roles in shaping the temperature sensitivity of soil microbial CH4 oxidation. Our findings highlight the importance of incorporating climate, soil resources, and methanotroph groups into biogeochemical models to more realistically predict forest soil CH4 sink under warming.

Crop rotation increases Tibetan barley yield and soil quality on the Tibetan Plateau

Tibetan barley (Hordeum vulgare) accounts for over 70% of the total food production in the Tibetan Plateau. However, continuous cropping of Tibetan barley causes soil degradation, reduces soil quality and causes yield decline. Here we explore the benefits of crop rotation with wheat and rape to improve crop yield and soil quality. We conducted 39 field experiments on the Tibetan Plateau, comparing short-term (≤5 years), 5–10 years and long-term (≥10 years) continuous cropping with rotation of Tibetan barley with wheat or rape. Results showed that Tibetan barley–wheat and Tibetan barley–rape rotations increased yields by 17% and 12%, respectively, while improving the soil quality index by 11% and 21%, compared with long-term continuous cropping. Both Tibetan barley rotations with wheat and rape improved soil quality and consequently yield, mainly by increasing soil microbial biomass nitrogen and microbial biomass carbon and decreasing pH. By contrast, long-term continuous cropping led to decreased soil organic matter, lower microbial biomass nitrogen and increased pH, contributing to yield decline. The benefits of rotations on crop yield and soil quality increased over time. Implementing crop rotation with wheat or rape thus offers a sustainable agricultural strategy for improving food security on the Tibetan Plateau.

Responses

Your email address will not be published. Required fields are marked *