{'fraction_non_alphanumeric': 0.031147977879938702, 'fraction_numerical': 0.03058165100939436, 'mean_word_length': 4.522259013980868, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The formation and growth of calcium sulfate crystals through oxidation of SO2 by O3 on size-resolved calcium carbonate

Calcium sulfate is a major constituent of atmospheric sulfate, with a typical rod-like morphology ranging from several hundred nanometers to approximately two micrometers observed in field studies. However, the chemical formation mechanism is still not well known. In this study, the kinetics and mechanism for the formation and growth of rod-like calcium sulfate crystals through oxidation of SO2 by O3 on size-resolved CaCO3 at different relative humidity (RH) were investigated using diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) and scanning electron microscopy (SEM). We found that the concentration and formation rate of sulfate decreased with the increasing diameter of CaCO3 particles, and thus smaller particles could enhance the formation of sulfate due to more reactive sites on smaller particles. The rod-like calcium sulfate crystals were formed only at RH above 60% and in the presence of reactant gases through the heterogeneous pathway. The liquid-like water layer formed by promotion of high RH in the presence of reactant gases could facilitate the formation and aggregation of calcium sulfate hydrates and thus promote the formation and growth of rod-like calcium sulfate crystals. This study provides a possible mechanism for the formation and growth of rod-like calcium sulfate crystals existing in the atmosphere.


Introduction
Haze, with a high frequency of occurrence in China, exerts adverse impacts on the global climate and human health such as inducing many respiratory diseases. 1,2 Atmospheric sulfate presented an explosive growth during haze episodes and contributed signicantly to the formation of severe haze. [3][4][5] Furthermore, as one of the most important constituents of aerosols in atmosphere, the sulfate can serve as a cloud condensation nuclei (CCN) and scatter solar radiation, and thus affects the climate and weather. [6][7][8][9] Sulfate also can impact on reactivity, acidity and hygroscopicity of aerosols by changing the mixing state, components and morphology of particles. [10][11][12][13] Sulfate originates from primary sources including sea spray aerosol and volcanic ash and from secondary sources by oxidation of SO 2 , carbonyl sulfur (COS) or dimethyl sulde (DMS). 14,15 SO 2 with concentration up to several tens ppb in polluted days, is an important gaseous pollutant largely released by fossil fuel combustion. 16,17 Sulfate is produced mainly from the oxidation of SO 2 by OH radical in the gasphase 15 and by H 2 O 2 , O 3 and NO 2 in liquid phase. 18,19 The catalyzation of transition metal ions (for instance, Fe 3+ and Mn 2+ ) on the oxidation of SO 2 is another important pathway to form sulfate. 20,21 Heterogeneous oxidation of SO 2 by oxidants such as O 3 and NO 2 on mineral dust, sea salt or soot is proposed as a signicant source of sulfate. [22][23][24][25] These processes lead to the formation of different types of sulfate, for example, sodium sulfate, ammonium sulfate and calcium sulfate.
Calcium sulfate is a main component of atmospheric sulfate. 26 The strong correlation between SO 4 2À and Ca 2+ within aerosol particles has been revealed in eld measurements. 27,28 Also, a large fraction of calcium-containing particles is found to be calcium sulfate component, 29 which is a major type of sulfate in ne particles and composes almost total sulfate in coarse particles. 30 Rod-like morphology ranging from several hundred nanometers to approximately two micrometers are identied as a typical morphology of calcium sulfate in eld studies. [31][32][33][34] Until now, little attention has focused on the mechanism for the formation of rod-like calcium sulfate crystals. Rod-like calcium sulfate with diameters of 10 to 500 nm could be generated by the adsorption of SO 2 on particles at high RH due to the role of adsorbed water proposed by eld studies. 35 Li et al. also found that the formation of secondary regular particles like gypsum was promoted by high humidity. 36 From above, it could be speculated that the formation and growth of rod-like calcium sulfate crystals is inuenced by humidity. Furthermore, Gustafsson et al. found that the amount of water adsorbed on the mineral dust aerosols decreased with the increasing particle diameters. 37 Takahashi et al. suggested that the ratio of gypsum in total calcium-containing particles is larger for smaller particles, due to the more thorough reaction of SO 2 in smaller particles. 38 However, the key factors for formation and growth of the rod-like calcium sulfate crystals are still unclear in the atmosphere. CaCO 3 is proved to be one of the most important precursors of calcium sulfate in the atmosphere. 39,40 CaCO 3 comprises up to 30% of the total mineral dust, 41 and could be transported to long distance and wide range, 42 and thus it is oen a naturally available precursor. As one of reactive components of atmospheric mineral dust, 43 CaCO 3 is liable to be converted into calcium sulfate through reacting with trace acidic gaseous pollutants in the atmosphere. 44 Several studies have focused on the formation of calcium sulfate transformed from CaCO 3 . It was found that RH and temperature could inuence the formation of sulfate on CaCO 3 particles. [45][46][47] Ma et al. showed that a deliquescent layer is crucial for the formation of gypsum. 48 However, few studies have been conducted with respect to the evolution of morphology for the sulfate formed on the CaCO 3 . In this study, we aim at providing a mechanistic understanding on how rod-like calcium sulfate crystals are formed in the atmosphere and investigating whether RH and particle diameters of reactant particles could impact on the formation process. O 3 is an efficient oxidant and could accelerate the formation of sulfate by oxidizing sulte on particles, 49,50 hence is chosen as a model oxidant in this study. Therefore, we have undertaken an experimental study to explore the formation and growth of calcium sulfate crystals by oxidation of SO 2 by O 3 on size-resolved CaCO 3 at different RH. Products and kinetics were investigated by diffuse reectance infrared Fourier transform spectroscopy (DRIFTS). The evolution of the sizes and morphology of calcium sulfate was observed by scanning electron microscopy (SEM).

Sources of chemicals
Three kinds of CaCO 3 particles with different diameters were purchased from commercial sources (InnoChem, Beijing InnoChem Science & Technology Co., Ltd.). The shapes and diameters were observed by SEM and analyzed by image analysis soware (NIH, Imagej). The smallest CaCO 3 particles were approximately square, with diameter of 65 AE 10 nm determined by edge length. The other two larger particles were spindle shape, and the diameters were calculated to be 270 AE 50 nm and 1000 AE 300 nm, respectively, by the mean value of the short and long axes of the spindles. 51 The three kinds of particles were denoted as 65 nm, 270 nm and 1000 nm CaCO 3 particles, respectively, with the corresponding Brunauer-Emmett-Teller (BET) surface area of the samples measured to be 23.742 m 2 g À1 , 13.461 m 2 g À1 and 5.301 m 2 g À1 (Autosorb-1-MP automatic equipment (Quanta Chrome Instrument Co.)). CaCO 3 particles were heated for 3 h for outgassing before the BET measurement. N 2 (>99.999%, Beijing Tailong Electronics Co., Ltd.) and O 2 (>99.998%, Orient Center Gas Science & Technology Co., Ltd.) were used to simulate the ambient air, both of which were desiccated by silica gel and molecular sieve to ensure that RH was less than 1%. SO 2 (500 ppm, Beijing Huayuan Gas Chemical Industry Co. Ltd.) and ozone (generated by irradiating a dry ow of pure O 2 with an ultraviolet lamp (Heraeus Ltd., ZSZ-8)) were used as the reactant gases. Before each experiment, pure O 2 was passed through an ultraviolet lamp compartment of the ozone generator (UVP, Model SOG-2) for 1 h to achieve the steady concentration of ozone. The ozone concentration measured by an ozone analyzer (Model 49i, Thermo, USA) was controlled by varying the ux of oxygen exposed to ultraviolet light. Ultrapure water (18.2 MU cm) used for humidifying N 2 was obtained from Water Purication System (Dura FV).

Gas supply system
The gas supply system was composed of ve Teon inlet lines. The rst and second line supplied SO 2 and a mixture of ozone and oxygen, respectively. The concentrations of both SO 2 and ozone were 4.9 Â 10 14 molecules cm À3 in this study. The third line was used to supply humid N 2 produced by bubbling a dry ow of pure N 2 through ultrapure water. The fourth and h line provided dry N 2 and O 2 for dilution to simulate the ambient air. The mixing of all gases with a total ow of 400 sccm before entering the reaction chamber resulted in a synthetic air (20% O 2 and 80% N 2 ) with additions of SO 2 , O 3 and H 2 O. Mass ow controllers (Beijing Sevenstar electronics Co., Ltd.) were used to adjust the ux of the ve lines to get the desired concentration of SO 2 and O 3 , as well as the expected RH. A commercial humidity and temperature sensor (HMT330, Vaisala) was used to measure the RH and temperature at the inow of the sample cell with a measurement accuracy of AE 1% RH and AE 0.2 K, respectively.

Measurement methods
DRIFTS was employed to record the reaction process in situ on CaCO 3 particles. A Nicolet FTIR spectrometer 6700 was equipped with a liquid-nitrogen-cooled narrow band mercurycadmium-telluride (MCT) detector, Praying Mantis accessory (Model DRP, Harrick Scientic Corp.), and a high temperature reaction chamber (Model HVC-DRP-3, Harrick Scientic Corp.), used for recording the DRIFTS spectra in the spectral range of 650 to 4000 cm À1 at a resolution of 4 cm À1 . 100 scans were averaged for each spectrum corresponding to a time resolution of 40 s in this study. DRIFTS connected with the gas supply system has been described in detail elsewhere. [52][53][54] The reaction chamber was cleaned with alcohol and then evacuated under high vacuum over 1 h. Samples of about 30 mg, 37 mg and 45 mg were placed into the stainless steel cup (10 mm in diameter, 0.5 mm in depth) for 65 nm, 270 nm and 1 mm CaCO 3 , respectively. Aerwards, the sample cup was put into the reaction chamber, with temperature controlled by an automatic temperature controller (ATC-024-2, Harrick Scientic Corp.) combined with circulating cooling water. The sample was kept at 673 K for 1 h to remove adsorbed species from the surface and then cooled to 298 K for 1 h to ensure a stable condition. 55 The spectra of unreacted particles at a specic RH were recorded as background aer synthetic air was introduced into the reaction chamber for about 1 h to achieve water adsorption equilibrium on the samples. SO 2 and O 3 with specic concentrations were stabilized through the bypass line for about 30 min and 1 h, respectively, and then were introduced into the reaction chamber with an average residence time of 2.5 s. Spectra of products were shown as positive bands while losses of surface species as negative bands. All the spectral data was automatically collected by series program in OMINC soware during a typical experiment time of 120 min. Each experiment at the same RH was repeated at least three times, and the reproducibility of results was within acceptable error probably caused by articial operation such as packing samples in the sample cup. 46 S-4800 SEM (Hitachi, Japan) operated at 10 kV was employed to determine the morphology features and particle diameters of CaCO 3 particles before and aer the reaction. Several hundred micrograms samples were sprayed onto the surface of silicon wafers coated by gold. The silicon wafer was then put into the SEM chamber for observation of the unreacted particles aer heating at 673 K for 1 h. Three silicon wafers sprayed with the three size-resolved CaCO 3 particles were placed into the sample cup and reacted in the DRIFTS chamber under the same conditions as samples directly in the sample cup. Aer reaction, the silicon wafers coated by gold were used to determine the morphology features of particles by SEM.

Results and discussion
The surface products, the kinetics and morphology features during the heterogeneous reactions were analyzed based on the measurements by DRIFTS and SEM. Then, a mechanism was proposed based on the result analysis. Fig. 1 represented the in situ DRIFTS spectra of surface products when the size-resolved samples were exposed to SO 2 and O 3 simultaneously for 120 min at different RHs. The adsorption bands were almost the same for three size-resolved CaCO 3 at the same RH. The peaks in the range from 990 to 1270 cm À1 were ascribed to the stretching mode (n) of adsorbed sulfate (SO 4 2À ). 56 Under dry condition (1.5% RH), peaks at 985 cm À1 belonged to n(SO 3 2À ) in calcium sulte. 57 The peaks at 839 cm À1

Characterization of surface products
and 1689 cm À1 were assigned to the bending mode (d) and n(C]O) of CO 3 2À in carbonic acid adsorbed on the particle surface, respectively, indicating the formation of H 2 CO 3 . 58 Hence, the products were CaSO 4 , CaSO 3  and n(OH) of adsorbed water on the particle surfaces. 47 Thus, the products were CaSO 4 and adsorbed H 2 CO 3 at 19% RH. At 53% RH, peaks at 839 cm À1 and 1689 cm À1 also disappeared. Peaks at 1196 cm À1 shied to 1190 cm À1 , which indicated the formation of new sulfate products. New bands at 1629, 3552, and 3610 cm À1 were assigned to the d(OH) and n(OH) of water in bassanite (CaSO 4 $0.5H 2 O). 59 Hence, bassanite was formed at 53% RH. At 90% RH, peaks at 995 cm À1 and 1136 cm À1 shied to 1006 cm À1 and 1142 cm À1 , respectively, which belonged to n(SO 4 2À ) in calcium sulfate hydrates, and thus suggested that more calcium sulfate hydrates were formed on the CaCO 3 particles. 60 Double peaks of d(OH) at 1629 and 1685 cm À1 and three peaks of n(OH) at 3248, 3406, and 3552 cm À1 belonged to absorbance of water in CaSO 4 $2H 2 O (gypsum). 61,62 Therefore, bassanite and gypsum were formed at 90% RH in this study.

The sulfate formation kinetics on size-resolved CaCO 3 particles
The concentration of the ions formed during the reactions showed non-linear dependence on integrated area calculated by Kubelka-Munk function, while showed linear relationship with integrated absorbance. 54 Hence, integrated absorbance from 1020 to 1270 cm À1 as a function of reaction time was used to follow the concentration of sulfate ions during the reaction of SO 2 with O 3 on the surface of 65 nm, 270 nm and 1 mm CaCO 3 particles at various RH, respectively. The amount of sulfate formed during the reaction was calculated by DRIFTS calibration curve (seen in Fig. S1 †). The DRIFTS calibration curve was made by a series of known weight ratios of CaSO 4 and CaCO 3 mixtures, which showed that the amount of sulfate had a good linear relationship with the DRIFTS signal and could be expressed by: Here f is the conversion factor. The calculated f value was 6.83 Â 10 18 ions g À1 ABU À1 (ABU, absorbance unit) for three kinds of CaCO 3 particles, consistent with the order of magnitude of the values in previous studies. 45,46 Additional experiments proved that particle diameter and RH had no inuence on f value. Fig. 2 represents sulfate ion concentrations as a function of particle diameter. The concentrations of sulfate ions decreased with increasing particle diameter of CaCO 3 . The differences of sulfate ion concentrations for three size-resolved CaCO 3 were almost the same. It could be explained by the fact that the smaller particles with higher specic surface area and more active sites could adsorb more water and SO 2 , and thus concentration of sulfate decreased with increasing particle diameter of CaCO 3 particles.
The formation rate of sulfate d[SO 4 2À ]/dt can be quantied through the slope of sulfate concentration as a function of reaction time. The formation process of sulfate on CaCO 3 surfaces exhibited two stages. At the beginning of the reaction, the rapid increase of sulfate concentration was denoted as the initial stage, aerwards, the formation rate of sulfate decreased gradually until reached a stable rate aer a certain time, and the stable stage appeared. 63 The formation rate of sulfate ions decreased with increasing particle diameter in the two stages (seen in Fig. 3), consistent with the variation of concentration of sulfate ions as a function of particle diameter. In initial stage, the differences of sulfate formation rate for three size-resolved CaCO 3 were almost unchanged at 1.5% RH, 19% RH and 53% RH, and decreased at 90% RH. During the initial stage, SO 2 had no time to diffuse into the CaCO 3 samples, as a result, the formation rate of sulfate was mainly determined by surface reactive sites and the amount of adsorbed SO 2 . Since the surface adsorbed water could occupy the reactive sites, the surface reactive sites decreased with increasing RH. 65 nm CaCO 3 particles with higher specic surface area could adsorb more water than 270 nm and 1 mm CaCO 3 particles at the same RH and thus the surface reactive sites decreased much faster. At 90% RH, aer the water adsorption equilibrium was achieved, the le surface reactive sites were not enough for the adsorption of SO 2 , hence, the formation rate decreased for 65 nm CaCO 3 compared with that at 19% and 53% RH. Therefore, the differences of sulfate formation rate in initial stage among the three kinds of CaCO 3 decreased at 90% RH. In stable stage, the differences of sulfate formation rate for three size-resolved CaCO 3 were larger under humid condition than under dry condition. It has been proved that the amount of surface adsorbed water increased with enhanced RH, which could promote the adsorption of SO 2 and improve the ion mobility. 64 Consequently, more SO 2 could diffuse into the bulk of the samples under humid condition than dry condition and thus accelerated the formation of sulfate in the stable stage. Considering more water adsorbed on the surface of smaller CaCO 3 particles, it could be understood that the differences of formation rate in stable stage were larger under humid condition compared with that under dry condition. The differences of sulfate ion concentrations seemed to be determined by formation rate both in initial stage and stable stage. The reactive uptake coefficient g, was dened as the ratio of the rate of the reactive gases collisions with the surface to the rate of total gases surface collisions (Z). The g for the heterogeneous oxidation of SO 2 was deduced from DRIFTS experiments. On the basis of the assumption that the rate of reactive SO 2 collision with the surface was equal to the sulfate formation rate d[SO 4 2À ]/dt, the reactive uptake coefficients of SO 2 can be calculated by the following equations.
Here, Z is the rate of total collisions between SO 2 and particles, c represents the mean molecular velocity of SO 2 , A surface is the effective surface area of CaCO 3 particles, and [SO 2 ] is the gas-phase concentration of SO 2 . R is the molar gas constant, T refers to the temperature, and M SO 2 represents the Fig. 1 In situ DRIFTS spectra when SO 2 and O 3 were exposed simultaneously to 65 nm, 270 nm and 1 mm CaCO 3 particles, respectively, at (a) dry condition, 1.5% RH, (b) 19% RH, (c) 53% RH, (d) 90% RH. Absorption lines are recorded at 10 (black), 30 (olive), 60 (red), 90 (magenta) and 120 min (blue), respectively. molecular weight of SO 2 . If the reaction probability is high, sulfate was formed only on the surface of the samples since SO 2 has no time to diffuse into the CaCO 3 samples, and thus the geometric surface area of the samples would be the effective surface area. If the reaction probability is low, sulfate formed was distributed evenly through the whole samples since SO 2 could diffuse into the sample entirely, and the effective surface area would be the BET surface area. 65 Hence, the calculated uptake coefficients of SO 2 on CaCO 3 samples based on the BET and geometric surface area were considered as the lower and upper limits, respectively. 49 Uptake coefficients were denoted as g 0 and g ss in initial stage and stable stage, respectively. In this study, g 0 using geometric surface areas were determined to be (6.71 AE 0.55) Â 10 À4 , (6.41 AE 0.45) Â 10 À4 and (5.06 AE 0.55) Â 10 À4 for 65 nm, 270 nm and 1 mm CaCO 3 particles, respectively. These values were in the range of values determined under the same condition by Wu et al. (2011) and Li et al. (2006). On the basis of BET surface areas, our uptake coefficients g 0 were (7.51 AE 0.61) Â 10 À8 , (1.00 AE 0.07) Â 10 À7 and (1.66 AE 0.18) Â 10 À7 for 65 nm, 270 nm and 1 mm CaCO 3 particles, respectively. The values for 65 nm and 270 nm CaCO 3 particles were lower than g 0 determined under the same condition in previous studies. 45,46 It could be attributed to the fact that the BET surface areas of 65 nm and 270 nm CaCO 3 particles were much larger than that in literatures. 45,46 g 0 and g ss varied with diameters of CaCO 3 particles for adsorption of SO 2 at various RHs were listed in Tables 1 and 2, respectively. Due to different amount of 65 nm, 270 nm and 1 mm CaCO 3 particles packed in the sample cup, the g based on geometric surface areas did not show regular change with the diameters of CaCO 3 particles. In contrast to the variation tendency of sulfate formation rate with particle diameters, the g using BET surface areas increased with increasing diameters of CaCO 3 particles, due to the great differences of the BET surface areas among three size-resolved CaCO 3 particles.

The variation in morphology features of particles
The morphology features of particles before and aer the reaction observed by the SEM were shown in Fig. 4 and 5. It is interesting that the morphologies of sulfate products were almost the same for size-resolved CaCO 3 particles at the same RH. At about 50% RH, some microcrystallites appeared on the surface of the samples at the reaction time of 120 min, which were consistent with the appearance of bassanite based on in situ DRIFTS spectra. At 90% RH, rod-like calcium sulfate crystals with several hundred nanometers in length and approximately 100 nm in width occurred as the reaction lasted 120 min, corresponding to the formation of bassanite and gypsum based on the IR results. Aerwards, more rods with length ranging from several hundred nanometers (reaction lasted 200 min) to over one micrometers (reaction lasted 400 min) were observed with the extended reaction time. These rods grew on and stretched out of the samples and covered the surface of particles. The length of rods was no more than 2 mm even aer the reaction time of 24 h under 90% RH condition (Fig. 5(2)) accompanied with appearance of some microcrystallites. According to the results of DRIFTS and SEM, when bassanite appeared without the presence of gypsum at 50% RH, microcrystallites were observed without rods. At higher RH, when gypsum was observed with certain amount, the rods also appeared. Hence, it could be speculated that microcrystallites and rods were largely composed of bassanite and gypsum, respectively. Fig. 2 The concentration of sulfate ions formed on particle surfaces during the reaction at 1.5% RH (olive), 19% RH (red), 53% RH (magenta) and 90% RH (blue) and reaction time of 120 min as a function of CaCO 3 particle diameter. ]/dt (a) in initial stage and (b) in stable stage at 1.5% RH (olive), 19% RH (red), 53% RH (magenta) and 90% RH (blue) as a function of CaCO 3 particle diameter.
To conrm the formation conditions of rod-like calcium sulfate crystals, further experiments were conducted. First, when the reaction lasted 24 h at 50% RH (Fig. 5(1)), only microcrystallites were observed. Second, the reaction time was 120 min at 50% RH, and then the RH was elevated to 90% for about 24 h in the absence of SO 2 and O 3 (Fig. 5(3)). Only microcrystallites appeared, and no rods could be observed. Third, at the reaction time of 120 min and 60% RH, rod-like calcium sulfate crystals with sizes of 50-100 nm in length and several tens nanometers in width appeared along with the microcrystallites similar to observed at 50% RH (Fig. S2 †). It indicated that even more microcrystallites were formed, they could not aggregate together to form rod-like crystals at 50% RH. Furthermore, microcrystallites formed at 50% RH also could not coagulate to form rod-like crystals even the RH increased to 90%. It could be speculated that presence of reaction gases at RH above 60% was necessary for the formation of rod-like calcium sulfate crystals. The comparisons of size of calcium sulfate crystals formed at 60% and 90% RH indicated that sizes of rod-like calcium sulfate crystals increased with increasing RH above 60%.

Proposed mechanisms
A mechanism for the formation and growth of rod-like calcium sulfate crystals by the heterogeneous reactions of SO 2 with O 3 on size-resolved CaCO 3 particles at different RHs was proposed. Under dry condition, CaSO 4 was formed on the surface of CaCO 3 particles upon adsorption of SO 2 and oxidation by O 3 . Under humid condition, adsorbed water on calcium carbonate surface could facilitate the formation of Ca(OH)(HCO 3 ), which was reactive towards SO 2 and thus promoted the adsorption of SO 2 . 64 Hence, the amount of adsorbed SO 2 increased with elevated RH and thus accelerated the formation of sulfate. In contrast, when SO 2 existed in the reaction systems, the particle surface became more hydrophilic upon reaction, 47 and thus promoted further adsorption of water. When the amount of adsorbed SO 2 and water were both enough, calcium sulfate hydrates could be formed immediately with the reaction. Hence, bassanite was formed at 50% RH and resulted in the Table 1 The reactive uptake coefficient in initial stage, g 0 when SO 2 and O 3 were exposed simultaneously to size-resolved CaCO 3 particles at different RH  Table 2 The reactive uptake coefficient in stable stage, g ss when SO 2 and O 3 were exposed simultaneously to size-resolved CaCO 3 particles at different RH  (1) Unreacted particles, (2 and 3) after SO 2 (4.9 Â 10 14 molecules cm À3 ) and O 3 (4.9 Â 10 14 molecules cm À3 ) were exposed simultaneously for 120 min at 50% RH and 90% RH, respectively.
formation of calcium sulfate microcrystallites on the surface of CaCO 3 particles. At higher RH, gypsum was also formed by the heterogeneous reaction process. Above 50% RH, the liquid-like water layer was formed by multilayer coverages of surface adsorbed water accompanied with the promotion effect of SO 2 , which could facilitate the formation and aggregation of calcium sulfate hydrates during the reactions, and thus led to the formation of rod-like calcium sulfate crystals. Therefore, it could be understood that the high RH (above 60%) and the presence of reactant gases is requirement for the formation and growth of rod-like calcium sulfate crystals. Upon the aggregation of sulfate on the particle surface, fresh reactive sites of CaCO 3 were exposed and thus further reinforced the formation of sulfate. The sizes and amount of rodlike calcium sulfate crystals increased with RH above 60% and increased with reaction time below 400 min. When the whole surface of CaCO 3 have been covered by calcium sulfate crystals as shown in Fig. 5(2a), the reactive sites were prevented from exposing to reactant gases, and thus further formation of sulfate was inhibited in the deeper position of the samples. As a result, the limiting size of rod-like calcium sulfate crystals was approximately 2 mm in length during the heterogeneous reactions.

Conclusions
This study provides experimental evidence for heterogeneous formation pathway of rod-like calcium sulfate crystals in the atmosphere. The key factors for the pathway to form rod-like calcium sulfate crystals are RH above 60% and the presence of reactant gases. We show that liquid-like water layer formed by promotion of high RH in the presence of reactant gases could promote the formation and growth of rod-like calcium sulfate crystals by facilitating the formation and aggregation of calcium sulfate hydrates. Reactant particle diameters could affect the formation of sulfate. Results here imply that concentration and formation rate of sulfate decreased with increasing diameter of CaCO 3 particles. Along with the longer suspending time in the air, CaCO 3 particles with smaller diameters could accordingly contribute greater to sulfate formation.
The rod-like shape is a typical morphology of calcium sulfate crystals observed in most eld studies, which is consistent with the morphology of calcium sulfate crystals acquired by laboratory heterogeneous reactions in this study. Thus, the heterogeneous reaction pathway is a possible source of rod-like calcium sulfate crystals in the atmosphere. Calcium sulfate crystals with different morphology also existed in the atmosphere. Guo et al. found that calcium sulfate derived from coal y ash were euhedral crystals or globular melt droplets. 66 Pan et al. showed that calcium sulfate emitted from coal-red power plant aer wet ue gas desulfurization were prismatic or clintheriform shape with size distribution centered at approximately 0.2 mm. 67 It could be speculated that the calcium sulfate observed in eld study may not be originated from the direct emission of coal-red power plant or coal y ash. Therefore, it is a feasible method to recognize different origins or formation pathways of sulfate by the morphology of the particles in the atmosphere. The present study has examined only the morphology of calcium sulfate formed by heterogeneous reactions. Further laboratory studies on characterizing morphology features of sulfate formed through other formation pathways are needed to better understand the origin of sulfate.

Conflicts of interest
There are no conicts to declare.