Advanced Processing of Multispectral Satellite Data for Detecting and Learning Knowledge-based Features of Planetary Surface Anomalies

abstract: The marked increase in the inflow of remotely sensed data from satellites have trans- formed the Earth and Space Sciences to a data rich domain creating a rich repository for domain experts to analyze. These observations shed light on a diverse array of disciplines ranging from monitoring Earth system components to planetary explo- ration by highlighting the expected trend and patterns in the data. However, the complexity of these patterns from local to global scales, coupled with the volume of this ever-growing repository necessitates advanced techniques to sequentially process the datasets to determine the underlying trends. Such techniques essentially model the observations to learn characteristic parameters of data-generating processes and highlight anomalous planetary surface observations to help domain scientists for making informed decisions. The primary challenge in defining such models arises due to the spatio-temporal variability of these processes. This dissertation introduces models of multispectral satellite observations that sequentially learn the expected trend from the data by extracting salient features of planetary surface observations. The main objectives are to learn the temporal variability for modeling dynamic processes and to build representations of features of interest that is learned over the lifespan of an instrument. The estimated model parameters are then exploited in detecting anomalies due to changes in land surface reflectance as well as novelties in planetary surface landforms. A model switching approach is proposed that allows the selection of the best matched representation given the observations that is designed to account for rate of time-variability in land surface. The estimated parameters are exploited to design a change detector, analyze the separability of change events, and form an expert-guided representation of planetary landforms for prioritizing the retrieval of scientifically relevant observations with both onboard and post-downlink applications.Dissertation/ThesisDoctoral Dissertation Computer Engineering 201

Vegetation responses to variations in climate: A combined ordinary differential equation and sequential Monte Carlo estimation approach

Vegetation responses to variation in climate are a current research priority in the context of accelerated shifts generated by climate change. However, the interactions between environmental and biological factors still represent one of the largest uncertainties in projections of future scenarios, since the relationship between drivers and ecosystem responses has a complex and nonlinear nature. We aimed to develop a model to study the vegetation’s primary productivity dynamic response to temporal variations in climatic conditions as measured by rainfall, temperature and radiation. Thus, we propose a new way to estimate the vegetation response to climate via a non-autonomous version of a classical growth curve, with a time-varying growth rate and carrying capacity parameters according to climate variables. With a Sequential Monte Carlo Estimation to account for complexities in the climate-vegetation relationship to minimize the number of parameters. The model was applied to six key sites identified in a previous study, consisting of different arid and semiarid rangelands from North Patagonia, Argentina. For each site, we selected the time series of MODIS NDVI, and climate data from ERA5 Copernicus hourly reanalysis from 2000 to 2021. After calculating the time series of the a posteriori distribution of parameters, we analyzed the explained capacity of the model in terms of the linear coefficient of determination and the parameters distribution variation. Results showed that most rangelands recorded changes in their sensitivity over time to climatic factors, but vegetation responses were heterogeneous and influenced by different drivers. Differences in this climate-vegetation relationship were recorded among different cases: (1) a marginal and decreasing sensitivity to temperature and radiation, respectively, but a high sensitivity to water availability; (2) high and increasing sensitivity to temperature and water availability, respectively; and (3) a case with an abrupt shift in vegetation dynamics driven by a progressively decreasing sensitivity to water availability, without any changes in the sensitivity either to temperature or radiation. Finally, we also found that the time scale, in which the ecosystem integrated the rainfall phenomenon in terms of the width of the window function used to convolve the rainfall series into a water availability variable, was also variable in time. This approach allows us to estimate the connection degree between ecosystem productivity and climatic variables. The capacity of the model to identify changes over time in the vegetation-climate relationship might inform decision-makers about ecological transitions and the differential impact of climatic drivers on ecosystems.Estación Experimental Agropecuaria BarilocheFil: Bruzzone, Octavio Augusto. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Bariloche; ArgentinaFil: Bruzzone, Octavio Augusto. Consejo Nacional de Investigaciones Cientificas y Tecnicas. Instituto de Investigaciones Forestales y Agropecuarias Bariloche; ArgentinaFil: Perri, Daiana Vanesa. Instituto Nacional de Tecnologia Agropecuaria (INTA). Estación Experimental Agropecuaria Bariloche. Área de Recursos Naturales; ArgentinaFil: Perri, Daiana Vanesa. Consejo Nacional de Investigaciones Cientificas y Tecnicas. Instituto de Investigaciones Forestales y Agropecuarias Bariloche; ArgentinaFil: Easdale, Marcos Horacio. Instituto Nacional de Tecnologia Agropecuaria (INTA). Estación Experimental Agropecuaria Bariloche. Área de Recursos Naturales; ArgentinaFil: Easdale, Marcos Horacio. Consejo Nacional de Investigaciones Cientificas y Tecnicas. Instituto de Investigaciones Forestales y Agropecuarias Bariloche; Argentin

Dynamical Approach for Real-Time Monitoring of Agricultural Crops

In this paper, a novel approach for exploiting multitemporal remote sensing data focused on real-time monitoring of agricultural crops is presented. The methodology is defined in a dynamical system context using state-space techniques, which enables the possibility of merging past temporal information with an update for each new acquisition. The dynamic system context allows us to exploit classical tools in this domain to perform the estimation of relevant variables. A general methodology is proposed, and a particular instance is defined in this study based on polarimetric radar data to track the phenological stages of a set of crops. A model generation from empirical data through principal component analysis is presented, and an extended Kalman filter is adapted to perform phenological stage estimation. Results employing quad-pol Radarsat-2 data over three different cereals are analyzed. The potential of this methodology to retrieve vegetation variables in real time is shown.This work was supported in part by the Spanish Ministry of Economy and Competitiveness (MINECO) and EU FEDER under Project TEC2011-28201-C02-02 and in part by the Generalitat Valenciana under Project ACOMP/2014/136

Assimilation of NDVI data in a land surface – Vegetation model for leaf area index predictions in a tree-grass ecosystem

Periodic observations of vegetation index, such as the normalized difference vegetation index (NDVI), can be used for data assimilation in heterogenous ecosystems. Indeed, the new Sentinel 2 Multispectral instrument and Landsat 8 Operational Land Imager sensor data are available at such high temporal and spatial resolutions that can be used to detect the patches of the main vegetation components (grass and trees) of heterogenous ecosystems, and capture their dynamics. We demonstrate the possibility to merge grass and tree NDVI observations and models, to optimally provide robust predictions of grass and tree leaf area index. The proposed assimilation approach assimilates NDVI data through the Ensemble Kalman filter (EnKF) and dynamically calibrates a key vegetation dynamic model parameter, the maintenance respiration coefficient (ma ). In the presence of large bias of the grass and tree ma base values, only the use of the proposed assimilation approach removes the large bias in the biomass balance, dynamically calibrating maintenance respiration coefficients, and corrects the model. The performance of a land surface - vegetation model was improved by assimilating observations of NDVI. The effective impact of the proposed assimilation approach on the evapotranspiration and CO2 uptake predictions in the heterogenous ecosystem is also demonstrated

Microwave Indices from Active and Passive Sensors for Remote Sensing Applications

Past research has comprehensively assessed the capabilities of satellite sensors operating at microwave frequencies, both active (SAR, scatterometers) and passive (radiometers), for the remote sensing of Earth’s surface. Besides brightness temperature and backscattering coefficient, microwave indices, defined as a combination of data collected at different frequencies and polarizations, revealed a good sensitivity to hydrological cycle parameters such as surface soil moisture, vegetation water content, and snow depth and its water equivalent. The differences between microwave backscattering and emission at more frequencies and polarizations have been well established in relation to these parameters, enabling operational retrieval algorithms based on microwave indices to be developed. This Special Issue aims at providing an overview of microwave signal capabilities in estimating the main land parameters of the hydrological cycle, e.g., soil moisture, vegetation water content, and snow water equivalent, on both local and global scales, with a particular focus on the applications of microwave indices

Désagrégation spatiale de températures Météosat par une méthode d'assimilation de données (lisseur particulaire) dans un modèle de surface continentale

Land surface temperature (LST) is one of the most important meteorologicalvariables giving access to water and energy budgets governing the Biosphere-Atmosphere continuum. To better monitor vegetation and energy states, we need hightemporal and spatial resolution measures of LST because its high variability in spaceand time.Despite the growing availability of Thermal Infra-Red (TIR) remote sensing LSTproducts, at different spatial and temporal resolutions, both high spatial resolution(HSR) and high temporal resolution (HTR) TIR data is still not possible because ofsatellite resolutions trade-off : the most frequent LST products being low spatial resolution(LSR) ones.It is therefore necessary to develop methods to estimate HSR/HTR LST from availableTIR LSR/HTR ones. This solution is known as "downscaling" and the presentthesis proposes a new approach for downscaling LST based on Data Assimilation (DA)methods. The basic idea is to constrain HSR/HTR LST dynamics, simulated by a dynamicalmodel, through the minimization of their respective aggregated LSTs discrepancytoward LSR observations, assuming that LST is homogeneous at the land cover typescale inside the LSR pixel.Our method uses a particle smoother DA method implemented in a land surfacemodel : SETHYS model (Suivie de l’Etat Hydrique de Sol). The proposed approach hasbeen firstly evaluated in a synthetic framework then validated using actual TIR LSTover a small area in South-East of France. Meteosat LST time series were downscaledfrom 5km to 90m and validated with ASTER HSR LST over one day. The encouragingresults conducted us to expand the study area and consider a larger assimilation periodof seven months. The downscaled Meteosat LSTs were quantitatively validated at1km of spatial resolution (SR) with MODIS data and qualitatively at 30m of SR withLandsat7 data. The results demonstrated good performances with downscaling errorsless than 2.5K at MODIS scale (1km of SR).La température des surfaces continentales (LST) est une variable météorologiquetrès importante car elle permet l’accès aux bilans d’énergie et d’eau ducontinuum Biosphère-Atmosphère. Sa haute variabilité spatio-temporelle nécessite desmesures à haute résolution spatiale (HRS) et temporelle (HRT) pour suivre au mieuxles états hydriques du sol et des végétations.La télédétection infrarouge thermique (IRT) permet d’estimer la LST à différentesrésolutions spatio-temporelles. Toutefois, les mesures les plus fréquentes sont souventà basse résolution spatiale (BRS). Il faut donc développer des méthodes pour estimerla LST à HRS à partir des mesures IRT à BRS/HRT. Cette solution est connue sous lenom de désagrégation et fait l’objet de cette thèse.Ainsi, une nouvelle approche de désagrégation basée sur l’assimilation de données(AD) est proposée. Il s’agit de contraindre la dynamique des LSTs HRS/HRT simuléespar un modèle en minimisant l’écart entre les LST agrégées et les données IRT àBRS/HRT, sous l’hypothèse d’homogénéité de la LST par type d’occupation des sols àl’échelle du pixel BRS. La méthode d’AD choisie est un lisseur particulaire qui a étéimplémenté dans le modèle de surface SETHYS (Suivi de l’Etat Hydrique du Sol).L’approche a été évaluée dans une première étape sur des données synthétiques etvalidée ensuite sur des données réelles de télédétection sur une petite région au Sud-Est de la France. Des séries de températures Météosat à 5 km de résolution spatialeont été désagrégées à 90m et validées sur une journée à l’aide de données ASTER.Les résultats encourageants nous ont conduit à élargir la région d’étude et la périoded’assimilation à sept mois. La désagrégation des produits Météosat a été validée quantitativementà 1km à l’aide de données MODIS et qualitativement à 30m à l’aide dedonnées Landsat7. Les résultats montrent de bonnes performances avec des erreursinférieures à 2.5K sur les températures désagrégées à 1km

Analyse de séries temporelles d’images à moyenne résolution spatiale : reconstruction de profils de LAI, démélangeage : application pour le suivi de la végétation sur des images MODIS

This PhD dissertation is concerned with time series analysis for medium spatial resolution (MSR) remote sensing images. The main advantage of MSR data is their high temporal rate which allows to monitor land use. However, two main problems arise with such data. First, because of cloud coverage and bad acquisition conditions, the resulting time series are often corrupted and not directly exploitable. Secondly, pixels in medium spatial resolution images are often “mixed” in the sense that the spectral response is a combination of the response of “pure” elements.These two problems are addressed in this PhD. First, we propose a data assimilation technique able to recover consistent time series of Leaf Area Index from corrupted MODIS sequences. To this end, a plant growth model, namely GreenLab, is used as a dynamical constraint. Second, we propose a new and efficient unmixing technique for time series. It is in particular based on the use of “elastic” kernels able to properly compare time series shifted in time or of various lengths.Experimental results are shown both on synthetic and real data and demonstrate the efficiency of the proposed methodologies.Cette thèse s’intéresse à l’analyse de séries temporelles d’images satellites à moyenne résolution spatiale. L’intérêt principal de telles données est leur haute répétitivité qui autorise des analyses de l’usage des sols. Cependant, deux problèmes principaux subsistent avec de telles données. En premier lieu, en raison de la couverture nuageuse, des mauvaises conditions d’acquisition, ..., ces données sont souvent très bruitées. Deuxièmement, les pixels associés à la moyenne résolution spatiale sont souvent “mixtes” dans la mesure où leur réponse spectrale est une combinaison de la réponse de plusieurs éléments “purs”. Ces deux problèmes sont abordés dans cette thèse. Premièrement, nous proposons une technique d’assimilation de données capable de recouvrer des séries temporelles cohérentes de LAI (Leaf Area Index) à partir de séquences d’images MODIS bruitées. Pour cela, le modèle de croissance de plantes GreenLab estutilisé. En second lieu, nous proposons une technique originale de démélangeage, qui s’appuie notamment sur des noyaux “élastiques” capables de gérer les spécificités des séries temporelles (séries de taille différentes, décalées dans le temps, ...)Les résultats expérimentaux, sur des données synthétiques et réelles, montrent de bonnes performances des méthodologies proposées

Sources of Atmospheric Fine Particles and Adsorbed Polycyclic Aromatic Hydrocarbons in Syracuse, New York

Land surface temperature (LST) images from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor have been widely utilized across scientific disciplines for a variety of purposes. The goal of this dissertation was to utilize MODIS LST for three spatial modeling applications within the conterminous United States (CONUS). These topics broadly encompassed agriculture and human health. The first manuscript compared the performance of all methods previously used to interpolate missing values in 8-day MODIS LST images. At low cloud cover (\u3c30%), the Spline spatial method outperformed all of the temporal and spatiotemporal methods by a wide margin, with median absolute errors (MAEs) ranging from 0.2°C-0.6°C. However, the Weiss spatiotemporal method generally performed best at greater cloud cover, with MAEs ranging from 0.3°C-1.2°C. Considering the distribution of cloud contamination and difficulty of implementing Weiss, using Spline under all conditions for simplicity would be sufficient. The second manuscript compared the corn yield predictive capability across the US Corn Belt of a novel killing degree day metric (LST KDD), computed with daily MODIS LST, and a traditional air temperature-based metric (Tair KDD). LST KDD was capable of predicting annual corn yield with considerably less error than Tair KDD (R2 /RMSE of 0.65/15.3 Bu/Acre vs. 0.56/17.2 Bu/Acre). The superior performance can be attributed to LST’s ability to better reflect evaporative cooling and water stress. Moreover, these findings suggest that long-term yield projections based on Tair and precipitation alone will contain error, especially for years of extreme drought. Finally, the third manuscript assessed the extent to which daily maximum heat index (HI) across the CONUS can be estimated by MODIS multispectral imagery in conjunction with land cover, topographic, and locational factors. The derived model was capable of estimating HI in 2012 with an acceptable level of error (R 2 = 0.83, RMSE = 4.4°F). LST and water vapor (WV) were, by far, the most important variables for estimation. Expanding this analytical framework to a more extensive study area (both temporally and spatially) would further validate these findings. Moreover, identifying an appropriate interpolation and downscaling approach for daily MODIS imagery would substantially increase the utility of the corn yield and HI models

Synergistic optical and microwave remote sensing approaches for soil moisture mapping at high resolution

Aplicat embargament des de la data de defensa fins al dia 1 d'octubre de 2022Soil moisture is an essential climate variable that plays a crucial role linking the Earth’s water, energy, and carbon cycles. It is responsible for the water exchange between the Earth’s surface and the atmosphere, and provides key information about soil evaporation, plant transpiration, and the allocation of precipitation into runoff, surface flow and infiltration. Therefore, an accurate estimation of soil moisture is needed to enhance our current climate and meteorological forecasting skills, and to improve our current understanding of the hydrological cycle and its extremes (e.g., droughts and floods). L-band Microwave passive and active sensors have been used during the last decades to estimate soil moisture, since there is a strong relationship between this variable and the soil dielectric properties. Currently, there are two operational L-band missions specifically devoted to globally measure soil moisture: the ESA’s Soil Moisture and the Ocean Salinity (SMOS), launched in November 2009; and the NASA’s Soil Moisture Active Passive (SMAP), launched in January 2015. The spatial resolution of the SMOS and SMAP radiometers, in the order of tens of kilometers (~40 km), is adequate for global applications. However, to fulfill the needs of a growing number of applications at local or regional scale, higher spatial detail (< 1 km) is required. To bridge this gap and improve the spatial resolution of the soil moisture maps, a variety of spatial enhancement or spatial (sub-pixel) disaggregation approaches have been proposed. This Ph.D. Thesis focuses on the study of the Earth’s surface soil moisture from remotely sensed observations. This work includes the implementation of several soil moisture retrieval techniques and the development, implementation, validation and comparison of different spatial enhancement or downscaling techniques, applied at local, regional, and continental scale. To meet these objectives, synergies between several active/passive microwave sensors (SMOS, SMAP and Sentinel-1) and optical/thermal sensors (MODIS) have been explored. The results are presented as follows: - Spatially consistent downscaling approach for SMOS using an adaptive moving window A passive microwave/optical downscaling algorithm for SMOS is proposed to obtain fine-scale soil moisture maps (1 km) from the native resolution (~40 km) of the instrument. This algorithm introduces the concept of a shape-adaptive window as a central improvement of the disaggregation technique presented by Piles et al. (2014), allowing its application at continental scales. - Assessment of multi-scale SMOS and SMAP soil moisture products across the Iberian Peninsula The temporal and spatial characteristics of SMOS and SMAP soil moisture products at coarse- and fine-scales are assessed in order to learn about their distinct features and the rationale behind them, tracing back to the physical assumptions they are based upon. - Impact of incidence angle diversity on soil moisture retrievals at coarse and fine scales An incidence angle (32.5°, 42.5° and 52.5°)-adaptive calibration of radiative transfer effective parameters single scattering albedo and soil roughness has been carried out, highlighting the importance of such parameterization to accurately estimate soil moisture at coarse-resolution. Then, these parameterizations are used to examine the potential application of a physically-based active-passive downscaling approach to upcoming microwave missions, namely CIMR, ROSE-L and Sentinel-1 Next Generation. Soil moisture maps obtained for the Iberian Peninsula at the three different angles, and at coarse and fine scales are inter-compared using in situ measurements and model data as benchmarks.La humedad del suelo es una variable climática esencial que juega un papel crucial en la relación de los ciclos del agua, la energía y el carbono de la Tierra. Es responsable del intercambio de agua entre la superficie de la Tierra y la atmósfera, y proporciona información crucial sobre la evaporación del suelo, la transpiración de las plantas y la distribución de la precipitación en escorrentía, flujo superficial e infiltración. Por lo tanto, es necesaria una estimación precisa de la humedad del suelo para mejorar las predicciones climáticas y meteorológicas, y comprender mejor el ciclo hidrológico y sus extremos (v.g., sequías e inundaciones). Los sensores pasivos y activos en banda L se han usado durante las últimas décadas para estimar la humedad del suelo debido a la relación directa que existe entre esta variable y las propiedades dieléctricas del suelo. Actualmente, hay dos misiones operativas en banda L específicamente dedicadas a medir la humedad del suelo a escala global: la misión Soil Moisture and Ocean Salinity (SMOS) de la ESA, lanzada en noviembre de 2009; y la misión Soil Moisture Active Passive (SMAP) de la NASA, lanzada en enero de 2015. La resolución espacial de los radiómetros SMOS y SMAP, del orden de unas decenas de kilómetros (~40 km), es adecuada para aplicaciones a escala global. Sin embargo, para satisfacer las necesidades de un número creciente de aplicaciones a escala local o regional, se requiere más detalle espacial (<1 km). Para solventar esta limitación y mejorar la resolución espacial de los mapas de humedad, se han propuesto diferentes técnicas de mejora o desagregación espacial. Esta Tesis se centra en el estudio de la humedad de la superficie terrestre a partir de datos obtenidos a través de teledetección. Este trabajo incluye la implementación de distintos algoritmos de recuperación de la humedad del suelo y el desarrollo, implementación, validación y comparación de distintas técnicas de desagregación, aplicadas a escala local, regional y continental. Para cumplir estos objetivos, se han explorado sinergias entre diferentes sensores de microondas activos/pasivos (SMOS, SMAP y Sentinel-1) y sensores ópticos/térmicos. Los resultados se presentan de la siguiente manera: - Técnica de desagregación espacialmente consistente, basada en una ventana móvil adaptativa, aplicada a los datos SMOS Se propone un algoritmo de desagregación del píxel basado en datos obtenidos de medidas radiométricas de microondas en banda L y datos ópticos, para mejorar la resolución espacial de los mapas de humedad del suelo desde la resolución nativa del instrumento (~40 km) hasta resoluciones de 1 km. El algoritmo introduce el concepto de una ventana de contorno adaptativo, como mejora principal sobre la técnica de desagregación presentada en Piles et al. (2014), permitiendo su implementación a escala continental. - Análisis multiescalar de productos de humedad del suelo SMAP y SMOS sobre la Península Ibérica Se han evaluado las características temporales y espaciales de distintos productos de humedad del suelo SMOS y SMAP, a baja y a alta resolución, para conocer sus características distintivas y comprender las razones de sus diferencias. Para ello, ha sido necesario rastrear los supuestos físicos en los que se basan. - Impacto del ángulo de incidencia en la recuperación de la humedad del suelo a baja y a alta resolución Se ha llevado a cabo una calibración adaptada al ángulo de incidencia (32.5°, 42.5° y 52.5°) de los parámetros efectivos, albedo de dispersión simple y rugosidad del suelo, descritos en el modelo de transferencia radiativa � − �, incidiendo en la importancia de esta parametrización para estimar la humedad del suelo de forma precisa a baja resolución. El resultado de las mismas se ha utilizado para estudiar la potencial aplicación de un algoritmo activo/pasivo de desagregación basado en la física para las próximas misiones de microondas, llamadas CIMR, ROSE-L y Sentinel-1 Next Generation. Los mapas de humedad recuperados a los tres ángulos de incidencia, tanto a baja como a alta resolución, se han obtenido para la Península Ibérica y se han comparado entre ellos usando como referencia mediciones de humedad in situ.Postprint (published version