Improving ecological forecasts using model and data constraints

Terrestrial ecosystems are essential to human well-being, but their future remains highly uncertain, as evidenced by the huge disparities in model projections of the land carbon sink. The existence of these disparities despite the recent explosion of novel data streams, including the TRY plant traits database, the Landsat archive, and global eddy covariance tower networks, suggests that these data streams are not being utilized to their full potential by the terrestrial ecosystem modeling community. Therefore, the overarching objective of my dissertation is to identify how these various data streams can be used to improve the precision of model predictions by constraining model parameters. In chapter 1, I use a hierarchical multivariate meta-analysis of the TRY database to assess the dependence of trait correlations on ecological scale and evaluate the utility of these correlations for constraining ecosystem model parameters. I find that global trait correlations are generally consistent within plant functional types, and leveraging the multivariate trait space is an effective way to constrain trait estimates for data-limited traits and plant functional types. My next two chapters assess the ability to measure traits using remote sensing by exploring the links between leaf traits and reflectance spectra. In chapter 2, I introduce a method for estimating traits from spectra via radiative transfer model inversion. I then use this approach to show that although the precise location, width, and quantity of spectral bands significantly affects trait retrieval accuracy, a wide range of sensor configurations are capable of providing trait information. In chapter 3, I apply this approach to a large database of leaf spectra to show that traits vary as much within as across species, and much more across species within a functional type than across functional types. Finally, in chapter 4, I synthesize the findings of the previous chapters to calibrate a vegetation model's representation of canopy radiative transfer against observed remotely-sensed surface reflectance. Although the calibration successfully constrained canopy structural parameters, I identify issues with model representations of wood and soil reflectance that inhibit its ability to accurately reproduce remote sensing observations

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

Multiscale remote sensing of plant physiology and carbon uptake

This study investigated the use of optical remote sensing for estimating leaf and canopy scale light use efficiency (LUE) and carbon exchange. In addition, a new leaf level model capable of predicting dynamic changes in apparent reflectance due to chlorophyll fluorescence was developed. A leaf level study was conducted to assess the applicability of passive remote sensing as a tool to measure the reduction, and the subsequent recovery, of photosynthetic efficiency during the weeks following transplantation. Spectral data were collected on newly planted saplings for a period of 8 weeks, as well as gas exchange measurements of LUE and PAM fluorescence measurements. A set of spectral indices, including the Photochemical Reflectance Index (PRI), were calculated from the reflectance measurements. A marked depression in photosynthetic rate occurred in the weeks after outplanting followed by a gradual increase, with recovery occurring in the later stages of the experimental period. As with photosynthetic rate, there was a marked trend in PRI values over the study period but no trend was observed in chlorophyll based indices. The study demonstrated that hyperspectral remote sensing has the potential to be a useful tool in the detection and monitoring of the dynamic effects of transplant shock. Relationships between hyperspectral reflectance indices, airborne carbon exchange measurements and satellite observations of ground cover were then explored across a heterogeneous Arctic landscape. Measurements were collected during August 2008, using the University of Edinburgh’s research aircraft, from an Arctic forest tundra zone in northern Finland as part of the Arctic Biosphere Atmosphere Coupling at Multiple Scales (ABACUS) study. Surface fluxes of CO2 were calculated using the eddy covariance method from airborne data that were collected from the same platform as hyperspectral reflectance measurements. Airborne CO2 fluxes were compared to MODIS vegetation indices. In addition, LUE was estimated from airborne flux data and compared to airborne measurements of PRI. There were no significant relationships between MODIS vegetation indices and airborne flux observations. There were weak to moderate (R2 = 0.4 in both cases) correlations between PRI and LUE and between PRI and incident radiation. A new coupled physiological radiative transfer model that predicts changes in the apparent reflectance of a leaf, due to chlorophyll fluorescence, was developed. The model relates a physically observable quantity, chlorophyll fluorescence, to the sub leaf level processes that cause the emission. An understanding of the dynamics of the processes that control fluorescence emission on multiple timescales should aid in the interpretation of this complex signal. A Markov Chain Monte Carlo (MCMC) algorithm was used to optimise biochemical model parameters by fitting model simulations of transient chlorophyll fluorescence to measured reflectance spectra. The model was then validated against an independent data set. The model was developed as a precursor to a full canopy scheme. To scale to the canopy and to use the model on trans-seasonal time scales, the effects of temperature and photoinhibition on the model biochemistry needs to be taken into account, and a full canopy radiative transfer scheme, such as FluorMOD, must be developed

Hyperspectral Analysis of Leaf Pigments and Nutritional Elements in Tallgrass Prairie Vegetation

Understanding the spatial distribution of forage quality is important to address critical research questions in grassland science. Due to its efficiency and accuracy, there has been a widespread interest in mapping the canopy vegetation characteristics using remote sensing methods. In this study, foliar chlorophylls, carotenoids, and nutritional elements across multiple tallgrass prairie functional groups were quantified at the leaf level using hyperspectral analysis in the region of 470–800 nm, which was expected to be a precursor to further remote sensing of canopy vegetation quality. A method of spectral standardization was developed using a form of the normalized difference, which proved feasible to reduce the interference from background effects in the leaf reflectance measurements. Chlorophylls and carotenoids were retrieved through inverting the physical model PROSPECT 5. The foliar nutritional elements were modeled empirically. Partial least squares regression was used to build the linkages between the high-dimensional spectral predictor variables and the foliar biochemical contents. Results showed that the retrieval of leaf biochemistry through hyperspectral analysis can be accurate and robust across different tallgrass prairie functional groups. In addition, correlations were found between the leaf pigments and nutritional elements. Results provided insight into the use of pigment-related vegetation indices as the proxy of plant nutrition quality

Understanding the temporal dimension of the red-edge spectral region for forest decline detection using high-resolution hyperspectral and Sentinel-2a imagery

The operational monitoring of forest decline requires the development of remote sensing methods that are sensitive to the spatiotemporal variations of pigment degradation and canopy defoliation. In this context, the red-edge spectral region (RESR) was proposed in the past due to its combined sensitivity to chlorophyll content and leaf area variation. In this study, the temporal dimension of the RESR was evaluated as a function of forest decline using a radiative transfer method with the PROSPECT and 3D FLIGHT models. These models were used to generate synthetic pine stands simulating decline and recovery processes over time and explore the temporal rate of change of the red-edge chlorophyll index (CI) as compared to the trajectories obtained for the structure-related Normalized Difference Vegetation Index (NDVI). The temporal trend method proposed here consisted of using synthetic spectra to calculate the theoretical boundaries of the subspace for healthy and declining pine trees in the temporal domain, defined by CItime=n/CItime=n+1 vs. NDVItime=n/NDVItime=n+1. Within these boundaries, trees undergoing decline and recovery processes showed different trajectories through this subspace. The method was then validated using three high-resolution airborne hyperspectral images acquired at 40 cm resolution and 260 spectral bands of 6.5 nm full-width half-maximum (FWHM) over a forest with widespread tree decline, along with field-based monitoring of chlorosis and defoliation (i.e., ‘decline’ status) in 663 trees between the years 2015 and 2016. The temporal rate of change of chlorophyll vs. structural indices, based on reflectance spectra extracted from the hyperspectral images, was different for trees undergoing decline, and aligned towards the decline baseline established using the radiative transfer models. By contrast, healthy trees over time aligned towards the theoretically obtained healthy baseline. The applicability of this temporal trend method to the red-edge bands of the MultiSpectral Imager (MSI) instrument on board Sentinel-2a for operational forest status monitoring was also explored by comparing the temporal rate of change of the Sentinel-2-derived CI over areas with declining and healthy trees. Results demonstrated that the Sentinel-2a red-edge region was sensitive to the temporal dimension of forest condition, as the relationships obtained for pixels in healthy condition deviated from those of pixels undergoing decline.JRC.D.1-Bio-econom

Remote sensing tools for monitoring grassland plant leaf traits and biodiversity

Rocchini, Duccio1This project has received funding from the European Union’s Horizon 2020 Research and Innovation Program under the Marie Skłodowska-Curie Grant No. 721995 (project Trustee).openGrasslands are one of the most important ecosystems on Earth, covering approximately onethird of the Earth’s surface. Grassland biodiversity is important as many services provided by such ecosystems are crucial for the human economy and well-being. Given the importance of grasslands ecosystems, in recent years research has been carried out on the potential to monitor them with novel remote sensing techniques. Improved detectors technology and novel sensors providing finescale hyperspectral imagery have been enabling new methods to monitor plant traits (PTs) and biodiversity. The aims of the work were to study different approaches to monitor key grassland PTs such as Leaf Area Index (LAI) and biodiversity-related traits. The thesis consists of 3 parts: 1) Evaluating the performance of remote sensing methods to estimate LAI in grassland ecosystems, 2) Estimating plant biodiversity by using the optical diversity approach in grassland ecosystems, and 3) Investigating the relationship between PTs variability with alpha and beta diversity for the applicability of the optical diversity approach in a subalpine grassland of the Italian Alps To evaluate the performance of remote sensing methods to estimate LAI, temporal and spatial observations of hyperspectral reflectance and LAI were analyzed at a grassland site in Monte Bondone, Italy (IT-MBo). In 2018, ground temporal observations of hyperspectral reflectance and LAI were carried out at a grassland site in Neustift, Austria (AT-NEU). To estimate biodiversity, in 2018 and 2019 a floristics survey was conducted to determine species composition and hyperspectral data were acquired at two grassland sites: IT-MBo and University of Padova’s Experimental Farm, Legnaro, Padua, Italy (IT-PD) respectively. Furthermore, in 2018, biochemistry analysis of the biomass samples collected from the grassland site IT-MBo was carried out to determine the foliar biochemical PTs variability. The results of the thesis demonstrated that the grassland spectral response across different spectral regions (Visible: VIS, red-edge: RE, Near-infrared: NIR) showed to be both site-specific and scale-dependent. In the first part of the thesis, the performance of spectral vegetation indices (SVIs) based on visible, red-edge (RE), and NIR bands alongside SVIs solely based or NIRshoulder bands (wavelengths 750 - 900 nm) was evaluated. A strong correlation (R2 > 0.8) was observed between grassland LAI and both RE and NIR-shoulder SVIs on a temporal basis, but not on a spatial basis. Using the PROSAIL Radiative Transfer Model (RTM), it was demonstrated that grassland structural heterogeneity strongly affects the ability to retrieve LAI, with high uncertainties due to structural and biochemical PTs co-variation. In the second part, the applicability of the spectral variability hypothesis (SVH) was questioned and highlighted the challenges to use high-resolution hyperspectral images to estimate biodiversity in complex grassland ecosystems. It was reported that the relationship between biodiversity (Shannon, Richness, Simpson, and Evenness) and optical diversity metrics (Coefficient of variation (CV) and Standard deviation (SD)) is not consistent across plant communities. The results of the second part suggested that biodiversity in terms of species richness could be estimated by optical diversity metrics with an R2 = 0.4 at the IT-PD site where the grassland plots were artificially established and are showing a lower structure and complexity from the natural grassland plant communities. On the other hand, in the natural ecosystems at IT-MBo, it was more difficult to estimate biodiversity indices, probably due to structural and biochemical PTs co-variation. The 18 effects of canopy non-vegetative elements (flowers and dead material), shadow pixels, and overexposed pixels on the relationship between optical diversity metrics and biodiversity indices were highlighted. In the third part, we examined the relationship between PTs variability (at both local and community scales, measured by standard deviation and by the Euclidean distances of the biochemical and biophysical PTs respectively) and taxonomic diversity (both α-diversity and βdiversity, measured by Shannon’s index and by Jaccard dissimilarity index of the species, families, and functional groups percent cover respectively) in Monte Bondone, Trentino province, Italy. The results of the study showed that the PTs variability metrics at alpha scale were not correlated with α-diversity. However, the results at the community scale (β-diversity) showed that some of the investigated biochemical and biophysical PTs variations metrics were associated with β-diversity. The SVH approach was also tested to estimate β-diversity and we found that spectral diversity calculated by spectral angular mapper (SAM) showed to be a better proxy of biodiversity in the same ecosystem where the spectral diversity failed to estimate alpha diversity, this leading to the conclusion that the link between functional and species diversity may be an indicator of the applicability of optical sampling methods to estimate biodiversity. The findings of the thesis highlighted that grassland structural heterogeneity strongly affects the ability to retrieve both LAI and biodiversity, with high uncertainties due to structural and biochemical PTs co-variation at complex grassland ecosystems. In this context, the uncertainties of satellite-based products (e.g., LAI) in monitoring grassland canopies characterized by either spatially or temporally varying structure need to be carefully taken into account. The results of the study highlighted that the poor performance of optical diversity proxies in estimating biodiversity in structurally heterogeneous grasslands might be due to the complex relationships between functional diversity and biodiversity, rather than the impossibility to detect functional diversity with spectral proxiesopenImran, H.A

Quantitative estimation of vegetation traits and temporal dynamics using 3-D radiative transfer models, high-resolution hyperspectral images and satellite imagery

Large-scale monitoring of vegetation dynamics by remote sensing is key to detecting early signs of vegetation decline. Spectral-based indicators of phys-iological plant traits (PTs) have the potential to quantify variations in pho-tosynthetic pigments, chlorophyll fluorescence emission, and structural changes of vegetation as a function of stress. However, the specific response of PTs to disease-induced decline in heterogeneous canopies remains largely unknown, which is critical for the early detection of irreversible damage at different scales. Four specific objectives were defined in this research: i) to assess the feasibility of modelling the incidence and severity of Phytophthora cinnamomi and Xylella fastidiosa based on PTs and biophysical properties of vegetation; ii) to assess non-visual early indicators, iii) to retrieve PT using radiative transfer models (RTM), high-resolution imagery and satellite observations; and iv) to establish the basis for scaling up PTs at different spatial resolutions using RTM for their retrieval in different vegetation co-vers. This thesis integrates different approaches combining field data, air- and space-borne imagery, and physical and empirical models that allow the retrieval of indicators and the evaluation of each component’s contribution to understanding temporal variations of disease-induced symptoms in heter-ogeneous canopies. Furthermore, the effects associated with the understory are introduced, showing not only their impact but also providing a compre-hensive model to account for it. Consequently, a new methodology has been established to detect vegetation health processes and the influence of biotic and abiotic factors, considering different components of the canopy and their impact on the aggregated signal. It is expected that, using the presented methods, existing remote sensors and future developments, the ability to detect and assess vegetation health globally will have a substantial impact not only on socio-economic factors, but also on the preservation of our eco-system as a whole