Indoor Air Pollution from Residential Stoves: Examining the Flooding of Particulate Matter into Homes during Real-World Use

: This study concerns the levels of particulate matter (PM 2.5 and PM 1 ) released by residential stoves inside the home during ‘real world’ use. Focusing on stoves that were certiﬁed by the UK’s Department of Environment, Food, and Rural Affairs (DEFRA), PM sensors were placed in the vicinity of 20 different stoves over four weeks, recording 260 uses. The participants completed a research diary in order to provide information on time lit, amount and type of fuel used, and duration of use, among other details. Multivariate statistical tools were used in order to analyse indoor PM concentrations, averages, intensities, and their relationship to aspects of stove management. The study has four core ﬁndings. First, the daily average indoor PM concentrations when a stove was used were higher for PM 2.5 by 66.24% and PM 1 by 69.49% than those of the non-use control group. Second, hourly peak averages are higher for PM 2.5 by 55.34% and for PM 1 by 57.09% than daily averages, showing that PM is ‘ﬂooding’ into indoor areas through normal use. Third, the peaks that are derived from these ’ﬂooding’ incidents are associated with the number of fuel pieces used and length of the burn period. This points to the opening of the stove door as a primary mechanism for introducing PM into the home. Finally, it demonstrates that the indoor air pollution being witnessed is not originating from outside the home. Taken together, the study demonstrates that people inside homes with a residential stove are at risk of exposure to high intensities of PM 2.5 and PM 1 within a short period of time through normal use. It is recommended that this risk be reﬂected in the testing and regulation of residential stoves.


Introduction
As a component of air pollution, particulate matter with an aerodynamic diameter that is equal to 2.5 µm or less (PM 2.5 ) has long been linked to adverse health effects. In terms of mortality, it causes seven-million deaths per year [1]. In terms of health effects, it causes inflammation and oxidative stress, which compromises pulmonary immunity and increases the susceptibility to infection [2]. As these particulates can move into every organ in the body, the illnesses that are associated with their presence range from lung cancer, bronchitis, and other respiratory infections, through to strokes, dementia, and Parkinson's disease [3]. Effects such as these are particularly pronounced for children, pregnancies, and the elderly [4]. While much research focuses on particulate emissions that are generated by industry and vehicles, in the United Kingdom (UK) the primary source for PM 2.5 is which the stoves adhere, if at all, tend not to be detailed. This makes it difficult to generalise findings to categories of stove that share fundamental design features. Where stove standards are described, those chosen tend to have been approved by regulators outside the UK. For instance, [17,21,22,31] have focused on stoves that are approved by environmental regulators, but these are limited to the USA and Canadian contexts. Taken together, this relationship between indoor emissions and UK-specific regulations that govern stove design and testing requires investigation.
Fourth, few of the existing studies examine Ultra Fine Particles (UFP), which are defined as particles with a diameter of less than 100 nm, or Particle Number Concentration (PNC), which is defined as the total number of particles measured per cubic centimeter in a given sample. Measuring PNC along with the regular mass concentration measurements of PM 2.5 is important because PNC and PM 2.5 are not representative of each other [32], with Pearson's r lying between 0.09-0.64 and high levels of PM 2.5 not necessarily causing high levels of PNC or vice versa. Therefore, measures that are taken to reduce or regulate PM 2.5 may be different to those that are needed to tackle the problem of increasing PNC. Indeed, Penttinen et al. [33] found a stronger negative association between PNC and peak expiratory flow (PEF) than PM 2.5 amongst asthmatic children. Therefore, UFP may pose a substantial health risk since PNC exposure increases remarkably in the smallest size fractions.
When considering these limitations, this study has four aims. First, it seeks to determine real-world indoor PM exposure from the use of residential heating stoves over 30 days. This period was chosen to increase the number of uses from which data could be derived without instructing participants to use their stoves, minimise intrusion into the research setting, and more accurately capture 'real-world' use. Second, it detects and identifies the existence of peak indoor PM 2.5 and PM 1 levels as a result of stove use. Third, it seeks to clarify whether the level of indoor air pollution is originating from indoor or outdoor sources. Finally, it seeks to determine the extent to which these emissions are coming from a specific category of stoves; those that are certified as a 'Smoke Exempt Appliance' by the UK's Department for Environment, Farming, and Rural Affairs (DEFRA). These stoves are modified in order to restrict incoming air and limit smoke produced from combustion, differentiating them from the older equipment of focus in Semmens et al. [17]. If a stove passes the official testing process [34], they are certified to be exempt from the Smoke Control Area regulations covering most of the UK's towns and cities. However, this testing is limited to measuring outdoor air pollution via flue emissions and heat output; none of the applicable standards that are required by DEFRA are concerned with indoor PM emissions from stoves (see PD 6434: 1969; BS 3841: Part 1: 1994; BS 3841: Part 2: 1994). Even the latest 'EcoDesign' standards, which call up EN 16510:2018, do not introduce testing for indoor emissions. Indeed, when taken together, the DEFRA testing regime rests on a baseline assumption that stoves do not pollute indoors, or only do so when a fault is present. The results of this study test the validity of that foundational assumption. Taken together, this work makes three core contributions: i It presents a framework in order to determine real-world indoor PM exposure from the use of residential heating stoves. ii It can detect and identify the existence of peak indoor PM 2.5 , PM 1 , and PNC levels as a result of stove use. iii It analyses the results in relation to the DEFRA regulations and determines the extent of these emissions from a specific category of stoves; those that are certified as a 'Smoke Exempt Appliance' by DEFRA.
In making these contributions, the study seeks to determine whether health risks are posed during normal operation and, in turn, whether DEFRA testing standards need modification in light of this reality.
The remainder of this paper is organised, as follows. Section 2 describes the experimental framework along with sensor calibration and evaluation in Section 2.2. Section 3 presents the findings and analysis, which is followed by the conclusion in Section 4.  1 The stoves were installed in a house but used under controlled conditions; 2 280 uses in total but 20 removed due to incomplete data.

Sampling Area and Study Design
Sheffield (53 • 23 ′ N 1 • 28 ′ W), the chosen study site that is shown in Figure 1a, is a geographically diverse city that is located in the county of South Yorkshire, England. Built on several hills, it is situated at an elevation of 29 m-500 m above sea level, covers a total area of 367.9 km 2 , and it has a growing population of 582,506 [36]. Sheffield has a temperate climate; July is considered to be the hottest month, with an average maximum temperature of 20.8 •C and January-February to be the coldest months. Air pollution in the city is primarily from road transport and industrial emissions and, to a lesser extent, fossil fuels run processes, such as energy supply and commercial or domestic heating systems [37]. Twenty households with solid fuel stoves were recruited between January and April 2020. An indoor and outdoor low-cost air quality monitor was installed in each of the houses. The indoor sensor was placed at a minimum of 3 m distance from the wood burner for safety, but in the same room. The outdoor unit was put in weatherproof casing and then attached to a window or drain pipe outside of the house (see Figure 1c). Data from each household were collected over a total period of four weeks. Data were recorded on days when stoves were used and left unused, thus providing two groups of data. The control group contained 10 users who had stoves and, over a 30 day period, used them around 30% of the time. Control group data were taken from 20 days of non-usage. In total, 10 out of the 20 participants were identified as the control users for the study.
Pollutants that were measured real-time for both indoors and outdoors were PM 10 , PM 2.5 , PM 1 , PNC (0.3 µm-1 µm), Nitrogen Dioxide (NO 2 ), Carbon Monoxide (CO), and Ammonia (NH 3 ). The meteorological parameters include temperature, Relative Humidity, and Atmospheric Pressure. The data were sampled every 145 s. Data for NO 2 , CO, and NH 3 were omitted for research purposes and only visualised as trend levels due to the lack of calibration instruments. For indoor air pollution levels, the focus of our analysis was PM 2.5 and PM 1 .
One participant from each household completed a survey prior to the measurement period and maintained a research diary throughout the study. Among other data, the research diary recorded stove usage timings, indicating when the stove was lit and when the last piece of fuel was added, type and total amount of fuel, and type and total amount of kindling used each time the stove was active. Any other activities carried out during stove use, such as cooking or lighting of candles, was also recorded. The air pollution level indoors was calculated between the time that the stove was lit until one hour after the last piece of fuel was added. This was done to allow for the complete combustion of the fuel that was fed to the stove. Each participant was provided with a tablet computer. This displayed a dashboard containing real-time information on indoor and outdoor pollution levels that were collated from their sensors. A state-of-the-art cloud-based dashboard was built for each participant, as shown in Figure 2. The data from the monitoring units were sent to the cloud based server that was hosted by the University of Sheffield, which was then displayed on the dashboard. The information refreshed by default every minute. The graph panels plotted over a period of 30 days displayed daily average, minimum, and maximum values of each pollutant. Real-time sensor readings were also made available in the form of dynamic gauges.

Sensor Validation and Correction: Accuracy, Evaluation and Limitations
The Urban Flows Observatory [38] at the University of Sheffield have developed Enviro+ (Figure 1b), an air quality measurement device, in collaboration with Pimoroni, which is a local electronics company. Enviro+ is a pHat, which is an add-on board that sits on top of raspberry pi Zero and is suitable for both indoor and outdoor air quality measurement. Sensors onboard this pHat include a BME280, which is a weather sensor monitoring temperature, pressure. and relative humidity, an LTR-559 light and proximity sensor, a MICS6814 analog gas sensor monitoring NO 2 , CO and NH 3 , ADS1015 analog to digital converter (ADC), a MEMS microphone for noise measurement, and a 0.96" colour LCD (160 × 80) for display. A connector for a particulate matter (PM) sensor is also available onboard, to which was connected the low-cost optical sensor PMS5003 (Plantower) Enviro+ with the connected PMS5003, which was used to conduct the particulate level measurement. Enviro+ with the connected PMS5003 was housed in a casing and installed outside the house for outdoor air pollution measurements ( Figure 1c).
All of the units were collocated with Sheffield City Council's Reference Air Quality Monitoring station at Lowfield four weeks prior to the study. The high end Palas Fidas 200 instrument installed at Lowfield Station by Sheffield City Council was used as a reference in order to correct the PMS5003 sensors PM 2.5 measurement.
The procedure for correction of the collocated sensors is discussed below:

1.
Raw data, including PM 2.5 , PM 10 , Temperature (T), and Relative Humidity (RH), were received every 160 s. This was converted to hourly averages in order to match the reference station data, because only hourly reference data are publicly available.

2.
The hour average was excluded if less than 90% of the measurements were available in that hour average.

3.
Humidity Correction: PM 2.5 concentrations can be relatively high from low-cost PM sensors at high RH levels. The hygroscopic growth of particles at high humidity, along with mist and fog particles, makes it detectable as particulates, as previously reported [39,40]. A Nephelometer, such as PMS5003, measures particulates based on light scattering principle. The particulates' refractive indices are dependent on relative humidity [41] and, thus, affects the sensor readings. While ambient temperature directly has a very limited role in sensors performance [40] (apart from extreme temperature), it affects the measurements indirectly. Jayaratne et al. [39] reports that, when the ambient temperature reaches the dew point temperature, the conditions become suitable for the formation of fog droplets in the air and fall within the detection size of such sensors. Figure 3 presents an example of the relationship between RH values and PM2.5 data from PMS5003 collocated. A Humidity-based bias correction approach was taken, as described here [42], while using the κ -Köhler theory [43]. The hygroscopic growth factor g (RH), as defined in Equation (1), where D dry is the diameter of the dry particle and D wet (RH) is the diameter of the particle at a given RH value.
RH dependence [44] was established while using Equation (2), as follows: where κ is a parameter that describes the degree of hygroscopicity of a particle and taken as 0.62, which is suitable for Sheffield [45]. Therefore, using Equations (1) and (2), hygroscopic growth factor g (RH) was calculated in order to obtain the humidity correction factor. Two additional PMS5003 have also been collocated at the same station permanently since 23rd April 2019 have been used to ensure correction factor accuracy. A conditional Quartile plot in Figure 4 below uses the corresponding values for both reference and low cost sensors, splitting the values into evenly spaced bins. For each low cost sensor value bin, the corresponding reference sensor values are identified and the median, 25/75t,h and 10/90 percentile (quantile) are calculated for that bin. The data are plotted in order to show how these values vary across all bins. The blue line shows the results for a perfect model i.e., zero error between low cost PMS5003 sensor and the reference Palas FIDAS 200 sensor. In the plot in Figure 4, the red line shows that the LCS tends to slightly over-report for PM 2.5 (NMB ≈ 0.2-0.3).

4.
Concentration Range Correction: a correction was applied based on the relationship between pollutant concentration range and sensor performance. Multivariate Linear regression model were used in order to establish the relationship. Palas Fidas 200: PM 2.5re f is used as the dependent variable and PMS5003 sensor data: PM 2.5lcs , T, and RH as predictors, as shown in Equation (3).
β 0 , β 1 and β 2 are calculated by training with the model generated. To note, β 3 is not used here, as it is obtained from the previous step.

5.
Evaluation of LCS: PMS5003 corrected data are evaluated by comparing to the Palas Fidas 200 values in the holdout data set. From the field evaluation through collocation between January-April 2020, PMS5003 showed high linear correlation with reference instrument with R 2 value 0.81 for the hourly averaged data. This is an improvement in accuracy when compared to the findings from previous studies on evaluating Plantower sensors [10,46] with R 2 values lying between 0.71-0.77 for PMS5003 without applying any correction factors. The inter-sensor comparison showed a high correlation, with an R 2 value between 0.98-0.99. Figure 5, below, shows the scatter plot between the reference and corrected PMS5003 sensor. Figure 6, below, also shows a consistently high linear correlation factor with an average R 2 value of 0.81 when analysed and split with relative humidity as the third variable.   Table 2, below also shows the R 2 value for different concentrations of PM 2.5 and compared to the Daily Air Quality Index (DAQI) bands and breakpoints for PM 2.5 , as set by DEFRA'S Air Quality Expert Group [47]. During the field evaluation, there was not enough data to evaluate the sensor for high and very high conditions (DAQI = 8-10).

Sensor Limitations
The reference station does not provide PM 1 and PNC data and, therefore, this data cannot be subjected to this correction. Further research is underway in order to evaluate sensor performance and evaluation in this specific regard. Finally, the study has not been able to account for UFP due to these sensors being unable to detect or measure particles below 300 nm. As such, the measured PNC has been limited to a size of 0.3 µm-1 µm.

Monitoring Outdoor Air Quality and Adjusting for Weather: A Generalized Boosted Regression Model
While the data that were collected from the indoor unit ( Figure 1b) were used to analyse the pollution emissions from the stoves, the first purpose of the outdoor unit (Figure 1c) was twofold. First, it was used for the general monitoring of outdoor PM levels. This allowed for the detection of any unusual levels of outdoor pollution that could impact the air quality indoors. Second, the sensors could indicate whether the outdoor air quality was also being influenced during stove use. While the outdoor sensors served the first purpose, achieving the second was complicated by multiple covariates, such as meteorological factors, local garden waste burning, neighbours using wood stoves, and traffic. Figure 7 plots the average weekly variation of outdoor PM 2.5 and PM 10 levels of the participants houses over the three-month period. Meteorology plays a crucial role in the estimation of levels of particulate matter. Therefore, when trying to understand the trends of outdoor pollution levels, it can be very challenging to determine whether a pollution episode is caused by local emissions or meteorology. Therefore, a Machine Learning (ML) based algorithm based on Generalised Boosted Regression Model [48] was used in order to explore and adjust for the non-linear relationships between the meteorological covariates and particulate matter PM 2.5 levels. The partial dependencies in Figure 8 show the relationship between PM 2.5 and the covariates that were used in the model while holding the value of other covariates at their mean levels. As can be seen, wind speed (16.1%) and wind direction (12.7%) play a crucial role in determining PM 2.5 levels; hence, its impact should be accounted for in order to better understand the air quality around the participating households. The popular R deweather and openair package [49] was used in creating the prediction model and plotting. The model is formed, as shown in Equation (4).
[PM 2.5 ] =RH +ū + φ + Tθ + t hour + t weekday + tJD (4) whereū is the mean hourly wind speed, φ is the mean hourly wind direction (degrees, clockwise from the north), and Tθ is the mean hourly temperature ( • C). Variables representing hour of the day, t hour , day of the week, t weekday , and day of the year, tJD were also considered for the model development. From Figures 7 and 8, it is evident that, during weekdays, the outdoor levels of PM 2.5 and PM 10 are higher than during the weekend. It can also be seen that the levels are considerably higher outside during the evening, which corresponds to the usage pattern of stoves by participants. This indicates that even DEFRA-certified solid fuel stoves could affect the local air quality outdoors. A more sophisticated source apportionment study is required in order to further investigate this. The high level of (PM 2.5 and PM 10 gradually decreases throughout the night, with the lowest levels being attained at around 5:30 am-6:00 am GMT (see Figure 7 hourly plot). Ten-fold cross validation [50] was used for evaluating the model performance and the model fitting results are shown in Figure 9.

Data Processing and Storage
The real-time sensor data are pushed to a cloud-based database over WiFi while using a Python script running on the Raspberry Pi Zero. The data are stored in a database designed and installed in a virtual server hosted by the University of Sheffield. These data are then made available to be accessed through an Application Programming Interface (API), called Enviro-API, developed for data retrieval and displayed on the dashboard. The API developed along with the database had two goals:

•
The ability to ingest a high volume of time series data with dynamic data from the sensors.
• The ability to return this time series data with basic querying parameters such as sensor ID and timestamps.
This allowed for us to create an end-to-end secure infrastructure for real-time sensor data collection, storage, and retrieval system for our study.

Data Analysis
Missing data have been treated. The usage days were only included if 90% of the hourly data were available. Data analyses were performed while using Excel, R, and Python programming languages.
The statistical significance of the results was calculated based on Welch's t-test (Moser & Stevens, 1992) whlie using the standard equations: and the degree of freedom of Welch t-test is calculated, as follows: • A and B represent the Control and Experimental group. , and strong (±0.5-1). Data were removed during such periods while stoves were lit in order to avoid data being influenced by emissions from cooking, burning candles, or incense sticks. The influence of outdoor air pollution on indoor emissions data was anticipated, but adjustment was unnecessary, due to the absence of notable outdoor pollution levels.

Study Limitations
The study exhibits several limitations that are associated with variability in the research setting due to its exploratory design and focus on real-world stove use. First, the study does not account for the impact of room size, seal, ventilation, and dwelling age on the duration of air pollution exposure witnessed. Nor does it relate the levels of air pollution to specific stages of the combustion cycle. Further study is needed in order to understand these aspects of indoor air pollution, requiring a sampling frame that is determined by more than the stove type and a research design that is appropriate for lab conditions. Second, despite using outdoor sensors to illustrate that the indoor air pollution is not coming from outside sources (see Section 3.2), further details on air pollution at the indoor-outdoor interface were beyond the design of this study. This is a characteristic of air pollution research more broadly, as reflected in the UK government's recent multi-million-pound call for research that is able to develop solutions to air pollution problems at the indoor/outdoor interface [51]. Relatedly, windspeed could influence the infiltration rate of outdoor air indoors, but, again, this was beyond the remit here. As such, further research into this relationship is recommended. Finally, the influence of sensor data on participant stove management practice has not been explored in detail. This will be drawn out more fully in a separate paper. Table 3 summarises the daily PM 2.5 and PM 1 mean, and hourly peak PM 2.5 and PM 1 mean from 20 households and 260 stove usages, along with the statistical analysis and distribution. Data on the average pieces of fuel per use (FP) and kindling per use (KP), along with the average duration of use, have also been presented. The hourly indoor mean PM 2.5 and PM 1 concentrations that were observed during stove usage ranged from 2.27 µg/m 3 and 1.11 µg/m 3 to 47.60 µg/m 3 and 36.15 µg/m 3 , respectively, with a high coefficient of variation 0.9 for PM 2.5 and 0.94 for PM 1 . The hourly PNC average that was observed indoors in the particle size range (0.3-1 µm diameter) was 2607 particles/0.1 litre (L) of air when each stove was used, but the hourly peak PNC average observed was 4345 particles/0.1 L with an hourly maximum of 9978 particles/0.1 L. The average number of fuel pieces (9.58 wooden logs) and kindling (8.37 pieces) used varied significantly between the households, with a coefficient of variation 0.69 and 0.67, respectively. The average duration of use was approximately 4 h, with most households using their stove between 6 pm and 10 pm.
The analysis in the three quartiles-(i) <25 percentile, (ii) >25 <75 percentile, and (iii) >75 percentile, representing low, medium, and peak concentrations, showed an increase for PM 2.5 (69.12%, 70.69%, 56.10%) and PM 1 (71.78%, 70.41%, 67.67%). The overall average concentrations were higher for PM 2.5 by 66.24% and PM 1 by 69.49% when used. Figure 10 density rug plots show the distribution and levels of PM 25 and PM 1 for users. Figure 11 compare the control group's indoor pollution levels with the experimental group. For reasons of visualisation, scaling the x-axis in the graph (see Figure 10) is limited to 60 µg/m 3 . Figure 10 reveals that the levels of PM that people are exposed to can vary, with a maximum peak average of 47.60 µg/m 3 for PM 2.5 and 36.15 µg/m 3 for PM 1 . While calculating the averages smooths the graph, these findings demonstrate that some users are exposed to maximum values of up to 160 µg/m 3 PM 2.5 . Control users experience much lower indoor particulate levels when their stoves are not lit when compared to users that do, as indicated by Figure 11.  Wood burner usage data from 20 households collected between January and April 2020.
In Figure 10, comparing the concentration levels between usage and non-usage days for the control group also illustrates an increase for PM 2.5 (58.24%, 76.60%, 76.22%) and PM 1 (56.92%, 80.51%, 78.56%) when stoves are used. The overall average concentrations were higher for PM 2.5 by 81.23% and PM 1 by 73.76%.   Figure 12, below, it is clear that indoor and outdoor values vary significantly between 10-45 µg/m 3 concentration levels. This variation is because the mean and hourly peak indoor PM lies within this range and, thus, the indoor levels are much higher than the corresponding outdoor levels. Further analysis of average indoor and outdoor PM 2.5 levels indicated a weak correlation (R 2 = 0.19) between them, which suggests that outdoor air quality is not a driving factor behind the high indoor pollution levels that were seen during stove usage.

Indoor Outdoor
While we acknowledge that indoor PM 2.5 levels can impact outdoor air quality, no measurements were taken from the chimney/flue. The air quality sensor outside the house indicates immediate outdoor air pollution levels and, thus, it is difficult to measure any leakage at the interface. Future research studies should focus on indoor air pollution and its influence on outdoor air quality in order to address this limitation of our study.

Hourly Peak PM Average Higher than Daily PM Average
The analysis of Table 3 shows hourly peak PM 2.5 and PM 1 is strongly correlated with daily mean PM 2.5 and PM 1 (r = 0.75). Statistical analysis shows that the hourly peak mean PM 2.5 (27.34  There exists high variation in exposure concentrations, concerning both short peaks and daily levels. This characteristic is related to the "real-world" nature of the study. The research diary tool provided data on not only the amount of fuel and kindling pieces used, but also their type. On average, participants used 9.58 pieces of solid fuel and 8.32 pieces of kindling per use. The number of fuel pieces used varied between a minimum of seven to a maximum of 40, while kindling varied between a minimum of one and a maximum of 32. All participants used dried and seasoned logs, but the sizes varied. There was also a diversity of kindling used, taking the form of firelighters, newspapers, balls of paper, twigs, sawdust, packing cardboard, greeting cards, and even empty egg boxes. Echoing the findings of existing studies [20,35,52]. This means that the same wood burner may emit different levels of indoor air pollution depending on the quantity and type of fuel and kindling used. While suggesting a link between indoor air pollution and fuel quantity, and type of fuel and kindling, following other studies in the next section, demonstrates that this is actually linked with the stove door being opened.
Epidemiology studies and policymaking are focused around hourly average concentration monitoring by regulatory air quality stations. This leads to the omission of short-term high exposure through the "flooding" of indoor spaces with PM 2.5 and PM 1 . Very few studies have reflected on short term peak concentration exposure. Lin et al.'s study [27] associated increased risk factors with hourly peak concentrations of PM 2.5 . Similarly, Delfino et al. [53] associated peak PM levels with Asthma attacks in children, but in outdoor environments. Therefore, the present study encourages future researchers to study the occurrences and effects of relatively short-term peak PM exposure on human health.

Hourly Peak Average PM Has a Moderate Correlation to the Pieces of Fuel Used
While Table 3 indicates a weak correlation between fuel pieces and mean PM 2.5 (r = 0.17), and with PM 1 (r = 0.15), comparing the hourly peak concentration of PM 2.5 (r= 0.44) and PM 1 (r = 0.43) exhibits a moderate correlation with fuel pieces. The scatter plots in Figures 13 and 14 chart the relation between peak hourly levels to the possible co-factors of fuel amount and duration of usage. In Figure 13a While correlation between fuel pieces and hourly mean concentration is weak, it is stronger when compared to the hourly peak concentration. Therefore, the findings suggest that the peak hourly concentrations are often higher by a minimum of 250% and a maximum of 400% when participants have refuelled their stove more than once during a usage compared to one refuel or none at all. As such, the findings indicate that the 'flooding' of indoor space occurs as a result of the stove door being opened for refuelling. This accords with several existing real-world [21,22,35], and lab-based [20] studies into stoves outside the UK. While the findings point to the opening of the stove door as the origin for indoor PM emissions, further lab-based research is required into how this might relate to duration, timings, and the point in the burn cycle at which the opening occurs. The hourly peak concentrations explain the shape of the rug plots, as seen in Figure 10. The shape of the curves exhibit a distinct broad frequency distribution in the lower PM concentration. This indicates that most of the sensor readings are lower during stove use, but there are also smaller spikes towards the right of x-axis, indicating sensor readings that correspond to higher levels of PM pollution. A 'leakage' would result in a more uniform shape, and, thus, the presence of the smaller spikes cannot be explained. This echoes Salthammer et al.'s findings [35] and provides further support for the theory of opening doors being the cause of the indoor air pollution seen rather than a leakage, which appears to be more common to open fires than 'closed' stoves (see [54]). The PM fraction gets dispersed quickly throughout the room due to its smaller size, reverting to lower hourly average concentrations. Table 3 also illustrates a non-linear relationship between the duration of use and mean PM 2.5 (r = 0.017). This is similar to PM 1 (r = 0.021), although, again, comparing the hourly peak concentrations of PM 2.5 (r= 0.4) and PM 1 (r = 0.38), it exhibits a moderate correlation with the duration of use. The scatter plots in Figure 14a,b also reflect this, with higher levels of peak values being continuously registered during the stove use.

Hourly Peak Averages Illustrate a Moderate Correlation with Duration of Use
Longer usage is associated with greater numbers of fuel pieces used. This result supports the explanation for the 'flooding' phenomenon observed, with higher short-term peak concentrations being seen during longer periods of use, because these periods are sustained by more refueling actions. This accords with [20], who also found the lighting and refueling aspects of stove management to form the main pollutant-generating phases of operation.

Conclusions
The present study aimed to understand the extent to which PM was emitted indoors and under real-world conditions by DEFRA-certified residential stoves. The findings indicate that real-world indoor PM exposure from these stoves is higher when lit as compared to the period in which they are not in use. When compared to periods of non-use, the overall average concentrations were higher for PM 2.5 by 66.24% and PM 1 by 69.49%. Peak hourly concentrations of PM were often found to be higher by 250-400% when the participants had refueled their stove more than once during a single usage. The findings also provide information on PNC, with an average hourly peak of 9978 particles/0.1 L emitted during a single usage. These 'flooding' events correlated with the opening of the stove door, which indicated that such incidents occurred as fuel was added. Data from outdoor sensors clarified that this was not originating from outdoors. On the basis of these results, it is recommended that DEFRA testing standards be modified in order to account for these normative health risks. The PM that is released into the home is not an aberration from normal use, but results directly from it. This is because real-world operation cannot occur without opening the stove door. It may be that with regulatory encouragement stove designs can be modified in a way that limits such instances. In the meantime, or in the event that appropriate modification cannot be achieved, it is also recommended that new residential stoves be accompanied by a health warning at the point of sale in order to indicate the normative health risks posed to users.