Chemical characterisation of water-soluble ions in atmospheric particulate matter on the east coast of Peninsular Malaysia

. Air quality on the east coast of Peninsular Malaysia is inﬂuenced by local anthropogenic and biogenic emissions as well as marine air masses from the South China Sea and aged emissions transported from highly polluted East Asian regions during the winter monsoon season. An atmospheric observation tower has been constructed on this coastline at the Bachok Marine Research Station. Daily PM 2 . 5 samples were collected from the top of the observation tower over a 3-week period, and ion chromatography


Introduction
The tropical Maritime Continent, a region in Southeast Asia between 10 • S-20 • N and 90-150 • E, is a complex distribution of islands and peninsulas and incorporates countries such as Malaysia, Indonesia, the Philippines and Papua New Guinea (Neale and Slingo, 2003). It lies within a tropical warm pool that extends eastwards from the Indian Ocean to the Western Pacific and is home to some of the warmest ocean temperatures in the world. Tropical regions such as the Maritime Continent are of central importance for the chemistry-climate system (Carpenter et al., 2010). For example, high photochemical activity in these regions means that global atmospheric lifetimes of key atmospheric species, such as methane and ozone, are controlled by destruction rates in the tropics (Lawrence et al., 2001;Bloss et al., 2005). In terms of ocean productivity, the observed decrease in primary productivity in low-latitude oceans has been linked to a reduced availability of nutrients for phytoplankton growth, caused by changes in upper-ocean temperature and stratification (Behrenfeld et al., 2006). Furthermore, the wind circulation system in the Maritime Continent is influenced by seasonal Asian monsoons, which are controlled by the natural oscillation of the intertropical convergence zone (ITCZ). During the Northern Hemisphere winter, a large anticyclone forms over Siberia each year, creating strong northeasterly monsoon winds in the South China Sea (northeast monsoon). These strong northeasterlies can transport air masses from rapidly developing East Asian countries (e.g. China, Japan, Vietnam, and North and South Korea) across the South China Sea to the Maritime Continent (Zhang et al., 1997;Garreaud, 2001;Oram et al., 2017). In addition, cold-surge events occur regularly throughout the winter monsoon season and last several days. Cold surges occur as a result of a southeasterly movement of the anticyclone and are characterised by cold air masses over Southern China and strengthening of the northeasterly monsoon winds in the South China Sea (Zhang et al., 1997). The transport of pollution from East Asia to the tropics during the monsoon season, particularly during cold surges, means that rural areas such as the east coast of Peninsular Malaysia are potentially at an elevated risk of the detrimental effects of poor air quality.
Tropical regions are highly important for atmospheric research, and whilst long-term atmospheric observations exist (Robinson et al., 2014;Pyle et al., 2011), there are fewer measurements than in the mid and high latitudes. The Bachok Marine Research Station (lat 6.00892, long. 102.42504) has been set up on the east coast of Peninsular Malaysia and is ideally located for studying the outflow of these highly industrialised regions and for investigating the interaction with cleaner air in the Southern Hemisphere. The research station, located approximately 30 km away from Kota Bharu, forms part of the Institute of Ocean and Earth Sciences at the University of Malaya (UM). An atmospheric observation tower facing the South China Sea has been constructed at the research station; this has been built for the specific purpose of monitoring long-range transported pollution, airsea exchange and coastal meteorology. The research station is working towards designation as a regional Global Atmospheric Watch (GAW) centre, which will be a valuable addition to the network of other global and regional GAW sites in the Maritime Continent, as shown in Fig. 1 (GAW-SIS, 2017).
In January and February 2014, an instrument demonstration campaign was carried out to assess the capabilities of the new research station. This was funded by the Natural Environment Research Council (NERC) and UM and involved several UK universities, as well as the National Cen- tre for Atmospheric Science (NCAS), UM and the Malaysian Meteorological Department (MMD). As part of this study, Dunmore et al. (2016) used a specialised multidimensional gas chromatography technique to accurately measure atmospheric mixing ratios of C 5 -C 13 volatile organic compounds (VOCs) with a wide range of functionalities. Furthermore, Dominick et al. (2015) characterised the particulate matter in Bachok by studying the influence of northeasterly winds on the patterns of particle mass and particle number concentration size distributions. Both studies highlighted the fact that the site is influenced by a mixture of local anthropogenic and biogenic emissions, clean marine air masses, and aged emissions transported from East Asia.
To extend upon these studies, it is important to investigate atmospheric aerosol composition in the Bachok region. A better understanding of aerosol chemical composition is essential as aerosols play an important role in atmospheric processes and climate change. For example, aerosols can modify the global radiation budget both directly, by scattering and absorbing solar radiation, and indirectly, by altering cloud properties and lifetime (Charlson et al., 1991). The strength of these direct and indirect effects depends partly on the particle concentration and size distribution, but also on the chemical composition.
There are a limited number of studies that focus on particulate matter composition on the east coast of Peninsular Malaysia, and to our knowledge the composition of ionic species has not been determined at any rural locations along this coastline. For example, Tahir et al. (2013) studied the composition of major elements and water-soluble ionic species in PM 2.5 and PM 10 samples on the east coast of Peninsular Malaysia, but the samples were collected at an urban coastal city, Kuala Terengganu. This study used principal component analysis to determine the main sources of both fine and coarse particles, which were found to be soil dust, marine aerosol, vehicle exhaust, secondary aerosol, road dust and biomass burning. In addition, Ismail et al. (2016) studied PM 10 concentrations in three major cities (Kota Bharu, Kuala Terengganu and Kuantan) on the east coast of Peninsular Malaysia between 2006 and 2012. The study showed that during the northeast monsoon, the air arriving at the sites had originated from China and the Philippines and travelled over the South China Sea. During the southwest monsoon, the air came from Indonesia via the Strait of Malacca. Over the 6-year period, it was found that the atmospheric PM 10 mass was directly proportional to the rate of urbanisation in each of the three cities.
In this study, measurements of water-soluble ions in atmospheric aerosol at a rural coastal location on the east coast of Peninsular Malaysia are presented. Analysis of temporal variation of different ionic species has been carried out and backward air mass trajectories have been used to determine the influence of air mass origin on aerosol composition. Correlation analyses of the ionic species and assessment of ratios between different ions have provided an insight into common sources and formation pathways of key atmospheric ions. This study contributes towards a better understanding of atmospheric composition in the rural Bachok region, which is vulnerable to the effects of poor air quality, largely as a result of rapid industrialisation in East Asia.

Sample collection and extraction
Thirty PM 2.5 samples were collected at the Bachok Marine Research Station (lat 6.00892, long. 102.42504) between 18 January and 7 February 2014. The samples were collected at the top of an atmospheric observation tower (18 m height) using a high-volume air sampler (Ecotech HiVol 3000, Victoria, Australia) operating at 1.13 m 3 min −1 over 24 h sampling intervals. The tower is located on the coastline of the South China Sea and is within 100 m of the shore. A 3-day intensive measurement period was in operation between midday (local time) on 30 January 2014 and midday on 2 February 2014, in which filters were collected every 4-8 h. The quartz fibre filters (20.3 cm × 25.4 cm) supplied by Whatman (Maidstone, UK) were pre-baked at 550 • C for a minimum of 12 h prior to sample collection. After sample collection, the filters were wrapped in aluminium foil and stored at −18 • C until analysis. To prepare the samples for analysis, 5.7 cm 2 of each sample was dissolved in 2 mL Milli-Q water and sonicated for 30 min at room temperature. The extract was filtered using a Millex-GP 33 mm diameter hydrophilic syringe filter with a pore size of 0.22 µm (Millipore UK Limited, Watford, UK) and made up to a final volume of 2.5 mL.

Eluents and analytical standards
Ultrapure Milli-Q water (18 M cm −1 ) from an ELGA Lab-Water purification system was used to prepare all the required eluents and analytical standards. A 20 mM solution of methanesulfonic acid was used as the eluent for cation exchange chromatography, and for anion exchange chromatography, a solution of 8 mM Na 2 CO 3 / 1 mM NaHCO 3 was prepared. Using a variety of salts and organic acids, individual analytical standards containing 500 ppm of each target ion (Cl − , NO − 2 , NO − 3 , PO 3− 4 , SO 2− 4 , CH 3 SO − 3 , C 2 O 2− 4 , Na + , NH + 4 , K + , Mg 2+ and Ca 2+ ) were prepared in Milli-Q water. The salts and organic acids were purchased from either Sigma-Aldrich Ltd. (Dorset, UK) or Fisher Scientific Ltd. (Loughborough, UK).

Chromatographic analysis
Chromatographic analysis was carried out using a Thermo Scientific Dionex ICS-1100 ion chromatography system equipped with an AS-DV autosampler. The column configuration used for anion exchange consisted of an Ion-Pac AG14A guard column (4 mm × 50 mm) and an Ion-Pac AS14A analytical column (4 mm × 250 mm). Cation exchange chromatography was performed using an IonPac CG12A guard column (4 mm × 50 mm) and an IonPac CS12A analytical column (4 mm × 250 mm). ASRS 300 and CSRS 300 self-regenerating suppressors (4 mm) were used for anion and cation exchange respectively. All columns and suppressors were supplied from Thermo Scientific Dionex. The run times for the anion and cation separations were 18 and 15 min respectively. The suppressor current was 45 mA for anion exchange mode and 59 mA for cation exchange mode. For all separations, the instrument was operated in isocratic mode at a flow rate of 1 mL min −1 and a column oven temperature of 30 • C. The injection volume was 100 µL and the data collection rate was 5 Hz. The system relied on a DS6 heated conductivity cell for ion detection, and all data were analysed using Thermo Scientific Chromeleon 7.1 Chromatography Data System software.

Method validation
Using isocratic elution methods for both cation and anion exchange chromatography, the target ions were successfully separated. Recovery tests were performed by spiking 5.7 cm 2 of quartz fibre filters (Whatman, Maidstone, UK) with 1 µg of each target ion (20 µL of a 50 ppm mixed ion solution). Prior to spiking, the filters were pre-baked at 550 • C for 6 h and wrapped in aluminium foil and stored at −18 • C until required. The spiked filters were dissolved in 2 mL of Milli-Q water and sonicated for 30 min at room temperature. The extract was filtered using a Millex-GP 33 mm diameter hy-drophilic syringe filter with a pore size of 0.22 µm (Millipore UK Limited, Watford, UK) and made up to a final volume of 2.5 mL. The recovery of the target ions from the filter papers ranged from 74.5 % to 98.2 % for the target anions, and 78.3 % to 87.3 % for the target cations, with the exception of Ca 2+ for which a recovery level of 123.3 % was calculated. Recovery tests were carried out in triplicate, and %RSD rec remained below 8 % for all the ions. The recovery of Ca 2+ should not have exceeded 100 %, and the result may be attributed to inconsistencies in the amount of Ca 2+ present on the blank filter, or due to Ca 2+ contamination during sample collection or storage. A 100 % recovery was assumed for Ca 2+ during the data analysis process. Further details of the individual recovery levels and associated errors can be found in the Supplement (Table S1). Procedural blanks were also carried out using quartz fibre filters (5.7 cm 2 ), and blank subtractions were applied to any target ions found in detectable amounts. The blank peak areas for each ion and average blank contribution to field samples over the entire sampling period are provided in the Supplement, Table S1. The main instrumental parameters of the IC system were evaluated and are also detailed in the Supplement (Table S2). Instrumental limits of detection (LODs) and limits of quantification (LOQs) were calculated according to the EPA protocol 40 CFR 136; multiplying the standard deviation (N = 10, 5 ng for cations, 25 ng for anions) by the t value (N = 10, 95 % confidence interval) gave the LOD, and multiplying the standard deviation by 10 gave the LOQ (EPA, 2017;Ripp, 1996). For anion exchange chromatography, LODs and LOQs were in the range 5.5-21.0 ng and 25.3-144.2 ng respectively. For cation exchange chromatography, the LODs ranged from 0.5 to 2.1 ng and the LOQs ranged from 2.5 to 6.1 ng. On average, the instrument precision (%RSD ins , n = 10) was 4.2 % for the target cations and 12.8 % for the target anions. Total errors were estimated by combining errors with the instrument and the recovery process and remained below 15.4 % for all ions except NO − 3 (22.6 %).

Additional measurements
Individual VOCs were measured using a combined heartcut and comprehensive two-dimensional gas chromatography system (GC − GC × GC); a detailed description of the instrument design is provided in a separate study (Dunmore et al., 2016). Measurements of NO and NO 2 were performed using a two channel 42i-TL commercial gas analyser (Thermo Scientific, MA, USA), and SO 2 measurements were made using a Thermo Scientific 43i SO 2 analyser. Meteorological data from the nearest meteorological station, the Sultan Ismail Petra Airport (lat 6.17208, long. 102.29288), were accessed from the Integrated Surface Database (NOAA, 2003). Hourly measurements of wind direction, wind speed, air temperature, dew point, atmospheric pressure and relative humidity were obtained. The 10-day backward air mass trajectories arriving at the sampling site were run every 3 h throughout the entire measurement period. A receptor height of 10 m was chosen to represent the measurements made on the sampling tower. The trajectories were computed using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model (Stein et al., 2015;Draxler, 1999;Draxler andHess, 1998, 1997), and the data were analysed using the openair package in RStudio (Carslaw and Ropkins, 2012;Carslaw, 2015).

ISORROPIA II model
Predictions of aerosol pH were made using the ISORROPIA-II thermodynamic equilibrium model. Although these models produce better results when gas phase measurements such as NH 3 and HNO 3 are available, it is possible to use the ion measurements obtained in this study to make a prediction of aerosol pH (Fountoukis and Nenes, 2007). Calculations were made in reverse mode, in which known quantities are temperature, relative humidity, and particle phase concentrations of NH + 4 , SO 2− 4 , Na + , Cl − , NO − 3 , Ca 2+ , K + and Mg 2+ . ISOR-ROPIA II assumes that the particles are internally mixed; this is a reasonable assumption for this study, as relative humidity was high (average = 77 %), and the aerosol arriving at the measurement site is often aged. The aerosol was assumed to be thermodynamically stable; i.e. the aerosol can exist as both solid and liquid, and salts are able to precipitate if the aqueous phase becomes saturated with respect to them. The ambient temperature and relative humidity data were taken from the measurements made nearby at the Sultan Ismail Petra Airport. It is likely that these measurements are representative of the Bachok research station, as further investigation of data from two other meteorological stations showed that temperature and relative humidity remain consistent along the coastline.

Bachok demonstration campaign
The filter samples were collected at the Bachok atmospheric observation tower. The Bachok district, located in the state of Kelantan, is a rural area and the primary economic activity comes from tobacco and kenaf plantations. Other agrarian activities in the wider Kelantan region include the production of rice and rubber as well as additional economic activities such as livestock rearing and fishing. Figure 2 shows the 7-day backward air mass trajectories arriving at the measurement site during the demonstration campaign, with each trajectory coloured by the mean altitude (m). For comparison with pollutant trends later in this study, the trajectories have been separated into two plots according to their mean height. As expected during the winter months, the strong anticyclone system known as the Siberian High led to the arrival of northeasterly onshore winds along the east coast of Peninsular Malaysia. Some of the air masses experienced a  significant continental influence from highly industrialised countries such as China, Japan, North and South Korea, and the island of Taiwan, whilst other air masses had a stronger marine influence from both the East China Sea and the South China Sea.
Although there was no reliable meteorological data recorded at the measurement site during the demonstration campaign, data from a nearby meteorological station were available. The station is located approximately 23 km away at the Sultan Ismail Petra Airport in Kota Bharu, as shown in Fig. S1 (Supplement). Whilst the meteorological data from the Airport will not be exactly representative of the measurement site, the patterns in wind direction are consistent with observations made by the field scientists during the campaign. In addition, during a study of the influence of northeast monsoon cold surges on air quality in Southeast Asia, Ashfold et al. (2017) used meteorological data from three locations that they believed to lie in the path of cold surges during the northeast monsoon, offering the best possibility of observing a cold-surge influence on air pollution. These sites were in Kota Bharu (lat 6.141, long. 102.247), Kuala Terengganu (lat 5.308, long. 103.118) and Kemaman (lat 4.271, long. 103.428); this provides further confirmation that data from the Sultan Ismail Petra Airport site are suitable for evaluating broad-scale transport at the Bachok measurement site. Figure S2 in the Supplement shows average hourly wind speed and wind direction conditions across the entire duration of the measurement campaign (18 January to 6 February 2014). On most days, gentle southwesterlies from the land (Peninsular Malaysia) were observed in the early hours of the morning to around 11:00 LT (Malaysia time). At this stage, a dramatic shift in wind direction occurred, Table 1. Mean and maximum ion concentrations measured throughout the measurement period. The average percent mass contribution of each ion to the total measured ions is included as well as the percent of samples in which each target ion is found (%Qt).

Ion
Mean ( and until around 19:00 LT a strong onshore breeze from the northeast (South China Sea) was observed. From 20:00 LT to the early hours of each morning, a calmer sea breeze predominantly from the east was seen. Hourly VOC measurements were conducted at the measurement site, and the development of a sea breeze at approximately 11:00 LT dramatically influenced the diurnal profiles of the measured species (Fig. S3, Supplement). In the morning, when the air being sampled was coming over the land, high levels of NO, NO 2 and anthropogenic VOCs such as toluene and C 10 aliphatics were observed; the main source of these species the was local burning of waste (Dunmore et al., 2016). When the sea breeze developed, the concentration of these species dropped significantly.

Aerosol composition and determination of non-sea-salt and sea salt components
The total concentration of measured water-soluble ions in the PM 2.5 during the campaign ranged from 8.06 to 27.0 µg m −3 , with an average concentration of 16.2 µg m −3 . Table 1 shows the mean and maximum water-soluble ion concentrations measured throughout the campaign, and Fig. 3 shows time series for all the ions measured in the aerosol.
As the composition of the water-soluble ions present in aerosol collected at the Bachok site was influenced by both marine and continental sources, it is useful to make an estimation of non-sea-salt (nss) and sea salt (ss) components, using Eqs. (1)-(4). Total Na + and Ca 2+ concentrations have been measured in this study, and the mean Ca 2+ /Na + ratio in the crust and mean Ca 2+ /Na + ratio in seawater have been estimated as 1.78 w/w and 0.038 w/w respectively (Bowen, 1979). Therefore it is possible to solve Eqs. (1)-(4) simultaneously for ssNa + , nssNa + , ssCa 2+ and nssCa 2+ (Boreddy and Kawamura, 2015). Furthermore, the resulting estimate of ssNa + , which can be used as a sea spray marker, can also be used to predict the contribution of nssSO 2− 4 and nssK + in the aerosol, as shown in Eqs. (5) and (6) respectively. Whilst there may be some uncertainty in the mean Ca 2+ /Na + crustal ratio, due to challenges associated with predicting the composition of the crust (Bowen, 1979), this approach provides more accuracy than simply using total Na + as a sea spray marker. Figure 4 shows a series of stacked bar charts to summarise the average mass composition of water-soluble ions in atmospheric aerosol and the distribution of non-sea-salt and sea salt components. The results show that the water-soluble ion fraction of the aerosol is dominated by SO 2− 4 , which on average made up 65.6 % of the total ion content by mass. NH + 4 and NO − 3 concentrations were significantly lower, with mean concentrations of 1.69 and 0.61 µg m −3 respectively. Na + and Cl − made up 11.1 % of the total ion content, and 72 % of the measured Na + was attributed to ssNa + . The average concentrations of nssK + and nssCa 2+ were 0.30 and 0.09 µg m −3 respectively; these ions can be used as tracers for biomass burning (nssK + ) and atmospheric dust (nssCa 2+ ). NO − 2 and CH 3 SO − 3 were the least abundant ions, with average concentrations of 0.05 and 0.08 µg m −3 . The two ions were only observed in a subset of samples; CH 3 SO − 3 was quantified in 67 % of samples and quantification of NO − 2 was only achieved in 23 % of samples.

Sources and formation of sulfate (SO 2− 4 )
The average SO 2− 4 concentration during the measurement campaign was 10.7 µg m −3 , with a maximum concentration of 20.8 µg m −3 recorded. The formation of SO 2− 4 in the particle phase occurs when emitted SO 2 is oxidised by OH in the gas phase, or by O 3 or H 2 O 2 in the aqueous phase (Fisher et al., 2011). The most dominant anthropogenic sources of SO 2 include fuel and industrial emissions, as well as open biomass burning.
By using ssNa + as a sea spray marker to determine nonsea-salt and sea salt components of the aerosol, it was found that on average 96 % of the measured SO 2− 4 was nssSO 2− 4 , and only 4 % of the SO 2− 4 was from sea salt. As a potential biogenic source of nssSO 2− 4 is DMS emissions from marine biota, and the main atmospheric source of MSA is the oxidation of DMS, it is possible to use the MSA − /nssSO 2− 4 ratio as a tracer to assess the contribution of biogenic sources to nssSO 2− 4 in the atmosphere (Legrand and Pasteur, 1998). In this study, MSA − concentrations ranged from 0.02 to 0.22 µg m −3 , with an average concentration of 0.08 µg m −3 . As a result, the MSA − /nssSO 2− 4 ratio ranged from 1.6 × 10 −3 to 2.2 × 10 −2 , with an average value of 8.0 × 10 −3 . These values are low compared to MSA − /nssSO 2− 4 ratios recorded at remote sites; for example, a MSA − /nssSO 2− 4 mean mass ratio of 7 × 10 −2 has been measured on Fanning Island and American Samoa (Savoie and Prospero, 1989).
The lower ratios recorded in Bachok suggest the majority of nssSO 2− 4 at the site originates from anthropogenic sources. As previously discussed, a study by Dunmore et al. (2016) revealed that levels of NO x and anthropogenic VOCs at the Bachok measurement site were significantly higher when the air being sampled had passed over nearby land and dropped significantly at around 11:00 LT when a sea breeze developed. This indicated that air quality in Bachok is influenced by local sources of pollution, such as vehicle emissions and burning domestic waste. Pollution rose plots (Fig. 5) show the relationship between wind direction and speed and SO 2 concentration at the Bachok measurement site. The left panel shows recorded SO 2 concentrations ≥ 5 ppb and for comparison, the right panel shows all SO 2 data recorded during the measurement period. SO 2 concentrations < 5 ppb have been excluded from the left panel plot to further investigate the meteorological conditions when the spikes in SO 2 concentration occur. In the lower right corner of each pollution rose, "mean" represents the mean SO 2 concentration and "calm" represents the fraction of data that cannot be attributed to a specific wind direction. During the elevated SO 2 periods (≥ 5 ppb), weak southwesterlies dominated and the average wind speed was 1.1 m s −1 . The mean SO 2 concentration was 24.7 ppb during these events. Over the entire measurement campaign, however, average wind speed was considerably higher, 2.8 m s −1 , and the air arrived predominantly from the east. In summary, the majority of higher SO 2 events were observed in calmer conditions when the air arriving at the site had passed over nearby land; these observations provide further evidence for the influence of local sources of pollution.
At the Bachok measurement site, no obvious relationship was observed between SO 2 and particulate SO 2− 4 concentration or between SO 2− 4 concentration and wind direction. This is likely to be because it takes time for SO 2 to oxidise to SO 2− 4 and that the SO 2− 4 fraction of the aerosol is more heavily influenced by long-range transport of aged emissions from East Asia. To investigate this further, the backward air mass trajectories were coloured by the concentration of Figure 5. Pollution rose plots to show the relationship between wind direction and SO 2 concentration at the Bachok measurement site. Panel (a) represents recorded SO 2 concentrations ≥ 5 ppb and panel (b) represents all SO 2 data recorded during the measurement period. Plot constructed using the openair package in RStudio (Carslaw and Ropkins, 2012;Carslaw, 2015). SO 2− 4 , as shown in Fig. 6. To allow for comparison with Fig. 1, trajectories with a mean height of 80 m or below are shown on the left panel and trajectories with a mean height greater than 80 m are shown on the right panel. At the lower altitudes (≤ 80 m), it is clear that the SO 2− 4 content of the aerosol is highest (ca. 15-20 µg m −3 ) when the site is influenced by continental air masses from industrialised regions of East Asia and lower when the air masses have a more significant marine influence. This trend is less clear at higher altitudes (right panel, mean height > 80 m), but nevertheless some of the lowest SO 2− 4 levels are observed (∼ 5 µg m −3 ) when marine air that has passed over the East and South China seas arrives at the site. With this information available, it is useful to perform cluster analysis on the back trajectories; this type of analysis groups air masses of similar geographic origin together, which provides more information on pollutant species with similar chemical histories. A distance matrix is used to create a specified number of clusters with the most different air mass pathways. Figure 7 shows the fivecluster solution to back trajectories calculated for the Bachok site during the measurement campaign.
Clusters 1 and 3 were associated with mean SO 2− 4 concentrations of 14.4 and 13.8 µg m −3 respectively. The mean altitude of cluster 1 was 72 m and the mean altitude of cluster 3 was 501 m. Full details of the air mass trajectories within each cluster are detailed in the Supplement (Table S3). These clusters contained air masses that had passed over several highly industrialised regions en route to Bachok, including cities such as Guangzhou (China) and Ho Chi Minh City (Vietnam). Average SO 2− 4 concentrations for clusters 2 and 4 were 8.4 and 8.3 µg m −3 , with mean altitudes of 169 and 37 m respectively. These air masses are likely to have experienced some continental influence from the east coast of China and the island of Taiwan but have a much more significant marine influence from passing over the South China Sea before arriving at Bachok. For air masses within clusters 1-4, the concentration of SO 2− 4 increases with decreas-ing altitude and increasing influence from industrialised regions. The mean altitude of air masses within cluster 5 was 1027 m, i.e. significantly higher than air masses in clusters 1-4. Furthermore, the air masses in cluster 5 arriving at the site were not a result of northeasterly monsoon winds from East Asian countries travelling across the South China Sea. Instead, there is evidence of a cyclonic weather system off the coast of the island of Borneo over the South China Sea. The air masses in cluster 5 arrived during the final 24 h period of the measurement campaign, and only one SO 2− 4 measurement is available (18.1 µg m −3 ). It is not possible to fully understand the nature of air masses from this different region without further measurements.
Despite potential differences in atmospheric lifetimes and behaviour, other pollutants aside from SO 2− 4 can provide further evidence that substantial amounts of industrial pollution from East Asia are undergoing atmospheric transport to tropical regions of the Western Pacific. For example, Oram et al. (2017) measured chlorine-containing, very short-lived substances (Cl-VSLSs) at the Bachok research station during the winter monsoon season in late January-early February 2014. Cl-VSLSs are ozone-depleting species with short atmospheric lifetimes, typically less than 6 months. Species include dichloromethane (CH 2 Cl 2 ) and 1,2-dichloroethane (CH 2 ClCH 2 Cl). A 7-day pollution or cold-surge event was reported between 19 and 26 January 2014, when significantly enhanced concentrations of Cl-VSLSs were observed. During this pollution episode, the measured samples were heavily impacted by emissions from the East Asian mainland, whilst this influence was less significant during the cleaner, non-polluted periods. In fact, the total median concentration of the four measured Cl-VSLSs was 546 ppt between 20 and 26 January and 243 ppt during the less polluted period (27 January to 5 February 2014). Oram et al. (2017) noted that even after the cold-surge event, the levels of Cl-VSLSs were still significantly higher than expected, indicating that this region of the South China Sea is widely im-  pacted by emissions from East Asia. Many other chemical pollutants, aside from short-lived chlorinated gases, will be present in these air masses from East Asia and will have a large impact on regional air quality. Oram et al. (2017) performed a Numerical Atmospheric dispersion Modelling Environment (NAME) trajectory analysis using carbon monoxide (CO) as a tracer of industrial emissions from regions north of 20 • N for six winter seasons (2009-2010 to 2014-2015). A strong correlation between CO and CH 2 Cl 2 (a measured Cl-VSLS) was observed during the pollution episode in late January 2014. Analysis of CO time series over the six winter seasons revealed that cold-surge events are likely to be repeated regularly each winter, demonstrating that pollution rapidly undergoes long-range transport across the South China Sea on a regular basis during the northeast monsoon.

Correlation of SO 2− 4 with NH + 4 and implications for aerosol acidity
Ammonium (NH + 4 ) was the second most abundant ion in the aerosol; on average it made up 10.4 % of the total ion content, and the mean and maximum concentrations were 1.69 and 4.73 µg m −3 respectively. The strong positive correlation between SO 2− 4 and NH + 4 was observed (R = 0.77, p < 0.001). A similar observation was reported by Keywood et al. (2003) during an investigation of the sources of particles contributing to haze in the Klang Valley, Malaysia. The strong relationship between these species is due to neutralisation of SO 2− 4 by NH + 4 . Local sources of NH 3 in the rural Bachok region are likely to come from agricultural practices such as animal husbandry, fertiliser use and agricultural waste burning. The average NH + 4 /SO 2− 4 molar ratio was 0.81, which indicated that there was insufficient gaseous NH 3 in the atmosphere to neutralise SO 2− 4 . Under an ammonia-poor regime, the uptake of SO 2− 4 is preferential to the uptake of NO − 3 because sulfuric acid has a lower vapour pressure than nitric acid (Seinfeld and Pandis, 2006). Although measurements of the total amounts of ammonia and sulfate in the gas, aqueous and solid phases would provide a better prediction of the aerosol acidity, the results presented in this study indicated that an ammonia-poor regime exists and that the aerosol is likely to be acidic. In these scenarios, the NH 3 partial pressure is low, and therefore the NH 3 -HNO 3 partial pressure product is also low, meaning that the concentrations of ammonium nitrate are low or zero (Seinfeld and Pandis, 2006). This hypothesis can be supported by the fact that NO − 3 concentrations in this study were very low, ranging from 0.005 to 1.52 µg m −3 . As a result, NO − 3 made up a significantly smaller fraction of the total ion content compared to SO 2− 4 ; the average percentage mass composition of the total ion content was 3.8 % for NO − 3 and 65.6 % for SO 2− 4 .
To estimate proton loading in atmospheric particles, the strong acidity approach can be used, as shown in Eq. (7). This approach assumes that any deficit in measured cation charge compared to measured anion charge can be attributed to H + . Total anion and total cation equivalents can be estimated using Eqs. (8)  As shown in Fig. 8, strong acidity values ranged from 0.03 to 0.19 µeq. m −3 , with an average value of 0.11 µeq. m −3 . The positive strong acidity values provide an initial indication that the aerosol is acidic and allows an estimate of the proton loading to be made (average H + = 0.11 µg m −3 ). To obtain error bars for the strong acidity predictions (Fig. 8), H + max and H + min were calculated according to Eqs. (10) and (11) respectively. For H + max the anions were adjusted up to within their uncertainties (i.e. +%RSD tot ), and the cations were adjusted down to within their uncertainties (i.e. −%RSD tot ). For H + min the anions were adjusted down and the cations were adjusted up (Murphy et al., 2017). The %RSD tot value was estimated by combining the error of the recovery process for each ion (%RSD rec , Table S1) and the error of the instrument for each ion (%RSD ins , In most cases, H + min remains above zero. However, between 2 and 6 February 2014, slightly negative H + min values between −1 × 10 −3 and −2 × 10 −2 µeq. m −3 were calculated, which are physically implausible (Murphy et al., 2017). These results highlight the possible sources of error associated with the strong acidity approach for estimating aerosol acidity. For example, Hennigan et al. (2015) report that organic acids (which are mostly excluded from this study, except for MSA and oxalic acid) can have an important influence on aerosol acidity, especially at relatively low acidities where organic acids dissociate and contribute to the ion balance. Furthermore, they can form salt complexes with inorganic species e.g. ammonium oxalate. Neglecting organic acids, as well as other atmospheric species such as HCO − 3 and basic amines, will lead to inaccuracies in the calculated H + .
Thermodynamic equilibrium models such as ISOR-ROPIA II can also be used to predict aerosol pH (Fountoukis and Nenes, 2007). The ambient temperature and relative humidity data recorded nearby at the Sultan Ismail Petra Airport were assumed to be representative of the Bachok research station, as further investigation of data from two other meteorological stations revealed that temperature and relative humidity remained consistent. A map to show the location of the three meteorological stations along the east coast, as well as the Bachok research station, can be found in the Supplement (Fig. S1). The two other stations are Narathiwat Airport (lat 6.520, long. 101.743) and Sultan Mahmud Airport (lat 5.383, long. 103.103). Average relative humidity between 18 January and 7 February 2014 was 77.7 % at Sultan Petra Ismail Airport, 75.4 % at Narathiwat Airport and 77.0 % at Sultan Mahmud Airport. Average temperatures recorded at the stations during this time were 24.8, 25.6 and 25.4 • C for Sultan Petra Ismail, Sultan Mahmud and Narathiwat airports respectively.
The ISORROPIA II thermodynamic model predictions of PM 2.5 pH are shown in Fig. 9. Particle pH was estimated with ISORROPIA II run in the reverse mode without gas phase species input and ranged from −0.97 to 1.12 during the measurement period, implying that the aerosol was highly acidic. The pH prediction for the aerosol collected between midday on 3 and 4 February 2014 was 7.06 and has been excluded from Fig. 9. The particle concentrations input on this day correspond to negative values of strong acidity, and therefore the model balances charge by assuming [OH − ] > [H + ]; this leads to a calculated pH of greater than 7. Murphy et al. (2017) reported that pH prediction is sensitive to strong acidity in the limit of strong acidity approaching zero and that the model can be drastically improved if gas phase NH 3 and HNO 3 measurements are included. The gas-to-particle partitioning of these species is sensitive to pH under conditions commonly encountered in the atmosphere; therefore gaseous NH 3 and HNO 3 measurements provide better constraint on the thermodynamic model. Acidic particles can have detrimental effects on human health, air quality, and the health of aquatic and terrestrial ecosystems (Hennigan et al., 2015). For example, Gwynn et al. (2000) performed a time-series analysis of acidic PM and daily mortality and morbidity in the Buffalo, New York, region; several significant pollutant-health effect associations were identified, the strongest being the correlation between atmospheric SO 2− 4 concentration and respiratory hospital admissions. Furthermore, particle acidity can influence various atmospheric chemical processes, including SO 2 oxidation, halogen chemistry, and the partitioning of ammonia, nitric acid, organic acids and isomeric epoxy diols from isoprene photooxidation (IEPOX) (Hennigan et al., 2015;Surratt et al., 2010;Lin et al., 2012). In summary, whilst some of the risks associated with aerosol acidity in Bachok origi-  nate from local sources of pollution, it is possible that people living in these rural areas are also exposed to an additional risk, as the region appears to be sensitive to the effects of industrialisation further afield in East Asia.

Sources and formation of oxalate, C 2 O 2− 4
Oxalic acid is the most abundant dicarboxylic acid in tropospheric aerosol (Sareen et al., 2016). This major watersoluble organic component can alter the hygroscopicity of aerosols and can either act as cloud condensation nuclei (CCN), or reduce the surface tension of particles to form CCN (Saxena and Hildemann, 1996;Novakov and Penner, 1993;Facchini et al., 1999;Kerminen, 2001). In this study, oxalate made up 2.6 % of the total measured water-soluble ion content. The average concentration was 0.42 µg m −3 , and throughout the measurement period the concentration ranged from 0.15 to 0.65 µg m −3 . Interestingly, such levels of oxalate in atmospheric aerosol are typical of urban environments, despite the fact that the Bachok research station is located in a rural coastal region. For example, Freitas et al. (2012) report average oxalate concentrations in PM 10 at an urban site and a rural site in Londrina city, Brazil, of 0.57 and 0.03 µg m −3 respectively. Other reported oxalate con-centrations in urban TSPs (total suspended particles) include measurements of 0.10-0.48 µg m −3 in Shanghai (Jiang et al., 2011) and 0.27 µg m −3 in Tokyo (Sempere and Kawamura, 1994).
To investigate possible oxalate sources and formation pathways, it is necessary to consider the correlation of oxalate with different atmospheric species. Jiang et al. (2011) report using NO − 2 as an indicator for vehicle emissions, nssSO 2− 4 and NO − 3 for secondary formation through different pathways, and K + for biomass burning. A study carried out by Huang et al. (2006) in the urban area of Shenzhen (Southern China) reported that, whilst good correlation of droplet oxalate with K + was observed (R 2 = 0.75, average diameter: 1.0 µm), there was poor correlation between oxalate and K + in the condensation mode (R 2 = 0.10, average diameter: 0.4 µm). This implied that whilst biomass burning was probably not an important primary source of condensation mode oxalate, it is likely that biomass burning particles act as effective CCN, promoting in-cloud sulfate and oxalate formation. In this study, nssK + levels remained below 0.7 µg m −3 ; this likely represents background biomass burning emissions over the region during the measurement period, as evidenced by fire hotspots from satellite data shown in Fig. S4 (Giglio et al., 2003). A strong correlation was ob-served between oxalate and nssK + (R = 0.78, p < 0.001), and the oxalate / nssK + ratio ranged from 0.68 to 4.64 during the Bachok measurement campaign. These ratios are significantly higher than those found in aerosol directly emitted from vegetation fires in the Amazon Basin (Yamasoe et al., 2000), suggesting that biomass burning in the wider region influenced the secondary formation of oxalate, rather than acting as a primary source. A similar hypothesis was proposed by Huang and Yu (2007), who measured ambient PM 2.5 in an urban environment in the Pearl River Delta region of China, and reported oxalate / K + ratios of 0.57 and 0.33 in summer and winter respectively. It is worth noting that whilst the higher ratios reported in both studies indicate that biomass burning is not a major primary source of oxalate, measurements of oxalate and K + in local biomass burning aerosols (rather than aerosol in the Amazon Basin) would provide a better indication of the source contribution by biomass burning.
There was no significant correlation observed between oxalate and NO − 3 , or between oxalate and NO − 2 . The lack of correlation with NO − 3 suggests that the two species do not have similar formation pathways and that vehicle emissions are not an important secondary source of oxalate. It is also unlikely that vehicle emissions contribute to the primary sources of oxalate, due to the lack of correlation between oxalate and NO − 2 . However, NO − 2 was only detected in 7 out of the 30 samples collected, so it is difficult to ascertain whether there is a relationship between these two species or not. A strong correlation between oxalate and nssSO 2− 4 (R = 0.69, p < 0.001) was observed, suggesting a common formation pathway of the two species. It is well known that SO 2− 4 forms via aqueous oxidation (Seinfeld and Pandis, 2006), and modelling studies also suggest that aqueous chemistry is a large contributor of oxalate formation globally (Myriokefalitakis et al., 2011). Furthermore, Carlton et al. (2006) report that whilst there are likely to be many sources of oxalate, oxidation of pyruvate in the aqueous phase is known to form oxalate at dilute (cloud-relevant) concentrations. Tan et al. (2012) also state that aqueous acetate oxidation is a key source of oxalate. A positive correlation between oxalate and NH + 4 was also observed (R = 0.73, p < 0.001). A similar observation was reported by Jiang et al. (2011) in a study of aerosol oxalate in Shanghai. Using size distribution data, they were able to propose that the correlation was due to the presence of ammonium oxalate in the aerosol. In this study there are no size distribution data available, and so it is important to consider the fact that the correlation may be linked to the influence of sulfate on both NH + 4 and oxalate in aerosol; NH + 4 partitions to the aerosol from gaseous NH 3 in an attempt to neutralise acidic sulfate particles, whilst oxalate exhibits similar formation pathways to SO 2− 4 .

Sea salt aerosol and factors affecting chloride depletion
On average, Na + and Cl − contributed 7.0 % and 4.1 % to the total measured water-soluble ion content respectively, and 72 % of the measured Na + was attributed to ssNa + . The concentration of ssNa + ranged from 0.24 to 2.35 µg m −3 , whilst the concentration of Cl − ranged from 0.003 to 2.38 µg m −3 . There was no correlation between nssNa + and Cl − (R = −0.01), but a strong positive correlation between ssNa + and Cl − was observed (R = 0.84, p < 0.001). During the measurement period, the Cl − /ssNa + molar ratio ranged from 0.003 to 1.10 with an average value of 0.40. A time series of Cl − /ssNa + molar ratio can be found in the Supplement (Fig. S5). Significant Cl − depletion was observed, as all of the ratios recorded were lower than that of bulk seawater, 1.18 (Boreddy and Kawamura, 2015). Figure 10 shows backward air mass trajectories arriving at the Bachok research station, coloured by the Cl − /ssNa + molar ratio.
A study of marine aerosol at the remote Chichi-jima island in the western North Pacific during 2001 and 2002 reported that mean Cl − /Na + molar ratios were highest (1.34) in September 2001 and lowest (0.30) in May 2002, and the mean ratio across the 2-year period was 1.10 (Boreddy et al., 2014). Boreddy et al. (2014) reported that the observed chloride depletion was likely due to acid displacement occurring as a result of atmospheric mixing of anthropogenic pollutants such as SO x and NO x . Acid displacement can occur when sea salt particles react with acids such as H 2 SO 4 , HNO 3 , oxalic acid (C 2 H 2 O 4 ) and methanesulfonic acid (CH 3 SO 3 H) in the atmosphere. Such processes are of atmospheric importance as they lead to the formation of gaseous HCl and potentially affect acid deposition conditions in the region. It is widely accepted that Cl − depletion through the volatilisation of HCl via acid displacement occurs particularly in polluted marine air masses (Newberg et al., 2005;Sturges and Shaw, 1993). In this study, the lowest Cl − /ssNa + ratios were observed when air masses arriving at the site had previously passed over highly industrialised countries such as China and Vietnam (Fig. 10). Higher Cl − /ssNa + ratios are found in marine air masses from less polluted regions. It is possible that aerosol, particularly in northern China, may be rich in chloride in the winter due to coal burning. However, this study suggests that the influence of other anthropogenic pollutants, such as HNO 3 and H 2 SO 4 , means that as the air is transported over the South China Sea to the Bachok region, significant Cl − depletion can occur prior to arrival at the measurement site. In a more recent study, Boreddy and Kawamura performed regression analysis between the Cl − /Na + mass ratio and various acidic species, including nssSO 2− 4 , NO − 3 , MSA − and oxalic acid. They found a moderate negative correlation between Cl − /Na + mass ratio and nssSO 2− 4 , and a moderate to weak negative correlation with NO − 3 ; this suggested that sulfate had a higher influence on chloride depletion than nitrate. Furthermore, whilst MSA − moderately correlated with the Cl − /Na + mass ratio in the summer, there was a significant negative correlation between oxalic acid and the mass ratio in the other three seasons, providing confirmation that oxalic acid plays an important role in chloride loss at Chichi-jima island. A similar regression analysis was carried out on the measurements obtained during the Bachok measurement campaign, and whilst a weak negative correlation was found between SO 2− 4 and the Cl − /ssNa + mass ratio (R = −0.41, p < 0.001), there was no correlation of the Cl − /ssNa + mass ratio with oxalate or MSA − . These results imply that whilst H 2 SO 4 played an important role in chloride depletion, methanesulfonic acid and oxalic acid may not have. Interestingly, a strong positive correlation was observed between NO − 3 and Cl − /ssNa + mass ratio (R = 0.82, p < 0.001). These measurements may be linked to each other through the important role of H 2 SO 4 in the atmosphere. Acid displacement, when sea salt reacts with H 2 SO 4 , leads to the removal of Cl − from the aerosol as gaseous HCl, and a partitioning of SO 2− 4 to the aerosol as Na 2 SO 4 . Furthermore, under an ammonia-poor regime (as observed in this study), H 2 SO 4 has a lower vapour pressure than HNO 3 , leading to the preferential formation of ammonium sulfate over ammonium nitrate when there is insufficient NH 3 available to fully neutralise sulfate and nitrate (Seinfeld and Pandis, 2006).

Using nssCa 2+ as a potential tracer for dust episodes
NssCa 2+ can be used as a tracer for atmospheric dust (Boreddy and Kawamura, 2015). The average concentration of nssCa 2+ during the measurement campaign was 0.07 µg m −3 , and the maximum concentration recorded was 0.28 µg m −3 . On most days, the nssCa 2+ concentration was below 0.10 µg m −3 , but between midday on 30 January 2014 and midnight on 31 January 2014 (local time), elevated nssCa 2+ levels were observed, and the average concentration during this period was 0.22 µg m −3 . The same trend was observed for PO 3− 4 ; the average concentration across the entire measurement period was 0.34 µg m −3 , but during this 12 h episode the PO 3− 4 concentration increased to 1.88 µg m −3 . The increase in ion concentration was less pronounced for nssNa + and Mg 2+ , but the concentrations were still 0.26 and 0.88 µg m −3 above the average for the whole measurement period respectively. Time-series plots for nssCa 2+ , PO 3− 4 , Mg 2+ and nssNa + can be found in the Supplement (Fig. S6).
Between midday on 30 January 2014 and midnight on 31 January 2014, the contribution of nssCa 2+ to the total measured water-soluble ion content was more significant (1.26 %) than the average contribution during the remainder of the measurement period (0.35 %). The Bachok research station is located in the outflow region of Asian dusts, and these measurements suggest that long-range transport of Asian dusts over the measurement site has occurred during this time. In fact, during this episode, the back trajectories arriving at the site can be traced back to the North China Plains and the Horqin Desert in eastern China; source attribution studies by Ginoux et al. (2012) have revealed that large anthropogenic dust sources are found in these regions. In the future, it is important that longer term measurements are carried out at the Bachok research station to provide confirmation that dust episodes in these regions of Asia are responsible for the elevated levels of nssCa 2+ and other associated ions. Dust is one of the most abundant types of aerosol in the atmosphere and can have important impacts on both air quality and climate; therefore it is important that seasonal and annual trends in water-soluble ions are studied in more detail at the Bachok research station. centred on the Bachok research station, it was possible to observe the diurnal wind pattern behaviour and assess where the air masses arriving at the site originated from. Air quality at this remote location is influenced by local anthropogenic and biogenic emissions, as well as marine air masses from the South China Sea and aged emissions transported from highly polluted East Asian regions during the winter monsoon season. In general, the site was influenced by southwesterlies coming from the land from the early hours of the morning until approximately 11:00 LT, and then a dramatic shift in wind direction occurred and a sea breeze was present for the remainder of the day. This shift was accompanied by a drop in the concentrations of NO x and anthropogenic VOCs (Dunmore et al., 2016).
Twelve atmospheric water-soluble ions were measured in this study, and SO 2− 4 was found to be the most dominant ion present, making up 66 % of the total measured ion content on average. The non-sea-salt and sea salt components of SO 2− 4 , Na + , K + and Ca 2+ were determined, and it was found that 96 % of the measured SO 2− 4 was non-sea-salt SO 2− 4 . Predictions of aerosol pH were made using the ISORROPIA II thermodynamic model and it was estimated that the aerosol was highly acidic, with pH values ranging from −0.97 to 1.12; such levels of acidity are likely to have a detrimental impact on human health and the health of ecosystems at this remote coastal location. A clear difference in aerosol composition was found between continental air masses originating from industrialised regions of East Asia and marine air masses predominantly influenced by the South China Sea. For example, elevated SO 2− 4 concentrations were observed when continental air masses that had passed over highly industrialised regions of East Asia arrived at the measurement site.
Correlation analyses amongst ionic species and assessment of ratios between different ions provided an insight into common sources and formation pathways of key atmospheric ions. Oxalate concentrations were recorded and found to be more comparable to measurements made at urban locations rather than rural ones. A strong correlation of C 2 O 2− 4 with SO 2− 4 suggested a common aqueous oxidation formation pathway. A strong correlation between C 2 O 2− 4 and K + coupled with high C 2 O 2− 4 /nssK + ratios indicated that biomass burning was an important secondary source of oxalate in the Bachok region, whereas a lack of correlation with NO − 2 and NO − 3 suggested that vehicular emissions were not an important source. The average Cl − /ssNa + molar ratio during the measurement campaign was 0.40, significantly lower than that of bulk seawater (1.18). Analysis of back trajectories revealed that chloride depletion was greater when the aerosol was more influenced by anthropogenic sources of pollution. Elevated levels of nssCa 2+ and other ions such as PO 3− 4 , Mg 2+ and nssNa + were observed between midday on 30 January 2014 and midnight on 31 January 2014. Assuming that nssCa 2+ can be used as a tracer for atmospheric dust, it was proposed that the increased concentrations were a result of air masses arriving at the site from the North China Plains and Horqin Desert. Longer term measurements are required to fully investigate the influence of Asian dusts at this remote coastal location.
To our knowledge, time-resolved measurements of watersoluble ions in PM 2.5 are virtually non-existent in rural locations on the east coast of Peninsular Malaysia. The data presented in this study have demonstrated the capabilities of the new atmospheric tower at the Bachok research station and has provided an initial insight into factors affecting aerosol composition on this coastline. In the future, it is important that longer term measurements are carried out, with increased time-resolved sampling and particle size fractionation, to provide a better understanding of the factors affecting aerosol composition at this measurement site. This remote location is susceptible to the effects of local, regional and international air pollution, and rapid industrialisation in East Asia is influencing air quality along the east coast of Peninsular Malaysia.
Author contributions. All authors contributed to the final version of this article. NJF analysed the aerosol samples and wrote the paper under the supervision of JFH. RED collected the aerosol samples, and BJB made the SO 2 measurements. MIM, MSMN, AAS and SMP coordinate and manage the University of Malaya BMRS. WTS coordinated the Bachok demonstration "International Opportunities Fund" campaign.