System Engineering Analyses for the Study of Future Multispectral Land Imaging Satellite Sensors for Vegetation Monitoring

Vegetation monitoring is one of the key applications of earth observing systems. Landsat data have spatial resolution of 30 meters, moderate temporal coverage, and reasonable spectral sampling to capture key vegetation features. These characteristics of Landsat make it a good candidate for generating vegetation monitoring products. Recently, the next satellite in the Landsat series has been under consideration and different concepts have been proposed. In this research, we studied the impact on vegetation monitoring of two proposed potential design concepts: a wider field-of-view (FOV) instrument and the addition of red-edge spectral band(s). Three aspects were studied in this thesis: First, inspired by the potential wider FOV design, the impacts of a detector relative spectral response (RSR) central wavelength shift effect at high angles of incidence (AOI) on the radiance signal were studied and quantified. Results indicate: 1) the RSR shift effect is band-dependent and more significant in the green, red and SWIR 2 bands; 2) At high AOI, the impact of the RSR shift effect will exceed sensor noise specifications in all bands except the SWIR 1 band; and 3) The RSR shift will cause SWIR2 band more to be sensitive to atmospheric conditions. Second, also inspired by the potential wider FOV design, the impacts of the potential new wider angular observations on vegetation monitoring scientific products were studied. Both crop classification and biophysical quantity retrieval applications were studied using the simulation code DIRSIG and the canopy radiative transfer model PROSAIL. It should be noted that the RSR shift effect was also considered. Results show that for single view observation based analysis, the higher view angular observations have limited influence on both applications. However, for situations where two different angular observations are available potentially from two platforms, up to 4% improvement for crop classification and 2.9% improvement for leaf chlorophyll content retrieval were found. Third, to quantify the benefits of a potential new design with red-edge band(s), the impact of adding red-edge spectral band(s) in future Landsat instruments on agroecosystem leaf area index (LAI) and canopy chlorophyll content (CCC) retrieval were studied using a real dataset. Three major retrieval approaches were tested, results show that a potential new spectral band located between the Landsat-8 Operational Land Imager (OLI) red and NIR bands slightly improved the retrieval accuracy (LAI: R2 of 0.787 vs. 0.810 for empirical vegetation index regression approach, 0.806 vs. 0.828 for look-up-table inversion approach, and 0.925 vs. 0.933 for machine learning approach; CCC: R2 of 0.853 vs. 0.875 for empirical vegetation index regression approach, 0.500 vs. 0.570 for look-up-table inversion approach, and 0.854 vs. 0.887 for machine learning approach). In general, for the potential wider FOV design, the RSR shift effect was found to cause noticable radiance signal difference that is higher than detector noise in all OLI bands except SWIR1 band, which is not observed in the current OLI design with its 15 degree FOV. Also both the new wider angular observations and potential red-edge band(s) were found to slightly improve the vegetation monitoring product accuracy. In the future, the RSR shift effect in other optical designs should be evaluated since this study assumed the angle reaching the filter array is the same as the angle reaching the sensor. In addition to improve the accuracy of the off angle imaging study, a 3D vegetation geometry model should be explored for vegetation monitoring related studies instead of the 2D PROSAIL model used in this thesis

Multiscale Imaging of Evapotranspiration

Evapotranspiration (ET; evaporation + transpiration) is central to a wide range of biological, chemical, and physical processes in the Earth system. Accurate remote sensing of ET is challenging due to the interrelated and generally scale dependent nature of the physical factors which contribute to the process. The evaporation of water from porous media like sands and soils is an important subset of the complete ET problem. Chapter 1 presents a laboratory investigation into this question, examining the effects of grain size and composition on the evolution of drying sands. The effects of composition are found to be 2-5x greater than the effects of grain size, indicating that differences in heating caused by differences in reflectance may dominate hydrologic differences caused by grain size variation. In order to relate the results of Chapter 1 to the satellite image archive, however, the question of information loss between hyperspectral (measurements at 100s of wavelength intervals) laboratory measurements and multispectral (≤ 12 wavelength intervals) satellite images must be addressed. Chapter 2 focuses on this question as applied to substrate materials such as sediment, soil, rock, and non-photosynthetic vegetation. The results indicate that the continuum that is resolved by multispectral sensors is sufficient to resolve the gradient between sand-rich and clay-rich soils, and that this gradient is also a dominant feature in hyperspectral mixing spaces where the actual absorptions can be resolved. Multispectral measurements can be converted to biogeophysically relevant quantities using spectral mixture analysis (SMA). However, retrospective multitemporal analysis first requires cross-sensor calibration of the mixture model. Chapter 3 presents this calibration, allowing multispectral image data to be used interchangeably throughout the Landsat 4-8 archive. In addition, a theoretical explanation is advanced for the observed superior scaling properties of SMA-derived fraction images over spectral indices. The physical quantities estimated by the spectral mixture model are then compared to simultaneously imaged surface temperature, as well as to the derived parameters of ET Fraction and Moisture Availability. SMA-derived vegetation abundance is found to produce substantially more informative ET maps, and SMA-derived substrate fraction is found to yield a surprisingly strong linear relationship with surface temperature. These results provide context for agricultural applications. Chapter 5 investigates the question of mapping and monitoring rice agricultural using optical and thermal satellite image time series. Thermal image time series are found to produce more accurate maps of rice presence/absence, but optical image time series are found to produce more accurate maps of rice crop timing. Chapter 6 takes a more global approach, investigating the spatial structure of agricultural networks for a diverse set of landscapes. Surprisingly consistent scaling relations are found. These relations are assessed in the context of a network-based approach to land cover analysis, with potential implications for the scale dependence of ET estimates. In sum, this thesis present a novel approach to improving ET estimation based on a synthesis of complementary laboratory measurements, satellite image analysis, and field observations. Alone, each of these independent sources of information provides novel insights. Viewed together, these insights form the basis of a more accurate and complete geophysical understanding of the ET phenomenon

Multi-Season Phenology Mapping of Nile Delta Croplands Using Time Series of Sentinel-2 and Landsat 8 Green LAI

Space-based cropland phenology monitoring substantially assists agricultural managing practices and plays an important role in crop yield predictions. Multitemporal satellite observations allow analyzing vegetation seasonal dynamics over large areas by using vegetation indices or by deriving biophysical variables. The Nile Delta represents about half of all agricultural lands of Egypt. In this region, intensifying farming systems are predominant and multi-cropping rotations schemes are increasing, requiring a high temporal and spatial resolution monitoring for capturing successive crop growth cycles. This study presents a workflow for cropland phenology characterization and mapping based on time series of green Leaf Area Index (LAI) generated from NASA’s Harmonized Landsat 8 (L8) and Sentinel-2 (S2) surface reflectance dataset from 2016 to 2019. LAI time series were processed for each satellite dataset, which were used separately and combined to identify seasonal dynamics for a selection of crop types (wheat, clover, maize and rice). For the combination of L8 with S2 LAI products, we proposed two time series smoothing and fitting methods: (1) the Savitzky–Golay (SG) filter and (2) the Gaussian Processes Regression (GPR) fitting function. Single-sensor and L8-S2 combined LAI time series were used for the calculation of key crop Land Surface Phenology (LSP) metrics (start of season, end of season, length of season), whereby the detection of cropland growing seasons was based on two established threshold methods, i.e., a seasonal or a relative amplitude value. Overall, the developed phenology extraction scheme enabled identifying up to two successive crop cycles within a year, with a superior performance observed for the seasonal than for the relative threshold method, in terms of consistency and cropland season detection capability. Differences between the time series collections were analyzed by comparing the phenology metrics per crop type and year. Results suggest that L8-S2 combined LAI data streams with GPR led to a more precise detection of the start and end of growing seasons for most crop types, reaching an overall detection of 74% over the total planted crops versus 69% with S2 and 63% with L8 alone. Finally, the phenology mapping allowed us to evaluate the spatial and temporal evolution of the croplands over the agroecosystem in the Nile Delta.E.A. was supported by the predoctoral scholarship, grant number ACIF/2019/187, funded by the Generalitat Valenciana and co-funded by the European Social Fund. J.V. and S.B. were supported by the European Research Council (ERC) under the ERC-2017-STG SENTIFLEX project, grant number 755617. J.V. was additionally supported by a Ramón y Cajal Contract (Spanish Ministry of Science, Innovation and Universities). S.B. was additionally supported by the Generalitat Valenciana SEJIGENT program (SEJIGENT/2021/001) and European Union—NextGenerationEU (ZAMBRANO 21-04)

Spatio-temporal modelling of bluetongue virus distribution in Northern Australia based on remotely sensed bioclimatic variables

The presence of Bluetongue virus (BTV) in Northern Australia poses an ongoing threat for animal health and although clinical disease has not been detected in livestock, it limits export of livestock from the infected areas. BTV presence is governed by variable environmental conditions, which influence vector and host habitats. The National Arbovirus Monitoring Program (NAMP) was established to determine the extent of virus activity and control the risk of infection spread. Groups of young cattle, previously unexposed to infection, are regularly tested to detect evidence of transmission. This approach is labour and cost intensive and difficult to operate in the remote areas of Northern Australia. The resulting data are therefore characterised by spatial and temporal gaps. The aim of this research is to assess the use of remotely sensed environmental and climatic data as a means of predicting the distribution of BTV seroprevalence throughout Northern Australia to complement conventional surveillance.Environmental factors relating to the viruses’ host and vector habitats and the transmission cycle of BTV have been identified based on the extensive review of virus ecology. Different data sources have been assessed to provide sufficient spatial and temporal coverage for the definition of spatio-temporal environmental variables that can be used to explain and predict the distribution of BTV. Following this assessment, satellite data products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Tropical Rainfall Measuring Mission (TRMM) were acquired for the Pilbara in Western Australia, and the Northern Territory. These were reprojected and processed into spatio-temporal variables for the period between the years 2000 and 2009. Due to uncertainty in the precision of the geographic location and timing of animals tested for seropositivity, summary statistics of bioclimatic variables were generated at the station (i.e. property) level for each year. Different combinations of these variables, including vegetation greenness and phenology, land surface temperature and precipitation were screened for correlation with BTV presence using a Generalised Additive Model approach. A final model was developed to predict the presence or absence of BTV seropositivity on the basis of statistical significance of the remotely sensed predictor variables, and informed by knowledge of virus ecological principles.The model, based on the maximum seasonal Normalised Difference Vegetation Index (NDVI), and mean and maximum land surface temperature variables provided excellent discriminatory ability and the basis for the generation of prediction maps of BTV seropositivity for the first eight years. Besides internal assessment, the model’s predictive capabilities were validated using monitoring data from the season 2008/09.It has been demonstrated that the predictions are useful in complementing complement NAMP surveillance by identifying areas at higher risk for seropositivity in cattle, which aids planning of livestock movement and further monitoring activities. Uncertainty in the model was attributed to the spatio-temporal inconsistency in the precision of the available serosurveillance data. The discriminatory ability of models of this type could be further improved by ensuring that exact location details and date of NAMP BTV test events are consistently recorded