Fusing simulated GEDI, ICESat-2 and NISAR data for regional aboveground biomass mapping

Accurate mapping of forest aboveground biomass (AGB) is critical for better understanding the role of forests in the global carbon cycle. NASA's current GEDI and ICESat-2 missions as well as the upcoming NISAR mission will collect synergistic data with different coverage and sensitivity to AGB. In this study, we present a multi-sensor data fusion approach leveraging the strength of each mission to produce wall-to-wall AGB maps that are more accurate and spatially comprehensive than what is achievable with any one sensor alone. Specifically, we calibrate a regional L-band radar AGB model using the sparse, simulated spaceborne lidar AGB estimates. We assess our data fusion framework using simulations of GEDI, ICESat-2 and NISAR data from airborne laser scanning (ALS) and UAVSAR data acquired over the temperate high AGB forest and complex terrain in Sonoma County, California, USA. For ICESat-2 and GEDI missions, we simulate two years of data coverage and AGB at footprint level are estimated using realistic AGB models. We compare the performance of our fusion framework when different combinations of the sparse simulated GEDI and ICEsat-2 AGB estimates are used to calibrate our regional L-band AGB models. In addition, we test our framework at Sonoma using (a) 1-ha square grid cells and (b) similarly sized irregularly shaped objects. We demonstrate that the estimated mean AGB across Sonoma is more accurately estimated using our fusion framework than using GEDI or ICESat-2 mission data alone, either with a regular grid or with irregular segments as mapping units. This research highlights methodological opportunities for fusing new and upcoming active remote sensing data streams toward improved AGB mapping through data fusion.</p

Forest Height Inversion by Combining Single-Baseline TanDEM-X InSAR Data with External DTM Data

Forest canopy height estimation is essential for forest management and biomass estimation. In this study, we aimed to evaluate the capacity of TanDEM-X interferometric synthetic aperture radar (InSAR) data to estimate canopy height with the assistance of an external digital terrain model (DTM). A ground-to-volume ratio estimation model was proposed so that the canopy height could be precisely estimated from the random-volume-over-ground (RVoG) model. We also refined the RVoG inversion process with the relationship between the estimated penetration depth (PD) and the phase center height (PCH). The proposed method was tested by TanDEM-X InSAR data acquired over relatively homogenous coniferous forests (Teruel test site) and coniferous as well as broadleaved forests (La Rioja test site) in Spain. Comparing the TanDEM-X-derived height with the LiDAR-derived height at plots of size 50 m × 50 m, the root-mean-square error (RMSE) was 1.71 m (R2 = 0.88) in coniferous forests of Teruel and 1.97 m (R2 = 0.90) in La Rioja. To demonstrate the advantage of the proposed method, existing methods based on ignoring ground scattering contribution, fixing extinction, and assisting with simulated spaceborne LiDAR data were compared. The impacts of penetration and terrain slope on the RVoG inversion were also evaluated. The results show that when a DTM is available, the proposed method has the optimal performance on forest height estimation.This work was supported in part by the National Natural Science Foundation of China under Grant 41820104005, Grant 42030112, and Grant 41904004, Hunan Natural Science Foundation under Grant 2021JJ30808, and in part by the Spanish Ministry of Science and Innovation, Agencia Estatal de Investigacion, under Projects PID2020-117303GB-C22/AEI/10.13039/501100011033 and PROWARM (PID2020-118444GA-I00/AEI/10.13039/501100011033)

The Global Ecosystem Dynamics Investigation: High-resolution laser ranging of the Earth’s forests and topography

Obtaining accurate and widespread measurements of the vertical structure of the Earths forests has been a longsought goal for the ecological community. Such observations are critical for accurately assessing the existing biomass of forests, and how changes in this biomass caused by human activities or variations in climate may impact atmospheric CO2 concentrations. Additionally, the three-dimensional structure of forests is a key component of habitat quality and biodiversity at local to regional scales. The Global Ecosystem Dynamics Investigation (GEDI) was launched to the International Space Station in late 2018 to provide high-quality measurements of forest vertical structure in temperate and tropical forests between 51.6 N & S latitude. The GEDI instrument is a geodetic-class laser altimeter/waveform lidar comprised of 3 lasers that produce 8 transects of structural information. Over its two-year nominal lifetime GEDI is anticipated to provide over 10 billion waveforms at a footprint resolution of 25 m. These data will be used to derive a variety of footprint and gridded products, including canopy height, canopy foliar profiles, Leaf Area Index (LAI), sub-canopy topography and biomass. Additionally, data from GEDI are used to demonstrate the efficacy of its measurements for prognostic ecosystem modeling, habit and biodiversity studies, and for fusion using radar and other remote sensing instruments. GEDI science and technology are unique: no other space-based mission has been created that is specifically optimized for retrieving vegetation vertical structure. As such, GEDI promises to advance our understanding of the importance of canopy vertical variations within an ecological paradigm based on structure, composition and function

The European Space Agency BIOMASS mission: Measuring forest above-ground biomass from space

The primary objective of the European Space Agency's 7th Earth Explorer mission, BIOMASS, is to determine the worldwide distribution of forest above-ground biomass (AGB) in order to reduce the major uncertainties in calculations of carbon stocks and fluxes associated with the terrestrial biosphere, including carbon fluxes associated with Land Use Change, forest degradation and forest regrowth. To meet this objective it will carry, for the first time in space, a fully polarimetric P-band synthetic aperture radar (SAR). Three main products will be provided: global maps of both AGB and forest height, with a spatial resolution of 200 m, and maps of severe forest disturbance at 50 m resolution (where “global” is to be understood as subject to Space Object tracking radar restrictions). After launch in 2022, there will be a 3-month commissioning phase, followed by a 14-month phase during which there will be global coverage by SAR tomography. In the succeeding interferometric phase, global polarimetric interferometry Pol-InSAR coverage will be achieved every 7 months up to the end of the 5-year mission. Both Pol-InSAR and TomoSAR will be used to eliminate scattering from the ground (both direct and double bounce backscatter) in forests. In dense tropical forests AGB can then be estimated from the remaining volume scattering using non-linear inversion of a backscattering model. Airborne campaigns in the tropics also indicate that AGB is highly correlated with the backscatter from around 30 m above the ground, as measured by tomography. In contrast, double bounce scattering appears to carry important information about the AGB of boreal forests, so ground cancellation may not be appropriate and the best approach for such forests remains to be finalized. Several methods to exploit these new data in carbon cycle calculations have already been demonstrated. In addition, major mutual gains will be made by combining BIOMASS data with data from other missions that will measure forest biomass, structure, height and change, including the NASA Global Ecosystem Dynamics Investigation lidar deployed on the International Space Station after its launch in December 2018, and the NASA-ISRO NISAR L- and S-band SAR, due for launch in 2022. More generally, space-based measurements of biomass are a core component of a carbon cycle observation and modelling strategy developed by the Group on Earth Observations. Secondary objectives of the mission include imaging of sub-surface geological structures in arid environments, generation of a true Digital Terrain Model without biases caused by forest cover, and measurement of glacier and icesheet velocities. In addition, the operations needed for ionospheric correction of the data will allow very sensitive estimates of ionospheric Total Electron Content and its changes along the dawn-dusk orbit of the mission

GEDI launches a new era of biomass inference from space

Accurate estimation of aboveground forest biomass stocks is required to assess the impacts of land use changes such as deforestation and subsequent regrowth on concentrations of atmospheric CO2. The Global Ecosystem Dynamics Investigation (GEDI) is a lidar mission launched by NASA to the International Space Station in 2018. GEDI was specifically designed to retrieve vegetation structure within a novel, theoretical sampling design that explicitly quantifies biomass and its uncertainty across a variety of spatial scales. In this paper we provide the estimates of pan-tropical and temperate biomass derived from two years of GEDI observations. We present estimates of mean biomass densities at 1 km resolution, as well as estimates aggregated to the national level for every country GEDI observes, and at the sub-national level for the United States. For all estimates we provide the standard error of the mean biomass. These data serve as a baseline for current biomass stocks and their future changes, and the mission's integrated use of formal statistical inference points the way towards the possibility of a new generation of powerful monitoring tools from space

The global tree carrying capacity (keynote)

editorial reviewe

FUSING GEDI LIDAR AND TANDEM-X INSAR OBSERVATIONS FOR IMPROVED FOREST STRUCTURE AND BIOMASS MAPPING

The upcoming NASA’s Global Ecosystem Dynamics Investigation (GEDI) mission presents an unprecedented opportunity to advance current global biomass estimates. However, gaps are expected between GEDI’s ground tracks, requiring the development of fusion-based methodologies to contiguously map forest biomass at satisfactory resolutions and accuracies. This dissertation is built on the complementary advantages of observations from GEDI and DLR’s TerraSAR-X/TanDEM-X (TDX)) Interferometric Synthetic Aperture Radar (InSAR) mission. To meet the goal of mapping forest structure and biomass contiguously and accurately, three types of fusion strategies have been investigated. First, a simulated GEDI-derived digital terrain model (DTM) was utilized to improve height estimation from TDX. Forest heights were initially derived from TDX coherence alone as a baseline using the widely used Random Volume over Ground (RVoG) scattering model. Here, assumptions about RVoG parameters – extinction coefficient (σ) and ground-to-volume amplitude ratio (µ) – were made. Using an external DTM derived from simulated GEDI lidar data, RVoG model was used to calculate spatially varied σ values and derived forest heights with better accuracy. TDX forest height estimation was further improved with the aid of simulated GEDI-derived DTM and canopy heights. The additional use of simulated GEDI canopy heights as RVoG input not just refined σ but also enabled the estimation of µ. Based on these parameters, forest heights were improved across three different forest types; biases were reduced from 1.7–3.8 m using only simulated GEDI DTMs to -0.9–1.1 m by using both simulated GEDI DTMs and canopy heights. Finally, wall-to-wall TDX heights were used to improve biomass estimates from simulated GEDI data over three contrasting forest types. When using simulated GEDI sampled observations alone, uncertainties were estimated statistically to be 9.0–19.9% at 1 km. These were improved to 5.2–11.7% at the same resolution by upscaling simulated GEDI footprint biomass with TDX heights. The GEDI/TDX data fusion also enabled the generation of biomass maps at a fine spatial resolution of 100 m, with uncertainties estimated to be 6.0–14.0%. Through the exploration of these fusion strategies, it has been demonstrated that a fusion-based mapping method could realize the generation of forest biomass products from GEDI with unprecedented resolutions and accuracies, while taking advantage of global seamless observations from TDX

Forest attributes mapping with SAR data in the romanian South-Eastern Carpathians requirements and outcomes

Esta tesis doctoral se centra en la estimación de variables forestales en la zona Sureste de los Cárpatos Rumanos a partir de imágenes de radar de apertura sintética. La investigación abarca parte del preprocesado de las imágenes, métodos de generación de mosaicos y la extracción de la cobertura de bosque, sus subtipos o su biomasa. La tesis se desarrolló en el Instituto Nacional de Investigación y Desarrollo Forestal Marín Dracea (INCDS) y la Universidad de Alcalá (UAH) gracias a varios proyectos: el proyecto EO-ROFORMON del INCDS (Prototyping an Earth-Observation based monitoring and forecasting system for the Romanian forests), y el proyecto EMAFOR de la UAH (Synthetic Aperture Radar (SAR) enabled Analysis Ready Data (ARD) cubes for efficient monitoring of agricultural and forested landscapes). El proyecto EO-ROFORMON fue financiado por la Autoridad Nacional para la Investigación Científica de Rumania y el Fondo Europeo de Desarrollo Regional. El proyecto EMAFOR fue financiado por la Comunidad Autónoma de Madrid (España). El objetivo de esta tesis es el desarrollo de algoritmos para la extracción de variables forestales de uso general como la cobertura, el tipo o la biomasa del bosque a partir de imagen de radar de apertura sintética. Para alcanzar dicho propósito se analizaron posibles fuentes de sesgo sistemático que podrían aparecer en zonas de montaña (ej., normalización topográfica, generación de mosaicos), y se aplicaron técnicas de aprendizaje de máquina para tareas de clasificación y regresión. La tesis contiene ocho secciones: una introducción, cinco publicaciones en revistas o actas de congresos indexados, una pendiente de publicación (quinto capítulo) y las conclusiones. La introducción contextualiza la importancia del bosque, cómo se recoge la información sobre su estado (ej., inventario forestal) y las iniciativas o marcos legislativos que requieren dicha información. A continuación, se describe cómo la teledetección puede complementar la información de inventario forestal, detallando el contexto histórico de las distintas tecnologías, su funcionamiento, y cómo pueden ser aplicadas para la extracción de información forestal. Por último, se describe la problemática y el monitoreo del bosque en Rumanía, detallando el objetivo de la tesis y su estructura. El primer capítulo analiza la influencia del modelo digital de elevaciones (MDE) en la calidad de la normalización topográfica, analizando tres MDE globales (SRTM, AW3D y TanDEM-X DEM) y uno nacional (PNOA-LiDAR). Los experimentos se basan en la comparación entre órbitas, con un MDE de referencia, y la variación del acierto en la clasificación dependiendo del MDE empleado para la normalización. Los resultados muestran una menor diferencia ente órbitas al utilizar un MDE con una mejor resolución (ej. TanDEM-X, PNOA-LIDAR), especialmente en el caso de zonas con fuertes pendientes o formas del terreno complejas, como pueden ser los valles. En zonas de alta montaña las imágenes de radar de apertura sintética (SAR) sufren frecuentes distorsiones. Estas distorsiones dependen de la geometría de adquisición, por lo que es posible combinar imágenes adquiridas desde varias órbitas para que la cobertura sea lo más completa posible. El segundo capítulo evalúa dos metodologías para la clasificación de usos del suelo utilizando datos de Sentinel-1 adquiridos desde varias órbitas. El primer método crea clasificaciones por órbita y las combina, mientras que el segundo genera un mosaico con datos de múltiples órbitas y lo clasifica. El acierto obtenido mediante combinación de clasificaciones es ligeramente mayor, mientras que la clasificación de mosaicos tiene importantes omisiones de las zonas boscosas debido a problemas en la normalización topográfica y a los efectos direccionales. El tercer capítulo se enfoca en separar la cobertura forestal de otras coberturas del suelo (urbano, vegetación baja, agua) analizando la utilidad de las variables basadas en la coherencia interferométrica. En él se realizan tres clasificaciones de máquina vector-soporte basadas en un conjunto concreto de variables. El primer conjunto contiene las estadísticas anuales de la retrodispersión (media y desviación típica anual), el segundo añade la coherencia a largo plazo (separación temporal mayor a un año), el tercero incluye las estadísticas de la coherencia a corto plazo (mínima separación temporal). Utilizar variables basadas en la coherencia aumenta el acierto de la clasificación hasta un 5% y reduce los errores de omisión de la cobertura forestal. El cuarto capítulo evalúa la posibilidad de detectar talas selectivas utilizando datos de Sentinel-1 y Sentinel-2. Sus resultados muestran que la detección resulta muy difícil debido a la saturación de los sensores y la confusión introducida por el efecto de la fenología. El quinto capítulo se centra en la clasificación de tipos de bosque basado en una serie temporal de datos Sentinel-1. Se basa en la creación de un conjunto de modelos que describen la relación entre la retrodispersión y el ángulo local de incidencia para un determinado tipo de bosque y fecha concreta. Para cada píxel se calcula el residuo respecto al modelo de cada uno de los tipos de bosque, acumulando dichos residuos a lo largo de la serie temporal. Hecho esto, cada píxel es asignado al tipo de bosque que acumula un menor residuo. Los resultados son prometedores, mostrando que frondosas y coníferas tienen un comportamiento distintivo, y que es posible separar ambos tipos de bosque con un alto grado de acierto. El sexto capítulo está dedicado a la estimación de biomasa utilizando datos Sentinel-1, ALOS PALSAR y regresión Random Forest. Se obtiene un error similar para ambos sensores a pesar de utilizar una banda diferente (band-C vs. -L), con poca reducción en el error cuando ambas bandas se utilizan conjuntamente. Sin embargo, el ajuste de un estimador adaptado a las condiciones locales de Rumanía sí ofreció una reducción de del error al ser comparado con las estimaciones globales de biomasa

Evaluating the potential of full-waveform lidar for mapping pan-tropical tree species richness

AIM: Mapping tree species richness across the tropics is of great interest for effective conservation and biodiversity management. In this study, we evaluated the potential of full‐waveform lidar data for mapping tree species richness across the tropics by relating measurements of vertical canopy structure, as a proxy for the occupation of vertical niche space, to tree species richness. LOCATION: Tropics. TIME PERIOD: Present. MAJOR TAXA STUDIED: Trees. METHODS: First, we evaluated the characteristics of vertical canopy structure across 15 study sites using (simulated) large‐footprint full‐waveform lidar data (22 m diameter) and related these findings to in‐situ tree species information. Then, we developed structure–richness models at the local (within 25–50 ha plots), regional (biogeographical regions) and pan‐tropical scale at three spatial resolutions (1.0, 0.25 and 0.0625 ha) using Poisson regression. RESULTS: The results showed a weak structure–richness relationship at the local scale. At the regional scale (within a biogeographical region) a stronger relationship between canopy structure and tree species richness across different tropical forest types was found, for example across Central Africa and in South America [R^{2} ranging from .44–.56, root mean squared difference as a percentage of the mean (RMSD%) ranging between 23–61%]. Modelling the relationship pan‐tropically, across four continents, 39% of the variation in tree species richness could be explained with canopy structure alone (R^{2} = .39 and RMSD% = 43%, 0.25‐ha resolution). MAIN CONCLUSIONS: Our results may serve as a basis for the future development of a set of structure–richness models to map high resolution tree species richness using vertical canopy structure information from the Global Ecosystem Dynamics Investigation (GEDI). The value of this effort would be enhanced by access to a larger set of field reference data for all tropical regions. Future research could also support the use of GEDI data in frameworks using environmental and spectral information for modelling tree species richness across the tropics

Characterizing tree species diversity in the tropics using full-waveform lidar data

Tree species diversity is of paramount value to maintain forest health and to ensure that forests are able to provide all vital functions, such as creating oxygen, that are needed for mankind to survive. Most of the world’s tree species grow in the tropical region, but many of them are threatened with extinction due to increasing natural and human-induced pressures on the environment. Mapping tree species diversity in the tropics is of high importance to enable effective conservation management of these highly diverse forests. This dissertation explores a new approach to mapping tree species diversity by using information on the vertical canopy structure derived from full-waveform lidar data. This approach is of particular interest in light of the recently launched Global Ecosystem Dynamics Investigation (GEDI), a full-waveform spaceborne lidar. First, successful derivation of vertical canopy structure metrics is ensured by comparing canopy profiles from airborne lidar data to those from terrestrial lidar. Then, the airborne canopy profiles were used to map five successional vegetation types in Lopé National Park in Gabon, Africa. Second, the relationship between vertical canopy structure and tree species richness was evaluated across four study sites in Gabon, which enabled mapping of tree species richness using canopy structure information from full-waveform lidar. Third, the relationship between canopy structure and tree species richness across the tropics was established using field and lidar data collected in 16 study sites across the tropics. Finally, it was evaluated how the methods and applications developed here could be adapted and used for mapping pan-tropical tree species diversity using future GEDI lidar data products