Tesis 699-150818 PDF
Tesis 699-150818 PDF
Tesis 699-150818 PDF
TESIS DOCTORAL:
Dirigida por:
I
Índice
Índice
Agradecimientos.......................................................................................... I
Resumen......................................................................................................... IX
Abstract.......................................................................................................... X
1. Introducción general............................................................................. 1
1.1. Aerosoles atmosféricos.......................................................................... 3
1.2. Estado del arte....................................................................................... 5
1.3. Motivación............................................................................................. 6
1.4. Objetivos de la tesis doctoral................................................................ 7
1.5. Estaciones de medida............................................................................ 8
1.5.1. Tamanrasset (Argelia).................................................................... 9
1.5.2. Tenerife (Islas Canarias, España)................................................... 10
1.5.2.1. Observatorio Atmosférico de Izaña (IZO)......................... 12
1.5.2.2. Observatorio de Santa Cruz de Tenerife (SCO)................ 14
1.6. Estructura de la memoria..................................................................... 15
2. Metodología.............................................................................................. 17
2.1. Gestión de la calidad de la instrumentación radiométrica solar...... 19
2.2. Medidas de aerosoles de la red AERONET........................................ 20
2.2.1. La red AERONET.......................................................................... 20
2.2.2. Instrumento estándar: Fotómetro solar Cimel 318......................... 21
2.2.3. Protocolo de calibración................................................................. 23
2.2.3.1. Calibración solar directa................................................... 24
2.2.3.2. Calibración de radiancia................................................... 25
2.2.4. Procesado de datos......................................................................... 26
2.2.4.1. Parámetros calculados a partir de las medidas de
radiación solar directa..................................................... 26
2.2.4.2. Parámetros calculados a partir de las medidas de
radiancia.......................................................................... 27
2.2.5. Distribución de los datos................................................................ 30
III
Índice
IV
Índice
V
Índice
VI
Índice
VII
Resumen
IX
Abstract
X
Capítulo 1
Introducción general
Introducción general
3
Capítulo 1
4
Introducción general
5
Capítulo 1
1.3. Motivación
La región subtropical del Hemisferio Norte, a pesar de ser una región muy amplia,
está poco estudiada desde el punto de vista atmosférico debido a que está ocupada en su
mayor parte por zonas desérticas y océanos. Sin embargo, como se ha mencionado
anteriormente, el desierto del Sáhara y la región subtropical del Atlántico nororiental
son zonas clave para el estudio del polvo mineral, su transporte y sus efectos a escala
global.
Conscientes de esta importancia, el programa de aerosoles en columna del Centro de
Investigación Atmosférica de Izaña (CIAI, http://izana.aemet.es/) perteneciente a la
Agencia Estatal de Meteorología (AEMET, http://www.aemet.es/) ha centrado sus
esfuerzos en la implantación y gestión de un fotómetro solar en Tamanrasset (Argelia),
en medio del desierto del Sáhara, destinado al estudio a largo plazo del polvo mineral
6
Introducción general
7
Capítulo 1
Figura 1.2. Localización de Tamanrasset (Argelia), en el corazón del desierto del Sáhara, y de Tenerife
(Islas Canarias, España), en la región subtropical el Atlántico nororiental.
8
Introducción general
9
Capítulo 1
Figura 1.3. Fotómetro solar Cimel instalado en la azotea del edificio principal del Centro Meteorológico
Regional de la ONM en Tamanrasset.
Tamanrasset es una ciudad libre de actividades industriales que está cerca de varias
fuentes de polvo mineral localizadas en el sur de Argelia, Mali, Libia y Chad. Su clima
está modulado por la influencia de las masas de aire húmedas del monzón en verano,
mientras que el invierno se caracteriza por ser seco, con ocasionales cirros y nubes
medias (Cuesta et al., 2008; Guirado et al., 2014).
10
Introducción general
Figura 1.4. a) Mapa topográfico de Tenerife (la equidistancia entre las curvas de nivel es de 100 metros
de altitud). b) Contorno vertical de Tenerife entre las estaciones de Izaña (IZO) y de Santa Cruz de
Tenerife (SCO) junto con la estratificación de la baja troposfera en la SENAR, que suele estar bien
establecida en la época estival. Fuentes: a) Rodríguez y Guerra, 2001 y b) Alastuey et al., 2005.
11
Capítulo 1
Figura 1.5. a) Vista general de la montaña de Izaña con el volcán Teide al fondo, b) edificio principal del
Observatorio Atmosférico de Izaña (IZO) y c) fotómetros solares Cimel CE-318 instalados en la azotea de
la torre del edificio principal.
IZO forma parte del sistema de VAG de la OMM cuya finalidad es la de conocer la
evolución a largo plazo de la composición química y características físicas de la
atmósfera con el fin de mejorar el conocimiento sobre los procesos atmosféricos y su
relación con las modificaciones del clima a diferentes escalas espacio-temporales. En
este contexto se enmarca su participación en la red VAG de Radiómetros de Filtros de
12
Introducción general
13
Capítulo 1
Figura 1.6. a) Vista general del Observatorio de Santa Cruz de Tenerife (SCO), cerca del puerto de la
ciudad, b) instalaciones del Observatorio en la azotea del edificio y c) fotómetro solar Cimel CE-318
instalado en SCO.
14
Introducción general
15
Capítulo 1
16
Capítulo 2
Metodología
19
Capítulo 2
20
Metodología
Figura 2.1. Distribución global de las estaciones que proporcionan observaciones de aerosoles en
columna en el marco de la red AERONET. Imagen tomada de la página web de AERONET
(http://aeronet.gsfc.nasa.gov) el 14 de marzo de 2015.
Para posibilitar la gestión de una red tan amplia se ha impuesto una estandarización
de aspectos tales como la instrumentación, el protocolo de calibración y el procesado de
datos.
21
Capítulo 2
Las medidas solares directas se realizan en todas las longitudes de onda disponibles
en el instrumento. Para cada una de estas longitudes de onda se determina el espesor
óptico de aerosoles (AOD, Aerosol Optical Depth), salvo para el canal de 940 nm que
se utiliza para la obtención del contenido en vapor de agua de la columna atmosférica
(PWV, Precipitable Water Vapor). Siguiendo la secuencia de medidas programada por
defecto en la caja de la electrónica, el instrumento comienza a medir por la mañana para
una masa óptica de 7 y termina para la misma masa óptica por la tarde. A partir de la
primera medida el resto de ellas se realizan cada 0.25 intervalos de masa óptica hasta la
masa óptica 2 a partir de la cual se realizan cada 15 minutos, repitiéndose la misma
secuencia de intervalos de medida para las masas ópticas de la tarde. Cada observación
consiste, en realidad, en tres medidas tomadas aproximadamente cada 30 segundos,
creando un triplete de medidas para cada longitud de onda, que facilita la detección de
medidas que deben ser desechadas, ya sea por fallo del instrumento o por la presencia
de nubes durante la medida.
En lo que respecta a las observaciones de radiancia, hay dos tipos de medidas (plano
principal y almucántar) y se realizan para cuatro longitudes de onda: 440, 670, 870 y
1020 nm. Como paso previo a las medidas de radiancia se realiza siempre una medida
solar directa que permite fijar la posición del sol (ángulos cenital y acimutal). El primer
tipo de medida de radiancia se realiza hasta nueve veces al día y consiste en seguir el
plano principal solar, es decir, fijado el ángulo acimutal solar se varía el ángulo cenital
desde 6º por debajo del sol hasta 140º en sentido opuesto (Figura 2.3a y Tabla 2.1). El
segundo tipo se realiza hasta seis veces al día y consiste en seguir el almucántar solar, es
decir, fijado el ángulo cenital solar se varía el ángulo azimutal 180º hacia la derecha
primero y, tras apuntar de nuevo al sol, 180º hacia la izquierda después (Figura 2.3b y
Tabla 2.1). El objetivo es obtener medidas de radiancia de la aureola solar y del cielo en
22
Metodología
Figura 2.3. Descripción visual de las medidas de (a) plano principal (PP) y (b) almucánar (ALM).
Imágenes adaptadas del manual del usuario del fotómetro solar CE-318
(http://support.cimel.fr/photo/pdf/man_ce318_us.pdf).
Tabla 2.1. Secuencia angular de las medidas de almucántar y plano principal. Los ángulos del almucántar
corresponden a posiciones acimutales respecto del Sol (0º) y se repiten a la derecha y a la izquierda del
mismo. Los ángulos del plano principal corresponden a posiciones cenitales respecto del Sol (0º) con
valores negativos indicando las medidas por debajo del mismo.
Tras la realización tanto de las medidas solares directas como de radiancia, el robot
está programado para que la cabeza sensora adopte una posición de reposo, muy
próxima al nadir. Además, el fotómetro está conectado a un sensor de humedad que al
activarse fuerza al fotómetro a adoptar la posición de reposo. De esta manera se evita
que se depositen partículas extrañas o gotas de agua en la óptica del fotómetro.
23
Capítulo 2
24
Metodología
I λ = I 0λ ·e −τ λ ·m (2.1)
25
Capítulo 2
del GSFC que es transferida a otros centros de calibración tales como el CIAI, el GOA
y el LOA. La transferencia de la calibración desde el GSFC se realiza mediante un
fotómetro maestro viajero que actualiza la calibración de las esferas integrantes
secundarias aproximadamente cada cuatro meses.
Cabe mencionar las deficiencias de este método de transferencia de la calibración.
Por un lado, la precisión del instrumento que transfiere la calibración es del mismo
orden que la del instrumento que pretende ser calibrado. Por otro lado, la longitud de
onda nominal no es la misma para todos los fotómetros, con el consiguiente error
asociado en la transferencia de la calibración. Por todo ello, la precisión de la
calibración en radiancia se estima en torno a un 5% (Holben et al, 1998).
Los parámetros derivados a partir de las medidas de radiación solar directa y que
han sido analizados en el presente trabajo se describen a continuación:
• El AOD se calcula utilizando la ley de Beer-Bouguer-Lambert (ecuación 2.1) a
partir de las medidas de extinción solar en cada longitud de onda. La atenuación
debida a la dispersión Rayleigh (τR) y la absorción del ozono (τO) y de otros gases
atmosféricos (τg) son estimadas para eliminar su contribución al espesor óptico total:
AOD ≡ τa = τ - τR - τO - τg. Tal y como se mencionó anteriormente la incertidumbre
asociada a las medidas de AOD es, aproximadamente, 0.01 excepto en el rango
ultravioleta (UV) que asciende a 0.02.
• El PWV se determina (en cm) a partir del canal de 940 nm. La contribución de la
absorción de los aerosoles y de la dispersión Rayleigh se estima mediante
extrapolación a partir de las medidas de los canales 670 y 870 nm que es sustraída
de las medidas en 940 nm. La incertidumbre asociada al cálculo del PWV se estima
en ±10% (Holben et al., 2001).
26
Metodología
τ a (λ ) = β ·λ−α (2.3)
donde τa(λ) es el AOD para una longitud de onda dada, β se conoce como el
coeficiente de turbiedad (corresponde al AOD para λ = 1µm) y α es el EA.
El EA es un indicador cualitativo del tamaño predominante de los aerosoles y puede
ser calculado a partir de diferentes longitudes de onda (Schuster et al., 2006). Para
estudios climatológicos se suele usar el calculado en el rango 440-870 nm, a partir
de las longitudes de onda en las que haya medidas de AOD dentro de ese rango.
Puede tomar valores entre 0 y 4, aunque para los aerosoles lo habitual es que varíe
entre 0 y 2.5. Los valores más altos (mayores que 1) indican el dominio del aerosol
fino mientras que los valores más bajos (menores que 1) indican el dominio del
aerosol grueso. Los valores próximos a 4 están relacionados con la dispersión
molecular (Rayleigh). A partir de la incertidumbre del AOD (0.01-0.02) se estima la
del EA entre 0.03 y 0.04 (Schuster et al., 2006).
• La fracción del modo fino (FMF, Fine Mode Fraction) se define:
FMF=AODfine/AOD (2.4)
donde AODfine es el espesor óptico del modo fino.
La FMF se computa, en 500 nm, a partir del AOD y del EA mediante el algoritmo
de deconvolución espectral (SDA, Spectral Deconvolution Algorithm) desarrollado
por O’Neill et al. (2001; 2003). Este algoritmo infiere las propiedades ópticas de los
modos fino y grueso de la distribución de partículas asumiendo que la derivada
espectral del EA del modo grueso es casi nula, y estableciendo una función
aproximada que relaciona el EA del modo fino y su derivada (O’Neill et al., 2001;
2003). La incertidumbre asociada a los parámetros derivados del SDA varía con la
inversa del AOD total, teniendo en cuenta no sólo los errores de las medidas sino
también la incertidumbre de los valores asumidos para la curvatura espectral de cada
modo (O’Neill et al., 2001).
27
Capítulo 2
• El radio efectivo (Ref) de las partículas (en µm) se obtiene a partir de la distribución
del número de partículas tanto para la distribución total como para la del modo fino
y la del modo grueso, a partir de la expresión:
R max
dN ( R )
∫ R3
d ln R
d ln R
Ref = R min
R max
(2.7)
dN ( R )
∫
R min
R2
d ln R
d ln R
Propiedades ópticas
• El índice de refracción complejo depende de la longitud de onda incidente λ. Se
calcula para cada una de las cuatro longitudes de onda en las que se realizan las
medidas de radiancia y da cuenta de la dispersión y de la absorción a través de su
parte real e imaginaria:
m(λ)=n(λ)-ik(λ) (2.8)
La parte imaginaria es un indicador de la absorción de las partículas con valores
menores cuanto menos absorbentes sean las partículas, siendo nula para partículas
no absorbentes.
• El albedo de dispersión simple (SSA, Single Scattering Albedo) es un indicador de
la capacidad de absorción de las partículas y se define como la fracción de luz
dispersada con respecto al total de luz extinguida al interactuar con las partículas:
SSA=σdis/σext (2.9)
siendo σdis y σext los coeficientes de dispersión y extinción, respectivamente, que se
definen:
28
Metodología
R2
σ dis (λ ) = π ∫ R 2 Qdis (m, λ , R) ⋅ N ( R) dR (2.10)
R1
R2
σ ext (λ ) = π ∫ R 2 Qext (m, λ , R ) ⋅ N ( R) dR (2.11)
R1
Figura 2.4. Parámetro de asimetría, g(λ), como indicador de la redistribución angular de la radiación
incidente. Se muestran ejemplos de casos en los que domina (a) la dispersión de la luz hacia delante, (b)
la dispersión isótropa y (c) la dispersión hacia atrás.
29
Capítulo 2
a) Es positivo cuando las partículas son del orden o mayores que la longitud de
onda incidente (teoría de dispersión de Mie), dominando la dispersión de la luz
hacia delante (Θ<90º, forward-scattering).
b) Es nulo en el caso en el que las partículas son mucho más pequeñas que la
longitud de onda incidente (teoría de dispersión de Rayleigh), dando lugar a una
dispersión isótropa.
c) Es negativo cuando las partículas tienen un índice de refracción muy alto,
dominando la retrodispersión, es decir, la dispersión de la luz hacia atrás
(90º<Θ<180º, backward-scattering).
Hay que puntualizar que el cálculo del SSA y del índice de refracción complejo en
el nivel 2.0 de AERONET se realiza solo para medidas de AOD(440 nm)>0.4. El
motivo es que la precisión de estos dos parámetros, en condiciones de baja carga de
aerosoles, decrece significativamente: hasta un 80-100% para la parte real del índice de
refracción, un 0.05-0.07 para la parte imaginaria del índice de refracción y un 0.05 para
el SSA (Dubovik et al., 2000, 2002).
En lo que respecta a la incertidumbre de las medidas, Dubovik et al. (2000) realizan
un detallado análisis de sensibilidad de los numerosos factores que pueden influir en el
error asociado a los productos del algoritmo de inversión. Dicho estudio analiza desde
la cobertura de diferentes ángulos de dispersión hasta errores aleatorios y sistemáticos
(p. ej., de calibración, de apuntamiento, del modelo de transferencia radiativa).
30
Metodología
combinado del EA, entre 440 y 870 nm, y de su curvatura espectral, representada por la
diferencia del exponente de Ångström δAE=EA(440-670 nm)-EA(670-870 nm). Estos
parámetros se relacionan, a su vez (Figura 2.5), con la FMF (en %) y el tamaño del
modo fino del aerosol (en µm) mediante puntos de referencia correspondientes a la
distribución bimodal de tamaños de partículas esféricas, que han sido determinados
mediante la teoría de Mie, en base al índice de refracción típico del aerosol
urbano/industrial (m=1.4-0.001i). No es de esperar que la suposición de partículas
esféricas afecte significativamente en los resultados (Gobbi et al., 2007). En lo que
respecta a la sensibilidad al índice de refracción del método gráfico (Figura 2.5), el nivel
de indeterminación es del orden del ±25% para el radio del modo fino del aerosol (Rf) y
del ±10% para la FMF, cuando se varía el índice de refracción entre m=1.33–0.000i
(típico de las gotas de agua) y m=1.53–0.003i (típico del polvo mineral). Para evitar
errores mayores que 30% en EA y en δAE, este método se aplica solo a valores de AOD
(a 675 nm) mayores que 0.15 Gobbi et al. (2007). Basart et al. (2009) aplicaron este
método gráfico para rastrear y caracterizar mezclas de polución y polvo mineral,
confirmando la solidez de esta metodología.
31
Capítulo 2
Tabla 2.2. Características principales del Lidar Micropulsado versión 3 (MPL-3) en operación en la
estación de Santa Cruz de Tenerife
Características MPL-3
Longitud de onda (nm) 523
Tipo de láser Nd:YLF
Frecuencia de repetición de los pulsos (Hz) 2500
Energía máxima por pulso (µJ) 7
Duración de cada pulso (ns) 10
Tiempo de integración
60
de la señal retrodispersada (segundos)
Resolución vertical (m) 75
Láser inofensivo para la vista (eye-safe) Sí
Sistema receptor Telescopio Matsukov-Cassegrain
Fabricante SES Inc., EE. UU.
32
Metodología
(Fernald et al., 1972; Klett, 1981, 1985; Fernald, 1984; Sasano and Nakane, 1984) que
determina los coeficientes de extinción a partir de la señal del Lidar combinada con
datos AERONET de AOD de la estación SCO. Además, hay que tener en cuenta que las
medidas de niveles inferiores a los 250 m sobre el nivel del suelo son descartadas
debido a limitaciones intrínsecas a la configuración del instrumento (overlap).
Figura 2.6. LIDAR instalado en Santa Cruz de Tenerife y vista del láser durante la noche.
2.4.2. CALIOP
El instrumento CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) es un
Lidar de retrodispersión elástica instalado a bordo del satélite CALIPSO (Cloud-
Aerosol Lidar and Infrared Pathfinder Satellite Observation) que fue lanzado en abril
de 2006 por la NASA en colaboración con el CNES (Centre National d'Etudes
Spatiales, Francia; www.cnes.fr). CALIOP emite luz linealmente polarizada en las
longitudes de onda a 532 y 1064 nm para proporcionar observaciones, con resolución
vertical y a escala global, de los aerosoles y las nubes (Hunt et al., 2009; Winker et al.,
2009). En este trabajo se han utilizado los perfiles de extinción CALIOP, de nivel 2 y
versión 3.01, en la longitud de onda de 532 nm calculados sobre Tamanrasset (en un
radio de 1.5º) con una resolución vertical de 60 metros (por debajo de los 20.2 km de
altura) y con una resolución horizontal de 5x5 km2. Los datos, que están disponibles en
la base de datos de la NASA (https://eosweb.larc.nasa.gov/cgi-
bin/searchTool.cgi?Dataset=CAL_IIR_L1-Prov-V1-10), han sido filtrados siguiendo la
metodología descrita por Tesche et al. (2013) con el fin de seleccionar sólo los perfiles
de máxima calidad.
33
Capítulo 2
Figura 2.7. Recreación artística del funcionamiento del sensor CALIOP a bordo del satélite CALIPSO.
Fuente: http://www.nasa.gov/mission_pages/calipso/mission/mission-objectives.html.
34
Metodología
La capacidad del MISR para obtener medidas en diferentes longitudes de onda y con
diferentes ángulos de observación le permite inferir propiedades de los aerosoles sobre
las superficies desérticas, altamente reflectantes (Kahn el al., 2005, 2010). Es por este
motivo por el que se ha elegido para determinar la distribución espacial de las posibles
fuentes de polvo mineral en el norte de África. Para ello se han usado los promedios
diarios de AOD (nivel 3) disponibles en el canal verde (555 nm) con una resolución
espacial de 0.5º×0.5º (aproximadamente 50x50 km2).
Figura 2.8. Recreación artística del funcionamiento del sensor MISR a bordo del satélite Terra,
representando los ángulos de visión de cada una de sus nueve cámaras que son sensibles en cuatro bandas
espectrales: azul, verde, rojo e infrarrojo cercano. Fuente: https://www-misr.jpl.nasa.gov/Mission/
(Imagen cortesía de Shigeru Suzuki y Eric M. De Jong, NASA JPL).
35
Capítulo 2
−1
N N
C ij = ∑ nijk ⋅ ∑ C k nijk (2.13)
k =1 k =1
donde nij es el número de horas que las trayectorias pasan en la celda (i,j).
36
Metodología
37
Capítulo 2
Figura 2.9. Recreación artística del satélite geoestacionario MSG que lleva a bordo el sensor SEVIRI.
Fuente: https://directory.eoportal.org/web/eoportal/satellite-missions/m/meteosat-second-generation.
38
Metodología
Figura 2.10. Recreación artística del funcionamiento del sensor MODIS a bordo del satélite Terra.
Fuente: http://www.ntsg.umt.edu/sites/ntsg.umt.edu/files/imce/EOS_AM1_scan.jpg.
39
Capítulo 3
43
Capítulo 3
Instrumento Características
Cabeza sensora: Ø80 mm
Fotómetro solar CE-318
Longitud cabeza-colimadores: 400 mm
Filtros (y anchura a altura mitad) en nm:
UV: 340 (2), 380 (2)
Visible: 440 (10), 500 (10), 675 (10), 870 (10)
Medida de vapor de agua: 940 (10)
IR: 1020 (10), 1640 (25)
Especrofotómetro Brewer MkIII
Dimensiones carcasa impermeable
espectrofotómetro: 700x460x340 mm.
Dimensiones tracker:
Chasis: 305x305x356 mm
Altura sin el trípode: 457 mm
Peso: 34 Kg
Piranómetro
Yankee Environmental System
(YES) UVB-1
Piránometro
Kipp&Zonen Series CMA y CMP
44
Laboratorio de calibración óptica del CIAI
45
Capítulo 3
46
Laboratorio de calibración óptica del CIAI
Figura 3.1. Plano del laboratorio de calibración óptica del Observatorio Atmosférico de Izaña en el que
se indican las zonas del mismo mencionadas en el texto.
47
Capítulo 3
48
Laboratorio de calibración óptica del CIAI
Figura 3.2. Esquema de diseño del sistema de calibración de irradiancia absoluta en posición horizontal.
Siguiendo este diseño, los componentes básicos del sistema han de ser los
siguientes:
• Un banco óptico sobre el que se asiente todo el sistema facilitando su alineamiento.
• Un soporte horizontal regulable, que se desplace sobre el banco óptico, para la
sujeción de los instrumentos a ser calibrados, teniendo en cuenta las diferencias
(dimensión y peso) de los mismos (piranómetros y radiómetros descritos en la Tabla
3.1).
• Un soporte, que se desplace sobre el banco óptico, para la sujeción de bafles
ignífugos, con diferentes aperturas, que minimicen la luz parásita que pueda llegar al
detector del instrumento a calibrar.
• Un soporte, que se desplace sobre el banco óptico, para la sujeción de lámparas
halógenas tipo FEL (código asignado por ANSI, American National Standards
Institute, http://www.ansi.org/) de 1000 W.
• Un soporte, que se desplace sobre el banco óptico, para la sujeción de un láser que
garantice el alineado del sistema.
El siguiente paso consiste en un análisis pormenorizado de las características
específicas que debe reunir el sistema para conseguir que el proceso de calibración sea
lo más sencillo posible pero ajustándose a todas las necesidades que requiere dicho
proceso:
• Las lámparas halógenas tipo FEL de 1000 W han de ser trazables al estándar de
irradiancia primario del Physikalisch-Technische Bundesanstalt (PTB,
http://www.ptb.de/index_en.html) y han de estar estacionadas (es decir, previo a su
uso en laboratorio, han debido estar en operación a su máxima intensidad durante el
49
Capítulo 3
50
Laboratorio de calibración óptica del CIAI
Figura 3.3. Esquema de diseño del sistema de calibración de irradiancia absoluta en posición vertical.
De acuerdo con este diseño, los componentes básicos del sistema han de ser los
siguientes:
• Un soporte vertical al que deben ir acoplados tres brazos horizontales de 40 cm cada
uno:
Un primer brazo horizontal para la sujeción del láser-guía de alineación vertical
(eje Z).
Un segundo brazo horizontal para la sujeción de lámparas halógenas tipo DXW
(código ANSI) de 1000 W.
Un tercer brazo de sujeción de bafles ignífugos, con diferentes aperturas, para
minimizar la luz parásita que pueda llegar al detector del instrumento a calibrar.
• Un láser-guía de alineación horizontal (eje X) que debe estar acoplado en el soporte
vertical del sistema.
• Bases regulables en altura para el alojamiento, soporte y nivelación de los
instrumentos a calibrar.
Una vez establecido el diseño básico del sistema se procede al análisis
pormenorizado de las características específicas que debe verificar para garantizar el
correcto proceso de calibración.
• Las lámparas halógenas tipo DXW de 1000 W han de estar estacionadas y ser
trazables al estándar de irradiancia primario del PTB.
• Debe disponer de un sistema de control y adquisición de la intensidad de la corriente
que circula por la lámpara para garantizar que, durante todo el proceso de
calibración, permanece en el rango establecido por el certificado de calibración de la
misma.
51
Capítulo 3
52
Laboratorio de calibración óptica del CIAI
En lo que respecta a las características específicas que deben reunir los componentes
del sistema son los siguientes:
• El certificado de calibración de la esfera integrante debe ser trazables al estándar de
radiancia primario del National Institute of Standards and Technology (NIST,
http://www.nist.gov/).
• Debe disponer de un sistema de control y adquisición de la intensidad de las
lámparas de la esfera para garantizar que, durante todo el proceso de calibración,
permanece en el rango establecido por el certificado de calibración de la misma.
53
Capítulo 3
Siguiendo este diseño, los componentes básicos del sistema han de ser los
siguientes:
• Un brazo principal de más de 1.0 m de altura con tres brazos horizontales acoplados
(de 40 cm cada uno):
Un primer brazo horizontal para la sujeción del láser-guía de alineación vertical
(eje N-S).
Un segundo brazo horizontal para la sujeción de lámparas halógenas tipo DXW
de 1000 W.
Un tercer brazo de sujeción de bafles ignífugos, con diferentes aperturas, para
minimizar la luz parásita que pueda llegar al detector del instrumento a calibrar.
• Un láser-guía de alineación horizontal que debe estar acoplado en el soporte vertical
del sistema.
• Un sistema mecánico angular que permita que el brazo gire 180º, en un plano
vertical E-W, paralelo al eje óptico vertical del radiómetro, considerando que el
radiómetro se encuentra en la posición sur.
A partir del diseño básico del sistema se procede al análisis pormenorizado de las
características específicas que debe verificar para garantizar el correcto proceso de
calibración.
• Las lámparas halógenas tipo DXW de 1000 W han de estar estacionadas.
• Debe disponer de un sistema de control y adquisición de la intensidad de la corriente
que circula por la lámpara para garantizar que, durante todo el proceso de
calibración, permanece constante.
54
Laboratorio de calibración óptica del CIAI
55
Capítulo 3
Figura 3.6. Esquema del diseño del sistema de calibración de respuesta espectral.
De acuerdo con este diseño, los componentes básicos del sistema han de ser los
siguientes:
• Un sistema de iluminación compuesto por lámparas que emitan en el rango
ultravioleta y en el visible (entre 290 y 700 nm).
• Un sistema dispersivo monocromador simple o doble. A la hora de seleccionar uno
u otro hay que tener en cuenta que el primero se ve afectado por luz parásita o stray
light y en el segundo se produce una disminución de la relación señal/ruido.
• Un módulo óptico de salida que permita acoplar a la salida del monocromador tanto
el sensor calibrado espectralmente como los instrumentos a calibrar.
Partiendo del diseño básico del sistema, se realiza un análisis pormenorizado de las
características específicas que debe reunir el mismo para permitir la calibración de los
distintos instrumentos a calibrar:
• Las lámparas del sistema de iluminación han de ser de Deuterio (200-400 nm) para
el rango ultravioleta y de Tungsteno (250-2500 nm) para el rango visible (Figura
3.7). Cada lámpara debe disponer de una fuente de alimentación independiente.
56
Laboratorio de calibración óptica del CIAI
Figura 3.7. (a) Esquema del diseño del sistema de iluminación de Optronic Laboratories OL 750-
20D/UV. (b) Camino óptico del haz de luz procedente de la lámpara de Deuterio. (c) Camino óptico del
haz de luz procedente de la lámpara de Tungsteno.
Figura 3.8. Esquema del diseño del monocromador doble OL 750-M-D. Los asteriscos rojos indican la
posición del soporte cinemático de las redes de difracción.
57
Capítulo 3
Figura 3.9. Esquema del diseño del módulo óptico de salida OL 750-660 en el que se ha indicado la
ubicación del puerto de entrada del haz de luz dispersada procedente del monocromador doble y los
puertos de salida (P. S.) en los que se pueden acoplar el sensor calibrado y el instrumento a calibrar.
58
Laboratorio de calibración óptica del CIAI
Figura 3.10. (a) Santiago González, técnico de la Universidad de Valladolid, durante la instalación, en el
laboratorio de calibración óptica del Observatorio Atmosférico de Izaña, de los sistemas de calibración
(b1) de respuesta angular, (c) de irradiancia absoluta en posición vertical, (d) de irradiancia absoluta en
posición horizontal y (e) de radiancia absoluta. Además, se muestra (b2) detalle del plato divisor y del
motor paso a paso del sistema mecánico angular.
59
Capítulo 3
60
Laboratorio de calibración óptica del CIAI
Figura 3.13. Bafle original. Figura 3.14. Prototipo de bafle para pruebas.
61
Capítulo 3
Para estas nuevas pruebas se decide usar un piranómetro YES UVB-1 de la casa
Yankee Environmental System, de menor precisión que el espectrofotómetro Brewer
pero más rápido y fácil de controlar, poseyendo las características adecuadas para el
estudio pretendido.
En primer lugar, se realizan una serie de medidas similares a las realizadas con el
espectrofotómetro. Es decir, se fijan las posiciones del difusor del piranómetro y de la
lámpara de calibración y se modifica la posición del bafle original. En este caso se sitúa
en tres posiciones diferentes a 52 cm, 77 cm y 94 cm respecto de la regla del sistema de
calibración vertical (Figura 3.15). Los resultados (Tabla 3.2) indican, por un lado, una
diferencia del 7.34% entre las medidas a 77 cm respecto de las de 55 cm y, por otro
lado, una diferencia del 4.33% entre las medidas a 94 cm respecto de las de 55 cm. Se
confirma de esta manera las variaciones excesivas en la señal medida por el instrumento
en función de la posición del bafle.
62
Laboratorio de calibración óptica del CIAI
A la luz de estos resultados Santiago González diseñó unos nuevos bafles para los
tres sistemas de calibración afectados. El nuevo diseño se caracteriza porque el borde
interno del agujero de cada bafle termina en punta para evitar reflexiones de la luz
procedente de la lámpara (Figura 3.17).
63
Capítulo 3
Figura 3.17. Diseño de los nuevos bafles con el borde interno del agujero terminado en punta para evitar
reflexiones.
Tabla 3.4. Voltajes promedio medidos con el piranómetro YES UVB-1 durante las pruebas de distancia
entre el bafle nuevo y el difusor.
Posición del bafle en la Voltaje promedio
regla del sistema de medido por el
calibración (cm) piranómetro (mV)
52 54.5281
77 54.6096
94 54.5261
La mejoría en los resultados obtenidos respecto de las pruebas realizada con los
antiguos bafles es notable. Las diferencias han pasado de 7.43%, en los casos extremos,
a diferencias como máximo de 0.15% con el nuevo diseño, que entra en el rango de
variación que cabría esperar bajo la influencia de la posición de los bafles. Con estos
resultados se da por concluida la validación de los nuevos bafles.
64
Laboratorio de calibración óptica del CIAI
65
Capítulo 3
Tabla 3.5. Diferentes configuraciones para la esfera de Izaña medidas con el fotómetro solar master
viajero Cimel #329 del GSFC.
Figura 3.18. Calibración de la esfera integrante de Izaña para diferentes configuraciones medidas, sin
filtro de densidad neutra, por el master viajero de AERONET.
66
Laboratorio de calibración óptica del CIAI
Figura 3.19. Calibración de la esfera integrante de Izaña para diferentes configuraciones medidas, con
filtro de densidad neutra, por el master viajero de AERONET.
A priori, se decide usar por defecto la configuración IZO 07, pero bajo diferentes
condiciones de sensibilidad del fotómetro o cambios en la esfera integrante (p. ej.,
cambios en sus lámparas), se pueden usar el resto de las configuraciones. A partir de
este punto ya se puede considerar que el sistema está totalmente operativo para la
realización de las calibraciones de los fotómetros solares Cimel.
67
Capítulo 3
68
Laboratorio de calibración óptica del CIAI
Figura 3.20. Diferentes pasos del proceso de calibración de irradiancia absoluta del espectrofotómetro
Brewer#201
Figura 3.21. Intensidad y voltaje de la lámpara de 1000 W del sistema de calibración de irradiancia
absoluta. Las líneas verticales de color negro indican el intervalo de tiempo en el que se realizó la
calibración del espectrofotómetro Brewer#201
69
Capítulo 3
Figura 3.22. Parámetro Hg step y temperatura del Brewer#201. El parámetro Hg step debe estar en el
intervalo de valores delimitado por el sombreado de color gris. Las líneas verticales de color negro
indican el intervalo de tiempo en el que se realizó la calibración.
70
Laboratorio de calibración óptica del CIAI
Figura 3.24. Espectros bruto (en cuentas) y calibrado (en W/m2/nm) medidos por el Brewer#201 el 21 de
marzo de 2010 a las 13:30 TMG
71
Capítulo 3
Figura 3.25. Evolución temporal de los factores de calibración, en cada longitud de onda, de la esfera
integrante de Izaña realizadas con el fotómetro solar master viajero de AERONET
Cabe destacar que tanto para la calibración de la esfera como para la calibración de
un fotómetro se realizan varias series de medidas, dando por válida una calibración sólo
si la desviación estándar entre ellas es menor del 1%, con el fin de garantizar la
repetitividad de la calibración.
En lo que respecta a las calibraciones del fotómetro solar Cimel#244 se obtienen
multiplicando la radiancia de la esfera (última calibración realizada con el master
viajero) por la señal medida por el instrumento en cada longitud de onda y en cada
canal. Es decir, se realiza una calibración para el canal de las medidas de aureola, en las
proximidades del sol, (Figura 3.26) y otra para el canal utilizado en las medidas de cielo
(Figura 3.27).
Figura 3.26. Evolución temporal de los factores de calibración, en cada longitud de onda, del fotómetro
solar Cimel#244 (master de Izaña) en el canal de aureola
72
Laboratorio de calibración óptica del CIAI
Figura 3.27. Evolución temporal de los factores de calibración, en cada longitud de onda, del fotómetro
solar Cimel#244 (master de Izaña) en el canal de cielo
Los ficheros con los factores de calibración calculados, para cada canal y para cada
longitud de onda, en unidades de µW/(cm2·sr·nm) son enviados a AERONET para que
sean aplicados, en tiempo real y hasta la siguiente calibración, a las medidas rutinarias
de almucántar y plano principal del fotómetro. La incertidumbre asociada a los factores
de calibración calculados se estima en torno a un 5% (Holben et al, 1998).
73
Capítulo 3
Figura 3.28. Descripción de la calibración de la respuesta angular del piranómetro YES UVB-1 #970839
Figura 3.29. Intensidad (color marrón) y voltaje (color verde) aplicados a la lámpara durante la
calibración tanto en el plano NS (panel superior) como en el plano EW (panel inferior)
74
Laboratorio de calibración óptica del CIAI
Tabla 3.6. Función respuesta angular normalizada del piranómetro YES UVB-1 #970839
Ángulo cenital (º) Plano N Plano S Plano E Plano W Promedio
0 1.00000 1.00000 1.00000 1.00000 1.00000
5 0.98822 0.99486 0.99322 0.99804 0.99359
10 0.97682 0.97799 0.98229 0.98491 0.98050
15 0.95112 0.96210 0.96070 0.95903 0.95824
20 0.91687 0.92152 0.92693 0.92266 0.92200
25 0.87066 0.88621 0.88124 0.87987 0.87949
30 0.81769 0.83317 0.83187 0.82876 0.82788
35 0.75995 0.77344 0.77236 0.77192 0.76942
40 0.69364 0.71003 0.70096 0.69812 0.70069
45 0.62260 0.63032 0.62897 0.62997 0.62797
50 0.54042 0.55155 0.55153 0.55116 0.54866
55 0.45157 0.46877 0.46303 0.46379 0.46179
60 0.36164 0.37109 0.37767 0.37845 0.37221
65 0.27306 0.28259 0.27955 0.28870 0.28097
70 0.19002 0.19499 0.19301 0.19863 0.19417
75 0.11756 0.12244 0.12307 0.12399 0.12176
80 0.05578 0.05632 0.04896 0.05965 0.05518
83 0.02184 0.01877 0.01794 0.02728 0.02146
85 0.01476 0.00954 0.00781 0.00896 0.01027
87 0.00589 0.00729 0.00528 0.00567 0.00603
90 0.00550 0.01051 0.01288 0.00955 0.00961
75
Capítulo 3
Figura 3.30. Función respuesta angular normalizada en los cuatro planos del piranómetro YES UVB-1
#970839 representado junto a la ley coseno
El error coseno, calculado como el ratio entre la respuesta angular y la ley coseno,
se muestra en la Figura 3.31. La incertidumbre de la calibración realizada se estima
entre un 1.2% para ángulos cenitales menores y el 10% para ángulos cenitales mayores.
Figura 3.31. Error coseno en los cuatro planos del piranómetro YES UVB-1 #970839
76
Laboratorio de calibración óptica del CIAI
Figura 3.32. Respuesta espectral relativa del piranómetro YES UVB-1 #970839
77
Capítulo 3
78
Laboratorio de calibración óptica del CIAI
79
Capítulo 3
de aerosoles y vapor de agua por la noche (Barreto et al., 2013a; 2013b). Por otro lado,
el CIAI está colaborando con la empresa Sieltec Canarias S.L. (http://www.sieltec.es/)
en el desarrollo de un nuevo instrumento robusto, de bajo coste y mantenimiento,
denominado Digital Sky Color Radiometer (DSCR), para la medida de espesor óptico
de aerosoles, que está siendo probado y caracterizado en el laboratorio.
Finalmente, cabe destacar el papel del laboratorio de calibración óptica en las
siguientes actividades del programa de aerosoles en columna:
• El CIAI contribuye al programa SDS-WAS con cuatro fotómetros solares instalados
en lugares estratégicos en Argelia, Marruecos, Egipto y Túnez, que son totalmente
calibrados en Izaña. En particular, los datos del fotómetro solar Cimel instalado en
Tamanrasset (Argelia) han sido utilizados en el Capítulo 4 de la presente memoria.
• IZO como banco de pruebas del WMO-CIMO ha participado, hasta ahora, en dos
tipos de actividades:
Cursos de entrenamiento y actividades de desarrollo de capacidades sobre
operación, calibración y evaluación de datos de los fotómetros solares a personal
de servicios meteorológicos de países del norte de África.
Colaboración con fabricantes de instrumentación de medida de aerosoles de bajo
coste, proporcionando soporte en el Observatorio de Izaña para la calibración de
dichos instrumentos, como por ejemplo el fotómetro solar Calitoo de TENUM
(http://www.calitoo.fr/).
80
Laboratorio de calibración óptica del CIAI
3.8. Apéndices
Documento creado por Carmen Guirado Fuentes en el marco del Convenio Específico de Colaboración
entre el Instituto Nacional de Meteorología y la Universidad de Valladolid para el establecimiento de
metodologías y sistemas de control de calidad para los programas de fotometría, radiometría, ozono y
aerosoles atmosféricos dentro del programa de Vigilancia Atmosférica Global de la Organización
Meteorológica Mundial
81
Capítulo 3
Documento creado por Carmen Guirado Fuentes en el marco del Convenio Específico de Colaboración
entre el Instituto Nacional de Meteorología y la Universidad de Valladolid para el establecimiento de
metodologías y sistemas de control de calidad para los programas de fotometría, radiometría, ozono y
aerosoles atmosféricos dentro del programa de Vigilancia Atmosférica Global de la Organización
Meteorológica Mundial
82
Laboratorio de calibración óptica del CIAI
1 Operador
2 Fecha
3 Fotómetro
4 Configuración de la esfera integrante
Intensidad de la fuente de alimentación de la lámpara
4A
interna (A)
Intensidad de la fuente de alimentación de la lámpara
4B
externa (A)
4C Valor del picoamperímetro (A)
4D Valor del micrómetro (graduación negra) (mm)
4E Distancia fotómetro-esfera (cm)
4F Distancia pantalla negra-colimadores (cm)
4G Versión de la configuración
5 Comprobación de la limpieza del interior de la esfera
6 Revisiones iniciales
Limpieza de las lentes y los colimadores del instrumento
6A
(si procediese)
6C Saturación de canales
6D Cambio de ganancias
6E Control remoto del instrumento
6F Valor de la batería interna (V)
6G Valor de la temperatura en el interior del instrumento (ºC)
6H Valor de la batería externa (V)
6I Valor de la humedad en el interior del instrumento
6J Envío de los últimos datos de campo al PC
6K Vaciado de la memoria del instrumento
6L Comprobación de la correcta fijación de los colimadores
Documento creado por Carmen Guirado Fuentes en el marco del Convenio Específico de Colaboración
entre el Instituto Nacional de Meteorología y la Universidad de Valladolid para el establecimiento de
metodologías y sistemas de control de calidad para los programas de fotometría, radiometría, ozono y
aerosoles atmosféricos dentro del programa de Vigilancia Atmosférica Global de la Organización
Meteorológica Mundial
83
Capítulo 3
Documento creado por Carmen Guirado Fuentes en el marco del Convenio Específico de Colaboración
entre el Instituto Nacional de Meteorología y la Universidad de Valladolid para el establecimiento de
metodologías y sistemas de control de calidad para los programas de fotometría, radiometría, ozono y
aerosoles atmosféricos dentro del programa de Vigilancia Atmosférica Global de la Organización
Meteorológica Mundial
84
Laboratorio de calibración óptica del CIAI
1 Operador
2 Fecha
3 Instrumento
Encendido de la fuente de alimentación del monocromador
4
doble: Hora (TMG)
Encendido de la fuente de alimentación de la/las lámpara/s que
5
se vaya a usar para la calibración: Hora (TMG)
6 Estabilización de la lámpara (15 minutos)
Selección y colocación de las rendijas y/o aperturas adecuadas
7
para la calibración
Selección de las redes de difracción adecuadas para la
8
calibración
9 Ejecución del software de control del sistema
9A Seguir los pasos del asistente de ejecución del programa
10 Acople al sistema del fotodiodo calibrado
Ejecución del modo Calibration scan para calibrar el sistema
11
con la configuración seleccionada
Desacople del fotodiodo calibrado y acople del instrumento a
12
calibrar
Ejecución del modo Measurement scan para calibrar el
13
instrumento
14 Desacople del instrumento a calibrar
Apagado de la fuente de alimentación de la/las lámpara/s:
15
Hora (TMG)
Apagado de la fuente de alimentación del monocromador
16
doble: Hora (TMG)
Documento creado por Carmen Guirado Fuentes en el marco del Convenio Específico de Colaboración
entre el Instituto Nacional de Meteorología y la Universidad de Valladolid para el establecimiento de
metodologías y sistemas de control de calidad para los programas de fotometría, radiometría, ozono y
aerosoles atmosféricos dentro del programa de Vigilancia Atmosférica Global de la Organización
Meteorológica Mundial
85
Capítulo 3
86
Laboratorio de calibración óptica del CIAI
87
Capítulo 3
88
Capítulo 4
91
Capítulo 4
Figura 4.1. Ejemplo de suciedad depositada en las lentes de un fotómetro solar Cimel que modifica el
camino óptico de la radiación solar hasta el detector.
Este falso ciclo diurno de los datos durante este período es aprovechado por el
método del KCICLO (Cachorro et al., 2004, 2008a) para corregir la serie de medidas.
Este procedimiento fue desarrollado con el fin de detectar, evaluar y corregir posibles
problemas en la calibración de los equipos. En particular, permite corregir el error de
calibración producido por la acumulación de suciedad en las lentes que modifica el
camino óptico de la radiación solar. Este error de calibración se traduce en un error
sistemático absoluto de las medidas de AOD (Romero y Cuevas, 2002). La magnitud de
este error sistemático es mayor a mediodía, ya que varía con la inversa de la masa de
aire que atraviesa la radiación solar, dando lugar a un patrón artificial de ciclo diurno,
que es característico (Figura 4.2) y que puede ser corregido con este método.
Figura 4.2. Ejemplo de patrón artificial de ciclo diurno en el AOD medido en Tamanrasset.
92
Caracterización de aerosoles en Tamanrasset
El método introduce una constante K que se define como el ratio entre la constante
de calibración aplicada (considerada ‘incorrecta’) y la constante de calibración ‘válida’.
K cuantifica el error del factor de calibración de tal modo que K=1 indica que la
constante de calibración es correcta, mientras que K>1 (o K<1) indica una
sobrestimación (o subestimación) de la constante de calibración aplicada y, por tanto,
una forma convexa (o cóncava) del ciclo diurno (Cachorro et al., 2004, 2008a).
Cachorro et al. (2008b) han estimado que las diferencias relativas entre una serie de
nivel 2.0 de AERONET y una corregida con el KCICLO son del 8.5% en AOD
(alrededor de un 0.01 en valor absoluto) y del 2.4% en EA. Este método ha sido
previamente utilizado por diversos autores para corregir series de AOD (p. ej., Toledano
et al., 2007; Barreto et al., 2014).
Sin embargo, la aplicación de este procedimiento in-situ de corrección de la
calibración requiere que haya un número suficiente de días estables y despejados para
cada período de tiempo concreto que quiera ser corregido. Para minimizar el error en la
aplicación del método, los días seleccionados para aplicarlo deberían cumplir los
siguientes requisitos (Cachorro et al., 2008a):
1) Rango de masas de aire mayor que 4, partiendo de una masa de aire mínima
mayor que 1.4 (el rango típico varía entre 1.7 y 6).
2) Espesor óptico de aerosoles menor que 0.12 para el canal de 440 nm con una
variabilidad, en el rango de masas especificado, menor que el 5%.
3) Número de medidas mayor que 12 para cada día.
4) Desviación estándar del ajuste para determinar la constante K menor que 0.01.
Por lo tanto, es recomendable que, para un período concreto, haya disponibles entre
un 5% y un 10% de días que cumplan los requisitos anteriores para poder aplicar con
éxito el método del KCICLO. Esto implica que este procedimiento no siempre es
aplicable con garantías en todos los lugares de medida ni en todos los períodos de
tiempo.
Para el presente estudio se dispuso, entre el 18 de noviembre de 2007 y el 20 de
junio de 2008, de un número suficiente de días (94 en total) que cumplieran los
requisitos para poder aplicar el método del KCICLO y completar la serie de AOD y de
EA en Tamanrasset. El método detectó dos períodos diferentes en cuanto a la desviación
en la constante de calibración, es decir, indicó la existencia de dos tipos de
contaminación en las lentes (con distinta cantidad de suciedad acumulada en cada
lente). Los valores promedio de la constante K para cada período se muestran en la
Tabla 4.1.
93
Capítulo 4
Tabla 4.1. Valores promedio de la constante K y su desviación estándar (σ) para cada canal de medida y
cada período de correccióna.
Una vez que se han determinado los valores promedio de la constante K para cada
período, el término lnK se divide por la masa óptica a la que se realizó cada medida y se
resta del correspondiente valor de AOD o EA que se desea corregir. Como ejemplo de
la corrección realizada, en la Figura 4.3 se muestran tanto la serie original como la serie
corregida de AOD y EA para los dos períodos de corrección. Se aprecia que, en general,
el ciclo diurno ficticio de los datos se ve ampliamente corregido, tanto para el AOD
como para el EA.
Figura 4.3. (a) EA en el rango 440-870 nm y (b) AOD a 440 nm representados, entre el 15 de marzo y el
1 de abril de 2008, tanto antes como después de la aplicación de la corrección de KCICLO. Se aplicaron
dos factores de corrección diferentes, antes y después del 23 de marzo.
94
Caracterización de aerosoles en Tamanrasset
Finalmente, cabe destacar que las medidas de AERONET de nivel 1.5 de PWV
entre octubre de 2006 y febrero de 2009 no se ven afectadas por ningún ciclo diurno
ficticio. El espesor óptico del vapor de agua y, por tanto, el PWV no se ven
prácticamente afectados por obstrucciones en el camino óptico debido a que el
algoritmo para su cálculo se basa en la sustracción de medidas experimentales (capítulo
2, apartado 2.2.4.1). Por lo tanto, la serie de PWV analizada está formada por los datos
de nivel 2.0 junto a los de nivel 1.5 cuando los primeros no estén disponibles.
Figura 4.4. Fotómetro solar manual tipo J-309 utilizado en Tamanrasset para medir AOD tres veces al
día (9, 12 y 15 UTC). Foto cortesía del Sr. Lahouari Zeudmi.
Las medidas realizadas, entre octubre de 2006 y febrero de 2009, en el canal de 500
nm a las 9, 12 y 15 UTC (Tiempo Universal Coordinado) son las que se van a usar para
comparar con las medidas del fotómetro solar Cimel. Para ello se han seleccionado las
medidas de AERONET en el canal de 440 nm que sean más próximas, en un intervalo
de ±15 minutos, a las del fotómetro manual. La comparación (Figura 4.5) se ha
realizado para tres series distintas de datos del fotómetro Cimel: 1) la serie de nivel 2.0
95
Capítulo 4
comprendida entre octubre de 2006 y febrero de 2009, 2) la serie de nivel 1.5 entre
noviembre de 2007 y junio de 2008 sin aplicar ningún tipo de corrección, y 3) la serie
de nivel 1.5 entre noviembre de 2007 y junio de 2008 tras aplicar la corrección de
KCICLO.
Figura 4.5. Correlación entre las medidas coincidentes (en un intervalo de ±15 minutos) de AOD a 440
nm de AERONET y a 500 nm del fotómetro manual de la NOAA. La serie de nivel 2.0 de AERONET
(verde claro) abarca el período entre octubre de 2006 y febrero de 2009. La serie de AERONET de nivel
1.5 abarca el período entre noviembre de 2007 y junio de 2008 tanto sin aplicar la corrección de KCICLO
(gris) como tras aplicarla (verde oscuro). Cada grupo de datos se representa con su correspondiente
regresión lineal. La línea continua negra corresponde a la línea de referencia 1:1.
96
Caracterización de aerosoles en Tamanrasset
Tabla 4.2. Resultados del ajuste lineal por mínimos cuadrados entre las medidas de AOD en 500 nm del
fotómetro manual de la NOAA y las de 440 nm del fotómetro Cimel para tres grupos de datos (nivel 1.5
de AERONET, tanto antes como después de aplicar el KCICLO, y nivel 2.0). Los parámetros que se
muestran son los siguientes: pendiente del ajuste, ordenada en el origen, coeficiente de correlación (R2),
raíz cuadrada del error cuadrático medio (RECM) y el número de observaciones coincidentes (±15
minutos) entre ambas series.
AERONET nivel AERONET nivel AERONET
1.5 pre-KCICLO 1.5 post-KCICLO nivel 2.0
Pendiente 1.15±0.02 1.07±0.01 1.02±0.01
Ord. en el origen 0.031±0.006 -0.014±0.004 0.001±0.001
2
R 0.968 0.981 0.983
RECM 0.044 0.031 0.024
Nº observaciones 450 450 1241
AOD(λ ) = β ⋅ λ− EA (4.1)
97
Capítulo 4
− EA ( 440 −870 nm )
500
AOD(500 nm) = AOD(440 nm) ⋅ (4.2)
440
Las diferencias relativas entre los canales de 500 nm y 440 nm son, en general,
menores de 0.01, excepto para valores de AOD superiores a 0.1 en que, en algunos
casos, dichas diferencias relativas pueden alcanzar 0.04 como máximo. Debido a estas
pequeñas diferencias, en lugar de analizar los datos estimados en 500 nm, se han elegido
los datos medidos en 440 nm como representativas para el análisis del AOD ya que,
además, el canal de 440 nm es uno de los que está presente en todos los modelos de
fotómetro solar CE-318.
Las estadísticas mensuales del AOD, EA y FMF se muestran en la Figura 4.6. El
AOD permanece relativamente estable y bajo, en torno a 0.1, entre noviembre y febrero,
con un mínimo absoluto de 0.07 en enero. En cambio, el AOD sobrepasa el valor de 0.3
entre abril y septiembre, alcanzando el máximo absoluto en junio con un valor de 0.43
(Figura 4.6a y Tabla 4.3).
Figura 4.6. Diagrama de cajas de los promedios mensuales correspondientes a los datos diarios de (a)
AOD en el canal de 440 nm, (b) EA en el rango 440-870 nm y (c) FMF en 500 nm durante el período en
estudio. Los círculos corresponden a los valores promedio, el área sombreada en gris indica el rango de
valores comprendidos entre el promedio más y menos la desviación estándar, las cajas corresponden a los
percentiles 25 y 75, y las barras verticales se extienden hasta el valor extremo dentro de 1.5 veces el rango
intercuartil.
98
Caracterización de aerosoles en Tamanrasset
Hay una mayor variabilidad del AOD (desviación estándar mayor de 0.30) en los
meses en los que se registran los valores más altos (entre abril y agosto, excepto mayo)
mientras que la menor variabilidad (desviación estándar en torno a 0.10) coincide con
las medidas más bajas de AOD (entre noviembre y enero). La presencia de valores muy
bajos de AOD entre noviembre y febrero ya había sido previamente detectada por
Cuesta et al. (2008), en su análisis centrado en el año 2006, y por Kim et al. (2011), en
un estudio de las propiedades ópticas del polvo mineral en diferentes estaciones, entre
ellas Tamanrasset.
Tabla 4.3. Promedios mensuales del AOD, el EA y el PWV en Tamanrasset entre octubre de 2006 y
febrero de 2009a.
El ciclo anual del EA y de la FMF es la opuesta del AOD (Figura 4.6). Entre
noviembre y febrero se registran los valores más altos de EA y FMF, alcanzando un
máximo en diciembre (0.72 y 0.57, respectivamente) y en enero (0.69 y 0.58,
respectivamente). Estos valores decrecen progresivamente hasta alcanzar un mínimo en
mayo (0.15 y 0.24, respectivamente). En agosto se observa un máximo secundario, con
valores de 0.44 para el EA y de 0.38 para la FMF, que está relacionado con un descenso
del modo grueso y un ligero incremento del modo fino. Este incremento se estudiará en
detalle en el apartado 4.5.
Respecto a los patrones estacionales mostrados en la Figura 4.6, Cuesta et al. (2008)
identificaron una marcada evolución estacional del contenido de aerosol atmosférico y
sus propiedades ópticas, que relacionaron con el régimen monzónico a lo largo del año
2006. Además, expusieron la influencia de las propiedades de las masas de aire
transportadas hasta Tamanrasset sobre la evolución de la altura de la capa límite
99
Capítulo 4
Figura 4.7. Series temporales de (a) la altura de la CLC calculada a partir de los sondeos lanzados a las
12 UTC en Tamanrasset y los valores promedio diarios de (b) PWV, (c) AOD a 440 nm y (d) EA entre
440 y 870 nm. Las líneas continuas corresponden a las medias móviles a 30 días.
100
Caracterización de aerosoles en Tamanrasset
noviembre y febrero). Cabe destacar que, en el intervalo de tiempo analizado, los meses
de marzo y octubre han actuado como meses de transición entre las dos estaciones
principales, sin identificarse claramente con ninguna de ellas. La estación húmeda está
modulada por una fuerte y gruesa CLC (Figura 4.7.a) y se ve afectada por vientos
húmedos del este relacionados con el régimen monzónico. Como consecuencia, en este
período, fuertes tormentas de polvo son frecuentes en Tamanrasset, cuando la profunda
CLC favorece la mezcla vertical de las capas de polvo que se elevan del suelo (Cuesta et
al., 2009). Los valores más altos de AOD y PWV (Figura 4.7) coinciden con la época en
la que la CLC está completamente desarrollada (entre 4 y 6 km s.n.m.). Por el contrario,
durante el resto del año, el flujo predominante de vientos secos del oeste conduce a una
CLC poco profunda (entre 1 y 2 km s.n.m.) y a registros menores de AOD y PWV.
Estos resultados coinciden con Cuesta et al. (2008), que describieron, para el año 2006,
una estación estival modulada por una CLC de 5 a 6 km de profundidad que evolucionó
a una menos profunda (entre 1.5 y 2 km) en el invierno.
En este contexto, se puede establecer que la estación seca y fría se caracteriza por
valores bajos de AOD (~0.09 en 440 nm) y de PWV (~0.51 cm) y valores relativamente
altos de EA (~0.62). En cuanto a la estación húmeda y cálida, se caracteriza por valores
más altos de AOD (~0.39), más bajos de EA (~0.28) y se doblan los valores de PWV
con respecto a la estación seca (~1.06 cm). En la Tabla 4.4 se muestran las estadísticas
correspondientes a estos dos períodos.
Tabla 4.4. Estadísticas del AOD, EA y PWV entre octubre de 2006 y febrero de 2009 en Tamanrasseta.
101
Capítulo 4
Las distribuciones de frecuencia del AOD y del EA para las estaciones seca y
húmeda se muestran en la Figura 4.8. El AOD muestra una distribución unimodal
asimétrica positiva para ambas estaciones. El valor modal de AOD para la estación
húmeda es 0.15 (con el 35% de las medidas por debajo de este valor), mientras que el
modo para la estación seca es más estrecho (con un 90% de las medidas por debajo de
0.15) y está centrado en 0.1. Estas cifras indican que Tamanrasset tiene una atmósfera
más limpia que aquellas estaciones situadas en el cercano Sahel, en las que tan solo el
15% de las medidas de AOD están por debajo 0.15 (Basart et al., 2009). Este hecho
podría ser parcialmente explicado por la altura a la que se encuentra esta estación.
En lo que respecta al EA, la distribución es bimodal tanto para la estación seca como
para la húmeda. La distribución para la estación seca es ligeramente bimodal (con
valores modales 0.4 y 0.7) y simétrica (la media y la mediana del EA es la misma, Tabla
4.4). Sin embargo, la distribución de la estación húmeda es simétrica positiva
mostrando, por un lado, un primer modo estrecho y centrado en 0.15 y, por otro lado, un
modo más ancho centrado en 0.4 (coincidiendo con el primer valor modal de la estación
seca).
Figura 4.8. Frecuencia relativa (FR) y frecuencia acumulada (FA) de (a) AOD a 440 nm y (b) EA entre
440 y 870 nm en Tamanrasset. Se muestran por separado los histogramas correspondientes a la estación
seca y a la estación húmeda de cada uno de los parámetros.
102
Caracterización de aerosoles en Tamanrasset
Figura 4.9. Diferencia del exponente de Ångström δAE=EA(440-670 nm)-EA(670-870 nm) en función
del EA y del AOD en la estación de Tamanrasset (10460 observaciones). Los eventos relacionados con
fuertes intrusiones de polvo mineral están indicados mediante un recuadro rojo mientras que la mezcla de
diferentes tipos de aerosol se indica mediante un recuadro verde.
Las características de los aerosoles en Tamanrasset (Figura 4.9) son similares a las
descritas por Basart et al. (2009) en otras zonas áridas y desérticas, como Banizombou
(Níger) o Saada (Marruecos). Se aprecian grandes variaciones del AOD respecto del EA
con un comportamiento, en general, inversamente proporcional entre ambos, lo cual
indica que las partículas mayores son las que producen una mayor extinción de la
radiación solar. Además, δAE es negativo o ligeramente positivo, lo cual indica que el
modo de un único tipo de partícula es el dominante en Tamanrasset. Las condiciones
típicas de puro polvo sahariano (recuadro rojo de la Figura 4.9) se caracterizan por
valores altos de extinción (AOD>0.7), con EA<0.3 y δAE<0, que corresponden a
FMF<40% y Rf~0.3µm. Cabe destacar que, para EA entre 0.6 y 1.1 se aprecian AODs
mayores de los esperados (casi hasta 0.4), lo cual ocurre en el 8.7% de los casos
(recuadro verde de la Figura 4.9). Estas medidas se caracterizan por valores variables de
δAE, FMF y Rf, que oscilan entre -0.3 y 0.2, 30% y 70%, y 0.10 y 0.20 µm,
respectivamente. Estas características pueden ser relacionadas con una mezcla de polvo
mineral y partículas más pequeñas de otro origen (Basart et al., 2009) y se observan sólo
103
Capítulo 4
104
Caracterización de aerosoles en Tamanrasset
La fracción en volumen diaria del modo fino (Vf/Vt) varía entre 0.03 y 0.46 (Figura
4.10b) indicando el predominio del modo grueso. Sin embargo, como ya se discutió en
los apartados 4.2.1 y 4.2.2 al analizar la FMF, tanto las partículas finas como las gruesas
dominan, dependiendo de la estación del año, en cuanto a su contribución al AOD total.
La relación entre estos dos parámetros cuantitativos del modo fino se muestra en la
Figura 4.10b. A diferencia de la estación seca, en la estación húmeda las partículas
gruesas dominan tanto en términos del espesor óptico como de la concentración en
volumen, ya que hay muy pocas medidas que verifiquen FMF>0.5 y Vf/Vt >0.25. La
mayor parte de estas últimas medidas están relacionadas con la presencia de aerosoles
finos en la estación húmeda que será analizada en el apartado 4.5.
Con respecto al radio efectivo total (Tabla 4.5), tiene un esperado comportamiento
estacional opuesto al EA, con un máximo en mayo (0.86 µm), un mínimo en
noviembre-diciembre (~0.58 µm) y un mínimo secundario en agosto (0.61 µm). El radio
efectivo del modo fino alcanza un máximo en la estación seca (~0.16 µm en enero y
diciembre), decreciendo hasta su valor más bajo en julio y agosto (~0.12 µm), lo cual
indica un comportamiento estacional muy próximo al opuesto del AOD. De forma
similar, el radio efectivo del modo grueso muestra los valores más altos en enero (1.92
µm) y parece mantenerse bastante estable durante la estación húmeda, oscilando entre
1.62 y 1.72 µm. Este último resultado, la disminución del radio efectivo del modo
grueso al aumentar sus concentraciones bajo condiciones de polvo desértico (Tabla 4.5),
ha sido previamente enunciado y relacionado con una distribución en volumen de
tamaño de partículas que sea prácticamente monomodal (p. ej. Prats et al., 2011 y las
referencias que contiene).
105
Capítulo 4
106
Caracterización de aerosoles en Tamanrasset
Tabla 4.6. Promedios estacionales del SSA, la parte real e imaginaria del índice de refracción (Ind. Ref.)
y parámetro de asimetría (Asim.) para los canales 440, 675, 870 y 1020 nma junto con el número de días
con observaciones disponibles (N).
107
Capítulo 4
Figura 4.11. Correlación entre las medidas de vapor de agua precipitable del Cimel y el radiosondeo
Vaisala RS92 en Tamanrasset para un total de 610 medidas coincidentes. La línea negra discontinua
corresponde a la regresión lineal.
108
Caracterización de aerosoles en Tamanrasset
109
Capítulo 4
Figura 4.12. Promedio del AOD en 555 nm del MISR en el período 2007-2008 durante (a) la estación
seca (entre noviembre y febrero) y (b) la estación húmeda (entre abril y septiembre). La posición
geográfica de Tamanrasset se indica mediante una estrella negra mientras que los círculos negros indican
estaciones AERONET que presentan series continuas de datos hasta la actualidad (2014). Diferentes
zonas fuente de polvo mineral, que se nombran en el texto, aparecen identificadas (líneas rojas/azules) y
numeradas en el mapa de la siguiente manera: 1, triángulo formado por los macizos montañosos Adrar
des Ifoghas, Hoggar y Aïr; 2, desierto libio oriental; 3, depresión del Bodélé; 4, zona del Sáhara
occidental; y 5, frontera entre Libia y Túnez. El promedio (línea negra) y la mediana (línea negra
discontinua) de los coeficientes de extinción en Tamanrasset del sensor CALIOP en 532 nm se muestran
en el período 2007-2008 para (c) la estación seca (43 perfiles disponibles) y (d) la estación húmeda (95
perfiles disponibles). El área sombreada en gris corresponde al rango entre los percentiles 20 y 80. La
línea roja indica el umbral para condiciones de atmósfera limpia (coeficiente de extinción < 0.05 km-1).
Los niveles 2600 y 5600 m s.n.m., que son significativos para el análisis CWT, se muestran mediante
líneas negras punteadas.
110
Caracterización de aerosoles en Tamanrasset
Figura 4.13. Estación seca (de noviembre a febrero) en Tamanrasset: Mapas resultantes de la aplicación
del método CWT a nivel de superficie para (a) el AOD y para (b) el EA. Retrotrayectorias HYSPLIT con
puntos de llegada a (c) 2600 m s.n.m. y (d) 5600 m s.n.m. Tamanrasset se encuentra localizada en el
punto de intersección de los cuatro cuadrantes.
Durante la estación seca, las retrotrayectorias de las masas de aire a 2600 y 5600 m
s.n.m muestran, en general, una procedencia del oeste (Figura 4.13c, d), moduladas por
la circulación general de la atmósfera ya que a estas alturas dominan las condiciones de
troposfera libre sobre el, relativamente bajo, tope de la CLC. Esta estación del año se
caracteriza por valores bajos de AOD y más bien altos del EA, asociados con
retrotrayectorias a nivel de superficie, procedentes del primer cuadrante (Figura 4.13a,
b).
111
Capítulo 4
112
Caracterización de aerosoles en Tamanrasset
Figura 4.14. Estación húmeda (de abril a septiembre) en Tamanrasset: Mapas resultantes de la aplicación
del método CWT para (a) el AOD a nivel de superficie, (b) el EA a nivel de superficie, (c) el AOD a 2600
m s.n.m. y (d) el EA a 2600 m s.n.m; (e) retrotrayectorias HYSPLIT con punto de llegada a 5600 m
s.n.m. Tamanrasset se encuentra localizada en el punto de intersección de los cuatro cuadrantes.
En la estación húmeda, el nivel de 5600 m s.n.m. está sobre el tope de la CLC que es
bastante alto en esta época del año (Figura 4.14e). Los resultados del análisis del
método CWT a nivel de superficie y a 2600 m s.n.m son muy similares, tanto para el
AOD como para el EA, lo cual sugiere una CLC bien mezclada. Estos resultados
(Figura 4.14a-d) muestran un recorrido curvo principal de las masas de aire con valores
relativamente altos del EA (partículas más pequeñas) desde el norte de Libia central que
pasa sobre las fuentes de polvo 1 y 2, al igual que en la estación seca. Se observa un
segundo recorrido curvo procedente de la frontera entre Libia y Túnez (zona fuente 5,
Figura 4.12) identificada previamente (p. ej. Caquineau et al., 2002) como fuente de
polvo mineral. Estas últimas masas de aire transportan, a nivel de superficie, partículas
113
Capítulo 4
más grandes (EA más bajos) hasta Tamanrasset. Finalmente, algunas trayectorias
proceden del oeste pasando sobre una extensa área (zona fuente 4, Figura 4.12) que
abarca el norte de Mali, el norte de Mauritania y el oeste de las montañas Hoggar
(Prospero, 2002; Brooks and Legrand, 2003, Alonso-Pérez et al., 2012). Se trata de una
compleja distribución de fuentes de polvo mineral que es particularmente activa entre
abril y septiembre.
En lo que respecta a la depresión del Bodélé (zona fuente 3, Figura 4.12)
identificada como una de las fuentes de polvo mineral más importantes del mundo (p.
ej. Goudi and Middleton, 2001; Prospero, 2002; Brooks and Legrand, 2003), el análisis
del método CWT muestra que es una fuente menor de polvo mineral para Tamanrasset.
114
Caracterización de aerosoles en Tamanrasset
Figura 4.15. Promedios diarios de AOD de AERONET (círculos verdes) y del modelo NMMB/BSC-
Dust (círculos negros) en el período 2007-2008.
115
Capítulo 4
Figura 4.16. Composición a 550 nm del algoritmo Deep Blue del sensor MODIS para (a) el AOD, en los
21 días de máximo AOD (Dmax) en Tamanrasset afectados por SCMs, y para la anomalía promedio del
AOD correspondiente a (b) los 21 días de máximo AOD y (c) el día anterior a cada evento (Dmax-1). La
ubicación de Tamanrasset está indicada por una estrella negra. Las retrotrayectorias HYSPLIT a dos días
con llegada a Tamanrasset a nivel de superficie (líneas negras) se muestran también en el panel (a).
El análisis del mapa promedio de AOD del MODIS Deep Blue (Figura 4.16a), para
los 21 días de máximo AOD en Tamanrasset debidos a SCMs, indica que hay diversas
regiones con valores altos de AOD, incluidos los alrededores de Tamanrasset. Sin
embargo, se aprecia que es sólo al sur de Tamanrasset donde existe una fuerte anomalía
positiva de AOD, con valores superiores a 0.20 (Figura 4.16b). La anomalía se produce
en una zona que coincide con la fuente de polvo mineral 1 de la Figura 4.12, y sus
alrededores, y es debida la presencia de SCMs modulados por el desplazamiento hacia
el norte de la discontinuidad intertropical. Las retrotrayectorias HYSPLIT (Figura
4.16a) muestran que el flujo de masas de aire dominantes en Tamanrasset proceden de
esta región con anomalía positiva del AOD. Al mismo tiempo, se observa (Figura
4.16b) una anomalía negativa en la zona este de Mali que es, probablemente, debida a
fenómenos de precipitación asociados a los SCMs ya que esta anomalía estaba situada
más al este, en Níger, en la Figura 4.16c correspondiente a la anomalía promedio del día
116
Caracterización de aerosoles en Tamanrasset
117
Capítulo 4
para cada una de las retrotrayectorias HYSPLIT a cinco días, en el período 2007-2008,
con llegada en Tamanrasset tanto a nivel de superficie como a 2600 m s.n.m. De entre
todas las retrotrayectorias con índice de permanencia superior al 10% en alguna de las
regiones predefinidas, se han seleccionado solo aquellas con un valor del EA superior a
0.7 (Tabla 4.7) como casos de estudio de transporte de aerosol fino antropogénico hasta
Tamanrasset.
Tabla 4.7. Características de los casos de estudio de transporte de aerosol fino antropogénico hasta
Tamanrasset.
118
Caracterización de aerosoles en Tamanrasset
Figura 4.17. Retrotrayectorias HYSPLIT, con punto de llegada en Tamanrasset (marcador azul), a nivel
de superficie (líneas rojas) y a 2600 m s.n.m. (líneas amarillas) para diversos días en estudio. El fondo del
mapa corresponde al producto Nighttime Lights del DMSP, que identifica las llamaradas procedentes de
la combustión de gas mediante el color verde.
119
Capítulo 4
120
Caracterización de aerosoles en Tamanrasset
121
Capítulo 5
Tabla 5.1. Fotómetros solares Cimel que han operado en la estación de Santa Cruz de Tenerife
etiquetados con el número asignado por la red AERONET e indicando el número de días que estuvieron
instalados.
En este período se dispone de 1964 días con observaciones válidas libres de nubes,
sobre un total de 2747 días, lo cual supone un 72% de los días disponibles (Figura 5.1).
Cabe mencionar la presencia de dos intervalos de tiempo en los que no hay datos
disponibles debido a la ausencia de fotómetro en esta estación: entre el 8 de mayo y el
14 de agosto de 2007 y entre el 14 de julio y el 25 de noviembre de 2010.
Figura 5.1. Serie temporal de promedios diarios y mensuales de AOD a 500 nm en la estación de Santa
Cruz de Tenerife.
125
Capítulo 5
5.1.2. Izaña
Las primeras observaciones AERONET en la estación de Izaña (IZO) se realizaron
en los meses de junio y julio de 1997, en el marco de la segunda campaña Aerosol
Characterization Experiment (ACE-2; Raes et al., 2000). Sin embargo, la calidad de
dichos datos es cuestionable ya que algunos de los canales de medida presentan un
ficticio ciclo diurno mientras que en otros canales la relación entre los espesores ópticos
no verifica la ley de Ångström. Al tratarse de una serie corta de medidas aislada
temporalmente del resto de las observaciones, se ha decidido no intentar corregir ni
tener en cuenta dichas medidas para el presente estudio.
La serie continuada de datos de nivel 2.0 de IZO, que está disponible a finales de
2014 en AERONET, comprende el intervalo entre el 26 de junio de 2004 y el 30 de
diciembre de 2012. En este período se dispone de 2479 días con observaciones válidas
libres de nubes, sobre un total de 3110 días, lo cual supone un 80% de los días
disponibles (Figura 5.2). El listado de los fotómetros utilizados para obtener dicha serie
y los períodos de tiempo correspondientes se muestran en la Tabla 5.2.
Hay que tener en cuenta que la estación de Izaña es un centro de calibración
internacional de fotómetros solares Cimel. Por este motivo, la mayoría de los
fotómetros usados para la serie de nivel 2.0 son fotómetros de referencia (maestros) que
estuvieron allí instalados para su calibración al mismo tiempo que otros fotómetros
maestros. Esto dificulta en gran manera la obtención de una serie de nivel 2.0 debido a
la gran variedad de instrumentos disponibles al mismo tiempo y la permanencia
relativamente corta, entre unas pocas semanas y unos pocos meses, de cada uno de ellos
en la estación (menor que el período de operación habitual en las estaciones en las que
se realizan medidas rutinarias, que es de un año). La obtención de esta serie por parte de
AERONET ha requerido una ardua labor de selección de los mejores datos disponibles
para componerla. La ventaja, en cambio, es que los fotómetros maestros están
sometidos a controles de calidad más estrictos y son más precisos (~0.005 en AOD) que
los fotómetros de campo destinados a medidas rutinarias (~0.01 en AOD).
126
Caracterización de aerosoles en Tenerife
Tabla 5.2. Fotómetros solares Cimel que han operado en la estación de Izaña, etiquetados con el número
asignado por la red AERONET e indicando el número de días que estuvieron instalados.
127
Capítulo 5
Figura 5.2. Serie temporal de promedios diarios y mensuales de AOD a 500 nm en la estación de Izaña.
128
Caracterización de aerosoles en Tenerife
Figura 5.3. Diagrama de cajas de los promedios mensuales correspondientes a los datos diarios de AOD
en el canal de 500 nm en (a) IZO y (d) SCO, EA en el rango 440-870 nm en (b) IZO y (e) SCO, y PWV
en (c) IZO y (f) SCO durante el período en estudio. Los círculos corresponden a los valores promedio, el
área sombreada en gris indica el rango de valores comprendidos entre el promedio más y menos la
desviación estándar, las cajas corresponden a los percentiles 25 y 75, y las barras verticales se extienden
hasta el valor extremo dentro de 1.5 veces el rango intercuartil.
Los ciclos anuales de AOD, EA y PWV en IZO y SCO comparten una serie de
rasgos comunes. En lo que respecta al AOD, entre noviembre y febrero se registran los
menores valores promedio. En IZO dichos valores son estables y muy bajos, propios de
129
Capítulo 5
una atmósfera muy limpia (~0.02). En SCO el AOD está en torno a 0.12 durante los
cuatro meses, con un mínimo en febrero (~0.11), coincidiendo con valores propios de
una atmósfera dominada por el aerosol marino (<0.15, Dubovik et al., 2002). Los
valores promedios máximos de AOD se producen en julio (~0.16 en IZO y ~0.27 en
SCO) y agosto (~0.14 en IZO y ~0.26 en SCO), coincidiendo con los valores mínimos
de EA tanto en julio (~0.60 en IZO y ~0.56 en SCO) como en agosto (~0.58 en IZO y
~0.54 en SCO). Cabe recordar que, en esta época del año, se dan las condiciones
meteorológicas que favorecen el transporte de polvo mineral desde África a Canarias
por encima de la capa de inversión térmica, como se describe en el capítulo 1, apartado
1.5.2. La alta variabilidad de las medidas de AOD (desviación estándar muy alta) en los
meses de julio y agosto, tanto en Santa Cruz como en Izaña, pone de manifiesto la
alternancia entre las correspondientes condiciones de fondo e intensas intrusiones de
polvo mineral. Los resultados obtenidos para los promedios mensuales de AOD son
similares a los presentados por García et al. (2014a) en su análisis de los datos de Izaña
entre marzo de 2009 y agosto de 2012. Sin embargo, el valor mínimo del EA obtenido
por ellos en julio y agosto es mayor (~0.64) que el del presente estudio. En este último,
al analizarse un período de tiempo mayor, el efecto de la variabilidad interanual en los
promedios mensuales se ve atenuado.
Aparte de los patrones previamente identificados, cabe destacar en marzo un
máximo secundario de AOD (~0.06 en IZO y ~0.19 en SCO) y una disminución de EA
(~1.06 en IZO y ~0.68 en SCO), coincidiendo con el transporte de polvo mineral en el
seno de la MBL (capítulo 1, apartado 1.5.2.) que, en ocasiones, puede alcanzar los 3 km
de altura (Cuevas et al., 2015), afectando a IZO. Por lo tanto, las intrusiones de polvo
mineral en esta época afectan más a la MBL que a la FT (Alonso-Pérez et al., 2007).
En lo que respecta al PWV, en febrero se registran, en ambas estaciones, los
promedios mensuales mínimos (~0.21 cm en IZO y ~1.48 cm en SCO). Durante la
primavera y el verano el contenido de PWV va aumentando en la atmósfera hasta
alcanzar su máximo en agosto en IZO (~0.65 cm) y en septiembre en SCO (~2.52 cm).
Estos resultados coinciden con los obtenidos por Romero et al. (2011) en su análisis de
los radiosondeos lanzados en la isla de Tenerife entre 1994 y 2011. Estos autores
indicaron que, en los meses invernales, la mayor cantidad de PWV se encuentra por
debajo de los 1500 m s.n.m mientras que en verano el incremento generalizado de PWV
se produce en distinta proporción entre este primer estrato y los estratos superiores,
siendo mayor en estos últimos. La explicación aportada por estos autores sobre el
incremento estival del PWV está relacionada, por un lado, con el aumento de la
temperatura superficial del mar, el correspondiente incremento de vapor de agua cerca
de la superficie marina y su mejor distribución a lo largo de la columna atmosférica. Por
otro lado, las masas de aire africanas que alcanzan Tenerife en niveles altos en verano
aportan una mayor humedad relativa que las masas de aire en condiciones de fondo
procedentes de la troposfera media y alta sobre el Atlántico Norte (Cuevas, 1995). En lo
que respecta a los valores máximos de PWV en agosto y septiembre, Romero et al.
130
Caracterización de aerosoles en Tenerife
131
Capítulo 5
Figura 5.4. Frecuencia relativa y frecuencia acumulada de AOD en el canal de 500 nm en (a) IZO y (c)
SCO, y EA en el rango 440-870 nm en (b) IZO y (d) SCO. En las representaciones de EA se han incluido
los valores mínimos y los valores modales asociados a los histogramas.
132
Caracterización de aerosoles en Tenerife
las utilizadas por otros autores. Por ejemplo, para IZO Basart et al. (2009) identificaron
la presencia de partículas gruesas para AOD(675 nm)>0.15 mientras que García et al.
(2014a) consideraron los valores AOD(500 nm)>0.1 y EA<0.75 como condiciones
representativas de la presencia de polvo mineral. En lo que respecta a SCO, Cuevas et
al. (2015) establecieron el criterio EA(440-670 nm)<0.75 para discriminar los días en
los que las masas de aire procedentes del Sáhara fueran ricas en polvo mineral y
Córdoba-Jabonero et al. (2011) seleccionaron los valores AOD(500 nm)>0.15 y
EA(440-670 nm)<0.5 para determinar condiciones de presencia de polvo mineral de
moderadas a altas.
A partir de los criterios experimentales establecidos se realiza una separación de las
observaciones de las estaciones IZO y SCO para cada uno de los casos a estudiar. Las
estadísticas sobre el número de días por mes a los que corresponde cada caso se muestra
en la Figura 5.5.
Figura 5.5. Diagrama de cajas de los promedios mensuales correspondientes al número de días por mes
en los que las condiciones atmosféricas se corresponden a: (a) FT y FT-SAL en IZO y (b) MBL y MBL-
SAL en SCO durante el período en estudio. Los círculos corresponden a los valores promedio, las cajas
corresponden a los percentiles 25 y 75, y las barras verticales se extienden hasta el valor extremo dentro
de 1.5 veces el rango intercuartil.
133
Capítulo 5
a la ausencia de datos en 2007 (apartado 5.1.1). En el resto de los años estudiados los
valores extremos en mayo oscilan entre 13 y 28 días al mes. Este hecho ya había sido
previamente observado por Cuevas et al. (2015) en su análisis de los perfiles verticales
de extinción del LIDAR de SCO, indicando que los meses de mayo y junio se
caracterizan por ser el período con una atmósfera más limpia en la que se producen muy
pocos episodios de transporte de polvo mineral. Entre noviembre y mayo los promedios,
ligeramente superiores para MBL-SAL que para FT-SAL (Figura 5.5), indican que las
intrusiones de polvo mineral afectan más a la MBL que a la FT, es decir, afectan más a
las capas bajas que a las capas altas de la baja troposfera.
Una vez que las observaciones diarias ya han sido identificadas con cada una de las
condiciones descritas anteriormente, se procede a la caracterización de las mismas.
134
Caracterización de aerosoles en Tenerife
Tabla 5.3. Estadísticas del AOD, EA y PWV entre julio de 2005 y diciembre de 2012 para FT, FT-SAL y
condiciones de mezcla en Izañaa.
Las masas de aire oceánicas de la MBL se caracterizan (Tabla 5.4) por valores
relativamente bajos de AOD (~0.09), partículas de tamaño intermedio (EA~0.95),
debido a la mezcla entre el aerosol fino y el marino (más grueso), y un alto contenido en
PWV (~1.84 cm). En el período estudiado, la MBL se ve afectada por la SAL en menos
de la mitad de los casos (441 días) que aquellos en los que la MBL se encuentra sin
perturbar salvo por la contaminación local (1063 días). Bajo las condiciones de MBL-
SAL (Tabla 5.4), el AOD asciende a 0.38, aumenta significativamente el tamaño de las
partículas (EA~0.29) y aumenta el contenido de PWV (~2.04 cm). Se observa que el
tamaño promedio de las partículas es mayor en IZO que en SCO bajo la influencia de la
SAL, principalmente debido a que la mezcla de polvo mineral con partículas finas es
mayor en SCO.
Los casos de mezcla tanto en IZO como en SCO requieren un análisis en mayor
profundidad que será realizado en próximos trabajos.
135
Capítulo 5
Tabla 5.4. Estadísticas del AOD, EA y PWV entre julio de 2005 y diciembre de 2012 para MBL, MBL-
SAL y condiciones de mezcla en Santa Cruz de Tenerifea.
136
Caracterización de aerosoles en Tenerife
Figura 5.6. Diagrama dispersión de AOD y EA en (a) IZO y (d) SCO, y particularizado para (b) FT y (e)
MBL. Diferencia del exponente de Ångström δAE=EA(440-670 nm)-EA(670-870 nm) en función del EA
y del AOD (código de colores) para (c) FT-SAL y (f) MBL-SAL.
137
Capítulo 5
138
Caracterización de aerosoles en Tenerife
139
Capítulo 5
Las condiciones de FT muestran una ligera bimodalidad con predominio del modo
fino (Vf/Vt~0.73). Presenta una distribución más estrecha del modo fino (valor modal
del radio geométrico en torno a 0.15 µm) y una distribución más ancha del modo
grueso, con un Ref~1.98 µm, y el punto de inflexión en torno a 0.76 µm. En cambio, en
condiciones de FT-SAL, el modo dominante es el grueso (Vf/Vt~0.12) y el punto de
inflexión decrece (~0.44 µm), dando lugar a unos valores de concentración en volumen
de ~0.126 µm3·µm-2 para el modo grueso frente a ~0.015 µm3·µm-2 para el modo fino.
En condiciones de MBL, la distribución de tamaños de partículas es bimodal con
dominio del modo grueso (Vf/Vt~0.28), debido a la presencia de las partículas de
aerosol marino, con concentraciones en volumen ~0.010 µm3·µm-2 y ~0.029 µm3·µm-2
para los modos fino y grueso respectivamente. Esta distribución de tamaños es similar a
la obtenida por Smirnov et al. (2002) para dos estaciones en el Océano Pacífico —Lanai
(Hawai) y Tahití— que son representativas de las condiciones de fondo para el aerosol
marino. El punto de inflexión de la distribución se sitúa en torno a 0.45 µm, con un
Ref~1.72 µm para el modo grueso y Ref~0.14 µm para el modo fino. Al irrumpir la SAL
en la MBL, la distribución de tamaños pasa a ser dominada por el modo grueso
(Vf/Vt~0.12) con un ligero desplazamiento del punto de inflexión hacia radios
geométricos en torno a 0.44 µm. La concentración en volumen del modo grueso
asciende hasta ~0.204 µm3·µm-2, con valores de ~0.026 µm3·µm-2 para el modo fino.
140
Caracterización de aerosoles en Tenerife
de polvo mineral (1.53± 0.05 según Dubovik et al. (2002) y las referencias que
contiene) aunque se aproxima a los valores obtenidos por este autor (1.48±0.05) bajo
condiciones desérticas para la estación de Cabo Verde (muy próxima a las Islas
Canarias).
Tabla 5.6. Promedios del SSA, la parte real e imaginaria del índice de refracción (Ind. Ref.) y parámetro
de asimetría (Asim.) para los canales 440, 675, 870 y 1020 nma junto con el número de observaciones
diarias (N) analizadas tanto para las condiciones de FT y FT-SAL en IZO como para las condiciones de
MBL y MBL-SAL en SCO, en el período entre julio de 2005 y diciembre de 2012.
141
Capítulo 5
entre 1.44 y 1.45. En los casos en los que la SAL se introduce en la MBL, se aprecia un
aumento del parámetro de asimetría (entre 0.73 y 0.75) asociado a la irrupción del
aerosol desértico con unos valores muy próximos a los que se observan para el caso FT-
SAL. La menor absorción de las partículas de polvo queda reflejada en la dependencia
espectral del SSA y del índice de refracción complejo, con valores más altos y más
bajos, respectivamente, en las longitudes de onda entre 675 y 1020 nm. Valores
similares en cuanto a la capacidad de absorción del polvo mineral han sido identificados
previamente (p. ej. Dubovik et al., 2002; Cachorro et al., 2010). Finalmente, el índice de
refracción real asciende ligeramente (entre 1.45 y 1.46) asociado al aumento en la
eficiencia de la dispersión de la radiación solar por parte del polvo mineral.
142
Caracterización de aerosoles en Tenerife
Figura 5.8. Promedio (línea azul) y mediana (línea roja) de los coeficientes de extinción en 532 nm del
lidar micropulsado de SCO entre 2007 y 2012 (salvo 2010) para (a) FT (16820 perfiles disponibles), (b)
FT-SAL (2654 perfiles disponibles), (c) MBL (12432 perfiles disponibles) y (d) MBL-SAL (7140 perfiles
disponibles). El área sombreada en gris claro corresponde al rango entre los percentiles 20 y 80. En los
paneles (a) y (b) se muestra, mediante una línea negra continua, la altura de 2400 m s.n.m., ya que los
criterios utilizados para determinar las condiciones de FT y FT-SAL, respectivamente, se han establecido
a partir de dicha altura.
Los perfiles verticales del coeficiente de extinción (Figura 5.8a) muestran una
atmósfera cada vez más limpia a partir de 2400 m s.n.m. por lo que se va a considerar
que el análisis del mapa resultante de la aplicación del método CWT para las
retrotrayectorias con llegada a esta altura es suficientemente representativo para
determinar el origen de las masa de aire en la FT (Figura 5.9).
Figura 5.9. Mapas resultantes de la aplicación del método CWT, en condiciones de FT, a 2400 m s.n.m.
para (a) el AOD y (b) el EA.
143
Capítulo 5
144
Caracterización de aerosoles en Tenerife
Figura 5.10. Mapas resultantes de la aplicación del método CWT, en condiciones de FT-SAL, a 2400 m
s.n.m. para (a) el AOD y (b) el EA y a 3500 m s.n.m. para (c) el AOD y (d) el EA.
Se aprecia que, respecto de las condiciones de FT (Figura 5.9), las masas de aire
procedentes del Atlántico al norte y noroeste de Canarias han desaparecido y, en
cambio, dominan las procedentes del noroeste de África. Las principales fuentes de
polvo identificadas (con valores asociados de AOD más altos, Figura 5.10a y c) son la
zona norte de Mali y Mauritania (que se corresponde con la región identificada como
zona fuente 4 en la Figura 4.12), Argelia central, el sur de Argelia (zona fuente 1, Figura
4.12) y la vertiente sur del extremo sur de la cordillera del Atlas (Marruecos). También
se identifica la zona de Sahel al oeste de Mauritania como fuente de polvo, aunque
menos frecuente. Estas zonas coinciden, a grandes rasgos, con las presentadas por otros
autores como fuentes de polvo que afectan a la región del SENAR (p. ej. Díaz et al.,
2006; Alonso-Pérez et al., 2012). En cuanto a la depresión del Bodélé, una de las zonas
fuente más importante del mundo, apenas llega a Tenerife polvo mineral procedente de
esta región.
Cabe mencionar la existencia de masas de aire que recorren el Atlántico, al suroeste
de Canarias, y que aportan valores similares de carga de aerosoles (AOD, Figura 5.10a
y c) y de tamaño de partículas (EA, Figura 5.10b y d) que las procedentes directamente
de África. Esta situación se debe a la presencia de polvo mineral emitido desde África
(Sáhara-Sahel) al sur de Canarias, que dichas masa de aire se encuentran a su paso
transportándolo hasta Tenerife. Díaz et al. (2006), aparte de las masas de aire
procedentes de la troposfera oceánica por encima de 2 km que relacionaron con las
condiciones de fondo de la FT (mencionadas en el apartado 5.3.5.1), detectaron que otra
145
Capítulo 5
parte de dichas masas de aire están asociadas a una importante presencia de partículas
de polvo. Estos autores identificaron estos episodios principalmente en primavera
debido a que, en esta época del año, la discontinuidad intertropical (ITD, Inter-Tropical
Discontinuity) se sitúa sobre la zona del Sahel, que se convierte en una importante
fuente de emisión del polvo a la atmósfera. En cambio, en verano la ITD se desplaza
hacia el norte, aumentando en la zona del Sáhara las emisiones de polvo mineral que es
transportado al oeste, alcanzando alturas de hasta 6 km, hacia la región subtropical
(Rodríguez et al., 2014). Al mismo tiempo, la lluvia asociada al monzón se desplaza
hacia el Sahel disminuyendo las emisiones de polvo en esta zona (Knippertz and Todd,
2010; Ashpole and Washington, 2013).
146
Caracterización de aerosoles en Tenerife
Figura 5.11. Mapas resultantes de la aplicación del método CWT, en condiciones de MBL, a nivel de
superficie para (a) el AOD y (b) el EA, a 600 m s.n.m. para (c) el AOD y (d) el EA, y a 1500 m s.n.m.
para (e) el AOD y (f) el EA.
En cuanto al valor del EA asociado a cada masa de aire, cabe recordar que el valor
promedio de esta magnitud para la MBL es ~0.95 (Tabla 5.4). En la Figura 5.11b, d y f
se aprecia una diferencia entre los valores de este parámetro dependiendo de la
procedencia de las masas de aire, con valores más altos del EA al noroeste de Canarias
y valores más bajos al noreste. En este caso, este resultado no se debe a la aplicación de
la función de peso al método CWT ya que también se aprecia en los mapas resultantes
previos a la aplicación de dicha función (no mostrado). La hipótesis barajada a este
respecto es que las masas de aire procedentes del oeste suelen tener un recorrido más
alto que las que llegan a Canarias desde el este, independientemente de la altura del
punto de llegada. Por lo tanto, las masas de aire del este aportan una mayor cantidad de
aerosol marino (grueso) que las que viajan a mayor altura. Esta hipótesis requiere un
análisis en mayor profundidad que será realizado en el futuro.
147
Capítulo 5
Los perfiles verticales del coeficiente de extinción (Figura 5.8d) muestran una serie
de rasgos que se procede a analizar. En niveles bajos se aprecia un aumento
considerable de la extinción con respecto a los perfiles de la MBL, asociado al
transporte de polvo sahariano en capas bajas de la atmósfera, principalmente en invierno
(capítulo 1, apartado 1.5.2). Por debajo de los 1500 m s.n.m. se produce un mínimo
local del coeficiente de extinción, indicando un descenso en el tope de la MBL debido a
que la SAL al desplazarse sobre la MBL hace que disminuya su espesor (Alastuey et al.,
2005). En torno a 2400 m s.n.m. se produce un máximo local, asociado a las intrusiones
en altura, propias del verano, que pueden llegar hasta los 6000 m s.n.m. como se
comentó en el apartado 5.3.5.2. Las amplias diferencias entre los percentiles 20 y 80
indican la gran variabilidad en los coeficientes de extinción debida a la combinación,
por un lado, de transporte de polvo en la MBL y condiciones limpias en capas
superiores y, por otro lado, transporte de polvo por encima de la MBL con condiciones
relativamente limpias en las capas inferiores. Por lo tanto, los mapas CWT que van a ser
analizados corresponden a los puntos de llegada de las retrotrayectorias a nivel de
superficie, a 600, 1500, 2400 y 3500 m s.n.m. (Figura 5.12).
A primera vista llama la atención el mapa a nivel de superficie en el cual, el origen
predominante de las masas de aire es el Atlántico, al noroeste de Canarias, pero está
asociado a valores muy altos del AOD (parámetro que se mide en toda la columna
atmosférica), que es la estación de llegada. Hay tres posibles situaciones que explican
este hecho:
1. Masas de aire a nivel de superficie que tienen su origen en el Atlántico, y son
predominantemente limpias, combinada con masas de aire en altura que
transportan polvo mineral y llevan asociadas valores altos de AOD. Estas son las
situaciones típicas de verano (capítulo 1, apartado 1.5.2) en las cuales el polvo
mineral afecta a la MBL de modo indirecto, por deposición gravitacional.
2. Los denominados “episodios de transición” en los cuales tras una intrusión de
polvo mineral, las masas de aire limpias procedentes del Atlántico empiezan a
reemplazar a las afectadas por el polvo pero aún quedan en la atmósfera
partículas de este tipo de aerosol (Alonso-Pérez et al., 2007).
3. La presencia de polvo mineral al norte de Canarias transportado por las masas de
aire atlánticas que se encuentran en su recorrido hasta Tenerife.
148
Caracterización de aerosoles en Tenerife
Figura 5.12. Mapas resultantes de la aplicación del método CWT, en condiciones de MBL-SAL, a nivel
de superficie para (a) el AOD y (b) el EA, a 600 m s.n.m. para (c) el AOD y (d) el EA, a 1500 m s.n.m.
para (e) el AOD y (f) el EA, a 2400 m s.n.m. para (g) el AOD y (h) el EA, y a 3500 m s.n.m. para (i) el
AOD y (j) el EA.
149
Capítulo 5
150
Caracterización de aerosoles en Tenerife
151
Capítulo 5
152
Conclusiones y líneas futuras
Conclusiones generales
153
Conclusiones y líneas futuras
154
Conclusiones y líneas futuras
155
Conclusiones y líneas futuras
La capa límite marina afectada por polvo mineral se caracteriza por un aumento
del AOD (~0.38), del tamaño de las partículas (EA~0.29) y del contenido de
PWV (~2.04 cm) respecto de las condiciones de fondo. Es la situación
dominante en SCO en julio y agosto (~43% de los días), debido a la mayor
frecuencia de masas saharianas en capas altas de la atmósfera (hasta unos 6 km
s.n.m.) siendo depositado gravitacionalmente. Su principal origen es el norte del
Sáhara Occidental y el oeste de Argelia, es decir, las regiones fuente más
próximas a Canarias. En el período invernal, especialmente en marzo y
noviembre, se produce un ligero repunte en la frecuencia de este transporte
(~18% de los días) que tiene lugar directamente en capas bajas de la atmósfera
(~2 km s.n.m.). En este caso, las principales regiones de origen de las masas de
aire son el sur del Atlas (Marruecos), Argelia central y la zona norte de Mali y
Mauritania.
156
Conclusiones y líneas futuras
157
Conclusiones y líneas futuras
158
Bibliografía
Bibliografía
Alastuey, A., Querol, X., Castillo, S., Escudero, M., Ávila, A., Cuevas, E., Torres, C.,
Romero, P. M., Expósito, F., García O., Díaz, J. P., Dingenen, R. V., Putaud, J. P.:
Characterisation of TSP and PM2.5 at Izana and Sta. Cruz de Tenerife (Canary
Islands, Spain) during a Saharan Dust Episode (July 2002), Atmos. Environ. 39(26),
4715–4728, doi:10.1016/j.atmosenv.2005.04.018, 2005.
Alonso-Pérez, S., Cuevas, E., Querol, X., Viana, M., and Guerra, J. C.: Impact of the
Saharan dust outbreaks on the ambient levels of total suspended particles (TSP) in
the Marine Boundary Layer (MBL) of the Subtropical Eastern North Atlantic
Ocean, Atmos. Environ., 41, 9468-9480, doi:10.1016/j.atmosenv.2007.08.049,
2007.
Alonso-Pérez, S., Cuevas, E., Querol, X.: Objective identification of synoptic
meteorological patterns favouring African dust intrusions into the Marine Boundary
Layer of the Subtropical Eastern North Atlantic Region, Meteorol. Atmos. Phys.,
113, 109–124, doi:10.1007/s00703-011-0150-z, 2011a.
Alonso-Pérez, S., Cuevas, E., Querol, X., Baldasano, J. M., Draxler, R., de Bustos, J. J.:
Trend changes of African airmass intrusions in the marine boundary layer over the
subtropical Eastern North Atlantic region in winter, Tellus 63B, 255-265, doi:
10.1111/j.1600-0889.2010.00524.x, 2011b.
Alonso-Pérez, S., Cuevas, E., Querol, X., Guerra, J. C., and Pérez, C.: African dust
source regions for observed dust outbreaks over the Subtropical Eastern North
Atlantic region above 25°N, J. Arid Environ., 78, 100-109,
doi:10.1016/j.jaridenv.2011.11.013, 2012.
Andrey, J., Cuevas, E., Parrondo, M. C., Alonso-Pérez, S., Redondas, A., and Gil-
Ojeda, M.,:Quantification of ozone reductions within the Saharan air layer through a
13-year climatologic analysis of ozone profiles, Atmos. Environ., 84, 28–34, 2014.
Ångström, A.: On the atmospheric transmission of sun radiation and on dust in the air,
Geogr. Ann., 11, 156–166, 1929.
Ångström, A.: The parameters of atmospheric turbidity, Tellus, 16, 64-76, 1964.
Ashpole, I., and Washington, R.: Intraseasonal variability and atmospheric controls on
daily dust occurrence frequency over the central and western Sahara during the
boreal summer, J. Geophys. Res. Atmos., 118, 12915–12926,
doi:10.1002/2013JD020267, 2013.
Barreto, A., Cuevas, E., Damiri, B., Guirado, C., Berkoff, T., Berjón, A. J., Hernández,
Y., Almansa, F., and Gil, M.: A new method for nocturnal aerosol measurements
159
Bibliografía
160
Bibliografía
Nauels, A., Xia, Y., Bex V., and Midgley P. M. (Eds.), Cambridge University Press,
Cambridge, United Kingdom and New York, NY, USA, 571-657, 2013.
Brooks, N., and Legrand, M.: Dust variability over northern Africa and rainfall in the
Sahel, in: Linking Climate Change to Land Surface Change, McLaren, S. J., and
Kniverton, D. R. (Eds.), Kluwer Academic Publishers, Dordrecht, Netherlands, 1-
25, 2003.
Cachorro, V. E., Romero, P. M., Toledano, C, Cuevas, E, and de Frutos, A. M.: The
fictitious diurnal cycle of aerosol optical depth: A new approach for in situ
calibration and correction of AOD data series, Geophys. Res. Lett., 31, L12106,
doi:10.1029/2004GL019651, 2004.
Cachorro, V. E., Toledano, C., Berjón, A., de Frutos, A. M., Torres, B., Sorribas, M.,
and Laulainen, N. S.: An ‘‘in situ’’ calibration correction procedure (KCICLO)
based on AOD diurnal cycle: Application to AERONET–El Arenosillo (Spain)
AOD data series, J. Geophys. Res., 113, D12205, doi:10.1029/2007JD009673,
2008a.
Cachorro, V. E., Toledano, C, Sorribas, M., Berjón, A., de Frutos, A. M., and
Laulainen, N.: An “in situ” calibration-correction procedure (KCICLO) based on
AOD diurnal cycle: Comparative results between AERONET and reprocessed
(KCICLO method) AOD-alpha data series at El Arenosillo, Spain, J. Geophys. Res.,
113, D02207, doi:10.1029/2007JD009001, 2008b.
Cachorro, V. E., Toledano, C., Antón, M., Berjón, A., de Frutos, A., Vilaplana, J. M.,
Arola, A., and Krotkov, N. A.: Comparison of UV irradiances from Aura/Ozone
Monitoring Instrument (OMI) with Brewer measurements at El Arenosillo (Spain) –
Part 2: Analysis of site aerosol influence, Atmos. Chem. Phys., 10, 11867-11880,
doi:10.5194/acp-10-11867-2010, 2010.
Campbell, J. R., Hlavka, D. L., Welton, E. J., Flynn, C. J., Turner, D. D., Spinhirne, J.
D., Scott, V. S., and Hwang, I. H.: Full-time, eye-safe cloud and aerosol lidar
observation at atmospheric radiation measurement program sites: instrument and
data processing, J. Atmos. Ocean. Tech., 19, 431–442, 2002.
Caquineau, S., Gaudichet, A., Gomes, L., and Legrand, M.: Mineralogy of Saharan dust
transported over northwestern tropical Atlantic Ocean in relation to source regions,
J. Geophys. Res., 107(D15), AAC 4-1–AAC 4-12, doi:10.1029/2000JD000247,
2002.
Carlson, T. N., Prospero, J. M.: The large-scale movement of Saharan air outbreaks
over the northern equatorial Atlantic, Journal of Applied Meteorology, 11, 283–297,
1972.
Córdoba-Jabonero, C., Sorribas, M., Guerrero-Rascado, J. L., Adame, J. A., Hernández,
Y., Lyamani, H., Cachorro, V., Gil, M., Alados-Arboledas, L., Cuevas, E., and de la
Morena, B.: Synergetic monitoring of Saharan dust plumes and potential impact on
161
Bibliografía
surface: a case study of dust transport from Canary Islands to Iberian Peninsula,
Atmos. Chem. Phys., 11, 3067-3091, doi:10.5194/acp-11-3067-2011, 2011.
Cuesta, J., Edouart, D., Mimouni, M., Flamant, P. H., Loth, C., Gibert, F., Marnas, F.,
Bouklila, A., Kharef, M., Ouchene, B., Kadi, M., and Flamant, C.: Multiplatform
observations of the seasonal evolution of the Saharan atmospheric boundary layer in
Tamanrasset, Algeria, in the framework of the African Monsoon Multidisciplinary
Analysis field campaign conducted in 2006, J. Geophys. Res., 113, D00C07,
doi:10.1029/2007JD009417, 2008.
Cuesta, J., Marsham, J. H., Parker, D. J., and Flamant, C.: Dynamical mechanisms
controlling the vertical redistribution of dust and the thermodynamic structure of the
West Saharan atmospheric boundary layer during summer, Atmosph. Sci. Lett., 10,
34–42, doi:10.1002/asl.207, 2009.
Cuesta, J., Lavaysse, C., Flamant, C., Mimouni, M., and Knippertz, P.: Northward
bursts of the West African monsoon leading to rainfall over the Hoggar Massif,
Algeria, Q. J. R. Meteorol. Soc., 136, 174–189, doi:10.1002/qj.439, 2010.
Cuevas, E.: Estudio del comportamiento del ozono troposférico en el Observatorio de
Izaña (Tenerife) y su relación con la dinámica atmosférica, Ph.D. Thesis,
Universidad Complutense de Madrid, 1995.
Cuevas, E., Pérez, C., Baldasano, J. M., Camino, C., Alonso-Pérez , S., and Basart, S.:
Meningitis linked to mineral dust transport in the Sahel. In: Long-term retrospective
analysis of dust and AI with meningitis epidemics data, MACC O-INT WP3., 2nd.
delivery report, 2011.
Cuevas, E., González, Y., Rodríguez, S., Guerra, J. C., Gómez-Peláez, A. J., Alonso-
Pérez, S., Bustos, J., and Milford, C.: Assessment of atmospheric processes driving
ozone variations in the subtropical North Atlantic free troposphere, Atmos. Chem.
Phys., 13, 1973-1998, doi:10.5194/acp-13-1973-2013, 2013.
Cuevas, E., Camino, C., Benedetti, A., Basart, S., Terradellas, E., Baldasano, J. M.,
Morcrette, J. J., Marticorena, B., Goloub, P., Mortier, A., Berjón, A., Hernández, Y.,
Gil-Ojeda, M., and Schulz, M.: The MACC-II 2007–2008 reanalysis: atmospheric
dust evaluation and characterization over northern Africa and the Middle East,
Atmos. Chem. Phys., 15, 3991-4024, doi:10.5194/acp-15-3991-2015, 2015.
D’Almeida, G. A.: A model for Saharan dust transport, Journal of Climate and Applied
Meteorology, 25, 903–916, 1986.
D’Almeida, G. A.: On the variability of desert aerosol radiative characteristics, J.
Geoph. Res., 92, 3017-3026, 1987.D’Almeida, G.A., and Schütz, L.: Number, mass
and volume distributions of mineral aerosols and soils of Sahara, J. Clim. Appl.
Meteorol., 22, 233–243, 1983.
162
Bibliografía
de Longueville, F., Ozer, P., Doumbia, S., and Henry, S., Int. J. Biometeorol., 57 (1), 1-
19, doi: 10.1007/s00484-012-0541-y, 2013.
Díaz, A. M., Díaz, J. P., Expósito, F. J., Hernández-Leal, P. A., Savoie, D., and Querol,
X.: Air masses and aerosols chemical components in the Free Troposphere at the
Subtropical Northeast Atlantic region, Journal of Atmospheric Chemistry, 53 (1),
63-90, doi:10.1007/s10874-006-2644-5, 2006.
Díaz, J. P., Expósito, F. J., Torres, C. J., Herrera, F., Prospero, J. M., Romero, M. C.:
Radiative properties of aerosols in African dust outbreaks using ground-based and
satellite data: applications to radiative forcing, Journal of Geophysical Research,
106(16), 18403-18416, 2001.
Diner, D. J., Beckert, J. C., Reilly, T. H., Bruegge, C. J., Conel, J. E., Kahn, R. A.,
Martonchik, R. A., Ackerman, T. P., Davies, R., Gerstl, S., Gordon, H. R., Muller,
J.-P., Myneni, R., Sellers, P., Pinty, B., and Verstraete, M.: Multi-angle Imaging
Spectro Radiometer (MISR) instrument description and experiment overview, IEEE
T. Geosci. Remote, 36, 1072–1087, 1998.
Draxler, R.R., and Hess, G.D.: An overview of the HYSPLIT_4 modelling system for
trajectories, dispersion, and deposition, Aust. Meteorol. Mag., 47, 295–308, 1998.
Dubovik, O., and King, M. D.: A flexible inversion algorithm for retrieval of aerosol
optical properties from Sun and sky radiance measurements, J. Geophys. Res., 105,
20673-20696, 2000.
Dubovik, O., Smirnov, A., Holben, B. N., King, M. D., Kaufman, Y. J., Eck, T. F., and
Slutsker, I.: Accuracy assessments of aerosol optical properties retrieved from
Aerosol Robotic Network (AERONET) Sun and sky radiance measurements, J.
Geophys. Res., 105, 9791–9806, 2000.
Dubovik, O., Holben, B. N., Eck, T. F., Smirnov, A., Kaufman, Y. J., King, M. D.,
Tanré, D., and Lutsker, I.: Variability of absorption and optical properties of key
aerosols types observed in worldwide locations, J. Atmos. Sci., 59, 590–608, 2002.
Dubovik, O., Sinyuk, A., Lapyonak, T., Holben, B.N., Mishchenko, M., Yang, P., Eck,
T.F., Volten, H., Muñoz, O., Veihelmann, B., van der Zande, W.J., Leon, J-F.,
Sorokin, M., and Slutsker, I.: Application of spheroid models to account for aerosol
particle nonsphericity in remote sensing of desert dust, J. Geophys. Res., 111,
D11208, doi:10.1029/2005JD006619, 2006.
Elvidge, C. D., Baugh, K. E., Kihn, E. A., Kroehl, H. W., and Davis, E. R.: Mapping of
city lights using DMSP Operational Linescan System data, Photogrammetric
Engineering and Remote Sensing, 63(6), 727–734, 1997.
Fernald, F. G., Herman, B. M., and Reagan, J. A.: Determination of aerosol height
distribution by lidar, J. Appl. Meteorol., 11, 482–489, 1972.
163
Bibliografía
164
Bibliografía
González, Y., Rodríguez, S., Guerra, J. C., Trujillo, J. L. and García, R.: Ultrafine
particles pollution in urban coastal air due to ship emissions, Atmospheric
Environment, 45 (28), 4907-4914, doi:10.1016/j.atmosenv.2011.06.002, 2011.
González, Y. and Rodríguez, S.: A comparative study on the ultrafine particle episodes
induced by vehicle exhaust, a crude oil refinery and ship emissions. Atmospheric
Research, 120-121, 43-54, doi:10.1016/j.atmosres.2012.08.001, 2013.
Goudi, A. S., Desert dust and human health disorders, Environ. Int., 63, 101-113, doi:
10.1016/j.envint.2013.10.011, 2014.
Goudie, A. S., and Middleton, N. J.: Saharan dust storms: nature and consequences,
Earth-science Reviews, 56(1), 179-204, doi:10.1016/S0012-8252(01)00067-8, 2001.
Guirado, C., Cuevas, E., Cachorro, V. E., Toledano, C., Alonso-Pérez, S., Bustos, J. J.,
Basart, S., Romero, P. M., Camino, C., Mimouni, M., Zeudmi, L., Goloub, P.,
Baldasano, J. M., and de Frutos, A. M.: Aerosol characterization at the Saharan
AERONET site Tamanrasset, Atmos. Chem. Phys., 14, 11753-11773,
doi:10.5194/acp-14-11753-2014, 2014.
Haustein, K., Pérez, C., Baldasano, J. M., Müller, D., Tesche, M., Schladitz, A.,
Freudenthaler, V., Heese, B., Esselborn, M., Weinzierl, B., Kandler, K., and von
Hoyningen-Huene, W.: Regional dust model performance during SAMUM 2006,
Geophys. Res. Lett., 36, L03812, doi:10.1029/2008GL036463, 2009.
Haustein, K., Pérez, C., Baldasano, J. M., Jorba, O., Basart, S., Miller, R. L., Janjic, Z.,
Black, T., Nickovic, S., Todd, M. C., Washington, R., Müller, D., Tesche, M.,
Weinzierl, B., Esselborn, M., and Schladitz, A.: Atmospheric dust modeling from
meso to global scales with the online NMMB/BSC-Dust model – Part 2:
Experimental campaigns in Northern Africa, Atmos. Chem. Phys., 12, 2933-2958,
doi:10.5194/acp-12-2933-2012, 2012.
Heffter, J. L.: Air Resources Laboratories Atmospheric Transport and Dispersion Model
(ARL-ATAD), Air Resources Laboratories, Silver Spring, Maryland, 1980.
Heinold, B., Knippertz, P., Marsham, J. H., Fiedler, S., Dixon, N. S., Schepanski, K.,
Laurent, B.,and Tegen, I.: The role of deep convection and nocturnal low-level jets
for dust emission in summertime West Africa: Estimates from convection permitting
simulations, J. Geophys. Res. Atmos., 118, 4385–4400, doi:10.1002/jgrd.50402,
2013.
Heintzenberg, J., Raes, F., and Schwartz, S.: Tropospheric aerosols, in: Atmospheric
Chemistry in a Changing World, Brasseur, G., Prinn, R., and Pszenny, A. P. (Eds.),
Springer, Berlin, Heidelberg, 125-156, 2003.
Holben, B. N., Eck, T. F., Slutsker, I., Tanré, D., Buis, J. P., Setzer, A., Vermote, E.,
Reagan, J. A., Kaufman, Y. J., Nakajima, T., Lavenu, F., Jankowiak, I., and
165
Bibliografía
Smirnov, A.: AERONET- A federated instrument network and data archive for
aerosol characterization, Rem. Sens. Env., 66(1), 1–16, 1998.
Holben, B. N., Tanre, D., Smirnov, A., Eck, T. F., Slutsker, I., Abuhassan, N.,
Newcomb, W.W., Schafer, J., Chatenet, B., Lavenue, F., Kaufman, Y. J., Castle, J.
V., Setzer, A., Markham, B., Clark, D., Frouin, R., Halthore, R., Karnieli, A.,
O’Neill, N. T., Pietras, C., Pinker, R. T., Voss, K., and Zibordi, G.: An emerging
ground-based aerosol climatology: Aerosol Optical Depth from AERONET,
Geophys. Res., 106(12), 12067–12097, 2001.
Hsu, N. C., Tsay, S.-C., King, M. D., and Herman, J. R.: Aerosol properties over bright-
reflecting source regions, IEEE T. Geosci. Remote Sens., 42, 557–569, 2004.
Hsu, N. C., Tsay, S.-C., King, M. D., and Herman, J. R.: Deep Blue retrievals of Asian
aerosol properties during ACE-Asia, IEEE T. Geosci. Remote Sens., 44, 3180–
3195, 2006.
Huneeus, N., Schulz, M., Balkanski, Y., Griesfeller, J., Prospero, J., Kinne, S., Bauer,
S., Boucher, O., Chin, M., Dentener, F., Diehl, T., Easter, R., Fillmore, D., Ghan, S.,
Ginoux, P., Grini, A., Horowitz, L., Koch, D., Krol, M. C., Landing, W., Liu, X.,
Mahowald, N., Miller, R., Morcrette, J.-J., Myhre, G., Penner, J., Perlwitz, J., Stier,
P., Takemura, T., and Zender, C. S.: Global dust model intercomparison in
AeroCom phase I, Atmos. Chem. Phys., 11, 7781-7816, doi:10.5194/acp-11-7781-
2011, 2011.
Hunt, W. H., Winker, D. M., Vaughan, M. A., Powell, K. A., Lucker, P. L., and
Weimer, C.: CALIPSO lidar description and performance assessment, J. Atmos.
Oceanic Technol., 26, 1214–1228, doi:10.1175/2009JTECHA1223.1, 2009.
Jakob, C., and Tselioudis, G.: Objective identification of cloud regimes in the Tropical
Western Pacific, Geophys. Res. Lett., 30(21), 2082, doi:10.1029/2003GL018367,
2003.
Jickells, T. D., An, Z. S., Andersen, K. K., Baker, A. R., Bergametti, G., Brooks, N.,
Cao, J. J., Boyd, P. W., Duce, R. A., Hunter, K. A., Kawahata, H., Kubilay, N., la
Roche, J. S. L. P., Mahowald, N., Prospero, J. M., Ridgwell, A. J., Tegen, I., and
Torres, R.: Global iron connections between desert dust, ocean biogeochemistry,
and climate, Science, 308, 67–71, 2005.
Jung, J., Lee, H., Kim, Y. J., Liu, X., Zhang, Y., Gu, J., and Fan, S.: Aerosol chemistry
and the effect of aerosol water content on visibility impairment and radiative forcing
in Guangzhou during the 2006 Pearl River Delta campaign, Journal of
Environmental Management, 90, 3231-3244, 2009.
Kahn, R. A., Gaitley, B. J., Martonchik, J. V., Diner, D. J., Crean, K. A., and Holben, B.
N.: Multiangle Imaging Spectroradiometer (MISR) global aerosol optical depth
validation based on 2 years of coincident Aerosol Robotic Network (AERONET)
observations, J. Geophys. Res., 110, D10S04, doi:10.1029/2004JD004706, 2005.
166
Bibliografía
Kahn, R. A., Gaitley, B. J., Garay, M. J., Diner, D. J., Eck, T. F., Smirnov, A., and
Holben, B. N.: Multiangle Imaging SpectroRadiometer global aerosol product
assessment by comparison with the Aerosol Robotic Network, J. Geophys. Res.,
115, D23209, doi:10.1029/2010JD014601, 2010.
Kaufman, Y. J., Tanré, D., Remer, L. A., Vermote, E. F., Chu, A., and Holben, B. N.:
Operational remote sensing of tropospheric aerosol over land from EOS moderate
resolution imaging spectroradiometer, J. Geophys. Res., 102, 17–51, 1997.
Kim, D., Chin, M., Yu, H., Eck, T. F., Sinyuk, A., Smirnov, A., and Holben, B. N.: Dust
optical properties over North Africa and Arabian Peninsula derived from the
AERONET dataset, Atmos. Chem. Phys., 11, 10733-10741, doi:10.5194/acp-11-
10733-2011, 2011.
Klett, J. D.: Stable analytic inversion solution for processing lidar returns, Appl. Optics,
20, 211–220, 1981.
Klett, J. D.: Lidar inversion with variable backscatter/extinction ratios, Appl. Optics, 24,
1638–1643, 1985.
Knippertz, P., and Todd, M. C.: The central west Saharan dust hot spot and its relation
to African easterly waves and extratropical disturbances, J. Geophys. Res., 115,
D12117, doi:10.1029/2009JD012819, 2010.
Knippertz, P., and Todd, M. C.: Mineral dust aerosols over the Sahara: Meteorological
controls on emission and transport and implications for modeling. Rev. Geophys.
50, RG1007. doi:10.1029/2011RG000362, 2012.
Kunz, G. and de Leeuw, G.: Inversion of lidar signals with the slope method, Appl.
Optics, 32, 3249–3256, 1993.
Li, J., Carlson, B. E., Dubovik, O., and Lacis, A. A.: Recent trends in aerosol optical
properties derived from AERONET measurements, Atmos. Chem. Phys., 14, 12271-
12289, doi:10.5194/acp-14-12271-2014, 2014.
Mahowald, N. M., Baker, A. R., Bergametti, G., Brooks, N., Duce, R. A., Jickells, T.
D., Kubilay, N., Prospero, J. M., and Tegen, I.: The atmospheric global dust cycle
and iron inputs to the ocean, Global Biogeochem. Cy., 19, GB4025,
doi:10.1029/2004GB002402, 2005.
Mahowald, N. M., Kloster, S., Engelstaedter, S., Moore, J. K., Mukhopadhyay, S., Mc-
Connell, J. R., Albani, S., Doney, S. C., Bhattacharya, A., Curran, M. A. J., Flanner,
M. G., Hoffman, F. M., Lawrence, D. M., Lindsay, K., Mayewski, P. A., Neff, J.,
Rothenberg, D., Thomas, E., Thornton, P. E., and Zender, C. S.: Observed 20th
century desert dust variability: impact on climate and biogeochemistry, Atmos.
Chem. Phys., 10, 10875–10893, doi:10.5194/acp-10-10875-2010, 2010.
Mallet, M., Dubovik, O., Nabat, P., Dulac, F., Kahn, R., Sciare, J., Paronis, D., and
Léon, J. F.: Absorption properties of Mediterranean aerosols obtained from multi-
167
Bibliografía
year ground-based remote sensing observations, Atmos. Chem. Phys., 13, 9195-
9210, doi:10.5194/acp-13-9195-2013, 2013.
Marsham, J. H., Parker, D. J., Grams, C. M., Taylor, C. M., and Haywood, J. M.: Uplift
of Saharan dust south of the intertropical discontinuity, J. Geophys. Res., 113,
D21102, doi:10.1029/2008JD009844, 2008.
Marsham, J. H., Knippertz, P., Dixon, N. S., Parker, D. J., and Lister, G. M. S.: The
importance of the representation of deep convection for modeled dust-generating
winds over West Africa during summer, Geophys. Res. Lett., 38, L16803,
doi:10.1029/2011GL048368, 2011.
Marsham, J. H., Hobby, M., Allen, C. J. T., Banks, J. R., Bart, M., Brooks, B. J.,
Cavazos-Guerra, C., Engelstaedter, S., Gascoyne, M., Lima, A. R., Martins, J. V.,
McQuaid, J. B., O’Leary, A., Ouchene, B., Ouladichir, A., Parker, D. J., Saci, A.,
Salah-Ferroudj, M., Todd, M. C., and Washington, R.: Meteorology and dust in the
central Sahara: Observations from Fennec supersite-1 during the June 2011
Intensive Observation Period, J. Geophys. Res. Atmos., 118, 4069–4089,
doi:10.1002/jgrd.50211, 2013.
Miloshevich, L. M., Vömel, H., Whilteman, D. N., and Leblanc, T.: Accuracy
assessment and correction of Vaisala RS92 radiosonde water vapor measurements,
J. Geophys. Res., 114, D11305, doi:10.1029/2008JD011565, 2009.
Morcrette, J.-J., Boucher, O., Jones, L., Salmond, D., Bechtold, P., Beljaars, A.,
Benedetti, A., Bonet, A., Kaiser, J. W., Razinger, M., Schulz, M., Serrar, S.,
Simmons, A. J., Sofiev, M., Suttie, M., Tompkins, A. M., and Untch, A.: Aerosol
analysis and forecast in the ECMWF integrated forecast system. Part I: Forward
modelling, J. Geophys. Res., 114, D06206, doi:10.1029/2008JD011235, 2009.
Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang, J., Koch, D.,
Lamarque, J.-F., Lee, D. Mendoza, B. Nakajima, T. Robock, A., Stephens, G.,
Takemura, T., and Zhang, H.: Anthropogenic and Natural Radiative Forcing. In:
Climate Change 2013: The Physical Science Basis. Contribution of Working Group
I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change,
Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M. Allen, S. K., Boschung, J.,
Nauels, A., Xia, Y., Bex V., and Midgley P. M. (Eds.), Cambridge University Press,
Cambridge, United Kingdom and New York, NY, USA, 659-740,2013
Naseema Beegum, S., Krishna Moorthy, K., Gogoi, Mukunda M., Suresh Babu, S., and
Pandey, S. K.: Multi-year investigations of aerosols from an island station, Port
Blair, in the Bay of Bengal: climatology and source impacts, Ann. Geophys., 30,
1113-1127, doi:10.5194/angeo-30-1113-2012, 2012.
Nickovic, S., Kallos, G., Papadopoulos, A., and Kakaliagou, O.: A model for prediction
of desert dust cycle in the atmosphere, J. Geophys. Res., 106, 18113–18130,
doi:10.1029/2000JD900794, 2001.
168
Bibliografía
169
Bibliografía
Prospero, J. M., Ginoux, P., Torres, O., Nicholson, S. E., and Gill, T. E.: Environmental
characterization of global sources of atmospheric soil dust identified with the
nimbus 7 total ozone mapping spectrometer (TOMS) absorbing aerosol product,
Rev. Geophys., 40(1), 1002, doi:10.1029/2000RG000095, 2002.
Prospero, J. M., Landing, W. M., and Schulz, M.: African dust deposition to Florida:
Temporal and spatial variability and comparisons to models, J. Geophys. Res., 115,
D13304, doi:10.1029/2009JD012773, 2010.
Raes, F., Van Dingenen, R., Cuevas, E., Van Velthoven, P. F. J., and Prospero J. M.:
Observations of aerosols in the free troposphere and marine boundary layer of the
subtropical Northeast Atlantic: Discussion of processes determining their size
distribution, J. Geophys. Res., 102(D17), 21315-21328, doi:10.1029/97JD01122,
1997.
Raes, F., Bates, T., McGovern, F., and Van Liedekerke, M.: The 2nd Aerosol
Characterization Experiment (ACE-2): general overview and main results, Tellus B,
52, 111–125, doi: 10.1034/j.1600-0889.2000.00124.x, 2000.
Reddy, P. J.: Instructions for J-Series Handheld Sunphotometer, NOAA, Boulder, 1986.
Redelsperger, J.-L., Thorncroft, C., Diedhiou, A., Lebel, T., Parker, D., and Polcher, J.:
African monsoon multidisciplinary analysis - An international research project and
field campaign, Bull. Amer. Meteor. Soc., 87(12), 1739-1746,
doi:http://dx.doi.org/10.1175/BAMS-87-12-1739, 2006.
Remer, L. A., Kaufman, Y. J., Tanr’e, D., Mattoo, S., Chu, D. A., Martins, J. V., Li, R.
R., Ichoku, C., Levy, R. C., Kleidman, R. G., Eck, T. F., Vermote, E., and Holben,
B. N.: The MODIS aerosol algorithm, products, and validation, J. Atmos. Sci., 62,
947–973, 2005.
Roberts, A. J.: Anomalously heavy rainfall and dust in the arid Sahara and northern
Sahel, In: Convective Episodes near the Intertropical Discontinuity in Summertime
West Africa: Representation in Models and Implications for Dust Uplift, PhD thesis,
University of Leeds, Leeds, UK, 2014.
Roberts, A. J, Knippertz, P., and Marsham, J. H.: The Formation of Convectively
Generated Dusty Episodes in the Sahara during Summer, DUST-2014, International
Conference on Atmospheric Dust, Castellaneta Marina, Italy, June 1-6, 2014.
Rodríguez, S.: Comparación de las variaciones de ozono superficial asociadas a
procesos de transporte sobre y bajo la inversión de temperatura subtropical en
Tenerife, Ph.D. Thesis, Universidad de La Laguna, Tenerife, Canarias, 1999.
Rodríguez, S., and Guerra, J. C.: Monitoring of ozone in a marine environment in
Tenerife (Canary Islands), Atmos. Environ., 135, 1829–1841, doi:10.1016/S1352-
2310(00)00550-1, 2001.
170
Bibliografía
Rodríguez, S., Alastuey, A., Viana, M. and Querol, X.: Estudios de contaminación por
material particulado en Canarias durante los años 2007 a 2010, informe elaborado
por el Instituto de Diagnóstico Ambiental y Estudios del Agua (IDÆA - CSIC), el
Centro de Investigación Atmosférica de Izaña (Agencia Estatal de Meteorología -
AEMET) y la Universidad de Huelva para la Consejería del Medio Ambiente y
Ordenación Territorial del Gobierno de Canarias, 2010.
Rodríguez, S., Alastuey, A., Alonso-Pérez, S., Querol, X., Cuevas, E., Abreu-Afonso,
J., Viana, M., Pérez, N., Pandolfi, M., and de la Rosa, J.: Transport of desert dust
mixed with North African industrial pollutants in the subtropical Saharan Air Layer,
Atmos. Chem. Phys., 11, 6663-6685, doi:10.5194/acp-11-6663-2011, 2011.
Rodríguez, S., Cuevas, E., Prospero, J. M., Alastuey, A., Querol, X., López-Solano, J.,
García, M. I., and Alonso-Pérez, S.: Modulation of Saharan dust export by the North
African dipole, Atmos. Chem. Phys. Discuss., 14, 26689-26719, doi:10.5194/acpd-
14-26689-2014, 2014.
Romero, P. M., and Cuevas, E.: Diurnal variation of the aerosol optical depth: artifact or
reality?, in: Proceeding of 3ª Asamblea Hispano Portuguesa de Geofísica y
Geodesia, Valencia, Spain, 4-8 February 2002, 2(S13), 1252-1256, 2002.
Sasano, Y. and Nakane, H.: Significance of the extinction/backscatter ratio and the
boundary value term in the solution for the two-component lidar equation, Appl.
Optics, 23, 11–13, 1984.
Schepanski, K., Tegen, I., Todd, M. C., Heinold, B., Bönisch, G., Laurent, B., and
Macke, A.: Meteorological processes forcing Saharan dust emission inferred from
MSG-SEVIRI observations of subdaily dust source activation and numerical
models, J. Geophys. Res., 114, D10201, doi:10.1029/2008JD010325, 2009.
Schmetz, J., Pili, P., Tjemkes, S., Just, D., Kerkmann, J., Rota, S., and Ratier, A.: An
introduction to Meteosat Second Generation (MSG), Bull. Am. Meteorol. Soc., 83,
977–992, 2002.
Schneider, M., Romero, P. M., Hase, F., Blumenstock, T., Cuevas, E., and Ramos, R.:
Continuous quality assessment of atmospheric water vapour measurement
techniques: FTIR, Cimel, MFRSR, GPS, and Vaisala RS92, Atmos. Meas. Tech., 3,
323-338, doi:10.5194/amt-3-323-2010, 2010.
Schulz, M., Prospero, J. M., Baker, A. R., Dentener, F., Ickes, L., Liss, P. S.,
Mahowald, N. M., Nickovic, S., Pérez García-Pando, C., Rodríguez, S., Sarin, M.,
Tegen, I., and Duce, R. A, Atmospheric Transport and Deposition of Mineral Dust
to the Ocean: Implications for Research Needs, Environmental Science &
Technology, 46 (19), 10390-10404, doi:10.1021/es300073u, 2012.
Schuster, G. L., Dubovik, O., and Holben, B. N.: Angstrom exponent and bimodal
aerosol size distributions, J. Geophys. Res., 111, D07207,
doi:10.1029/2005JD006328, 2006.
171
Bibliografía
Schuster, G. L., Vaughan, M., MacDonnell, D., Su, W., Winker, D., Dubovik, O.,
Lapyonok, T., and Trepte, C.: Comparison of CALIPSO aerosol optical depth
retrievals to AERONET measurements, and a climatology for the lidar ratio of dust,
Atmos. Chem. Phys., 12, 7431–7452, doi:10.5194/acp-12-7431-2012, 2012.
Smirnov, A., Holben, B. N., Kaufman, Y. J., Dubovik, O., Eck, T. F., Slutsker, I.,
Pietras, C., Halthore, R.: Optical properties of atmospheric aerosol in maritime
environments, J. Atmos. Sci., 59, 501–523, 2002.
Sorribas, M: Caracterización y calibración de instrumentos de medida de radiación
solar, Diploma de Estudios Avanzados, Universidad de Valladolid, 2003.
Stefanski, R. and Sivakumar, M. V. K.: Impacts of Sand and Dust Storms on
Agriculture and Potential Agricultural Applications of a SDSWS, IOP C. Ser. Earth
Env., 7, 012016, doi:10.1088/1755-1307/7/1/012016, 2009.
Su, L., and Toon, O. B.: Saharan and Asian dust: similarities and differences determined
by CALIPSO, AERONET, and a coupled climate-aerosol microphysical model,
Atmos. Chem. Phys., 11, 3263-3280, doi:10.5194/acp-11-3263-2011, 2011.
Tanré, D., Kaufman, Y. J., Herman, M., and Mattoo, S.: Remote sensing of aerosol
properties over oceans using the MODIS/EOS spectral radiances, J. Geophys. Res.,
102, 16971–16988, doi:10.1029/96JD03437, 1997.
Tegen, I., Schepanski, K., and Heinold, B.: Comparing two years of Saharan dust source
activation obtained by regional modelling and satellite observations, Atmos. Chem.
Phys., 13, 2381–2390, doi:10.5194/acp-13-2381-2013, 2013.
Tesche, M., Wandinger, U., Ansmann, A., Althausen, D., Müller, D., and Omar A. H.:
Ground-based validation of CALIPSO observations of dust and smoke in the Cape
Verde region, J. Geophys. Res. Atmos., 118, 2889–2902, doi:10.1002/jgrd.50248,
2013.
Toledano, C., Cachorro, V. E., Berjon, A., de Frutos, A. M., Sorribas, M., de la Morena,
B. A., and Goloub, P.: Aerosol optical depth and Ångström exponent climatology at
El Arenosillo AERONET site (Huelva, Spain), Q. J. R. Meteorol. Soc., 133, 795–
807, doi:10.1002/qj.54, 2007.
Toledano, C., Wiegner, M., Garhammer, M., Seefeldner, M., Gasteiger, J., Müller, D.,
and Koepke, P.: Spectral aerosol optical depth characterization of desert dust during
SAMUM 2006, Tellus B, 61(1), 216-228, doi:10.1111/j.1600-0889.2008.00382.x,
2009.
Torres, C., Cuevas, E., Guerra, J. C., Carreño, V.: Characterización de las masas de aire
en la región subtropical, in: Proceedings of the V Simposio Nacional de Predicción,
Instituto Nacional de Meteorología, Madrid, ISBN:84-8320-192-5, 2001.
Viana, M., Querol, X., Alastuey, A., Cuevas, E., and Rodríguez, S.: Influence of
African dust on the levels of atmospheric particulates in the Canary Islands air
172
Bibliografía
173
Anexo 1
Sección Especial: 37th AMASOM / Special Section: 37th AMASOM
Radiation and Atmospheric Components
Optical calibration facility at the Izaña Atmospheric Research Center
Laboratorio de calibración óptica del Centro de Investigación Atmosférica
de Izaña
C. Guirado(1,2,*), R. Ramos(2), Á. de Frutos(1,S), A. Berjón(1,3,S), A. Redondas(2), C. López(4),
V. Cachorro(1,S), E. Cuevas(2), R. González(1), S. González(1), M. Hernández(4)
1. Atmospheric Optics Group, Valladolid University (GOA‐UVA), Spain.
2. Izaña Atmospheric Research Center, Meteorological State Agency of Spain, AEMET, Spain.
3. Laboratoire d’Optique Atmosphérique, Université Lille, France.
4. Sieltec Canarias SL, Spain.
Email: [email protected]
(*) S: miembro de SEDOPTICA / SEDOPTICA member
Recibido / Received: 30/01/2011. Aceptado / Accepted: 30/08/2011.
ABSTRACT:
During the last years a new optical calibration facility has been developed and deployed at the Izaña
Observatory for the calibration and characterization of the radiation measuring instruments within
research activities. These activities require a traceable Quality Assurance & Quality Control system.
This new facility allows the absolute, spectral and cosine response calibration. At present there are
six set‐ups running at the laboratory for radiance calibration, angular response determination,
spectral response characterization, slit function determination and absolute irradiance calibration
(vertical and horizontal set‐ups). Each method and procedure has required the development of the
corresponding protocol. These systems, as well as some examples of instrument calibrations, are
presented in this work.
Keywords: Calibration Facility, Quality Assurance, Quality Control, Radiation, Izaña.
RESUMEN:
Durante los últimos años ha sido implementado un nuevo laboratorio de calibración óptica en el
Observatorio de Izaña para la caracterización de los instrumentos de medida de radiación, en el
marco de actividades de investigación que requieren un estricto control de calidad de las medidas
realizadas. Esta nueva instalación dispone de diversos dispositivos para la calibración absoluta,
espectral y de respuesta coseno de los instrumentos radiométricos mencionados anteriormente. En
este momento hay seis sistemas funcionando en el laboratorio que permiten la calibración en
radiancia, el cálculo de la respuesta angular, la caracterización de la respuesta espectral, la
determinación de la función rendija y, dos de ellos, la calibración absoluta en irradiancia (tanto en
posición horizontal como vertical). Cada procedimiento o método de calibración ha requerido la
realización del correspondiente protocolo. En este artículo se realiza una descripción de cada
dispositivo y se muestran algunos ejemplos de calibración instrumental.
Palabras clave: Laboratorio de Calibración, Control de Calidad, Radiación, Izaña.
REFERENCIAS Y ENLACES / REFERENCES AND LINKS
[1]. http://www.pmodwrc.ch/euvc/euvc.php?topic=calibration_services
[2.] H. Slaper, P. den Outer, Analysis of spectral wavelength shifts, spectral anomalies, spectral resolution
and calculation of weighted UV‐doses by means of the software package SHICrivm, User Guide version
3.075, July 2002
2. Calibration set‐ups
2.a. Absolute irradiance calibration
There are two different systems for the absolute
irradiance calibration.
One of them has been designed for the
calibration of the instruments in a horizontal Fig. 2. Vertical absolute irradiance calibration set‐up.
alignment of the system is slightly more difficult. possible but only ±90º of rotation is used for the
Therefore the horizontal set‐up is used to calibration while the instrument is illuminated
calibrate the smaller instruments (pyranometers by the uniform and parallel light beam of the
and multifilter radiometers) while the larger lamp.
ones (Brewer and Bentham spectrometers) are
calibrated using the vertical set‐up.
The basis of the absolute irradiance scale
consists on a set of DXW‐type 1000 W lamps (for
vertical set‐up) and 1000 W FEL lamps (for
horizontal set‐up) traceable to the primary
irradiance standard of the Physikalisch‐
Technische Bundesanstalt (PTB). The accuracy
of the current intensity of the lamp during the
calibration and the distance between the optical
entrance of the instrument and the lamp are
tightly controlled.
2.b. Absolute radiance calibration
This system has been developed for the
instruments which perform sky radiance
measurements as the Cimel sun photometers of
the AERONET‐PHOTONS‐RIMA networks (Fig. Fig. 3. Radiance calibration set‐up with a cimel sun‐
photometer.
3).
The light source for the radiance calibration
comes from an integrating sphere wich provides
a spatially integrated radiant flux. Its calibration
certificate is traceable to the National Institute of
Standards and Technology (NIST) and three
times a year a recalibration of the sphere is
accomplished by comparison to the travel
Master Cimel sun photometer of the National
Aeronautics and Space Administration (NASA).
The accuracy of the calibration is dependent on
the stability of the radiance emission and the
alignment between the sphere and the detector
to be calibrated.
2.c. Angular response calibration
The irradiance measurements require the
detector to weight incoming radiation with the
cosine of the incoming angle relative to normal
incidence [1]. This set‐up (Fig. 4) is used to
quantify the deviations from an ideal cosine
response of radiometer and spectrophotometer
measurements.
The relative angular response function is
measured rotating the mechanical arm where
the seasoned DXW‐type 1000 W lamp is located.
The rotation is controlled by a stepper motor Fig. 4. Angular response function set‐up with a YES UVB‐1
with 0.01º of precision. Rotation over 360º is pyranometer.
Fig. 6. Slit function determination set‐up. Fig. 7. Radiance calibration of the Cimel#380.
TABLE I
Aureole radiance calibration factors
Time 1020 nm 870 nm 670 nm 440 nm
9:48:31 0.019529 0.014046 0.019755 0.033614
9:58:14 0.019497 0.014035 0.019741 0.033593
14:39:04 0.019396 0.014002 0.019711 0.033542
14:51:16 0.019361 0.01399 0.019695 0.033506
Mean 0.019446 0.014018 0.019726 0.033564
St Dev(%) 0.411559 0.189308 0.138867 0.145846
TABLE II
Sky radiance calibration factors
Time 1020 nm 870 nm 670 nm 440 nm
9:48:31 0.004445 0.002957 0.004017 0.008079
9:58:14 0.004437 0.002955 0.004014 0.00807
14:39:04 0.004407 0.002947 0.004002 0.008046
14:51:16 0.004398 0.002946 0.003998 0.008038
Mean 0.004422 0.002951 0.004008 0.008058
St Dev(%) 0.514844 0.188404 0.229026 0.240712
Finally the ultraviolet response was
calculated (Fig. 9) to apply it to the UV radiation
measurements of this instrument.
8.001
8
7.999
13:30 14:00 14:30 15:00 15:30 16:00
Fig. 10. Several steps of the irradiance calibration of the
Voltage (V)
112.8
Brewer#201.
112.75
112.7
112.65
13:30 14:00 14:30 15:00
March 10th, 2010
15:30 16:00 4. Conclusions
Fig. 8. Intensity and voltage applied to the lamp (vertical A new optical calibration facility is available at
black lines show the interval time of calibration). the Izaña Observatory. At present, six set‐ups
have been developed for the absolute, spectral
Brewer#201: Irradiance calibration and cosine response calibration of the radiation
9000 measuring instruments installed at IZO. This
8750 facility allows a better Quality Assurance &
Counts / ( W / m2 )
8500
Quality Control of the radiation measurements.
8250
8000 The systems are used for periodical calibrations
7750 of radiometers, spectroradiometers and
7500 sunphotometers involved in research activities
7250
7000
within GAW, NDACC, BSRN, AERONET and
2865 2975 3085 3195 3305 3415 3525 3635 RBCC‐E programmes.
Wavelength (A)
Fig. 9. UV response of the Brewer#201. The “break” at about
3500 Å is due to a change in the position of the slitmask to
measure the larger wavelengths.
Acknowledgements
This work was developed within the Specific
Agreement of Collaboration between the
University of Valladolid and the CIAI‐AEMET
“Establish methodologies and quality assurance
systems for programs of photometry,
radiometry, atmospheric ozone and aerosols
within the atmospheric monitoring program of
the World Meteorological Organization.
Financial supports from the Spanish MICIIN (ref.
CGL2008‐05939‐CO3‐00/CLI and CGL 2009‐
09740) and from the GR‐220 Project of the Junta
de Castilla y León are gratefully acknowledged.
Barcelona, Spain
4 Universidad Europea de Canarias, Laureate International Universities, La Orotava, Spain
5 Earth Sciences Department, Barcelona Supercomputing Center-Centro Nacional de Supercomputación,
Received: 28 March 2014 – Published in Atmos. Chem. Phys. Discuss.: 24 June 2014
Revised: 1 September 2014 – Accepted: 25 September 2014 – Published: 11 November 2014
Abstract. More than 2 years of columnar atmospheric tify potential sources of air masses arriving at Tamanrasset
aerosol measurements (2006–2009) at the Tamanrasset site at several heights for each season. Microphysical and optical
(22.79◦ N, 5.53◦ E, 1377 m a.s.l.), in the heart of the Sahara, properties and precipitable water vapour were also investi-
are analysed. Aerosol Robotic Network (AERONET) level gated.
2.0 data were used. The KCICLO (K is the name of a con-
stant and ciclo means cycle in Spanish) method was applied
to a part of the level 1.5 data series to improve the quality
of the results. The annual variability of aerosol optical depth 1 Introduction
(AOD) and Ångström exponent (AE) has been found to be
strongly linked to the convective boundary layer (CBL) ther- The regional characterization of mineral dust, particularly
modynamic features. The dry-cool season (autumn and win- close to source areas, has become a valuable tool for re-
ter) is characterized by a shallow CBL and very low mean searchers from different fields. It will lead to reducing some
turbidity (AOD ∼ 0.09 at 440 nm, AE ∼ 0.62). The wet-hot existing uncertainties about direct radiative forcing by atmo-
season (spring and summer) is dominated by high turbidity spheric aerosols (Forster et al., 2007), and to a better un-
of coarse dust particles (AE ∼ 0.28, AOD ∼ 0.39 at 440 nm) derstanding about the potential impact of aerosol on human
and a deep CBL. The aerosol-type characterization shows health and air quality (e.g. De Longueville et al., 2010; Perez
desert mineral dust as the prevailing aerosol. Both pure Saha- et al., 2012).
ran dust and very clear sky conditions are observed depend- The Sahara and its margins are the largest and most contin-
ing on the season. However, several case studies indicate an uous dust sources in the world. Several satellite and ground-
anthropogenic fine mode contribution from the industrial ar- based observation analyses have led to identifying the base of
eas in Libya and Algeria. The concentration weighted trajec- the Ahaggar and Tibesti mountains and the Bodélé Depres-
tory (CWT) source apportionment method was used to iden- sion as the major sources in this area (Goudi and Middle-
ton, 2001; Prospero et al., 2002; Ginoux et al., 2012). During
recent years, several field campaigns in different locations nent (AE) time series. Specific characterizations have been
focused on the analysis of Saharan dust features (Todd et made for the first time: annual evolution and seasonal fea-
al., 2013, and references therein). In particular, Tamanras- tures of precipitable water vapour (PWV), fine mode frac-
set (main city in the Hoggar, also known as Ahaggar Moun- tion (FMF), and aerosol microphysics, as well as an identi-
tains, in Algeria) hosted a specific soil and aerosol sam- fication of potential source regions. The paper is structured
pling analysis at the beginning of the 1980s (D’Almeida as follows: measurement site, data sets and tools used are de-
and Schütz, 1983), the African turbidity monitoring network scribed in Sect. 2. In Sect. 3.1 the main aerosol and PWV
(1980–1984) for climate modelling purposes (D’Almeida, seasonal features are analysed, an aerosol-type classifica-
1986, 1987), and the more recent African Monsoon Multi- tion is performed and microphysical and optical properties
disciplinary Analysis (AMMA) campaign (Redelsperger et are discussed. In Sect. 3.2 the concentration weighted tra-
al., 2006). jectory method is used to identify potential source regions.
During AMMA’s intensive observation periods in 2006, In Sect. 3.3 the transport of anthropogenic fine aerosols to
Tamanrasset was a fully equipped ground-based station for Tamanrasset is discussed. In Sect. 4 the main concluding
aerosol and radiation measurements. This campaign has pro- points are provided.
vided comprehensive analysis of several features at Taman-
rasset and the Hoggar Mountains (e.g. Flamant et al., 2007;
Bou Karam et al., 2008; Cuesta et al., 2008, 2009, 2010). In 2 Methodology
addition, aerosol observations carried out at Tamanrasset in
2.1 Unique characteristics of the Tamanrasset site
2006 have been part of selected aerosol data sets used for sev-
eral model validations (e.g. Toledano et al., 2009; Haustein On 30 September 2006, a sun photometer was installed on
et al., 2009, 2012; Su and Toon, 2011). In spite of these stud- the roof of the main building of the Regional Meteorological
ies, limited aerosol observations, mainly confined to shorter Center (Direction Météo Régional Sud, Office National de
period campaigns, are available for this area which is strate- la Météorologie, Algeria) at Tamanrasset (22.79◦ N, 5.53◦ E,
gically located in the heart of the Sahara. 1377 m a.s.l.) in southern Algeria. Tamanrasset is free from
Consequently, Tamanrasset was considered to be a key industrial activities and is representative of pure desert dust
place to initiate the Saharan Air Layer Analysis and Mon- aerosols (Guirado et al., 2011). It is near dust sources located
itoring (SALAM) project as part of the Global Atmo- in Mali, southern Algeria, Libya and Chad, on the north-
spheric Watch (GAW) twinning cooperation program be- ern edge of the zonal dust pathway identified by the MISR
tween l’Office Nationale de la Météorologie (ONM, Alge- (Multi-angle Imaging SpectroRadiometer, onboard NASA’s
ria) and the Meteorological State Agency of Spain (AEMET, Terra satellite) AOD retrieval (Fig. 8). Moreover, this ge-
formerly INM) through the Izaña GAW station (Canary Is- ographical location is very significant since ground-based
lands, Spain). In the framework of this project, at the end measurements of atmospheric constituents from continental
of September 2006, a Cimel sun photometer was set up at Africa are very limited, especially in the surrounding area
Tamanrasset and integrated into the Aerosol Robotic Net- of Tamanrasset. This station is involved in several interna-
work (AERONET). In 2010 the station was incorporated tional measurement programs such as the Global Climate
into the World Meteorological Organization (WMO) Sand Observing System (GCOS) – Upper-Air Network (GUAN),
and Dust Storm Warning Advisory and Assessment System the Baseline Surface Radiation Network (BSRN), and the
(SDS-WAS) Regional Center for Northern Africa, Middle GAW program of the WMO.
East and Europe (http://sds-was.aemet.es/) for near-real time The climate of the region is modulated by the influence of
and long-term dust model evaluation. The new aerosol data the monsoon during summer and the westerly winds during
set from Tamanrasset has been used for a preliminary char- the rest of the year (Cuesta et al., 2008). In July and August
acterization of aerosol properties (Guirado et al., 2011), for easterly winds, moist air masses and scarce rainfall are the
space-based remote sensing evaluation (e.g. Schuster et al., prevailing weather conditions. In September the influence of
2012), and for model validation (e.g. Tegen et al., 2013). Re- the westerly winds appears at high altitude and draws suc-
garding dust optical properties, Kim et al. (2011) provide an cessively closer to the ground until the end of autumn. This
analysis of the single scattering albedo, asymmetry param- system is maintained, although wind strengths vary, during
eter, real refractive index, and imaginary refractive index at the winter and even springtime until June when the influence
several stations, including Tamanrasset from 2006 to 2009. of the easterly winds starts in layers close to the ground (Du-
The present work focuses on a detailed characterization bief, 1979). The winter season is characterized by dry con-
of aerosol properties at the Tamanrasset site. Very prelim- ditions and occasional midlevel and cirrus clouds (Cuesta et
inary results, briefly shown by Guirado et al. (2011), have al., 2008).
been carefully revised and extended taking into account cor-
rected data. The KCICLO (K is the name of a constant and
ciclo means cycle in Spanish) method has been used to cor-
rect the aerosol optical depth (AOD) and the Ångström expo-
2.2 Cimel sun photometer data set is difficult, and where dust events lead to dirt accumulating
on the optics. Intense dust storms cause a rapid accumula-
2.2.1 AERONET data tion of dust while less intense dust transport leads to a more
progressive build up of dust on the optics. These inconsisten-
The Cimel sun photometer (model CE-318 operating at 340, cies are difficult to correct and deteriorate the quality of the
380, 440, 500, 670, 870, 940 and 1020 nm nominal wave- measurements.
lengths) installed at Tamanrasset is one of the standard in-
struments in AERONET. Data acquisition protocols, cali-
2.2.2 KCICLO correction
bration procedures and data processing methods have been
extensively described (Holben et al., 1998; Dubovik et al.,
2000; Smirnov et al., 2000; O’Neill et al., 2003). Solar ex- The analysis of the AOD period from 18 November 2007
tinction measurements are used to compute AOD at each to 20 June 2008 reveals a systematic and strong AOD and
wavelength, except for the 940 nm channel, used to retrieve AE fictitious diurnal cycle, most likely caused by dirt on
PWV (Eck et al., 1999). AE, which is a measure of the AOD the sun photometer’s front windows (Guirado et al., 2011).
spectral dependence with the wavelength of incident light, is Other possible causes, such as the effect of temperature on
a qualitative indicator of the aerosol’s predominant particle the detector and incorrect sun-pointing, were analysed and
size and it can be computed for two or more wavelengths discarded. Measurements corresponding to 340 and 500 nm
(Schuster et al., 2006). For climatological studies, linear did not achieve level 2.0 in the entire analysed period due to
fit determination of AE in the 440–870 nm range is com- the significant degradation of these filters.
puted for three or four nominal wavelengths (440 nm, 500 nm The KCICLO method is used to detect, evaluate and cor-
when available, 670 nm, and 870 nm). AERONET FMF from rect possible calibration problems, after discarding a real at-
the spectral deconvolution algorithm (SDA; O’Neill et al., mospheric effect or instrument malfunctions (Cachorro et al.,
2003) has also been included in the present analysis. Further- 2004, 2008a). Particularly, the obstruction in the optical path,
more, several aerosol microphysical and optical properties due to dirt on the sun photometer’s front windows, leads to
retrieved from the AERONET inversion algorithm (Dubovik a distinct and artificial diurnal cycle pattern that can be cor-
and King, 2000; Dubovik et al., 2006) are discussed. Par- rected using the KCICLO method. This fictitious diurnal cy-
ticularly, particle size distribution, volume concentration, ef- cle is due to the systematic absolute error in the AOD mea-
fective radius, as well as single scattering albedo, asymme- surements as a consequence of calibration errors (Romero
try factor, and complex refractive index are analysed because and Cuevas, 2002): the magnitude of this absolute error is
they are closely related to aerosol radiative effects. greatest at midday because it varies as the inverse of the so-
The AOD uncertainty is approximately 0.01–0.02 (spec- lar air mass (Cachorro et al., 2008a). This method introduces
trally dependent, with higher errors in the UV) and it al- a constant K defined as the ratio between “incorrect” cur-
ters the AE by 0.03–0.04 (Eck et al., 1999; Schuster et al., rent and true calibration constants. K quantifies calibration
2006). The PWV uncertainty is estimated to be ± 10 % (Hol- factor error in such a way that K = 1 corresponds to a cor-
ben et al., 2001). The amplitude of the errors of the derived rect calibration constant and K > 1 (K < 1) will result in an
parameters from SDA retrieval varies as the inverse of the to- overestimation (underestimation) of the current calibration
tal AOD. In addition to measurement errors, there are errors constant and a convex (concave) curve shape in the diurnal
in the AOD retrieval due to the uncertainty in the assumed cycle (Cachorro et al., 2004, 2008a). The relative differences
values of the spectral curvature in each mode (O’Neill et of AOD between AERONET level 2.0 and KCICLO data se-
al., 2001), which are most critical in coarse mode-dominated ries are estimated to be 8.5 % (or about 0.01 in absolute AOD
conditions. Dubovik et al. (2002) provide a detailed descrip- values) and 2.4 % for AE (Cachorro et al., 2008b).
tion of the expected error in aerosol size distribution, com- The application of this in situ correction-calibration proce-
plex refractive index, and single scattering albedo. dure requires a sufficient number of clear-sky and stable days
At present, AERONET level 2.0 data from Tamanrasset is for a given period to be corrected. The selected days must
only available from October 2006 to February 2009, except fulfil a set of requirements about air mass range (higher than
from 18 November 2007 to 20 June 2008. Data for the period 4 and typically between 1.7 and 6), turbidity (AOD (440 nm)
February 2009–October 2012 will likely never be promoted < 0.12 and variability lower than 5 % in the specified air mass
to level 2.0, and level 1.5 data in this period do not have the range), number of data points (at least 12 per day), and stan-
sufficient quality to be properly corrected with the KCICLO dard deviation of the fit to quantify the calibration factor er-
method that will be discussed below. Data after November ror (lower than 0.01) (Cachorro et al., 2008a). Therefore, the
2012 are expected to achieve AERONET level 2.0 and might successful application of the KCICLO method over a given
be incorporated in the future to perform a relatively longer- period is associated with a sufficient number of days (5–
term analysis. Long AOD data series fulfilling the highest 10 %) fulfilling all the above-mentioned requirements. As a
quality criteria are difficult to obtain in remote stations such consequence, the application of the method is not always fea-
as Tamanrasset, in which the annual exchange of instruments sible at all stations or at all periods of time. The KCICLO
method has been previously used to correct the AOD data Following the data processing and quality control proce-
series (e.g. Toledano et al., 2007; Barreto el al., 2014). dures described above, AERONET level 2.0 and KCICLO-
At Tamanrasset, a sufficient number of days (94) from corrected level 1.5 data (AOD and AE) were used for aerosol
18 November 2007 to 20 June 2008 were available to prop- characterization. Due to the degradation of the 500 nm fil-
erly apply the KCICLO method and complete the AOD/AE ter, AOD measurements at 440 nm were selected for anal-
data set. This method confirmed a calibration shift between ysis. However, since AOD at 500 nm is more suitable for
November 2007 and June 2008. Only two different correc- satellite and modelling comparisons, it was estimated from
tion periods, i.e. two different types of lens contamination AOD (440 nm) and AE (440–670–870 nm) applying the
(amount of dirt and lenses affected), were detected and the Ångström power law (Ångström, 1929). With regard to the
corresponding mean K values (Table 1) were computed. Part PWV record, AERONET level 1.5 measurements were not
of the original and the corrected AOD and AE data for both affected by any fictitious diurnal cycle. The water vapour
periods is shown in Fig. 1. Note that the fictitious diurnal cy- optical depth, and consequently the PWV product, is not
cle is largely reduced both in the AOD and the derived AE. strongly affected by obstructions in the optical path because
Additionally, it was possible to apply an external qual- the calculation algorithm is based on a subtraction of exper-
ity control of the KCICLO correction. Since 1995, in the imental measurements (Schmid et al., 2001). Therefore, the
framework of the GAW program, a J-309 handheld sun pho- analysed PWV data series comprised AERONET levels 2.0
tometer (Reddy, 1986) supplied by the National Oceanic and and 1.5 when level 2.0 was not available. Limitations and
Table 2. Least squares linear fit results (dimensionless) between the NOAA handheld sun photometer AOD measurements and three Cimel
sun photometer AOD data sets (AERONET level 1.5 before and after KCICLO correction, and AERONET level 2.0). The parameters are the
following: slope of the regression, Y intercept, correlation coefficient (R 2 ), root-mean-square error (RMSE), and number of observations.
special features regarding the analysed AERONET inversion cgi?Dataset=CAL_IIR_L1-Prov-V1-10), have been filtered
retrievals for single scattering albedo and complex refractive following the methodology of Tesche et al. (2013).
index will be discussed in Sect. 3.1.4. All the analysed daily,
monthly, and seasonal averages have been calculated from
2.3.3 Concentration weighted trajectory method
the corresponding sun photometer single measurements.
2.3 Ancillary data The concentration weighted trajectory (CWT) source appor-
tionment method (Seibert et al., 1994; Hsu et al., 2003) was
2.3.1 Meteorological radiosonde data used to identify pathways of aerosol-laden air masses for
the period 2006–2009 in the dry season (from November
A GCOS-GUAN meteorological radiosonde (Vaisala RS92) to February) and the wet season (from April to September).
is launched twice a day (at 00:00 UTC and 12:00 UTC) The resulting information about air mass pathways was com-
at Tamanrasset airport: data available at the University bined with the information about aerosol source regions re-
of Wyoming web site (http://weather.uwyo.edu/upperair/ ported by several authors (Sect. 3.2, Fig. 8) to detect poten-
sounding.html). Radiosonde data at 12:00 UTC were used tial sources affecting Tamanrasset. This method combines
for calculation of the convective boundary layer (CBL) top data measured at the receptor site with air mass back tra-
altitude from 2006 to 2009. The criteria used to account jectories. Although this method was originally designed and
for the overshooting thermals are 1θ/1z ≥ 0.0025 K m−1 widely used for weighting trajectories with concentrations
and θtop − θbase ≥ 1 K, where 1θ/1z is the potential tem- measured at a receptor site, we used AERONET daily AOD
perature lapse rate and θtop and θbase refer to the top and and AE observations at Tamanrasset to identify aerosol con-
base of the layer, respectively (Heffter, 1980; Cuesta et al., tent and type, respectively. A similar approach to connect
2008). Additionally, PWV retrieved from radiosonde was distinct sources with different aerosol types has been previ-
compared with corresponding AERONET PWV as will be ously performed by other authors (e.g. Naseema Beegum et
shown in Sect. 3.1.5. Estimated PWV precision of the ra- al., 2012). A weighted AOD or AE value is assigned to each
diosonde RS92 is around 5 % but for very dry conditions it is grid cell by averaging the values associated with the trajecto-
about 10–20 % (Miloshevich et al., 2009). ries crossing that grid cell:
1.00 80 < nij
0.70 20 < nij ≤ 80
Wij = , (2)
0.42 10 < nij ≤ 20
0.05 nij ≤ 10
Table 3. Monthly means of AOD, AE, and PWV for the period October 2006 to February 2009 at Tamanrasseta .
Month AOD (440 nm)b AE (440–670–870 nm)b PWV (cm) No. of days
January 0.07 (0.08) 0.69 (0.25) 0.37 (0.16) 93
February 0.12 (0.15) 0.49 (0.23) 0.48 (0.23) 66
March 0.23 (0.22) 0.31 (0.17) 0.57 (0.37) 62
April 0.40 (0.39) 0.19 (0.11) 0.64 (0.29) 60
May 0.37 (0.22) 0.15 (0.08) 0.99 (0.29) 62
June 0.43 (0.34) 0.17 (0.14) 0.97 (0.26) 60
July 0.39 (0.32) 0.32 (0.20) 1.15 (0.24) 62
August 0.41 (0.34) 0.44 (0.33) 1.39 (0.45) 62
September 0.33 (0.24) 0.36 (0.20) 1.22 (0.32) 61
October 0.20 (0.14) 0.41 (0.22) 1.01 (0.28) 93
November 0.10 (0.06) 0.54 (0.21) 0.68 (0.24) 90
December 0.09 (0.12) 0.72 (0.25) 0.49 (0.26) 93
a Corresponding standard deviations are shown in brackets.
b Dimensionless.
Figure 7. (a) Monthly means of aerosol particle size distribution (µm3 µm−2 ) at Tamanrasset for the period 2006–2009. Same colours are
used for the dry season (blue), the wet season (orange), and the transition months (grey). (b) Scatter plot of FMF (dimensionless) and fine
mode volume fraction (Vf /Vt ) (dimensionless) for the dry and wet seasons (157 and 183 coincident observations, respectively).
Table 5. Monthly means of volume particle concentration (VolCon) of total, fine and coarse modes, fine mode volume fraction (Vf /Vt ), and
effective radius (Reff ) for the period October 2006–February 2009 at Tamanrasseta .
3.1.4 Aerosol optical properties and imaginary parts of refractive index, respectively, and
0.05 for SSA (Dubovik et al., 2000, 2002). Therefore, no in-
formation of these parameters is available in the AERONET
AERONET level 2.0 retrievals for single scattering albedo
database for the dry season at Tamanrasset. Regarding the
(SSA) and complex refractive index are limited to measure-
wet season, dust optical properties (from March to October)
ments of AOD (440 nm) > 0.4. The reason is that the accuracy
are reported by Kim et al. (2011). To perform an analysis
of these two parameters significantly decreases under lower
for the dry season, we have filtered level 1.5 data following
aerosol loading conditions: 80–100 % and 0.05–0.07 for real
Figure 8. Averaged MISR AOD at 555 nm (dimensionless blue/red scale) for the period 2007–2008 during (a) the dry season (from November
to February) and (b) the wet season (from April to September). Geographical location of Tamanrasset (black star) in the Hoggar Mountains
(Algeria) and present (2014) continuous monitoring AERONET stations (black dots) are indicated. Several potential dust sources, discussed
in the text, have been identified (solid red/blue lines) and numbered as follows: 1, triangle formed by the Adrar des Ifoghas, Hoggar and
Aïr massifs; 2, eastern Libyan desert; 3, Bodélé Depression; 4, western Sahara region; and 5, Libya–Tunisia border. Mean (black solid line)
and median (black dashed line) CALIOP extinction coefficients at 532 nm (km−1 ) are displayed for the period 2007–2008 during (c) the
dry season (43 available profiles) and (d) the wet season (95 available profiles) over Tamanrasset. Grey-shaded area shows the range of
values between the 20th and 80th percentiles. The red line marks the threshold of pristine conditions (extinction coefficient <0.05 km−1 ).
Significant height levels, except the ground level, for the CWT analysis (2600 and 5600 m a.g.l.) are marked (black dotted lines).
the same AERONET criteria but applying a smaller thresh- than expected for dust conditions. Deviations ranging ± 0.05
old to AOD (above 0.1, instead of 0.4, at 440 nm). A similar or more from 1.53 have been previously reported by Dubovik
approach has been previously considered by other authors to et al. (2002) and references therein.
investigate the role of fine aerosols on the absorption of solar
radiation (e.g. Mallet et al., 2013). 3.1.5 Precipitable water vapour
The imaginary part of refractive index, the SSA, and the
asymmetry parameter are relatively constant in the interval
The observed PWV atmospheric content shows an annual cy-
675–1020 nm during both the dry and the wet seasons (Table
cle quite similar to that of the CBL (Fig. 4a, b). The lowest
6). However, the SSA is lower at 440 nm whereas both the
multiannual monthly mean of PWV (Table 3) is observed in
imaginary part (absorption) of the refractive index, as well
January (0.37 ± 0.16 cm) showing a low year to year vari-
as the asymmetry parameter, are higher at 440 nm. These
ability and increases during winter and spring peaking in
spectral patterns of SSA and complex refractive index are
August (1.39 ± 0.45 cm) under the monsoon regime. PWV
consistent with dust-dominant conditions (e.g. Dubovik et
retrieved from radiosondes launched at 12:00 UTC has been
al., 2002). For the dry season this spectral dependency is
compared (not shown) with the corresponding AERONET
smoothed and the asymmetry parameter is slightly lower due
PWV (average of the measurements taken from 12:00 UTC
to the decrease in the coarse mode dominance. It is worth not-
to 13:00 UTC), with a good correlation (0.94) for an over-
ing that the real part of the refractive index ranges between
all number of 610 coincident measurements. The slope of
1.43 and 1.46 during the wet season. These values are lower
the least-squares regression line is 1.14 and the RMSE is
Table 6. Seasonal means (dimensionless) of SSA, real and imaginary parts of the refractive index (Ref. index), and asymmetry parameter
(Asym.) at 440, 675, 870 and 1020 nma . Number of daily available observations (N) is also indicated.
1.15 mm. These results are in good agreement with sim- variability of the cluster classification obtained for each sea-
ilar comparisons, such as that performed by Schneider et son (dry and wet) and for each altitude. For this reason the
al. (2010), who reported a correlation of 0.96 between 675 CWT method was applied to AOD and AE parameters.
AERONET and radiosonde coincident measurements (1 hour Air mass back trajectories at 2600 and 5600 m a.g.l. show
as temporal coincidence criterion) at the Izaña station. a westerly component in the dry season (Fig. 9c, d), driven
by the general circulation, since these levels correspond to
3.2 Potential source regions the free troposphere over the relatively low CBL top. The
dry season is characterized by low AOD and rather high AE
3.2.1 Concentration weighted trajectory analysis associated with short air mass back trajectories at ground
level from the first quadrant (Fig. 9a, b). Dust source re-
Recently, several comprehensive reviews of potential dust gions identified as 1 and 2 in Fig. 8 might potentially affect
sources in Northern Africa (e.g. Formenti et al., 2011; Gi- Tamanrasset in this season. The region located in the triangle
noux et al., 2012) have been provided. However, our goal is formed by Adrar des Ifoghas, Hoggar Mountains and the Aïr
to identify the potential dust sources affecting the Tamanras- Massif (dust source 1; Fig. 8), has been previously identified
set station. This study has been performed through the anal- (D’Almeida, 1986; Prospero, 2002; Schepanski et al., 2009;
ysis of primary air mass pathways and their relationship with Alonso-Pérez et al., 2012) as a Saharan dust source formed
AERONET AOD and AE measurements at Tamanrasset. by a drainage system of ephemeral rivers and streams. This
CALIOP aerosol extinction profiles at 532 nm (Fig. 8c, d) source is sensitive to the effects of mesoscale winds inten-
have been analysed to link aerosol extinctions and air mass sified by the orography (Ginoux et al., 2012). A second po-
pathways at certain heights. The 20th percentile of the extinc- tential dust source (dust source 2; Fig. 8) extends from the
tion in the wet season (Fig. 8d) has been selected as a thresh- northwest side of the Tibesti Mountains in Chad over the
old of pristine conditions. The features of the CBL top identi- eastern Libyan Desert (D’Almeida, 1986; Caquineau et al.,
fied from CALIOP agree quite well with those obtained from 2002; Prospero, 2002; Ginoux et al., 2012). This source is
the radiosondes. Taking into account the averaged CALIOP formed by a large basin with sand seas and the northern part
profiles, HYSPLIT back trajectories at several heights have is marked with a chain of wadis (and associated complexes
been calculated for each day of the period 2007–2008. The of salt/dry lakes). It is active during much of the year but it is
end-point heights of the back trajectories have been selected especially intense in May–June.
according to the CBL top height during both the dry and In relation to AE (Fig. 9b), the highest values (smaller par-
the wet seasons. The three selected height levels provide in- ticles) are found around Tamanrasset. It could be the result
formation about air mass transport near the surface (ground of a mixture of desert dust and local pollution produced by
level), at an intermediate layer (2600 m a.g.l.), which is just cooking and heating bonfires with firewood, common in this
above the CBL top in the dry season and within the CBL region, which is not well dispersed by the low-level atmo-
during the wet season, as well as at 5600 m a.g.l., above the spheric circulation. The potential influence of biomass burn-
CBL (free troposphere) all year long (Figs. 4a, 8c, d). A first ing from the Sahel region to Tamanrasset during the dry sea-
cluster analysis was performed using the k-means clustering son can be considered as nonexistent according to the CWT
algorithm following the Jakob and Tselioudis (2003) proce- analysis.
dure. However, no conclusive results were found due to the
Figure 9. CWT maps at ground level for (a) AOD and (b) AE, and HYSPLIT back trajectories ending at (c) 2600 and (d) 5600 m a.g.l. during
the dry season (from November to February). Refer to dimensionless white/red scales for colour description of AOD and AE. Tamanrasset
is located at the intersection of the four quadrants.
In the wet season, only the 5600 m a.g.l. level is over the 3.2.2 Mesoscale convective system analysis
top of a high CBL typical of summertime, showing mainly
westerly trajectories (Fig. 10e). CWT analysis for AOD and
Mesoscale weather systems (dry boundary layer convection,
AE at ground and 2600 m a.g.l. levels give similar results
“haboob” dust storms, nocturnal low-level jets, and southerly
suggesting a well-mixed CBL in this season. AOD and AE
monsoon flow) influence dust emission, transport, and depo-
CWT plots at ground level and 2600 m a.g.l. (Fig. 10a–d)
sition over the central western Sahara (Marsham et al., 2008,
show a curved dust pathway with relatively high values of
2013; Knippertz and Todd, 2010, 2012; Ashpole and Wash-
AE (smaller particles) from northern central Libya passing
ington, 2013). During 2006, Cuesta et al. (2008) observed
over dust sources 1 and 2 as occurs in the dry season. A
several summertime dust transport events over Tamanras-
second curved dust pathway from the Libya–Tunisia bor-
set associated with MCSs (mesoscale convective systems).
der (Caquineau et al., 2002) (dust source 5; Fig. 8) is ob-
However, MCSs cannot be well captured by global meteo-
served to transport larger particles (low AE) to Tamanrasset
rological models or regional dust models (Marsham et al.,
at ground level. A few air mass trajectories originate from the
2011; Heinold et al., 2013) nor by the HYSPLIT back-
west passing over the large dust source 4 (Fig. 8) located in
trajectory parameterization. Consequently, we performed an
northern Mali, northern Mauritania and the western flanks of
additional analysis to identify the influence of MCSs on dust
the Hoggar Mountains (Prospero, 2002; Brooks and Legrand,
transport events over Tamanrasset.
2003, Alonso-Pérez et al., 2012). It is a complex distribution
We have analysed 21 episodes of MCSs that have been
of dust sources with extensive dune systems that are particu-
selected through comparison between observed AERONET
larly active from April to September.
AOD and NMMB/BSC-Dust model AOD over Taman-
Regarding one of the most significant dust sources in the
rasset in the period 2007–2008 (Fig. 11). The simula-
world, the Bodélé Depression (Goudi and Middleton, 2001;
tion of the Non-hydrostatic Multiscale Model (NMMB)
Prospero, 2002; Brooks and Legrand, 2003) (dust source 3,
Barcelona Supercomputing Center (BSC) v1 is generated
Fig. 8), CWT analysis shows that it is a minor dust source
using the National Center for Environmental Prediction
affecting Tamanrasset.
(NCEP) reanalysis-2 (1◦ grid) and initial and boundary con-
ditions from the Global Land Data Assimilation System
(GLDAS). The resolution is set at 0.5◦ in the horizontal
Figure 10. CWT maps for AOD and AE at (a and b) ground level and (c and d) 2600 m a.g.l., and (e) HYSPLIT back trajectories ending at
5600 m a.g.l., during the wet season (from April to September). Refer to dimensionless white/red scales for colour description of AOD and
AE. Tamanrasset is located at the intersection of the four quadrants.
and to 40 hybrid sigma-pressure model layers in the verti- Forecasts (ECMWF) ERA-Interim reanalysis data from the
cal. A detailed description of the model is provided by Pérez IFS-Cy31r model analysis. Satellite information and mete-
et al. (2011). The NMMB/BSC-Dust model properly repro- orological data were jointly computed and visualized with
duces dust transport associated with synoptic-scale meteoro- McIDAS (Man computer Interactive Data Access System)
logical processes observed during most of the year (Fig. 11). software.
However, from June to September, although the AOD trend is Once we have identified and confirmed all the MCS events
well reproduced, the model is not capable of capturing strong impacting Tamanrasset, the Moderate Resolution Imaging
and fast dust outbreaks associated with MCSs. The summer- Spectroradiometers (MODIS) Deep Blue 550 nm AOD re-
time observation-model AOD discrepancies have been used trieval has been used in a similar approach to Roberts (2014)
to identify the potential MCSs affecting Tamanrasset. The and Roberts et al. (2014). The advantage of the MODIS
convective origin of each event has been evaluated by using Deep Blue aerosol retrieval algorithm regarding other satel-
high temporal and spatial RGB (red, green, blue) dust com- lite products over bright surfaces in the visible (such as
posites from the Meteosat Second Generation-Spinning En- deserts) is that the former employs radiances from the blue
hanced Visible and Infrared Imager (MSG-SEVIRI) sensor channels where the surface reflectance is relatively low (Hsu
combined with European Centre for Medium-range Weather et al., 2004; 2006). The MODIS Deep Blue composite AOD
Figure 12. Composite MODIS Deep Blue 550 nm (a) AOD and AOD-averaged anomaly corresponding (b) to the 21 days of maximum
(Dmax ) AOD at Tamanrasset during MCS events, and (c) to the previous day for each maximum AOD event (Dmax -1). The Tamanrasset
station is marked with a black star. Two-day HYSPLIT back trajectories arriving at Tamanrasset at ground level (black solid lines) are also
displayed in panel (a).
nitrate aerosols. These results agree with the Rodríguez et months. The AOD shows the same but opposite seasonal
al. (2011) observations. cycle to the AE and FMF. AOD remains stable around 0.1
from November to February but exceeds 0.3 from April to
September, reaching an absolute maximum of 0.43 in June.
4 Summary and conclusions The maximum AE and FMF observations are reached in De-
cember (0.72 and 0.57, respectively) and January (0.69 and
Tamanrasset is a strategic site for aerosol research placed in
0.58, respectively), tending to decrease until May (minimum
the heart of the Sahara. An aerosol characterization at this
of 0.15 and 0.25, respectively). Minimum PWV is recorded
site has been provided based on more than 2 years (October
in January (0.37 ± 0.16 cm) whereas maximum values are
2006–February 2009) of AERONET level 2.0 and KCICLO-
reached in August (1.39 ± 0.45 cm) linked to the monsoon
corrected Cimel sun photometer measurements. The top of
regime.
the CBL over Tamanrasset has been characterized by both
Coarse mode (modal radius around 2.24 µm) prevails over
radiosonde data and CALIOP extinction vertical profiles. A
the fine mode (modal radius around 0.10 µm) showing lower
strong seasonal cycle linked to the CBL is observed. The
volume concentrations during the dry season and maxima in
dry-cool season (November–February) is characterized by a
July. Spectral patterns of SSA and complex refractive index
shallow CBL, low AOD (∼ 0.09 at 440 nm), moderate-low
also indicate coarse mode dominance conditions. However,
AE values (∼ 0.62) and low PWV (∼ 0.51 cm). The wet-hot
FMF and fine mode volume fraction show values correspond-
season (April–September) is characterized by a deep CBL,
ing to fine mode dominance in terms of optical depth, and
higher AOD (∼ 0.39 at 440 nm), low AE (∼ 0.28) and higher
coarse mode dominance in terms of volume concentration
PWV (∼ 1.06 cm) and it is affected by strong and frequent
during the dry season. In addition, around 60 % of the AOD
dust storms. March and October are considered transition
Figure 13. HYSPLIT back trajectories arriving at Tamanrasset (blue pin) at ground level (red lines) and 2600 m a.g.l. (yellow lines) are
displayed for several case studies. DMSP night-time lights (shown as background) identify gas flares (green colour).
measurements are below 0.15, showing a cleaner atmosphere rasset. Dust uplift sources associated with summertime MCS
than sites located in the Sahel. The AE dry season distribu- events located south of Tamanrasset have been also identi-
tion is slightly bimodal (0.4 and 0.7 modal values) and sym- fied.
metrical, indicating a similar frequency for the two differ- Evidence of the arrival of fine particles to Tamanrasset dur-
ent particle populations (desert dust and background condi- ing summer has been detected and nine events of polluted air
tions). During the wet season AE shows two clear modes, a masses coming from urban/industrial areas in Libya and Al-
narrowed first mode centred at 0.15 (high dust events) and geria have been shown.
a smaller but wider mode centred at 0.4 (background condi-
tions).
The aerosol-type characterization at Tamanrasset indicates Acknowledgements. The AERONET sun photometer at Tamanras-
desert mineral dust is the prevailing aerosol. Higher extinc- set has been calibrated within AERONET-EUROPE TNA supported
by the PHOTONS and RIMA networks and partially financed by
tions linked to larger particles and dominance of the one-
the European Community – Research Infrastructure Action under
particle mode have been identified by the observation of
the Seventh Framework Programme (FP7/2007-2013) “Capacities”
large variations of AOD, AE almost inversely proportional specific programme for Integrating Activities, ACTRIS grant
to AOD, and the Ångström exponent difference (δAE), neg- agreement no. 262254. The authors gratefully acknowledge the
ative or slightly positive. Moreover, typical pure Saharan NOAA Air Resources Laboratory (ARL) for the provision of the
dust conditions have been observed, i.e. high-extinction val- HYSPLIT transport and dispersion model and READY website
ues (AOD > 0.7) with AE < 0.3 and δAE < 0 corresponding (http://ready.arl.noaa.gov) used in this publication. We also thank
to FMF < 40 % and fine mode radius (Rf ) of around 0.3 µm. Google Earth™ . Financial support from the Spanish MINECO
However, an anthropogenic fine mode contribution has been (projects CGL2011-23413, CGL2012-33576 and CGL2012-37505)
found mixed with mineral dust (8.7% of total cases), i.e. is also gratefully acknowledged. We are particularly grateful to
AOD of up to 0.4 for AE values ranging between 0.6 and 1.1 the Tamanrasset Global Atmospheric Watch (GAW) station’s staff
(l’Office National de la Météorologie, Algeria) for supporting the
and δAE, FMF and Rf between −0.3 and 0.2, 30 and 70 %,
measurement program. J. M. Baldasano and S. Basart acknowledge
and 0.10 and 0.20 µm, respectively. the Supercomputación y eCiencia project (CSD2007-0050) from
Potential sources of the air masses arriving at Tamanras- the Consolider-Ingenio 2010 and Severo Ochoa (SEV-2011-00067)
set have been identified: the triangle formed by Adrar des programs of the Spanish Government. We also acknowledge our
Ifoghas, Hoggar Mountains and Aïr Massif; a complex dis- colleague Dr. Celia Milford for improving the English language of
tribution of dust sources including northern Mali, northern the manuscript.
Mauritania and the western flanks of the Hoggar Mountains;
and the eastern Libyan Desert. However, the Bodélé Depres- Edited by: W. Lahoz
sion has been found to be a minor potential source at Taman-
References Caquineau, S., Gaudichet, A., Gomes, L., and Legrand, M.: Min-
eralogy of Saharan dust transported over northwestern tropical
Alonso-Pérez, S., Cuevas, E., Querol, X., Viana, M., and Atlantic Ocean in relation to source regions, J. Geophys. Res.,
Guerra, J. C.: Impact of the Saharan dust outbreaks on 107, AAC 4-1–AAC 4-12, doi:10.1029/2000JD000247, 2002.
the ambient levels of total suspended particles (TSP) in Cuesta, J., Edouart, D., Mimouni, M., Flamant, P. H., Loth, C.,
the Marine Boundary Layer (MBL) of the Subtropical East- Gibert, F., Marnas, F., Bouklila, A., Kharef, M., Ouchene, B.,
ern North Atlantic Ocean, Atmos. Environ., 41, 9468–9480, Kadi, M., and Flamant, C.: Multiplatform observations of the
doi:10.1016/j.atmosenv.2007.08.049, 2007. seasonal evolution of the Saharan atmospheric boundary layer in
Alonso-Pérez, S., Cuevas, E., Querol, X., Guerra, J. C., Tamanrasset, Algeria, in the framework of the African Monsoon
and Pérez, C.: African dust source regions for observed Multidisciplinary Analysis field campaign conducted in 2006,
dust outbreaks over the Subtropical Eastern North At- J. Geophys. Res., 113, D00C07, doi:10.1029/2007JD009417,
lantic region above 25◦ N, J. Arid Environ., 78, 100–109, 2008.
doi:10.1016/j.jaridenv.2011.11.013, 2012. Cuesta, J., Marsham, J. H., Parker, D. J., and Flamant, C.: Dynami-
Ångström, A.: On the atmospheric transmission of sun radiation and cal mechanisms controlling the vertical redistribution of dust and
on dust in the air, Geogr. Ann., 2, 156–166, 1929. the thermodynamic structure of the West Saharan atmospheric
Ashbaugh, L. L., Malm, W. C., and Sadeh, W. D.: A residence time boundary layer during summer, Atmos. Sci. Lett., 10, 34–42,
probability analysis of sulfur concentrations at Grand Canyon doi:10.1002/asl.207, 2009.
National Park, Atmos. Environ., 19, 1263–1270, 1985. Cuesta, J., Lavaysse, C., Flamant, C., Mimouni, M., and Knip-
Ashpole, I. and Washington, R.: Intraseasonal variability and atmo- pertz, P.: Northward bursts of the West African monsoon leading
spheric controls on daily dust occurrence frequency over the cen- to rainfall over the Hoggar Massif, Algeria, Q. J. Roy. Meteor.
tral and western Sahara during the boreal summer, J. Geophys. Soc., 136, 174–189, doi:10.1002/qj.439, 2010.
Res. Atmos., 118, 12915–12926, doi:10.1002/2013JD020267, D’Almeida, G. A.: A model for Saharan dust transport, J. Clim.
2013. Appl. Meteorol., 25, 903–916, 1986.
Barreto, A., Cuevas, E., Pallé, P., Romero, P.M., Almansa, F., and D’Almeida, G. A.: On the variability of desert aerosol radiative
Wehrli, C.: Recovering Long-term Aerosol Optical Depth Series characteristics, J. Geophys. Res., 92, 3017–3026, 1987.
(1976–2012) from an Astronomical Potassium-based Resonance D’Almeida, G. A. and Schütz, L.: Number, mass and volume dis-
Scattering Spectrometer, Atmos. Meas. Tech. Discuss., 7, 4093– tributions of mineral aerosols and soils of Sahara, J. Clim. Appl.
4121, doi:10.5194/amtd-7-4093-2014, 2014. Meteorol., 22, 233–243, 1983.
Basart, S., Pérez, C., Cuevas, E., Baldasano, J. M., and Gobbi, G. P.: De Longueville, F., Hountondji, Y. C., Henry, S., and Ozer, P.: What
Aerosol characterization in Northern Africa, Northeastern At- do we know about effects of desert dust on air quality and human
lantic, Mediterranean Basin and Middle East from direct-sun health in West Africa compared to other regions?, Sci. Total En-
AERONET observations, Atmos. Chem. Phys., 9, 8265–8282, viron., 409, 1–8, doi:10.1016/j.scitotenv.2010.09.025, 2010.
doi:10.5194/acp-9-8265-2009, 2009. Draxler, R. R. and Hess, G. D.: An overview of the HYSPLIT_4
Bou Karam, D., Flamant, C., Knippertz, P., Reitebuch, O., Pelon, J., modelling system for trajectories, dispersion, and deposition,
Chong, M., and Dabas, A.: Dust emissions over the Sahel asso- Aust. Meteorol. Mag., 47, 295–308, 1998.
ciated with the West African monsoon intertropical discontinuity Dubief, J.: Review of the North African climate with particular em-
region: a representative case-study, Q. J. Roy. Meteor. Soc., 134, phasis on the production of eolian dust in the Sahel Zone and in
621–634, doi:10.1002/qj.244, 2008. the Sahara, in: Saharan Dust: Mobilization, Transport, Deposi-
Brooks, N. and Legrand, M.: Dust variability over northern Africa tion, edited by: Morales, C., John Wiley and Sons Ltd., Hoboken,
and rainfall in the Sahel, in: Linking Climate Change to Land NJ, 27–48, 1979.
Surface Change, edited by: McLaren, S. J. and Kniverton, D. R., Dubovik, O. and King, M. D.: A flexible inversion algorithm for
Kluwer Academic Publishers, Dordrecht, Netherlands, 1–25, retrieval of aerosol optical properties from sun and sky radiance
2003. measurements, J. Geophys. Res., 105, 20673–20696, 2000.
Cachorro, V. E., Romero, P. M., Toledano, C., Cuevas, E., and Dubovik, O., Smirnov, A., Holben, B. N., King, M. D., Kauf-
de Frutos, A. M.: The fictitious diurnal cycle of aerosol op- man, Y. J., Eck, T. F., and Slutsker, I.: Accuracy assessments of
tical depth: a new approach for in situ calibration and cor- aerosol optical properties retrieved from Aerosol Robotic Net-
rection of AOD data series, Geophys. Res. Lett., 31, L12106, work (AERONET) Sun and sky radiance measurements, J. Geo-
doi:10.1029/2004GL019651, 2004. phys. Res., 105, 9791–9806, 2000.
Cachorro, V. E., Toledano, C., Berjón, A., de Frutos, A. M., Dubovik, O., Holben, B. N., Eck, T. F., Smirnov, A., Kaufman, Y. J.,
Torres, B., Sorribas, M., and Laulainen, N. S.: An “in King, M. D., Tanré, D., and Lutsker, I.: Variability of absorption
situ” calibration correction procedure (KCICLO) based on and optical properties of key aerosols types observed in world-
AOD diurnal cycle: application to AERONET–El Arenosillo wide locations, J. Atmos. Sci., 59, 590–608, 2002.
(Spain) AOD data series, J. Geophys. Res., 113, D12205, Dubovik, O., Sinyuk, A., Lapyonak, T., Holben, B. N.,
doi:10.1029/2007JD009673, 2008a. Mishchenko, M., Yang, P., Eck, T. F., Volten, H., Muñoz, O.,
Cachorro, V. E., Toledano, C, Sorribas, M., Berjón, A., de Fru- Veihelmann, B., van der Zande, W. J., Leon, J.-F., Sorokin, M.,
tos, A. M., and Laulainen, N.: An “in situ” calibration-correction and Slutsker, I.: Application of spheroid models to account for
procedure (KCICLO) based on AOD diurnal cycle: comparative aerosol particle nonsphericity in remote sensing of desert dust, J.
results between AERONET and reprocessed (KCICLO method) Geophys. Res., 111, D11208, doi:10.1029/2005JD006619, 2006.
AOD-alpha data series at El Arenosillo, Spain, J. Geophys. Res.,
113, D02207, doi:10.1029/2007JD009001, 2008b.
Eck, T. F., Holben, B. N., Reid, J. S., Dubovik, O., Smirnov, A., Heffter, J. L.: Air Resources Laboratories Atmospheric Transport
O’Neill, N. T., Slutsker, I., and Kinne, S.: Wavelength depen- and Dispersion Model (ARL-ATAD), Air Resources Laborato-
dence of the optical depth of biomass burning, urban, and desert ries, Silver Spring, Maryland, 1980.
dust aerosol, J. Geophys. Res., 104, 31333–31350, 1999. Heinold, B., Knippertz, P., Marsham, J. H., Fiedler, S., Dixon,
Elvidge, C. D., Baugh, K. E., Kihn, E. A., Kroehl, H. W., and N. S., Schepanski, K., Laurent, B., and Tegen, I.: The role of
Davis, E. R.: Mapping of city lights using DMSP Operational deep convection and nocturnal low-level jets for dust emission
Linescan System data, Photogramm. Eng. Rem. S., 63, 727–734, in summertime West Africa: Estimates from convection per-
1997. mitting simulations, J. Geophys. Res. Atmos., 118, 4385–4400,
Flamant, C., Chaboureau, J.-P., Parker, D. J., Taylor, C. M., Cam- doi:10.1002/jgrd.50402, 2013.
mas, J.-P., Bock, O., Timouk, F., and Pelon, J.: Airborne obser- Holben, B. N., Eck, T. F., Slutsker, I., Tanré, D., Buis, J. P., Setzer,
vations of the impact of a convective system on the planetary A., Vermote, E., Reagan, J. A., Kaufman, Y. J., Nakajima, T.,
boundary layer thermodynamics and aerosol distribution in the Lavenu, F., Jankowiak, I., and Smirnov, A.: AERONET – a fed-
inter-tropical discontinuity region of the West African Monsoon, erated instrument network and data archive for aerosol character-
Q. J. Roy. Meteor. Soc., 133, 1175–1189, doi:10.1002/qj.97, ization, Remote Sens. Environ., 66, 1–16, 1998.
2007. Holben, B. N., Tanré, D., Smirnov, A., Eck, T. F., Slutsker, I.,
Formenti, P., Schütz, L., Balkanski, Y., Desboeufs, K., Ebert, M., Abuhassan, N., Newcomb, W. W., Schafer, J. S., Chatenet,
Kandler, K., Petzold, A., Scheuvens, D., Weinbruch, S., and B., Lavenu, F., Kaufman, Y. J., Vande Castle, J., Setzer, A.,
Zhang, D.: Recent progress in understanding physical and chem- Markham, B., Clark, D., Frouin, R., Halthore, R., Karneli, A.,
ical properties of African and Asian mineral dust, Atmos. Chem. O’Neill, N. T., Pietras, C., Pinker, R. T., Voss, K. and Zibordi,
Phys., 11, 8231–8256, doi:10.5194/acp-11-8231-2011, 2011. G.: An emerging ground-based aerosol climatology: aerosol opti-
Forster, P., Ramaswamy, V., Artaxo, P., Berntsen, T., Betts, R., Fa- cal depth from AERONET, J. Geophys. Res., 106, 12067–12097,
hey, D. W., Haywood, J., Lean, J., Lowe, D. C., Myhre, G., doi:10.1029/2001JD900014, 2001.
Nganga, J., Prinn, R., Raga, G., Schulz, M., and Dorland, R. V.: Hsu, Y. K., Holsen, T. M., and Hopke, P. K.: Comparison of hy-
Changes in atmospheric constituents and in radiative forcing, in: brid receptor models to locate PCB sources in Chicago, Atmos.
Climate Change 2007: The Physical Science Basis, Contribution Environ., 37, 545–562, 2003.
of Working Group I to the Fourth Assessment Report of the Inter- Hsu, N. C., Tsay, S.-C., King, M. D., and Herman, J. R.: Aerosol
governmental Panel on Climate Change, edited by: Solomon, S., properties over bright-reflecting source regions, IEEE T. Geosci.
Qin, D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B., Tig- Remote Sens., 42, 557–569, 2004.
nor, M., and Miller, H. L., Cambridge University Press, Cam- Hsu, N. C., Tsay, S.-C., King, M. D., and Herman, J. R.: Deep Blue
bridge, UK and New York, NY, USA, 129–234, 2007. retrievals of Asian aerosol properties during ACE-Asia, IEEE T.
Ginoux, P., Prospero, J. M., Gill, T. E., Hsu, N. C., and Geosci. Remote Sens., 44, 3180–3195, 2006.
Zhao, M.: Global-scale attribution of anthropogenic and nat- Hunt, W. H., Winker, D. M., Vaughan, M. A., Powell, K. A., Lucker,
ural dust sources and their emission rates based on MODIS P. L., and Weimer, C.: CALIPSO lidar description and perfor-
Deep Blue aerosol products, Rev. Geophys., 50, RG3005, mance assessment, J. Atmos. Oceanic Technol., 26, 1214–1228,
doi:10.1029/2012RG000388, 2012. doi:10.1175/2009JTECHA1223.1, 2009.
Gobbi, G. P., Kaufman, Y. J., Koren, I., and Eck, T. F.: Classification Jakob, C. and Tselioudis, G.: Objective identification of cloud
of aerosol properties derived from AERONET direct sun data, regimes in the Tropical Western Pacific, Geophys. Res. Lett., 30,
Atmos. Chem. Phys., 7, 453–458, doi:10.5194/acp-7-453-2007, 2082, doi:10.1029/2003GL018367, 2003.
2007. Kim, D., Chin, M., Yu, H., Eck, T. F., Sinyuk, A., Smirnov, A.,
Goudie, A. S. and Middleton, N. J.: Saharan dust storms: nature and and Holben, B. N.: Dust optical properties over North Africa
consequences, Earth-Sci. Rev., 56, 179–204, doi:10.1016/S0012- and Arabian Peninsula derived from the AERONET dataset,
8252(01)00067-8, 2001. Atmos. Chem. Phys., 11, 10733–10741, doi:10.5194/acp-11-
Guirado, C., Cuevas, E., Cachorro, V. E., Mimouni, M., Zeudmi, L., 10733-2011, 2011.
Toledano, C., Alonso-Pérez, S., Basart, S., Blarel, L., Goloub, P., Knippertz, P. and Todd, M. C.: The central west Saharan dust
and Baldasano, J. M.: Preliminary characterization of columnar hot spot and its relation to African easterly waves and
aerosols properties (AOD–AE) at the Saharan Tamanrasset (Al- extratropical disturbances, J. Geophys. Res., 115, D12117,
geria) station, Óptica Pura y Aplicada, 44, 635–639, 2011. doi:10.1029/2009JD012819, 2010.
Haustein, K., Pérez, C., Baldasano, J. M., Müller, D., Tesche, M., Knippertz, P. and Todd, M. C.: Mineral dust aerosols over the
Schladitz, A., Freudenthaler, V., Heese, B., Esselborn, M., Sahara: Meteorological controls on emission and transport
Weinzierl, B., Kandler, K., and von Hoyningen-Huene, W.: Re- and implications for modeling, Rev. Geophys. 50, RG1007,
gional dust model performance during SAMUM 2006, Geophys. doi:10.1029/2011RG000362, 2012.
Res. Lett., 36, L03812, doi:10.1029/2008GL036463, 2009. Mallet, M., Dubovik, O., Nabat, P., Dulac, F., Kahn, R., Sciare, J.,
Haustein, K., Pérez, C., Baldasano, J. M., Jorba, O., Basart, S., Paronis, D., and Léon, J. F.: Absorption properties of Mediter-
Miller, R. L., Janjic, Z., Black, T., Nickovic, S., Todd, M. C., ranean aerosols obtained from multi-year ground-based remote
Washington, R., Müller, D., Tesche, M., Weinzierl, B., Esselborn, sensing observations, Atmos. Chem. Phys., 13, 9195–9210,
M., and Schladitz, A.: Atmospheric dust modeling from meso doi:10.5194/acp-13-9195-2013, 2013.
to global scales with the online NMMB/BSC-Dust model – Part Marsham, J. H., Parker, D. J., Grams, C. M., Taylor, C. M.,
2: Experimental campaigns in Northern Africa, Atmos. Chem. and Haywood, J. M.: Uplift of Saharan dust south of the
Phys., 12, 2933–2958, doi:10.5194/acp-12-2933-2012, 2012.
intertropical discontinuity, J. Geophys. Res., 113, D21102, Redelsperger, J.-L., Thorncroft, C., Diedhiou, A., Lebel, T.,
doi:10.1029/2008JD009844, 2008. Parker, D., and Polcher, J.: African monsoon multidisciplinary
Marsham, J. H., Knippertz, P., Dixon, N. S., Parker, D. J., analysis – an international research project and field campaign,
and Lister, G. M. S.: The importance of the representation B. Am. Meteorol. Soc., 87, 1739–1746, doi:10.1175/BAMS-87-
of deep convection for modeled dust-generating winds over 12-1739, 2006.
West Africa during summer, Geophys. Res. Lett., 38, L16803, Roberts, A. J.: Anomalously heavy rainfall and dust in the arid Sa-
doi:10.1029/2011GL048368, 2011. hara and northern Sahel, in: Convective Episodes near the In-
Marsham, J. H., Hobby, M., Allen, C. J. T., Banks, J. R., Bart, tertropical Discontinuity in Summertime West Africa: Represen-
M., Brooks, B. J., Cavazos-Guerra, C., Engelstaedter, S., Gas- tation in Models and Implications for Dust Uplift, PhD thesis,
coyne, M., Lima, A.R., Martins, J. V., McQuaid, J. B., O’Leary, University of Leeds, Leeds, UK, 2014.
A., Ouchene, B., Ouladichir, A., Parker, D. J., Saci, A., Salah- Roberts, A. J., Knippertz, P., and Marsham, J. H.: The Formation
Ferroudj, M., Todd, M. C., and Washington, R.: Meteorology and of Convectively Generated Dusty Episodes in the Sahara during
dust in the central Sahara: Observations from Fennec supersite-1 Summer, DUST-2014, International Conference on Atmospheric
during the June 2011 Intensive Observation Period, J. Geophys. Dust, Castellaneta Marina, Italy, June 1–6, 2014.
Res. Atmos., 118, 4069–4089, doi:10.1002/jgrd.50211, 2013. Rodríguez, S., Alastuey, A., Alonso-Pérez, S., Querol, X.,
Miloshevich, L. M., Vömel, H., Whilteman, D. N., and Leblanc, T.: Cuevas, E., Abreu-Afonso, J., Viana, M., Pérez, N., Pandolfi, M.,
Accuracy assessment and correction of Vaisala RS92 radiosonde and de la Rosa, J.: Transport of desert dust mixed with North
water vapor measurements, J. Geophys. Res., 114, D11305, African industrial pollutants in the subtropical Saharan Air
doi:10.1029/2008JD011565, 2009. Layer, Atmos. Chem. Phys., 11, 6663–6685, doi:10.5194/acp-
Naseema Beegum, S., Krishna Moorthy, K., Gogoi, Mukunda M., 11-6663-2011, 2011.
Suresh Babu, S., and Pandey, S. K.: Multi-year investigations of Romero, P. M. and Cuevas, E.: Diurnal variation of the aerosol op-
aerosols from an island station, Port Blair, in the Bay of Bengal: tical depth: artifact or reality?, in: Proceeding of 3a Asamblea
climatology and source impacts, Ann. Geophys., 30, 1113–1127, Hispano Portuguesa de Geofísica y Geodesia, Valencia, Spain,
doi:10.5194/angeo-30-1113-2012, 2012. 4–8 February 2002, 2, 1252–1256, 2002.
O’Neill, N. T., Dubovik, O., and Eck, T. F.: Modified Ångström ex- Ryder, C. L., Highwood, E. J., Rosenberg, P. D., Trembath, J.,
ponent for the characterization of submicrometer aerosols, Appl. Brooke, J. K., Bart, M., Dean, A., Crosier, J., Dorsey, J., Brind-
Optics, 40, 2368–2375, 2001. ley, H., Banks, J., Marsham, J. H., McQuaid, J. B., Sodemann, H.,
O’Neill, N. T., Eck, T. F., Smirnov, A., Holben, B. N., and Thulasir- and Washington, R.: Optical properties of Saharan dust aerosol
aman, S.: Spectral discrimination of coarse and fine mode optical and contribution from the coarse mode as measured during the
depth, J. Geophys. Res., 108, 4559, doi:10.1029/2002JD002975, Fennec 2011 aircraft campaign, Atmos. Chem. Phys., 13, 303–
2003. 325, doi:10.5194/acp-13-303-2013, 2013.
Pérez, C., Haustein, K., Janjic, Z., Jorba, O., Huneeus, N., Bal- Schepanski, K., Tegen, I., Todd, M. C., Heinold, B., Bönisch, G.,
dasano, J. M., Black, T., Basart, S., Nickovic, S., Miller, R. Laurent, B., and Macke, A.: Meteorological processes forcing
L., Perlwitz, J. P., Schulz, M., and Thomson, M.: Atmospheric Saharan dust emission inferred from MSG-SEVIRI observations
dust modeling from meso to global scales with the online of subdaily dust source activation and numerical models, J. Geo-
NMMB/BSC-Dust model – Part 1: Model description, annual phys. Res., 114, D10201, doi:10.1029/2008JD010325, 2009.
simulations and evaluation, Atmos. Chem. Phys., 11, 13001– Schmid, B., Michalsky, J. J., Slater, D. W., Barnard, J. C.,
13027, doi:10.5194/acp-11-13001-2011, 2011. Halthore, R. N., Liljegren, J. C., Holben, B. N., Eck, T. F.,
Perez, L., Tobías, A., Querol, X., Pey, J., Alastuey, A., Díaz, J., and Livingston, J. M., Russell, P. B., Ingold, T., and Slutsker, I.:
Sunyer, J.: Saharan dust, particulate matter and cause-specific Comparison of columnar water-vapor measurements from so-
mortality: a case–crossover study in Barcelona (Spain), Environ. lar transmittance methods, Appl. Optics, 40, 1886–1896,
Int., 48, 150–155, doi:10.1016/j.envint.2012.07.001, 2012. doi:10.1364/AO.40.001886, 2001.
Polissar, A. V., Hopke, P. K., Paatero, P., Kaufmann, Y. J., Schneider, M., Romero, P. M., Hase, F., Blumenstock, T.,
Hall, D. K., Bodhaine, B. A., Dutton, E. G., and Harris, J. M.: Cuevas, E., and Ramos, R.: Continuous quality assessment of at-
The aerosol at Barrow, Alaska: long-term trends and source lo- mospheric water vapour measurement techniques: FTIR, Cimel,
cations, Atmos. Environ., 33, 2441–2458, 1999. MFRSR, GPS, and Vaisala RS92, Atmos. Meas. Tech., 3, 323–
Prats, N., Cachorro, V. E., Berjón, A., Toledano, C., and De Fru- 338, doi:10.5194/amt-3-323-2010, 2010.
tos, A. M.: Column-integrated aerosol microphysical proper- Schuster, G. L., Dubovik, O., and Holben, B. N.: Angstrom expo-
ties from AERONET Sun photometer over southwestern Spain, nent and bimodal aerosol size distributions, J. Geophys. Res.,
Atmos. Chem. Phys., 11, 12535–12547, doi:10.5194/acp-11- 111, D07207, doi:10.1029/2005JD006328, 2006.
12535-2011, 2011. Schuster, G. L., Vaughan, M., MacDonnell, D., Su, W., Winker, D.,
Prospero, J. M., Ginoux, P., Torres, O., Nicholson, S. E., and Dubovik, O., Lapyonok, T., and Trepte, C.: Comparison of
Gill, T. E.: Environmental characterization of global sources of CALIPSO aerosol optical depth retrievals to AERONET mea-
atmospheric soil dust identified with the nimbus 7 total ozone surements, and a climatology for the lidar ratio of dust, Atmos.
mapping spectrometer (TOMS) absorbing aerosol product, Rev. Chem. Phys., 12, 7431–7452, doi:10.5194/acp-12-7431-2012,
Geophys., 40, 1002, doi:10.1029/2000RG000095, 2002. 2012.
Reddy, P. J.: Instructions for J-Series Handheld Sunphotometer,
NOAA, Boulder, 1986.
Seibert, P., Kromp-Kolb, H., Baltensperger, U., Jost, D. T., Todd, M. C., Allen, C. J. T., Bart, M., Bechir, M., Bentefouet, J.,
Schwikowski, M., Kasper, A., and Puxbaum, H.: Trajectory anal- Brooks, B. J., Cavazos-Guerra, C., Clovis, T., Deyane, S., Dieh,
ysis of aerosol measurements at High Alpine Sites, in: Transport M., Engelstaedter, S., Flamant, C., Garcia–Carreras, L., Gan-
and Transformation of Pollutants in the Troposphere, edited by: dega, A., Gascoyne, M., Hobby, M., Kocha, C., Lavaysse, C.,
Borrell, P. M., Borrell, P., Cvitas, T., and Seiler, W., Academic Marsham, J. H., Martins, J. V., McQuaid, J. B., Ngamini, J. B.,
Publishing, Den Haag, 689–693, 1994. Parker, D. J., Podvin, T., Rocha-Lima, A., Traore, S., Wang, Y.,
Smirnov, A., Holben, B. N., Eck, T. F., Dubovik, O., and and Washington, R.: Meteorological and dust aerosol conditions
Slutsker, I.: Cloud screening and quality control algorithms for over the western Saharan region observed at Fennec Supersite-2
the AERONET data base, Remote Sens. Environ., 73, 337–349, during the intensive observation period in June 2011, J. Geophys.
2000. Res.-Atmos., 118, 8426–8447, doi:10.1002/jgrd.50470, 2013.
Su, L. and Toon, O. B.: Saharan and Asian dust: similarities and Toledano, C., Cachorro, V. E., Berjon, A., de Frutos, A. M., Sor-
differences determined by CALIPSO, AERONET, and a coupled ribas, M., de la Morena, B. A., and Goloub, P.: Aerosol opti-
climate-aerosol microphysical model, Atmos. Chem. Phys., 11, cal depth and Ångström exponent climatology at El Arenosillo
3263–3280, doi:10.5194/acp-11-3263-2011, 2011. AERONET site (Huelva, Spain), Q. J. Roy. Meteor. Soc., 133,
Tegen, I., Schepanski, K., and Heinold, B.: Comparing two years 795–807, doi:10.1002/qj.54, 2007.
of Saharan dust source activation obtained by regional modelling Toledano, C., Wiegner, M., Garhammer, M., Seefeldner, M.,
and satellite observations, Atmos. Chem. Phys., 13, 2381–2390, Gasteiger, J., Müller, D., and Koepke, P.: Spectral aerosol op-
doi:10.5194/acp-13-2381-2013, 2013. tical depth characterization of desert dust during SAMUM 2006,
Tesche, M., Wandinger, U., Ansmann, A., Althausen, D., Tellus B, 61, 216–228, doi:10.1111/j.1600-0889.2008.00382.x,
Müller, D., and Omar, A. H.: Ground-based validation of 2009.
CALIPSO observations of dust and smoke in the Cape Winker, D. M., Vaughan, M. A., Omar, A., Hu, Y., Pow-
Verde region, J. Geophys. Res. Atmos., 118, 2889–2902, ell, K. A., Liu, Z., Hunt, W. H., and Young, S. A.:
doi:10.1002/jgrd.50248, 2013. Overview of the CALIPSO mission and CALIOP data pro-
cessing algorithms, J. Atmos. Ocean. Tech., 26, 2310–2323,
doi:10.1175/2009JTECHA1281.1, 2009.
Mateos, D., Cachorro, V. E., Toledano, C., Burgos, M. A., Bennouna, Y., Torres, B.,
Fuertes, D., González, R., Guirado, C., Calle, A., and de Frutos, A. M.: Columnar
and surface aerosol load over the Iberian Peninsula establishing annual cycles,
trends, and relationships in five geographical sectors, Sci. Tot. Env., 518, 378–
392,doi:10.1016/j.scitotenv.2015.03.002, 2015.
Schneider, M., González, Y., Dyroff, C., Christner, E., Wiegele, A., Barthlott, S.,
García, O. E., Sepúlveda, E., Hase, F., Andrey, J., Blumenstock, T., Guirado, C.,
Ramos, R., and Rodríguez, S.: Empirical validation and proof of added value of
MUSICA's tropospheric δD remote sensing products, Atmos. Meas. Tech., 8, 483-
503, doi:10.5194/amt-8-483-2015, 2015.
Barreto, A., Cuevas, E., Pallé, P., Romero, P. M., Guirado, C., Wehrli, C. J., and
Almansa, F.: Recovering long-term aerosol optical depth series (1976–2012) from
an astronomical potassium-based resonance scattering spectrometer, Atmos. Meas.
Tech., 7, 4103-4116, doi:10.5194/amt-7-4103-2014, 2014.
Guirado, C., Cuevas, E., Cachorro, V. E., Toledano, C., Alonso-Pérez, S., Bustos, J. J.,
Basart, S., Romero, P. M., Camino, C., Mimouni, M., Zeudmi, L., Goloub, P.,
Baldasano, J. M., and de Frutos, A. M.: Aerosol characterization at the Saharan
AERONET site Tamanrasset, Atmos. Chem. Phys., 14, 11753-11773,
doi:10.5194/acp-14-11753-2014, 2014.
Barreto, A., Cuevas, E., Damiri, B., Guirado, C., Berkoff, T., Berjón, A. J., Hernández,
Y., Almansa, F., and Gil, M.: A new method for nocturnal aerosol measurements
with a lunar photometer prototype, Atmos. Meas. Tech., 6, 585-598,
doi:10.5194/amt-6-585-2013, 2013.
Guirado, C., Ramos, R., de Frutos, A., Berjón, A., Redondas, A., López, C., Cachorro,
V., Cuevas, E., González, R., González, S., and Hernández, M.: Optical calibration
facility at the Izaña Atmospheric Research Center, Óptica Pura y Aplicada, 45(1),
57-62, 2012.
Guirado C., Cuevas, E., Cachorro, V. E., Mimouni, M., Zeudmi, L., Toledano, C.,
Alonso-Pérez, S., Basart, S., Blarel, L., Goloub, P., and Baldasano, J. M.:
Preliminary characterization of columnar aerosols properties (AOD-AE) at the
Saharan Tamanrasset (Algeria) station, Óptica Pura y Aplicada, 44(4), 635-639,
2011.
Toledano, C., Cachorro, V. E., Berjón, A., de Frutos, A. M., Fuertes, D., González, R.,
Torres, B., Rodrigo, R., Bennouna, Y., Martín, L., and C. Guirado: RIMA-
AERONET network: long-term monitoring of aerosol properties, Óptica Pura y
Aplicada, 44(4), 629-633, 2011.
209
Anexo 4
Experiencia en formación
En el marco del proyecto “Sand and Dust Storm Warning Advisory and Assessment
System (SDS-WAS)-África” de la Organización Meteorológica Mundial se han
realizado las siguientes actividades formativas:
• Tres cursos de 15 horas de duración sobre instalación, operación, mantenimiento,
detección y resolución de problemas, y procesado de datos del fotómetro solar
Cimel, impartidos en el CIAI a los operadores de los fotómetros solares Cimel de
las estaciones AERONET de Túnez, El Cairo y Tamanrasset entre septiembre y
diciembre de 2014. Estos cursos forman parte, además, de las actividades
213
realizadas en el OAI como Campo de pruebas CIMO-OMM de instrumentos de
teledetección de aerosoles y vapor de agua.
• Curso de 12 horas de duración sobre instalación, operación, mantenimiento,
detección y resolución de problemas, y procesado de datos del fotómetro solar
Cimel, dentro del “WMO Course on Cimel sunphotometers for Northern Africa
stations within the SDS-Africa Project” impartido en el OAI en junio de 2013.
• Curso intensivo de instalación, operación, mantenimiento, resolución de
problemas y procesado de datos del fotómetro solar Cimel, impartido al personal
de la División Técnica y de Equipamiento de la Dirección Meteorológica de
Marruecos en la estación de Ouarzazate en marzo de 2013.
• Curso intensivo de instalación, operación, mantenimiento y resolución de
problemas del fotómetro solar Cimel, impartido al personal del Departamento de
Calidad del Aire de la Autoridad Meteorológica Egipcia en la estación de El Cairo
en mayo de 2010.
214
Anexo 5
En la presente memoria los acrónimos que son más reconocidos por sus siglas en
lengua inglesa se han mantenido como tales:
ACTRIS Aerosols, Clouds, and Trace gases Research InfraStructure Network
AEMET: Agencia Estatal de Meteorología
AERONET Aerosol Robotic Network (GSFC)
AMMA African Monsoon Multidisciplinary Analysis: Análisis Multidisciplinar del
Monzón Africano
ANSI American National Standards Institute
AOD Aerosol Optical Depth: Espesor Óptico de Aerosoles
BSRN Baseline Surface Radiation Network
CALIOP Cloud-Aerosol Lidar with Orthogonal Polarization
CALIPSO Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation
CE Cimel Electronique
CIAI: Centro de Investigación Atmosférica de Izaña (AEMET)
CIMO Commission for Instruments and Methods of Observations: Comisión para
Instrumentos y Métodos de Observación (OMM)
CLC: Capa límite convectiva
CNES Centre National d'Etudes Spatiales
CWT Concentration Weighted Trajectory
DMSP Defense Meteorological Satellite Program
DSCR Digital Sky Color Radiometer
EA: Exponente de Ångström
ECMWF European Centre for Medium-range Weather Forecasts
EOS Earth Observing System
EUMETSAT European Organisation for the Exploitation of Meteorological Satellites
FMF Fine Mode Fraction: Fracción del modo fino"
FOV Field of View: Campo de iluminación media
FT Free Troposphere: Troposfera libre
GCOS-GUAN Global Climate Observing System Upper-Air Network
GDAS Global Data Assimilation System
GOA: Grupo de Óptica Atmosférica (UVA)
GSFC Goddard Space Flight Center (NASA)
HYSPLIT Hybrid Single-Particle Lagrangian Integrated Trajectory
217
INTA: Instituto Nacional de Técnica Aeroespacial
IPCC Intergovernamental Panel for Climate Change: Panel Intergubernamental del
Cambio Climático
ITD Inter-Tropical Discontinuity: Discontinuidad Intertropical
IZO: Observatorio Atmosférico de Izaña (CIAI)
LIDAR Light Detection and Ranging
LOA Laboratoire d' Optique Atmosphérique (Universidad de Lille)
MBL Marine Boundary Layer: Capa límite marina
MISR Multi-angle Imaging SpectroRadiometer
MODIS Moderate resolution Imaging Spectrometer
MPLNET Micro Pulse Lidar Network (NASA)
MSG-SEVIRI Meteosat Second Generation-Spinning Enhanced Visible and Infrared
Imager
NASA National Aeronautics and Space Administration
NCEP National Center for Environmental Prediction
NDACC Network for the Detection of Atmospheric Composition Change: Red para la
detección del cambio en la composición atmosférica (NOAA)
NIST National Institute of Standards and Technology
NMMB/BSC Non-hydrostatic Multiscale Model/Barcelona Supercomputing Center
NOAA National Oceanic and Atmospheric Administration
OLS Operational Linescan System
OMM: Organización Meteorológica Mundial
ONM Office National de la Météorologie: Oficina Nacional de Meteorología argelina
PFR Precision Filter Radiometers: Radiómetros de Filtros de Precisión
PHOTONS Photométrie pour le Traitement Opérationnel de Normalisation Satellitaire
(Universidad de Lille)
PTB Physikalisch-Technische Bundesanstalt
PWV Precipitable Water Vapor: Vapor de Agua Precipitable
QA/QC Quality Assurance/Quality Control: Control de aseguramiento y control de
calidad
RBCC-E Regional Brewer Calibration Centre - Europe: Centro de Calibración
Regional Brewer para Europa
RGB Red, Green, Blue: Rojo, Verde y Azul
RIMA: Red Ibérica de Medida fotométrica de Aerosoles
SAL Saharan Air Layer: Capa de aire sahariano
SALAM Saharan Air Layer Analysis and Monitoring
SCMs: Sistemas Convectivos Mesoescalares
218
SCO: Observatorio de Santa Cruz de Tenerife (CIAI)
SDA Spectral Deconvolution Algorithm: Algoritmo de deconvolución espectral
SDS-WAS Sand and Dust Storm – Warning Advisory and Assessment System: Sistema
de Evaluación y Avisos de Tormentas de Polvo y Arena (OMM)
SENAR Subtropical Eastern North Atlantic Region: Región subtropical del Atlántico
nororiental
SSA Single Scattering Albedo: Albedo de dispersión simple
TIL Temperature Inversion Layer: Capa de inversión térmica
UTC Universal Time Coordinated: Tiempo Universal Coordinado
UV: Ultravioleta
UVA: Universidad de Valladolid
VAG: Vigilancia Atmosférica Global (OMM)
WRC World Radiation Centre: Centro de Radiación Mundial
219
Lista de símbolos
220