Boron Isotopic Analysis of Representative Atmospheric Aerosols Derived From Long-Range Transported/Local Emission on an Islet Offshore NE Taiwan

To identify pollutant origins and their potential sources from either long-range transported (LRT) or local emissions, the Boron (B) concentrations and the B isotopic compositions (δ11B) in the water-soluble fraction of aerosols were determined. These aerosols were collected from an islet offshore in northeastern Taiwan, Peng Chia Yu (PCY), a non-residential islet, from January 1998 to March 2000. This islet may be influenced by various pollutants or chemical transport during monsoon periods, but suffers minimal human perturbation locally. The B in the specimens falls to 0.3–1.63 ng m−3 during the SW monsoon seasons when compared to 0.46–2.56 ng m−3 in the NE monsoon. However, the δ11B results show no clear variations in both monsoon seasons (10.7–24.3‰), regardless of differences in air mass origin. A two end-member mixing scenario is proposed to explain our observations using the obtained δ11B and 1/[B] results. The ocean endmember is characterized by high 1/[B] and high δ11B; while other endmember is from continental endmember or anthropogenic contributions that are characterized by of 1/[B] and lowδ11B (δ11B < 10‰). Based on these chemical and isotopic results, we found aerosol emissions in northern Taiwan are characterized as low B with nearly constant δ11B, when compared with long-range transported continental endmembers derived from nearby regions. This study provides preliminary B and δ11B levels in aerosols derived from LRT/local emissions and discusses potential monsoonal effects on aerosols offshore of NE Taiwan.


INTRODUCTION
Air pollution in Taiwan derives primarily from local emissions and long-range transported (LRT; Lin, 2001;Hsu et al., 2006). Identifying the sources of pollutants is the first step in developing greater air quality control program in the environment. The emissions from local sources are generated mainly by the domestic combustion of fossil fuels. In Southeast Asia, aerosol transportation from China has rapidly increased over the last few decades, because of a marked increase in coal consumption (Fogg and Duce, 1985;Park and Schlesinger, 2002;Zhao and Liu, 2010). The sources of pollutants varies seasonally due to differing weather patterns and wind directions (Cheng et al., 2010;Cheng et al., 2011). As a nonresidential islet, Peng Chia Yu, offers limited environmental perturbations that make it is an excellent location for monitoring through airborne in pollutants southeastern Asia.
Water-soluble ions are major components (up to 60-70%) of the atmospheric aerosols, which include sulfate, nitrate, ammonium and alkali cations produced from various anthropogenic processes (Ali-Mohamed, 1991). The chemical and isotopic composition of water-soluble ions in aerosols provide valuable information for a better understanding of their chemical characteristics, pollutant sources, and transportation mechanisms (Hoffman et al., 1977;Zhang et al., 2011). Stable isotopes are powerful tools for investigating sources of atmospheric aerosols. For example, Pb isotopic fingerprints are tracers used for anthropogenic sources, but not for natural sources. B is a trace constituent in the environment that can be divided into a natural or an anthropogenic source in the atmosphere. Significant mass differences between 10 B and 11 B demonstrates and leads to larger δ 11 B variations from ∼60‰ rainwater and aerosols (Park and Schlesinger, 2002;Schlesinger and Vengosh, 2016). Its variations in aerosols are much larger than other natural isotopic systems that make it a powerful source for a future study on identification study.
B exists in as both particulate and gaseous forms in the atmosphere with total concentrations ranging from a few ng m −3 to more than 100 ng m −3 (Fogg and Duce, 1985;Miyata et al., 2000). Early works found B sources in the atmosphere were derived from multiple origins, including the evaporation of boric acid from seawater and sea spray by large fluxes (1.44 Tg B/year), erosion from terrestrial materials (0.017 Tg/year), emissions from the combustion of biomass (0.26-0.43 Tg B/year), and coals (0.20 Tg B/year) (Park and Schlesinger, 2002). B has two stable isotopes, 10 B and 11 B, with relative abundance 19.9 and 80.1%, respectively. In recent studies, B and B isotopes were used as tracers for rainwater and aerosols to identify the sources and transport routes of air pollutants (Sakata et al., 2010;Zhao and Liu, 2010;Sakata et al., 2013;Roux et al., 2017). B in the atmosphere was further divided into a natural or an anthropogenic source (Sakata et al., 2013;Roux et al., 2017;Sakata et al., 2017). Sea salt and the evaporation of boric acid from seawater are primary natural emissions characterized with high levels of δ 11 B close to seawater that of (δ 11 B 39.6‰). Most terrestrial materials (e.g., rocks, evaporates, crude oil, and water) had δ 11 B values that ranged from −30 to 45‰ and were mostly greater than 0‰ (Williams and Hervig, 2004). Most of the anthropogenic contributions of B in the atmosphere were derived from the combustion of coals and characterized by much lower δ 11 B values (δ 11 B −70 to 8‰) (Park and Schlesinger, 2002;Zhao and Liu, 2010). As mentioned above, sea salt aerosols and the evaporation of boric acid from seawater are primary natural emissions of B and the anthropogenic contributions of B to the atmosphere was derived from the combustion of coal. The B content vapor was condensed to aerosol and transmitted to PCY by NE or summer monsoon.
There has been insufficient research on B isotope studies in rainwater and aerosols because of limited amounts of B in the atmosphere. To date, B isotopes have only been used to trace the sources of particle emissions, including Shizuoka, Matsuura and Nakanoto in Japan (Sakata et al., 2013;Sakata et al., 2017) and in Lorraine and Paris in France (Chetelat et al., 2009;Roux et al., 2017).
In this study, we extracted and analyzed the elemental concentrations and δ 11 B in water-soluble fractions of aerosol specimens in PCY. This represents a preliminary attempt to analyze aerosols with B isotopic compositions for improving air-quality monitoring in Taiwan/Asia. Our objectives are to provide δ 11 B data in PCY, to study representative LRT/local emissions in NE Taiwan, and to identify any potential sources using B and δ 11 B data.

MATERIALS AND METHODS
Aerosol samples were collected at a weather station at PCY under the Taiwan Central Weather Bureau (25°37′46″ N, 122°04′17″ E; Figure 1), from January 1998 to March 2000. The Meteorologic data for PCY is influenced alternatively by the Pacific subtropical high-pressure system in the summer and the NE monsoon in the winter. Aerosol sources changes caused by wind patterns alternate seasonally. To study source variations with seasonal monsoons, we choose samples from the summer monsoon (June to August) and the northeastern monsoon periods (September to March). Additionally, we selected samples from typhoon or dust storm events for analysis. A total of 29 samples were selected for this study.
The PM2.5 fractions of aerosol were collected with the flow rate of 1.05 m 3 min −1 on acetate membrane filters (Whatman Grade 41 Quantitative Filter Pape, Thermo Fisher, Scientific ™ , 20 × 25 cm 2 ). The sampler was placed on the rooftop of PCY weather station, which include an air sampler coupled with suction motor (GBM-2000H, General Metal Works Inc.), stainless steel grid (GFH-2100, General Metal Works Inc.), and continuous flow system (G105, General Metal Works Inc.). Those filters were kept in desiccator in the refrigerator at low temperature to prevent potential contaminations after sampling and put it in the sealed zipper bags to prevent element lost through evaporation.
To extract the water-soluble elements, aerosol filters were cut and placed in a clean Teflon vial, and then treated with 6 ml Milli-Q water for 90 min at room temperature under ultrasonic irradiation. Then the aerosol filters were removed from the Teflon vials after soaking in Milli-Q water at room temperature overnight. These solutions were passed through a 0.45 µm-membrane filter then acidified by adding 0.3 M HNO 3 conditions that were used for the determination of element concentrations.
The concentration of major elements sodium (Na) and sulfur (S) were measured by inductively coupled plasma optical emission spectrometry (ICP-OES, Thermo-Fisher, iCAP) at the Earth Dynamic System Research Center, National Cheng Kung University (EDSRC, NCKU). The concentration of B and lead (Pb) were measured by Inductively Coupled Plasma Quadrupole Mass Spectrometry (ICP-Q-MS Thermo-Fisher, iCAP-Q) at EDSRC, NCKU. The analytical precision of the dissolved elements is better than 2%. B in the sample solution was purified by using a microsublimation technique to separate B from all other matrices in natural water or low-B marine biogenic carbonates (Wang et al., 2010). This technique shows additional advantages for low procedural blanks (∼8 pg) and is easy to operate, i.e., it is an ideal way to efficiently purify low-B-content samples, such as those in aerosols and rainwater. 50 µl of acidified watersoluble ion samples were deposited at the center of the cap of a 5 ml conical vial. A minimum of 50 ng of B is needed for analysis; as a result, some samples with low B content were combined in multiple vials. These vials were tightly closed in an upside-down position on a homemade thermostatic hotplate with a Teflon coated graphite rack. They were then heated for 18 h under 98 ± 0.1°C. Each sample was diluted with 0.1 N HNO 3 to a final volume of 1 ml for the MC-ICP-MS measurements.
All isotopic measurements were carried out using a Multi-Collector Inductively Coupled Plasma Mass Spectrometry (MC-ICP-MS, Thermo Fisher Scientific, Neptune) at EDSRC, NCKU. The B isotopic compositions are presented in delta notation relative to a Standard Reference Material (SRM) 951 boric acid standards, prepared by National Institute of Standards and Technology (NIST) and using the equation below: The reproducibility (2SE) of δ 11 B for NIST SRM 951 in this laboratory was approximately of 0.4‰. The evaporation of seawater is an important source of B in aerosols, by calculating the Enrichment Factors (EFs) using B/Na, the values can then be used to evaluate the total marine contribution in the data (Chabas and Lefevre, 2000). The EF values were calculated using the equation below: where B and Na are the concentrations of water-soluble B and Na in samples and seawater, respectively (B/Na 4.2 × 10 -4 ; Li, 1991). The air trajectory analysis allows for the evaluation of the pathways of air mass during different seasons. The Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model was used.  The S and Pb concentrations in aerosols ranged from 0.04 to 3.3 μg m −3 and from 0.04 to 18.81 ng m −3 , respectively. By mapping the two chemical values with the seasonal variation chart ( Figure 2B), three samples (A18, A19 and A26) had relatively high Pb concentrations (14.7-18.8 ng m −3 ) that corresponded to NE monsoon the seasons. There were no significant changes in S and Pb during typhoons, dust storms and potential volcanic emissions events. Furthermore, the concentrations of S and Pb in aerosols during the NE monsoon period (with the exception of dust and typhoon sample) correlate positively (R 2 0.64, p < 0.001; Figure 3) and demonstrate a similar trend over time.

Seasonal δ 11 B Variations
The δ 11 B in aerosols during the NE monsoon period ranged from 10.7 to 24.0‰ ( Figure 2C), where sample A41, with the heaviest δ 11 B, was characterized by a relatively low B concentration (0.58 ng m −3 ). On the other hand, the sample with the lightest δ 11 B is characterized by the highest B concentration. The δ 11 B in dust storm samples (A50 and A51) were characterized by more negative values, 14.1 and 10.9‰, respectively ( Figures 2C, 4).
The δ 11 B in aerosols range from 13.1 to 24.3‰ ( Figure 2C) in the summer monsoon period and sample A37 had the heaviest δ 11 B. The typhoon sample (A15) and volcanic emissions (A13) were characterized relatively negative δ 11 B ( Figures 2C, 5).

Sea Salt and Seawater Evaporation
The B and the EF seawater showed similar seasonal changes ( Figure 2; Table 1), the EFs values were higher during the NE monsoon than in summer monsoon seasons (0.7-4.8). In this study, dissolved Na was selected as reference element for EF seawater calculation. EFs higher than 1 indicated B in sample relative enrichment compared to seawater. Relatively high EF seawater (1.2-9.7 with the exception of the least sample A41) during the NE monsoons implies that severe contamination was caused by anthropogenic emissions. During summer monsoon, the relatively low EF seawater (0.7-4.8) is less than 5 and indicate a significant influence from seawater. Some samples of EF seawater valued less than 1, which was likely due to B were lost through the evaporation.
B isotopes can be used as tracers for aerosols to identify the contributions and transport routes for natural and anthropogenic sources. Although seawater evaporation is a primary source of B in aerosols, the marine endmember of B isotope in aerosols has yet to be clarified from the literature. The δ 11 B in the seawater and coastal rainwater ranged up to +45‰ (Chetelat et al., 2005). Previous studies utilized the general δ 11 B values between 39.6 and 45‰ to estimate the contributions of each source in the aerosols (Zhao and Liu, 2010;Sakata et al., 2013). In Sakata et al. (2013), they examined the B isotopic compositions in aerosols at Matsuura and Nakanoto in Japan, which showed that the highest δ 11 B value (21.7‰) in their study was different from the seawater values (39.6‰), even when the levels of B in this aerosol sample were mostly derived from seawater. All of the δ 11 B in aerosols published were markedly lower than the seawater values.
The δ 11 B in seawater and coastal rainwater should not refer directly to the δ 11 B for marine endmembers in aerosols. Significant negative isotopic fractionation (Δ particle-liquid ) may occur during the transformation from liquid B to particulate B, due to the δ 11 B values of particulate B were expected to be 0-20‰ lower than gaseous B, and the liquid vapor isotopic fractionation (Δ liquid-vapor ) is positive (Sakata et al., 2017). Roux et al. (2017) examined the evolution of B and B isotopic compositions of atmospheric particulate deposition and total dissolved deposition between 2012 and 2013. The δ 11 B in particulate depositions was much lower than those from the dissolved depositions. As mentioned above, there were significant isotopic fractionation between rainwater and aerosol. The marine endmembers of B isotope within aerosols (particles) should be clarified in a future study.
We suggest the highest δ 11 B in marine aerosols can be potentially represented by the δ 11 B value for marine endmembers in aerosols. Aerosols and rainwater derived from seawater evaporation should display the highest δ 11 B, as seawater has the highest δ 11 B source in the atmosphere. As the δ 11 B and  Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 697844 serves as a common in coastal rainwater increased to 45‰, the fractionation factors associated with the evaporation of seawater (α seawater-vapor 1.025) and the vapor-condensation processes (α liquid-vapor 1.031) are shown giving the first rains at 45‰ (Chetelat et al., 2005;Rose-Koga et al., 2006). To our knowledge, we reported the highest δ 11 B (24.3‰, A37) in aerosols in this study. Therefore, the δ 11 B marine aerosols is 24.3‰ and the rainaerosol isotopic fractionation is 20.7‰ (δ 11 B rain − δ 11 B aerosol ). Sakata et al. (2017) supported this hypothesis because they showed gaseous B was transformed into particulate B through the reaction of condensed B(OH) 3 with chemical constituents on particles to that precipitate borates. The isotopic fractionation was about 20‰ between B(OH) 3 and B(OH) 4 − , with B (OH) 4 − concentrating 10 B (Kakihana et al., 1977). Our results (δ 11 B rain -− δ 11 B aerosol δ 11 B liquid − δ 11 B particle 20.7‰) also supported this hypothesis.

Anthropogenic Emissions
The S and Pb concentrations in aerosols are often indicators of anthropogenic activities. Combining data with B isotopes can provide more concrete evidence for its evolutionary processes. There are two specimens (A18, A26) with the lightest B isotopic values (without the influences of typhoons) and also holds the highest Pb and S concentration levels (Figures 6, 7). The air mass from the LRT during the NE monsoon season was affected by anthropogenic activities in this region ( Figure 8A). This Frontiers in Environmental Science | www.frontiersin.org June 2021 | Volume 9 | Article 697844 indicated that B had common origins or similar transportation pathways as those of S and Pb. The marked Pb concentrations in aerosols were measured in the NE monsoon season and enhanced by LRT from China (Hsu et al., 2005). The dissolved Pb in PCY derived from Taiwan and LRT from Shanghai varied largely from 8.9 to 53.5% and 46.5 to 91.1%, respectively (Cheng et al., 2011). Sakata et al. (2013) measured the atmospheric concentrations and isotopic composition of B and S in aerosol samples collected at two sites along the coast of the Japan Sea. They found that aerosols during the winter when air flows from northern China were dominant and characterized by negative δ 11 B values, which accurately reflected the conditions of δ 11 B for the coal use in this region.
To clarify possible anthropogenic sources of LRT from China in the NE monsoon season, we explored further the study of δ 11 B in aerosols from Zhao and Liu (2010) as measured with precipitation in Guiyang, China was important. The B extracted from the precipitation of Guiyang showed four B sources, which were marine, biomass, and organic matter, along with coal combustion at a regional level. This suggests that the δ 11 B values for aerosol biomass or organic matter and coal combustion were −5 and 8‰, respectively, and can be used to calculate the contribution of different end-members.
The air mass backward trajectory results ( Figure 8A) indicated a predominant contribution by LRT from China during the NE monsoon periods (Cheng et al., 2010). Most B anthropogenic contributions in the atmosphere were derived from the combustion of coals (Park and Schlesinger, 2002;Zhao and Liu, 2010). The results showed a significant correlation (R 2 0.57, p < 0.001, Figure 9) between δ 11 B and 1/ [B]. Figure 9 shows that it    intercepted at approximately 11.6‰, with two lightest δ 11 B values of 10.7 and 12.6‰. These values were rather close to coal combustion in Guiyang, China, 8‰ (Zhao and Liu, 2010). Hence, we assumed that the boron content and isotopic composition in the winter monsoon season composed of a combination of evaporation seawater and coal burning.

Potential Volcanic Emissions
In retrospect, the B concentrations in aerosol during the summer monsoon period show a narrow range (0.3-1.29 ng m −3 ), with the exception of there being the highest B sample ( Figure 5). The highest B sample (A13) was four times greater than the summer average B and collected from July 13, 1998 to July 15, 1998. On July 11, 1998, the eruptions of the Mount Merapi volcano in Indonesia, which lasted for several weeks. The air mass backward trajectory results indicate the volcano ashes and gases were delivered across the equator, the South China Sea, and finally arrived at PCY (Cheng et al., 2010). Cheng et al. (2010) reported the anomalous strontium isotopic ratio ( 87 Sr/ 86 Sr) during this period, which was explained by the addition of the less radiogenic 87 Sr/ 86 Sr substances derived from the Merapi volcano eruption in July 1998. During this period, the anomalous B and δ 11 B in aerosols was likely affected by the Merapi volcano eruption. In the B sourcesink model, the flux of atmospheric B from volcanic eruption was only 0.02 Tg B/year, which was minimal to the overall atmospheric B flux (Park and Schlesinger, 2002). However, due to the high B concentration in the volcanic gases (0.2-104 µg m −3 ; Fogg and Duce, 1985), the atmospheric B of the region will significantly increase in short periods of time during eruption. The anomalous B and δ 11 B can be explained by additions of the substances derived from the Merapi volcano eruption in July 1998. The aerosol specimen analyzed clearly indicates the influence of volcanic gases with the value of δ 11 B 14.4‰.

Dust Storm and Typhoon Events
Two samples (A50 and A51) taken during dust storm events and characterized by high [B] and low δ 11 B (Figures 4, 8C) during the NE monsoon. These samples had relatively low S and Pb levels when compared to anthropogenic sources (Figures 6, 7), which and indicated that B was not derived from coal burning. Dust storms arise when strong winds blow loose sand and dirt from a dry surface. The sandy loess in northeast China had a low δ 11 B (−12.1 to 9.4‰; Lei et al., 2021), which an important source of B aerosols during storm events. There were no significant change in B and δ 11 B during typhoon events. Sample A19 was affected by typhoon ZEB in 1998, which was the strongest typhoon in the western North Pacific that year. The NE monsoon wind will be enhanced for a short time prior to a typhoon making landfall. The unchanged wind direction and atmospheric sources may be the main reason for a lack of significant changes in B and δ 11 B during winter typhoon events.

Rayleigh-like Evolution δ 11 B Fractionation
As no significant differences within B concentration level and air mass backward trajectory results ( Figure 8B) during the summer monsoon, we can assume that the changes to B isotopic ratios were not caused by mixture of marine and high B content endmembers (i.e., biofuel, coal). Apart from the effects of endmember mixing on the atmospheric B isotope, during the sea salt evaporation, condensation, or deposition processes, the isotopic values change and were subsequently affected the B isotope aerosol values. The combination of low B concentration with low δ 11 B was observed early on and explained by the evaporation of gaseous boric acid from the seawater and subsequent integration in rainwater of isotopically fractionated B (Roux et al., 2017). Boron isotopes in rainwater were then modeled by Rayleigh-like a distillation processes (Rose-Koga et al., 2006;Roux et al., 2017):  With R and R 0 as the isotopic ratio of the rainwater and the original atmospheric gaseous B reservoir, respectively. f is the mass fraction of residual gaseous 10 B and α is the boron isotopic fractionation factor between the gaseous reservoir and α liquid-vapor 1.031 (Rose-Koga et al., 2006). The large variations in the 11 B/ 10 B ratio of atmospheric B, either as vapor or as precipitation, are related to the seawater boron reservoir via a Rayleigh-like distillation process (Rose-Koga et al., 2006). As mentioned above, we suggest the fractionation factor between liquid and particle is constant (α liquid-particle 1.021). B isotopes in aerosols should also be affected by a Rayleigh-like distillation process. Using the empirical α value (α liquid-vapor 1.031), the Rayleigh distillation model of marine aerosols and water/vapor exchange processes can explain most of the B data during the summer monsoon period (Figure 7).
The results show a significant correlation (R 2 0.47, p < 0.001; Figure 10A) between δ 11 B and mean temperature for the summer samples. The mean temperature data for the collected day was retrieved from NOAA GSOD (https://www.ncdc.noaa.gov/). Palmer et al. (1987) demonstrated the fractionation factor between the dissolved and adsorbed species (α liquid-particle ) remained constant at 5°C and at 40°C. Rose-Koga et al. (2006) demonstrated the fractionation factors between vapor and seawater (α liquid-vapor ) also indicated no significant differences between 15°C (1.031) and 30°C (1.026). The variation isotope fractionation with temperature during evaporation and condensation can be ignored. Condensation of marine-derived gaseous B would be the dominant process during the summer when the higher temperatures promote the evaporation of B (Roux et al., 2017). This also enhanced the Rayleigh distillation process, strengthened B isotope fractionation, and decreased the B and δ 11 B in the aerosol.
There was no significant correlation (R 2 0.03, p > 0.001; Figure 10B) between δ 11 B and temperature for the NE monsoon (with the exception of dust and typhoon sample). The δ 11 B fractionation in Rayleigh distillation were enhanced in temperature increased, this process should be significant for the NE monsoon sample with large temperature range (12.5-27.3°C; Table 2). From our perspective, the δ 11 B fractionation during phase transformation was significantly lower than two end-member mixing during the NE monsoon months. The B content in coal is much higher than seawater. The B and δ 11 B in aerosol will be strongly affected by coal emission, if their derived from coal emission and seawater evaporation. The mixing of two end-member has played the most important role for δ 11 B during the NE monsoon in the winter.

Estimate Relative Contributions of Each Source in Peng Chia Yu Aerosols
As mentioned above, a two endmember mixing model (seawater evaporation and the coal burning) can be used qualitatively for the B and δ 11 B variations in aerosols during the NE monsoon period, which was supported by the plot of δ 11 B vs. 1/[B] (Figure 9). The contributions of the sources to total B accumulation can be estimated from an isotopic mass balance perspective as follows: δ 11 B f δ 11 B marin aerosol + 1 − f δ 11 B coal burning where the δ 11 B marine aerosol and δ 11 B coal burning is the δ 11 B in aerosol derived from seawater evaporation (δ 11 B 24.3‰) and coal burning (δ 11 B 8‰), respectively. f is the fraction of B derived from seawater evaporation. During the NE monsoon with the exception of dust storm period, the B average contributions derived from seawater evaporation and coal burning from LRT was 57% (17-98%) and 43% (2-83%), respectively ( Table 2). Samples taken during dust storm and volcanic emissions events were inadvisable for estimation of the estimate contributions for the sources from this formula, i.e., these samples contain dust  2 | δ 11 B, analytical errors (2SE, n 10), temperature and the estimated relative contributions (%) from seawater evaporation, long-range transported (LRT) in selected aerosols at PCY. Samples were taken during NE monsoon (NE), summer monsoon (S), typhoon, dust storm and volcanic eruption periods. The mean temperature data for the collected day was retrieved from NOAA GSOD (https://www.ncdc.noaa.gov/). or volcanic gas with unknown δ 11 B. As we mentioned above, Rayleigh-like evolution would be the dominant B fractionation process during the summer monsoon. The B contents in the summer monsoon samples demonstrated no significant differences, which implied that these samples were not mixtures of high B content endmembers (i.e. coal). From our perspective, most B came from seawater evaporation that occurred in the summer monsoon period.

CONCLUSION
In this study, we presented various B sources, and the characteristics of B and δ 11 B, in atmospheric aerosols at PCY. Our major observations are summarized below: 1) Seawater evaporation and coal burning are the major sources, characterized by high δ 11 B (δ 11 B 24.3‰) and low δ 11 B (δ 11 B 8‰), respectively. Climatic events including volcanic emissions and dust storms significantly influenced the B and δ 11 B over short periods of time during the events.
2) The Rayleigh distillation processes would be a dominant cause changes in of B and δ 11 B in the summer. The mixing of seawater evaporation and coal burning from China has played an important role during the NE monsoon in the winter.
3) The B in aerosols were derived mostly from seawater in the summer. Average B contributions were derived from marine and LRT coal burning at 57 and 43%, respectively, during the NE monsoon.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
C-FY, Y-PL, and T-YK conceived and planned the experiments. C-HC, Y-PL, and T-YK carried out the experiments and sample preparation. C-FY, Y-PL, and T-YK contributed to the interpretation of the results. C-FY and Y-PL took the lead in writing the manuscript. C-HC has helped to improve the English in the manuscript. All authors provided critical feedback and helped shape the research, analysis, and manuscript.

FUNDING
This study is supported by grants from the Ministry of Science and Technology, Taiwan to C-FY (MOST 108-2116-M-006-011).