Abstract
Exposure to household air pollution is a leading cause of ill-health globally. The Household Air Pollution Intervention Network (HAPIN) randomized controlled trial evaluated the impact of a free liquefied petroleum gas stove and fuel intervention on birth outcomes and maternal and child health. As part of HAPIN, an extensive exposure assessment was conducted. Here, we report on PM2.5 and CO exposures of young children (≤ 15 months old) reconstructed using a Bluetooth-beacon based time-activity monitoring system coupled with microenvironmental pollutant monitors. Median (IQR) exposures to PM2.5 were 65.1 (33 – 128.2) µg/m3 in the control group and 22.9 (17.2 – 35.3) µg/m3 in the intervention group; for CO, median (IQR) exposures were 1.1 (0.3 – 2.9) ppm and 0.2 (0 – 0.7) ppm for control and intervention group, respectively. Exposure reductions were stable over time and consistent with previous findings for the children’s mothers. In the intervention group, 75% of children’s reconstructed exposures were below the WHO interim target guideline value of 35 µg/m3, while 26% were below the standard in the control group. Our findings suggest that an LPG fuel and stove intervention can substantially reduce children’s exposure to household air pollution.
Introduction
Between 3 and 4 billion people cook with solid fuels (coal, wood, dung, or crop residues)1, 2, resulting in high exposures to household air pollution (HAP), including particulate matter with an aerodynamic diameter of 2.5 µm or less (PM2.5), carbon monoxide (CO), and black carbon (BC), among others. In 2019, 2.3 million deaths were attributed to HAP exposure, approximately 460,000 of which occurred in children under the age of 5.3
HAP-reducing interventions – including improved or advanced biomass stoves and clean fuel stoves – have had varying exposure reduction effectiveness to date.4–6 Trials evaluating the impact of interventions on health have often insufficiently measured changes in exposures7 or failed to reduce exposure to levels expected to result in improved health.8, 9 Among those that do measure exposure, many have reported little or no change in exposure, arguably due to mixed use of intervention and traditional stoves or poor-performing interventions that do not reduce exposure in real-world settings.10, 11
The multi-country Household Air Pollution Intervention Network (HAPIN) clean stove and fuel randomized controlled trial (RCT) sought to address these shortcomings. First, it provided a liquefied petroleum gas LPG stove and continuous clean fuel supply to pregnant women, with behavioral reinforcement, leading to high fidelity and adherence to the intervention during pregnancy12 and after childbirth.13 Second, we performed extensive personal exposure assessments among participants throughout the trial. Among expectant mothers, we found significantly lower exposures to HAP among the intervention group than in the control group: PM2.5 exposures were 66% lower (23.9 vs. 70.7 µg/m3), BC exposures were 71% lower (2.8 vs. 9.6 µg/m3), and CO exposures were 83% lower (0.2 vs. 1.1 ppm).14
Exposures are often measured for primary cooks, typically women; less is known about exposures to infants and young children,15 despite evidence from the ambient air pollution literature that early-in-life exposure can have respiratory16 and cardiovascular17, 18 impacts that may set health trajectories into adulthood. The scarcity of exposure data on children is due partly to challenges of wearability, compliance, and induced behavior change; alternatively, exposures can be reconstructed using microenvironmental measurements coupled with time-activity data either from sensors19, 20 or self-reported questionnaires.21 Here, we report the impact of HAPIN’s LPG intervention on estimated PM2.5 and CO exposures for HAPIN children using indirect exposure reconstruction methods validated during HAPIN’s formative phase.20
Materials and Methods
Study population
The HAPIN RCT of an LPG fuel and stove intervention occurred in Guatemala, India, Peru, and Rwanda. The approach and sites have been described in detail.22–24 Briefly, we selected rural sites in each country with low background air pollution concentrations, few additional air pollution sources, and where most households used biomass stoves for cooking.14, 25, 26 The trial assessed intervention health effects among pregnant women (n ∼ 3200), their resulting newborn children (n ∼ 3200), and non-pregnant adult women living in the same household (n ∼ 418), split evenly between intervention and control arms. Primary health outcomes are birth weight,27, 28 incidence of severe pediatric pneumonia,29 stunted growth in the children,22 and blood pressure in the non-pregnant adult women.22
The study protocol has been reviewed and approved by institutional review boards (IRBs) and Ethics Committees at Emory University (00089799), Johns Hopkins University (00007403), Sri Ramachandra Institute of Higher Education and Research (IEC-N1/16/JUL/54/49), the Indian Council of Medical Research – Health Ministry Screening Committee (5/8/4-30/(Env)/Indo-US/2016-NCD-I), Universidad del Valle de Guatemala (146-08-2016), Guatemalan Ministry of Health National Ethics Committee (11-2016), Asociación Benefica PRISMA (CE2981.17), the London School of Hygiene and Tropical Medicine (11664-5), the Rwandan National Ethics Committee (No.357/RNEC/2018), and Washington University in St. Louis (201611159). The study has been registered with ClinicalTrials.gov (Identifier NCT02944682).
Recruitment
Pregnant women were identified and enrolled through partnerships with local clinics and community health workers. Eligible women were aged 18-35 years old at 9 to <20 weeks’ gestation, with a viable singleton pregnancy (confirmed by ultrasound); used primarily biomass fuel for cooking; and agreed to participate via informed consent. Exclusion criteria included tobacco use, plans to move outside the study area, and plans to switch to primary use of clean fuels.
Intervention
Participants were randomly assigned (one-to-one) to receive an LPG stove, continuous fuel delivery, and regular behavioral messaging versus continued use of a biomass-burning stove. In India and Peru, stratified randomization was used to ensure balance between distinct geographical regions within each country. No additional stratification was used in Rwanda and Guatemala, where the study areas were deemed homogenous.
The HAPIN intervention was refined during formative research and described in detail previously.22 Stoves and fuel were procured at each IRC based on local cooking practices and equipment and fuel availability. The exact stove deployed varied by country; all stoves had at least two burners and included additional components for preparation of traditional foods (e.g., a flat griddle for cooking tortillas in Guatemala; a roasting device in Rwanda). Intervention households received the stove and a continuous supply of LPG at no cost throughout follow-up. Behavioral support included (i) a pledge by intervention homes to use the LPG stove for all cooking throughout the trial, (ii) safety training, (iii) tailored messaging to encourage exclusive use of LPG and discourage use of traditional stoves, and (iv) behavioral reinforcements upon detection of traditional stove use.
Monitoring plan and microenvironmental PM2.5 and CO measurements
We attempted collection of three 24-hour PM2.5 and CO measurements in children in the first year of life (<3 months, ∼6 months, and ∼12 months). We used a validated indirect approach to assess the exposure of children; using standard personal exposure monitors is impractical due to their size and participant preferences.22 This approach pairs microenvironmental pollutant sampling and participant proximity sensing.23 Details of the indirect exposure assessment process and approach are described below. Microenvironmental PM2.5 concentrations were measured using the Enhanced Children’s MicroPEM (ECM, RTI International, Durham, NC USA), the same device used for personal PM2.5 exposure monitoring in the HAPIN trial. The ECM measures PM2.5 concentrations every 10 seconds using a nephelometer and simultaneously collects integrated gravimetric samples on 15 mm polytetrafluoroethylene filters (Measurement Technology Laboratories).23 Gravimetrically corrected nephelometric PM2.5 concentrations are used in infants’ indirect PM2.5 exposure estimates. Continuous microenvironmental CO concentrations at 1-minute intervals were logged using the Lascar EL-USB-300 (Lascar Electronics), a large pen-sized device with a sensing range between 0 and 300 ppm. Details on instrumentation, sampling strategy, sample analysis, data process, and quality have been described previously.14, 23
Indirect exposure measurements
Infants’ exposure is estimated using an indirect microenvironmental approach linking objectively determined locations (microenvironments) and assessed PM2.5/CO concentrations in these same areas at the times when the infant was deemed present in the microenvironment with the monitor. We used a Bluetooth®-based beacon system for time-location monitoring among infants. The beacon system consists of 1) emitters (Model O, Roximity Inc. Denver, CO, USA or EM Microelectronic, La Tène, Neuchâtel, Switzerland) and 2) loggers (Berkeley Air Monitoring Group, Berkeley, CA, USA) to monitor the location of infants within a home. Monitored microenvironments include the most commonly occupied rooms (i.e., the main living/sleeping areas, the kitchen), and mothers who consented to be a mobile microenvironment and wore an exposure monitor and beacon logger.23
The emitter (hereafter referred to as a beacon) is a coin-size device that constantly emits Bluetooth signals.20 Infants wore two Beacons each on their clothing during the sampling period. The Beacon logger is a smartphone-sized device that receives and logs Bluetooth signals emitted from the Beacon.20 The beacon logger records the beacon’s unique media access control (MAC) address and the received signal strength indicator (RSSI) every 20 seconds. The RSSI value is proportional to the distance between the Beacon and the Beacon logger and is thus used to determine the infant’s microenvironmental location.20
In each microenvironment, we co-located a fixed-position logger with a PM2.5 and a CO monitor. The time-microenvironment and corresponding area concentrations are assigned at 5-minute intervals as the Beacon logger recorded the strongest average RSSI from the two Beacons worn by infants.20 Reconstructed personal PM2.5/CO exposures for infants were estimated by integrating corresponding microenvironmental concentrations with time spent in the respective locations within the 24-hour monitoring period.
Statistical analysis
All analyses were performed in R (version 4). We calculated descriptive statistics for valid infant indirect PM2.5 and CO exposures by study group (control vs. intervention), study visit (<3 months, ∼6 months, and ∼12 months), and IRC. We evaluated the overall and stove/fuel stratified Spearman’s correlations 1) of measurements of the same pollutants collected at each visit and 2) between PM2.5 and CO exposure at each measurement point. We assessed the differences in PM2.5 and CO exposures between the controls and interventions by study visit (at <3 months, ∼6 months, and ∼12 months) using non-parametric tests (Wilcoxon Rank Sum, Kruskal-Wallis, and Dunn’s tests). Moreover, we evaluated the proportion of children’s exposures that were less than or equal to the WHO Annual Interim Target 1 (WHO IT-1) level of 35 μg/m3 for PM2.5 and the WHO 24-hour guideline value of 4 mg/m3 (∼3.5 ppm) for CO. Additionally, we visually examined the trends in PM2.5 and CO exposures by plotting exposures against the age of infants.
Since each child had up to three exposure measurements, we used mixed-effects models to assess the impact of the intervention on personal exposures. We included a random intercept to account for correlations among repeated measurements on the same child. Exposures are natural log-transformed in the models given their right-skewed distribution. As secondary analyses, we also modeled non-transformed exposures to assess the absolute changes in exposures. The model (Eq. 1) assesses the difference in mean PM2.5/CO exposures between the control and intervention groups during the first year of life.
where yij is the PM2.5 or CO exposure for child i at visit j (<3 months, ∼6 months, or ∼12 months visit), Study Armj is the child’s assigned study group (reference: control), X$ through X$# are indicator variables for the study sites within which the randomization took place (1 in Guatemala, 2 in India, 6 in Peru, and 1 in Rwanda), bi is the random intercept for child i, eij represents the model residual, assumed to be normally distributed. The parameter of interest (Study Arm) was exponentiated, subtracted from one, and multiplied by 100 to obtain the percent reduction in personal exposure due to intervention. In addition, we calculated the intraclass correlation coefficient (ICC) using mixed-effects models with a random intercept for each child to evaluate the within and between individual variability.
Results and Discussion
Household characteristics. 3195 pregnant women randomized to the control (N = 1605) or intervention (N = 1590) groups yielded 3061 live births (control: 1525; intervention: 1536). Among those, 75 exited the study post-birth, and 2986 children remained in the trial (control: 1489; intervention: 1497). Table 1 summarizes selected baseline household and child characteristics by study group and IRC. The household characteristics of children who remained in the trial were similar between the control and intervention groups at baseline trial-wide and within each IRC. About 80% of households had access to electricity at home, which was also reported as the primary light source across the four IRCs. The primary fuel type varied in each IRC, with wood being predominant in Guatemala and India. Cow dung was reported as the primary fuel in Peru. In Rwanda, most households used wood, but a quarter reported charcoal as the primary household energy source. Overall, 47% of the households owned only one stove (mainly in India and Rwanda), and 48% used two stoves (mainly in Guatemala and Peru). Trial-wide only 10% of the households reported having one or more active smokers at home, mainly from India (34%). The mean (SD) gestational age at delivery was 39.3 (1.5) weeks for the control group and 39.4 (1.5) weeks for the intervention group. 4.8% (N = 71) and 5.1% (N = 77) of the births were classified as preterm (defined as gestational age at birth < 37 weeks) in the control and intervention groups, respectively.
Indirect exposure measurements, data completeness, compliance, and quality assurance and control
HAPIN field staff made 7227 exposure monitoring visits for 2986 children during their first year of life. 61% and 63% of the children had at least one valid indirect PM2.5 and CO exposure sample, respectively. The proportion of (a) enrolled children with valid indirect measurements and (b) exposure visits with valid indirect measurements vary by study visit and IRC (Table S1-S4). Generally, the percentage of valid samples went down over time; the <3 month visit had the most valid samples, and the ∼12 month visit had the least. The moderate to low percentage of valid samples is partially due to increased equipment failure as time passed, field measurements interrupted by the COVID-19 pandemic lockdown, and the challenge and complexity of successfully merging multiple real-time data sources in the indirect microenvironmental approach.
Child exposure summary
24-hour child indirect exposures to PM2.5 and CO are summarized by study arm and visit in Table 2 and in Figure 1. Exposure visits occurred as planned, with the first visit occurring on average 3.2 (SD 0.43) months after birth, the second visit occurring on average 6.3 (SD 0.48) months after birth, and the final visit occurring 12.1 (0.49) months after birth (SI Figure S1).
Trial-wide and within IRCs, we observed stable reductions in PM2.5 between post-birth measurement rounds. The overall median (IQR) child PM2.5 and CO exposures are 65.1 (33.0 — 128.2) µg/m3 and 1.1 (0.3 — 2.9) ppm for the control group, and 22.9 (17.2 — 35.3) µg/m3 and 0.2 (0.0 — 0.7) ppm for the intervention group. Median child exposures to PM2.5 in the intervention arm were 65% lower than in the control arm during the <3 months visit (67 vs. 23 µg/m3), 62% lower at the ∼6 months visit (60 vs. 23 µg/m3), and 69% lower at the ∼12 months visit (72 vs. 22 µg/m3). CO exposures were similarly reduced: 85% at <3 months (1.1 vs. 0.2 ppm), 82% at ∼6 months (1.1 vs. 0.2 ppm), and 88% at ∼12 months (1.0 vs 0.1 ppm). Across all rounds, 26% of control samples and 75% of intervention samples were below the WHO-IT1 target guideline value of 35 µg/m3. Similarly, 95% of intervention samples and 79% of control samples were below the WHO annual guideline value for CO (3.5 ppm). Both descriptive statistics of PM2.5 and CO exposures and the proportion of measurements below the WHO guidelines were similar across the three post-birth visits.
Exposures over time
We plotted children’s exposures to PM2.5 by age (the number of days between the exposure measurement date and the child’s birthday, Figure 2). We observe consistent and distinct separation of exposures between the control and intervention children in the first year of life. Control and intervention samples were significantly different at each round for each pollutant (p < 0.001). For control households, there was no significant difference between PM2.5 or CO exposures by study visits. Exposures were 109 (SD 125), 112 (SD 147), and 124 (SD 162) for <3 months, ∼6 months, and ∼12 months visits. Similar trends were observed for CO: 2.4 (SD 4.0), 2.6 (SD 4.6), and 2.7 (SD 5.1) ppm at < 3 months, ∼6 months, and ∼ 2 months. There were no significant differences between PM2.5 or CO exposures by study round for intervention households. Exposures were 39 (SD 120), 35 (SD 69), and 39 (SD 88) for visits at <3 months, ∼6 months, and ∼12 months. Similar trends were observed for CO: 0.8 (SD 2.3), 0.7 (SD 1.5), and 0.9 (SD 2.6) ppm at <3 months, ∼6 months, and ∼12 months.
Correlations between pollutants and between study rounds
Correlations between pollutants were moderate for the control arm (Spearman’s ρ 0.45) and poor for the intervention arm (Spearman’s ρ 0.08). Overall correlations between study rounds (depicted in SI Figure S2) were weak for CO (Spearman’s ρ 0.37) and were moderate for PM (Spearman’s ρ 0.58); for both pollutants, they were stronger in the control arm than in the intervention arm.
Modeling the effect of the LPG stove and fuel intervention on personal exposures
Modeled impacts of the HAPIN LPG fuel and stove intervention on children’s exposure indicated significant reductions for PM2.5 and CO (Figure 3) of 62% (95% CI 59, 64) and 80 (95% CI 77, 83), respectively. In models with untransformed outcomes, these percent reductions translate to predicted absolute reductions of 79.9 µg/m3 (95% CI 69.7, 90.1) and 1.8 ppm (95% CI 1.5, 2.1) for PM2.5 and CO, respectively.
Discussion
Our findings add to evidence of the effectiveness of a well-used and highly adhered to LPG intervention for lowering exposures to household air pollution, measured here as PM2.5 and CO.12, 13 Our findings are consistent with previous work done as part of the HAPIN trial to characterize exposures during pregnancy; among the same households, we find remarkable consistency between exposure reductions during pregnancy and for young children, as reported here. As we reported previously for pregnant women, these exposures are, to the best of our knowledge, among the most substantial for a household energy intervention to date. Our findings of reduced exposures are consistent and stable across our four low-and middle-income settings.
Comparisons with measurements made during HAPIN are especially relevant, as they provide evidence of stable reductions throughout different phases and measurement regimens of the trial. Prior to the birth of the child, maternal exposure to HAP was measured once at baseline, prior to randomization, and at most twice during pregnancy. Findings from those initial measurements of the intervention’s effectiveness indicated similar reductions in PM2.5 and CO exposures (61% and 81%, respectively, as compared to 62% and 80% in the current study, Figure 4). Taken in total, these findings suggest the HAPIN intervention reduced exposures in a significant and consistent pattern throughout pregnancy and the child’s first year of life.
Several studies have measured personal exposure to HAP among children under five. Early studies in Guatemala,30, 31 Gambia,32–34 and South Africa35 measured CO using passive diffusion tubes attached to children’s clothing. Recent studies used personal monitors to measure children’s exposures to PM2.5 in Malawi,36 the United States,37 Rwanda,38 Ghana,39 Vietnam,40 and CO in Ghana,8, 17 Peru,41 Rwanda,42 or both Uganda and Ethiopia.43 Among identified studies, four measured PM2.5/CO exposures using personal monitors in children under 1-year old. In these studies, sampling instruments were all attached to the customized clothing of the caregiver.8, 17, 36, 41, 42 Several direct measurements were conducted among children between ages 1 and 5 with lightweight personal monitors in a backpack/shoulder bag or a vest pocket.37, 38, 40 Studies also used microenvironment measurements plus time-activity logs and statistical modeling to estimate the exposures of children.32, 39, 43
Despite these different approaches, reported values are uniformly high among LMIC households using solid fuels. Kirby et al. (2019) reported high child PM2.5 exposure in Rwanda, where the median (IQR) exposure was 154 (85–267) µg/m3 in the intervention households using an improved biomass stove and 161 (91–285) µg/m3 in control households.38 In Malawi, the geometric mean (GM) (95% CI) of PM2.5 exposure among children at six months was 53.9 (49.4 - 58.7) in a cleaner-burning biomass-fueled cookstove intervention group and 49.0 (45.0 - 53.3) in the control group using open fires.36 Similarly, the GM (GSD) of young infants’ exposure to PM2.5 in households using biomass was 80.2 (1.34) and 97.0 (1.89) in Uganda and Ethiopia, respectively.43 Clean fuels resulted in lower exposures; in Vietnam, the reported mean (SD) child exposure was 36 (24) µg/m3 in households using electricity for cooking and 44 (38) µg/m3 in homes using LPG.40
Personal CO exposure of children under five old followed a similar pattern. The 48-hour median (IQR) personal exposure observed in Ghana was 1) 0.48 (0.17–1.23) ppm in traditional cookstove households, 2) 0.51 (0.17–1.10) ppm in improved biomass cookstove households, and 3) 0.39 (0.12–0.94) ppm in LPG stove households.8 In Rwanda, the estimated child CO exposure from burning solid fuels was much higher, with a 24-hour mean of 6.8 ppm.42 The reported 24-hour GM was 1.4 in Peru,41 0.64 in Uganda, and 2.4 in Ethiopia.43
Our study has a number of strengths. It contributes to this growing body of evidence of HAP exposures among young children. It is among a small number that has used sensor-based measurements to assess microenvironmental locations and coupled this time-activity data with microenvironmental PM2.5 and CO measures, enabling exposure reconstruction. We also acknowledge a number of limitations. While we highlight the reduced exposures for both expecting mothers and their young children resulting from the HAPIN intervention, we acknowledge fairly large, unexplained residual exposures. Future work should focus on identifying the drivers of this residual exposure, which may include outdoor exposures from trash burning, from unmonitored biomass stoves, or from other household energy sources (animal food preparation, heating water for bathing, etc). Second, we acknowledge the unique circumstances of the trial, in which LPG was supplied for free, and during which households were provided behavioral nudges to continue its use, neither of which will continue post-trial. Third, despite the scale of the exposure assessment exercise reported here, we only reconstructed periodic snapshots of exposure that may not represent the ‘true’ participant experience. However, we note the remarkable stability of both the control arm and intervention exposures post-randomization (Figure 4) for mothers and children. Next, microenvironmental reconstructions required the use of numerous data streams linked by unique beacons, beacon loggers, and household IDs. Loss of any single data stream – from data processing errors, file naming errors, or mis-entry during digital record taking while deploying samples – resulted in lost samples. Future work will describe efforts to model both mother and children’s exposure using limited subsets of the data collected during HAPIN and leveraging repeat measurements in the same household, where available. Finally, we acknowledge that the study children’s behavioral patterns at the three measurement timepoints likely are not representative of each other nor of later periods of childhood and adolescence. Nonetheless, the findings provide evidence (and rationale) for these types of measurements and for further investigation of exposures during critical periods of child development.
While the HAPIN study is an efficacy trial under which LPG was provided continuously and for free, our findings, coupled with previous findings for pregnant women14 and forthcoming findings for non-pregnant adult women from the same households (in preparation), provide evidence that adoption and sustained, near-exclusive use of clean fuels can reduce exposure to household air pollution for multiple household members, including vulnerable young children. Given the potentially deleterious impacts of exposure across the lifespan, additional measurements are warranted for children and other household members whose experience of pollution is not currently well characterized.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Acknowledgements
The HAPIN trial is funded by the U.S. National Institutes of Health (cooperative agreement 1UM1HL134590) in collaboration with the Bill & Melinda Gates Foundation [OPP1131279].
The investigators would like to thank the members of the advisory committee - Patrick Breysse, Donna Spiegelman, and Joel Kaufman - for their valuable insight and guidance throughout the implementation of the trial. We also wish to acknowledge all research staff and study participants for their dedication to and participation in this important trial.
A multidisciplinary, independent Data and Safety Monitoring Board (DSMB) appointed by the National Heart, Lung, and Blood Institute (NHLBI) monitors the quality of the data and protects the safety of patients enrolled in the HAPIN trial. NHLBI DSMB: Nancy R. Cook, Stephen Hecht, Catherine Karr (Chair), Joseph Millum, Nalini Sathiakumar, Paul K. Whelton, and Gail Weinmann and Thomas Croxton (Executive Secretaries). Program Coordination: Gail Rodgers, Bill & Melinda Gates Foundation; Claudia L. Thompson, National Institute of Environmental Health Sciences; Mark J. Parascandola, National Cancer Institute; Marion Koso-Thomas, Eunice Kennedy Shriver National Institute of Child Health and Human Development; Joshua P. Rosenthal, Fogarty International Center; Concepcion R. Nierras, NIH Office of Strategic Coordination Common Fund; Katherine Kavounis, Dong-Yun Kim, Antonello Punturieri, and Barry S. Schmetter, NHLBI.
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the US National Institutes of Health or Department of Health and Human Services.
Footnotes
↵# (co-first)
Synopsis
Household air pollution (HAP) exposure assessment rarely includes estimates for infants and young children. We contrast HAP exposures between children in households with a clean fuel intervention and those using biomass-fueled stoves.