Efectos a corto plazo de la contaminación atmosférica sobre la mortalidad: Resultados del proyecto EMECAM en Valencia, 1994-96 (página 2)

MATERIAL Y MÉTODOS

Ámbito

La ciudad de Valencia tiene en la actualidad una población aproximada de 750.000 habitantes (censo de 1996). Su situación a orillas del mar Mediterráneo le confiere un clima de tipo mesotermal, con inviernos templados y veranos secos y calurosos1.

Mortalidad

Para el proyecto EMECAM se seleccionaron diferentes indicadores de mortalidad, todos ellos obtenidos del registro de mortalidad de la Comunidad Valenciana. Los casos seleccionados se refieren a aquellas personas con residencia en la ciudad y cuya defunción ocurrió en la misma. Los indicadores fueron a) mortalidad por todas las causas, menos las externas, en todas las edades; b) mortalidad por todas las causas, menos las externas, en personas de 70 ó más años, c) mortalidad por enfermedades respiratorias y d) mortalidad por enfermedades del aparato circulatorio. Se dispuso de las series diarias desde el 1 de enero de 1994 hasta el 31 de diciembre de 1996.

Contaminación atmosférica

Valencia dispone en la actualidad de dos redes de vigilancia de la contaminación atmosférica, una red manual con 22 captadores y otra automática con 5, ambas gestionadas por el Ayuntamiento de la ciudad. En la red manual se monitorizan los niveles diarios de humos negros (British smoke, determinados por el método reflectométrico) y de SO2 (método de la Thorina). La red automática esta formada por 5 captadores y mide los niveles de SO2 (fluorescencia ultravioleta), NO2 (quimioluminiscencia), CO (absorción infrarroja) y O3 (absorción ultravioleta). De la red manual se obtuvieron los niveles medios diarios (24 h) de humos negros, descartando los niveles de SO2 por el método de la Thorina, por disponer de la medida automatizada. De la red automática se calcularon tanto los niveles medios diarios como el valor máximo horario (1hora) de SO2, NO2 y CO. Para el ozono se calcularon los niveles máximos de 8 horas (franja horaria de 8 horas con el máximo nivel de O3). Siguiendo los criterios de proyecto EMECAM4 se seleccionaron solamente aquellos captadores que tuvieran al menos un 75% de lecturas válidas para toda la serie. Con este criterio se incluyeron, finalmente, 7 captadores manuales y 4 estaciones automáticas (figura 1). Las correlaciones entre los diferentes captadores varió entre 0,30 y 0,62 en la red manual y entre 0,39 y 0,76 en la automática.

Figura 1 Mapa de los Distritos Municipales de la ciudad de Valencia en el que se muestra la distribución de los captadores manuales (círculos) y estaciones automáticas (triángulos) de la Red de Vigilancia de la Contaminación Atmosférica seleccionados en el estudio. La estrella representa el emplazamiento del Centro Meteorológico Zonal de Valencia.

Variables meteorológicas

Las dos variables seleccionadas en el proyecto EMECAM fueron la temperatura media (media de la temperatura mínima y máxima de cada día) y la humedad relativa diaria (media de la humedad a las 7, 13 y 18 horas). En el caso de esta última no se dispuso de las lecturas correspondientes a las 0 horas, por lo que los niveles calculados pueden infraestimar los reales. Ambas variables meteorológicas se obtuvieron en el Centro Meteorológico Zonal de Valencia, instalado dentro del casco urbano de la ciudad (figura 1).

Otras variables

La incidencia de gripe se obtuvo del sistema de Enfermedades de Declaración Obligatoria (EDO), asignando a cada día la séptima parte del número de casos declarados durante la semana. La cobertura del sistema de declaración en Valencia se ha mostrado como aceptable5.

Entre los efectos calendario incluidos en los modelos basales, aparte de los días de la semana, se comprobó la influencia de los días de fiesta (diferentes al domingo) o de eventos especiales como la huelga de médicos, que en nuestra ciudad se desarrolló entre el 16 de mayo y el 20 de junio de 1995. Tanto para los días de fiesta como para la huelga se creó una variable indicadora (dummy) que tomaba el valor 1 en los días en los que se produjo el acontecimiento y 0 en el resto de la serie.

Análisis estadístico

De acuerdo con la metodología del proyecto EMECAM4 se construyeron modelos de regresión autoregresiva de Poisson, controlando por los diferentes factores de confusión (estacionalidad, tendencia, calendario, variables meteorológicas e incidencia de gripe).

RESULTADOS

En la tabla 1 se resumen los valores de los diferentes indicadores de mortalidad. La mediana del numero diario de muertes por todas las causas menos las externas fue 16 (rango de 5 a 34 defunciones). De éstas, como era de esperar, una gran parte afectaron a personas con 70 ó más años. Por causas específicas, fueron más numerosas las muertes por enfermedades del aparato circulatorio que por respiratorias. La estacionalidad de las series de mortalidad mostró un mayor número de casos durante el semestre frío que en el cálido.

Tabla 1 Estadísticos descriptivos de los indicadores de mortalidad (casos diarios), variables meteorológicas, incidencia de gripe y contaminantes. Valencia, periodo 1994-1996

Los niveles de contaminación atmosférica en Valencia para el periodo de análisis se situaron en general dentro de los estándares de calidad del aire establecidos por la OMS en 19876 (Tabla 1). Como excepción, los niveles de humos negros durante el semestre más frío superaron el valor guía de 125 µg/m3 en siete días (0,6 % de la serie) y el NO2 24h (promedio de 24 horas) el umbral de 150 µg/m3 un día (0,09%). La estacionalidad de las series de humos negros, SO2 y CO, fue similar a las de mortalidad, con niveles más altos durante la época más fría del año. Por el contrario, el ozono presentó una estacionalidad opuesta con niveles mayores durante los meses más cálidos. El NO2 no presentó una estacionalidad identificable.

Según se muestra en la tabla 2 los niveles de humos negros, SO2 y CO presentaron una alta correlación lineal positiva entre ellos. El O3 se correlacionó negativamente con el resto de contaminantes y el NO2 mostró una correlación lineal bastante pobre con el resto de las series. La temperatura se relacionó positivamente con el ozono y negativamente con los humos negros, SO2 y CO.

Tabla 2 Coeficientes de correlación de Pearson entre las variables meteorológicas y contaminantes atmosféricos

Los residuos de los modelos de regresión de Poisson construidos para cada indicador de mortalidad no presentaron autocorrelaciones significativas excepto en la serie de mortalidad por todas las causas, evidenciando algún problema de especificación, que se solucionó introduciendo términos autorregresivos7.

En la tabla 3 se muestran los resultados para cada contaminante, expresados como el riesgo relativo (RR) de mortalidad para un incremento en 10 µg/m3 (1 mg/m3 para CO). Los efectos más significativos correspondieron a humos negros, SO2 1h, NO2 1h y CO 24h, para los indicadores de mortalidad por todas las causas con mayor número de casos (todas las edades y mayores de 70 años). Los retardos seleccionados fueron aquellos con una estimación más precisa, tal y como se puede comprobar en los gráficos de correlaciones cruzadas (CCF) de la figura 2. Los coeficientes de correlación de los retardos negativos (parte izquierda del gráfico) indican la influencia del contaminante en el día indicado (0 el mismo día, -1 el día anterior y así sucesivamente) sobre la mortalidad diaria. La parte derecha (retardos positivos) indicaría la hipotética influencia de la mortalidad sobre la contaminación, por lo que deben de permanecer dentro de las bandas de no significación estadística.

Figura 2 Funciones de correlaciones cruzadas (CCF) entre la mortalidad por todas las causas y los promedios diarios de humos negros (a), y de monóxido de carbono (b), en Valencia, periodo 1994-1996, después de ajustar por el modelo basal correspondiente.

Tabla 3 Asociación entre contaminación atmosférica y mortalidad en Valencia, periodo 1994-1996.

En el análisis estratificado por semestres, los resultados fueron similares a los de todo el periodo. Sin embargo, para la mortalidad por enfermedades circulatorias se observó una modificación de efecto significativa del CO con respecto al semestre de estudio. El CO 24h, retardo 5, presentó un efecto positivo durante la época más cálida (RR 1 mg: 1,068; IC95%: 1,007 – 1,133) y que no se encontró durante el resto del año (RR 1 mg: 0,990; IC95% 0,959 – 1,022). Este fenómeno de interacción significativa se observó para otros retardos (1 y 3) del mismo contaminante.

En el caso de las enfermedades respiratorias se observó de forma aislada un fenómeno de interacción en sentido contrario al descrito en el párrafo anterior para los humos negros en el retardo 0. El efecto fue negativo y significativo en el semestre cálido (RR 10 µg/m3: 0,937; IC95% 0,880 – 0,998) y positivo en el más frío (RR 10 µg/m3: 1,021; IC95% 0,990 – 1,054). Esta interacción no se evidenció para el resto de retardos de humos negros ni para los demás contaminantes analizados.

En la figura 3 se representan, para la mortalidad por todas las causas y en mayores de 70 años, los RR correspondientes a incrementos equivalentes (diferencia entre percentiles 90 y 10: P90-10) para diferentes contaminantes. Los incrementos esperados en la mortalidad por todas las causas fueron de un 6% para todos los contaminantes seleccionados salvo el SO2, con un 3,9%. Este incremento se situó en torno a un 7% en la mortalidad de mayores de 70 años.

Figura 3 Asociación de la mortalidad por todas las causas (a) y en mayores de 70 años (b) con la contaminación atmosférica en Valencia, periodo 1994-1996. El efecto representado corresponde al riesgo relativo de mortalidad (e intervalo de confianza del 95%) asociado a un incremento en los niveles de contaminante equivalente a la diferencia entre los percentiles 90 y 10 (P90-10).

Los retardos para cada contaminante (entre paréntesis) son los mismos que los expuestos en la tabla 3.

 CONCLUSIONES

1. Los resultados muestran una asociación significativa de la mortalidad diaria en la ciudad de Valencia, especialmente con la contaminación por partículas, aunque pueden estar implicados otros contaminantes como CO, NO2 y SO2.
2. Las asociaciones entre los niveles de contaminación y el número de defunciones fueron más consistentes para los indicadores de mortalidad total y en mayores de 70 años. Para la mortalidad por enfermedades respiratorias, al igual que para las circulatorias, los hallazgos fueron menos sólidos, no hallándose un efecto significativo para ningún contaminante cuando se analizaba todo el periodo.
3. Los efectos de los diferentes contaminantes son muy similares cuando se calculan para incrementos equivalentes en sus niveles (P90-10), situándose en torno al 6% para la mortalidad por todas las causas y al 7% en mayores de 70 años. Esto podría ser debido, en parte, a la correlación entre contaminantes.
4. Con relación al estudio anterior sobre la mortalidad, realizado para el periodo 1991-19931, la asociación entre humos negros y mortalidad total y en mayores de 70 años es, aunque igualmente significativa, de magnitud ligeramente superior.

BIBLIOGRAFÍA

1. Ballester F, Corella D, Pérez-Hoyos S, Hervás A. Air pollution and mortality in Valencia, Spain: a study using the APHEA methodology. J Epidemiol Community Health 1996; 50: 527-33.

2. A committee of the Environmental and Occupational Health Assembly of the American Thoracic Society. Health effects of outdoor air pollution. Am J Respir Crit Care Med 1996; 153:3-50.

3. Sanz JM. La contaminación atmosférica. Unidades temáticas ambientales de la Secretaría de Estado para las Políticas del Agua y el Medio Ambiente. Madrid: Ministerio de Obras Públicas y Transportes; 1991.

4. Ballester F, Sáez M, Alonso ME, Taracido M, Ordóñez JM, Aguinaga I, et al. El proyecto EMECAM: Estudio multicéntrico sobre la relación entre la contaminación atmosférica y la mortalidad. Rev Esp Salud Pub ().1999; 73:165-175

5. Tello O, Cano R, Hernández G. Situación actual y análisis de la vigilancia epidemiológica en España. En: Vigilancia Epidemiológica: ¿Hacia que modelo vamos?. Madrid: Centro Nacional de Epidemiología; 1992.p. 91-9.

6. World Health Organization. Air quality guidelines for Europe. WHO Regional Publications European Series Nº 23. Copenhagen: WHO; 1987.

7. Sáez M, Pérez-Hoyos S, Tobías A, Saurina C, Barceló MA, Ballester F. Métodos de series temporales en los estudios epidemiológicos sobre contaminación atmosférica. Rev Esp Salud Pub  1999; 73: 133-143.

José María Tenías Burillo (1,2), Santiago Pérez-Hoyos (1), Rosa Molina Quilis (3), Julián González-Aracil (1), Ferran Ballester Díez (1). (1) Institut Valencià d’Estudis en Salut Pública (IVESP), D. G. Salud Pública. Conselleria de Sanidad. Generalitat Valenciana (2) Servei de Medicina Preventiva. Hospital Lluís Alcanyís, Xàtiva (3) Centro de Salud Pública de Alzira. D. G. Salud Pública. Conselleria de Sanidad. Generalitat Valenciana

Correspondencia: Ferrán Ballester Díez. Institut Valencià d’Estudis en Salut Pública. C/Joan de Garay, 21. 46017 Valencia