Temperaturas mínimas estivales en el sureste de la península ibérica a partir de termografías satelitales

Revista de Estudios Andaluces (REA) | Núm. 42, julio (2021) | e-ISSN: 2340-2776. https://dx.doi.org/10.12795/rea.2021.i42 | © Editorial Universidad de Sevilla 2021 https://editorial.us.es/es/revistas/revista-de-estudios-andaluces https://revistascientificas.us.es/index.php/REA CC BY-NC-ND 4.0 8 8 https://dx.doi.org/10.12795/rea.2021.i42.01 Temperaturas mínimas estivales en el sureste de la península ibérica a partir de termografías satelitales Minimum summer temperatures in the southeast of the iberian peninsula through satellite thermographies David Espín-Sánchez [email protected] 0000-0003-4807-5450 Departamento de Geografía, Facultad de Letras, Universidad de Murcia. C/ Santo Cristo. 30001 Murcia, España. INFO ARTÍCULO RESUMEN Recibido: 11-12-2020 Revisado: 15-04-2021 Aceptado: 22-04-2021 PALABRAS CLAVE Gradiente altitudinal Inversión térmica Polo frío Termografía Canícula Temperaturas mínimas Sureste peninsular Teledetección La posición latitudinal del sureste de la península ibérica genera temperaturas medias elevadas (entre 25 y 30oC) durante la época estival, si bien, su compleja orografía y el importante gradiente altitudinal, genera importantes diferencias en las temperaturas mínimas, influenciadas especialmente por intensos procesos de inversiones térmicas nocturnas (ITN) en valles de interior. El incremento de nuevas estaciones meteorológicas y el aumento de la densidad de las redes de observación, unida a la implantación de una red propia de termómetros registradores en polos fríos (CAP) de las provincias de Murcia, Albacete o Jaén, lleva a un mayor conocimiento de la distribución de temperaturas nocturnas superficiales. El principal objetivo de la investigación es la generación de cartografía (termografías) que refleje la distribución de las temperaturas nocturnas en el sureste de la Península Ibérica, realizada a través del tratamiento de productos satelitales (VIIRS LST Band I5 del satélite SUOMI-NPP. Unido al uso de las imágenes satelitales (termografías) es posible analizar detalladamente los contrastes térmicos existentes durante el periodo canicular a través de 147 estaciones meteorológicas y termómetros registradores (15 de julio al 15 de agosto) de los años 2017 a 2020. La verificación de las termografías utilizadas (producto Band I5) muestra un notable ajuste estadístico (R2 de 0,87), con un Bias y RMSE de -0,2 y 1,83oC respectivamente. La cartografía es realizada para las diferencias provincias que integran el sureste de la península ibérica (Alicante, Albacete, Murcia, Jaén, Granada y Almería). Las diferencias térmicas diarias aportadas por las estaciones meteorológicas superficiales arrojan contrastes de temperaturas mínimas de hasta 23,4oC entre valles elevados de interior y ámbitos litorales. KEYWORDS ABSTRACT Lapse rate Temperature inversion Cold air pools Thermography Hottest days of summer Sureste peninsular Teledetection The latitudinal position of the southeast of the Iberian Peninsula generates high average temperatures (between 25 and 30oC) during the summer season, although its complex orography and the important lapse rate (LR) give rise to important differences in minimum temperatures, influenced by intense processes of nocturnal temperature inversion (ITN) in inland valleys. The increase of new meteorological stations and the rise in the density of observation networks, together with the implementation of its own network in cold air pools (CAP) in the provinces of Murcia, Albacete or Jaén, leads to a greater understanding of the distribution of surface temperatures. The main objective of the research is the generation of cartography (thermographs) that reflects the distribution of night temperatures in the southeast of the Iberian Peninsula, carried out through the treatment of satellite products (VIIRS LST Band I5 of the SUOMI-NPP satellite). Together with the use of satellite images (thermography), it is possible to analyze in detail the thermal contrasts that exist during the hottest days of summer through 147 weather stations and recording thermometers DOI: https://dx.doi.org/10.12795/rea.2021.i42.01 Formato de cita / Citation: Espín-Sánchez, D. (2021). Minimum summer temperatures in the southeast of the iberian peninsula through satellite thermographies. Revista de Estudios Andaluces, 42, 8-33. https://dx.doi.org/10.12795/rea.2021.i42.01 Correspondencia autores: [email protected] (David Espín-Sánchez). David Espín-Sánchez / REA N. 42 (2021) 8 - 33 9 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 (July 15 to August 15) of the year 2017 to 2020. Verification of the thermographs used (Band I5 product) shows a notable statistical adjustment (R2 of 0.87), with a Bias and RMSE of -0.2 and 1.83oC respectively. The cartography is carried out for the differences in the provinces that make up the SE of the Iberian Peninsula (Alicante, Albacete, Murcia, Jaén, Granada and Almería), with daily thermal differences in minimum temperatures of up to 23.4oC between inland valleys and coastal areas. 1. INTRODUCCIÓN Una inversión térmica nocturna (ITN) supone un cambio térmico en el gradiente altitudinal, con temperaturas más bajas por debajo del estrato de inversión, y valores más elevados en un nivel superior. La deposición de aire más frío y denso en los fondos de valle y depresiones se lleva a cabo durante horas nocturnas de estabilidad atmosférica. El viento sinóptico débil y ausencia de nubosidad durante las madrugadas estables en zonas orográficamente deprimidas o de valle, provoca una pérdida de radiación y calor en la superficie del suelo (Whiteman et al., 2004), destruyéndose por convección después del amanecer (Arduini et al., 2020). Se generan las denominadas piscinas de aire frío (cold air pools) o CAP, una capa de aire estable y estratificada que se encuentra confinada en el fondo de un valle o cuenca. El proceso de enfriamiento local genera un estrato superficial de aire frío confinado y desacoplado del flujo sinóptico más cálido de niveles altitudinalmente superior (Sheridan, 2019). Un gran número de investigaciones (Sauberer & Dirmhirn, 1954, 1956; Litschauer, 1962, Whiteman et al., 2004, 2008; Eisenbach et al., 2003; Pospichal et al., 2003) han puesto de manifiesto la importante variabilidad espacial en las temperaturas mínimas sobre orografía irregular y compleja, y el acopio máximo de frío sobre plataformas kársticas y campos de dolinas, determinando, que es ahí, donde se registran las temperaturas más bajas del planeta, como los -56oC en Peter Sink (Utah, EE.UU) o -50oC en Gstettneralm (Austria). La acumulación de aire frío en los fondos de la depresión se ve reforzada cuando actúan flujos catabáticos, potenciándose el proceso a razón del incremento del confinamiento de las tierras bajas (Gustavsson et al., 1998; Lindkvist et al., 2000; Lundquist et al., 2008). Así, como resultado de la influencia de estas variables climáticas y topográficas, la intensidad y el espesor de los fenómenos de inversión térmica se alteran de acuerdo con la situación sinóptica, la persistencia del balance energético negativo y la configuración topográfica del valle (McKee & O’Neal, 1989; Whiteman, 2000; Clements et al., 2003; Daly et al., 2009). Los extraordinarios contrastes térmicos entre los fondos de valle y ámbitos geográficamente accidentados lleva a la necesidad de poseer información espacial con un creciente interés científico en busca de nuevas metodologías satelitales para rellenar los vacíos climáticos existentes. Los recientes avances en los instrumentos de teledetección y potencia de cálculo, junto con mejoras en parametrizaciones y algoritmos numéricos, dan lugar a análisis detallados de posibles acumulaciones de aire frío en fondos de valle. Se demuestra, por tanto, que las ITN tienen un efecto muy notable sobre el clima, con formación persistente de estratos fríos superficiales (Zhang et al., 2011). La previsión en invierno de la acumulación y de aire frío persistentes se ve obstaculizada por la falta actual de comprensión en los mecanismos involucrados, y es considerado uno de los problemas más difíciles de pronosticar en cuencas y valles En años recientes, los importantes avances en teledetección, han conseguido el desarrollo de nuevos algoritmos y ecuaciones para la obtención de temperatura del aire a 1,5 m de altura (Ta) a través de la temperatura de la superficie a 2 m (Ts), obtenida a partir de imágenes de sensores remotos satelitales (Vancutsem et al., 2010; Zhu et al., 2013). La Ta es un parámetro de gran importancia para una amplia gama de campos, como la agricultura (evapotranspiración de cultivos, predicción del rendimiento de cultivos, hidrología, ecología, medio ambiente y cambio climático). Los productos LST de diferentes sensores (MODIS, SUOMI NPP, NOAA20, AVHRR…) se han utilizado con éxito en diversos ámbitos científicos, por ejemplo, para la estimación de la evapotranspiración (Anderson et al., 2011), estimación de la temperatura del aire (Vancutsem et al., 2010) y el monitoreo de islas de calor urbano (Rajasekar & Weng, 2009). El uso de productos LST puede ayudar al desarrollo de registros de datos climáticos LST de resolución moderada a largo plazo (Yu et al., 2008). David Espín-Sánchez / REA N. 42 (2021) 8 - 33 10 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 La gran resolución espacial y temporal (diaria), datos en abierto, y el fácil acceso a las imágenes satelitales, ha generado un importante incremento de los estudios que emplean la Ta (Noi et al., 2016; Sismanidis et al., 2016; Zhou et al., 2017), especialmente utilizando información del instrumento Visible Infrared Imaging Radiometer Suite (VIIRS) a bordo del satélite Suomi National Polar-Orbiting Partnership (Suomi NPP) (Hillger et al., 2013). También son muy utilizados los productos generados por el sensor Moderate Resolutions Imaging Spectroradiometer (MODIS6), que viaja a bordo de las plataformas satelitales Terra (1999) y Aqua (2002) (Wan et al., 2015). Se han determinado correlaciones estadísticas muy notables (R2 > 0,90) y un sesgo mínimo (RMSE y MD < 7oC) entre los productos que estiman Ta y la temperatura del aire a 2 m registradas por estaciones meteorológicas superficiales (T), demostrando ser lo suficientemente precisos para la estimación de temperaturas mínimas superficiales multitud de ámbitos planetarios (Liu et al., 2015; Hulley et al., 2017; Niclòs et al., 2018). La interpolación espacial ha sido ampliamente aplicada para la cartografía de temperaturas mínimas (Thornton et al., 1997; Jarvis & Stuart, 2001), donde la mayoría de los modelos de interpolación espacial incorpora una corrección de gradiente vertical con el fin de adaptar la diferencia de elevación y evitar los importantes sesgos generados (Dodson & Marks, 1997; Nalder & Wein, 1998). El gradiente varía temporalmente a escalas estacionales, diarias e incluso horarias. Además, la ubicación de la estación de referencia puede afectar al gradiente vertical (Bolstad et al., 1998) causando un sesgo positivo para elevaciones más bajas. Por tanto, aunque el método de interpolación ha sido ampliamente aplicado para la estimación de temperaturas mínimas durante fenómenos de ITN, su limitación espacial y vertical limita su uso. Para paliar las limitaciones de las técnicas de interpolación espacial de las últimas décadas (Dodson & Marks, 1997), debido a la distribución desigual y escasa de los puntos de observación (Bolstad et al., 1998), se utilizan productos satelitales donde se recupera la temperatura superficial (LST). Para reflejar la distribución de temperaturas nocturnas en el sureste de la península ibérica, se genera cartografía (termografías) a través de productos satelitales, siendo uno de los principales objetivos de la investigación. Se analizan, también los contrastes térmicos en la época estadísticamente más cálida del año (periodo canicular de 2017 a 2020), con el apoyo de datos observacionales en estaciones meteorológicas y termómetros registradores, que son utilizados también para el segundo gran objetivo de la investigación, el de la validación del produc- to satelital en el área de estudio 2. DATOS Y MÉTODOS 2.1. Datos El análisis de temperaturas mínimas superficiales se sustenta a través de la utilización de temperaturas mínimas diarias y horarias pertenecientes a 147 estaciones meteorológicas automáticas (EMA) (mapa 1) distribuidas espacialmente en el sureste de la península ibérica (provincias de Alicante, Albacete, Murcia, Almería, Granada y Jaén). Un total de 92 estaciones pertenecen al Sistema de Información Agroclimática para el Regadío (SIAR) del Ministerio de Agricultura, Pesca y Alimentación; 32 a la Red de estaciones agroclimáticas de la Consejería de Agricultura, Ganadería, Pesca y Desarrollo Sostenible de la Junta de Andalucía; 8 a la Agencia Estatal de Meteorología (AEMET), y finalmente 15 termómetros registradores pertenecientes a la red propia del autor (Espín-Sánchez et al., 2018b). La utilización de la red propia de datalogger (cuadro 1) se justifica por la utilización de enclaves con temperaturas medias de las mínimas más bajas que las aportadas por las redes oficiales, en un esfuerzo de llegar a áreas pobladas de alta montaña y zonas remotas. Los termómetros registradores (dataloggers) utilizados son en su mayoría de marca Elitech RC-5 de 25.000 registros (con certificación de calibrado del CNAS, China National Accreditation Service), dos de ellos son Microlite 32K (32.000 registros), y finalmente, ubicados en los dos emplazamientos con las temperaturas medias más frías, y bajo condiciones invernales más severas, se utilizan dos registradores HOBO MX25 de 50.000 registros. En la mayoría de observatorios se utiliza una garita meteorológica artesanal compuesta por 11 platos blancos, en 3 de ellas la garita meteorológica homologada es tipo Davis (7714), y finalmente, para los termómetros ubicados en los parajes de Navalasno y Nava del Pol- David Espín-Sánchez / REA N. 42 (2021) 8 - 33 11 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 vo se ubican en el interior de dos garitas meteorológicas tipo pagoda, cedidas por AEMET Murcia (imagen 1). Los registradores utilizados poseen un rango temporal diezminutal, con una precisión de ±0,3˚C y una resolución térmica de 0,1°C. Los ubicados en Monterilla y Navalasno son calibrados por la Delegación Territorial de AEMET de Murcia, en el Centro Meteorológico Territorial de Guadalupe. Los responsables del mantenimiento del instrumental determinan una desviación de +0,2o en las temperaturas mínimas. Las temperaturas mínimas diarias y de rango temporal diezminutal han sido analizadas exclusivamente para cuatro periodos caniculares (del 15 de julio al 15 de agosto) para los años 2017, 2018, 2019 y 2020. Desde el punto de vista estadístico, el periodo canicular posee la temperatura media más elevada del año (Mosiño & García, 1966). La realización de la cartografía de temperaturas nocturnas superficiales (termografías) ha sido elaborada a través del producto satelital VIIRS LST, obtenido del instrumento VIIRS, perteneciente al satélite SUOMI-NPP (https://worldview.earthdata.nasa.gov/). Los productos satelitales del instrumento VIIRS se dividen en 22 bandas con una cobertura espectral de 0,4 a 12 μm. Entre las 22 bandas hay 5 canales de imágenes de alta resolución (bandas I), 16 canales de resolución moderada (bandas M) y banda de 1 día / noche (DNB). (Cao et al., 2013). Se utiliza específicamente el subproducto Band I5 (una capa de temperatura de brillo calcula a partir de radiancias calibradas de VIIRS VNP02), con una resolución espacial de 375 m, y una resolución temporal diaria (Cao et al., 2013). El sesgo de medición y la precisión especificada para el producto es de 1,4 y 0,5 K respectivamente, que deben cumplirse cuando la máscara de nube indica una alta confianza en condiciones claras (Minnett et al., 2014). Posee una órbita sincrónica solar a una altitud de 829 km y con unos tiempos de paso sobre la vertical del área de estudio entre las 1.40 y 2.40 h de la madrugada (Niclòs et al., 2018). Para la validación estadística del producto VIIRS LST (Band I5) se utilizan las temperaturas superficiales de los observatorios meteorológicos anteriormente descritos, tomadas a las 2.00 h AM debido a su mayor proximidad con el paso del satélite SUOMI NPP sobre la vertical del área de estudio. Mapa 1. Estaciones meteorológicas utilizadas en el área de estudio. Fuente: Elaboración propia. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 12 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Cuadro 1. Características de los registradores de temperatura (dataloggers) de la red propia. Paraje Provincia Altitud Protección Sensor Nava del Polvo Jaén 1836 Pagoda AEMET HOBO MX25 Navalasno Jaén 1770 Pagoda AEMET Microlite 32K Monterilla Jaén 1660 David 7714 HOBO MX25 Viboreros Albacete 1555 Garita casera Elitech RC-5 Pico de la Selva Murcia 1521 Garita Davis Elitech RC-5 Vega de Pontones Jaén 1338 Garita casera Elitech RC-5 Nevazo de Arriba Murcia 1181 Garita casera Microlite 32K La Rogativa Murcia 1152 Garita casera Elitech RC-5 Cañada Grande Murcia 1087 Garita casera Elitech RC-5 La Risca Murcia 1047 David 7714 Elitech RC-5 Gollaín Murcia 1021 Garita casera Elitech RC-5 Bco. del Agua Murcia 906 Garita casera Elitech RC-5 Bco. del Cenajo Murcia 865 Garita casera Elitech RC-5 Benablón Murcia 777 Garita casera Elitech RC-5 Avilés Murcia 762 Garita casera Elitech RC-5 Fuente: Elaboración propia. Imagen 1. Sensor térmico (datalogger) localizado en el poljé colgado de Navalasno (Jaén), a una altitud de 1770 m. Garita meteorológica de tipo Pagoda, cedida por AEMET Murcia, con sensor microlite 32K en su interior. Fuente: Autor. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 13 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 2.2. Métodos El producto VIIRS LST (Band I5) utiliza el algoritmo de ventana dividida (Yu et al., 2005) recurriendo a la temperatura de brillo (Ti) medidas en dos bandas espectrales VIIRS: M15 colocada a 10,76 mm (i = 15) y M16 colocada a 12,01 mm (i = 16). Se utilizan diferentes conjuntos de coeficientes (aj, con j = 0–4) para el día y la noche y para 17 tipos de superficie de los mapas de clasificación del Programa Internacional de Geosfera-Biosfera (IGBP) (Guillevic et al., 2014; Liu et al., 2015). La línea de base es un algoritmo de ventana dividida (SW) que aplica temperaturas de brillo medidas en las bandas VIIRS M15 y M16 centradas en longitudes de onda de 10,8 μm y 12,0 μm, respectivamente (Sun & Pinker, 2003; Schroeder & Giglio, 2017). El algoritmo opcional de ventana de doble división (DSW) aplica dos bandas infrarrojas de onda corta adicionales M12 y M13 centradas en longitudes de onda de 3,75 μm y 4,0 μm, respectivamente. Desde el 11 de agosto el producto opera utilizando un algoritmo de ventana dividida única, que es insensible a la radiación solar (Yu et al., 2005. Se utiliza el algoritmo nocturno Nighttime DSW (Ecuación 1): LSTi = b0 (i) + b1 (i)T11 + + b2 (T11 – T12) + b1. (i) (secθ-1) + b4 (i)T3.75 + b5 (i)T4.0 + b6 (i)T3.752 + b7 (i)T4.02 + b8 (i) (T11 – T12)2 (1) donde i es el índice de 17 tipos de superficies IGBP; aj (i), bj (i) y cj (i) son los coeficientes de regresión del algoritmo en los que j representa la posición secuencial del término en la ecuación; T3.75, T4.0, T11 y T12 son las temperaturas de brillo de las bandas VIIRS M12, M13, M15 y M16, respectivamente; θ es el ángulo cenital del satélite y φ es el ángulo cenital solar. El algoritmo DSW consta de versiones diurnas y nocturnas, y se aplicó como algoritmo de referencia antes del 10 de agosto de 2012. La elaboración de la cartografía final (termografías) para las diferentes provincias analizadas es realizada con el promedio de 50 mapas diarios a través del software de licencia gratuita QGIS 3.10.7. Se trata de 50 madrugadas estables (con procesos de inversión térmica nocturna), sin cobertura nubosa (0%), y complementadas con el apoyo del análisis de datos de temperaturas mínimas superficiales durante los periodos caniculares de observación (2017, 2018, 2019 y 2020). La validación estadística de Band I5 es efectuada mediante un modelo lineal de regresión con las temperaturas registradas en las estaciones meteorológicas. Es llevada a cabo para diferentes umbrales altitudinales y térmicos mediante la función “Extract Multi Values to point” del software ArcGis 10.3. A través de la exportación de ficheros de texto plano (txt) se obtiene información de la altitud y la temperatura de la termografía promedio generada para el pixel de cada una de las 147 estaciones meteorológicas utilizadas Se utiliza el test estadístico de Shapiro Wilk para contrastar la normalidad (distribución normal o no) de un conjunto de datos (Shapiro, 1998). Si la muestra es normal se aplica el test de correlación lineal de Pearson (Eq. 2): W = 1/n2 (Σh j=1 ain(xn-i+1-xi)2 (2) donde: s2 es la varianza muestral, y los coeficientes ain suelen aparecer tabulados en los manuales. Por el contrario, si la muestra analizada presenta asimetría y no presenta una distribución normal, se aplica el test no paramétrico del coeficiente de correlación de Spearman. Dicho test examina la fuerza y la dirección de la relación monótona entre dos variables continuas u ordinales. El valor del coeficiente de correlación puede variar de –1 a +1. Cuanto mayor sea el valor absoluto del coeficiente, más fuerte será la relación entre las variables. Para calcular ρ, los datos son ordenados y reemplazados por su respectivo orden (Eq. 3): rR = 1 - (6Σi di2) n(n2–1) donde n= número de puntos de datos de las dos variables y di= diferencia de rango del elemento “n” El procedimiento aplicado minimiza simultáneamente dos funciones de costo definidas como el sesgo (BIAS) y el error cuadrático medio (RMSE) entre la LST simulada de Band I5 y la temperatura observada. El error cuadrático medio o root mean square deviation (RMSE), que mide la cantidad de error existente entre dos conjuntos de datos, es decir, un valor predicho y un valor observado o conocido (Eq. 4) (Chai & David Espín-Sánchez / REA N. 42 (2021) 8 - 33 14 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Draxler, 2014). El RMSE, por tanto, es ideal para la validación de las estimaciones de temperatura del produc- to Band I5 (predicho), y las temperaturas registradas en estaciones meteorológicas (observado o conocido). Este estadístico ha sido utilizado en numerosas validaciones de productos como MODIS, Band I5… (Guillevic et al., 2012; Li et al., 2014; Duan et al., 2019). RMSE = n donde: Σ es el sumatorio, Pi es el valor predicho para la ith observación en el conjunto de datos, Oi es el valor observado para la ith observación en el conjunto de datos, y n es el tamaño de la muestra. El Bias es una medición de la “exactitud” del sistema de medición y representa al error sistemático del sistema. Es la contribución al error total debido a los efectos combinados de todas las fuentes de variación, conocidas o desconocidas (Voinov et al., 1993). Este estudio estadístico se resuelve a través de una prueba de hipótesis t Student de una muestra y una constante, con las siguientes hipótesis. Cuando el p value es mayor a 0.05 se acepta la hipótesis nula y podemos concluir que el sistema de medición no presenta ningún sesgo. Bias (^θS,θ) = ES Dm [^θs] – θ donde Es denota el valor esperado sobre la distribución estadística, ^θ sirve como estimador de θ basado en cualquier dato observado. 3. RESULTADOS 3.1. Validación estadística del producto VIIRS LST (BandI5) El producto VIIRS LST (Band I5) es sometido a una verificación estadística con las temperaturas superficiales registradas en las estaciones meteorológicas de superficie. Las correlaciones estadísticas para el conjunto de la muestra son relevantes (R2 de Spearman de 0,87), con un nivel de confianza de 0,99 (gráfico 1). El Bias, y el RMSE de las diferencias entre Band I5 y las estaciones meteorológicas es de -0,2 y 1,83 ºC respectivamente para el conjunto del área de estudio. RMSE es una forma útil para analizar qué tan bien un modelo de regresión (o cualquier modelo que produzca valores predichos) es capaz de “ajustarse” a un conjunto de datos. En esta ocasión los resultados arrojan una gran similitud entre ambos conjuntos (cuadro 2). Las correlaciones estadísticas más importantes se dan en zonas de elevada altitud, aquellos ámbitos con temperaturas medias más bajas (R2 de 0,76). Las diferencias térmicas son las más importantes de los rangos altitudinales analizados, motivados principalmente por los procesos de inversión térmica nocturna que se localizan en valles de poca extensión territorial. El análisis nocturno efectuado arroja excelentes resultados, pues existe una mejor predicción para temperaturas mínimas (Tmin) que para las temperaturas máximas (Tmax) en el análisis de productos LST (Yang et al., 2017). En general, y a excepción de la franja altitudinal superior a 1000 m, la gran mayoría del área de estudio presenta un Bias y RMSE no superior a -0,51 y 1,83oC. La franja altitudinal más elevada (>1000 m) posee una menor muestra estadística (102) ante el menor número de estaciones meteorológicos en ámbitos montañosos (gráfico 2d.). Aun así, obtiene una correlación estadística importante (R2 0,77) y un Bias no demasiado importante (-1,28oC). La diferencia térmica entre el producto Band I5 y las temperaturas de las estaciones meteorológicas utilizadas es de -0,2oC. Es decir, las imágenes de satélite infraestiman muy ligeramente las temperaturas mínimas registradas en superficie, especialmente allí donde los procesos de inversión térmica nocturna son más intensos, generalmente en zonas altitudinalmente más elevadas. Se trata, por tanto, de una aproximación realmente notable, constituyendo una herramienta esencial para el estudio de la distribución de las temperaturas nocturnas. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 15 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Las correlaciones estadísticas en el resto de intervalos altitudinales son similares, entre 0,80 y 0,75, con diferencias térmicas poco importantes (-0,5 a -0,3oC). Entre 250 y 1000 m. de altitud, los resultados estadísticos muestran resultados muy importantes (Cuadro 2). Sin embargo, los observatorios con altitudes inferiores a 250 m (con los valores más templados) registran una correlación algo inferior (0,54), pero corroborada estadísticamente con un Bias y RMSE de 0,41 y 1,63oC (gráfico 2d). En este rango altitudinal, la diferencia de temperatura es de 0,4oC, donde los valores estimados por el producto Band I5 infraestiman la temperatura superficial registradas por los observatorios, seguramente motivado por el efecto de las islas de calor urbana (ICU), al coincidir con el área más densamente pobladas. Cuadro 2. BIAS, RMSE, correlación estadística y diferencia media calculadas entre VIIRS LST (Band I5) y las temperaturas registradas en estaciones meteorológicas de superficie. Área de estudio Altitud >1000 m Altitud 500 – 1000 m Altitud 250 – 500 m Altitud < 250 m N (nº muestra) 799 102 265 150 282 T. Shapiro-Wilk (p value < 0.01) 2,2e-16 0.03063 0.0006263 0.02751 0.3362 Tipo de test Spearman Pearson Spearman Pearson Pearson Rho / R20.872 0.766 0.749 0.840 0.563 Bias -0.22 -1.28 -0.31 -0.51 0.41 RMSE 1.83 2.84 1.73 1.48 1.63 Diferencia media -0.22 -1.27 -0.31 -0.51 0.40 Fuente: Elaboración propia. Gráfico 1. BIAS, RMSE y correlación estadística calculadas entre VIIRS LST (Band I5) y las temperaturas registradas en estaciones meteorológicas de superficie (TS). Fuente: Elaboración propia. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 16 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Gráfico 2. BIAS, RMSE y correlación estadística calculadas entre VIIRS LST (Band I5) y las temperaturas registradas en estaciones meteorológicas de superficie (TS) en diferentes franjas altitudinales: a) >1000 m y b) 500-1000 m, c) 250-500 m y d) <250 m. Fuente: Elaboración propia. 3.2. Distribución espacial de las temperaturas nocturnas estivales 3.2.1. Provincia de Jaén La distribución espacial de los valores térmicos nocturnos en la provincia de Jaén viene motivada por la complejidad de la orografía. Los píxeles más elevados según el producto Band I5 aparecen en el casco urbano de Andújar (21,0oC), Sierra de Andújar (20,7oC), Jaén capital (20,8oC) y el núcleo urbano de Linares (20,8oC), donde los valores de las jornadas analizadas superan el umbral de los 20oC (mapa 2). La diferencia de los cascos urbanos frente a las zonas rurales limítrofes es importante, con valores que descienden a 17,0oC en la periferia de los cascos urbanos. Por su parte, atendiendo a los registros de las EMAs utilizadas, los valores más elevados se registran en Jaén capital (21,9oC), con temperaturas de hasta 29,8oC (realmente notables). Otras cabeceras municipales como Bedmar o Jódar, situadas en resaltes topográficas y fueras de zonas de ITN, registran valores mínimos superiores a 28,0oC (cuadro 3). Las temperaturas más bajas de la provincia (Band I5) aparecen en el Calar de Hernán Pelea (Sierra de Segura), con un valor mínimo de 8,9oC en el paraje del Campo del Espino, considerado en recientes publicaciones como una de las zonas más frías de la Demarcación Hidrográfica del Segura (DHS) y del sur peninsular (Espín Sánchez et al., 2018b). Las temperaturas inferiores a 10,0oC aparecen en las inmediaciones de Santiago de la Espada (Vega río Zumeta), en el río Madera o en Pontones (zona occidental y de mayor altitud media de la provincia). Las temperaturas en el Calar de Magina se sitúan en 10,4oC, o valores de 11,9oC al sur de la provincia, en las inmediaciones de los municipios de Alcalá la Real y Guadahortuna – Huelma. Se trata principalmente de calares, depresiones o valles situados a altitudes superiores a 1300 m de altitud. Zonas de intensos procesos de ITN incluso en época estival. Las temperaturas registradas en estos sectores, a través de termómetros registradores, se sitúan en temperaturas medias de las mínimas de 8,4oC en la Nava del Polvo, 8,6oC en Monterilla, o 12,4oC en el núcleo poblacional de Pontones, donde los registros mínimos durante los últimos periodos caniculares son excepcionales en la mitad sur de la península ibérica, llegando a descender a 5,3oC en Pontones (28/07/2018). David Espín-Sánchez / REA N. 42 (2021) 8 - 33 23 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Mapa 5. Termografía elaborada a través del producto VIIRS LST (Band I5) para la provincia de Albacete. Estaciones meteorológicas utilizadas y promedio de las temperaturas mínimas durante el periodo de canícula (15 julio – 15 agos- to) en el periodo temporal (2017-2020). Recuadro de color azul representa los ámbitos con los píxeles más fríos, mientras, que, los recuadros negros hacen referencia a las temperaturas más elevadas. Fuente: Elaboración propia. Cuadro 6. Cuadro estadístico de las temperaturas mínimas diarias registradas en las estaciones meteorológicas de la provincia de Albacete durante las canículas (15 julio – 15 agosto) de 2017 a 2020. Estación meteorológica Temperatura media de las mínimas (oC) Temperatura mínima absoluta (oC) Temperatura máxima de las mínimas (oC) Pozo Cañada 17,8 13,5 21,7 Caudete 17,7 13,5 22,2 Albacete 17,6 13,2 22,5 Juanaco 17,5 11,8 22,4 Ontur 17,2 12,4 21,5 Motilleja 16,6 12,4 20,3 Almansa 16,3 12,0 21,1 La Gineta 16,2 12,0 20,0 Cañadas de Nerpio 13,5 8,7 21,9 Viboreros 12,7 4,9 19,9 Fuente: Elaboración propia. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 24 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 3.2.5. Provincia de Murcia La provincia de Murcia se caracteriza por presentar dos ámbitos geográficos con los promedios de temperatura mínimas más bajas, al oeste de la Región (comarca del Noroeste), y la zona más septentrional (Altiplano de Murcia). El resto de la provincia registra temperaturas, por lo general, superiores a 18/19oC, destacando especialmente algunos ámbitos donde las temperaturas de los píxeles (Band I5) superan los 20oC. Destaca el área metropolitana de Murcia, donde se registran temperaturas de 21,4oC en el centro o 21,0oC en el Centro Meteorológico Territorial (CMT) de AEMET. Otros municipios próximos a la capital registran valores nocturnos elevados, como Alcantarilla (20,7oC), Molina de Segura (20,6oC) o Archena y Ricote (20,3oC). Los píxeles más elevados aparecen en el sector litoral más meridional, con valores más elevados en La Manga (22,9oC), Cabo Cope (22,3oC) o Cabo Tiñoso y Cartagena (22,0oC) (mapa 6). Las temperaturas medias de las mínimas de las EMAs son especialmente elevadas en la ciudad de Murcia y zonas litorales. El CMT de AEMET (Murcia) obtiene una temperatura de 23,0oC, una de las más elevadas de las estaciones meteorológicas analizadas. En los cuatro periodos caniculares observados, la temperatura mínima absoluta es de 20,0oC, registrando 127 madrugadas tropicales sobre 128 jornadas posibles (cuadro 7). La Pilica (22,6oC) o La Alberca (22,4oC), ubicados en zonas próximas a la capital de Murcia, en terrenos orográficamente accidentados, registran temperaturas medias de las mínimas superiores a 22,0oC, al igual que ocurre en Los Belones (22,6oC) o en Tébar (22,1oC). Las temperaturas mínimas son realmente contrastadas con el sector más occidental de la provincia, donde destacan valles relativamente elevados (>1000 m) muy propenso a los procesos de ITN (Mapa 6). El píxel estimado (Band I5) más bajo durante el periodo canicular aparece al sur de El Entredicho (11,5oC), en el macizo de Revolcadores (11,7oC), o en el altiplano de Caravaca de la Cruz, en las pedanías de El Moralejo y La Junquera (11,8 a 12,0oC). Se trata en general, del área con las temperaturas más bajas de la provincia, donde destacan también algunos ámbitos habitados como el Campo de San Juan (12,0oC), Valle de la Rogativa (12,1oC) o Calar de la Santa (12,5oC). Las temperaturas arrojadas por los termómetros registradores en el sector más occidental de la provincia son similares, destacando los valles de Gollaín, La Rogativa, Nevazo o La Risca, con promedios entre 13,5 y 13,7oC. En estos ámbitos se registran temperaturas mínimas absolutas ligeramente inferiores a 7 ºC en pleno periodo canicular (6,9oC en La Risca el 17/07/2018), y ausencia de noches tropicales en Benablón o Gollaín, donde las temperaturas mínimas más elevadas no superan 19,3 y 18,2oC respectivamente. Gráfico 6. Diagrama de cajas o boxplot de las temperaturas mínimas registradas en cuatro estaciones meteorológicas de la provincia de Albacete durante cuatro periodos caniculares (2017-2020). Fuente: Elaboración propia. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 25 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Las diferencias entre los extremos térmicos son notables, con temperaturas mínimas absolutas promedios que rozan 10oC entre el CMT de Murcia y La Risca. El gráfico 7 muestra la importante diferencia entre el P25 de las temperaturas mínimas del CMT de Murcia (21.7oC) y el P95 de La Risca (14.8oC). Mapa 6. Termografía elaborada a través del producto VIIRS LST (Band I5) para la provincia de Murcia. Estaciones meteorológicas utilizadas y promedio de las temperaturas mínimas durante el periodo de canícula (15 julio – 15 agosto) en el periodo temporal (2017-2020). Recuadro de color azul representa los ámbitos con los píxeles más fríos, mientras, que, los recuadros negros hacen referencia a las temperaturas más elevadas. Fuente: Elaboración propia. Cuadro 7. Cuadro estadístico de las temperaturas mínimas diarias registradas en las estaciones meteorológicas de la provincia de Murcia durante las canículas (15 julio – 15 agosto) de 2017 a 2020. Estación meteorológica Temperatura media de las mínimas (oC) Temperatura mínima absoluta (oC) Temperatura máxima de las mínimas (oC) CMT Murcia 23,0 20,0 26,5 Los Belones 22,6 18,1 26,2 La Pilica 22,6 18,6 27,1 La Alberca 22,4 19,8 26,2 Tébar 22,1 19,0 25,5 Benablón 14,5 10,0 19,3 La Risca 13,7 6,9 20,5 El Nevazo 13,7 7,2 22,1 La Rogativa 13,6 9,5 19,1 Gollaín 13,5 8,3 18,2 Fuente: Elaboración propia. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 26 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 3.2.6. Provincia de Alicante La provincia de Alicante posee una orografía similar a Murcia, lo que determina un litoral cálido durante las noches estivales, y un interior de mayor altitud media y multitud de valles, que le otorgan un carácter menos tórrido. Los píxeles más elevados según el producto Band I5 se sitúan en el litoral más meridional, como Cabo Roig (Orihuela) con 22,1oC, Torrevieja (21,9oC), puerto de Alicante (21,8oC) o Santa Pola (21,8oC). También el litoral más septentrional como Denia registra valores medios de 21,3oC. Finalmente, algunos ámbitos prelitorales como Crevillente y Elche registran valores medios de 21,1oC (mapa 7). Las temperaturas registradas por las EMAs son similares a las estimadas por las termografías, destacando los valores medios en Crevillente (22,7oC), Alicante (22,6oC) o Elche (22,4oC), donde las temperaturas mínimas absolutas difícilmente descienden de 20oC (Cuadro 8). Las temperaturas mínimas absolutas más elevadas llegan a registrar más de 27oC en Crevillente, Alicante o Villajoyosa. El interior orográficamente más accidentado es el que registra las temperaturas nocturnas estivales más bajas de la provincia. Destaca, como pixel más frío (Band I5), Banyeres de Mariola (13,7oC), mientras que, otros emplazamientos como el sur de Alcoy (13,8 ºC), sur de Banyeres (13,8oC) o el sur de Villena (14,1oC) registran valores frescos en las noches estivales. Las estaciones meteorológicas automáticas utilizadas no están ubicadas en los sectores geográficos anteriormente descritas, por lo que las temperaturas medias de las mínimas son algo más elevadas (>17oC). Destacan los municipios de Planes (18,3oC), Pinoso (18,2oC) o la más baja en Villena (17,2oC). Las temperaturas mínimas absolutas, lejos de ser tan bajas como en valles de alta y media montaña, registran temperaturas mínimas absolutas de 11,7oC en Planes o 12,6oC en Villena (cuadro 8). Las temperaturas nocturnas no son tan contrastadas como en las provincias de Murcia o Jaén, donde se han utilizado datalogger ubicados en polos fríos relevantes. Aun así, las diferencias medias alcanzan los 5,5 ºC de media entre Villena y Crevillente, con P25 de 16oC en Villena y 21.5oC en Alicante (gráfico 8). Gráfico 7. Diagrama de cajas o boxplot de las temperaturas mínimas registradas en cuatro estaciones meteorológicas de la provincia de Murcia durante cuatro periodos caniculares (2017-2020). Fuente: Elaboración propia. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 27 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Mapa 7. Termografía elaborada a través del producto VIIRS LST (Band I5) para la provincia de Alicante. Estaciones meteorológicas utilizadas y promedio de las temperaturas mínimas durante el periodo de canícula (15 julio – 15 agosto) en el periodo temporal (2017-2020). Recuadro de color azul representa los ámbitos con los píxeles más fríos, mientras, que, los recuadros negros hacen referencia a las temperaturas más elevadas. Fuente: Elaboración propia. Cuadro 8. Cuadro estadístico de las temperaturas mínimas diarias registradas en las estaciones meteorológicas de la provincia de Alicante durante las canículas (15 julio – 15 agosto) de 2017 a 2020. Estación meteorológica Temperatura media de las mínimas (oC) Temperatura mínima absoluta (oC) Temperatura máxima de las mínimas (oC) Crevillente 22,7 19,4 27,3 Alicante 22,6 17,5 27,2 Elche EFA 22,4 19,6 25,9 Villajoyosa 22,3 17,8 27,4 Altea 22,2 17,3 27,1 Agost 19,4 15,1 23,6 Planes 18,3 11,7 22,2 Pinoso 18,2 14,3 22,1 Camp de Mirra 18,1 13,4 21,6 Villena 17,2 12,6 20,9 Fuente: Elaboración propia. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 28 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Gráfico 8. Diagrama de cajas o boxplot de las temperaturas mínimas registradas en cuatro estaciones meteorológicas de la provincia de Alicante durante cuatro periodos caniculares (2017-2020). Fuente: Elaboración propia. 4. DISCUSIÓN Y CONCLUSIONES La necesidad de cartografiar los fenómenos de inversión térmica nocturna (ITN) y la distribución de las temperaturas en un territorio, ha llevado al desarrollo y proliferación de las imágenes de satélite a través de las termografías. Debido a la complejidad de la distribución espacial de temperaturas nocturnas durante madrugadas estables (conformación de CAP en fondos de valle), y a la imposibilidad de mantener una red de estaciones meteorológicas lo suficientemente densa para identificar bien el fenómeno, se utilizan las termografías nocturnas o distribución térmica a través de bandas satelitales. La validación a la que han sido sometidos diferentes productos satelitales (entre los que se encuentran Band I5) muestran que la temperatura mínima del aire puede recuperarse de forma óptima con los productos de superficie terrestre LST (Zhu et al., 2013). En ese sentido, los excelentes resultados estadísticos en la validación del producto Band I5 (0,87) en los últimos años, a excepción de las ligeras discrepancias en funciones espectrales del sensor, supone una excelente relación entre las radiaciones del sensor medidas por VIIRS y la estimación de temperaturas mínimas, resultados similares a los arrojados en el oeste de EE.UU (Guillevic et al., 2014). El producto Band I5, tras los excelentes ajustes estadísticos obtenidos en diferentes investigaciones, ha sido especialmente utilizado en los últimos años para la estimación de la temperatura superficial en diferentes regiones del planeta. Los resultados arrojados en la presente investigación, con una correlación estadística de R2 de 0,87 y diferencia térmica de -0,22oC se encuentran en la horquilla de las principales investigaciones realizadas, parecidas, por ejemplo, a las de Neteler (2010), quien determina correlaciones entre 0,70 y 0,98, y diferencias térmicas de 0,5 ºC en los Alpes centro-orientales (noreste de Italia). Los resultados también indican una precisión muy razonable de los productos actuales de VIIRS LST (Band I5), con sesgos promedio de 0,36 K y –0,58 K y RMSE promedio de 2.7 K y 1.5 K durante el día y la noche respectivamente (Niclòs et al., 2018; Liu et al., 2015). Se trata de resultados muy ajustados a los obtenidos en el área de estudio, con un RMSE promedio de 1,83 y Bias de -0,22oC. Es preciso considerar que la mayoría de la población se ubica en torno a las capitales de provincia, ejes litorales y prelitorales, donde, por lo general, las temperaturas mínimas estivales son más elevadas. Las David Espín-Sánchez / REA N. 42 (2021) 8 - 33 29 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 capitales de provincia de Granada, Jaén o Albacete, ubicadas en zonas de interior, se ven fuertemente influenciadas por las ICUs nocturnas, con un importante incremento térmico respecto a su ámbito periurbano. Se sitúan aquí la gran mayoría de estaciones meteorológicas, por lo que los sectores con los valores promedios más elevados se encuentran bien representadas. La principal dificultad radica en la observación de áreas despobladas de interior, valles montañosos poco habitados y media/alta montaña. Aunque no ha sido posible la instalación de sensores térmicos registradores en la totalidad de la provincia, se han utilizado datalogger ubicados en los principales polos fríos de las provincias de Murcia, Jaén, Albacete y Almería (Espín Sánchez et al., 2018b). El uso de las redes propias de termómetros registradores (datalogger) ha sido ampliamente utilizada en ámbitos montañosos desprovistos de series meteorológicas. En la península, se han implementado principalmente en Picos de Europa (Cordillera Cantábrica) o Pirineos o en la plataforma kárstica del Pale di San Martino (Italia) (Renon 2011; Jiménez et al., 2015; Miró et al., 2017; Pagès et al., 2017; Iglesias González et al., 2018), válidos para multitud de investigaciones donde el análisis de temperaturas en ambientes microclimáticos es clave para el estudio de vegetación, gradientes térmicos, formación y disipación de piscinas de aire frío (CAP) y validación de imágenes de satélite (termografías). En los últimos años, la teledetección se ha convertido en una herramienta metodológica de primer nivel para muchos campos de la climatología. La temperatura de la superficie terrestre (LST) se considera una de las formas más adecuadas para recuperar la temperatura de la superficie del aire (Ta) (Phan et al., 2019). Se ha utilizado, durante los últimos años, para evaluar el potencial TVX (Temperatura-Vegetación Index) aplicado a las imágenes térmicas MODIS con el fin de estimar las temperaturas mínimas y máximas diarias del aire en el Tibet (Pepin et al., 2019); o para el estudio de las islas de calor urbana (ICU), por ejemplo, en Madrid (Allende et al., 2018). El análisis efectuado en el área de estudio durante la época canicular muestras notables resultados en las provincias analizadas, con una diferencia absoluta de las temperaturas medias de las mínimas de 15,1oC entre la Nava del Polvo (8,4oC) y Los Belones (23,5oC) separados por 182,2 km en línea recta. Y es que, en los últimos cuatro periodos caniculares (2016-2020) y un total de 128 jornadas de análisis, el Área Metropolitana de Murcia, próxima al centro de la ciudad (CMT AEMET), registra el 100% de las madrugadas tropicales, influenciada, sin duda, por la ICU nocturna. Espín et al., (2018) determinan notables diferencias térmicas en la Vega Media del río Segura (Murcia), concretamente entre La Contraparada y el CMT de AEMET (en resalte topográfico), con un desnivel altitudinal de 4 m y 4,2 km de distancia. La tasa de cambio (lapse rate o LR) media es de 80,2oC/100 m, con valores absolutos de 192,9oC/100 m el 06/11/2013 (temperatura mínima absoluta de 16,0oC en el CMT AEMET y 8,3oC en la Contraparada). En el mismo ámbito temporal, las temperaturas se mantienen próximas a la helada, con valores mínimos absolutos de 1,0oC en la Sierra de la Cabrilla (Nava del Polvo, Jaén). La diferencia diaria más importante en temperaturas mínimas absolutas se registra el 29/07/2018, con 23,4oC, sin duda un valor realmente destacado. Durante la citada jornada se registra una temperatura mínima de 1,2oC en el paraje de Monterilla (Jaén), en el Calar de Hernán Perea, bajo un intenso proceso de ITN. En el mismo ámbito anteriormente descrito (Los Belones) la temperatura mínima no desciende de 24,6oC. En el sondeo atmosférico de Murcia (00 h UTC) se registra una temperatura de 21,7oC en el geopotencial de 850 hPa, por lo que la capacidad de inversión (CINV) es realmente importante en Hernán Perea (Espín Sánchez et al., 2018). En general, el aire seco y templado en el nivel gepotencial de 850 hPa favorece las intensas inversiones nocturnas. (Vihma et al., 2011), e incrementa los contrastes térmicos entre zonas propensas a ITN y resaltes orográficamente accidentados. En definitiva, los intensos procesos de ITN en valles de elevada altitud de zonas del interior sureste peninsular, generan importantes diferencias térmicas con las zonas litorales, prelitorales, y medianías orográficamente accidentadas, donde las temperaturas mínimas estivales no suelen descender en la mayoría de ocasiones de 20 ºC. La utilización de termómetros registradores en éstas zonas montañosas ha sido capaz de descubrir temperaturas mínimas inusualmente bajas en pleno periodo canicular en calares, ámbitos remotos no habitados, o valles de montaña, con temperaturas mínimas próximas a 0oC. Incluso en zonas habitadas como La Risca o Benablón (Murcia), Cañadas de Nerpio (Albacete), Cañada Grande (Almería), Guadix (Granada) o Pontones (Jaén), las temperaturas mínimas descienden de 10oC, por lo que se trata de una recurrencia importante en sólo 4 periodos caniculares de observación. David Espín-Sánchez / REA N. 42 (2021) 8 - 33 30 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Agradecimientos Especial agradecimiento a la Agencia Estatal de Meteorología (AEMET) en su Delegación Territorial de Murcia, por la cesión de dos garitas oficiales de tipo pagoda, y por la calibración de dos termómetros registradores. Declaración responsable y conflicto de intereses El autor declara que no existe ningún conflicto de interés con relación a la publicación de este artículo. REFERENCIAS Allende, F. A., Fernández, F. G., Rasilla, D.A. & Alcaide, J.M. (2018). Isla de calor nocturna estival y confort térmico en Madrid: avance para un planeamiento térmico en áreas urbanas. Ciudad y Territorio Estudios Territoriales (CyTET), 50(195), 101-120. https://doi.org/10.37230/CyTET.2021.M21 Anderson, M.C., Kustas, W.P., Norman, J.M., Hain, C.R., Mecikalski, J.R., Schultz, L. & Gao, F. (2011). Mapping daily evapotranspiration at field to continental scales using geostationary and polar orbiting satellite imagery. Hydrology and Earth System Sciences, 15(1), 223-239. https://doi.org/10.5194/hess-15-223-2011 Arduini, G., Chemel, C. & Staquet, C. (2020). Local and non-local controls on a persistent cold-air pool in the Arve River Valley. Quarterly Journal of the Royal Meteorological Society. https://doi.org/10.1002/qj.3776 Bolstad, P.V., Swift, L., Collins, F. & Régnière, J. (1998). Measured and predicted air temperatures at basin to regional scales in the southern Appalachian Mountains. Agricultural and Forest Meteorology, 91(3-4), 161-176. https://doi. org/10.1016/S0168-1923(98)00076-8 Cao, C., Xiong, X., Wolfe, R., De Luccia, F., Liu, Q., Blonski, S.,... & Hillger, D. (2013). Visible infrared imaging radiometer suite (VIIRS) sensor data record (SDR) user’s guide. NOAA Technical Report NESDIS: College Park, MD, USA. Chai, T. & Draxler, R.R. (2014). Root mean square error (RMSE) or mean absolute error (MAE). Geoscientific Model Development Discussions, 7(1), 1525-1534. https://doi.org/10.5194/gmd-7-1247-2014 Clements, C.B., Whiteman, C.D. & Horel J.D. (2003). Cold-air-pool structure and evolution in a mountain basin: Peter Sinks, Utah. Journal of Applied Meteorology, 42, 752-768. https://doi.org/10.1175/1520-0450(2003)042<0752:CSAEIA>2.0.CO;2 Daly, CH. Conklin, D.R. & Unsworth, M.H. (2009). Local atmospheric decoupling in complex topography alters climate change impacts. International Journal of Climatology, 30(12), pp. 1857-1864. https://doi.org/10.1002/joc.2007 Dodson, R. & Marks, D. (1997). Daily air temperature interpolated at high spatial resolution over a large mountainous region. Climate Research, 8(1), 1-20. https://doi.org/10.3354/cr008001 Duan, S.B., Li, Z.L., Li, H., Göttsche, F.M., Wu, H., Zhao, W. & Coll, C. (2019). Validation of Collection 6 MODIS land surface temperature product using in situ measurements. Remote sensing of environment, 225, 16-29. https://doi.org/10.1016/j.rse.2019.02.020 Eisenbach, S., Pospichal, B., Whiteman, C.D., Steinacker, R. & Dorninger, M. (2003). Classification of cold air pool events in the Gstettneralm, a sinkhole in the Eastern Alps. [Comunicación en congreso]. Extended Abstracts, Int. Conf. on Alpine Meteorology and MAP-Meeting, Brig, Switzerland, MeteoSwiss, Publication 66, 157–160. Espín Sánchez, D. Conesa García C. & Castejón Porcel, G. (2018a). Temperature inversions due to warm air advections at low levels: significant thermal constrast in the Vega Media of the Segura River (Southeast Spain)”. In: Advances in Environmental Research Publisher: Nova Science Publishers (New York). http://doi.org/10.3989/estgeogr.201416 Espín Sánchez, D., Ruiz Álvarez, V., Martí Talavera, J. & García Marín, R. (2018b). Estudio preliminar de las inversiones térmicas en el sureste de la Península Ibérica: El caso de los campos de Hernán Perea. Pirineos, 173, e036. https:// doi.org/10.3989/pirineos.2018.173003 Guillevic, P.C., Privette, J.L., Coudert, B., Palecki, M.A., Demarty, J., Ottle, C. & Augustine, J.A. (2012). Land Surface Temperature product validation using NOAA’s surface climate observation networks—Scaling methodology for the Visible Infrared Imager Radiometer Suite (VIIRS). Remote Sensing of Environment, 124, 282-298. https://doi.org/10.1016/j.rse.2012.05.004 David Espín-Sánchez / REA N. 42 (2021) 8 - 33 31 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Guillevic, P.C., Biard, J.C., Hulley, G.C., Privette, J. L., Hook, S.J., Olioso, A.,... & Csiszar, I. (2014). Validation of Land Surface Temperature products derived from the Visible Infrared Imaging Radiometer Suite (VIIRS) using ground-based and heritage satellite measurements. Remote Sensing of Environment, 154, 19-37. https://doi.org/10.1016/j.rse.2014.08.013 Gustavsson, T.; Karlsson, M.; Bogren J. & Lindqvist, S. (1998). Development of temperature patterns during clear nights. Journal of Applied Meteorology, 37, pp. 559-571. https://doi.org/10.1175/1520-0450(1998)037<0559:DOTPDC>2.0.CO;2 Hillger D. (2013). First-light imagery from Suomi NPP VIIRS. Bull. Amer. Meteorol. Society, 94, 1019–1029. https://doi. org/10.1175/BAMS-D-12-00097.1 Hulley, G.C., Malakar, N.K., Islam, T. & Freepartner, R.J. (2017). NASA’s MODIS and VIIRS Land Surface Temperature and Emissivity Products: A Long-Term and Consistent Earth System Data Record. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 11(2), 522-535. https://10.1109/JSTARS.2017.2779330 Iglesias González, M., Acuña, J.L, García H., Rodríguez, A., Pajares, S., Rodríguez, J., Ruiz-Verdú, A. & Jesús Pérez H. (2018): Proyecto Jous: temperaturas mínimas absolutas en la cordillera cantábrica y su relación con las piscinas de aire frío. [Comunicación en congreso]. XXXV Jornadas científicas de la AME - 19º encuentro hispano-luso de meteorología. https:// doi.org/10.30859/ameJrCn35 Jarvis, C.H. & Stuart, N. (2001). A comparison among strategies for interpolating maximum and minimum daily air temperatures. Part II: The interaction between number of guiding variables and the type of interpolation method. Journal of Applied Meteorology, 40(6), 1075-1084. https://doi.org/10.1175/1520-0450(2001)040<1075:ACASFI>2.0.CO;2 Jiménez, M.A., Ruiz, A. & Cuxart, J. (2015). Estimation of cold pool areas and chilling hours through satellite-derived surface temperatures. Agricultural and Forest Meteorology, 207, 58-68. https://doi.org/10.1016/j.agrformet.2015.03.017 Li, H., Sun, D., Yu, Y., Wang, H., Liu, Y., Liu, Q.,... & Cao, B. (2014). Evaluation of the VIIRS and MODIS LST products in an arid area of Northwest China. Remote Sensing of Environment, 142, 111-121. https://doi.org/10.1016/j.rse.2013.11.014 Lindkvist, L.; Gustavsson, T. & Bogren, J. (2000). A frost assessment method for mountainous areas. Agricultural and Forest Meteorology, 102, 51-67. https://doi.org/10.1016/S0168-1923(99)00087-8 Litschauer, D., (1962). Untersuchung der Entwicklung von Kaltluftseen in Dolinen- und Beckenlagen (Investigation of the development of cool air pools in sinkholes and basins) [Tesis Doctoral]. University of Vienna. Liu, Y., Yu, Y., Yu, P., Göttsche, F. & Trigo, I. (2015). Quality assessment of S-NPP VIIRS land surface temperature product. Remote Sensing, 7(9), 12215-12241. https://doi.org/10.3390/rs70912215 Lundquist J.D. & Huggett B. (2008). Evergreen trees as inexpensive radiation shields for temperature sensors. Water Resour. Res. 44: W00D04. https://doi.org/10.1029/2008WR006979 McKee, T.B. & O’Neal, R.O. (1989). The role of valley geometry and energy budget in the formation of nocturnal valley winds. Journal of Applied Meteorology, 28, 445-456. https://doi.org/10.1175/1520-0450(1989)028<0445:TROVGA>2.0.CO;2 Minnett, P.J., Evans, R.H., Podestá, G.P. & Kilpatrick, K.A. (2014). Sea-Surface Temperature from Suomi-NPP VIIRS: Algorithm development and uncertainty estimation. [Comunicación en congreso] Ocean Sensing and Monitoring VI. International Society for Optics and Photonics (9111, p. 91110C). Miró, J.R., Peña, J. C., Pepin, N., Sairouni, A. & Aran, M. (2017). Key features of cold-air pool episodes in the northeast of the Iberian Peninsula (Cerdanya, eastern Pyrenees). International Journal of Climatology, 37. https://doi.org/10.1002/joc.5236 Mosiño, A.P. & García, E. (1966). Evaluación de la sequía intraestival en la República Mexicana. Proc. Conf. Reg. Latinoamericana Unión Geogr. Int, 3, 500-516. Nalder, I.A. & Wein, R.W. (1998). Spatial interpolation of climatic normals: test of a new method in the Canadian boreal forest. Agricultural and forest meteorology, 92(4), 211-225. https://doi.org/10.1016/S0168-1923(98)00102-6 Neteler, M. (2010). Estimating daily land surface temperatures in mountainous environments by reconstructed MODIS LST data. Remote sensing, 2(1), 333-351. https://doi.org/10.3390/rs1020333 Niclòs, R., Pérez-Planells, L., Coll, C., Valiente, J.A. & Valor, E. (2018). Evaluation of the S-NPP VIIRS land surface temperature product using ground data acquired by an autonomous system at a rice paddy. ISPRS Journal of Photogrammetry and Remote Sensing, 135, 1-12. https://doi.org/10.1016/j.isprsjprs.2017.10.017 Noi, P., Kappas, M. & Degener, J. (2016). Estimating daily maximum and minimum land air surface temperature using MODIS land surface temperature data and ground truth data in Northern Vietnam. Remote Sensing, 8(12), 1002. https://doi.org/10.3390/rs8121002 Pagès, M., Pepin, N. & Miró, J.R. (2017). Measurement and modelling of temperature cold pools in the Cerdanya valley (Pyrenees), Spain. Meteorological Applications, 2(24), 290-302 https://doi.org/10.1002/met.1630 David Espín-Sánchez / REA N. 42 (2021) 8 - 33 32 © Editorial Universidad de Sevilla 2021 | Sevilla, España| CC BY-NC-ND 4.0 | e-ISSN: 2340-2776 | doi: https://dx.doi.org/10.12795/rea.2021.i42.01 Pepin, N., Deng, H., Zhang, H., Zhang, F., Kang, S. & Yao, T. (2019). An examination of temperature trends at high elevations across the Tibetan Plateau: The use of MODIS LST to understand patterns of elevation-dependent warming. Journal of Geophysical Research: Atmospheres, 124(11), 5738-5756. https://doi.org/10.1029/2018JD029798 Phan, T.N., Kappas, M., Nguyen, K.T., Tran, T.P., Tran, Q.V. & Emam, A.R. (2019). Evaluation of MODIS land surface temperature products for daily air surface temperature estimation in northwest Vietnam. International Journal of Remote Sensing, 40(14), 5544-5562. https://doi.org/10.1080/01431161.2019.1580789 Pospichal, B., Eisenbach, S., Whiteman, C.D., Steinacker, R. & Dorninger, M. (2003). Observations of the cold air outflow from a basin cold pool through a low pass. [Comunicación en congreso]. In: Ext. Abstr., Vol. A, Intl Conf. Alpine Meteor. and the MAP-Meeting, 153-156. Rajasekar, U. & Weng, Q. (2009). Urban heat island monitoring and analysis using a non-parametric model: A case study of Indianapolis. ISPRS Journal of Photogrammetry and Remote Sensing, 64(1), 86-96. https://doi.org/10.1016/j.isprsjprs.2008.05.002 Renon, B. (2011). Le fabbriche naturali del freddo. Dipartimento Regionale per la Sicurezza del Territorio. ARPAV Sauberer, F. & Dirmhirn, I. (1954). Über die Entstehung der extremen Temperaturminima in der Doline Gstettner-Alm”. Archiv für Meteorologie, Geophysik und Bioklimatologie, Serie B, 5(3-4), 307-326. https://doi.org/10.1007/BF02242757 Sauberer, F. & Dirmhirn, I. (1956). Weitere Untersuchungen über die kaltluftansammungen in der Doline Gstettner-Alm bei Lunz im Niederösterreich (Further investigations of the cold air buildup in the Gstettner-Alm doline near Lunz in lower Austria). Wetter Leben, 8, 187-196. Schroeder, W. & Giglio, L. (2017). VIIRS/NPP Thermal Anomalies/Fire 6-Min L2 Swath 750m V001. NASA EOSDIS Land Processes DAAC. https://doi.org/10.5067/VIIRS/VNP14.001 Shapiro, S.S. (1998). Distribution assessment, Statistics (17): Applications (Eds.N. Balakrishnan and C. R. Rao), pp. 475–494. Elsevier. Sheridan, P.F. (2019). Synoptic-flow interaction with valley cold-air pools and effects on cold-air pool persistence: Influence of valley size and atmospheric stability. Quarterly Journal of the Royal Meteorological Society, 145(721), 1636-1659. https://doi.org/10.1002/qj.3517 Sismanidis, P., Keramitsoglou, I., Bechtel, B. & Kiranoudis, C. (2016). Improving the downscaling of diurnal land surface temperatures using the annual cycle parameters as disaggregation kernels. Remote Sensing, 9(1), 23. https://doi. org/10.3390/rs9010023 Sun, D. & Pinker, R.T. (2003). Estimation of land surface temperature from a Geostationary Operational Environmental Satellite (GOES-8). Journal of Geophysical Research: Atmospheres, 108(D11). https://doi.org/10.1029/2002JD002422 Thornton, P.E., Running, S.W. & White, M.A. (1997). Generating surfaces of daily meteorological variables over large regions of complex terrain. Journal of Hydrology, 190(3-4), 214-251. https://doi.org/10.1016/S0022-1694(96)03128-9 Vancutsem, C., Ceccato, P., Dinku, T. & Connor, S.J. (2010). Evaluation of MODIS land surface temperature data to estimate air temperature in different ecosystems over Africa. Remote Sensing of Environment, 114(2), 449-465. https://doi. org/10.1016/j.rse.2009.10.002 Vihma, T., Kilpeläinen, T., Manninen, M., Sjöblom, A., Jakobson, E., Palo, T. & Maturilli, M. (2011). Characteristics of temperature and humidity inversions and low-level jets over Svalbard fjords in spring. Advances in Meteorology, 2011. https://doi.org/10.1155/2011/486807 Voinov, Vassily G., Nikulin & Mikhail S. (1993). Unbiased estimators and their applications. Univariate case. Dordrect: Kluwer Academic Publishers. Wan, Z., Hook, S. & Hulley, G. (2015). MYD11A1 MODIS/Aqua Land Surface Temperature/Emissivity Daily L3 Global 1km SIN Grid V006. NASA EOSDIS Land Processes DAAC. https://10.5067/MODIS/MYD11A1.006 Whiteman, C.D. (2000). Mountain Meteorology: Fundamentals and Applications. Oxford University Press. Whiteman, C.D., Pospichal, B., Eisenbach, S., Weihs, P., Clements, C.B., Steinacker, R. & Dorninger, M., (2004). Inversion breakup in small Rocky Mountain and Alpine basins. Journal of Applied Meteorology, 43(8), 1069-1082. https://doi. org/10.1175/1520-0450(2004)043<1069:IBISRM>2.0.CO;2 Whiteman, C.D., Hoch, S.W., Hahnenberge R.M., Muschinski, A., Hohreiter, V., Behn, M. & Clements, C.B., (2008). METCRAX 2006: Meteorological experiments in arizona’s meteor crater. Bulletin of the American Meteorological Society, 89(11), 1665-1680. https://doi.org/10.1175/2008BAMS2574.1