R. Niclòs et al.inversión <strong>de</strong> la ecuación [1], requiere conocer laemisividad efectiva para cada superficie. Estaemisividad fue obtenida aplicando el Método <strong>de</strong>Cobertura Vegetal (MCV, Valor y Caselles, 1996y 2005). La proporción <strong>de</strong> vegetación se estimóconociendo las dimensiones <strong>de</strong> la estructura vegetal,asociada a un mo<strong>de</strong>lo <strong>de</strong> cubos <strong>de</strong> acuerdocon el MCV, y se comprobó a partir <strong>de</strong> ortofotos.En la zona A la proporción <strong>de</strong> vegetaciónpara la visión nadir era <strong>de</strong> 48,8%, aumentando a60% para un ángulo <strong>de</strong> 20° y 95% para 55°. Losvalores <strong>de</strong> emisividad <strong>de</strong> la vegetación, matorralcompuesto predominantemente por romero(Rosmarinus officinalis L.) y aliaga (Ulex parviflorusL.), y <strong>de</strong>l suelo (70% <strong>de</strong> Luvisol crómicoy 30% <strong>de</strong> Leptosol lítico) se obtuvieron conel Método <strong>de</strong> la Caja (Rubio et al., 1997; Rubioet al., 2003). La emisividad efectiva resultante<strong>de</strong> aplicar el MCV para esta zona es <strong>de</strong>0,986 ± 0,007 para un ángulo <strong>de</strong> observación <strong>de</strong>20° y <strong>de</strong> 0,989 ± 0,008 para 55°. Para la zona B,había una proporción <strong>de</strong> vegetación <strong>de</strong>l 100%durante los meses <strong>de</strong> verano (ver Fig. 2b), conuna emisividad <strong>de</strong> 0,985 ± 0,005.ANÁLISIS DE DATOSY RESULTADOST s -T a (K)151050–5T s -T a = 0,01308 I – 4,6r 2 = 0,73 error ajuste = 0,9 K0 200 400 600 800 1.000I (W m –2 )Figura 4. Depen<strong>de</strong>ncia <strong>de</strong> T s -T a con I en condiciones <strong>de</strong>cielos <strong>de</strong>spejados (zona A).Con el fin <strong>de</strong> establecer correlaciones entre laT a y la T s obtenida con la ecuación [1], en primerlugar nos centramos en los datos medidosdurante el día, cuando la irradiancia solar tieneun efecto importante. Para ello, consi<strong>de</strong>ramosun valor límite <strong>de</strong> 20° en el ángulo <strong>de</strong> elevaciónsolar, ya que por <strong>de</strong>bajo <strong>de</strong> este valor la imprecisiónen la medida <strong>de</strong> la irradiancia solar es elevada.Por otro lado, también distinguimos los datosmedidos con y sin nubosidad. Para ello,simulamos la irradiancia solar, I, que mediríanuestro piranómetro con el mo<strong>de</strong>lo SBDART 2.4(Ricchiazzi et al., 1998), el cual tiene en cuentala elevación <strong>de</strong>l terreno, y calculamos el cocienteentre el valor instantáneo medido en campo yel simulado. Cuando este cociente es próximo ala unidad (> 0,8) significa que el dato fue medidoen condiciones <strong>de</strong> nubosidad nula o <strong>de</strong>spreciable,y cuando difiere que la medida fue realizadacon cielos nubosos. A<strong>de</strong>más, para estadiscriminación también se consi<strong>de</strong>ró un valor límiteen la medida radiométrica <strong>de</strong> la temperaturaatmosférica, que fue <strong>de</strong> unos 253 K. En cualquiercaso este valor <strong>de</strong>pen<strong>de</strong> <strong>de</strong>l IFOV <strong>de</strong>l instrumentoutilizado y también ligeramente <strong>de</strong> lalatitud geográfica y la estación en que se realizanlas medidas.Los primeros resultados <strong>de</strong> este análisis muestranque la diferencia <strong>de</strong> temperaturas superficie-airetiene una <strong>de</strong>pen<strong>de</strong>ncia con la I y la velocidad<strong>de</strong>l viento en superficie, U. La Figura 4muestra dicha <strong>de</strong>pen<strong>de</strong>ncia con la I para los datosmedidos con cielos <strong>de</strong>spejados en la zona A,a<strong>de</strong>más <strong>de</strong> un ajuste lineal, con el que obtenemosun coeficiente <strong>de</strong> estimación razonable <strong>de</strong>0,7 y un error estándar <strong>de</strong> estimación para elajuste <strong>de</strong> ± 0,9 K.Así, la T a pue<strong>de</strong> obtenerse a partir <strong>de</strong> la T susando una función lineal <strong>de</strong> I (W m –2 ):T a = T s – ( a I + b ) [2]don<strong>de</strong> a = 0,01308 ± 0,00005 Km 2 W –1 y b =–4,6 ± 0,9 K. Si ajustamos una ecuación <strong>de</strong> estetipo separando los datos por rangos <strong>de</strong> velocida<strong>de</strong>s<strong>de</strong> viento en superficie, los coeficientes a yb muestran una ligera <strong>de</strong>pen<strong>de</strong>ncia con esta velocidad,pero los errores obtenidos indican quela precisión en la <strong>de</strong>terminación <strong>de</strong> T a no mejoraríaapreciablemente. El uso <strong>de</strong> la ecuación [2]para obtener mapas <strong>de</strong> T a a partir <strong>de</strong> imágenes<strong>de</strong> satélite <strong>de</strong> T s requiere a su vez disponer <strong>de</strong> informaciónespacial <strong>de</strong> irradiancia solar, que podríaobtenerse mediante un mo<strong>de</strong>lo <strong>de</strong> transferenciaradiativa teniendo en cuenta la topografía<strong>de</strong>l terreno. Dado que a escala regional el relieve<strong>de</strong>l terreno pue<strong>de</strong> ser variable, en futuros estudiospreten<strong>de</strong>mos comprobar si el ajuste pro-40 Revista <strong>de</strong> Tele<strong>de</strong>tección. ISSN: 1988-8740. 2010. 34: 36-43
Análisis <strong>de</strong> correlaciones entre la temperatura <strong>de</strong>l aire y la temperatura <strong>de</strong> las superficies vegetadas...151515101010T s -T a (K)5Ts-Ta (K)5T s -T a (K)50000 ms –1