ORIGINAL RESEARCH article

Front. Public Health, 02 September 2025

Sec. Environmental Health and Exposome

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1593285

Association between daily variations in the levels of atmospheric O3, PM2.5, and NO2 and the frequency of hospital visits due to respiratory diseases and hypertension in Mexico City using Generalized Additive Mixed Models

  • 1. Instituto de Ciencias de la Atmósfera y Cambio Climático, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico

  • 2. Department of Soil Sciences, Colegio de Postgraduados, Montecillo, Mexico

  • 3. Laboratorio Nacional de Ciencias de la Sostenibilidad, Instituto de Ecología, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico

  • 4. Facultad de Medicina, Universidad Anáhuac México Sur, Mexico City, Mexico

  • 5. Department of Geography and Environmental Sciences, Northumbria University, Newcastle upon Tyne, United Kingdom

  • 6. Unidad de Investigación Médica en Epidemiología Clínica, UMAE Hospital de Pediatría Siglo XXI, Instituto Mexicano del Seguro Social (IMSS), Mexico City, Mexico

  • 7. Dirección de Monitoreo de Calidad del Aire, Secretaria del Medio Ambiente, Gobierno de la Ciudad de Mexico, Mexico City, Mexico

  • 8. Coordinación de Vigilancia Epidemiológica, Instituto Mexicano del Seguro Social, Mexico City, Mexico

  • 9. Unidad de Educación, Investigación y Políticas de Salud, Instituto Mexicano del Seguro Social, Ciudad de Mexico, Mexico

  • 10. Population Health Sciences Institute, Newcastle University, Newcastle upon Tyne, United Kingdom

  • 11. Laboratorio de Ecología, Unidad de Biología de la Conservación, Parque Científico y Tecnológico de Yucatán, Universidad Nacional Autónoma de México (UNAM), Mérida, Mexico

  • 12. Análisis y Modelación de Calidad del Aire, Secretaría de Medio ambiente de la Ciudad de México, Mexico City, Mexico

  • 13. Department of Mathematics, Physics and Electrical Engineering, Northumbria University, Newcastle upon Tyne, United Kingdom

  • 14. Genómica del Cáncer, Instituto Nacional de Medicina Genómica, Mexico City, Mexico

Abstract

Background:

Environmental pollution is a significant public health issue in Mexico City. Patients with respiratory or cardiovascular diseases such as asthma or hypertension often experience exacerbations triggered by environmental factors. This retrospective time series epidemiological study analyzed the association between daily air pollution levels and the numbers of medical visits for exacerbations of cardiorespiratory diseases.

Methods:

Records from primary, secondary, and tertiary hospitals of the Mexican Institute of Social Security (IMSS) were reviewed between 2017 and 2019. Air quality data, including PM2.5, PM10, O3, and NO2 concentrations were collected hourly from air quality monitoring stations at fixed sites. To fit the models and

take into account temporal autocorrelation and the complex non-linear relationships between pollutants and cardiorespiratory diseases, Generalized Additive Mixed Models (GAMM) were applied.

Results:

PM2.5, O3, and NO2 exposure showed a strong association with an increase in visits for upper respiratory diseases. Age was a relevant factor: individuals aged between 19–35 years were the most affected, whereas those aged 66–100 years were least affected. An effect on the weekdays was identified, with peaks in visits on Mondays and decreases on Saturdays. Also, seasonally, results showed an increase in October and declines in August. Regarding hypertensive diseases, only PM2.5 exhibited an association. The models confirmed the link between air pollution levels and respiratory disease exacerbations, highlighting key public health implications.

Conclusion:

This study provides strong evidence that short term exposure to elevated concentrations of atmospheric pollutants, specifically PM2.5, O3, and NO2, are significantly associated with an increase in medical visits for upper respiratory diseases, asthma, and hypertension in Mexico City. By employing Generalized Additive Mixed Models and analyzing health records from IMSS-affiliated hospitals between 2017 and 2019, we observed that the effects of air pollution varied by disease, age group, season, and day of the week.

Introduction

PM2.5 consists of a complex mixture of chemical compounds, including heavy metals, polycyclic aromatic hydrocarbons, organic carbon, trace elements, and secondary ions such as sulfates and nitrates. A study conducted across multiple regions in Asia revealed that a significant proportion of PM2.5 mass originated from organic matter, primarily associated with biomass burning and from secondary inorganic aerosols, predominantly sulfates. The main contributors identified included vehicular emissions, biomass combustion, inorganic ions, and mineral dust. Additionally, elevated levels of nitrates and heavy metals were detected. These pollutants were particularly prevalent in samples collected from urban, agricultural, and industrial areas. Key contributors include traffic emissions, construction activities, open burning, and industrial processes. This evidence reinforces the predominant role of anthropogenic sources in PM2.5 formation. Furthermore, the study demonstrated that meteorological variables, including temperature, humidity, and wind speed, significantly influence both the composition and ambient concentration of PM2.5 ().

Although PM2.5 is a key contributor to the pathophysiological mechanisms underlying various respiratory and cardiovascular diseases, it rarely exists alone. Instead, it frequently interacts with other atmospheric pollutants, leading to complex chemical dynamics. For instance, carbon monoxide (CO) and sulfur dioxide (SO2) can undergo multiple oxidation processes in the atmosphere, resulting in the formation of hydroxyl radicals, ozone (O3), and sulfate ions (SO42−). These reactive species can modify the chemical composition of the air and influence PM2.5 concentrations. Evidence shows that when CO levels exceed 3.5 μg/m3 and SO2 surpasses 6.0 μg/m3, there is a significant increase in PM2.5 concentrations (). Furthermore, a separate study analyzing the interaction of these pollutants found that nitrogen dioxide (NO2) alone was significantly associated with respiratory diseases. However, synergistic effects between NO2 and PM2.5 were linked to a higher incidence of circulatory diseases, while combined exposure to NO2 and SO2 was associated with increased respiratory morbidity ().

International studies have extensively explored the link between respiratory and cardiopulmonary diseases and air pollution. This research consistently shows that exposure to pollutants like particulate matter (PM2.5 and PM10), nitrogen dioxide (NO2), and ozone (O3) contributes to the development and exacerbation of respiratory conditions, including asthma, chronic obstructive pulmonary disease (COPD), and lung cancer. Beyond respiratory diseases, air pollution is strongly linked to cardiopulmonary conditions, such as ischemic heart disease, heart failure, arrhythmias, and stroke. Air pollution affects cardiopulmonary health through oxidative stress, inflammation, endothelial dysfunction, and autonomic imbalance. All these factors contribute to the progression of atherosclerosis, thrombosis, and, ultimately, cardiovascular events ().

Impact on cardiovascular health

Air pollution is actually ranked as the 4th highest risk factor for cardiovascular mortality, attributing more deaths to it than high LDL cholesterol, high body-mass index, physical inactivity, or alcohol use (, ). From a pathophysiological standpoint, to PM2.5 exposure is associated with abnormal activation of the hemostatic system, which can lead to thrombus formation (). PM2.5 also promotes vascular inflammation and endothelial cell injury contributing to coronary and carotid atherosclerosis. Furthermore, it can trigger atherosclerotic plaque rupture by activating metalloproteinases (). A study demonstrated that a 10 μg/m3 increase in PM2.5 was associated with an 11% increase in cardiovascular mortality (). Another cohort study evaluated the long-term effects of PM2.5 exposure and found a positive association between long-term exposure to PM2.5 and stroke incidence ().

Although the physiopathology of asthma is not well known, it is suggested that air pollution promotes the release of epithelial cytokines, driving TH2 responses, which induce oxidative stress and produce proinflammatory cytokines, disrupting of the airway epithelial barrier (). Among environmental factors, fine particulate matter (PM2.5) stands out, given its potential to trigger inflammatory signaling and induce oxidative stress (). Long-term air pollution exposure has been associated with increased hospital admissions, reduced lung function, and a heightened risk of premature death. Large-scale cohort studies across North America and Europe to Asia consistently show that individuals living in areas with higher pollution levels face an elevated risk of developing these diseases. For instance, research in Europe and the United States have found a clear correlation between elevated levels of PM2.5 and increased mortality from cardiovascular and respiratory diseases. Moreover, studies in developing countries, where air pollution levels are often higher due to rapid industrialization and urbanization, have further underscored air pollution's global health burden.

These findings highlight the urgent need for policies to reduce emissions and improve air quality to protect public health (, , , , ).

In Mexico City, the location and the orography also contribute to worsened air pollution. For example, the highest O3 concentrations are more frequent in the south of the urban area. Conversely, the PM2.5 maximum concentrations are generally observed in the city's north. NO2 levels peak in the north and center, correlating with heavier vehicular traffic ().

This study used Generalized Additive Mixed Models (GAMMs) to fit the frequency of respiratory and cardiovascular hospital visits to air pollutant concentrations in Mexico City. GAMMs consider temporal autocorrelation and the complex non-linear relationships between pollutants and upper respiratory and cardiovascular diseases.

Generalized Additive Mixed Models (GAMM) are an advanced extension of Generalized Additive Models (GAM) and Generalized Linear Mixed Models (GLMM) (). Their primary purpose is to analyze complex relationships in data, especially in time-series studies, such as the present study, where assumptions of linearity or independence may not be met (, ).

We employed GAMM to analyze the association between daily atmospheric pollutant levels and the frequency of hospital visits. The choice of GAMM is due to their design to capture temporal autocorrelation and the non-linear relationships that exist among pollutants, meteorological factors, and health outcomes. Furthermore, GAMM possess the ability to handle heterogeneity in terrain and meteorological conditions and can include random effects ().

Smoothing functions or splines are an advantage of GAMM, as they allow for the modeling of these non-linear relationships without assuming a rigid parametric form (such as a straight line or a simple polynomial) (, ).

Materials and methods

Study area

The Metropolitan Zone of the Valley of Mexico is located on one of the central plateaus of the country (Figure 1), with an altitude ranging between 2,223 and 3,510 meters above sea level. It is relevant that the sectors and neighborhoods inhabited at the highest altitude are more than 3,000 meters above the sea level in the south of the city because this creates complex and irregular meteorological behavior. Figure 1 shows an elevation profile from the north to the south of the city. It is possible to observe the high-altitude variation (approximately 1,500 meters in just 60 km), and a plateau can be observed at the south of the urban territory (center of Mexico City). The southern part of Mexico City (53% of the geographical demarcation) is designated as a conservation zone, primarily consisting of mountainous terrain. The climate is temperate, with rainfall in summer (from 600 mm per year in the northern plains to more than 1,500 mm in the southern mountains), with a wet season in summer, from May to October, with a maximum from June to September; the rest of the year is dry, warm from March to May and cold from November to February. During the winter, there can be frost and even occasional snowfall in the upper parts of the surrounding mountains. The months of the dry seasons of the year (warm-dry and cold-dry) are the ones with the highest levels of air pollution. Additionally, the prevailing winds tend to be weak, especially during the dry seasons of the year, with little capacity to carry pollutants, a condition that, together with the direction of the prevailing winds (from the north in winter and northeast in summer), causes that in their displacement many times they are not able to overcome the mountains located in the opposite position, which surround the city on three of its flanks (south, east and west). These situations configure frequent thermal inversions, especially in winter, which aggravate the atmospheric pollution conditions of this urban agglomeration.

Figure 1

In the Air Quality and Health Annual, the Environmental Authorities in Mexico City reported that during December 2023, the concentration of PM2.5 ranged between 60 μg/m3 and 108 μg/m3, showing that the values were above the Mexican National Ambient Air Quality Standards (NAAQS) (). In addition, the concentrations of O3 in the downtown area of the city recorded in that period were maximum average concentrations of 111 ppb and minimums of 14 ppb, whereas NO2 had a maximum of 51 ppb and a minimum of 17 ppb (). Moreover, the Mexican National Standards (NOMs) establish an average concentration of 24-hr of 45 μg/m3 for PM2.5, 70 ppb for O3, and 210 ppm for NO2 ().

Retrospective time series analysis

Hospital visits data

Data on disease-related hospital visits were obtained from the Mexican Institute of Social Security (IMSS) Registration Department. This data encompassed all public hospitals and clinics (primary, secondary, and tertiary care levels) affiliated with IMSS, spanning from 2017 to 2019. Starting in 2017, the IMSS implemented a system capable of recording daily medical care, which is why data from earlier years, reported monthly, could not be included in this analysis. The dataset included the number of hospital admissions according to asthma, upper respiratory diseases, hypertension, and myocardial cardiovascular diseases. Diseases were classified using the International Classification of Diseases (ICD-10) codes in two main categories: acute respiratory diseases (comprising upper respiratory tract infection, bronchitis, bronchiolitis, and/or pneumonia) and asthma. Demographic information, including age and gender, is included. We focused our analysis on children under 5 years old and adults aged 19 years and older. These age groups were selected due to their heightened vulnerability and higher incidence of specific conditions: acute respiratory diseases in children, and acute myocardial infarctions, hypertensive crises, and chronic obstructive pulmonary diseases in adults.

Air pollution data

Hourly air quality data for the period 2017–2019 for PM2.5, O3, and NO2 were obtained from ten fixed air quality monitoring stations in Mexico City and the metropolitan area. Only those stations that provided data for at least 80% of the days and the time period covered were included. Humidity and temperature data were obtained from the Mexico City meteorological network stations. The air quality data was generated by the Automatic Network of Atmospheric Monitoring System of Mexico City, established by Mexico City's Environmental Agency (SEDEMA, according to its Spanish acronym) (). The pollutant averaging times used for the study correspond to current air quality standards in Mexico established by the Mexican government and outlined in the official Mexican regulations (“Norma Oficial Mexicana”) (NOM-025-SSA1-2014, NOM-020-SSA1-2014, NOM-023-SSA1-1994) () which are, at the same time, congruent with health studies and are listed below:

  • O3: Daily maximum of 8 h moving average (ppb).

  • PM2.5: 24 hours average (μg/m3).

  • NO2: Daily maximum of 1-h concentration (ppb).

Data were included in the analysis if they met an 80% completeness threshold. This ensured that data for each pollutant were calculated for at least 17 h per day. To maximize the number of observations for the analysis, missing data were imputed using a moving average method that incorporated both previous and succeeding available data points. The maximum number of consecutively missing values permitted for imputation was set to six to avoid unreliable estimations over extended periods (see Supplementary Figure 1).

Pollutant concentrations were assigned to hospitals based on the spatial representativeness of each air quality monitoring station. Representativeness areas were defined in a study by SEDEMA () which evaluated land use, emission sources, and environmental conditions surrounding each station It follows the criteria recommended by the U.S. Environmental Protection Agency (EPA) in Title 40 of the Code of Federal Regulations (CFR 40). For each hospital, pollutant levels were assigned from the nearest monitoring station whose representativeness buffer included the hospital.

To avoid misclassification from distant monitoring stations, a maximum distance threshold was applied. Only hospitals located within 5 km of an eligible station were assigned pollutant data. This strategy is justified by the findings of the SEDEMA study. Where the mean representativeness radius across all pollutants was 5,006 meters, and 5,346 meters specifically for PM2.5, O3, and NO2. Thus, a 5 km threshold was selected as a limit to ensure that assigned concentrations reflected local exposure conditions while preserving adequate spatial coverage. The hospital's location and the spatial representativeness of monitoring stations are shown in Figure 2.

Figure 2

The spatial scale of the representativeness of an atmospheric monitoring station defines the geographic area whose air quality conditions can be approximated by the measurements from that station. This scale varies depending on the type of pollutant, the characteristics of the emission source, the topography, and the weather patterns.

The representativeness of the stations is urban and extends up to a 5 km radius, representing population exposure, general and background conditions, pollutant transport, and their impact on wellbeing.

A 5 km radius around an air quality monitoring station is commonly used as a guideline to assess the representativeness of the station's data. This limit is based on the usual spatial variation of air pollutants and the practical constraints of monitoring networks. It indicates that within this radius, the air quality is likely similar enough for the station's measurements to reflect the area's conditions.

Statistical analysis

Generalized Linear Mixed Models (GLMM) with Poisson regression were used to fit the time series of selected medical events. The method was adapted from that described by Tadano et al. (). Poisson regression was chosen because the outcome variable, daily counts of medical events, is a non-negative integer count variable that follows a Poisson distribution. Poisson distribution is commonly used in time series studies of environmental epidemiology to model the occurrence of events, such as hospital admissions or emergency visits over time. These time series were related to daily concentrations of air pollutants (derived from the hourly observations), controlling for temporal trends, weather indicators, day of the week, and epidemics. Likewise, the average daily (2017–2019 analysis) relative humidity and temperature and the season of the year were controlled. In addition, for the daily analysis (2017–2019), the day of the week and holidays were controlled.

The relative risks of a hospital admission associated with an increase of 10 units (μg m-3) of PM10 and an increase of 5 units of PM2.5 and NO2 were calculated. In addition, the relative risks for an increase in the interquartile range of contaminants were calculated. For the daily analysis, delay times from 0 to 7 days were considered. The Akaike Information Criterion (AIC) was used to identify the best-fitting model. Analyses were performed separately for adults (over 18 years) and for children (under 5 years).

GAMMs helped to fit the hospital ad into account temporal autocorrelation and the complex non-linear relationships between pollutants and upper respiratory diseases. GAMMs models can handle the complex variations in Mexico City that are caused mainly by the orography shown in Figure 1. GAMMs extend the capabilities of generalized additive models, including random effects, in addition to the parametric additive terms. Hence, the splines can be replaced by random effects as Equation 1 shows ().

Where:

git) Is the link function (log) connecting the mean response μit(expected hospital admissions for locationiat timet) to the predictors

= is a linear parametric term with parameter

= is an additive term with functions α(1 ), …, αm

contains the random effects.

In addition, the use of different distributions can be complemented by modeling residual autocorrelation through Auto Regressive Moving Average (ARMA) functions ().

For model validation and selection, the Akaike Information Criteria (AIC) was used, where its lowest values indicate better model fitting. The AIC is defined in Equation 2:

Where K is the number of estimated parameters in the model and is the maximum value of likelihood function of that model.

Results

The monthly behavior of visits had a similar trend among hospitals, with a peak in 2018 where seasonal behavior can be observed. However, part of the information about ten hospitals was missing. In this case, there was no imputation because of a lack of randomness in those missing points, which consist of continuous periods of time (Figure 3).

Figure 3

The models that exhibited the best fit were the ones that accounted for all three pollutants as factors that influence the frequency of disease-related visits. The results showed a significant association between the number of visits for upper respiratory diseases and the three pollutants under consideration. Specifically, O3 and NO2 had a positive effect on the outcome measure, while PM2.5 exhibited a non-linear impact requiring a smoother. Similar behavior was found in relative humidity, with different percentages having different effects on the number of visits. From a range between 0 to 12%, it had a negative partial effect. Above 12%, the effect is positive, and then it inverts the effect from 17%. These findings are supported by the results presented in Table 1 and Figure 4.

Table 1

VariableEstimate(Pr > t) < 0.05Visits
O34.80E-04*+
NO22.16E-04*+
Age 66–1003.30E+00*+
Age 51–653.62E+00*+
Age 36–503.89E+00*+
Age 0–54.53E+00*+
Age 19–354.66E+00*+
Monday2.13E-01*+
Tuesday9.14E-02*+
Wednesday6.49E-02*+
Thursday1.07E-02NS*+
Friday5.30E-03NS+
Saturday−2.52E-02*-
August−2.54E-01*-
NS
July−2.16E-01*-
June−1.48E-01*-
May−9.36E-02-
November1.23E-01*+
September1.23E-01*+
March3.57E-01*+
January4.26E-01*+
December4.33E-01*+
February4.70E-01*+
October4.82E-01*+
Lag day14.00E-02*+
Lag day 29.97E-03*+

Result of the best fitting model for upper respiratory diseases.

The sign in the visits column means increment (+) or reduction (–) in the number of visits. NS, Statistically non-significant (p > 0.05), R2: 0.69. Smoother: s(rh), s(pm25). *p-value < 0.05.

Figure 4

Regarding age, individuals aged between 19–35 and 66–100 years were most and least affected, respectively. The analysis also showed that Monday was associated with an increase in the number of visits, whereas Saturday was associated with a reduction.

Regarding seasonal effects, October had the most substantial increase in the number of visits, whereas August had a negative impact. Regarding day lags, a significant positive effect was found on the first 2 days of exposition.

Regarding the ACF and PACF plots, the residuals demonstrated minimal correlation, with values falling inside the autocorrelation thresholds. This indicates that the correlations are likely due to random chances and are not meaningful. This was accomplished by incorporating p = 1 and q = 1 terms into the correlation matrix within the GAMM model.

Regarding fitting, the Poisson distribution was more suitable for the data. However, there is still information that the model couldn't fit.

Relative risk for O3 and NO2 is defined as:

Which means that if there is an interval of 50 ppb and then there is an increment to 60 ppb for O3, the visits will increase in 10. In the same way, an increase in NO2 concentration from 60–70 ppb for NO2, the visits due to upper respiratory diseases will be increased in 10.

On the other hand, PM2.5 could not be interpreted only with one relative risk. This can be observed in the smoothers for this pollutant (Figure 5). For instance, in 37 μg/m3 and 63 μg/m3, there are the maximum and minimum relative risks respectively:

Figure 5

The smoother plots provide a visual representation of the non-linear relationship between PM2.5 concentrations and the relative risk of upper respiratory diseases. These plots are generated using GAMMs, which allow for the inclusion of random effects and parametric additive terms. This modeling approach is particularly useful in capturing the complex interactions between pollutants and health effects.

The non-linear pattern observed in the smoother plots highlights the importance of considering varying pollutant levels rather than relying on a single RR value. For example, at lower concentrations of PM2.5, the RR may increase, indicating a higher risk of health issues. However, as the concentration increases, the RR may decrease, suggesting a different impact on health effects.

Regarding hypertensive diseases, the effect of pollutants was only found for PM2.5 as smoother. However, the partial effect was nearly zero (Figure 5). A positive effect was observed in older people about age. Table 2 shows that, in older people, visits to hospitals due to hypertensive diseases were more common. People aged 19 to 35 years were the least affected by hypertensive disease. Regarding the day of the week, all days have a positive and similar effect on the visits, except Saturday, which had a low and adverse effect on that variable. Regarding the year's seasonality, only January, May, and June had statistically significant effects, negatively influencing the visits. The impact of lag was found for day 1. The model could fit the data with an ARMA correlation matrix of p = 1 and q = 1 for AR and MA components.

Table 2

VariableEstimate(Pr > t) < 0.05Visits
Age 19–351.81*+
Age 36–502.87*+
Age 51–653.40*+
Age 66–1003.83*+
Monday0.66*+
Tuesday0.66*+
Wednesday0.68*+
Thursday0.64*+
Friday0.65*+
Saturday0.03*-
August−0.02NS-
July−0.05NS-
June−0.06*-
May−0.07*-
November−0.02NS-
September−0.01NS-
March−0.02NS-
January−0.08NS-
December−0.18*-
February−0.02NS-
January−0.08*-
October0.09NS+
Lag day 10.38*-

As a result of the best model for hypertensive diseases, the sign in the visits column means an increment (+) or reduction (-) in the number of visits.

NS, Statistically non-significant (p > 0.05), R2: 0.87. Smoother: s(pm25). *p-value < 0.05.

The relative risk for PM2.5 remains constant throughout the entire range of data examined in the study. This is expressed as:

The group most affected by asthma was those aged 36–50 years, in contrast to the least affected group of hypertensive diseases being those aged 66–100 years. All days of the week positively affected the number of visits, while PM2.5 and O3 had significant non-linear effects, exhibiting a cyclic pattern. This behavior refers to the recurring patterns observed in the effects of PM2.5. These patterns can be visualized using smoother plots shown in Figures 5, 6, which show how the RR varies with different concentrations of the pollutant. Although Saturday showed a negative effect, it was not substantial. Seasonality effects were only observed in January and February, negatively impacting the number of visits. Regarding lags, all days had adverse effects, except for days 4 and 10 which showed positive effects. Autocorrelation in the residuals was negligible based on ACF and PACF plots, indicating that no significant temporal information remained. This was verified using ARMA terms (q = 1, p = 1) incorporated into the correlation matrix of the GAMM.

Figure 6

For O3 and PM2.5 relative risks. They can be only interpreted in terms of the partial effects in their smoother (Table 3). Figure 6 shows the following minimum and maximum relative risks for PM2.5 in 87 and 37 respectively:

In the case of O3 the maximum and minimum relative risks are in:

Table 3

VariableEstimate(Pr > t) < 0.05Visits
Age 0–51.94*+
Age 19–352.29*+
Age 36–502.37*+
Age 51–652.16*+
Age 66–1001.66*+
Monday0.52*+
Tuesday0.43*+
Wednesday0.46*+
Thursday0.41*+
Friday0.30*+
Saturday0.03-
August0.11-
July0.20-
June0.07-
May0.06+
November0.05+
September0.27+
March0.06-
December0.10-
February0.26*-
January0.25*-
October0.24+
Lag day 11.52*-
Lag day 20.26-
Lag day 31.64*-
Lag day 41.05*+
Lag day 50.67*-
Lag day 61.23*-
Lag day 70.84*-
Lag day 81.09*-
Lag day 90.31*-
Lag day 100.76*+

As a result of the best model for asthma, the sign in the visits column means an increment (+) or reduction (–) in the number of visits.

NS, Statistically non-significant (p > 0.05), R2: 0.55. Smoother: s(O3). *p-value < 0.05.

Discussion

This study demonstrated a relationship between different environmental pollution components and increased medical consultations for respiratory diseases, asthma, and hypertension. Specifically O3 and NO2 were positively associated with the number of respiratory disease consultations, with the most affected age groups being 19–35 and over 66. For hypertension consultations, PM2.5 showed a more significant effect than O3 and NO2, and a stronger correlation was observed with increasing age. Asthma consultations followed a cyclical pattern with both PM2.5 and O3 across all age groups, though the correlation was weaker in the over 66 group.

The present study was conducted using data from the Mexican Social Security Institute population, because of the available and reliable data regarding medical care demand. During the study period, the eligible population for medical care through Social Security was almost half of Mexico City's total population (). This suggests that our findings may be extrapolated to at least 50% of Mexico City's population and therefore to other cities with similar demographic and healthcare characteristics. By including consultations from all levels of medical care, namely in family medicine units, general hospitals, and high specialty hospitals, our results may be applicable to all levels of care.

The data recording is highly reliable because it is conducted under a social security institution. However, it requires individuals seeking medical attention to be registered and enrolled in the IMSS (Mexican Social Security Institute). This study collected information exclusively from the IMSS database. Within IMSS, qualified personnel register medical consultations in the database, strictly adhering to the criteria set by the International Classification of Diseases, further ensuring data quality.

The measurement of environmental pollution in Mexico City is standardized. As a result, the operation of monitoring equipment data collection, and information recording undergo constant review and validation. This standardization provides robust quality assurance, thereby ensuring reliable data on exposure to various environmental pollutants throughout Mexico City.

Respiratory diseases

The strongest association was observed with O3 and NO2, indicating that these pollutants significantly impact. While PM2.5 showed a weaker association, it was still statistically significant.

Respiratory diseases, particularly bronchiolitis, are a leading cause of hospital admission in children under 2 years old. Although poor air quality affects everyone, children are disproportionately impacted due to their physiological and immune system immaturity. Their smaller airway caliber, higher metabolism, and greater respiratory rate compared to adults, lead to larger inhaled air volumes. This amplifies the effect of pollutants, overwhelming their limited capacity to neutralize and eliminate these environmental contaminants (, ). Consistent with this, studies have linked NO2, primarily from vehicle emissions in urban areas, to an increased incidence of lower respiratory tract infections and a higher risk of hospital admission for bronchiolitis (, ). Furthermore, Yamazaki et al., observed a positive association between PM2.5 and O3 concentrations and a reduced lung function in children ().

For upper respiratory illnesses, our study found a positive correlation between these pathologies and elevated PM2.5, NO2, and O3 concentrations. These findings align with previous research; for instance, one study also reported a positive association between PM2.5 and the monthly case counts in outpatient consultations (). Another article similarly linked daily PM2.5 concentrations to hospital admissions for upper respiratory tract infections in adults (). It is believed that PM2.5 leads to inflammation and cell damage by producing free radicals and reactive oxygen species. Such inflammation and respiratory system damage could increase susceptibility to upper respiratory diseases ().

Regarding O3 and NO2, a 2022 time-series study identified a positive association between short-term exposure to these pollutants and emergency department visits for respiratory diseases and upper respiratory infections (). As O3 and NO2 are the two main oxidative gaseous pollutants in the atmosphere; the mechanisms and pathways through which these pollutants could affect the respiratory system include systemic oxidative injuries and airway inflammation, leading to bronchoconstriction and systemic inflammatory responses ().

Hypertension

In the present study, PM2.5 particles showed the strongest association with hypertension.

Our findings align with research demonstrating that PM2.5 is linked to an increment in clinic visits for hypertensive diseases ().

Further support comes from investigations in Guangzhou, China, which identified a significant risk correlation between ambient PM2.5 exposure and outpatient clinic visits and hospitalizations for hypertension (, ). In that study, PM2.5 was the only pollutant to show a substantial effect. A similar pattern emerged from a research involving older adults, revealing a significant association with higher systolic blood pressure and both non-significant and significant results for diastolic blood pressure ().

Several pathophysiological mechanisms can explain the association between high PM2.5 concentrations and the prevalence of hypertension. Key among these are systemic inflammatory and oxidative responses, which progressively impair vascular function and lead to arterial remodeling (). Other proposed mechanisms include: an altered autonomic nervous system that mediate arterial vasoconstriction, systemic lipid peroxidation leading to the production of pro- inflammatory cytokines or an overstimulated sympathetic nervous system. This can enhance renin- angiotensin activity, improving myocardial contractility, and finally increasing blood pressure (, ).

Asthma

Our study revealed stronger associations between asthma and both PM2.5 and O3 compared to previous research using the same Mexico City data (). A notable distinction lies in our examination of various age groups where we found associations across all demographics. The strength of this association was weaker in individuals over 66 years old. This nuanced finding was not reported in our earlier article ().

These findings are consistent with prior studies, which also linked increased ozone and PM2.5 levels to hospital visits for asthma exacerbations (, ).

Overall, our results align with the broader literature reporting an association between O3 and PM2.5 and acute asthma exacerbations. It is also widely recognized that air pollutant-induced asthma exacerbations are more prevalent in children than in adults ().

From the physiopathological perspective, air pollutants primarily trigger an immune response mediated by oxidative stress. PM2.5 contains biological and organic compounds with oxidative stress effects, which cause damage to the airway's mucosal barrier and cellular epithelium. This triggers a type 2 inflammatory response leading to the release of proinflammatory cytokines and an increased flow of inflammatory cells. This activates various signaling pathways such as NF-κB or PI3K/Akt, driving inflammatory T2 cell activity and, consequently, airway hyperreactivity (, ).

Regarding O3, its mechanism generates an accumulation of reactive oxygen species (ROS) through lipid peroxidation of pulmonary surfactant factor and lipid membranes. This ROS accumulation activates the release of IL-1β, IL-6, IL-23, IL- 33, TNF-α, and TSLP, causing a proinflammatory cascade in the respiratory mucosa. Subsequent activation of the RORγt gene leads to the transcription of IL-17A and IL-22. IL-1, IL-17A, and IL-22 then induce neutrophil activation, prompting these active neutrophils to release more ROS. This cascade results in a high flow of inflammatory cells, representing a critical point in severe asthma exacerbations ().

Public health implications

In LATAM, 120 million people were reported to be exposed to pollution levels that exceed WHO guidelines, and consequently, 7 million deaths. Mexico, as a county of the global south, experiences high levels of contamination due to strong industrialization, caused by fossil-fuel generated energy, vehicular pollution, thermal power generation, industrial energy production, open burning of municipal or agricultural waste, and especially due to poor management by the government. In LATAM, Mexico is the country with the most deaths attributable to air pollution from burning fossil fuels (). The global burden of disease reported an estimated 4.58 million deaths attributable to exposure to PM2.5 in 2017, in a study conducted in 195 countries and territories worldwide ().

In Mexico, this problem has been addressed with the instauration of programs and policies. Exposure to ambient PM2.5 contributes to a high pollution linked death rate. The improvement of the Vehicle Verification Program, the control of ostensibly polluting vehicles in Mexico City, has reportedly improved the air quality since 2018, resulting in 9% decrease in PM2.5 particles, 4% in PM10 and nitrogen dioxide, 3% in sulfur dioxide, and 2% in carbon monoxide (). The Government of Mexico started a program “PROAIRE” in 2017. It aims to reduce by 2030, 75% of SO2; 58% of PM2.5; 46% of VOCs; and 32% of NO2 from total emissions () through strategies like: air quality monitoring program, atmospheric environmental contingency programs, sentinel health units' program for evaluating health effects from air pollutants and an inventory of atmospheric pollutant emissions ().

In other countries, policies have been implemented to control air pollution, thereby reducing the public health implications. In Global North countries like Europe and North America, the push for the manufacturing and use of electric vehicles and low-emission fuels has helped decrease emissions from burning fuel. Continuous improvement in public transport and the establishment of bike lanes and sidewalks have also improved air quality in urban areas (). China stands as an example that improved air pollution regulation policies lead to significant improvements in public health. In 2013, the State Council implemented the Air Pollution Prevention and Control Action Plan (APPCAP), which was, at the time, the strictest air pollution regulation policy to date. This initiative resulted in 47,240 fewer deaths and 710,020 fewer years of life lost in 74 major cities in 2017 compared to 2013. Furthermore, APPCAP led to a 5.6% decrease in mortality attributable to PM2.5 in Beijing in 2018 compared to 2014 (). These results highlight the importance of understanding the health implications of specific air pollution components and the necessity of implementing effective actions to diminish them. To protect public health, action-oriented research should inform and drive targeted interventions aimed at reducing population exposure to harmful air pollutants. Our findings offer critical insights into exposure-risk thresholds, which can be translated into evidence-based public health policies. Specifically, our exposure-response analyses demonstrate that relatively small increases in PM 2.5 and O3 concentrations are associated with measurable adverse health impacts. These findings underscore the need for comprehensive air quality interventions, including emissions controls, the promotion of clean transportation, and the expansion of urban green spaces. In addition, efforts to improve access to healthcare services can help mitigate the health burden of air pollution and address disparities in exposure among vulnerable urban populations ().

Currently, mobile applications make it easy to ascertain the levels of environmental pollutants in the air. This capability would allow decision-makers to guide the population to take greater precautions when pollutant levels rise to figures that can be associated with an exacerbation of cardiovascular and respiratory diseases.

The study has some limitations. It relied solely on data from individuals enrolled in the IMSS, which may exclude informal workers or those covered by other healthcare systems, thus limiting generalizability to the entire population.

These findings underscore the importance of designing educational strategies targeted at populations with asthma and hypertension to prevent acute episodes that could escalate into severe events, especially during periods of heightened pollutant levels. The most vulnerable populations to respiratory illnesses should also remain vigilant when pollutant levels exceed the thresholds established by Mexican regulations.

It is also critical to implement public health policies and educational interventions aimed at reducing exposure to harmful air pollutants. In several of the most polluted countries in the world, particularly in China, and now increasingly in Mexico, there is growing awareness of the negative effects of pollution, not only from an environmental standpoint but also as a major public health concern. This heightened awareness has led to greater action, with policies aimed at reducing atmospheric pollutant levels. However, as demonstrated in this study, these efforts remain insufficient, given the persistently high concentrations of environmental pollutants, often exceeding the levels recommended by the WHO, and their impact on the public healthcare system. There is still significant work to be done in Mexico to meet the targets set for 2030.

Strengthening environmental regulations to maintain pollutant levels below national and international thresholds could significantly reduce the burden of preventable cardiorespiratory diseases. Future research should continue exploring these associations and assess the long-term effects of repeated exposure to atmospheric pollutants. It could also assess the impact of recent air quality regulations. Expanding the analysis to other Mexican cities or Latin American urban areas would also strengthen the regional understanding of pollution related health outcomes.

Conclusion

There is a well-established global association between air pollution and both respiratory and cardiovascular diseases. This study provides strong evidence that short term exposure to elevated concentrations of atmospheric pollutants, specifically PM2.5, O3, and NO2, are significantly associated with an increase in medical visits for upper respiratory diseases, asthma, and hypertension in Mexico City. By employing Generalized Additive Mixed Models and analyzing health records from IMSS-affiliated hospitals between 2017 and 2019, we observed that the effects of air pollution varied by disease, age group, season, and day of the week. Respiratory diseases showed the strongest associations with O3 and NO2, particularly affecting individuals aged 19–35 and children under 5 years of age. Hypertension was predominantly linked to PM2.5, with older adults being the most affected. Asthma was related to PM2.5 and O3, affecting all age groups. A key strength of this study is its ability to provide insight into the public health burden of air pollution within an urban context with unique topographic and meteorological characteristics like Mexico City.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the Ethics and National Committee of Scientific Research approved this study with number: R-2020-785-045. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

EV: Validation, Writing – review & editing, Conceptualization, Investigation, Supervision, Methodology, Resources, Writing – original draft, Data curation. AR-O: Validation, Writing – review & editing, Software, Methodology, Formal analysis, Data curation, Writing – original draft. YM: Supervision, Validation, Writing – review & editing, Methodology, Data curation, Investigation, Writing – original draft. RR-C: Methodology, Supervision, Validation, Investigation, Data curation, Writing – review & editing, Writing – original draft. LB: Writing – original draft, Investigation, Data curation, Validation, Supervision, Methodology, Writing – review & editing. RF: Writing – review & editing, Supervision, Methodology, Writing – original draft, Validation, Investigation, Data curation. AN: Investigation, Validation, Methodology, Data curation, Supervision, Writing – original draft, Writing – review & editing. JN-E: Methodology, Data curation, Supervision, Investigation, Validation, Writing – review & editing, Writing – original draft. JE: Methodology, Supervision, Data curation, Investigation, Validation, Writing – review & editing, Writing – original draft. MJ-P: Data curation, Validation, Investigation, Methodology, Supervision, Writing – review & editing, Writing – original draft. NT-R: Methodology, Supervision, Writing – review & editing, Investigation, Writing – original draft, Data curation, Validation. RR-P: Methodology, Writing – review & editing, Data curation, Investigation, Supervision, Validation, Writing – original draft. DD-R: Validation, Writing – review & editing, Formal analysis, Software, Writing – original draft, Data curation, Methodology. RM: Validation, Writing – review & editing, Data curation, Methodology, Writing – original draft, Supervision, Investigation. JG-B: Validation, Methodology, Data curation, Investigation, Writing – review & editing, Supervision, Writing – original draft. CR-M: Writing – review & editing, Supervision, Methodology, Writing – original draft, Data curation, Investigation, Validation. ÁF-C: Supervision, Methodology, Writing – review & editing, Validation, Investigation, Data curation, Writing – original draft. LH: Validation, Data curation, Methodology, Writing – review & editing, Investigation, Supervision, Writing – original draft. AP: Methodology, Supervision, Validation, Data curation, Writing – review & editing, Investigation, Writing – original draft. JM-A: Validation, Methodology, Investigation, Conceptualization, Supervision, Writing – review & editing, Writing – original draft, Data curation.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Secretaria de Educacion, Ciencia, Tecnologia e Innovacion funded project (Award No. SECTEI/203/2019, Mexico) and the British Council (Award No. 527635178). National Institutes of Health (Administrative Supplement 3R01CA262263 to Juan Manuel Mejía-Arangure, 2U24ES028524 sub-award number 00011320 to Juan Manuel Mejía-Arangure), the Consejo Nacional de Humanidades, Ciencias y Tecnologías (CONAHCYT, CF-2023-G-1399 to Juan Manuel Mejía-Arangure) and Dirección General de Políticas de Investigación en Salud (DGPIS), Financiamiento de Proyectos de Investigación para la Salud (FPIS 2024) FPIS2024-INMEGEN-7593. The funders had a role in the study design, in the collection and analysis of data, in the writing of the report, and in the decision to submit the article for publication. All researchers have independence from funders.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1593285/full#supplementary-material

References

Summary

Keywords

hypertension, respiratory diseases, particulate matter 2.5 μm (PM) 2.5, Mexico City (CDMX), air pollution, asthma, children, aging

Citation

Vega E, Ruiz-Olivares A, Miquelajauregui Y, Rentería-Campos R, Bramwell L, Feytmans RBM, Namdeo A, Núñez-Enríquez JC, Entwistle JA, Jaimes-Palomera M, Torres-Rojas NM, Rascón-Pacheco RA, Duarte-Rodríguez DA, McNally RJQ, García-Burgos J, Resendiz-Martinez CG, Fragoso-Chino Á, Hayes L, Portas A and Mejia-Arangure JM (2025) Association between daily variations in the levels of atmospheric O3, PM2.5, and NO2 and the frequency of hospital visits due to respiratory diseases and hypertension in Mexico City using Generalized Additive Mixed Models. Front. Public Health 13:1593285. doi: 10.3389/fpubh.2025.1593285

Received

13 March 2025

Accepted

24 July 2025

Published

02 September 2025

Volume

13 - 2025

Edited by

Haider Khwaja, University at Albany, United States

Reviewed by

Worradorn Phairuang, Chiang Mai University, Thailand

Katalin Hubai, University of Pannonia, Hungary

Updates

Copyright

*Correspondence: Elizabeth Vega Juan Manuel Mejia-Arangure

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics