Analysis of atmospheric ammonia over South and East Asia based on the MOZART-4 model and its comparison with satellite and surface observations

Limited availability of atmospheric ammonia (NH3) observations limits our understanding of controls on its spatial and temporal variability and its interactions with the ecosystem. Here we used the Model for Ozone and Related chemical Tracers version 4 (MOZART-4) global chemistry transport model and the Hemispheric Transport of Air Pollution version 2 (HTAP-v2) emission inventory to simulate global NH3 distribution for the year 2010. We presented a first comparison of the model with monthly averaged satellite distributions and limited ground-based observations available across South Asia. The MOZART-4 simulations over South Asia and East Asia were evaluated with the NH3 retrievals obtained from the Infrared Atmospheric Sounding Interferometer (IASI) satellite and 69 ground-based monitoring stations for air quality across South Asia and 32 groundbased monitoring stations from the Nationwide Nitrogen Deposition Monitoring Network (NNDMN) of China. We identified the northern region of India (Indo-Gangetic Plain, IGP) as a hotspot for NH3 in Asia, both using the model and satellite observations. In general, a close agreement was found between yearly averaged NH3 total columns simulated by the model and IASI satellite measurements over the IGP, South Asia (r = 0.81), and the North China Plain (NCP), East Asia (r = 0.90). However, the MOZART-4-simulated NH3 column was substantially higher over South Asia than East Asia, as compared with the IASI retrievals, which show smaller differences. Model-simulated surface NH3 concentrations indicated smaller concentrations in all seasons than surface NH3 measured by the ground-based observations over South and East Asia, although uncertainties remain in the available surface NH3 measurements. Overall, the comparison of East Asia and South Asia using both MOZART-4 model and satellite observations showed smaller NH3 columns in East Asia compared with South Asia for comparable emissions, indicating rapid dissipation of NH3 due to secondary aerosol formation, which can be explained by larger emissions of acidic precursor gases in East Asia. Published by Copernicus Publications on behalf of the European Geosciences Union. 6390 P. V. Pawar et al.: Analysis of atmospheric ammonia over South and East Asia based on the MOZART-4

Abstract. Limited availability of atmospheric ammonia (NH 3 ) observations limits our understanding of controls on its spatial and temporal variability and its interactions with the ecosystem. Here we used the Model for Ozone and Related chemical Tracers version 4 (MOZART-4) global chemistry transport model and the Hemispheric Transport of Air Pollution version 2 (HTAP-v2) emission inventory to simulate global NH 3 distribution for the year 2010. We presented a first comparison of the model with monthly averaged satellite distributions and limited ground-based observations available across South Asia. The MOZART-4 simulations over South Asia and East Asia were evaluated with the NH 3 retrievals obtained from the Infrared Atmospheric Sounding Interferometer (IASI) satellite and 69 ground-based monitoring stations for air quality across South Asia and 32 groundbased monitoring stations from the Nationwide Nitrogen Deposition Monitoring Network (NNDMN) of China. We identified the northern region of India (Indo-Gangetic Plain, IGP) as a hotspot for NH 3 in Asia, both using the model and satel-lite observations. In general, a close agreement was found between yearly averaged NH 3 total columns simulated by the model and IASI satellite measurements over the IGP, South Asia (r = 0.81), and the North China Plain (NCP), East Asia (r = 0.90). However, the MOZART-4-simulated NH 3 column was substantially higher over South Asia than East Asia, as compared with the IASI retrievals, which show smaller differences. Model-simulated surface NH 3 concentrations indicated smaller concentrations in all seasons than surface NH 3 measured by the ground-based observations over South and East Asia, although uncertainties remain in the available surface NH 3 measurements. Overall, the comparison of East Asia and South Asia using both MOZART-4 model and satellite observations showed smaller NH 3 columns in East Asia compared with South Asia for comparable emissions, indicating rapid dissipation of NH 3 due to secondary aerosol formation, which can be explained by larger emissions of acidic precursor gases in East Asia.

Introduction
Gaseous pollution due to various forms of nitrogen emissions plays an important role in environmental processes. Specifically, ammonia (NH 3 ) emitted from various agricultural activities, such as the use of synthetic fertilisers and animal farming, together with nitrogen oxides (NO x ), is one of the largest sources of reactive nitrogen (N r ) emission to the atmosphere. Ammonia has great environmental implications due to its substantial influence on the global nitrogen cycle and its associated air pollution, on the ecosystem and on public health (Behera et al., 2013;Liu et al., 2017b;Zhou et al., 2016). Emission estimates provided by the latest Emission Database for Global Atmospheric Research (EDGAR v4.3.2) emission inventory suggest that globally about 59 Tg of NH 3 was emitted in the atmosphere in 2012, of which direct soil emissions contributed about 56 %, manure management (on farm) contributed about 19 % and agricultural burning contributed about 1.5 %, while the biomass burning contribution is not included in the emission estimate. Furthermore, due to a lack of observed emission factors and the high uncertainty of agricultural statistics, the uncertainty of NH 3 is the largest among all other pollutants in EDGAR v4.3.2 (Crippa et al., 2018). Ammonia is a key precursor of aerosol formation, as the reactions in the atmosphere lead to an increase in different forms of sulfates and nitrates that contribute to secondary aerosol formation (Pinder et al., 2007(Pinder et al., , 2008. India and China together accounted for an estimated 64 % of the total amount of NH 3 emissions in South Asia during 2000-2014 (Xu et al., 2018). Emissions of NO x and NH 3 are increasing substantially over South Asia (Sutton et al., 2017a, b), which contributes to an increase in particulate mass loading, visibility degradation, acidification and eutrophication (Behera et al., 2013;Ghude et al., 2008Ghude et al., , 2013Ghude et al., , 2016. Asia is responsible for the largest share of global NH 3 emissions (Janssens-Maenhout et al., 2012). A further increase in NH 3 emission will increase its negative impacts and societal cost (Sutton et al., 2017b).
In India, around 50 % of total NH 3 emissions is estimated to be from the fertiliser application and the remainder from livestock and other NH 3 sources (Aneja et al., 2011;Behera et al., 2013). However, there are large uncertainties in emissions of ammonia, its deposition to the surface, its chemistry and its transport (Sutton et al., 2013;Zhu et al., 2015). Urea is mostly used as a fertiliser (Fertiliser Association of India Annual Report 2018-2019, 2018) and alone contributes more than 90 % of total fertiliser used for agricultural activities (Sharma et al., 2008). India is currently the second largest consumer of fertilisers after China, and fertiliser usage is bound to increase with further intensification of agriculture. The fertiliser input of India is expected to have doubled by 2050 (Alexandratos and Bruinsma, 2012).
Recent studies based on Infrared Atmospheric Sounding Interferometer (IASI) satellite measurements show very high concentrations of NH 3 over the Indo-Gangetic Plain (IGP) and the North China Plain (NCP), which were mainly related to agricultural (Van Damme et al., 2014a, b, 2015b and industrial activity Van Damme et al., 2018). The seasonality was shown to be more pronounced in the Northern Hemisphere, with peak columns in spring and summer seasons (Van Damme et al., 2014a). Van Damme et al. (2015a) attempted first to validate IASI NH 3 measurements using existing independent ground-based and airborne datasets. This study does not include comparison of ground-based NH 3 datasets with IASI measurements particularly over South Asia (India) due to the limited availability of NH 3 measurements. Liu et al. (2017a) estimated the groundbased NH 3 concentrations over East Asia, combining IASI NH 3 columns and NH 3 profiles from the Model for Ozone and Related chemical Tracers (MOZART-4) and validated the data with 44 sites of the Chinese Nationwide Nitrogen Deposition Monitoring Network (NNDMN). In one of the recent studies over South Asia, interannual variability of atmospheric NH 3 using IASI observations revealed large seasonal variability in atmospheric NH 3 concentrations, which were equivalent to the highest number of urea fertiliser plants. This study highlights the importance of the role of agriculture statistics and fertiliser consumption/application in determining ammonia concentration in South Asia (Kuttippurath et al., 2020). The available global ammonia emission inventory does not include comprehensive bottom-up NH 3 emissions for South Asia compared to East Asia that are suitable for input to atmospheric models by taking into consideration actual statistical data of various NH 3 sources, such as livestock excreta, fertiliser application, agricultural soil, nitrogen-fixing plants, crop residue compost, biomass burning, urine from rural populations, chemical industry, waste disposal and traffic, which are currently missing (Behera et al., 2013;Huang et al., 2012;Janssens-Maenhout et al., 2015;Li et al., 2017;Zhang et al., 2010). Han et al. (2020) suggested that an updated emission inventory as per the source activity is essential for South Asia to reduce the uncertainties in simulated NH 3 over this region. A recent study by Wang et al. (2020) examined the NH 3 column observed over the IGP during summer using a regional model driven with the MIX emission inventory. The study suggested that high agriculture activity and high summer temperature contribute to high NH 3 emission fluxes over the IGP, which lead to large total columns. A summertime increase in NH 3 concentration at the surface over certain sites in the IGP regions is also observed from the ground-based monitoring network (Datta et al., 2012;Mandal et al., 2013;Saraswati et al., 2019;Sharma et al., 2012Sharma et al., , 2014b. In this study, we examined the spatio-temporal variability of atmospheric NH 3 over Asia (South and East Asia) and focus on two hotspot regions of ammonia, the Indo-Gangetic Plain (IGP) and the North China Plain (NCP). The approach for this study is a combination of simulations using chemical transport modelling, satellite observations and in situ ammonia measurements over South Asia (69 stations) and East Asia (32 stations). The analysis applies the Model for Ozone and Related chemical tracers (MOZART-4) driven by a priori ammonia emissions based on the Hemispheric Transport of Air Pollution version 2 (HTAP-v2) emission inventory. It applies HTAP-v2 data for emissions to produce estimated total columns of NH 3 and aerosol species for the year 2010 over Asia. Model simulations were evaluated and compared with NH 3 data from IASI (over South and East Asia) and selected ground-based observations (noted above). In addition to the regional comparison, we examine why certain emission hotspot regions in East Asia show lower NH 3 total columns compared with similar hotspot regions in South Asia, when analysed with both model and satellite observations.
2 Data and methodology

MOZART-4 model
The global chemical transport model MOZART-4 has been employed in this study to conduct a year-long (2010) simulation of atmospheric trace gases and aerosols over Asia using the updated HTAP-v2 emission inventory (Janssens-Maenhout et al., 2015). These simulations were performed earlier to meet the objectives of the Task Force on Hemispheric Transport of Air Pollution phase 2 multi-model experiments (Surendran et al., 2015(Surendran et al., , 2016. The model domain covers the entire globe at a horizontal grid resolution of 1.9 • × 2.5 • and 56 vertical levels from the surface up to 1 hPa. The model has approximately 10 levels in the boundary layer (below 850 hPa). MOZART-4 takes into account surface emissions, convection, advection, boundary layer transport, photochemistry, and wet and dry deposition. The model simulations were driven by the input meteorological dataset of 1.9 • ×2.5 • resolution from Modern Era Retrospective Analysis for Research and Applications (MERRA) of the Goddard Earth Observing System Data Assimilation System (GEOS-DAS). Model simulations were performed for the complete year of 2010 (1 January to 31 December 2010), and its outputs were saved every 6 h (4 time steps each day) with a spin-up time of 6 months (1 July 2009 to 31 December 2009). MOZART-4 includes 157 gas-phase reactions, 85 gas-phase species, 39 photolysis compounds and 12 bulk aerosol compounds (Emmons et al., 2010). Dry deposition of gases and aerosols was calculated online according to the parameterisation of Wesely (1989), and wet deposition of soluble gases was calculated as described by the method of Emmons et al. (2010). Land use cover (LUC) maps used in MOZART-4 are based on the Advanced Very High Resolution Radiometer (AVHRR) and Moderate Resolution Imaging Spectroradiometer (MODIS) data based on the NCAR Community Land Model (CLM) (Oleson et al., 2010). MOZART-4 represents the land surface as a hierarchy of sub-grid types: glacier, lake, wetland, urban land and vegetated land. The vegetated land is further divided into a mosaic of plant function types (PFTs). These same maps are used for the dry deposition calculations (Emmons et al., 2010;Oleson et al., 2010;Lawrence and Chase, 2007). In MOZART-4 the tropospheric aerosol component is built on the extended work of Tie et al. (2001Tie et al. ( , 2005. The online Fast Tropospheric Ultraviolet Visible (FTUV) scheme, based on the TUV model (Tie et al., 2003), is used for the calculation of photolysis rates in MOZART-4. For long-lived species like CH 4 and H 2 , surface boundary conditions are constrained by observations from NOAA/ESRL/GMD (Dlugokencky et al., 2005(Dlugokencky et al., , 2008Novelli, 1999), and N 2 O concentrations are set to the values as described in the Intergovernmental Panel on Climate Change 2000 report (IPCC, 2000). Biogenic emissions of isoprene and monoterpenes are calculated online using the Model of Emissions of Gases and Aerosols from Nature (MEGAN) (Guenther et al., 2006), using the implementation described by Pfister et al. (2008). Surface moisture flux and all relevant physical parameters are used to calculate water vapour (H 2 O) online. Biomass burning emissions of a wide range of gaseous components, including NH 3 , SO 2 and individual volatile organic compounds, were provided by the Global Fire Emission Database (GFED-v3), determined by scaling the GFED CO 2 emissions by the emission factors provided at 1.9 • ×2.5 • grid resolution (Emmons et al., 2010). In MOZART-4, the ammonium nitrate distribution is determined from NH 3 emissions and the parameterisation of gas-aerosol partitioning using the equilibrium simplified aerosol model (EQSAM) by Metzger et al. (2002), which is a set of approximations to the equilibrium constant calculation (Seinfeld et al., 1998), based on the level of sulfate present. In Metzger et al. (2002), cations other than NH + 4 , e.g. sodium (Na + ), potassium (K + ), calcium (Ca 2+ ), and magnesium (Mg 2+ ), as well as organic acids, have been neglected for the gas-aerosol partitioning calculations. Metzger et al. (2006) found that the NH 3 / NH + 4 ratio (calculated by accounting for ammonium-sulfate-nitrate-sodium-chloride-water system (updated EQSAM2 parameterisation considering mineral elements and organic acids)) was 15 % higher than that calculated from the parameterisation similar to EQSAM. Ammonia has a stronger affinity towards the neutralisation of sulfuric acid (H 2 SO 4 ) than nitric acid (HNO 3 ), whereas the formation of ammonium chloride (NH 4 Cl(s) or (aq)) in atmosphere is unstable and can dissociate reversibly to NH 3 and HCl. These aerosols in both dry and aqueous phase evaporate faster than the corresponding ammonium nitrate (NH 4 NO 3 ) aerosols (Seinfeld and Pandis, 2012). In the current modelling setup, NH 3 / NH + 4 partitioning is mainly controlled by sulfate and subsequently by nitrate. A recent study (Acharja et al., 2020) based on the analysis of water-soluble inorganic chemical ions of PM 1 , PM 2.5 and atmospheric trace gases over the IGP revealed that NH + 4 was one of the dominant ions and collectively with Cl − , NO − 3 and SO − 4 constituted more than 95 % of the measured ionic mass in both PM 1 and PM 2.5 . Remaining ionic species (i.e. Na + , K + , Ca 2+ and Mg 2+ ) formed constituted only about 3 % of the total measured ions. Although major mineral cations (i.e. Na + , K + , Ca 2+ and Mg 2+ ) contribute actively to the neutralisation reaction, but their concentration in the IGP was found to be very low. Also over the NCP, mineral cations contributed less than 5 % in both PM 1 and PM 2.5 (Dao et al., 2014). Furthermore, a recent study by Xu et al. (2017) over East Asia revealed that NH + 4 was the predominant neutralising cation with the highest neutralisation factor (NF) (above 1), whereas Na + , K + , Ca 2+ and Mg 2+ had relatively low contributions (below 0.2). Therefore, consideration of mineral cations and organic acids in the NH 3 / NH + 4 partitioning might be limited and will not have a significant impact on the results of this study.

Emission inventory (HTAP-v2)
The HTAP-v2 bottom-up database is used in this study as an input for anthropogenic emissions of NH 3 for the year 2010 (Janssens-Maenhout et al., 2015). The HTAP-v2 dataset is embedded with the activity data as per harmonised emission factors, international standards and gridded emissions with global proxy data. It includes important point sources providing high spatial resolution and emission grid maps with global coverage. This dataset consists of monthly mean NH 3 emission maps with 0.1 • × 0.1 • grid resolution for the year 2010. The HTAP-v2 dataset is compiled using various regional gridded emission inventories by the Environmental Protection Agency (EPA) for the United States, Environment Canada for Canada, the European Monitoring Evaluation Programme (EMEP) and the Netherlands Organisation for Applied Scientific Research for Europe and Model Intercomparison Study for Asia (MICS Asia) for China, India and other Asian countries. The Emission Database for Global Atmospheric Research (EDGARv4.3) is used for the rest of the world (mainly South America, Africa, Russia and Oceania). The MICS Asia dataset incorporated into the HTAP-v2 dataset includes an anthropogenic emission inventory developed in 2010 , which incorporates several local emission inventories, including the Multi-resolution Emission Inventory for China (MEIC), the NH 3 emission inventory from Peking University (Huang et al., 2012) and the Regional Emission inventory in ASia version 2.1 (REAS2.1) (Kurokawa et al., 2013) for areas where local emission data are not available. A detailed description of HTAP-v2 datasets can be found in Janssens- Maenhout et al. (2015).
For this study, we used emissions from five important sectors, such as agricultural, residential (heating/cooling of buildings and equipment/lighting of buildings and waste treatment), energy (power industry), transport (ground transport) and industries (manufacturing, mining, metal, cement, chemical, solvent industry) for the year 2010. Aircraft and international shipping are not considered for NH 3 emissions in the HTAP-v2 bottom-up database. These emissions also include natural emissions such as soil from the Community Earth System Model (CESM) and biomass burning from the Global Fire Emission Database (GFED-v3) (Randerson et al., 2013). All these emissions are re-gridded to 1.9 • × 2.5 • to match the model resolution.
The spatial distribution of the total NH 3 emissions over Asian region is shown in Fig. 1. It shows the highest emissions over both South and East Asia, especially over the IGP and NCP region (shown with the black box in Fig. 1). The agricultural sector is the main contributor to NH 3 emissions, including management of manure and agricultural soils (application of nitrogen fertilisers, including animal waste). It also includes emissions from livestock and crop cultivation, excluding emissions from agricultural waste burning and savannah burning (Janssens-Maenhout et al., 2015). Minor contributions from the residential sector are also observed for the Asian countries due to biomass combustion and coal burning, which are also included in the emissions. Spatial proxies such as population density, road networks and land use information have been used to allocate the area of emission sources. For the REAS2 emission inventory over India, the agricultural sector follows the spatial proxy of the total population . The use of this approach is expected to be the main source of spatial uncertainty in the estimated NH 3 emissions to the extent that total human population is only approximately correlated with the spatial distribution of fertiliser use and livestock numbers. Seasonal variation of average NH 3 emission over the IGP and NCP region for anthropogenic (HTAP-v2), biomass burning (GFED-v3) and soil emission (CESM) is shown in Fig. 2. Anthropogenic NH 3 emissions do not show any strong seasonal variability over the IGP region; however, over the NCP region, NH 3 emissions show strong seasonality with peak emissions during May-September. It can be seen that the magnitude of peak emissions is 2 times more over the NCP region than the IGP region. On the other hand, seasonality in biomass burning NH 3 emissions is strong over the IGP region, which shows the highest emissions in the spring season (MAM). Also, contribution of NH 3 emissions from the IGP region is significantly higher compared to the NCP region during peak burning season, but the magnitude of biomass burning emission is 6 times lower compared to the magnitude of anthropogenic emissions.

Satellite NH 3 observations
The NH 3 total column data used in this study are derived from the IASI space-borne remote sensing instrument on board MetOp-A, which was launched in 2006 in a polar sun-synchronous orbit. The IASI operates in the thermal infrared spectral range (645-2760 cm −1 ) with a mean local solar overpass time of 09:30 and 21:30 . It covers the globe twice a day, and each observation is composed of 4 pixels with a circular footprint of 12 km diameter at nadir and an elliptical footprint at the end of the swath (20×39 km). IASI is a suitable tool for evaluation of regional and global models due to its relatively high spatial and temporal sampling, and retrieval algorithms have been continuously improved (Whitburn et al., 2016). The NH 3 total column retrievals show satisfactory agreement with monthly averaged integrated ground-based measurements with Fouriertransform infrared (FTIR) column data (Van Damme et al., 2015a). IASI measurements are also found to be consistent with other NH 3 satellite products (Clarisse et al., 2010;Someya et al., 2020;Viatte et al., 2020). In the present study, we have used the ANNI-NH 3 -v2.2R-I dataset for the year 2010, which relies on ECMWF ERA-Interim meteorological input data, along with surface temperature retrieved from a dedicated network (Van Damme et al., 2017). An improved retrieval scheme for IASI spectra relies on the calculation of a dimensionless hyperspectral range index, which is successively converted to the total column and allows for a better identification of weak point sources of atmospheric NH 3 (Van Damme et al., 2017;Whitburn et al., 2016). More details about the IASI satellite and the NH 3 data product are given in Clerbaux et al. (2009), Van Damme et al. (2017 and Whitburn et al. (2016). We have considered the daily NH 3 cloud-free satellite total column data and compared them with the modelled daily NH 3 total column averaging paired observations across the months, seasons and year. We have used only morning overpasses at 09:30 measurements, as the relative errors due to the lower thermal contrast are larger for the night-time measurements (21:30 overpass).
For consistency with satellite retrievals, first the model output (11:30 LT) at each day close to satellite overpass time (09:30 LT) is interpolated in space to the location of valid satellite retrievals. Since the IASI retrieval algorithm only provides total columns, in the second step, we made an unweighted average distribution of the daily paired data to obtain a monthly, seasonal and annual mean value of satellite and model total NH 3 columns at each horizontal resolution of the model (1.9 • × 2.5 • ).

Ground-based observations
To evaluate the model's performance in South Asia, we used hourly NH 3 measurements from the air quality monitoring station (AQMS) network operated by the Central Pollution Control Board (CPCB) across India. The CPCB follows a national programme for sampling of ambient air quality as well as weather parameter measurements. An automatic analyser (continuous) method is adopted at each monitoring location. NH 3 is measured by the chemiluminescence method as NO x following the oxidation of NH 3 to NO x . In this approach, NH 3 is determined from the difference between NO x concentration with and without inclusion of NH 3 oxidation (CPCB, 2011). The quality assurance and control process followed for these air quality monitoring instruments is given by the CPCB (2014, 2020). Surface observations of NH 3 are taken from 69 different stations in South Asia. Most of the NH 3 monitoring stations from India used in the current study are situated in the cities representing the urban environment. Sampling of ambient NH 3 is done through a sampling inlet of 1 m above the roof top of the AQMS container with a height of 2.5 m (Technical specifications for CAAQM station, 2019). The details of these monitoring locations are given in Table S1 (in the Supplement), and the geographical locations are shown in Fig. 3. Out of these stations, 35 locations in Delhi, six in Bangalore, four in Hyderabad and two in Jaipur are averaged to get a single value for the same geographical location, and the remaining 22 locations are considered independently, representing 26 respective cities. Hourly NH 3 concentrations (in µg m −3 ) used in the study are for the duration of 2016 to 2019. The quality control and assurance method, followed by the Central Pollution Control Board (CPCB) for these air quality monitoring stations, is given in the CPCB (2011,2020). The calibration procedures for the NH 3 analyser conform to United States Environmental Protection Agency (USEPA) methodologies and include daily calibration checks, biweekly precision checks and linearity checks every 6 weeks. All analysers undergo full calibration every 6 weeks. For details on the calibration procedure, refer to Technical specifications for CAAQM station (2019) and CPCB (2020). Furthermore, we take the following steps to assure the quality of NH 3 observations from the CPCB network stations. For data quality, we rejected all the observations values below the lowest detection limit of the instrument (1 µg m −3 ) (Technical specifications for CAAQM station, 2019) because most of the sites are situated in the urban environment. For cities where more than one monitoring station is available, we rejected all the observations above 250 µg m −3 at a given site if other sites in the network do not show values outside this range. This step aims to eliminate any short-term local influence that cannot be captured in the models and to retain the regional-scale variability. Second, we removed single peaks characterised by a change of more than 100 µg m −3 in just 1 h for all the data in CPCB monitoring stations. This step filters random fluctuations in the observations. Third, we removed some very high NH 3 values that appeared in the time series right after the missing values. For any given day, we removed the sites from the consideration that either experience instrument malfunction or appear to be very heavily influenced by strong local sources. In order to verify the data quality of the CPCB monitoring site, we have inter-compared the NH 3 measurement at the CPCB monitoring station (R.K. Puram) in Delhi with the NH 3 measurements at Indira Gandhi International (IGI) Airport taken during Winter Fog Experiment (WiFEX) (Ghude et al., 2017) using the Measurement of Aerosols and Gases (MARGA) instrument during the winter season of 2017-2018. More details on the NH 3 measurements using MARGA are available in Acharja et al. (2020). Both sites were situated in the same area of Delhi (less than 1 km apart). Our inter-comparison shows that NH 3 measured at the CPCB monitoring station using the chemiluminescence method are slightly (on an average 9.8 µg m −3 ) on the higher side compared to NH 3 measured by ion chromatography (IC) using MARGA (Fig. S1 in the Supplement). The differences that were observed could partly be related to the different NH 3 measurement techniques and partly to the locations of the two monitoring sites, which were not at exactly the same location. Evidently, the difference of 9.8 µg m −3 indicates that the NH 3 measurements from the CPCB do not suffer from the calibration issue. However, rigorous validation is required in the future with more datasets. Given the presence of relatively high NO x concentrations, especially at urban locations, it is recognised that the measurement of NH 3 by difference (i.e. between NO x and NO x plus oxidised NH 3 ) is a potentially significant source of uncertainty. Future measurement inter-comparisons are planned (rescheduled from 2020 to 2021 because of COVID-19) to allow the chemiluminescence method as used in the Indian network to be compared with a range of other NH 3 measurement methods (Móring et al., 2021; the Global Challenges Research Fund (GCRF).
To further evaluate the model's performance over East Asia, we used monthly mean NH 3 measurements from the 32 stations of the Nationwide Nitrogen Deposition Monitoring Network (NNDMN) of China, operated by China Agricultural University. The details of these monitoring locations are given in Table S2 (in the Supplement), and the geographical locations are shown in Fig. 3. Monthly mean NH 3 concentrations (in µg m −3 ) used in the study are for the duration of 2010 to 2015. Ambient concentrations of gaseous NH 3 were measured using an active Denuder for Long-Term Atmospheric sampling (DELTA) system. More detail about the data product is given by Xu et al. (2019). To compare the model with the observation, simulated NH 3 from the model is compared with the surface-based observations by using bilinear interpolation of model output to the geographical location and elevation of the observational sites.
3 Results and discussion

Annual mean NH 3 total columns over South Asia
Yearly averaged 2010 distribution of NH 3 total columns over Asia simulated by MOZART-4 model and also retrieved with IASI instrument is shown in Fig. 4a and b. The total NH 3 columns simulated by the model show higher tropospheric vertical column densities (TVCDs) of about 0.5-7 × 10 16 molecules cm −2 over the IGP region of India compared to any other regions of Asia. This to an extent justifies the larger range of NH 3 column values for the South Asian model domain, with both more polluted and cleaner conditions. These high TVCDs values coincide with the high nitrogen fertiliser and livestock numbers, as scaled according to human population density in Fig. 1. Spatial differences between model-simulated data and satellite data for NH 3 total column distribution are shown in Fig. 4c. On a quantitative level, the MOZART-4 model is found to overestimate the NH 3 total column compared with IASI by 1-4 × 10 16 molecules cm −2 over South Asia, especially over north-east India and Bangladesh. Conversely, the MOZART-4 model underestimates NH 3 in comparison with IASI over the arid region of north-western India (state of Rajasthan adjacent to Pakistan) and centring on Pakistan. There are several possible reasons for the spatial differences shown in Fig. 4c, including (a) uncertainties in the mapped NH 3 emissions data (e.g. between Afghanistan, Bangladesh, India and Pakistan, due to different relationships between human population and livestock/fertiliser activities); (b) uncertainties related to turbulent mixing and dispersion (this may affect both the simulations in MOZART-4 and the assumed vertical profiles for the IASI retrievals); and (c) uncertainties related to precipitation scavenging of ammonia and ammonium, noting that the eastern part of the IGP is substantially wetter than the western part.
According to Fig. 1, the magnitude of NH 3 emissions over the NCP is similar to the IGP. By contrast, much smaller TVCDs of the NH 3 columns are estimated by MOZART-4 and IASI over the NCP compared with the IGP. The MOZART-4 and IASI estimates are found to be in close agreement, with slightly smaller values estimated by MOZART-4. The possible reasons for the difference in NH 3 concentrations in the IGP and the NCP are discussed in Sect. 3.4. The relationship between modelled and IASIretrieved NH 3 total columns is further analysed in terms of scatter plots in Fig. 5a and b, over the IGP region of South Asia (20-32 • N, 70-95 • E) and the NCP region of East Asia (30-40 • N, 110-120 • E) (rectangular areas shown in Fig. 1). Correlation coefficients (r) between model-and satelliteobserved annual mean total columns over the IGP and the NCP are found to be 0.81 and 0.90 respectively for 2010. This indicates that spatial variability in simulated NH 3 by the model and satellite observation is in closer agreement, both over the IGP and the NCP region. The model-simulated annual mean total NH 3 columns give larger values over the IGP region (normalised mean bias (NMB) = 38 %) as well as over the entire South Asia region (NMB = 44 %), whereas over the NCP region (NMB = −35 %) and the entire East Asia region (NMB = −32 %), the model gives values which are smaller than IASI. Other statistical indicators are summarised in Table 1. Larger estimates of NH 3 columns from an atmospheric chemistry transport model (CTM) compared with IASI were also found in an earlier study for South Asia .
The overall higher value of the model-simulated NH 3 over South Asia compared with IASI could be due to the combination of the uncertainties in both approaches. This includes uncertainties in emissions from the HTAP-v2 datasets used for the model simulations, inaccurate modelling of the chemistry in MOZART-4, errors in dry and wet deposition schemes used in the model and biases inherent to infrared satellite remote sensing. For IASI, firstly, only cloud-free satellite  scenes are processed, which could result in missing some of the NH 3 values during cloudy periods and biomass burning events. Secondly, NH 3 vertical columns retrieved from the IASI observations are actually sampled around 09:30 local time, while the MOZART-4-simulated model output close to overpass time (11:30 LTC) was used. Finally, the retrieval of NH 3 from infrared satellites is sensitive to inaccuracies in the temperature profile, and biases in the IASI L2 temperature profiles can result in biases in the retrieved NH 3 (Whitburn et al., 2016). The HTAP-v2 dataset uses proxy values for agricultural activities (i.e. distributed by human population) instead of actual values for field fertiliser application and livestock excretion over South Asia. This could also result in additional uncertainty of NH 3 emissions from agricultural activities. Further work is ongoing to integrate NH 3 emissions inventories for different countries in South Asia based on national datasets, which should allow the emission-related uncertainties to be reduced in future. Similarly, slight underestimation over East Asia might originate from the countryspecific emission inventory used for China (Huang et al., 2012) in the MOSAIC HTAP-v2 emission inventory and the limitations discussed above. The application of any equilibrium model (EQM) in global atmospheric studies is associated with considerable uncertainties. In MOZART-4 chemistry, the ammonium nitrate distribution is determined from NH 3 emissions and the parameterisation of gas-aerosol par-   Metzger et al. (2002), based on the level of sulfate present. The emission fluxes of SO 2 and NO x in the HTAP-v2 dataset also have large uncertainties over the IGP (Jena et al., 2015b;Wang et al., 2020), which can introduce additional uncertainty in NH 3 / NH + 4 gas-aerosol partitioning. In MOZART-4 chemistry, uncertainty can also be associated with the dry and wet deposition scheme, which can result in overestimation (Emmons et al., 2010). Figure 6 shows the model (left) and IASI satellite (middle) seasonal distributions of NH 3 total columns over Asia. These seasons are represented as 3-month periods: winter, December-January-February (DJF, first row); spring, March-April-May (MAM, second row); summer, June-July-August (JJA, third row); and autumn, September-October-November (SON, fourth row). It can be seen in Fig. 6 that there is larger seasonal variation in IASI NH 3 total columns, while MOZART-4 presents limited seasonality in South Asia compared to better seasonal variation estimated in East Asia, as shown by both IASI and the MOZART-4 model. In general, during autumn, spring, summer and winter seasons, MOZART-4 shows higher NH 3 total columns compared with IASI estimates over most of South Asia. However, this difference is more pronounced during autumn (SON) and winter (DJF) seasons (Fig. 6 right). We have seen that (Fig. 2) anthropogenic emission of NH 3 is nearly the same in all months, and biomass burning has a peak during MAM over South Asia in the MOZART-4 model, whereas seasonality is better represented in NH 3 emission for East Asia.

Seasonal variability of NH 3 total columns
Major drivers in anthropogenic NH 3 seasonal variation include differences in management and timing of fertiliser, which are not well represented in the emission over South Asia (Janssens-Maenhout et al., 2012). This can be expected to have a direct effect on NH 3 total columns over South Asia. It is recognised that NH 3 emission can be strongly affected by both short-term meteorological variation and longer term climatic differences (Sutton et al., 2013). This means that NH 3 emissions may be expected to increase in warm summer conditions than in winter (Battye and Barrows, 2004). However, the magnitude of these emissions is expected to be smaller in comparison with anthropogenic emissions and may not contribute significantly to larger summertime NH 3 columns observed from IASI retrievals over South Asia and East Asia than MOZART-4. Additional drivers in NH 3 seasonal variation include meteorological variation. For example, strong subsidence, lower temperature and lighter winds over South Asia in the autumn and winter months prevent venting of low-altitude pollution to the higher altitudes. This means that emitted air pollutants tend to accumulate close to the source region in wintertime conditions (Ghude et al., 2010(Ghude et al., , 2011. Considering the comparison of IGP with NCP, Figure 6. Seasonal NH 3 total column distribution (×10 16 molecules cm −2 ) in 2010 simulated by  and measured by the IASI satellite (middle). On the right, spatial differences between MOZART-4 and IASI during (top to bottom) winter (DJF), spring (MAM), summer (JJA) and autumn (SON) seasons are shown. accumulation of pollutants in the boundary layer is more pronounced over the IGP region due to flat land topography, and it is higher during winter than the autumn months (Surendran et al., 2016). We saw that simulated mean planetary boundary layer height (PBLH) is lower (approximately 400 m; Fig. S2 in the Supplement), and winds are lighter in winter months, compared to summer months, over South Asia, and particularly over the IGP region (Surendran et al., 2016). Figure 7a and b show the time-height distribution of NH 3 and mean PBLH averaged over the IGP region, respectively. It can be seen that during winter months higher atmospheric stability prevents mixing of boundary layer NH 3 to the free troposphere over the IGP (Fig. 7a), which is reflected in the higher wintertime values of MOZART-4 NH 3 columns. Similarly, a higher NH 3 / NH + 4 ratio (Fig. S3 in the Supplement) and lower dry and wet deposition (Figs. S4 and S5 in the Supplement) of NH 3 over the IGP in winter months enhances the accumulation of NH 3 in the boundary layer compared to summer months. On the other hand, much less NH 3 gets detected by the satellite at the higher altitudes, where detection sensitivity of the satellite is higher than that at the surface (Clarisse et al., 2010). The limited sensitivity of IASI measurements to detect boundary layer NH 3 (Van Damme et al., 2014a) could be one of the reasons for large differences (1-4×10 16 molecules cm −2 ) between MOZART-4 and IASI in winter seasons. Also, the sowing of wheat crops over the IGP involves a higher rate of fertiliser application during peak winter months (Sharma et al., 2014a) that releases a significant quantity of NH 3 into the atmosphere. However, this seasonality is largely missing in the emissions (Fig. 2a), indicating that higher MOZART-4 NH 3 is largely driven by the wintertime meteorology over this region.
It is interesting to note from Fig. 6 (right) that during spring the difference between modelled and observed column NH 3 is smaller over the IGP region compared with the winter season. Heating of the land mass due to large solar incidence suppresses the wintertime subsidence over the IGP and leads to a deeper boundary layer during spring and early summer. It can be seen that ( Fig. 7b and Fig. S2 in the Supplement) the average PBLH is about 1100 m and 600 m deeper during spring and summer compared to winter over the IGP. During this season, significant transport of the boundary layer pollution in the mid-troposphere and upper troposphere due to enhanced convective activities and large-scale vertical motion can be noticed in Fig. 7a and is consistent with the earlier studies over this region (Lal et al., 2014;Surendran et al., 2016). Vertical motion associated with the convective activities is expected to redistribute the NH 3 concentration in the column, which leads to more NH 3 at higher altitudes, where detection sensitivity of the satellite is higher than that at the surface (Clarisse et al., 2010). As a result, more NH 3 gets detected by the satellite and we see less difference between observations and the model over the IGP. This may also partly explain the higher IASI estimates of the NH 3 column for summertime prior to the monsoon season. However, this hypothesis needs to be tested with higher sensitivity experiments as part of future work. During the spring season, MOZART-4 reflects a widespread NH 3 total column from the entire Indian land mass, and IASI observations do capture the increase in the NH 3 total column, at least for the seasonal mean cycle (Fig. 8a). This seasonal maximum in the NH 3 total column identified both in IASI and MOZART-4 over South Asia can be explained by two factors: the meteorology factor and biomass burning emissions. Volatilisation of NH 3 is enhanced with an increase in temperature (Sutton et al., 2013); hence higher temperatures during these drier periods over the IGP partly enhance NH 3 emission to the environment, which is also evident from the soil NH 3 emissions in Fig. 2e and f. However, the magnitude of these emissions is expected to be smaller in comparison with anthropogenic emissions. In the Indian region, emissions from biomass burning (crop residue burning) peak in March to May (Jena et al., 2015a), and emission of NH 3 from biomass burning is at a maximum during this period ( Fig. 2c and d). However, MOZART-4 estimates smaller NH 3 total columns compared with IASI over Myanmar, Laos and Thailand during the period March-May (Fig. 6 right). This period is estimated to be associated with large-scale forest fires (and open crop burning) (Chan, 2017;Wu et al., 2018;Zheng et al., 2017), the effect of which appears to be underestimated in the MOZART-4 simulations. This suggests that the Global Fire Emissions Database (GFED-v3) underestimates the fire emissions over this region, agreeing with  and Huang et al. (2013). During the monsoon (JJA) (Fig. 6  right) and spring (MAM) season, IASI NH 3 total columns are larger than the MOZART-4 estimates over the north-western arid region of South Asia, where monsoon rainfall is lowest (less than 30 cm). On the other hand, NH 3 columns estimated by IASI are lower in the north-western IGP than the MOZART-4 simulations. Figure 8 shows the comparison between IASI and modelled monthly time series of NH 3 total columns over the IGP (20-32 • N, 70-95 • E) and the NCP (30-40 • N, 110-120 • E), respectively (rectangular areas shown in Fig. 1). We found a better consistency between modelled and measured seasonal NH 3 total column over the NCP than the IGP. Monthly NH 3 columns over the IGP show bimodal distribution in the model. However, IASI does not show such bimodal variation. Seasonal statistics show large normalised mean bias (38 %) and poor correlation (r = 0.41) between the model and IASI. The bimodal distribution in NH 3 total columns is partly driven by the biomass burning emissions, which show a major peak in spring and another small peak in autumn ( Fig. 2c and d), and partly by the meteorology, as discussed in the previous section. During monsoon months (JJA), when the whole of South Asia receives significant rainfall, model simulations present lower NH 3 total columns, which is not seen in the IASI observations or in the surface observations (Figs. 8a and 9b) over the IGP. The reason for this discrepancy may be related to the flat NH 3 emission over South Asia (Fig. 2). Usually a large amount of fertilisation application is expected during the warm months of June and July in the IGP, which is not represented in the HTAP-v2 emissions, and therefore lower values in the model during monsoon months are mostly driven by the model meteorology. Lower values observed during the monsoon season in general are attributed to increased wet scavenging of NH 3 due to monsoon rain (Fig. S5 left in the Supplement) and influx of cleaner marine air from the Bay of Bengal and Arabian Sea through south-easterly and south-westerly wind (Ghude et al., 2008). On the other hand, monthly variation in IASI NH 3 total columns over East Asia is found to be captured well by the model (Fig. 8b) and seems to follow the variation observed in the anthropogenic NH 3 emission (Fig. 2), except for the month of July where IASI estimates substantially higher NH 3 total columns than the model. The reason for this peak in the IASI data for July may be related to urea fertiliser application in warm July conditions (see temporal course of the enhanced vegetation index; Li et al., 2014), which seems not to be represented well in the HTAP-v2 emissions. The overall statistics show slight good correlation (r = 0.61) between observed and simulated NH 3 columns and negative normalised mean bias (NMB = −41 %).

Comparison between surface NH 3 measurements and simulated NH 3 concentrations in South and East Asia
To evaluate modelled surface NH 3 concentrations in South Asia, we have used NH 3 surface measurements from 69 monitoring locations over India for the years from 2016 to 2019. As 2010 data were not available, we make the hypothesis that measurement from 2016-2019 can be considered as representative of what was measured in 2010. Out of these stations, 35 locations in Delhi, six in Bangalore, four in Hyderabad and two in Jaipur are averaged to get a single value for the same geographical location, and the remaining 22 locations are considered independently, representing 26 respective cities. Due to the lack of ground-based measurements performed in 2010, the following comparison will mainly be qualitative, although it is estimated that the main spatial features of Indian agriculture and NH 3 emissions will be consistent between 2010 and 2016-2019. As per the RCP8.5 scenario (Kumar et al., 2018), NH 3 emission from South Asia is expected to have increased by less than 20 % from 2010 to 2020. Assuming a linear relationship between emission and surface concentration, it is expected that NH 3 concentrations could be higher by about 10 %-15 % in 2016 to 2019. It is interesting to note that the correlation between annual and monthly mean MOZART-4-simulated and surfacemeasured NH 3 concentration (r = 0.82 and r = 0.62) is better than the comparison between MOZART-4 and IASI for South Asia (Fig. 9). However, MOZART-4 has systematically smaller estimated NH 3 concentrations compared with the ground-based measurement network (NMB = −47 %). It should be noted that most of the monitoring stations are situated in urban regions (cities) of India and therefore represent the urban environment, which may have higher NH 3 concentrations locally due to traffic and human activities (Sharma et al., 2014b). Since the MOZART-4 model is run at relatively coarse (1.9 • × 2.5 • ) grid resolution, the emissions may not capture the true variability in emissions at city scale. These surface NH 3 sites are influenced by local emissions that are therefore not resolved by the MOZART-4 model. Therefore, when comparing coarse-scale models to observations, the model may have difficulties in resolving local-scale effects (Surendran et al., 2015). Until the planned further evaluation of the chemiluminescence monitoring method for ammonia (measured by difference with NO x ) is evaluated (as noted in Sect. 2.4), it is not possible to be certain of the extent to which possible uncertainties in the measurement method contribute to the differences shown in Fig. 9b. While noting these uncertainties, it is worth noting that the ground-based NH 3 observation network confirms the occurrence of higher ground-level NH 3 concentrations in autumn and winter, as simulated using MOZART-4 using the HTAP-v2 emissions inventory (Fig. 9b).
Comparison of Figs. 8a and 9b shows that the time course of ground-level NH 3 concentrations (as estimated by MOZART-4) is significantly different than the time course of the total NH 3 column (as also estimated by MOZART-4), whereas the total column is largest in the summer (reflective of deeper atmospheric mixing and recirculation), and the ground-level concentrations are largest during winter. Although it is not easy to use the IASI data to infer groundlevel NH 3 concentrations, the stronger summer maximum of IASI (Fig. 8a) compared with MOZART-4 suggests that IASI would be in less close agreement with the ground-based measurement network than MOZART-4 (Fig. 9b). While recognising uncertainties in this interpretation, the key point is that large NH 3 columns estimated by IASI for May-July are not reflected in the ground-based NH 3 measurements from the Indian monitoring network. Figure 10 shows the comparison between monthly mean (from 2010 to 2015 observations) NH 3 surface measurements from 32 monitoring locations over China and modelled surface NH 3 concentrations from the same location over East Asia. Similar to South Asia, MOZART-4 has systematically smaller estimated NH 3 concentrations compared with the ground-based measurement network (NMB = −44 %) over East Asia. Figure 10b shows that maximum NH 3 concentration occurred in summer (JJA), demonstrating agreement with IASI measurements. Other statistical indicators are summarised in Table 2. Furthermore, high NH 3 concentration from ground-based measurements during JJA is consistent with the higher HTAP-v2 emissions (Fig. 2) (Huang et al., 2012) and higher NH 4 NO 3 concentration (Fig. S6 in the Supplement). Higher concentration of NH 4 NO 3 can also lead to higher NH 3 concentrations, especially during summer due to its semi-volatile and unstable character at higher temperatures, as observed in East Asia. This implies that the NH 3 emissions may play a vital role in determining the seasonal pattern of the ground NH 3 con- centrations. The summer peak may originate from fertiliser application, livestock emissions and volatilisation of NH 3 , which is enhanced at higher temperatures (Liu et al., 2017a).
3.4 Why were NH 3 total columns low over high NH 3 emission over East Asia compared to the high NH 3 emission over South Asia?
Fine-scale details of the NH 3 emissions over Asia in Figs. 1 and 2 clearly revealed larger emission values in areas where there is intensive agricultural management. This is the case especially in the NCP and the IGP (Fig. 1, shown with box). Earlier emission estimates suggest that fertiliser application and livestock contribute 2.6 Tg per year (yr −1 ) and 1.7 Tg yr −1 NH 3 emissions respectively from South Asia (Aneja et al., 2011). Over South Asia, urea accounts for   Fig. 1) account for ∼ 45 % of the NH 3 emitted from fertilisation in East Asia and South Asia (Huang et al., 2012). We find that satellite observations show larger NH 3 columns over the IGP than over similar higher emission regions of the NCP. However, in addition, we also find that the MOZART-4 model is able to capture these contrasting columnar NH 3 levels between the IGP and the NCP. This indicates that the difference between the IGP and the NCP is unrelated to differences between the mosaic of emissions over South Asia and East Asia in HTAP-v2 and similarly not related to uncertainties in satellite retrievals. Instead, the analysis from MOZART-4 demonstrates that the difference can be explained by differences in atmospheric chemistry between the two regions, linked to higher SO 2 and NO x emissions in the NCP than in the IGP. A recent study by  shows that emission fluxes of SO 2 and NO x over the IGP are only one-quarter of those over the NCP.
As ammonia is a highly alkaline gas with an atmospheric lifetime usually of few hours (and rarely a few days) (Dammers et al., 2019), it readily reacts with acid present in the atmosphere to form aerosols, which are eventually deposited to the earth's surface by either dry or wet deposition processes (Figs. S4 and S5 in the Supplement). In the atmosphere, ammonia therefore reacts rapidly with atmospheric sulfuric acid (H 2 SO 4 ), nitric acids (HNO 3 ) and hydrochloric acid (HCl) to contribute to ambient levels of fine particles, forming ammonium sulfate, ammonium nitrate and ammonium chloride, shown in Reactions (R1) and (R2): In the atmosphere, the ammonium ion (NH + 4 ) as an aerosol is estimated to have a lifetime of about 1-15 d (Aneja et al., 1998), though this is obviously dependent on the amount of atmospheric acids (Seinfeld and Pandis, 2012). In addition to the large fertiliser application and livestock management activities which are characteristic of both the IGP and the NCP, industrial and transportation activities are higher over the NCP (China), which also result in higher emission of NO x and SO 2 over the NCP compared with the IGP (Zhao et al., 2013). Ammonia has greater affinity towards oxides of sulfur; hence it first reacts to form ammonium sulfate, and then the remaining ammonia further reacts to form ammonium nitrate (Seinfeld et al., 1998). The differences in the secondary aerosol formation over the NCP and the IGP are compared by considering the MOZART-4 model estimates of volume mixing ratio (VMR) in parts per billion (×10 9 ppb) of total sulfate, ammonium, ammonium nitrate at the surface and the total column of NO x (Fig. 11). Although vertical profiles of the aerosol components are small, there are strong vertical gradients in NO x concentrations, and for this reason we consider the comparison with the total NO x column more reflective of overall NO x chemistry than the ground-level NO x VMR. Figure 11 shows that the total sulfate VMR (Fig. 11a) and the NO x total column (Fig. 11c) are significantly higher over the NCP region than the IGP. Similarly, the total ammonium VMR (Fig. 11b) is significantly larger over the NCP than the IGP, indicating how a higher fraction of the gaseous ammo- Figure 11. MOZART-4-simulated spatial distribution of annual averaged (a) total sulfate aerosol (×10 9 ppb), (b) total ammonium aerosol (×10 9 ppb), (c) NO x total columns (×10 16 molecules cm −2 ) and (d) total ammonium nitrate aerosol (×10 9 ppb) over Asia.
nia is transformed to form ammonium over the NCP region. In addition, Fig. 11d shows higher estimated levels of ammonium nitrate in MOZART-4 over the NCP, reflective of the higher NO x emissions in this region. As a consequence of the different SO 2 and NO x sources, gaseous NH 3 is more quickly removed from the atmosphere over East Asia, with a residence time of approximately 6 h ( Fig. S7 in the Supplement) (higher values indicates lower mean residence time), which is reflected in the higher VMR of ammonium, sulfate and ammonium nitrate (Fig. 11a, b and d). It can be seen that the NH 3 / NH + 4 ratio denotes lower values of 0-1 (Fig. S3 in the Supplement) over East Asia than South Asia, suggesting that NH + 4 partitioning is more over East Asia. As a result, the NH 3 total columns over the NCP are much smaller than over the IGP, even though the magnitude of NH 3 emission fluxes is greater over the NCP than the IGP. This difference indicates that the high NH 3 loading over the IGP is partly coming from the low gas-particle partitioning of NH 3 caused by low SO 2 and NO x emission over South Asia. In contrast, high SO 2 and NO x emissions promote the conversion of gaseous NH 3 into particulate ammonium in the NCP. However, a rapid decline of acidic (SO 2 ) emissions over China after 2000, which may not be reflected correctly in HTAP_v2 (Mortier et al., 2020;Tong et al., 2020;Zheng et al., 2018), will lead to higher NH 3 loading due to less partitioning of NH 3 .

Conclusions
In this work, we have compared NH 3 total columns simulated by the MOZART-4 model with IASI NH 3 satellite observations over South and East Asia. The annual mean distribution reveals a consistent spatial pattern between MOZART-4 and IASI, but MOZART-4 tends to show larger NH 3 columns over South Asia than IASI, particularly over the Indo-Gangetic Plain (IGP), whereas it is in close agreement over East Asia (including the North China Plain, NCP), with the exception of a July peak seen in the IASI dataset, which may be related to the specific timing of fertiliser-related NH 3 emissions. Comparison of the seasonally and monthly resolved IASI total column with the MOZART-4 simulations shows inconsistencies in the spatial and temporal pattern over South Asia. This inconsistency is due to the uncertainties in the emission estimate, which does not include the seasonality pattern in HTAP-v2 over South Asia, as well as uncertainties in the processing of the IASI data. Both the MOZART-4 results and IASI estimates involve assumptions that could considerably affect the comparison between total columns of NH 3 .
In a comparison with estimates from a ground-based NH 3 monitoring network for both South and East Asia, our results showed that MOZART-4 systematically gives smaller NH 3 concentration estimates than the monitoring network. The NH 3 measurement sites used in present study mostly represent urban locations, and the model may not be able to capture the actual concentration at point locations due to the coarser grid resolution over India. In addition, further assessment is needed to demonstrate the reliability of the NH 3 measurement technique used in the monitoring network, where NH 3 is measured by difference with NO x concentrations, which may be uncertain in urban areas with high NO x concentrations.
Despite the high NH 3 emission over both South and East Asia, a larger NH 3 total column is observed over South Asia in both the IASI and MOZART-4 estimates. This difference is explained by the MOZART-4 simulation, which treats the full atmospheric chemistry interaction with SO 2 and NO x emissions, leading to aerosol formation. The MOZART-4 model showed a higher sulfate volume mixing ratio and NO x total column over East Asia, especially in the NCP, which is reflected in the ammonium aerosol volume mixing ratio (VMR) over East Asia. This suggests that the formation of ammonium aerosols (dominated by ammonium, sulfate and ammonium nitrate) is quicker over East Asia than in South Asia, leading to lower NH 3 total columns in East Asia.
To examine the present findings, future studies should investigate the effect of changing emissions of NO x and SO 2 on NH 3 columns, for example by using perturbation of these emissions through counterfactual modelling scenarios. The comparison between model simulations using MOZART-4, satellite-derived estimates from IASI and ground-based monitoring of NH 3 concentrations has highlighted the known uncertainties in emissions, satellite retrievals and measurements at point locations. In order to reduce the uncertainties in ammonia emission, it would be key to create an NH 3 emission inventory specifically over South Asia, which is now currently under development as part of the GCRF South Asian Nitrogen Hub. This includes work to improve the bottom-up NH 3 emission inventory, taking into account primary agricultural statistics on fertiliser use and animal number distributions. There is also potential for top-down (inverse) constraints on NH 3 and NO x emissions by taking inference from the model, satellite and ground-based evidence. Here it is essential to recognise the need for more ground-based observational sites to measure NH 3 air concentrations in rural areas, where agriculture activity is predominant. Such measurements at present are currently very few in South Asia. Coarser global models fail to resolve the local-scale emissions; hence higher resolution regional models with advance chemistry are also needed to resolve the sources and chemical processes on urban and rural scales.
Author contributions. All authors contributed to the research. SDG designed the research; PVP conducted the research; PVP and SDG wrote the paper; CJ and DS performed the MOZART model simulations; AM and MAS formulated the research; MVD, LC, and PFC performed the IASI experiments; and SK, DML, GG, XL, WX, JJ, and TKA contributed to writing.
Review statement. This paper was edited by Frank Dentener and reviewed by three anonymous referees.