Asian Journal of atmospheric environment
Asian Journal of atmospheric environment Asian Journal of atmospheric environment
Asian Journal of atmospheric environment
  Aims and Scope Type of Manuscripts Best Practices Contact Information  
  Editor-in-Chief Associate Editors Editorial Advisory Board  

Current Issue

Asian Journal of Atmospheric Environment - Vol. 16 , No. 4

[ Technical Information ]
Asian Journal of Atmospheric Environment - Vol. 14, No. 2
Abbreviation: Asian J. Atmos. Environ
ISSN: 1976-6912 (Print) 2287-1160 (Online)
Print publication date 30 Jun 2020
Received 25 Sep 2019 Revised 30 Jan 2020 Accepted 07 Apr 2020
DOI: https://doi.org/10.5572/ajae.2020.14.2.155

Estimating Mass Concentration Using a Low-cost Portable Particle Counter Based on Full-year Observations: Issues to Obtain Reliable Atmospheric PM2.5 Data
Sayako Ueda1), 2), * ; Kazuo Osada1) ; Makiko Yamagami3) ; Fumikazu Ikemori3) ; Kunihiro Hisatsune3)
1)Graduate School of Environmental Studies, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
2)Japan Society for the Promotion of Science, Kojimachi, Chiyoda-ku, Tokyo 102-0083, Japan
3)Nagoya City Institute for Environmental Sciences, 5-16-8 Toyoda, Minami-ku, Nagoya, Aichi 457-0841, Japan

Correspondence to : *Tel: +81-52-789-4306 E-mail: ueda.sayako@d.mbox.nagoya-u.ac.jp


Copyright © 2020 by Asian Association for Atmospheric Environment
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Funding Information ▼

Abstract

Expanding the use of a recently introduced low-cost particle monitor (DC1700 Dylos Air Quality Monitor) for sensing atmospheric PM2.5 requires comparison with data obtained using a certified method for PM2.5 based on appropriate atmospheric observations. Full-year measurements of atmospheric aerosols were taken in Nagoya, Japan during March 2017-March 2018 using the DC1700 to measure the particle number concentrations of >0.5 and >2.5 μm diameter particles and to measure the PM2.5 mass concentration (Mdry, PM2.5) using an automated β attenuation mass monitor (PM712). The number-size distribution was measured using an optical particle counter (KC01D). The dried mass concentration of 0.5-2.5 μm particles (Mdry, 0.5-2.5) was estimated from the ambient relative humidity and the DC1700 number concentration. The values of Mdry, 0.5-2.5 were invariably less than those of Mdry, PM2.5. The coefficient of determination and slope of Mdry, 0.5-2.5 to Mdry, PM2.5 for the year were, respectively, 0.68 and 0.40. Slope values changed seasonally from 0.24 in July and August 2017 to 0.55 in May and April 2017. Light absorbing particles, smaller-fine particles, and the estimation method of Mdry, 0.5-2.5 were inferred as causes of the difference between Mdry, 0.5-2.5 and Mdry, PM2.5. Especially, we estimated a large contribution (ca. 54% underestimation of Mdry, 0.5-2.5 into Mdry, PM2.5) of particles smaller than the minimum detection diameter of DC1700. The seasonal variation of Mdry, 0.5-2.5/Mdry, PM2.5 was related to the volume fraction of particles smaller than 0.5 μm. Good correlation of Mdry, 0.5-2.5 to Mdry, PM2.5 suggests that data obtained using DC1700 with a correction factor are useful as a rough proxy of atmospheric PM2.5 within a season. However, precise estimation of PM2.5 from the DC1700 number concentrations should include appropriate corrections of the size distribution, not only hygroscopicity.


Keywords: Atmospheric aerosols, Dylos, Low-cost monitor, Optical sensor, PM2.5

1. INTRODUCTION

Atmospheric aerosol particles have been studied extensively because of their important role in assessing air quality and Earth climate. Fine particulate matter has been recognized as adversely affecting human health, leading to premature mortality in many people (Zhang et al., 2017; Lelieveld et al., 2015). Among fine particles, those with aerodynamic diameter less than 2.5 μm are defined as PM2.5. Actually, PM2.5 mass concentrations have been monitored widely using various instruments based on tapered-element oscillating microbalances, beta attenuation, and a hybrid of beta attenuation and light scattering, in conjunction with an impactor or a cyclone at the inlet (EPA, 2013; Kulkarni et al., 2011).

Recently, low-cost monitors of PM2.5 have received increasing attention for measuring temporal and spatial concentration variations in indoor and outdoor environments (Rai et al., 2017; Jovašević-Stojanović et al., 2015; Kumar et al., 2015). Particularly developed have been instruments with optical sensors that detect light scattered from particles. Their use has expanded rapidly because of their low cost, compact size, and high time resolution (Liu et al., 2017; Nakayama et al., 2017; Zikova et al., 2017; Jiao et al., 2016; Austin et al., 2015). Although some points of caution have been indicated, such as the need for correction under high humidity conditions ( Jayaratne et al., 2018; Zheng et al., 2018; Han et al., 2017), many reports have described good correlation between output values from these devices and PM2.5 mass concentrations measured using conventional methods ( Johnson et al., 2018; Kelly et al., 2017; Nakayama et al., 2017; Rai et al., 2017; Zikova et al., 2017). Such low-cost monitors, because of their beneficial characteristics of mobility and time-resolution, offer great potential for use in many situations such as investigation of streets and building environments.

Results obtained using the two-channel (>0.5 μm and >2.5 μm) optical particle counters (DC1100 and DC1700; Dylos Corp.) adopted for this study have been compared widely with PM2.5 mass concentrations measured in atmospheric, laboratory, and closed environments (Rai et al., 2017; Jones et al., 2016; Manikonda et al., 2016; Sousan et al., 2016; Dacunto et al., 2015; Semple et al., 2013). Earlier reports have described good correlation between PM2.5 mass concentration and outputs from DC series. The device is sold at a reasonable price. Therefore, we expect that the device can be used easily for measuring PM2.5 at various sites to provide multipoint and high-time-resolution data. Nevertheless, the light scattering intensity differs depending on the optical properties of aerosol species such as dust, sea salt, and diesel fumes (Sousan et al., 2016). Atmospheric aerosols have various size distributions and compositions, which differ by time and location. Such aerosol physicochemical properties can affect their optical properties. Therefore, detection sensitivities of optical sensors vary according to time and location. Using data from a particle counter as a proxy of atmospheric PM2.5 therefore requires knowledge of the reliability and limitations of measurements. Although one report has described comparison with PM2.5 mass concentration (Han et al., 2017), the observation periods used for that study were less than 1 month. Moreover, they covered few seasons. A shorter observation duration makes it difficult to detect variation according to various parameters. Therefore, realistic evaluation of optical responses is difficult. For another low-cost sensor PMS, Sayahi et al. (2019) and Zheng et al. (2018) evaluated the performance for measuring PM2.5 mass concentrations based on long-term observations. They reported seasonal differences in sensor performance. Zheng et al. (2018) compared the sensor performance to meteorological parameters, especially relative humidity. However, the relation of the performance with seasonal/event changes of atmospheric aerosol properties such as light absorbing aerosols and size distribution was not explained. For these reasons, data obtained using a low-cost monitor should be compared for at least one year to those obtained from a certified PM2.5 monitor with aerosol properties.

For this study, simultaneous measurements obtained using one optical particle counter (DC1700; Dylos Corp.) were used to compile particle number concentrations of >0.5 and >2.5 μm diameter particles for comparison with PM2.5 data obtained based on β attenuation method during one year in Nagoya, Japan. Most earlier studies using data obtained using a DC1700 had evaluated the correlation coefficient and regression formula for PM2.5 mass concentrations directly from the number concentration (Rai et al., 2017; Manikonda et al., 2016; Sousan et al., 2016; Dacunto et al., 2015; Semple et al., 2013). However, the degree of undercounting or overcounting by DC1700 was less apparent when using the earlier method. In the present study, to evaluate underestimation and overestimation of PM2.5 and temporal variation more quantitatively, the mass concentration of particles counted by the DC1700 was estimated and compared to the PM2.5 mass concentration. Several studies have demonstrated the necessity of considering aerosol hygroscopicity. Therefore, we clarified the long-term difference originated in other factors by estimation of the mass concentration after humidity correction. Special attention is devoted to factors affecting uncertainty derived from temporal changes in the number-size distribution.


2. MEASUREMENTS AND DATA TREATMENTS
2. 1 Observation Sites and Instruments

Table 1 presents locations and instrumentation used for this study. Simultaneous observations were made using the DC1700 device and standard monitor of PM2.5 at the Nagoya City Institute for Environmental Sciences (NCIES) located in Nagoya city, central Japan, during March 2017-March 2018. Fig. 1 portrays a location map of the observation sites in Nagoya, a major city with about 2.3 million residents. Yamagami et al. (2019) reported details of NCIES and PM2.5 measurements at the site. The air mass at the site is often affected by long-range transport from continental eastern Asia in addition to domestic emissions. Particle diameters detected using the DC1700 show separation at >0.5 μm and >2.5 μm. The DC1700 was put in a weather shield box fitted with air ducts passing outside air easily. The weather shield box was installed outdoors under eaves to prevent sunlight and rain exposure. Hourly PM2.5 mass concentrations were monitored using a β attenuation mass monitor (PM712; Kimoto Electric Co. Ltd.) installed in an outdoor case on the flat roof of the NCIES building. The PM712 is a PM2.5 automatic measuring instrument that uses the standard method accepted in Japan. The temperature and relative humidity of the sample air and ambient air were measured respectively downstream of a PTFE tape filter and outside of the outdoor case. The PM712 also measured the mass concentration of optical black carbon (OBC, MOBC) based on attenuation of near-IR scattering on a sampling spot of PM2.5. Daily samples of PM2.5 were collected using a pair of FRM-2000 samplers (Rupprecht and Patashnick Co. Inc., Albany, NY, USA) with PTFE filters (TK15-G3M; Pall Corp., Port Washington, NY, USA) for ion components and quartz fiber filters (2500QATUP; Pall Corp.) for carbon at NCIES. These samples were collected from 10:00 a.m. through 9:30 a.m. the next day. Details of FRM samplers at NCIES were reported by Ueda et al. (2016) and by Yamagami et al. (2019). Although original data obtained using the DC1700 were recorded every minute, their hourly averages were used for comparison with PM2.5 mass concentrations.

Table 1. 
Locations and instruments discussed.
Station Nagoya City Institute for
Environmental Sciences (NCIES)
Nagoya University (NU) National air monitoring station
in Nagoya at Kanokoden (KNK)
Location 35.10°N, 136.92°E,
0 m a.s.l.
35.16°N, 136.97°E,
49 m a.s.l.
35.18°N, 136.98°E,
59 m a.s.l.
Direction and
distance
from NCIES
- Northeast, about 8 km Northeast, about 9 km
Instrument DC1700
 Number concentrations of
 >0.5, >2.5 μm
DC1700
 Number concentrations of
 0.5, >2.5 μm
ACSA
 Mass concentration of PM2.5
PM712
 Mass concentration of PM2.5
 OBC mass concentration of PM2.5
 RH and Temp. of ambient and sample air
KC01D
 Number concentrations of
 >0.3, >0.5, >1.0, >2.0, >5.0 μm
FRM-2000
 Mass concentration of PM2.5
Period Mar. 2017-Mar. 2018 Mar. 2017-Mar. 2018 Apr. 2017-Mar. 2018


Fig. 1. 
Map showing locations of observation sites (Nagoya City Institute for Environmental Sciences (NCIES), Nagoya University (NU), National air monitoring station in Nagoya at Kanokoden (KNK), and the local meteorological observatory (Met. Obs.)) in Nagoya, Japan.

The number-size distribution of aerosol particles, which was measured continually at Nagoya University (NU), was also used to evaluate data at NCIES. The distance between NCIES and NU is about 8 km (Fig. 1). The number-size distributions were measured using an optical particle counter (KC01D; Rion Co. Ltd.) with a laser diode. The KC01D measured the number concentrations for particles larger than 0.3, 0.5, 1.0, 2.0, and 5.0 μm diameter. The irradiated volume v of the laser was 5.0×10-4 cm3. Sample air was introduced into the KC01D located in the observation room of the seventh floor after drying in a diffusion dryer to less than 20% relative humidity (RH). The length of tubing from the inlet to KC01D was about 3 m. The sample air flow rate was 0.5 L min-1. A diffusion dryer used to keep the sample RH below 20% consisted of a 20-cm-long, 7-cm-diameter acrylic pipe and two desiccants mainly consisting of magnesium chloride (Nisso Dry-M; Nisso Fine Co., Ltd.). The KC01D was replaced by another instrument of the same type for maintenance during 12 May 2017 and 8 December 2017 because of air flow clogging. The use of KC01D data was limited to evaluation of the relative change of particle size distribution. Data of the optical particle counter needed to be corrected for the loss for coincidence counting error, depending on v and the particle number concentration. For comparison to data at the NCIES, measurements using the DC1700 in a weather shield box were also conducted at NU. In addition, PM2.5 mass concentrations were monitored at the national air monitoring station in Nagoya at Kanokoden (KNK), which is near NU (within 2 km). Those data were used to evaluate regional differences of PM2.5 mass concentrations. From March 2017, the PM2.5 mass concentrations were measured hourly (ACSA-14; Kimoto Electric Co. Ltd.) using the same method as that used for PM712.

Hourly precipitation and temperature data were obtained from a local meteorological observatory at Nagoya. The distance from NCIES to KNK and the meteorological observatory is about 9 km (Fig. 1). The meteorological observatory is also within 2 km of NU.

2. 2 PM2.5 Mass Concentration by PM712 and FRM

Mass concentrations of PM2.5 used for routine air monitoring such as the Federal Reference Method (FRM), which is widely regarded as reliable, were measured under the condition of 35±5% RH around the sample. PM712 provides the mass concentration of PM2.5 equivalent to the value of FRM after correcting the water contents of ambient particles (Kimoto Electric Co. Ltd., 2007).

Mtotal,PM2.5Mdry,PM2.5=Mdry,PM2.5+Mw,PM2.5Mdry,PM2.5=1+a×expRHsample100×b(1) 

Therein, Mtotal, PM2.5, Mdry, PM2.5, and Mw, PM2.5 respectively represent the mass concentrations of PM2.5 at the relative humidity of the sample air (RHsample), the mass concentration of dried PM2.5, and the mass concentration of water contained in PM2.5 at RHsample. In addition, a and b (respectively 0.010 and 6.000) are coefficients ascertained from correlations of daily data between the FRM and the PM712 obtained in Japan. A correction curve based on equation (1) for RH<90% is similar to a theoretical equation, assuming hygroscopicity of ammonium sulfate and ammonium nitrate (Iwamoto et al., 2018; Snider et al., 2016).

For equivalence checking of PM2.5 mass concentration in this study, the PM2.5 mass concentration found using PM712 at NCIES was compared with that estimated from daily sample of PM2.5 by FRM. Fig. 2 presents scatter plots of daily PM2.5 mass concentrations obtained using FRM and PM712 during the observation period. Daily PM2.5 mass concentrations of PM712 were averaged for 24 hr of data obtained during 10:00 a.m. to 10:00 a.m of the next day, using hourly data after correction by equation (1). Although plots for PM2.5 mass concentration before humidity correction (Mtotal, PM2.5) show dispersion, the relation between PM2.5 mass concentration by FRM and Mdry, PM2.5 showed good correlation of R2=0.95 and slope=1.01. The greatest outlier plot (arrowed plot in Fig. 2) was that of data from a Kosa (dust) event of 7 May 2017. However, such cases were rarely observed. Most of the PM712 data showed a nearly 1:1 relation with those of FRM, suggesting that the daily average of mass concentration of PM2.5 by PM712 after humidity correction using equation (1) was usually reliable to that by FRM. For this study, the hourly PM2.5 mass concentration by PM712 was referenced to compare data of DC1700. The manufacturer states that averaged and standard deviations of hourly blanks of PM712 are within ±3.5 μg m-3 for more than 15 data. A recent report has described that PM2.5 mass concentration by some instruments using β attenuation can have different bias between night and daytime (Hasegewa et al., 2018). However, their test for 12-hr samples using FRM showed good correspondence with the PM2.5 mass concentration of PM712 in many cases, except at some times during daytime in winter.


Fig. 2. 
Scatter plots of mass concentrations of PM2.5 by PM712 and FRM at NCIES for the year of Mar. 2017 through Mar. 2018.

2. 3 Estimation of Mass Concentration from DC1700 Data

Based on the difference between particle number concentrations of >0.5 μm and >2.5 μm diameter obtained using DC1700, the number concentrations (N0.5-2.5) for 0.5 μm to 2.5 μm particles are calculable. The mass concentration (M0.5-2.5) for sizes of 0.5 μm and 2.5 μm was estimated from N0.5-2.5, given as

M0.5-2.5=ρπ6D3N0.5-2.5(2) 

where ρ and D respectively stand for the particle density and particle diameter. The D value was calculated as the geometric average between 0.5 μm and 2.5 μm. The ρ value was assumed for density (1.8 g/cm3) of ammonium sulfate as the most abundant compound in PM2.5 at Nagoya (Yamagami et al., 2019; Ueda et al., 2016; Ikemori et al., 2015). According to their reports, sulfate, organic carbon, ammonium, elemental carbon, and nitrate were the major components, comprising more than 65% of the PM2.5 mass concentration for all seasons. Chloride, sodium, potassium, magnesium, and calcium were less than 5% in PM2.5. Among the other densities of materials considered from major components, the respective densities of sulfuric acid and ammonium nitrate are 1.84 and 1.73 g/cm3, which are similar values to those found for ammonium sulfate. The density of carbonaceous particles (approx. 1.5 g/cm3; Slowik et al., 2004) is slightly lower than that of ammonium sulfate. Therefore, the value of Mdry, 0.5-2.5 using the density of ammonium sulfate might be overestimated by as much as 17% under the condition that all particles are carbonaceous particles.

The value of M0.5-2.5 from DC1700 must be corrected for the hygroscopic increase of particle mass because measurements are conducted under ambient conditions without humidity control. For this study, the dried mass concentration (Mdry, 0.5-2.5) obtained using the DC1700 was estimated, given as

Mdry,0.5-2.5=M0.5-2.51+a×expRHambient100×b-1(3) 

where RHambient represents the relative humidity of ambient air. This correction uses the same method for PM712 as that shown in equation (1). Instead of RHsample of equation (1), RHambient was applied to the correction for DC1700 because the length to the optical sensor in DC 1700 is short.


3. RESULTS AND DISCUSSION
3. 1 Temporal Variation of Mdry, 0.5-2.5 and Mdry, PM2.5

Fig. 3 portrays temporal variations in March 2017 as an example of well-synchronized variation of Mdry, PM2.5 and Mdry, 0.5-2.5: (a) the number concentration of aerosol particles measured using DC1700 at NCIES and NU; (b) Mdry, 0.5-2.5, Mdry, PM2.5, and MPM2.5 by FRM and Mdry, 0.5-2.5/Mdry, PM2.5 at NCIES; (c) MOBC at NCIES; (d) the size-segregated volume concentration (0.3-5.0 μm) by KC01D at NU; (e) the relative humidity of ambient air measured at outside of PM712 (RHambient) and downstream of the sample filter of PM712 (RHsample); and (f) the temperature and precipitation amounts recorded at the local meteorological observatory. The volume concentrations were calculated from the number concentrations for the respective size ranges of the KC01D, given as shown below.

VDi-Di+1=π6DiDi+132Ni-Ni+1(4) 

Fig. 3. 
Temporal variations are shown: (a) number concentrations found using DC1700 at NCIES and NU; (b) mass concentrations of 0.5-2.5 μm diameter particles found using DC1700 (Mdry, 0.5-2.5), mass concentration of PM2.5 by PM712 (Mdry, PM2.5), and by FRM (MPM2.5 FRM), and the ratio of Mdry, 0.5-2.5 and Mdry, PM2.5 at NCIES; (c) mass concentration of OBC (MOBC) found using PM712; (d) volume concentrations of aerosols with 0.3-5 μm diameter by OPC at NU; (e) relative humidity of ambient air measured at outside of PM712 (RHambient) and downstream of the sample filter of PM712 (RHsample) at NCIES; and (f) temperature and hourly precipitation at a local meteorological observatory for March 2017.

Therein, Di and Ni respectively stand for the minimum particle diameter and number concentration of channel i counting for particles larger than Di of KC01D.

As portrayed in Fig. 3a, the temporal variations of the number concentrations obtained using DC1700 at NU and NCIES were almost equal. Such synchronous variation was observed throughout the observation period. Fig. 4 presents a scatter plot of the number concentrations (>0.5 μm and >2.5 μm) obtained by DC1700 at NCIES and NU and scatter plots of PM2.5 mass concentrations at NCIES and KNK for the observation periods. Their respective coefficients of determination (R2) and slope values of linear regression were 0.78-0.87 and 0.91-1.13. The good degree of correlation suggests that the air masses of these locations can be regarded as almost equivalent. As portrayed in Fig. 3b, the value of Mdry, PM2.5 was almost equal to that of MPM2.5 by FRM. The temporal variation of Mdry, 0.5-2.5 correlated well with that of Mdry, PM2.5 (Fig. 3b) and the volume concentration (0.3-5.0 μm) obtained by KC01D (Fig. 3d). However, the values of Mdry, 0.5-2.5 were always less than that of Mdry, PM2.5: about 50% of Mdry, PM2.5 in March 2017. Also, RHambient was usually less than 80%, as depicted in Fig. 3e. The period of RHambient higher than 80% in the period was recorded during precipitation events (Fig. 3f). During such events, Mdry, PM2.5 tended to be lower. Higher Mdry, PM2.5 (>20 μg m-3) events were observed mostly under conditions of less than 70% RHambient. Although ambient and sample air humidity changed diurnally, Mdry, 0.5-2.5/Mdry, PM2.5 (Fig. 3b) is less dependent on humidity. The values of MOBC (Fig. 3c) were usually less than 10% of Mdry, PM2.5.


Fig. 4. 
Scatter plots of number concentrations of >0.5 μm and 2.5 μm measured using DC1700 at NCIES and NU, and of the mass concentrations of PM2.5 at NCIES and KNK.

3. 2 Relation between Mdry, 0.5-2.5 and Mdry, PM2.5

Scatter plots of Mdry, PM2.5 and Mdry, 0.5-2.5 at NCIES for two months (a, March-April 2017; g, March 2018) and for all data in one plot (h) are presented in Fig. 5. The Mdry, 0.5-2.5 values are shown below the 1 : 1 line to Mdry, PM2.5 values. The coefficient of determination (R2) and slope values of linear regression for all data in this study were, respectively, 0.68 and 0.40. Although the observation values for the whole period (Fig. 5h) showed large dispersion, the relations within two months (as shown in Figs. 5a, b, e and g) were found to have stronger correlation. For the scatter plot showing data of January and February (Fig. 5f ), greater variation was apparent at higher concentrations during particular events in February (red dots in Figs. 5f and h), as discussed in section 3.4. The respective maximum and minimum values of R2 were 0.90 found for March-April 2017 (Fig. 5a) and 0.62 found for July-August (Fig. 5c) 2017. The respective maximum and minimum values of slopes were also 0.55 found for March-April 2017 and 0.24 found for July-August 2017. The slope values of linear regression differed by season. Fig. 6 portrays monthly box plots for ratios of Mdry, 0.5-2.5/Mdry, PM2.5 at NCIES. For March-April 2017 and December 2017-March 2018, the median values of Mdry, 0.5-2.5/Mdry, PM2.5 were 0.4-0.5. However, the median values were less than 0.4 for May-November 2017. They were especially low (approx. 0.25) in summer ( July-August). The ratio of Mdry, 0.5-2.5/Mdry, PM2.5 decreased once in summer but recovered in winter and the following spring.


Fig. 5. 
Scatter plots of mass concentrations of 0.5-2.5 μm in diameter by DC1700 versus mass concentration of PM2.5 by PM712 for two months (a-g) and for a year from Mar. 2017 through Mar. 2018. Periods of events B (18:00 19-8:00 20 February 2018), and C (21:00 254-9:00 26 February 2018) and D (16:00-20:00 28 February 2018) are denoted by red dots.


Fig. 6. 
Monthly box plots for ratio of mass concentration of 0.5-2.5 μm by DC1700 (Mdry, 0.5-2.5) to mass concentration of PM2.5 by PM712. The lower boundary of the box shows the 25th percentile, the line within the box represents the median, and the upper boundary of the box shows the 75th percentile. Whiskers above and below the box respectively show the 90th and 10th percentiles. Numbers shown above boxes are the numbers of observed data.

Han et al. (2017) also measured ambient aerosols using DC1700 in Houston, Texas. Correlation between the number concentration of 0.5-2.5 μm and PM2.5 mass concentration found using standard instruments for 12 days was a similar level (R2=0.778). Seasonal changes in detection for another low-cost sensor to measure atmospheric PM2.5 also reported from results of several studies. Sayahi et al. (2019) measured mass concentrations of particulate matter using PMS 1003 and 5003 at Salt Lake City, Utah, in the USA. Correlation to 24-hr averaged PM2.5 mass concentration by FRM was good (R2>0.88) in winter, but poorer (R2 of 0.18-0.32) in spring. Nakayama et al. (2017) evaluated a new palm-sized optical PM2.5 sensor (PM2.5 sensor; Panasonic Corp.) at four urban and suburban sites in Japan. Their slope to standard PM2.5 mass concentration tended to be higher in winter than in summer, which was inferred as attributable to small particles by photochemical formation. They also reported that the sensor tended to overestimate the PM2.5 mass concentration by hygroscopic growth of particles.

3. 3 Cause of Difference between Mdry, 0.5-2.5 and Mdry, PM2.5 and Seasonal Variation of Mdry, 0.5-2.5/Mdry, PM2.5

The Mdry, 0.5-2.5 found in this study was corrected for hygroscopic growth. Therefore, differences between Mdry, 0.5-2.5 and Mdry, PM2.5 reflect the contribution of other factors. The Mdry, 0.5-2.5 values were always less than those of Mdry, PM2.5. Sousan et al. (2016) reported tests of the detection efficiency of DC1700 using monodispersed fine salt (0.1, 0.2 and 0.3 μm) and large oleic acid (1.3, 2, 3 and 5 μm) particles. The actual detection efficiency for the DC1700 in the size range of 0.5-2.5 μm was 52% for 1.3 μm particles. The difference of actual detection efficiency from manufacturer specifications can affect the underestimation of Mdry, 0.5-2.5 to Mdry, PM2.5. Moreover, the difference of Mdry, 0.5-2.5 to Mdry, PM2.5 in this study tended to vary seasonally. This finding suggests that the difference can change along with the seasonal change of atmospheric aerosols. As factors affecting the seasonal change of Mdry, 0.5-2.5/Mdry, PM2.5, we can infer the existence of light-absorbing particles and smaller fine particles below the detection limit, and infer a difference of size distribution in the estimation equations to that of actual atmospheric aerosols. Effects on estimation of the mass concentration of PM2.5 using DC1700 and the seasonal change by each factor are explained below.

3. 3. 1 Existence of Light-absorbing Particles

Light-absorbing particles are detected inappropriately as being of smaller size because of the lower intensity of light scattering for the same physical size of non-light-absorbing particles. According to Sousan et al. (2016), who described laboratory tests for DC 1700, the detected diesel exhaust particles were one or more orders fewer than those detected using non-light absorbing salts having the same mass concentration. In the present study, the average fraction of OBC mass concentrations to PM2.5 was 7.8±9.5% (average±standard deviation for this observation period). The light absorption of aerosols might partly affect underestimation at such a level. Fig. 7a, b and c present monthly box plots for Mdry, PM2.5, MOBC and ratios of MOBC/Mdry, PM2.5 at NCIES. Median values of MOBC/Mdry, PM2.5 were less than 0.2 for all data (Fig. 7c). The ratios were slightly high in October with slight variation. In that year, Mdry, PM2.5 of October was lowest (Fig. 7a) because of the arrival of multiple typhoons. MOBC of October was of a comparable level to the other months (Fig. 7b). Therefore, the mass fraction of black carbon was apparently high. However, comparison to Mdry, 0.5-2.5/Mdry, PM2.5 in Fig. 6 shows that the seasonal variation of MOBC/Mdry, PM2.5 was less and was poorly matched.


Fig. 7. 
Monthly box plots for mass concentrations of PM2.5 (a) and OBC (b) and the ratio of mass concentration of OBC to PM2.5 (c) by PM712 at NCIES and the ratio of volume concentrations of 0.5-2 μm to 0.3-2 μm (d) measured by OPC at NU. The lower boundary of the box shows the 25th percentile, the line within the box represents the median, and the upper boundary of the box shows the 75th percentile. Whiskers above and below the box respectively show the 90th and 10th percentiles. The numbers shown above boxes are the numbers of observed data.

3. 3. 2 Existence of Smaller Fine Particles below the Detection Limit

The minimum size detectable by DC1700 is 0.5 μm as the optical diameter. Earlier studies have also demonstrated the significant contribution of particles with diameters smaller than the detectable size limit of the DC 1700 sensor (Zikkova et al., 2017; Manikoda et al., 2016; Sousan et al., 2016). Fig. 7d presents monthly box plots for the volume ratios of 0.5-2.0 μm to 0.3-2.0 μm (V0.5-2.0 μm/V0.3-2.0 μm) by KC01D at NU. The median values were maximum (0.5) in December and minimum (0.3) in July. The seasonal variation of high values in winter and low values in summer were close to those of Mdry, 0.5-2.5/Mdry, PM2.5 in Fig. 6. The volume fraction of particles smaller than 0.5 μm in PM2.5 might affect the seasonal difference between Mdry, 0.5-2.5 and Mdry, PM2.5. Fig. 8 portrays the monthly averaged number-size distribution measured using the KC01D device at NU and their volume-size distribution for December and July 2017. For the volume-size distribution, the area of each bar for a size range corresponds to the volume concentration for the size range. Aerosol measurements were taken using the DC1700 device under ambient relative humidity. Therefore, some hygroscopic particles having less than 0.5 μm dry diameter can be counted if the particles become larger than 0.5 μm under ambient conditions. The dry size (Ddry, 0.5) of the minimum detectable size (0.5 μm) by the DC1700 in ambient conditions was estimated using the hygroscopicity of ammonium sulfate (Snider et al., 2016). Actually, Ddry, min was 0.39±0.06 μm, on average, during the observation periods. Similarly, the dry size (Ddry, 2.5) of the upper detection size (2.5 μm) obtained using the DC1700 device was estimated as 1.9±0.3 μm. The Ddry, 0.5 and Ddry, 2.5 are shown as dashed and solid red lines in Fig. 8b.


Fig. 8. 
Monthly averaged number-size distribution (a) and volume-size distributions (b) of aerosol particles by KC01D for July and December 2017. Red dash-dotted and solid lines of (b) respectively represent the yearly averaged dried size of 0.5 μm and 2.5 μm particles estimated by assuming ammonium sulfate.

PM2.5 particles by PM712 were classified as less than 2.5 μm under ambient conditions by a virtual impactor. Therefore, the size range of dried PM2.5 was also estimated as less than Ddry, 2.5 (ca. 2 μm). Although particles less than Ddry, 0.5 can affect underestimation of Mdry, 0.5-2.5 to Mdry, PM2.5, the estimated volume concentration for sizes of 0.3-2 μm in Fig. 8b was a maximum at 0.3-0.5 μm. The Ddry, 0.5 was in 0.3-0.5 μm. According to parameters for atmospheric volume-size distributions derived from averages of measurements reported by Whitby (1978), the respective mode peaks of the accumulation and coarse mode particles were 0.32 μm and 5.7 μm for urban aerosols. The volume fraction of particles smaller than 0.5 μm to particles smaller than 2.5 μm can be estimated as 62% based on Whitby’s parameters. Considering size changes attributable to hygroscopic growth, the volume fraction of particles smaller than 0.39 μm (Ddry, min) among particles smaller than 2 μm was 54%. The result suggests that the contribution of particles smaller than the detection limit especially accounts for the underestimation of Mdry, 0.5-2.5 to Mdry, PM2.5.

Although the volume concentration of particles smaller than 0.3 μm was unknown for this study, the volume ratio of 0.3-0.5 μm to 0.5-2 μm of July 2017 was minimum (Fig. 7d). The volume fraction of particles smaller than the detection limit (Ddry, 0.5) was expected to be high in July 2017; consequently, the values of Mdry, 0.5-2.5/Mdry, PM2.5 were small. Given a volume concentration of less than 0.3 μm as V<0.3 μm, the ratio of Mdry, 0.5-2.5 to Mdry, PM2.5 was almost equivalent to the ratio of V0.5-2.0 μm to (V<0.3 μm+V0.3-2.0 μm). The ratio of Mdry, 0.5-2.5 to Mdry, PM2.5 is expected to be somewhat larger because the DC1700 counted hygroscopic particles of less than 0.5 μm. In consideration of general and observed volume-size distribution, we assumed simple values of V<0.3 μm: the values of V<0.3 μm are V0.3-2.0 μm and 0.5·V0.3-2.0 μm, respectively, when the accumulation mode diameters are small (as in July) and large (as in December). Then, the ratios of V0.5-2.0 μm to (V<0.3 μm+V0.3-2.0 μm) were estimated respectively as 0.15 and 0.33. Consequently, the size change can become a factor of difference of 0.18 by this assumption, which was comparable to a difference of Mdry, 0.5-2.5/Mdry, PM2.5 between July and December (ca. 0.2) in Fig. 6. As implied by this estimate and by Whitby’s size distribution, the seasonal variation of V0.5-2.0 μm/V0.3-2.0 μm suggests that a particle volume smaller than the detection limit of DC1700 varies with the season. The contributions of smaller particles might be greater in summer because the formation and growth of secondary particles is more active in summer (Okada, 1985). These results suggest that simple conversion of number data from Dylos to PM2.5, which is independent of the size distribution, would engender large bias, with variation derived from seasonal changes in the size distribution.

3. 3. 3 Difference of Size Distribution in the Estimation Equations to Actual Atmospheric Aerosols

From estimation of Mdry, 0.5-2.5 using equations (2) and (3), the geometric average between 0.5 and 2.5 μm was used as D in this study. In the equations, a dV/dlogD was used for the size range. However, the equation is correct in theory when dV/dlogD is linear to logD because excess volume for a size range larger than the geometric average between Ddry, 0.5 and Ddry, 2.5 is equal to the shortage volume for a smaller size range. For the observed volume-size distribution as presented in Fig. 8b, the volume concentration around Ddry, 0.5 tended to be somewhat high, especially in July 2017, rather than showing a linear change. Therefore, the assumption of D in this study (i.e. a geometric average between 0.5 and 2.5 μm) can engender some overestimation of the actual Mdry, 0.5-2.5. Considering that point, the mass concentration of dried particles counted by DC1700 would be smaller than the estimated Mdry, 0.5-2.5, especially when the volume concentration of approximately 0.5 μm was high. Therefore, the difference of size distribution between assumption and actual aerosols is regarded as reducing the seasonal difference of Mdry, 0.5-2.5/Mdry, PM2.5, rather than the factor of the difference.

3. 4 Example of Short-term Event Effects on a Particular Size Distribution

Fig. 9 presents temporal variations of Mdry, 0.5-2.5, Mdry, PM2.5 at NCIES, Mdry, PM2.5 at KNK, and other parameters by KC01D, along with meteorological data for February 2018. As described earlier, Mdry, 0.5-2.5/Mdry, PM2.5 was usually about 0.4. However, Mdry, 0.5-2.5/Mdry, PM2.5 were occasionally higher than 0.7, with high mass concentration events such as B (18:00 19-8:00 20 February 2018), C (21:00 24-9:00 26 February 2018), and D (16:00-20:00 28 February 2018). As depicted in Fig. 5f and h with data colored in red, these high ratios of Mdry, 0.5-2.5/Mdry, PM2.5 were rarely observed during the remainder of the year. For the other high mass concentration event A (8 February 2018) in Fig. 9a, the ratio of Mdry, 0.5-2.5/Mdry, PM2.5 was at a normal level (ca. 0.5). Values of MPM2.5 by FRM and Mdry, PM2.5 at KNK during events A-D also corresponded well to Mdry, PM2.5 at NCIES (Fig. 9a). These increases of Mdry, PM2.5 in events A-D occurred during the dry condition (RHambient<70%) without precipitation (Fig. 9a, e and f).


Fig. 9. 
Temporal variations are shown: (a) mass concentration of 0.5-2.5 μm in diameter by DC1700 (Mdry, 0.5-2.5) and mass concentrations of PM2.5 by PM712 and FRM (Mdry, PM2.5 at NCIES and MPM2.5 FRM) at NCIES and mass concentrations of PM2.5 by ASCA-14 at Kanokoden (Mdry, PM2.5 at KNK); (b) mass concentration of OBC; (c) volume concentrations of aerosols with 0.3-5 μm diameter by OPC at NU; (d) volume fractions of 0.3-0.5 μm, 0.5-1 μm and 1-2 μm particles to 0.3-2 μm particles; (e) relative humidity of ambient air (RHambient) and downstream of sample filter of PM712 (RHsample) at NCIES; and (f) temperature and hourly precipitation at local meteorological observatory for February 2018. Periods of high-concentration events A, B, C, and D are highlighted by green and yellow bars.

As estimated in section 3.3, particles smaller than the detection limit can account for the large fraction of underestimation of Mdry, 0.5-2.5 to Mdry, PM2.5 among some factors. Therefore, size distribution differences might affect the ratio of Mdry, 0.5-2.5 to Mdry, PM2.5. In events C-D, the fraction of 0.5-1.0 μm particles (Fig. 9d) increased with higher Mdry, PM2.5. In event A (8 February 2018), an increase of smaller (0.3-0.5 μm) particles contributed to the increase of Mdry, PM2.5. For comparison, volume-size distributions of events A and C are presented in Fig. 10. As the figure shows, the volume concentration of 0.5-1.0 μm particles of event C was clearly greater than that of event A. Therefore, the volume fraction of particles larger than Ddry, 0.5 was apparently higher for event C. The precipitous change observed for the difference between Mdry, 0.5-2.5 and Mdry, PM2.5 was regarded as mainly reflecting rapid changes in the volume-size distribution.


Fig. 10. 
Volume-size distributions of aerosol particles by KC01D for events A and C. Red dash-dotted and solid lines respectively show the dried sizes of 0.5 μm and 2.5 μm particles estimated by assuming ammonium sulfate.

As shown in Fig. 5h, drastic events of Mdry, 0.5-2.5/Mdry, PM2.5 change, such as events B, C, and D, rarely occurred throughout the year. Nevertheless, these events demonstrated that a particle fraction that is undetectable by DC1700 can change drastically in a single day, not only on the time-scale of seasons. Therefore, to use DC1700 data as a proxy of atmospheric PM2.5, attention must be devoted to temporal changes of the number-size distribution.


4. SUMMARY AND CONCLUSIONS

To assess the usability of a low-cost optical particle counter for measurements of atmospheric PM2.5, measurements were taken using DC1700, PM712, and FRM for a year in Nagoya city, Japan. The daily average of Mdry, PM2.5 by PM712 was well correlated with the mass concentration of PM2.5 by FRM. The value of Mdry, 0.5-2.5 measured using DC1700 was always less than Mdry, PM2.5 measured by PM712, but they showed good correlation (slope=0.40, R2=0.68) for the year. The slope values of Mdry, 0.5-2.5 to Mdry, PM2.5 tended to change seasonally: the minimum was 0.24 of July-August 2017; the maximum was 0.55 of May-April 2017.

As causes of difference between Mdry, 0.5-2.5 and Mdry, PM2.5, we considered some factors in addition to the detection efficiency of instruments evaluated by an earlier study: light absorbing particles, smaller fine particles with less than the minimum detection diameter, and the method of estimating Mdry, 0.5-2.5. The contributions of light absorption of aerosols and particles less than the minimum detection diameter were inferred respectively as ca. 8% (on average during this observation) and ca. 54% (in general). The seasonal variation of Mdry, 0.5-2.5/Mdry, PM2.5 apparently reflects the seasonal variation of the number-size distribution better than that of light absorbing particles. Three short-term events observed during the year changed the ratio of Mdry, 0.5-2.5 to Mdry, PM2.5 to high values with the increase of Mdry, PM2.5. The sporadic event affecting Mdry, 0.5-2.5 was attributed to the increase of the detectable particle fraction with the change in the number-size distribution.

The mass concentration estimated from DC1700 showed some degree of correlation with PM2.5 for each of two months. Our results suggest that data conversion using a simple correction factor can provide a rough value of PM2.5 within a season. However, as this study demonstrated, the mass fraction of particles of undetectable size found using DC1700 was large in atmospheric PM2.5. This outcome can engender underestimation or overestimation of PM2.5 from DC1700 data according to differences in the size distribution. Consequently, to estimate PM2.5 more reliably, appropriate correction of temporal changes in the size distribution should be included not only to adjust for hygroscopicity.


Acknowledgments

This work was performed with the support of a Grantin-Aid from the Environment Research and Technology Development Fund, Grant No. 5-1604, provided by Environmental Restoration and Conservation Agency of Japan.


References
1. Austin, E., Novosselov, I., Setom, E., Yost, M.G. (2015) Laboratory Evaluation of the Shinyei PPD42NS Low-Cost Particulate Matter Sensor. PLOS ONE, 10(9).
2. Dacunto, P.J., Klepeis, N.E., Cheng, K.-C., Acevedo-Bolton, V., Jiang, R.-T., Repace, J.L., Ott, W.R., Hildemann, L.M. (2015) Determining PM2.5 Calibration Curves for a Low-Cost Particle Monitor: Common Indoor Residential Aerosols. Environmental Science: Processes & Impacts, 17, 1959-1966.
3. EPA (2013) National Ambient Air Quality Standards for Particulate Matter; Final Rule, Federal Register of Environmental Protection Agency, USA, 78, 3086-3287.
4. Han, I., Symanski, E., Stock, T.H. (2017) Feasibility of using low-cost portable particle monitors for measurement of fine and coarse particulate matter in urban ambient air. Journal of the Air & Waste Management Association, 67(3), 330-340.
5. Hasegawa, S., Yamagami, M., Suzuki, Y., Kumagai, K., Nishimura, R. (2018) Verification of measured values by PM2.5 automatic measuring instruments using the standard method. Journal of Environmental Laboratories Association, 43(1), 40-46 (in Japanese). http://tenbou.nies.go.jp/science/institute/region/journal/JELA_4301040_2018.pdf
6. Ikemori, F., Honjyo, K., Yamagami, M., Nakamura, T. (2015) Influence of contemporary carbon originating from the 2003 Siberian forest fire on organic carbon in PM2.5 in Nagoya, Japan. Science of Total Environment, 530-531, 403-441.
7. Iwamoto, Y., Sekine, H., Saito, S., Miura, K., Nishikawa, M., Nagano, K., Osada, K. (2018) Continuous measurement of hygroscopic characteristics of PM2.5 using an optical particle counter - Including case analysis of high concentration events in December 2016. Earosoru Kenkyu, 33(4), 238-247.
8. Jayaratne, R., Liu, X., Thai, P., Dunbabin, M., Morawska, L. (2018) The influence of humidity on the performance of a low-cost air particle mass sensor and the effect of atmospheric fog. Atmospheric Measurement Techniques, 11(8), 4883-4890.
9. Jiao, W., Hagler, G., Williams, R., Sharpe, R., Brown, R., Garver, D., Judge, R., Caudill, M., Rickard, J., Davis, M., Weinstock, L., Zimmer-Dauphinee, S., Buckley, K. (2016) Community Air Sensor Network (CAIRSENSE) project: evaluation of low-cost sensor performance in a suburban environment in the southeastern United States. Atmospheric Measurement Techniques, 9(11), 5281-5292.
10. Johnson, K., Bergin, M.H., Russell, A.G., Hagler, G.S.W. (2018) Field Test of Several Low-Cost Particulate Matter Sensors in High and Low Concentration Urban Environments. Aerosol and Air Quality Research, 18, 565-578.
11. Jones, S., Renée Anthony, T., Sousan, S., Altmaier, R., Park, J.H., Peters, T.M. (2016) Evaluation of a Low-Cost Aerosol Sensor to Assess Dust Concentrations in a Swine Building. The Annal of Occupational Hygiene, 60(5), 597-607.
12. Jovašević-Stojanović, M., Bartonova, A., Topalović, D., Lazović, I., Pokrić, B., Ristovski, Z. (2015) On the use of small and cheaper sensors and devices for indicative citizen-based monitoring of respirable particulate matter. Environmental Pollution, 206, 696-704.
13. Kelly, K.E., Whitaker, J., Petty, A., Widmer, C., Dybwad, A., Sleeth, D., Martin, R., Butterfield, A. (2017) Ambient and laboratory evaluation of a low-cost particulate matter sensor. Environmental Pollution, 221, 491-500.
14. Kulkarni, P., Baron, P.A., Willeke, K. (2011) Aerosol Measurement: Principles, Techniques, and Applications, Third ed., John Wiley & Sons, Inc. Hoboken, NJ, 883pp.
15. Kumar, P., Morawska, L., Martani, C., Biskos, G., Neophytou, M., Di Sabatino, S., Bell, M., Norford, L., Britter, R. (2015) The rise of low-cost sensing for managing air pollution in cities. Environment International, 75, 199-205.
16. Lelieveld, J., Evans, J.S., Fnais, M., Giannadaki, D., Pozzer, A. (2015) The Contribution of Outdoor Air Pollution Sources to Premature Mortality on a Global Scale. Nature, 525, 367-371. https://www.nature.com/articles/nature15371
17. Liu, D., Zhang, Q., Jiang, J., Chen, D.-R. (2017) Performance calibration of low-cost and portable particular matter (PM) sensors. Journal of Aerosol Science, 112, 1-10.
18. Manikonda, A., Zíková, N., Hopke, P.K., Ferro, A.R. (2016) Laboratory assessment of low-cost PM monitors. Journal of Aerosol Science, 102, 29-40.
19. Nakayama, T., Matsumi, Y., Kawahito, K., Watabe, Y. (2017) Development and evaluation of a palm-sized optical PM2.5 sensor. Aerosol Science and Technology, 52(1), 2-12.
20. Rai, A.C., Kumar, P., Pilla, F., Skouloudis, A.N., Sabatino, S.D., Ratti, C., Yasar, A., Rickerby, D. (2017) End-user perspective of low-cost sensors for outdoor air pollution monitoring. Science of The Total Environment, 607-608, 691-705.
21. Sayahi, T., Butterfield, A., Kelly, K.E. (2019) Long-term field evaluation of the Plantower PMS low-cost particulate matter sensors. Environmental Pollution, 245, 932-940.
22. Semple, S., Ibrahim, A.E., Apsley, A., Steiner, M., Turner, S. (2013) Using a new, low-cost air quality sensor to quantify second-hand smoke (SHS) levels in homes. Tobacco Control, 24(2).
23. Slowik, J.G., Stainken, K., Davidovits, P., Williams, L.R., Jayne, J.T., Kolb, C.E., Worsnop, D.R., Rudich, Y., DeCarlo, P.F., Jimenez, J.L. (2004) Particle Morphology and Density Characterization by Combined Mobility and Aerodynamic Diameter Measurements. Part 2: Application to Combustion-Generated Soot Aerosols as a Function of Fuel Equivalence Ratio. Aerosol Science and Technology, 38(12), 1206-1222.
24. Snider, G., Weagle, C.L., Murdymootoo, K.K., Ring, A., Ritchie, Y., Stone, E., Walsh, A., Akoshile, C., Anh, N.X., Balasubramanian, R., Brook, J., Qonitan, F.D., Dong, J., Griffith, D., He, K., Holben, B.N., Kahn, R., Lagrosas, N., Lestari, P., Ma, Z., Misra, A., Norford, L.K., Quel, E.J., Salam, A., Schichtel, B., Segev, L., Teipathi, S., Wang, C., Yu, C., Zhang, Q., Zhang, Y., Brauer, M., Cohen, A., Gibson, M.D., Liu, Y., Martins, J.V., Rudich, Y., Martin, R.V. (2016) Variation in global chemical composition of PM2.5: emerging results from SPARTAN. Atmospheric Chemistry and Physics, 16, 9629-9653.
25. Sousan, S., Koehler, K., Thomas, G., Park, J.H., Hillman, M., Halterman, A., Peters, T.M. (2016) Inter-comparison of lowcost sensors for measuring the mass concentration of occupational aerosols. Aerosol Science and Technology, 50(5), 462-473.
26. Ueda, K., Yamagami, M., Ikemori, F., Hisatsune, K., Nitta, H. (2016) Associations Between Fine Particulate Matter Components and Daily Mortality in Nagoya, Japan. Journal of Epidemiology, 26(5), 249-257.
27. Yamagami, M., Ikemori, F., Nakashima, H., Hisatsune, K., Osada, K. (2019) Decreasing trend of elemental carbon concentration with changes in major sources at Mega city Nagoya, Central Japan. Atmospheric Environment, 199, 155-163.
28. Zhang, Q., Jiang, X., Tong, D., Davis, S.J., Zhao, H., Geng, G., Feng, T., Zheng, B., Lu, Z., Streets, D.G., Ni, R., Brauer, M., van Donkelaar, A., Martin, R.V., Huo, H., Liu, Z., Pan, D., Kan, H., Yan, Y., Lin, J., He, K., Guan, D. (2017) Transboundary Health Impacts of Transported Global Air Pollution and International Trade. Nature, 543, 705-709. https://www.nature.com/articles/nature21712
29. Zheng, T., Bergin, M.H., Johnson, K.K., Tripathi, S.N., Shirodkar, S., Landis, M.S., Sutaria, R., Carlson, D.E. (2018) Field evaluation of low-cost particulate matter sensors in high-and low-concentration environments. Atmospheric Measurement Technology, 11, 4823-4846.
30. Zikova, N., Masiol, M., Chalupa, D.C., Rich, D.Q., Ferro, A.R., Hopke, P.K. (2017) Estimating hourly concentrations of PM2.5 across a metropolitan area using low-cost particle monitors. Sensors, 17(8), 1922.