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

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

Uncertainties and Perspectives on Forest Height Estimates by Sentinel-1 Interferometry

Forest height is a key parameter in forestry. SAR interferometry (InSAR) techniques have been extensively adopted to retrieve digital elevation models (DEM) to give a representation of the continuous variation of the Earth’s topography, including forests. Unfortunately, InSAR has been proven to fail over vegetation due to low coherence values; therefore, all phase unwrapping algorithms tend to avoid these areas, making InSAR-derived DEM over vegetation unreliable. In this work, a sensitivity analysis was performed with the aim of properly initializing the relevant operational parameters (baseline and multilooking factor) to maximize the theoretical accuracy of the height difference between the forest and reference point. Some scenarios were proposed to test the resulting “optimal values”, as estimated at the previous step. A simple model was additionally proposed and calibrated, aimed at predicting the optimal baseline value (and therefore image pair selection) for height uncertainty minimization. All our analyses were conducted using free available data from the Copernicus Sentinel-1 mission to support the operational transfer into the forest sector. Finally, the potential uncertainty affecting resulting height measures was quantified, showing that a value lower than 5 m can be expected once all user-dependent parameters (i.e., baseline, multilooking factor, temporal baseline) are properly tuned

Estimation of change in forest variables using synthetic aperture radar

Large scale mapping of changes in forest variables is needed for both environmental monitoring, planning of climate actions and sustainable forest management. Remote sensing can be used in conjunction with field data to produce wall-to-wall estimates that are practically impossible to produce using traditional field surveys. Synthetic aperture radar (SAR) can observe the forest independent of sunlight, clouds, snow, or rain, providing reliable high frequency coverage. Its wavelength determines the interaction with the forest, where longer wavelengths interact with larger structures of the trees, and shorter wavelengths interact mainly with the top part of the canopy, meaning that it can be chosen to fit specific applications. This thesis contains five studies conducted on the Remningstorp test site in southern Sweden. Studies I – III predicted above ground biomass (AGB) change using long wavelength polarimetric P- (in I) and L-band (in I – III) SAR data. The differences between the bands were small in terms of prediction quality, and the HV polarization, just as for AGB state prediction, was the polarization channel most correlated with AGB change. A moisture correction for L-band data was proposed and evaluated, and it was found that certain polarimetric measures were better for predicting AGB change than all of the polarization channels together. Study IV assessed the detectability of silvicultural treatments in short wavelength TanDEM-X interferometric phase heights. In line with earlier studies, only clear cuts were unambiguously distinguishable. Study V predicted site index and stand age by fitting height development curves to time series of TanDEM-X data. Site index and age were unbiasedly predicted for untreated plots, and the RMSE would likely decrease with longer time series. When stand age was known, SI was predicted with an RMSE comparable to that of the field based measurements. In conclusion, this thesis underscores SAR data's potential for generalizable methods for estimation of forest variable changes

Temporal Characteristics of Boreal Forest Radar Measurements

Radar observations of forests are sensitive to seasonal changes, meteorological variables and variations in soil and tree water content. These phenomena cause temporal variations in radar measurements, limiting the accuracy of tree height and biomass estimates using radar data. The temporal characteristics of radar measurements of forests, especially boreal forests, are not well understood. To fill this knowledge gap, a tower-based radar experiment was established for studying temporal variations in radar measurements of a boreal forest site in southern Sweden. The work in this thesis involves the design and implementation of the experiment and the analysis of data acquired. The instrument allowed radar signatures from the forest to be monitored over timescales ranging from less than a second to years. A purpose-built, 50 m high tower was equipped with 30 antennas for tomographic imaging at microwave frequencies of P-band (420-450 MHz), L-band (1240-1375 MHz) and C-band (5250-5570 MHz) for multiple polarisation combinations. Parallel measurements using a 20-port vector network analyser resulted in significantly shorter measurement times and better tomographic image quality than previous tower-based radars. A new method was developed for suppressing mutual antenna coupling without affecting the range resolution. Algorithms were developed for compensating for phase errors using an array radar and for correcting for pixel-variant impulse responses in tomographic images. Time series results showed large freeze/thaw backscatter variations due to freezing moisture in trees. P-band canopy backscatter variations of up to 10 dB occurred near instantaneously as the air temperature crossed 0⁰C, with ground backscatter responding over longer timescales. During nonfrozen conditions, the canopy backscatter was very stable with time. Evidence of backscatter variations due to tree water content were observed during hot summer periods only. A high vapour pressure deficit and strong winds increased the rate of transpiration fast enough to reduce the tree water content, which was visible as 0.5-2 dB backscatter drops during the day. Ground backscatter for cross-polarised observations increased during strong winds due to bending tree stems. Significant temporal decorrelation was only seen at P-band during freezing, thawing and strong winds. Suitable conditions for repeat-pass L-band interferometry were only seen during the summer. C-band temporal coherence was high over timescales of seconds and occasionally for several hours for night-time observations during the summer. Decorrelation coinciding with high transpiration rates was observed at L- and C-band, suggesting sensitivity to tree water dynamics.The observations from this experiment are important for understanding, modelling and mitigating temporal variations in radar observables in forest parameter estimation algorithms. The results also are also useful in the design of spaceborne synthetic aperture radar missions with interferometric and tomographic capabilities. The results motivate the implementation of single-pass interferometric synthetic aperture radars for forest applications at P-, L- and C-band

First Demonstration of Space-Borne Polarization Coherence Tomography for Characterizing Hyrcanian Forest Structural Diversity

Structural diversity is recognized as a complementary aspect of biological diversity and plays a fundamental role in forest management, conservation, and restoration. Hence, the assessment of structural diversity has become a major effort in the primary international processes, dealing with biodiversity and sustainable forest management. Because of prohibitive costs associated with the ground measurements of forest structure, despite their high accuracy, space-borne polarization coherence tomography (PCT) can introduce an alternative approach given its ability to provide a vertical reflectivity profile and spatiotemporal resolutions related to detecting forest structural changes. In this study, for the first time ever, the potential of space-borne PCT was evaluated in a broad-leaved Hyrcanian forest of Iran over 308 circular sample plots with an area of 0.1 ha. Two aspects of horizontal structure diversity, including standard deviation of diameter at breast height (σdbh) and the number of trees (N), were predicted as important characteristics in wood production and biomass estimation. In addition, the performance of prediction algorithms, including multiple linear regression (MLR), k-nearest neighbors (k-NN), random forest (RF), and support vector regression (SVR) were compared. We addressed the issue of temporal decorrelation in space-borne PCT utilizing the single-pass TanDEM-X interferometer. The data were acquired in standard DEM mode with single polarization of HH. Consequently, airborne laser scanning (ALS) was used to estimate initial values of height hv and ground phase φ0. The Fourier–Legendre series was used to approximate the relative reflectivity profile of each pixel. To link the relative reflectivity profile averaged within each plot with corresponding ground measurements of σdbh and N, thirteen geometrical and physical parameters were defined (P1−P13). Leave-one-out cross validation (LOOCV) showed a better performance of k-NN than the other algorithms in predicting σdbh and N. It resulted in a relative root mean square error (rRMSE) of 32.80%, mean absolute error (MAE) of 4.69 cm, and R2* of 0.25 for σdbh, whereas only 22% of the variation in N was explained using the PCT algorithm with an rRMSE of 41.56%. This study revealed promising results utilizing TanDEM-X data even though the accuracy is still limited. Hence, an entire assessment of the used framework in characterizing the reflectivity profile and the possible effect of the scale is necessary for future studies

First Demonstration of Space-Borne Polarization Coherence Tomography for Characterizing Hyrcanian Forest Structural Diversity

Structural diversity is recognized as a complementary aspect of biological diversity and plays a fundamental role in forest management, conservation, and restoration. Hence, the assessment of structural diversity has become a major effort in the primary international processes, dealing with biodiversity and sustainable forest management. Because of prohibitive costs associated with the ground measurements of forest structure, despite their high accuracy, space-borne polarization coherence tomography (PCT) can introduce an alternative approach given its ability to provide a vertical reflectivity profile and spatiotemporal resolutions related to detecting forest structural changes. In this study, for the first time ever, the potential of space-borne PCT was evaluated in a broad-leaved Hyrcanian forest of Iran over 308 circular sample plots with an area of 0.1 ha. Two aspects of horizontal structure diversity, including standard deviation of diameter at breast height (σdbh) and the number of trees (N), were predicted as important characteristics in wood production and biomass estimation. In addition, the performance of prediction algorithms, including multiple linear regression (MLR), k-nearest neighbors (k-NN), random forest (RF), and support vector regression (SVR) were compared. We addressed the issue of temporal decorrelation in space-borne PCT utilizing the single-pass TanDEM-X interferometer. The data were acquired in standard DEM mode with single polarization of HH. Consequently, airborne laser scanning (ALS) was used to estimate initial values of height hv and ground phase φ0. The Fourier–Legendre series was used to approximate the relative reflectivity profile of each pixel. To link the relative reflectivity profile averaged within each plot with corresponding ground measurements of σdbh and N, thirteen geometrical and physical parameters were defined (P1−P13). Leave-one-out cross validation (LOOCV) showed a better performance of k-NN than the other algorithms in predicting σdbh and N. It resulted in a relative root mean square error (rRMSE) of 32.80%, mean absolute error (MAE) of 4.69 cm, and R2* of 0.25 for σdbh, whereas only 22% of the variation in N was explained using the PCT algorithm with an rRMSE of 41.56%. This study revealed promising results utilizing TanDEM-X data even though the accuracy is still limited. Hence, an entire assessment of the used framework in characterizing the reflectivity profile and the possible effect of the scale is necessary for future studies

Temporal survey of P- A nd L-band polarimetric backscatter in boreal forests

Environmental conditions and seasonal variations affect the backscattered radar signal from a forest. This potentially causes errors in a biomass retrieval scheme using data from the synthetic aperture radar (SAR) data. A better understanding of these effects and the electromagnetic scattering mechanisms in forests is required to improve biomass estimation algorithms for current and upcoming P- A nd L-band SAR missions. In this paper, temporal changes in HH-, VV-, and HV-polarized P- A nd L-band radar backscatter and temporal coherence from a boreal forest site are analyzed in relation to environmental parameters. The radar data were collected from a stand of mature Norway spruce ( Picea abies (L.) Karst.) with an above-ground biomass of approximately 250 tons/ha at intervals of 5 min from January to August 2017 using the BorealScat tower-based scatterometer. It was observed that subzero temperatures during the winters cause large variations (4 to 10 dB) in P- A nd L-band backscatter, for which the HH/VV backscatter ratio offered some mitigation. High wind speeds were also seen to cause deviations in the average backscatter at P-band due to decreased double-bounce scattering. Severe temporal decorrelation was observed at L-band over timescales of days or more, whereas the P-band temporal coherence remained high (> 0.9) for at least a month neglecting windy periods. Temporal coherence at P-band was highest during night times when wind speeds are low

The SAR Handbook: Comprehensive Methodologies for Forest Monitoring and Biomass Estimation

This Synthetic Aperture Radar (SAR) handbook of applied methods for forest monitoring and biomass estimation has been developed by SERVIR in collaboration with SilvaCarbon to address pressing needs in the development of operational forest monitoring services. Despite the existence of SAR technology with all-weather capability for over 30 years, the applied use of this technology for operational purposes has proven difficult. This handbook seeks to provide understandable, easy-to-assimilate technical material to remote sensing specialists that may not have expertise on SAR but are interested in leveraging SAR technology in the forestry sector

Estimation of biophysical parameters in boreal forests from ERS and JERS SAR interferometry

The thesis describes investigations concerning the evaluation of ERS and JERS SAR images and repeat-pass interferometric SAR images for the retrieval of biophysical parameters in boreal forests. The availability of extensive data sets of images over several test sites located in Sweden, Finland and Siberia has allowed analysis of temporal dynamics of ERS and JERS backscatter and coherence, and of ERS interferometric phase. Modelling of backscatter, coherence and InSAR phase has been performed by means of the Water Cloud Model (WCM) and the Interferometric Water Cloud Model (IWCM); sensitivity analysis and implications for the retrieval of forest biophysical parameters have been thoroughly discussed. Model inversion has been carried out for stem volume retrieval using ERS coherence, ERS backscatter and JERS backscatter, whereas for tree height estimation the ERS interferometric phase has been used. Multi-temporal combination of ERS coherence images, and to a lesser extent of JERS backscatter images, can provide stem volume estimates comparable to stand-wise ground-based measurements. Since the information content of the interferometric phase is strongly degraded by phase noise and uncorrected atmospheric artefacts, the retrieved tree height shows large errors