Model Inter-Comparison Study for Asia (MICS-Asia) phase III: multimodel comparison of reactive nitrogen deposition over China

Atmospheric nitrogen deposition in China has attracted public attention in recent years due to the increasing anthropogenic emission of reactive nitrogen (Nr) and its impacts on the terrestrial and aquatic ecosystems. However, limited long-term and multisite measurements have restrained the understanding of the mechanism of the Nr deposition and the chemical transport model (CTM) improvement. In this study, the performance of the simulated wet and dry deposition for different Nr species, i.e., particulate NO−3 and NH + 4 , gaseous NOx, HNO3 and NH3 have been conducted using the framework of Model Inter-Comparison Study for Asia (MICS-Asia) phase III. A total of nine models, including five Weather Research and Forecasting models coupled with the Community Multiscale Air Quality (WRF-CMAQ) models, two self-developed regional models, a global model and a Regional Atmospheric Modeling System coupled with the Community Multiscale Air Quality (RAMS-CMAQ) model have been selected for the comparison. For wet deposition, observation data from 83 measurement sites from the East Asia Acid Deposition Monitoring Network (EANET), Chinese Ecosystem Research Network (CERN), China Agricultural University Deposition Network (CAUDN), National Acid Deposition Monitoring Network (NADMN) and Department of Ecological Environment (DEE) of China have been collected and normalized for comparison with model results. In general, most models show the consistent spatial and temporal variation of both oxidized N (Nox) and reduced N (Nrd) wet deposition in China, with the normalized mean error (NME) at around 50 %, which is lower than the value of 70 % based on EANET observation over Asia. Both the ratio of wet or dry deposition to the total inorganic N (TIN) deposition and the ratios of TIN to their emissions have shown consistent results with the Nationwide Nitrogen Deposition Monitoring Network (NNDMN) estimates. The performance of ensemble results (ENMs) was further assessed with satellite measurements. In different regions of China, the results show that the simulated Nox wet deposition was overestimated in northeastern China (NE) but underestimated in the south of China, namely southeastern (SE) and southwestern (SW) China, while the Nrd wet deposition was underestimated in all regions by all models. The deposition of Nox has larger uncertainties than the Nrd, especially in northern China (NC), indicating the chemical reaction process is one of the most important factors affecting the model performance. Compared to the critical load (CL) value, the Nr deposition in NC, SE and SW reached or exceeded reported CL values and resulted in serious ecological impacts. The control of Nrd in NC and SW and Nox in SE would be an effective mitigation measure for TIN deposition in these regions. The Nr deposition in the Tibetan Plateau (TP) with a high ratio of TIN / emission (∼ 3.0), indicates a significant transmission from outside. Efforts to reduce these transmissions ought to be paramount due the climatic importance of the Tibetan region to the sensitive ecosystems throughout China.

B. Ge et al.: MICS-Asia III Abstract. Atmospheric nitrogen deposition in China has attracted public attention in recent years due to the increasing anthropogenic emission of reactive nitrogen (N r ) and its impacts on the terrestrial and aquatic ecosystems. However, limited long-term and multisite measurements have restrained the understanding of the mechanism of the N r deposition and the chemical transport model (CTM) improvement. In this study, the performance of the simulated wet and dry deposition for different N r species, i.e., particulate NO − 3 and NH + 4 , gaseous NO x , HNO 3 and NH 3 have been conducted using the framework of Model Inter-Comparison Study for Asia (MICS-Asia) phase III. A total of nine models, including five Weather Research and Forecasting models coupled with the Community Multiscale Air Quality (WRF-CMAQ) models, two self-developed regional models, a global model and a Regional Atmospheric Modeling System coupled with the Community Multiscale Air Quality (RAMS-CMAQ) model have been selected for the comparison. For wet deposition, observation data from 83 measurement sites from the East Asia Acid Deposition Monitoring Network (EANET), Chinese Ecosystem Research Network (CERN), China Agricultural University Deposition Network (CAUDN), National Acid Deposition Monitoring Network (NADMN) and Department of Ecological Environment (DEE) of China have been collected and normalized for comparison with model results. In general, most models show the consistent spatial and temporal variation of both oxidized N (N ox ) and reduced N (N rd ) wet deposition in China, with the normalized mean error (NME) at around 50 %, which is lower than the value of 70 % based on EANET observation over Asia. Both the ratio of wet or dry deposition to the total inorganic N (TIN) deposition and the ratios of TIN to their emissions have shown consistent results with the Nationwide Nitrogen Deposition Monitoring Network (NNDMN) estimates. The performance of ensemble results (ENMs) was further assessed with satellite measurements. In different regions of China, the results show that the simulated N ox wet deposition was overestimated in northeastern China (NE) but underestimated in the south of China, namely southeastern (SE) and southwestern (SW) China, while the N rd wet deposition was underestimated in all regions by all models. The deposition of N ox has larger uncertainties than the N rd , especially in northern China (NC), indicating the chemical reaction process is one of the most important factors affecting the model performance. Compared to the critical load (CL) value, the N r deposition in NC, SE and SW reached or exceeded reported CL values and resulted in serious ecological impacts. The control of N rd in NC and SW and N ox in SE would be an effective mitigation measure for TIN deposition in these regions. The N r deposition in the Tibetan Plateau (TP) with a high ratio of TIN / emission (∼ 3.0), indicates a significant transmission from outside. Efforts to reduce these transmissions ought to be paramount due the climatic importance of the Tibetan region to the sensitive ecosystems throughout China.

Introduction
Atmospheric nitrogen (N) deposition is defined as N-related gases, and particles are deposited with precipitation occurring (wet deposition) and without precipitation occurring (dry deposition; Clark and Kremer, 2005). These deposits on the Earth's surface are either close to the sources or in remote regions (e.g., chemical transformation and long-range transport of oxidized and reduced N; hereafter N ox and N rd ) located far from human activities and labeled as the N-limited areas (Phoenix et al., 2006;Holtgrieve et al., 2011). Evidence shows that the effects of reactive N (N r = N ox + N rd ) deposition on the environment are numerous, including decreased biological diversity, increased soil acidification and lake eutrophication (Clark and Tilman, 2008;Janssens et al., 2010;Holtgrieve et al., 2011;Phoenix et al., 2006;Galloway et al., 2004). Different human activities disturb the natural N cycle in serious ways (Galloway et al., 2004); for example, using artificial fertilizers to increase crop production  or excessively relying on fossil fuels for industrial production. N r production increased from approximately 15 Tg N yr −1 in 1860 to 187 Tg N yr −1 in 2005, and more than 50 % of that N r has reportedly been deposited onto the ground (Nicolas and Galloway, 2008). In the past two decades, high rates of N r deposition were widely documented in developed countries, such as the USA (Fenn et al., 1998) and Europe (Dise and Wright, 1995). Great efforts have been made to fight against these negative effects in the USA, and the N ox deposition was decreased dramatically in recent years (Li et al., 2016). However, the growing human demand for food and energy on a global scale has resulted in increased emissions of N r into the environment , particularly in large developing countries like China and India (Chen et al., 2019a;Liu et al., 2013).
A nationwide estimate of long-term N deposition in China, based on the bulk measurements and summaries from reported references in 270 sites by Liu et al. (2013), showed an increasing rate of 0.41 kg N ha −1 per year from 1980 to 2010. In contrast to the increasing importance of N rd deposition, due to apparently decreasing N ox resulting from the air quality control policies in the USA in the past decades (Li et al., 2016), the ratio of N rd / N ox recorded from bulk and/or wet deposition decreased from 5 : 1 in 1980 to 2 : 1 in 2010. This suggests a more and more important role of N ox in China (Liu et al., 2013). The ratio in highly developed regions, such as the North China Plain, was even lower than 1 : 1 in recent years . However, very limited long-term observations in China challenge our capacity to understand and control the increase in N r deposition. The published long-term N deposition monitoring network, which includes the East Asia Acid Deposition Monitoring Network (EANET; http://www.eanet.asia, last access: 6 September 2020), the Nationwide Nitrogen Deposition Monitoring Network (NNDMN), established in 2010 by the China Agriculture University (CAU; Xu et al., 2015), the Chinese Ecosystem Research Network (CERN) in the North China Plain, established by the Chinese Academy of Science , the Acid Rain Monitoring Network, run by the China Meteorological Administration (CMA-ARMN; Tang et al., 2007Tang et al., , 2010Ge et al., 2011), and the National Acid Deposition Monitoring Network (NADMN; Li et al., 2019b) have been identified with many shortcomings. Monitoring sites are widely spread over a large geographical area, and therefore data records, due to the high cost of the measurement and unstable financial support, are incomplete. Chemical transport model (CTM) simulation is another option to offset these drawbacks and also to quantify the longrange transport of deposition in a global or regional map (Seinfeld and Pandis, 2006). It is important to know the accuracy of the CTM before it is employed for investigating the spatial and temporal variation of the depositions.  and Mann et al. (2014) explain that different parameters in CTMs can result in large variations and the inaccuracies. The multimodel ensemble mean (ENM) shows better performance than any single model (Carmichael et al., 2002;Hayami et al., 2008;Holloway et al., 2008;. Additionally, to better localize the applications of CTM, comprehensive evaluations of the strengths and weaknesses of current CTMs for simulating the acid deposition and their precursors in a unified framework, with certain regulated rules and the same inputs to models, must be undertaken. The Model Inter-Comparison Study for Asia (MICS-Asia) provides an opportunity to investigate the CTM application with different models in Asia. MICS-Asia was first employed in 1998, with the target of long-range transport and deposition in SO 2− 4 in the first stage (MICS-Asia phase I; Carmichael et al., 2002) and sulfur, nitrogen and ozone in the second stage (MICS-Asia phase II; Carmichael et al., 2008). The findings and methodologies developed in the previous intercomparison studies contributed to common understandings of the performance and uncertainties of CTM applications in East Asia Carmichael et al., 2008;Han et al., 2008;Wang et al., 2008). The comprehensive multimodel intercomparison study on acid deposition in China is becoming an urgent issue as the high emissions in China are causing acid deposition in neighboring countries (Lin et al., 2008;Kajino et al., 2011Kajino et al., , 2013Itahashi et al., 2018). In this study, 1 year of simulated N r depositions, i.e., N ox and N rd in both wet and dry depositions, using the framework of MICS-Asia phase III, have been compared with each other and validated by the observed wet deposition from EANET, NNDMN, CERN and by the Department of Ecological Environment (DEE; formerly known as the Environmental Protection Administration -EPA) over the whole of China. The ENM results were also compared to the vertical column density (VCD) from satellite and emission inventories. Finally, the uncertainties of the sources of N r depositions and their ecological impacts have been quantified. The results from this study will not only provide an impor-tant reference for establishing a suitable N deposition model, but the localized application of CTMs in China will also be tested.

Description of the participant models
In phase III of MICS-Asia, 14 chemical transport models (CTM; M1-M14) were used to compare and evaluate current multiscale air quality models (called topic 1 in MICS-Asia III). The same number index was used to measure aerosols and ozone levels, which was reported by Chen et al. (2019b) and Li et al. (2019a). However, the fully coupled online Weather Research and Forecasting model with chemistry (WRF-Chem), which has been indexed as M7-M10, was not included in the deposition comparison part in the overview of model intercomparison and evaluation for acid deposition in Asia . Briefly, the Weather Research and Forecasting model coupled with the Community Multiscale Air Quality (WRF-CMAQ) has been numbered M1-M6, with a different version of version 5.0.2 for M1 and M2, version 5.0.1 for M3 and version 4.7.1 for M4-M6. M11 and M12 are the independent models developed in Japan and China, namely a nonhydrostatic model coupled with chemistry (NHM-Chem; Kajino et al., 2019) and the Nested Air Quality Prediction Model System (NAQPMS), respectively. A global 3D chemical transport model (namely the Goddard Earth Observing System global chemical transport (GEOS-Chem) version 9.1.3), M13, was also used as the long-range transport and future change prediction model in MICS-Asia III. The last, M14, was represented as the Regional Atmospheric Modeling System coupled with CMAQ (RAMS-CMAQ). It should be noted that the last two models, M13 and M14, were not driven by the "standard" meteorological fields from the WRF version 3.4.1 model. Basic information about the configuration of each model is summarized in Table 1. More detailed descriptions can also be found in previous studies Chen et al., 2019b;Li et al., 2019a).

Model inputs and simulation domain
As mentioned by Chen et al. (2019b), the same ("standard") meteorological fields, emission inventories and boundary conditions have been prepared for the CTM intercomparison in MICS-Asia III to reduce the uncertainties from model inputs. However, some models, such as M13 and M14, imported "nonstandard" inputs due to their specific characteristics. The "standard" meteorological inputs were simulated by WRF version 3.4.1, with the initial and lateral boundary conditions from the National Centers for Environmental Prediction (NCEP) Final Analysis (FNL) data. A 4D data assimilation (FDDA) nudging was adopted every 6 h to improve the accuracy of the meteorological parameters simu-  Pielke et al. (1992). For the emission inputs, all the participant models were using the same emission inventory, which included the MIX anthropogenic emissions over Asia developed for MICS-Asia phase III (Li et al., 2017), the biogenic emissions calculated by the Model of Emissions of Gases and Aerosols from Nature (MEGAN) version 2.04 (Guenther et al., 2006), and the biomass burning emissions from Global Fire Emission Database ( MICS-Asia phase III provided two sets of lateral boundary conditions derived from GEOS-Chem  and the CHASER chemical transport model (Sudo et al., 2002), respectively. The boundary conditions from GEOS-Chem were run with 2.5 • × 2 • resolution and 47 vertical layers, while those from CHASER were run with 2.8 • × 2.8 • and 32 vertical layers. M4, M5, M6, M11 and M12 used the output from CHASER as the boundary conditions, and M1, M13 and M14 were from GEOS-Chem. Only M2 used the default boundary condition field provided in CMAQ.
The standard simulation domain covers the region of East Asia (15.4 • S-58.3 • N, 48.5-160.2 • E) with 180 × 170 grids at 45 km horizontal resolution. M1-M6, M11 and M12 followed standard simulation domains, while M13 and M14 employed different modeling domains with 0.5 • latitude × 0.667 • longitude and 64 × 64 km, respectively. In this study, the analyzed region was only focused in China, and all participant models covered it. Therefore, simulated reactive N deposition in each model can be analyzed and compared to show the performance of the participant models. All model outputs of N depositions have been classified as oxidized N (N ox = gHNO 3 + gNO x + pNO − 3 , including gaseous nitrate acid, NO x and particulate nitrate) and reduced N (N rd = gNH 3 + pNH + 4 , including gaseous ammonia and particulate ammonium) for comparison. Several previous studies have reported that the contributions of the other N species, e.g., PAN and isoprene nitrates to total N dry deposition, are less important than the inorganic N, e.g., HNO 3 and NO x (Yuan et al., 2018;Wolfe et al., 2011). Thus, these organic N species were not included in this study.

Observation data
China has a large land mass -almost 5500 km from south to north (3.5-53.3 • N) and 5200 km from west to east (75.5-135 • E). It goes from coastal to inland and from a tropical climate to a frigid zone. Only eight sites located in Guangdong, Fujian, Sichuan and Shanxi in EANET were insufficient to show the real performance of CTMs in China. Besides the eight EANET sites, data from 83 sites recorded daily, weekly or yearly from the CERN , NNDMN (Xu et al., 2015;Liu et al., 2013) and DEE in Guangdong, Liaoning and Xinjiang province and Shanghai were employed in this study to compare the simulated wet deposition in MICS-Asia III in China. Figure 1 shows the location of the 83 measurement sites and the divided regions of the whole of China. There were 50 urban sites and 33 rural sites. More detailed information on each measurement site can also be found in Table S1 in the Supplement.
The daily wet deposition was measured by a wet-only sampler to collect precipitation samples during the rainfall events in EANET. Analysis methods for NO − 3 and NH + 4 were based on ion chromatography and checked by the ion balance and conductivity agreement. Detailed descriptions can be found in the source document (EANET, 2010). Daily rainwater samples at 10 sites located in the North China Plain were collected using a custom wet-dry automatic collector (APS-2B; Xianglan Scientific Instruments Co., Ltd., Changsha, China) in CERN. Inorganic N, including NO − 3 and NH + 4 , in the precipitation samples was measured using an ion chromatography system (model ICS-90; Dionex Corporation, Sunnyvale, CA, USA) and the standard laboratory procedure of the State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry (LAPC; Wang et al., 2012). The detection limit (DL) of N for this instrument was 5 µg L −1 . A detailed description can be found in the research of Pan et al. (2012). The wet and/or bulk NO − 3 and NH + 4 deposition data from 25 sites of China Agricultural University Deposition Network (CAUDN), which was renamed as NNDMN in China in 2010, has been collected and reanalyzed as yearly data (Xu et al., 2015;Liu et al., 2013). At all monitoring sites, precipitation samples were collected using precipitation gauges (SDM6; Tianjin Weather Equipment Inc., China) located beside the DEnuder for Long-Term Atmospheric sampling (DELTA) systems (ca. 2m). After their collection, the samples were analyzed in CAU's laboratory, based on the standard laboratory procedure of CAU (Xu et al., 2015). Routine NO − 3 and NH + 4 wet depositions collected in each rainfall event at 40 sites were provided by the DEEs of Guangdong, Liaoning and Xinjiang provinces and Shanghai city. The analytical process was the same used in the laboratory procedure of the China National Environmental Monitoring Center (CNEMC).
The temporal resolutions of the wet deposition data are different from each other, i.e., daily in EANET and CERN, rainfall event collection in DEE and yearly in NNDMN. For the convenience of comparison, all data from daily or rainfall event collecting samples at each type of measurement site have been summarized and interpolated as monthly wet deposition data to compare with the monthly simulation in MICS-Asia III in this study, except for the yearly data provided by NNDMN. The VCD of NO 2 from SCIAMACHY (http://www.temis.nl/airpollution, last access: 6 Septem-ber 2020) and NH 3 from the Infrared Atmospheric Sounding Interferometer (IASI; https://iasi.aeris-data.fr/nh3/, last access: 6 September 2020) have also been used for comparison with the total deposition and the emissions.

Validation of wet deposition
3.1.1 Yearly comparison and monthly variation of oxidized N Yearly simulated wet deposition of N ox has been evaluated by observed nitrate wet deposition in 83 sites over China. Table 2 lists the statistical parameters of the simulated wet deposition of N ox compared with the observed data in rural and urban sites in China. For all sites, M1, M5 and M11 overestimated the yearly wet deposition of N ox with a normalized mean bias (NMB) of +30.3 %, +55.4 % and +67.2 %, respectively. M6, M12 and M13 simulated almost comparable results with a NMB of −6.8 %, +0.6 % and +0.1 %, respectively. M2, M4 and M14 underestimated the yearly wet deposition of N ox , with a NMB of −38.7 %, −10.7 % and −47.4 %, respectively. The NME was around 50 % with the highest being 82.2 % in M11, which is lower than that reported over East Asia with the value of 70 %, based on EANET observations by Itahashi et al. (2020). However, the correlation coefficient R was around 0.2 − 0.3 (n = 83), which is lower than the East Asia comparison based on the EANET data (0.3 − 0.4, n = 54; Itahashi et al., 2020). In order to eliminate influences from rainfall uncertainties (R = 0.82), the volume weighted mean (VWM) concentration of N ox in precipitation has also been evaluated. In contrast to the low R value of the yearly wet deposition of N ox , the correlation R increased to almost 0.5 for the VWM concentrations.
To judge the agreement between simulation and observation, the percentages within a factor of 2 (FAC2) have been calculated in this study. Approximately 50 % of model results corresponded within the percentages within FAC2. M1 and M13 showed better agreement with 60 % and 59 % within FAC2, while M2 and M14 showed only 36 % and 45 % agreement within FAC2. The ground surface measurement sites were divided into 49 urban sites and 34 rural sites, according to their location. Overall, all the models showed better performance with the R in 0.2−0.4 and FAC2 in 50 %−60 % in urban sites than that with R in 0.05−0.3 and FAC2 in 40 %−50 % in rural sites. This difference may not be due to the uncertainties in rainfall simulation, as the simulated VWM concentration of N ox in precipitation may eliminate the rainfall uncertainties and also shows better agreement in urban locations than that in rural sites (Table 2). Figure 2 shows the percentile boxplot of the yearly wet deposition of N ox simulated by nine participant models in five regions of China (i.e., northern China -NC; northeastern    Fig. S1 in the Supplement. The model results show large intraregion or interregion uncertainties, especially in NC, NE and SE. The highest wet deposition of N ox simulated by M11 is almost 3-4 times the lowest result simulated by M14 in the above regions (Fig. 2). Specifically, two models simulate 30 %-50 % higher N ox wet deposition, while four models are 20 %-40 % lower compared to the averaged observations in NC, with the averaged value of 6.5 kg N ha −1 a −1 . For the wet deposition of N ox in SE and SW and TP, most of the participant models are more than 50 % underestimated, with the largest underestimation being 75 % from M14, even though the precipitation in this region is overestimated. Additionally, the divergence of observed N ox wet deposition between different sites in NC, SE and SW, which was shown as the length of the red box in Fig. 2a, d and e, is significantly larger than the multimodels' results. The scattered distribution of the measurement sites in these regions is responsible for the large divergence in ob-servations. However, most of the participant models failed to capture the large difference, indicating that the coarse grid in MICS-Asia III (45 km) is not suitable for the performance of detailed characterization at a local scale. A global assessment of the ensemble simulated wet depositions in the Task Force on Hemispheric Transport of Atmospheric Pollutants (TF HTAP) by Vet et al. (2014) also indicated the underestimation of the models in a number of sites in North America, Europe, Central Africa and part of East Asia. The underprediction in Europe was found to be due to the large underpredictions of precipitation depth, while the reason for the error in East Asia is still unknown. However, most of the models overestimate the wet deposition of N ox in NE. Several models, including M1, M5 and M11, simulate more than 10 kg N ha −1 a −1 N ox wet deposition, almost double the observed value of 5 kg N ha −1 a −1 . Both the multimodels and the observations show very low values of 3-4 kg N ha −1 a −1 N ox wet deposition in NW, where the precipitation depth was very low compared to the other regions of China ( Fig. S1 in the Supplement).
Regarding to the comparison over the whole of East Asia reported in the overview of acid deposition in MICS-Asia III , similar overestimations were found in M5 and M11, while there was underestimation in M2, M4 and M14. It should be noted that the EANET sites are mostly located around Japan, Korean Peninsula and southeastern Asia, and only eight sites are located in China. The similar performances between the validations in East Asia and China indicated that the general underestimation (overestimation) of M2, M4 and M14 (M5 and M11) were reliable in these models. For the rest of models, different results were found between China and East Asia, i.e., the simulated N ox wet deposition in M1 was significantly overestimated in China ( Fig. 6 in Itahashi et al., 2020) but comparable with the observations over the rest of East Asia. Due to the absence of the observations for atmospheric NO 2 / NO − 3 , we cannot validate their model performances directly. Instead, another companion paper (Chen et al., 2019b) reported that most of the models overestimated NO − 3 concentrations based on 14 sites in China, with most sites located in NC (Fig. S5 in Chen et al., 2019b). In summary, the relationship between the atmospheric concentration of NO − 3 and the wet deposition in NC was not obvious, which is also same outcome as was found in East Asia .
Further evaluations of the temporal variations in both the urban and rural sites of NC and NE are displayed in Fig. 3. Generally, all of the models and observations found high levels of depositions in spring and summer and low values in winter in the two regions. High depositions were due to large precipitation depths in the rainy season. However, this was not always true in some urban sites. For example, higher depositions of N ox were observed in May and June, with lower rainfall volume, than in July and August, with higher rainfall, in the urban sites of NC. Similar examples were found at urban sites in NE. Previous studies confirmed that there is a decreasing trend in the variations of chemical components in precipitation as rainfall varies (Aikawa and Hiraki, 2009;Aikawa et al., 2014;Xu et al., 2017). If the rainfall lasted long enough, or if the rainfall volume was large enough, the concentrations of chemical components in precipitation remained at low levels and were attributed to the effects of the in-cloud scavenging process. That is, the large rainfall in an event may not cause the high level of monthly wet depositions due to the low level of in-cloud deposition as compared to the wet depositions by several different precipitation events, especially in polluted regions in urban sites. Unfortunately, only the monthly data of wet depositions and precipitation have been compared in this MICS-Asia III. Detailed comparisons of rainfall events should be considered in the future.

Yearly comparison and monthly variation of reduced N
Simulated wet deposition of N rd in MICS-Asia III has been evaluated using multisource observations from the same sites referred to in N ox . It is shown in Table 3 that all of the models underestimated the N rd wet depositions, with the negative NMB in both urban and rural sites. Although little difference between rural and urban sites was found in M11 and M14, a better performance in rural areas was manifested by lower NMB and higher FAC2 in rural sites than the urban sites in most of models (−13.6 % to −23.2 % versus −37.3 % to −45.6 % for NMB and 55.9 %-70.6 % versus 42.9 %-55.1 % for FAC2, except M11 and M14). The underestimation of the simulated N rd wet depositions was also found in the VWM concentration of N rd in precipitation with similar NMB and FAC2. However, compared with the wet deposition, the correlation between the simulated and observed N rd VWM concentration in precipitation was significant, with the R increasing from ∼ 0.3 to ∼ 0.8, which was similar to that shown in N ox . This indicates that the simulated VWM concentration of N rd in precipitation by MICS-Asia III has a better performance in the spatial variation than the simulation of N rd wet deposition over China. The underestimation of N rd wet deposition was also found over the whole of East Asia and reported in the overview of acid deposition in MICS-Asia III . This implies that the current CTM models might underestimate the prediction of N rd wet deposition not only in China but also in the whole of East Asia. The close correlations between the atmospheric concentration of NH + 4 and wet deposition of N rd , with overestimation in the atmosphere but underestimation in precipitation, were found over all of East Asia . Consistent relationships in NC were also found between the results of Chen et al. (2019b; overestimated NH + 4 concentration) and this study (underestimated N rd wet deposition). Bae et al. (2012) reported that the below-cloud scavenging process was important in the simulation of N rd wet deposition, which was not explicitly sep- arated into in-cloud and below-cloud scavenging but computes it as a whole in the CMAQ model. Note that the wet scavenging processes in most of the models (including M11 and M12) of MICS-Asia III were similar to those treated in CMAQ module, except for M13 (Table 1). It is too simple to accurately simulate wet deposition with the absence of an accurate below-cloud wet scavenging simulation. This would be one reason for the underestimation of N rd wet deposition, especially considering the high concentration of gaseous ammonia in the surface layer of NC Kong et al., 2019). Besides, the underestimation NH 3 emissions in China would also lead to the underestimation of N rd wet deposition. According to the improved inversion of NH 3 emissions in China by Kong et al. (2019), the significant underestimation of NH 3 emissions was found, especially in the North China Plain (NCP). For the whole of China, the prior emission and the inversion emission of NH 3 are 10.3 and 13.1 Tg yr −1 , respectively.
Specifically, the performance of N rd wet deposition prediction in MICS-Asia III has also been validated in five regions through the percentile boxplot in Fig. 4. Site-by-site validation of both the yearly wet deposition and the VWM concentration of N rd simulated by each model are displayed in Fig. S2 in the Supplement. While different from that found in N ox , almost similar behavior prediction has been found in some models, i.e., CMAQ models in M1-M6, except M12, which was driven by a different meteorological model. Other regional models and a global model show significantly different percentile distributions in all regions. Overall, both the medium and mean values of N rd wet deposition were underestimated in NC, SE and SW and TP, while they are found to be similar in NE and NW. The underestimation in NC was largely due to the underprediction in summertime, not only in urban sites (Fig. 3e) but also in rural sites (Fig. 3f). Unfortunately, we cannot obtain convincing temporal variations in SE and SW due to the scarcity of monthly data in these two regions (there are only one or two sites in each region). In NE, most models predicted similar temporal variations of N rd wet deposition, especially for the high depositions in the summer months. Figure 5 shows the map of the distribution of yearly N ox wet deposition simulated by each participant model, the ENM results and the observed results over China. Most models show similar spatial patterns, with high levels of depositions in central to eastern China and low levels in western China. However, the threshold value in the hotspot areas (Fig. 5, from light yellow to orange and red) varies significantly among the models, and the average is much higher than the N r deposition threshold value of 10 kg N ha −1 compared to the temperate ecosystems suggested by Bleeker et al. (2011). For example, M1, M5 and M11 simulated very high wet depositions of N ox (almost reaching at 20 kg N ha −1 ) in the middle Yangtze River and Yangtze River Delta (YRD), the basin of the Sichuan Province, the south of NC, and the Liaoning Province located in NE. In contrast, M2 and M14 fail to show the relative hotspot of N wox in such areas, and M4, M6, M12 and M13 show an obscure hotspot with a small value of 10 kg N ha −1 . The significant differences not only exist between different models but also in the same model (CMAQ), i.e., M1, M2, M4, M5 and M6. Since most models were driven by the meteoroidal field and standard emission input, except M13 (Geos-Chem) and M14 (RAMS-CMAQ), the differences in simulated N ox wet deposition should come from the CTMs themselves, e.g., the diffusion and convection process, the oxidation and chemical transformation, and the wet scavenging and deposition processes. The comparison of the long-lifetime species of CO (Kong et al., 2020) and weak chemical activity species of BC (Chen et al., 2019b) revealed that the model uncertainties are less than other factors, i.e., O 3  and NO − 3 (Chen et al., 2019b), which have strong chemical activity and short lifetimes in the atmosphere. These results indicate that the chemical reaction process, rather than the diffusion and convection process, is one of the most important factors affecting the model uncertainties in MICS-Asia III.  Figure 6 shows the map distribution of the reduced N (N rd ) wet deposition over China. All of the models show a similar spatial pattern, with high values in central and eastern China but low levels of deposition in NW and northwestern NE. Compared with the N ox , few differences in the simulated N rd wet deposition were found among the nine models, except M11 which predicted significantly lower values. N ox wet deposition of five agriculturally dominant provinces -namely Shandong, Henan, Hubei, Hunan and Anhui -is higher than the threshold value of 10 kg N ha −1 , according to the simulated results by most models. Unfortunately, the small number of observations in these areas makes it harder to validate their findings. Evidence shows the high level of N rd wet deposition over the threshold based on the observations in Hebei, YRD and the Pearl River Delta (PRD). Almost all of the models underpredicted the results in these areas. Liu et al. (2013) reported on the important contribution of N rd to the total N deposition in China based on the long-term national scale of observed nitrogen deposition data. In the predominantly agricultural areas, ammonia emission is the main contributor to N rd deposition (Liu et al., 2011;Kang et al., 2016).

Coefficient of variations for N depositions in MICS-Asia III
Besides the wet deposition of oxidized and reduced N, dry deposition was also an important process for the total deposition in China (Liu et al., 2013;Pan et al., 2012). The coefficient of variation (hereinafter CV), defined as the standard deviation divided by the mean value of all selected model results with a large value denoting lower consistency among the models, is applied for the model comparison of simulated reactive N depositions both for dry and wet deposition processes in MICS-Asia III. Figure 7 shows the distribution of CV for each type of simulated reactive N deposition. Since the low level of mean values of the deposition are more likely to be associated with a higher CV, the gridded CV was only calculated in the area with simulated depositions higher than 0.5 kg N ha −1 (hereafter analyzed value) in this study. As shown in Fig. 7, the spatial distribution of CV only covers eastern China, southern China and northeastern China, which indicates that the quarterly and yearly fluxes of reactive N deposition in these regions were higher than the analyzed value. For the annual case, the CV value of N rd was lower when compared with N ox both for dry and wet depositions. This means that the multimodel simulations are more consistent in N rd depositions than in N ox depositions. Specifically, the N rd in wet depositions have the lowest CV values, followed by N rd in dry depositions and then the N ox in wet and dry depositions. This suggests that the simulated wet depositions of N rd have fewer uncertainties than other types of reactive N depositions. More complicated patterns were shown in seasonal variations of each type of deposition. The simulated N ox for dry deposition in Fig. 7a showed larger uncertainties in southern China (south of 30 • N, with the CV > 0.4) than that in northern China (north of 30 • N, with the CV < 0.3) in all seasons except summer. Similar spatial and temporal patterns of the CV values were found in N rd dry deposition. It is worth noting that the large CV values, with the range of 0.4-0.6, were exhibited in central China (i.e., Henan, Hebei and Shandong provinces) during summer and autumn in spite of the high flux of N rd dry depositions in these regions (map distributions of N ox and N rd dry depositions, simulated in nine participant models, are displayed in Figs. S3 and S4 in the Supplement). This suggests, importantly, that the uncertainties of the physical and chemical processes in the participant models -including gas-particle equilibrium , dry deposition parameter scheme (Zhang et al., 2003), transportation and the chemical reaction with other acidifying substances  -in the regions of high emissions originating from agricultural activities in growing seasons may lead to significant deviations of simulated N rd dry depositions.
For wet deposition of N ox , large uncertainties were located in southern China in summer and autumn, with the CV values higher than 0.6 compared to the CV values lower than 0.4 in other regions (Fig. 7c). This high value of CV was not found in the summertime of simulated N rd wet deposition (Fig. 7d). Due to the high portion of summertime flux in the total annual wet deposition, high CV values in N ox contributed the most important part of the significantly larger annual CV value than that shown in the N rd case. Due to the same rainfall input for the wet deposition in the framework of MICS-Asia III, except for models 13 and 14, the different CV values for N rd and N ox in same region (i.e., lower CV values of N rd wet depositions in NC, SE and central China) can be attributed to their precursor concentrations in the air mass and the different wet scavenging processes (Seinfeld and Pandis, 2006). This will be discussed in the following section.

Comparison of precursors in the air mass
As is well known, the depositions from both the dry and wet parts of a certain substance were found to have origi-nated from its precursor in the air mass. The uncertainties of the nitrogen-related species in the air mass simulated during MICS-Asia III were therefore an important index for estimating the performance of deposition simulations. It should be noted that only concentrations of most of the related species at the surface layer have been submitted to MICS-Asia III, except for NO 2 vertical column density data (VCD). According to the comparison of CV between the NO 2 concentration at the surface layer and VCD data (Fig. S5 in the Supplement), it was shown that there is a similar spatial pattern over the whole of China. This indicates that the simulated concentration on the surface layer can partly represent the difference of the whole column among the multimodels, especially for providing a broad overview in China. Thus, only concentration data at the surface layer have been used in this study. Figure 8 shows the distribution of CV for gaseous NO x , particulate NO − 3 , gaseous NH 3 and particulate NH + 4 in the air mass simulated by the nine participant models during four seasons and the annual mean values. There were significant seasonal variations among the spatial patterns of the CV for each type of the N-related air pollutant. It is interesting to note that not only the seasonal variations but also the spatial patterns of the simulated precursors' CVs were reasonably consistent with those previously shown in the deposition part (Fig. 7). For example, high CV values were found in the simulation of particulate NO − 3 in southern China during summer, reaching or even exceeding 0.8 in SE China (Fig. 8b). The high CV values were also found in summertime N ox wet depositions (Fig. 7c). As the most important precursor of N ox wet deposition , the consistent distribution of CV between the precursor and the deposition is reasonable. Different to the particulate NO − 3 , very low CV values were shown in particulate NH + 4 during summer, leading to less deviation of simulated N rd wet deposition than N ox . However, uncertainties in precursors cannot explain everything. For example, the high CV values of N ox wet deposition in south China correspond to the low CV values of NO − 3 in autumn. Some other factors, such as the scavenging process, might be responsible for the unknown uncertainties. Xu et al. (2017Xu et al. ( , 2019 first compared the below-cloud wet scavenging coefficients based on the different estimation methods and found a large difference, even at magnitude level, between each method. Thus, a more detailed comparison -such as incloud and below-cloud wet scavenging coefficients in each participant model -should be carried out in the next phase of MICS-Asia. For N ox dry depositions, the anomalies of deposition and NO x concentration in the air are shown in Figs. S6 and S7 in the Supplement. It shows some variations among the models, i.e., of both higher deposition and concentration in M1, M5, M11 and M13 and lower deposition in M2, M4, M6, M12 and M14. This has also been proved in the distribution of CV values in NO x (Fig. 8a) and N ox dry depositions (Fig. 7a) during autumn and winter. As the most important precursor of N rd dry deposition, gaseous NH 3 also shows large CV values in central China during summertime (> 0.6). There were also significantly high CV values in the south of the Yangtze River during the autumn and winter period (0.7-  and NO x emission from MICS-Asia (d); ENM results for dry deposition (e) and wet deposition (f) of N rd , VCD of NH 3 from IASI (g) and NH 3 emission from MICS-Asia (h). Figure 10. ENM results of each process of (a) N deposition flux, and the fraction of (b) total inorganic N (TIN) in MICS-Asia III. The labels of w_N, w_A, d_NO2, d_NH3, d_HNO3, d_ammonium and d_nitrate represent the wet deposition of nitrate, wet deposition of ammonium, dry deposition of NH 3 , dry deposition of HNO 3 , dry deposition of ammonium and dry deposition of nitrate, respectively. that the ensemble mean value of simulated N ox wet deposition has large uncertainties, while N rd wet deposition was underpredicted by a factor of 2 in MICS-Asia III. Compared to each single model, the ensemble mean shows a higher R value than most single models. However, due to a lack of direct observations of dry deposition, validation for dry and total deposition of reactive N cannot be achieved. Instead, the column densities from satellites and emission spatial distribution were employed to address the reasonability of the ensemble mean of four types of reactive N depositions simulated in nine models. As displayed in Fig. 9, dry depositions of N ox and N rd are concentrated in NC, the YRD and Henan province, which corresponds to the distribution of their emissions and VCDs, respectively. Meanwhile, wet depositions of N ox and N rd are centered in central China provinces, such as Hubei and Hunan, and the Chengdu region. There were especially high wet depositions of N rd in the southwest of the Hubei province and northeast of Chengdu city, where high values of emissions and the VCDs for NH 3 were absent. The reasons that these regions were loaded with high wet depositions were mainly due to the high volume of rainfall (for more details, see Fig. S12 in the Supplement) and the longrange transport of acidic substances (Ge et al., 2011).

Contributions to the total inorganic N depositions and their potential effects
Total inorganic N (TIN) deposition, which includes the reduced and oxidized forms of inorganic N deposition from both wet and dry processes, has been calculated to estimate its ecosystem effects in this study as it was measured, in most cases, before Liu et al., 2013). Figures 10  and 11 show the pathway of each type of N deposition to the TIN from a spatial distribution view and six regions' statisti-cal results, respectively. The ENM dry depositions of gaseous HNO 3 and NH 3 were the two most important contributors to the TIN, both of which contributed to 18 % of the TIN over the whole country, while the wet deposition of NO − 3 and NH + 4 were another two main components, with percentages of 23 % and 28 % (Table 4), respectively. Consistent with that reported in the global assessment under HTAP (Vet et al., 2014) and in the Nationwide Nitrogen Deposition Monitoring Network (NNDMN) estimation (Xu et al., 2015), the N rd in China dominated the TIN deposition with an average percentage of 52 % for the ensemble results, which is slightly lower when compared with 60 % and 58 % in the two previous works. The overall contribution of wet and dry deposition to TIN was almost half and half, which is consistent with that reported in NNDMN by Xu et al. (2015). Considering the total emissions, the depositions in all of China took about 67 %, 65 % and 66 % in the 2010 emissions of NH 3 , NO x and total N (NH 3 -N + NO x -N), respectively. It is interesting to show that the relationship between the gridded average N rd deposition and the N ox deposition with their relevant emissions in six regions (shown in Fig. 12, with the slope -0.56 and r 2 -0.97 for N rd and the slope -0.47 and r 2 -0.88 for N ox ) were consistent with those reported by Xu et al. (2015;slope -0.51 and r 2 -0.89 for N rd and slope -0.48 and r 2 -0.81 for N ox ). Even the increasing order of the regions, from the lowest in the TP to the highest in NC, was the same as the previous measurement study. This implies that the spatial distribution and the relationships between deposition and emission are comparable with those measured in the NNDMN. Pan et al. (2013) also compared the correlations between the observed depositions and emissions and attributed the inconsistent distribution between them in the NCP to the uncertainties of the emission. However, the patterns of depositions were also influenced by the regional transport, in addition to the emissions. In this study, the significant positive correlations of the simulated N ox (N rd ) depositions with the corresponding NO x (NH 3 ) emissions reflect the control role of the relative emissions to the spatial distribution of the depositions. Although most regions were located below a 1 : 1 ratio of deposition to emission (Fig. 12), a few regions, such as TP and NE, were close to or above the 1 : 1 ratio, thus implicating the impacts of transport on deposition among the regions.
For regions, the area-average deposition of TIN was as high as 29.2 kg N ha −1 and 27 kg N ha −1 in NC and SE, followed by 15 kg N ha −1 and 10.1 kg N ha −1 in SW and NE, respectively. The TIN in NW and TP was as low as 3.1 and 2.7 kg N ha −1 . In the two highest regions of NC and SE, the deposition of TIN was similar, but the pathways to them were different. The N rd deposition (53 %) and the dry deposition (54 %) contributed to more than half of the TIN in NC, while the N ox deposition (55 %) dominated the TIN in SE. Considering the lower ratio of NO x / NH 3 emissions in SE (21.4/21.6, 0.99) than in NC (30.4/24.4, 1.25), higher contribution of N ox to TIN in SE indicated a higher nitrogen oxidant ratio (i.e., the ratio of oxidation from NO 2 to NO − 3 ) than in NC. Our companion paper (Tan et al., 2020) also revealed the higher nitrogen oxidation ratio in SE as 0.4-0.6, as compared with that in NC as 0.2-0.4. For more oxidant N species, i.e., HNO 3 and NO − 3 , both dry and wet depositions were higher in SE than those shown in NC (5.8 versus 4.9 for dry deposition of gaseous HNO 3 and 6.9 versus 6.3 for wet deposition of particulate NO − 3 ). While for less oxidant N and the reduced N, all types of depositions -such as the dry deposition of gaseous NO x , gaseous NH 3 and the particulate NH + 4 -were less in SE than NC, except for the wet deposition of particulate NH + 4 , due to the much higher volume of rainfall in SE (Fig. S5 in the Supplement). Overall, the oxidant N made the emitted NO x easier to scavenge in SE, with the ratio of N ox deposition / NO x emission reaching 70 %, while the reduced N is more likely to be scavenged from its emission with the ratio of 64 % in NC. The total ratio of TIN/emission in NC and SE was 53 % and 63 %, respectively. Compared to the critical load (CL; Duan et al., 2001;Zhao et al., 2009;Liu et al., 2011), which is a judgment of the deposited N effects on the ecosystem, the two regions almost reached and even exceeded to the CL value (Table 4) in some cases, indicating serious ecological impacts of the N deposition in NC and SE. More attention should be given to  Liu et al. (2011). controlling N-related species, especially the N rd in NC and N ox in SE. In the less developed economic and social area of SW, due to the high emission of NH 3 , 60 % of the TIN was contributed by N rd deposition. The ratio of NO x / NH 3 emissions reached 0.49 as more NH 3 was emitted from agricultural activity than NO x from fossil fuel consumption. The ratio of wet deposition/TIN was 55 %, which was lower than the HTAP comparison during 2000 (60 %-70 %; Vet et al., 2014) but higher than the results of NNDMN (45 %; Xu et al., 2015). Although socially undeveloped, the TIN deposition was almost as high as the CL value according to Zhao et al. (2009). Besides the high emission of NH 3 , the high ratio of N ox deposition / NO x emission of up to 94 % reflects the importance of N ox from high emission areas, such as SE and NC, and should attract our attention in this region. Although the N deposition in TP was not as high as CL valuewhich was the lowest in all regions of China, with the value of 2.7 kg N ha −1 -the N ecological impacts cannot be ne-glected due to its sensitive ecosystem (Shen et al., 2019) and the important climatic influence it has on all of China. Considering the high ratio of TIN / emissions, which were larger than 1 : 1 -with 3 : 1 for TIN, 2.71 : 1 for N rd and 4 : 1 for N ox -the imports from outside the region were responsible for the N deposition in TP.

Conclusions
Reactive N depositions over China simulated in the framework of MICS-Asia III have been compared within each participant model. Wet deposition was also validated by multisource observations, i.e., recorded data from EANET, CAS, NNDMN and EPA in the Guangdong and Liaoning provinces. Most models show the consistent spatial and temporal variations of both N ox and N rd wet deposition in China with the NME around 50 %, which is lower than the value of 70 % based on EANET observations over Asia. The coefficient of variation (CV) was applied for model comparison of dry deposition and the related precursor's concentration in the air mass. Consistency of both the spatial and temporal variation of CV in deposition and the concentration in air mass indicates that performance of the precursors' simulation was highly correlated with their deposition.
A large deposition of ensemble simulations of N rd deposition in eastern China corresponds with a low level of VCD from satellite measurements, while the case of N ox was just the contrary. The total emissions of NO x and NH 3 were similar, at 8 kg N ha −1 , in China. This indicates that the allocation of both NO x and NH 3 from the deposition to the ground and the amount staying in the atmosphere were conserved from their emissions into the air, which also suggests that there were reasonable simulations of deposition for N ox and N rd in MICS-Asia III.
Wet deposition of nitrate and ammonium and the dry deposition of gaseous NH 3 and HNO 3 were the important pathways to TIN deposition, with the percentages of 18 %, 18 %, 23 % and 28 % for ensemble results, respectively. The gridded average N rd in China dominated the TIN deposition, with the average percentage found to be 52 %, which is slightly lower than the reported 60 % and 58 % in HTAP and NNDMN measurements. The contribution of wet and dry deposition to TIN was almost half and half and consistent with that reported in NNDMN. Even the ratio of TIN/emissions was similar with the NNDMN, indicating that the spatial distribution and the relationships of deposition and emissions are comparable with those measured in the NNDMN.
For different regions of China, the simulated N ox wet deposition was overestimated in NE but underestimated in SE and SW, while large uncertainties were shown in NC. Two models simulated 30 %-50 % higher N ox wet deposition, and four models were 20 %-40 % lower compared to observations in NC. The large divergences not only exist between different models but also in the same (CMAQ) model, i.e., M1-M6. For the simulation of N rd wet deposition, all the models underpredicted in all regions, with the largest underestimation in NC and SE. While different from N ox , an almost similar behavior prediction of the less oxidative species, such as the N rd wet deposition, has been found in CMAQ models, indicating the chemical reaction process is one of the most important factors affecting the model uncertainties in MICS-Asia III. Compared to the CL value, the reactive N deposition in NC, SE and SW reached or exceeded the reported CL value and indicates serious ecological impacts. The control of N rd in NC and SW and N ox in SE would be effective in mitigating the TIN deposition in these regions. For the lowest reactive N deposition in TP, however, the N ecological impacts cannot be neglected since it has a sensitive ecosystem, and it has an important climatic influence on all of China, especially considering the high ratio of TIN/emissions, which was mainly caused by outside sources. The joint prevention and control of air pollution in China should be carefully considered and implemented in the future.
Data availability. The observed data are available, on request, for scientific research purposes from Baozhu Ge at the Institute of Atmospheric Physics, Chinese Academy of Sciences (gebz@mail.iap.ac.cn).
Author contributions. BG designed the whole structure of this work, and prepared the manuscript with contributions from all coauthors. BG, SI and KS led the deposition analysis group in MICS-Asia III. DX, JW, FF and QT helped with the data processing. JSF, XW, KY, TN, JL, MK, HL and MZ performed the model simulations and contributed by submitting their simulated deposition results. ZW performed the meteorological model simulation and ex-amined the model performance. ML, JW, JK and QW prepared the emission inventory data. YP and XL supported the observation data in China. ZW was involved in the scientific interpretation and discussion.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Regional assessment of air pollution and climate change over East and Southeast Asia: results from MICS-Asia Phase III". It is not associated with a conference.