Source: http://slideplayer.es/slide/106816/
Timestamp: 2017-03-27 11:28:31+00:00

Document:
M. Pilar Utrillas Esteban Universitat de Valencia, ppt descargar
M. Pilar Utrillas Esteban Universitat de Valencia, 2002
Publicada porBelén Mascorro
Presentación del tema: "M. Pilar Utrillas Esteban Universitat de Valencia, 2002"— Transcripción de la presentación:
M. Pilar Utrillas Esteban Universitat de Valencia, 2002Características de aerosoles y nubes, de interés para la resolución de la Ecuación de Transferencia Radiativa José A. Martínez Lozano M. Pilar Utrillas Esteban Universitat de Valencia, 2002
Estructura Fundamentos teóricos Caracterización de aerosolesCaracterización de nubes Algunos resultados de campañas experimentales
Fundamentos Teóricos FUNDAMENTOS TEÓRICOS
Dispersión de la radiación por partículas atmosféricas Fundamentos Teóricos Fundamentos Teóricos Dispersión de la radiación por partículas atmosféricas Transferencia radiativa en una atmósfera dispersiva
Fundamentos Teóricos. DispersiónExtinción de la radiación Absorción : dL= - sa Ldx ( sa: coeficiente de absorción) Dispersión: dL= - ss Ldx (ss: coeficiente de dispersión) Espesor óptico :
Fundamentos Teóricos. DispersiónDependencia angular de la dispersión f(q): función angular de dispersión Función de fase:	Normalización:
Fundamentos Teóricos. DispersiónDispersión en la atmósfera Moléculas Partículas (aerosoles) Aerosol:	Toda partícula sólida o líquida, de tamaño mucho mayor que una molécula, en suspensión en la atmósfera, que sigue el movimiento de las masas de aire.
Fundamentos Teóricos. DispersiónClasificación de los constituyentes atmosféricos en cuanto a la dispersión Parámetro de tamaño: x = 2pr/l x < 0.6p/n dispersión por moléculas 0.6p/n < x < dispersión por partículas 5 < x reflexión
Fundamentos Teóricos. DispersiónDispersión molecular Simetría en la función de fase Fuerte dependencia con la longitud de onda (l-4) Formulación teórica: Rayleigh
Fundamentos Teóricos. DispersiónDispersión por partículas Asimetría de la función de fase. La dispersión hacia delante aumenta con el tamaño de las partículas. Prácticamente independiente de la longitud de onda Dependencia muy compleja respecto al ángulo de observación Formulación teórica: Mie x=1 x=10 x=40
Fundamentos Teóricos. DispersiónTeoría de Mie. Bases Físicas Interacción de una OEM primaria con una partícula que se comporta como un multipolo emitiendo una OEM secundaria. Desfases espaciales y temporales entre las ondas secundarias parciales emitidas, que originan interferencias. Estas interferencias dependen de: longitud de onda de la OEM incidente (l) tamaño e índice de refracción de la partícula (r,n) dirección angular de observación (q)
Fundamentos Teóricos. DispersiónTeoría de Mie. Formulación matemática La función angular de dispersión puede expresarse como i1 e i2 : funciones de distribución de intensidades de luz polarizada en las direcciones perpendicular y paralela al plano de observación. Para una partícula esférica e isótropa
Fundamentos Teóricos. DispersiónFormulación matemática de la teoría de Mie Yn y en funciones de Riccati-Bessel: Yn(mx)= z Jn (z) en (x) = z Jn (z) - i zYn (z) Jn y Yn funciones esféricas de Bessel.
Fundamentos Teóricos. DispersiónParámetros de la dispersión de Mie Factor de eficiencia de extinción de Mie Factor de eficiencia de dispersión de Mie Truncamiento de la serie: N = 2x + 3
Fundamentos Teóricos. DispersiónParámetros de la dispersión de Mie Coeficiente de dispersión Una partícula Dispersión:	sspart = pr2 Qdis(r,l,n) Extinción:	separt = pr2 Qext(r,l,n)
Fundamentos Teóricos. DispersiónParámetros de la dispersión de Mie Coeficiente de dispersión N partículas (coeficiente de volumen) n(r) función de distribución de tamaños Dispersión: Extinción:
Fundamentos Teóricos. DispersiónParámetros de la dispersión de Mie Espesor óptico. Dirección cenital. Dirección distinta a la cenital k al * = kal mal mal : masa óptica relativa
Fundamentos Teóricos. DispersiónParámetros de la dispersión de Mie Función de fase Parámetro de asimetría Representa la fracción de radiación dispersada hacia delante 19
Fundamentos Teóricos. DispersiónParámetros de la dispersión de Mie Albedo de dispersión simple Indica la capacidad de absorción de los aerosoles.
Fundamentos Teóricos. Transferencia RadiativaDispersión de la radiación por partículas atmosféricas Transferencia radiativa en una atmósfera dispersiva
Fundamentos Teóricos. Transferencia RadiativaLa Ecuación de Transferencia Radiativa (ETR) Forma general de la ecuación general de transferencia radiativa:
Fundamentos Teóricos. Transferencia RadiativaLa Ecuación de Transferencia Radiativa Función fuente: Función fuente de dispersión Función fuente de emisión
Fundamentos Teóricos. Transferencia RadiativaLa Ecuación de Transferencia Radiativa Aproximación de atmósfera plano-paralela: Las variaciones horizontales de las variables implicadas son despreciables en comparación con sus variaciones verticales. ETR en coordenadas (z,m,f) (m=cosq)
Fundamentos Teóricos. Transferencia RadiativaLa Ecuación de Transferencia Radiativa Radiación terrestre (rango espectral de onda larga) Atmósfera absortiva-emisiva (no dispersiva) Radiación solar (rango espectral de onda corta) Atmósfera dispersiva no emisiva
Fundamentos Teóricos. Transferencia RadiativaRadiación solar (rango espectral de onda corta). Atmósfera dispersiva no emisiva. Función fuente:
Fundamentos Teóricos. Transferencia RadiativaRadiación solar (rango espectral de onda corta). Atmósfera dispersiva no emisiva Considerando la radiancia global como la suma de las radiancias directa y difusa, e introduciendo el espesor óptico (kl) en lugar de la altura (z):
Fundamentos Teóricos. Transferencia RadiativaParámetros susceptibles de ser determinados experimentalmente para introducir en la ETR Espesor óptico: kl Albedo de dispersión simple: wl (kl) Función de fase Pl (kl) Parámetro imprescindible en cualquier caso: Espesores ópticos de los distintos constituyentes atmosféricos
Caracterización aerosolesCARACTERIZACIÓN DE AEROSOLES
Caracterización aerosoles- Parámetros físicos y radiativos - Distribución de tamaños - Parametrización del espesor óptico - Caracterización de aerosoles a partir de medidas experimentales
Caracterización aerosolesParámetros para determinar su influencia en la ETR Parámetros Físicos: Forma Tamaño Cantidad total Índice de refracción Parámetros Radiativos: Coeficientes de extinción Espesor óptico Función de fase Albedo de dispersión simple
Caracterización aerosolesParámetros Físicos Forma Normalmente se considera esférica, lo que permite aplicar la Teoría de Mie Esta aproximación es válida para partículas líquidas y partículas sólidas higroscópicas: sulfatos, sales marinas, etc. No es válida, por ejemplo, para aerosoles de tipo desértico ni cristales de hielo.
Caracterización aerosolesParámetros Físicos Tamaño Aerosoles esféricos: se caracterizan por su radio (r) Imposible estudiar un único aerosol en la atmósfera: se estudia el conjunto de todos los aerosoles en una columna a lo largo de la atmósfera, por lo que el radio no es único. Es necesario definir una función de distribución de tamaños: n(r) dr Representa el número de partículas por unidad de volumen con radios comprendidos entre r y r+dr
Caracterización aerosolesParámetros Físicos Cantidad total Se representa por el número total de aerosoles por volumen unidad de aire: N De forma que n(r) = dN(r)/dr (función de distribución de tamaños) Dimensiones	[N(r)] = L-3 (partículas por unidad de volumen) [n(r)] = L-4
Caracterización aerosolesParámetros Físicos Índice de refracción n = nR - inI ( parte imaginaria distinta de cero: aerosoles absorbentes) Parámetros físicos más relevantes Función de distribución de tamaños Todos los aerosoles que presentan los mismos valores de índice de refracción y distribución de tamaños son iguales a efectos de su comportamiento radiativo, e imposibles de distinguir a partir de medidas de radiación a nivel de suelo.
Caracterización aerosolesParámetros radiativos Los relativos a la teoría de Mie Todos ellos son función de: Parámetro de tamaño: x = 2pr/l Indice de refracción: n = nR - inI Por tanto son funciones de tres variables (r, l,n) En algunos casos presentan además dependencia angular respecto a la dirección de observación (q).
Caracterización aerosolesParámetros radiativos Factor de extinción: Qext(r, l,n) Su uso es debido a que normalmente para los aerosoles no es posible distinguir (experimentalmente) entre los procesos de: Dispersión Qs(r, l,n) Absorción Qa(r, l,n)
Caracterización aerosolesParámetros radiativos Factor de extinción: Qext(r, l,n) (para un única partícula) A partir de los factores de extinción se definen, para una partícula, las secciones eficaces (o coeficientes moleculares): Dispersión Ps(l) = pr2 Qs(r, l,n) Absorción Pa(l) = pr2 Qa(r, l,n) Extinción Pext (l) = pr2 Qext (r, l,n) Para simplificar la notación suele emplearse simplemente: Pa(l) = pr2 Q(r, l,n) Pa(l): sección eficaz de extinción para aerosoles Q: coeficiente de extinción de Mie para aerosoles
Caracterización aerosolesParámetros radiativos Factor de extinción para un conjunto de partículas aerosoles Distribución de tamaños Coeficiente de extinción de volumen:
Caracterización aerosolesParámetros radiativos Espesor óptico Dirección cenital: donde representa la función de distribución de tamaños en una columna vertical.
Caracterización aerosolesParámetros radiativos Espesor óptico Si en lugar de la dirección vertical consideramos cualquier otra dirección: ka(l)* = ka(l) ma(l) siendo ma(l) la masa óptica relativa de aerosoles. El espesor óptico de aerosoles es el parámetro más simple de obtener experimentalmente.
Caracterización aerosolesParámetros radiativos Parámetro de asimetría Representa la fracción de dispersión hacia delante respecto a la dispersión total donde g(r, l,n) es el parámetro de asimetría para una partícula
Caracterización aerosolesParámetros radiativos Función de fase Proporciona la distribución angular de la radiación dispersada por el conjunto de partículas donde p(r, l,n,q) es la función de fase para una partícula
Caracterización aerosolesParámetros radiativos Albedo de dispersión simple (single scattering albedo) Indica la capacidad de absorción de los aerosoles. Para aerosoles no absorbentes wO(l) = 1 Es fundamental para estudiar: - forzamiento radiativo de los aerosoles - posible impacto sobre le cambio climático. Actualmente es el parámetro óptico cuya determinación experimental presenta mayores incertidumbres (del orden del 50%).
Caracterización aerosolesParámetros radiativos Consideración final ka(l) depende de Qext(r, l,n) g(l) y P(q,l) dependen de Qs(r, l,n) wO(l) depende de Qs(r, l,n) y Qext(r, l,n) El espesor óptico de aerosoles puede determinarse a partir de medidas de extinción de la radiación (medidas de irradiancia directa). El resto de parámetros radiativos requieren medidas de irradiancia difusa procedente del cielo (radiancia) en diferentes ángulos de observación.
Caracterización aerosoles. Distribución de tamañosActualmente se admite que existen 3 tipos de aerosoles. Esta clasificación no se realiza tanto en base a su tamaño como a su origen: Partículas Aitken (nuclei mode) Partículas finas (accumulation mode) Partículas gruesas (coarse particle mode) Partículas Aitken: Tamaño muy pequeño. r < 0.05 mm. Originadas en la propia atmósfera por nucleaciones moleculares homogéneas. Partículas finas: Tamaño del orden de r < 1 mm. Se originan por condensación y acumulación (antiguo smoke). Partículas gruesas: Tamaño r > 1 mm. Origen mecánico de fragmentación (antiguo dust).
Caracterización aerosoles. Distribución de tamaños¿Cómo establecer una distribución de tamaños? ¿Se pueden englobar los tres tipos en una única distribución? ¿Es necesario considerar tres distribuciones independientes?. Preguntas a responder A) ¿Es necesario considerar dos modos de acumulación? B) ¿Es necesario considerar una discontinuidad entre los modos de acumulación y el de origen mecánico?
Caracterización aerosoles. Distribución de tamañosRespuesta pregunta A La respuesta debe buscarse en base a las ecuaciones dinámicas de nucleación-acumulación. Se demuestra que el tiempo necesario para que los núcleos pasen al modo de acumulación mediante procesos de auto-coagulación es mucho mayor que el tiempo de residencia media en la atmósfera. Respuesta pregunta B A partir del análisis teórico de la velocidad de aumento de volumen, se observa que ésta decrece rápidamente por encima de los 0.5 mm. Para radios de aproximadamente 1 mm es del orden de 10-9 mm3/cm3s. Si queremos hacer crecer aerosoles de tipo fino hasta los tamaños de tipo grueso a estas velocidades: varias semanas. Supera con creces los tiempos de residencia media en la atmósfera.
Caracterización aerosoles. Distribución de tamañosCondiciones que deberá cumplir una función de distribución de tamaños que sea válida para los diferentes tipos de aerosoles presentes en la atmósfera: a) Permitir ajustar distribuciones de tamaño en un rango amplio: de mm a 100 mm. b) Ser independiente del procedimiento seguido en la normalización. Se debe poder ajustar número, superficie, volumen o masa de aerosoles. c) Estar establecida sobre bases físicas.
Caracterización aerosoles. Distribución de tamañosSe puede demostrar que si en una distribución, la variable a distribuir está afectada por mecanismos de crecimiento, en los cuales el efecto de los mismos es proporcional al valor de la variable, la distribución debe ser log-normal. En base a los razonamientos anteriores, lo más adecuado será utilizar una función multimodal (con tres modos) que sea la superposición de tres log-normales: Distribución lognormal trimodal
Caracterización aerosoles. Distribución de tamañosReflexiones finales Si esto parece tan obvio - ¿porqué se siguen utilizando otras funciones de distribución como la potencial (Junge) o la distribución gamma? - ¿porqué se siguen proponiendo distribuciones monomodales?
Caracterización aerosoles. Distribución de tamañosUn enfoque inicial poco afortunado de repercusiones históricas Primera función de distribución propuesta (Junge, 1953): n(r) = C r-n-1 En forma más general n(r) = A rB	(función potencial) A partir de ella dN = A rB dr = A rB+1 dr/r = A rB+1 dlnr dN/ dlnr = A rB+1 ln (dN/ dlnr) = lnA + (B+1) lnr
Caracterización aerosoles. Distribución de tamañosSi para un distribución de partículas obtengo una gráfica de la forma de la figura, aparentemente es una distribución potencial monomodal. ln (dN/ dlnr) = lnA + (B+1) lnr
Caracterización aerosoles. Distribución de tamañosAhora bien, si en lugar de la representación de la forma anterior (en N) elijo una en volumen dV/dlnr = f(lnr) Aparecen claramente dos modos
Caracterización aerosoles. Distribución de tamañosrM: radio de la moda (radio para el cual la función adquiere el valor máximo) <r>: radio medio rEF: radio efectivo El radio efectivo se introduce debido a que en la dispersión de Mie todos los parámetros son proporcionales a pr2
Caracterización aerosoles. Distribución de tamañosArea eficaz de la distribución: Volumen eficaz de la distribución: Importancia del radio efectivo	rEF = (3/4) V/A Permite estudiar, con limitaciones empíricas, problemas de dispersión múltiple en términos de dispersión simple equivalente.
Caracterización aerosoles. Parametrización del espesor ópticoAngstrom (1929) kal = bl-a (l en mm) Sentido físico b = tal (l = 1 mm). Proporciona una medida de la cantidad de aerosoles presentes en la atmósfera y se conoce como coeficiente de turbiedad. a está relacionado con la distribución de tamaño de los aerosoles. De acuerdo con Angstron b y a son independientes de la longitud de onda.
Caracterización aerosoles. Parametrización del espesor ópticoDado que estos coeficientes no dependen de l, es suficiente disponer de medidas en un par de longitudes de onda para determinarlos. La experiencia demuestra que si se toman únicamente dos longitudes de onda los errores que se cometen son muy elevados, y además los resultados dependen del par de longitudes de onda elegido. Es mucho más preciso utilizar varias longitudes de onda, realizando un ajuste a partir de las medidas disponibles en todas las longitudes de onda ln kal = - a lnl + ln b
Caracterización aerosoles. Parametrización del espesor óptico
Caracterización aerosoles. Parametrización del espesor ópticoAngstrom se basa en la distribución de tamaños de Junge, y esta función de distribución no es válida para cualquier tipo de aerosoles ni para todo el rango espectral (actualmente se da por valida la lognormal multimodal). Algunos autores consideran que bajo determinadas condiciones de turbiedad Angstrom no es válida porque: no tiene en cuenta la absorción por aerosoles la mayoría de las distribuciones de tamaño son multimodales las partículas gigantes tienden a dispersar la luz de manera independiente respecto a la longitud de onda. No se ha propuesto hasta ahora una alternativa válida para la fórmula de Angstrom.
Caracterización aerosoles. Parametrización del espesor ópticoInterés de a Actualmente se intenta utilizar como un parámetro simple para identificar aerosoles a partir de las funciones de distribución de tamaños. Su interés radica en que es un parámetro relativamente fácil de determinar desde satélite, particularmente con los modernos sensores hiperespectrales.
Caracterización aerosoles. Métodos experimentalesCaracterización de aerosoles a partir de medidas experimentales Técnicas actuales: Contadores de partículas Lidar Medidas de radiación solar a nivel de suelo
Caracterización aerosoles. Métodos experimentalesCaracterización de aerosoles a partir de medidas experimentales Medidas espectrales de radiación solar a nivel de suelo Irradiancia directa (extinción) Permiten determinar el espesor óptico de aerosoles espectral. A partir de él se puede obtener, mediante métodos de inversión, la función de distribución de tamaños. Irradiancia difusa direccional (radiancia de cielo) Permite determinar el espesor óptico de dispersión y los coeficientes de dispersión de Mie. Sirven como datos de entrada en algoritmos para el cálculo de: Parámetro de asimetría Función de fase Albedo simple de aerosoles
Caracterización aerosoles. Métodos experimentalesCaracterización de aerosoles a partir de medidas experimentales Medidas espectrales de radiación solar a nivel de suelo Instrumentación Obtención del espesor óptico espectral Obtención de la distribución de tamaños
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiación solar directa Espectroradiómetro provisto de un tubo limitador de radiancia (colimador) Los colimadores usados habitualmente tienen un ángulo de visión (Field of view, FOV) no superior a 5º, de forma que no sea necesaria la corrección de la radiación circumsolar. Instrumentación disponible en Valencia: Espectroradiómetro Licor 1800 Espectroradiómetro Optronic 754
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiación solar directa Espectrorradiómetro Li-cor 1800 Simple monocromador con red de difracción holográficas de 800 líneas/mm Rango de medida nm Resolución 6 nm Tubo limitador de radiancia: 4.7° de abertura de campo de visión Calibración: Li-Cor Optical Radiation Calibrator
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiación solar directa Espectrorradiómetro Li-cor 1800 Microcomputador interno Coseno receptor Monocromador Filtro rotatorio Detector de silicio
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiación solar directa Espectrorradiómetro Li-cor 1800 0.00 0.20 0.40 0.60 0.80 1.00 1.20 1.40 300 400 500 600 700 800 900 1000 1100 Global Directa Difusa
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiación solar directa Espectrorradiómetro Li-cor 1800
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiación solar directa Espectrorradiómetro Optronic 754 Rango: nm Ancho de banda: 2 nm Mínimo paso banda: 0.05 nm Doble monocromador: 1200 lineas/mm Receptor: esfera integradora ( PTFE) Detector: fotodiodo silicio estabilizado en temperatura por Peltier
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiancia de cielo Optronic + Telescopio	Cimel
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas de radiancia espectral en varias longitudes de onda para diferentes ángulos acimutales en el plano del almucántar del sol (Optronic) Valencia, día 30/05/00, masa óptica 1.05.
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiancia de cielo Cimel Tipos de medida - Radiancia de cielo Colimador de cielo + alta ganancia - Irradiancia directa de sol Colimador solar + baja ganancia - Aureola Colimador solar + alta ganancia
Métodos experimentales. InstrumentaciónCaracterización de aerosoles a partir de medidas experimentales Instrumentación Medidas espectrales de radiancia de cielo Cimel Medida en el Almucantar (ALM) Se ejecuta cada hora Consiste en 4 barridos independientes Filtros empleados: 440, 670, 870, 1020nm Duración: 60s aprox. cada barrido Protocolo masa óptica: m=4, 3, 2, 1.7 Paso angular variable con la posición
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Medidas espectrales de radiación solar a nivel de suelo Instrumentación Obtención del espesor óptico espectral Obtención de la distribución de tamaños
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Determinación del espesor óptico Medidas de irradiancia directarealizadas con el Licor 1800
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Determinación del espesor óptico La atenuación de la radiación a través de un medio puede evaluarse a partir de la ley de Bouger-Lambert-Beer. Aplicada a una radiación monocromática a través de la atmósfera toma la forma Inl = (1/r2)I0l exp(- klm) Inl: Irradiancia medida a nivel del suelo I0l: Irradiancia en el límite de la atmósfera r : radio vector relativo Tierra-Sol kl : espesor óptico de la atmósfera (total) m : masa óptica relativa
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Determinación del espesor óptico A partir de los valores experimentales de Inl, y de la curva espectral de I0l, es posible obtener el espesor óptico total de la atmósfera, kl. El espesor óptico total puede expresarse como una suma de las contribuciones de todos los procesos que tienen lugar en la atmósfera kl = kRl+ kGl+ kal kRl : espesor óptico de dispersión molecular (Rayleigh) kal : espesor óptico de extinción por aerosoles (engloba absorción y dispersión) kGl : espesor óptico de absorción por diferentes gases atmosféricos
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Determinación del espesor óptico El término correspondiente a los gases atmosféricos es habitual separarlo en tres contribuciones debidas a los espesores óptcos de absorción de: ozono vapor de agua mezcla homogénea de gases Para evitar en lo posible el aplicar correcciones para la absorción por todos los gases presentes en la atmósfera se eligen bandas espectrales en las que la absorción sea mínima. En la ventana del visible ( nm) – ozono (banda de Chapuis) NO2 (solo importante en zonas urbanas altamente contaminadas) En estas condiciones podemos expresar kal = kl - kRl - kol
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Determinación del espesor óptico La dispersión de Rayleigh es un fenómeno perfectamente conocido. Recientemente se ha propuesto una aproximación que mejora los resultados cuando la transmisividad de la atmósfera es muy alta (días muy claros) y Rayleigh adquiere una mayor importancia relativa. (Bodahine et al., On Rayleigh optical depth calculations, J. Atmos. Oceanic Technol., 16, , 1999)
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Determinación del espesor óptico Para aplicar la corrección debida al ozono es necesario conocer la cantidad de este elemento presente en una columna atmosférica. Alternativas: Ozonosondeo Absorción diferencial (Brewer) TOMS
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Determinación del espesor óptico Medida del Ozono Método seguido en Valencia: Microtops II
Métodos experimentales. Espesor ópticoCaracterización de aerosoles a partir de medidas experimentales Determinación del espesor óptico Curvas de espesor óptico de aerosoles obtenidas por el método descrito
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Medidas espectrales de radiación solar a nivel de suelo Instrumentación Obtención del espesor óptico espectral Obtención de la distribución de tamaños
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Espesor óptico de aerosoles en la teoría de Mie que se puede expresar como donde nC(r) es la distribución de tamaños en una columna vertical de aire.
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Determinar nC(r) a partir de ka(l,n), teniendo en cuenta que se dispone de un valor de ka para cada una de las longitudes de onda del intervalo espectral de medidas. La solución de este tipo de ecuaciones se conoce como problema inverso, y normalmente la solución no es única (problemas "ill-posed“) La solución depende del método utilizado y de las condiciones iniciales impuestas . En este tipo de ecuaciones: - nC(r): función desconocida a determinar - Q(r, l,n): kernel de la función
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Métodos utilizados en Valencia. - King (1978, 1982) Medidas de extinción de irradiancia solar directa - Nakajima (1985 y siguientes) Medidas de extinción de irradiancia solar directa y de difusa de cielo
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Algoritmo de King Se utiliza una aproximación numérica. La integral se sustituye por un sumatorio sobre intervalos del radio, cada uno de los cuales a su vez está formado por varios subintervalos. La función n(r) se sustituye por el producto de dos funciones n(r) = h(r) f(r) donde h(r) varia rápidamente con r, mientras que f(r) es prácticamente constante en los distintos intervalos de r, con lo que la integral toma la forma
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Algoritmo de King Como f(r) se considera constante en cada uno de los intervalos discretos definidos, es posible escribir la integral como un sistema de ecuaciones lineales de la forma g = Af + e gi = kal e es un vector cuyas componentes ei representan las desviaciones entre los valores experimentales gi y los valores estimados Aij fi.
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Algoritmo de King Para evitar las inestabilidades que aparecen en la solución cuando se intenta resolver la ecuación directamente, se introduce una condición de ligadura. Existen distintas posibles condiciones de ligadura, aunque la más utilizada consiste en minimizar la segunda derivada de la curva. Las soluciones se obtienen por un proceso iterativo
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Algoritmo de King En el método de inversión los errores se calculan a través de la matriz de covariancia, dada por las expresiones S = s2(gH + AT A)-1 S = (gH + AT Se-1A)-1 Una vez obtenido el error que afecta a f, el error de n(r) se determina mediante propagación de errores.
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Algoritmo de King Algunos resultados para Valencia
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Algoritmo de King Inconvenientes del algoritmo Algunas inversiones divergen al aplicar las iteraciones El método es sensible al intervalo de radios al que se aplica (rMIN, rMAX) Para aplicar el método es necesario fijar el tamaño de los intervalos de radios. El método es sensible al número de intervalos de radio utilizados, y a si estos son o no regulares. El método se aplica a un número finito de longitudes de onda. Es sensible al número de longitudes de onda utilizadas y a como se distribuyen en el rango espectral de medidas experimentales.
Métodos experimentales. Distribución de tamañosCaracterización de aerosoles a partir de medidas experimentales Determinación de la función de distribución de tamaños Algoritmo de King Actualmente se está trabajando en determinar: número óptimo de intervalos de r. número óptimo de valores de l. distribución de estos valores en el espectro de r y l. La determinación de la bondad de los resultados se determina a partir de los valores reconstruidos del espesor óptico de aerosoles ka(l,n) EXP ® nC(r) ® ka(l,n) EST ka(l,n) EXP «ka(l,n) EST
Caracterización nubesCARACTERIZACIÓN DE NUBES
Caracterización nubesNubes y temperatura de equilibrio del sistema Tierra-Atmósfera Distribución de tamaños La resolución de la ETR en presencia de nubes Parametrización de las propiedades radiativas de las nubes
Caracterización nubesNubes y temperatura de equilibrio del sistema Tierra-Atmósfera Efectos de las nubes sobre el equilibrio radiativo Tierra-Atmósfera: Aumento del albedo del sistema Temperatura inferior a la del suelo Consecuencia: Disminución de la temperatura del sistema
Caracterización nubesNubes y temperatura de equilibrio del sistema Tierra-Atmósfera Situación actual: Cobertura nubosa 51 % Albedo planetario: 0,31-0,33 Temperatura del sistema: 253 K Situación hipotética: Cobertura 100 % (tipo Venus) Albedo planetario: 0,71 Temperatura del sistema: 205 K Situación hipotética: Cobertura 0 % (tipo Marte): Albedo planetario: 0,17 Temperatura del sistema: 267 K Efecto a nivel del suelo Cobertura nubosa: 51 %	Temperatura: 288 K Cobertura nubosa: 0 %	Temperatura: 304 K
Caracterización nubesNubes y temperatura de equilibrio del sistema Tierra-Atmósfera Efectos generales sobre la radiación Radiación solar (OC): Absorción nula Gran dispersión No es necesario reformular la ETR (atmósfera dispersiva no absortiva) Radiación terrestre (OL): Absorción en las nubes de gotas de agua Dispersión en las nubes de cristales de hielo Es necesario reformular la ETR (atmósfera absortiva y dispersiva) Microondas: Transparentes Rango adecuado para comunicaciones sin interferencias en la señal
Caracterización nubesDistribución de tamaños Nube: Conjunto de aerosoles monomodales Gotas de agua (esféricos) Cristales de hielo (no esféricos) Nubes de agua: caracterizadas por el contenido en agua líquida (LWC) Para una gota LWC = (4/3)pr3rw Para una nube completa 101
Caracterización nubesDistribución de tamaños Función n(r) utilizada habitualmente: G modificada Como ejemplo Para cúmulos: n(r) = 2,37 r6 exp(-1,5r) siendo rm = 4,0 mm
Caracterización nubesDistribución de tamaños Nubes de hielo: caracterizadas por el contenido en agua en forma de hielo (IWC) Este parámetro se puede obtener a partir de la temperatura de la nube a partir de la relación empírica: Ln (IWC) = -7,6 + exp [-0,244x10-3(t-20)2,45] Los cristales se caracterizan en tamaño por L: dimensión máxima del cristal (longitud)
Caracterización nubesDistribución de tamaños Distribución de tamaños nubes de hielo: n(L) donde B1 y B2 son ctes a determinar, IWC se expresa en g/m3 y n(L) se expresa en m-3/mm N100 y N son las concentraciones de cristales con L igual a 100 mm y 1000 mm
Caracterización nubesLa resolución de la ETR en presencia de nubes Implica incluir término de dispersión en la resolución de la ETR para onda larga Partiendo de la aproximación plano-paralela Se consideran separadamente: Nubes de agua (+ vapor de agua): absorción, sa Nubes de hielo: dispersión, ss
Caracterización nubesLa resolución de la ETR en presencia de nubes La ETR toma la forma Que se puede expresar como: Siendo se = sa + ss Introduciendo el albedo de dispersión simple, se tiene
Caracterización nubesLa resolución de la ETR en presencia de nubes La función fuente Jndis está asociada a la dispersión múltiple, y por lo tanto a la función de fase P(m,m’). En el rango IR es suficiente considerar la contribución independiente del azimut Para la función fuente Jnem se utiliza la aproximación habitual del cuerpo negro.
Caracterización nubesLa resolución de la ETR en presencia de nubes Para resolver la ETR se utiliza normalmente algún método de ordenadas discretas (DISORT o similares) Permite discretizar el problema en azimut, obteniéndose expresiones a resolver del tipo: La solución final es de la forma donde kj y fj son, respectivamente, los vectores y valores propios
Caracterización nubesLa resolución de la ETR en presencia de nubes Consideramos la nube rodeada por dos capas atmosféricas (encima y debajo de ella) que a su vez absorben y emiten radiación Podemos determinar la la contribución de la nube a la radiación en los límites superior e inferior de la atmósfera. Límite superior: Límite inferior (suelo):
Caracterización nubesParametrización de las propiedades radiativas de las nubes Consideramos de forma independiente los rangos: radiación solar (OC) radiación terrestre (OL) Además se consideran separadamente: nubes de gotas de agua nubes de cristales de hielo
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de agua Todas las características radiativas de la dispersión dependen (como en el caso de aerosoles) de la distribución de tamaños n(r) Se define el radio efectivo como: (promedio de la distribución de tamaños respecto a la sección geométrica eficaz)
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de agua Relación entre ref y LWC Se define LWP como el contenido de LWC en una columna vertical Se define el espesor óptico de la forma habitual
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de agua Dispersión de Mie Cuando (de hecho se puede considerar cte para r > 10 mm) Dividiendo LWP y k obtenemos A partir de LWP y k se puede desarrollar una parametrización sencilla de la ETR utilizando la aproximación de los dos flujos.
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de agua Alternativa a la resolución de la ETR Se pueden obtener propiedades radiativas de banda ancha con una precisión aceptable definiendo la reflectividad y la absortividad de la forma Reflectancia Los resultados muestran que depende significativamente del ref Ejemplo. LWP de 100 g/m2 Una nube con ref pequeño (del orden de 2 mm) muestra una reflectancia doble que una nube de ref de 16 mm
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de agua Absortancia Tiene comportamiento diferente según el valor del LWP. Para LWP > 100 g/m2 las gotas de mayor radio absorben más que las gotas de radio pequeño. Para LWP < 100 g/m2 las gotas de mayor radio absorben menos que las gotas de radio pequeño. Razón: las gotas grandes tienen un área de sección eficaz más pequeña que las gotas pequeñas.
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de hielo Parámetro de partida: Es necesario estudiar cada tipo de cristal separadamente. Uno de los más sencillos y abundantes: hexagonales Siendo L la longitud máxima del cristal y D longitud transversal del cristal (aproximadamente 2r)
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de hielo Para una columna vertical definimos IWP = Dz IWC Sección eficaz de extinción (cristales hexagonales) Tamaño (o grosor) efectivo medio de la distribución de cristales
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de hielo Espesor óptico (cristales hexagonales). Se demuestra que k = IWP (a + b/ De) siendo a y b coeficientes empíricos a determinar. Por tanto k = f(IWP, De) que a su vez son función de la temperatura
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de hielo Una posible parametrización de las nubes de hielo puede escribirse en la forma siendo ani coeficientes empíricos Como en el caso de las nubes de agua, esta parametrización permite resolver la ETR de forma sencilla a partir de la aproximación de los dos flujos.
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de hielo Alternativa a la resolución de la ETR Rflectancia y absortancia en banda ancha Reflectancia Para una De determinada, la reflectancia aumenta cuando aumenta IWP Para un IWP determinado, la reflectancia aumenta cuando disminuye el De (debido a que la sección eficaz geométrica es mayor)
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OC. Nubes de hielo Alternativa a la resolución de la ETR Rflectancia y absortancia en banda ancha Absortancia Comportamiento totalmente análogo al de las nubes de agua Ejemplo. Para un valor de IWP de 100 g/m2 Valores de R de 70, 55, 40 y 30% para valores de De de 25, 50, 75 y 100 mm, respectivamente.
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OL. Emisividad Para simplificar el problema se considera despreciable la contribución de la dispersión en este rango. Por lo tanto: atmósfera absortiva sin dispersión. Validez de la aproximación: Partículas de nubes pequeñas (hasta 10 mm) Longitud de onda incidente l > 10 mm
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OL. Emisividad Emisividad monocromática de la nube siendo sn el coeficiente de absorción de masa para partículas de nubes. El flujo debido a la emisividad viene dado por: siendo el factor de difusividad (se asume que la emisividad se aproxima a 1 cuando sn o w son muy grandes)
Caracterización nubesParametrización de las propiedades radiativas de las nubes Rango de OL. Emisividad El valor de sn varía relativamente poco con n. El espectro de absorción del agua líquida y sólida presenta un aspecto mucho más continuo que el del vapor de agua. Esto permite aplicar la aproximación del cuerpo gris:
Resultados experimentalesALGUNOS RESULTADOS DE CAMPAÑAS EXPERIMENTALES (Grupo de Radiación Solar de Valencia) AREA DE VALENCIA. VERANO 1998 AREA DE BARRAX. VERANOS 1998, 1999 Y 2000
Campaña: Area de ValenciaAREA DE VALENCIA. VERANO 1998 Objetivos de la campaña: Análisis de la posible influencia de los aerosoles urbanos en las áreas circundantes (limitación del radio de estudio: 50 km). Análisis de la variabilidad de los aerosoles en función del régimen de vientos propios de la zona. Diferencias entre los vientos de origen marítimo (Levante) y terrestre (Poniente).
Campaña: Area de ValenciaSitios de observación experimental: Centro urbano. Expuesto a los aerosoles emitidos en el área urbana y a los provenientes del mar. Parque natural de La Albufera. Localizado 20 km al sur del centro urbano. Muy próximo al mar, libre de edificios, y presumiblemente fuera de los límites de influencia del núcleo urbano, debido a los vientos provenientes del mar. Chulilla. Localizado a una distancia de 50 km en dirección Noroeste, a una altura sobre el nivel del mar de 600 m. Ningún tipo de actividad industrial. Economía basada en la agricultura y el turismo rural.
Campaña: Area de ValenciaLocalización de los puntos de medida
Campaña: Area de ValenciaExtensión temporal: Días comprendidos entre el 22 de Junio y el 10 de Julio de 1998. Medidas simultáneas en las tres localidades cada media hora, desde las 7:00 hasta las 16:00 hora solar. Sólo las medidas correspondientes a días totalmente despejados en las tres localidades se seleccionaron para el análisis de datos: 6 días.
Campaña: Area de ValenciaCondiciones sinópticas durante la campaña: Cinco de los días seleccionados: Condiciones típicas de verano en el área de Valencia. Combinación de altas presiones debidas al anticiclón de las Azores y una baja térmica relativa originada en el norte de África. Da lugar a vientos ligeros de componente este o sureste, con humedad relativa del 50%, que puede originar nieblas matinales y ocasionalmente formación de nubes de evolución diurna a últimas horas de la tarde. Uno de los días seleccionados: Viento de dirección noroeste. Se caracteriza básicamente por valores de humedad relativa muy bajos en relación a los valores medios anuales de la zona considerada.
Campaña: Area de ValenciaMedidas experimentales Irradiancia espectral solar directa en incidencia normal Instrumentación Espectrorradiómetros Li-cor 1800 provistos de colimadores 133
Campaña: Area de ValenciaIntercomparación de los tres instrumentos antes de la campaña Día 18/06/98 (masa óptica 1.4)
Campaña: Area de ValenciaDecisión previa: Establecer la posible influencia en las medidas de la ubicación del espectrorradiómetro a utilizar en el área urbana.
Campaña: Area de ValenciaUbicación urbana elegida: Terraza de la E.T.S.I. Industriales, de la UPV horizonte perfectamente libre de obstáculos disponibilidad de infraestructura
Campaña: Area de ValenciaResultados Correspondientes a dos días: 3 de Julio (viento terrestre), día que hemos catalogado de atípico 6 de Julio (viento marítimo), representativo de las condiciones típicas de verano en el área considerada En la Figura se observa que durante el día 3 de Julio la irradiancia directa integrada es alrededor de un 20% mayor que durante el día 6 de Julio.
Campaña: Area de ValenciaResultados. Espesor óptico Curvas espectrales de espesor óptico de aerosoles Espesor óptico en el visible para La Albufera Espesor óptico en el visible para Chulilla
Campaña: Area de ValenciaResultados. Espesor óptico Espesor óptico de aerosoles para 550 nm. Evolución diaria. Viento marítimo	Viento terrestre
Campaña: Area de ValenciaResultados. Espesor óptico El espesor óptico de aerosoles para viento marítimo es bastante más elevado que para viento terrestre. La forma de las curvas espectrales es muy diferente para ambos condiciones, con una pendiente mucho más pronunciada en condiciones de viento marítimo. Presencia de aerosoles de distinto tipo en ambos días.
Campaña: Area de ValenciaResultados. Coeficiente a de Angstrom. Ajuste directo de los valores experimentales del espesor óptico de aerosoles.
Campaña: Area de ValenciaResultados. Coeficiente a de Angstrom. Dependencia de los coeficientes de Angstrom con el rango de longitudes de onda considerado.
Campaña: Area de ValenciaResultados. Coeficiente a de Angstrom. En el rango superior de longitudes de onda las correlaciones son menores. El coeficiente de correlación decrece rápidamente a medida que el intervalo utilizado para obtener a incluye longitudes de onda mayores. Limitamos el análisis de las curvas espectrales al rango nm
Campaña: Area de ValenciaResultados. Coeficiente a de Angstrom. Conclusiones Los valores de a en condiciones de viento terrestre indican la presencia de aerosoles de menor tamaño. Este resultado conduce a una aparente contradicción. El presunto origen continental de los aerosoles presentes bajo estas condiciones debería conducir a resultados opuestos.
Campaña: Area de ValenciaResultados. Retrotrayectorias. Día 3 de julio (viento terrestre) Retrotrayectoria correspondiente a 700 HPa
Campaña: Area de ValenciaResultados. Retrotrayectorias. Día 6 de julio (viento marítimo) Retrotrayectoria correspondiente a 700 HPa
Campaña: Area de ValenciaResultados. Distribución de tamaños Día 3 de julio (viento terrestre)
Campaña: Area de ValenciaResultados. Distribución de tamaños Día 6 de julio (viento marítimo)
Campaña: Area de ValenciaResultados. Distribución de tamaños Carácter multimodal de las distribuciones Se pone claramente de manifiesto si en lugar de representar la función de probabilidad de densidad de partículas se representa la correspondiente a volumen.
Campaña: Area de ValenciaResultados. Distribución de tamaños Carácter multimodal de las distribuciones Las distribuciones responden a la forma lognormal Función de distribución lognormal multimodal. El ajuste de la misma permite determinar la posición de los modos en el espectro de los radios
Campaña: Area de ValenciaResultados. Distribución de tamaños Carácter multimodal de las distribuciones Chulilla. Día 6 de julio a la 12:00 N rM(mm) Modo 1 0.001 Modo2 0.087
BARRAX. VERANOS 1998, 1999 y 2000 Campañas DAISEXCampañas DAISEX (Digital Airborne Imaging Spectrometer Experiment) Financiadas por: ESA (European Space Agency) Programa marco: Land Surface Processes and Interactions Mission (LSPIM) Earth Observation Preparatory Program.
Campañas DAISEX Campañas DAISEX Área de estudio: Área agrícola de 10x10 km2 en Barrax (Albacete) Marco temporal: Días seleccionados de los veranos de 1998, 1999, 2000
Campañas DAISEX Resultados campaña Espesor óptico Espesor óptico de aerosoles. Banda nm. Coefientes de ajuste Día 04/06/99 a las 12:00 Rango a b r -0.01 0.11 0.01 1.05 0.05 0.67 -0.46 0.15 0.43
Campañas DAISEX Resultados campaña Retrotrayectoria 700 HPa
Campañas DAISEX DAISEX 1999 Test del método gracias a la disponibilidad del avión ARAT Marco: Programa STAAARTE II (CEE)
Campañas DAISEX DAISEX 1999 Test del método gracias a la disponibilidad del avión ARAT Comparación con medidas en vuelo de: Lidar Nefelómetro
Campañas DAISEX DAISEX 1999 LIDAR Señal retrodifundida por las moléculas y aerosoles de la atmósfera, de acuerdo con su concentración y propiedades físicas. La señal proviene de un láser de Nd-Yag a 532 y 1064 nm. Se emite hacia arriba (cenit) y hacia abajo (nadir) y es recogida con un telescopio, procesada y almacenada.
Campañas DAISEX DAISEX 1999 LIDAR
Campañas DAISEX DAISEX 1999 NEFELOMETRO Mide la radiación difundida por los aerosoles en tres longitudes de onda: 450 nm, 550 nm y 700 nm. Toma muestras de aire, en condiciones de presión, temperatura y humedad controladas, y las hace pasar la luz emitida por una lámpara. Un sistema de lentes y filtros recoge esa luz, que es amplificada por un fotomultiplicador.
Campañas DAISEX DAISEX 1999 NEFELOMETRO
Campañas DAISEX DAISEX 1999 PERFILES ATMOSFÉRICOS Vapor de agua Ozono
Campañas DAISEX DAISEX 1999 Resultados Espesor óptico para 550 nm Nefelómetro: ± 0.020 Lidar: ± 0.030 Irradiancia directa: ±
Descargar ppt "M. Pilar Utrillas Esteban Universitat de Valencia, 2002"

References: resolución 
 Resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución