Skip to main content



Globally, nearly 3 billion people rely on solid fuels for cooking and heating, the vast majority residing in low- and middle-income countries (LMICs). The resulting household air pollution (HAP) is a leading environmental risk factor, accounting for an estimated 1.6 million premature deaths annually. Previous interventions of cleaner stoves have often failed to reduce exposure to levels that produce meaningful health improvements. There have been no multicountry field trials with liquefied petroleum gas (LPG) stoves, likely the cleanest scalable intervention.


This paper describes the design and methods of an ongoing randomized controlled trial (RCT) of LPG stove and fuel distribution in 3,200 households in 4 LMICs (India, Guatemala, Peru, and Rwanda).


We are enrolling 800 pregnant women at each of the 4 international research centers from households using biomass fuels. We are randomly assigning households to receive LPG stoves, an 18-month supply of free LPG, and behavioral reinforcements to the control arm. The mother is being followed along with her child until the child is 1 year old. Older adult women (40 to <80years of age) living in the same households are also enrolled and followed during the same period. Primary health outcomes are low birth weight, severe pneumonia incidence, stunting in the child, and high blood pressure (BP) in the older adult woman. Secondary health outcomes are also being assessed. We are assessing stove and fuel use, conducting repeated personal and kitchen exposure assessments of fine particulate matter with aerodynamic diameter 2.5μm (PM2.5), carbon monoxide (CO), and black carbon (BC), and collecting dried blood spots (DBS) and urinary samples for biomarker analysis. Enrollment and data collection began in May 2018 and will continue through August 2021. The trial is registered with (NCT02944682).


This study will provide evidence to inform national and global policies on scaling up LPG stove use among vulnerable populations.



Globally, nearly 3 billion people rely on solid fuels (wood, dung, coal, charcoal, or agricultural crop waste) for cooking and heating (Bonjour et al. 2013). These fuels are often burned in inefficient and poorly ventilated combustion devices (e.g., open fires, traditional stoves). The resulting household air pollution (HAP) accounts for an estimated 1.6 million premature deaths per year and 59.5 million disability-adjusted life-years (GBD 2018). Despite progress in recent years, this largely preventable exposure remains a leading risk factor for morbidity and mortality worldwide. Poor populations in low- and middle-income countries (LMICs) bear most of this burden (GBD 2018).
Several studies have documented association between HAP and multiple diseases or health conditions including chronic lung disease, lung cancer, cancers of the aerodigestive tract, cervical cancer in adults, pediatric acute lower respiratory infections (ALRI), or pneumonia, low birth weight, stillbirth, preterm birth, childhood stunting (short length for age), tuberculosis, and impaired cognitive development (Bruce et al. 2015; Smith et al. 2014; Gordon et al. 2014; Quansah et al. 2017; Thakur et al. 2018). However, except for pediatric ALRI, current methods used for the global HAP-attributable disease burden calculations are focused primarily on cardiovascular disease (CVD) and chronic respiratory diseases in adults using integrated exposure–response functions drawn substantially from estimates of effects of other sources of air pollution (Burnett et al. 2014). Additionally, adverse birth outcomes such as low birth weight and preterm birth, along with childhood development and growth, are not included in the global burden of disease estimates. Therefore, the current estimate for the burden of disease related to HAP is uncertain and likely underestimated.
The state of the science illustrates several compelling reasons to undertake a multicountry randomized controlled trial (RCT) for reducing HAP using a clean fuel intervention. We will address several knowledge gaps in this study. First, liquefied petroleum gas (LPG) is currently the most widely available clean fuel in LMICs (IEA 2017), but to date, no LPG trials have been conducted that demonstrate significantly improved health outcomes among children and adults. Currently available cleaner-burning biomass combustion stoves are unlikely to achieve or sustain health-relevant exposure reductions (Clark et al. 2013; Bruce et al. 2015; Anenberg et al. 2013; Sambandam et al. 2015). This trial will provide needed evidence from a randomized LPG stove intervention to support policy formulation on national levels. Second, focusing on a combination of child (birth weight, pneumonia, and stunting) and adult cardiovascular [blood pressure (BP)] primary outcomes is strategic for public health goals in LMICs, as these conditions contribute the most to HAP-associated health burdens in these settings (GBD 2018). Although HAP has been identified as a risk factor for these outcomes, intervention efforts directed at birth outcomes, child health, and noncommunicable disease risk have been included in a few RCTs aimed at reducing HAP. Furthermore, we will evaluate biomarkers and other indicators that are known to predict noncommunicable disease occurrence and/or severity across the lifespan. Third, establishing exposure–response relationships across a range of personal exposures to HAP is needed to close critical gaps in our current understanding of exposure to disease relationships (Steenland et al. 2018). Exposure–response relationships are also critical for transferability of trial results across settings and for benchmarking future intervention efforts. Finally, most stove intervention studies have failed to adequately investigate and address the behaviors necessary to overcome concurrent use of polluting stoves (a practice known as stacking) to ensure consistent and sustained use of cleaner stoves and displacement of polluting ones (Rosenthal et al. 2017).
This paper summarizes the rationale, study design, and methods of the Household Air Pollution Intervention Network (HAPIN) trial, a recently launched RCT that seeks to provide the evidence necessary for policy makers to determine what health benefits that can be achieved by implementing a scalable (in many areas) intervention aimed directly at reducing HAP in LMICs. The trial represents the first multicountry RCT to assess the effect of a stove intervention LPG on exposure to HAP and on a broad range of maternal, child, and adult health outcomes.

Study Aims

The study has three specific aims. The first is to determine the effect of a randomized LPG stove and fuel intervention on health in four diverse biomass-using LMIC populations across the world using a common protocol. We hypothesize that compared to pregnant women (18 to <35years of age) in control households (n=1,600), those who receive LPG stoves and fuel (n=1,600) will have offspring with increased birth weight, reduced severe pneumonia incidence, and improved growth [less stunting, defined as length-for-age z-score less than 2 standard deviations (SD) below the median z-score based on World Health Organization (WHO) child growth standards] up to 12 months of age. We also hypothesize that compared to control households, older adult women (40 to <80years of age) living in households that receive LPG stoves and fuel will have reduced BP. In addition to these primary outcomes, the study will assess multiple secondary outcomes on mothers, infants, and older adult women.
The second aim is to evaluate the exposure–response associations for HAP and health outcomes in four diverse LMIC populations. Using repeated 24-h personal and indirect measurements of exposure to fine particulate matter with aerodynamic diameter 2.5μm (PM2.5), carbon monoxide (CO), and black carbon (BC), we will characterize the exposure–response associations for all primary and secondary outcomes (assessing potential nonlinearity) while adjusting for confounders. While evidence for an overall effect of the intervention will be available from the first aim, analysis of exposure–response is critical for quantitative risk assessment and policy determinations of acceptable levels of HAP regardless of the cooking technologies and/or fuels in use.
The third aim is to evaluate the extent to which biomarkers of exposure and health effects, including targeted and exploratory (e.g., metabolomics) analyses, are associated with intervention status or exposure. We hypothesize that participants residing in households that receive LPG stoves and fuel or have lower levels of exposure to HAP will have lower levels of carcinogenic polycyclic aromatic hydrocarbons and volatile organic compounds, such as urinary 1-OH-pyrene, 2-naphthol, 9-phenanthrene, as well as of chronic disease indicators, such as inflammatory, endothelial, inflammatory, oxidative stress, and glycemic control/diabetes biomarkers [e.g., C-reactive protein (CRP), endothelin-1, E-selectin, interleukin 6, and hemoglobin A1c (HbA1c)] when compared to participants in control households.
By comparing an LPG stove and continuous free fuel intervention with standard cooking practices (typically traditional solid biomass in these settings), the HAPIN trial will provide an estimate of a potentially achievable level of HAP reduction and the associated impact on select maternal, infant, and adult health outcomes. Establishing exposure–response relationships in LMIC settings will allow for estimates of the range of improvement that can be expected across real-world conditions where clean stoves and fuels are often combined with traditional biomass stoves.


Study Overview

The study is an RCT of an LPG stove and continuous fuel distribution intervention and promotion of its exclusive use among 3,200 households in four LMICs (India, Guatemala, Peru, and Rwanda). Following an 18-month period of planning, piloting, and formative research, the study began recruiting participants in May 2018 and is expected to complete enrollment in February 2020; follow-up data collection will continue through August 2021. In each country, eligible pregnant women are recruited and their households randomly assigned to intervention and control groups on a 1:1 ratio, and they are followed for 18 months until their newborn child is 1 year old. Intervention households receive a free LPG stove and free unlimited supply of LPG for the 18-month follow-up period. Control group households do not receive an LPG stove and fuel during the study period, and it is anticipated that they will continue cooking with solid biomass fuels during the trial. After enrollment, assessments will be made on a regular schedule over the course of the pregnancy (baseline, 24–28 wk gestation, 32–36 wk gestation), at 3 months of age, 6 months of age, and 12 months of age for the child, and at the same time points for the older adult woman in the household (Table 1). Control group compensation is summarized below and described elsewhere (Quinn et al. 2019).
Table 1 Schedule of exposure and outcome assessment.
Child age (study time point)<20wk gestation (baseline)24–28 wk gestation (1–3 months postrandomization)32–36 wk gestation/birth (3–5 months postrandomization)3monthsold (9months postrandomization)6monthsold (12months postrandomization)9monthsold (15months postrandomization)12monthsold (18months postrandomization)
Personal exposure
 Urinary biomarkersX XX XX X XX XX    XX
Primary outcomes
 Birth weight       Xa             
 Severe pneumoniab                     
 Stunting                   X 
 Blood pressure  X  X  X  X  X     X
Secondary outcomes
 Maternal blood pressureX  X  X  X  X     X  
 Fetal growth X  X  X             
 Child linear growth (continuous)       Xa  X  X  X  X 
 Preterm birth       Xa             
 Child development          X  X  X  X 
 WHO severe pneumoniaa                     
 BART  X                 X
 CIMT  X                 X
 SGRQ  X                 X
 SF-36  X                 X
 Expenditures/time useX X                 X
 Chronic disease biomarkers  X  X  X  X  X     X
 Metabolomics/microRNAc             XX      
 Hemoglobin (anemia)X  X  X      X     Xa 
 SociodemographicsX X                 X
 Clinical historyX XX XX XXXX XX XX XX
 Weight/height/BMIX XX XX X  X  X     X
 Diet/food securityX X               X X
 IYCF          X  X  X  X 
Biospecimens collected
 UrineX XX XX X XX XX    XX
 Dried blood spotsX XX XXXaX XX XX    XX
Note: BART, brachial artery reactivity testing; BC, black carbon; BMI, body mass index; C, child; CIMT, carotid intima-media thickness; CO, carbon monoxide; IYCF, infant and young child feeding practices; NM, new mother; O, older adult woman in household; PM2.5, fine particulate matter with aerodynamic diameter 2.5μm; PW, pregnant woman; SF-36, Short Form 36 survey; SGRQ, St. George Respiratory Questionnaire.
Measured at birth.
Recorded whenever children present to HAPIN health facilities with respiratory symptoms.
Metabolomics/microRNA biomarker discovery in 100/site for the older adult woman and child.

Study Settings and Formative Research

The study is being conducted across four LMIC settings in which large portions of the population use solid biomass as the primary fuel type. To increase generalizability, the settings were purposefully selected to represent a diversity of characteristics expected to influence intervention effects, including altitude, population density, cooking practices, baseline pollution levels, and sources of pollution other than cooking (Table 2). Other factors that may influence intervention effects, such as fuel types, dwelling characteristics, and socioeconomic conditions, are being measured and recorded. Within each country, candidate sites were selected after evaluation in formative research over 12 months. The formative research consisted of four phases: a) initial scoping to identify potentially suitable sites and developing contextually grounded behavior change strategies to promote intervention adoption, b) pilot intervention to determine the HAP exposure contrast that might be expected from the intervention, c) pilot assessment to test trial procedures and methods, and d) respiratory rate/pulse oximetry assessment to define context-specific tachypnea and oxyhemoglobin saturation thresholds in the study sites.
Table 2 Summary of key characteristics of international research centers by country based on sampling, government information, or published studies.
Tamil Nadu, selected blocks:
Kalrayan Hills (Villupuram District); Kezhvelur, Keelayur, Vedaranyam, and Thalaignayiru (Nagapattinam District)
Jalapa municipality:
Regions: Santa Maria Xalapan, Ladinos Pardos
Department of Puno:
Provinces of Puno, San Roman, Azángaro, Huancané, El Collao, and Chicuito
Eastern Province:
Kayonza District:
Selected sectors: Kabare, Kabarondo, Murama
Altitude (m above sea level)10–165871–2,6773,8251,300–1,700
Population density (perkm2)892–1,09823317.6274
Cooking practices (stoves, fuel, location)Traditional plastered clay/mud stoves fueled with biomass, with 90% cooking indoors. Use of chimneys is negligible.Chimney stoves and open fires; 97% wood use; cooking mainly indoorsRural households use traditional biomass (typically dung)-burning stoves daily for cookingTraditional three-stone fires (62.6%) or Rondereza (34.6%) fueled with wood (89.9%) or charcoal (8.1%), cooking indoors (71.9%)
Baseline 24-h PM2.5 concentration (μg/m3)Mean:
Kitchen area: 210
Personal: 155
Median (Q1–Q3):
Personal: 115 (80–265)
Kitchen: 349 (244–524)
Bedroom: 25 (17–120)
Patio: 30 (16–50)
Median indoor: 130Mean 24 h:
Personal: 329
Kitchen: 437
Outdoor: 44
Other sources of air pollutionIncense, mosquito coils, etc. (55% of households), trash burning (38% of households)Trash burning and seasonal burning crop residue (February–April)Low ambient air pollutionMinimal unless near road or in more urban/dense area.
Prevalence of smoking in the homeWomen rarely smoke. Other smokers reported in 24% of households but smoking inside house reported in less than 1% of households.Women rarely smoke (<1% ever smoked). Secondhand smoke in 20% of homes, but limited to 1–2 cigarettes/daySelf-reported daily smoking is 0.2%.Province level:
Never: 76.4%
Daily: 21.9%
Note: Data are from Guatemala: Fujisada H et al. 2012; Government of the Republic of Guatemala 2019; Johnson M et al. 2018. Peru: Pollard SL et al. 2014; Jaganath D et al. 2015; Hollada J et al. 2017. India: Balakrishnan et al. 2013, 2018. Rwanda: National Institute of Statistics of Rwanda, 2010, 2012; unpublished data from Kirby et al. 2016; HAPIN formative research results (unpublished). PM2.5, fine particulate matter with aerodynamic diameter 2.5.

Eligibility Criteria, Screening, and Recruitment

Study teams led by experienced local investigators work in collaboration with clinics and community health workers in each country to identify candidate pregnant women. To be eligible to participate in the study, a pregnant woman must meet the following inclusion criteria: confirmed pregnancy (human chorionic gonadotropin–positive blood or urine test); 18 to <35years of age (confirmed by government-issued ID, whenever possible), cooks primarily with biomass stoves, lives in the study area, 9 to <20wk gestation with a viable singleton pregnancy with normal fetal heart rate confirmed by ultrasound, continued pregnancy at the time of randomization (via self-report), and agrees to participate with informed consent. Eligible pregnant women are excluded if they currently smoke cigarettes or other tobacco products, plan to move permanently outside the study area in the next 12 months, use a clean fuel stove predominantly, or are likely to use LPG or another clean fuel predominantly in the near future. Ultrasound measurements are conducted by trained personnel (who are also additionally certified centrally) in a clinic or home setting to determine eligibility and assess fetal growth using a portable ultrasound [Edge (Edge Ultrasound System), Sonosite/Fujifilm Edge (FUJIFILM SonoSite Inc.)].
Across the trial locations, up to 800 older adult women 40 to <80years of age (confirmed by government-issued ID whenever possible) who reside in the same households as an enrolled pregnant woman are being recruited (one per household), provided they do not fall within the following exclusion criteria: currently smoking cigarettes or other tobacco products, pregnant (via self-report), or planning to permanently move out of current household in the next 12 months.

Baseline Surveys and Assessments; Randomization

Following recruitment and informed consent, a baseline visit is made to the household by a trained fieldworker to conduct surveys and other assessments. This baseline visit includes a survey that covers a range of topics like cooking behaviors, household composition, socioeconomic and demographic information, housing characteristics, and pregnancy-related information. Pregnant women are also surveyed about their medical and gynecological history, including medication use. Separate questionnaires assess physical activity, dietary diversity, household food insecurity, and household expenditures.
The baseline visit includes assessments of health, biomarkers, and exposure for both the pregnant woman and older adult woman. For pregnant women, a trained fieldworker or nurse measures resting BP (model HEM-907XL; Omron®) in triplicate and maternal weight (seca 876/874 scales; Seca) and height (seca 213 stadiometer; Seca) in duplicate. Among older adult women, respiratory symptoms are assessed using the Saint George’s Respiratory Questionnaire (SGRQ) (, and quality-of-life measures are assessed using the Short Form 36 (SF-36) ( Further, carotid intima-media thickness (CIMT) and, in Peru, only due to its exploratory nature, brachial artery reactivity testing (BART) are measured among older adult women using the ultrasound devices described above. Common carotid artery (CCA) ultrasound will be performed as a marker of atherosclerosis using a high-resolution linear transducer to image the distal 1-cm CCA region (just proximal to the bifurcation) (Stein et al. 2008). CIMT will be obtained from end-diastolic B-mode images as the average of the posterior wall segments from both right and left CCAs, measured using an automated system with an edge detection algorithm and manual override capacity (100 separate dimensional measurements are obtained from the 1-cm segment and averaged to obtain mean and maximal CIMT values). BART will be performed with a portable ultrasound to assess endothelial function using a high-resolution linear transducer to image the brachial artery (BA) above the antecubital fossa to measure rest diameter (Corretti et al. 1995). The BA will be occluded for 5 min by a BP cuff inflated to suprasystolic levels in the forearm; after 5 min of occlusion, the BP will be released, and the BA diameter will be imaged every 30 s for 2.5 min after cuff deflation to calculate percent BA dilation after hyperemia.
Urine samples (first morning void) and dried blood spots (DBS) via finger prick are obtained from all participating pregnant women and older adult women. Finally, pregnant women and the older adult women are monitored for 24-h personal and household level exposure to HAP (PM2.5, CO, and BC) using the procedures described below. Baseline data on ambient temperature and humidity in the home and primary cooking area are also collected by the data loggers of the PM2.5 samplers, and dimensions of the kitchen are measured by the surveyor.
After the baseline surveys and assessments are completed, households are randomly assigned to intervention or control arms, stratified by country. In India and Peru, additional stratified randomization is used to ensure a balance between discrete geographical regions within the study area. In Rwanda and Guatemala, the study areas are deemed homogenous so that such further stratification is not necessary.


The intervention consists of an LPG stove, a continuous supply of LPG fuel delivered to the homes for 18 months, as well as education and behavioral messages (described in future papers) to promote safe, exclusive use of the LPG stove for cooking. Stoves are procured in each country and vary (footprint, base, burner size, nob position, and griddle) based on local cooking practices; however, all include at least two burners and meet applicable safety requirements. The intervention (stove and fuel) is provided free of charge to all intervention households after baseline measurements are conducted. On each visit to provide additional fuel cylinders in Rwanda, Guatemala, and Peru, stove condition is examined, any necessary repairs performed, and the weight of LPG cylinders measured and recorded in order to help monitor use and anticipate the need for additional refills. In India, per national governmental regulations, the public sector oil marketing company is responsible for stove installation, cylinder refills, and repairs; study staff facilitate those activities for intervention participants. The rate of LPG usage is monitored by calculating average kilograms of LPG used per household member per day (using data on fuel cylinder weights and the number of days between installation and exchange of each cylinder).

Control Compensation

Control households receive compensation designed to meet three aims. First, it must comply with applicable ethics requirements for treatment of controls. Second, we are compensating control participants for the burden associated with this study, with the view of minimizing losses to follow-up. Third, we offset the economic advantage to intervention households accorded by the provision of free stoves and fuel. While the details vary across the four countries, compensation was designed based on a uniform set of trial-wide principles that address the above aims, with details informed by focus group discussions in the communities selected for the intervention. Controls receive either an LPG stove and a supply of fuel at the end of the trial or preferred alternatives of comparable value during or at the end of the trial. Details concerning the development of the compensation strategy are provided elsewhere (Quinn et al. 2019).

Primary Outcomes

Our primary health outcomes are birth weight, severe pneumonia in the first 12 months of life, stunting at 12 months of age, and BP in the older adult woman over 18 months of follow-up.
Birth weight is measured in duplicate to the nearest gram within 24 h of birth by a trained fieldworker or nurse using a routinely calibrated seca 334 mobile digital baby scale (Seca). If the first two measurements differ by more than 10g, a third measurement is taken. Newborns are weighed naked or in a preweighed blanket, typically at the health facility where infants are delivered.
For the purposes of the HAPIN trial, a case of severe pneumonia is adapted from the revised WHO classification of childhood pneumonia (WHO 2014b) and is consistent with expert opinion (Goodman et al. 2019; WHO 2014a). This classification includes two independent algorithms for severe pneumonia, which we have enriched with more objective imaging and pulse oximetry criteria (Simkovich et al., in press). The first algorithm requires the presence of cough and/or difficulty breathing, at least one general danger sign, and primary end point pneumonia on either lung ultrasound or chest radiograph imaging. Lung ultrasound is our preferred imaging modality; however, if logistical reasons preclude lung ultrasound, then we will consider chest radiograph images interpreted according to WHO methodology. The second severe pneumonia algorithm requires the presence of cough and/or difficulty breathing and hypoxemia, measured noninvasively by pulse oximetry. We will define hypoxemia based on the physiologic threshold of 92% for altitudes <2,500m above sea level and 86% for altitudes 2,500m above sea level.
BP is measured in the right arm of an older adult woman in triplicate (with at least 2 min resting between repeat measurements) using an automatic BP monitor (model HEM-907XL; Omron®). The first measurement is ignored to reduce white coat hypertension. The average of the second and third measurements will be used for analysis. Measurements are taken after receiving assurances that the older adult woman has not smoked, consumed alcohol or a caffeinated beverage (coffee, tea, or cola), or cooked using biomass in the 30 min prior to the measurement. She is asked to sit in a chair in a quiet room for 5 min with legs uncrossed, back supported by a chair, and arm supported by a table prior to commencing the measurements.

Secondary Outcomes

Secondary outcomes are maternal BP, preterm birth, fetal growth, infant linear growth (as a continuous outcome), infant development using the long-form version of Caregiver Reported Early Development Instruments (CREDI) (McCoy et al. 2017), WHO severe pneumonia among children <12months old, WHO pneumonia (both severe and nonsevere) confirmed by ultrasound in Guatemala, burns, CIMT, BART (Peru only), SGRQ (respiratory symptoms), SF-36 (quality of life), household expenditures for fuel and healthcare, household time/activity, and chronic disease biomarkers in the older adult woman. Detailed methods on these outcomes will be included in the papers reporting the results.

Exposure Assessment

Two approaches to measure personal HAP exposure are used. For pregnant women (baseline, 24–28 wk gestation, and 32–36 wk gestation) and older adult women (baseline and five additional times during follow-up) (Table 1), exposure is measured through instrumentation placed on the participant. Pollutants measured include PM2.5, CO, and BC. PM2.5 concentrations are assessed using both gravimetric and real-time methods, and the PM2.5 filters are being assessed for BC with SootScan™ Model OT21 Transmissometers (Magee Scientific). The primary instrument for this purpose is the Enhanced Children’s MicroPEM™ (ECM) developed by RTI International (Johnson et al. 2020). Real-time PM2.5 concentrations are estimated with the ECM’s nephelometer. Real-time CO concentrations are measured personally for all pregnant women and older adult women (and for infants in Guatemala) and area monitoring locations with EL-USB-300 CO monitors (Lascar Electronics).
For infants for whom the ECM is impractical to use due to their size or participant preference, a microenvironmental approach is employed at 3, 6, and 12 months of age (Balakrishnan et al. 2004; Xu et al. 2018; Zuk et al. 2007), whereby area sampling via ECMs is used in main living/sleeping areas, including the kitchen. The child is outfitted with a coin-sized location beacon (Roximity) linked to receivers in these same areas to objectively assess location. Personal exposures for the child are estimated by integrating corresponding area concentrations with time spent in the respective locations along with time spent near the mother, who is instrumented with an ECM and receiver whenever feasible. Details of this approach can be found in Liao et al. (2019), which presents results from piloting the indirect assessment with proximity beacons for the HAPIN trial in Guatemala.
To track changes in ambient PM2.5 concentrations and broadly characterize potential for regional air quality impacting participant exposure, ambient measurements are being made in at least two locations at each international research sites (IRCs) field site. Real-time ambient concentrations of PM2.5 are measured using the E-Sampler (Met One Instruments) or comparable systems with an integrated 24-h filter-based measurement.
A full description of the exposure assessment procedures is described elsewhere (Johnson et al. 2020).

Stove Use Monitoring

Stove use is monitored throughout the 18-month follow-up period in both intervention and control households using a combination of observations, reports, and instruments. Use of the LPG stoves in intervention arm households is also tracked via data on the distribution and use of LPG cylinders. Additionally, all homes in the intervention arm, and a subset of 80 homes in the control arm of each country, are equipped with stove use monitors (SUMS), which are temperature data loggers placed on traditional stoves (Ruiz-Mercado et al. 2012; Pillarisetti et al. 2014). SUMS data by household, supplemented with observations made by local study teams, are reviewed weekly to identify any continued traditional stove usage. Where it occurs, local teams, following trial-wide written standard operation procedures (SOPs), visit each home to address barriers to exclusive LPG stove use that may be responsible and reinforce behavioral messages developed in relation to local cooking needs to minimize future events. For example, barriers related to capabilities and skills are addressed by how-to training, whereas barriers related to motivation are targeted with an appeal to emotions such as trust and security and conscious decision-making. Barriers related to opportunity and context are addressed with suggestions for adapting to contextual realities. The 80 control homes per country comprise the 20% random subsample in which additional exposure measurements are taking place (as described above). The SUMS are Geocene Temperature Loggers™, which attach high-temperature thermocouples to the various stoves used in a home.


Biomarker analyses are an integral part of the trial, drawing on existing biomarker research and developing and validating novel biomarkers associated with HAP exposure and health. The Biomarker Center is based at Emory University (Emory), where samples from three international research centers are analyzed. Samples from India are principally analyzed at a laboratory based in India, where extensive infrastructure and capacity already exists, and that laboratory collaborates closely with Emory in methods and data analysis. Using urine and DBS collected at regular intervals from the pregnant woman, older adult woman, and child (Table 1), the Biomarker Center will: a) assess repeated measures, including targeted biomarkers of exposure, biomarkers of tobacco smoke exposure, biomarkers of effect that are predictive of clinical outcomes, and exploratory analyses that include metabolomics, mRNA, miRNA, and DNA methylation; b) perform HAP-specific biomarker development and validation; c) create and operate a biospecimen repository; and d) create a data set for deposit in the National Institutes of Health (NIH) BioLINCC repository and ensure that study procedures are compatible with that program’s requirements for biospecimen collection, labeling, and storage. Biomarker sampling methodologies are described in detail elsewhere (Barr et al. 2020; Hu et al. 2000; Liu et al. 2014).
Prior studies link exposure to HAP with cardiometabolic biomarkers and metabolomic profiles as measured in DBS (Clark et al. 2009), a method that overcomes collection, transportation, and storage limitations of venipuncture sampling. The ability to evaluate biomarkers of exposure in young children is an innovative tool for HAP research. Biomarkers of exposure (e.g., polycyclic aromatic hydrocarbons, volatile organic chemicals, and levoglucosan, a marker of wood smoke) are measured at multiple time points in all study participants; correlations will be assessed with measured 24-h pollutant concentrations. In the older adult woman, measurements include a suite of biomarkers of endothelial function, inflammation, oxidative stress, glycemic control/diabetes (HbA1c), a marker with specific relevance to lung cancer (P53 tumor-associated antigen antibodies), and enzyme induction (cytochrome P450). These biomarkers were chosen to inform mechanistic pathways and/or because of their capabilities to predict future disease risk. Biomarker discovery (metabolomics and miRNA evaluations) (Espín-Pérez et al. 2014) will be conducted using DBS collected in a subsample of children (n100/site) and older adult women (n100/site).
We are also conducting a substudy focused on cancer biomarkers, conducted in collaboration with the National Cancer Institute, among older adult women from the Peru and Guatemala research centers. These centers were chosen because early data indicated that they have a higher percentage of households with older women than the other sites. Among these older women, buccal cells, oral rinse, nasal turbinate, and peripheral (venous) blood are collected at baseline and a year later (in both intervention and control arms) and analyzed to provide additional data on the effects of biomass emission exposure on biomarkers related to cancer. Sample size for the cancer substudy was based on available budgetary resources and in the absence of any prior knowledge of the effect size for change in the prevalence of biomarkers.

Follow-Up Assessments

Each participating household is followed from enrollment until the index child reaches his/her first birthday. The data collection schedule during this period is summarized in Table 1.

Trial Management

The trial is led by a steering committee composed of the study’s multiple principal investigators (MPIs), the lead investigators of each of the international research centers, the directors of the biomarker core, one NIH project officer, one NIH scientific officer, and one representative of the Bill and Melinda Gates Foundation (BMGF). The steering committee obtains guidance on particular areas of expertise from five study cores (behavior and economics, exposure, clinical and imaging, biomarkers, and data management) and two working groups (pneumonia; anthropometry and nutrition). Day-to-day management of the trial is led by the trial coordinating center, based at Emory. The steering committee is supported by an external advisory committee and the pneumonia working group by an outside expert group. A data safety monitoring board (DSMB) appointed by the National Heart, Lung, and Blood Institute (NHLBI) is responsible for safeguarding the interests of study participants, assessing the safety and efficacy of study procedures, ensuring data quality, and monitoring the overall conduct of the study.

Data Management and Reporting

Field staff collect data on password-protected tablets, then upload the data (daily) to a secure REDCap™ (research electronic data capture) server [Health Insurance Portability and Accountability Act (HIPAA) and Federal Information Security Management Act (FISMA) compliant] hosted by Emory. Every day after uploads, the tablets are refreshed, and all data are deleted from the mobile device. Emory has a daily offsite backup on the REDCap™ server.
We also implemented algorithms for real-time data quality checking wherever possible, which triggers an error or warning message as the staff collect the data, even if they are offline. For example, if a staff member enters 20kg for a baby’s weight, the tablet shows a warning message that the value is more than the maximum allowed, and she/he can correct the entry.
The data management core (DMC), jointly with data managers in each country, is responsible for collecting survey and other field data on tablet-based platforms that incorporate GIS (geographical information systems) positioning and QR readers, allowing data collection teams to scan barcodes placed on LPG stoves, fuel tanks, exposure monitoring equipment, and biospecimen collection material. The DMC developed and maintains SOPs for all data management activities, which detail specific elements of data storage, entry, quality, transfer, processing, and security. The main study database is maintained in REDCap™ (for forms), LIMS™ (laboratory information management system for laboratory data), and TRICE™ (for clinical imaging). These systems are integrated using LabKey™.
Data reports with all relevant data and metadata for each of the analytic aims of the study are made available regularly to members of the study steering committee, the DSMB, and the relevant study investigators after consultation of exact data needs. However, until results are analyzed and reported by the DMC, investigators are blinded to study arm assignment except for the purpose of monitoring compliance with the intervention. Nevertheless, because of the nature of the intervention, fieldworkers are assumed to have knowledge of study arm assignment for data collected at the household level.
Final study data associated with HAPIN publications will be archived by the DMC to allow data sharing as per NIH requirements following a registration and review process by the DMC. The DMC will manage data requests and document those receiving data.

Protocol Standardization and Data Quality Control

Several activities have been undertaken to ensure standardization of field protocols within and between the four countries involved in this trial. These include the use of SOPs, step-by-step instruction sheets for each case report form (CRF), and a step-by-step visit schedule, all of which were translated into local languages and reviewed during in-person trainings at the field sites. In addition, during these in-person trainings, the representative of each core or work group directly observed field teams conduct practice visits and provided real-time feedback. Finally, field supervisors conducted weekly direct observations of enumerators and completion of a supervisor checklist that included key procedures to observe for each measurement, as determined by their respective core or work group.
With specific regard to CRFs, these were labeled as “M” for CRFs collecting information on the pregnant women, “C” for those collecting information on the child, “H” for those collecting information on the household, and “A” for those collecting information on older adult women. While the CRF content was translated into the local language, the CRF codes (e.g., “M10” for the mother’s demographic information) and question numbers were not changed to facilitate communication. In addition, bilingual color-coded CRFs were made available to all study staff in order to address any ambiguity on original questions and translations.
Data quality was monitored by each core or work group, as they had the required content knowledge to identify issues and propose appropriate solutions. For example, for the primary outcomes of birth weight and stunting (length-for-age z-score), the nutrition work group received the raw REDCap™ data download for the CRFs containing those measurements every 2 months from the central data management center at Emory. They then ran a standard report analysis code to generate a standard set of parameters established by a group of nutrition experts to assess data quality. These included summarizing missing data; descriptive statistics (mean, standard deviation, and minimum and maximum for continuous variables and prevalence of low birth weight and stunting, which are compared to DHS estimates); flagging outliers (defined as greater than ±3SD) and informing field teams so that they can identify potential errors in recording; tracking individual child growth trajectories and flagging children that lose weight or get shorter as they age, as these may be potential errors in recording and evaluating digit preference, the number of repeat measurements that are identical, and the time between repeat measurements as well as the time between birth and the birth weight measurement (should be within 24 h).

Data Analysis

Analysis of primary outcomes is described below; the analyses for secondary outcomes will be analogous. For aims 1 and 3, comparisons are made between the study arms to which participants were assigned, regardless of their actual adherence to the intended condition. Two-sided tests for each primary outcome will be performed separately at an αlevel of 0.0125 (using a conservative Bonferroni correction, i.e., α level of 0.05/4, the number of primary outcomes). The primary analyses will include indicator variables to account for the stratification of random sampling with different geographical areas in India and Peru, as well as individual sites.
Analysis for a continuous measure like birth weight, absent confounding because of the randomization, essentially is a t-test for mean birth weight between the two arms. For the two binary primary outcomes, pneumonia incidence over a 1-y follow-up period and stunting at 12 months of age, we will compare the two groups using their relative risk and accounting for randomization strata. BP analyses in older adult women will exclude women who self-report the use of BP medication at baseline (anticipated to be 5%). Repeated measurements of BP of the older adult woman will be analyzed using a random intercept model to examine differences in mean BP over the follow-up period, controlling for baseline BP. Women who begin BP medication during follow-up will be censored at that time in this analysis. Using longitudinal BP measurements, we will also consider whether the intervention is associated with the rate of change in BP changes over time in a sensitivity analysis by examining the interaction term between time and intervention.
In a supplemental analysis, women beginning BP medication after enrollment will be included via the method “j” described by Tobin et al. (2005), which is based on the assumption that treated women have BPs that are right censored in a normal distribution, with the censoring occurring above the level at which hypertensives are typically treated.
For exposure–response analyses, the general regression model is given by
where Yi is the primary outcome of interest, g() is the appropriate link function (identity for birth weight and BP, logistic for stunting and pneumonia), Xi is the continuous exposure of interest, and Zi is the vector of confounders. Potential confounders include, but are not limited to, age, sex, body mass index, socioeconomic status, tobacco use, secondhand smoke, physical activity, dietary intake, and season. These potential confounders have been identified as being associated with exposure or outcomes in prior studies. However, we expect that exact list of confounders to differ by outcome. For each outcome, we determine potential confounders using previous literature and directed acyclic graph methods. We will determine final models by evaluating change in effect estimates for meaningful changes. Exploration of models will consider the consideration of possible intermediate variables and effect modification. While we do not anticipate a great amount of missing data, we will evaluate whether the missing data patterns are differential between intervention and control groups, as well as evaluate the temporal pattern using Kaplan-Meier curves. For all analyses, we will assume missing data to be missing at random (MAR) a priori. The MAR mechanism specifies that the complete data distribution can be modeled using the observed data.
For repeated measurements of BP, an individual-level random intercept will be included, and baseline BP will be included in the model. We will consider PM2.5 mass to be our primary exposure measurement; we will also evaluate CO and BC in secondary analyses. For birth weight, a gestational exposure level will be obtained by averaging the two or three 24-h average measurements available during pregnancy. For birth weight, analyses will be restricted to full-term births, or we will take into account gestational age by using as the outcome z-scores for birth weight adjusted for gestational age [derived from INTERGROWTH tables (]. For stunting and pneumonia, we will consider average gestational and average first-year-of-life exposures (average of the 24-h measurements). For BP, we will consider time-varying exposures (24-h averages), available at the time of BP measurement.

Study Power

Study power for primary outcomes is given in Table 3, which summarizes the minimal detectable difference in mean (for birth weight and BP) and minimal detectable relative risk (for stunting and pneumonia) associated with an 80% power and a type I error rate α level of 0.0125, assuming a 10% attrition during follow-up. For BP, we show power calculations for 200 older adult women per arm, although we expect about 240 per arm, given that initial recruitment data suggest that there will be older adult women in about 15% of households. Table 3 shows we have good power to detect a smaller difference than has been previously reported in the literature.
Table 3 Minimum detectable effects, with 80% power and α=0.0125.
 Key parameter for estimating poweraSample size per armMinimal detectable effectPrevious studies’ estimate (95% CI)
Change in mean birth weight (g)σ21,6005489 (27, 204)b
Change in mean blood pressure (mm Hg)σ22002.583.7 (0.6, 8.1)c
Relative risk for stuntingp11,6000.810.79 (0.70, 0.89)d
Relative risk for HAPIN severe pneumoniap11,6000.640.67 (0.45, 0.98)e
Note: For a fuller description of the literature and the choice of three of these effect estimates (birth weight, blood pressure, infant pneumonia), see Steenland et al. (2018). CI, confidence interval; HAPIN, Household Air Pollution Intervention Network.
Our key population parameters were the variance of continuous measures (σ2) or the incidence rates among controls (p1) for relative risks. We took our key parameter estimates from the previous studies listed in the last column, with the exception of the control rate for pneumonia, which was estimated from our early data from the trial (Smith et al. 2011).
Thompson et al. (2011) is a substudy of a randomized trial of improved cookstoves in Guatemala (RESPIRE) in which those using stoves were compared to those not using stoves after adjusting for confounders.
McCracken et al. (2007) is a study of 120 older women (>38 years of age) from the same randomized trial (RESPIRE), which adjusted for confounders.
This comprehensive review article of the health effects of HAP exposure by Bruce et al. (2013) is based on two observational studies for moderate stunting [z-score<2standarddeviations(SD)].
Smith et al. (2011), from the RESPIRE randomized trial, provides the estimated relative risk (RR) of 0.67, the RR for clinician-diagnosed severe pneumonia for children under 18 months of age. Our estimated background rate of 0.09 for controls is based on our observed severe pneumonia rate in both arms (treatment and control) in the HAPIN study (using the HAPIN severe pneumonia definition), with 20% of the child person time observed, and assuming the RESPIRE RR of 0.67 for intervention vs. controls.


HAP is a leading cause of morbidity and mortality, especially among young children in LMICs (GBD 2018). Large-scale programs have succeeded in replacing traditional stoves with improved biomass stoves that reduce fuel consumption (Pope et al. 2017). However, deployment of these stoves has not generally achieved major reductions of fine PM, nor have they achieved meaningful improvements in health (Sambandam et al. 2015). Moreover, continued use of traditional stoves, along with improved solid fuel stoves (i.e., stacking), has reduced the potential contribution of such stoves to improved health. Systematic reviews have linked cleaner cooking fuels with improvements in health (Quansah et al. 2017; Thakur et al. 2018; WHO 2014c; Pope et al. 2017; Bruce et al. 2013), but evidence is still not strong, suffers from few well-designed studies, and varies by end point. Intervention studies of chimney stoves or improved solid biomass stoves have provided mixed evidence of health benefits (Sambandam et al. 2015; McCracken et al. 2007; Thompson et al. 2011; Hartinger et al. 2013; Tielsch et al. 2016; Johnson et al. 2013; Mortimer et al. 2017; Lee et al. 2019; Alexander et al. 2018). Evidence from observational studies is moderate for birth weight, pneumonia, and BP and still sparse for stunting. Recent trials have shown mixed results from clean cooking fuels, including ethanol and LPG (Alexander et al. 2018; Olopade et al. 2017; Alexander et al. 2017). Taken together, by addressing multiple outcomes at different ages and studying several LMIC contexts, the HAPIN trial is well positioned to fill critical scientific gaps with direct relevance to national policies.
This study addresses many of these gaps while also advancing knowledge in new areas. For example, the prevalence of low birth weight continues to be high in many LMICs, and preterm birth remains the leading cause of death among children under 5 years of age globally (Liu et al. 2015). While substantial literature exists on the effect of active tobacco smoking and ambient air pollution on low birth weight and preterm birth, few studies have examined the effect of these in LMICs or specifically the effect of HAP exposures on these outcomes (Thakur et al. 2018; Sambandam et al. 2015; Balakrishnan et al. 2013, 2018; Alexander et al. 2017; Liu et al. 2015; Amegah et al. 2014; Wylie et al. 2014). Only three studies have quantified the association between HAP and child developmental outcomes (Suter et al. 2018; Dix-Cooper et al. 2012; Munroe and Gauvain 2012), and only two quantified the association between HAP and stunting (Wylie et al. 2017; Kim et al. 2017), but none have done so using a longitudinal approach in children under 1 years of age. Our case definition of pneumonia is designed to balance the need for sensitivity and specificity (Goodman et al. 2019).
Similarly, the CVD burden from HAP is estimated from integrated exposure–response models (Burnett et al. 2014) that are not informed by direct measurements. Demonstrating potential improvements in cardiovascular biomarkers and BP through the LPG stove intervention can provide strategic information for the prevention of CVD itself. This is especially important for LMICs, where the health burden from noncommunicable diseases is increasing rapidly (WHO 2017). Finally, the inclusion of molecular biomarkers of exposure and early biological effect for a range of chronic health end points, including cancer, fills a critical need to establish biomarker-based approaches to assess long-term health impacts from HAP (IEA 2017; Caravedo et al. 2016).
While providing policy-relevant results, our research will incorporate important technical and training innovations, including expert clinical outcomes and imaging support for and training in pulmonary (point-of-care ultrasound for severe pneumonia), noninvasive markers of cardiovascular health (ultrasound assessment of endothelial function and CIMT), and obstetrics–gynecology (fetal growth). We will analyze targeted biomarkers of exposure, susceptibility, and effect, providing the most comprehensive analysis of relevant biomarkers in a HAP study to date. We will conduct discovery metabolomics and miRNA, mRNA, and DNA methylation analyses to evaluate perturbations occurring in response to HAP exposure to discover novel HAP biomarkers. We are implementing comprehensive stove-use monitoring, building on previous research (Ruiz-Mercado et al. 2012; Pillarisetti et al. 2014). These data, in conjunction with information gained in formative research on barriers to adoption and sustained, exclusive use of LPG stoves, will be used to maximize adherence and will help inform efforts to minimize stacking when stoves are delivered.
Our study will include a more comprehensive exposure assessment of PM2.5 and CO in RCTs than has been attempted to date, with multiple measurements, modeling of exposures, and assessments of stove stacking via SUMs. We expect to be able to estimate longer-term exposure averages, which can be used in exposure–response analyses. We anticipate a wide range of intervention to control and cross-cohort exposure contrasts across and within our four countries, making exposure–response analyses more powerful.
Moreover, our study will be the first HAP intervention trial to be conducted in multiple country settings using a common protocol, collecting a rich set of data with the potential to help define the generalizability criteria of the findings. While this efficacy study includes elements that are unlikely to be operationalized at scale (free fuel, purposefully selected settings with low to moderate ambient pollution, and frequent behavioral reinforcement to minimize stacking with traditional stoves), it will yield important comparable information about the most scalable clean fuel intervention across different world settings. As governments worldwide undertake efforts to reduce reliance on solid biomass fuels, the trial will provide data-supported evidence of the immediate and long-term health benefits that may be achieved by expanding access to and promoting the exclusive use of clean cooking fuels by those relying on solid biomass.


The HAPIN trial is funded by the NIH (cooperative agreement 1UM1HL134590) in collaboration with the Bill & Melinda Gates Foundation (OPP1131279). Participating NIH organizations include the NHLBI; National Institute of Environmental Health Sciences (NIEHS); National Cancer Institute (NCI); National Institute of Child Health and Human Development (NICHD); Fogarty International Center (FIC); and the NIH Common Fund.
The study protocol has been reviewed and approved by institutional review boards or ethics committees at Emory University (00089799), Johns Hopkins University (00007403), Sri Ramachandra Institute of Higher Education and Research (IEC-N1/16/JUL/54/49) and 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/11-2016) and Guatemalan Ministry of Health National Ethics Committee (11-2016), Asociacion Benefica PRISMA (CE2981.17), London School of Hygiene and Tropical Medicine (11664-5), Rwandan National Ethics Committee (No. 357/RNEC/2018), and Washington University in St. Louis (201611159). The study has been registered with (NCT02944682).
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: N.R. Cook, S. Hecht, C. Karr, J. Millum, N. Sathiakumar (Chair), P.K. Whelton, G.G. Weinmann (Executive Secretary). Program Coordination: G. Rodgers, Bill & Melinda Gates Foundation; C.L. Thompson, National Institute of Environmental Health Science; M.J. Parascandola, National Cancer Institute; D.M. Krotoski and M. Koso-Thomas, Eunice Kennedy Shriver National Institute of Child Health and Human Development; J.P. Rosenthal, Fogarty International Center; C.R. Nierras, NIH Office of Strategic Coordination Common Fund; K. Kavounis, D.Y. Kim, A. Punturieri, and B.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 U.S. National Institutes of Health or Department of Health and Human Services.

Article Notes

The authors declare they have no actual or potential competing financial interests.


Alexander DA, Northcross A, Karrison T, Morhasson-Bello O, Wilson N, Atalabi OM, et al. 2018. Pregnancy outcomes and ethanol cookstove intervention: a randomized-controlled trial in Ibadan, Nigeria. Environ Int 111:152–163.,
Alexander D, Northcross A, Wilson N, Dutta A, Pandya R, Ibigbami T, et al. 2017. Randomized controlled ethanol cookstove intervention and blood pressure in pregnant Nigerian women. Am J Respir Crit Care Med 195(12):1629–1639.,
Amegah AK, Quansah R, Jaakkola J. 2014. Household air pollution from solid fuel use and risk of adverse pregnancy outcomes: a systematic review and meta-analysis of the empirical evidence. PLoS One 9(12):e113920.,
Anenberg SC, Balakrishnan K, Jetter J, Masera O, Mehta S, Moss J, et al. 2013. Cleaner cooking solutions to achieve health, climate, and economic cobenefits. Environ Sci Technol 47(9):3944–3952.,
Balakrishnan K, Ghosh S, Ganguli B, Sambandam S, Bruce N, Barnes DF, et al. 2013. State and national household concentrations of PM2.5 from solid cookfuel use: results from measurements and modeling in India for estimation of the global burden of disease. Environ Health 12(1):77.,
Balakrishnan K, Ghosh S, Thangavel G, Sambandam S, Mukhopadhyay K, Puttaswamy N, et al. 2018. Exposures to fine particulate matter (PM2.5) and birthweight in a rural-urban, mother-child cohort in Tamil Nadu, India. Environ Res 161:524–531.,
Balakrishnan K, Sambandam S, Ramaswamy P, Mehta S, Smith KR. 2004. Exposure assessment for respirable particulates associated with household fuel use in rural districts of Andhra Pradesh, India. J Expo Sci Environ Epidemiol 14(suppl 1):S14–S25.,
Barr DB, Puttaswamy N, Jaacks LM, Steenland K, Rajkumar S, Gupton S, et al. 2020. Design and rationale of the biomarker center of the household air pollution intervention network (HAPIN) trial. Environ Health Perspect 128(4):047010,
Bonjour S, Adair-Rohani H, Wolf J, Bruce NG, Mehta S, Prüss-Ustün A, et al. 2013. Solid fuel use for household cooking: country and regional estimates for 1980–2010. Environ Health Perspect 121(7):784–790.,
Bruce NG, Dherani MK, Das JK, Balakrishnan K, Adair-Rohani H, Bhutta ZA, et al. 2013. Control of household air pollution for child survival: estimates for intervention impacts. BMC Public Health 13(suppl 3):S8.,
Bruce N, Pope D, Rehfuess E, Balakrishnan K, Adair-Rohani H, Dora C. 2015. WHO indoor air quality guidelines on household fuel combustion: strategy implications of new evidence on interventions and exposure–risk functions. Atmos Environ 106:451–457,
Burnett RT, Pope CA 3rd, Ezzati M, Olives C, Lim SS, Mehta S, et al. 2014. An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure. Environ Health Perspect 122(4):397–403.,
Caravedo MA, Herrera PM, Mongilardi N, de Ferrari A, Davila-Roman VG, Gilman RH, et al. 2016. Chronic exposure to biomass fuel smoke and markers of endothelial inflammation. Indoor Air 26(5):768–775.,
Clark ML, Peel JL, Balakrishnan K, Breysse PN, Chillrud SN, Naeher LP, et al. 2013. Health and household air pollution from solid fuel use: the need for improved exposure assessment. Environ Health Perspect 121(10):1120–1128.,
Clark ML, Peel JL, Burch JB, Nelson TL, Robinson MM, Conway S, et al. 2009. Impact of improved cookstoves on indoor air pollution and adverse health effects among Honduran women. Int J Environ Health Res 19(5):357–368.,
Corretti MC, Plotnick GD, Vogel RA. 1995. Technical aspects of evaluating brachial artery vasodilatation using high frequency ultrasound. Am J Physiol 286(4 Pt 2):H1397–H1408.,
Dix-Cooper L, Eskenazi B, Romero C, Balmes J, Smith KR. 2012. Neurodevelopmental performance among school age children in rural Guatemala is associated with prenatal and postnatal exposure to carbon monoxide, a marker for exposure to woodsmoke. Neurotoxicology 33(2):246–254.,
Espín-Pérez A, Krauskopf J, de Kok TM, Kleinjans JC. 2014. ‘OMICS-based’ biomarkers for environmental health studies. Curr Envir Health Rpt 1(4):353–362,
Fujisada H, Urai M, Uwasaki A. 2012. Technical methodology for ASTER global DEM. IEEE Trans Geosci Remote Sens 50(10):3725–3736,
GBD (Global Burden of Disease) Risk Factor Collaborators. 2018. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 392(10159):1923–1994.,
Goodman D, Crocker ME, Pervaiz F, McCollum ED, Steenland K, Simkovich SM, et al. 2019. Challenges in the diagnosis of paediatric pneumonia in intervention field trials: recommendations from a pneumonia field trial working group. Lancet Respir Med 7(12):1068–1083.,
Gordon SB, Bruce NG, Grigg J, Hibberd PL, Kurmi OP, Lam KB, et al. 2014. Respiratory risks from household air pollution in low and middle income countries. Lancet Respir Med 2(10):823–860.,
Government of the Republic of Guatemala. Results of the 2018 National Census Guatemala City, Guatemala 2019. Available from: [accessed 18 September 2019].
Hartinger SM, Commodore AA, Hattendorf J, Lanata CF, Gil AI, Verastegui H, et al. 2013. Chimney stoves modestly improved indoor air quality measurements compared with traditional open fire stoves: results from a small-scale intervention study in rural Peru. Indoor Air 23(4):342–352.,
Hollada J, Williams KN, Miele CH, Danz D, Harvey SA, Checkley W. 2017. Perceptions of improved biomass and liquefied petroleum gas stoves in Puno, Peru: implications for promoting sustained and exclusive adoption of clean cooking technologies. Int J Environ Res Public Health 14(2):182.,
Hu YA, Barr DB, Akland G, Melnyk L, Needham L, Pellizzari ED, et al. 2000. Collecting urine samples from young children using cotton gauze for pesticide studies. J Expo Anal Environ Epidemiol 10(6 Pt 2):703–709.,
IEA (International Energy Agency). 2017. Energy Access Outlook 2017: From Poverty to Prosperity. Paris, France: International Energy Agency,
Jaganath D, Miranda JJ, Gilman RH, Wise RA, Diette GB, et al. 2015. Prevalence of chronic obstructive pulmonary disease and variation in risk factors across four geographically diverse resource-limited settings in Peru. Respir Res 16(1):40.,
Johnson M, Piedrahita R, Garland C, Pillarisetti A, Sambandam S, et al. 2012. Exposures to PM2.5 associated with LPG stove and fuel interventions in four countries: pilot results from the HAPIN Trial. Presented at: ISEE Conference; August 26–30, 2018; Ottawa, Canada.
Johnson MA, Pilco V, Torres R, Joshi S, Shrestha RM, Yagnaraman M, et al. 2013. Impacts on household fuel consumption from biomass stove programs in India, Nepal, and Peru. Energ Sustain Dev 17(5):403–411,
Johnson M, Steenland K, Piedrahita R, Clark M, Pillarisetti A, Balakrishnan K, et al. 2020. Air pollutant exposure and stove use assessment methods for the Household Air Pollution Intervention Network (HAPIN) trial. Environ Health Perspect 128(4):047009,
Kim R, Mejia-Guevara I, Corsi DJ, Aguayo VM, Subramanian SV. 2017. Relative importance of 13 correlates of child stunting in South Asia: insights from nationally representative data from Afghanistan, Bangladesh, India, Nepal, and Pakistan. Soc Sci Med 187:144–154.,
Kirby MA, Nagel CL, Rosa G, Iyakaremye L, Zambrano LD, Clasen TF 2016. Faecal contamination of household drinking water in Rwanda: A national cross-sectional study. Sci Total Environ 571:426–434,
Lee AG, Kaali S, Quinn A, Delimini R, Burkart K, Opoku-Mensah J, et al. 2019. Prenatal household air pollution is associated with impaired infant lung function with sex-specific effects: evidence from GRAPHS, a cluster randomized cookstove intervention trial. Am J Respir Crit Care Med 199(6):738–746.,
Liao J, McCraken J, Piedrahita R, Thompson L, Millinedo E, Canuz E, et al. 2019. The use of Bluetooth low energy beacon systems to assess indirect personal exposure to household air pollution. J Expo Sci Environ Epidemiol. In press.,
Liu L, Oza S, Hogan D, Perin J, Rudan I, Lawn JE, et al. 2015. Global, regional, and national causes of child mortality in 2000–13, with projections to inform post-2015 priorities: an updated systematic analysis. Lancet 385(9966):430–440.,
Liu L, Xia T, Zhang X, Barr DB, Alamdar A, Zhang J, et al. 2014. Biomonitoring of infant exposure to phenolic endocrine disruptors using urine expressed from disposable gel diapers. Anal Bioanal Chem 406(20):5049–5054.,
McCoy DC, Sudfeld CR, Bellinger DC, Muhihi A, Ashery G, Weary TE, et al. 2017. Development and validation of an early childhood development scale for use in low-resourced settings. Popul Health Metr 15(1):3.,
McCracken JP, Smith KR, Díaz A, Mittleman MA, Schwartz J. 2007. Chimney stove intervention to reduce long-term wood smoke exposure lowers blood pressure among Guatemalan women. Environ Health Perspect 115(7):996–1001.,
Mortimer K, Ndamala CB, Naunje AW, Malava J, Katundu C, Weston W, et al. 2017. A cleaner burning biomass-fuelled cookstove intervention to prevent pneumonia in children under 5 years old in rural Malawi (the Cooking and Pneumonia Study): a cluster randomised controlled trial. Lancet 389(10065):167–175.,
Munroe RL, Gauvain M. 2012. Exposure to open-fire cooking and cognitive performance in children. Int J Environ Health Res 22(2):156–164.,
NISR (National Institute of Statistics of Rwanda). Rwanda DHS 2010.
NISR. Ministry of Finance and Economic Planning (MINECOFIN) [Rwanda], 2012. Rwanda Fourth Population and Housing Census.
Olopade CO, Frank E, Bartlett E, Alexander D, Dutta A, Ibigbami T, et al. 2017. Effect of a clean stove intervention on inflammatory biomarkers in pregnant women in Ibadan, Nigeria: a randomized controlled study. Environ Int 98:181–190.,
Pillarisetti A, Vaswani M, Jack D, Balakrishnan K, Bates MN, Arora NK, et al. 2014. Patterns of stove usage after introduction of an advanced cookstove: the long-term application of household sensors. Environ Sci Technol 48(24):14525–14533.,
Pollard SL, Williams DL, Breysse PN, Baron PA, Grajeda LM, et al. 2014. A cross-sectional study of determinants of indoor environmental exposures in households with and without chronic exposure to biomass fuel smoke. Environ Health 13(1):21,
Pope D, Bruce N, Dherani M, Jagoe K, Rehfuess E. 2017. Real-life effectiveness of ‘improved’ stoves and clean fuels in reducing PM2.5 and CO: systematic review and meta-analysis. Environ Int 101:7–18.,
Quansah R, Semple S, Ochieng CA, Juvekar S, Armah FA, Luginaah I, et al. 2017. Effectiveness of interventions to reduce household air pollution and/or improve health in homes using solid fuel in low-and-middle income countries: a systematic review and meta-analysis. Environ Int 103:73–90.,
Quinn AK, Williams K, Thompson LM, Rosa G, Díaz-Artiga A, Thangavel G, et al. 2019. Compensating control participants when the intervention is of significant value: experience in Guatemala, India, Peru and Rwanda. BMJ Glob Health 4(4):e001567.,
Rosenthal J, Balakrishnan K, Bruce N, Chambers D, Graham J, Jack D, et al. 2017. Implementation science to accelerate clean cooking for public health. Environ Health Perspect 125(1):A3–A7. Jan 1.,
Ruiz-Mercado I, Canuz E, Smith KR. 2012. Temperature dataloggers as stove use monitors (SUMs): field methods and signal analysis. Biomass Bioenergy 47:459–468.,
Sambandam S, Balakrishnan K, Ghosh S, Sadasivam A, Madhav S, Ramasamy R, et al. 2015. Can currently available advanced combustion biomass cook-stoves provide health relevant exposure reductions? Results from initial assessment of select commercial models in India. Ecohealth 12(1):25–41.,
Simkovich S, Underhill L, Kirby M, Goodman D, Crocker M, Hossen S, et al. Design and conduct of facility-based surveillance for severe childhood pneumonia in the Household Air Pollution Intervention Network (HAPIN) trial. ERJ Open Research. In press.
Smith KR, Bruce N, Balakrishnan K, Adair-Rohani H, Balmes J, Chafe Z, et al. 2014. Millions dead: how do we know and what does it mean? Methods used in the comparative risk assessment of household air pollution. Annu Rev Public Health 35:185–206.,
Smith KR, McCracken JP, Weber MW, Hubbard A, Jenny A, Thompson LM, et al. 2011. Effect of reduction in household air pollution on childhood pneumonia in Guatemala (RESPIRE): a randomised controlled trial. Lancet 378(9804):1717–1726.,
Steenland K, Pillarisetti A, Kirby M, Peel J, Clark M, Checkley W, et al. 2018. Modeling the potential health benefits of lower household air pollution after a hypothetical liquified petroleum gas (LPG) cookstove intervention. Environ Int 111:71–79.,
Stein JH, Korcarz CE, Hurst RT, Lonn E, Kendall CB, Mohler ER, et al. 2008. Use of carotid ultrasound to identify subclinical vascular disease and evaluate cardiovascular disease risk: a consensus statement from the American Society of Echocardiography Carotid Intima-Media Thickness Task Force. Endorsed by the Society for Vascular Medicine. J Am Soc Echocardiogr 21(2):93–111.,
Suter M, Karr C, John-Stewart G, Gómez L, Moraa H, Nyatika D, et al. 2018. Implications of combined exposure to household air pollution and HIV on neurocognition in children. Int J Environ Res Public Health 15(1):163.,
Thakur M, Nuyts PAW, Boudewijns EA, Flores Kim J, Faber T, Babu GR, et al. 2018. Impact of improved cookstoves on women’s and child health in low and middle income countries: a systematic review and meta-analysis. Thorax 73(11):1026–1040.,
Thompson LM, Bruce N, Eskenazi B, Diaz A, Pope D, Smith KR. 2011. Impact of reduced maternal exposures to wood smoke from an introduced chimney stove on newborn birth weight in rural Guatemala. Environ Health Perspect 119(10):1489–1494.,
Tielsch JM, Katz J, Khatry SK, Shrestha L, Breysse P, Zeger S, et al. 2016. Effect of an improved biomass stove on acute lower respiratory infections in young children in rural Nepal: a cluster-randomised, step-wedge trial. Lancet Glob Health 4:S19,
Tobin MD, Sheehan NA, Scurrah KJ, Burton PR. 2005. Adjusting for treatment effects in studies of quantitative traits: antihypertensive therapy and systolic blood pressure. Stat Med 24(19):2911–2935.,
WHO (World Health Organization). 2014a. Integrated Management of Childhood Illness Chart Booklet. Geneva, Switzerland: World Health Organization.
WHO. 2014b. Revised WHO Classification and Treatment of Pneumonia in Children at Health Facilities: Evidence Summaries. Geneva, Switzerland: World Health Organization.
WHO. 2014c. WHO guidelines for indoor air quality: household fuel combustion. Geneva, Switzerland: World Health Organization.
WHO. 2017. Noncommunicable Diseases Progress Monitor 2017. License: CC BY-NC-SA 3.0 IGO. Geneva, Switzerland: World Health Organization.
Wylie BJ, Coull BA, Hamer DH, Singh MP, Jack D, Yeboah-Antwi K, et al. 2014. Impact of biomass fuels on pregnancy outcomes in central East India. Environ Health 13(1):1.,
Wylie BJ, Kishashu Y, Matechi E, Zhou Z, Coull B, Abioye AI, et al. 2017. Maternal exposure to carbon monoxide and fine particulate matter during pregnancy in an urban Tanzanian cohort. Indoor Air 27(1):136–146.,
Xu H, Li Y, Guinot B, Wang J, He K, Ho KF, et al. 2018. Personal exposure of PM2.5 emitted from solid fuels combustion for household heating and cooking in rural Guanzhong Plain, northwestern China. Atmos Environ 185:196–206,
Zuk M, Rojas L, Blanco S, Serrano P, Cruz J, Angeles F, et al. 2007. The impact of improved wood-burning stoves on fine particulate matter concentrations in rural Mexican homes. J Expo Sci Environ Epidemiol 17(3):224–232.,

Information & Authors


Published In

Environmental Health Perspectives
Volume 128Issue 4April 2020
PubMed: 32347766


Received: 21 October 2019
Revision received: 7 February 2020
Accepted: 11 February 2020
Published online: 29 April 2020



Thomas Clasen
Department of Environmental Health, Rollins School of Public Health, Emory University, Atlanta, Georgia, USA
William Checkley
Division of Pulmonary and Critical Care, School of Medicine, Johns Hopkins University, Baltimore, Maryland, USA
Jennifer L. Peel
Department of Environmental and Radiological Health Sciences, Colorado State University, Fort Collins, Colorado, USA
Kalpana Balakrishnan
Department of Environmental Health Engineering, ICMR Center for Advanced Research on Air Quality, Climate and Health, Sri Ramachandra Institute for Higher Education and Research (Deemed University), Chennai, Tamil Nadu, India
John P. McCracken
Center for Health Studies, Universidad del Valle de Guatemala, Guatemala City, Guatemala
Ghislaine Rosa
Department of Disease Control, Faculty of Infections and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
Lisa M. Thompson
Nell Hodgson Woodruff School of Nursing, Emory University, Atlanta, Georgia, USA
Dana Boyd Barr
Department of Environmental Health, Rollins School of Public Health, Emory University, Atlanta, Georgia, USA
Maggie L. Clark
Department of Environmental and Radiological Health Sciences, Colorado State University, Fort Collins, Colorado, USA
Michael A. Johnson
Berkeley Air Monitoring Group, Berkeley, California, USA
Lance A. Waller
Department of Biostatistics and Bioinformatics, Rollins School of Public Health, Emory University, Atlanta, Georgia, USA
Lindsay M. Jaacks
Department of Global Health and Population, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, USA
Kyle Steenland
Department of Environmental Health, Rollins School of Public Health, Emory University, Atlanta, Georgia, USA
J. Jaime Miranda
CRONICAS Center of Excellence in Chronic Diseases, Universidad Peruana Cayetano Heredia, Lima, Peru
Howard H. Chang
Department of Biostatistics and Bioinformatics, Rollins School of Public Health, Emory University, Atlanta, Georgia, USA
Dong-Yun Kim
Office of Biostatistics Research, National Heart Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland, USA
Eric D. McCollum
Eudowood Division of Respiratory Sciences, Department of Pediatrics, School of Medicine, Johns Hopkins University, Baltimore, Maryland, USA
Victor G. Davila-Roman
Cardiovascular Imaging and Clinical Research Core Laboratory, Cardiovascular Division, Washington University School of Medicine, St. Louis, Missouri, USA
Aris Papageorghiou
Nuffield Department of Women’s and Reproductive Health, University of Oxford, Oxford, UK
Joshua P. Rosenthal
Division of Epidemiology and Population Studies, Fogarty International Center, National Institutes of Health, Bethesda, Maryland, USA


Address correspondence to Thomas Clasen, Department of Environmental Health, Rollins School of Public Health, Emory University, 1518 Clifton Rd. NE, Atlanta, GA 30322 USA. Telephone: (404) 727-3480. Email: [email protected]

Metrics & Citations


About Article Metrics


Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click DOWNLOAD.

Cited by

  • Household Air Pollution Interventions to Improve Health in Low- and Middle-Income Countries: An Official American Thoracic Society Research Statement, American Journal of Respiratory and Critical Care Medicine, 10.1164/rccm.202402-0398ST, 209, 8, (909-927), (2024).
  • Impact of indoor environment on children’s pulmonary health, Expert Review of Respiratory Medicine, 10.1080/17476348.2024.2307561, 17, 12, (1249-1259), (2024).
  • Liquefied Petroleum Gas or Biomass Cooking and Severe Infant Pneumonia, New England Journal of Medicine, 10.1056/NEJMoa2305681, 390, 1, (32-43), (2024).
  • Effects of Cooking with Liquefied Petroleum Gas or Biomass on Stunting in Infants, New England Journal of Medicine, 10.1056/NEJMoa2302687, 390, 1, (44-54), (2024).
  • Effects of a liquefied petroleum gas stove intervention on stillbirth, congenital anomalies and neonatal mortality: A multi-country household air pollution intervention network trial, Environmental Pollution, 10.1016/j.envpol.2024.123414, 345, (123414), (2024).
  • Estimating personal exposures to household air pollution and plastic garbage burning among adolescent girls in Jalapa, Guatemala, Chemosphere, 10.1016/j.chemosphere.2023.140705, 348, (140705), (2024).
  • Human Milk Composition Is Associated with Maternal Body Mass Index in a Cross-Sectional, Untargeted Metabolomics Analysis of Human Milk from Guatemalan Mothers, Current Developments in Nutrition, 10.1016/j.cdnut.2024.102144, (102144), (2024).
  • Cooking with liquefied petroleum gas or biomass and fetal growth outcomes: a multi-country randomised controlled trial, The Lancet Global Health, 10.1016/S2214-109X(24)00033-0, 12, 5, (e815-e825), (2024).
  • Lung Function Impairment Among Women Cooking with Biomass Fuels Compared to LPG Fuel in Rural Tamil Nadu, Biomedical and Pharmacology Journal, 10.13005/bpj/2728, 16, 3, (1505-1511), (2023).
  • Heat Exposure among Adult Women in Rural Tamil Nadu, India, Environmental Science & Technology, 10.1021/acs.est.3c03461, 58, 1, (315-322), (2023).
  • See more

View Options

View options


View PDF

Get Access

Restore your content access

Enter your email address to restore your content access:

Note: This functionality works only for purchases done as a guest. If you already have an account, log in to access the content to which you are entitled.







Copy the content Link

Share on social media