Skip to main content

The effects of indoor environmental exposures on pediatric asthma: a discrete event simulation model



In the United States, asthma is the most common chronic disease of childhood across all socioeconomic classes and is the most frequent cause of hospitalization among children. Asthma exacerbations have been associated with exposure to residential indoor environmental stressors such as allergens and air pollutants as well as numerous additional factors. Simulation modeling is a valuable tool that can be used to evaluate interventions for complex multifactorial diseases such as asthma but in spite of its flexibility and applicability, modeling applications in either environmental exposures or asthma have been limited to date.


We designed a discrete event simulation model to study the effect of environmental factors on asthma exacerbations in school-age children living in low-income multi-family housing. Model outcomes include asthma symptoms, medication use, hospitalizations, and emergency room visits. Environmental factors were linked to percent predicted forced expiratory volume in 1 second (FEV1%), which in turn was linked to risk equations for each outcome. Exposures affecting FEV1% included indoor and outdoor sources of NO2 and PM2.5, cockroach allergen, and dampness as a proxy for mold.


Model design parameters and equations are described in detail. We evaluated the model by simulating 50,000 children over 10 years and showed that pollutant concentrations and health outcome rates are comparable to values reported in the literature. In an application example, we simulated what would happen if the kitchen and bathroom exhaust fans were improved for the entire cohort, and showed reductions in pollutant concentrations and healthcare utilization rates.


We describe the design and evaluation of a discrete event simulation model of pediatric asthma for children living in low-income multi-family housing. Our model simulates the effect of environmental factors (combustion pollutants and allergens), medication compliance, seasonality, and medical history on asthma outcomes (symptom-days, medication use, hospitalizations, and emergency room visits). The model can be used to evaluate building interventions and green building construction practices on pollutant concentrations, energy savings, and asthma healthcare utilization costs, and demonstrates the value of a simulation approach for studying complex diseases such as asthma.

Peer Review reports


In the US, asthma is among the most common chronic diseases of childhood across all socioeconomic classes and is the most frequent cause of hospitalization among children after birth [1]. Higher asthma prevalence has been documented in low-income inner-city children in many cities [25]. A number of studies have documented relationships between asthma exacerbation and exposure to indoor environmental stressors found in residential settings, such as allergens (e.g., dust mites, cockroach allergens), air pollutants (e.g., ozone (O3), nitrogen dioxide (NO2), fine particulate matter (PM2.5)), and environmental tobacco smoke (ETS) [6, 7]. Asthma exacerbations and related events are clearly influenced by numerous additional factors, including but not limited to access to health care, medication compliance, and rhinovirus and other infectious agents [810].

Because of the multi-factorial nature of asthma exacerbations, it can be challenging to design optimal intervention strategies. Studies have demonstrated public health benefits of residential interventions such as using community health workers to provide asthma education [11, 12], conducting integrated pest management [13, 14], or multi-factorial indoor interventions [15]. It is difficult to make generalized conclusions about intervention efficacy, because of significant differences in context and risk factors for various populations, as well as the possibility that intensive interventions using community health workers have significant social support components that provide benefits beyond improvements in the physical environment [16]. A recent systematic review by the Centers for Disease Control and Prevention (CDC) found that multi-trigger, multi-component, home based environmental interventions for children were effective at improving asthma quality of life and productivity, but did not specify which components were essential elements [17].

Indoor environmental interventions are also prone to complex tradeoffs among pollutants, as interventions that influence ventilation can have opposing effects on indoor and outdoor sources, and interventions that address pests can lead to increased pesticide exposures if not designed and implemented appropriately [18]. For example, higher NO2 concentrations tend to be measured indoors in homes with gas stoves [19, 20], indicating that improved venting of gas stoves or increased ventilation in general will reduce NO2 concentrations. However, increasing general ventilation increases indoor NO2 concentrations from outdoor sources, particularly in urban settings with high traffic.

Simulation modeling can be a valuable tool for evaluating intervention strategies across a range of outcomes, especially in the presence of significant tradeoffs. In this context, simulation modeling refers to a systems science approach involving modeling of a complex system that evolves over time given changes in state variables that occur at defined points in time [21]. Such models have been used for a number of health policy analyses such as evaluating alternative interventions for schizophrenia [22], malaria [23] diabetes [24], and breast cancer [25, 26]. In spite of the flexibility of a simulation modeling approach and its applicability to a multi-exposure setting for a complex multi-factorial disease such as asthma, modeling applications in either environmental exposures or asthma have been limited to date.

Prior modeling work in asthma has been limited to Markov state-transition models. These models have focused on either the adult population or interventions involving medication adherence alone [2729]. While these studies provided valuable insight about medication cost-effectiveness, Markov state-transition models are limited given their inability to track individuals, take account of multiple individual attributes, and capture interactions and non-linear effects. By definition, Markov models are memory-less, and the probability of transitioning between states does not depend on prior history within the simulation [30]. In contrast, asthma has been characterized as a classic example of a dynamic and non-linear disease with numerous influential factors [31], for which alternative modeling approaches may be more informative.

In this article, we provide the conceptual design and parameterization for a discrete event simulation model of pediatric asthma, focusing on the efficacy of indoor environmental interventions in low-income multi-family dwellings. These residences are characterized by having smaller living spaces and high occupant density, as well as an increased prevalence of peeling paint, indoor water leaks, and structural deficiencies that provide points of entry for cockroaches and other pests [18]. The model begins with a baseline population of high-risk children – characterized by demographic, residential, and behavioral factors. For each household, we modeled indoor environmental concentrations of multiple contaminants for fifty thousand simulated households across a variety of conditions. Based on literature syntheses, these concentrations are used to predict the percent predicted forced expiratory volume in one second (FEV1%), which in turn determines the probability of asthma exacerbations and health care utilization. To establish model credibility, we compared outputs with the published literature. As model evaluation is an iterative process, we will continue this as new data come to light. Finally, we conduct a simulation experiment with a simple hypothetical intervention to illustrate the insights available through this model approach, and discuss the strengths and weaknesses of our model structure.



We designed the discrete event simulation to model children’s exposure over time to environmental factors and their risk of having an asthma event on a daily basis. We modeled health outcomes by first estimating children’s exposure to indoor combustion pollutants, allergens, and other risk factors and using these factors to predict daily changes in lung function (FEV1%), as done elsewhere in the context of policy models for asthma medication [2729]. Children were assigned a baseline FEV1% which was affected by daily random variation as well as external conditions including pollutant exposures. The daily risk of asthma outcomes (symptoms, oral steroid use, emergency room (ER) visits, hospitalizations) was calculated based on changes in FEV1% (Figure 1). The model parameters, described and justified in detail below, are summarized in Additional file 1.

Figure 1
figure 1

Schematic of discrete event simulation model showing an apartment within a multi-family building from CONTAM software [[76]], and relationships between model inputs.

The model, built in R (R 2.12.1, The R foundation of Statistical Computing), generates an ensemble or cohort of children and their associated households that have characteristics typical of Boston public housing. We used the stream package for random number generation. We used the Scientific Computer Facilities at Boston University to simulate fifty thousand children individually over 10 years, or until they turned 18 years of age. Individual results were aggregated to obtain results for the cohort.

Study population and housing characteristics

The target simulated population was children living in low-income multi-family housing consistent with public housing residents, a population with high asthma prevalence and severity [3234] – a community survey in two Boston Housing Authority developments found pediatric asthma prevalence of 22% [35], compared to 8.5% in the general population [36]. Inputs describing demographic and housing characteristics were drawn from studies in Boston public housing and similar settings and are presented in Table 1. In the case of age, gender, apartment level, and type of smoker (light versus heavy), values were uniformly distributed across the simulated population. The remaining characteristics obtained from studies in Boston public housing were either derived from a community-based survey administered to 53 households in one development [32], or from an asthma intervention study enrolling 78 children from 61 households in 3 developments [34]. We simulated fifty thousand children to ensure an ability to detect changes in relatively infrequent asthma events (like hospitalizations) associated with changes in a single environmental risk factor.

Table 1 Baseline occupant and household characteristics of a simulated population of low-income asthmatic children

While we lacked sufficient symptom and severity data to apply the National Heart, Blood and Lung Institute (NHLBI) classification guidelines for managing asthma [8], we wanted to characterize asthma severity to determine the medications that each child would typically be prescribed. We therefore used the FEV1% cutoffs that correspond with severity classification for persistent asthmatics (> 80% for mild, 60-80% for moderate, and < 60% for severe), and used these to determine the prescribed medications. Intermittent asthmatics were excluded from the simulation because of the limited environmental literature on intermittent asthmatics – inclusion criteria for most epidemiology and environmental studies require having persistent asthma. Because of our focus on persistent asthmatics above the age of 5, many of whom are likely sensitized to one or more allergens, we assumed that asthma remained throughout childhood (i.e., no remission). Time-varying characteristics included age, indoor and outdoor NO2 and PM2.5, outdoor temperature, indoor and outdoor relative humidity (RH), cockroach allergen, mold exposure (defined by “damp” housing), daily random variation in baseline FEV1%, and changes in FEV1% due to all risk factors.

Lung function (FEV1% predicted)

Lung function was modeled using forced expiratory volume at 1 second expiration (FEV1), a measure that is obtained by spirometry. FEV1 is generally thought to reflect airway obstruction, as observed in asthma and chronic obstructive pulmonary disease (COPD). The values of FEV1 are generally translated into a percentage of predicted FEV1 based on age, race, and gender (FEV1%), and are thought to be normal if ranging from 80 to 120% of predicted, given variations in effort by the patient or administration of the test itself. There is sufficient evidence of the relationship between FEV1% and both environmental exposures and asthma outcomes, and FEV1% is a key predictor of asthma morbidity within evidence-based national guidelines for asthma care [8]. Further, this metric has previously been used in simulation modeling analyses of asthma interventions [2729]. Other measures of lung function, including the ratio of FEV1 to forced vital capacity (FVC), peak expiratory flow (PEF), and forced expiratory flow from the 25th percentile to the 75th percentile of FVC (FEF25-75%), either had an insufficient literature examining the relationship with key environmental exposures and/or were not robust measures of asthma severity. Although some investigators have shown weak correlations between FEV1% and asthma events [37, 38], we considered it to be the most robust and interpretable predictor available.

The distribution of baseline FEV1% was derived using spirometry data for children with persistent asthma who participated in the Healthy Public Housing Initiative (HPHI), an intervention study focused on pediatric asthmatics living in multi-family public housing developments in Boston [13, 14, 19, 39]. For FEV1% measurements taken prior to HPHI environmental interventions, mean baseline FEV1% was 88.4% (SD = 11.6%). However, this distribution reflects the likely cumulative effect of several environmental and socioeconomic factors present within the HPHI cohort, including indoor environmental influences. Thus, if we assumed that this distribution corresponded with a ‘no-exposure’ scenario and subsequently observed further reductions in FEV1% given the simulated exposures in our study, we would be systematically biased in our characterization of asthma severity. To address this issue, we reset the mean baseline value but kept the relative FEV1% distribution. We first ran our simulation model with the HPHI FEV1% distribution and determined the mean decrease in FEV1% attributable to all environmental exposures in comparison with a no-exposure scenario (decrease = 19.8%). We then approximated the no-exposure scenario by increasing all baseline FEV1% values by 19.8%, with a maximum allowed FEV1% of 120% to avoid findings which might be inconsistent with normal pulmonary function test results. We subsequently simulate the influence of multiple environmental exposures and medication use on FEV1% and evaluate the resulting distribution, as described in more detail below.

To characterize day-to-day variation in FEV1% within individuals, we constructed a distribution of daily variability using data collected from the HPHI study. We averaged the variability of individual repeated spirometry measures for 49 children pre-intervention over multiple testing sessions in the year, with each session including at least one week of twice-daily spirometry. A total of 110 testing sessions were recorded, with an average of 10 spirometry measurements included per child per testing session, and each pre-intervention spirometry session including at least one week of twice-daily spirometry. The mean variation between measurements within sessions was 10%, which includes random variation and variation due to environmental effects. Our objective is to include only random variation, and there are no formal analytical approaches available to determine how much of the 10% was random versus due to environmental effects. For the simulation, we assumed daily random variation of 5%, half of that observed in HPHI. Thus, daily FEV1% was assigned by drawing a value from a normal distribution, with a mean value of the baseline FEV1%, and a standard deviation of 5%.

Long term decline in FEV1% was determined based on coefficients presented by O’Byrne et al. (Table 2, [40]), where children began the simulation with a typical yearly decline in FEV1% (−0.8% for 5–10 year olds and 0% for 11–17 year olds). These values were relevant to both those who took corticosteroids and those who did not. In theory, all persistent asthmatics should be prescribed corticosteroids or other controller medications according to NHLBI clinical guidelines, but studies have found gaps related to health care access, compliance, miscommunication between providers and patients, misunderstanding of asthma as a chronic vs. acute disease, and other factors. For our study, we therefore use outputs from O’Byrne et al. to represent asthmatics with appropriately prescribed medications (termed “compliant”) and those without (termed “non-compliant”). Long term decline coefficients were modified by the occurrence of a severe asthma-related event (SARE), based on the definition of SARE and coefficients presented by O’Byrne et al. (Table 2, [40]). In our simulation a SARE was defined as a hospitalization or ER visit. In 5–10 year olds, when a SARE occurred, lung function decreased by −2.1%/year in non-compliant children and remained at −0.8%/year in compliant children. In 11–17 year olds, when a SARE occurred, lung function decreased by −1%/year in non-compliant children and −0.3%/year in compliant children. These yearly rates were converted to daily rates and applied to the daily FEV1% on the day after the SARE occurred. Children were allowed to return to pre-SARE FEV1% decline rates if they remained SARE-free for 3 years, the time period used in the O’Byrne study.

Table 2 Asthma medication prescription and usage as a function of FEV1%, based on NHLBI guidelines given asthma severity classification[8]

Following this approach, mean FEV1% at the end of a 10-year simulation for 50,000 asthmatic children was 84.9% (SD = 11.4%, range: 51-118%), indicating that we have reasonably represented the underlying distribution of spirometry values from HPHI, with deviations related to modeled variability in exposures not explicitly considered in the original distribution.

Asthma medication and compliance

Medication usage is both a key outcome variable influenced by the frequency of asthma events and a potential modifier of the effect of environmental exposures on health outcomes. The probability of a child having a prescription for controller medications was determined from HPHI data. A prior publication [41] gave the percentage of children in each severity category (mild, moderate, or severe persistent) who reported being prescribed at least one controller medication. However, as our model does not have sufficient information to formally classify severity, we reanalyzed raw data from HPHI to estimate the relationship between FEV1% and the probability of being prescribed a controller medication. We used SAS (Proc Logit, version 9.1, SAS Institute Inc., Cary, NC) to calculate the odds of being prescribed asthma medication, and converted the odds ratio to a probability estimate. The resulting probability equation was:

P med = exp 2 . 228 - 2 .854*FEV1% 1 + exp 2 . 228 - 2 .854*FEV1% ,

where P med is the probability of being prescribed and reporting taking a controller medication, and FEV1% is the baseline lung function value. Of note, we are interested in both children who are not prescribed controller medications and those who are prescribed the medications but are non-adherent. Lacking data to quantify each of these components, we assume that the self-report from HPHI represents a combination of the two. Therefore, if the caregiver reported that a child was prescribed and was taking a controller medication, it was assumed the child was “compliant” with the medications and its protective effects. A child who should have, but was not prescribed a controller, was “non-compliant” according to clinical guidelines. This may have resulted in overestimating the use of controller medication, given the literature documenting failure to adhere to medication use even with a prescription [42], but reliable compliance data were not available from the HPHI study population and literature values would not be well aligned with the available data.

Children were evaluated at the end of each simulated year to determine changes in their asthma medication prescription, approximating adjustments that would happen during a yearly physical exam. At every year anniversary we compared each child’s average FEV1% during the past year to his or her FEV1% at the beginning of the year. Although severity classification and changes in medication are based on many components beyond FEV1%, we simplified this step by using the standard ranges of FEV1% associated with each severity classification [8], and changing their asthma medication prescription accordingly (Table 2).

Linking indoor pollutant exposures and FEV1%

The model includes four contaminants that can potentially affect a child’s lung function and asthma status. We modeled two combustion pollutants – NO2 from gas stoves or outdoors, and PM2.5 from cooking, smoking, or outdoors – and two allergens (cockroach and dampness as a proxy for mold). Other common pollutants associated with asthma exacerbations such as ozone, mouse, cat, dog, and dust mite allergen were not included because we either lacked a critical mass of literature linking the exposure with FEV1%, or an ability to readily model indoor concentrations.

For each of the contaminants included in the study, we conducted a systematic literature review using PubMed. We initially utilized broad search criteria described below, and subsequently reviewed each article for its interpretability and applicability to our discrete event simulation model. For example, for NO2 and PM2.5, studies needed to utilize exposure metrics that could be constructed from our indoor air quality modeling (described below), and needed to include asthmatic children who were reasonably representative of our simulated population. For cockroach allergens, articles were included if they measured allergen levels in the household dust of study participants, given the difficulty in linking physical environmental measurements to measures such as IgE expression. For studies of mold and/or moisture, articles were included if they included exposure characterization that would be feasible within our exposure modeling, which focused on moisture/mold characterization given RH outputs.

For PM2.5 and NO2, we conducted a joint literature review given a desire to avoid double-counting of effects (i.e., with single-pollutant models not controlling for other pollutants). PubMed search terms included [“lung function” AND(PM2.5 OR particulate matter OR nitrogen dioxide OR NO2) AND(asthma OR children)], [("childhood asthma" OR “pediatric asthma”) AND("particulate matter" OR “nitrogen dioxide”)], and [asthma AND("particulate matter" or “nitrogen dioxide”) AND FEV]. These searches collectively yielded 413 publications, though with multiple duplicates.

Of these publications, 17 were sufficiently relevant to our study to merit more in-depth evaluation [4359]. Studies were excluded for a variety of reasons, including a lack of focus on short-term changes in lung function due to short-term changes in air pollution. We eliminated 12 of these studies given issues such as the use of FEV1 rather than FEV1% as the outcome measure without sufficient data to allow for conversion [45, 46, 50, 5355, 5759], season-specific characterization [56], or missing quantification of values necessary for our study [47, 51]. For the remaining five studies, two [43, 44] were panel studies in Los Angeles, two [49, 60] were studies of the same panel of children in Windsor, Ontario, and one [52] was based on children within the Inner City Asthma Study. We did not formally pool these available estimates because of the small number of applicable independent studies, as well as the significant heterogeneity in exposure metrics and statistical approaches. Instead, we chose the one study by O’Connor et al. [52] that provided multi-pollutant estimates among the full study population of children with exposure metrics that were available from our indoor air quality simulations. To ensure that this choice of study did not significantly bias our model, we conducted inverse-variance weighted pooling of the other four studies using the most comparable estimates available, and compared the results to those derived strictly from O’Connor et al. In spite of some limitations in the four studies and key differences in methods and assumptions, the pooled concentration-response unction was reasonably similar to the values from O’Connor et al. (20% higher for PM2.5 and 34% lower for NO2), providing reassurance that our choice was reasonable.

More specifically, we used coefficients reported in Table III of O’Connor et al. to estimate the effect of NO2 and PM2.5 on FEV1%. They reported a change in FEV1% at the 90th percentile of NO2 and PM2.5 concentration relative to the 10th percentile for exposure to a 5-day average outdoor pollutant concentration, adjusted for site, month, site-by-month interaction, temperature, intervention group and multiple pollutants (PM2.5, NO2, and O3). We converted the coefficients to a percent change in FEV1% per unit increase in pollutant, and divided by an infiltration factor to convert outdoor concentrations to equivalent indoor pollutant concentrations. This method was applied elsewhere to align epidemiologic effect estimates based on outdoor measurements with an indoor exposure model [61]. The infiltration factor for NO2 was 0.58 (average of factors reported by [20, 6264]), and 0.72 for PM2.5 (average of factors reported by [6466]). Coefficients were a −0.093% (SE = 0.030) change in FEV1% per 1 ppb increase in NO2, and −0.077% (SE = 0.032) change in FEV1% per 1 μg/m3 increase in PM2.5.

For cockroach allergen exposure, we found 36 results using the PubMed search term [“pediatric asthma AND cockroach”] and 75 results using the search term [“childhood asthma AND cockroach”], but none of these articles met our inclusion criteria. Examination of other articles known to the authors did not provide adequate information for our application, largely due to a lack of quantitative exposure measures in residential dust, which would be necessary to model the marginal benefits of intervention strategies. One option was to use studies of IgE expression in children [67]. However, this would require us to dynamically model the association between allergen concentrations in dust and serum IgE, which would be quite complex and uncertain. The other option was to use studies that used a more relevant exposure metric but focused on adults. We opted for the latter approach, and selected an individual study with all relevant attributes but conducted in adults (asthmatics and non-asthmatics).

In this study, Weiss et al. found that log-transformed dust concentrations of Bla g 1 and Bla g 2 were both significantly associated with longitudinal FEV1 decline (ΔFEV1), with multiple linear regression coefficients of −194.14 mL/year and −94.83 mL/year respectively, after adjusting for age, initial FEV1, ever smoking, and Der p 1, Der f 1, and Fel d 1 allergen concentrations [68]. The study did not report functions for asthmatics only, so we used values for the entire population, noting that the relationship between dust concentrations and FEV1 was not appreciably different for the non-asthmatic population than the population as a whole. We converted change in FEV1 (ΔFEV1) to change in FEV1% by dividing ΔFEV1 by FEV1 predicted, where FEV1 predicted was calculated using the NHANES equation below [69], using the average age and height reported in Table 1 of the Weiss study.

FEV 1 predicted = 0.554 0.013 * Age 0.0002 * Age 2 0.0001 * Height 2 ,

where age = 57.5 years, and height = 174.42 cm [68]. This conversion would ideally be based on individual height and age, but lacking this information, using the average values provides a reasonable approximation given the relatively narrow age and height ranges within the study (57.50 +/− 2.58 and 174.42 +/− 2.33, respectively). The resulting yearly decreases in FEV1% per unit increase in log transformed allergen concentration (log10 U/g) were −0.055% (SE = 0.013) and −0.027% (SE = 0.007) for Bla g 1 and Bla g 2, respectively.

Articles on mold and moisture were found in PubMed using the phrases [“childhood asthma AND mold”], [“childhood asthma AND moisture”], [“pediatric asthma AND mold”], and [“pediatric asthma AND moisture”], with 88, 13, 44, and 1 articles found, respectively. Five articles [7074] matched the inclusion criteria, but only two measured our lung function of interest, FEV1, and only one paper was a study of asthmatics [73]. If a house was characterized as damp, we used the coefficient reported in this paper to reduce FEV1% by 10.6% (SE = 4.95, 95% CI: 1.0 – 20.3). The coefficient reported in this paper was adjusted for unemployment, household income (above/below £200), respondent smoker, other smoker in house, and pet ownership.

Indoor pollutant exposure

For NO2 and PM2.5, daily 24-hour average exposures were estimated with regression models developed using the multi-zone simulation software output from CONTAM2.4c (NIST, Gaithersburg, MD,, an approach described in more detail elsewhere [75]. Briefly, within CONTAM, we selected the building most typical of Boston public housing and other low-income multi-family dwellings in Boston –a building 4 stories, 1940–1969 construction, and naturally ventilated [76]. A family of 2 adults and 2 children were simulated living in each 703 square foot apartment, which included a bedroom, bathroom, living room, and kitchen (Figure 1). Sources of NO2 included the gas stove used for cooking, the gas oven used for supplemental heat in the winter, and outdoors. Sources of PM2.5 included environmental tobacco smoke, cooking, and outdoors.

Because CONTAM could not be directly linked with the discrete event simulation model across all 50,000 children on an hourly basis for 10 years, we instead constructed regression models to explain variability in CONTAM outputs from a series of runs across key factors known to influence indoor concentrations. Regression predictors for pollutant/source combinations included terms consistent with a one-compartment box model for indoor concentrations, building characteristics (floor level, air exchange rate), occupant behaviors (kitchen and bathroom exhaust fan usage), and meteorological conditions (season, outdoor RH). Regression models had good predictive power (R2 from 0.89 to 0.98 across models) and physical interpretability, and outputs corresponded well with literature values [75]. Based on these equations, the 24-hour concentration was updated daily in the simulation model. See Additional file 1 for equations.

Regression models for RH were calculated in a similar way using CONTAM outputs, and were used to determine the likelihood of mold growth and dampness, described below. Indoor sources of RH included occupant behaviors such as breathing, showering, cooking, and dishwashing. Our models indicated a strong effect of season, air exchange rate, and outdoor RH and temperature [75]. We constructed exposure measures to be homologous with epidemiological evidence linking dampness/mold with FEV1%, described above. Mold growth was calculated using a set of differential equations developed to model mold growth on wooden material [77] which estimate a daily mold index (M) (Table 3). The index is a function of critical RH necessary for mold growth, temperature, surface quality, time, current RH, wood species, and some constant coefficients. Research suggests that models for pine sapwood most closely approximate materials with high nutrient content, as is the case with modern building materials [78], thus we applied the pine sapwood parameters to model mold growth in our buildings.

Table 3 Description of mold index developed to describe mold growth in wood[77]

All homes were assigned a mold index value of 0 at the beginning of the simulation. Each day a change in mold index (dM/dt) was calculated based on the parameters described, where if a critical RH was reached, the index would increase, and if the RH was below a critical value, the index would decrease. Once level 2 was reached, a house was categorized as “damp”, and had an effect on FEV1% as described above. This effect was reversible, that is, if the index dropped below 2 on subsequent days, the “damp” classification was eliminated, as was the effect on FEV1%. If level 4 was reached, then the house was permanently classified as “damp”, that is, even if the index dropped below 2, the house retained its “damp” value and effect on FEV1%. See Additional file 1 for equations.

For cockroach allergen, simulations of daily concentrations were infeasible, as models such as CONTAM are not applicable. Instead, distributions for Bla g 1 and Bla g 2 were calculated from HPHI data, where Peters et al. measured cockroach allergen in air, dust in bedrooms and dust in kitchens [39]. As the Weiss et al. epidemiological study used to derive the concentration-response function chose the highest concentration among kitchen, living room, and bedroom as their exposure metric, we utilized the kitchen measurements from Peters et al. 2007, which were consistently highest in HPHI. Peters presented multivariate models where cockroach allergen concentrations depended on having holes in walls and ceilings, and the degree of cleanliness and clutter in the apartment, termed housekeeping practices. We derived three distributions of cockroach allergen from the HPHI data, where homes were categorized as: 1.“with holes and below average housekeeping”, 2.“with holes and average or above average housekeeping”, and 3.“without holes and average or above average housekeeping”. Holes were defined as “open cracks or holes thicker than a dime found in the inside walls or ceilings”. Bla g 1 geometric mean concentrations were 143.5 (GSD = 3.6), 42.7 (GSD = 6.2), and 8.2 (GSD = 14.6) U/g, for each category respectively, and Bla g 2 geometric mean concentrations were 691.4 (GSD = 8.6), 117.3 (GSD = 9.0), and 21.9 (GSD = 12.5) U/g, respectively. Because the concentrations were skewed, daily cockroach allergen concentrations were drawn from the log transformed data, truncated at one SD.

Daily NO2, PM2.5, and cockroach allergen were multiplied by a factor of 0.7 to account for the time children spent inside their homes. The 0.7 factor was calculated using data in Table 15–3 of the EPA Exposure Factors Handbook [79], which lists average times spent on indoor and outdoor activities for children 6 to 8 years old. Indoor activities included sleeping, personal care, household work, eating, studying, playing, TV, and reading. We calculated the average time spent indoors, weighted for weekend and weekday differences. For the remaining time, we assumed that outdoor NO2 and PM2.5 reasonably represented exposures (lacking data on other microenvironments) and that there were no cockroach allergen exposures.

Outcomes: asthma exacerbations and health care utilization

Asthma outcomes were computed daily for each child based on changes in FEV1% associated with environmental exposures, derived from a prior model of the association between FEV1% and asthma symptoms or serious asthma events reported by Fuhlbrigge et al. [80]. The Fuhlbrigge study reported on 407 mild to moderate asthmatic children enrolled in the placebo branch of the Childhood Asthma Management Program (CAMP), aged 5–12 years, and recruited from 8 centers around the US [81]. In that asthma study, serious asthma event referred to ER visits, hospitalizations and oral steroid bursts, which were evaluated separately from asthma symptoms. We developed functions to approximate continuous associations between FEV1% and both asthma symptoms and serious asthma events, given only categorical models in the original publication, and we utilize other data streams to determine all outcomes of interest (which were not independently reported in the original publication).

First, the frequency of asthma symptoms was characterized in Fuhlbrigge et al. (listed in that article’s Figure 1), which shows the number of episode-free days per 4-month period across four categories of FEV1% (<60%, 60-79%, 80-99%, ≥ 100%). An episode-free day was defined as “a day with an asthma diary asthma score of 0, and no report of night awakening, morning and evening peak flow >80% personal best, no albuterol use for symptoms or prednisone use, absence from school as a result of asthma, or physician contact as a result of asthma”. We focused on the number of days with symptoms to be better aligned with our model structure. To convert this into a continuous function of FEV1%, we used the estimated midpoint of each FEV1% (50%, 70%, 90%, and 110%) category and fit the following polynomial expression:

P symptom_‐day = 2 . 95 FEV 1 % 3 6 . 93 FEV 1 % 2 + 4 . 68 FEV 1 % 0.27

where Psymptom_day is the daily probability of having a day with asthma symptoms as defined above. While the polynomial expression is clearly over determined given four categories of FEV1%, the intent was not to build a predictive model or a model for out-of-sample characterization, but simply to capture previously observed trends. The equation is valid for values of FEV1% between 0.5 and 1.2, a range of values consistent with our simulated study population.

A similar process was used to fit an equation predicting “serious asthma events”, defined in Fuhlbrigge et al. as oral steroid use, hospitalization, or emergency room visit. Table 3 of Fuhlbrigge et al. provides a multivariate regression model including the influence of FEV1% (again in four categories) as well as night awakenings and previous hospitalizations. To convert the reported odds ratios into a probability of a serious asthma event based on a continuous FEV1% scale, we first determined the baseline rate of serious asthma events and converted it to a probability of a serious asthma event. Fuhlbrigge et al. reported that their study population had a baseline rate of 0.26 serious asthma events per 4 month period, or approximately 0.0022 events per day (probability of 0.0022). Distributing this rate on a population-weighted basis following odds ratios and population numbers in Table 1 of Fuhlbrigge et al. yields daily event probabilities of 0.0068, 0.0032, 0.0022, and 0.0017 in the four FEV1% categories of decreasing severity. Fitting a polynomial expression to these values leads to a resulting equation of:

P serious event = 0.045 FEV 1 % 3 + 0.1277 FEV 1 % 2 0.1224 FEV 1 % + 0.0417

Where Pserious event is the daily probability of having a serious asthma event. As previously, the equation is valid for values of FEV1% between 0.5 and 1.2. Pserious event was multiplied by a season factor which accounted for residual variation due to seasonally dependent variables such as respiratory virus infections and outdoor allergen concentrations [9, 8285]. The season factors were 1.11, 0.60, 1.23, and 1.05 for spring, summer, fall, and winter respectively, and were estimated from monthly ER utilization raw data collected in the Quality Improvement study between 2006 and 2008, using Boston Medical Center Health Net Plan (BMCHP) data [86]. Similar trends were observed in a 23-year survey of asthma hospitalizations in Canada (Figure 1A,[87]).

Based on data published in the Fuhlbrigge study, if a child had a hospitalization due to asthma in the previous 12 months , their probability of having a serious asthma event increased (Table 3, [80]). We calculated this multiplicative factor following the same process described above, with the resulting polynomial equation:

O hospit = 45 . 7 FEV 1 % 3 + 129 . 7 FEV 1 % 2 124 . 4 FEV 1 % + 42 . 4

where Ohospit is the increased odds of having a serious asthma event given an asthma hospitalization in the last 12 months, and was equal to 1 if no hospitalization had occurred.

While the Pserious event equation provides a robust expression for serious asthma events, data are not provided in Fuhlbrigge et al. to determine the relative distribution among these events, a critical component for policy analysis. To approximate this distribution, we used asthma statistics reported by the Centers for Disease Control and Prevention. We assumed that these average hospitalization and ER visit rates apply to the population from the Fuhlbrigge study, and that the odds ratios from Table 3 apply to each individual serious asthma related event, allowing us to create distributions of hospitalization rates and emergency room visits proportional to those calculated above for asthma events. The remaining events are presumed to be oral steroid bursts.

For hospitalization, we estimated a baseline hospitalization rate due to asthma of 0.023 hospitalizations per year per asthmatic child, combining the National Hospital Discharge Survey [88] which reports an asthma pediatric hospitalization discharge rate of 24.8 per year per 10,000 children under 15 years of age and the National Health Interview Survey [89] which reports that 10.7% of children between 5 and 11 years have asthma. This was similar to the approach used by the US Environmental Protection Agency to characterize baseline asthma hospitalization rates for regulatory analyses [90]. For this model we used more recent data.

We constructed a polynomial equation to predict the daily probability of hospitalization based on FEV1% using the approach described above, with the resulting equation:

P hosp = 0.0013 FEV 1 % 3 + 0.0037 FEV 1 % 2 0.0036 FEV 1 % + 0.0012

where Phosp includes direct hospitalizations and transfers from the ER to the hospital.

For ER visits, we used an estimated emergency department visit rate of 0.1 per year per asthmatic [91] to build a similar equation, where the daily probability of going to the ER is:

P ER = 0.0057 FEV 1 % 3 + 0.0162 FEV 1 % 2 0.0155 FEV 1 % + 0.0053

Because 8% of ER visits result in hospitalization and are already accounted for in Phosp, we multiplied PER by 0.92 so as not to overestimate ER visits [92].

Oral steroid bursts were estimated by subtracting Phosp and PER from Pserious event.


Model outcomes and evaluation

Using the population baseline characteristics we simulated 50,000 children over 10 years. Average yearly health outcomes stratified by FEV1% are presented in Table 4. As expected, children with lower lung function had a higher incidence of days with asthma episodes and severe asthma events. Average health outcome rates across the total simulated population align closely with our baseline model inputs drawn from the literature: 0.023 hospitalizations/year [88, 89], 0.1 ER visits/year [91], 0.78 serious events/year [80]. The percentage of days with asthma symptoms per year was 41%, averaged across all simulated children. At the end of the simulation, 67% of the children had FEV1% > 80%, 32% had FEV1% between 60% and 80%, and 1% had FEV1% < 60%. This distribution is consistent with what has been observed in the field [93].

Table 4 Healthcare outcomes from a baseline simulation of 50,000 asthmatic children over 10 years

As a measure of disease progression, we found that 13.2% of children dropped to a more severe asthma classification over the 10 year simulation, solely based on a comparison of average FEV1% over the previous year to FEV1% asthma severity classification guidelines. This change in severity classification would then simulate the changes in asthma medication prescription that might occur at a yearly physical checkup. Because available evidence indicates a declining or unchanging lung function for asthmatic children over time, all else being equal, children would not move to a less severe classification within our simulation. Average FEV1% decline for compliant (assigned to get a controller medication) versus non-compliant (not assigned to get a controller medication) children was 0.18% and 0.37% respectively. Of the children who were non-compliant, 15% dropped to a more severe asthma classification over the 10 year simulation, compared to 12.5% of the compliant children.

For indoor pollutants, 24 hour average indoor NO2 was normally distributed with a mean concentration of 54 ppb (SD = 23 ppb, range 3–140 ppb). In the literature, reported values from field studies in Boston public housing include 43 ppb (SD = 20 ppb) measured in the kitchen, and 36 ppb (SD = 17 ppb) measured in the living room [19]. Another study measured an average of 19.6 ppb (SD = 11.0 ppb, range 5.7-61 ppb) [64]. The concentrations we simulated were above those reported in field studies, likely because of the high percentage of simulated homes with gas stoves and non-operational kitchen exhaust fans. All homes in the highest NO2 quartile (mean = 80 ppb) had gas stoves, 99.95% had non-operating kitchen exhaust fans, and 42% had low air exchange rates (i.e. leaky category 1, as defined in Table 1). In contrast, in the lowest NO2 quartile (mean = 22 ppb), only 56% of homes owned a gas stove, 51% had a non-operating fan and 16% had low air exchange rates.

Indoor PM2.5 was lognormally distributed, with a mean of 55 μg/m3 (SD = 34 μg/m3, range 14–394 μg/m3). In the literature, the mean PM2.5 concentrations in homes where smoking was rarely reported was 20.3 μg/m3 (SD = 12.5, range 6.77–74.9 μg/m3)[64]. This did not include ETS, estimated to increase indoor PM2.5 concentrations between 7 and 49 μg/m3[65, 9497]. The PTEAM study reported indoor PM2.5 concentrations of 48 μg/m3 averaged over 178 smoker and non-smoker homes [65]. In the simulation, in the highest PM2.5 quartile (mean = 101 μg/m3), 82% of homes had light or heavy smokers, 94% had non-operating exhaust fans, and 44% had low air exchange rates. In contrast, in the lowest PM2.5 quartile (mean = 25 μg/m3), only 1% of homes had smokers, 66% had non-operating exhaust fans, and 5% had low air exchange rates.

Cockroach allergen was distributed according to housekeeping practices and the presence of holes in walls and ceilings of the home. Bla g 1 and Bla g 2 median concentrations were 180 U/g (range 90–402 U/g) and 1270 U/g (range 245–2870 U/g) for below average housekeeping and presence of holes, 66.7 U/g (range 12.9-166 U/g) and 220 U/g (range 18–476 U/g) for above average housekeeping and presence of holes, and 20.3 U/g (range 0.8-51.8 U/g) and 49.1 U/g (range 5.9-92.7 U/g) for above average housekeeping and no holes. Peters et al. reported median kitchen floor Bla g 1 concentrations of approximately 48 U/g (Figure 1), [98]) measured in a longitudinal component of HPHI, and median kitchen floor concentrations of 61.8 U/g and 198 U/g for Bla g 1 and Bla g 2 in a cross-sectional component of the same study [39]. Gergen et al. reported geometric mean Bla g 1 values of 68.7 U/g in kitchens of inner-city homes [99], and Arbes et al. [100] found a geometric mean concentration of Bla g 1 in kitchen floors to be 287 U/g.

Approximately 19% of homes were classified as “damp” at the end of a 10 year simulation, comparable to prevalence values observed in studies of Boston public housing, where Hynes et al. reported 20% of homes had observed mold growth [34], and Brugge et al. reported that 43% of residents had smelled or seen mold in their homes [32].

As shown in Figure 2, average daily pollutant exposure decreased FEV1% by 3.4%, 4.0%, 1.9%, 8.4%, and 5.5% for PM2.5, NO2, mold, Bla g 1 and Bla g 2 respectively, relative to a no-exposure scenario, and are similar to the literature values we used as inputs. Cockroach allergen had the highest impact on FEV1%, followed by NO2 and PM2.5. In the 19% of children living in a damp (i.e. moldy) home, FEV1% decreased by 9.8%, higher than any other pollutant.

Figure 2
figure 2

Average daily decrease in FEV1% for 50,000 simulated children over 10 years, relative to a no-exposure scenario.

Application example, fixing exhaust fans

We conducted a simulation experiment with a simple hypothetical intervention by simulating what would happen if the kitchen and bathroom exhaust fans were improved for the entire cohort. To evaluate this we set the exhaust fans in the kitchen and bathroom to be operational for 100% of the children, versus 13% in the baseline scenario, ran the simulation for 10 years, and compared results to the baseline 10 year simulation. In this scenario, average daily NO2 concentrations were reduced to 28 ppb (−48%), and PM2.5 concentrations were reduced to 35 μg/m3 (−36%). Cockroach allergen concentrations remained the same, and the percentage of homes classified as damp at the end of the simulation was 18.7%, compared to 19.0% in the baseline simulation.

In terms of health outcomes, mean daily FEV1% was 87.8% (range: 52-119%), 3% higher than the baseline scenario. Only 10.8% of children dropped to a more severe asthma category compared to 13.2% at baseline. Days with asthma symptoms, days with serious asthma events, ER visits, hospitalizations, and oral steroid use dropped 2%, 7.5%, 9.1%, 5.5%, and 8.0% respectively across all asthmatics.

For validation and evaluation purposes, we compared the magnitude of our results to relevant intervention studies as a plausibility check for our model. Unfortunately very few field studies have measured the effect of building interventions on asthma outcomes and none are directly comparable. Morgan et al. reported an 11-15% reduction in ER visits after an environmental intervention trial combining education with reducing allergen and ETS concentrations [15]. A research group in New Zealand reported a 46% drop in NO2 when they replaced a majority of unflued gas heaters with non-polluting heaters, with a corresponding reduction in asthma wheezing of 29% (non-significant) and a 45% reduction in sleep disturbed by wheeze (significant)[101]. Kattan et al. published on the Inner City Asthma study, showing that remediation reducing allergens to which children were specifically allergic (dust mites, passive smoking, cockroaches, pets, and rodents) decreased environmental allergens, and corresponded to a 19% decrease in unscheduled clinic visits. Since none of these interventions are directly comparable to our simulated interventions, all we can conclude is that our numbers are plausible. Further analysis comparing other interventions will be the subject of future applications of our simulation model. One approach we will take to reduce uncertainties is to focus on the marginal changes in exposures and outcomes when evaluating the effect of building interventions, rather than focusing on the absolute number of outcomes at baseline and post-intervention.



One limitation in our set of assumptions was the inability to model in more detail the cumulative effects of pollutants on FEV1%. We assumed simple additivity of effects on FEV1%, but given different mechanisms of action, the true combined effect may be different (with potential synergistic or antagonistic relationships). There is insufficient evidence in the current published literature to allow us to model these relationships. Similarly, we did not have adequate information to model the effect of non-chemical factors such as socioeconomic status, stress, or race, some of which are known to influence medication compliance or to be related to asthma outcomes. As more data become available on these relationships the model can be expanded and refined.

The equations linking FEV1% to asthma events were developed based on data from the CAMP study, which includes children with mild to moderate asthma from 8 centers across the US [81]. Thus there are some limitations in interpreting the results for the severe asthmatics as classified in our model, as well as likely differences in relationships between the CAMP children and our simulated low-income children living in multi-family housing. It is possible that using the CAMP data underestimates the severity of the response compared to lower-income populations. The same is true for our baseline estimates of asthma outcomes, which were based on national databases or the CAMP study.

The model as parameterized here is specific to multi-family housing with resident characteristics of inner-city Boston residents. Because the model relies heavily on data from the HPHI project, which involved a small number of households from a few housing developments in Boston, the generalizability of our numerical results is constrained to housing developments with similar construction (i.e. multifamily buildings) and heating systems. High rise buildings and buildings with different heating systems would require separate analysis. However, we know that this setting is relevant to many low-income urban populations where asthma prevalence is elevated and associated housing-related risk factors are also common.

The asthma medication module is an oversimplified version of what happens in the real world with respect to asthma severity classification, asthma medication prescription, and adherence. We classified asthmatics solely based on FEV1%, for the purpose of determining medications, although severity classification is far more complex [8], and medication prescription and adherence are influenced by many factors, including asthma status, access to health-care, income status, emotional status, and social experience [102]. Also, the model does not incorporate daily symptoms and level of asthma controls, and how this influences medication use, or cases where symptoms worsen over time prior to becoming a severe attack. Future work should include developing more complex models of asthma medication use which takes into account these many factors.


In order to extend this framework to other settings and populations, changes in the current model can be implemented by modifying any of the parameters presented in Table 1 and elsewhere (e.g., Additional file 1), which include resident characteristics and behavior. That said, our findings are likely robust to a number of basic demographic assumptions. Even the number of residents in the unit would have a small influence on our findings given our model structure, within a reasonable range of values, so the quantitative conclusions readily extend to one-adult households.

Changes in building type or building characteristics can be implemented by selecting a different building plan from the 209 buildings available in CONTAM [76] and repeating the modeling process described, making the modeling approach generalizable.

Potential applications of this model include evaluating the effect of building construction and public housing policy changes on pollutants and asthma, evaluating green building practices on indoor air quality and health of residents in remodeled homes, and conducting cost-benefit analyses comparing energy savings to cost of healthcare utilization.


We developed a discrete event simulation model of pediatric asthma that can be used to evaluate the effect of building interventions on multiple air pollutants and allergens, as well as on healthcare (asthma attacks, hospitalizations, ER visits) and asthma medicine utilization. This work presents a novel framework linking environmental exposures to FEV1% and pediatric asthma outcomes, can be expanded to any housing type, and can be refined as more data becomes available regarding the different relationships. The model can help determine health-optimal strategies as buildings are renovated or constructed, and can be used to consider a variety of environmental and non-environmental interventions targeting pediatric asthma. The outputs are directly relevant to policy discussions in low-income urban communities.



Boston Medical Center Health Net Plan


Childhood Asthma Management Program


Centers for Disease Control and Prevention


chronic obstructive pulmonary disease


emergency room


environmental tobacco smoke


forced expiratory flow from the 25th percentile to the 75th percentile of FVC


forced expiratory volume in 1 second


percent predicted forced expiratory volume in 1 second


forced vital capacity


Healthy Public Housing Initiative


National Heart, Lung and Blood Institute

NO2 :

Nitrogen dioxide

O3 :



Peak expiratory flow


Parts per billion


Relative humidity


Severe asthma-related event


Standard error.


  1. Merrill CT, Elixhauser A: Hospitalization in the United States, 2002: HCUP Fact Book No. 6. 2005, Rockville, MD: Agency for Healthcare Research and Quality

    Google Scholar 

  2. ARC: Asthma in New England. Part II: children. 2004, Dorchester, MA: New England Asthma Regional Council

    Google Scholar 

  3. Crain EF, Weiss KB, Bijur PE, Hersh M, Westbrook L, Stein RE: An estimate of the prevalence of asthma and wheezing among inner-city children. Pediatrics. 1994, 94: 356-362.

    CAS  Google Scholar 

  4. Joseph CL, Foxman B, Leickly FE, Peterson E, Ownby D: Prevalence of possible undiagnosed asthma and associated morbidity among urban schoolchildren. J Pediatr. 1996, 129: 735-742.

    CAS  Google Scholar 

  5. Spiegel J, Love AS, Wood PR, Griffith M, Taylor KR, Williams SG, Redd SC: The Inner-City Asthma Intervention: description of a community-based implementation of an evidence-based approach to asthma management. Ann Allergy Asthma Immunol. 2006, 97: S6-S10.

    Google Scholar 

  6. Sharma HP, Hansel NN, Matsui E, Diette GB, Eggleston P, Breysse P: Indoor environmental influences on children's asthma. Pediatr Clin North Am. 2007, 54: 103-120. ix

    Google Scholar 

  7. IOM: Clearing the air: asthma and indoor exposures. 2000, Washington DC: National Academy Press

    Google Scholar 

  8. NHLBI: National Asthma Education and Prevention Program Expert Panel Report 3: Guidelines for the diagnosis and management of asthma. 2007, Bethesda, MD: National Heart Lung and Blood Institute, National Institutes of Health

    Google Scholar 

  9. Johnston NW, Sears MR: Asthma exacerbations. 1: epidemiology. Thorax. 2006, 61: 722-728.

    CAS  Google Scholar 

  10. Spahn JD, Cherniack R, Paull K, Gelfand EW: Is forced expiratory volume in one second the best measure of severity in childhood asthma?. Am J Respir Crit Care Med. 2004, 169: 784-786.

    Google Scholar 

  11. Krieger J, Takaro TK, Song L, Beaudet N, Edwards K: A randomized controlled trial of asthma self-management support comparing clinic-based nurses and in-home community health workers: the Seattle-King County Healthy Homes II Project. Arch Pediatr Adolesc Med. 2009, 163: 141-149.

    Google Scholar 

  12. Postma J, Karr C, Kieckhefer G: Community health workers and environmental interventions for children with asthma: a systematic review. J Asthma. 2009, 46: 564-576.

    Google Scholar 

  13. Clougherty JE, Levy JI, Hynes HP, Spengler JD: A longitudinal analysis of the efficacy of environmental interventions on asthma-related quality of life and symptoms among children in urban public housing. J Asthma. 2006, 43: 335-343.

    CAS  Google Scholar 

  14. Levy JI, Brugge D, Peters JL, Clougherty JE, Saddler SS: A community-based participatory research study of multifaceted in-home environmental interventions for pediatric asthmatics in public housing. Soc Sci Med. 2006, 63: 2191-2203.

    Google Scholar 

  15. Morgan WJ, Crain EF, Gruchalla RS, O'Connor GT, Kattan M, Evans R, Stout J, Malindzak G, Smartt E, Plaut M: Results of a home-based environmental intervention among urban children with asthma. N Engl J Med. 2004, 351: 1068-1080.

    CAS  Google Scholar 

  16. Clougherty JE, Kubzansky LD, Spengler JD, Levy JI: Ancillary benefits for caregivers of children with asthma participating in an environmental intervention study to alleviate asthma symptoms. J Urban Health. 2009, 86: 214-229.

    Google Scholar 

  17. Crocker DD, Kinyota S, Dumitru GG, Ligon CB, Herman EJ, Ferdinands JM, Hopkins DP, Lawrence BM, Sipe TA: Effectiveness of home-based, multi-trigger, multicomponent interventions with an environmental focus for reducing asthma morbidity a community guide systematic review. Am J Prev Med. 2011, 41: S5-S32.

    Google Scholar 

  18. Adamkiewicz G, Zota AR, Fabian MP, Chahine T, Julien R, Spengler JD, Levy JI: Moving environmental justice indoors: understanding structural influences on residential exposure patterns in low-income communities. Am J Public Health. 2011, 101: S238-S245.

    Google Scholar 

  19. Zota A, Adamkiewicz G, Levy JI, Spengler JD: Ventilation in public housing: implications for indoor nitrogen dioxide concentrations. Indoor Air. 2005, 15: 393-401.

    CAS  Google Scholar 

  20. Lee K, Levy JI, Yanagisawa Y, Spengler JD, Billick IH: The Boston residential nitrogen dioxide characterization study: classification and prediction of indoor NO2 exposure. J Air Waste Manag Assoc. 1998, 48: 736-742.

    CAS  Google Scholar 

  21. Law AM, Kelton WD: Simulation modeling and analysis. 2000, Boston, MA: McGraw-Hill Higher Education

    Google Scholar 

  22. Heeg BM, Damen J, Buskens E, Caleo S, de Charro F, van Hout BA: Modelling approaches: the case of schizophrenia. Pharmacoeconomics. 2008, 26: 633-648.

    Google Scholar 

  23. McKenzie FE, Wong RC, Bossert WH: Discrete-event simulation models of plasmodium falciparum malaria. Simulation. 1998, 71: 250-261.

    Google Scholar 

  24. Chen J, Alemao E, Yin D, Cook J: Development of a diabetes treatment simulation model: with application to assessing alternative treatment intensification strategies on survival and diabetes-related complications. Diabetes Obes Metab. 2008, 10 (Suppl 1): 33-42.

    Google Scholar 

  25. Tosteson AN, Stout NK, Fryback DG, Acharyya S, Herman BA, Hannah LG, Pisano ED: Cost-effectiveness of digital mammography breast cancer screening. Ann Intern Med. 2008, 148: 1-10.

    Google Scholar 

  26. Fryback DG, Stout NK, Rosenberg MA, Trentham-Dietz A, Kuruchittham V, Remington PL: The wisconsin breast cancer epidemiology simulation model. J Natl Cancer Inst Monogr. 2006, 36: 37-47.

    Google Scholar 

  27. Fuhlbrigge AL, Bae SJ, Weiss ST, Kuntz KM, Paltiel AD: Cost-effectiveness of inhaled steroids in asthma: impact of effect on bone mineral density. J Allergy Clin Immunol. 2006, 117: 359-366.

    CAS  Google Scholar 

  28. Wu AC, Paltiel AD, Kuntz KM, Weiss ST, Fuhlbrigge AL: Cost-effectiveness of omalizumab in adults with severe asthma: results from the Asthma Policy Model. J Allergy Clin Immunol. 2007, 120: 1146-1152.

    CAS  Google Scholar 

  29. Paltiel AD, Fuhlbrigge AL, Kitch BT, Liljas B, Weiss ST, Neumann PJ, Kuntz KM: Cost-effectiveness of inhaled corticosteroids in adults with mild-to-moderate asthma: results from the asthma policy model. J Allergy Clin Immunol. 2001, 108: 39-46.

    CAS  Google Scholar 

  30. Le Lay A, Despiegel N, Francois C, Duru G: Can discrete event simulation be of use in modelling major depression?. Cost Eff Resour Alloc. 2006, 4: 19-

    Google Scholar 

  31. Frey U: Predicting asthma control and exacerbations: chronic asthma as a complex dynamic model. Curr Opin Allergy Clin Immunol. 2007, 7: 223-230.

    Google Scholar 

  32. Brugge D, Rice PW, Terry P, Howard L, Best J: Housing conditions and respiratory health in a Boston public housing community. New Solut. 2001, 11: 149-164.

    CAS  Google Scholar 

  33. Digenis-Bury EC, Brooks DR, Chen L, Ostrem M, Horsburgh CR: Use of a population-based survey to describe the health of Boston public housing residents. Am J Public Health. 2008, 98: 85-91.

    Google Scholar 

  34. Hynes H, Brugge D, Watts J, Lally J: Public health and the physical environment in Boston Public Housing: a community-based survey and action agenda. Planning Practice & Research. 2000, 15: 21-49.

    Google Scholar 

  35. Brugge D, Melly S, Finkelman A, Russell M, Bradeen L, Perez R: A community-based participatory survey of public housing conditions and associations between renovations and possible building-related symptoms. Appl Environ Sci Pub Health. 2003, I: 89-101.

    Google Scholar 

  36. Moorman JE, Rudd RA, Johnson CA, King M, Minor P, Bailey C, Scalia MR, Akinbami LJ: National surveillance for asthma–United States, 1980–2004. MMWR Surveill Summ. 2007, 56: 1-54.

    Google Scholar 

  37. Paull K, Covar R, Jain N, Gelfand EW, Spahn JD: Do NHLBI lung function criteria apply to children? A cross-sectional evaluation of childhood asthma at National Jewish Medical and Research Center, 1999–2002. Pediatr Pulmonol. 2005, 39: 311-317.

    Google Scholar 

  38. Wildhaber JH, Sznitman J, Harpes P, Straub D, Moller A, Basek P, Sennhauser FH: Correlation of spirometry and symptom scores in childhood asthma and the usefulness of curvature assessment in expiratory flow-volume curves. Respir Care. 2007, 52: 1744-1752.

    Google Scholar 

  39. Peters JL, Levy JI, Rogers CA, Burge HA, Spengler JD: Determinants of allergen concentrations in apartments of asthmatic children living in public housing. J Urban Health. 2007, 84: 185-197.

    Google Scholar 

  40. O'Byrne PM, Pedersen S, Lamm CJ, Tan WC, Busse WW: Severe exacerbations and decline in lung function in asthma. Am J Respir Crit Care Med. 2009, 179: 19-24.

    Google Scholar 

  41. Levy JI, Welker-Hood LK, Clougherty JE, Dodson RE, Steinbach S, Hynes HP: Lung function, asthma symptoms, and quality of life for children in public housing in Boston: a case-series analysis. Environ Health. 2004, 3: 13-

    Google Scholar 

  42. Janson SL, Earnest G, Wong KP, Blanc PD: Predictors of asthma medication nonadherence. Heart Lung. 2008, 37: 211-218.

    Google Scholar 

  43. Delfino RJ, Quintana PJ, Floro J, Gastanaga VM, Samimi BS, Kleinman MT, Liu LJ, Bufalino C, Wu CF, McLaren CE: Association of FEV1 in asthmatic children with personal and microenvironmental exposure to airborne particulate matter. Environ Health Perspect. 2004, 112: 932-941.

    CAS  Google Scholar 

  44. Delfino RJ, Staimer N, Tjoa T, Gillen D, Kleinman MT, Sioutas C, Cooper D: Personal and ambient air pollution exposures and lung function decrements in children with asthma. Environ Health Perspect. 2008, 116: 550-558.

    Google Scholar 

  45. Lewis TC, Robins TG, Dvonch JT, Keeler GJ, Yip FY, Mentz GB, Lin X, Parker EA, Israel BA, Gonzalez L, Hill Y: Air pollution-associated changes in lung function among asthmatic children in Detroit. Environ Health Perspect. 2005, 113: 1068-1075.

    Google Scholar 

  46. Koenig JQ, Larson TV, Hanley QS, Rebolledo V, Dumler K, Checkoway H, Wang SZ, Lin D, Pierson WE: Pulmonary function changes in children associated with fine particulate matter. Environ Res. 1993, 63: 26-38.

    CAS  Google Scholar 

  47. Svendsen ER, Yeatts KB, Peden D, Orton S, Alexis NE, Creason J, Williams R, Neas L: Circulating neutrophil CD14 expression and the inverse association of ambient particulate matter on lung function in asthmatic children. Ann Allergy Asthma Immunol. 2007, 99: 244-253.

    Google Scholar 

  48. Liu L, Poon R, Chen L, Frescura AM, Montuschi P, Ciabattoni G, Wheeler A, Dales R: Acute effects of air pollution on pulmonary function, airway inflammation, and oxidative stress in asthmatic children. Environ Health Perspect. 2009, 117: 668-674.

    CAS  Google Scholar 

  49. Dales R, Chen L, Frescura AM, Liu L, Villeneuve PJ: Acute effects of outdoor air pollution on forced expiratory volume in 1 s: a panel study of schoolchildren with asthma. Eur Respir J. 2009, 34: 316-323.

    CAS  Google Scholar 

  50. Barraza-Villarreal A, Sunyer J, Hernandez-Cadena L, Escamilla-Nunez MC, Sienra-Monge JJ, Ramirez-Aguilar M, Cortez-Lugo M, Holguin F, Diaz-Sanchez D, Olin AC, Romieu I: Air pollution, airway inflammation, and lung function in a cohort study of Mexico City schoolchildren. Environ Health Perspect. 2008, 116: 832-838.

    Google Scholar 

  51. Epton MJ, Dawson RD, Brooks WM, Kingham S, Aberkane T, Cavanagh JA, Frampton CM, Hewitt T, Cook JM, McLeod S: The effect of ambient air pollution on respiratory health of school children: a panel study. Environ Health. 2008, 7: 16-

    Google Scholar 

  52. O'Connor GT, Neas L, Vaughn B, Kattan M, Mitchell H, Crain EF, Evans R, Gruchalla R, Morgan W, Stout J: Acute respiratory health effects of air pollution on children with asthma in US inner cities. J Allergy Clin Immunol. 2008, 121: 1133-1139. e1131

    Google Scholar 

  53. Allen RW, Mar T, Koenig J, Liu LJ, Gould T, Simpson C, Larson T: Changes in lung function and airway inflammation among asthmatic children residing in a woodsmoke-impacted urban area. Inhal Toxicol. 2008, 20: 423-433.

    CAS  Google Scholar 

  54. Holguin F, Flores S, Ross Z, Cortez M, Molina M, Molina L, Rincon C, Jerrett M, Berhane K, Granados A, Romieu I: Traffic-related exposures, airway function, inflammation, and respiratory symptoms in children. Am J Respir Crit Care Med. 2007, 176: 1236-1242.

    CAS  Google Scholar 

  55. Rosenlund M, Forastiere F, Porta D, De Sario M, Badaloni C, Perucci CA: Traffic-related air pollution in relation to respiratory symptoms, allergic sensitisation and lung function in schoolchildren. Thorax. 2009, 64: 573-580.

    CAS  Google Scholar 

  56. Nitschke M, Pilotto LS, Attewell RG, Smith BJ, Pisaniello D, Martin J, Ruffin RE, Hiller JE: A cohort study of indoor nitrogen dioxide and house dust mite exposure in asthmatic children. J Occup Environ Med. 2006, 48: 462-469.

    CAS  Google Scholar 

  57. Peters JM, Avol E, Gauderman WJ, Linn WS, Navidi W, London SJ, Margolis H, Rappaport E, Vora H, Gong H, Thomas DC: A study of twelve Southern California communities with differing levels and types of air pollution. II. Effects on pulmonary function. Am J Respir Crit Care Med. 1999, 159: 768-775.

    CAS  Google Scholar 

  58. Trenga CA, Sullivan JH, Schildcrout JS, Shepherd KP, Shapiro GG, Liu LJ, Kaufman JD, Koenig JQ: Effect of particulate air pollution on lung function in adult and pediatric subjects in a Seattle panel study. Chest. 2006, 129: 1614-1622.

    Google Scholar 

  59. Romieu I, Sienra-Monge JJ, Ramirez-Aguilar M, Tellez-Rojo MM, Moreno-Macias H, Reyes-Ruiz NI, del Rio-Navarro BE, Ruiz-Navarro MX, Hatch G, Slade R, Hernandez-Avila M: Antioxidant supplementation and lung functions among children with asthma exposed to high levels of air pollutants. Am J Respir Crit Care Med. 2002, 166: 703-709.

    Google Scholar 

  60. Liu LJ, Box M, Kalman D, Kaufman J, Koenig J, Larson T, Lumley T, Sheppard L, Wallace L: Exposure assessment of particulate matter for susceptible populations in Seattle. Environ Health Perspect. 2003, 111: 909-918.

    Google Scholar 

  61. MacIntosh DL, Minegishi T, Kaufman M, Baker BJ, Allen JG, Levy JI, Myatt TA: The benefits of whole-house in-duct air cleaning in reducing exposures to fine particulate matter of outdoor origin: A modeling analysis. J Expo Sci Environ Epidemiol. 2010, 20: 213-224.

    CAS  Google Scholar 

  62. Monn C, Fuchs A, Hogger D, Junker M, Kogelschatz D, Roth N, Wanner HU: Particulate matter less than 10 microns (PM10) and fine particles less than 2.5 microns (PM2.5): relationships between indoor, outdoor and personal concentrations. Sci Total Environ. 1997, 208: 15-21.

    CAS  Google Scholar 

  63. Levy JI, Lee K, Spengler JD, Yanagisawa Y: Impact of residential nitrogen dioxide exposure on personal exposure: an international study. J Air Waste Manag Assoc. 1998, 48: 553-560.

    CAS  Google Scholar 

  64. Baxter LK, Clougherty JE, Laden F, Levy JI: Predictors of concentrations of nitrogen dioxide, fine particulate matter, and particle constituents inside of lower socioeconomic status urban homes. J Expo Sci Environ Epidemiol. 2007, 17: 433-444.

    CAS  Google Scholar 

  65. Ozkaynak H, Xue J, Spengler J, Wallace L, Pellizzari E, Jenkins P: Personal exposure to airborne particles and metals: results from the Particle TEAM study in Riverside, California. J Expo Anal Environ Epidemiol. 1996, 6: 57-78.

    CAS  Google Scholar 

  66. Long CM, Suh HH, Catalano PJ, Koutrakis P: Using time- and size-resolved particulate data to quantify indoor penetration and deposition behavior. Environ Sci Technol. 2001, 35: 2089-2099.

    CAS  Google Scholar 

  67. Lin YC, Su HJ, Hsiue TR, Lee CH, Chen CW, Guo YL: Levels of house dust mite-specific IgE and cockroach-specific IgE and their association with lower pulmonary function in Taiwanese children. Chest. 2002, 121: 347-353.

    CAS  Google Scholar 

  68. Weiss ST, O'Connor GT, DeMolles D, Platts-Mills T, Sparrow D: Indoor allergens and longitudinal FEV1 decline in older adults: the Normative Aging Study. J Allergy Clin Immunol. 1998, 101: 720-725.

    CAS  Google Scholar 

  69. Hankinson JL, Odencrantz JR, Fedan KB: Spirometric reference values from a sample of the general U.S. population. Am J Respir Crit Care Med. 1999, 159: 179-187.

    CAS  Google Scholar 

  70. Cuijpers CE, Swaen GM, Wesseling G, Sturmans F, Wouters EF: Adverse effects of the indoor environment on respiratory health in primary school children. Environ Res. 1995, 68: 11-23.

    CAS  Google Scholar 

  71. Norback D, Bjornsson E, Janson C, Palmgren U, Boman G: Current asthma and biochemical signs of inflammation in relation to building dampness in dwellings. Int J Tuberc Lung Dis. 1999, 3: 368-376.

    CAS  Google Scholar 

  72. Strachan DP, Sanders CH: Damp housing and childhood asthma; respiratory effects of indoor air temperature and relative humidity. J Epidemiol Community Health. 1989, 43: 7-14.

    CAS  Google Scholar 

  73. Williamson IJ, Martin CJ, McGill G, Monie RD, Fennerty AG: Damp housing and asthma: a case–control study. Thorax. 1997, 52: 229-234.

    CAS  Google Scholar 

  74. Have WH, Berg NJ, Palenz J, Aalderen WMC, Bindels PJE: Residential exposure to mould and dampness is associated with adverse respiratory health. Clin Exp Allergy. 2007, 37: 1827-1832.

    Google Scholar 

  75. Fabian P, Adamkiewicz G, Levy J: Simulating indoor concentrations of NO2 and PM2.5 in multi-family housing for use in health-based intervention modeling. Indoor Air. 2011, 22 (1): 12-23.

    Google Scholar 

  76. Persily A, Musser A, Leber D: A collection of homes to represent the U.S. housing stock. 2006, Washington DC: National Institute of Standards and Technology

    Google Scholar 

  77. Hukka A, Viitanen H: A mathematical model of mold growth on wooden material. Wood Science and Technology. 1999, 33: 475-485.

    CAS  Google Scholar 

  78. Viitanen H, Ojanen T: Improved model to predict mold growth in building materials. 2007, ASHRAE

    Google Scholar 

  79. EPA: U.S. EPA. Exposure Factors Handbook (External Review Draft) 2009 Update. 2009, Washington, DC: U.S. Environmental Protection Agency

    Google Scholar 

  80. Fuhlbrigge AL, Weiss ST, Kuntz KM, Paltiel AD: Forced expiratory volume in 1 second percentage improves the classification of severity among children with asthma. Pediatrics. 2006, 118: e347-355.

    Google Scholar 

  81. Fabian P, Adamkiewicz G, Levy J: The childhood asthma management program (CAMP): design, rationale, and methods. childhood asthma management program research group. Control Clin Trials. 1999, 20: 91-120.

    Google Scholar 

  82. Johnston NW, Johnston SL, Duncan JM, Greene JM, Kebadze T, Keith PK, Roy M, Waserman S, Sears MR: The September epidemic of asthma exacerbations in children: A search for etiology. J Allergy Clin Immunol. 2005, 115: 132-138.

    Google Scholar 

  83. Johnston SL, Pattemore PK, Sanderson G, Smith S, Campbell MJ, Josephs LK, Cunningham A, Robinson BS, Myint SH, Ward ME: The relationship between upper respiratory infections and hospital admissions for asthma: a time-trend analysis. Am J Respir Crit Care Med. 1996, 154: 654-660.

    CAS  Google Scholar 

  84. Johnston SL, Pattemore PK, Sanderson G, Smith S, Lampe F, Josephs L, Symington P, O'Toole S, Myint SH, Tyrrell DA: Community study of role of viral infections in exacerbations of asthma in 9–11 year old children. BMJ (Clinical research ed. 1995, 310: 1225-1229.

    CAS  Google Scholar 

  85. Jariwala SP, Kurada S, Moday H, Thanjan A, Bastone L, Khananashvili M, Fodeman J, Hudes G, Rosenstreich D: Association between tree pollen counts and asthma ED Visits in a high-density urban center. J Asthma. 2011, 48: 442-448.

    CAS  Google Scholar 

  86. Sandel M: 2006-2008 Boston Medical Center Health Net Plan data Quality Improvement study, personal communication with P. Fabian. 2011, Boston

    Google Scholar 

  87. Johnston NW, Johnston SL, Norman GR, Dai J, Sears MR: The September epidemic of asthma hospitalization: school children as disease vectors. J Allergy Clin Immunol. 2006, 117: 557-562.

    Google Scholar 

  88. CDC: National Hospital Discharge Survey. 2007, Atlanta: Centers for Disease Control

    Google Scholar 

  89. CDC: Summary Health Statistics for US Children, National Health Interview Survey. 2009, Atlanta: Centers for Disease Control

    Google Scholar 

  90. US Environmental Protection Agency: The Benefits and Costs of the Clean Air Act: 1990 to 2020. 2011, Washington, DC: Office of Air and Radiation

    Google Scholar 

  91. Akinbami LJ: The State of Childhood Asthma, United States, 1980–2005. 2006, Atlanta: CDC

    Google Scholar 

  92. DPHMA: The burden of asthma in Massachusetts. 2009, Boston: Department of Public Health

    Google Scholar 

  93. Fuhlbrigge AL, Kitch BT, Paltiel AD, Kuntz KM, Neumann PJ, Dockery DW, Weiss ST: FEV(1) is associated with risk of asthma attacks in a pediatric population. J Allergy Clin Immunol. 2001, 107: 61-67.

    CAS  Google Scholar 

  94. Dockery DW, Spengler JD: Personal exposure to respirable particulates and sulfates. J Air Pollut Control Assoc. 1981, 31: 153-159.

    CAS  Google Scholar 

  95. Chahine T, Schultz B, Zartarian V, Subramanian SV, Spengler J, Hammitt JK, Levy JI: Modeling geographic and demographic variability in residential concentrations of environmental tobacco smoke using national datasets. J Exp Sci and Env Epi. 2011, 21: 646-655.

    CAS  Google Scholar 

  96. Myatt TA, Minegishi T, Allen JG, Macintosh DL: Control of asthma triggers in indoor air with air cleaners: a modeling analysis. Environ Health. 2008, 7: 43-

    Google Scholar 

  97. Klepeis NE, Nazaroff WW: Modeling residential exposure to secondhand tobacco smoke. Atm Environ. 2006, 40: 4393-4407.

    CAS  Google Scholar 

  98. Peters JL, Levy JI, Muilenberg ML, Coull BA, Spengler JD: Efficacy of integrated pest management in reducing cockroach allergen concentrations in urban public housing. J Asthma. 2007, 44: 455-460.

    Google Scholar 

  99. Gergen PJ, Mortimer KM, Eggleston PA, Rosenstreich D, Mitchell H, Ownby D, Kattan M, Baker D, Wright EC, Slavin R, Malveaux F: Results of the national cooperative inner-city asthma study (NCICAS) environmental intervention to reduce cockroach allergen exposure in inner-city homes. J Allergy Clin Immunol. 1999, 103: 501-506.

    CAS  Google Scholar 

  100. Arbes SJ, Sever M, Mehta J, Gore JC, Schal C, Vaughn B, Mitchell H, Zeldin DC: Abatement of cockroach allergens (Bla g 1 and Bla g 2) in low-income, urban housing: month 12 continuation results. J Allergy Clin Immunol. 2004, 113: 109-114.

    Google Scholar 

  101. Howden-Chapman P, Pierse N, Nicholls S, Gillespie-Bennett J, Viggers H, Cunningham M, Phipps R, Boulic M, Fjallstrom P, Free S: Effects of improved home heating on asthma in community dwelling children: randomised controlled trial. BMJ (Clinical research ed. 2008, 337: a1411-

    Google Scholar 

  102. Kattan M, Mitchell H, Eggleston P, Gergen P, Crain E, Redline S, Weiss K, Evans R, Kaslow R, Kercsmar C: Characteristics of inner-city children with asthma: the national cooperative inner-city asthma study. Pediatr Pulmonol. 1997, 24: 253-262.

    CAS  Google Scholar 

Download references


The authors would like to thank Kadin Tseng from the Boston University Office of Information Technology and Daniel Kamalic, Research Computing Manager at the Boston University College of Engineering for their help running the models at the Scientific Computer Facility.

The project described was supported by Award number R21ES017522 from the National Institute of Environmental Health Sciences. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Environmental Health Sciences or the National Institutes of Health.

Author information

Authors and Affiliations


Corresponding author

Correspondence to M Patricia Fabian.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

NKS, MS, JIL, GA, and MPF participated in the design of the study, and defined parameters used in the model. AG researched literature and helped define parameters used in the model. CR developed the initial discrete event modeling code. MPF wrote the modeling program and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Fabian, M.P., Stout, N.K., Adamkiewicz, G. et al. The effects of indoor environmental exposures on pediatric asthma: a discrete event simulation model. Environ Health 11, 66 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: