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
|
|
|
|
|
| Semestre cálido* | Semestre frío* |
| Media | Mínimo | P10 | P90 | Máximo | Media | Media |
Mortalidad por todas | 16,1 | 5 | 11 | 22 | 34 | 14,8 | 17,4 |
Mortalidad en mayores | 11,1 | 3 | 7 | 16 | 24 | 10,0 | 12,3 |
Mortalidad por | 1,5 | 0 | 0 | 3 | 7 | 1,2 | 1,8 |
Mortalidad por | 6,2 | 0 | 3 | 10 | 18 | 5,5 | 7,0 |
Temperatura media (°C) | 18,6 | 6,5 | 11,7 | 26,6 | 33,8 | 22,9 | 14,3 |
Humedad relativa (%) | 63,2 | 23,0 | 43,0 | 79,0 | 96,0 | 65,3 | 61,0 |
Gripe (casos/día) | 172,5 | 3,0 | 13,0 | 398,6 | 660,4 | 68,6 | 277,9 |
Humos negros (µg/m3) | 44,2 | 8,0 | 22,3 | 71,0 | 136,2 | 37,2 | 51,3 |
SO2 24 h | 25,6 | 4,4 | 11,3 | 40,3 | 68,4 | 21,6 | 29,7 |
SO2 1h | 54,9 | 8,8 | 23,2 | 92,0 | 157,4 | 48,2 | 61,7 |
NO2 24 h | 66,8 | 11,9 | 31,8 | 104,5 | 155,1 | 66,8 | 66,9 |
NO2 1h | 116,6 | 31,1 | 58,4 | 176,0 | 304,9 | 116,3 | 116,8 |
CO 24 h (mg/m3) | 2,75 | 0,25 | 1,50 | 4,17 | 6,71 | 2,23 | 3,28 |
CO 1h (mg/m3) | 6,18 | 0,60 | 3,18 | 9,65 | 17,85 | 5,05 | 7,31 |
O3 8h (µg/m3) | 45,5 | 10,2 | 19,5 | 71,6 | 125,0 | 54,2 | 36,6 |
24 h promedio diario; 1h |
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
| Temperatura media | Humedad | Humos negros | SO2 24 h | NO2 1h | CO 24h | O3 8h |
Temperatura media | 1 | 0,07 | -0,41 | -0,42 | -0,10 | -0,55 | 0,45 |
Humedad |
| 1 | 0,07 | -0,18 | -0,10 | -0,02* | -0,10 |
Humos negros |
|
| 1 | 0,63 | 0,33 | 0,64 | -0,57 |
SO2 24 h |
|
|
| 1 | 0,22 | 0,74 | -0,35 |
NO2 1h |
|
|
|
| 1 | 0,03* | -0,26 |
CO 24h |
|
|
|
|
| 1 | -0,42 |
O3 8h |
|
|
|
|
|
| 1 |
Temperatura media: |
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.
| Mortalidad por todas las causas | Mortalidad en mayores de 70 años | Mortalidad por enfermedades respiratorias | Mortalidad por enfermedades circulatorias | ||||
| Ret. | RR (IC 95%)* | Ret. | RR (IC 95%)* | Ret. | RR (IC 95%)* | Ret. | RR (IC 95%)* |
Humos negros | 1 | 1,0127 (1,0029 1,0225) | 1 | 1,0168 (1,0050 1,0288) | 3 | 0,9811 (0,9519 1,0111) | 1 | 1,0095 (0,9949 1,0244) |
SO2 24 h | 3 | 0,9910 (0,9743 1,0083) | 0 | 1,0166 (0,9947 1,0389) | 0 | 1,0309 (0,9785 1,0860) | 2 | 0,9822 (0,9569 1,0081) |
SO2 1h | 0 | 1,0055 (0,9990 1,0121) | 0 | 1,0102 (1,0021 1,0184) | 0 | 1,0118 (0,9916 1,0325) | 4 | 0,9943 (0,9843 1,0044) |
NO2 24 h | 1 | 1,0076 (0,9976 1,0176) | 1 | 1,0116 (0,9997 1,0237) | 4 | 0,9867 (0,9562 1,0181) | 1 | 1,0039 (0,9941 1,0138) |
NO2 1h | 1 | 1,0048 (0,9997 1,0100) | 5 | 1,0068 (1,0007 1,0129) | 1 | 1,0068 (0,9892 1,0247) | 5 | 1,0022 (0,9963 1,0081) |
CO 24 h | 1 | 1,0243 (1,0033 1,0458) | 1 | 1,0243 (0,9988 1,0505) | 1 | 1,0310 (0,9680 1,0982) | 5 | 1,0065 (0,9783 1,0355) |
CO 1h | 3 | 1,0071 (0,9998 1,0145) | 5 | 1,0064 (0,9977 1,0152) | 4 | 0,9940 (0,9720 1,0165) | 5 | 1,0035 (0,9928 1,0143) |
O3 8h | 2 | 1,0118 (0,9985 1,0252) | 2 | 1,0147 (0,9989 1,0307) | 5 | 0,9716 (0,9299 1,0152) | 5 | 1,0127 (0,9929 1,0330) |
* Riesgo relativo de |
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
- 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. - 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. - 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. - 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à dEstudis 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à dEstudis en Salut Pública. C/Joan de Garay,
21. 46017 Valencia
Página anterior | Volver al principio del trabajo | Página siguiente |