Skip to main content

Association of PM2.5 concentration with health center outpatient visits for respiratory diseases of children under 5 years old in Lima, Peru

  • The Correction to this article has been published in Environmental Health 2020 19:11



Lima is one of the more polluted cities in Latin America. High levels of PM2.5 have been shown to increase health center outpatient visits of respiratory diseases.


Health center outpatient visits for children < 5 years for childhood respiratory disease (acute lower respiratory infections (ALRI), pneumonia and acute bronchiolitis/asthma) from 498 public clinics in Lima were available on a weekly basis from 2011 to 2015 from Peru’s Ministry of Health (MINSA). The association between the average weekly concentrations of PM2.5 was evaluated in relation to the number of weekly health center outpatient visits for children. Weekly PM2.5 values were estimated using a recently developed model that combined data observed from ground monitors, with data from space satellite and meteorology. Ground monitoring data came from 10 fixed stations of the Peruvian National Service of Meteorology and Hydrology (SENAMHI) and from 6 mobile stations located in San Juan de Miraflores by Johns Hopkins University. We conducted a time-series analysis using a negative binomial model.


We found a significant association between exposure to PM2.5 and all three types of respiratory diseases, across all age groups. For an interquartile increase in PM2.5, we found an increase of 6% for acute lower respiratory infections, an increase of 16–19% for pneumonia, and an increase of 10% for acute bronchiolitis / asthma.


Higher emissions of environmental pollutants such as PM2,5 could be a trigger for the increase of health center outpatients visits for respiratory diseases (ALRI, pneumonia and asthma), which are themselves risk factors for mortality for children in Lima province, Peru.

Peer Review reports


A World Health Organization (WHO) report regarding global outdoor air pollution in 2014 noted that Lima, the capital of Peru, was one of more polluted cities in the Americas [1]. Average PM2.5 levels during 2014–2015 were 26 μg/m3 [2]. Ambient air pollution has been associated with respiratory diseases in children [3]. Environmental air pollution is one of the causes of mortality and morbidity, due to cardiovascular diseases, acute respiratory infections, pneumonias and acute bronchiolitis/asthma [4, 5]. Respiratory diseases are among the leading causes of death in the world [6].

According to the WHO, acute lower respiratory infections (ALRI) cause the death of 4.3 million children under 5 years old, which represents 30% of the total annual deaths of children in this age group [6]. Also ALRI are the leading cause of premature death in Peru, approximately 222 child deaths per year for every 100,000 live births [6].

PM2.5 are fine particles with aerodynamic diameter of 2.5 μm or less, which are emitted from a large variety of sources including automotive vehicles, industry, power generation and engine combustion [7, 8]. Some prior studies have found an association between PM2.5 and increased health center outpatient visits for respiratory diseases [9]. It has also been shown that PM2.5 levels have a greater effect on admission to outpatient health clinics due to respiratory problems than do PM10 levels [4]. Childhood respiratory diseases may be partially preventable with better control of environmental contaminants, such as PM2.5 [10].

We evaluated the association between the average weekly concentration of PM2.5, and the number of weekly health center outpatients visits for ALRI, pneumonia and acute bronchiolitis/asthma in children under 5 years old in public health centers in Lima province, Peru between 2011 and 2015.


We conducted a time-series analysis, with data related to health center outpatient visits due to childhood respiratory diseases, including acute lower respiratory inflammatory, pneumonia, asthma. Children with respiratory diseases were divided according to respiratory diseases and age groups: ALRI (acute lower respiratory disease, J00-J11.1 ICD-10 (International Classification of Diseases 10th)) from 0 to 2 months, 2 to 11 months and 1 to 5 years; pneumonia (J12.0-J18.9 ICD-10) from 2 to 11 months and 1 to 4 years; and acute bronchiolitis/asthma (J21.0-J21.9 / J44.8–46.0 ICD-10)(acute bronchiolitis (J21) and asthma (J44.8-J46) are similar diseases and may have similar causes). Data were grouped by age, with categories varying slightly by disease category. The age and disease groupings are those used by Peru’s Ministry of Health (MINSA), and we received the data with these groupings.

Our study included all child outpatients who lived in Lima province, considering all the districts (43 districts) of the Lima province. In the Lima province, there are 498 MINSA health centers which receive outpatients. The number of health center outpatient visits was obtained for January to December of 2011–2015 (n = 3,099,438); cases were available only grouped by week, so weekly cases were our unit of analysis. The outcome data were obtained from the National Center for Epidemiology, Prevention and Disease Control, a part of the Peruvian Ministry of Health (MINSA).

PM2.5 values were estimated using a model that combined data observed from ground monitors, with space satellite and meteorological data, to estimate daily PM2.5 at a 1 km2 resolution. The ground monitoring was carried out at 10 fixed stations of the National Service of Meteorology and Hydrology of Peru (SENAMHI) and from 6 mobile stations located in San Juan de Miraflores, collected by Johns Hopkins University [7]. This model was a good predictor of the observed ground monitoring data (R-square = 0.70).

Daily PM2.5 estimates by district, weighted by population density, were averaged to get weekly means for each district where a patient lived. District of residence was available from health center records. The weekly mean temperature, weekly mean relative air humidity, seasons (summer, autumn, fall and spring), years, indicator variables for districts, all of which may act as confounders, were included as variables in the models; the meteorological data were obtained from SENAHMI [6]. The analysis includes 11,050 observations (52 weekly visits * 5 years * 43 districts), with a 1.2% loss of observations due to missing data.

The statistical analysis was first done via a negative binomial model (the Poisson model was over-dispersed) with week-long lags of 0 to 3 weeks. Goodness of fit was evaluated using the Akaike Information Criterion (AIC). The rate ratio for health center outpatient visits was estimated for an increase in the interquartile range (IQR) of PM2.5, the increase from the 25th to the 75th percentiles (which was 7.1 μg/m3 during the years considered). Also, we analyzed PM2.5 as a categorical variable in quintiles (Q): 1stQ: < 15.64 μg/m3; 2ndQ: 15.64–17.48 μg/m3; 3rdQ:17.49–19.71 μg/m3; 4thQ: 19.72–25.08 μg/m3 and 5thQ: 25.09–48.62 μg/m3.

The inclusion of a dummy variable for district in the model essentially stratified on district, so overall results are a weighted average of strata-specific estimates. The use of a dummy variable in the model for district is a more efficient way to take district level effects into account, rather than to create an observation for each district for each week/year/season, which would increase the size of data set 43 times, and require a meta-analysis to combined district-specific estimates. We do present some district-specific results for the large districts, in Additional file 1.

Analysis was conducted using the statistical software Stata 12 (Stata, Inc., Texas, USA) and Excel (Microsoft Office Excel 2007; Microsoft Corporation). A p-value < 0.05 was considered as significant. The data were obtained through agreements with the MINSA and SENAMHI. As this is a record based study with no contact with patients, and the data were anonymized, it was not submitted to an Internal Review Board (Ethics Committee).


During the study period, there were 3,099,438 cases of health center outpatient visits for respiratory diseases of children under 5 years old in Lima province, Peru between 2011 and 2015. Table 1 gives the mean weekly number of visits for the various outcomes by respiratory diseases according age group, as well as the average weekly PM2.5 across all years, and the average weekly temperature, and relative humidity values across all years. The population-weighted average PM2.5 estimated for Lima province, across all districts and years, was 20.5 μg/m3 (SE = 6.3) (Table 1).

Table 1 Average weekly values across the study period for outcome and predictor variables, Lima-Peru, 2011–2015

In 20% of the weeks considered, the weekly mean concentrations of PM2.5 presented values above the limit daily recommended by the WHO (25 μg/m3) and almost 100% presented values above the annual recommended limit (10 μg/m3) [11]. In relation to the Peruvian Ministry of the Environmental (MINAM) recommendations, 20% of the weeks considered were above the annual limit (25 μg/m3), although none surpassed the daily limit (50 μg/m3) (Fig. 1).

Fig. 1

Modeling of PM2.5 maximum concentration daily in Lima province, Peru, 2011–2015. = MINAM Air Quality Guidelines 24 h standard (50 μg/m3)  — = WHO Air Quality Guidelines 24 h standard (25 μg/m3) and MINAM Air Quality Guidelines annual (25 μg/m3 =WHO Air Quality Guidelines annual (10 μg/m3) Units on x axis are weekly concentration of PM2.5 (μg/m3)

Figure 1 shows the distribution daily of concentrations of PM2.5, in Lima province, Peru, 2011–2015.

Lag (0) fit better than the other lags according to the AIC for all 3 types of respiratory diseases), and was used in the regressions.

The rate ratios (RR) for outpatient visits for an IQR increase in PM2.5 are shown in Table 2. We found a significant association between exposures to PM2.5 for all three types of respiratory diseases, and across all age groups. All ALRI increases 6% per IQR, while asthma increases 10% and pneumonia increases 17%.

Table 2 Rate ratios for respiratory diseases associated with each interquartile range increase in PM2.5 concentration in different age group in Lima province, Peru, 2011–2015a

Table 3 shows analyses by quintile of PM2.5. We found consistent monotonic increases in rate ratios with increasing PM2.5, with increases of more than 10% in rates of clinic visits in the highest quintile (25.1–48.6 μg/m3) in all disease groups.

Table 3 Relationship between respiratory diseases with PM2.5 quintiles in children under 5 years in Lima-Perua


Lima province has 8565.213 inhabitants according to the 2017 National Census [12]. The districts of Lima province have a large childhood population exposed to levels of air pollution well above the recommended annual WHO Air Quality Guidelines for PM2.5 (10 μg/m3) [11].The association between PM2.5 concentrations and health center outpatients visits for respiratory diseases in our study was evident at concentrations which are in the range of the MINAM’s annual permissible level in Lima (25 μg/m3) [13, 14]. Our data suggest that it is necessary to take corrective measures to reduce the effects of environmental pollutants, such as PM2.5.

We have studied health center outpatient visits, which are generally less severe than visits to emergency rooms, and we have focused on children under 5. Our study found a significant association between an increase of 7.1 μg/m3 in PM2.5 with increased outpatient visits for respiratory disease, with effects for all types of respiratory disease (ALRI, pneumonia and acute bronchiolitis/asthma).

Most studies in Lima, Peru about contamination by PM2.5 are restricted to describing exposure, without any association with disease [2, 8, 14]. On exception is the paper by Tapia et al. (2019), which found that for each interquartile range (IQR) increase in PM2.5, respiratory disease, emergency room visits in Lima increased 3% (95% CI: 1.02–1.04%) for those under 18 years [15]. In our study we found a 6% (95% CI: 1.05–1.08%) increase in respiratory clinic visits children under 5 years. It is possible that children under 5 are more susceptible to the effects of air pollution. On the other hand, a recent meta-analysis of 16 time-series studies of hospital admissions in children younger than 5 years, found a respiratory disease increase of 2.7% (95 CI = 0.9–7.7%) per 10 μg / m3 increase in PM2. 5 [16], also lower than our estimate.

Our study has some limitations. We don’t know how accurate our local estimates of PM2.5 are, although the model [7] had a very good correlation with daily observed levels on the ground (r-square = 0.70). Another possible limitation is that, although data from health center outpatients visits were obtained from an official source (MINAM), they may contain diagnostic errors, and they were available only for epidemiological week, rather than daily. Also we had no data on disease severity, although it is likely that cases were not severe ones, as they did not go to the emergency room.


Higher emissions of environmental pollutants such as PM2,5 could be a trigger for the increase of health center outpatients visits for respiratory diseases (ALRI, pneumonia and asthma), which are themselves risk factors for mortality for children in Lima province, Peru.

Availability of data and materials

Data from this study is available upon request.

Change history

  • 30 January 2020

    The original version of this article [1], published on 15 January 2020, contained incorrect name of the co- author. In this Correction the affected part of the article is shown.



Akaike information criterion


Acute lower respiratory infections


Confidence interval


International Classification of Diseases 10th


Interquartile range


Peruvian Ministry of the Environmental


Peru’s Ministry of Health






Rate ratio


Standard error


Peruvian National Service of Meteorology and Hydrology


World Health Organization


  1. 1.

    Hansel NN, Romero KM, Pollard SL, Bose S, Psoter KJ, Underhill L, et al. Ambient Air Pollution Adversely Impacts Various Domains of Asthma Morbidity among Peruvian Children. Ann Am Thorac Soc. 2018;16(3):1–32 Available from: [cited 2019 Mar 29].

    Google Scholar 

  2. 2.

    Silva J, Rojas J, Norabuena M, et al. Particulate matter levels in a South American megacity: the metropolitan area of Lima-Callao, Peru. Environ Monit Assess. 2017;189:635.

  3. 3.

    Thi N, Nhung T, Amini H, Schindler C, Künzli N. Short-term association between ambient air pollution and pneumonia in children : a systematic review and meta-analysis of time-series and. Environ Pollut. 2017;230:1000–8.

    CAS  Article  Google Scholar 

  4. 4.

    Koutrakis P, Sax SN, Sarnat JA, Coull B, Demokritou P, Oyola P, et al. Journal of the Air &amp; Waste Management Association Analysis of PM 10 , PM 2.5 , and PM 2.5–10 Concentrations in Santiago, Chile, from 1989 to 2001 Analysis of PM 10 , PM 2.5 , and PM 2.5–10 Concentrations in Santiago, Chile, from 1989 to 2001. J Air Waste Manage Assoc. 2005;55(3):342–51 Available from: [cited 2019 Apr 25].

    CAS  Article  Google Scholar 

  5. 5.

    Gouveia N, Junger WL, Romieu I, Cifuentes LA, de Leon AP, Vera J, et al. Effects of air pollution on infant and children respiratory mortality in four large Latin-American cities. Environ Pollut. 2018;232:385–91 Available from: 1. [cited 2019 Apr 16].

    CAS  Article  Google Scholar 

  6. 6.

    Foro de las Sociedades Respiratorias Internacionales. El impacto mundial de la Enfermedad Respiratoria- Segunda Edición. 2017. Available from: [cited 2019 May 16]

  7. 7.

    Vu BN, Sánchez O, Bi J, Xiao Q, Hansel NN, Checkley W, et al. Developing an advanced PM 2.5 exposure model in Lima, Peru. Remote Sens. 2019;11(6):1–18.

    Article  Google Scholar 

  8. 8.

    Tapia V, Carbajal L, Vásquez V, Espinoza R, Vásquez-Velásquez C, Steenland K, et al. Reordenamiento vehicular y contaminación ambiental por material particulado (2,5 y 10), dióxido de azufre y dióxido de nitrógeno en Lima Metropolitana, Perú. Rev Peru Med Exp Salud Publica. 2018;35(2):190 Available from: [cited 2019 Apr 16].

    Article  Google Scholar 

  9. 9.

    Slama A, Śliwczyński A, Woźnica J, et al. Impact of air pollution on hospital admissions with a focus on respiratory diseases: a time-series multi-city analysis. Environ Sci Pollut Res. 2019;26:16998–17009.

    CAS  Article  Google Scholar 

  10. 10.

    Barraza F, Lambert F, Jorquera H, Villalobos AM, Gallardo L. Temporal evolution of main ambient PM2.5 sources in Santiago, Chile, from 1998 to 2012. Atmos Chem Phys. 2017;17(16):10093–107.

    CAS  Article  Google Scholar 

  11. 11.

    World Health Organization. Guías de calidad del aire de la OMS relativas al material particulado, el ozono, el dióxido de nitrógeno y el dióxido de azufre. Actualización mundial 2005 [Internet]. 2006. Available from:;jsessionid=7B0591E84285605728DE66380E4F3DB6?sequence=1.

  12. 12.

    Instituto Nacional de Estadística e Informática. Perú. Resultados Definitivos [Internet]. Vol. Tomo V, Instituto Nacional de Estadística e Informática. 2018. Available from:

  13. 13.

    Ministerio del Ambiente-Perú. Decreto Supremo No 003-2017-MINAM [Internet]. Normas Legales El Peruano. 2017. Available from:

  14. 14.

    Pacsi Valdivia SA. Analisis temporal y espacial de la calidad del aire determinado por material particulado PM 10 y PM 2,5 en Lima Metropolitana. An Científicos. 2016;77(2):273–83 Available from:

    Article  Google Scholar 

  15. 15.

    Tapia V, Steenland K, Sarnat SE, et al. Time-series analysis of ambient PM2.5 and cardiorespiratory emergency room visits in Lima, Peru during 2010–2016. J Expo Sci Environ Epidemiol. 2019.

  16. 16.

    Requia WJ, Adams MD, Arain A, Papatheodorou S, Koutrakis P, Mahmoud M. Global Association of Air Pollution and Cardiorespiratory Diseases: A Systematic Review, Meta-Analysis, and Investigation of Modifier Variables. Am J Public Heal. 2018;108(S2):S123–30 Available from: [cited 2019 Jun 16].

    Article  Google Scholar 

Download references


Research reported in this publication was supported by the NIH Fogarty International Center, National Institutes of Environmental Health Sciences (NIEHS) National Cancer Institute, National Institute for Occupational Safety and Health, and the NIH under Award Number U01 TW0101 07.

The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.


The present study was founded by the National Institutes of Health (Fogarty Program) [Grant U01TW010107, 1/2 Regional GEOHealth Hub centered in Peru].

Author information




The study was conceived by GFG and VV; VT, VV and LO participated in data collection; JED, VT, BV and KS participated in data analysis. JED and KS wrote the first draft of the manuscript. All authors reviewed the final version and agreed on its publication.

Corresponding author

Correspondence to Jennifer Estefanía Davila Cordova.

Ethics declarations

Ethics approval and consent to participate

The study has the approval (Certificate 84–03-16) of the University Cayetano Heredia ethics committee.

Consent for publication

All authors agree with the publication of the research.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original version of this article was revised: incorrect name of the co-author has been updated.

Additional file

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Davila Cordova, J., Tapia Aguirre, V., Vasquez Apestegui, V. et al. Association of PM2.5 concentration with health center outpatient visits for respiratory diseases of children under 5 years old in Lima, Peru. Environ Health 19, 7 (2020).

Download citation


  • Health center outpatient visits
  • PM2.5
  • Respiratory diseases
  • Air pollution
  • Children