Here, we examined multitemporal behavior of fully polarimetric SAR (PolSAR) parameters at L-band in relation to the stem volume of boreal forests. The PolSAR parameters were evaluated in terms of their temporal consistency, inter-dependence and suitability for forest stem volume estimation across several seasonal conditions (frozen, thaw and unfrozen). The satellite SAR data were represented by a time series of PolSAR images acquired during several seasons in the years 2006 to 2009 by the ALOS PALSAR sensor. The study area was in central Finland, and represented a managed area in typical boreal mixed forest land. Utility of different PolSAR parameters, their temporal stability and cross-correlations were studied along with reference stand-level stem volume data from forest inventory. Further, two polarimetric parameters, cross-polarization backscatter and co-polarization coherence, were chosen for further investigation and stem volume retrieval. A relationship between forest stem volume and PolSAR parameters was established using the kNN regression approach. Ways of optimally combining PolSAR images were evaluated as well. For a single scene, best results were observed with polarimetric coherence (RMSE ≈ 38.8 m³/ha) for scene acquired in frozen conditions. An RMSE of 40.8 m³/ha (42.9%, R² = 0.66) was achieved for cross-polarization backscatter in the best case. Cross-polarization backscatter was a better predictor than polarimetric coherence for few summer scenes. Multitemporal aggregation of selected PolSAR scenes improved estimates for both studied PolSAR parameters. Stronger improvement was observed for coherence with RMSE down to 34 m³/ha (35.8%, R² = 0.77) compared to 38.8–51.6 m³/ha (40.8–54.3%) from separate scenes. Finally, the accuracy statistics reached RMSE of 32.2 m³/ha (34%, R² = 0.79) when multitemporal HHVV coherence was combined with multitemporal HV-backscatter.

The applicability of optical and synthetic aperture radar (SAR) data for land cover classification to support REDD+ (Reducing Emissions from Deforestation and Forest Degradation) MRV (measuring, reporting and verification) services was tested on a tropical to sub-tropical test site. The 100 km by 100 km test site was situated in the State of Chiapas in Mexico. Land cover classifications were computed using RapidEye and Landsat TM optical satellite images and ALOS PALSAR L-band and Envisat ASAR C-band images. Identical sample plot data from Kompsat-2 imagery of one-metre spatial resolution were used for the accuracy assessment. The overall accuracy for forest and non-forest classification varied between 95% for the RapidEye classification and 74% for the Envisat ASAR classification. For more detailed land cover classification, the accuracies varied between 89% and 70%, respectively. A combination of Landsat TM and ALOS PALSAR data sets provided only 1% improvement in the overall accuracy. The biases were small in most classifications, varying from practically zero for the Landsat TM based classification to a 7% overestimation of forest area in the Envisat ASAR classification. Considering the pros and cons of the data types, we recommend optical data of 10 m spatial resolution as the primary data source for REDD MRV purposes. The results with L-band SAR data were nearly as accurate as the optical data but considering the present maturity of the imaging systems and image analysis methods, the L-band SAR is recommended as a secondary data source. The C-band SAR clearly has poorer potential than the L-band but it is applicable in stratification for a statistical sampling when other image types are unavailable.

Due to the high cost of traditional forest plot measurements, the availability of up-to-date in situ forest inventory data has been a bottleneck for remote sensing image analysis in support of the important global forest biomass mapping. Capitalizing on the proliferation of smartphones, citizen science is a promising approach to increase spatial and temporal coverages of in situ forest observations in a cost-effective way. Digital cameras can be used as a relascope device to measure basal area, a forest density variable that is closely related to biomass. In this paper, we present the Relasphone mobile application with extensive accuracy assessment in two mixed forest sites from different biomes. Basal area measurements in Finland (boreal zone) were in good agreement with reference forest inventory plot data on pine ( $R^2=0.75$ , $RMSE=5.33$ m ${}^2$ /ha), spruce ( $R^2=0.75,RMSE=6.73$ m ${}^2$ /ha) and birch ( $R^2=0.71,RMSE=4.98$ m ${}^2$ /ha), with total relative $RMSE(\%)=29.66\%$ . In Durango, Mexico (temperate zone), Relasphone stem volume measurements were best for pine ( $R^2=0.88$ , $RMSE=32.46$ m ${}^3$ /ha) and total stem volume ( $R^2=0.87$ , $RMSE=35.21$ m ${}^3$ /ha). Relasphone data were then successfully utilized as the only reference data in combination with optical satellite images to produce biomass maps. The Relasphone concept has been validated for future use by citizens in other locations.

Deep learning (DL) models are gaining popularity in forest variable prediction using Earth observation (EO) images. However, in practical forest inventories, reference datasets are often represented by plot- or stand-level measurements, while high-quality representative wall-to-wall reference data for end-to-end training of DL models are rarely available. Transfer learning facilitates expansion of the use of deep learning models into areas with sub-optimal training data by allowing pretraining of the model in areas where high-quality teaching data are available. In this study, we perform a “model transfer” (or domain adaptation) of a pretrained DL model into a target area using plot-level measurements and compare performance versus other machine learning models. We use an earlier developed UNet based model (SeUNet) to demonstrate the approach on two distinct taiga sites with varying forest structure and composition. The examined SeUNet model uses multi-source EO data to predict forest height. Here, EO data are represented by a combination of Copernicus Sentinel-1 C-band SAR and Sentinel-2 multispectral images, ALOS-2 PALSAR-2 SAR mosaics and TanDEM-X bistatic interferometric radar data. The training study site is located in Finnish Lapland, while the target site is located in Southern Finland. By leveraging transfer learning, the SeUNet prediction achieved root mean squared error (RMSE) of $2.70$ m and R² of 0.882, considerably more accurate than traditional benchmark methods. We expect such forest-specific DL model transfer can be suitable also for other forest variables and other EO data sources that are sensitive to forest structure.

A method to detect abrupt land cover changes using hierarchical clustering of multi-temporal satellite imagery was developed. The Autochange method outputs the pre-change land cover class, the change magnitude, and the change type. Pre-change land cover information is transferred to post-change imagery based on classes derived by unsupervised clustering, enabling using data from different instruments for pre- and post-change. The change magnitude and change types are computed by unsupervised clustering of the post-change image within each cluster, and by comparing the mean intensity values of the lower level clusters with their parent cluster means. A computational approach to determine the change magnitude threshold for the abrupt change was developed. The method was demonstrated with three summer image pairs Sentinel-2/Sentinel-2, Landsat 8/Sentinel-2, and Sentinel-2/ALOS 2 PALSAR in a study area of 12,372 km² in southern Finland for the detection of forest clear cuts and tested with independent data. The Sentinel-2 classification produced an omission error of 5.6% for the cut class and 0.4% for the uncut class. Commission errors were 4.9% for the cut class and 0.4% for the uncut class. For the Landsat 8/Sentinel-2 classifications the equivalent figures were 20.8%, 0.2%, 3.4%, and 1.6% and for the Sentinel-2/ALOS PALSAR classification 16.7%, 1.4%, 17.8%, and 1.3%, respectively. The Autochange algorithm and its software implementation was considered applicable for the mapping of abrupt land cover changes using multi-temporal satellite data. It allowed mixing of images even from the optical and synthetic aperture radar (SAR) sensors in the same change analysis.

Dense time series of stripmap RADARSAT-2 data acquired in the Multilook Fine mode were used for detecting and mapping the extent of selective logging operations in the tropical forest area in the northern part of the Republic of the Congo. Due to limited radiometric sensitivity to forest biomass variation at C-band, basic multitemporal change detection approach was supplemented by spatial texture analysis to separate disturbed forest from intact. The developed technique primarily uses multi-temporal aggregation of orthorectified synthetic aperture radar (SAR) imagery that are acquired before and after the logging operations. The actual change analysis is based on textural features of the log-ratio image calculated using two SAR temporal composites compiled of SAR scenes acquired before and after the logging operations. Multitemporal aggregation and filtering of SAR scenes decreased speckle and made the extracted textural features more prominent. The overall detection accuracy was around 80%, with some underestimation of the area of forest disturbance compared to reference based on optical data. The user’s accuracy for disturbed forest varied from 76.7% to 94.9% depending on the accuracy assessment approach. We conclude that change detection utilizing RADARSAT-2 time series represents a useful instrument to locate areas of selective logging in tropical forests.