Next Article in Journal
Psychosocial Resources for Hedonic Balance, Life Satisfaction and Happiness in the Elderly: A Path Analysis
Previous Article in Journal
Ecological Effects of Heavy Metal Pollution on Soil Microbial Community Structure and Diversity on Both Sides of a River around a Mining Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Liquefied Petroleum Gas Cookstoves and Wood Cooking Fires on PM2.5 Trends in Brick Workers’ Homes in Nepal

1
Department of Public Health, Brigham Young University, Provo, UT 84602, USA
2
Department of Plant and Wildlife Sciences, Brigham Young University, Provo, UT 84602, USA
3
Marriott School of Business, Marketing & Global Supply Chain, Brigham Young University, Provo, UT 84602, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(16), 5681; https://doi.org/10.3390/ijerph17165681
Submission received: 14 June 2020 / Revised: 4 August 2020 / Accepted: 4 August 2020 / Published: 6 August 2020
(This article belongs to the Section Environmental Health)

Abstract

:
Prior studies document a high prevalence of respiratory symptoms among brick workers in Nepal, which may be partially caused by non-occupational exposure to fine particulate matter (PM2.5) from cooking. In this study, we compared PM2.5 levels and 24 h trends in brick workers’ homes that used wood or liquefied petroleum gas (LPG) cooking fuel. PM2.5 filter-based and real-time nephelometer data were collected for approximately 24 h in homes and outdoors. PM2.5 was significantly associated with fuel type and location (p < 0.0001). Pairwise comparisons found significant differences between gas, indoor (geometric mean (GM): 79.32 μg/m3), and wood, indoor (GM: 541.14 μg/m3; p = 0.0002), and between wood, indoor, and outdoor (GM: 48.38 μg/m3; p = 0.0006) but not between gas, indoor, and outdoor (p = 0.56). For wood fuel homes, exposure peaks coincided with mealtimes. For LPG fuel homes, indoor levels may be explained by infiltration of ambient air pollution. In both wood and LPG fuel homes, PM2.5 levels exceeded the 24 h limit (25.0 µg/m3) proposed by the World Health Organization. Our findings suggest that increasing the adoption of LPG cookstoves and decreasing ambient air pollution in the Kathmandu valley will significantly lower daily PM2.5 exposures of brick workers and their families.

1. Introduction

Household air pollution from indoor burning of solid biomass and gaseous fuels is responsible for approximately 3.8 million deaths annually worldwide [1]. Indoor burning of solid biomass fuels, including wood, on open fires produces high levels of air pollutants, including inhalable (PM10) and fine (PM2.5) particulate matter, carbon monoxide (CO), nitrogen oxides (NOx), and an array of chemicals that are toxic or carcinogenic to humans [2,3,4,5,6]. Chronic exposure to household air pollution is associated with acute lower respiratory tract infections, respiratory illness, impaired immune function, and low birth weight in children, and chronic obstructive pulmonary disease (COPD), respiratory illnesses, tuberculosis (TB), impaired immune function, cardiovascular disease, cataracts, and lung cancer in adults [5,6,7]. Efforts to reduce exposure to household air pollution by providing improved stoves to exposed populations have been met with mixed results, largely due to cultural barriers, cost, and local social and environmental factors [8,9,10,11]. Thus, understanding population-specific exposure patterns, fuel use, and cooking and heating practices is a necessary first step in developing tailored, culturally appropriate interventions to reduce exposure.
Brick manufacturing is a major industry in Nepal and a potentially large source of household air pollution exposure within the country. Estimates of the number of operating brick kilns in Nepal range from 429 to 750 [12,13,14], with a workforce estimate of up to 400,000 seasonal workers annually [13,15]. In the Kathmandu valley alone, there are over 100 brick kilns that employ approximately 30,000 workers in total [16,17,18]. Brick workers in Nepal experience high rates of respiratory symptoms, including chronic cough, chronic phlegm, chronic bronchitis, wheezing, asthma, dyspnea, and chest tightness [19]. These respiratory symptoms may be partially due to occupational exposures, particularly for dusty tasks such as moving red bricks from the kiln after they are fired [20]. Brick workers’ exposures to silica and other hazardous dusts are well documented globally [21,22,23]. In Nepal, it is common for brick workers to live in impoverished conditions, often on-site at the brick kiln. Brick workers and their families typically live in rudimentary huts composed of brick walls and a tin roof, with an average floor area of less than 8.0 m2 [24]. During mealtimes, some brick workers and their families rely on wood fuel for cooking, which may emit dangerous levels of particulate matter and other pollutants inside the home. Thus, household air pollution exposure during non-working hours may compound brick workers’ occupational exposures and may partially explain the high rates of respiratory symptoms in this population.
To date, there are few studies on indoor air pollution exposure among brick workers and their families globally or specifically in Nepal. One previous study evaluated indoor PM2.5 exposure in brick workers’ homes in Bhaktapur [24]. However, exposure monitoring in that study occurred during the middle of the day for approximately seven hours, during which time most of the residents were not home. As a result, PM2.5 levels were not measured during peak cooking hours, and results showed no significant difference between indoor and outdoor concentrations [24]. The two primary cooking methods used among brick workers and their families in Nepal are cook stoves using liquefied petroleum gas (LPG) and open fires using wood. Given the relatively low cost and global availability of LPG, the World Health Organization (WHO) has stated that studies of the impact of LPG on reducing household air pollution exposure are needed as soon as possible [25]. Therefore, we aimed in this study to characterize indoor PM2.5 levels in LPG vs. wood fuel homes of brick workers and their families over a 24 h period. Understanding exposure patterns related to fuel use may help to guide future interventions to reduce PM2.5 exposure and respiratory symptoms in this vulnerable population of workers and their families.

2. Materials and Methods

2.1. Study Design

We used a cross-sectional study design to compare indoor PM2.5 levels, air temperature (°C), and relative humidity (RH; %) in brick workers’ homes according to the primary fuel used for cooking. Homes were classified as either LPG cookstove or woodfire homes. Figure 1 shows the typical construction of homes sampled in this study. Figure 2 illustrates LPG and wood fire cooking arrangements observed in this study. Homes (N = 19) were selected by convenience sampling at a single brick kiln in Bhaktapur, Nepal. All study homes were located at the worksite, within an approximately 300 m radius of the brick kiln. Samples were collected from 30 April to 3 May 2019. Measurements were collected by placing the monitoring devices on a tripod positioned at 1.2 m from the floor, which we estimated to be the approximate breathing zone height of an adult crouching in the home. For comparison, we simultaneously collected outdoor PM2.5, air temperature, and RH measurements away from observed emission sources and within 200 m of the study homes. A brief housing survey was also administered to an adult occupant of the home by means of an interpreter. Brigham Young University’s Institutional Review Board (IRB) reviewed this study. This study did not meet the definition of human subjects research outlined in 45 CFR 46 [26] because the unit of study was the home rather than the individual and, therefore, it was determined that the study did not require IRB approval.

2.2. Indoor PM2.5 Measurement

All PM2.5 samples were collected using RTI International’s MicroPEMs V 3.2A (RTI International, Research Triangle Park, NC, USA). The MicroPEMs allow for both gravimetric (filter-based) sampling and real-time datalogging using an on-board nephelometer. The MicroPEM nephelometers were set to record PM2.5 concentrations every 10 s over the sampling period. Gravimetric PM2.5 samples were collected on 3.0 μm PTFE 25 mm filters (Zefon International, Ocala, FL, USA). We pre-weighed filters using a Mettler Toledo (Mettler Toledo, Columbus, OH, USA) model XP2U microbalance in a temperature and humidity-controlled room. Filters were conditioned in the room for 24 h prior to weighing. We used Teflon™ coated tweezers (Mettler Toledo) to handle filters. We used a Haug (Haug North America, Williamsville, NY, USA) U-bar deionizer to remove static electricity from each filter prior to weighing. We weighed and deionized each filter 3 times. The mean of these weights was recorded as the pre-weight value. Filters were then placed in SKC filter keepers (SKC, Inc., Eighty Four, PA, USA) until use. The same procedure was followed after returning from Nepal for determining post-weights.
Prior to sample collection each day, we used Docking Station software version 2.0 (RTI International, Research Triangle Park, NC, USA) to interface with the MicroPEMs to set the instrument date/time, flow rate, and nephelometer offset. Nephelometer offsets were adjusted by placing a high efficiency particulate air (HEPA) filter in-line before the flowmeter to provide particle free air. The offset for the 780 nm IR laser was adjusted until the particle count was zero. A TSI model 4140 mass flowmeter (TSI, Shoreview, MN, USA) was used to set the MicroPEM flowrate at 0.50 L/min. Sintered stainless steel impaction plates were cleaned and oiled by RTI International prior to data collection. MicroPEMs were attached to tripods by means of a pouch made of low-dander ripstop nylon. After sample collection, nephelometer data were downloaded from the MicroPEMs using Docking Station. MicroPEM flow rates were checked after sampling, and all instruments were found to have minimal drift of less than ± 5% from 0.50 L/min and were thus determined to be valid samples [27]. Finally, the memory was cleared from each instrument before preparing the instrument for the next sampling event. For all homes sampled, the median sampling time was 21.21 h (interquartile range = 2.04 h).

2.3. Indoor Temperature and RH Measurement

Air temperature and RH data were collected in each home using Extech SD500 dataloggers (Extech Instruments, Nashua, NH, USA). Prior to data collection each day, we installed fresh batteries, cleared the SD card, set the date and time, and set the logging interval to record measurements every five minutes. Following the sampling period, data were downloaded from the instruments and saved as Excel spreadsheets. For all homes, the median sampling time was 21.42 h (interquartile range = 1.83 h).

2.4. Outdoor PM2.5, Temperature, and RH Measurement

For comparison with the indoor samples, we collected daily outdoor PM2.5, temperature, and RH measurements in a centralized location at the brick kiln, away from any observed pollution sources. Outdoor samples were collected under a covered pavilion on the kiln property. The pavilion roof was approximately 8 m high, with open sides, allowing for natural air flow. An RTI MicroPEM and Extech SD500 datalogger were positioned at 1.2 m from the ground on a tripod under the pavilion. For outdoor samples across all days, the median sampling time for the MicroPEMs was 21.98 h, and the mean and standard deviation (SD) for the Extech SD500 dataloggers were 20.62 and 1.31 h, respectively. We handled all instruments and data with the same methods as those described for the indoor samples.

2.5. Housing Questionnaire

Housing factors were assessed using an extant 14-item questionnaire [24] that was modified for use in this study. Specifically, we limited the number of questions to five, including number of people living in the home, number of children living in the home, primary fuel used for cooking, smokers in the home, and, if so, the number of smokers. Fuel and cooking device were verified by visual inspection of the home. While questionnaires were administered, study personnel measured the living area (m2) of the home. Occupant density was calculated as the number of people living in the home divided by the living area. Study personnel filled out the questionnaires by interviewing an adult occupant of each home with the help of an interpreter.

2.6. Statistical Analyses

Prior to conducting statistical analyses, we corrected the MicroPEM nephelometer readings based on the gravimetric results from the on-board filter, as recommended by RTI International [28]. The correction factor was calculated as the filter concentration divided by the mean nephelometer concentration over the sampling period. We then multiplied all of the MicroPEM nephelometer PM2.5 data by the correction factor.
We conducted all statistical analyses using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA) and used a significance level of α = 0.05. We calculated frequencies and percentages for categorical brick kiln home characteristics and the arithmetic mean, standard deviation, minimum, first quartile, median, third quartile, and maximum for continuous home characteristics and for mean (over the sampling period) PM2.5 concentration, RH, and temperature.
We used simple linear regression models to estimate unadjusted associations, with 95% confidence intervals, between brick kiln home characteristics and three primary outcome variables: mean (over the sampling period) PM2.5 concentration, RH, and temperature. The distribution of mean PM2.5 concentration was right-skewed, so we natural logarithm transformed the mean PM2.5 concentrations before including that variable as the outcome in linear regression models and then exponentiated the regression coefficients. We evaluated pairwise differences in the outcome variables among categories of fuel type and location and used the Tukey method to adjust p-values for multiple comparisons. Additionally, we made box plots that showed mean PM2.5 concentration, RH, and temperature according to fuel type and location. We also fit some multiple linear regression models that contained two home characteristics as independent variables.
To assess the diurnal variation of PM2.5 concentration, RH, and temperature, we made line graphs of mean (over fuel type and location categories) PM2.5 concentration, RH, and temperature over time.

3. Results

Samples were collected over a three-day period, during which we sampled 19 homes and collected daily outdoor measurements from a central location at a single brick kiln in Bhaktapur, Nepal. The median home area was 10.25 m2, mean number of people in the home was 3.42, and median occupant density was 29.70 people per 100 m2 (Table 1). Fifty-three percent of homes had 0–1 child in the home, while 2–3 children were present in 47% of homes. Smokers were present in 53% of homes and the mean number of smokers in the home was 0.76. Sixty-seven percent of homes used LPG, while 33% used wood as the primary fuel for cooking.
The median average (over the sampling period) PM2.5 concentration for all samples combined was 118.46 μg/m3 (Table 2). The mean average RH was 57.87% and the mean average temperature was 24.58 °C for all samples combined.
For the mean (over the sampling period) PM2.5 concentration, there were significant associations with smokers in the home (yes: geometric mean (GM) = 301.35 μg/m3; no: GM = 61.65 μg/m3; p = 0.0005), number of smokers in the home (exponentiated regression coefficient = 1.74; p = 0.03), and fuel type and location (p < 0.0001) (Table 3). Pairwise comparisons found significant differences between gas, indoor (GM = 79.32 μg/m3), and wood, indoor (GM = 541.14 μg/m3; p = 0.0002), and between wood, indoor, and outdoor (GM = 48.38 μg/m3; p = 0.0006) but not between gas, indoor, and outdoor (p = 0.56). For mean RH, there was a significant association with fuel type and location (p < 0.0001). Pairwise comparisons found significant differences between gas, indoor (mean = 54.40%), and outdoor (mean = 67.97%; p < 0.0001) and between wood, indoor (mean = 58.51%), and outdoor (p = 0.006) but not between gas, indoor, and wood, indoor (p = 0.15). Smokers in the home (yes: mean = 25.63 °C; no: mean = 24.58 °C; p = 0.05), number of smokers in the home (regression coefficient = 0.55 °C; p = 0.04), and fuel type and location (p < 0.0001) were significantly associated with mean temperature. Pairwise comparisons found significant differences between gas, indoor (mean = 24.67 °C) and wood, indoor (mean = 26.20 °C; p = 0.02), between gas, indoor, and outdoor (mean = 22.11 °C; p = 0.001), and between wood, indoor, and outdoor (p < 0.0001). No other home characteristic was significantly associated with mean PM2.5 concentration, RH, or temperature.
Box plots for mean PM2.5 concentration, RH, and temperature according to fuel type and location showed patterns that were similar to those just described (e.g., the distribution of mean PM2.5 concentration was higher for wood, indoor than for gas, indoor, and outdoor) (Figure 3).
Including smokers in the home and fuel type and location as independent variables in the multiple linear regression model for mean PM2.5 concentration did not change the significance of associations (smokers in the home: p = 0.04; fuel type and location: p = 0.03). However, fuel type and location (p = 0.008) remained significantly associated with mean PM2.5 concentration, but number of smokers in the home did not (p = 0.87), when number of smokers in the home and fuel type and location were included as independent variables in the multiple linear regression model for mean PM2.5 concentration. Fuel type and location (p = 0.03) remained significantly associated with mean temperature, but smokers in the home did not (p = 0.63), when smokers in the home and fuel type and location were included as independent variables in the multiple linear regression model for mean temperature. Similarly, fuel type and location (p = 0.04) remained significantly associated with mean temperature, but number of smokers in the home did not (p = 0.69), when number of smokers in the home and fuel type and location were included as independent variables in the multiple linear regression model for mean temperature.
The line graph of mean (over fuel type and location categories) PM2.5 concentration over time showed increases in mean PM2.5 concentrations primarily at mealtimes, particularly in homes that burned wood indoors (Figure 4). Mean RH decreased until mid-afternoon, then increased until the morning, and then decreased again, particularly outdoors. The pattern for mean temperature was opposite to the pattern for mean RH.

4. Discussion

Indoor cooking and heating with biomass fuels on open fires, as opposed to cleaner biomass burning devices such as chimney cookstoves and gasifiers, is recognized as a major source of PM2.5 exposure globally, particularly among populations living below the poverty level [25]. Prior studies in Nepal report high levels of household air pollutants due to biomass fuel use among various populations [29,30,31]. However, there is little data specific to brick kiln workers, who represent a uniquely vulnerable and relatively large workforce in the Kathmandu valley. Brick kiln workers are at increased risk for respiratory diseases due to crowded housing [24], low socio-economic status, and occupational exposures associated with brick manufacturing [19,32]. Adding to these risks, the results presented in this study suggest that brick kiln workers and their families are potentially exposed to dangerously high levels of PM2.5 during non-working hours. Specifically, GM PM2.5 levels in homes using wood and LPG as the primary cooking fuels were 21.6 and 3.2 times the WHO’s recommended 24 h limit of 25.0 µg/m3 [25], respectively. Considering the high prevalence of respiratory symptoms in this population, household air pollution exposures should be considered as a potentially significant environmental risk factor, warranting additional research on the topic.
Our findings indicate that PM2.5 concentration was significantly associated with device-fuel type (open wood fire vs. LPG cookstove). This finding differs from PM2.5 levels reported by Thygerson et al. (2019) showing that differences by fuel type were not significant [24]. Of note, however, air samples collected by Thygerson et al. were collected during the middle of the workday and did not capture exposures during morning and evening mealtimes. Thus, their results may mostly reflect the infiltration of outdoor air pollution into homes. One advantage of this study is that we collected samples across an approximately 24 h period, capturing mealtime cooking. We also simultaneously sampled PM2.5 using optical and gravimetric methods, allowing us to identify trends in concentration over time. Using this sampling strategy, we were able to show that increases in mean PM2.5 concentration largely appear to correspond with mealtimes, particularly for wood fuel homes.
Although PM2.5 levels in LPG homes in this study were still significantly higher than the WHO’s guidelines, they were much lower than levels in wood fuel homes. A sizeable proportion of the indoor PM2.5 pollution in LPG homes likely originates from ambient air pollution infiltrating into homes. Study homes were mostly made of un-mortared brick walls, with gaps between the bricks, the roof and walls, and at the door opening. Thygerson et al. (2019) compared PM2.5 pollution in brick workers’ homes at four kilns in Bhaktapur, Nepal, and found that the overall PM2.5 concentration and the elemental and carbon components did not differ between indoor and outdoor samples. Indeed, the indoor/outdoor (I/O) ratio for PM2.5 was 0.98, suggesting that gaps in building materials for these homes allow for significant infiltration of outdoor air. Most brick-making tasks occur outside, where workers are exposed to daily ambient air pollution while simultaneously being exposed to occupational dust and silica [20,32]. Based on our results, brick kiln workers do not appear to have a recovery period from hazardous inhalation exposure, as indoor PM2.5 levels remained above the WHO’s guidelines overnight, even in LPG homes. In addition to providing LPG fuel and stoves for brick workers (as opposed to wood or other biofuels), efforts to decrease PM2.5 exposure in this population must also focus on lowering ambient air pollution levels in the Kathmandu valley, which have increased significantly in recent years [30,33,34,35].
We are not aware of any formal assessment of household fuel use among brick workers in Nepal. Anecdotally, it appeared that more homes used LPG fuel than wood for cooking at the kiln included in this study. However, approximately 60% of the homes in the study by Thygerson et al. (2019) used wood fuel [24], suggesting that fuel use patterns may vary by kiln. The WHO recommends that, prior to policy planning, household energy use (fuels used for cooking, heating, and lighting) should be assessed for specific populations [25]. Based on findings from this study and that of Thygerson et al., we agree that a fuel use assessment in brick workers’ housing is needed to inform future intervention measures. LPG gas stoves may already be adopted by a relatively large proportion of brick workers and may provide a cleaner alternative to wood fuel. However, based on the WHO’s guidelines, a formal evaluation is needed to determine the cultural acceptability, sustainability, cost, supply, and potential exposures and safety to home occupants of transitioning to LPG cooking fuel [25].
In the absence of large fossil fuel sources, Nepal is obliged to import LPG fuel primarily from India [36]. The demand for LPG fuel, particularly in urban areas such as the Kathmandu valley, is expected to increase by 3.9%–4.9% per year through 2040 [37], even while the cost of LPG has increased by 8% annually [38]. Currently, the government of Nepal subsidizes the cost of LPG, but there is a movement to discontinue these subsidies and to divert that money to the development of hydroelectric power [36]. Thus, efforts to transition brick workers from open wood fires to LPG cookstoves will require monitoring and reporting to identify ongoing sociopolitical impacts, such as increasing LPG fuel cost, that may create barriers to use. Given the vast, undeveloped hydroelectric potential in Nepal [36], low-cost electric cookstoves may provide a long-term solution to decrease indoor air pollution exposure among brick workers. However, there may also be a significant time lag for the development of hydroelectric power generation and distribution of sufficient electricity throughout the Kathmandu valley. This may result in ongoing, long-term exposure to household air pollution among brick workers that could, in the interim, be lowered by the temporary use of LPG cookstoves. In addition, cooking stacking (using multiple fuel/device options) may be an obstacle to cleaner energy use in some populations [39]. Providing brick workers with more than one clean fuel-device option may help to prevent reversion to pollution prone options such as open wood fires. A clearer understanding of fuel use patterns, and barriers to adopting LPG fuel, is an essential first step to guide interventions to decrease household air pollution exposure among brick workers in Nepal, with electric cookstoves or clean cooking stacking options as potential long-term solutions.
A previous study conducted on the prevalence of respiratory illnesses among brick workers in Nepal shows that smoking was a significant predictor (p < 0.001) of chronic bronchitis and chronic wheezing [19]. Of the world’s 1.1 billion global smokers, around 80% live in low- and middle-income countries [40]. In Nepal, 27.4% of adult males and 5.5% of adult females smoke tobacco, and the prevalence of smoking is higher among those with lower levels of education [41]. In our study, approximately half of the homes sampled had at least one smoker living in the home. This is similar to the findings of Thygerson et al. (2019). In their study, 67% of homes had a least one smoker, many of which reported smoking indoors [24]. We also found that PM2.5 levels were significantly higher in homes with smokers compared to homes with no smokers. Thus, we recognize that multi-faceted interventions that include smoking cessation are necessary in order to decrease respiratory illness among brick workers in Nepal. The primary national intervention strategies to limit tobacco use in Nepal are prescribed in the Tobacco Product (Control and Regulatory) Act, which was signed into law in 2011. This law requires the size of warning labels on tobacco products to cover at least 90% of the package area, prohibits the distribution and advertisement of tobacco products to children, levies a tax on tobacco products, and bans smoking in public areas, among other provisions [42,43]. Worksite smoking cessation interventions have shown some success in developed countries [44,45]. In addition to the national intervention measures, public health agencies and/or non-governmental organizations (NGOs) could work directly with individual brick kilns to implement worksite smoking cessation programs that include financial incentives and other evidence-based strategies [44,45].
Our findings regarding temperature and RH may be cause for further concern regarding brick workers’ respiratory and overall health. The average RH in brick workers’ homes was approximately 58%. As evenings dropped to cooler temperatures, RH in both gas and wood burning homes increased to levels between 55% and 75% and remained within this range for several hours, often until almost noon on the next day. This is a concern because house dust mites (HDMs) grow at RH levels between 55% and 75% [46,47]. HDM allergens are among the most clinically significant antigen exposures affecting humans [48,49,50] and exposure to them is associated with the development of asthma and other allergic diseases [48,49,50]. In light of the already present risks of ambient and indoor air pollution, the possibility of HDMs within brick workers’ bedding presents an added health concern. We recommend further study investigating the presence of HDMs as a possible contributing factor to the high levels of respiratory problems among brick kiln workers in Nepal [19,32].
To our knowledge, this is the first study to characterize indoor PM2.5 levels in brick workers’ homes in Nepal for an approximately 24 h period. One strength of our study design is that we simultaneously collected filter-based (gravimetric) and optical particle counts, allowing us to identify PM2.5 trends over time. This was advantageous in that it allowed us not only to identify mean differences between wood and LPG fuel homes but also to see the influence that cooking has on indoor pollution levels, which was not studied in previous research by Thygerson et al., 2019 [24].

Limitations

This study was limited to a relatively small number of homes at a single brick kiln in Bhaktapur, Nepal. Thus, our study may have been underpowered to detect some associations (i.e., type II error) between brick kiln home characteristics and the three outcomes and our point estimates had wide 95% CI. It is also possible that our small sample size could have led to spurious results (i.e., type I error). In addition, we used a convenience sample of homes and our findings may not be generalizable to other kilns in the Kathmandu valley or to other clay brick kilns worldwide. In addition, we were limited to sampling brick workers’ homes over one 24 h period during a single season, which also may have contributed to lower statistical power and wider 95% CI than if we had been able to collect multiple 24 h samples from each home. Longer-term sampling or repeated sampling of the same homes across seasons will likely provide a much better understanding of PM2.5, temperature, and RH levels in brick workers’ homes. Unmeasured confounding by temporal factors could have affected associations between brick kiln home characteristics and the three outcomes. We did not collect information regarding technologies (e.g., open fires vs. cookstoves that had venting/chimneys) and, therefore, could not estimate associations between use of technologies and PM2.5, RH, and temperature levels or the combined effect of use of fuels and technologies on these outcomes. Finally, this study was limited to environmental measurements. Considering the relatively high levels of indoor PM2.5 in both LPG and wood fuel homes, studies of respiratory and other health effects associated with brick workers’ housing are needed.

5. Conclusions

In this study, indoor PM2.5 levels for both wood and LPG fuel homes exceeded the WHO’s recommended 24 h limit of 25.0 µg/m3, but the exceedance was most pronounced for wood fuel homes. Peaks in PM2.5 appeared to coincide largely with mealtimes in wood fuel homes. For LPG homes, mealtime peaks were significantly less pronounced, and overall levels may be significantly influenced by the infiltration of ambient air pollution. Based on our findings and those of a previous study [24], LPG cookstoves appear to be used by a relatively large proportion of brick workers, suggesting that they may be a culturally appropriate alternative to wood fuel cooking. Future studies should focus on understanding barriers to LPG stove use, the potential for the adoption of clean cooking stacking, as well as interventions to decrease the use of wood fuel for cooking in this population of workers.

Author Contributions

Conceptualization, J.D.J., J.D.B.; Methodology, J.D.J., J.D.B.; Investigation, J.D.J., M.E.H., H.B.J., L.A.J.; Formal analysis, J.D.B., M.E.H.; Project administration and supervision, J.D.J., J.D.B.; Resources, J.D.J., J.D.B.; Validation, J.D.J., J.D.B., M.E.H.; Visualization, J.D.B., M.E.H., J.D.J.; Writing—original draft, J.D.J., J.D.B., M.E.H., H.B.J., L.A.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Our special thanks to the brick kiln manager and staff who supported this study, to Sirish, our interpreter, and to the brick workers and their families for allowing us into their homes for sample collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization: Household Air Pollution and Health. Fact Sheet. Available online: http://www.who.int/en/news-room/fact-sheets/detail/household-air-pollution-and-health (accessed on 10 June 2020).
  2. Ezzati, M.; Kammen, D.M. Household Energy, Indoor Air Pollution and Health in Devleoping Countries: Knowledge Base for Effective Interventions. Annu. Rev. Energy Environ. 2002, 27, 233–270. [Google Scholar] [CrossRef] [Green Version]
  3. Zhang, J.; Smith, K.R. Indoor Air Pollution: A Global Health Concern. Br. Med. Bull. 2003, 68, 209–225. [Google Scholar] [CrossRef] [Green Version]
  4. Naeher, L.P.; Brauer, M.; Lipsett, M.; Zelikoff, J.T.; Simpson, C.D.; Koenig, J.Q.; Smith, K.R. Woodsmoke Health Effects: A Review. Inhal. Toxicol. 2007, 19, 67–106. [Google Scholar] [CrossRef]
  5. Fullerton, D.G.; Bruce, N.; Gordon, S.B. Indoor Air Pollution from Biomass Fuel Smoke is a Major Health Concern in the Developing World. Trans. R. Soc. Trop. Med. Hyg. 2008, 102, 843–851. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zhang, J.J.; Smith, K.R. Household Air Pollution from Coal and Biomass Fuels in China: Measurements, Health Impacts, and Interventions. Environ. Health Perspect. 2007, 115, 848–855. [Google Scholar] [CrossRef]
  7. Kim, K.H.; Jahan, S.A.; Kabir, E. A Review of Diseases Associated with Household Air Pollution Due to the Use of Biomass Fuels. J. Hazard. Mater. 2011, 192, 425–431. [Google Scholar] [CrossRef]
  8. Gill, J. Improved Stoves in Developing Countries—A Critique. Energy Policy 1987, 15, 135–144. [Google Scholar] [CrossRef]
  9. Barnes, D.F. What Makes People Cook with Improved Biomass Stoves? A Comparative International Review of Stove Programs; World Bank: Washington, DC, USA, 1994; pp. 13–19. [Google Scholar]
  10. Clark, M.L.; Peel, J.L.; Balakrishnan, K.; Breysse, P.N.; Chillrud, S.N.; Naeher, L.P.; Rodes, C.E.; Vette, A.F.; Balbus, J.M. Health and Household Air Pollution from Solid Fuel Use: The Need for Improved Exposure Assessment. Environ. Health Perspect. 2013, 121, 1120–1128. [Google Scholar] [CrossRef]
  11. Baldwin, S.; Geller, H.; Dutt, G.; Ravindranath, N.H. Improved Woodbunning Cookstoves: Signs of Success. Ambio 1985, 14, 280–287. [Google Scholar]
  12. Gyawali, K.; Sharma, S.; Sharma, R.K. A Rapid Assessment of Children in the Brick Industry; World Education: Kathmandu, Nepal, 2012; Available online: http://3vvxza37b4sa2ci1ug2o2qwg-wpengine.netdna-ssl.com/wp-content/uploads/2014/02/nbnp-ra-brick-industry.pdf (accessed on 12 June 2020).
  13. Larmar, S.; O’Leary, P.; Chui, C.; Benfer, K.; Zug, S.; Jordan, L.P. Hazardous Child Labor in Nepal: The Case of Brick Kilns. Child Abus. Negl. 2017, 72, 312–325. [Google Scholar] [CrossRef]
  14. MinErgy. A Roadmap for Efficient Brick Production in Nepal. 2012. Available online: http://www.repic.ch/files/6113/7544/1756/SB_AGHH_Minergy_Road_Map.pdf (accessed on 12 June 2020).
  15. GEFONT. Nepal: Labour under the Chimney. A Study on the Brick Kilns of Nepal; GEFONT: Kathmandu, Nepal, 2007; p. 14. Available online: https://www.gefont.org/assets/upload/downloads/brick_klin_report_2008.pdf (accessed on 12 June 2020).
  16. Haack, B.N.; Khatiwada, G. Rice and Bricks: Environmental Issues and Mapping of the Unusual Crop Rotation Pattern in the Kathmandu Valley, Nepal. Environ. Manag. 2007, 39, 774–782. [Google Scholar] [CrossRef] [PubMed]
  17. ENPHO. A Study on Status of Brick Industry in the Kathmandu Valley; Environmental and Public Health Organization: Kathmandu, Nepal, 2001. [Google Scholar]
  18. Raut, A.K. Brick Kilns in Kathmandu Valley: Current Status, Environmental Impacts and Future Options. Himal. J. Sci. 2003, 1, 59–61. [Google Scholar] [CrossRef] [Green Version]
  19. Sanjel, S.; Khanal, S.N.; Thygerson, S.M.; Carter, W.S.; Johnston, J.D.; Joshi, S.K. Respiratory Symptoms and Illnesses Related to the Concentration of Airborne Particulate MATTER among Brick Kiln Workers in Kathmandu Valley, Nepal. Ann. Occup. Environ. Med. 2017, 29. [Google Scholar] [CrossRef] [Green Version]
  20. Sanjel, S.; Khanal, S.N.; Thygerson, S.M.; Carter, W.; Johnston, J.D.; Joshi, S.K. Exposure to Respirable Silica among Clay Brick Workers in Kathmandu Valley, Nepal. Arch. Environ. Occup. Health 2017, 73, 347–350. [Google Scholar] [CrossRef]
  21. Myers, J.E.; Lewis, P.; Hofmeyr, W. Respiratory health of brickworkers in Cape Town, South Africa. Background, Aims and Dust Exposure Determinations. Scand. J. Work Environ. Health 1989, 15, 180–187. [Google Scholar] [CrossRef] [PubMed]
  22. Wiecek, E.; Gościcki, J.; Indulski, J.; Stroszejn-Mrowca, G. Dust and Occupational Diseases in Brickyards. Med. Pr. 1983, 34, 35–45. [Google Scholar] [PubMed]
  23. Love, R.G.; Waclawski, E.R.; Maclaren, W.M.; Wetherill, G.Z.; Groat, S.K.; Porteous, R.H.; Soutar, C.A. Risks of Respiratory Disease in the Heavy Clay Industry. Occup. Environ. Med. 1999, 56, 124–133. [Google Scholar] [CrossRef] [Green Version]
  24. Thygerson, S.M.; Beard, J.D.; House, M.J.; Smith, R.L.; Burbidge, H.C.; Andrus, K.N.; Weber, F.X.; Chartier, R.; Johnston, J.D. Air-Quality Assessment of On-Site Brick-Kiln Worker Housing in Bhaktapur, Nepal: Chemical Speciation of Indoor and Outdoor PM 2.5 Pollution. Int. J. Environ. Res. Public Health. 2019, 16, 4114. [Google Scholar] [CrossRef] [Green Version]
  25. World Health Organization. WHO Indoor Air Quality Guidelines: Household Fuel Combustion; World Health Organization: Geneva, Switzerland, 2014; Available online: https://www.who.int/airpollution/guidelines/household-fuel-combustion/IAQ_HHFC_guidelines.pdf (accessed on 12 June 2020).
  26. Code of Federal Regulations Title 45, Part 46.102. Protection of Human Subjects; Definitions for Purposes of this Policy. Available online: https://www.ecfr.gov/cgi-bin/retrieveECFR?gp=&SID=83cd09e1c0f5c6937cd9d7513160fc3f&pitd=20180719&n=pt45.1.46&r=PART&ty=HTML#se45.1.46_1102 (accessed on 12 June 2020).
  27. Waldron, P.F. Principles and instrumentation for calibrating air sampling equipment. In The Occupational Environment–Its Evaluation and Control; DiNardi, R., Ed.; AIHA Press: Fairfax, VA, USA, 1997; pp. 155–175. [Google Scholar]
  28. Rodes, C.E. Standard Operating Procedure MicroPEM v 3.2A [Gates Foundation Version] Aerosol Exposure SOP M-300: Setup, Preparation, Calibration, Deployment, Data Processing and Servicing; RTI International: Research Triangle Park, NC, USA, 2013; pp. 1–46. [Google Scholar]
  29. Singh, A.; Tuladhar, B.; Bajracharya, K.; Pillarisetti, A. Assessment of Effectiveness of Improved Cook Stoves in Reducing Indoor Air Pollution and Improving Health in Nepal. Energy Sustain. Dev. 2012, 16, 406–414. [Google Scholar] [CrossRef]
  30. Gurung, A.; Bell, M.L. The State of Scientific Evidence on Air Pollution and Human Health in Nepal. Environ. Res. 2013, 124, 54–64. [Google Scholar] [CrossRef]
  31. Shrestha, I.L.; Shrestha, S.L. Indoor Air Pollution from Biomass Fuels and Respiratory Health of the Exposed Population in Nepalese Households. Int. J. Occup. Environ. Health 2005, 11, 150–160. [Google Scholar] [CrossRef] [PubMed]
  32. Sanjel, S.; Thygerson, S.M.; Khanal, S.N.; Joshi, S.K. Environmental and Occupational Pollutants and their Effects on Health among Brick Kiln Workers. Open J. Saf. Sci. Technol. 2016, 6, 81–98. [Google Scholar] [CrossRef] [Green Version]
  33. Saud, B.; Paudel, G. The Threat of Ambient Air Pollution in Kathmandu, Nepal. J. Environ. Public Health 2018, 2018. [Google Scholar] [CrossRef] [Green Version]
  34. Shakya, K.M.; Rupakheti, M.; Shahi, A.; Maskey, R.; Pradhan, B.; Panday, A.; Puppala, S.P.; Lawrence, M.; Peltier, R.E. Near-Road Sampling of PM2.5, BC, and Fine-Particle Chemical Components in Kathmandu Valley, Nepal. Atmos. Chem. Phys. 2017, 17. [Google Scholar] [CrossRef] [Green Version]
  35. Aryal, R.K.; Lee, B.K.; Karki, R.; Gurung, A.; Baral, B.; Byeon, S.H. Dynamics of PM2.5 Concentrations in Kathmandu Valley, Nepal. J. Hazard. Mater. 2009, 168. [Google Scholar] [CrossRef]
  36. Bhandari, R.; Pandit, S. Electricity as a Cooking Means in Nepal—A Modelling Tool Approach. Sustainability 2018, 10, 2841. [Google Scholar] [CrossRef] [Green Version]
  37. Malla, S. Household Energy Consumption Patterns and its Environmental Implications: Assessment of Energy Access and Poverty in Nepal. Energy Policy 2013, 61, 990–1002. [Google Scholar] [CrossRef]
  38. Parajuli, R.; Østergaard, P.A.; Dalgaard, T.; Pokharel, G.R. Energy Consumption Projection of Nepal: An Economic Approach. Renew. Energy 2014, 63, 432–444. [Google Scholar] [CrossRef]
  39. Medina, P.; Berrueta, V.; Cinco, L.; Ruiz-Garcia, V.; Edwards, R.; Olaya, B.; Schilmann, A.; Masera, O. Understanding Household Energy Transitions: From Evaluating singLe Cookstoves to “Clean Stacking” Alternatives. Atmosphere 2019, 10, 693. [Google Scholar] [CrossRef] [Green Version]
  40. Perez-Warnisher, M.T.; de Miguel, M.D.P.C.; Seijo, L.M. Tobacco Use Worldwide: Legislative Efforts to Curb Consumption. Ann. Glob. Health 2019, 85, 9. [Google Scholar] [CrossRef] [Green Version]
  41. Shrestha, N.; Mehata, S.; Pradhan, P.M.S.; Joshi, D.; Mishra, S.R. A nationally Representative Study on Socio-Demographic and Geographic Correlates, and Trends in Tobacco Use in Nepal. Sci. Rep. 2019, 9, 2682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Bhatta, D.N.; Crosbie, E.; Bialous, S.A.; Glantz, S. Defending Comprehensive Tobacco Control Policy Implementation in Nepal from Tobacco Industry Interference (2011–2018). Nicotine Tob. Res. 2020. [Google Scholar] [CrossRef] [PubMed]
  43. Chang, C.P.; Siwakoti, B.; Sapkota, A.; Gautam, D.K.; Lee, Y.C.A.; Monroe, M.; Hashibe, M. Tobacco Smoking, Chewing Habits, Alcohol Drinking and the Risk of Head and Neck Cancer in Nepal. Int. J. Cancer 2019, 147, 866–875. [Google Scholar] [CrossRef]
  44. van den Brand, F.A.; Nagelhout, G.E.; Winkens, B.; Chavannes, N.H.; van Schayck, O.C. Effect of a Workplace-Based Group Training Programme Combined with Financial Incentives on Smoking Cessation: A Cluster-Randomised Controlled Trial. Lancet Public Health 2018, 3, e536–e544. [Google Scholar] [CrossRef] [Green Version]
  45. Cahill, K.; Lancaster, T. Workplace Interventions for Smoking Cessation. Cochrane Database Syst. Rev. 2014, 2. [Google Scholar] [CrossRef]
  46. Arlian, L.G.; Veselic, M.M. Reevaluation of the Humidity Requirements of the House Dust Mite Dermatophagoides farinae (Acari: Pyroglyphidae). J. Med. Entomol. 1981, 18, 351–352. [Google Scholar] [CrossRef]
  47. Arlian, L.G.; Platts-Mills, T.A.E. The Biology of Dust Mites and the Remediation of Mite Allergens in Allergic Disease. J. Allergy Clin. Immunol. 2001, 107, S406–S413. [Google Scholar] [CrossRef] [PubMed]
  48. Sporik, R.; Holgate, S.T.; Plattsmills, T.A.E.; Cogswell, J.J. Exposure to House-Dust Mite Allergen (Der-P-I) and the Development of Asthma in Childhood—A Prospective Study. N. Engl. J. Med. 1990, 323, 502–507. [Google Scholar] [CrossRef]
  49. Portnoy, J.; Miller, J.D.; Williams, P.B.; Chew, G.L.; Miller, J.D.; Zaitoun, F.; Phipatanakul, W.; Kennedy, K.; Barnes, C.; Grimes, C.; et al. Environmental Assessment and Exposure Control of Dust Mites: A Practice Parameter. Ann. Allergy Asthma Immunol. 2013, 111, 465–507. [Google Scholar] [CrossRef] [Green Version]
  50. Lodge, C.J.; Lowe, A.J.; Gurrin, L.C.; Hill, D.J.; Hosking, C.S.; Khalafzai, R.U.; Hopper, J.L.; Matheson, M.C.; Abramson, M.J.; Allen, K.J.; et al. House Dust Mite Sensitization in Toddlers Predicts Current Wheeze at Age 12 Years. J. Allergy Clin. Immunol. 2011, 128, 782–788. [Google Scholar] [CrossRef]
Figure 1. Typical construction of brick workers’ homes sampled in this study. Construction consisted of dirt floors and un-mortared brick walls. Tin roofs were constructed by placing bamboo or wooden poles across the top of the walls, with sheet metal placed on the poles. Rocks, bricks, soil, and other heavy objects were used to hold the sheet metal in place. Indoor samples were collected by attaching PM2.5 and temperature/relative humidity instruments to tripods as shown. Tripods with the attached instruments were placed inside the home during the entire sampling period. Abbreviation: PM2.5, particulate matter with an aerodynamic diameter less than 2.5 μm.
Figure 1. Typical construction of brick workers’ homes sampled in this study. Construction consisted of dirt floors and un-mortared brick walls. Tin roofs were constructed by placing bamboo or wooden poles across the top of the walls, with sheet metal placed on the poles. Rocks, bricks, soil, and other heavy objects were used to hold the sheet metal in place. Indoor samples were collected by attaching PM2.5 and temperature/relative humidity instruments to tripods as shown. Tripods with the attached instruments were placed inside the home during the entire sampling period. Abbreviation: PM2.5, particulate matter with an aerodynamic diameter less than 2.5 μm.
Ijerph 17 05681 g001
Figure 2. Illustrative examples of the two cooking practices observed in brick workers’ homes in Bhaktapur, Nepal: (a) indoor wood fire and (b) LPG (liquefied petroleum gas) cook stove.
Figure 2. Illustrative examples of the two cooking practices observed in brick workers’ homes in Bhaktapur, Nepal: (a) indoor wood fire and (b) LPG (liquefied petroleum gas) cook stove.
Ijerph 17 05681 g002
Figure 3. Box plots of the mean of samples inside or outside a homes by fuel type and location at a brick kiln in Bhaktapur, Nepal, May 2019: (a) PM2.5, (b) relative humidity, (c) temperature. Reading from the top to the bottom, the horizontal lines on each box plot represent the maximum, third quartile, median, first quartile, and minimum, and the black triangle represents the mean. Abbreviation: PM2.5, particulate matter with an aerodynamic diameter less than 2.5 μm. a Nineteen samples were from inside and four samples were from outside the homes (for PM2.5, 17 samples were from inside and three samples were from outside the homes).
Figure 3. Box plots of the mean of samples inside or outside a homes by fuel type and location at a brick kiln in Bhaktapur, Nepal, May 2019: (a) PM2.5, (b) relative humidity, (c) temperature. Reading from the top to the bottom, the horizontal lines on each box plot represent the maximum, third quartile, median, first quartile, and minimum, and the black triangle represents the mean. Abbreviation: PM2.5, particulate matter with an aerodynamic diameter less than 2.5 μm. a Nineteen samples were from inside and four samples were from outside the homes (for PM2.5, 17 samples were from inside and three samples were from outside the homes).
Ijerph 17 05681 g003aIjerph 17 05681 g003b
Figure 4. Line graphs of the mean of samples inside or outside a homes over time by fuel type and location at a brick kiln in Bhaktapur, Nepal, May 2019: (a) PM2.5, (b) subset of PM2.5 (i.e., from 0 to 2000 μg/m3), (c) relative humidity, (d) temperature. Abbreviation: PM2.5, particulate matter with an aerodynamic diameter less than 2.5 μm. a Nineteen samples were from inside and four samples were from outside the homes (for PM2.5, 17 samples were from inside and three samples were from outside the homes).
Figure 4. Line graphs of the mean of samples inside or outside a homes over time by fuel type and location at a brick kiln in Bhaktapur, Nepal, May 2019: (a) PM2.5, (b) subset of PM2.5 (i.e., from 0 to 2000 μg/m3), (c) relative humidity, (d) temperature. Abbreviation: PM2.5, particulate matter with an aerodynamic diameter less than 2.5 μm. a Nineteen samples were from inside and four samples were from outside the homes (for PM2.5, 17 samples were from inside and three samples were from outside the homes).
Ijerph 17 05681 g004aIjerph 17 05681 g004b
Table 1. Characteristics of homes at a brick kiln in Bhaktapur, Nepal, May 2019.
Table 1. Characteristics of homes at a brick kiln in Bhaktapur, Nepal, May 2019.
CharacteristicHomes, n (%)Missing, nMeanSDMinQ1MedianQ3Max
Total19 (100)
Home area, m2 11.305.635.418.0310.2511.6031.40
Number of people in home 3.421.431.002.004.004.007.00
Occupant density, number of people/100 m2 33.2515.699.5524.9129.7039.0373.96
Number of children in home
 0-110 (53)
 2-39 (47)
Smokers in home 2
 No8 (47)
 Yes9 (53)
Number of smokers in home 20.761.030.000.001.001.004.00
Fuel type 1
 Gas12 (67)
 Wood6 (33)
Abbreviations: Max, maximum; Min, minimum; Q1, first quartile; Q3, third quartile; SD, standard deviation.
Table 2. Summary statistics for the mean of samples inside or outside a homes at a brick kiln in Bhaktapur, Nepal, May 2019.
Table 2. Summary statistics for the mean of samples inside or outside a homes at a brick kiln in Bhaktapur, Nepal, May 2019.
VariableSamples, nMissing, nMeanSDMinQ1MedianQ3Max
PM2.5, μg/m3203272.13389.5419.3758.73118.46278.491384.44
Relative humidity, %23 57.876.3746.5252.6658.3960.4773.70
Temperature, °C23 24.581.6620.6623.7224.8825.8926.87
Abbreviations: Max, maximum; Min, minimum; PM2.5, particulate matter with an aerodynamic diameter less than 2.5 μm; Q1, first quartile; Q3, third quartile; SD, standard deviation. a Nineteen samples were from inside and four samples were from outside the homes (for PM2.5, 17 samples were from inside and three samples were from outside the homes).
Table 3. Associations between home characteristics and the mean of samples inside or outsidea homes at a brick kiln in Bhaktapur, Nepal, May 2019.
Table 3. Associations between home characteristics and the mean of samples inside or outsidea homes at a brick kiln in Bhaktapur, Nepal, May 2019.
CharacteristicPM2.5, μg/m3Relative Humidity, %Temperature, °C
GM b95% CI bp-Value bMean c95% CI cp-Value cMean c95% CI cp-Value c
Home area, m20.97 d0.87, 1.08 d0.540.01 e−0.37, 0.40 e0.95−0.03 e−0.13, 0.07 e0.58
Number of people in home0.88 d0.58, 1.35 d0.54−1.04 e−2.47, 0.39 e0.140.01 e−0.39, 0.41 e0.96
Occupant density, number of people/100 m21.01 d0.97, 1.05 d0.48−0.04 e−0.18, 0.10 e0.540.02 e−0.02, 0.05 e0.30
Number of children in home
 0–1160.5770.82, 364.06 57.1354.40, 59.87 25.2624.50, 26.02
 2–3150.1856.46, 399.500.9154.2151.32, 57.100.1424.9324.13, 25.730.54
Smokers in home
 No61.6536.05, 105.43 54.1650.86, 57.45 24.5823.82, 25.34
 Yes301.35176.20, 515.370.000556.9153.80, 60.020.2125.6324.91, 26.350.05
Number of smokers in home1.74 d1.08, 2.81d0.031.19 e−1.10, 3.48 e0.290.55 e0.04, 1.06 e0.04
Fuel type and location
 Gas, indoor79.3249.82, 126.28 54.4051.88, 56.92 24.6724.04, 25.29
 Wood, indoor541.14288.31, 1015.68 58.5154.94, 62.07 26.2025.32, 27.08
 Outdoor48.3819.86, 117.87<0.0001 f67.9763.60, 72.33<0.0001 g22.1121.04, 23.19<0.0001 h
Abbreviations: CI, confidence interval; GM, geometric mean; PM2.5, particulate matter with an aerodynamic diameter less than 2.5 μm. a Nineteen samples were from inside and four samples were from outside the homes (for PM2.5, 17 samples were from inside and three samples were from outside the homes). Except for analyses involving fuel type and location, only indoor samples were used for analyses shown in this table. b Estimated via linear regression models of the natural logarithm transformed values. c Estimated via linear regression models of the original values. d Exponentiated regression coefficient and 95% CI (i.e., GM PM2.5 concentration ratio for a specified change in the independent variable or exp(β)−1 = percent change in GM PM2.5 concentration for a specified change in the independent variable). e Regression coefficient and 95% CI (i.e., change in relative humidity or temperature for a specified change in the independent variable). f Using the Tukey method to adjust for multiple comparisons, p-values for tests of pairwise differences among fuel type and location categories were as follows: gas, indoor vs. wood, indoor: 0.0002; gas, indoor vs. outdoor: 0.56; wood, indoor vs. outdoor: 0.0006. g Using the Tukey method to adjust for multiple comparisons, p-values for tests of pairwise differences among fuel type and location categories were as follows: gas, indoor vs. wood, indoor: 0.15; gas, indoor vs. outdoor: <0.0001; wood, indoor vs. outdoor: 0.006. h Using the Tukey method to adjust for multiple comparisons, p-values for tests of pairwise differences among fuel type and location categories were as follows: gas, indoor vs. wood, indoor: 0.02; gas, indoor vs. outdoor: 0.001; wood, indoor vs. outdoor: <0.0001.

Share and Cite

MDPI and ACS Style

Johnston, J.D.; Hawks, M.E.; Johnston, H.B.; Johnson, L.A.; Beard, J.D. Comparison of Liquefied Petroleum Gas Cookstoves and Wood Cooking Fires on PM2.5 Trends in Brick Workers’ Homes in Nepal. Int. J. Environ. Res. Public Health 2020, 17, 5681. https://doi.org/10.3390/ijerph17165681

AMA Style

Johnston JD, Hawks ME, Johnston HB, Johnson LA, Beard JD. Comparison of Liquefied Petroleum Gas Cookstoves and Wood Cooking Fires on PM2.5 Trends in Brick Workers’ Homes in Nepal. International Journal of Environmental Research and Public Health. 2020; 17(16):5681. https://doi.org/10.3390/ijerph17165681

Chicago/Turabian Style

Johnston, James D., Megan E. Hawks, Haley B. Johnston, Laurel A. Johnson, and John D. Beard. 2020. "Comparison of Liquefied Petroleum Gas Cookstoves and Wood Cooking Fires on PM2.5 Trends in Brick Workers’ Homes in Nepal" International Journal of Environmental Research and Public Health 17, no. 16: 5681. https://doi.org/10.3390/ijerph17165681

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop