American Journal of Epidemiology Advance Access originally published online on February 17, 2007
American Journal of Epidemiology 2007 165(8):927-935; doi:10.1093/aje/kwk088
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ORIGINAL CONTRIBUTIONS |
Air Pollution Sources and Childhood Asthma Attacks in Cataño, Puerto Rico
1 Office of Surveillance and Biometrics, Center for Devices and Radiological Health, Food and Drug Administration, US Department of Health and Human Services, Rockville, MD
2 Department of Biostatistics, Bloomberg School of Public Health, Johns Hopkins University, Baltimore, MD
3 Center for Imaging Science, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD
4 Department of Psychiatry, School of Medicine, Johns Hopkins University, Baltimore, MD
5 Division of Environmental and Occupational Health, School of Medicine, Johns Hopkins University, Baltimore, MD
6 Department of Epidemiology, Bloomberg School of Public Health, Johns Hopkins University, Baltimore, MD
7 District of Columbia Cancer Registry, District of Columbia Health Department, Washington, DC
Correspondence to Dr. Nilsa I. Loyo-Berríos, Center for Devices and Radiological Health, Food and Drug Administration, US Department of Health and Human Services, 1350 Piccard Drive, HFZ-541, Rockville, MD 20850 (nilsa.loyo-berrrios{at}fda.hhs.gov).
Received for publication May 9, 2006. Accepted for publication October 9, 2006.
| ABSTRACT |
|---|
|
|
|---|
Asthma prevalence in the Cataño Air Basin of Puerto Rico is 27% for children aged 1314 years and 45% for children aged 56 years. There is concern that these rates are related to air pollution. The authors conducted a nested case-control study to evaluate whether proximity to air pollution point sources was associated with increased risk of asthma attacks. For 19972001, 1,382 asthma-related medical visits (International Classification of Diseases, Ninth Revision, codes 493 and 493.9) in children under 17 were identified through health insurance claims. Controls were children with no asthma attacks who were randomly selected from enrollees in two health insurance companies by incidence density sampling (1:5) and matched to cases on gender, age, insurance company, and event date. The distance from a point source to the subject's residence area represented a surrogate exposure measurement. Odds ratios for a 1-km decrease in distance were obtained by conditional logistic regression. Risk of asthma attack was associated with residing near a grain mill (odds ratio (OR) = 1.35), petroleum refinery (OR = 1.44), asphalt plant (OR = 1.23), or power plant (OR = 1.28) (all p's < 0.05). Residence near major air emissions sources (>100 tons/year) increased asthma attack risk by 108% (p < 0.05). These results showed that proximity to some air pollution sources is associated with increased risks of asthma attacks.
air pollution; asthma; child
Abbreviations: CI, confidence interval; OR, odds ratio
| INTRODUCTION |
|---|
|
|
|---|
Asthma is a common disease that affects people of all ages. It is a chronic inflammatory disorder of the lungs characterized by episodic and reversible symptoms of airflow obstruction and nonspecific airway hyperactivity (1, 2). In the United States, the self- or proxy-reported 12-month prevalence of asthma increased 73.9 percent from 1980 to 1996, when approximately 14.6 million persons (54.6 per 1,000 population) reported that they had had an asthma attack during the previous 12 months (3). In 1997, the National Health Interview Survey was redesigned, and the new estimate for 12-month prevalence was 40.7 per 1,000 population. However, an insufficient number of years with the new measures have passed to determine whether the trends in asthma are increasing or decreasing (3).
Ambient air pollution has been related to an increased risk of asthma attacks. Temporal associations have been shown between increases in the number of asthma attacks (emergency room visits/hospital admissions) and increases in levels of ambient air pollution in different parts of the world (412). Other studies have shown spatial associations, with asthma attacks being more common in areas that are closer to industrial parks, air pollution point sources, and high-traffic-volume roads (1317).
Ambient air pollution and its health effects have not been studied in detail in Puerto Rico. The Cataño Air Basin, located west of San Juan Bay, is known for its air pollution problems and high asthma prevalences among children aged 67 years (45 percent) and 1314 years (27 percent) (18). The area has been identified by the federal and local governments as a nonattainment area with regard to Environmental Protection Agency air quality standards for total suspended particulates and particulate matter less than 10 microns in diameter. The Cataño Air Basin includes the Cataño municipality (30,071 total population) and portions of four nearby municipalities: San Juan (434,374 total population), Guaynabo (100,053 total population), Bayamón (224,044 total population), and Toa Baja (94,085 total population) (19). The combined total population living in the Cataño Air Basin has been estimated at approximately 300,000 (20). Within its boundaries, there are numerous industrial point sources of air pollution, including a rum distillery, electric power plants, petroleum refineries, a sewage incinerator, sewage treatment plants, grain mills, and cement plants. Additionally, mobile sources include motor vehicle emissions from two highways and emissions from activities in a port zone. These air pollution sources are mostly concentrated in the Cataño municipality and its surrounding areas, generating concern about the observed high asthma prevalence rates being associated with the ambient air pollution.
Therefore, we conducted a study to evaluate whether there was an association between proximity to air pollution point sources and risk of asthma attacks in the Cataño municipality. We hypothesized that proximity to the air pollution point sources increased the risk of asthma attacks.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Data sources
Two health insurance companies that operate in Cataño (referred to as company A and company B) provided health and demographic data on the study participants. Data on emergency room visits, hospital admissions, and visits to physicians' offices were included in the study. Each company provides services to different sectors of the population (figure 1); company A provides services to people who can afford its private health insurance, and company B contracts with the government of Puerto Rico to provide services to persons living below the US poverty level.
|
Air pollution data and the locations of air pollution point sources were obtained from the Environmental Protection Agency (Envirofacts Data Warehouse (21) and AirData (22)). The National Climatic Data Center provided data on weather conditions (unpublished data), and census block maps were obtained from the US Census Bureau (24).
Study area and population
Puerto Rico is located in the Caribbean, and its weather conditions do not vary greatly across seasons of the year. The average temperature is 27°C (
80°F), and humidity averages 77 percent (National Climatic Data Center, unpublished data). The Cataño municipality is located in northern Puerto Rico, west of San Juan Bay; it is a small municipality with an area of 5.16 km2 (24).
The study population consisted of children younger than 17 years who were residents of the Cataño municipality from January 1997 through December 2001. By the year 2000, they represented 30 percent of Cataño's total population. The member rosters of the two health insurance companies included 69 percent of the child population through the year 2000 (figure 1).
Study design
We conducted a nested case-control study to evaluate the study hypothesis. This study design can be used when a cohort is available; it combines the features and advantages of both cohort and case-control designs, and because the cases and controls are selected from the same source cohort, selection bias is minimized (25). The study period consisted of 5 years of follow-up, starting in January 1997 and continuing through December 2001, and the two health insurance companies provided data on the cohort of children to be followed during that period. This cohort was dynamic; new children could enroll or children could leave the cohort, either because they became older than the age criteria or because their parents changed their health insurance company.
All of the cases observed during the follow-up period were included in the nested case-control study, and incidence density sampling was used to select the controls. Use of this sampling technique allowed for the selection of controls who were at risk and available in the cohort at the time an asthma attack was observed.
Definitions of cases, controls, and exposure
Cases.
An asthma attack was counted as one case. Multiple attacks in the same child were considered separate cases. An incident attack was defined as an attack that occurred on any day for which the previous day was free of events. Asthma attacks were identified through insurance claims that had an asthma diagnosis made by a physician on the basis of symptoms and were coded from the International Classification of Diseases, Ninth Revision, as 493.0 (extrinsic asthma) or 493.9 (asthma, unspecified). Because this study focused on childhood asthma, the code used to indicate late onset of asthma (493.1) was not included. A total of 1,382 cases were identified during the study period.
Controls.
Three control groups were used in the analysis. The first control group (control group 1) was defined as children who were free of an asthma attack on the day a case was observed. This group included asthmatic children and children who could have had other respiratory conditions (e.g., any respiratory condition other than asthma, such as "catarrh, catarrhal (inflammation)" or "croup") on the day of the case. The second control group (control group 2) was defined similarly to the first group, except that children with other respiratory conditions on the day of the case attack were not included. Finally, the third control group (control group 3) consisted of children who never had an asthma attack or other respiratory condition during the 5 years of follow-up. Controls were matched to the cases on gender, age, insurance company, and date of the asthma event. Incidence density sampling was used to select five controls per case. A total of 6,910 controls were included in each control group.
Exposure.
The physical distance from the air pollution point sources to the center of the census block in which the child's residence was located was used as a surrogate measure for average exposure to ambient air pollution. A global positioning system was used to obtain the latitudes and longitudes of the center of each census block. Then, distance was calculated with R software (26) using the Haversine Formula (equations 15), which takes into account the curvature of the earth (27).
![]() | (1) |
|
| (2) |
![]() | (3) |
|
| (4) |
|
| (5) |
Statistical procedures
The geographic distribution of the point sources was graphed using census maps in MapInfo Professional, version 7 (24).
Conditional logistic regression models were used to assess the association between distance from the residence to point sources and risk of asthma attack. R software, version 1.6.2, was used to conduct the analyses (26). The risk of asthma attack was estimated for a 1-km decrease in the distance to air pollution point sources. Three measurements of distance were used.
First, we used data from the Envirofacts Data Warehouse (21). The Envirofacts Data Warehouse classifies air pollution point sources as being major (
100 tons/year) or minor (<100 tons/year) sources of air emissions (based on the total air emissions per year). That classification was used to estimate the mean distance from the residence to major or minor air emissions sources for each child.
Second, we also classified point sources according to types of industries that have been related to respiratory problems (power plants, grain mills, refineries, distilleries, asphalt plants, and quarries). We conducted a subanalysis for this exposure measurement using distance categories.
Third, we developed a corrected distance measurement to take wind direction into account (equation 6) (28).
![]() | (6) |
|
| RESULTS |
|---|
|
|
|---|
Exploratory analysis
Figure 3 presents the geographic distribution of the point sources, as well as their classification by air emission type. There were 67 sources identified in the area, from which 20 were located within the municipality and 43 in the surrounding areas. There were no latitude and longitude data for four industries. As reported in the Envirofacts database (21), 17 sources were classified as major air emissions sources and 44 were identified as minor air emissions sources. Information on emission type was not available for six industries; those are represented as minor sources on the map and were excluded from risk estimation. Most of the major air emissions sources were located southeast of town. Two power plants, seven grain mills, five petroleum refineries, two asphalt plants, one quarry, and one rum distillery were identified in the area.
|
Multivariate analysis
Odds ratios were estimated for a 1-km decrease in the mean distance to major or minor air pollution sources. Table 1 presents the results. When sources were evaluated independently as major or minor air emissions sources, a 1-km decrease in the distance to major sources was associated with a 69 percent increase in risk (odds ratio (OR) = 1.69, 95 percent confidence interval (CI): 1.50, 1.91), and minor sources were associated with a risk increase of 47 percent (OR = 1.47, 95 percent CI: 1.32, 1.63). However, after controlling for type of air emissions, minor sources were no longer important, whereas major sources presented a 108 percent increased risk (OR = 2.08, 95 percent CI: 1.61, 2.70) (table 1). The observed associations were similar with the comparisons for control groups 2 and 3.
|
Table 2 presents the results for mean distance by type of industry. With the exception of power plant 1 (OR = 0.76, 95 percent CI: 0.72, 0.81) and the rum distillery (OR = 0.78, 95 percent CI: 0.73, 0.84), all of the sources were associated with increased risks of asthma attack; odds ratios ranged from 1.28 to 1.44, and the 95 percent confidence intervals excluded 1. The associations were similar for the comparisons with control groups 2 and 3. No effect modification was observed by wind speed or direction, temperature, humidity, or air pollution (data not shown). Figure 4 presents the results for the subanalysis by distance category and type of industry. Children living within 2,000 m of the power plants and the rum distillery did not have increased risk. Interestingly, the risk estimates for the two power plants differed; the plant located north of town (power plant 1) never posed increased risk, whereas for the plant located southeast of town (power plant 2), the risk increased more rapidly and was higher for children living within 2,0013,000 m. Furthermore, for all of the sources located to the south and southeast, proximity was associated with increased risks of asthma attack, and the risks decreased with distance.
|
|
Table 3 presents the results for the corrected distance models. The association observed for the grain mills, the petroleum refineries, and the asphalt plants remained after wind correction (tables 2 and 3). However, the corrected distance model showed a significant 12 percent increased risk for the power plants (OR = 1.12, 95 percent CI: 1.03, 1.22) after wind correction. Similar results were observed for the comparisons with control groups 2 and 3.
|
| DISCUSSION |
|---|
|
|
|---|
The study hypothesis postulated that proximity to industrial point sources of air pollution could increase the risk of asthma attacks. Various measurements of distance were used as indicators of average exposure to ambient air pollution. Results appeared to be consistent: Proximity to some point sources was related to increased risk of asthma attacks. The highest risk estimates were observed for point sources located in the south and southeast.
The results of our study are supported by reports in the scientific literature suggesting that proximity to air pollution point sources is related to increased risk of respiratory health problems (13, 1517). A study conducted in northeastern England showed that persons living closest to air pollution sources had higher risks of having had more than 12 asthma attacks in the previous year (17). In another study in West Virginia, investigators found a higher risk of respiratory symptoms among children attending schools in a valley where industrial sources were located (13).
In the present study, we not only evaluated distance as a risk factor but also identified some other aspects that can influence the spatial association between air pollution and asthma attack risk. These factors are the amount of emissions and the location of the point sources. The data showed that major air emissions sources strongly influenced the risk estimate, with a 108 percent increase in risk. These data suggest that if those 17 sources identified as major air polluters by weight of output were able to decrease their annual emissions, the risk of asthma attacks in the Cataño Air Basin could be minimized. Additionally, the location of the sources proved to be another important factor determining risk. The sources that were located north of town did not present increased risks. This was not the case for sources located to the south and southeast, all of which posed increased risks. Furthermore, this analysis identified two sources of the same type (power plant 1 and power plant 2) that presented opposite risks. The only difference between these two sources was their location; the power plant located in the southeast presented a risk that increased by 28 percent for each 1-km decrease in distance, whereas the power plant located in the north did not present increased risk.
Wind direction could be the reason for the high risk observed in the southeast. In this area, the wind mostly blows from the northeast, but there are times when it blows from the east. It is possible that emissions from the sources located to the southeast are carried towards the population living in the southeastern part of town. One of the strengths of this study is that a measurement of distance was developed to account for the wind direction on the day of the event (an asthma attack) (figure 2). This model assumes that same-day exposure will trigger an attack. After wind correction, the results changed only for the power plants, which presented a 12 percent increase in the risk associated with a 1-km decrease in distance.
Other factors that may also influence the risk are weather conditions (temperature and humidity), parents' occupations, and where children go to school. Weather conditions in Puerto Rico remain similar across all seasons of the year. Therefore, weather factors were considered to be relatively constant over the 5 years of study. In fact, analyses stratified by temperature and humidity did not show effect modification for the observed association between distance and increased risk of asthma attack. On the other hand, parents' occupations may affect the risk estimate if their working environments involve exposure to contaminants that can trigger asthma attacks and if those contaminants are brought home on the parents' clothing. This study did not include data on parents' occupations; further research examining indoor air quality and related factors is warranted. Additionally, where children go to school is an important factor to consider, since children spend most of their time at school. This study considered residence area rather than school location, and information on school location was not available. However, the study area is a small municipality, the school system in Puerto Rico assigns children to schools by residence area, and the norm is for children to go to schools within their residence area. Very rarely will a student go to a school that is either outside of his/her residential area or outside the municipality. Therefore, in this study, we considered residence area an appropriate indirect measurement for average exposure to outdoor air pollution.
Another aspect to consider when interpreting these results is that health insurance claims data capture only those children who need to seek health services. In this study, it was possible that children who had an asthma attack but had medications at home and did not seek health services for their attack were not captured as cases. Data on medications were not available for this study, and missing cases or misclassification of cases as controls may have resulted in lower odds ratios. We addressed this issue by using three different control groups: Control groups 1 and 2 included asthmatic children with a history of attacks, and control group 3 excluded asthmatic children. Control groups 1 and 2 should have included children with a history of asthma attacks and medications at home. The risk estimates for these two groups will tend to have been lower if cases were missed or if case children were misclassified as controls because of self-medication. Note that self-medication may not fully treat severe attacks, so some children with medications at home may still have needed to seek health services and therefore would have been captured as cases. Note also that the observed odds ratios showed increased risk, and they could have been underestimated. Future research considering medication use is warranted.
Another limiting aspect of these data was the diagnosis of asthma for children younger than 1 year of age. Asthma diagnosis for this age group is less specific, which could lead to misclassification of controls as cases. In this study, only 3.9 percent of the cases were younger than 1 year (data not shown); therefore, we would expect the effect of such misclassification to have been minimal.
This study had several other limitations. Correlation between locations of point sources prevented the identification of a specific source as the sole cause of increased asthma attack risk. When distance to an independent source is used, the estimated risk carries the risk of nearby sources. In this study, we addressed this limitation by classifying point sources as major or minor air emission sources, irrespective of their location. After controlling for emission type, increased risk was associated with proximity to point sources with major air emissions. Most of these major polluters were located south and southeast of town.
Finally, there was a lack of data on smokestack height, temperature of emissions, and emissions rate, which could have been used to determine how the emissions were distributed in the area and affecting the community. Nevertheless, this study showed that when these data are not available, valid estimates of risk can be obtained. Not only can this method be used in communities with similar health and air pollution problems, it can also be used to examine other health outcomes and exposures.
In future research in Cataño, investigators should consider evaluating how emissions are distributed and how they affect the community, perhaps using the Gaussian Plume model (29). Special attention should be given to the two power plants. Evaluations of the point sources located south and southeast of town are recommended, especially those classified as major air emissions sources. Minimizing emissions from these sources could reduce the occurrence of asthma attacks in this population.
| ACKNOWLEDGMENTS |
|---|
This research was conducted to fulfill the requirements for a doctoral degree in epidemiology and was funded by the Department of Epidemiology at the Johns Hopkins Bloomberg School of Public Health (Baltimore, Maryland). Dr. Nilsa I. Loyo-Berríos currently works for the Food and Drug Administration and prepared this paper under the Professional Development Program.
The authors thank Dr. Luis Marini, Dr. Benjamin Santiago, Carlos Morel, and Dr. Alejandro Franco, from the health insurance companies that provided the data on asthma. They also acknowledge Dr. Erika Avila-Tang for her help with data abstraction.
Conflict of interest: none declared.
| References |
|---|
|
|
|---|
- Kemp JP. (1998) Comprehensive asthma management: guidelines for clinicians. J Asthma 35:60120.[ISI][Medline]
- Mannino DM, Homa DM, Pertowski CA, et al. (1998) Surveillance for asthmaUnited States, 1960 1995. MMWR CDC Surveill Summ 47:127.[Medline]
- Mannino DM, Homa DM, Akinbami LJ, et al. (2002) Surveillance for asthmaUnited States, 1980 1999. MMWR Surveill Summ 51:113.[Medline]
- Castellsague J, Sunyer J, Saez M, et al. (1995) Short-term association between air pollution and emergency room visits for asthma in Barcelona. Thorax 50:10516.[Abstract]
- Norris G, YoungPong SN, Koenig JQ, et al. (1999) An association between fine particles and asthma emergency department visits for children in Seattle. Environ Health Perspect 107:48993.[ISI][Medline]
- Schouten JP, Vonk JM, de Graaf A. (1996) Short term effects of air pollution on emergency hospital admissions for respiratory disease: results of the APHEA project in two major cities in The Netherlands, 1977 89. J Epidemiol Community Health 50:suppl 1, s229.
- Burnett RT, Dales RE, Raizenne ME, et al. (1994) Effects of low ambient levels of ozone and sulfates on the frequency of respiratory admissions to Ontario hospitals. Environ Res 65:17294.[Medline]
- Burnett RT, Dales R, Krewski D, et al. (1995) Associations between ambient particulate sulfate and admissions to Ontario hospitals for cardiac and respiratory diseases. Am J Epidemiol 142:1522.
[Abstract/Free Full Text] - Burnett RT, Smith-Doiron M, Stieb D, et al. (1999) Effects of particulate and gaseous air pollution on cardiorespiratory hospitalizations. Arch Environ Health 54:1309.[ISI][Medline]
- Burnett RT, Smith-Doiron M, Stieb D, et al. (2001) Association between ozone and hospitalization for acute respiratory diseases in children less than 2 years of age. Am J Epidemiol 153:44452.
[Abstract/Free Full Text] - Romieu I, Meneses F, Ruiz S, et al. (1996) Effects of air pollution on the respiratory health of asthmatic children living in Mexico City. Am J Respir Crit Care Med 154:3007.[Abstract]
- Yu O, Sheppard L, Lumley T, et al. (2000) Effects of ambient air pollution on symptoms of asthma in Seattle-area children enrolled in the CAMP study. Environ Health Perspect 108:120914.[ISI][Medline]
- Ware JH, Spengler JD, Neas LM, et al. (1993) Respiratory and irritant health effects of ambient volatile organic compounds: The Kanawha County Health Study. Am J Epidemiol 137:1287301.
[Abstract/Free Full Text] - Anto JM and Sunyer J. (1986) A point-source asthma outbreak. Lancet 1:9003.[Medline]
- Anto JM and Sunyer J. (1990) Epidemiologic studies of asthma epidemics in Barcelona. Chest 98:suppl, 185S90S.
- English PN, Raymond PN, Scalf R, et al. (1999) Examining associations between childhood asthma and traffic flow using a geographic information system. Environ Health Perspect 107:7617.[ISI][Medline]
- Bhopal RS, Moffatt S, Pless-Mulloli T, et al. (1998) Does living near a constellation of petrochemical, steel, and other industries impair health? Occup Environ Med 55:81222.[Abstract]
- Environmental Health Secretariat, Puerto Rico Department of Health. (1996) Prevalence of respiratory and allergic problems among school-aged children in the Cataño Air Basin: final report. (Puerto Rico Department of Health, San Juan, Puerto Rico).
- Bureau of the Census, US Department of Commerce. Fact sheet: United States, 2000. In: Population Finder. (Database). (2000) (Bureau of the Census, Washington, DC) (http://factfinder.census.gov/servlet/SAFFFacts?_submenuId=factsheet_1&_sse=on).
- Environmental Health Secretariat, Puerto Rico Department of Health. Cataño workshop final report: The Cataño Air Basin Asthma Prevalence Study, 1996. (1997) (Puerto Rico Department of Health, San Juan, Puerto Rico).
- Environmental Protection Agency. (2000) Locational Reference Tables (LRT). In: Envirofacts Data Warehouse. (Database). (Environmental Protection Agency, Washington, DC) (http://www.epa.gov/enviro/html/locational/lrt/).
- Environmental Protection Agency. (2002) AirData: Access to Air Pollution Data. (Database). (Environmental Protection Agency, Washington, DC) (http://www.epa.gov/air/data/index.html).
- Bureau of the Census, US Department of Commerce. (2002) Census 2000 block maps. (Bureau of the Census, Washington, DC) (http://www.census.gov/geo/www/maps/descriptwindows/census_blockmaps.htm).
- MapInfo Corporation. MapInfo Professional, version 7. (2002) (MapInfo Corporation, Troy, NY) (http://www.mapinfo.com/).
- Szklo M and Nieto FJ. (2000) Epidemiology: beyond the basics. (Aspen Publishers, Inc, Gaithersburg, MD).
- R Core Group. R software. (19982003) (R Foundation for Statistical Computing, c/o Vienna University of Economics, Vienna, Austria) (http://www.r-project.org/).
- Sinnott RW. (1984) Virtues of the Haversine. Sky Telescope 68:159.
- Gellert W, Gottwald S, Hellwich M, et al. (1989) The VNR concise encyclopedia of mathematics. 2nd ed (Van Nostrand Reinhold Company, New York, NY).
- Curtiss PS and Rabl A. (1996) Impacts of air pollution: general relationships and site dependence. Atmos Environ 30:333147.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||






), Cataño municipality, Puerto Rico, 19972001