Resumen de Estimación de flujos de agua entre suelo, vegetación y atmósfera mediante teledetección

Gorka Mendiguren González

• español

  En la frontera entre la superficie terrestre y la atmósfera se producen numerosos procesos físicos relacionados con el ciclo hidrológico. Cuando se producen precipitaciones en forma de lluvia, y el agua alcanza la superficie terrestre, una parte llega al suelo y otra parte puede ser interceptada por la vegetación. La fracción que llega al suelo se infiltra en la zona no saturada donde se almacena, lo humedece, disuelve los elementos que son absorbidos posteriormente por la vegetación y modifica las propiedades físicas del suelo. Para que la vegetación pueda desarrollarse es necesario que la planta abra los estomas, absorba CO2 y realice la fotosíntesis. Durante este proceso se produce una pérdida de agua a través de la hoja, que si es lo suficientemente grande puede llegar a hacer que la planta marchite si no es capaz de reponerla del suelo. El agua del suelo es devuelta a la atmósfera posteriormente mediante la evaporación y la transpiración de las plantas.

  La primera parte del trabajo se ha centrado en la estimación de parámetros biofísicos y estructurales de la vegetación, concretamente los relacionados con el contenido de agua. Para ello se han empleado numerosos datos recogidos en campo a lo largo de dos años fenológicos completos y se relacionaron con las medidas espectrales a dos escalas diferentes, campo y sensor MODIS (500 m). El contenido de agua se calculó usando tres métricas diferentes calculadas a partir de la misma muestra, el Contenido de Humedad de la Vegetación (FMC), el Espesor Equivalente de Agua (EWT) y el Contenido de Agua del Dosel (CWC). Además se usaron dos estimaciones a partir de Modelos de Transferencia Radiativa (RTM) para la obtención del FMC y CWC que fueron comparados con las obtenidas a partir de los modelos empíricos creados a partir los índices espectrales. Otras variables relacionadas como el contenido de materia seca (Dm) y el índice de área foliar (LAI) fueron también evaluadas usando índices de vegetación. Entre los resultados destacables de este estudio se encuentran en primer lugar los relacionados con el protocolo de recogida de datos en campo. En este estudio se obtuvieron evidencias de que las diferencias temporales a la hora de recoger datos en campo son más importantes que las diferencias espaciales en este ecosistema. Además se demostró la necesidad de mostrar consistencia en el protocolo de muestreo: tamaño de la muestra, hora de recogida de las muestras, etc. y en la importancia de evitar, en lo posible la toma de decisiones, generalmente subjetivas, por parte de los operadores de campo. Otro resultado destacable ha sido demostrar la existencia de una alta variabilidad del Dm a lo largo del año. Esto indica que asumir, como sugieren algunos autores, un valor constante de Dm para la estimación del espesor equivalente de agua a partir del contenido de humedad de la vegetación no es una opción viable en este ecosistema. De los índices de vegetación que fueron comparados en el estudio, el que presentó menores correlaciones fue el Índice de Vegetación Resistente a la Atmósfera (VARI). Se observaron algunas diferencias en el comportamiento de los modelos empíricos obtenidos con MODIS y las producidas a partir de medidas espectrales de campo, obteniendo resultados algo mejores en el caso de MODIS. Este hecho posiblemente sea debido a que las adquisiciones de del sensor MODIS presentan diferentes ángulos de observación, lo que reduce la proporción de suelo captada por el sensor y por lo tanto capturando una mayor proporción del dosel. La comparación entre los modelos empíricos y las estimaciones a partir de RTM demostró que en este caso los modelos empíricos aún mejoran las estimaciones de los modelos físicos desarrollados en zonas similares para estimar el contenido de humedad de la vegetación. La segunda parte del trabajo se ha centrado en la estimación del contenido de humedad del suelo combinando datos ópticos y térmicos mediante el cálculo del Índice de Temperatura y Sequedad de la Vegetación (TVDI) cuya obtención se basa en la técnica del triángulo. Se han investigado diferentes factores que afectan a la definición del triángulo, y cómo estos afectan los valores del TVDI y a su relación final con el contenido de humedad del suelo. En este trabajo se introdujo una modificación al cálculo del TVDI en la que se sustituyó el Índice de Vegetación de Diferencia Normalizada (NDVI) por el Índice de Diferencia Infraroja Normalizada (NDII). Esta modificación se tradujo en una mejora en el comportamiento de los modelos empíricos para estimar el contenido de humedad del suelo.

  Finalmente en la tesis se investiga el comportamiento de la EF en la zona de estudio y su estimación a partir de teledetección. El principal motivo del empleo de la EF es que ha sido ampliamente utilizada para estimar la evapotranspiración diaria, asumiendo que la EF es constante a lo largo del día. A partir de las medidas recogidas por una torre de flujos se han evaluado las variaciones diarias y se han validado las estimaciones de EF calculadas a partir de imágenes Landsat. Se ha usado una nueva versión modificada de la técnica del triángulo en la que se ha introducido el índice de área foliar adaptado a la escala Landsat a partir del producto MODIS (de 1 km a 30 m) como sustituto del índice de vegetación. Además se muestra un innovador método basado en las estadísticas propias del triángulo para la selección de las fechas a incluir en el análisis estadístico. La validación de las estimaciones de EF se ha llevado a cabo de dos maneras diferentes: usando las contribuciones de todos los pixeles incluidos en la zona de influencia de la torre; y utilizando el valor del único pixel correspondiente a la localización de la torre, mostrando ambas aproximaciones escasas diferencias en cuanto a resultados. Además se han comparado las EF diarias y la correspondiente a la hora de la pasada de Landsat sobre la zona de estudio. En este caso se observaron mayores diferencias, lo cual indica que el supuesto de una EF constante a lo largo del día ha de ser tomada con ciertas precauciones si el objetivo final es el cálculo de la evapotranspiración diaria.

• English

  In the boundary layer between vegetation and soil with the atmosphere different processes within the hydrologic cycle take place. When rainfall occurs, part of the water reaches the inland surface, part is infiltrated in the soil, and part is intercepted by vegetation. The fraction reaching the soil infiltrates in the unsaturated zone of the soil making it wet, dissolving the elements that will be absorbed by the plants, and changing the physical properties of the soil. Vegetation needs to open the stomata during the absorption of CO2 during photosynthesis. In this process loss of water by vegetation trough the leaves is produced. If this is loss is large, the wilting point of the plant can be reached, and therefore there is a need of water uptake which is done by the roots. Water returns to the atmosphere trough evaporation and plant transpiration. Evapotranspiration (ET), which can be defined as the process that combines the evaporation from the soil and the transpiration from the plants transferring water vapor to the atmosphere, is one the most relevant processes as it has strong influence on local climate, weather and many other biophysical processes. Finally, ET is an energy driven process that relies on land-surface energy fluxes like net radiation (Rn), sensible heat flux (H) and ground heat flux (G).

  It is possible to collect in situ different variables from vegetation, soil and atmosphere related to water fluxes in the boundary layer. However, these variables are difficult to measure due to its spatial and temporal variability. Remote sensing platforms on the other hand offer an alternative to in situ measurements.

  Nowadays, remote sensing offers a great variety of sensors with different spatial and temporal resolutions, which makes this tool very attractive to monitor these processes, although there is a need for accounting for the existing differences between remote sensing estimates and in situ measurements. This issue requires rigorous validations of remote sensing estimates as well as in situ sampling protocols.

  In this thesis, different variables related to water fluxes between vegetation, soil and atmosphere have been analyzed. The study focused on vegetation and soil water content, biophysical parameters of vegetation and in the behavior of Evaporative Fraction (EF), is defined as the ratio between ET and available energy (Rn-G). All data used were collected in a dehesa ecosystem at Las Majadas del Tiétar, in the province of Cáceres in Spain. In the study site there is an eddy covariance flux tower inside the FLUXNET network that has been operation since 2003.

  The first part of this work was focused on the estimation of biophysical and structural variables of the vegetation, being the study mainly oriented to estimate vegetation water content. A large data set collected over two full phenological years was used to relate different metrics of vegetation water content with spectral measurements at two different scales, using proximal sensing data and using MODIS surface reflectance data (500 m). Three different water content metrics simultaneously obtained from the same sample - Fuel Moisture Content (FMC), Equivalent Water Thickness (EWT) and Canopy Water Content (CWC) - were related to spectral vegetation indices (VI) calculated from MODIS and proximal sensing data. In addition estimates of FMC and CWC by inversion techniques of Radiative Transfer Model (RTM) were tested and compared against the empirical models. Dry Matter content (Dm) and Leaf Area Index (LAI) were also evaluated empirically with VI.

  As result from this study, several findings were revealed. The first one is related to the field data collection: results showed that temporal changes in FMC, EWT and CWC are more critical than their spatial variation within the MODIS pixel.

  The field protocols should hence be adapted in order sample more frequently rather than conducting extensive spatial samplings. Field data analysis also highlighted the need of being consistent with the size of the sample collected, and avoiding as much as possible subjective operator interpretations when collecting the data. Another relevant result from this study is that due to the high seasonal Dm variability, a constant annual value would not be recommended to predict EWT from FMC. Furthermore, Vegetation Atmospheric Resistant Index (VARI) provided the worst results in all cases. The empirical estimators differed between sensors, slightly better for MODIS than proximal sensing, probably due to differences in view angles, and as result the proportions of canopy observed by the sensor increases reducing the soil effect. These empirical methods still exceed RTM inversions developed for other sites to predict FMC and CWC.

  The second part focused on estimation of soil moisture (SM) combining optical and thermal sensors and evaluating the Temperature Vegetation Dryness Index (TVDI) calculated using the triangle method. Some of the factors involved in the triangle parameterization were studied to better understand how they affect the TVDI values and the implications in the final model performance for estimating SM.

  A modification was introduced in the TVDI calculation and the Normalized Difference Vegetation Index (NDVI) was substituted by the Normalized Difference Infrared Index (NDII). This modification translated in a better performance of empirical methods to estimate SM using this technique.

  Finally, this study investigates the behavior of EF on the site and the estimate using remote sensing. The reason why EF is investigated is that EF has been extensively used to retrieve daily ET, assuming that EF remains constant during daytime. Using the EF measured by the Eddy Covariance (EC), we have investigate how EF varies temporary and validated the EF calculated from Landsat using a modified triangle approach in which the VI was substituted by MODIS LAI downscaled to Landsat (from 1km to 30m). A novel method for selecting the highest quality days for the analysis based on the statistics of the triangle is also presented.

  The validation of EF estimates from remote sensing was carried out using either the pixel contribution based on the information of the EC footprint, or the single pixel located in the EC tower, showing both cases very similar results. Furthermore, daily EF and instantaneous EF at the time of the satellite overpass were also compared with the EF calculated from Landsat. Results of this comparison showed larger differences indicating that for this ecosystem the EF self-preservation assumption should be carefully taken into account if daily ET has to be obtained from EF.