Improving Satellite Leaf Area Index Estimation Based On Various Integration Methods

Leaf Area Index (LAI) is an important land surface biophysical variable that is used to characterize vegetation amount and activity. Current satellite LAI products, however, do not satisfy the requirements of the modeling community due to their large uncertainties and frequent missing values. Each LAI product is currently generated from only one satellite sensor data. There is an urgent need for advanced methods to integrate multiple LAI products to improve the product's accuracy and integrality for various applications. To meet this need, this study proposes four methods, including the Optimal Interpolation (OI), Bayesian Maximum Entropy (BME), Multi-Resolution Tree (MRT) and Empirical Orthogonal Function (EOF), to integrate multiple LAI products. Three LAI products have been considered in this study: Moderate Resolution Imaging Spectroradiometer (MODIS), Multi-angle Imaging SpectroRadiometer (MISR) and Carbon cYcle and Change in Land Observational Products from an Ensemble of Satellites (CYCLOPES) LAI. As the basis of data integration, this dissertation first validates and intercompares MODIS and CYCLOPES LAI products and also evaluates their geometric accuracies. The CYCLOPES LAI product has smoother temporal profiles and fewer spatial variations, but tends to produce spurious large errors in winter. The Locally Adjusted Cubic-spline Capping algorithm is revised to smooth multiple years' average and variance. Although OI, BME and MRT based methods have been used in other fields, this is the first research to employ them in integrating multiple LAI products. This dissertation also presents a new integration method based on EOF to solve the problem of large data volume and inconsistent temporal resolution of different datasets. High resolution LAI reference maps generated with ground measurements are used to validate these algorithms. Validation results show that all of these four methods can fill data gaps and reduce the errors of the existing LAI products. The data gaps are filled with information from adjacent pixels and background. These algorithms remove the spurious large temporal and spatial variation of the original LAI products. The combination of multiple satellite products significantly reduces bias. OI and BME can reduce the RMSE from 1.0 (MODIS) to 0.7 and reduce the bias from +0.3 (MODIS) and -0.2 (CYCLOPES) to -0.1. MRT can produce similar results with OI but with significantly improved efficiency. EOF also generates the results with the RMSE of 0.7 but zero bias. Limited ground measurement data hardly prove which methods outperform the others. OI and BME theoretically produce statistically optimal results. BME relaxes OI's linear and Gaussian assumption and explicitly considers data error, but bears a much higher computational burden. MRT has improved efficiency but needs strict assumptions on the scale transfer function. EOF requires simpler model identification, while it is more "empirical" than "statistical". The original contributions of this study mainly include: 1) a new application of several different integration methods to incorporate multiple satellite LAI products to reduce uncertainties and improve integrality, 2) an enhancement of the Locally Adjusted Cubic-spline Capping by revising the end condition, 3) a novel comprehensive comparison of MODIS C5 LAI product with other satellite products, 4) the development of a new LAI normalization scheme by assuming the linear relationship between measurement error and LAI natural variance to account for the inconsistency between products, and finally, 5) the creation of a new data integration method based on EOF

Potential of using remote sensing techniques for global assessment of water footprint of crops

Remote sensing has long been a useful tool in global applications, since it provides physically-based, worldwide, and consistent spatial information. This paper discusses the potential of using these techniques in the research field of water management, particularly for ‘Water Footprint’ (WF) studies. The WF of a crop is defined as the volume of water consumed for its production, where green and blue WF stand for rain and irrigation water usage, respectively. In this paper evapotranspiration, precipitation, water storage, runoff and land use are identified as key variables to potentially be estimated by remote sensing and used for WF assessment. A mass water balance is proposed to calculate the volume of irrigation applied, and green and blue WF are obtained from the green and blue evapotranspiration components. The source of remote sensing data is described and a simplified example is included, which uses evapotranspiration estimates from the geostationary satellite Meteosat 9 and precipitation estimates obtained with the Climatic Prediction Center Morphing Technique (CMORPH). The combination of data in this approach brings several limitations with respect to discrepancies in spatial and temporal resolution and data availability, which are discussed in detail. This work provides new tools for global WF assessment and represents an innovative approach to global irrigation mapping, enabling the estimation of green and blue water use

Estimating the fraction of absorbed photosynthetically active radiation from multiple satellite data

The fraction of absorbed photosynthetically active radiation (FAPAR) is a critical input parameter in many climate and ecological models. The accuracy of satellite FAPAR products directly influences estimates of ecosystem productivity and carbon stocks. The targeted accuracy of FAPAR products is 10%, or 0.05, for many applications. This study evaluates satellite FAPAR products, presents a new FAPAR estimation model and develops data fusion schemes to improve the FAPAR accuracy. Five global FAPAR products, namely MODIS, MISR, MERIS, SeaWiFS, and GEOV1 were intercompared over different land covers and directly validated with ground measurements at VAlidation of Land European Remote sensing Instruments (VALERI) and AmeriFlux sites. Intercomparison results show that MODIS, MISR, and GEOV1 agree well with each other and so do MERIS and SeaWiFS, but the difference between these two groups can be as large as 0.1. The differences between the products are consistent throughout the year over most of the land cover types, except over the forests, because of the different assumptions in the retrieval algorithms and the differences between green and total FAPAR products over forests. Direct validation results show that the five FAPAR products have an uncertainty of 0.14 when validating with total FAPAR measurements, and 0.09 when validating with green FAPAR measurements. Overall, current FAPAR products are close to, but have not fulfilled, the accuracy requirement, and further improvements are still needed. A new FAPAR estimation model was developed based on the radiative transfer for horizontally homogeneous continuous canopy to improve the FAPAR accuracy. A spatially explicit parameterization of leaf canopy and soil background reflectance was derived from a thirteen years of MODIS albedo database. The new algorithm requires the input of leaf area index (LAI), which was estimated by a hybrid geometric optic-radiative transfer model suitable for both continuous and discrete vegetation canopies in this study. The FAPAR estimates by the new model was intercompared with reference satellite FAPAR products and validated with field measurements at the VALERI and AmeriFlux experimental sites. The validation results showed that the FAPAR estimates by the new method had slightly better performance than the MODIS and the MISR FAPAR products when using corresponding satellite LAI product values as input. The FAPAR estimates can be further improved with the LAI estimates from the presented model as input. The improvements are apparent at grasslands and forests with an 8% reduction of uncertainty. The new model can successfully identify the growing seasons and produce smooth time series curves of estimated FAPAR over years. The root mean square error (RMSE) was reduced from 0.16 to 0.11 for MODIS and from 0.18 to 0.1 for MISR overall. Application of the presented model at a regional scale generated consistent FAPAR maps at 30 m, 500 m, and 1100 m spatial resolutions from the Landsat, MODIS, and MISR data. As an alternative method to improve FAPAR accuracy, in addition to developing FAPAR estimation models, two data fusion schemes were applied to integrate multiple satellite FAPAR products at two scales: optimal interpolation at the site scale and multiple resolution tree at the regional scale. These two fusion schemes removed the bias and resulted in a 20% increase in the R2 and a 3% reduction in the RMSE as compared with the average of the individual FAPAR products. The regional scale fusion filled in the missing values and provided spatially consistent FAPAR distributions at different resolutions. The original contribution of this study is that multiple FAPAR products have been assessed with a comprehensive set of measurements from two field experiments at the global scale. This study improved the accuracy of FAPAR using a new model and local pixel based soil background and leaf canopy albedos. High FAPAR accuracy was achieved through integration at both the temporal and spatial domains. The improved accuracy of FAPAR values from this study by 5% would help to decrease an equal amount of uncertainty in the estimation of gross and net primary production and carbon fluxes

Validation of a forage production index (FPI) derived from MODIS fcover time-series using high-resolution satellite imagery: methodology, results and opportunities

An index-based insurance solution was developed to estimate and monitor near real-time forage production using the indicator Forage Production Index (FPI) as a surrogate of the grassland production. The FPI corresponds to the integral of the fraction of green vegetation cover derived from moderate spatial resolution time series images and was calculated at the 6 km x 6 km scale. An upscaled approach based on direct validation was used that compared FPI with field-collected biomass data and high spatial resolution (HR) time series images. The experimental site was located in the Lot and Aveyron departments of southwestern France. Data collected included biomass ground measurements from grassland plots at 28 farms for the years 2012, 2013 and 2014 and HR images covering the Lot department in 2013 (n = 26) and 2014 (n = 22). Direct comparison with ground-measured yield led to good accuracy (R-2 = 0.71 and RMSE = 14.5%). With indirect comparison, the relationship was still strong (R-2 ranging from 0.78 to 0.93) and informative. These results highlight the effect of disaggregation, the grassland sampling rate, and irregularity of image acquisition in the HR time series. In advance of Sentinel-2, this study provides valuable information on the strengths and weaknesses of a potential index-based insurance product from HR time series images

Remote Sensing of Biophysical Parameters

Vegetation plays an essential role in the study of the environment through plant respiration and photosynthesis. Therefore, the assessment of the current vegetation status is critical to modeling terrestrial ecosystems and energy cycles. Canopy structure (LAI, fCover, plant height, biomass, leaf angle distribution) and biochemical parameters (leaf pigmentation and water content) have been employed to assess vegetation status and its dynamics at scales ranging from kilometric to decametric spatial resolutions thanks to methods based on remote sensing (RS) data.Optical RS retrieval methods are based on the radiative transfer processes of sunlight in vegetation, determining the amount of radiation that is measured by passive sensors in the visible and infrared channels. The increased availability of active RS (radar and LiDAR) data has fostered their use in many applications for the analysis of land surface properties and processes, thanks to their insensitivity to weather conditions and the ability to exploit rich structural and texture information. Optical and radar data fusion and multi-sensor integration approaches are pressing topics, which could fully exploit the information conveyed by both the optical and microwave parts of the electromagnetic spectrum.This Special Issue reprint reviews the state of the art in biophysical parameters retrieval and its usage in a wide variety of applications (e.g., ecology, carbon cycle, agriculture, forestry and food security)

An Earth Observation Land Data Assimilation System (EO-LDAS)

Current methods for estimating vegetation parameters are generally sub-optimal in the way they exploit information and do not generally consider uncertainties. We look forward to a future where operational data assimilation schemes improve estimates by tracking land surface processes and exploiting multiple types of observations. Data assimilation schemes seek to combine observations and models in a statistically optimal way taking into account uncertainty in both, but have not yet been much exploited in this area. The EO-LDAS scheme and prototype, developed under ESA funding, is designed to exploit the anticipated wealth of data that will be available under GMES missions, such as the Sentinel family of satellites, to provide improved mapping of land surface biophysical parameters. This paper describes the EO-LDAS implementation, and explores some of its core functionality. EO-LDAS is a weak constraint variational data assimilation system. The prototype provides a mechanism for constraint based on a prior estimate of the state vector, a linear dynamic model, and Earth Observation data (top-of-canopy reflectance here). The observation operator is a non-linear optical radiative transfer model for a vegetation canopy with a soil lower boundary, operating over the range 400 to 2500 nm. Adjoint codes for all model and operator components are provided in the prototype by automatic differentiation of the computer codes. In this paper, EO-LDAS is applied to the problem of daily estimation of six of the parameters controlling the radiative transfer operator over the course of a year (> 2000 state vector elements). Zero and first order process model constraints are implemented and explored as the dynamic model. The assimilation estimates all state vector elements simultaneously. This is performed in the context of a typical Sentinel-2 MSI operating scenario, using synthetic MSI observations simulated with the observation operator, with uncertainties typical of those achieved by optical sensors supposed for the data. The experiments consider a baseline state vector estimation case where dynamic constraints are applied, and assess the impact of dynamic constraints on the a posteriori uncertainties. The results demonstrate that reductions in uncertainty by a factor of up to two might be obtained by applying the sorts of dynamic constraints used here. The hyperparameter (dynamic model uncertainty) required to control the assimilation are estimated by a cross-validation exercise. The result of the assimilation is seen to be robust to missing observations with quite large data gaps

Assimilation de données satellitaires pour le suivi des ressources en eau dans la zone Euro-Méditerranée

Une estimation plus précise de l'état des variables des surfaces terrestres est requise afin d'améliorer notre capacité à comprendre, suivre et prévoir le cycle hydrologique terrestre dans diverses régions du monde. En particulier, les zones méditerranéennes sont souvent caractérisées par un déficit en eau du sol affectant la croissance de la végétation. Les dernières simulations du GIEC (Groupe d'Experts Intergouvernemental sur l'Evolution du Climat) indiquent qu'une augmentation de la fréquence des sécheresses et des vagues de chaleur dans la région Euro-Méditerranée est probable. Il est donc crucial d'améliorer les outils et l'utilisation des observations permettant de caractériser la dynamique des processus des surfaces terrestres de cette région. Les modèles des surfaces terrestres ou LSMs (Land Surface Models) ont été développés dans le but de représenter ces processus à diverses échelles spatiales. Ils sont habituellement forçés par des données horaires de variables atmosphériques en point de grille, telles que la température et l'humidité de l'air, le rayonnement solaire et les précipitations. Alors que les LSMs sont des outils efficaces pour suivre de façon continue les conditions de surface, ils présentent encore des défauts provoqués par les erreurs dans les données de forçages, dans les valeurs des paramètres du modèle, par l'absence de représentation de certains processus, et par la mauvaise représentation des processus dans certaines régions et certaines saisons. Il est aussi possible de suivre les conditions de surface depuis l'espace et la modélisation des variables des surfaces terrestres peut être améliorée grâce à l'intégration dynamique de ces observations dans les LSMs. La télédétection spatiale micro-ondes à basse fréquence est particulièrement utile dans le contexte du suivi de ces variables à l'échelle globale ou continentale. Elle a l'avantage de pouvoir fournir des observations par tout-temps, de jour comme de nuit. Plusieurs produits utiles pour le suivi de la végétation et du cycle hydrologique sont déjà disponibles. Ils sont issus de radars en bande C tels que ASCAT (Advanced Scatterometer) ou Sentinel-1. L'assimilation de ces données dans un LSM permet leur intégration de façon cohérente avec la représentation des processus. Les résultats obtenus à partir de l'intégration de données satellitaires fournissent une estimation de l'état des variables des surfaces terrestres qui sont généralement de meilleure qualité que les simulations sans assimilation de données et que les données satellitaires elles-mêmes. L'objectif principal de ce travail de thèse a été d'améliorer la représentation des variables des surfaces terrestres reliées aux cycles de l'eau et du carbone dans le modèle ISBA grâce à l'assimilation d'observations de rétrodiffusion radar (sigma°) provenant de l'instrument ASCAT. Un opérateur d'observation capable de représenter les sigma° ASCAT à partir de variables simulées par le modèle ISBA a été développé. Une version du WCM (water cloud model) a été mise en œuvre avec succès sur la zone Euro-Méditerranée. Les valeurs simulées ont été comparées avec les observations satellitaires. Une quantification plus détaillée de l'impact de divers facteurs sur le signal a été faite sur le sud-ouest de la France. L'étude de l'impact de la tempête Klaus sur la forêt des Landes a montré que le WCM est capable de représenter un changement brutal de biomasse de la végétation. Le WCM est peu efficace sur les zones karstiques et sur les surfaces agricoles produisant du blé. Dans ce dernier cas, le problème semble provenir d'un décalage temporel entre l'épaisseur optique micro-ondes de la végétation et l'indice de surface foliaire de la végétation. Enfin, l'assimilation directe des sigma° ASCAT a été évaluée sur le sud-ouest de la France.More accurate estimates of land surface conditions are important for enhancing our ability to understand, monitor, and predict key variables of the terrestrial water cycle in various parts of the globe. In particular, the Mediterranean area is frequently characterized by a marked impact of the soil water deficit on vegetation growth. The latest IPCC (Intergovernmental Panel on Climate Change) simulations indicate that occurrence of droughts and warm spells in the Euro-Mediterranean region are likely to increase. It is therefore crucial to improve the ways of understanding, observing and simulating the dynamics of the land surface processes in the Euro-Mediterranean region. Land surface models (LSMs) have been developed for the purpose of representing the land surface processes at various spatial scales. They are usually forced by hourly gridded atmospheric variables such as air temperature, air humidity, solar radiation, precipitation, and are used to simulate land surface states and fluxes. While LSMs can provide a continuous monitoring of land surface conditions, they still show discrepancies due to forcing and parameter errors, missing processes and inadequate model physics for particular areas or seasons. It is also possible to observe the land surface conditions from space. The modelling of land surface variables can be improved through the dynamical integration of these observations into LSMs. Remote sensing observations are particularly useful in this context because they are able to address global and continental scales. Low frequency microwave remote sensing has advantages because it can provide regular observations in all-weather conditions and at either daytime or night-time. A number of satellite-derived products relevant to the hydrological and vegetation cycles are already available from C-band radars such as the Advanced Scatterometer (ASCAT) or Sentinel-1. Assimilating these data into LSMs permits their integration in the process representation in a consistent way. The results obtained from assimilating satellites products provide land surface variables estimates that are generally superior to the model estimates or satellite observations alone. The main objective of this thesis was to improve the representation of land surface variables linked to the terrestrial water and carbon cycles in the ISBA LSM through the assimilation of ASCAT backscatter (sigma°) observations. An observation operator capable of representing the ASCAT sigma° from the ISBA simulated variables was developed. A version of the water cloud model (WCM) was successfully implemented over the Euro-Mediterranean area. The simulated values were compared with those observed from space. A more detailed quantification of the influence of various factors on the signal was made over southwestern France. Focusing on the Klaus storm event in the Landes forest, it was shown that the WCM was able to represent abrupt changes in vegetation biomass. It was also found that the WCM had shortcomings over karstic areas and over wheat croplands. It was shown that the latter was related to a discrepancy between the seasonal cycle of microwave vegetation optical depth (VOD) and leaf area index (LAI). Finally, the direct assimilation of ASCAT sigma° observations was assessed over southwestern France

Development and Extrapolation of a General Light Use Efficiency Model for the Gross Primary Production

The global carbon cycle is one of the large biogeochemical cycles spanning all living and non-living compartments of the Earth system. Against the background of accelerating global change, the scientific community is highly interested in analyzing and understanding the dynamics of the global carbon cycle and its complex feedback mechanism with the terrestrial biosphere. The international network FLUXNET was established to serve this aim with measurement towers around the globe. The overarching objective of this thesis is to exploit the powerful combination of carbon flux measurements and satellite remote sensing in order to develop a simple but robust model for the gross primary production (GPP) of vegetation stands. Measurement data from FLUXNET sites as well as remote sensing data from the NASA sensor MODIS are exploited in a data-based model development approach. The well-established concept of light use efficiency is chosen as modeling framework. As a result, a novel gross primary production model is established to quantify the carbon uptake of forests and grasslands across a broad range of climate zones. Furthermore, an extrapolation scheme is derived, with which the model parameters calibrated at FLUXNET sites can be regionalized to pave the way for spatially continuous model applications

Improving data-oriented light use efficiency models of gross primary productivity with remotely sensed spectral indices

Several studies sustained the possibility that a photochemical reflectance index (PRI) directly obtained from satellite data can be used as a proxy for ecosystem light use efficiency (LUE) in diagnostic models of gross primary productivity. This modelling approach would avoid the complications that are involved in using me- teorological data as constraints for a fixed maximum LUE. However, no unifying model predicting LUE across climate zones and time based on MODIS PRI has been published to date. This study evaluates the effectiveness with which MODIS-based PRI can be used to estimate ecosystem light use efficiency at study sites of different plant functional types and vegetation densities. The objective is to examine if known limitations such as dependence on viewing and illumination geometry can be overcome and a single PRI-based model of LUE (i.e. based on the same reference band) can be applied under a wide range of conditions. Fur- thermore, this thesis examines the effect of using different fraction of absorbed photosynthetically active radiation (faPAR) products on the in-situ LUE used as ground truth and thus on the whole evaluation exercise. The conclusion of this study is that estimating LUE at site-level based on PRI reduces uncertainty compared to the approaches relying on a maximum LUE reduced by minimum temperature and vapour pressure deficit. Despite the advantages of using PRI to estimate LUE at site-level, a universally applicable light use efficiency model based on MODIS PRI could not be established. Models that were optimised for a pool of data from several sites did not perform well