EFECTOS A CORTO PLAZO DE LA CONTAMINACIÓN ATMOSFÉRICA SOBRE LA MORTALIDAD: RESULTADOS DEL PROYECTO EMECAM EN DOS CIUDADES DE ASTURIAS
Alvaro Cañada Martínez (1), Jesús Vicente García González (2), Valentín Rodríguez Suárez (1), Federico Fernández Noval (3), Carmen Fernández Rodríguez (4), Ismael Huerta González (1).
(1) Dirección Regional de Salud Pública. Consejería de Servicios Sociales del Principado de Asturias
(2) Gerencia de Atención Primaria. Area Sanitaria V (Gijón). Insalud
(3) Unidad de Salud Pública Area Sanitaria V (Gijón). Consejería de Servicios Sociales del Principado de Asturias
(4) Dirección Regional de Medio Ambiente. Consejería de Fomento del Principado de Asturias
Correspondencia:
Alvaro Cañada Martínez. ]]>
Sección de Análisis e Investigación Epidemiológica. Dirección Regional de Salud Pública.
c/ General Elorza, 32
33001 Oviedo
Tel: 985 10 65 00 (ext 4124) Fax: 985 10 65 20
Correo electrónico: alvarocm@salud.asturias.org
| RESUMEN Fundamento: Los estudios realizados hasta la fecha sobre la posibilidad de que los contaminantes atmosféricos puedan tener efectos sobre la salud no han alcanzado resultados homogéneos. El estudio pretende estimar las medidas de asociación entre los niveles diarios de los contaminantes y la mortalidad a corto plazo en las 2 ciudades más pobladas de Asturias (Gijón y Oviedo), así como contribuir a aumentar el poder estadístico y la representatividad del proyecto EMECAM, en el que se enmarca. Métodos: Estudio ecológico de series temporales. Estimación de niveles de exposición colectivos a partir de los valores recogidos en las estaciones de control de la contaminación. Modelado de las series de mortalidad, incluyendo variables de control, mediante regresión de Poisson. Estimación de riesgos relativos de cada contaminante para la mortalidad, controlando la autocorrelación seriada. Resultados: Durante el período 1993-1996 la contaminación por partículas en suspensión y por CO fué superior en Gijón y la debida a SO2 y NO2 fué mayor en Oviedo. En ambas ciudades los niveles pueden considerarse bajos y dentro de lo admitido por la legislación. La mayoría de las estimaciones de riesgos relativos ha estado próxima al efecto nulo, aunque se han observado relaciones significativas tanto positivas (especialmente para NO2) como negativas. Las asociaciones significativas halladas no se han demostrado consistentes en ambas ciudades para los períodos estudiados. Conclusiones: Con los resultados hallados en este estudio no podemos concluir una asociación clara entre los contaminantes estudiados (partículas, SO2, NO2, CO) y la mortalidad a corto plazo, al menos a los niveles detectados en Gijón y Oviedo. ]]> Palabras clave: Contaminación atmosférica. Mortalidad. Series temporales | ABSTRACT The Short-term Impact of Air Pollution on the Mortality. Results of the EMECAM Project in two cities in Asturias Spain Background: The studies conducted to date regarding the possibility that air pollutants are capable of having impact on human health have not led to homogenous findings. This study is aimed at estimating the degrees of relationship between the daily levels of the pollutants and the death rate on a short-terms basis in the two most populated cities in Asturias (Gijón and Oviedo), as well as contributing to increasing the statistical importance and the representative nature of the EMECAM Project, within which this study is comprised. Methods: Ecological time series study, Estimate of degrees of group exposure based on the readings taken at the pollution control stations. Modeling of the death rate series, including control variables, by means of Poisson regression. Estimating risks related to each pollutant for the death rate, controlling the series-based autocorrelation. Findings: Throughout the 1993-1996 period, the pollution by means of particles in suspension and CO was greater in Gijón, that involving SO2 and NO2 having been greater in Oviedo. In these two cities, the levels can be considered to be low and to fall within what is considered admissible under the laws currently in impact. Most of the relative risk forecasts neared the zero impact point, although significant positive (especially for NO2) as well as negative relationships have been found to exist. The significant relationships found were not proven to be consistent in these two cities for the periods studied. Conclusions: Based on the findings of this study, the conclusion cannot be drawn that a clear-cut relationship exists between the pollutants studied (particles, SO2, NO2, CO) and the death rate on a short-term basis, at least at the levels detected in Gijón and Oviedo. Key words: Air pollution. Mortality. Time series. |
INTRODUCCIÓN
Desde mediados del presente siglo se ha intensificado el estudio de los efectos de la contaminación atmosférica sobre la salud1. En los últimos años se están investigando los efectos a corto plazo sobre la salud de niveles de contaminación moderados2-11.
]]> En Asturias se han elaborado varios informes sobre la situación de la contaminación atmosférica y sobre las relaciones entre calidad del aire y salud12, pero hasta ahora no se había abordado específicamente la posible relación a corto plazo entre contaminación y mortalidad.La existencia en Asturias de una Red de Vigilancia de la Calidad del Aire y de un Registro de Mortalidad permite una aproximación al problema a través de investigaciones de correlación entre ambas variables. El objetivo del presente estudio es investigar la posible relación a corto plazo entre cada uno de los contaminantes y la mortalidad, tomando como unidad de análisis cada día del período estudiado, y realizar estimaciones de las medidas de asociación. Asimismo, se pretende colaborar a aumentar la potencia estadística y la representatividad geográfica del proyecto EMECAM, aportando los datos para el análisis conjunto de las 14 ciudades.
Los requisitos establecidos en el protocolo de estudio en cuanto a población, número de estaciones de control de la contaminación y período de tiempo mínimo, han llevado a seleccionar para este trabajo las ciudades de Oviedo y Gijón entre los años 1993 y 1996. Se valoró inicialmente la posibilidad de incluir otras ciudades asturianas, algunas de ellas con importante desarrollo industrial, pero el número de habitantes y, por tanto, de fallecimientos, era demasiado reducido para los criterios antedichos.
MATERIAL Y MÉTODOS
Se trata de un estudio ecológico de series temporales. Se han analizado las defunciones ocurridas entre el 1 de Enero de 1993 y el 31 de Diciembre de 1996, en Oviedo y Gijón, de personas residentes en estas mismas ciudades. Según el padrón de 1996, Gijón tenía 264.381 habitantes y Oviedo 200.049.
Las variables analizadas, todas ellas recogidas con desagregación diaria, fueron las siguientes: mortalidad por causas, niveles de contaminación atmosférica por dióxido de azufre (SO2), dióxido de nitrógeno (NO2), monóxido de carbono (CO), partículas en suspensión medidas como humos negros (PSUS) y partículas totales en suspensión (PST). Además, se han utilizado como variables de control las siguientes: temperatura media, humedad relativa, incidencia de gripe y variables temporales (tendencia, cambios cíclicos a través de términos sinusoidales, día de la semana, semestre, año, días festivos y otros días "excepcionales"), basándose en el conocimiento o la sospecha de su relación con la variable estudiada (mortalidad) a través de diferentes investigaciones anteriores.
Las fuentes de información han sido las siguientes:
La mortalidad (variable resultado) se ha obtenido del Registro de Mortalidad de la Dirección Regional de Salud Pública del Principado de Asturias. Los datos de contaminación proceden de la Red de Vigilancia de la Contaminación Atmosférica, gestionada por la Dirección Regional de Medio Ambiente del Principado de Asturias, que integra la información de su propia red automática y la que recibe de las redes manuales municipales. La información sobre temperatura y humedad procede del Instituto Nacional de Meteorología (INM). Los datos de incidencia de gripe han sido aportados por el Sistema de Vigilancia Epidemiológica de la Dirección Regional de Salud Pública.
La mortalidad se analizó en 4 categorías (no mutuamente excluyentes): todas las causas, excepto las externas (Clasificación Internacional de Enfermedades, 9ª revisión 13 (CIE-9) 001-799), todas las causas excepto las externas en personas de 70 y más años de edad, causas respiratorias (CIE-9 460-519) y causas cardiovasculares (CIE-9 390-459).
]]> La red manual de medición de la contaminación presenta valores diarios. La red automática presenta los datos integrados cada 15 minutos. Se han calculado los valores horarios como las medias de los 4 valores recogidos cada hora y se han obtenido los valores medios diarios y, para el SO2 y el NO2, los valores horarios máximos de cada día.La estimación de la exposición diaria de la población a cada contaminante se valoró como el promedio de las mediciones de cada una de las estaciones.
La temperatura media diaria ofrecida por el INM se calcula como la media simple entre los valores horarios mínimo y máximo. La humedad relativa como el promedio de 3 determinaciones diarias: a las 7, 13 y 18 horas solares.
Para el período a estudio se disponía inicialmente de información de 3 estaciones automáticas y 13 manuales en Oviedo y 4 estaciones automáticas y 10 manuales en Gijón. Tras la aplicación de los criterios de exclusión establecidos en el protocolo EMECAM, se eliminaron algunas de ellas para determinados contaminantes, por lo que finalmente se trabajó con 3 estaciones automáticas en ambas ciudades, además de las siguientes estaciones manuales: 10 en Gijón y 9 en Oviedo para partículas en suspensión y 9 en Gijón y 7 en Oviedo para SO2.
Las estaciones están distribuídas geográficamente de modo que recogen datos de contaminación de toda el área urbana en ambas ciudades (figura 1).
Figura 1
Ubicación de las estaciones incluídas en el estudio
]]> RESULTADOS
El promedio diario de fallecimientos por todas las causas (excluyendo las externas) entre los residentes en el municipio, que ocurrieron en el propio municipio, fué de 6,3 en Gijón y 4,5 en Oviedo.
Las correlaciones entre las diferentes estaciones captadoras que medían un mismo contaminante y que fueron finalmente incluidas, son siempre mayores de 0,5, excepto en el caso del NO2 en Gijón (tanto en el promedio diario como en el máximo horario), que presenta correlaciones muy bajas. Las correlaciones entre los contaminantes (series promedio definitivas) se presentan en la tabla 1. Se han comprobado también las correlaciones entre las mediciones de SO2 en los casos en los que hay una estación manual y otra automática ubicadas en el mismo lugar o muy próximas (2 en Oviedo y 3 en Gijón), hallándose valores de los coeficientes que oscilan entre 0,34 y 0,54.
Tabla 1
Correlaciones entre contaminantes en Gijón y Oviedo
| GIJÓN | PSUS | ]]> PST | SO2 man | SO2 auto | SO2 máx | NO2 | NO2 máx |
| PST | 0.63 | ||||||
| SO2 (manual) | 0.61 | ]]> 0.43 | |||||
| SO2 (automát) | 0.55 | 0.64 | 0.52 | ||||
| SO2 máximo | 0.35 | 0.50 | 0.38 | 0.86 | |||
| ]]> NO2 | 0.38 | 0.45 | 0.28 | 0.38 | 0.29 | ||
| NO2 máximo | 0.46 | 0.54 | 0.33 | ]]> 0.43 | 0.36 | 0.85 | |
| CO | 0.58 | 0.62 | 0.34 | 0.61 | 0.47 | 0.37 | ]]> 0.51 |
| OVIEDO | PSUS | PST | SO2 man | SO2 auto | SO2 máx | ]]> NO2 | NO2 máx |
| PST | 0.55 | ||||||
| SO2 (manual) | 0.75 | 0.46 | |||||
| SO2 (automát) | 0.69 | 0.68 | ]]> 0.60 | ||||
| SO2 máximo | 0.55 | 0.63 | 0.51 | 0.90 | |||
| NO2 | 0.58 | 0.61 | 0.52 | ]]> 0.72 | 0.63 | ||
| NO2 máximo | 0.44 | 0.51 | 0.42 | 0.56 | 0.52 | 0.84 | |
| CO | ]]> 0.80 | 0.61 | 0.57 | 0.73 | 0.62 | 0.67 | 0.56 |
En la tabla 2 presentamos los estadísticos descriptivos de los niveles de contaminación y de las variables de mortalidad y de control.
Tabla 2 ]]> Estadísticos descriptivos de las variables estudiadas, en unidades/día.
| MEDIA | P10 | MEDIANA | P90 | MEDIA SEMEST FRÍO | MEDIA SEMEST CÁLIDO | |||||||
| GIJÓN | ]]> OVIEDO | GIJÓN | OVIEDO | GIJÓN | OVIEDO | GIJÓN | OVIEDO | GIJÓN | OVIEDO | GIJÓN | ]]> OVIEDO | |
| MORT TODAS NO EXT (n) | 6.3 | 4.5 | 3 | 2 | 6 | 4 | 10 | 7 | ]]> 7.1 | 4.9 | 5.6 | 4.1 |
| MORT TODAS NO EXT 70+ (n) | 4.7 | 3.3 | 2 | 1 | 4 | ]]> 3 | 8 | 6 | 5.3 | 3.6 | 4.1 | 3.0 |
| MORT RESPIRATORIA (n) | 0.7 | 0.5 | ]]> 0 | 0 | 0 | 0 | 2 | 1 | 0.9 | 0.6 | 0.5 | 0.4 |
| ]]> MORT CARDIOCIRCULAT (n) | 2.2 | 1.5 | 0 | 0 | 2 | 1 | 4 | 3 | 2.5 | ]]> 1.7 | 2.0 | 1.4 |
| PSUS (m g/m3) | 52.0 | 28.9 | 25.3 | 11.7 | 45.1 | 20.8 | ]]> 86.0 | 59.5 | 66.9 | 40.3 | 37.8 | 17.6 |
| SO2 manual (m g/m3) | 24.6 | 25.4 | 11.2 | ]]> 11.9 | 19.3 | 21.6 | 44.8 | 44.0 | 31.7 | 32.5 | 17.8 | 17.6 |
| SO2 auto (m g/m3) | ]]> 34.1 | 44.5 | 11.3 | 18.3 | 28.7 | 37.7 | 63.9 | 81.9 | 42.4 | 55.4 | ]]> 26.0 | 33.7 |
| SO2 máximo horario (m g/m3) | 87.8 | 103.5 | 27.1 | 37.1 | 75.0 | 88.3 | 168.6 | ]]> 192.6 | 99.1 | 123.7 | 76.8 | 83.6 |
| PST (m g/m3) | 82.9 | 79.1 | 51.4 | 52.7 | ]]> 75.0 | 73.0 | 127.0 | 116.0 | 92.7 | 83.6 | 73.3 | 74.6 |
| NO2 (m g/m3) | 45.1 | ]]> 50.4 | 28.2 | 34.0 | 43.3 | 49.3 | 62.3 | 67.7 | 49.1 | 56.2 | 41.1 | ]]> 44.7 |
| NO2 máximo horario (m g/m3) | 77.6 | 83.7 | 52.6 | 58.0 | 75.9 | 79.7 | 105.0 | 111.7 | ]]> 84.0 | 92.6 | 71.7 | 75.0 |
| CO (mg/m3) | 1.9 | 1.5 | 1.0 | 0.8 | 1.7 | ]]> 1.3 | 3.1 | 2.6 | 2.2 | 1.9 | 1.6 | 1.2 |
| TEMPERATURA MEDIA (ºC) | 13.8 | 13.2 | ]]> 8.3 | 7.4 | 13.7 | 13.2 | 19.5 | 19.1 | 10.6 | 10.0 | 17.0 | 16.4 |
| ]]> HUMEDAD RELATIVA (%) | 79.4 | 78.0 | 67.0 | 63.0 | 81.0 | 79.0 | 90.0 | 92.0 | 78.6 | ]]> 76.3 | 80.3 | 79.6 |
| GRIPE (n) | 59.4 | 39.1 | 8.1 | 6.1 | 42.4 | 26.7 | ]]> 131.1 | 81.1 | 95.1 | 61.3 | 24.3 | 17.3 |
| * Semestre Frío: Noviembre-Abril. Semestre Cálido: Mayo-Octubre | ||||||||||||
Las estimaciones del riesgo relativo de cada contaminante respecto a cada causa de mortalidad, obtenidas tras el análisis mediante regresión de Poisson, se presentan en las tablas 4 y 5 y en la figura 3. Para la mortalidad total se encuentra una asociación positiva y estadísticamente significativa en el NO2 en Gijón, mientras que en Oviedo, aunque la estimación puntual del riesgo relativo es positiva, el intervalo de confianza engloba la unidad. Se detectan asociaciones negativas estadísticamente significativas para las partículas en Gijón y para el SO2 en Oviedo (aunque la significación estadística es inconsistente cuando comparamos las 2 técnicas de determinación, manual y automática). Por último, para la mortalidad cardiocirculatoria la mayoría de las asociaciones encontradas son positivas y, en el caso de Gijón, estadísticamente significativas, especialmente con el NO2.
]]> Tabla 4| GIJÓN | MORT. TOTAL (excl. externa) | MORT. 70+ (total excl. externa) | MORT. CARDIOCIRCULATORIA | MORT. RESPIRATORIA | ||||||||
| ]]> ret* | RR | ICRR 95% | ret* | RR | ICRR 95% | ret* | RR | ICRR 95% | ret* | ]]> RR | ICRR 95% | |
| PSUS | 4 | 0,9869 | 0,9761 - 0,9979 | 4 | 0,9848 | 0,9722 - 0,9974 | 5 | ]]> 1,0214 | 1,0031 - 1,0401 | 4 | 0,9488 | 0,9197 - 0,9789 |
| SO2 manual | 5 | 1,0120 | 0,9957 - 1,0285 | 5 | ]]> 1,0091 | 0,9902 - 1,0283 | 5 | 1,0227 | 0,9966 - 1,0495 | 4 | 0,9751 | 0,9332 - 1,0188 |
| SO2 auto | 0 | ]]> 0,9946 | 0,9843 -1,0051 | 0 | 0,9896 | 0,9778 - 1,0015 | 5 | 1,0197 | 1,0024 - 1,0372 | 2 | 0,9791 | ]]> 0,9486 - 1,0104 |
| SO2 máx | 5 | 1,0031 | 0,9993 - 1,0070 | 0 | 0,9957 | 0,9912 - 1,0001 | 5 | 1,0072 | ]]> 1,0008 - 1,0136 | 2 | 0,9940 | 0,9821 - 1,0061 |
| PST | 0 | 1,0060 | 0,9985 - 1,0136 | 0 | 1,0062 | ]]> 0,9975 - 1,0150 | 5 | 1,0137 | 1,0010 - 1,0267 | 2 | 0,9784 | 0,9568 - 1,0005 |
| NO2 | 5 | 1,0173 | ]]> 1,0011 - 1,0338 | 5 | 1,0164 | 0,9975 - 1,0358 | 5 | 1,0486 | 1,0220 - 1,0758 | 2 | 0,9437 | 0,8979 - 0,9919 |
| ]]> NO2 máx | 1 | 1,0079 | 0,9976 - 1,0183 | 1 | 1,0047 | 0,9927 - 1,0168 | 4 | 1,0215 | 1,0039 - 1,0395 | ]]> 2 | 0,9747 | 0,9451 - 1,0051 |
| CO | 1 | 0,9819 | 0,9547 - 1,0098 | 1 | 0,9758 | 0,9444 - 1,0083 | ]]> 1 | 0,9666 | 0,9211 - 1,0143 | 5 | 0,9556 | 0,8830 - 1,0341 |
Tabla 5
Estimaciones de la asociación entre contaminantes y mortalidad en Oviedo: riesgo relativo para cada 10 m g/m3 (1 mg/m3 para CO) e intervalos de confianza.
Se especifica el retardo (ret*) de contaminante con asociación más significativa.
| OVIEDO | MORT. TOTAL (excl. externa) | MORT. 70+ (total excl. externa) | MORT. CARDIOCIRCULATORIA | MORT. RESPIRATORIA | ||||||||
| ret* | RR | ICRR 95% | ]]> ret* | RR | ICRR 95% | ret* | RR | ICRR 95% | ret* | RR | ICRR 95% | |
| PSUS | ]]> 5 | 0,9883 | 0,9725 - 1,0044 | 5 | 0,9803 | 0,9609 - 1,0001 | 2 | 1,0153 | 0,9926 - 1,0384 | 1 | ]]> 1,0262 | 0,9864 - 1,0677 |
| SO2 manual | 2 | 0,9711 | 0,9499 - 0,9929 | 2 | 0,9717 | 0,9441 - 1,0001 | 5 | ]]> 0,9862 | 0,9533 - 1,0202 | 1 | 1,0538 | 0,9952 - 1,1159 |
| SO2 auto | 5 | 0,9909 | 0,9795 - 1,0025 | 5 | ]]> 0,9813 | 0,9679 - 0,9949 | 4 | 1,0077 | 0,9900 - 1,0256 | 5 | 0,9866 | 0,9554 - 1,0188 |
| SO2 máx | 2 | ]]> 0,9965 | 0,9920 - 1,0010 | 5 | 0,9946 | 0,9895 - 0,9998 | 2 | 0,9959 | 0,9889 - 1,0030 | 3 | 1,0051 | ]]> 0,9924 - 1,0180 |
| PST | 0 | 0,9937 | 0,9835 - 1,0041 | 5 | 0,9896 | 0,9771 - 1,0022 | 2 | 1,0136 | ]]> 0,9965 - 1,0309 | 2 | 1,0181 | 0,9879 - 1,0493 |
| NO2 | 1 | 1,0128 | 0,9900 - 1,0361 | 5 | 0,9838 | ]]> 0,9574 - 1,0110 | 3 | 1,0335 | 0,9971 - 1,0713 | 5 | 0,9557 | 0,8939 - 1,0218 |
| NO2 máx | 3 | 1,0039 | ]]> 0,9916 - 1,0163 | 5 | 0,9898 | 0,9756 - 1,0042 | 2 | 1,0165 | 0,9976 - 1,0357 | 2 | 1,0302 | 0,9963 - 1,0654 |
| ]]> CO | 5 | 0,9810 | 0,9378 - 1,0262 | 5 | 0,9598 | 0,9094 - 1,0129 | 3 | 1,0513 | 0,9880 - 1,1187 | ]]> 2 | 1,0687 | 0,9572 - 1,1932 |
DISCUSIÓN
Los niveles de contaminación encontrados en Oviedo y Gijón pueden considerarse globalmente bajos y, ocasionalmente, medios. En Gijón, en el período a estudio, los niveles de partículas en suspensión y de CO son mayores que en Oviedo. Por contra Oviedo presenta niveles superiores de NO2 y SO2. Estas diferencias pueden explicarse, en parte, por la ubicación geográfica de ambas ciudades: Gijón en la costa (facilitando el viento la dispersión de contaminantes), Oviedo en el interior y próxima a una cadena montañosa (más afectada por las situaciones de inversión térmica). También influirían los diversos orígenes de la contaminación: en Oviedo fundamentalmente urbana (tráfico y calefacciones), en Gijón mixta (urbana-industrial). Las mediciones automáticas de SO2 dan valores superiores que las manuales, especialmente en Oviedo. Una posible explicación a este hecho estaría en el diferente número y ubicación de las estaciones pero, incluso en aquéllas que están muy próximas, las correlaciones son sólo moderadas y los niveles son mayores en la medición automática. Aunque en caso de discrepancia suelen considerarse más fiables las determinaciones de las estaciones automáticas, los niveles establecidos por la legislación para el SO2 están basados en técnicas manuales.
En el análisis simple de Poisson, incluyendo sólo una variable de contaminación y una de mortalidad, todos los contaminantes presentan asociación positiva significativa (<0,10) para la mortalidad total en ambas ciudades, excepto el SO2 máx en Gijón y las PST en Oviedo.
Los modelos basales definitivos ajustados para la mortalidad total con el objeto de tener en cuenta el efecto de las variables de control, incluyeron, tanto en Gijón como en Oviedo, términos sinusoidales, gripe y temperatura (las variables originales o alguno de sus retardos o funciones cuadráticas), junto con las variables correspondientes a días de la semana, que se incluían en todos los casos siguiendo el protocolo. Además en Gijón se mantuvieron en este modelo final la humedad relativa y la tendencia (tabla 3).
]]> Tabla 3.| MORTALIDAD | GIJÓN | OVIEDO |
| TOTAL (excl. externa) | sen/cos de 1º y 6º orden retardos 3º y 6º de la T2 humedad relativa y hum2 tendencia y tend2 retardo 4º de la gripe | sen/cos de 1º y 5º orden retardo 5º de la temperatura retardo 2º de la T2 retardos 3º y 5º de la gripe |
| ]]> TOTAL 70+ AÑOS (excl. externa) | sen/cos de 1º y 6º orden T2 y su 6º retardo retardo 6º de la humedad relativa y hum2 tendencia y tend2 retardo 4º de la gripe | sen/cos de 1º orden retardo 7º de la temperatura retardo 7º de la T2 años retardos 4º y 5º de la gripe |
| CARDIOCIRCULATORIA | sen/cos de 1º, 5º y 6º orden retardo 6º de la temperatura humedad relativa, hum2 y su 5º retardo años | sen/cos de 4º orden T2 humedad relativa y hum2 |
| RESPIRATORIA | sen/cos de 1º y 5º orden retardo 4º de la humedad relativa retardos 1º y 3º de la hum2 retardo 7º de la gripe | retardos 3º y 4º de la temperatura gripe |
| Notas: sen/cos: pareja de términos sinusoidales; T2: función cuadrática de la temperatura media; hum2 : función cuadrática de la humedad relativa; tend2 : función cuadrática de la tendencia ]]> Siguiendo el protocolo EMECAM se ha incluido siempre la variable día de la semana, como una variable "dummy". | ||
En el modelado de las series de mortalidad se ha conseguido globalmente un buen ajuste, aunque no podemos descartar totalmente la existencia de factores no controlados que puedan afectar a las estimaciones de riesgo.
En la figura 2 podemos comprobar, a modo de ejemplo, la calidad del ajuste del modelo basal para la mortalidad total en Gijón y en Oviedo, mediante las series secuenciales de muertes observadas, muertes predichas por el modelo y residuos. Asimismo, se ha comprobado gráficamente la posible autocorrelación de los residuos y su eliminación tras la inclusión de las series retardadas de mortalidad en el modelo final, a través del uso de funciones de autocorrelación y autocorrelación parcial.
Figura 2
Series de residuos, muertes observadas y esperadas en los modelos basales de mortalidad total.
NOTA: Se ha sumado una constante (10 a los observados y 5 a los esperados) para evitar el solapamiento de las series en el gráfico.
Las estimaciones de riesgos para cada contaminante presentan cierta similitud entre Gijón y Oviedo, como podemos comprobar en la figura 3. Los efectos más significativos no se encuentran siempre con el mismo retardo en ambas ciudades. En Gijón los efectos retardados 4 y 5 días son los que aparecen más frecuentemente como significativos. En Oviedo hay una mayor variabilidad en la máxima significación estadística de los efectos retardados. Debemos tener en cuenta que la probabilidad de asociación espúrea, probablemente, sea mayor cuanto más avanzado sea el retardo.
Figura 3
Estimaciones de riesgo relativo e intervalos de confianza para cada contaminante y cada mortalidad, en ambas ciudades
]]>
NOTA: Están representados los retardos más significativos de cada contaminante
En el análisis estratificado por semestres comprobamos que las estimaciones de riesgo difieren en el semestre frío (más próximas al efecto global) y el semestre cálido. Entre las posibles explicaciones estaría la diferente exposición "efectiva" a los contaminantes en ambos períodos, teniendo en cuenta que el tiempo que se permanece en el exterior es mayor en épocas cálidas que en las frías, aunque no podemos descartar otros factores.
La mayor parte de las estimaciones de riesgo tienen intervalos de confianza que incluyen o están muy próximas al efecto nulo. En esta valoración deben tenerse en cuenta dos consideraciones: el número de habitantes (y por tanto el número de fallecimientos diario), tanto en Gijón como en Oviedo son moderados y los niveles de contaminación en ambas ciudades no pueden calificarse de excesivos. Por estos motivos era de esperar que los efectos detectables, si los hubiera, sean pequeños. A pesar de ello y teniendo en cuenta que prácticamente toda la población está expuesta al factor de riesgo contaminación, incluso la detección de riesgos relativos muy pequeños debe considerarse importante desde el punto de vista de la salud pública. En cualquier caso, la inclusión de estos datos en el análisis conjunto del estudio EMECAM, con un volumen muy importante de población y con representatividad de todo tipo de ciudades españolas, aumentaría la capacidad de detectar efectos significativos si los hubiera.
Se han ensayado algunas alternativas al modelado propuesto en el protocolo, en cuanto a los criterios para introducir y mantener variables de control, con objeto de probar la estabilidad de los resultados. Aunque en la mayoría de los casos no varía el signo de las asociaciones, sí se observan pérdidas de significación estadística, con inclusión del valor "no efecto" en el intervalo de confianza.
Teniendo en cuenta todo lo anterior, especialmente las inconsistencias en el sentido de la asociación y en el orden del retardo más significativo para un mismo contaminante respecto a las diferentes causas de mortalidad y la variabilidad de las estimaciones de riesgo obtenidas, con nuestros datos no podemos concluir la existencia de una clara asociación entre los niveles de contaminantes estudiados y la mortalidad.
Un estudio de series temporales como este tiene, al menos, dos importantes ventajas: por un lado los factores de confusión que no varían de un día a otro (por ejemplo la proporción de fumadores en la población o el nivel socioeconómico de los individuos) no afectarían a las estimaciones; por otra parte, es más importante la estimación de las variaciones diarias de la contaminación que un conocimiento exacto de la exposición.
]]> Como todos los estudios ecológicos, la presente investigación tiene también una serie de limitaciones, la mayoría ligadas a la dificultad en la estimación de la exposición a los factores de riesgo. Por un lado, existe el riesgo de detectar asociaciones atribuibles a un sesgo de agregación o falacia ecológica. Por otra parte, aunque estamos atribuyendo una exposición común a toda la población, las diferencias entre los sensores indican que la exposición no es tan homogénea las correlaciones moderadas indicarían que los incrementos y reducciones en la exposición no son totalmente paralelos entre las distintas estaciones, por lo que la búsqueda de correlaciones con la mortalidad (medida de modo conjunto para todo el municipio) podría verse limitada. Una pequeña parte de la población del municipio (el 5% en Gijón y el 9% en Oviedo) reside fuera del casco urbano y no estaría sometida a los mismos niveles de contaminación. El hospital que atiende un mayor volumen de población en Gijón está fuera del casco urbano, y una proporción importante de fallecimientos se produce tras una estancia en el mismo, muchas veces de una duración mayor que el intervalo máximo entre exposición y efecto (5 días) utilizado en este trabajo. Respecto a la variable resultado, debemos tener en cuenta que la mortalidad es una variable sólo moderadamente sensible para detectar los efectos sobre la salud de la contaminación y que las agrupaciones de causas de mortalidad incluyen causas relativamente heterogéneas, que muy probablemente tengan diferente grado de asociación con los niveles de contaminación.
BIBLIOGRAFÍA
1. Ware JH, Thibodeau LA, Speizer FE, Colome S, Ferris BG. Assessment of the health effects of atmospheric sulphur oxides and particulate matter: evidence from observational studies. Environ Health Perspect 1981; 41: 255-76 [ Links ]
2. Derriennic F, Richardson S, Mollie A, Lellouch J. Short-term effects of sulphur dioxide pollution on mortality in two French cities. International J Epidemiol 1989; 18: 186-97 [ Links ]
3. Schwartz J, Marcus A. Mortality and air pollution in London: a time series analysis. Am J Epidemiol 1990; 131: 185-94 [ Links ]
4. Sunyer J, Antó JM, Murillo C, Sáez M. Effects of urban air pollution on emergency room admissions for chronic obstructive pulmonary disease. Am J Epidemiol 1991; 134: 277-86 [ Links ]
5. Touloumi G, Pocock SJ, Katsouyanni K, Trichopoulos D. Short-term effects of air pollution on daily mortality in Athens: a time series analysis. Int J Epidemiol 1994; 23: 957-67 [ Links ]
6. Momas Y, Pirard P, Quenel P, Medina S, Le Moullec Y, Ferry R et al. Pollution atmosphérique urbaine et mortalité: une synthèse des études épidémiologiques publiées entre 1980 et 1991. Rev Epidémiol Santé Publ 1993; 41: 30-43 [ Links ]
7. Ballester Díaz F, Merino Egea C, Pérez Hoyos S. La asociación entre contaminación atmosférica y mortalidad: una revisión de los estudios epidemiológicos recientes. Rev Esp Salud Pública 1995; 69: 177-88 [ Links ]
8. Katsouyanni K, Zmirou D, Spix C, Sunyer J, Schouten JP, Pönka A et al. Short-term effects of air pollution on health: a European approach using epidemiological time-series data. The APHEA project: background, objectives, design. Eur Respir J 1995; 8: 1030-8 [ Links ]
9. Katsouyanni K, Schwartz J, Spix C, Touloumi G, Zmirou D, Zanobetti A et al. Short term effects of air pollution on health: a European approach using epidemiologic time series data: the APHEA protocol. J Epidemiol Community Health 1996; 50 Supl 1: S12-18 [ Links ]
10. Touloumi G, Katsouyanni K, Zmirou D, Schwartz J, Spix C, Ponce de León A. et al. Short-term effects of ambient oxidant exposure on mortality: a combined analysis within the APHEA project. Am J Epidemiol 1997; 146: 177-85 [ Links ]
11. 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 [ Links ]
12. Grupo de trabajo del Plan de Salud para Asturias. Calidad del Aire y Salud. En: Medio ambiente físico y salud. Documentos monográficos de apoyo al Plan de Salud para Asturias. Oviedo: Consejería de Sanidad y Servicios Sociales. 1993. p. 17-80 [ Links ]
13. Organización Mundial de la Salud. Manual de la Clasificación Estadística Internacional de Enfermedades, Traumatismos y Causas de Defunción. 9ª revisión. Publicación núm 353. Washington : Organización Panamericana de la Salud;1978 [ Links ] ]]>