Utilization of bistatic TanDEM-X data to derive land cover information

Forests have significance as carbon sink in climate change. Therefore, it is of high importance to track land use changes as well as to estimate the state as carbon sink. This is useful for sustainable forest management, land use planning, carbon modelling, and support to implement international initiatives like REDD+ (Reducing Emissions from Deforestation and Degradation). A combination of field measurements and remote sensing seems most suitable to monitor forests. Radar sensors are considered as high potential due to the weather and daytime independence. TanDEM-X is a interferometric SAR (synthetic aperture radar) mission in space and can be used for land use monitoring as well as estimation of biophysical parameters. TanDEM-X is a X-band system resulting in low penetration depth into the forest canopy. Interferometric information can be useful, whereas the low penetration can be considered as an advantage. The interferometric height is assumable as canopy height, which is correlated with forest biomass. Furthermore, the interferometric coherence is mainly governed by volume decorrelation, whereas temporal decorrelation is minimized. This information can be valuable for quantitative estimations and land use monitoring. The interferometric coherence improved results in comparison to land use classifications without coherence of about 10% (75% vs. 85%). Especially the differentiation between forest classes profited from coherence. The coherence correlated with aboveground biomass in a R² of about 0.5 and resulted in a root mean square error (RSME) of 14%. The interferometric height achieved an even higher correlation with the biomass (R²=0.68) resulting in cross-validated RMSE of 7.5%. These results indicated that TanDEM-X can be considered as valuable and consistent data source for forest monitoring. Especially interferometric information seemed suitable for biomass estimation

Estimating mangrove aboveground biomass from airborne LiDAR data: a case study from the Zambezi River delta

Mangroves are ecologically and economically important forested wetlands with the highest carbon (C) density of all terrestrial ecosystems. Because of their exceptionally large C stocks and importance as a coastal buffer, their protection and restoration has been proposed as effective mitigation strategy for climate change. The inclusion of mangroves in mitigation strategies requires the quantification of C stocks (both above and belowground) and changes to accurately calculate emissions and sequestration. A growing number of countries are becoming interested in using mitigation initiatives, such as REDD+, in these unique coastal forests. However, it is not yet clear how methods to measure C traditionally used for other ecosystems can be modified to estimate biomass in mangroves with the precision and accuracy needed for these initiatives. Airborne lidar (ALS) data has often been proposed as the most accurate way for larger-scale assessments but, the application of ALS for coastal wetlands is scarce, primarily due to a lack of contemporaneous ALS and field measurements. Here, we evaluated the variability in field and lidar-based estimates of aboveground biomass (AGB) through the combination of different local and regional allometric models and standardized height metrics that are comparable across spatial resolutions and sensor types. The end result being a simplified approach for accurately estimating mangrove AGB at large-scales and determining the uncertainty by combining multiple allometric models. We then quantified wall-to-wall aboveground biomass stocks of a tall mangrove forest in the Zambezi Delta, Mozambique. Our results indicate that the Lidar H100 height metric correlates well with AGB estimates, with R2 between 0.80 and 0.88 and RMSE of 33% or less. When comparing lidar H100 AGB derived from three allometric models, mean AGB values range from 192 Mg. ha-1 up to 252 Mg. ha-1. We suggest the best model to predict AGB was based on the East Africa specific allometry and a power based regression that used Lidar H100 as the height input with a R2 of 0.85 and a RMSE of 122 Mg.ha-1 or 33%. The total AGB of the lidar inventoried mangrove area (6654 ha) was 1,350,902 Mg with a mean AGB 203 Mg. ha-1. Because the allometry suggested here was developed using standardized height metrics, it is recommended that the models can generate AGB estimates using other remote sensing instruments that are more readily accessible over other mangrove ecosystems on a large scale, and as part of future carbon monitoring efforts in mangroves

Forest disturbance and recovery: A general review in the context of spaceborne remote sensing of impacts on aboveground biomass and canopy structure

Abrupt forest disturbances generating gaps \u3e0.001 km2 impact roughly 0.4–0.7 million km2a−1. Fire, windstorms, logging, and shifting cultivation are dominant disturbances; minor contributors are land conversion, flooding, landslides, and avalanches. All can have substantial impacts on canopy biomass and structure. Quantifying disturbance location, extent, severity, and the fate of disturbed biomass will improve carbon budget estimates and lead to better initialization, parameterization, and/or testing of forest carbon cycle models. Spaceborne remote sensing maps large-scale forest disturbance occurrence, location, and extent, particularly with moderate- and fine-scale resolution passive optical/near-infrared (NIR) instruments. High-resolution remote sensing (e.g., ∼1 m passive optical/NIR, or small footprint lidar) can map crown geometry and gaps, but has rarely been systematically applied to study small-scale disturbance and natural mortality gap dynamics over large regions. Reducing uncertainty in disturbance and recovery impacts on global forest carbon balance requires quantification of (1) predisturbance forest biomass; (2) disturbance impact on standing biomass and its fate; and (3) rate of biomass accumulation during recovery. Active remote sensing data (e.g., lidar, radar) are more directly indicative of canopy biomass and many structural properties than passive instrument data; a new generation of instruments designed to generate global coverage/sampling of canopy biomass and structure can improve our ability to quantify the carbon balance of Earth\u27s forests. Generating a high-quality quantitative assessment of disturbance impacts on canopy biomass and structure with spaceborne remote sensing requires comprehensive, well designed, and well coordinated field programs collecting high-quality ground-based data and linkages to dynamical models that can use this information

The NASA AfriSAR campaign: Airborne SAR and lidar measurements of tropical forest structure and biomass in support of current and future space missions

International audienceIn 2015 and 2016, the AfriSAR campaign was carried out as a collaborative effort among international space and National Park agencies (ESA, NASA, ONERA, DLR, ANPN and AGEOS) in support of the upcoming ESA BIOMASS, NASA-ISRO Synthetic Aperture Radar (NISAR) and NASA Global Ecosystem Dynamics Initiative (GEDI) missions. The NASA contribution to the campaign was conducted in 2016 with the NASA LVIS (Land Vegetation and Ice Sensor) Lidar, the NASA L-band UAVSAR (Uninhabited Aerial Vehicle Synthetic Aperture Radar). A central motivation for the AfriSAR deployment was the common AGBD estimation requirement for the three future spaceborne missions, the lack of sufficient airborne and ground calibration data covering the full range of ABGD in tropical forest systems, and the intercomparison and fusion of the technologies. During the campaign, over 7000 km2 of waveform Lidar data from LVIS and 30,000 km2 of UAVSAR data were collected over 10 key sites and transects. In addition, field measurements of forest structure and biomass were collected in sixteen 1-hectare sized plots. The campaign produced gridded Lidar canopy structure products, gridded aboveground biomass and associated uncertainties, Lidar based vegetation canopy cover profile products, Polarimetric Interferometric SAR and Tomographic SAR products and field measurements. Our results showcase the types of data products and scientific results expected from the spaceborne Lidar and SAR missions; we also expect that the AfriSAR campaign data will facilitate further analysis and use of waveform lidar and multiple baseline polarimetric SAR datasets for carbon cycle, biodiversity, water resources and more applications by the greater scientific community

Estimation of Forest Biomass From Two-Level Model Inversion of Single-Pass InSAR Data

A model for aboveground biomass estimation from single-pass interferometric synthetic aperture radar (InSAR) data is presented. Forest height and canopy density estimates Delta h and eta(0), respectively, obtained from two-level model (TLM) inversion, are used as biomass predictors. Eighteen bistatic VV-polarized TanDEM-X (TDM) acquisitions are used, made over two Swedish test sites in the summers of 2011, 2012, and 2013 (nominal incidence angle: 41 degrees; height-of-ambiguity: 32-63 m). Remningstorp features a hemiboreal forest in southern Sweden, with flat topography and where 32 circular plots have been sampled between 2010 and 2011 (area: 0.5 ha; biomass: 42-242 t/ha; height: 14-32 m). Krycklan features a boreal forest in northern Sweden, 720-km north-northeast from Remningstorp, with significant topography and where 31 stands have been sampled in 2008 (area: 2.4-26.3 ha; biomass: 23-183 t/ha; height: 7-21 m). A high-resolution digital terrain model has been used as ground reference during InSAR processing. For the aforementioned plots and stands and if the same acquisition is used for model training and validation, the new model explains 65%-89% of the observed variance, with root-mean-square error (RMSE) of 12%-19% (median: 15%). By fixing two of the three model parameters, accurate biomass estimation can also be done when different acquisitions or different test sites are used for model training and validation, with RMSE of 12%-56% (median: 17%). Compared with a simple scaling model computing biomass from the phase center elevation above ground, the proposed model shows significantly better performance in Remningstorp, as it accounts for the large canopy density variations caused by active management. In Krycklan, the two models show similar performance

Biomass Representation in Synthetic Aperture Radar Interferometry Data Sets

This work makes an attempt to explain the origin, features and potential applications of the elevation bias of the synthetic aperture radar interferometry (InSAR) datasets over areas covered by vegetation. The rapid development of radar-based remote sensing methods, such as synthetic aperture radar (SAR) and InSAR, has provided an alternative to the photogrammetry and LiDAR for determining the third dimension of topographic surfaces. The InSAR method has proved to be so effective and productive that it allowed, within eleven days of the space shuttle mission, for acquisition of data to develop a three-dimensional model of almost the entire land surface of our planet. This mission is known as the Shuttle Radar Topography Mission (SRTM). Scientists across the geosciences were able to access the great benefits of uniformity, high resolution and the most precise digital elevation model (DEM) of the Earth like never before for their a wide variety of scientific and practical inquiries. Unfortunately, InSAR elevations misrepresent the surface of the Earth in places where there is substantial vegetation cover. This is a systematic error of unknown, yet limited (by the vertical extension of vegetation) magnitude. Up to now, only a limited number of attempts to model this error source have been made. However, none offer a robust remedy, but rather partial or case-based solutions. More work in this area of research is needed as the number of airborne and space-based InSAR elevation models has been steadily increasing over the last few years, despite strong competition from LiDAR and optical methods. From another perspective, however, this elevation bias, termed here as the “biomass impenetrability”, creates a great opportunity to learn about the biomass. This may be achieved due to the fact that the impenetrability can be considered a collective response to a few factors originating in 3D space that encompass the outermost boundaries of vegetation. The biomass, presence in InSAR datasets or simply the biomass impenetrability, is the focus of this research. The report, presented in a sequence of sections, gradually introduces terminology, physical and mathematical fundamentals commonly used in describing the propagation of electromagnetic waves, including the Maxwell equations. The synthetic aperture radar (SAR) and InSAR as active remote sensing methods are summarised. In subsequent steps, the major InSAR data sources and data acquisition systems, past and present, are outlined. Various examples of the InSAR datasets, including the SRTM C- and X-band elevation products and INTERMAP Inc. IFSAR digital terrain/surface models (DTM/DSM), representing diverse test sites in the world are used to demonstrate the presence and/or magnitude of the biomass impenetrability in the context of different types of vegetation – usually forest. Also, results of investigations carried out by selected researchers on the elevation bias in InSAR datasets and their attempts at mathematical modelling are reviewed. In recent years, a few researchers have suggested that the magnitude of the biomass impenetrability is linked to gaps in the vegetation cover. Based on these hints, a mathematical model of the tree and the forest has been developed. Three types of gaps were identified; gaps in the landscape-scale forest areas (Type 1), e.g. forest fire scares and logging areas; a gap between three trees forming a triangle (Type 2), e.g. depending on the shape of tree crowns; and gaps within a tree itself (Type 3). Experiments have demonstrated that Type 1 gaps follow the power-law density distribution function. One of the most useful features of the power-law distributed phenomena is their scale-independent property. This property was also used to model Type 3 gaps (within the tree crown) by assuming that these gaps follow the same distribution as the Type 1 gaps. A hypothesis was formulated regarding the penetration depth of the radar waves within the canopy. It claims that the depth of penetration is simply related to the quantisation level of the radar backscattered signal. A higher level of bits per pixels allows for capturing weaker signals arriving from the lower levels of the tree crown. Assuming certain generic and simplified shapes of tree crowns including cone, paraboloid, sphere and spherical cap, it was possible to model analytically Type 2 gaps. The Monte Carlo simulation method was used to investigate relationships between the impenetrability and various configurations of a modelled forest. One of the most important findings is that impenetrability is largely explainable by the gaps between trees. A much less important role is played by the penetrability into the crown cover. Another important finding is that the impenetrability strongly correlates with the vegetation density. Using this feature, a method for vegetation density mapping called the mean maximum impenetrability (MMI) method is proposed. Unlike the traditional methods of forest inventories, the MMI method allows for a much more realistic inventory of vegetation cover, because it is able to capture an in situ or current situation on the ground, but not for areas that are nominally classified as a “forest-to-be”. The MMI method also allows for the mapping of landscape variation in the forest or vegetation density, which is a novel and exciting feature of the new 3D remote sensing (3DRS) technique. Besides the inventory-type applications, the MMI method can be used as a forest change detection method. For maximum effectiveness of the MMI method, an object-based change detection approach is preferred. A minimum requirement for the MMI method is a time-lapsed reference dataset in the form, for example, of an existing forest map of the area of interest, or a vegetation density map prepared using InSAR datasets. Preliminary tests aimed at finding a degree of correlation between the impenetrability and other types of passive and active remote sensing data sources, including TerraSAR-X, NDVI and PALSAR, proved that the method most sensitive to vegetation density was the Japanese PALSAR - L-band SAR system. Unfortunately, PALSAR backscattered signals become very noisy for impenetrability below 15 m. This means that PALSAR has severe limitations for low loadings of the biomass per unit area. The proposed applications of the InSAR data will remain indispensable wherever cloud cover obscures the sky in a persistent manner, which makes suitable optical data acquisition extremely time-consuming or nearly impossible. A limitation of the MMI method is due to the fact that the impenetrability is calculated using a reference DTM, which must be available beforehand. In many countries around the world, appropriate quality DTMs are still unavailable. A possible solution to this obstacle is to use a DEM that was derived using P-band InSAR elevations or LiDAR. It must be noted, however, that in many cases, two InSAR datasets separated by time of the same area are sufficient for forest change detection or similar applications

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

An Effective Method for InSAR Mapping of Tropical Forest Degradation in Hilly Areas

Current satellite remote sensing methods struggle to detect and map forest degradation, which is a critical issue as it is likely a major and growing source of carbon emissions and biodiveristy loss. TanDEM-X InSAR phase height is a promising variable for measuring forest disturbances, as it is closely related to the mean canopy height, and thus should decrease if canopy trees are removed. However, previous research has focused on relatively flat terrains, despite the fact that much of the world's remaining tropical forests are found in hilly areas, and this inevitably introduces artifacts in sideways imaging systems. In this paper, we find a relationship between InSAR phase height and aboveground biomass change in four selectively logged plots in a hilly region of central Gabon. We show that minimising multilooking prior to the calculation of InSAR phase height on a pixel-by-pixel basis. This shows that TanDEM-X InSAR can measure the magnitude of degradation, and that topographic effects can be mitigated if data from multiple SAR viewing geometries are available