Contaminación atmosférica por partículas y salud en Valencia,

1994-1996

F. Ballester / C. Iñiguez / S. Pérez-Hoyos / J.M. Tenías
Unitat d'Epidemiologia i Estadística. Escola Valenciana d'Estudis per a la Salut (EVES). Conselleria de Sanitat. 
Generalitat Valenciana.

Correspondencia: Ferran Ballester. Escola Valenciana d'Estudis per a la Salut (EVES). Joan de Garay, 21. 46017 Valencia.
Correo electrónico: ballester_fer@gva.es

Este trabajo se ha beneficiado de las ayudas a la investigación del Fondo de Investigación Sanitaria (FIS: 95/0050-02; 97/0051-01 y 00/0010-1), y de la Conselleria de Sanitat de la Generalitat Valenciana (96/014 y 00/020).

(Particulate air pollution and health in Valencia [Spain] 1994-1996)

Resumen
Objetivo: Estimar la asociación a corto plazo entre los valores de contaminación atmosférica por partículas en suspensión y diferentes indicadores de morbilidad y mortalidad en la ciudad de Valencia.
Métodos: Estudio ecológico de series temporales en el que se relacionan las variaciones diarias en los valores de contaminación atmosférica por humos negros e indicadores de mortalidad, ingresos hospitalarios y visitas a urgencias en la ciudad de Valencia, durante el período 1994-1996. La magnitud de la asociación se estimó mediante regresión autorregresiva de Poisson utilizando modelos aditivos generalizados. Se examinó la forma de la relación, el efecto retardado de la contaminación y la posible modificación de efecto por otros contaminantes o períodos del año.
Resultados: Existe una asociación significativa entre los valores de partículas en suspensión y varios indicadores analizados en Valencia. La forma de dicha relación es lineal. Un incremento de 10 μg/m³ en los valores diarios de humos negros se asoció con un incremento del 1,8% (intervalo de confianza [IC] del 95%, 0,9-2,7%) en el número de defunciones del día siguiente. El mismo incremento en la contaminación se asoció con un aumento del 1,5% (0,1-2,8%) en las defunciones por causa circulatoria, en un 1,3% (0,0-2,6%) en los ingresos por causas circulatorias y en un 5,4 (0,6-10,4%) en el número de visitas a urgencias por asma. Esta asociación no está confundida por los valores de otros contaminantes. No se encontró asociación de los humos negros con la mortalidad por enfermedades respiratorias, los ingresos hospitalarios por enfermedades cerebrovasculares ni las visitas a urgencias por enfermedad pulmonar obstructiva crónica (EPOC).
Conclusiones: Este estudio proporciona evidencia de que los valores de partículas en suspensión en Valencia se asocian con las urgencias por asma, los ingresos por enfermedades circulatorias y las defunciones diarias. La asociación encontrada es consistente con los resultados de otros estudios realizados y es coherente con conocimientos recientes sobre los mecanismos fisiopatológicos.
Palabras clave: Contaminación atmosférica. Partículas en suspensión. Mortalidad. Ingresos hospitalarios. Visitas a urgencias. Enfermedades circulatorias. Asma. EPOC. Series temporales.
]]> Summary
Objective: To estimate the short-term association between levels of air pollution due to suspended particulates and several indicators of morbidity and mortality in the city of Valencia, Spain.
Methods: We performed an ecological time-series study. Daily levels of air pollution from black smoke were related to indicators of mortality, hospital admissions, and visits to emergency departments in the city of Valencia between 1994 and 1996. The magnitude of the association was estimated through Poisson autoregression using generalized additive models. The form of the relationship, the delayed effect of pollution, and the possible modification of this effect by other pollutants or periods of the year were assessed.
Results: A significant association was found between levels of suspended particulates and several of the indicators analyzed in Valencia. The form of this relationship was linear. An increment of 10 μg/m³ in the daily levels of black smoke was associated with an increase of 1.8% (95% confidence interval: 0.9-2.7%) in the number of deaths on the following day. The same increment in pollutant levels was associated with an increase of 1.5% (0.1-2.8%) in deaths from all cardiovascular causes, an increase of 1.3% (0.0-2.6%) in admissions for all cardiovascular diseases and in an increase of 5.4% (0.6-10.4%) in the number of emergency visits for asthma. This association was not confounded by levels of other pollutants. Black smoke was not associated with mortality from respiratory diseases, with hospital admissions for cerebrovascular diseases, or with emergency visits for chronic obstructive pulmonary disease.
Conclusions: This study provides evidence that levels of suspended particles in Valencia are associated with emergency visits for asthma, hospital admissions for cardiovascular diseases, and daily deaths. The association found is consistent with the results of previous studies and fits with recent knowledge of psychopathological mechanisms. 
Key words: Air pollution. Suspended particles. Mortality. Hospital admissions. Emergency visits. Cardiovascular diseases. Asthma. COPD. Time series.

Introducción

En los últimos años una larga serie de trabajos ha puesto de manifiesto la relación entre contaminación atmosférica y diferentes indicadores de salud, incluso con valores de contaminación considerados moderados o bajos. Los estudios llevados a cabo en numerosas ciudades como las incluidas en el proyecto APHEA en Europa^1-4, en el estudio NMMAPS en los EE.UU.⁵ y en otras partes del mundo^6,7 ponen de manifiesto que la contaminación atmosférica sigue representando un riesgo para la salud de la población. En España los resultados de estudios en ciudades como Barcelona^8-10, Madrid^11,12, Valencia^13-15 y el conjunto de las ciudades participantes en el proyecto EMECAM^16-18 indican, de igual modo, la existencia de una relación significativa entre contaminación atmosférica y salud.

La mayoría de los estudios incluyen las partículas en suspensión, y con gran consistencia los resultados muestran una asociación positiva con los indicadores de salud¹⁹. La contaminación por partículas se mide utilizando distintos tipos de indicadores: en Europa, tradicionalmente se ha utilizado el total de partículas en suspensión (TSP) o los humos negros como indicadores de la contaminación por partículas. Este último indicador refleja la concentración de partículas carbonosas en el aire menores de 4,5 micras de diámetro²⁰. En la nueva directiva europea sobre calidad del aire²¹ se establece que el indicador que se debe medir debe ser las partículas con diámetro aerodinámico menor de 10 micras, conocidas como PM₁₀. Por otro lado, los últimos trabajos epidemiológicos y toxicológicos sugieren que son las partículas finas²² (con diámetro inferior a 2,5 micras [PM_2,5]) o las ultrafinas²³ (menos de 0,1 micras de diámetro) las que poseen mayor capacidad tóxica.

Para poder proporcionar estimadores del posible impacto en salud, los estudios epidemiológicos que analizan la relación a corto plazo entre contaminación atmosférica y salud deben tener la calidad y validez suficientes. Para ello son necesarios un diseño adecuado, la obtención de datos fiables (de las variables respuesta, exposición y control) y un método de análisis que tenga en cuenta los factores implicados en la relación contaminación-salud, así como su relación funcional. En este sentido, es importante valorar la relación dosis-respuesta entre el contaminante y el indicador en salud, la estructura de su relación dinámica (es decir, la secuencia temporal de la asociación), así como la posible confusión y modificación de los estimadores obtenidos debidas a otros contaminantes o a variaciones estacionales²⁴.

El objetivo principal de este estudio es estimar la asociación a corto plazo entre los valores de contaminación atmosférica por partículas en suspensión y diferentes indicadores de morbilidad y mortalidad en la ciudad de Valencia, valorando la consistencia y la coherencia de los resultados. Al mismo tiempo se pretende revisar los aspectos metodológicos fundamentales en este tipo de estudios.

Material y métodos

Se trata de un estudio de series temporales en el que se relacionan las variaciones diarias en los valores de contaminación atmosférica y de diferentes indicadores de salud en Valencia para el período 1994-1996.

Datos

En la tabla 1 se exponen los indicadores de salud incluidos en el estudio, las fuentes de información, las poblaciones de referencia y el número diario de casos correspondiente. Para hacer comparables las magnitudes de los diferentes indicadores se presentan también como la tasa de incidencia diaria para una población hipotética de 1.000.000 habitantes. Del Registro de Mortalidad se obtuvo la información del número diario de defunciones en Valencia por todas las causas excepto las externas (Clasificación Internacional de Enfermedades, 9.^a edición [CIE-9]: 1-799), mortalidad por causas del aparato circulatorio (CIE-9: 390-459) y mortalidad por causas respiratorias (CIE-9: 460-519). De la revisión de todas las historias clínicas de urgencias médicas en personas mayores de 14 años del Hospital Clínico Universitario de Valencia durante los años 1994 y 1995 se obtuvo el número diario de casos diagnosticados de asma o de enfermedad pulmonar obstructiva crónica (EPOC). El procedimiento de recogida de datos de las visitas a urgencias y su fiabilidad se han descrito anteriormente²⁶. Por último, el número diario de ingresos hospitalarios por enfermedades del aparato circulatorio (CIE-9: 390-459), de enfermedades cardíacas (CIE-9: 410-414, 427, 428) y de enfermedades cerebrovasculares (CIE-9: 430-438) en los dos principales hospitales de la ciudad (Hospital Clínico Universitario y Hospital Universitario La Fe) se obtuvo de las bases de datos de los registros del Conjunto Mínimo de Bases de Datos hospitalarios (CMBD). Únicamente se incluyeron los ingresos urgentes, a fin de examinar los efectos agudos de la contaminación. Por otro lado, para los tres tipos de indicadores se excluyeron los casos de personas no residentes en la ciudad.

Los valores de contaminación atmosférica se recogieron de la red de vigilancia de la ciudad. Las estaciones de medición seleccionadas para el estudio cumplieron los criterios siguientes: estar orientadas a monitorizar los valores de inmisión de la población, no haber sido sometidas a cambios de emplazamiento y presentar datos válidos en más del 75% de los días de estudio. Se incluyeron seis captadores manuales (para los datos de humos negros) y cinco estaciones automáticas (para el resto de contaminantes) (fig. 1). La correlación entre las diferentes estaciones varió de 0,30 a 0,62 para la red manual (humos negros) y de 0,46 a 0,78 en la red automática (contaminantes gaseosos). En la tabla 2 se exponen los descriptivos correspondientes a las variables de contaminación y a las meteorológicas, así como el método de determinación analítica para las primeras. No se dispuso de información de PM₁₀ durante el período de estudio y, aunque sí se medía el total de partículas en suspensión, los datos disponibles eran menos completos que los de los humos negros.

Figura 1. Mapa de la ciudad de Valencia con indicación de los distritos municipales. 
Ubicación de las estaciones manuales (círculo) y de las estaciones automáticas (triángulo)  ]]> incluidas en el estudio. La estación meteorológica se indica con una estrella. 
Localización del Hospital Clínico Universitario (HC) y del Hospital Universitario La Fe (HF).

En la tabla 3 se puede observar las correlaciones entre estas variables. Los humos negros, el SO₂ y el CO presentan en Valencia un patrón estacional similar con valores más elevados en invierno y más bajos en verano. El ozono presenta un patrón estacional distinto con una correlación positiva con la temperatura. Por último, el NO₂ no presentó un patrón claro, pero sí una correlación positiva, aunque baja, con humos negros y SO₂.

Los datos de las variables meteorológicas se obtuvieron de la estación meteorológica de Valencia, situada en los Jardines de Viveros de la ciudad (fig. 1). Este emplazamiento se prefirió al del aeropuerto ya que refleja de manera más cercana la situación meteorológica en la ciudad. Los casos semanales de gripe se obtuvieron de los registros del sistema de declaración de enfermedades (sistema EDO).

Figura 2. Representación gráfica del ajuste no paramétrico de la relación entre los valores de humos negros y cada 
uno de los indicadores de salud. En cada gráfico se indica el resultado del contraste de la prueba de linealidad correspondiente.

Método

El modelo teórico sobre el que se sustentan las bases para el análisis de las series de indicadores en salud asume algunas hipótesis. La más importante es que las variaciones diarias en el número de defunciones o ingresos hospitalarios pueden explicarse, en parte, por factores que varían en el tiempo (p. ej., fenómenos meteorológicos o incidencia de gripe)²⁷. En consecuencia, estas variables deben ser consideradas para controlar la posible confusión en la estimación del efecto entre los indicadores de salud y los contaminantes atmosféricos. Otra asunción importante es que el número de acontecimientos diarios (defunciones, urgencias, ingresos hospitalarios) sigue un proceso de Poisson no estacionario y con sobredispersión. En consecuencia, la magnitud de la asociación se estima mediante una autorregresión de Poisson. Por último, se entiende que la relación entre los indicadores de salud y las otras variables no tiene por que ser lineal; por ello, para el análisis se utilizan modelos aditivos generalizados (GAM) que permiten la incorporación de funciones de suavizado no paramétricas, de modo que el efecto de las variables se puede controlar sin prefijar la forma de su relación con la variable respuesta²⁸.

El método de análisis empleado se estructuró en dos fases: en la primera se construyó un modelo basal para cada indicador de salud que se pretendía analizar, y en él se incluyeron las principales variables confusoras (estacionalidad y tendencia a largo plazo, periodicidad semanal, temperatura y humedad relativa, incidencia de gripe y acontecimientos no habituales [huelga de médicos en mayo 1995]); en la segunda fase se introdujeron los contaminantes a estudio en el modelo basal. En nuestro caso la variable de humos negros se incorporó primero en forma no paramétrica y, mediante un test de bondad de ajuste que compara las razones de verosimilitudes de los dos modelos, se comprobó su adaptación a la forma lineal. En caso de linealidad se obtuvo el estimador paramétrico de la asociación. Por último se comprobó la correlación seriada hasta el sexto orden y, en caso de existir, se controló mediante la introducción de los términos autorregresivos necesarios. Modelo y procedimiento han sido descritos anteriormente²⁹.

Dado que el efecto de las partículas sobre la salud puede no ser inmediato sino ocurrir algún día después de la exposición³⁰, se realizó una exploración de la relación dinámica evaluando el efecto de las partículas desde el mismo día de la exposición (retardo 0) hasta 5 días después. Por razones de síntesis y comparabilidad, se decidió utilizar promedios de la concentración de contaminantes en lugar de retardos. Así se evita también el riesgo de presentar únicamente los valores más significativos (usualmente positivos) cuando se han llevado a cabo comparaciones múltiples. Los promedios se seleccionaron basándose en los resultados propios y en los seleccionados en otros trabajos de la bibliografía^30,31.

Además de las partículas en suspensión otros contaminantes pueden causar un efecto para la salud²⁰. Aunque cada uno se mide (y se regula) por separado, en la vida real respiramos una mezcla de partículas y gases distinta en cada ambiente. Por ello, cuando se evalúa el efecto de uno de los contaminantes en una ciudad se debe tener en cuenta el resto de contaminantes. En nuestro estudio se valoró tanto la confusión como la modificación de efecto que las variaciones de los contaminantes gaseosos hubieran podido causar sobre los estimadores de las partículas.

En primer lugar, se valoró la confusión mediante la comprobación del cambio producido en los coeficientes de las partículas al introducir, alternativamente, cada uno de los otros cuatro contaminantes. En segundo lugar, se comprobó la posible modificación de efecto. Para ello, se utilizó un modelo de interacción entre el indicador de partículas y las cuatro categorías determinadas por los cuartiles del mismo indicador del contaminante gaseoso.

Varios estudios^31,32 han demostrado que el efecto de las partículas puede ser distinto según la época del año. Con objeto de evaluar esta posible modificación del efecto se utilizó un modelo que incluía la concentración de partículas, la variable calor (variable indicadora que distingue el semestre cálido [mayo a octubre] del frío) y la interacción de ambas, para disponer de la estimación del efecto del contaminante en el período frío. Análogamente, se utilizó la variable frío en la construcción del modelo dirigido a captar el efecto en el período cálido. Los resultados de la asociación se presentan como el incremento porcentual esperado en el número diario del indicador de salud correspondiente por incremento en 10 μ g/m³ en los valores de humos negros.

En la figura 2 se ilustra la relación entre cada uno de los indicadores de salud analizados y los valores de humos negros. Aunque algunas de las curvas parecen sugerir formas de la relación no lineales, los resultados de la prueba de contraste indican que en dicha relación no se puede rechazar la asunción de linealidad.

En la figura 3 se expone la secuencia de la relación dinámica entre los valores de humos negros y los diferentes indicadores de salud. El efecto de las partículas sobre la mortalidad parece ser mayor en los primeros 2 días tras la exposición. Por ello, para la mortalidad se eligió el promedio de los retardos 0 y 1 como indicador de la exposición a humos negros. En el caso de las urgencias por asma la secuencia del efecto presenta un patrón similar, por lo que se seleccionó el mismo promedio. Para la EPOC no se observa un efecto, por lo que en orden a la comparabilidad se eligió también el promedio 01. Por último, la relación con las enfermedades del aparato circulatorio presenta una secuencia algo más retardada, por lo que se seleccionó el promedio de los retardos 0, 1 y 2 como indicador de la exposición a humos negros.

Figura 3. Relación dinámica en el tiempo entre contaminación atmosférica por humos negros y los indicadores de salud, 
desde el día simultáneo (retardo 0) al quinto día después de la exposición. Todos los resultados se exponen 
como el incremento porcentual (y su intervalo de confianza del 95%) en el número de casos correspondiente 
a cada indicador de salud (defunciones, visitas a urgencias, ingresos hospitalarios) asociado con un 
incremento de 10 μ g/m ³ en los valores diarios de cada uno de los retardos de humos negros.

]]> En la tabla 4 se resume el resultado de la asociación entre los indicadores de salud y las concentraciones de partículas. Se encontró una asociación significativa con la mortalidad total, la mortalidad por enfermedades del aparato circulatorio y las urgencias por asma. La asociación es cercana a la significación estadística en el caso de los ingresos hospitalarios por enfermedades circulatorias y enfermedades cardíacas. No se encontró asociación con la mortalidad por enfermedades respiratorias, los ingresos hospitalarios por enfermedades cerebrovasculares ni las visitas a urgencias por EPOC.

En la figura 4 se presenta el resultado del análisis para controlar la confusión de los otros contaminantes. En general, los coeficientes de la asociación con humos negros no se vieron alterados sustancialmente al incluir cada uno de los otros contaminantes en los modelos correspondientes. Es de destacar, sin embargo, el cambio en la magnitud de la asociación de humos negros con mortalidad por todas las causas y por enfermedades del aparato circulatorio al incluir la variable con las concentraciones de SO₂ en el modelo. También llama la atención la variación de los estimadores de la asociación con ingresos por enfermedades circulatorias. Por un lado, el efecto de la inclusión del SO₂ sobre los coeficientes es distinto del observado con mortalidad, es decir, en este caso la magnitud de la asociación disminuye. Por otro lado, el control por NO₂ y CO conlleva el aumento de la magnitud del estimador central y se alcanza la significación estadística de la asociación.

Figura 4. Efecto de la inclusión de los contaminantes gaseosos en la estimación del riesgo de los humos negros sobre cada uno 
de los indicadores de salud estudiados. El riesgo está expresado como el incremento porcentual (y su intervalo de 
confianza del 95%) en el número de casos correspondiente a cada indicador de salud asociado con un incremento 
de 10 μ g/m ³ en los valores de humos negros del promedio de los retardos 0 y 1 (para mortalidad e ingresos por 
asma y por EPOC) y en los valores del promedio de los retardos 0, 1 y 2 para los ingresos hospitalarios.

]]> Las estimaciones de la asociación de las partículas demostraron ser robustas al análisis, considerando los valores de los otros contaminantes. En la figura 5 se ilustra la posible modificación de efecto tanto para defunciones como para ingresos hospitalarios por enfermedades del aparato circulatorio, indicadores que en el análisis anterior habían demostrado ser más sensibles a la inclusión de un segundo contaminante. Aunque no se encontró interacción estadísticamente significativa se observa cierta tendencia a obtener estimadores de mayor magnitud a cocentraciones más elevadas de NO₂ y de CO.

Figura 5. Análisis estratificado de la asociación de humos negros con la mortalidad (a) y los ingresos hospitalarios (b), ambas por causas del aparato circulatorio, según cuartiles de los contaminantes gaseosos. El riesgo está expresado como el incremento porcentual (y su intervalo de confianza del 95%) en el número de casos correspondiente a cada indicador de salud asociado 
a un incremento de 10 μ g/m ³ en los valores de humos negros del promedio de los retardos 0 y 1 (para mortalidad) 
y en los valores del promedio de los retardos 0, 1 y 2 para los ingresos hospitalarios.
 Los contaminantes gaseosos se introdujeron en los mismos retardos que los humos negros.

Por último, en la figura 6 se observan los estimadores para los dos semestres analizados (frío y cálido) junto con el estimador correspondiente a todo el período de estudio. Los coeficientes de la asociación con mortalidad son ligeramente mayores para los meses fríos, en especial para la mortalidad respiratoria y urgencias por asma y EPOC. Por contra, existe un efecto diferente para los ingresos por enfermedades del aparato circulatorio, y la asociación es mayor en los meses más calurosos.

Figura 6. Asociación entre contaminación atmosférica por humos negros y los indicadores de salud, para todo el período 
y para los semestres frío (noviembre a abril) y cálido (mayo a octubre). El riesgo está expresado como el incremento  ]]> porcentual (y su intervalo de confianza del 95%) en el número de casos correspondiente a cada indicador de salud 
asociado con un incremento de 10 μ g/m ³ en los valores de humos negros del promedio de los retardos 0 y 1 (para mortalidad 
e ingresos por asma y por EPOC) y en los valores del promedio de los retardos 0, 1 y 2 para los ingresos hospitalarios.

Discusión

Los resultados de este estudio demuestran la relación a corto plazo de la contaminación atmosférica por partículas y distintos indicadores de mortalidad y morbilidad en la ciudad de Valencia durante el período 1994-1996.

El método de análisis utilizado en este estudio introduce algunas variaciones, como el uso de modelos aditivos generalizados, sobre análisis anteriores en Valencia^13-15. Los resultados obtenidos con los dos métodos proporcionan estimadores del efecto consistentes, aunque con el método utilizado aquí se obtienen resultados más precisos²⁹; además permite examinar y controlar formas de relación no lineales entre variables. La asociación lineal de humos negros con los indicadores de salud encontrada en nuestro estudio es consistente con lo encontrado en estudios realizados en los EE.UU.^24,33, e incluso en España³⁴. Ello implica que no se puede asumir un valor umbral por debajo del cual no exista riesgo de efectos adversos para la población.

Los resultados del análisis de los modelos de dos contaminantes sugieren que las asociaciones encontradas no fueron debidas a confusión por otros contaminantes. Sin embargo, debido a la elevada correlación entre humos negros, SO₂ y CO, los estimadores pueden verse afectados por problemas de multicolinealidad. Por ello creemos recomendable no introducir más de dos contaminantes a la vez en los modelos y comparar, posteriormente, los cambios en los coeficientes del contaminante sujeto a estudio.

Por otro lado, aunque la correlación con el NO₂ es más baja, los resultados indican que con valores más elevados de NO₂ (y también de CO) el efecto de las partículas es mayor. Este hallazgo concuerda con lo encontrado en el estudio APHEA2³⁵ y podría indicar que cuando los valores de NO₂ y CO son más elevados las partículas en suspensión proceden en mayor medida de los combustibles del tráfico de vehículos.

Para la mortalidad por causas respiratorias y las urgencias por asma, sin embargo, se ha encontrado cierta diferencia entre los estimadores por estaciones, con valores mayores (aunque no significativos) para los meses de invierno. Un posible argumento a favor de esta diferencia de efecto sería que la exposición a partículas podría incrementar la susceptibilidad de la población de padecer infecciones respiratorias³⁸. En un análisis realizado en Filadelfia se encontró que las defunciones por causa cardiovascular en los días de contaminación más alta presentaban infecciones respiratorias como causas contribuyentes, mientras que en los días de contaminación más baja esto no ocurría³⁹.

La magnitud de las asociaciones encontradas varía según los indicadores estudiados, y es parecida para las defunciones y los ingresos hospitalarios, y sensiblemente mayor para las urgencias por asma. La asociación no fue significativa en la mortalidad por causas respiratorias, las urgencias por EPOC ni los ingresos por enfermedades cerebrovasculares. Sin embargo, en otros estudios realizados en Valencia se ha encontrado relación significativa entre las urgencias por EPOC con otros contaminantes, como el O₃ y el CO⁴⁰.

En general, los resultados obtenidos en Valencia se encuentran en el espectro de los obtenidos en otros estudios (tabla 5). A la hora de comparar nuestros resultados con los de la bibliografía hay que considerar tres aspectos importantes: por un lado, la mayoría de los estudios provienen de Norteamérica, con características ambientales, sociales y de sistemas de información sanitaria distintas de las nuestras; por otro lado, el indicador de contaminación por partículas rara vez se corresponde con el que nosotros hemos dispuesto (humos negros) excepto en el estudio APHEA³⁵; por último, los resultados de la tabla 5 provienen de la combinación de resultados individuales de un número considerable de estudios, por lo que su precisión será mayor que la obtenida en Valencia.

Para la mortalidad los estimadores del efecto en Valencia son algo más elevados que los obtenidos en los EE.UU. y se asemejan más a los obtenidos en el estudio italiano³¹. Tanto en Italia como en el proyecto APHEA se ha encontrado cierto patrón norte-sur con un efecto de la contaminación mayor en las ciudades con temperatura media más alta. En concreto, las ciudades incluidas en el proyecto APHEA que se situaron por encima del percentil 75 de la distribución de la temperatura media presentaron un incremento del riesgo del 0,8 y el 0,7%, para PM₁₀ y humos negros, respectivamente, frente al 0,6% global³⁵. En el estudio APHEA no se encontró asociación significativa de partículas (PM₁₀) con la mortalidad respiratoria, al igual que sucedió en Valencia con los humos negros. En cuanto a los indicadores de morbilidad, los resultados para enfermedades cardiovasculares en Valencia concuerdan con los obtenidos en los otros estudios⁴¹. Por otro lado, en el estudio APHEA2 no se encontró una asociación significativa entre humos negros e ingresos por enfermedades respiratorias⁴², similar a lo encontrado en Valencia.

A la hora de interpretar los resultados obtenidos hay que tener en cuenta que el estudio presenta algunas limitaciones: el corto período analizado (2 años para las series de urgencias), junto con el número escaso de algunos indicadores de salud, resulta en estimadores poco precisos; por otro lado, al tratarse de riesgos pequeños el error aleatorio podría tener una influencia relativa importante en los estimadores obtenidos. Por ello es recomendable llevar a cabo estudios multicéntricos que permitan reducir este problema y explorar los factores que pudieran explicar la heterogeneidad en los resultados.

Las limitaciones del uso de los datos sobre ingresos en estudios epidemiológicos se han planteado anteriormente, tanto de manera general⁴³ como específica, para los estudios de contaminación atmosférica⁴⁴. El análisis de la calidad de las bases de datos del CMBD en la Comunidad Valenciana ha evidenciado pocos errores en la gestión de los datos, excepto en la residencia⁴³. Este problema podría producir un sesgo de mala clasificación, que sería no diferencial, ya que no parece lógico que los errores en los datos sobre residencia del paciente (u otros campos) se relacionen con las variaciones diarias de partículas en suspensión. Por otro lado, un estudio llevado a cabo en Canadá⁴⁵ puso de manifiesto un buen grado de concordancia (kappa: 0,84) en la clasificación de las enfermedades cardíacas.

Otro motivo de crítica para el uso de los datos de visitas y/o ingresos hospitalarios es el posible efecto de los reingresos en los estimadores obtenidos. En otro reciente estudio en Canadá⁴⁴ se ha investigado este tema construyendo dos tipos de series de ingresos: una para pacientes con un solo ingreso durante el período de estudio y otra con los pacientes con múltiples ingresos. Los estimadores de efecto no difirieron o fueron más pequeños en el caso de pacientes con más de un ingreso. Esta evaluación sugeriría que el efecto de incluir las visitas repetidas sería, en caso de ocurrir, hacia la nulidad. Por otro lado, se debería tener en cuenta que el número de ingresos hospitalarios en un día concreto depende de la disponibilidad del sistema sanitario, y en situaciones epidémicas, por ejemplo, puede haber una disponibilidad limitada de camas hospitalarias. En este sentido, se ha argumentado que dado que el análisis efectuado se basa en las variaciones diarias esto no debería producir graves distorsiones³¹. En cualquier caso, la puesta en marcha reciente de este sistema de información (CMBD), sus objetivos orientados a la gestión sanitaria y las posibles fuentes de error en los datos recomiendan precaución a la hora de interpretar los resultados en estudios epidemiológicos en los que se hayan utilizado.

Un último aspecto que cabe considerar se relaciona con el diseño utilizado. Una crítica común a los estudios de series temporales es su naturaleza ecológica. Sin embargo, cuando el interés se centra en valorar la relación entre las variaciones en el tiempo y los distintos contaminantes, en salud a escala poblacional, el diseño de series temporales tiene validez, además de ser eficiente y de utilidad para la actuación en salud pública. La limitación fundamental de este diseño es que las estimaciones que se obtienen únicamente nos permiten medir los acontecimientos desencadenados por la contaminación a corto plazo. Por ello, para valorar de forma adecuada el impacto, a corto y largo plazo, de la contaminación sobre la salud (p. ej., años de vida perdidos) se necesitan estudios de cohortes.

En los últimos años se han propuesto varios mecanismos fisiopatológicos⁵¹. Una de las principales hipótesis es que las partículas inducen una activación de algunos mediadores que provocan un incremento de la coagulabilidad sanguínea^52,53. Otros mecanismos estudiados están relacionados con el control autonómico cardíaco en los que se ha demostrado una asociación significativa de las partículas con un aumento de la frecuencia cardíaca y con la disminución de su variabilidad^54-56. De forma reciente también se ha encontrado una asociación de las partículas con la hipertensión arterial⁵⁷. Por último, en estudios experimentales en humanos se han encontrado cambios inflamatorios y un incremento de valores de fibrinógeno, leucocitos y plaquetas en la zona alveolar asociados con partículas^58,59.

Dos cuestiones relacionadas con la medida de la exposición, tal como se ha utilizado en este trabajo, merecen mención. Una de ellas hace referencia a una limitación general de los estudios que utilizan estaciones medidoras del ambiente exterior. Debido a las diferencias individuales entre lo medido por las estaciones y por monitores personales, se ha puesto en duda la asociación encontrada en estudios de series temporales. Sin embargo, cuando se han comparado las variaciones diarias de las partículas entre los dos tipos de medidas, se ha encontrado una correlación alta⁶¹. Por otro lado, el examen estadístico de las implicaciones del error de medida de la exposición en este tipo de estudios sugiere que el posible error derivaría en un sesgo que, en todo caso, subestimaría la asociación encontrada⁶¹.

La otra cuestión se refiere al contaminante utilizado en este estudio: los humos negros. A pesar de haber sido sustituido por las PM₁₀ en los estándares de la Unión Europea²¹, hay estudios que indican que son un buen marcador de los efectos de las partículas sobre la salud^62,63. Además, los humos negros han demostrado ser un mejor marcador de la exposición a partículas derivadas del tráfico que las PM₁₀ e incluso que las PM_2,5⁶⁴. Lo anterior, junto con la existencia de series históricas de medidas fiables de humos negros, nos lleva a recomendar, desde el punto de vista de la epidemiología y la salud pública, que se mantengan las mediciones de humos negros. Junto con estos dos indicadores sería de gran interés que se fuera extendiendo el número de estaciones que determinen la concentración de partículas finas (PM_2,5) para las cuales ya existe una norma europea contrastada⁶⁵. Por último, dada la mezcla compleja que representa la contaminación atmosférica urbana, con fuentes comunes y alta correlación entre contaminantes, es difícil atribuir los efectos a un único contaminante. Más bien las partículas deben considerarse un marcador de los valores de contaminación en una ciudad.

En conclusión, este estudio proporciona evidencia de que los valores de partículas en suspensión en Valencia se asocian de forma estadísticamente significativa con las urgencias por asma, los ingresos por enfermedades circulatorias y las defunciones diarias. Esta asociación no está confundida por los valores de otros contaminantes. La asociación encontrada es consistente con los resultados de otros estudios realizados y es coherente con conocimientos recientes sobre los mecanismos fisiopatológicos. Por otro lado, la forma de la relación lineal indica que las medidas encaminadas a reducir los valores medios de contaminación atmosférica serían más eficientes que aquellas encaminadas a evitar unos pocos días con valores altos³⁴. A pesar de que los riesgos observados son pequeños, las implicaciones para la salud pública son importantes debido a que la exposición afecta a toda la población⁶⁶. Entre 1985 y 1998 la contaminación atmosférica por humos negros en la ciudad de Valencia se ha reducido prácticamente en un 50%⁶⁷, aunque como indican los resultados de este estudio aún representa un riesgo para la salud de la población. Además, otros contaminantes como el NO₂, el ozono o las partículas más finas podrían presentar una situación más desfavorable. Todo ello hace que sigan siendo necesarias estrategias y acciones encaminadas a que mejore la calidad ambiental y, por tanto, la salud de los ciudadanos.

Agradecimientos

Queremos expresar nuestro agradecimiento a los profesionales del Laboratorio de Medio Ambiente y del Servicio de Salud Pública del Ayuntamiento de Valencia; a los compañeros de los proyectos EMECAM-EMECAS, APHEA y APHEIS; a Julián González, Toni Merelles, Rosa Molina, Loreto Suay, Fernando Gómez, María Paz Rodríguez, Marisa Rivera, Voro Peiró, Julián Librero, Vicente Escoms y a todas las personas que nos han brindado su ayuda.

2. Zmirou D, Schwartz J, Saéz M, Zanobetti A, Wojtyniak B, Touloumi G, et al. Time-series analysis of air pollution and cause-specific mortality. Epidemiology 1998;9:495-503.        [ Links ]

3. Sunyer J, Spix C, Quenel P, Ponce A, Barumandzadeh T, Touloumi G, et al. Urban air pollution and emergency admissions for asthma in four European cities: the APHEA Project. Thorax 1997;52:760-5.        [ Links ]

4. Anderson HR, Spix C, Medina S, Schouten JP, Castellsague J, Rossi G, et al. Air pollution and daily admissions for chronic obstructive pulmonary disease in 6 European cities: results from the APHEA project. Eur Respir J 1997;10:1064-71.        [ Links ]

5. Samet JM, Zeger SL, Dominici F, Curriero F, Coursac I, Dockery DW, et al. The National Morbidity, Mortality, and Air Pollution Study (NMMAPS). Part 2. Morbidity and mortality from air pollution in the United States. 2000. Health Effects Institute.        [ Links ]

6. Romieu I, Meneses F, Ruiz S, Sierra Monge JJ, Huerta J, White MC, et al. Effects of air pollution on the respiratory health of asthmatic children living in Mexico City. Am J Respir Crit Care Med 1996;154:300-7.        [ Links ]

7. Xu X, Gao J, Dockery DW, Chen Y. Air pollution and daily mortality in residential areas of Beijing, China. Arch Environ Health 1994;49:216-22.        [ Links ]

8. Sunyer J, Sáez M, Murillo C, Castellsague J, Martínez F, Antó JM. Air pollution and emergency room admissions for chronic obstructive pulmonary disease: a 5-year study. Am J Epidemiol 1993;137:701-5.        [ Links ]

9. Castellsague J, Sunyer J, Sáez M, Antó JM. Short-term association between air pollution and emergency room visits for asthma in Barcelona. Thorax 1995;50:1051-6.        [ Links ]

10. Tobías A, Sunyer J, Castellsague J, Sáez M, Antó JM. Impacto de la contaminación atmosférica sobre la mortalidad y las urgencias por enfermedad pulmonar obstructiva crónica y asma en Barcelona. Gac Sanit 1998;12:219-26.        [ Links ]

11. Díaz J, Alberdi JC, Montero JC, Mirón IJ. Asociación entre la contaminación atmosférica por dióxido de azufre y partículas totales en suspensión y la mortalidad en la ciudad de Madrid (1986-1992). Gac Sanit 1998;12:203-11.        [ Links ]

12. Díaz J, García R, Ribera P, Alberdi JC. Modeling of air pollution and its relationship with mortality and morbidity in Madrid, Spain. Int Arch Occup Environ Health 1999;72:366-76.        [ Links ]

13. 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 ]

14. Tenías JM, Ballester F, Rivera ML. Asssociation between hospital emergency visits for asthma and air pollution in Valencia, Spain. Occup Environ Med 1998;55:541-7.        [ Links ]

15. Ballester F, Tenías JM, Pérez-Hoyos S. Air pollution and emergency hospital admissions for cardiovascular diseases in Valencia, Spain. J Epidemiol Community Health 2001;55:57-65.        [ Links ]

16. Grupo EMECAM. El proyecto EMECAM: estudio multicéntrico español sobre la relación entre la contaminación atmosférica y la mortalidad. Rev Esp Salud Pública 1999;73:105-314.        [ Links ]

17. Ballester F, Saez M, Perez-Hoyos S, Iñíguez C, Gandarillas A, Tobías, et al. The EMECAM project: a multi-center study on air pollution and mortality in Spain. Combined results for particulates and for sulphur dioxide. Occup Environ Med 2002;59:300-8.        [ Links ]

18. Saez M, Ballester F, Barceló MA, Pérez-Hoyos S, Bellido J, Tenías JM, et al. A combined analysis of the short-term effects of photochemical air pollutants on mortality within the EMECAM project. Environ Health Perspect 2002;110:221-8.        [ Links ]

19. Pope CA, Dockery DW. Epidemiology of particles effects. In: Holgate ST, Samet JM, Koren H, Maynard RL, editors. Air pollution and health. San Diego: Academic Press, 1999; p. 673-705.        [ Links ]

20. A Committee of the Environmental and Occupational Health Assembly of the American Thoracic Society. Health effects of outdoor air pollution (Part 1). Am J Respir Crit Care Med 1996;153:3-50.        [ Links ]

21. Council Directive 1999/30/EC of 22 April 1999 relating to limite values for sulphur dioxide, nitrogen dioxide and oxides of nitrogen, particulate matter and lead in ambient air. Official Journal of the European Communities 1999;L163.        [ Links ]

22. Pope CA. Epidemiology of fine particulate air pollution and human health: biologic mechanisms and who's at risk? Environ Health Perspect 2000;108(Suppl 4):713-23.        [ Links ]

23. Donaldson K, Stone V, Clouter A, Renwick L, MacNee W. Ultrafine particles. Occup Environ Med 2001.        [ Links ]

24. Schwartz J. Assessing confounding, effect modification, and thresholds in the association between ambient particles and daily deaths. Environ Health Perspect 2000;108:563-8.        [ Links ]

25. Bates DV. Health indices of the adverse effects of air pollution the question of coherence. Environ Res 1992;59: 336-49.        [ Links ]

26. Merelles A, Rivera ML, Ballester F, Antó JM, Oterino D. Identificación de urgencias por asma y EPOC en el servicio de urgencias de un hospital. Estudio de concordancia. Gac Sanit 1999;13:34-41.        [ Links ]

27. Schwartz J, Spix C, Touloumi G, Bacharova L, Barumandzadeh T, Le Tertre A, et al. Methodological issues in studies of air pollution and daily counts of deaths or hospital admissions. J Epidemiol Community Health 1996;50(Suppl 1): S3-11.        [ Links ]

28. Hastie TJ, Tisbhirani RJ, Hastie TJ, Tisbhirani RJ, editors. Generalized additive models. London: Chapman and Hall, 1990.        [ Links ]

29. Iñíguez C, Pérez-Hoyos S, Ballester F, Saez M, Gómez F, Tenías JM. Comparación de dos metodologías estadísticas para el análisis de los efectos a corto plazo de la contaminación en la salud [resumen]. Gac Sanit 2001;15(Supl 2):31.        [ Links ]

30. Schwartz J. The distributed lag between air pollution and daily deaths. Epidemiology 2000;11:320-6.        [ Links ]

31. Biggeri A, Bellini P, Terracini B, editors. Meta-analysis of the Italian studies on short-term effects of air pollution. Epidemiologia Prevenzione 2001;25.        [ Links ]

32. Michelozzi P, Forastiere F, Fusco D, Perucci CA, Ostro B, Ancona C, et al. Air pollution and daily mortality in Rome, Italy. Occup Environ Med 1998;55:605-10.        [ Links ]

33. Daniels MJ, Dominici F, Samet JM, Zeger S. Estimating particulate matter-mortality dose-response curves and threshold levels: an analysis of daily time-teries for the 20 Largest US cities. Am J Epidemiol 2000;152:397-406.        [ Links ]

34. Schwartz J, Ballester F, Saez M, Bellido J, Cambra K, Arribas F, et al. The concentration-response relation between air pollution and daily deaths. Environ Health Perspect 2001;109: 1001-6.        [ Links ]

35. Katsouyanni K, Touloumi G, Samoli E, Gryparis A, Le Tertre A, Monopolis Y, et al. Confounding and effect modification in the short-term effects of ambient particles on total mortality: results from 29 European cities within the APHEA2 Project. Epidemiology 2001;12:521-31.        [ Links ]

36. Pekkanen J, Brunner EJ, Anderson HR, Tiittanen P, Atkinson RW. Daily concentrations of air pollution and plasma fibrinogen in London. Occup Environ Med 2000;57:818-22.        [ Links ]

37. Katsouyanni K. Health effects of air pollution in southern Europe: are there interacting factors? Environ Health Perspect 1995;103(Suppl 2):23-7.        [ Links ]

38. Schwartz J. Is there harvesting in the association of airborne particles with daily deaths and hospital admissions? Epidemiology 2001;12:55-61.        [ Links ]

39. Schwartz J. What are people dying of on high air pollution days? Environmental Research 1994;64:26-35.        [ Links ]

40. Tenías JM, Ballester F, Pérez-Hoyos S, Rivera ML. Hospital emergency room admissions for chronic obstructive pulmonary disease and air pollution in Valencia, Spain. Arch Environ Health 2002;57:41-7.        [ Links ]

41. Le Tertre A, Medina S, Samoli E, Forsberg B, Michelozzi P, Boumghar A, et al. Short term effects of particulate air pollution on cardiovascular disease in eight European cities: a quantitative summary. J Epidemiol Community Health 2002;56: 773-9.        [ Links ]

42. Atkinson RW, Ross AH, Sunyer J, Ayres J, Baccini M, Vonk JM, et al. Acute effects of particulate air pollution on respiratory admissions. Results from APHEA 2 project. Am J Respir Crit Care Med 2001;164:1860-6.        [ Links ]

43. Librero J, Ordiñana R, Peiró S. Análisis automatizado de la calidad del conjunto mínimo de datos básicos. Implicaciones para los sistemas de ajuste de riesgos. Gac Sanit 1998;12: 9-21.        [ Links ]

44. Lippmann M, Ito K, Nadas A, Burnnet RT. Association of particulate matter components with daily mortality and morbidity in urban populations. Health Effects Institute, 2000.        [ Links ]

45. Stieb DM, Beveridge RC, Rowe BH, Walter S, Judek S. Assessing diagnostic classification in an emergency departament: implications for daily time series studies of air pollution. Am J Epidemiol 1998;148:666-70.        [ Links ]

46. Pope CA, Bates DV, Raizenne ME. Health effects of particulate air pollution: time for reassessment? Environ Health Perspect 1995;103:472-80.        [ Links ]

47. Bates DV. Lines that connect: assessing the causality inference in the case of particulate pollution. Environ Health Perspect 2000;108:91-2.        [ Links ]

48. Moolgavkar SH, Luebeck EG. A critical review of the evidence on particulate air pollution and mortality. Epidemiology 1996;7:420-8.        [ Links ]

49. Gamble JF, Lewis RJ. Health and respirable particulate (PM10) air pollution: a causal or statistical association? Environ Health Perspect 1996;104:838-50.        [ Links ]

50. Pope CA, Dockery DW, Kanner RE, Villegas GM, Schwartz J. Oxigen saturation, pulse rate, and particulate air pollution. Am J Respir Crit Care Med 1999;159:365-72.        [ Links ]

51. Seaton A, MacNee W, Donaldson K, Godden D. Particulate air pollution and acute health effects. Lancet 1995;345: 176-8.        [ Links ]

52. Peters A, Döring A, Wichmann HE, Koenig W. Increased plasma viscosity during an air pollution episode: a link to mortality? Lancet 1997;349:1582-7.        [ Links ]

53. Schwartz J. Air pollution and blood markers of cardiovascular risk. Environ Health Perspect 2001;109(Suppl 3):405-9.        [ Links ]

54. Pope CA, Verrier RL, Lovett EG, Larson AC, Raizenne ME, Kanner RE, et al. Heart rate variability associated with particulate air pollution. Am Heart J 1999;138:890-9.        [ Links ]

55. Liao D, Creason J, Shy C, Williams R, Watts R, Zweidinger R. Daily variation of particulate air pollution and poor cardiac autonomic control in the elderly. Environ Health Perspect 1999;107:521-5.        [ Links ]

56. Peters A, Perz S, Doring A, Stieber J, Koening W, Wichmann HE. Increases in heart rate during an air pollution episode. Am J Epidemiol 1999;150:1094-8.        [ Links ]

57. Ibald-Mulli A, Stieber J, Wichmann HE, Koenig W, Peters A. Effects of air pollution on blood pressure: a population-based approach. Am J Public Health 2001;91:571-7.        [ Links ]

58. Salvi S, Blomberg A, Rudell B, Kelly F, Sandstrom T, Holgate ST, et al. Acute inflammatory responses in the airways and peripheral blood after short-term exposure to diesel exhaust in healthy human volunteers. Am J Respir Crit Care Med 1999;159:702-9.        [ Links ]

59. Ghio AJ, Kim C, Devlin RB. Concentrated ambient air particles induce mild pulmonary inflammation in healthy human volunteers. Am J Respir Crit Care Med 2000;162:981-8.        [ Links ]

60. Janssen NA, Hoek G, Brunekreef B, Harssema H, Mensink I, Zuidhof A. Personal sampling of particles in adults: relation among personal, indoor, and outdoor air concentrations. Am J Epidemiol 1998;147:537-47.        [ Links ]

61. Zeger S, Thomas D, Dominici F, Samet JM, Schwartz J, Dockery DW, et al. Exposure measurement error in time-series studies of air pollution: concepts and consequences. Environ Health Perspect 2000;108:419-26.        [ Links ]

62. Brunekreef B, Janssen NA, De Hartog J, Harssema H, Knape M, Van Vliet P. Air pollution from truck traffic and lung function in children living near motorways. Epidemiology 1997;8: 298-303.        [ Links ]

63. Schwartz J, Marcus A. Mortality and air pollution in London: a time series analysis. Am J Epidemiol 1990;131: 185-94.        [ Links ]

64. Roemer W, Van Wijnen JH. Differences among black smoke, PM10, and PM1.0 levels at urban measurement sites. Environ Health Perspect 2001;109:151-4.        [ Links ]

65. Ramos MC, Rubio L, Martín D, Hernández P, García dos Santos R, Fernández R. Desarrollo de la norma europea (EN) para la determinación de las partículas PM2,5 (respirables) en aire ambiente. Método de captación por un filtro y posterior análisis gravimétrico [abstract]. Rev Salud Ambient 2001;1:102.        [ Links ]

66. Künzli N, Kaiser J, Medina S, Studnicka M, Chanel O, Filliger P, et al. Public Health impact of outdoor and traffic-related air pollution: a european assessment. Lancet 2000;356: 795-801.        [ Links ]

67. Fernández-Patier R, Ballester F. Objetivo 21: mejorar la calidad del aire. En: Álvarez-Dardet C, Peiró S, editores. Informe SESPAS 2000. Barcelona: Doyma, 2000; p. 279-89.        [ Links ] ]]> 1 1997 314 1658-63 2 1998 9 495-503 3 1997 52 760-5 4 1997 10 1064-71 5 2000 6 1996 154 300-7 7 1994 49 216-22 8 1993 137 701-5 9 1995 50 1051-6 10 1998 12 219-26 11 1998 12 203-11 12 1999 72 366-76 13 1996 50 527-33 14 1998 55 541-7 15 2001 55 57-65 16 Grupo EMECAM 1999 73 105-314 17 2002 59 300-8 18 2002 110 221-8 19 1999 673-705 20 A Committee of the Environmental and Occupational Health Assembly of the American Thoracic Society 1996 153 3-50 21 1999 L163 L163 22 2000 108 ^s4 ^s4 4 713-23 23 2001 24 2000 108 563-8 25 1992 59 336-49 26 1999 13 34-41 27 1996 50 ^s1 ^s1 1 S3-11 28 1990 29 2001 15 ^s2 ^s2 2 31 30 2000 11 320-6 31 2001 25 32 1998 55 605-10 33 2000 152 397-406 34 2001 109 1001-6 35 2001 12 521-31 36 2000 57 818-22 37 1995 103 ^s2 ^s2 2 23-7 38 2001 12 55-61 39 1994 64 26-35 40 2002 57 41-7 41 2002 56 773-9 42 2001 164 1860-6 43 1998 12 9-21 44 2000 45 1998 148 666-70 46 1995 103 472-80 47 2000 108 91-2 48 1996 7 420-8 49 1996 104 838-50 50 1999 159 365-72 51 1995 345 176-8 52 1997 349 1582-7 53 2001 109 ^s3 ^s3 3 405-9 54 1999 138 890-9 55 1999 107 521-5 56 1999 150 1094-8 57 2001 91 571-7 58 1999 159 702-9 59 2000 162 981-8 60 1998 147 537-47 61 2000 108 419-26 62 1997 8 298-303 63 1990 131 185-94 64 2001 109 151-4 65 2001 1 102 66 2000 356 795-801 67 2000 279-89