Chemical Characteristics of Water-Soluble Inorganic Ions in Different Types of Asian Dust in Wajima, a Background Site in Japan

: Two Asian dust (AD) events were observed in March 2021 (AD1: 16 March 2021 00:00 UTC~17 March 2021 12:00 UTC and AD2: 28 March 2021 00:00 UTC~31 March 2021 12:00 UTC). To determine the chemical characteristics of water-soluble inorganic ions (WSIIs) in different types of Asian dust, the total suspended particulates (TSP) were collected at Kanazawa University Wajima Air Monitoring Station (KUWAMS), a background site in Japan from 27 February to 4 March, 2021. Based on the lidar observations and the backwards trajectory analysis results, AD events were divided into two types: ADN (aerosols were mainly mineral dust) and ADP (aerosols were mixtures of spherical particles). During ADs, the concentrations of the TSP and WSII increased, with the highest TSP concentration in ADN (38.6 µ g/m 3 ) and the highest WSII concentration in ADP (5.82 µ g/m 3 ). The increase in (cations)/(anions) during AD indicates that the input of AD aerosol buffered the aerosol acidity. Additionally, a signiﬁcant increase in Cl depletion, along with ADN events, was found (Cl depletion = 73.8%). To comprehensively analyse the different types of ADs on WSIIs, we refer to the previous data from 2010 to 2015 at KUWAMS. As a result, the increased Cl depletion was caused by the heterogeneous reaction of HNO 3 with sea salt when the air mass passed over the Japanese Sea. Additionally, the chemical form of SO 42 − was highly dependent on the source and pathway, while SO 42 − mainly came from natural soil dust in ADN and from anthropogenic emissions in ADP. The enhancement of secondary NO 3 − was observed in AD via the heterogeneous hydrolysis of N 2 O 5 .


Introduction
Mineral dust plays a nonnegligible role in increasing global and regional aerosol loading [1][2][3]. According to the Intergovernmental Panel on Climate Change, the global quantity of dust released into the atmosphere is estimated to be approximately 2150 Tg/y, accounting for approximately 40% of the annual emissions of global aerosols from desert areas, such as the Gobi Desert, Taklimakan Desert, Tengger Desert, Loess Plateau, etc. [4][5][6][7][8]. TSP samples were collected starting at 9:00 a.m. by using a high-volume air sampler (AH-600, Sibata Sci. Tech. Ltd., Saitama, Japan) with a quartz fibre filter (8 inch × 10 inch, 2500QAT-UP, Pallflex Products, Putnam, CT, USA) at a flow rate of 700 L/min continuously from 27 February to 3 April, 2021. The filters were replaced once a week. After sampling, the filters were dried in a desiccator, kept in the dark, weighed and kept in the refrigerator (−20 °C) until use.
The pre-treatment methods were detailed in our previous study [37]. Briefly, each filter was cut into small pieces and placed in centrifuge tubes. After adding 15 mL ultrapure water, the compounds were sealed and ultrasonicated for 20 min. The extract was filtered through a 0.45-μm micropore filter (Millex ® -HP, Merck KGaA, Darmstadt, Germany) to remove particulate impurities. The solution was divided into two parts for the analysis of cations and anions. An ion chromatography system (883 Basic IC plus, Metrohm, Herisau, Switzerland), consisting of a conductivity detector, a suppressor, guard columns (Metrosep C 4 S-guard/4.0 for cations and Metrosep A Supp 5 S-Guard/4.0 for anions, Metrohm, Herisau, Switzerland) and separation columns (Metrosep C 4-250/4.0 for cations and Metrosep A Supp 5-250/4.0 for anions, Metrohm, Herisau, Switzerland), was used to analyse the above WSIIs.

Backwards Trajectory Analysis
The Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT-4, Windows-based version, NOAA Air Resources Laboratory) was used to identify the potential source and transport pathway of the air mass during the AD periods. Three-day backwards trajectories of arrival at KUWAMS (37.4° N, 136.9° E) during the two AD events were used. Trajectories were set at two different arrival heights of 500 m and 2000 m due to the different heights of the AD layers over KUWAMS ( Figure S1). The global data assimilation system (GDAS 1°) of days in 2021 provided by the National Centers for Environmental Prediction (NCEP) was used. For the first AD, the starting time was set as 16 March 00:00 UTC, and for the second AD, the starting time was 28 March 00:00 UTC.
The back trajectories were computed every 2 h in two AD events. Each trajectory is defined as a cluster. Every two clusters (trajectories) were combined as trajectory pairs, and the cluster spatial variance was calculated. The sum of the spatial variances of all clusters was defined as the total spatial variance (TSV). Cluster combinations of many different combinations perform cluster iteration. Most of the time, the clustering usually increases at a roughly constant rate. When it increases rapidly again, it indicates that the
The pre-treatment methods were detailed in our previous study [37]. Briefly, each filter was cut into small pieces and placed in centrifuge tubes. After adding 15 mL ultrapure water, the compounds were sealed and ultrasonicated for 20 min. The extract was filtered through a 0.45-µm micropore filter (Millex ® -HP, Merck KGaA, Darmstadt, Germany) to remove particulate impurities. The solution was divided into two parts for the analysis of cations and anions. An ion chromatography system (883 Basic IC plus, Metrohm, Herisau, Switzerland), consisting of a conductivity detector, a suppressor, guard columns (Metrosep C 4 S-guard/4.0 for cations and Metrosep A Supp 5 S-Guard/4.0 for anions, Metrohm, Herisau, Switzerland) and separation columns (Metrosep C 4-250/4.0 for cations and Metrosep A Supp 5-250/4.0 for anions, Metrohm, Herisau, Switzerland), was used to analyse the above WSIIs.

Backwards Trajectory Analysis
The Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT-4, Windows-based version, NOAA Air Resources Laboratory) was used to identify the potential source and transport pathway of the air mass during the AD periods. Three-day backwards trajectories of arrival at KUWAMS (37.4 • N, 136.9 • E) during the two AD events were used. Trajectories were set at two different arrival heights of 500 m and 2000 m due to the different heights of the AD layers over KUWAMS ( Figure S1). The global data assimilation system (GDAS 1 • ) of days in 2021 provided by the National Centers for Environmental Prediction (NCEP) was used. For the first AD, the starting time was set as 16 March 00:00 UTC, and for the second AD, the starting time was 28 March 00:00 UTC.
The back trajectories were computed every 2 h in two AD events. Each trajectory is defined as a cluster. Every two clusters (trajectories) were combined as trajectory pairs, and the cluster spatial variance was calculated. The sum of the spatial variances of all clusters was defined as the total spatial variance (TSV). Cluster combinations of many different combinations perform cluster iteration. Most of the time, the clustering usually increases at a roughly constant rate. When it increases rapidly again, it indicates that the clusters being combined are not very similar. Based on the case of a significant increase in TSV, the appropriate number of clusters in both ADs was 3.

Non-Sea Salt Analysis
The non-sea salt concentration is usually calculated by using the total Na + as a sea spray maker. However, based on our previous study on natural aeolian dust particles, the concentration of NaCl in the Gobi Desert was high [38]. Therefore, the Na + collected at KUWAMS was also from continental dust (nss-Na + ), and sea spray and sea salt Na + (ss-Na + ) can be calculated using the four-equation system reported below. The mean Ca 2+ /Na + ratio in the crust ((Na + /Ca 2+ ) = 1.78 w/w) and the mean Ca 2+ /Na + ratio in sea water ((Ca 2+ /Na + ) = 0.038 w/w) were referred [39].

Quality Assurance and Quality Control
The detection limits of Na + , Ca 2+ , Mg 2+ , K + , NH 4 + , Cl − , SO 4 2− and NO 3 − were 0.28, 0.13, 0.30, 0.43, 0.14, 0.23, 0.02 and 0.08 mg/L, respectively. The standard curves were linear, and the fitting degrees (R 2 ) of the standard curves were all greater than 0.99. After every 20 samples, a random replicate check was performed, and the relative standard deviation (RSD) of each ion was less than 5% in the reproducibility testing. The test results of the standard materials (Multication Standard Solution III 137-1461 and Anion Mixture Standard Solution 1019-24011, FUJIFILM, Osaka, Japan) indicated that the recoveries ranged from 89.9% to 102.6%. All WSIIs were undetectable in the blank samples [27].

Optical Properties and Time-Height Evolution of AD
The optical properties and vertical distribution of AD aerosols were further identified by referring to data from the Toyama lidar obverse station (TLO) [41]. The linear distance between TLO and KUWAMS is 80.13 km (Figure 1).
In this study, two AD events were observed at KUWAMS from 16 March 00:00 UTC to 17 March 12:00 UTC (AD1) and from 28 March 00:00 UTC to 31 March 12:00 UTC (AD2) in 2021 ( Figure S1). The dust layer can be identified by the depolarisation ratio at 532 nm ( Figure S1). Since aerosols are assumed to be a mixture of two components: one component is mineral dust, mainly dust particles, and the other component includes spherical particles such as anthropogenic particles and sea salt droplets. The observed total particle extinction (δ p532 ) should be divided into the dust extinction coefficient and the spherical particle extinction [42]. During AD1, the aerosol was regarded as a mixture of dust particles and spherical particles, and AD2 was mainly contributed by AD particles. Additionally, AD1 and AD2 had different vertical distributions of the extinction coefficient. During AD1, dust extinction was distributed from 1 km above ground level (AGL), and the extinction coefficient was approximately 0.15-0.35/km. AD2 shows strong dust extinction as the extinction coefficient reached 0.5/km lower than 2 km and approximately 0.2 in height from 2 km to 5 km ( Figure S1). A previous study proved that dust aerosols at different altitudes may originate from different source regions [43]. Therefore, a backwards trajectory analysis during two AD events was carried out in different settings. For AD1, the starting time was set at 16 March 00:00, 2021, at 500 m, with a total run time of 72 h. As shown in Figure 2a, the major air mass originated from the North-eastern China region at 1 km and was then transported at a low altitude (<2.5 km) over North-eastern China [44]. As shown in Figure S1, AD2 had a higher dust layer around 5 km at TLO, and the strong dust extinction was distributed around 2 km. Therefore, we hypothesised that the air mass with AD2 during outbreaking at the source region might be higher. To figure out the source of AD over KUWAMS, we set the starting height at 2 km. The starting time was 28 March 00:00, 2021, with a total running time of 72 h. As shown in Figure 2b, the dust aerosol was mainly from the Gobi Desert in Mongolia at a higher altitude (>3 km) and then passed over Northern eastern China.
particle extinction [42]. During AD1, the aerosol was regarded as a mixture of dust particles and spherical particles, and AD2 was mainly contributed by AD particles. Additionally, AD1 and AD2 had different vertical distributions of the extinction coefficient. During AD1, dust extinction was distributed from 1 km above ground level (AGL), and the extinction coefficient was approximately 0.15-0.35/km. AD2 shows strong dust extinction as the extinction coefficient reached 0.5/km lower than 2 km and approximately 0.2 in height from 2 km to 5 km ( Figure S1).
A previous study proved that dust aerosols at different altitudes may originate from different source regions [43]. Therefore, a backwards trajectory analysis during two AD events was carried out in different settings. For AD1, the starting time was set at 16 March 00:00, 2021, at 500 m, with a total run time of 72 h. As shown in Figure 2a, the major air mass originated from the North-eastern China region at 1 km and was then transported at a low altitude (<2.5 km) over North-eastern China [44]. As shown in Figure S1, AD2 had a higher dust layer around 5 km at TLO, and the strong dust extinction was distributed around 2 km. Therefore, we hypothesised that the air mass with AD2 during outbreaking at the source region might be higher. To figure out the source of AD over KU-WAMS, we set the starting height at 2 km. The starting time was 28 March 00:00, 2021, with a total running time of 72 h. As shown in Figure 2b, the dust aerosol was mainly from the Gobi Desert in Mongolia at a higher altitude (>3 km) and then passed over Northern eastern China.

Concentration of TSP
Since our sampling was weekly, we used dust period 1 (D1) and dust period 2 (D2) to represent AD1 and AD2, respectively. Figure 3a shows the concentration variations of TSP in the sampling period. The TSP concentration increased due to the input of AD aerosols, with the highest TSP concentration in D2 (38.5 μg/m 3 ), followed by D1 (18.7 μg/m 3 ) and non-AD (NAD; 14.7 ± 3.75 μg/m 3 ).

Concentration of TSP
Since our sampling was weekly, we used dust period 1 (D1) and dust period 2 (D2) to represent AD1 and AD2, respectively. Figure 3a shows the concentration variations of TSP in the sampling period. The TSP concentration increased due to the input of AD aerosols, with the highest TSP concentration in D2 (38.5 µg/m 3 ), followed by D1 (18.7 µg/m 3 ) and non-AD (NAD; 14.7 ± 3.75 µg/m 3 ). We also compared the particle concentrations at Ulaanbaatar (Mongolia) and several cities in Northern China to determine the atmospheric behaviour during ADs [45, 46]. As shown in Figure S2, the first AD broke out on 14 March in Ulaanbaatar and then quickly settled in Hohhot with the PM2.5 concentration increased significantly. However, the PM10 concentration in Ulaanbaatar and the PM2.5 concentration increased at the same time on We also compared the particle concentrations at Ulaanbaatar (Mongolia) and several cities in Northern China to determine the atmospheric behaviour during ADs [45,46]. As shown in Figure S2, the first AD broke out on 14 March in Ulaanbaatar and then quickly settled in Hohhot with the PM 2.5 concentration increased significantly. However, the PM 10 concentration in Ulaanbaatar and the PM 2.5 concentration increased at the same time on 27 March and were lower than those in the first AD, and then, a significant increase in the TSP concentration at KUWAMS was observed. Combined with the sea-surface pressure and backwards trajectory analysis during the two ADs, AD2 broke out at a higher pressure and then transported at a faster speed in the free troposphere, resulting in lower AD particle loss during the long-range transportation ( Figure S3) [47]. Therefore, the TSP concentration in AD2 was the highest. The TSP concentration in AD1 did not show an obvious increase, because most AD particles had already been removed from the air mass. Figure 3b shows the concentration variation of WSIIs during the sampling period. The WSII concentration was highest in D1 (5.82 µg/m 3 ), followed by NAD (5.19 ± 2.27 µg/m 3 ), and lowest in D2 (4.51 µg/m 3 ). This result suggests that the WSII concentration was determined by the source and long-range transport pathway of AD. During D1, the air mass was from North-eastern China, which has high anthropogenic pollutants due to the large consumption of coal, fossil fuels and biomass [48][49][50]. During long-range transportation, the air mass was able to mix with these anthropogenic pollutants, therefore resulting in higher concentration in D1. However, in D2, as we discussed before, the air mass rarely mixed with other types of aerosols. The WSII concentration in D2 was diluted by the additional significantly increased AD particles.

Concentration of WSIIs
Since NaCl exists in both the ocean and the Gobi Desert. In order to better illustrate the effect of AD, the effect of sea salt WSIIs was excluded in this study. As shown in Table S1, in the NAD and D1 periods, the ss-Na + concentration (NAD = 0.71 ± 0.39 µg/m 3 , D1 = 0.75 µg/m 3 ) was higher than the nss-Na + concentration (NAD = 0.17 ± 0.14 µg/m 3 , D1 = 0.31 µg/m 3 ), indicating that Na + was more contributed to by sea salt. However, in D2, nss-Na + was the predominant contributor, with a ss-Na + concentration of nearly 0 µg/m 3 . This result suggests that the Na + concentration will increase when the AD air mass mainly contains soil particles that originate from the Gobi Desert. As shown in Figure S7a, we also found that Cl depletion in D1 (73.8%) and D2 (32.4%) was higher than that in NAD (26.5%), suggesting that the reaction between Cl − and acidic components was enhanced. As a marker of dust particles, the nss-Ca 2+ concentration was higher in D1 (0.18 µg/m 3 ) and D2 (0.37 µg/m 3 ) than in NAD (0.10 ± 0.08 µg/m 3 ), reflecting the intensity of AD. The nss-K + concentration and nss-Mg 2+ concentration also increased as AD occurred.
As shown in Figure 3b, during the sampling period, secondary WSIIs (NO 3 − , NH 4 + and SO 4 2+ ) were the major components, suggesting the strong effects of anthropogenic pollution from the Asian continent, which is consistent with previous studies [28,29]. However, as shown in Table S1, the secondary WSII concentration did not necessarily increase with AD occurrence, since secondary WSIIs may also be affected by chemical reactions in the atmosphere. Although the concentration of NH 4 + was lower than those of nss-SO 4 2− or NO 3 − , NH 4 + played an important role in neutralising aerosol acidity [51].

Aerosol Acidity and Neutralisation Factor
The acidity of atmospheric aerosols is determined by the cation (CA) and anion (AN) micro-equivalent concentrations (µeq/m 3 ). The ratio of CA/AN = 1 represents neutrality, CA/AN > 1 represents alkalinity and CA/AN < 1 represents acidity. To better evaluate the influence of AD, the non-sea salt concentration was used to calculate the CA and AN.
Atmosphere 2022, 13, 1210 As shown in Figure 4a, the CA/AN ratios in D1 (0.43) and D2 (1.01) were higher than that in NAD (0.37 ± 0.19). This result indicates that the acidity of aerosols will be buffered due to the occurrence of AD. The difference between D1 and D2 can be explained by the transport pathway of the air mass. In D1, the air mass passed over Northern China and mixed with anthropogenic pollutants. The high SO 4 2− and NO 3 − emissions in Northern China resulted in the aerosol still being acidic. However, in D2, AD aerosols mainly contained natural AD particles, with less SO 4 2− and NO 3 − and higher Ca 2+ in the Gobi Desert, resulting in the acidity of aerosols tending to be neutralised. The acidic components in aerosols can be neutralised by alkaline components such as NH 3 and CaCO 3 . Therefore, the neutralisation factor (NF) of NH 4 + and Ca 2+ was calculated to evaluate the contribution of alkaline components. The equation is as follows [52]: influence of AD, the non-sea salt concentration was used to calculate the CA and AN.
As shown in Figure 4a, the CA/AN ratios in D1 (0.43) and D2 (1.01) were higher than that in NAD (0. 37 ± 0.19). This result indicates that the acidity of aerosols will be buffered due to the occurrence of AD. The difference between D1 and D2 can be explained by the transport pathway of the air mass. In D1, the air mass passed over Northern China and mixed with anthropogenic pollutants. The high SO4 2− and NO3 − emissions in Northern China resulted in the aerosol still being acidic. However, in D2, AD aerosols mainly contained natural AD particles, with less SO4 2− and NO3 − and higher Ca 2+ in the Gobi Desert, resulting in the acidity of aerosols tending to be neutralised. The acidic components in aerosols can be neutralised by alkaline components such as NH3 and CaCO3. Therefore, the neutralisation factor (NF) of NH4 + and Ca 2+ was calculated to evaluate the contribution of alkaline components. The equation is as follows [52]: As shown in Figure 4b, NH4 + (NFNH4+ = 65%) was the predominant neutralisation factor in the NAD period. During D1 and D2, nss-Ca 2+ played a more important role due to the enhanced input of AD particles (D1: NFnss-Ca2+ = 62%, D2: NFnss-Ca2+ = 66%).  As shown in Figure 4b, NH 4 + (NF NH4+ = 65%) was the predominant neutralisation factor in the NAD period. During D1 and D2, nss-Ca 2+ played a more important role due to the enhanced input of AD particles (D1: NF nss-Ca2+ = 62%, D2: NF nss-Ca2+ = 66%).

Long-Term Observation of AD
Due to the uncertainty and contingency of the sampling, it is not comprehensive to evaluate the effects of different types of ADs on aerosols only through two AD events. Therefore, in this study, we refer to the data previously observed at KUWAMS from 2010 to 2015 to systematically illustrate the effect of AD [28]. According to a previous study, Japan suffers from different atmospheric pressure patterns in the cold season and warm season: in the cold season, because of the strong Siberian High, air pollutants from the Asian continent tend to be transported to Japan; however, in the warm season, Japan is affected by air masses from the Pacific, since the Pacific High is strong in the summer [53]. To reduce the impact of different pressure patterns on the direction of the monsoon, we focus on the data in the cold season, also called the East Asian monsoon period, to illustrate. From 2010 to 2015, there were a total of 27 AD events observed by the Toyama lidar station during the East Asian monsoon period (Table S2 and Figure S1). As shown in Figure 5, ADs were generated from the Gobi Desert in Mongolia and North-eastern China and then transported at different altitudes, matching the observations of the two ADs in March 2021. Based on the lidar observation of the dust extinction coefficient and the spherical particle extinction in each AD event, we divided all ADs into two types: ADN (aerosols were mainly mineral dust) and ADP (aerosols were a mixture of spherical particles). There focus on the data in the cold season, also called the East Asian monsoon period, to illustrate. From 2010 to 2015, there were a total of 27 AD events observed by the Toyama lidar station during the East Asian monsoon period (Table S2 and Figure S1). As shown in Figure 5, ADs were generated from the Gobi Desert in Mongolia and North-eastern China and then transported at different altitudes, matching the observations of the two ADs in March 2021. Based on the lidar observation of the dust extinction coefficient and the spherical particle extinction in each AD event, we divided all ADs into two types: ADN (aerosols were mainly mineral dust) and ADP (aerosols were a mixture of spherical particles). There were 14 ADNs and 13 ADPs from 2010 to 2015. Two ADs observed in March 2021 might be classified as ADPs (AD1) and ADN (AD2).

General Description of ADs
As shown in Figure S4a, during AD events, the TSP concentration increased significantly (p < 0.001), and the highest TSP concentration appeared in ADN (40.1 ± 19.7 μg/m 3 ), followed by ADP (33.9 ± 16.7 μg/m 3 ) and NAD (19.9 ± 10.5 μg/m 3 ), which matched our discussion on ADs in 2021. This result indicates that AD particles tended to be deposited at KUWAMS if the AD plumes travelled at heights above 3 km, although TSP is thought to be quickly removed when it is passing the ocean region [54]. The WSII concentration variations were significantly correlated with TSP (p < 0.001), with WSII concentrations in ADN (13.5 ± 5.31 μg/m 3 ) being higher than those in ADP (10.8 ± 4.21 μg/m 3 ) and NAD (9.65 ± 3.31 μg/m 3 ) ( Figure S4b).
The increase in the nss-Ca 2+ concentration was the result of the input of AD particles, as discussed previously. Mg 2+ may come from both ocean and continental sources [55]. Through the linear fit regression between Mg 2+ and ss-Na + and nss-Ca 2+ , Mg 2+ had a better correlation with ss-Na + in the NAD period (R 2 = 0.673) and correlated with nss-Ca 2+ in the ADN period (R 2 = 0.554) ( Figure S5). This result indicates that Mg 2+ always tends to originate from ocean sources, since KUWAMS is located near the Japanese Sea; however, when 16% 43% 41% Figure 5. Cluster mean backward trajectory analysis during AD events from 2010 to 2015, and the different colours represent different transport pathways.

General Description of ADs
As shown in Figure S4a, during AD events, the TSP concentration increased significantly (p < 0.001), and the highest TSP concentration appeared in ADN (40.1 ± 19.7 µg/m 3 ), followed by ADP (33.9 ± 16.7 µg/m 3 ) and NAD (19.9 ± 10.5 µg/m 3 ), which matched our discussion on ADs in 2021. This result indicates that AD particles tended to be deposited at KUWAMS if the AD plumes travelled at heights above 3 km, although TSP is thought to be quickly removed when it is passing the ocean region [54]. The WSII concentration variations were significantly correlated with TSP (p < 0.001), with WSII concentrations in ADN (13.5 ± 5.31 µg/m 3 ) being higher than those in ADP (10.8 ± 4.21 µg/m 3 ) and NAD (9.65 ± 3.31 µg/m 3 ) ( Figure S4b).
The increase in the nss-Ca 2+ concentration was the result of the input of AD particles, as discussed previously. Mg 2+ may come from both ocean and continental sources [55]. Through the linear fit regression between Mg 2+ and ss-Na + and nss-Ca 2+ , Mg 2+ had a better correlation with ss-Na + in the NAD period (R 2 = 0.673) and correlated with nss-Ca 2+ in the ADN period (R 2 = 0.554) ( Figure S5). This result indicates that Mg 2+ always tends to originate from ocean sources, since KUWAMS is located near the Japanese Sea; however, when AD events occur, especially for AD originating from the desert region and transported at a higher altitude, the input of AD will cause an increase in Mg 2+ concentration. There are two sources of K + : anthropogenic sources (biomass burning) and natural sources (soil dust particles) [56,57]. Through linear fit regression, K + was better correlated with nss-SO 4 2− (ADN: R 2 = 0.757, ADP: R 2 = 0.717, NAD: R 2 = 0.414) than with nss-Ca 2+ (ADN: R 2 = 0.443, ADP: R 2 = 0.625, NAD: R 2 = 0.399) ( Figure S6a,b). However, since there is large anthropogenic emissions in Mongolia, to better evaluate the level of anthropogenic emissions, we further divided nss-SO 4 2− into the crustal fraction of sulphate (Cr-SO 4 2− ) and anthropogenic fraction of sulphate (an-SO 4 2− ) (Equations (12) and (13)) [58].
The observations in March 2021 showed that the acidity of aerosols was buffered due to AD. This result can also be found in long-term observations, as the CA/AN ratio was higher in ADP (0.62 ± 0.14) and ADN (0.62 ± 0.23) than in NAD (0.42 ± 0.19) ( Figure S4c). The results of the nonparametric tests also showed that CA/AN had significantly different distributions in different categories, with a p-value less than 0.001. For the neutralisation factor, NH 4 + contributed more than 80% to the NAD period. During AD, the effect of Ca 2+ on neutralisation increased, with Ca 2+ contributing more than 50% to the neutralising acid during ADN ( Figure S4d). The results of long-term observation further proved that AD would neutralise aerosol acidity due to the increased effect of Ca 2+ .

Cl Depletion
From the observation in March 2021, Cl depletion increased, along with the AD events occurring. Through long-term observation, Cl depletion was also enhanced simultaneously with AD, with the average Cl depletion showing the highest ADP (62% ± 19%), followed by NAD (52% ± 30%) and ADN (47% ± 35%), although the nonparametric test showed no significant difference (p = 0.493) ( Figure S7b).
Previous studies have demonstrated that Cl − in sea salt aerosols can be removed by reactions with acidic components, such as H 2 SO 4 and HNO 3 . These processes can produce reactive Cl compounds to further affect the oxidation capacity of aerosols [59].

HNO 3 (g) + NaCl (s) = HCl (g) + NaNO 3 (s)
H 2 SO 4 (g) + 2NaCl (s) = 2HCl (g) +Na 2 SO 4 (s) (15) where g represents the gas phases, and s represents the solid phase.  [60]. Additionally, a previous study implied that oxidation of sulphate can proceed rapidly only when the pH is maintained above 6-6.5 [61]. However, in this study, aerosols tend to be acidic. Therefore, the oxidation of SO 2 may be suppressed to further reduce the reaction between sea salt and sulphate.

Secondary Formation of Sulphate and Nitrate
In this study, we used the following equation to calculate the secondary formation of nss-SO4 2− , NO3 − and NH4 + during ADP and ADN [62].

Secondary Formation of Sulphate and Nitrate
In this study, we used the following equation to calculate the secondary formation of nss-SO 4 2− , NO 3 − and NH 4 + during ADP and ADN [62]. (16) where dA is the secondarily formed nss-SO 4 2− , NO 3 − and NH 4 + in sample i; C i A is the concentration of nss-SO 4

Secondary Formation of Sulphate and Nitrate
In this study, we used the following equation to calculate the secondary formation of nss-SO4 2− , NO3 − and NH4 + during ADP and ADN [62].
where dA is the secondarily formed nss-SO4 2− , NO3 − and NH4 + in sample i; C is the concentration of nss-SO4 2− , NO3 − and NH4 + in sample i; is the average concentration of nss-SO4 2− , NO3 − and NH4 + in NAD (C = 4.32, C = 0.74, C = 0.65); C is the average concentration of TSP in NAD (C =19.76) and C is the concentration of TSP in sample i. As a result, there was no significant amount of nss-SO4 2− and NH4 + that were secondarily formed during ADP and ADN, while the secondary formation of NO3 − was found in ADP (dNO3 − = 0.18 ± 1.02) and ADN (dNO3 − = 0.02 ± 0.88) (Figure 7).  SO 4 2− is mainly generated by oxidation reactions with SO 2 gas. Some SO 4 2− exists in sea salt particles in the form of Na 2 SO 4 , while CaSO 4 exists in soil dust particles [32,63]. Through the correlation analysis, nss-SO 4 2− was extremely significantly positively correlated with nss-Ca 2+ (R 2 = 0.625, p < 0.001) and NO 3 − (R 2 = 0.433, p < 0.001) but had a significantly negative correlation with the change in ss-Na + (R 2 = −0.612, p < 0.001). As shown in Figure S8, nss-SO 4 2− was highly correlated with NO 3 − in ADP (R 2 = 0.428) and better correlated with nss-Ca 2+ in ADN (R 2 = 0.536). These results suggest that the source of nss-SO 4 2− was highly affected by the transport pathway. In ADP, air masses mixed with atmospheric pollutants were generated from Northern China. Due to the high consumption of coal for heating in the cold season, the SO 4 2− concentration was high [64]. When AD aerosols passed over North-eastern China, nss-SO 4 2− mixed with AD particles to further increase the nss-SO 4 2− concentration downwind. However, in ADN, AD particles rarely mixed with other aerosols; therefore, the soil dust particles were the main components, such as CaSO 4 . NH 4 + may originate from the reaction of NH 3 gas with acidic gases, such as H 2 SO 4 , HNO 3 and HCl [65]. Therefore, NH 4 + always acts as a neutralisation factor, existing in the form of (NH 4 ) 2 SO 4 and NH 4 NO 3 . NR was calculated based on the equivalent concentrations of SO 4 2− , NO 3 − and NH 4 + and used to evaluate the neutralisation level of NH4 + .
Neutralisation ratio (NR) As shown in Figure S9, both in the 2021 observations (NAD: NR = 0.20 ± 0.09, D1: NR = 0.10 and D2: NR = 0.16) and long-term observations from 2010 to 2015 (NAD: NR = 0.30 ± 0.13, ADP: NR = 0.30 ± 0.11 and ADN: NR = 0.20 ± 0.11), the NR values were all less than 1, suggesting that sulphate and nitrate cannot be fully neutralised by NH 4 + . This result illustrates that the atmospheric environment was ammonium-deficient. We used excess NH 4 + concentration to evaluate the degree of ammonium-deficient environment by the following equation [66,67]. As a result, the excess NH 4 + concentration in all aerosol outflows was less than 0 ( Figure S10)" Previous studies have shown that, under NH 4 + -rich conditions, NO 3 − is formed by gas-phase reactions, while NO 3 − is formed by heterogeneous reactions of HNO 3 in liquid phase under NH 4 + -poor conditions [68]: cantly negative correlation with the change in ss-Na + (R 2 = −0.612, p < 0.001). As shown in Figure S8, nss-SO4 2− was highly correlated with NO3 − in ADP (R 2 = 0.428) and better correlated with nss-Ca 2+ in ADN (R 2 = 0.536). These results suggest that the source of nss-SO4 2− was highly affected by the transport pathway. In ADP, air masses mixed with atmospheric pollutants were generated from Northern China. Due to the high consumption of coal for heating in the cold season, the SO4 2− concentration was high [64]. When AD aerosols passed over North-eastern China, nss-SO4 2− mixed with AD particles to further increase the nss-SO4 2− concentration downwind. However, in ADN, AD particles rarely mixed with other aerosols; therefore, the soil dust particles were the main components, such as CaSO4. NH4 + may originate from the reaction of NH3 gas with acidic gases, such as H2SO4, HNO3 and HCl [65]. Therefore, NH4 + always acts as a neutralisation factor, existing in the form of (NH4)2SO4 and NH4NO3. NR was calculated based on the equivalent concentrations of SO4 2− , NO3 − and NH4 + and used to evaluate the neutralisation level of NH4 + .
Neutralisation ratio (NR) = NH (neq/m ) SO (neq/m ) + NO (neq/m ) As shown in Figure S9, both in the 2021 observations (NAD: NR = 0.20 ± 0.09, D1: NR = 0.10 and D2: NR = 0.16) and long-term observations from 2010 to 2015 (NAD:NR = 0.30 ± 0.13, ADP:NR = 0.30 ± 0.11 and ADN:NR = 0.20 ± 0.11), the NR values were all less than 1, suggesting that sulphate and nitrate cannot be fully neutralised by NH4 + . This result illustrates that the atmospheric environment was ammonium-deficient. We used excess NH4 + concentration to evaluate the degree of ammonium-deficient environment by the following equation [66,67]. As a result, the excess NH4 + concentration in all aerosol outflows was less than 0 ( Figure S10)" Previous studies have shown that, under NH4 + -rich conditions, NO3 − is formed by gas-phase reactions, while NO3 − is formed by heterogeneous reactions of HNO3 in liquid phase under NH4 + -poor conditions [68]: N2O5 (aq) + H2O (aq) ⇄ 2NO3 − (aq) + 2H + (aq) (20) As shown in Figure 8, the NO3 − concentration increased while the excess NH4 + concentration decreased, suggesting that the NO3 − concentration in this study cannot be explained by the homogeneous gas-phase reaction between NH3 and nitric acid. Nitrate can also be formed by heterogeneous hydrolysis of N2O5 on the moist aerosol surface during night-time, which depends on the concentration of N2O5 and water in the aerosol [69]. According to previous laboratory and field studies, the heterogeneous hydrolysis of N2O5 can be promoted by the enhanced hygroscopicity, surface area and acidity of the pre-existing particles [70][71][72]. The N2O5 uptake coefficients will increase if aerosols tend to be acidic [73]. Additionally, water plays an important role in the heterogeneous hydrolysis of N2O5, since water can promote the uptake of N2O5 on the aerosol surface [72]. However, the increase in water can decrease the aerosol acidity to further reduce the N2O5 uptake. In AD, the increased AD particles provided a surface area for heterogeneous reactions between the gases HNO3, NO2 and SO2 and particulate (NH4)2SO4. The acidic aerosols can promote the adsorption of N2O5 on the surface of particles as well. Previous studies have shown that the hygroscopicity of chemical species in dust particles depends on their transport pathways [74]. In ADP, hygroscopic species were formed by the interaction of AD particles with anthropogenic pollutants, such as Ca(NO3)2 and (NH4)2SO4, and during long-range transportation, with the aging of dust particles, the hygroscopes increased.
Through the correlation analysis, nss-SO4 2− was extremely significantly positively correlated with nss-Ca 2+ (R 2 = 0.625, p < 0.001) and NO3 − (R 2 = 0.433, p < 0.001) but had a significantly negative correlation with the change in ss-Na + (R 2 = −0.612, p < 0.001). As shown in Figure S8, nss-SO4 2− was highly correlated with NO3 − in ADP (R 2 = 0.428) and better correlated with nss-Ca 2+ in ADN (R 2 = 0.536). These results suggest that the source of nss-SO4 2− was highly affected by the transport pathway. In ADP, air masses mixed with atmospheric pollutants were generated from Northern China. Due to the high consumption of coal for heating in the cold season, the SO4 2− concentration was high [64]. When AD aerosols passed over North-eastern China, nss-SO4 2− mixed with AD particles to further increase the nss-SO4 2− concentration downwind. However, in ADN, AD particles rarely mixed with other aerosols; therefore, the soil dust particles were the main components, such as CaSO4. NH4 + may originate from the reaction of NH3 gas with acidic gases, such as H2SO4, HNO3 and HCl [65]. Therefore, NH4 + always acts as a neutralisation factor, existing in the form of (NH4)2SO4 and NH4NO3. NR was calculated based on the equivalent concentrations of SO4 2− , NO3 − and NH4 + and used to evaluate the neutralisation level of NH4 + .
Neutralisation ratio (NR) = NH (neq/m ) SO (neq/m ) + NO (neq/m ) (17) As shown in Figure S9, both in the 2021 observations (NAD: NR = 0.20 ± 0.09, D1: NR = 0.10 and D2: NR = 0.16) and long-term observations from 2010 to 2015 (NAD:NR = 0.30 ± 0.13, ADP:NR = 0.30 ± 0.11 and ADN:NR = 0.20 ± 0.11), the NR values were all less than 1, suggesting that sulphate and nitrate cannot be fully neutralised by NH4 + . This result illustrates that the atmospheric environment was ammonium-deficient. We used excess NH4 + concentration to evaluate the degree of ammonium-deficient environment by the following equation [66,67]. As a result, the excess NH4 + concentration in all aerosol outflows was less than 0 ( Figure S10)" Previous studies have shown that, under NH4 + -rich conditions, NO3 − is formed by gas-phase reactions, while NO3 − is formed by heterogeneous reactions of HNO3 in liquid phase under NH4 + -poor conditions [68]: N2O5 (aq) + H2O (aq) ⇄ 2NO3 − (aq) + 2H + (aq) (20) As shown in Figure 8, the NO3 − concentration increased while the excess NH4 + concentration decreased, suggesting that the NO3 − concentration in this study cannot be explained by the homogeneous gas-phase reaction between NH3 and nitric acid. Nitrate can also be formed by heterogeneous hydrolysis of N2O5 on the moist aerosol surface during night-time, which depends on the concentration of N2O5 and water in the aerosol [69]. According to previous laboratory and field studies, the heterogeneous hydrolysis of N2O5 can be promoted by the enhanced hygroscopicity, surface area and acidity of the pre-existing particles [70][71][72]. The N2O5 uptake coefficients will increase if aerosols tend to be acidic [73]. Additionally, water plays an important role in the heterogeneous hydrolysis of N2O5, since water can promote the uptake of N2O5 on the aerosol surface [72]. However, the increase in water can decrease the aerosol acidity to further reduce the N2O5 uptake. In AD, the increased AD particles provided a surface area for heterogeneous reactions between the gases HNO3, NO2 and SO2 and particulate (NH4)2SO4. The acidic aerosols can promote the adsorption of N2O5 on the surface of particles as well. Previous studies have shown that the hygroscopicity of chemical species in dust particles depends on their transport pathways [74]. In ADP, hygroscopic species were formed by the interaction of AD particles with anthropogenic pollutants, such as Ca(NO3)2 and (NH4)2SO4, and during long-range transportation, with the aging of dust particles, the hygroscopes increased.
As shown in Figure 8, the NO 3 − concentration increased while the excess NH 4 + concentration decreased, suggesting that the NO 3 − concentration in this study cannot be explained by the homogeneous gas-phase reaction between NH 3 and nitric acid. Nitrate can also be formed by heterogeneous hydrolysis of N 2 O 5 on the moist aerosol surface during night-time, which depends on the concentration of N 2 O 5 and water in the aerosol [69]. According to previous laboratory and field studies, the heterogeneous hydrolysis of N 2 O 5 can be promoted by the enhanced hygroscopicity, surface area and acidity of the preexisting particles [70][71][72]. The N 2 O 5 uptake coefficients will increase if aerosols tend to be acidic [73]. Additionally, water plays an important role in the heterogeneous hydrolysis of N 2 O 5 , since water can promote the uptake of N 2 O 5 on the aerosol surface [72]. However, the increase in water can decrease the aerosol acidity to further reduce the N 2 O 5 uptake. In AD, the increased AD particles provided a surface area for heterogeneous reactions between the gases HNO 3 , NO 2 and SO 2 and particulate (NH 4 ) 2 SO 4 . The acidic aerosols can promote the adsorption of N 2 O 5 on the surface of particles as well. Previous studies have shown that the hygroscopicity of chemical species in dust particles depends on their transport pathways [74]. In ADP, hygroscopic species were formed by the interaction of AD particles with anthropogenic pollutants, such as Ca(NO 3 ) 2 and (NH 4 ) 2 SO 4 , and during long-range transportation, with the aging of dust particles, the hygroscopes increased. Therefore, we indicate that the secondary generation of NO 3 − observed in ADP is due to the heterogeneous reaction of HNO 3 .
Atmosphere 2022, 13, x FOR PEER REVIEW 12 of 16 Therefore, we indicate that the secondary generation of NO3 − observed in ADP is due to the heterogeneous reaction of HNO3.

Conclusions
In this study, the effect of Asian dust (AD) on water-soluble inorganic ions (WSIIs) in total suspended particles (TSPs) was investigated at KUWAMS, a remote background site in Japan, from 27 February to 3 April 2021. To provide a comprehensive illustration, data from long-term observations of WSIIs at KUWAMS from 2010 to 2015 were referred During the East Asian monsoon period, the air mass mainly flowed out of the Asian continent. Based on the lidar observations and backwards trajectory analysis, the observed ADs were classified into two types: ADs mainly containing natural AD particles and ADs mixed with other anthropogenic pollutants. The TSP and WSII concentrations increased due to AD outbreaks, especially when AD was transported at a higher altitude with faster speed. The AD marker nss-Ca 2+ had a concentration variation similar to that of TSP. AD buffered the acidity of aerosols with the increase in Ca 2+ as a neutralisation factor. A higher Cl depletion was observed when the AD particles were mainly natural particles, which was caused by the heterogeneous reaction of HNO3 with sea salt. SO4 2− from soil dust existed in the form of CaSO4 in the natural AD particles, while it came from combustion sources when AD passed through polluted areas. When AD particles were transported with anthropogenic pollutants, the secondary formation of NO3 − was enhanced, indicating the effect of AD on secondary aerosols, which was attributed to the heterogeneous hydrolysis of N2O5. Based on long-term observations, this study systematically analysed the chemical composition changes of water-soluble inorganic ions during Asian dust, which can more comprehensively describe the impact of Asian dust on the changes of aerosol in downwind areas and reduce the specificity between Asian dusts.

Supplementary Materials:
The following supporting information can be downloaded at www.mdpi.com/xxx/s1: Figure S1. Lidar observation of the attenuated backscatter coefficient (532 nm) during AD events in 2021 and from 2010 to 2015. Each AD event is marked in red blocks, and the sampling period is marked with a black dotted line. Figure S2 Figure S5. Linear fit regression of Mg 2+ with nss-Ca 2+ (a) and ss-Na + (b) during long-term observation from 2010 to 2015. Figure S6. Linear fit regression of K + with nss-Ca 2+ (a), nss-SO4 2− (b), cr-SO4 2− (c)

Conclusions
In this study, the effect of Asian dust (AD) on water-soluble inorganic ions (WSIIs) in total suspended particles (TSPs) was investigated at KUWAMS, a remote background site in Japan, from 27 February to 3 April 2021. To provide a comprehensive illustration, data from long-term observations of WSIIs at KUWAMS from 2010 to 2015 were referred During the East Asian monsoon period, the air mass mainly flowed out of the Asian continent. Based on the lidar observations and backwards trajectory analysis, the observed ADs were classified into two types: ADs mainly containing natural AD particles and ADs mixed with other anthropogenic pollutants. The TSP and WSII concentrations increased due to AD outbreaks, especially when AD was transported at a higher altitude with faster speed. The AD marker nss-Ca 2+ had a concentration variation similar to that of TSP. AD buffered the acidity of aerosols with the increase in Ca 2+ as a neutralisation factor. A higher Cl depletion was observed when the AD particles were mainly natural particles, which was caused by the heterogeneous reaction of HNO 3 with sea salt. SO 4 2− from soil dust existed in the form of CaSO 4 in the natural AD particles, while it came from combustion sources when AD passed through polluted areas. When AD particles were transported with anthropogenic pollutants, the secondary formation of NO 3 − was enhanced, indicating the effect of AD on secondary aerosols, which was attributed to the heterogeneous hydrolysis of N 2 O 5 . Based on long-term observations, this study systematically analysed the chemical composition changes of water-soluble inorganic ions during Asian dust, which can more comprehensively describe the impact of Asian dust on the changes of aerosol in downwind areas and reduce the specificity between Asian dusts.
Supplementary Materials: The following supporting information can be downloaded at https:// www.mdpi.com/article/10.3390/atmos13081210/s1: Figure S1: Lidar observation of the attenuated backscatter coefficient (532 nm) during AD events in 2021 and from 2010 to 2015. Each AD event is marked in red blocks, and the sampling period is marked with a black dotted line. Figure S2 Figure S9: Neutralisation ratio in different aerosol outflows during sampling period in 2021 (a) and long-term observation from 2010 to 2015 (b). Figure S10: Concentration of (excess NH 4 + ) with (NO 3 − ) during long-term observation from 2010 to 2015. Table S1: Concentration of water-soluble inorganic ions (WSIIs) and non-sea salt WSIIs during the sampling period in 2021. Table S2: Time periods and type of AD from 2010 to 2015.
Author Contributions: N.T. designed the sampling work, the WSIIs analyses for particulate matters; N.T., L.Z., H.Z., X.Z. and Y.W. collected samples and conducted monitoring work; P.B. analysed the WSIIs and performed a statistical analysis and S.N., B.C. and N.T. gave recommendations to the paper. All authors have read and agreed to the published version of the manuscript.