Evaluation of the Uncertainty in Satellite-Based Crop State Variable Retrievals Due to Site and Growth Stage Specific Factors and Their Potential in Coupling with Crop Growth Models

Coupling crop growth models and remote sensing provides the potential to improve our understanding of the genotype x environment x management (G X E X M) variability of crop growth on a global scale. Unfortunately, the uncertainty in the relationship between the satellite measurements and the crop state variables across different sites and growth stages makes it diffcult to perform the coupling. In this study, we evaluate the effects of this uncertainty with MODIS data at the Mead, Nebraska Ameriflux sites (US-Ne1, US-Ne2, and US-Ne3) and accurate, collocated Hybrid-Maize (HM) simulations of leaf area index (LAI) and canopy light use effciency (LUECanopy). The simulations are used to both explore the sensitivity of the satellite-estimated genotype X management (G X M) parameters to the satellite retrieval regression coeffcients and to quantify the amount of uncertainty attributable to site and growth stage specific factors. Additional ground-truth datasets of LAI and LUECanopy are used to validate the analysis. The results show that uncertainty in the LAI/satellite measurement regression coeffcients lead to large uncertainty in the G X Mparameters retrievable from satellites. In addition to traditional leave-one-site-out regression analysis, the regression coeffcient uncertainty is assessed by evaluating the retrieval performance of the temporal change in LAI and LUECanopy. The weekly change in LAI is shown to be retrievable with a correlation coeffcient absolute value (|r|) of 0.70 and root-mean square error (RMSE) value of 0.4, which is significantly better than the performance expected if the uncertainty was caused by random error rather than secondary effects caused by site and growth stage specific factors (an expected |r| value of 0.36 and RMSE value of 1.46 assuming random error). As a result, this study highlights the importance of accounting for site and growth stage specific factors in remote sensing retrievals for future work developing methods coupling remote sensing with crop growth models

Improving Retrievals of Crop Vegetation Parameters from Remote Sensing Data

Agricultural systems are difficult to model because crop growth is driven by the strongly nonlinear interaction of Genotype x Environment x Management (G x E x M) factors. Due to the nonlinearity in the interaction of these factors, the amount of data necessary to develop and utilize models to accurately predict the performance of agricultural systems at an operational scale is large. Satellite remote sensing provides the potential to vastly increase the amount of data available for modelling agricultural systems as a result of its high revisit time and spatial coverage. Unfortunately, there have been significant difficulties in deploying remote sensing for many agricultural modelling applications because of the uncertainty involved in the retrievals. In this dissertation, we show that collecting farmer-provided agro-managment information has the potential to reduce the uncertainty in the retrieval products obtained from remote sensing observations. Specifically, both field-scale and regional-scale analysis are used to show that secondary factor variability is a very significant cause of uncertainty in both crop growth modelling and agricultural remote sensing that needs to be addressed through increased data collection. In order to address this need for increased data availability, a method is developed that allows geolocated crop growth model simulations to be used to train satellite-based crop state variable retrievals, which is then validated at regional scale. The method developed provides a general robust methodology to create a large-scale platform that would allow farmers to share data with government agencies and universities to improve crop state variable retrievals and crop growth modelling and provide farmers, government, industry, and researchers with insights and predictive capability into crop growth at both field and regional scales

Utilizing Collocated Crop Growth Model Simulations to Train Agronomic Satellite Retrieval Algorithms

Due to its worldwide coverage and high revisit time, satellite-based remote sensing provides the ability to monitor in-season crop state variables and yields globally. In this study, we presented a novel approach to training agronomic satellite retrieval algorithms by utilizing collocated crop growth model simulations and solar-reflective satellite measurements. Specifically, we showed that bidirectional long short-term memory networks (BLSTMs) can be trained to predict the in-season state variables and yields of Agricultural Production Systems sIMulator (APSIM) maize crop growth model simulations from collocated Moderate Resolution Imaging Spectroradiometer (MODIS) 500-m satellite measurements over the United States Corn Belt at a regional scale. We evaluated the performance of the BLSTMs through both k-fold cross validation and comparison to regional scale ground-truth yields and phenology. Using k-fold cross validation, we showed that three distinct in-season maize state variables (leaf area index, aboveground biomass, and specific leaf area) can be retrieved with cross-validated R2 values ranging from 0.4 to 0.8 for significant portions of the season. Several other plant, soil, and phenological in-season state variables were also evaluated in the study for their retrievability via k-fold cross validation. In addition, by comparing to survey-based United State Department of Agriculture (USDA) ground truth data, we showed that the BLSTMs are able to predict actual county-level yields with R2 values between 0.45 and 0.6 and actual state-level phenological dates (emergence, silking, and maturity) with R2 values between 0.75 and 0.85. We believe that a potential application of this methodology is to develop satellite products to monitor in-season field-scale crop growth on a global scale by reproducing the methodology with field-scale crop growth model simulations (utilizing farmer recorded field-scale agromanagement data) and collocated high-resolution satellite data (fused with moderate-resolution satellite data)

Unlocking the benefits of spaceborne imaging spectroscopy for sustainable agriculture

With the Environmental Mapping and Analysis Program (EnMAP) mission, launched on April 1st 2022, new opportunities unfold for precision farming and agricultural monitoring. The recurring acquisition of spectrometric imagery from space, contiguously resolving the electromagnetic spectrum in the optical domain (400—2500 nm) within close narrow bands, provides unprecedented data about the interaction of radiation with biophysical and biochemical crop constituents. These interactions manifest in spectral reflectance, carrying important information about crop status and health. This information may be incorporated in agricultural management systems to support necessary efforts to maximize yields against the backdrop of an increased food demand by a growing world population. At the same time, it enables the effective optimization of fertilization and pest control to minimize environmental impacts of agriculture. Deriving biophysical and biochemical crop traits from hyperspectral reflectance thereby always relies on a model. These models are categorized into (1) parametric, (2) nonparametric, (3) physically-based, and (4) hybrid retrieval schemes. Parametric methods define an explicit parameterized expression, relating a number of spectral bands or derivates thereof with a crop trait of interest. Nonparametric methods comprise linear techniques, such as principal component analysis (PCA) which addresses collinearity issues between adjacent bands and enables compression of full spectral information into dimensionality reduced, maximal informative principal components (PCs). Nonparametric nonlinear methods, i.e., machine learning (ML) algorithms apply nonlinear transformations to imaging spectroscopy data and are therefore capable of capturing nonlinear relationships within the contained spectral features. Physically-based methods represent an umbrella term for radiative transfer models (RTMs) and related retrieval schemes, such as look-up-table (LUT) inversion. A simple, easily invertible and specific RTM is the Beer-Lambert law which may be used to directly infer plant water content. The most widely used general and invertible RTM is the one-dimensional canopy RTM PROSAIL, which is coupling the Leaf Optical Properties Spectra model PROSPECT and the canopy reflectance model 4SAIL: Scattering by Arbitrarily Inclined Leaves. Hybrid methods make use of synthetic data sets created by RTMs to calibrate parametric methods or to train nonparametric ML algorithms. Due to the ill-posed nature of RTM inversion, potentially unrealistic and redundant samples in a LUT need to be removed by either implementing physiological constraints or by applying active learning (AL) heuristics. This cumulative thesis presents three different hybrid approaches, demonstrated within three scientific research papers, to derive agricultural relevant crop traits from spectrometric imagery. In paper I the Beer-Lambert law is applied to directly infer the thickness of the optically active water layer (i.e., EWT) from the liquid water absorption feature at 970 nm. The model is calibrated with 50,000 PROSPECT spectra and validated over in situ data. Due to separate water content measurements of leaves, stalks, and fruits during the Munich-North-Isar (MNI) campaigns, findings indicate that depending on the crop type and its structure, different parts of the canopy are observed with optical sensors. For winter wheat, correlation between measured and modelled water content was most promising for ears and leaves, reaching coefficients of determination (R2) up to 0.72 and relative RMSE (rRMSE) of 26%, and in the case of corn for the leaf fraction only (R2 = 0.86, rRMSE = 23%). These results led to the general recommendation to collect destructive area-based plant organ specific EWT measurements instead of the common practice to upscale leaf-based EWT measurements to canopy water content (CWC) by multiplication of the leaf area index (LAI). The developed and calibrated plant water retrieval (PWR) model proved to be transferable in space and time and is ready to be applied to upcoming EnMAP data and any other hyperspectral imagery. In paper II the parametric concept of spectral integral ratios (SIR) is introduced to retrieve leaf chlorophyll a and b content (Cab), leaf carotenoid content (Ccx) and leaf water content (Cw) simultaneously from imaging spectroscopy data in the wavelength range 460—1100 nm. The SIR concept is based on automatic separation of respective absorption features through local peak and intercept analysis between log-transformed reflectance and convex hulls. The approach was validated over a physiologically constrained PROSAIL simulated database, considering natural Ccx-Cab relations and green peak locations. Validation on airborne spectrometric HyMAP data achieved satisfactory results for Cab (R2 = 0.84; RMSE = 9.06 µg cm-2) and CWC (R2 = 0.70; RMSE = 0.05 cm). Retrieved Ccx values were reasonable according to Cab-Ccx-dependence plausibility analysis. Mapping of the SIR results as multiband images (3-segment SIR) allows for an intuitive visualization of dominant absorptions with respect to the three considered biochemical variables. Hence, the presented SIR algorithm allows for computationally efficient and RTM supported robust retrievals of the two most important vegetation pigments as well as of water content and is applicable on satellite imaging spectroscopy data. In paper III a hybrid workflow is presented, combining RTM with ML for inferring crop carbon content (Carea) and aboveground dry and fresh biomass (AGBdry, AGBfresh). The concept involves the establishment of a PROSAIL training database, dimensionality reduction using PCA, optimization in the sampling domain using AL against the 4-year MNI campaign dataset, and training of Gaussian process regression (GPR) ML algorithms. Internal validation of the GPR-Carea and GPR-AGB models achieved R2 of 0.80 for Carea, and R2 of 0.80 and 0.71 for AGBdry and AGBfresh, respectively. Validation with an independent dataset, comprising airborne AVIRIS NG imagery (spectrally resampled to EnMAP) and in situ measurements, successfully demonstrated mapping capabilities for both bare and green fields and generated reliable estimates over winter wheat fields at low associated model uncertainties (< 40%). Overall, the proposed carbon and biomass models demonstrate a promising path toward the inference of these crucial variables over cultivated areas from upcoming spaceborne hyperspectral acquisitions, such as from EnMAP. As conclusions, the following important findings arise regarding parametric and nonparametric hybrid methods as well as in view of the importance of in situ data collection. (1) Uncertainties within the RTM PROSAIL should always be considered. A possible reduction of these uncertainties is thereby opposed to the invertibility of the model and its intended simplicity. (2) Both physiological constraints and AL heuristics should be applied to reduce unrealistic parameter combinations in a PROSAIL calibration or training database. (3) State-of-the-art hybrid ML approaches with the ability to provide uncertainty intervals are anticipated as most promising approach for solving inference problems from hyperspectral Earth observation data due to their synergistic use of RTMs and the high flexibility, accuracy and consistency of nonlinear nonparametric methods. (4) Parametric hybrid approaches, due to their algorithmic transparency, enable deeper insights into fundamental physical limitations of optical remote sensing as compared to ML approaches. (5) Integration-based indices that make full use of available hyperspectral information may serve as physics-aware dimensionality reduced input for ML algorithms to either improve estimations or to serve as endmember for crop type discrimination when additional time series information is available. (6) The validation of quantitative model-based estimations is crucial to evaluate and improve their performance in terms of the underlying assumptions, model parameterizations, and input data. (7) In the face of soon-to-be-available EnMAP data, collection of in situ data for validation of retrieval methods should aim at high variability of measured crop types, high temporal variability over the whole growing season, as well as include area- and biomass-based destructive measurements instead of LAI-upscaled leaf measurements. Provided the perfect functionality of the payload instruments, the success of the EnMAP mission and the here presented methods depend critically on a low-noise, accurate atmospherically corrected reflectance product. High-level outputs of the retrieval methods presented in this thesis may be incorporated into agricultural decision support systems for fertilization and irrigation planning, yield estimation, or estimation of the soil carbon sequestration potential to enable a sustainable intensive agriculture in the future.Mit der am 1. April 2022 gestarteten Satellitenmission Environmental Mapping and Analysis Program (EnMAP) eröffnen sich neue Möglichkeiten für die Präzisionslandwirtschaft und das landwirtschaftliche Monitoring. Die wiederkehrende Erfassung spektrometrischer Bilder aus dem Weltraum, welche das elektromagnetische Spektrum im optischen Bereich (400—2500 nm) innerhalb von engen, schmalen Bändern zusammenhängend auflösen, liefert nie dagewesene Daten über die Interaktionen von Strahlung und biophysikalischen und biochemischen Pflanzenbestandteilen. Diese Wechselwirkungen manifestieren sich in der spektralen Reflektanz, die wichtige Informationen über den Zustand und die Gesundheit der Pflanzen enthält. Vor dem Hintergrund einer steigenden Nachfrage nach Nahrungsmitteln durch eine wachsende Weltbevölkerung können diese Informationen in landwirtschaftliche Managementsysteme einfließen, um eine notwendige Ertragsmaximierung zu unterstützen. Gleichzeitig können sie eine effiziente Optimierung der Düngung und Schädlingsbekämpfung ermöglichen, um die Umweltauswirkungen der Landwirtschaft zu minimieren. Die Ableitung biophysikalischer und biochemischer Pflanzeneigenschaften aus hyperspektralen Reflektanzdaten ist dabei immer von einem Modell abhängig. Diese Modelle werden in (1) parametrische, (2) nichtparametrische, (3) physikalisch basierte und (4) hybride Ableitungsmethoden kategorisiert. Parametrische Methoden definieren einen expliziten parametrisierten Ausdruck, der eine Reihe von Spektralkanälen oder deren Ableitungen mit einem Pflanzenmerkmal von Interesse in Beziehung setzt. Nichtparametrische Methoden umfassen lineare Techniken wie die Hauptkomponentenanalyse (PCA). Diese adressieren Kollinearitätsprobleme zwischen benachbarten Kanälen und komprimieren die gesamte Spektralinformation in dimensionsreduzierte, maximal informative Hauptkomponenten (PCs). Nichtparametrische nichtlineare Methoden, d. h. Algorithmen des maschinellen Lernens (ML), wenden nichtlineare Transformationen auf bildgebende Spektroskopiedaten an und sind daher in der Lage, nichtlineare Beziehungen innerhalb der enthaltenen spektralen Merkmale zu erfassen. Physikalisch basierte Methoden sind ein Oberbegriff für Strahlungstransfermodelle (RTM) und damit verbundene Ableitungsschemata, d. h. Invertierungsverfahren wie z. B. die Invertierung mittels Look-up-Table (LUT). Ein einfaches, leicht invertierbares und spezifisches RTM stellt das Lambert-Beer'sche Gesetz dar, das zur direkten Ableitung des Wassergehalts von Pflanzen verwendet werden kann. Das am weitesten verbreitete, allgemeine und invertierbare RTM ist das eindimensionale Bestandsmodell PROSAIL, eine Kopplung des Blattmodells Leaf Optical Properties Spectra (PROSPECT) mit dem Bestandsreflexionsmodell 4SAIL (Scattering by Arbitrarily Inclined Leaves). Bei hybriden Methoden werden von RTMs generierte, synthetische Datenbanken entweder zur Kalibrierung parametrischer Methoden oder zum Training nichtparametrischer ML-Algorithmen verwendet. Aufgrund der Äquifinalitätsproblematik bei der RTM-Invertierung, müssen potenziell unrealistische und redundante Simulationen in einer solchen Datenbank durch die Implementierung natürlicher physiologischer Beschränkungen oder durch die Anwendung von Active Learning (AL) Heuristiken entfernt werden. In dieser kumulativen Dissertation werden drei verschiedene hybride Ansätze zur Ableitung landwirtschaftlich relevanter Pflanzenmerkmale aus spektrometrischen Bilddaten vorgestellt, die anhand von drei wissenschaftlichen Publikationen demonstriert werden. In Paper I wird das Lambert-Beer'sche Gesetz angewandt, um die Dicke der optisch aktiven Wasserschicht (bzw. EWT) direkt aus dem Absorptionsmerkmal von flüssigem Wasser bei 970 nm abzuleiten. Das Modell wird mit 50.000 PROSPECT-Spektren kalibriert und anhand von In-situ-Daten validiert. Aufgrund separater Messungen des Wassergehalts von Blättern, Stängeln und Früchten während der München-Nord-Isar (MNI)-Kampagnen, zeigen die Ergebnisse, dass je nach Kulturart und -struktur, unterschiedliche Teile des Bestandes mit optischen Sensoren beobachtet werden können. Bei Winterweizen wurde die höchste Korrelation zwischen gemessenem und modelliertem Wassergehalt für Ähren und Blätter erzielt und sie erreichte Bestimmtheitsmaße (R2) von bis zu 0,72 bei einem relativen RMSE (rRMSE) von 26%, bei Mais entsprechend nur für die Blattfraktion (R2 = 0,86, rRMSE = 23%). Diese Ergebnisse führten zu der allgemeinen Empfehlung, Kompartiment-spezifische EWT-Bestandsmessungen zu erheben, anstatt der üblichen Praxis, blattbasierte EWT-Messungen durch Multiplikation mit dem Blattflächenindex (LAI) auf den Bestandswassergehalt (CWC) hochzurechnen. Das entwickelte und kalibrierte Modell zur Ableitung des Pflanzenwassergehalts (PWR) erwies sich als räumlich und zeitlich übertragbar und kann auf bald verfügbare EnMAP-Daten und andere hyperspektrale Bilddaten angewendet werden. In Paper II wird das parametrische Konzept der spektralen Integralratios (SIR) eingeführt, um den Chlorophyll a- und b-Gehalt (Cab), den Karotinoidgehalt (Ccx) und den Wassergehalt (Cw) simultan aus bildgebenden Spektroskopiedaten im Wellenlängenbereich 460-1100 nm zu ermitteln. Das SIR-Konzept basiert auf der automatischen Separierung der jeweiligen Absorptionsmerkmale durch lokale Maxima- und Schnittpunkt-Analyse zwischen log-transformierter Reflektanz und konvexen Hüllen. Der Ansatz wurde anhand einer physiologisch eingeschränkten PROSAIL-Datenbank unter Berücksichtigung natürlicher Ccx-Cab-Beziehungen und Positionen der Maxima im grünen Wellenlängenbereich validiert. Die Validierung mit flugzeuggestützten spektrometrischen HyMAP-Daten ergab zufriedenstellende Ergebnisse für Cab (R2 = 0,84; RMSE = 9,06 µg cm-2) und CWC (R2 = 0,70; RMSE = 0,05 cm). Die ermittelten Ccx-Werte wurden anhand einer Plausibilitätsanalyse entsprechend der Cab-Ccx-Abhängigkeit als sinnvoll bewertet. Die Darstellung der SIR-Ergebnisse als mehrkanalige Bilder (3 segment SIR) ermöglicht zudem eine auf die drei betrachteten biochemischen Variablen bezogene, intuitive Visualisierung der dominanten Absorptionen. Der vorgestellte SIR-Algorithmus ermöglicht somit wenig rechenintensive und RTM-gestützte robuste Ableitungen der beiden wichtigsten Pigmente sowie des Wassergehalts und kann in auf jegliche zukünftig verfügbare Hyperspektraldaten angewendet werden. In Paper III wird ein hybrider Ansatz vorgestellt, der RTM mit ML kombiniert, um den Kohlenstoffgehalt (Carea) sowie die oberirdische trockene und frische Biomasse (AGBdry, AGBfresh) abzuschätzen. Das Konzept umfasst die Erstellung einer PROSAIL-Trainingsdatenbank, die Dimensionsreduzierung mittels PCA, die Reduzierung der Stichprobenanzahl mittels AL anhand des vier Jahre umspannenden MNI-Kampagnendatensatzes und das Training von Gaussian Process Regression (GPR) ML-Algorithmen. Die interne Validierung der GPR-Carea und GPR-AGB-Modelle ergab einen R2 von 0,80 für Carea und einen R2 von 0,80 bzw. 0,71 für AGBdry und AGBfresh. Die Validierung auf einem unabhängigen Datensatz, der flugzeuggestützte AVIRIS-NG-Bilder (spektral auf EnMAP umgerechnet) und In-situ-Messungen umfasste, zeigte erfolgreich die Kartierungsfähigkeiten sowohl für offene Böden als auch für grüne Felder und führte zu zuverlässigen Schätzungen auf Winterweizenfeldern bei geringen Modellunsicherheiten (< 40%). Insgesamt zeigen die vorgeschlagenen Kohlenstoff- und Biomassemodelle einen vielversprechenden Ansatz auf, der zur Ableitung dieser wichtigen Variablen über Anbauflächen aus künftigen weltraumgestützten Hyperspektralaufnahmen wie jenen von EnMAP genutzt werden kann. Als Schlussfolgerungen ergeben sich die folgenden wichtigen Erkenntnisse in Bezug auf parametrische und nichtparametrische Hybridmethoden sowie bezogen auf die Bedeutung der In-situ-Datenerfassung. (1) Unsicherheiten innerhalb des RTM PROSAIL sollten immer berücksichtigt werden. Eine mögliche Verringerung dieser Unsicherheiten steht dabei der Invertierbarkeit des Modells und dessen beabsichtigter Einfachheit entgegen. (2) Sowohl physiologische Einschränkungen als auch AL-Heuristiken sollten angewendet werden, um unrealistische Parameterkombinationen in einer PROSAIL-Kalibrierungs- oder Trainingsdatenbank zu reduzieren. (3) Modernste ML-Ansätze mit der Fähigkeit, Unsicherheitsintervalle bereitzustellen, werden als vielversprechendster Ansatz für die Lösung von Inferenzproblemen aus hyperspektralen Erdbeobachtungsdaten aufgrund ihrer synergetischen Nutzung von RTMs und der hohen Flexibilität, Genauigkeit und Konsistenz nichtlinearer nichtparametrischer Methoden angesehen. (4) Parametrische hybride Ansätze ermöglichen aufgrund ihrer algorithmischen Transparenz im Vergleich zu ML-Ansätzen tiefere Einblicke in die grundlegenden physikalischen Grenzen der optischen Fernerkundung. (5) Integralbasierte Indizes, die die verfügbare hyperspektrale Information voll ausschöpfen, können als physikalisch-basierte dimensionsreduzierte Inputs für ML-Algorithmen dienen, um entweder Schätzungen zu verbessern oder um als Eingangsdaten die verbesserte Unterscheidung von Kulturpflanzen zu ermöglichen, sobald zusätzliche Zeitreiheninformationen verfügbar sind. (6) Die Validierung quantitativer modellbasierter Schätzungen ist von entscheidender Bedeutung für die Bewertung und Verbesserung ihrer Leistungsfähigkeit in Bezug auf die zugrunde liegenden Annahmen, Modellparametrisierungen und Eingabedaten. (7) Angesichts der bald verfügbaren EnMAP-Daten sollte die Erhebung von In-situ-Daten zur Validierung von Ableitungsmethoden auf eine hohe Variabilität der gemessenen Pflanzentypen und eine hohe zeitliche Variabilität über die gesamte Vegetationsperiode abzielen sowie flächen- und biomassebasierte destruktive Messungen anstelle von LAI-skalierten Blattmessungen umfassen. Unter der Voraussetzung, dass die Messinstrumente perfekt funktionieren, hängt der Erfolg der EnMAP-Mission und der hier vorgestellten Methoden entscheidend von einem rauscharmen, präzise atmosphärisch korrigierten Reflektanzprodukt ab. Die Ergebnisse der in dieser Arbeit vorgestellten Methoden können in landwirtschaftliche Entscheidungsunterstützungssysteme für die Dünge- oder Bewässerungsplanung, die Ertragsabschätzung oder die Schätzung des Potenzials der Kohlenstoffbindung im Boden integriert werden, um eine nachhaltige Intensivlandwirtschaft in der Zukunft zu ermöglichen

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)

Earth Observations and Integrative Models in Support of Food and Water Security

Global food production depends upon many factors that Earth observing satellites routinely measure about water, energy, weather, and ecosystems. Increasingly sophisticated, publicly-available satellite data products can improve efficiencies in resource management and provide earlier indication of environmental disruption. Satellite remote sensing provides a consistent, long-term record that can be used effectively to detect large-scale features over time, such as a developing drought. Accuracy and capabilities have increased along with the range of Earth observations and derived products that can support food security decisions with actionable information. This paper highlights major capabilities facilitated by satellite observations and physical models that have been developed and validated using remotely-sensed observations. Although we primarily focus on variables relevant to agriculture, we also include a brief description of the growing use of Earth observations in support of aquaculture and fisheries

Impacts Of Market Forces And Agricultural Practices On Land Surfaces Related To Scientific Modeling And Remote Sensing Applications

Due to the rotational needs of crops, homogenous crop fields, and external influences such as market and policy changes, crop production generates significant changes to the landscape on annual and semi-annual basis. In this study we looked at two aspects of this change. In the first aspect of the study, we attempt to account for market and policy driven producer’s decision making through a new model constructed by pairing an economics model with the ALMANAC crop simulation model via a two-way coupling. This coupled model approach integrated farmer’s land-use choices based on relative economic returns and produced dynamic land use probabilities for ALMANAC simulations through a feedback loop. The coupled model approach was inter-compared with static crop modeling through a historic acreage approach, and comparable accuracies were found from both modeling efforts for the 2014 growing season. Furthermore, as a proof-concept effort, the method was applied to evaluate the impact of two scenarios on crop simulations: major crops (maize, soybean, and wheat) intensification through price increases (e.g. market change), as well as incentivized grassland conservation (e.g. policy change). The results of this sensitivity study suggest that the coupled system has the capability of integrating economic factors into traditional crop simulation, allowing for insight into the impacts of changes in markets and policies on agricultural landscapes and crop yields. In the second aspect of this study, changes to surface albedo driven by these landscape changes are investigated. Using collocated Moderate Resolution Imaging Spectroradiometer (MODIS) derived Bidirectional Reflectance Distribution Function (BRDF) with the Cropland Data Layer (CDL), we computed the daily albedo of homogenous agricultural fields across the United States for 55 crop types by wavelength, sky-type, day of year, crop, and hardiness zone over a four-year period (2015-2018). This study suggests that cropland spectral albedo is complicated by large variations over the course of the growing season, which can result in changes in reflectivity up to a factor of 2 at most wavelengths. This change was found to be unique per crop type, but predictable year-to-year for individual crops within specific regions, so generating a lookup table that incorporates these factors for use in remote sensing and atmospheric modeling applications is viable for albedo estimation. Additionally, impacts of crop types on broadband albedo were studied and found to be less conspicuous than the individual wavelength counterpart, but still significant over cropland. The results were used to evaluate the accuracy of a common method of albedo estimation, where NDVI is used as a proxy for albedo over cropland, and the NDVI method was found to have some significant limitations dependent on wavelength and day of year. Finally, a database of surface albedo variations as a function of growing period is constructed for 55 crops common to croplands across the United States. The constructed database can be used to aid both satellite remote sensing applications and long-term weather modeling efforts by providing a method for parameter adjustments based on crop driven albedo changes, including changes in cropland composition related to commodity markets and other external factors

Enhancing regional estimates of evapotranspiration with earth observation data

Food security and food sustainability are high on the global policy agenda. Reliable information on crop water use and terrestrial water stress are important to ensure an optimal use of available water resources and for enhancing crop production. Remote sensing provides a feasible avenue to estimate regional evapotranspiration (ET), which can be employed to assess terrestrial water stress. However, the heterogeneity of land surfaces and the accumulated errors from various inputs often result in substantial biases in most global or regional ET models across different landscapes. Reducing uncertainties in available ET products or remote sensing (RS)-based models and obtaining regional ET estimates with improved accuracy is important for effectively using ET to support agricultural monitoring and water resources managements. This thesis first compared different Priestly-Taylor (PT)-based methods that use three Earth observation-based alternatives - apparent thermal inertia (ATI), microwave soil moisture (SM), and optical spectral indices based on shortwave infrared (SWIR) to assess soil evaporation over cropland and grassland regions. Using FLUXNET data as ET reference, the results illustrated that the incorporation of the SWIR-based soil moisture divergence index (SMDI) and microwave-based SM into monthly soil evaporation led to 6% and 5% increase in explained ET variances and reduced RMSE by 23.2% and 13.1% for cropland and grassland, respectively, as compared to PT-JPL using atmospheric reanalysis data only. The results suggested that a combination of optical SWIR and microwave SM has good potential to improve the PT-JPL model accuracy for agricultural landscapes. Based on the performance of different PT-based methods, ET estimates derived from the revised PT method were used to assess water budgets across 53 catchments in central-western Europe with a humid climate and were compared with three additional ET data sources (MOD16, GLEAM, and PT-JPL). Surprisingly, all RS-based ET estimates significantly diverged from water balance-based ET (ETWB) in 45 humid catchments, whereas most previous studies that focussed on arid catchments or on the global scale found significantly less divergence. Using ET retrievals from the Budyko framework and upscaled ET from FLUXCOM as references, the closure errors of water budgets were sensitive to errors arising from precipitation data in humid regions and the water balance approach was found to overestimate ET during heavy rainfall events. Instead, the Budyko framework that describes the partitioning of precipitation to ET as a functional balance between atmospheric water supply (precipitation, P) and demand (potential evapotranspiration, PET) had good correlation with ET ensemble from multiple data sources and upscaled ET from FLUXCOM product. 161 Summary The results indicated that errors from precipitation and terrestrial water storage anomalies introduce large uncertainties in ETWB, thereby complicating water balance validation in humid regions across multiple timesteps. To improve the application of ETWB for benchmarking ETEB in humid regions, high-quality input data should be used or – like the Budyko framework – energy constraints should be considered. The thesis then proceeds to explore causes for the notable deviations between observed and Budyko-predicted water balances in certain catchments. The results revealed that for humid catchments, topography and seasonal cumulative moisture surplus can explain the spatial distributions of Budyko scatter with r higher than 0.65, whereas soil properties and vegetation indices explained little of the variance (r≤0.30). Temporally, the interannual variability of Budyko scatter was negatively correlated with annual average vegetation indices, particularly for catchments with relatively low vegetation cover. This thesis provides valuable insights to the interpretation of the Budyko framework and offers possible solutions to improve its performance to predict the spatiotemporal variability of water balances. Lastly, to address the deviations from the predictive Budyko curve, additional controls of hydrological partitioning were introduced to correct Budyko scatter between catchments and between years. The results illustrated that the use of catchment climatic seasonality properties and topography attributes is effective in reproducing the Budyko parameter (w) with an r of 0.76 and RMSE of 0.49 for all 45 catchments in central-western Europe. After the correction of temporal Budyko scatter using interannual variability of vegetation information and the fraction of precipitation falling as snow, the performance of the modified Budyko-type equation improves with respect to the original Budyko framework, in comparison to ETWB at catchment scale (∆r of 0.26 and ∆RMSE of 19.19 mm/yr). When compared with the gridded ET ensemble using energy balance, the enhanced Budyko framework is generally effective to reproduce the spatial distribution of ET with good similarity, even in ungauged regions. Overall, the revised Budyko framework shows improved performance in predicting water balances and can be applied to assess crop water use and terrestrial water stress at regional scale, particularly in ungauged areas. Overall, this thesis contributes significantly to the enhancement of regional ET estimation using Earth observation. It proposes a novel blended parameterization for soil moisture constraints in the modified PT-JPL model, which is capable of capturing the soil evaporation more accurately within agroecosystems. Meanwhile, this thesis proposes a new water balance-based validation method that uses the Budyko framework integrated with environmental parameters. By developing improved RS-based models and water balance-based validation methods, this thesis provides valuable insights into the complexities of ET 162 Summary estimation at the regional scale. These findings are expected to advance the application of ET in decision-making regarding the management of agriculture and water resources

Surface Soil Moisture Retrievals from Remote Sensing:Current Status, Products &amp; Future Trends

Advances in Earth Observation (EO) technology, particularly over the last two decades, have shown that soil moisture content (SMC) can be measured to some degree or other by all regions of the electromagnetic spectrum, and a variety of techniques have been proposed to facilitate this purpose. In this review we provide a synthesis of the efforts made during the last 20 years or so towards the estimation of surface SMC exploiting EO imagery, with a particular emphasis on retrievals from microwave sensors. Rather than replicating previous overview works, we provide a comprehensive and critical exploration of all the major approaches employed for retrieving SMC in a range of different global ecosystems. In this framework, we consider the newest techniques developed within optical and thermal infrared remote sensing, active and passive microwave domains, as well as assimilation or synergistic approaches. Future trends and prospects of EO for the accurate determination of SMC from space are subject to key challenges, some of which are identified and discussed within. It is evident from this review that there is potential for more accurate estimation of SMC exploiting EO technology, particularly so, by exploring the use of synergistic approaches between a variety of EO instruments. Given the importance of SMC in Earth’s land surface interactions and to a large range of applications, one can appreciate that its accurate estimation is critical in addressing key scientific and practical challenges in today’s world such as food security, sustainable planning and management of water resources. The launch of new, more sophisticated satellites strengthens the development of innovative research approaches and scientific inventions that will result in a range of pioneering and ground-breaking advancements in the retrievals of soil moisture from space

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