the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
![](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-avatar-thumb150.png)
Spatial analysis of PM2.5 using a concentration similarity index applied to air quality sensor networks
Rósín Byrne
John C. Wenger
Stig Hellebust
Air quality sensor (AQS) networks are useful for mapping PM2.5 (particles with a diameter of 2.5 µm or smaller) in urban environments, but quantitative assessment of the observed spatial and temporal variation is currently underdeveloped. This study introduces a new metric – the concentration similarity index (CSI) – to facilitate a quantitative and time-averaged comparison of the concentration–time profiles of PM2.5 measured by each sensor within an air quality sensor network. Following development on a dataset with minimal unexplained variation and robust tests, the CSI function is used to represent an unbiased and fair depiction of the air quality variation within an area covered by a monitoring network. The measurement data is used to derive a CSI value for every combination of sensor pairs in the network, yielding valuable information on spatial variation in PM2.5. This new method is applied to two separate AQS networks, in Dungarvan and in the city of Cork, Ireland. In Dungarvan there was a lower mean CSI value (, ), indicating lower overall similarity between locations in the network. In both networks, the average diurnal plots for each sensor exhibit an evening peak in PM2.5 concentration due to emissions from residential solid-fuel burning; however, there is considerable variation in the size of this peak. Clustering techniques applied to the CSI matrices identify two different location types in each network; locations in central or residential areas that experience more pollution from solid-fuel burning and locations on the edge of the urban areas that experience cleaner air. The difference in mean PM2.5 between these two location types was 6 µg m−3 in Dungarvan and 2 µg m−3 in Cork. Furthermore, the examination of winter and summer months (January and May) indicates that higher PM2.5 levels during periods of increased residential solid-fuel burning act as a major driver for greater differences (lower similarity indices) between locations in both networks, with differences in mean seasonal CSI values exceeding 0.25 and differences in mean seasonal PM2.5 exceeding 7 µg m−3. These findings underscore the importance of including wintertime PM data in analyses, as the differences between locations is enhanced during periods when solid-fuel burning activities are at a peak. Additionally, the CSI method facilitates the assessment of the representativeness of the PM2.5 measured at regulatory air quality monitoring locations with respect to population exposure, showing here that location type is more important than physical proximity in terms of similarity and spatial representativeness assessments. Applying the CSI in this manner can allow for the placement of monitoring infrastructure to be optimised. The results indicate that the population exposure to PM2.5 in Dungarvan is moderately represented () by the current regulatory monitoring location, and the regulatory monitoring location assessed in Cork represented the city-wide PM2.5 levels well ().
- Article
(7383 KB) - Full-text XML
-
Supplement
(1093 KB) - BibTeX
- EndNote
Air pollution affects the environment and quality of life and is a major cause of premature death and disease (Cesaroni et al., 2013; Lelieveld et al., 2015; Pedersen et al., 2013; Raaschou-Nielsen et al., 2013). The category of air pollutant with the largest impact on human mortality and health is fine particulate matter, i.e. atmospheric particles with an aerodynamic diameter of 2.5 µm or less (PM2.5) (Pope et al., 2020; Pope and Dockery, 2012; Samoli et al., 2013). In many regions around the world, air quality monitoring and management have become critical endeavours to mitigate the detrimental effects of air pollution, and especially PM2.5, on citizens and the environment.
Over the years, technological advances have provided valuable tools to enhance our understanding of air pollution, and low-cost air quality sensors (AQSs) are emerging as promising instruments for collecting real-time air quality data at an improved spatial and temporal resolutions (Kumar et al., 2015; Munir et al., 2019). When used in networks, air quality sensors offer immense potential for enhancing and supplementing regulatory monitoring and assessment (Malings et al., 2020). However, further work needs to be carried out to assess the effectiveness of sensor networks and how to make best use of the data for gaining further insights into air pollution within a locality because the data quality obtained with such low-cost devices does not meet the standards for regulatory monitoring. Careful consideration must be given to the quality of the data provided by sensors, and the requirement for calibration must be assessed (Diez et al., 2022). Recent studies have shown that the performance and calibration of a PM2.5 sensor is dependent on the type of sensor and often on the measurement location, suggesting the need for site-specific and individual calibrations to correct for the absolute level of PM2.5 (Kaur and Kelly, 2023; Sayahi et al., 2019; Wang et al., 2015; Zamora et al., 2020). When these factors are considered and accounted for, AQS networks offer an unprecedented opportunity to gain further insights into the complex dynamics of air pollution in localised areas, such as urban environments, industrial zones, and residential neighbourhoods (Crawford et al., 2021; Frederickson et al., 2022; Heimann et al., 2015; Hodoli et al., 2023; O'Regan et al., 2022).
Information on the spatial variation in air quality is important because air pollution is not homogenous and can exhibit significant variations across different areas even on a local scale (Frederickson et al., 2022, 2023; Kassomenos et al., 2014; Wang et al., 2018). The variability in air pollution can be influenced by a multitude of factors such as traffic patterns, industrial activities, meteorological conditions, and local topography. Consequently, relying on single monitoring locations or limited data resolution can provide an incomplete picture and inadequate understanding of local air quality in a certain area (Li et al., 2019). Understanding these variations is crucial for targeted interventions and policy decisions aimed at improving air quality and safeguarding public health. Spatial analysis facilitated by sensor networks allows for a more accurate and nuanced understanding of how air quality, and therefore exposure to pollution, varies across a population centre.
In a recent study, we used data collected by a PM2.5 sensor network in the city of Cork, Ireland, to estimate the contribution of local pollution sources as separate and distinct from regional or transported air pollution (Byrne et al., 2023). The results highlighted the very localised nature of PM2.5 caused by residential solid-fuel burning during winter, which is a significant problem in many towns and cities in Ireland and elsewhere (Dall'Osto et al., 2013; Kourtchev et al., 2011; Lin et al., 2018, 2019; Ovadnevaite et al., 2021; Wenger et al., 2020; Zhang et al., 2021).
In this work, we propose a new approach for assessing the spatial profile of air quality using an AQS network. The method yields a time-averaged concentration similarity index (CSI) for quantitative assessment of the similarity between the complete data series produced by different sensors within the network. The CSI is built on the premise that sensors exposed to similar ambient conditions and pollutant sources will produce comparable PM2.5 temporal trends. Conversely, sensors subject to different conditions might display divergent PM2.5 concentration trends. The motivation for the development of an assessment method based on the temporal variation over an extended period is the realisation that the annual average is often an incomplete representation of true population exposure, which is experienced from hour to hour and day to day. If hourly or daily PM2.5 variability is high, it is not always adequate to merely compare annual averages of PM2.5 in different locations to compare PM2.5 exposure experienced by the local populations. While the annual average and hourly/daily values are often well correlated, numerous studies have found positive associations between short-term exposure to particulate matter and increased morbidity and mortality due to respiratory and cardiovascular diseases (Fajersztajn et al., 2017; Orellano et al., 2020; Weinmayr et al., 2010). This method aims to translate this idea into a quantifiable metric by calculating the time-averaged degree of similarity between two sensor datasets. After method development and testing, the CSI analysis is applied to an AQS network in the town of Dungarvan in Ireland to identify areas that may be experiencing persistently elevated or very localised PM2.5 pollution compared to others. Clustering techniques are used to group sensors based on the similarity of their PM2.5 measurements. The CSI method is also retrospectively applied to sensor network data collected in the city of Cork to investigate the transferability of the method between sensor networks and to explore any differences between the locations.
2.1 Data collection, preprocessing, and calibration
The collection, preprocessing, and calibration of the data collected by the PM2.5 sensor networks in Dungarvan and the city of Cork were carried out using the Julia programming language and the openair package written for the R programming language (Bezanson et al., 2017; Carslaw and Ropkins, 2012). Since low-cost AQSs are not of regulatory standard, great care needs to be taken with quality assessment and quality control of the data. In particular, the degree to which changes or differences in PM2.5 measurements between devices can be trusted needs to be considered. The methodology proposed here addresses these inherent issues to deliver an approach for assessing the spatial representativeness of any monitoring location and to facilitate comparison of different environment types, regardless of geographical distance to the location.
2.1.1 Dungarvan PM2.5 sensor network
The Dungarvan sensor network consisted of 18 solar powered Clarity Node-S devices (Clarity Movement Co., USA), which utilise the Plantower PM6003 sensor to measure PM2.5 within the range of 1–1000 µg m−3 and at a resolution of 1 µg m−3 (Clarity Movement Co., 2023; Node-S technical sheet, 2023). By default, the Node-S devices take measurements every 15 min, allowing sufficient data upload and battery sleep time in between sampling periods. However, this can be adjusted to higher or lower frequencies. The highest sampling interval achievable during winter without significantly affecting the battery performance was 8 min.
The Clarity Node-S devices were typically attached to street light poles between 2 and 4 m above the ground. The sensors were positioned in a range of different environments including urban background, residential, coastal, and roadside locations (Fig. S1 in the Supplement). Many of these locations were a mix of the different environments. The majority of devices were operational from 1 November 2022 to 31 May 2023; however three devices (AP7, AY9N, AY93) with the Clarity Wind Module were only deployed from 12 January 2023. Measurements were taken over a continuous period covering different meteorological seasons (mainly winter and spring/early summer), thus ensuring that temporal variations in PM2.5 concentrations were captured comprehensively.
Prior to and after deployment in Dungarvan, the Clarity Node-S devices were co-located on the roof of the Ellen Hutchins Building, University College Cork (51.895136, −8.516146), to compare their performance. Details of the three co-location periods are outlined in Table S1 in the Supplement. Although some devices were not available for all three co-location periods, the three periods combined provide a comparison between the sensors across different seasons. This co-location dataset enabled the CSI method to be developed on measurements that in theory should be equal, and the function could then be modified if necessary to allow for sensor behaviour, uncertainties, errors, and potential limitations.
The raw sensor data from the co-location periods and field deployment underwent a series of preprocessing steps to mitigate potential sources of error in the measurement and to ensure data quality and consistency. Data points outside of the operational range of the sensors (>1000 µg m−3) were identified and removed, although instances of these were minimal. The 8 min data were averaged to produce hourly measurements. Missing data points could potentially affect the temporal continuity of the data; however, the data coverage was overall very good for the co-location and measurement campaign periods. On average, the devices had an hourly measurement coverage of 87 % for the field measurement campaign. This corresponds to an average of 4443 hourly measurements per device for the campaign period.
Assessing the consistency of measurements across the sensor network was paramount. Although the PM2.5 readings were very well correlated when the devices were co-located (Table S2 in the Supplement), a data harmonisation procedure was performed to ensure the uniformity of sensor measurements, which is a prerequisite for the subsequent development of the concentration similarity index. Since there was no reference-grade PM2.5 data available during the co-location periods, the PM2.5 concentrations from each sensor were scaled to a common reference point, represented by the mean of all data points across the whole co-location dataset (Fig. S2 in the Supplement). The data series for each sensor was then individually compared with the calculated mean dataset and subsequently harmonised to the common reference point using a simple linear regression approach. The equations resulting from this harmonisation procedure were applied to the measurements collected from all devices during the subsequent field measurement campaign. While this procedure did not convert the measured PM2.5 to reference-equivalent concentrations, it minimised sensor output variability and facilitated a more equitable comparison between sensor measurements (Table S2).
2.1.2 Cork PM2.5 sensor network
For brevity when referring to the sensor network, the term “Cork” will herein refer exclusively to the city area of Cork. The Cork sensor network consisted of 16 PurpleAir PA-II-SD units that each contain two Plantower PMS5003 sensors to measure PM2.5 within the effective range of 0–500 µg m−3, with a maximum range of 1000 µg m−3, and at a resolution of 1 µg m−3 (PMS5003 series data manual, 2022). In this study, data recorded by the devices in the network for the periods from 1 January 2021 to 31 May 2021 and from 1 September 2021 to 31 December 2021 were collated and analysed. However, four devices were found to have limited data capture for the specified periods (<50 %) and were therefore omitted from the analysis. The 12 sensors used in this analysis had an average data capture of 85 % for the specified periods; their locations are shown in Fig. S3 in the Supplement.
Due to logistical constraints, it was not possible to co-locate all of the PurpleAir devices together to assess variability in PM2.5 concentrations. However, low inter-sensor and inter-unit variability was exhibited by four co-located PurpleAir devices in our previous study on the Cork network, where all inter-sensor and inter-unit comparisons yielded R2 values greater than 0.98 (Byrne et al., 2023). Moreover, PurpleAir PM2.5 measurements were highly correlated (R2=0.92) with hourly values of PM2.5 concentrations obtained using a Met-One (USA) beta-attenuation monitor (BAM-1020). The comparison yielded a low offset (0.3 µg m−3), although reference measurements tended to be lower than the sensor measurements (slope=0.57). A co-location dataset was then used to derive calibration factors incorporating the effects of temperature and relative humidity. The data processing procedures for obtaining the PM2.5 concentrations reported here are identical to those reported by Byrne et al. (2023).
The Cork dataset spans a similar measurement period to the Dungarvan dataset to allow for comparable results due to the known seasonality of PM2.5 pollution in Ireland (Ovadnevaite et al., 2021). Although the year 2021 included some periods of COVID-19 pandemic restrictions, such measures mainly affected NO2 concentrations and were not shown to have a significant impact on PM levels in Ireland (Environmental Protection Agency, EPA, 2020).
2.1.3 Meteorological measurements
Meteorological data was analysed in each location. For the city of Cork, data collected at Cork Airport by Ireland's National Meteorological Service, Met Éireann, was accessed from the website https://www.met.ie (last access: 25 April 2024). The airport weather station is located approximately 5.5 km from the Cork city centre.
There is no weather station located nearby Dungarvan that provides hourly measurements; however three of the Clarity Node-S devices were fitted with Clarity Wind Modules (AP7, AY93, AY9N), which provide high-time-resolution measurements of wind speed and direction (Clarity Movement Co., USA). Due to technical difficulties, device AY9N did not capture wind direction measurements, however its wind speed is included. The Wind Module contains a solid-state two-axis ultrasonic anemometer, which provides wind speed measurements with a range of 0–60.00 m s−1 and a resolution of 0.01 m s−1, along with wind direction at a resolution of 0.1° over a range of 0–359.9° (Wind Module technical sheet, 2024). These measurements have not been validated against reference meteorological data; however, they are included for indicative purposes.
2.2 Development of the concentration similarity index
The concentration similarity index (CSI) derived here quantifies the degree of likeness between PM2.5 concentration profiles from two sensors for a defined period of time and forms the basis for assessing the spatial disparities in PM2.5 measurements within sensor networks. The methodology proposed was developed through multiple iterations in order to adjust and improve the procedure. An overview of the development is described, showing the evolution towards the final method.
2.2.1 Original function application
The first phase of development was based directly on the work carried out by Piersanti et al. (2015), who used a concentration similarity function to assess the spatial representativeness of PM2.5 and O3 monitoring stations in the Italian air quality monitoring network. Using modelled hourly air pollutant data covering Italy with a 4 km×4 km grid cell resolution, Piersanti et al. (2015) produced maps showing how representative certain sites in the Italian monitoring infrastructure were. The application proposed here compares point measurement to point measurement as opposed to comparing modelled grid cell data; however, the underlying principle of comparing two concentration–time profiles to produce a single indication of similarity between them still applies. The function value fsite(x,y) used by Piersanti et al. (2015) to assess the spatial coverage of point measurements is given in Eq. (1):
where
and where represents the surface concentration from the modelled data in a grid point at time ti, represents the modelled data of a specific site of interest at time ti, and Nt is the total number of time steps. The study defined a modelled grid cell at the site of interest as representative of a surrounding grid cell area if the condition is true.
In the first step of our approach, this function was applied to the hourly average PM2.5 data obtained from the co-located Clarity Node-S units by comparing two sensor data series at a time. The concentration at the point of interest and surrounding grid cell concentration inputs were substituted for sensor PM concentration values from any given sensor A and sensor B pair, C(A,ti) and C(B,ti). Over a total of 1565 co-located hours, the mean number of comparable data points per C(A,ti), C(B,ti) pair was 654 due to devices being present at different stages during the co-location periods (Table S1).
Table 1Function values, fsite(x,y), for hourly averaged PM2.5 measured by a range of co-located Clarity Node-S devices. Device labels in the columns were set as , and device labels in the rows were set as .
![](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-t01-thumb.png)
In theory, the function value comparing two sensor data series would be 1, given that the measurements were collected in the same location and were known to all represent the same air parcel at each point in time. However, it was found that the function was not comprehensive enough to allow for an acceptable comparison of the sensor data. The results showed discrepancies between some device pairs because the function value deviated significantly from 1 in many cases (Table 1) and was as low as 0.51 in some cases, with an overall mean of 0.82.
2.2.2 Function parameter optimisation and introduction of PM limit
Analysis of the results obtained from direct application of the original function showed that the conditions set out by it were too strict to apply to the sensor data, given the variations that can occur in AQS measurements. The areas of the entire sensor networks discussed here could be within the original single grid cell size analysed by Piersanti et al. (2015). Therefore, overall pollution dynamics would vary significantly, in part because of hyper-local effects, and pollution averaging effects would be more pronounced when assessing larger areas. Moreover, the high hourly PM2.5 variation and very localised effects exhibited in a typical Irish winter PM2.5 profile are not suited to the original function (Byrne et al., 2023). While the original application contains a mathematical function examining the difference between two pollutant concentrations and is independent of specifications regarding area size and pollution dynamics, the threshold values can be adapted to reflect the specific application of the function. A second threshold value, a PM mass concentration limit, PMlim, was introduced to the function, with different relative concentration limits for the upper and lower PM values, Clim, upper and Clim, lower, respectively. Treating larger and smaller PM2.5 values differently when assessing the similarity between two data series is useful for capturing the nuanced relationships and patterns in the data. It allows for the real-world significance of the data to be reflected, acknowledging the varying implications of PM2.5 measurements based on the magnitude. Higher PM2.5 values can indicate a pollution episode or specific local pollution sources, while lower values can represent background levels. Therefore, treating lower PM2.5 values with more leniency in the similarity assessment recognises that minor fluctuations in low hourly concentrations might not be as concerning as similar deviations in higher concentrations and the health-related considerations associated with these high concentrations.
Another potential advantage of the PM limit concerns the varying degrees of accuracy of the AQS measurements. Allowing the leeway introduced here in assessing the similarity of lesser measurement values considers potential measurement uncertainties with these devices. However, it is important to note that this approach is not accommodating sensor limitations at the expense of accuracy, but rather it is a strategy to ensure that the assessment remains faithful to the underlying air quality dynamics while accounting for the potential deficiencies in measurement equipment.
When the function is applied to a pair of sensors, the resulting CSI can differ slightly depending on which sensor was classified as or , or sensor A or sensor B, in Eq. (1) when computing the difference at each time step. Due to the nature of the function, the denominator value of the relative difference calculation, the concentration of sensor A at a given time step, is what makes the difference. To counteract this and to avoid the possibility of large discrepancies between the CSI values for a sensor pair depending on which sensor is taken as A or B, the function was modified to use the geometric mean, or the square root of the product, of C(A,ti) and C(B,ti) as the denominator. This ensured symmetry in the function so that the CSI values were identical regardless of which sensor was classified as A or B in a sensor pair.
Equation (2) shows the next form of the concentration similarity function (function notation has been modified to be more suitable for this application).
where
and where C(A,ti) and C(B,ti) are the PM2.5 measurements from devices A and B at time, ti. Clim, lower and Clim, upper are the threshold values defining the acceptable level of difference between two concentrations, and PMlim is the PM mass concentration threshold value.
2.2.3 Development and testing of the modified equation
The PM limit and associated concentration similarity limits introduced were chosen by iteratively testing the similarity function on the co-location data using different limits. Each co-located sensor pair was tested with different PMlim values (5, 10, 15, 20 µg m−3) and with Clim values ranging from 0.1 to 2.0 in steps of 0.1 for both the upper and lower limits. This produced a Clim vs. CSI comparison for each A–B pair for data above and below the corresponding PMlim value. The Clim value for each sensor comparison, which gave a minimum CSI value of 0.95, was recorded with the overall mean of these Clim values above and below each PMlim value taken forward. The mean Clim pair values were then applied to the co-location measurements with the respective PMlim values to give final CSI values for each sensor pair, highlighting how the PM2.5 concentration profile of each sensor compares to that of all the other sensors. The highest mean CSI value for all co-located A–B pairs was found for PMlim=15 µg m−3, Clim, upper=0.2, and Clim, lower=0.7. When applying these new limits, all sensor pairs gave CSI>0.85, with 99 % of pairs above 0.90 and with an overall CSI mean of 0.98. These final limits enabled a good comparison for the hourly co-located AQS measurements (Table 2).
Table 2Concentration similarity indices for hourly averaged PM2.5 measured by a range of co-located Clarity Node-S devices. PMlim=15 µg m−3, Clim, upper=0.2, and Clim, lower=0.7.
![](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-t02-thumb.png)
The CSI function was also applied to data obtained from the four co-located PurpleAir devices in order to make sure that the function was applicable across the two AQS types. The data was harmonised by following the same procedure as the Clarity device data, through scaling each data point from each sensor to the mean data series of all four sensors. Although this co-location period was shorter than that of the Clarity dataset used for the function development, it still allowed for the CSI to be calculated from around 250 common data points per sensor pair. All device pairs reported a CSI close to 1.0, with a mean CSI of 0.99 (Table S3 in the Supplement).
The function described in Eq. (2) was further tested by comparing the sensors to numerous sets of synthetic data created from each sensor's measurements to assess the impact of a range of scenarios. Comparing a sensor dataset to itself establishes a baseline for the comparison where the CSI is 1 and any subsequent adjustments to the data to create the synthetic data can be explored, resulting in a new CSI. The first scenario investigated changes in CSI when outliers are present in the data. To explore this, the sensor data was changed so a certain number of data points could be considered outliers (n=1, 10, 500, 1000). An outlier data point was created by increasing a value by 100 µg m−3 in order to ensure discrepancy between it and the original value. The function was then tested in a scenario where the data was scaled linearly so the mean remained constant, but the variance of the data was increased, and it was also tested in a scenario where the entire dataset was offset by 5, 10, 15, and 20 µg m−3. The final test scenario involved the introduction of noise to the dataset, representing impactful variations in that data. Gaussian noise with various values of standard deviation was added to the data. The CSI results for the synthetic data tests were also compared to the results when the R2 was found between any two given datasets. Low variations were found during all synthetic data analyses, with the resulting CSI values having standard deviations ≤0.05 across the individual devices for each test. As an example, the effects of these tests on the CSI results for A6P are shown in Table 3, where 4406 data points were included in the calculations.
Table 3Influence of data outliers and other factors on CSI determined in test scenarios with device A6P.
![](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-t03-thumb.png)
It is clear that in the case of the linearly scaled data with higher variability but the same overall mean, the CSI is impacted (CSI=0.52) when the standard deviation is increased by just a factor of 1.5, indicating that such a dataset is dissimilar to the original. In comparison, the R2 is not an accurate reflection of the change, as it does not deviate from 1. Offsetting the data by different degrees also shows a major change in the CSI (CSI<0.55). However, this is not reflected well in the R2 values, which do not deviate from 1. The CSI method is quite robust with respect to outliers, whereas the R2 is more sensitive (R2=0.94) when 10 outliers are introduced to the dataset, which is approximately 0.2 % of the total data points. The R2 is significantly reduced (R2=0.45) with 500 outliers (∼11 % of the total data points), whereas the CSI is only slightly impacted (CSI=0.89). As the method yields a time-averaged result, low numbers of outliers do not hugely affect the index for a given sensor pair. So, two datasets that are generally similar but where one experiences some outliers will be deemed similar by the method. The R2 also shows a more limited response when larger amounts of Gaussian noise are added, resulting in a value of 0.96 when the standard deviation of the noise is 4 µg m−3, while the CSI is adjusted to 0.7. From a health-impact and exposure point of view, increased variation and higher offset represent very different exposure scenarios, whereas a large difference in the occasional hourly average in an otherwise similar exposure regime does not. The CSI offers a more appropriate comparison between hourly measurements collected at two locations.
2.3 Application to sensor networks and analysis of spatial trends
The CSI methodology developed above was subsequently applied to the Dungarvan and Cork sensor networks to evaluate the similarity and spatial variations in PM2.5. A systematic pairwise comparison approach was employed, wherein each sensor was individually compared to every other sensor within the network. Hierarchical clustering and fuzzy c-means (FCM) clustering were both performed on the CSI results to identify groupings based on each sensor's relationship to other sensors in the network that can then be reflected spatially.
With both clustering techniques, the quality of cluster assignments can be assessed with various evaluation metrics to choose the optimal number of clusters. As the “true” cluster classifications are not known here, validation must be performed using the clustering algorithm itself. To assess the quality of the hierarchical clustering assignments, the silhouette metric was used along with the Calinski–Harabasz index to assess the FCM assignments (Caliński and Harabasz, 1974; Rousseeuw, 1987). The silhouette score, ranging from −1 to +1, can be calculated for each member of a cluster and then the mean silhouette score from all members indicates an overall assignment quality for members of that cluster, with a high score closer to 1 indicating higher-quality clusters and a low or negative score indicating poorer cluster assignments. The Calinski–Harabasz index also quantifies the quality of cluster assignments, with higher scores indicating better quality. The metrics were used to test for the optimal number of clusters for each algorithm.
3.1 Dungarvan PM2.5 sensor network
Analysis of the harmonised data obtained from the sensors in the Dungarvan PM2.5 network was conducted to determine CSI values and assess the spatial variation in air pollution across the town. Although the PM2.5 concentrations are not as accurate as those collected by reference instrumentation, any relative differences between the sensors and between individual sensor data trends can be regarded as genuine due to the low inter-sensor variation observed after data harmonisation procedures, where the standard deviation of the mean PM2.5 co-located measurements was 1.7 µg m−3.
The temporal and spatial trends of PM2.5 across the Dungarvan sensor network are reflected in the average diurnal plots obtained for each sensor, Fig. 1. These diurnal profiles all show large evening peaks in PM2.5, which are typical of towns and cities in Ireland affected by residential solid-fuel burning during winter evenings (Dall'Osto et al., 2014; Healy et al., 2010; Wenger et al., 2020). However, there are clear disparities in some of the average evening peak values between the sensors. One group of sensors has maximum values above 35 µg m−3 (A3, A4, A8H, A9, AQ, A7, AW6, AQV), while the sensors with maxima below 35 µg m−3 can be further divided into three smaller groups. Sensors labelled AJ3, AWF, and AZ all have a maximum PM2.5 concentration around 30 µg m−3; sensors AY9N, AY93, ARF, A8Z, AYG, and A6P all have maxima in the 20–26 µg m−3 range, while AP7 has a significantly lower evening peak than all other devices.
![https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f01](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f01-thumb.png)
Figure 1Diurnal profiles for hourly averaged measurements of PM2.5 in the Dungarvan sensor network (September 2022 to May 2023).
Most sensors exhibited the diurnal maximum around the same time of day, between 18:00 and 20:00 LT; however AP7 and ARF showed a slightly delayed peak from 20:00 to 22:00 LT. AP7 had the lowest peak concentration and did not exhibit the sharp rise and subsequent decrease associated with evening solid-fuel burning that the other sensors showed. AP7 was located on the southwestern edge of the town, and since the predominant wind direction is southwesterly, it did not measure as much local pollution as other locations in the eastern part of the network.
Summary statistics obtained for the 18 sensors in the Dungarvan network are listed in Table 4. Unsurprisingly, most of the devices with diurnal maxima >35 µg m−3 have the highest mean, median, and maximum values. Out of this subset of devices, AQV has the lowest overall mean (15 µg m−3) but still has a relatively high standard deviation (22 µg m−3), indicating that the PM2.5 values tend to vary widely but are lower on average. This could be indicative of fluctuating particle concentrations, consistent with intermittent pollution sources such as residential solid-fuel burning.
Table 4Summary statistics of hourly average PM2.5 concentrations obtained for all sensors in the Dungarvan sensor network (September 2022 to May 2023).
![](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-t04-thumb.png)
The wind speed and direction recorded at sites AP7 and AY93 showed some variation (Fig. S4a and b in the Supplement); however, wind speeds measured at all three sites showed a moderate correlation with all R2 values above 0.65. The measured wind direction at the AP7 and AY93 sites reported a moderate correlation (R2=0.63). Both sites measured winds emanating from a broad range of directions. Both locations reported generally southerly winds 53 % of the time and southwesterly winds 30 % of the time. The temporal variations in wind speed measured at the three sites are detailed in Fig. S5 in the Supplement. Little diurnal variation is seen between devices AY93 and AY9N; however, it is clear that AP7 tended to experience slightly lower wind speeds than AY93 and AY9N during the measurement campaign. Nevertheless, this difference did not exceed 1 m s−1 in any of the temporal variation assessments, and all three sites reported the same overall trends in wind speed. The variations in wind measurements between the sites indicate some slight local meteorological differences; however, the overall meteorological field is not likely to differ greatly between the three sites.
3.1.1 Concentration similarity index
The matrix of CSI values obtained for the Dungarvan sensor network is shown in Table 5. The results can be analysed in a number of ways. Firstly, the indices for one sensor can be used to assess how similar or dissimilar the measurements are to all other sensors in the network, thus providing information on the spatial representativeness of that particular location. Secondly, the indices of all sensors can be looked at together to elucidate any potential relationship between sensor measurement locations.
Table 5Concentration similarity indices for the hourly averaged PM2.5 concentrations measured by Clarity Node-S devices in the Dungarvan sensor network.
![](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-t05-thumb.png)
The minimum CSI value (0.85) determined during the co-location deployment can act as the lower limit for when two sensor locations can be considered very similar. The reported CSI values for Dungarvan sensors ranged from 0.48 (ARF vs. A7) to 0.79 (AYG vs. AWF) with a mean of 0.61, indicating a significant difference in air quality representation between locations across the town. The device with the lowest mean of its CSI values with respect to the other locations was A4 (0.55), and although device ARF was only slightly above this (0.57), it reported a larger range of CSI values, including the lowest of the entire dataset. AJ3, AQV, and AYG all shared the highest mean CSI values (0.66).
To further investigate the effect of solid-fuel burning on local air quality, the CSI function was applied to data from two isolated months – January and May 2023. The purpose of this assessment was to evaluate the extent to which residential solid-fuel burning dictates the CSI between two sensors, given that one month (January) will have higher PM2.5 levels, with measurements heavily influenced by solid-fuel burning, and the other will not (May). For both months, all sensors had data capture above 65 %, and the mean capture was 94 % for January and 92 % for May. The January mean CSI from all comparisons was 0.51, and the May mean CSI was 0.84 (Tables S4 and S5 in the Supplement). The large discrepancy between the mean CSI for January and May is most likely due to the higher variation typically seen in wintertime PM2.5 (sJanuary=25 µg m−3, sMay=9 µg m−3) due to residential solid-fuel burning (Fig. 1). This highlights the importance of seasonality when assessing the spatial representativeness of monitoring network locations.
3.1.2 Clustering
Clustering techniques were employed on the CSI matrix to uncover any inherent spatial relationships between different locations in the network. Hierarchical clustering produced a dendrogram showing the hierarchical relationship between the sensor locations and was used to identify clusters (Fig. 2). The highest mean silhouette score was found with two clusters (Fig. S6 in the Supplement). However, it was not a high silhouette score (0.19), indicating that the quality of the cluster assignments was low. The highest Calinski–Harabasz index corresponded to the assignment of members to two clusters when applying the FCM clustering (Fig. S7 in the Supplement).
![https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f02](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f02-thumb.png)
Figure 2Dendrogram output from hierarchical clustering of the CSI data from the Dungarvan sensor network.
From both the dendrogram (Fig. 2) and the FCM membership weights (Fig. 3), it is clear that devices A4 through AQ are grouped together in one cluster (Cluster 1), and devices AQV to AP7 are grouped in another cluster (Cluster 2). This split is very similar to the easily visualised groupings shown in the diurnal profile maxima (Fig. 1), with the only difference being device AQV. The devices in Cluster 1 are also those with the highest mean PM2.5 for the measurement period. The mean CSI for each sensor mostly corresponds to the cluster assignments, with Cluster 1 devices having a mean CSI equal to or below 0.6 and all devices in Cluster 2 having a mean CSI above 0.6 except for device ARF. Interestingly, this grouping also appears to have spatial importance too, as shown in Fig. 4. Cluster 2 devices are mainly located around the edge of the town and generally experience cleaner air ( µg m−3, sPM2.5=27 µg m−3), while Cluster 1 devices are located in central and residential areas ( µg m−3, sPM2.5=17 µg m−3), which are more polluted during winter months.
![https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f04](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f04-thumb.jpg)
Figure 4Dungarvan AQS locations with two cluster groups indicated. Cluster 1 devices (red triangle markers) are mainly located in central and residential areas, while Cluster 2 devices (blue cross markers) are mainly located on the edge of the town. (map obtained from Esri, DigitalGlobe, GeoEye, i-cubed, USDA FSA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community).
3.2 Cork PM2.5 sensor network
The same approach as above was used to analyse the data collected by the Cork AQS network. In this case, the corrected measurements are indicative of the actual PM2.5 experienced in each location. The diurnal plots for each sensor in the Cork network are similar to those observed in Dungarvan, with a sizeable evening peak in PM2.5 concentrations (19:00–21:00 LT) due to emissions from residential solid-fuel burning. Again, there is considerable variation in the peak concentration of PM2.5 (Fig. 5). Device MTU showed the lowest diurnal average maximum of 9 µg m−3. This device is located on the western side of the city and has few upwind pollution sources contributing to air pollution at the location, as the prevailing wind direction is from the southwest. Devices CCC12 and CCC9 both showed the highest diurnal average maximum, 17 µg m−3. CCC12 is located northeast of the city and so likely experiences urban PM2.5 sources upwind from it or has strong localised sources. Similarly, CCC9 is located to the east of the city in a residential area. Table 6 contains summary statistics for each of the sensors in the Cork network. Some devices had very high PM2.5 maxima, e.g. 201 µg m−3 for CCC11, which were more than double the maxima of other devices, e.g. CCC8, which had the lowest overall maximum of 47 µg m−3. Device MTU had the lowest diurnal maximum value, indicating that this location was the least affected by local emissions from solid-fuel burning. However, it measured a significant overall PM2.5 maximum of 99 µg m−3 and significant spikes in pollution were occasionally observed, likely due to meteorological conditions or specific localised effects. When looking at all of the parameters listed in Table 6, CCC11 stands out. This sensor has the highest maximum hourly average PM2.5 concentration in the network, but the standard deviation (8 µg m−3) is in the middle of the range, indicating that the location had relatively stable PM2.5 levels throughout the measurement period with less variation than other devices, but it was still susceptible to occasional spikes in PM2.5.
![https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f05](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f05-thumb.png)
Figure 5Diurnal PM2.5 profiles for all AQSs in the Cork network (January to May and September to December 2021).
Table 6Summary statistics of hourly averaged PM2.5 obtained for all sensors in the Cork network (January to May and September to December 2021).
![](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-t06-thumb.png)
The meteorological data retrieved from Cork Airport, which is approximately 4–11 km from each device in the Cork sensor network, was investigated for the measurement period in 2021. While data obtained from the airport site indicates the meteorological conditions on a synoptic scale, the local weather experienced at individual locations within the city are additionally shaped by factors such as street canyon effects and local topography. Consequently, the wind direction measured at the airport site cannot be assumed to mirror that of all devices in the network. Wind speeds measured at the airport generally surpass those within the city as it is situated at a higher elevation than the city. However, the broader regional wind patterns are expected to exert a predominant influence on the overall meteorological conditions across the city, and therefore the relationship with meteorological conditions and local PM2.5 levels can be investigated. The Cork Airport site recorded southerly winds 59 % of the time and southwesterly winds 39 % of the time (Fig. S8 in the Supplement).
3.2.1 Concentration similarity index
The matrix of CSI values obtained for the Cork sensor network is shown in Table 7. The values range from 0.52 (CCC12 vs. MTU and CCC9 vs. MTU) to 0.85 (CCC2 vs. CCC11) with a mean of 0.71. The high maximum CSI indicates a high degree of similarity between those locations in the network, and overall, Cork locations show a higher degree of similarity compared to those in Dungarvan.
Table 7Concentration similarity indices for the hourly averaged PM2.5 concentrations measured by PurpleAir devices in the Cork AQS network.
![](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-t07-thumb.png)
The isolated CSI results for the months of January and May 2021 were also assessed for the city of Cork. The average data coverage during both periods was 92 %. The mean CSI value in January (0.55) was considerably lower than that observed in May (0.82), Tables S6 and S7 in the Supplement. This result is similar to that found for the Dungarvan network, again indicating that the large difference in mean scores between the two months can be attributed to higher wintertime PM2.5 variation due to residential solid-fuel burning (sJanuary=15 µg m−3, sMay=3 µg m−3).
3.2.2 Clustering
The two clustering algorithms were applied to investigate the CSI results of the Cork network. The silhouette scores for each number of assigned clusters (2 to 5) were low, with two clusters showing the highest mean score (Fig. S9 in the Supplement). Similarly, with the FCM analysis, two clusters showed the highest score with the Calinski–Harabasz indices (Fig. S10 in the Supplement).
The dendrogram produced from the hierarchical clustering and the membership weights for two clusters from FCM clustering are shown in Figs. 6 and 7, respectively. It is clear that devices MTU, CCC5, and CCC8 are all grouped together in one branch, Cluster 2, with the remainder of the devices in the other branch. The one assignment difference between the two clustering methods is CCC3, which has a higher membership weight towards Cluster 2 with the FCM method but does not branch with that cluster in the dendrogram. However, its membership weight is close to 0.5. CCC1 also shows a membership weight close to 0.5, however it is showing a higher weight towards Cluster 1, as per the hierarchical clustering results. Devices in Cluster 2, except for CCC3, all have the lowest mean CSI values.
![https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f06](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f06-thumb.png)
Figure 6Dendrogram output from hierarchical clustering of the CSI data from the Cork sensor network.
Similar to the Dungarvan results, there appears to be a spatial component to the cluster groupings, with devices in Cluster 2 being mainly on the western side of the city, Fig. 8. However, the contrast in cluster PM2.5 mean values is not as stark in the Cork clusters as in those in Dungarvan. Cluster 1 had a mean PM2.5 of 9 µg m−3, while Cluster 2 had a mean PM2.5 of 7 µg m−3. Interestingly, device CCC7, located in a commuter town on the western side of the city boundary, is grouped in Cluster 1 along with devices mainly in urban residential type sites instead of being grouped with other devices on the western edge of the city. This indicates that it has a more comparable CSI profile to the urban residential sites than the locations closer to it, further emphasising the importance of location type over physical proximity.
![https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f08](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f08-thumb.jpg)
Figure 8Cork AQS locations with two cluster groups indicated. Cluster 1 devices (red triangle markers) are located in the city centre and east/northeast, while Cluster 2 devices (yellow circle markers) are mainly located on the western side of the city. (map obtained from Esri, DigitalGlobe, GeoEye, i-cubed, USDA FSA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community).
3.3 Application of the CSI to assess representativeness of air quality monitoring locations
One key benefit of the CSI metric for AQS networks is that one sensor can be singled out and its overall degree of similarity to measurements from other locations can be determined. This analysis can be used to assess the spatial representativeness of a given location in the AQS network by quantitatively exploring how similar its PM2.5 profile is to other locations. If a network sensor is co-located with a reference instrument, then the CSI values for that sensor can be used to provide a measure of the representativeness of the designated monitoring location and how well it informs the assessment of population exposure to air pollution.
In Dungarvan, the device A6P was co-located with a PM2.5 instrument (Osiris, Turnkey) deployed as part of the national air quality monitoring network. The instrument is not a reference instrument but is certified to provide indicative measurements of PM2.5 (National Ambient Air Quality Monitoring Network, 2023; Osiris, 2024). A6P had a mean CSI of 0.63, the fifth-highest of the mean CSI values across all devices. The similarity indices for A6P are included in Table 5 and represented spatially in Fig. 9. All CSI values are below the minimum threshold of 0.85 for two Clarity S-node devices in the Dungarvan network to be considered very similar. The most similar devices are found to the northeast of this location, AJ3 and AYG. Interestingly, the similarity of PM profiles does not decrease with increasing distance from A6P. Devices on the furthest western (AZ, A8Z, AP7) and eastern (AWF, AY9N) edges of the town are within 0.6 to 0.7, yet devices A4, A7, AQ, and A9 are all at or below 0.6 despite being physically closer to A6P. This suggests that the location type as opposed to physical proximity is more important when it comes to assessing the similarity of locations within Irish towns, as A4, A7, AQ, and A9 are all fully surrounded by residential areas, whereas the other devices mentioned are in more open areas.
![https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f09](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f09-thumb.jpg)
Figure 9Dungarvan AQS locations with CSI results indicated in coloured circles (blue is the lowest CSI; yellow is the highest CSI) and the A6P location indicated by red pin marker. (map obtained from Esri, DigitalGlobe, GeoEye, i-cubed, USDA FSA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community).
One of the devices in the Cork sensor network, UCC, was co-located alongside a reference instrument (BAM-1020) at the national air quality monitoring location on the UCC campus. The CSI values for the device labelled UCC are shown in Fig. 10, showing how similar the measurements at this site are compared to the rest of the locations in the sensor network. The CSI scale on the map has been adjusted for these values. Similar to the Dungarvan case, there are devices that show high similarity (CCC4, CCC2, CCC12, CCC1) with UCC that are not located nearby.
![https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f10](https://amt.copernicus.org/articles/17/5129/2024/amt-17-5129-2024-f10-thumb.jpg)
Figure 10Cork AQS locations with UCC CSI results indicated in coloured circles (blue is the lowest CSI; yellow is the highest CSI) and the UCC location indicated by red pin marker. (map obtained from Esri, DigitalGlobe, GeoEye, i-cubed, USDA FSA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community).
A robust framework for comparing data series from individual air quality sensors in a network has been established, and a new metric, the concentration similarity index (CSI), has been developed, optimised, and tested on a co-location dataset. The CSI allows one to consider the monitoring network in terms of the similarity of the concentration–time profile of PM2.5 at one location to profiles at the other locations in the same network. The harmonised dataset with minimal unexplained inter-sensor variation underpins the development of the CSI method, along with robust tests to ensure that the function represents an unbiased and fair depiction of the inter-sensor relationships after deployment in a monitoring network.
The CSI method has been used to analyse data generated by PM2.5 sensor networks in two locations in Ireland: the coastal town of Dungarvan and the city of Cork. Clustering techniques are applied to the CSI matrix, and comparable similarity trends between locations drive the distinctions made with the clustering algorithm. The resulting groupings can provide several insights into the PM2.5 profile at each location, including the likelihood of similarity in pollution sources, spatial patterns, and temporal trends. An interesting contrast in the CSI results from the two monitoring networks was obtained from the clustering analysis. In Dungarvan, the locations generated clusters that were reflected well when comparing the individual diurnal profiles and specifically the diurnal maximum values, indicating that this factor has a major influence when relating the concentration–time profiles at each location to one another in this network. However, for the Cork network results, this was not as apparent. The clusters were not aligned based on diurnal peaks, but rather the differentiating factor was more nuanced. Both clusters contained locations with a mix of higher and lower diurnal maxima and overall maxima. However, both network groupings reflect the fact that devices may report dissimilar CSI results to other devices located nearby and that considering location specifications or types, such as residential areas, is more important than physical proximity when it comes to understanding and quantifying the similarities between locations.
The CSI function was also applied to two separate months in the network datasets, with January chosen to represent a period of higher PM2.5 levels due to solid-fuel burning emissions and May chosen to represent a period with lower PM2.5 concentrations due to reduced solid-fuel burning. In both locations, the mean CSI for the network comparisons was higher in May than in January, indicating that higher PM2.5 levels are a major driver of lower similarity indices between sensor locations. Combining this with the findings of our previous study, we provide further evidence that high levels of localised PM2.5 cause distinct disparities in exposure to poor air quality in different locations. Furthermore, to properly assess the burden of PM2.5 experienced by a population and to accurately compare the measurements at two locations, the wintertime PM data must be included in the assessment.
The similarity of PM2.5 measured at designated sites in the national air quality monitoring network compared to the rest of the locations in the sensor networks was analysed to give an estimation of the representativeness of the air pollution measured at the designated monitoring site. The national monitoring site location in Dungarvan was shown to be moderately representative of the other AQS network locations in the town, with CSI values ranging from 0.53 to 0.72. The CSI values for the Cork comparison ranged from 0.62 to 0.84, also showing a fair representation of the air pollution experienced in the rest of the network. The CSI function was also tested via synthetic datasets, which showed that a positive offset of just 5 µg m−3 resulted in almost halving the CSI, which was a lower CSI than most of the sensor comparisons in both network locations. So, while a CSI of 0.85 was used as a limit for two sensor measurement sets being very similar, CSI values between 0.6 and 0.7 are still moderately similar. In general, the CSI values in the city of Cork for the reference site comparison were higher (mean=0.75) than that of Dungarvan (mean=0.63), indicating less similarity between the reference site and devices in the Dungarvan network compared to those in the city of Cork.
While the function was developed and tested on multiple sensor pairs and further validated with additional co-located pairs, validation with co-located PM2.5 measurements of the PMlim, Clim, upper, and Clim, lower parameters for specific applications is recommended to ensure that the index represents the dataset accurately. Co-location assessments are also recommended to ensure minimal inter-sensor variation. Nonetheless, the differentiation between higher and lower PM values in the concentration similarity assessment is a strategic choice that acknowledges the complexity of PM2.5 data, the varying significance of concentration levels, and the limitations of sensors. It allows for a more accurate representation of similarities, while considering real world implications and measurement uncertainties, and minimises the potential biases that could arise from an indiscriminate approach, thus ensuring an impartial and unbiased evaluation.
The analysis and application of the CSI function displays the potential for AQS networks to be used in conjunction with a regulatory monitoring system. This study has shown the potential for sensor networks to assess the need for more regulatory monitoring in an area and to identify locations that are being poorly represented by the current system. Furthermore, the CSI method can be used to optimise a sensor network by carrying out short-term sensor deployments and identifying areas of similarity or dissimilarity and thus assessing where the best locations for sensors are, based on the similarity in exposure to air pollution.
All raw data are available upon request.
The supplement related to this article is available online at: https://doi.org/10.5194/amt-17-5129-2024-supplement.
RB and SH conceptualised the project and the methodology. RB carried out the formal data analysis, investigation, and visualisation and wrote the paper with supervision and contributions from SH and JCW.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
The authors acknowledge Cork City Council, especially Kevin Ryan, for developing and maintaining the air quality sensor network in the city of Cork and Waterford City and County Council for supporting the Dungarvan measurement campaign. In particular, we would like to thank Paul Flynn, who facilitated the physical deployment of the air quality units in Dungarvan.
This research has been supported by the Environmental Protection Agency Ireland and the EU LIFE programme through LIFE EMERALD (grant no. LIFE19 GIE/IE/001101).
This paper was edited by Albert Presto and reviewed by Matthew Johnson and one anonymous referee.
Bezanson, J., Edelman, A., Karpinski, S., and Shah, V. B.: Julia: A fresh approach to numerical computing, SIAM Rev., 59, 65–98, https://doi.org/10.1137/141000671, 2017.
Byrne, R., Ryan, K., Venables, D. S., Wenger, J. C., and Hellebust, S.: Highly local sources and large spatial variations in PM2.5 across a city: evidence from a city-wide sensor network in Cork, Ireland, Environmental Science: Atmospheres, 3, 919–930, https://doi.org/10.1039/D2EA00177B, 2023.
Caliñski, T. and Harabasz, J.: A Dendrite Method For Cluster Analysis, Commun. Stat., 3, 1–27, https://doi.org/10.1080/03610927408827101, 1974.
Carslaw, D. C. and Ropkins, K.: openair – An R package for air quality data analysis, Environ. Model. Softw., 27–28, 52–61, https://doi.org/10.1016/j.envsoft.2011.09.008, 2012.
Cesaroni, G., Badaloni, C., Gariazzo, C., Stafoggia, M., Sozzi, R., Davoli, M., and Forastiere, F.: Long-term exposure to urban air pollution and mortality in a cohort of more than a million adults in Rome, Environ. Health Persp., 121, 324–331, https://doi.org/10.1289/EHP.1205862, 2013.
Clarity Movement Co.: https://www.clarity.io/, last access: 29 August 2023.
Crawford, B., Hagan, D. H., Grossman, I., Cole, E., Holland, L., Heald, C. L., and Kroll, J. H.: Mapping pollution exposure and chemistry during an extreme air quality event (the 2018 Kīlauea eruption) using a low-cost sensor network, P. Natl. Acad. Sci. USA, 118, e2025540118, https://doi.org/10.1073/pnas.2025540118, 2021.
Dall'Osto, M., Ovadnevaite, J., Ceburnis, D., Martin, D., Healy, R. M., O'Connor, I. P., Kourtchev, I., Sodeau, J. R., Wenger, J. C., and O'Dowd, C.: Characterization of urban aerosol in Cork city (Ireland) using aerosol mass spectrometry, Atmos. Chem. Phys., 13, 4997–5015, https://doi.org/10.5194/acp-13-4997-2013, 2013.
Dall'Osto, M., Hellebust, S., Healy, R. M., Connor, I. P., Kourtchev, I., Sodeau, J. R., Ovadnevaite, J., Ceburnis, D., O'Dowd, C. D., and Wenger, J. C.: Apportionment of urban aerosol sources in Cork (Ireland) by synergistic measurement techniques, Sci. Total Environ., 493, 197–208, https://doi.org/10.1016/J.SCITOTENV.2014.05.027, 2014.
Diez, S., Lacy, S. E., Bannan, T. J., Flynn, M., Gardiner, T., Harrison, D., Marsden, N., Martin, N. A., Read, K., and Edwards, P. M.: Air pollution measurement errors: is your data fit for purpose?, Atmos. Meas. Tech., 15, 4091–4105, https://doi.org/10.5194/amt-15-4091-2022, 2022.
Environmental Protection Agency (EPA): Air Quality in Ireland 2020, https://www.epa.ie/publications/monitoring--assessment/air/air-quality-in-ireland-2020.php (last access: 17 August 2023), 2020.
Fajersztajn, L., Saldiva, P., Pereira, L. A. A., Leite, V. F., and Buehler, A. M.: Short-term effects of fine particulate matter pollution on daily health events in Latin America: a systematic review and meta-analysis, Int. J. Public Health, 62, 729–738, https://doi.org/10.1007/S00038-017-0960-Y, 2017.
Frederickson, L. B., Sidaraviciute, R., Schmidt, J. A., Hertel, O., and Johnson, M. S.: Are dense networks of low-cost nodes really useful for monitoring air pollution? A case study in Staffordshire, Atmos. Chem. Phys., 22, 13949–13965, https://doi.org/10.5194/acp-22-13949-2022, 2022.
Frederickson, L. B., Russell, H. S., Fessa, D., Khan, J., Schmidt, J. A., Johnson, M. S., and Hertel, O.: Hyperlocal air pollution in an urban environment – measured with low-cost sensors, Urban Clim., 52, 101684, https://doi.org/10.1016/J.UCLIM.2023.101684, 2023.
Healy, R. M., Hellebust, S., Kourtchev, I., Allanic, A., O'Connor, I. P., Bell, J. M., Healy, D. A., Sodeau, J. R., and Wenger, J. C.: Source apportionment of PM2.5 in Cork Harbour, Ireland using a combination of single particle mass spectrometry and quantitative semi-continuous measurements, Atmos. Chem. Phys., 10, 9593–9613, https://doi.org/10.5194/acp-10-9593-2010, 2010.
Heimann, I., Bright, V. B., McLeod, M. W., Mead, M. I., Popoola, O. A. M., Stewart, G. B., and Jones, R. L.: Source attribution of air pollution by spatial scale separation using high spatial density networks of low cost air quality sensors, Atmos. Environ., 113, 10–19, https://doi.org/10.1016/J.ATMOSENV.2015.04.057, 2015.
Hodoli, C. G., Coulon, F., and Mead, M. I.: Source identification with high-temporal resolution data from low-cost sensors using bivariate polar plots in urban areas of Ghana, Environ. Pollut., 317, 120448, https://doi.org/10.1016/J.ENVPOL.2022.120448, 2023.
Kassomenos, P. A., Vardoulakis, S., Chaloulakou, A., Paschalidou, A. K., Grivas, G., Borge, R., and Lumbreras, J.: Study of PM10 and PM2.5 levels in three European cities: Analysis of intra and inter urban variations, Atmos. Environ., 87, 153–163, https://doi.org/10.1016/J.ATMOSENV.2014.01.004, 2014.
Kaur, K. and Kelly, K. E.: Performance evaluation of the Alphasense OPC-N3 and Plantower PMS5003 sensor in measuring dust events in the Salt Lake Valley, Utah, Atmos. Meas. Tech., 16, 2455–2470, https://doi.org/10.5194/amt-16-2455-2023, 2023.
Kourtchev, I., Hellebust, S., Bell, J. M., O'Connor, I. P., Healy, R. M., Allanic, A., Healy, D., Wenger, J. C., and Sodeau, J. R.: The use of polar organic compounds to estimate the contribution of domestic solid fuel combustion and biogenic sources to ambient levels of organic carbon and PM2.5 in Cork Harbour, Ireland, Sci. Total Environ., 409, 2143–2155, https://doi.org/10.1016/J.SCITOTENV.2011.02.027, 2011.
Kumar, P., Morawska, L., Martani, C., Biskos, G., Neophytou, M., Di Sabatino, S., Bell, M., Norford, L., and Britter, R.: The rise of low-cost sensing for managing air pollution in cities, Environ. Int., 75, 199–205, https://doi.org/10.1016/J.ENVINT.2014.11.019, 2015.
Lelieveld, J., Evans, J. S., Fnais, M., Giannadaki, D., and Pozzer, A.: The contribution of outdoor air pollution sources to premature mortality on a global scale, Nature, 525, 367–371, https://doi.org/10.1038/nature15371, 2015.
Li, H. Z., Gu, P., Ye, Q., Zimmerman, N., Robinson, E. S., Subramanian, R., Apte, J. S., Robinson, A. L., and Presto, A. A.: Spatially dense air pollutant sampling: Implications of spatial variability on the representativeness of stationary air pollutant monitors, Atmos. Environ. X, 2, 100012, https://doi.org/10.1016/J.AEAOA.2019.100012, 2019.
Lin, C., Huang, R. J., Ceburnis, D., Buckley, P., Preissler, J., Wenger, J., Rinaldi, M., Facchini, M. C., O'Dowd, C., and Ovadnevaite, J.: Extreme air pollution from residential solid fuel burning, Nat. Sustain., 1, 512–517, https://doi.org/10.1038/s41893-018-0125-x, 2018.
Lin, C., Ceburnis, D., Huang, R.-J., Xu, W., Spohn, T., Martin, D., Buckley, P., Wenger, J., Hellebust, S., Rinaldi, M., Facchini, M. C., O'Dowd, C., and Ovadnevaite, J.: Wintertime aerosol dominated by solid-fuel-burning emissions across Ireland: insight into the spatial and chemical variation in submicron aerosol, Atmos. Chem. Phys., 19, 14091–14106, https://doi.org/10.5194/acp-19-14091-2019, 2019.
Malings, C., Tanzer, R., Hauryliuk, A., Saha, P. K., Robinson, A. L., Presto, A. A., and Subramanian, R.: Fine particle mass monitoring with low-cost sensors: Corrections and long-term performance evaluation, Aerosol Sci. Tech., 54, 160–174, https://doi.org/10.1080/02786826.2019.1623863, 2020.
Munir, S., Mayfield, M., Coca, D., Jubb, S. A., and Osammor, O.: Analysing the performance of low-cost air quality sensors, their drivers, relative benefits and calibration in cities-a case study in Sheffield, Environ. Monit. Assess., 191, 94, https://doi.org/10.1007/S10661-019-7231-8, 2019.
National Ambient Air Quality Monitoring Network: https://airquality.ie/, last access: 26 October 2023.
Node-S technical sheet: https://click.clarity.io/hubfs/Marketing%20Assets%20-%20PDFs/Product%20and%20Specification%20Sheets/Node-S%20Specifications%20Sheet.pdf, last access: 29 August 2023.
O'Regan, A. C., Byrne, R., Hellebust, S., and Nyhan, M. M.: Associations between Google Street View-derived urban greenspace metrics and air pollution measured using a distributed sensor network, Sustain. Cities Soc., 87, 104221, https://doi.org/10.1016/J.SCS.2022.104221, 2022.
Orellano, P., Reynoso, J., Quaranta, N., Bardach, A., and Ciapponi, A.: Short-term exposure to particulate matter (PM10 and PM2.5), nitrogen dioxide (NO2), and ozone (O3) and all-cause and cause-specific mortality: Systematic review and meta-analysis, Environ. Int., 142, 105876, https://doi.org/10.1016/J.ENVINT.2020.105876, 2020.
Osiris: https://turnkey-instruments.com/product/osiris/, last access: 26 February 2024.
Ovadnevaite, J., Lin, C., Rinaldi, M., Ceburnis, D., Buckley, P., Coleman, L., Facchini, M. C., Wenger, J., and O'Dowd, C.: Air Pollution Sources in Ireland, Environmental Protection Agency, Ireland, ISBN 978-1-80009-007-1, 2021.
Pedersen, M., Giorgis-Allemand, L., Bernard, C., Aguilera, I., Andersen, A. M. N., Ballester, F., Beelen, R. M. J., Chatzi, L., Cirach, M., Danileviciute, A., Dedele, A., van Eijsden, M., Estarlich, M., Fernández-Somoano, A., Fernández, M. F., Forastiere, F., Gehring, U., Grazuleviciene, R., Gruzieva, O., Heude, B., Hoek, G., Hoogh, K. de, van den Hooven, E. H., Håberg, S. E., Jaddoe, V. W. V., Klümper, C., Korek, M., Krämer, U., Lerchundi, A., Lepeule, J., Nafstad, P., Nystad, W., Patelarou, E., Porta, D., Postma, D., Raaschou-Nielsen, O., Rudnai, P., Sunyer, J., Stephanou, E., Sørensen, M., Thiering, E., Tuffnell, D., Varró, M. J., Vrijkotte, T. G. M., Wijga, A., Wilhelm, M., Wright, J., Nieuwenhuijsen, M. J., Pershagen, G., Brunekreef, B., Kogevinas, M., and Slama, R.: Ambient air pollution and low birthweight: A European cohort study (ESCAPE), Lancet Resp. Med., 1, 695–704, https://doi.org/10.1016/S2213-2600(13)70192-9, 2013.
Piersanti, A., Vitali, L., Righini, G., Cremona, G., and Ciancarella, L.: Spatial representativeness of air quality monitoring stations: A grid model based approach, Atmos. Pollut. Res., 6, 953–960, https://doi.org/10.1016/J.APR.2015.04.005, 2015.
PMS5003 series data manual: https://aqicn.org/air/sensor/spec/pms5003-english-v2.3.pdf, last access: 2 February 2022.
Pope, C. A. and Dockery, D. W.: Health Effects of Fine Particulate Air Pollution: Lines that Connect, J. Air Waste Manage., 56, 709–742, https://doi.org/10.1080/10473289.2006.10464485, 2012.
Pope, C. A., Coleman, N., Pond, Z. A., and Burnett, R. T.: Fine particulate air pollution and human mortality: 25+ years of cohort studies, Environ. Res., 183, 108924, https://doi.org/10.1016/J.ENVRES.2019.108924, 2020.
Raaschou-Nielsen, O., Andersen, Z. J., Beelen, R., Samoli, E., Stafoggia, M., Weinmayr, G., Hoffmann, B., Fischer, P., Nieuwenhuijsen, M. J., Brunekreef, B., Xun, W. W., Katsouyanni, K., Dimakopoulou, K., Sommar, J., Forsberg, B., Modig, L., Oudin, A., Oftedal, B., Schwarze, P. E., Nafstad, P., De Faire, U., Pedersen, N. L., Östenson, C. G., Fratiglioni, L., Penell, J., Korek, M., Pershagen, G., Eriksen, K. T., Sørensen, M., Tjønneland, A., Ellermann, T., Eeftens, M., Peeters, P. H., Meliefste, K., Wang, M., Bueno-de-Mesquita, B., Key, T. J., de Hoogh, K., Concin, H., Nagel, G., Vilier, A., Grioni, S., Krogh, V., Tsai, M. Y., Ricceri, F., Sacerdote, C., Galassi, C., Migliore, E., Ranzi, A., Cesaroni, G., Badaloni, C., Forastiere, F., Tamayo, I., Amiano, P., Dorronsoro, M., Trichopoulou, A., Bamia, C., Vineis, P., and Hoek, G.: Air pollution and lung cancer incidence in 17 European cohorts: Prospective analyses from the European Study of Cohorts for Air Pollution Effects (ESCAPE), Lancet Oncol., 14, 813–822, https://doi.org/10.1016/S1470-2045(13)70279-1, 2013.
Rousseeuw, P. J.: Silhouettes: A graphical aid to the interpretation and validation of cluster analysis, J. Comput. Appl. Math., 20, 53–65, https://doi.org/10.1016/0377-0427(87)90125-7, 1987.
Samoli, E., Stafoggia, M., Rodopoulou, S., Ostro, B., Declercq, C., Alessandrini, E., Díaz, J., Karanasiou, A., Kelessis, A. G., Tertre, A. Le, Pandolfi, P., Randi, G., Scarinzi, C., Zauli-Sajani, S., Katsouyanni, K., Forastiere, F., Alessandrini, E., Angelini, P., Berti, G., Bisanti, L., Cadum, E., Catrambone, M., Chiusolo, M., Davoli, M., de' Donato, F., Demaria, M., Gandini, M., Grosa, M., Faustini, A., Ferrari, S., Forastiere, F., Pandolfi, P., Pelosini, R., Perrino, C., Pietrodangelo, A., Pizzi, L., Poluzzi, V., Priod, G., Randi, G., Ranzi, A., Rowinski, M., Scarinzi, C., Stivanello, E., Zauli-Sajani, S., Dimakopoulou, K., Elefteriadis, K., Katsouyanni, K., G.Kelessis, A., Maggos, T., Michalopoulos, N., Pateraki, S., Petrakakis, M., Sypsa, V., Agis, D., Alguacil, J., Artiñano, B., Barrera-Gómez, J., Basagaña, X., de la Rosa, J., Diaz, J., Fernandez, R., Jacquemin, B., Linares, C., Ostro, B., Pérez, N., Pey, J., Querol, X., Sanchez, A., Sunyer, J., Tobias, A., Bidondo, M., Declercq, C., Le Tertre, A., Lozano, P., Medina, S., Pascal, L., and Pascal, M.: Associations between fine and coarse particles and mortality in Mediterranean cities: Results from the MED-PARTICLES project, Environ. Health Persp., 121, 932–938, https://doi.org/10.1289/EHP.1206124, 2013.
Sayahi, T., Butterfield, A., and Kelly, K. E.: Long-term field evaluation of the Plantower PMS low-cost particulate matter sensors, Environ. Pollut., 245, 932–940, https://doi.org/10.1016/J.ENVPOL.2018.11.065, 2019.
Wang, Y., Li, J., Jing, H., Zhang, Q., Jiang, J., and Biswas, P.: Laboratory Evaluation and Calibration of Three Low-Cost Particle Sensors for Particulate Matter Measurement, Aerosol Sci. Tech., 49, 1063–1077, https://doi.org/10.1080/02786826.2015.1100710, 2015.
Wang, Z., Zhong, S., He, H., Peng, Z. R., and Cai, M.: Fine-scale variations in PM2.5 and black carbon concentrations and corresponding influential factors at an urban road intersection, Build. Environ., 141, 215–225, https://doi.org/10.1016/J.BUILDENV.2018.04.042, 2018.
Weinmayr, G., Romeo, E., de Sario, M., Weiland, S. K., and Forastiere, F.: Short-Term effects of PM10 and NO2 on respiratory health among children with asthma or asthma-like symptoms: A systematic review and Meta-Analysis, Environ. Health Persp., 118, 449–457, 2010.
Wenger, J., Arndt, J., Buckley, P., Hellebust, S., Mcgillicuddy, E., O'Connor, I., Sodeau, J., and Wilson, E.: Source Apportionment of Particulate Matter in Urban and Rural Residential Areas of Ireland (SAPPHIRE), Environmental Protection Agency, Ireland, ISBN 978-1-84095-905-5, https://www.epa.ie/publications/research/air/research-318.php (last access: 26 October 2023), 2020.
Wind module technical sheet: https://click.clarity.io/hubfs/Marketing%20Assets%20-%20PDFs/Product%20and%20Specification%20Sheets/2%20Pager%20Flyer%20%2B%20Specifications%20%E2%80%94%C2%A0Wind%20&%20Met%20Module.pdf, last access: 29 April 2024.
Zamora, M. L., Rice, J., and Koehler, K.: One year evaluation of three low-cost PM2.5 monitors, Atmos. Environ., 235, 117615, https://doi.org/10.1016/j.atmosenv.2020.117615, 2020.
Zhang, Y., Shi, Z., Wang, Y., Liu, L., Zhang, J., Li, J., Xia, Y., Ding, X., Liu, D., Kong, S., Niu, H., Fu, P., Zhang, X., and Li, W.: Fine particles from village air in northern China in winter: Large contribution of primary organic aerosols from residential solid fuel burning, Environ. Pollut., 272, 116420, https://doi.org/10.1016/J.ENVPOL.2020.116420, 2021.