Search

Accurate estimation of forest terrain and canopy height is crucial for timely understanding of forest growth. Gao Fen-7 (GF-7) Satellite is China’s first sub-meter-level three-dimensional (3D) mapping satellite for civilian use, which was equipped with a two-line-array stereo mapping camera and a laser altimeter system that can provide stereo images and full waveform LiDAR data simultaneously. Most of the existing studies have concentrated on evaluating the accuracy of GF-7 for topographic survey in bare land, but few have in-depth studied its ability to measure forest terrain elevation and canopy height. The purpose of this study is to evaluate the potential of GF-7 LiDAR and stereo image for forest terrain and height measurement. The Airborne Laser Scanning (ALS) data were utilized to generate reference terrain and forest vertical information. The validation test was conducted in Pu’er City, Yunnan Province of China, and encouraging results have obtained. The GF-7 LiDAR data obtained the accuracy of forest terrain elevation with RMSE of 8.01 m when 21 available laser footprints were used for results verification; meanwhile, when it was used to calculate the forest height, R2 of 0.84 and RMSE of 3.2 m were obtained although only seven effective footprints were used for result verification. The canopy height values obtained from GF-7 stereo images have also been proven to have high accuracy with the resolution of 20 m × 20 m compared with ALS data (R2 = 0.88, RMSE = 2.98 m). When the results were verified at the forest sub-compartment scale that taking into account the forest types, further higher accuracy (R2 = 0.96, RMSE = 1.23 m) was obtained. These results show that GF-7 has considerable application potential in forest resources monitoring.

The role of tropical forests in the global carbon budget remains controversial, as carbon emissions from deforestation are highly uncertain. This high uncertainty arises from the use of either fixed forest carbon stock density or maps generated from satellite-based optical reflectance with limited sensitivity to biomass to generate accurate estimates of emissions from deforestation. New space missions aiming to accurately map the carbon stock density rely on direct measurements of the spatial structures of forests using lidar and radar. We found that lost forests are special cases, and their spatial structures can be directly measured by combining archived data acquired before and after deforestation by space missions principally aimed at measuring topography. Thus, using biomass mapping, we obtained new estimates of carbon loss from deforestation ahead of forthcoming space missions. Here, using a high-resolution map of forest loss and the synergy of radar and lidar to estimate the aboveground biomass density of forests, we found that deforestation in the 2000s in Latin America, one of the severely deforested regions, mainly occurred in forests with a significantly lower carbon stock density than typical mature forests. Deforestation areas with carbon stock densities lower than 20.0, 50.0, and 100.0 Mg C/ha accounted for 42.1%, 62.0%, and 83.3% of the entire deforested area, respectively. The average carbon stock density of lost forests was only 49.13 Mg C/ha, which challenges the current knowledge on the carbon stock density of lost forests (with a default value 100 Mg C/ha according to the Intergovernmental Panel on Climate Change Tier 1 estimates, or approximately 112 Mg C/ha used in other studies). This is demonstrated over both the entire region and the footprints of the spaceborne lidar. Consequently, our estimate of carbon loss from deforestation in Latin America in the 2000s was 253.0 ± 21.5 Tg C/year, which was considerably less than existing remote-sensing-based estimates, namely 400–600 Tg C/year. This indicates that forests in Latin America were most likely not a net carbon source in the 2000s compared to established carbon sinks. In previous studies, considerable effort has been devoted to rectify the underestimation of carbon sinks; thus, the overestimation of carbon emissions should be given sufficient consideration in global carbon budgets. Our results also provide solid evidence for the necessity of renewing knowledge on the role of tropical forests in the global carbon budget in the future using observations from new space missions.

With the advent of very high-resolution (VHR) imaging satellites, it is possible to measure the heights of forest stands or even individual trees more accurately. However, the accurate geometric processing of VHR images depends on ground control points (GCPs). Collecting GCPs through fieldwork is time-consuming and labor-intensive, which presents great challenges for regional applications in remote or mountainous regions, particularly for international applications. This study proposes a promising approach that leverages GF-7 VHR stereoscopic images and Google Earth’s multi-temporal historical imagery to accurately extract forest heights without the need for fieldworks. Firstly, an algorithm is proposed to collect GCPs using Multi-temporal Averaging of historical imagery provided by Google Earth (GE), known as MAGE. Digital surface model (DSM) is then derived using GF-7 stereoscopic imagery and MAGE GCPs in Switzerland. Forest heights are finally extracted by subtracting ground surface elevations from GF-7 DSM. Results show that absolute coordinate errors of MAGE GCPs are less than 2.0 m. The root mean square error (RMSE) of forest heights extracted from GF-7 DSM, derived using the original geolocation model, is 12.3 m, and the determination coefficient (R2) of linear estimation model is 0.72. When the geolocation model is optimized using MAGE GCPs, the RMSE is reduced to 1.5 m and the R2 increases to 0.95. These results not only demonstrate the effectiveness of MAGE GCPs but, more importantly, also reveal the significance of precise geometric processing of VHR stereoscopic imagery in forest height estimations.

Canopy cover is an important parameter affecting forest succession, carbon fluxes, and wildlife habitats. Several global maps with different spatial resolutions have been produced based on satellite images, but facing the deficiency of reliable references for accuracy assessments. The rapid development of unmanned aerial vehicle (UAV) equipped with consumer-grade camera enables the acquisition of high-resolution images at low cost, which provides the research community a promising tool to collect reference data. However, it is still a challenge to distinguish tree crowns and understory green vegetation based on the UAV-based true color images (RGB) due to the limited spectral information. In addition, the canopy height model (CHM) derived from photogrammetric point clouds has also been used to identify tree crowns but limited by the unavailability of understory terrain elevations. This study proposed a simple method to distinguish tree crowns and understories based on UAV visible images, which was referred to as BAMOS for convenience. The central idea of the BAMOS was the synergy of spectral information from digital orthophoto map (DOM) and structural information from digital surface model (DSM). Samples of canopy covers were produced by applying the BAMOS method on the UAV images collected at 77 sites with a size of about 1.0 km2 across Daxing’anling forested area in northeast of China. Results showed that canopy cover extracted by the BAMOS method was highly correlated to visually interpreted ones with correlation coefficient (r) of 0.96 and root mean square error (RMSE) of 5.7%. Then, the UAV-based canopy covers served as references for assessment of satellite-based maps, including MOD44B Version 6 Vegetation Continuous Fields (MODIS VCF), maps developed by the Global Land Cover Facility (GLCF) and by the Global Land Analysis and Discovery laboratory (GLAD). Results showed that both GLAD and GLCF canopy covers could capture the dominant spatial patterns, but GLAD canopy cover tended to miss scattered trees in highly heterogeneous areas, and GLCF failed to capture non-tree areas. Most important of all, obvious underestimations with RMSE about 20% were easily observed in all satellite-based maps, although the temporal inconsistency with references might have some contributions.

Directly mapping impervious surface area (ISA) at national and global scales using nighttime light data is a challenge due to the complexity of land surface components and the impacts of unbalanced economic conditions. Previous research mainly used the coarse spatial resolution Defense Meteorological Satellite Program’s Operational Linescan System (DMSP OLS) and Moderate Resolution Imaging Spectroradiometer (MODIS), normalized difference vegetation index (NDVI) data for ISA mapping; the improved spatial resolution and data quality in the Suomi National Polar-orbiting Partnership, Visible Infrared Imaging Radiometer Suite’s Day/Night Band (VIIRS DNB) and in Proba-V data provide a new opportunity to accurately map ISA distribution at the national scale, which has not been explored yet. This research aimed to develop a new index – modified impervious surface index (MISI) – based on VIIRS DNB and Proba-V data to improve ISA estimation and to compare the results with those from the combination of VIIRS DNB and MODIS NDVI data. Landsat data were used to develop ISA data for the typical sites for use as reference data. Regression analysis was used to establish the ISA estimation model in which the dependent variable was from the Landsat data and the independent variable was from the MISI, as well as the previously used Large-scale Impervious Surface Index (LISI). The results indicate that the major error is from the very small or very large proportion of ISA in a unit; improvement of spatial resolution through use of higher spatial resolution nighttime light data (e.g., VIIRS DNB) or NDVI (e.g., Proba-V NDVI) data is an effective approach to improve ISA estimation. Although different indices for the combination of nighttime light and NDVI data have been used, the MISI is especially valuable for reducing the estimation errors for the regions with a small or large ISA proportion.

Waveform broadening effects of large-footprint lidar caused by terrain slopes are still a great challenge limiting the estimation accuracy of forest aboveground biomass (AGB) over mountainous areas. Slope-adaptive metrics of waveforms were proposed in our previous studies. However, its validation was limited by the unavailability of enough reference data. This study made full validation of slope-adaptive metrics using data acquired by the Global Ecosystem Dynamics Investigation (GEDI) mission, meanwhile exploring GEDI waveforms on estimations of forest AGB. Three types of waveform metrics were employed, including slope-adaptive metrics (RHT), typical height metrics relative to ground peaks (RH), and waveform parameters (WP). In addition to terrain slopes, two other factors were also explored including the geolocation issue and signal start and ending points of waveforms. Results showed that footprint geolocations in the first version GEDI data products were shifted to the left forward of nominal geolocations with a distance of about 24 m~30 m and were substantially corrected in the second version; the fourth and fifth groups of signal start and ending points of waveforms had worse performance than the rest of the four groups because they used the maximum and minimum signal thresholds, respectively. Taking airborne laser scanner (ALS) data as reference, the root mean square error (RMSE) of terrain slopes extracted from the digital elevation model of the shuttle radar topography mission (SRTM DEM) was about 3°. The coefficients of determination (R2) of estimation models of forest AGB based on RH metrics were improved from 0.48 to 0.68 with RMSE decreased from 19.7 Mg/ha to 15.4 Mg/ha by the second version geolocations. The RHT and WP metrics gave the best and the worst estimation accuracy, respectively. RHT further improved R2 to 0.77 and decreased RMSE to 13.0 Mg/ha using terrain slopes extracted from SRTM DEM with a resolution of 1 arc second. R2 of estimation models based on RHT was finally improved to 0.8 with RMSE decreased to 11.7 Mg/ha using exact terrain slopes from ALS data. This study demonstrated the great potential of slope-adaptive metrics of GEDI waveforms on estimations of forest aboveground biomass over mountainous areas.

Sentinel-1 provides an extraordinary opportunity to explore the temporal behavior of backscattering of C-band synthetic aperture radar (SAR) due to its unique capability of successive observations every 12 days. This study reported new findings on the critical temporal features of Sentinel-1 backscattering over mountainous forested areas in northeast China and their application in regional forest mapping. Two interesting phenomena were discovered through the analysis of 450 scenes of images acquired by Sentinel-1A or Sentinel-1B over an area of 318,898.62 km2. The first phenomenon was that the dates of the largest drops of backscattering coefficients over forest and non-forest plots were different during the transition from autumn to winter. The largest drop of non-forest plots occurred around the date of the minimum temperature decreasing below 0 °C, while that of forest plots occurred around the date of the maximum temperature decreasing below 0 °C. The second phenomenon was that at the dates where these two drops occurred, the magnitude of the drop was negatively correlated with the forest canopy coverage for the first date and positively correlated for the second date. Based on these two phenomena, two methods for the forest mapping, referred to as the direct method and the indirect method, were proposed using only three dates of Sentinel-1 images, i.e., Date1: before the minimum temperature decreased below 0 °C, Date2: after the minimum temperature decreased below 0 °C but before the maximum temperature decreased below 0 °C, and Date3: after the maximum temperature decreased below 0 °C. The results showed that the overall accuracy of the forest map produced by the direct method was 93.60%, while that by the indirect method was 93.80%. Their accuracies were comparable with those of forest maps derived from publicly released land cover maps, which was approximately 94.42% for the best one. This study proposed a new way to do large-scale forest mapping in annually frozen regions using as few Sentinel-1 SAR images as possible.

Accurate estimates of forest aboveground biomass (AGB) after anthropogenic disturbance could reduce uncertainties in the carbon budget of terrestrial ecosystems and provide critical information to policy makers. Yet, the loss of carbon due to forest disturbance and the gain from post-disturbance recovery have not been sufficiently assessed. In this study, a sensitivity analysis was first conducted to investigate: (1) the influence of incidence angle and soil moisture on Synthetic Aperture Radar (SAR) backscatter; (2) the feasibility of cross-image normalization between multi-temporal and multi-sensor SAR data; and (3) the possibility of applying normalized backscatter data to detect forest biomass changes. An empirical model was used to reduce incidence angle effects, followed by cross-image normalization procedure to lessen soil moisture effect. Changes in forest biomass at medium spatial resolution (100 m) were mapped using both spaceborne and airborne SAR data. Results indicate that (1) the effect of incidence angle on SAR backscatter could be reduced to less than 1 dB by the correction model for airborne SAR data; (2) over 50% of the changes in SAR backscatter due to soil moisture could be eliminated by the cross-image normalization procedure; and (3) forest biomass changes greater than 100 Mg·ha−1 or above 50% of 150 Mg·ha−1 are detectable using cross-normalized SAR data.

Forest aboveground biomass (AGB) was mapped throughout China using large footprint LiDAR waveform data from the Geoscience Laser Altimeter System (GLAS) onboard NASA’s Ice, Cloud, and land Elevation Satellite (ICESat), Moderate Resolution Imaging Spectro-radiometer (MODIS) imagery and forest inventory data. The entire land of China was divided into seven zones according to the geographic characteristics of the forests. The forest AGB prediction models were separately developed for different forest types in each of the seven forest zones at GLAS footprint level from GLAS waveform parameters and biomass derived from height and diameter at breast height (DBH) field observation. Some waveform parameters used in the prediction models were able to reduce the effects of slope on biomass estimation. The models of GLAS-based biomass estimates were developed by using GLAS footprints with slopes less than 20° and slopes ≥ 20°, respectively. Then, all GLAS footprint biomass and MODIS data were used to establish Random Forest regression models for extrapolating footprint AGB to a nationwide scale. The total amount of estimated AGB in Chinese forests around 2006 was about 12,622 Mt vs. 12,617 Mt derived from the seventh national forest resource inventory data. Nearly half of all provinces showed a relative error (%) of less than 20%, and 80% of total provinces had relative errors less than 50%.

The Global Digital Elevation Model produced from stereo images of Advanced Spaceborne Thermal Emission and Reflection Radiometer data (ASTER GDEM) covers land surfaces between latitudes of 83°N and 83°S. The Phased Array type L-band Synthetic Aperture Radar (PALSAR) onboard Advanced Land Observing Satellite (ALOS) collected many SAR images since it was launched on 24 January 2006. The combination of ALOS/PALSAR interferometric data and ASTER GDEM should provide the penetration depth of SAR data assuming ASTER GDEM was the elevation of vegetation canopy top. It would be correlated with forest biomass because penetration depth could be affected by forest density and forest canopy height. Their combination held great promises for the forest biomass mapping over large area. The feasibility of forest biomass mapping through the data synthesis of ALOS/PALSAR InSAR data and ASTER GDEM was investigated in this study. A procedure for the extraction of penetration depth was firstly proposed. Then three models were built for biomass estimation: (I) model only using backscattering coefficients of ALOS/PALSAR data; (II) model only using penetration depth; (III) model using both of them. The biomass estimated from Lidar data was taken as reference data to evaluate the three different models. The results showed that the combination of backscattering coefficients and penetration depth gave the best accuracy. The forest disturbance has to be considered in forest biomass estimation because of the long time span of ASTER data for generating ASTER GDEM. The spatial homogeneity could be used to improve estimation accuracy.

Spaceborne stereoscopic systems have been growing in recent years, and the point cloud extracted from spaceborne stereo imagery has been used to measure forest spatial structures. These systems work on different viewing angles and image spatial resolutions, which are two critical factors determining the quality of the derived point cloud. In addition, the complex terrain is also a great challenge for the regional mapping of forest spatial structures using spaceborne stereo imagery. Although several theoretical models for simulating multi-view spectral features of forest canopies have been developed, there is hardly any report of a stereoscopic analysis using these models due to the limited size of the simulated forest scenes and the lack of a geometric sensory model (i.e., physical relationship between two-dimensional image coordinates and three-dimensional georeferenced coordinates). The stereoscopic features (i.e., parallax) are, as important as the spectral features contained in the multi-view images of a targeted area, the basis for the extraction of a point cloud. In this study, a new model, referred to as LandStereo model, has been proposed, which is capable of simulating the stereoscopic features of forest canopies over mountainous areas at landscape scales. The model comprised five parts, including defining the mountainous forest landscapes, setting the sun-senor observation geometry, simulating images, generating ground control points, and building geometric sensor models. The LandStereo model was validated over three different scenes, including flat forest landscapes, bare mountain landscapes, and mountainous forest landscapes. The results clearly demonstrated that the LandStereo model worked well on simulating stereoscopic features of both terrains and forest canopies at landscape scales. The extracted height of a forest canopy top from simulated stereo imagery was highly correlated to the truth (R2 = 0.96 and RMSE = 0.99 m) over the flat terrains and (R2 = 0.92 and RMSE = 1.15 m) over the mountainous areas. The LandStereo model provided a powerful tool to further our understanding of the relationships between forest spatial structures and point cloud extracted from stereo imagery acquired from different view angles and spatial resolutions under complex terrain conditions.

Applications of stereo imagery acquired by cameras onboard unmanned aerial vehicles (UAVs) as practical forest inventory tools are hindered by the unavailability of ground surface elevation. It is still a challenging issue to remove the elevation of ground surface in leaf-on stereo imagery to extract forest canopy height without the help of lidar data. This study proposed a method for the extraction of forest canopy height through the synthesis of UAV stereo imagery of leaf-on and leaf-off, and further demonstrated that the extracted forest canopy height could be used for the inventory of deciduous forest aboveground biomass (AGB). The points cloud of the leaf-on and leaf-off stereo imagery was firstly extracted by an algorithm of structure from motion (SFM) using the same ground control points (GCP). The digital surface model (DSM) was produced by rasterizing the point cloud of UAV leaf-on. The point cloud of UAV leaf-off was processed by iterative median filtering to remove vegetation points, and the digital terrain model (DTM) was generated by the rasterization of the filtered point cloud. The mean canopy height model (MCHM) was derived from the DSM subtracted by the DTM (i.e., DSM-DTM). Forest AGB maps were generated using models developed based on the MCHM and sampling plots of forest AGB and were evaluated by those of lidar. Results showed that forest AGB maps from UAV stereo imagery were highly correlated with those from lidar data with R2 higher than 0.94 and RMSE lower than 10.0 Mg/ha (i.e., relative RMSE 18.8%). These results demonstrated that UAV stereo imagery could be used as a practical inventory tool for deciduous forest AGB.

As the largest among terrestrial ecosystems, forests are vital to maintaining ecosystem services and regulating regional climate. The area and spatial distribution of trees in densely forested areas have been focused on in the past few decades, while sparse forests in agricultural zones, so-called agroforests or trees outside forests (TOF), have usually been ignored or missed in existing forest mapping efforts, despite their important role in regulating agricultural ecosystems. We combined Landsat and PALSAR data to map forests in a typical agricultural zone in the North China Plain. The resultant map, based on PALSAR and Landsat (PL) data, was also compared with five existing medium resolution (30–100 m) forest maps from PALSAR (JAXA forest map) and Landsat: NLCD-China, GlobeLand30, ChinaCover, and FROM-GLC. The results show that the PL-based forest map has the highest accuracy (overall accuracy of 95 ± 1% with a 95% confidence interval, and Kappa coefficient of 0.86) compared to those forest maps based on single Landsat or PALSAR data in the North China Plain (overall accuracy ranging from 85 ± 2% to 92 ± 1%). All forest maps revealed higher accuracy in densely forested mountainous areas, while the PL-based and JAXA forest maps showed higher accuracy in the plain, as the higher omission errors existed in only the Landsat-based forest maps. Moreover, we found that the PL-based forest map can capture more patched forest information in low forest density areas. This means that the radar data have advantages in capturing forests in the typical agricultural zones, which tend to be missing in published Landsat-based only forest maps. Given the significance of agroforests in regulating ecosystem services of the agricultural ecosystem and improving carbon stock estimation, this study implies that the integration of PALSAR and Landsat data can provide promising agroforest estimates in future forest inventory efforts, targeting a comprehensive understanding of ecosystem services of agroforests and a more accurate carbon budget inventory.

Estimating clumping indices is important for determining the leaf area index (LAI) of forest canopies. The spatial distribution of the clumping index is vital for LAI estimation. However, the neglect of woody tissue can result in biased clumping index estimates when indirectly deriving them from the gap probability and LAI observations. It is difficult to effectively and automatically extract woody tissue from digital hemispherical photos. In this study, a method for the automatic detection of trunks from Terrestrial Laser Scanning (TLS) data was used. Between-crown and within-crown gaps from TLS data were separated to calculate the clumping index. Subsequently, we analyzed the gap probability, clumping index, and LAI estimates based on TLS and HemiView data in consideration of woody tissue (trunks). Although the clumping index estimated from TLS had better agreement (R2 = 0.761) than that from HemiView, the change of angular distribution of the clumping index affected by the trunks from TLS data was more obvious than with the HemiView data. Finally, the exclusion of the trunks led to a reduction in the average LAI by ~19.6% and 8.9%, respectively, for the two methods. These results also showed that the detection of woody tissue was more helpful for the estimation of clumping index distribution. Moreover, the angular distribution of the clumping index is more important for the LAI estimate than the average clumping index value. We concluded that woody tissue should be detected for the clumping index estimate from TLS data, and 3D information could be used for estimating the angular distribution of the clumping index, which is essential for highly accurate LAI field measurements.

Mapping the magnitude and spatial distribution of forest aboveground biomass (AGB, in Mg·ha−1) is crucial to improve our understanding of the terrestrial carbon cycle. Landsat/TM (Thematic Mapper) and ICESat/GLAS (Ice, Cloud, and land Elevation Satellite, Geoscience Laser Altimeter System) data were integrated to estimate the AGB in the Changbai Mountain area. Firstly, four forest types were delineated according to TM data classification. Secondly, different models for prediction of the AGB at the GLAS footprint level were developed from GLAS waveform metrics and the AGB was derived from field observations using multiple stepwise regression. Lastly, GLAS-derived AGB, in combination with vegetation indices, leaf area index (LAI), canopy closure, and digital elevation model (DEM), were used to drive a data fusion model based on the random forest approach for extrapolating the GLAS footprint AGB to a continuous AGB map. The classification result showed that the Changbai Mountain region was characterized as forest-rich in altitudinal vegetation zones. The contribution of remote sensing variables in modeling the AGB was evaluated. Vegetation index metrics account for large amount of contribution in AGB ranges