Estudio estadístico de la correlación entre contaminantes atmosféricos y variables meteorológicas en la zona norte de Chiapas, México

Statistical study of the correlation between atmospheric pollutants and meteorological variables in northern Chiapas, Mexico

S Ramos–Herrera^1*, R Bautista–Margulis, A Valdez–Manzanilla

¹ División Académica de Ciencias Biológicas. UJAT km. 0.5 Carr. Villahermosa – Cárdenas. Villahermosa 86039 Tabasco, México. (SRH). *Correo electrónico: sergio.ramos@dacbiol.ujat.mx

RESUMEN

El análisis de regresión múltiple es un método estadístico empleado en muchas áreas del conocimiento. En este estudio, dicho análisis se aplicó a los datos de concentraciones de cuatro contaminantes atmosféricos (SO₂, NO₂, H₂S y PM₁₀), monitoreados en tres estaciones que se ubican en la Zona Norte del estado de Chiapas. El periodo que abarcó el estudio fue de enero 2001 a febrero 2005. El objetivo fue proponer funciones de regresión para describir la concentración en función del tiempo y/o las variables meteorológicas. Se empleó un análisis de regresión lineal múltiple paso a paso en la selección de variables regresoras. Las más importantes fueron la temperatura, la humedad relativa y la dirección del viento. Se obtuvieron funciones de regresión de la concentración anual, mensual y diaria de estos contaminantes. Se obtuvo una regresión lineal simple para explicar la concentración anual de SO₂, alcanzando un coeficiente de determinación de 0.927. Los modelos de la concentración mensual alcanzaron un coeficiente de determinación entre 0.417 y 0.846; mientras que para los de la concentración diaria, este coeficiente varió de 0.285 a 0.581. De este estudio paramétrico, se concluyó que las variables meteorológicas describieron adecuadamente la concentración anual y mensual, pero no la concentración diaria.

Palabras clave: Correlación estadística, regresión lineal, contaminantes atmosféricos.

ABSTRACT

The multiple regression analysis is a statistical method used in many academic areas. In this study, this analysis was applied to the concentration data of four atmospheric pollutants (SO₂, NO₂, H₂S y PM₁₀), monitored in three stations located in northern Chiapas. The study was carried out from January 2001 to February 2005. The purpose of the study was to propose regression functions to describe the concentration as a function of time and/or meteorological variables. A stepwise multiple linear regression analysis was used in the selection of regressive variables, the most important of which included the temperature, relative humidity and wind direction. Regression functions of the annual, monthly and daily concentrations of these pollutants were obtained. A simple linear regression was obtained to explain the annual concentration of SO₂, with a determination coefficient of 0.927. The models for the monthly concentration obtained determination coefficients of 0.417 to 0.846,whereas the coefficients for the daily concentration model varied from 0.285 to 0.581. This parametric study made it possible to conclude that the meteorological variables described the annual and monthly concentrations adequately, but not the daily concentrations.

Key words: Statistic correlation, linear regression, atmospheric pollutants.

INTRODUCCIÓN

La concentración de los contaminantes en la atmósfera se afecta por variables meteorológicas como, la temperatura, la humedad relativa y la radiación solar, que controlan a su vez la velocidad con la que se realizan las reacciones químicas atmosféricas. Por otro lado, la dispersión de estos contaminantes es influenciada por la velocidad y dirección de los vientos, así como la estabilidad que predomina en la atmósfera. A nivel regional se ha puesto atención al problema de la contaminación atmosférica, debido principalmente a las actividades petroleras que se realizan en el suroeste de México. A través de una red de monitoreo automatizada conocida como Sistema Automatizado de la Región Sur (SAMARS) de Petróleos Mexicanos (PEMEX) se monitorean los principales contaminantes atmosféricos (excepto el ozono) y variables meteorológicas. Esta red opera desde 1999 y cuenta con seis estaciones de monitoreo ubicadas en las periferias de las instalaciones petroleras (Baterías y Compresoras). La información recopilada a la fecha se ha empleado principalmente en la evaluación de la calidad del aire en la periferia de estas instalaciones, en la calibración de modelos de dispersión de los contaminantes, y en la distribución espacial de estos. Por ejemplo, Valdés et al. (2004) realizaron un análisis de la distribución espacial y estacional del dióxido de azufre. El análisis se llevó a cabo con datos de once meses (Noviembre de 2002 a Septiembre de 2003) en cinco estaciones de monitoreo atmosférico del SAMARS distribuidas en la Zona Norte de Chiapas (ZNC) y centro de Tabasco. De la distribución espacial se concluyó que la máxima concentración de SO₂ se presentó en la ZNC, concretamente en la zona de Reforma Chiapas, y que el máximo estacional ocurrió en el otoño. Sin embargo, se recomendó continuar con este tipo de estudios, abarcando un periodo mayor de análisis que incluyera además, variables meteorológicas como la velocidad, dirección del viento y la estabilidad atmosférica.

Existen estudios que muestran que los contaminantes atmosféricos y las variables meteorológicas se correlacionan. En un estudio realizado en la Zona Metropolitana de la Ciudad de México (ZMCM), Bravo et al. (2000) mostraron que la radiación solar es más útil en la modelación del máximo diario de ozono que la temperatura. Mediante un análisis descriptivo comprobaron que aunque ambas variables meteorológicas se correlacionan con el máximo diario de ozono, la temperatura máxima se presentó casi dos horas después de la hora del máximo de ozono en el día; mientras que la radiación solar máxima se presentó entre 1.69 y 2.35 h antes de la hora del máximo de ozono. Sin embargo, no propuso un modelo de los niveles diarios de ozono en función de la radiación solar máxima.

También existe correlación entre la humedad relativa (HR), la concentración de partículas suspendidas totales (PST) y las partículas de diámetro aerodinámico menor a 10 μm (PM₁₀), como lo muestran los resultados de una investigación realizada en el suroeste de la Ciudad de México, en una zona de tipo residencial. Ahí se encontró que la correlación es de tipo inversa y que es mayor con las PST (r = –0.828) que con las PM₁₀ (r = –0.749). También se estableció que existe diferencia significativa entre la concentración en temporada de secas y de lluvias, siendo mayor la concentración en la temporada de secas. Finalmente, se puso de manifiesto la importancia de estudiar el comportamiento anual de la contaminación ya que es importante en las exposiciones por periodos prolongados (Amador–Muñoz et al. 2001). La determinación de los patrones de distribución diurna y anual de los contaminantes ha permitido concluir si estos presentan ciclos diarios, tendencias a la disminución a largo plazo y comportamiento estacional. Para ello se ha analizado el comportamiento de estadísticas como la media aritmética, la mediana y el percentil 98 de los valores diarios, mensuales y anuales del contaminante Davydova–Belitskaya et al. 1999; Varó & Carratalá 2002).

Debido a la naturaleza aleatoria de los niveles atmosféricos de los contaminantes, se ha intentado explicar su comportamiento temporal y en el proceso se ha recurrido a los modelos de regresión. Por ejemplo, el comportamiento de los máximos y promedios semanales de ozono en la Zona Metropolitana de la Ciudad de México (ZMCM) con base en datos de ocho estaciones de monitoreo atmosférico de la RAMA (en el periodo enero de 1995 a diciembre de 1999) se explicó mediante el siguiente modelo de regresión no lineal:

Donde x_i, representa el tiempo transcurrido en semanas, y es la única variable regresora.

El modelo integró una componente de tendencia lineal, una componente cíclica y un componente aleatorio, ε_i, con las que se modeló la estacionalidad anual y la tendencia lineal creciente o decreciente del ozono. Aunque el modelo ajustado fue bastante bueno, se sugirió incluir variables meteorológicas en el análisis, con el propósito de eliminar su efecto en la variación del ozono y así evaluar mejor el impacto que las estrategias de control tienen sobre los datos de ozono observados (Jiménez, 2002).

En estudios de corte epidemiológico los niveles diarios de contaminación (por O₃, PM₁₀, SO₂, humos negros) se han empleado como variables explicativas de variables de salud (cómo el número de consultas a urgencias por enfermedades respiratorias agudas o del número de defunciones diarias por todas las causas de muerte o por enfermedades del aparato respiratorio y del aparato circulatorio). Aunque el objetivo principal ha sido estudiar la asociación entre estas variables, en el proceso de análisis, las variables meteorológicas (como la temperatura, humedad relativa y velocidad del viento) se han considerado como variables confusoras cuando se emplean métodos de series temporales en el análisis de la una variación temporal de la variable salud y contaminación diaria (Sáez et al. 1999). También han intervenido como variables predictoras en la estimación de los niveles faltantes de los contaminantes mediante la aplicación de técnicas de regresión lineal (Hernández–Cadena et al. 2000; Arrivas–Monzón et al. 2001).

MATERIALES Y MÉTODOS

El análisis se llevó a cabo en la ZNC (Figura 1) en una región ubicada entre los meridianos 92° 49' y 93° 34' y los paralelos 18° 34' y 17° 33'. Esta zona presenta un clima cálido húmedo con lluvias en verano; la temperatura promedio anual oscila alrededor de los 24 ^°C y la máxima precipitación anual del estado es de 3 977.5 mm (Anónimo 2003).

La base de datos se obtuvo de tres estaciones (Reforma, Giraldas y Artesa) de monitoreo de la red SAMARS. La distancia entre las estaciones de Reforma y Artesa es de 10.16 km, entre Artesa y Giraldas es de 25.94 km y entre Reforma y Giraldas es de 33.42 km. La estación Reforma, se ubica en la zona industrial de la ciudad de Reforma, Chiapas; la estación Artesa, se ubica en el poblado de Santa Teresa en un entorno rural; y la estación Giraldas, se ubica en el Campo San Manuel, una zona de tipo industrial, perteneciente al municipio de Huimanguillo, Tabasco (Tabla 1). De las tres estaciones, la primera y la última cuentan además, con una estación meteorológica. El estudio abarcó el periodo de enero 2001 a febrero de 2005. A partir de esta última fecha las estaciones quedaron fuera de operación y así continúan en la actualidad. Para lograr el objetivo de esta investigación, se recopilaron los valores horarios de los siguientes contaminantes y variables meteorológicas: dióxido de azufre (SO₂), dióxido de nitrógeno (NO₂), ácido sulfhídrico (H2S), partículas con diámetro aerodinámico menor a 10 micras (PM₁₀), velocidad del viento (VV), dirección del viento (DV), temperatura (T), humedad relativa (HR), presión barométrica (P) y radiación solar (RS). De los resultados del trabajo de Vázquez (2006) se obtuvieron los valores diarios de la altura de la capa de mezclado (L) en la región en el 2003. Las concentraciones se manejaron ppb para los gases y en μgm^–3 para las partículas. La temperatura en °C, la presión en mb, la radiación solar en Wm^–2 y la altura de la capa de mezclado en m.

Las medias aritméticas anuales, mensuales y valores diarios de los contaminantes y meteorológicos se calcularon con el 50% o más de los valores horarios correspondientes. Se tuvieron cuatro medias anuales, 50 medias mensuales y 1 520 valores diarios por variable.

Particularmente los valores diarios de concentración de NO2 y H2S se calcularon de la siguiente forma: para el NO2 el valor diario se tomó como el máximo valor horario del día; mientras que para el H2S, el valor diario se tomó como el máximo promedio móvil de ocho horas del día. En los cálculos de la media de velocidad y dirección del viento se usó el enfoque escalar y el enfoque vectorial (Anónimo 2000), generándose medias escalares (VVe, DVe) y medias vectoriales (VVv, DVv) para estos parámetros. Se realizó un análisis de correlación entre los contaminantes y las variables meteorológicas; primero entre los valores anuales, luego entre los valores mensuales y finalmente entre los valores diarios. Dicho análisis dio un indicio de qué variables pudieran servir como variables regresoras de la concentración. El modelo de regresión propuesto fue el siguiente:

Para ajustar el modelo, se empleó el paquete estadístico SPSS y el método de selección de variables STEPWISE (Ferrán 2001). Se analizaron los residuales para verificar el grado en el que se cumplieron las hipótesis de homogeneidad de varianzas y normalidad de los residuales (Montgomery et al. 2006). El modelo incluyó a la dirección del viento como una variable de rango y como una variable categórica, para ver cual funcionaba mejor. Como variable categórica, la dirección del viento incluyó las 16 direcciones siguientes: N, NNE, NE, ENE, E, ESE, SE, SSE, S, SSO, SO, OSO, O, ONO, NO, NNO. Esto se logró definiendo la variable indicadora binaria Z_i, que toma el valor de uno si la dirección del viento es igual a i, y cero si la dirección del viento es diferente a i; donde i es igual a las direcciones anteriores del viento menos una. También se analizó la correlación entre la concentración y la dirección del viento. Por último, se estudió la correlación diaria entre las tres estaciones y se ajustaron rectas de regresión lineal simple en los casos donde la correlación fue significativa y alta. Para medir las correlaciones, se empleó el coeficiente de correlación de Spearman. El modelo no incluyó términos para el comportamiento estacional de las variables meteorológicas debido a la necesidad de utilizar mayor información de campo y, por consiguiente, incrementar la complejidad del análisis.

RESULTADOS

Comportamiento anual del SO₂

El SO2 monitoreado en Reforma fue el único contaminante con suficiente información horaria (mayor al 50 %) para estimar las medias anuales 2001–2004. Los resultados (Figura 2) muestran una tendencia a la disminución de la media anual con el tiempo en años. Dicha tendencia se representó por el modelo de regresión simple: [SO₂] = 14.10 –2.362t, donde t es el tiempo transcurrido en años. De acuerdo a este, para el 2005 (t = 5 años) se espera una media anual de 2.29 ppb. No obstante la media observada fue de 5 ppb, quedando dentro de los límites del intervalo de predicción para este año.

Comportamiento mensual y diario de los contaminantes

El comportamiento de la concentración mensual se muestra en las Figuras 3, 4 y 5 para las estaciones de Reforma, Artesa y Giraldas, respectivamente. Ninguna de las estaciones registró el 100% de las medias mensuales dentro del periodo de análisis (50 meses). Los casos críticos se dieron en Artesa y Giraldas y para los contaminantes NO₂y H₂S, donde el porcentaje alcanzado fue menor a 50%.

Los niveles diarios de los contaminantes durante el periodo de estudio presentaron mayor variabilidad que los valores diarios de las variables meteorológicas. Para mejorar la homogeneidad de varianzas y la normalidad de los valores diarios, los niveles diarios se transformaron mediante el logaritmo natural. El comportamiento diario de SO₂, NO₂, H₂S y PM₁₀(Figuras 6, 7, 8 y 9) para cada estación de monitoreo se representó mediante gráficos del logaritmo natural de la concentración diaria del contaminante en función de los días transcurridos desde enero de 2005.

De acuerdo a la NOM–022–SSA1–1993 (Anónimo 1993a), el límite máximo permisible de SO₂en un promedio de 24 h es de 0.13 ppm ó 130 ppb; para NO₂es de 0.21 ppm ó 210 ppb en un promedio de una hora según lo establecido en la NOM–023–SSA1–1993 (Anónimo 1993b); y para PM₁₀ fue de 150 μgm^–3 hasta antes de septiembre de 2005 y después de esta fecha se estableció en 120 μgm^–3 acorde a la NOM–025–SSA1–1993 (Anónimo 1993c). En cuanto al H₂S, se han recomendado concentraciones máximas de exposición aguda de 42 μgm^–3 < 1 h y exposición crónica de 10 μgm^–3 (Anónimo 2008). En base a estos límites, los máximos alcanzados en algún día dentro del periodo de estudio (Tabla 2) que excedieron la norma fueron el SO₂en las tres estaciones, el NO₂en Reforma y Giraldas y las PM₁₀ en Artesa y Giraldas. Los máximos de H₂S excedieron los valores recomendados por exposición aguda o crónica en todas las estaciones.

El análisis de correlaciones mensuales entre los contaminantes y las variables meteorológicas se hizo solo en los casos de las estaciones de Reforma y Giraldas, pues Artesa no cuenta con estación meteorológica. La concentración mensual se correlacionó negativamente con el tiempo t, a excepción de las PM₁₀ observadas en Reforma cuya correlación fue positiva aunque baja y no significativa. La correlación de NO₂ y PM₁₀ con la humedad relativa (HR) fue negativa; y positiva con la temperatura (T) y la radiación solar (RS) en ambas estaciones de monitoreo. La presión (P) se correlacionó negativamente con el NO2 en las dos estaciones. La velocidad escalar del viento (VVe) se correlacionó positivamente solo con el NO₂ y las PM₁₀ en ambas estaciones, mientras que la dirección escalar del viento (DVe) se correlacionó significativamente y de forma negativa con las PM₁₀ de Reforma. La altura de la capa de mezcla (L) se correlacionó negativamente con el NO₂ y las PM₁₀ (Tabla 3). En todos los resultados obtenidos la correlación significativa al 5% se indicó con un asterisco y al 2.5 % se indicó con dos asteriscos como un superíndice en el valor.

Las concentraciones diarias en Reforma y Giraldas se correlacionaron negativamente con el tiempo (Tabla 4). La humedad relativa se correlacionó negativamente más con las PM₁₀ que con el NO₂. La temperatura se correlacionó positivamente con los cuatro contaminantes en Reforma pero solo con el NO₂ y las PM₁₀ en Giraldas. De modo general, la radiación solar se correlacionó positivamente mientras que la presión lo hizo negativamente con los contaminantes en las dos estaciones. A diferencia de los valores mensuales, los valores diarios de la contaminación se correlacionaron significativamente con la velocidad y dirección escalar del viento, aunque no para todos los contaminantes y ambas estaciones. Algo similar ocurrió con la altura diaria de la capa de mezclado, la correlación fue negativa y significativa en las dos estaciones aunque más baja que la correlación a escala mensual (Tabla 4). Esto tiene sentido ya que al aumentar la altura de la capa de mezcla, aumenta el volumen de aire disponible para dilución lo cual disminuye la concentración de los contaminantes.

En el análisis del efecto de la dirección categórica del viento en la concentración se observó la existencia de una correlación mediana a escala mensual y baja a escala diaria, que en algunos fue significativa y en otros no. Por ejemplo, a nivel mensual, el SO₂ en Reforma se correlacionó (0.390) positivamente con los vientos del noreste, mientras que las PM₁₀ lo hicieron de modo inverso (–0.377) con los vientos del sureste. En Giraldas el NO₂ se correlacionó positivamente con los vientos del sur. A nivel diario se encontró que la dependencia entre la concentración de los contaminantes y la dirección del viento es variada (Tabla 5). El SO₂, NO₂ y H₂S en Reforma se correlacionaron negativamente con el viento del este, que es el viento local que predomina tanto en Reforma como en Giraldas (Tabla 6). En Giraldas, solo el H₂S se correlacionó con el viento dominante. También se observó correlación significativa con los segundos vientos prevalentes, que resultaron ser en ambos casos los vientos del este noreste.

Correlación entre estaciones

Los gráficos de dispersión de los valores diarios de SO₂, NO₂y H₂S no parecen mostrar ningún patrón de asociación (Figuras 11, 12, 13 y 14) entre las estaciones. Sin embargo, el coeficiente de correlación fue significativo y varió entre 0.1 y 0.5. La asociación más fuerte se dio para el caso de las PM₁₀ (Figura 12) donde el coeficiente de correlación resultó significativo y mayor a 0.6.

Funciones de regresión

La Tabla 7 presenta las ecuaciones de regresión que describen la concentración mensual de los contaminantes para las estaciones de Reforma y Giraldas. La Tabla 8 presenta las funciones de regresión ajustadas al logaritmo natural de la concentración diaria y en ambos casos R²_ajus es el coeficiente de determinación ajustado por el número de observaciones y el número de variables independientes. Fue necesario transformar las concentraciones diarias a logaritmo de las concentraciones para cumplir moderadamente con las hipótesis de homogeneidad de varianza y normalidad de los residuales de los modelos.

La Tabla 9 muestra las funciones de regresión ajustadas en el análisis de correlación entre estaciones, donde R² es el coeficiente de determinación simple. Como puede notarse son modelos de regresión simple que explican una buena parte de la variabilidad observada entre variables del mismo tipo.

DISCUSIÓN

En los últimos años se han instalado plantas recuperadoras de azufre y quemadores ecológicos en las instalaciones petroleras, con el objetivo de reducir las tasas de emisión de contaminantes a la atmósfera, entre ellos el SO₂. Esto pudiera explicar la tendencia en la disminución de la concentración media anual de SO₂. Un resultado semejante encontró Varó & Carratalá (2002) en donde el periodo en el que se dan las disminuciones de humos negros y SO₂ en ALCOY, coinciden con la puesta en funcionamiento de una red de gas natural en las industrias y que reduce las emisiones debido a actividades industriales.

El grado de asociación entre los contaminantes y las variables meteorológicas a escala mensual (que oscila entre 0.373 y 0.772) es significativamente mayor que las que reporta Correa et al. (2002) quienes analizaron la relación entre los máximos diarios de ozono y los promedios diarios de las variables meteorológicas (menor de 0.3). Las variables que resultaron mejores descriptoras de la concentración mensual fueron la temperatura, la humedad relativa y la dirección de los vientos. En los modelos ajustados no se incluyeron simultáneamente estas variables, ya que están correlacionadas fuertemente (mayor a 0.6) y dejan modelos con problemas fuertes de multicolinealidad. Aunque la correlación entre la dirección categórica de los vientos y la concentración mensual de los contaminantes es baja (< 0.3), esta variable sí contribuyó a la mejora de los modelos, lo cual se manifestó en un incremento significativo en el coeficiente de correlación múltiple corregido. En Reforma, este coeficiente osciló entre 0.548 y 0.704 y en Giraldas osciló entre 0.533 y 0.846. Estos valores son aceptables; por ejemplo Correa (2004) ajustó ecuaciones de regresión para explicar la concentración media anual de los contaminantes en la ZMVM en un periodo de seis años con un coeficiente de determinación ajustado que varía entre 0.3 y 0.9 típicamente.

A nivel diario, la humedad relativa y la temperatura resultaron inadecuadas en la descripción del comportamiento de los contaminantes SO₂ y H₂S. Para mejorar los modelos, en el caso del SO₂, pudiera incluirse como variable regresora, la tasa de emisión desde las fuentes, que como ya se comentó en los últimos años ha ido disminuyendo paulatinamente. Para mejorar el modelos del H₂S quizás sea necesario incluir al SO₂ como una variable regresora, ya que de los resultados han mostrado una fuerte correlación entre estos contaminante (hasta de 0.5). Estas mismas variables resultaron adecuadas en la descripción de la contaminación diaria por NO₂ y PM₁₀. Los modelos ajustados alcanzaron coeficientes de determinación múltiple entre 0.2 y 0.6. Sin embargo, hay que tomar en cuenta que los modelos de regresión se ajustaron al logaritmo de las concentraciones diarias y al retornar a las concentraciones originales, la incertidumbre del modelo cambia.

Otro de los resultados interesantes en este estudio, es la existencia de moderada y fuerte correlación en las mediciones de PM₁₀, temperatura, humedad relativa, radiación solar y presión entre las estaciones de monitoreo. Los coeficientes de correlación múltiple de las rectas de regresión ajustadas a estas variables con relación a lo observado en Reforma varían entre 0.631 y 0.974; mientras que para las PM₁₀ oscilan entre 0.328 y 0.820. Esta aparente homogeneidad puede deberse a que las condiciones ambientales en esas zonas son aparentemente similares, ya que ambas casetas se ubican dentro de áreas Talleres de PEMEX. Lo mismo puede decirse de la presión y la radiación solar; ya que en el caso de la presión barométrica, se necesitan diferencias significativas en altitud para obtener lecturas que difieran entre sí. Quizás es menos obvio en el caso de la radiación solar y la temperatura.

De manera genérica se puede decir que la recta de regresión ajustada a la media anual de SO₂ es adecuada. La tendencia fuertemente negativa de la concentración en este contaminante con el tiempo puede deberse a la disminución en las emisiones de éste, presumiblemente debido a la puesta en marcha de plantas recuperadoras de azufre en las instalaciones petroleras.

Las ecuaciones de regresión ajustadas para explicar las concentraciones mensuales de los contaminantes tienen una bondad de ajuste bastante buena. Puede pensarse, que las variables meteorológicas mensuales como la humedad relativa y la temperatura son adecuadas en la descripción de la tendencia mensual en la contaminación atmosférica de la ZNC. Se recomienda en el futuro simular la distribución espacial y temporal de estos contaminantes usando modelos de dispersión y de pronóstico meteorológico (MM5), que tomen en cuenta la circulación a nivel regional, los vientos dominantes y sus variaciones estacionales. El análisis entre estaciones meteorológicas, indicó que probablemente están registrando información muy similar por lo que la falta de información en una estación pudiera complementarse por la otra.

AGRADECIMIENTOS

Se agradece a la Gerencia de Seguridad Industrial y Protección Ambiental de PEMEX Exploración y Producción Región Sur en el estado de Tabasco por la información proporcionada para el desarrollo de este estudio.

Anónimo (1993a) NOM–022–SSA1–1993; Salud ambiental. Criterio para evaluar la calidad del aire ambiental con respecto al bióxido de azufre (SO₂). Valor normado para la concentración de bióxido de azufre (SO₂) en el aire ambiente como medida de protección a la salud de la población. 4pp.         [ Links ]

Anónimo (1993b) NOM–023–SSA1–1993; Salud ambiental. Criterio para evaluar la calidad del aire ambiental con respecto al bióxido de nitrógeno (NO₂). Valor normado para la concentración de bióxido de nitrógeno (NO₂) en el aire ambiente como medida de protección a la salud de la población. 3pp.         [ Links ]

Anónimo (1993c) NOM–025–SSA1–1993; Salud ambiental. Criterios para evaluar la calidad del aire ambiente con respecto a material particulado. 11p.         [ Links ]

Anónimo (2000) Meteorological monitoring guidance for regulatory modeling applications. Environmental Protection Agency Oficina de Aire y Radiación Oficina de Planeación de Calidad del Aire y Estándares. EPA–454/R–99–005. Cap. 6. 7 pp.         [ Links ]

Anónimo (2003) Anuario estadístico edición 2003 Chiapas. Instituto Nacional de Estadística Geografía e Informática (INEGI). 606 pp.         [ Links ]

Anónimo (2008).Office of Environmental Health Hazard Assessment (OEHHA) and the Air Reosurces Borrad (ARB) for use 2588 in facility helath risk assessments conducted for the AB 2588 Air Toxics Hot Spots Program.         [ Links ]

Arrivas–Monzón F, Rabanaque MJ, Martos MC, Abad JM, Alcalá–Nalvaíz T, Navarro–Elipe M (2001) Efectos de la contaminación atmosférica sobre la mortalidad diaria en la ciudad de Zaragoza, España, 1991–1995. Salud Pública de México 43(4) 289–297.         [ Links ]

Amador–Muñoz O, Delgado–Rodríguez A, Villalobos–Pietrini R, Munive–Colín Z, Ortíz–Marttello R, Díaz–Gonzáles G, Bravo–Cabrera JL, Gómez–Arroyo S (2001) Partículas suspendidas, hidrocarburos aromáticos policíclicos y mutagenicidad en el Suroeste de la Ciudad de México. Contaminación Ambiental 17(4):193–204.         [ Links ]

Bravo JL, Nava MM, Muhlia A (2000) Relaciones entre la magnitud del valor máximo de ozono, la radiación solar y la temperatura ambiente en la Zona Metropolitana de la Ciudad de México. Contaminación Ambiental 16(2): 45–44.         [ Links ]

Correa AG, Muñoz RC, Cardoso GC (2002) Análisis de ozono y óxidos de nitrógeno mediante técnicas de series de tiempo y regresión lineal. Secretaría del Medio Ambiente Dirección General de Prevención y Control de la Contaminación México, D.F. 14pp.         [ Links ]

Correa AG (2004) Contaminantes atmosféricos en la zona metropolitana de la ciudad de México. Instituto Politécnico Nacional Universidad Autónoma Metropolitana México 280p.         [ Links ]

Davydova–Belitskaya V, Skiba YN, Bulgakov SN, Martínez A (1999) Modelación matemática de los niveles de contaminación en la ciudad de Guadalajara, Jalisco, México. Parte I, microclima y monitoreo de la contaminación. Revista Internacional de Contaminación Ambiental. 15 (2): 103–111.         [ Links ]

Ferrán MA (2001) SPSS para Windows análisis estadístico. España MCGRAW–HILL. 421 pp.         [ Links ]

Hernández–Cadena L, Téllez–Rojo MM, Sanin–Aguirre LH, Lacasaña–Navarro M, Campos A, Romieu I (2000) Relación entre consultas a urgencias por enfermedad respiratoria y contaminación atmosférica en Ciudad Juárez, Chihuahua. Salud Pública de México. 42(4): 288–297.         [ Links ]

Jiménez JCH (2002) Análisis de valores extremos de la concentración de ozono en la ciudad de México Tesis profesional. UJAT División Académica de Ciencias Básicas Cunduacán Tabasco 95p.         [ Links ]

Montgomery DC, Peck EA, Vining GG (2006) Introducción al análisis de regresión lineal (3a. ed.) México D.F. CECSA. 588 pp.         [ Links ]

Sáez M, Pérez–Hoyos S, Tobías A, Saurina C, Barceló MA, Ballester F (1999) Métodos de series temporales en los estudios epidemiológicos sobre contaminación atmosférica. Rev Esp Salud Pública. 73(2): 133–143.         [ Links ]

Varó PG, Carratalá AG (2002) Evolución de los niveles de inmisión de contaminación atmosférica en una ciudad industrial (ALCOY) desde 1989 a 2000. Revista Salud Ambiental, 2(1): 8–15.         [ Links ]

Valdés AM, Bautista RM, Ramos SH, Magaña EV (2004) Distribución espacial y estacional de las inmisiones de bióxido de azufre (SO₂) en la zona norte de Chiapas y centro de Tabasco. KUXULKAB' Revista de Divulgación. División Académica de Ciencias Biológicas. Universidad Juárez Autónoma de Tabasco. 10(19): 33–36.         [ Links ]

Vázquez CA (2006) Aplicación de un modelo de dispersión para determinar las concentraciones de monóxido de carbono (CO) emitidos por motores diesel de un pozo petrolero. Tesis profesional Universidad Juárez Autónoma de Tabasco División Académica de Ciencias Biológicas Villahermosa Tabasco 106p.         [ Links ]