Irrigated grassland monitoring using a time series of terraSAR-X and COSMO-skyMed X-Band SAR Data

[Departement_IRSTEA]Territoires [TR1_IRSTEA]SYNERGIE [Axe_IRSTEA]TETIS-ATTOSInternational audienceThe objective of this study was to analyze the sensitivity of radar signals in the X-band in irrigated grassland conditions. The backscattered radar signals were analyzed according to soil moisture and vegetation parameters using linear regression models. A time series of radar (TerraSAR-X and COSMO-SkyMed) and optical (SPOT and LANDSAT) images was acquired at a high temporal frequency in 2013 over a small agricultural region in southeastern France. Ground measurements were conducted simultaneously with the satellite data acquisitions during several grassland growing cycles to monitor the evolution of the soil and vegetation characteristics. The comparison between the Normalized Difference Vegetation Index (NDVI) computed from optical images and the in situ Leaf Area Index (LAI) showed a logarithmic relationship with a greater scattering for the dates corresponding to vegetation well developed before the harvest. The correlation between the NDVI and the vegetation parameters (LAI, vegetation height, biomass, and vegetation water content) was high at the beginning of the growth cycle. This correlation became insensitive at a certain threshold corresponding to high vegetation (LAI ~2.5 m2/m2). Results showed that the radar signal depends on variations in soil moisture, with a higher sensitivity to soil moisture for biomass lower than 1 kg/m². HH and HV polarizations had approximately similar sensitivities to soil moisture. The penetration depth of the radar wave in the X-band was high, even for dense and high vegetation; flooded areas were visible in the images with higher detection potential in HH polarization than in HV polarization, even for vegetation heights reaching 1 m. Lower sensitivity was observed at the X-band between the radar signal and the vegetation parameters with very limited potential of the X-band to monitor grassland growth. These results showed that it is possible to track gravity irrigation and soil moisture variations from SAR X-band images acquired at high spatial resolution (an incidence angle near 30°)

Coupling SAR C-band and optical data for soil moisture and leaf area index retrieval over irrigated grasslands

International audienceThe objective of this study was to develop an approach for estimating soil moisture and vegetation parameters in irrigated grasslands by coupling C-band polarimetric Synthetic Aperture Radar (SAR) and optical data. A huge dataset of satellite images acquired from RADARSAT-2 and LANDSAT-7/8, and in situ measurements were used to assess the relevance of several inversion configurations. A neural network (NN) inversion technique was used. The approach for this study was to use RADARSAT-2 and LANDSAT-7/8 images to investigate the potential for the combined use of new data from the new SAR sensor SENTINEL-1 and the new optical sensors LANDSAT-8 and SENTINEL-2. First, the impact of SAR polarization (mono-, dual- and full-polarizations configurations) and the Normalized Difference Vegetation Index (NDVI) calculated from optical data for the estimation error of soil moisture and vegetation parameters was studied. Next, the effect of some polarimetric parameters (Shannon entropy and Pauli components) on the inversion technique was also analyzed. Finally, configurations using in situ measurements of the fraction of absorbed photosynthetically active radiation (FAPAR) and the fraction of green vegetation cover (FCover) were also tested.The results showed that HH polarization is the SAR polarization most relevant to soil moisture estimates. An RMSE for soil moisture estimates of approximately 6 vol.% was obtained even for dense grassland cover. The use of in situ FAPAR and FCover only improved the estimate of the leaf area index (LAI) with an RMSE of approximately 0.37 m²/m². The use of polarimetric parameters did not improve the estimate of soil moisture and vegetation parameters. Good results were obtained for the biomass (BIO) and vegetation water content (VWC) estimates for BIO and VWC values lower than 2 and 1.5 kg/m², respectively (RMSE is of 0.38 kg/m² for BIO and 0.32 kg/m² for VWC). In addition, a high under-estimate was observed for BIO and VWC higher than 2 and 1.5 kg/m², respectively (a bias of -0.65 kg/m² on BIO estimates and -0.49 kg/m² on VWC estimates). Finally, the estimation of vegetation height (VEH) was carried out with an RMSE of 13.45 cm

Automatic methods for crop classification by merging satellite radar (sentinel 1) and optical (sentinel 2) . data and artificial intelligence analysis

Land use and land cover maps can support our understanding of coupled human- environment systems and provide important information for environmental modelling and water resource management. Satellite data are a valuable source for land use and land cover mapping. However, cloud-free or weather independent data are necessary to map cloud-prone regions. Merging radar with optical images would increase the accuracy of the study. Agricultural land cover is characterized by strong variations within relatively short time intervals. These dynamics are challenging for land cover classifications on the one hand, but deliver crucial information that can be used to improve the machine learning classifier’s performance on the other hand. A parcel-based map of the main crop classes of the Netherlands was produced implementing a script on GEE and using Copernicus data. The machine-learning model used is a Random Forest Classifier. This was done by combining time series of radar and multispectral images from Sentinel 1 and Sentinel 2 satellites, respectively. The results show the potential of providing useful information delivered by entirely open source data and uses a cloud computing-based approach. The algorithm combines the two satellites data of one year in a multibands image to feed in the classifier. Standard deviation and several vegetation indexes were added in order to have more variables for each 15-day-median image composite. The process paid particular attention to time variability of mean values of each field. This will provide useful information both for understanding differences among crops and variability over the phenology of the plant. The accuracy assessment demonstrates that several crop types (i.e. corn, tulip) can be better classified with both radar and optical images while others (i.e. sugar beet, barley) have an increased accuracy with only radar. The overall accuracy of RFC with optical and radar is 76% while it is 74% if only radar is used

Remote Sensing of Complex Permittivity and Penetration Depth of Soils Using P-Band SAR Polarimetry

A P-band SAR moisture estimation method is introduced for complex soil permittivity and penetration depth estimation using fully polarimetric P-band SAR signals. This method combines eigen- and model-based decomposition techniques for separation of the total backscattering signal into three scattering components (soil, dihedral, and volume). The incorporation of a soil scattering model allows for the first time the estimation of complex soil permittivity and permittivity-based penetration depth. The proposed method needs no prior assumptions on land cover characteristics and is applicable to a variety of vegetation types. The technique is demonstrated for airborne P-band SAR measurements acquired during the AirMOSS campaign (2012–2015). The estimated complex permittivity agrees well with climate and soil conditions at different monitoring sites. Based on frequency and permittivity, P-band penetration depths vary from 5 cm to 35 cm. This value range is in accordance with previous studies in the literature. Comparison of the results is challenging due to the sparsity of vertical soil in situ sampling. It was found that the disagreement between in situ measurements and SAR-based estimates originates from the discrepancy between the in situ measuring depth of the top-soil layer (0–5 cm) and the median penetration depth of the P-band waves (24.5–27 cm)

Assessment of multi - date sentinel - 1 polarizations and GLCM texture features capacity for onion and sunflower classification in an irrigated valley: An object level approach

The objective of this work is to evaluate the capacity of the C-band Synthetic Aperture Radar (SAR) time series imagery, acquired by the European satellite Sentinel-1 (S1), for the agriculture crop classification and its reliability to differentiate onion from sunflower, among others. The work then focused on classifying land cover in intensively cultivated agricultural regions. The study was developed in the Bonaerense Valley of the Colorado River (BVCR), Buenos Aires Province in Argentina, backed up by the field truth of 1634 field samples. In addition to the onion and sunflower crops, there are other crops present in the study area such as cereals, alfalfa, potatoes and maize, which are considered as the image background in the classification process. The field samples database was used for training and supporting a supervised classification with two machine learning algorithms—Random Forest (RF) and Support Vector Machine (SVM)—obtaining high levels of accuracy in each case. Different S1 SAR time-series features were used to assess the performance of S1 crop classification in terms of polarization VH+VV, Grey Level Co-occurrence Matrix (GLCM) image texture and a combination of both. The analysis of SAR data and their features was carried out at OBIA lot level (Object Based Image Analysis) showing an optimal strategy to counteract the effect of the residual and inherent speckle noise of the radar signal. In the process of differentiating onion and sunflower crops, the analysis of the VH+VV stack with the SVM algorithm delivered the best statistical classification results in terms of Overall Accuracy (OA) and Kappa Index, (Kp) when other crops (image background) were not considered (OA = 95.35%, Kp = 0.89). Certainly, the GLCM texture analysis derived from the S1 SAR images is a valuable source of information for obtaining very good classification results. When differentiating sunflower from onion considering also other crops present in the BVCR, the GLCM stack proved to be the most suitable dataset analyzed in this work (OA = 89.98%, Kp = 0.66 for SVM algorithm). This working methodology is applicable to other irrigated valleys in Argentina dedicated to intensive crops. There are also variables inherent to each lot, soil, crop and agricultural producer that differ according to the study area and that should be considered for each case in the future.Fil: Caballero, Gabriel. Universidad Blas Pascal. Centro de Investigación y Desarrollo Aplicado en Informática y Telecomunicaciones (CIADE-IT); ArgentinaFil: Platzech, Gabriel. INVAP. Government & Security Division; ArgentinaFil: Pezzola, Alejandro. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Casella, Alejandra. Instituto Nacional de Tecnología Agropecuaria (INTA). Instituto de Clima y Agua; ArgentinaFil: Winschel, Cristina Ines. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Hilario Ascasubi; ArgentinaFil: Silva, Samanta. Ministerio de Desarrollo Agrario (Buenos Aires, provincia). Colorado River Development Corporation (CORFO); ArgentinaFil: Ludueña, Emilia. INGTRADUCCIONES; ArgentinaFil: Pasqualotto, Nieves. Universidad de Valencia. Image Processing Laboratory (IPL); EspañaFil: Delegido, Jesús. Universidad de Valencia. Image Processing Laboratory (IPL); Españ

Satellite remote sensing priorities for better assimilation in crop growth models : winter wheat LAI and grassland mowing dates case studies

In a context of markets globalization, early, reliable and timely estimations of crop yields at regional to global scale are clearly needed for managing large agricultural lands, determining food pricing and trading policies but also for early warning of harvest shortfalls. Crop growth models are often used to estimate yields within the growing season. The uncertainties associated with these models contribute to the uncertainty of crop yield estimations and forecasts. Satellite remote sensing, through its ability to provide synoptic information on growth conditions over large geographic extents and in near real-time, is a key data source used to reduce uncertainties in biophysical models through data assimilation methods. This thesis aims at assessing possible improvements for the assimilation of remotely-sensed biophysical variables in crop growth models and to estimate their related errors reduction on modelled yield estimates. Assimilated observations are winter wheat leaf area index (LAI) and grassland mowing dates derived respectively from optical (MODIS) and microwave (ERS-2) data. These observations have been assimilated in WOFOST and LINGRA growth models. Observing System Simulation Experiments (OSSE) allowed assessing errors reduction on yield estimates after assimilation for different situations of accuracy and frequency of remotely-sensed estimates and for different assimilation strategies, indicating expected improvements with the currently available and forthcoming sensors; it supports also guidelines for future satellite missions. A main finding is the fact that yield estimate improvements are only possible for assimilation strategies able to correct the possible phenological discrepancies between the remotely-sensed and the simulated data. These phenological shifts arise mainly from uncertainties on the parameters and initial states driving the phenological stages in the models but are also due, in some situations, to lack of pixel purity because of the medium resolution of sensors such as MODIS. This thesis also identifies some of the main sources of uncertainties and assesses their impact on assimilation performances. The impact of water content and biomass on SAR backscattering of grasslands is specifically assessed. The backscattering of grasslands increases with the increases of water content and decreases with the biomass in dry conditions. Based on these results, methodologies to classify grasslands according to land use (mowing or grazing) and to detect mowings are designed and demonstrated. A good classification accuracy is observed (overall accuracy around 80%). Results for mowings detection are a bit lower as half of the mowings are correctly identified but the methodology can be considered as promising in particular in the perspective of very dense SAR time series as currently acquired operationally by Sentinel-1.(AGRO - Sciences agronomiques et ingénierie biologique) -- UCL, 201

Determination of Soil Moisture and Vegetation Parameters from Spaceborne C-Band SAR on Agricultural Areas

Soil moisture is an important factor influencing hydrological and meteorological exchange processes at the land surface. As ground measurements of soil moisture cannot provide spatial-ly distributed information, remote sensing of soil moisture using Synthetic Aperture Radar (SAR) offers an alternative. To derive soil moisture from vegetated areas with SAR, the influ-ence of vegetation parameters on SAR backscatter must be considered, though. The first part of the study analyses the potential to use a qualitative soil moisture index from ERS-SAR with high spatial resolution that can be used without ground truth soil moisture and vegetation data. The index ranges from low to high soil moisture instead of giving absolute soil moisture values. The method is applied to agricultural areas in the catchment of the river Rur in Germany. The soil moisture index represents wetting and drying tendencies well when compared to precipitation records and behaves like in-situ soil moisture regarding its variabil-ity. The analysis of spatial patterns from the soil moisture index by using semivariograms re-veals that differences in management that result for example in differences in evapotranspira-tion from one to the next agricultural field, are the only influence on spatial patterns of soil moisture in the Rur catchment. This study confirms the applicability of a high-resolution soil moisture index for monitoring soil moisture changes and to analyze spatial soil moisture pat-terns. The soil moisture index could be used as input to hydrological models and could substi-tute antecedent precipitation, which needs precipitation stations, as a proxy to soil moisture. The second part of the study examines the capability of dual-polarimetric C-Band SAR data with high incidence angles from the Sentinel-1 satellites to derive soil moisture and vegetation parameters quantitatively. A processing scheme for Sentinel-1 Level-1 data is presented to produce images of different SAR observables that are compared to extensive ground meas-urements of soil moisture and vegetation parameters. It shows that soil moisture retrieval is feasible from bare soil and maize with an RMSE of 7 Vol%. From other land use types, dif-ferent vegetation parameters could be retrieved with an error of around 25 % of their range, in median. Neither soil moisture nor vegetation parameters could be derived from grassland and triticale due to the influence of the thatch layer and the missing of a clear row structure. Both grassland and triticale are in contrast to the other crops not sown in rows on our research fields. The analysis has shown that the incidence angle is of main importance for the capability of C-band SAR to derive soil moisture and that the availability of at least one co- and cross-polarized channel is important for the quantitative retrieval of land surface parameters. The dual-pol H2α parameters were not meaningful for soil moisture and vegetation parameter re-trieval in this study

Modelling of Multi-Frequency Microwave Backscatter and Emission of Land Surface by a Community Land Active Passive Microwave Radiative Transfer Modelling Platform (CLAP)

Emission and backscattering signals of land surfaces at different frequencies have distinctive responses to soil and vegetation physical states. The use of multi-frequency combined active and passive microwave signals provides complementary information to better understand and interpret the observed signals in relation to surface states and the underlying physical processes. Such a capability also improves our ability to retrieve surface parameters and states such as soil moisture, freeze-thaw dynamics and vegetation biomass and vegetation water content (VWC) for ecosystem monitoring. We present here a prototype Community Land Active Passive Microwave Radiative Transfer Modelling platform (CLAP) for simulating both backscatter (&sigma;0) and emission (TB) signals of land surfaces, in which the CLAP is backboned by an air-to-soil transition model (ATS) (accounting for surface dielectric roughness) integrated with the Advanced Integral Equation Model (AIEM) for modelling soil surface scattering, and the Tor Vergata model for modelling vegetation scattering and the interaction between vegetation and soil parts. The CLAP was used to simulate both ground-based and space-borne multi-frequency microwave measurements collected at the Maqu observatory on the eastern Tibetan plateau. The ground-based systems include a scatterometer system (1&ndash;10 GHz) and an L-band microwave radiometer. The space-borne measurements are obtained from the X-band and C-band Advanced Microwave Scanning Radiometer 2 (AMSR2) radiation observations. The impacts of different vegetation properties (i.e., structure, water and temperature dynamics) and soil conditions (i.e., different moisture and temperature profiles) on the microwave signals were investigated by CLAP simulation for understanding factors that can account for diurnal variations of the observed signals. The results show that the dynamic VWC partially accounts for the diurnal variation of the observed signal at the low frequencies (i.e., S- and L-bands), while the diurnal variation of the observed signals at high frequencies (i.e., X- and C-bands) is more due to vegetation temperature changing, which implies the necessity to first disentangle the impact of vegetation temperature for the use of high frequency microwave signals. The model derived vegetation optical depth &tau; differs in terms of frequencies and different model parameterizations, while its diurnal variation depends on the diurnal variation of VWC regardless of frequency. After normalizing &tau; at multi-frequency by wavenumber, difference is still observed among different frequencies. This indicates that &tau; is indeed frequency-dependent, and &tau; for each frequency is suggested to be applied in the retrieval of soil and vegetation parameters. Moreover, &tau; at different frequencies (e.g., X-band and L-band) cannot be simply combined for constructing accurate long time series microwave-based vegetation product. To this purpose, it is suggested to investigate the role of the leaf water potential in regulating plant water use and its impact on the normalized &tau; at multi-frequency. Overall, the CLAP is expected to improve our capability for understanding and applying current and future multi-frequency space-borne microwave systems (e.g. those from ROSE-L and CIMR) for vegetation monitoring.</p

Mapping intra- and inter-annual dynamics in wetlands with multispectral, thermal and SAR time series

Kartierung der intra- und interannuellen Dynamik von Feuchtgebieten mit multispektralen, thermischen und SAR-Zeitreihen Die Analyse der aktuellen räumlichen Verbreitung und der zeitlichen Entwicklung von Feuchtgebieten stellt eine äußerst komplexe Aufgabe dar, welche durch die Saisonalität, die schwierige Zugänglichkeit und die besonderen Eigenschaften als Ökoton bedingt ist. Erdbeobachtungssysteme sind somit das am besten geeignete Werkzeug, um zeitliche und räumliche Muster von Feuchtgebieten auf globaler Ebene zu beobachten (saisonale Veränderungen und Langzeit-Trends) und um den Einfluss der menschlichen Aktivitäten auf ihre physischen und biologischen Eigenschaften zu untersuchen. Zur Kartierung von raum-zeitlichen Mustern wurden Zeitreihen von Radar- (Sentinel-1), Multispektral- (Sentinel-2) und Thermal-Satellitendaten (MODIS) in fünf Untersuchungsgebieten, mit für Feuchtgebiete unterschiedlichen typischen Charakteristika, untersucht. In Kapitel 1 werden die Problematik in Bezug auf die Definition von Feuchtgebieten erläutert und allgemeine Degradations-Trends beschrieben. Die Kapitel 2 und 3 behandeln einen Algorithmus, der Veränderungen mithilfe von SAR-Zeitreihen feststellt, sowie die Vorteile des Cloud-Computings für das operationelle Monitoring saisonaler Muster und die Erkennung kurzfristig auftretender Veränderungen. In den Kapiteln 4 und 5 werden die zwei Hauptursachen für den Verlust von Feuchtgebieten betrachtet: der Staudammbau und die Ausdehnung landwirtschaftlicher Flächen. In Kapitel 4 werden dichte Zeitreihen multispektraler (Sentinel-2) und SAR-Daten (Sentinel-1) verwendet, um die Feuchtgebiete Albaniens – eines Landes in dem konträre Pläne zum Ausbau seines Wasserkraftpotentials und dem Schutz intakter Flussökosysteme zu Spannungen führen – landesweit zu kartieren. Die synergetischen Vorteile, die sich durch die Fusionierung von multispektralen und SAR-Daten für die Klassifikation ergeben, werden dabei herausgestellt. Kapitel 5 veranschaulicht, dass die Kilombero-Überschwemmungsebene in Tansania ein großes und bedeutendes Feuchtgebiet ist, das in den vergangenen Jahren infolge der weitgehend unkontrollierten Ausbreitung landwirtschaftlicher Flächen in seiner Ausdehnung und seiner Ökologie stark beeinträchtigt wurde. Um die Auswirkungen der Landnutzungsänderungen des Feuchtgebietes während der vergangenen 18 Jahre zu analysieren, wurden eine Zeitreihe (2000 bis 2017) thermaler Daten (MODIS) analysiert. Die drei für die Zeitreihenanalyse angewandten Modelle zeigen, wie landwirtschaftliche Praktiken die Landoberflächentemperatur in den landwirtschaftlich genutzten Gebieten sowie in den angrenzenden natürlichen Feuchtgebieten erhöht haben.Due to wetlands’ seasonality, their difficult access and ecotone character, determining their actual extension and trends over time is a complex task. Earth Observation systems are the most appropriate tool to monitor their spatio-temporal patterns (seasonal changes and long term trends) at global scales, and to study the effects that human activities have in their physical and biological properties. In this work I use time series of radar (Sentinel-1), multispectral (Sentinel-2) and thermal (MODIS) imagery to map the spatio-temporal patterns in 5 wetlands of different characteristics. First, I introduce in chapter 1 the problematic of wetlands’ definitions and their degradation trends. I continue with a brief introduction on remote sensing, time series analysis, and their applications on wetlands’ research and management. In chapters 2 and 3 I implement an algorithm for change detection of time series of Sentinel-1 images and demonstrate the advantages of cloud computation for operational monitoring. In chapters 4 and 5 I address two of the main causes of wetland degradation: dam building and agricultural expansion. In chapter 4 I use dense time series of Sentinel-1 and Sentinel-2 images map all the wetlands of Albania; a country struggling between developing its large hydropower potential or preserving its intact and valuable river ecosystems. I evaluate the synergic advantages of fusing multispectral and radar imagery in combination with knowledge-based rules to produce classification of higher thematic and spatial resolutions. In chapter 5 I present how the Kilombero Floodplain, in Tanzania, has been degraded during the last years due to uncontrolled farmland expansion. I use a time series of thermal imagery (MODIS) from 2000 until 2017 to analyze the effect of land use changes on the wetland. I compare three models for time series analysis and reveal how farming practices have increased the surface temperature of the farmed area, as well as in adjacent natural wetlands.Mapeo de las dinámicas inter- e intra-anuales en humedales con series temporales de imágenes multiespectrales, termales y de radar Debido a la estacionalidad de los humedales, su difícil acceso y sus características de ecotono, determinar su actual extensión y sus tendencias a lo largo del tiempo es una tarea compleja. Los sistemas de observación terrestres son la herramienta más apropiada para monitorear sus patrones espacio-temporales (estacionalidad y tendencias a largo plazo) a escalas globales, y para estudiar los efectos que las actividades humanas causan en sus propiedades físicas y biológicas. En esta tesis uso series temporales de imágenes radar (Sentinel-1), multiespectrales (Sentinel-2) y termales (MODIS) para mapear los patrones espacio-temporales de 5 humedales de diferentes características. En el capítulo 1 describo los retos que derivan de las diferentes definiciones que existen de los humedales. También presento las tendencias globales de degradación que la mayoría de los humedales continúan experimentando en los últimos años. Continúo con una breve introducción de los sistemas de teledetección remota, análisis de series temporales, y sus aplicaciones a la investigación y gestión de los humedales. En los capítulos 2 y 3 implemento un algoritmo de detección de cambios para series temporales de imágenes radar, y muestro las ventajas de usar sistemas de computación en la nube para monitorear cambios en la cobertura del suelo a corto plazo. En los capítulos 4 y 5 trato con dos de las causas más comunes de degradación de humedales: la construcción de presas y la expansión de la agricultura. En el capítulo 4 uso series temporales de imágenes multiespectrales (Sentinel-2) y radar (Sentinel-1) para mapear todos los humedales Albania; un país que se debate entre desarrollar su potencial hidroenergético o preservar sus valiosos e intactos ecosistemas de rivera. Mediante la fusión de imágenes radar y multiespectrales y el uso de reglas de decisión genero un mapa de suficiente resolución espacial y temática para que pueda ser usado por sectores interesados y gestores. En el capítulo 5 presento como las llanuras inundables de Kilombero, en Tanzania, han sido degradadas durante los últimos años debido a la expansión incontrolada de la agricultura. Usando series temporales de imágenes termales (MODIS) desde 2000 hasta 2017 y mapas de cambios de usos del suelo, determino los efectos que estos cambios han tenido en el humedal. Comparo 3 modelos diferentes de análisis de series temporales y muestro cómo la expansión de la agricultura ha incrementado la temperatura superficial terrestre, no solo de la zona cultivada, sino también de zonas adyacentes aún naturales

Monitoring the Sustainable Intensification of Arable Agriculture:the Potential Role of Earth Observation

Sustainable intensification (SI) has been proposed as a possible solution to the conflicting problems of meeting projected increases in food demand and preserving environmental quality. SI would provide necessary production increases while simultaneously reducing or eliminating environmental degradation, without taking land from competing demands. An important component of achieving these aims is the development of suitable methods for assessing the temporal variability of both the intensification and sustainability of agriculture. Current assessments rely on traditional data collection methods that produce data of limited spatial and temporal resolution. Earth Observation (EO) provides a readily accessible, long-term dataset with global coverage at various spatial and temporal resolutions. In this paper we demonstrate how EO could significantly contribute to SI assessments, providing opportunities to quantify agricultural intensity and environmental sustainability. We review an extensive body of research on EO-based methods to assess multiple indicators of both agricultural intensity and environmental sustainability. To date these techniques have not been combined to assess SI; here we identify the opportunities and initial steps required to achieve this. In this context, we propose the development of a set of essential sustainable intensification variables (ESIVs) that could be derived from EO data