Your new experience awaits. Try the new design now and help us make it even better

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


Elizabeth Vega
Elizabeth Vega1*Alejandro Ruiz-OlivaresAlejandro Ruiz-Olivares2Yosune MiquelajaureguiYosune Miquelajauregui3Regina Rentería-CamposRegina Rentería-Campos4Lindsay BramwellLindsay Bramwell5Romain Bernard M. FeytmansRomain Bernard M. Feytmans4Anil NamdeoAnil Namdeo5Juan Carlos Núez-EnríquezJuan Carlos Núñez-Enríquez6Jane A. EntwistleJane A. Entwistle5Mnica Jaimes-PalomeraMónica Jaimes-Palomera7Nancy Minerva Torres-RojasNancy Minerva Torres-Rojas8R. Alberto Rascn-PachecoR. Alberto Rascón-Pacheco9David A. Duarte-RodríguezDavid A. Duarte-Rodríguez6Richard J. Q. McNallyRichard J. Q. McNally10Jimena García-BurgosJimena García-Burgos11Cinthia Gabriela Resendiz-MartinezCinthia Gabriela Resendiz-Martinez12ngel Fragoso-ChinoÁngel Fragoso-Chino12Louise HayesLouise Hayes10Antonio PortasAntonio Portas13Juan Manuel Mejia-Arangure
Juan Manuel Mejia-Arangure14*
  • 1Instituto de Ciencias de la Atmósfera y Cambio Climático, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico
  • 2Department of Soil Sciences, Colegio de Postgraduados, Montecillo, Mexico
  • 3Laboratorio Nacional de Ciencias de la Sostenibilidad, Instituto de Ecología, Universidad Nacional Autónoma de México (UNAM), Mexico City, Mexico
  • 4Facultad de Medicina, Universidad Anáhuac México Sur, Mexico City, Mexico
  • 5Department of Geography and Environmental Sciences, Northumbria University, Newcastle upon Tyne, United Kingdom
  • 6Unidad 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
  • 7Dirección de Monitoreo de Calidad del Aire, Secretaria del Medio Ambiente, Gobierno de la Ciudad de Mexico, Mexico City, Mexico
  • 8Coordinación de Vigilancia Epidemiológica, Instituto Mexicano del Seguro Social, Mexico City, Mexico
  • 9Unidad de Educación, Investigación y Políticas de Salud, Instituto Mexicano del Seguro Social, Ciudad de Mexico, Mexico
  • 10Population Health Sciences Institute, Newcastle University, Newcastle upon Tyne, United Kingdom
  • 11Laboratorio 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
  • 12Análisis y Modelación de Calidad del Aire, Secretaría de Medio ambiente de la Ciudad de México, Mexico City, Mexico
  • 13Department of Mathematics, Physics and Electrical Engineering, Northumbria University, Newcastle upon Tyne, United Kingdom
  • 14Genómica del Cáncer, Instituto Nacional de Medicina Genómica, Mexico City, Mexico

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 (1).

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 (2). 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 (3).

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 (46).

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 (7, 8). From a pathophysiological standpoint, to PM2.5 exposure is associated with abnormal activation of the hemostatic system, which can lead to thrombus formation (7). 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 (9). A study demonstrated that a 10 μg/m3 increase in PM2.5 was associated with an 11% increase in cardiovascular mortality (3). 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 (10).

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 (11). Among environmental factors, fine particulate matter (PM2.5) stands out, given its potential to trigger inflammatory signaling and induce oxidative stress (12). 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 (7, 8, 10, 12, 13).

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 (14).

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) (15). 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 (15, 16).

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 (15).

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) (15, 16).

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
Map of Mexico with an inset showing a detailed area outlined in red. A smaller graph displays elevation data over a distance, indicating a rise from approximately 2,500 to 3,000 meters. The inset map highlights the region with red crosses, indicating specific points of interest. Scale bars are present in both the main and inset maps.

Figure 1. Mexico City location and landscape. The black polygon is the Country, and the red polygon is Mexico City. The plot is the elevation profile of the red dashed line drawn in the city (north to south).

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) (7). 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 (7). 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 (6).

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) (14). 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) (1719) 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 (14) 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
Figure 2 displaying three maps (PM2.5, O3, and NO2) of Mexico City, where the circles indicate the areas of representativeness for each monitoring station. The black dots represent the hospitals from which the data is sourced.

Figure 2. Spatial representativeness of the PM2.5, O3, and NO2 monitoring stations (circles) and the hospitals (black dots) considered in the study.

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. (20). 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 (21).

g(μit)=xitTβ + j=1mαj( uitj)+ zitTbi= ηitpar+ ηitadd+ ηitrand    (1)

Where:

git) Is the link function (log) connecting the mean response μit (expected hospital admissions for location i at time t) to the predictors

xitTβ = ηitpar is a linear parametric term with parameter βT= β0+β1,, βp

j=1mαj( uitj) = ηitadd is an additive term with functions α(1 ), …, αm

zitTbi=ηitrand 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 (21).

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:

AIC=2k-2 ln ln(L^)    (2)

Where K is the number of estimated parameters in the model and L^ 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
Grid of line charts showing daily visit trends from 2017 to 2019. Each chart shows variable patterns with shaded areas indicating uncertainty. Axes labeled with years and “Visits per day.

Figure 3. ” Temporal behavior by month of the 49 hospitals during the period of analysis. The black line is the average daily visits, the red ribbon is the standard deviation, and the blue ribbon is the minimum and maximum.

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
www.frontiersin.org

Table 1. Result of the best fitting model for upper respiratory diseases.

Figure 4
(A) Composite image showing six plots: a QQ plot of residuals indicating a deviation from normality, a scatter plot of residuals versus linear predictor, a histogram of residuals, a scatter plot of observed versus fitted values, and plots of partial effects for variables “pm25” and “rh” with confidence intervals shaded. (B) Two side-by-side plots showing ACF and PACF analyses. The ACF plot on the left displays spikes at various lags, indicating autocorrelation. The PACF plot on the right shows partial autocorrelations with decreasing values. Both plots have dashed horizontal lines marking significance levels.

Figure 4. (A) Model fitting evaluation of upper respiratory diseases (left) and PM2.5 and RH smoothers. (B) Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) for upper respiratory diseases.

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:

RRO3=exp (4.8E-04) =1.0 visits per O3 unitRRNO2=exp (2.16E-04) = 1.0 visits per NO2 unit

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
(A) Four-panel plot displaying diagnostic graphs and a partial effect plot. Top left: QQ plot of residuals showing deviations from the line of normality. Top middle: Scatter plot of residuals versus linear predictor, indicating randomness. Bottom left: Histogram of residuals with a normal-like distribution. Bottom middle: Scatter plot of observed versus fitted values showing clustering along the line of equality. Right: Partial effect plot of PM2.5 with a smooth curve and gray confidence interval bands, showing slight positive and negative trends across the PM2.5 range. (B) Two bar graphs showing Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) with lag values on the x-axis and correlation values on the y-axis. The ACF graph shows decreasing bars with notable spikes around lags one to ten. The PACF graph displays a similar pattern with prominent spikes at early lags, then lower values as lag increases. Dashed horizontal lines indicate significance thresholds.

Figure 5. (A) Model fitting evaluation of hypertensive diseases. (B) Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) for hypertensive diseases.

RRPM2.5(37 μg/m3 )=exp (0.02) =1.02 visits per PM2.5 unitRRPM2.5(63 μg/m3)=exp(-0.08)=0.92 visits per PM2.5 unit

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
www.frontiersin.org

Table 2. 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.

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

RRPM2.5=exp (0) =1.0 visits per PM2.5 unit

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
(A) Composite image showcasing six plots related to residuals in statistical analysis. Top left: QQ plot of deviance residuals versus theoretical quantiles, showing a near-linear relation. Top center: Scatterplot of deviance residuals against a linear predictor line with Poisson family. Top right: Partial effect of O3 with confidence interval shaded. Bottom left: Histogram of deviance residuals showing normal distribution. Bottom center: Scatterplot of observed versus fitted values. Bottom right: Partial effect of PM25 with confidence interval shaded. Plots analyze residual distributions and predictors' effects. (B) Two bar graphs showing Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) with lag values on the x-axis and correlation values on the y-axis. The ACF graph shows decreasing bars with notable spikes around lags one to ten. The PACF graph displays a similar pattern with prominent spikes at early lags, then lower values as lag increases. Dashed horizontal lines indicate significance thresholds.

Figure 6. (A) Model fitting evaluation of asthma. (B) Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) for asthma.

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:

RRPM2.5(37 ug/m3 )=exp(0.025)=1.02 visits per PM2.5 unitRRPM2.5(87 ug/m3)=exp(-0.025)=0.97 visits per PM2.5 unit

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

RRO3(27 ug/m3)=exp(0.007)=1.0 visits per O3 unitRRO3(87 ug/m3)=exp(-0.045)=0.95 visits per O3 unit
Table 3
www.frontiersin.org

Table 3. 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.

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 (22). 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 (23, 24). 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 (23, 24). Furthermore, Yamazaki et al., observed a positive association between PM2.5 and O3 concentrations and a reduced lung function in children (25).

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 (26). Another article similarly linked daily PM2.5 concentrations to hospital admissions for upper respiratory tract infections in adults (27). 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 (26).

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 (28). 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 (28).

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 (29).

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 (25, 26). 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 (30).

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 (4). 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 (29, 31).

Asthma

Our study revealed stronger associations between asthma and both PM2.5 and O3 compared to previous research using the same Mexico City data (32). 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 (32).

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

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 (35).

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 (36, 37).

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 (38).

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 (3). 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 (1).

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 (39). 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 (39) 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 (40).

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 (3). 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 (41). 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 (42).

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.

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

1. Bu X, Xie Z, Liu J, Wei L, Wang X, Chen M, et al. Global PM2.5-attributable health burden from 1990 to 2017: Estimates from the Global Burden of disease study 2017. Environ Res. (2021) 197:111123. doi: 10.1016/j.envres.2021.111123

PubMed Abstract | Crossref Full Text | Google Scholar

2. Wei Q, Chen Y, Zhang H, Jia Z, Yang J, Niu B. Simulation and prediction of PM25 concentrations and analysis of driving factors using interpretable tree-based models in Shanghai, China. Environ Res. (2025) 270:121003. doi: 10.1016/j.envres.2025.121003

PubMed Abstract | Crossref Full Text | Google Scholar

3. Mangones SC, Cuéllar-Álvarez Y, Rojas-Roa NY, Osses M. Addressing urban transport-related air pollution in Latin America: Insights and policy directions. Latin Am Transp Stud. (2025) 3:100033. doi: 10.1016/j.latran.2025.100033

Crossref Full Text | Google Scholar

4. Liu C, Chen R, Sera F, Vicedo-Cabrera AM, Guo Y, Tong S, et al. Interactive effects of ambient fine particulate matter and ozone on daily mortality in 372 cities: two stage time series analysis. BMJ. (2023) 383:e075203. doi: 10.1136/bmj-2023-075203

PubMed Abstract | Crossref Full Text | Google Scholar

5. Chen ZY, Petetin H, Méndez Turrubiates RF, Achebak H, Pérez García-Pando C, Ballester J. Population exposure to multiple air pollutants and its compound episodes in Europe. Nat Commun. (2024) 15:2094. doi: 10.1038/s41467-024-46103-3

PubMed Abstract | Crossref Full Text | Google Scholar

6. CDMX. Gaceta oficial de la Ciudad de México. (2018). Available online at: http://www.aire.cdmx.gob.mx/descargas/monitoreo/normatividad/NADF-009-AIRE-2017.pdf (Accessed March 20, 2024).

PubMed Abstract | Google Scholar

7. Gobierno. Dirección de Monitoreo Atmosférico. (2023). Available online at: http://www.aire.cdmx.gob.mx/default.php (Accessed March 20, 2024).

PubMed Abstract | Google Scholar

8. Thangavel P, Park D, Lee YC. Recent insights into particulate matter (PM25)-mediated toxicity in humans: an overview. Int J Environ Res Public Health. (2022) 19:7511. doi: 10.3390/ijerph19127511

PubMed Abstract | Crossref Full Text | Google Scholar

9. Henning RJ. Particulate matter air pollution is a significant risk factor for cardiovascular disease. Curr Probl Cardiol. (2024) 49:102094. doi: 10.1016/j.cpcardiol.2023.102094

PubMed Abstract | Crossref Full Text | Google Scholar

10. Caldeira D, Franco F, Bravo Baptista S, Cabral S, Cachulo MDC, Dores H, et al. Air pollution and cardiovascular diseases: a position paper. Rev Port Cardiol. (2022) 41:709–17. doi: 10.1016/j.repc.2022.05.006

PubMed Abstract | Crossref Full Text | Google Scholar

11. Zhou X, Sampath V, Nadeau KC. Effect of air pollution on asthma. Ann Allergy Asthma Immunol. (2024) 132:426–32. doi: 10.1016/j.anai.2024.01.017

PubMed Abstract | Crossref Full Text | Google Scholar

12. Wang X, Guo Y, Cai M, Qian ZM, Zhang S, Zhang Z, et al. Constituents of fine particulate matter and asthma in 6 low- and middle-income countries. J Allergy Clin Immunol. (2022) 150:214–222.e5. doi: 10.1016/j.jaci.2021.12.779

PubMed Abstract | Crossref Full Text | Google Scholar

13. Fajersztajn L, Saldiva P, Pereira LAA, Leite VF, Buehler AM. Short-term effects of fine particulate matter pollution on daily health events in Latin America: a systematic review and meta-analysis. Int J Public Health. (2017) 62:729–38. doi: 10.1007/s00038-017-0960-y

PubMed Abstract | Crossref Full Text | Google Scholar

14. SEDEMA (Secretaría del Medio Ambiente de la Ciudad de México). Caracterización, evaluación y análisis del entorno físico y de la representatividad en las estaciones del Sistema de Monitoreo Atmosférico de la Ciudad de México. (2021). Available online at: http://www.aire.cdmx.gob.mx/descargas/publicaciones/simat-entornos.pdf (Accessed March 20, 2024).

Google Scholar

15. Fang X, Li R, Kan H, Bottai M, Fang F, Cao Y. Bayesian model averaging method for evaluating associations between air pollution and respiratory mortality: a time-series study. BMJ Open. (2016) 6:e011487. doi: 10.1136/bmjopen-2016-011487

PubMed Abstract | Crossref Full Text | Google Scholar

16. Cengiz MA, Terzi Y. Using of generalized additive model for model selection in multiple poisson regression for air pollution data. Sci Res Essay. (2009) 4:867–71. Available online at: https://academicjournals-org.webpkgcache.com/doc/-/s/academicjournals.org/article/article1380535000Terzi%20and%20Cengiz.pdf

Google Scholar

17. DOF - Diario Oficial de la Federación. NORMA Oficial Mexicana NOM-020-SSA1-2014, Salud ambiental. Valor límite permisible para la concentración de ozono (O3) en el aire ambiente y criterios para su evaluación. (2014). Available online at: https://dof.gob.mx/nota_detalle.php?codigo=5356801&fecha=19%2F08%2F2014#gsc.tab=0 (Accessed August 21, 2024).

Google Scholar

18. DOF - Diario Oficial de la Federación. NORMA Oficial Mexicana NOM-025-SSA1-2014, Salud ambiental. Valores límite permisibles para la concentración de partículas suspendidas PM10 y PM2.5 en el aire ambiente y criterios para su evaluación. (2014). Available online at: https://www.dof.gob.mx/nota_detalle.php?codigo=5357042&fecha=20%2F08%2F2014#gsc.tab=0 (Accessed August 21, 2024).

Google Scholar

19. DOF - Diario Oficial de la Federación. NORMA Oficial Mexicana NOM-023-SSA1-1993, Salud ambiental. Criterio para evaluar la calidad del aire ambiente con respecto al bióxido de nitrógeno (NO2). Valor normado para la concentración de bióxido de nitrógeno (NO2) en el aire ambiente como medida de protección a la salud de la población. (1994). Available online at: https://dof.gob.mx/nota_detalle.php?codigo=4780125&fecha=23%2F12%2F1994#gsc.tab=0 (Accessed August 21, 2024).

Google Scholar

20. Tadano YdS, Ugaya C, Franco AT. Methodology to assess air pollution impact on the population's health using the Poisson regression method. Ambiente Soc. (2009) 12:1–17. doi: 10.5772/33385

Crossref Full Text | Google Scholar

21. Groll A, Tutz G. Regularization for generalized additive mixed models by likelihood-based boosting. Methods Inf Med. 51:168–77. doi: 10.3414/ME11-02-0021

PubMed Abstract | Crossref Full Text | Google Scholar

22. Población con afiliación a servicios de salud por entidad federativa según institución. (2020). Available online at: https://www.inegi.org.mx/app/tabulados/interactivos/?pxq=Derechohabiencia_Derechohabiencia_02_822ebcc5-ef41-40c1-9901-22e397025c64 (Accessed September 5, 2024).

Google Scholar

23. Ortega-García JA, Martínez-Hernández I, Boldo E, Cárceles-Álvarez A, Solano-Navarro C, Ramis R, et al. Contaminación atmosférica urbana e ingresos hospitalarios por asma y enfermedades respiratorias agudas en la ciudad de Murcia (España) [Urban air pollution and hospital admissions for asthma and acute respiratory disease in Murcia city (Spain)]. An Pediatr (Engl Ed). (2020) 93:95–102. doi: 10.1016/j.anpedi.2020.01.012

PubMed Abstract | Crossref Full Text | Google Scholar

24. Martín Martín R, Sánchez Bayle M. Impacto de la contaminación ambiental en las consultas pediátricas de Atención Primaria: estudio ecológico [Impact of air pollution in paediatric consultations in Primary Health Care: Ecological study]. An Pediatr (Engl Ed). (2018) 89:80–5. doi: 10.1016/j.anpedi.2017.06.013

PubMed Abstract | Crossref Full Text | Google Scholar

25. Yamazaki S, Shima M, Ando M, Nitta H, Watanabe H, Nishimuta T. Effect of hourly concentration of particulate matter on peak expiratory flow in hospitalized children: a panel study. Environ Health. (2011) 10:15. doi: 10.1186/1476-069X-10-15

PubMed Abstract | Crossref Full Text | Google Scholar

26. Choo ELW, Janhavi A, Koo JR, Yim SHL, Dickens BL, Lim JT. Association between ambient air pollutants and upper respiratory tract infection and pneumonia disease burden in Thailand from 2000 to 2022: a high frequency ecological analysis. BMC Infect Dis. (2023) 23:379. doi: 10.1186/s12879-023-08185-0

PubMed Abstract | Crossref Full Text | Google Scholar

27. Zhang D, Tian Y, Zhang Y, Cao Y, Wang Q, Hu Y. Fine particulate air pollution and hospital utilization for upper respiratory tract infections in Beijing, China. Int J Environ Res Public Health. (2019) 16:533. doi: 10.3390/ijerph16040533

PubMed Abstract | Crossref Full Text | Google Scholar

28. Fu Y, Zhang W, Li Y, Li H, Deng F, Ma Q. Association and interaction of O3 and NO2 with emergency room visits for respiratory diseases in Beijing, China: a time-series study. BMC Public Health. (2022) 22:2265. doi: 10.1186/s12889-022-14473-2

PubMed Abstract | Crossref Full Text | Google Scholar

29. Wu R, Kang N, Zhang C, Song Y, Liao W, Hong Y, et al. Long-term exposure to PM25 and its components is associated with elevated blood pressure and hypertension prevalence: Evidence from rural adults. J Adv Res. (2024) 60:173–81. doi: 10.1016/j.jare.2023.07.012

PubMed Abstract | Crossref Full Text | Google Scholar

30. Lin YC, Cai Y, Huang HY, Liang D, Li J, Tang Y, et al. Air pollution and blood pressure in the elderly: evidence from a panel study in Nanjing, China. Heliyon. (2022) 8:e10539. doi: 10.1016/j.heliyon.2022.e10539

PubMed Abstract | Crossref Full Text | Google Scholar

31. Li G, Su W, Zhong Q, Hu M, He J, Lu H, et al. Individual PM25 component exposure model, elevated blood pressure and hypertension in middle-aged and older adults: a nationwide cohort study from 125 cities in China. Environ Res. (2022) 215:114360. doi: 10.1016/j.envres.2022.114360

PubMed Abstract | Crossref Full Text | Google Scholar

32. Hayes L, Mejia-Arangure JM, Errington A, Bramwell L, Vega E, Nunez-Enriquez JC, et al. Relationship between air quality and asthma-related emergency hospital admissions in Mexico City 2017-2019. Thorax. (2023) 79:43–9. doi: 10.1136/thorax-2022-219262

PubMed Abstract | Crossref Full Text | Google Scholar

33. Rosenquist NA, Metcalf WJ, Ryu SY, Rutledge A, Coppes MJ, Grzymski JJ, et al. Acute associations between PM25 and ozone concentrations and asthma exacerbations among patients with and without allergic comorbidities. J Expo Sci Environ Epidemiol. (2020) 30:795–804. doi: 10.1038/s41370-020-0213-7

PubMed Abstract | Crossref Full Text | Google Scholar

34. Kelchtermans J, Mentch F, Hakonarson H. Ambient air pollution sensitivity and severity of pediatric asthma. J Expo Sci Environ Epidemiol. (2024) 34:853–60. doi: 10.1038/s41370-023-00573-7

PubMed Abstract | Crossref Full Text | Google Scholar

35. Dondi A, Carbone C, Manieri E, Zama D, Del Bono C, Betti L, et al. Outdoor air pollution and childhood respiratory disease: the role of oxidative stress. Int J Mol Sci. (2023) 24:4345. doi: 10.3390/ijms24054345

PubMed Abstract | Crossref Full Text | Google Scholar

36. Liu K, Hua S. Song L. PM25 Exposure and asthma development: the key role of oxidative stress. Oxid Med Cell Longev. (2022) 2022:3618806. doi: 10.1155/2022/3618806

PubMed Abstract | Crossref Full Text | Google Scholar

37. Shen H, Zheng R, Du M, Christiani DC. Environmental pollutants exposure-derived extracellular vesicles: crucial players in respiratory disorders. Thorax. (2024) 79:680–91. doi: 10.1136/thorax-2023-221302

PubMed Abstract | Crossref Full Text | Google Scholar

38. Enweasor C, Flayer CH, Haczku A. Ozone-induced oxidative stress, neutrophilic airway inflammation, and glucocorticoid resistance in asthma. Front Immunol. (2021) 12:631092. doi: 10.3389/fimmu.2021.631092

PubMed Abstract | Crossref Full Text | Google Scholar

39. medidas que contribuyen a la mejora en la calidad del aire de la megalópolis de la ZMVM. (2018). Available online at: http://www.gob.mx/comisionambiental/articulos/38-medidas-que-contribuyen-a-la-mejora-en-la-calidad-del-aire-de-la-megalopolis-de-la-zmvm (Accessed June 28, 2025).

Google Scholar

40. Programa de Mejoramiento de Calidad de Aire. (2016). Available online at: http://www.gob.mx/comisionambiental/acciones-y-programas/programa-de-mejoramiento-de-calidad-de-aire (Accessed June 28, 2025).

Google Scholar

41. Wang Y, Zhang Z, Hao Z, Eriksson T. Environmental regulation and mental well-being: Evidence from China's air pollution prevention and control action plan. Soc Sci Med. (2025) 365:117584. doi: 10.1016/j.socscimed.2024.117584

PubMed Abstract | Crossref Full Text | Google Scholar

42. García-Burgos J, Miquelajauregui Y, Vega E, Namdeo A, Ruíz-Olivares A, Mejía-Arangure JM, et al. Exploring the spatial distribution of air pollution and its association with socioeconomic status indicators in Mexico City. Sustainability. (2022) 14:15320. doi: 10.3390/su142215320

Crossref Full Text | Google Scholar

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.

Edited by:

Haider Khwaja, University at Albany, United States

Reviewed by:

Worradorn Phairuang, Chiang Mai University, Thailand
Katalin Hubai, University of Pannonia, Hungary

Copyright © 2025 Vega, Ruiz-Olivares, Miquelajauregui, Rentería-Campos, Bramwell, Feytmans, Namdeo, Núñez-Enríquez, Entwistle, Jaimes-Palomera, Torres-Rojas, Rascón-Pacheco, Duarte-Rodríguez, McNally, García-Burgos, Resendiz-Martinez, Fragoso-Chino, Hayes, Portas and Mejia-Arangure. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Elizabeth Vega, ZXZlZ2FAYXRtb3NmZXJhLnVuYW0ubXg=; Juan Manuel Mejia-Arangure, am1lamlhQGlubWVnZW4uZ29iLm14

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.