Detecting Harvest Events in Plantation Forest Using Sentinel-1 and -2 Data via Google Earth Engine

: South Africa is reported to experience timber shortages as a result of growing timber demands and pulp production, coupled with the government’s reluctance to grant new forestry permits. Rampant timber theft in the country makes these circumstances worse. The emergence of cloud-based platforms, such as Google Earth Engine (GEE), has greatly improved the accessibility and usability of high spatial and temporal Sentinel-1 and -2 data, especially in data-poor countries that lack high-performance computing systems for forest monitoring. Here, we demonstrate the potential of these resources for forest harvest detection. The results showed that Sentinel-1 data are e ﬃ cient in detecting clear-cut events; both VH and VV backscatter signals decline sharply in accordance with clear-cutting and increase again when forest biomass increases. When correlated with highly responsive NDII, the VH and VV signals reached the best accuracies of 0.79 and 0.83, whereas the SWIR1 achieved –0.91. A Random Forest (RF) algorithm based on Sentinel-2 data also achieved over 90% accuracies for classifying harvested and forested areas. Overall, our study presents a cost-e ﬀ ective method for mapping clear-cut events in an economically important forestry area of South Africa while using GEE resources.


Introduction
Commercial forestry is a primary and widespread land-use activity along the eastern seaboard of South Africa, with nearly 1.2 million hectares being devoted to timber production [1]. The establishment of these plantations has been repeatedly hampered by disturbances, such as forest fires [2], the outbreak of insect pests and disease [3], invasive alien plants [4], and drought [5,6]. Additionally, forest land-use modification that arises from harvesting operations has become the main source of contemporary forest disturbance [7]. Because of the increase in the demand for timber and pulp products, combined with government's reluctance to grant new forestry permits, South Africa is currently facing timber shortages [8], which are, in part, further heightened by rampant timber theft [9]. Similar concerns about the effects of illegal logging and trade in timber prompted the European Commission in order to adopt the EU Action Plan for Forest Law Enforcement, Governance and Trade (FLEGT) in 2003-an effort that is intended to improve the legitimacy of the production and sale of timber [10]. This plan entails a variety of measures to timber producing countries, such as South Africa, which include promoting public procurement policies, safeguards for financing and investment, addressing the problem of conflict timber, and so on [11]. In this regard, the mapping of forest change provides insight into timber (6-12 days) when compared with previous SAR constellations [38]. Sentinel-1 has already displayed a high ability for monitoring forest applications, such as forest classification [39], forest cover mapping [40], burnt forested areas [41], forest logging [42], and deforestation [30]. In this paper, we further evaluate the utility of the weather-insensitive Sentinel-1 C-band, together with Sentinel-2 multispectral imagery, for monitoring forest harvesting events in KwaMbonambi, South Africa. Despite their potential for forest monitoring, only few studies have utilized these resources in order to detect clear-cut harvest events and, even are limited by exploring the GEE environment. Here, we demonstrate how to reconstruct high-density Sentinel-1 C-band time series in the process classify harvest events in planted forests while using Sentinel-2 within a cloud-based GEE platform and, consequently, stress the value of these resources for this purpose. We then compare the performance of Sentinel-1 with Landsat-8 data in detecting harvest events, as recommended by [43]. Finally, we illustrate how our results serve to identify harvest events in planted forests and how they could efficiently help monitor related disturbances, such as unauthorized harvesting and deforestation in different forest types.

Study Area
The study was performed in the KwaMbonambi plantation forests, which are located along the eastern seaboard of South Africa, 30 km northeast of Richards Bay ( Figure 1). They are dominated by even-aged Eucalyptus plantations of almost 40,000 ha, which are managed by Sappi, a large South African pulp and paper company. This area is representative of frequently harvested forests and it provides an ideal testing location. The growing stock mainly comprises 6-14-year-old Eucalyptus grandis (E. grandis), E. grandis × Eucalyptus camaldulensis (E. gxc), and E. grandis × Eucalyptus urophylla (E. gxu) hybrid clones. The stands are fairly uniform regarding canopy cover with tree density established at 1667 trees ha -1 . It is typical of the northern KwaZulu-Natal forestry region, which is subject to notable forest disturbances, such as fire and insect pest outbreaks. A subtropical climate characterizes the area, and the mean annual temperature is 22 • C [44]. The area receives an annual rainfall averaging 1200 mm, which is highly seasonal and peaks between November and February [45]; the potential evapotranspiration is 1772 mm per year [46]. The landscape of KwaMbonambi is flat and it consists of Quaternary alluvial sediments of clay sands of aeolian deposits [47] and soil with varying levels of organic matter [48], at mean elevation of 74 m above sea level. The high penetrability of the soils permits rapid leaching of soil nutrients due to the high rainfall in this region [47]. These conditions are favorable for fast-growing Eucalyptus plantations [49].
We used the vector polygons of known clear-cut compartments from Sappi, comprising felled date, areal extent, stand age classes, and forest types, in order to validate the clear-cut mapping results. These compartments were used to test whether high-density, readily available Sentinel-1 and -2 data can serve to detect harvest events with accuracies that are comparable with those acquired through forest inventories. Table 1 specifies the characteristics of the selected stands against which the competence of Sentinel-1 and -2 data were tested to detect clear-cut harvesting. We studied similar forest species (E. gxu) to compare, because Dostálová et al. [40] attributed differences in microwave backscatter variations to forest type and structure, as can be seen from the table.
where Ik(υ) is the radar intensity of output image k at pixel position υ, Ii(υ) is the radar intensity of the input image i, 〈 ( )〉 is the local average intensity of the input image i (window size of 7 × 7), and N is the number of images.  We also used Sentinel-2 Level-1C product representing Top of Atmosphere (TOA) reflectance for classification, again obtained and processed from the GEE environment. The spectrally sensitive Sentinel-2 is being increasingly utilized for various applications because of its high spatial (10-m) and temporal (3-5 days) resolutions.
We followed Nagler et al. [52] by using the Landsat-8-derived (L1 terrain corrected at TOA) normalized difference vegetation index (NDVI) and the normalized difference infrared index (NDII) values as references in order to evaluate the performance of Sentinel-1. In the region studied, Landsat data are impaired by data gaps due to persistent cloud cover. We then employed the spline interpolation to fit a cubic spline to each missing NDVI and NDII value over the time series, while

Satellite Data
All 110 available Sentinel-1, C-band SAR dense time series images, covering the KwaMbonambi area between June 2015 to December 2018, were used in order to detect clear-cut harvest events. The images were acquired and preprocessed by the GEE team while using the Sentinel-1 toolbox and available at no charge from the Google Earth Engine (GEE) cloud environment as image collection. The images were collected in Interferometric Wide Swath (IW, 250 km swath width) providing VV (Vertical transmit-Vertical receive) and VH (Vertical transmit-Horizontal receive) polarizations. This 10-m resolution, all-weather imagery with a 6-12-day revisit cycle was already processed at high-resolution Level-1 GRD level, which includes radiometric, geometric corrections and orthorectification [50]. We further applied the speckle filter using the JavaScript code editor in the GEE environment to reduce the speckle effect. As noted by Mermoz et al. [51], this filter produces images with reduced speckle effects, multi-temporal and multi-polarized (VH and VV) products, and are expressed, as follows: where I k (υ) is the radar intensity of output image k at pixel position υ, I i (υ) is the radar intensity of the input image i, I i (υ) is the local average intensity of the input image i (window size of 7 × 7), and N is the number of images. We also used Sentinel-2 Level-1C product representing Top of Atmosphere (TOA) reflectance for classification, again obtained and processed from the GEE environment. The spectrally sensitive Sentinel-2 is being increasingly utilized for various applications because of its high spatial (10-m) and temporal (3-5 days) resolutions.
We followed Nagler et al. [52] by using the Landsat-8-derived (L1 terrain corrected at TOA) normalized difference vegetation index (NDVI) and the normalized difference infrared index (NDII) values as references in order to evaluate the performance of Sentinel-1. In the region studied, Landsat data are impaired by data gaps due to persistent cloud cover. We then employed the spline interpolation to fit a cubic spline to each missing NDVI and NDII value over the time series, while using the R statistical package 'spline'. Cubic spline interpolation is widely used in order to generate a complete time series [53]. Because of different acquisition dates, both NDVI and NDII as well as Sentinel-1 profiles were averaged to represent monthly intervals, so that observations correspond with each other. Moreover, the use of SAR requires the removal of speckle and noise, which are efficiently reduced with monthly composites [41]. The vegetation indices were computed from top-of-atmosphere reflectance, using the following equations:

Background and the Use of the GEE Platform
The remote sensing of forests has reached the stage that enables scientists to expend their efforts on information generation, rather than data preparation [54]. This has been achieved through the emergence of a cloud-based GEE platform, a computing environment that provides easy access to a multi-petabyte ready-to-use satellite data archive [55]. It affords users superior computing power for processing and analyzing images without computational burden, and coupled with a high-performance, intrinsically parallel computation service. The GEE archive contains the complete Landsat data set as well as Sentinel-1 and Sentinel-2 images, and it also includes climate and other geophysical data sets [55]. The GEE can visualize the data, plot it in time series graphs, and the data can be downloaded for external processing.

GEE-Based Random Forest Classification
The GEE offers a variety of advanced machine learning classifiers for pixel-based classification that can be used for multi-temporal land use mapping [56]. Here, we implemented a Random Forest (RF; [57]) classifier for harvested and forested pixels' separation, and more information can be obtained from (https://developers.google.com/earth-engine/classification). RF has gained wide popularity, because of its high accuracy and ability to process complex datasets and produce satisfactory results with large numbers of input classification bands and training points [57]. We employed RF based on Sentinel-2 imagery, because it is ideal to classify very small plantations [58]. Certainly, some spectral bands are more important than others for the classification, and it has been established that, when more spectral bands are included, the accuracy is improved until a certain threshold is reached [59]. Thus, all of the spectral bands in the Sentinel-2 imageries were selected to train the classifier ( Table 2).
Reference data for training and validation comprised harvest records from Sappi management database, in the form of ArcMap shapefile polygons. The shapefiles contained fell dates, areal extent for each compartment, etc. Within GEE, Sentinel-2 TrueColor imagery corresponding to the fell dates were filtered and 60 points for each harvested and forested compartment were sampled throughout the studied area. Up to now, no literature stipulates the minimum number of training samples for machine learning algorithms, and because harvested compartments are clearly discernible from forested ones, we believed that 60 samples are ideal for this kind of classification. In order to produce robust classification results, the classification of Sentinel-2 images was run 1000 times by randomly selecting 50% of reference data from each class for training, and testing against the other 50% [58]. In classification, the accuracy assessment is the necessary step in evaluating the performance of algorithm used. We computed producer accuracy (PA), which is calculated by dividing the number of correctly classified pixels in each category by the number of "known" pixels to be of that category, the user's accuracy (UA), which is computed by dividing the number of correctly classified pixels in each category by the total number of pixels that were classified in that category, and the overall accuracy (OA) that is calculated by dividing the sum of correctly classified pixels by the total number of sampled pixels.

Statistical Analysis
A Taylor diagram was exploited to establish the relationship between spectral indices and the Sentinel-1 backscatter signals. This is a single diagram that synchronously combines three statistical performance metrics: the Pearson correlation coefficient, root mean square difference (RMSD), and the standard deviation in a series of points on a polar plot [60]. The RMSD is a measure of the average error that is produced by a model and it is widely used to gauge overall model performance. However, it does not reveal the types or sources of the error, which would assist greatly in refining the models [61]. The RMSD is calculated as: Figure 2 illustrates the comparisons between VH and VV backscatter signals and the field inventory record. It can be seen that both VH and VV bands detect the clear-cut events over the same area and time period. However, the VH appears fairly constant around −13 dB and −15 dB in uncut stands with a clear decline in backscatter signals over clear-cut stands reaching −21 dB. On the other hand, the VV bands seem to be more sensitive to various forest changes, as illustrated by the spread of weaker backscattered signals that were observed across the plantation, even over uncut stands. Despite this observation, the patterns of clear-cut stands are distinguishable from the VV backscatter, because they are dominant over clear-cut areas. Xu et al. [62] also noted that the VV band suffers from the complex effects of changes in the scattering mechanism. The variation in VV backscatter in 2016 is believed to be associated with the exceptional drought conditions reported by Xulu et al. [6]. Chatziantoniou et al. [63] noted that extreme drought affects the VV signal, because it reduces backscatter, which is sensitive to moisture conditions. The catastrophic character of the 2015-2016 drought led to widespread tree dieback in the Zululand forestry region [64]. This outcome appears to be illustrated by the extensive reduced VV backscatter over the uncut forest stands, as demonstrated in Figure 2. to be illustrated by the extensive reduced VV backscatter over the uncut forest stands, as demonstrated in Figure 2. The classification of harvested and forested compartments while using RF based on Sentinel-2 imagery successfully produced overall accuracies as high as 99% for all sampled years ( Figure 2). The producer's and user's accuracies for both classes were consistently high, with each class attaining over 90%. Our results are consistent with Nomura and Mitchard [58], in their classification of a plantations in a complex forested landscape. While both of the Sentinel-1 bands seem to be reactive to clear-cutting, and despite the small differences between them, the VV bands are seen to be slightly sensitive than the VH band. This observation is also consistent with greater correlation of VV band to highly responsive NDII (shown in the next sections).

Spatial Patterns of VH and VV Backscatter over Clear-Cut Stands and the Classification
We performed a time series analysis of VH and VV polarization indices and Landsat-derived NDVI and NDII (as reference) for each year separately over known harvested stands in the KwaMbonambi plantations in order to evaluate the effectiveness of Sentinel-1 backscattering to characterize clear-cut events. These sensors record different, yet related, forest attributes. In particular, VH and VV radar backscattering is highly correlated with fresh forest biomass, whereas NDVI is correlated with photosynthetic activity [31]. Given that clear-cutting is an extreme and The classification of harvested and forested compartments while using RF based on Sentinel-2 imagery successfully produced overall accuracies as high as 99% for all sampled years ( Figure 2). The producer's and user's accuracies for both classes were consistently high, with each class attaining over 90%. Our results are consistent with Nomura and Mitchard [58], in their classification of a plantations in a complex forested landscape.
While both of the Sentinel-1 bands seem to be reactive to clear-cutting, and despite the small differences between them, the VV bands are seen to be slightly sensitive than the VH band. This observation is also consistent with greater correlation of VV band to highly responsive NDII (shown in the next sections).
We performed a time series analysis of VH and VV polarization indices and Landsat-derived NDVI and NDII (as reference) for each year separately over known harvested stands in the KwaMbonambi plantations in order to evaluate the effectiveness of Sentinel-1 backscattering to characterize clear-cut events. These sensors record different, yet related, forest attributes. In particular, VH and VV radar backscattering is highly correlated with fresh forest biomass, whereas NDVI is correlated with photosynthetic activity [31]. Given that clear-cutting is an extreme and abrupt forest change, and that sensors may yield signals with variations as a response to other forest dynamics, this study only considers the most significant breakpoints. We decided to only select three harvested compartments for each year to exemplify the behavior of Sentinel-1 and spectral bands when harvest events transpire, and the analysis is presented below. Figure 3 shows the temporal profiles of Sentinel-1 VH and VV backscatter (a), NDVI and NDII (b), and the shortwave infrared (SWIR1) (c) for the compartment harvested in 2016. The VH and VV polarized backscatter exhibited similar results for this compartment, generally with a virtually uniform pattern throughout the study period with the exception of a notable change in the middle of 2016. An abrupt reduction of 5 dB for VH (−13.7 to −18 dB) and 4 dB for VV (−8.4 to −11.6 dB) backscattering in June 2016 corresponds with a clear-cut harvest event, as confirmed by Sappi's official inventory records. During this period, it appears, as expected from cleared forest, the VH polarization revealed much less backscattering than the VV signal. This observation is consistent with Nguyen et al. [65], who reported that VH polarized backscatter is generally less than VV, and that the backscattering coefficients for both polarizations progressively increase with tree growth over time. There was an obvious increase in the backscatter signal for both polarizations, corresponding with tree canopy increase after September 2016, reaching a peak of −8 dB (VV) and −13 dB (VH) around March 2017, when the forests reached maturity stage, as illustrated in Figure 3a. An equivalent finding was made by Macelloni et al. [66], who observed an increase in VH cross-polarized backscatter with increasing leaf area index (LAI) over planted forests in the Italian landscape. In a related study, Huang et al. [67] found that VH polarization is more efficient than VV for characterizing forest canopy dynamics, because it is less influenced by soil moisture. Nguyen et al. [66] attributed the reduced backscatter in the VV signal to the lower moisture content of trees during the harvest period; conversely, the VH signal seemed to be the least affected by this factor. Moreover, we agree with Karjalainen et al. [68] that VH polarization is better than VV to characterize the actual biomass variations in forests.

Temporal Profile over Harvested Compartment in 2016
The NDII and NDVI values resembled the Sentinel-1 backscatter signals and declined during clear-cutting (Figure 3b). The NDII declined from 1.7 to −0.13 and recovered to 0.39 a few months later. The NDVI also decreased, from 0.54 to 0.32, after a clear-cut event and rose to 0.77 during the forest growth phase.
It is important to note the increase of the SWIR1 signal during clear-cutting. At this time, the SWIR1 rose from 0.15 to its highest maximum of 0.28 (Figure 3c). Indeed, clear-cutting is well-known to increase the reflectance in the SWIR1 band [69]. This is particularly so, because the SWIR1 band is sensitive to vegetation density [70]. Thus, as forest canopy complexity increases with time, reflectance in the SWIR1 decreases [71]. Our results are in harmony with this observation-the SWIR1 signal declines as the forest canopy matures. In a related study, Schroeder et al. [72] classified forest disturbances while using Landsat-derived vegetation indices and found the SWIR1 channel to be the most effective for distinguishing clear-cut harvests.
Forests 2020, 11, x FOR PEER REVIEW 9 of 17 as shown in Figure 4c. An initial increase in SWIR1 from 0.11 to 0.36 in April was followed by a further increase to 0.37 during a clear-cutting period.   (Figure 4b). The increase in SWIR1 corresponds with reductions in the SAR backscatter signals and vegetation indices, as shown in Figure 4c. An initial increase in SWIR1 from 0.11 to 0.36 in April was followed by a further increase to 0.37 during a clear-cutting period.

Temporal Profile over Harvested Compartment in 2018
The pattern of VH and VV backscatter signals for the compartment that was harvested in 2018 was similar to that for 2016 and 2017 harvested compartments. The VH backscatter dropped from −13.81 dB in January 2018 to its lowest minimum (−19.82 dB) in July; the VV backscatter also declined, from −8.33 to −12.54 dB (Figure 5a). The NDVI value decreased from 0.72 to 0.21, whereas the NDII declined from 0.48 to −0.19 in the corresponding periods ( Figure 5b). As expected, the SWIR1 increased from 0.11 to 0.34 over the time of the clear-cutting (Figure 5c).

Temporal Profile over Harvested Compartment in 2018
The pattern of VH and VV backscatter signals for the compartment that was harvested in 2018 was similar to that for 2016 and 2017 harvested compartments. The VH backscatter dropped from −13.81 dB in January 2018 to its lowest minimum (−19.82 dB) in July; the VV backscatter also declined, from −8.33 to −12.54 dB (Figure 5a). The NDVI value decreased from 0.72 to 0.21, whereas the NDII declined from 0.48 to −0.19 in the corresponding periods ( Figure 5b). As expected, the SWIR1 increased from 0.11 to 0.34 over the time of the clear-cutting (Figure 5c).   Figure 5 seems to indicate that Sentinel-1 can be an efficient weatherinsensitive source from which clear-cutting activities are detectable. It is noteworthy that the Sentinel-1 backscatter appear to be responsive to various disturbances, but the dominant backscatter is obvious during clear-cutting. Similarly, the vegetation indices (NDVI and NDII) and the SWIR1 band showed a corresponding pattern to the Sentinel-1 backscatter signals, and these were useful in confirming the results. We now explain the relationship of these signals having described the temporal profiles of data derived from Sentinel-1 and Landsat-8.

Correlation between NDII Sentinel-1 and Other Vegetation Indices
Given the greater sensitivity of NDII to clear-cutting events, we decided to use it as a reference variable in order to illustrate its relationship with microwave indices (VH and VV) as well as other vegetation criteria. For this purpose, Taylor diagrams [60] were computed in order to provide a graphical representation of how closely a pattern corresponds to the reference signal. The Taylor Figures 3-5 seems to indicate that Sentinel-1 can be an efficient weather-insensitive source from which clear-cutting activities are detectable. It is noteworthy that the Sentinel-1 backscatter appear to be responsive to various disturbances, but the dominant backscatter is obvious during clear-cutting. Similarly, the vegetation indices (NDVI and NDII) and the SWIR1 band showed a corresponding pattern to the Sentinel-1 backscatter signals, and these were useful in confirming the results. We now explain the relationship of these signals having described the temporal profiles of data derived from Sentinel-1 and Landsat-8.

Correlation between NDII Sentinel-1 and Other Vegetation Indices
Given the greater sensitivity of NDII to clear-cutting events, we decided to use it as a reference variable in order to illustrate its relationship with microwave indices (VH and VV) as well as other vegetation criteria. For this purpose, Taylor diagrams [60] were computed in order to provide a graphical representation of how closely a pattern corresponds to the reference signal. The Taylor diagram is explanatory diagram that presents the visualization of the comparative strength of independent parameters to the actual target variable. In the Taylor diagram, the two different statistical metrics (i.e., correlation coefficients (R 2 ) and standard deviations of each model) are used to quantify the comparability between the calculated data (e.g., models) and actual data. The distance from the reference point is a measure of the centered RMSD. The Taylor diagram is designed in such a way that it can graphically indicate which calculated data (or models) is most realistic. Figure 6 shows the corresponding Taylor diagram for the three different compartments, harvested in 2016, 2017, and 2018. The VV polarization was better correlated with NDII (compartment harvested in 2016, R 2 = 0.78; compartment harvested in 2017, R 2 = 0.81; compartment harvested in 2018, R 2 = 0.82) than VH backscatter (2016 harvested compartment, R 2 = 0.67; 2017 harvested compartment, R 2 = 0.72; 2018 harvested compartment, R 2 = 0.79). The values that correspond to these polarizations were grouped around the normalized standard deviation between 1.25 and 1.81, with VV showing less deviation from NDII than VH backscatter. The VV backscatter also had a smaller RMSD value (~1) than the VH signal (~2). Similar to those of Nagler et al. [52], these results reveal the good quality of both products that are derived from independent data sources. diagram is explanatory diagram that presents the visualization of the comparative strength of independent parameters to the actual target variable. In the Taylor diagram, the two different statistical metrics (i.e., correlation coefficients (R 2 ) and standard deviations of each model) are used to quantify the comparability between the calculated data (e.g., models) and actual data. The distance from the reference point is a measure of the centered RMSD. The Taylor diagram is designed in such a way that it can graphically indicate which calculated data (or models) is most realistic. Figure 6 shows the corresponding Taylor [52], these results reveal the good quality of both products that are derived from independent data sources.  As expected, the NDVI and SWIR1 results were more highly correlated with NDII than the Sentinel-1 bands; this was because the former were partly derived from similar bands to those of It is also important to note that the increasing water deficits increase the spectral reflectance in the SWIR1 region, and this largely influences the change in the vegetation values of indices that exploits this band, such as NDII [73]. The greater sensitivity of Sentinel-1 to moisture variations, particularly the VV band, could explain its higher correlation with the NDII.

Conclusions
Our study demonstrated the impressive capability of Sentinel-1 and -2 in order to detect and map clear-cut events in a complex commercial forestry area. The results showed high RF classification accuracies (99%) for harvested compartments based on Sentinel-2 data, which was consistent with Sentinel-1 backscatter. Our results also exhibited a similar temporal profile of VH and VV backscattering indices to those of NDVI and NDII; the respective signals were consistently relatively strong during the mature forest stage, but abruptly decreased after clear-cutting. In contrast, the SWIR1 band exhibited the opposite pattern with a sharp increase coinciding with clear-cutting, confirming reduced leaf water content prevailing at this time. When correlated with the highly responsive NDII, the VH and VV signals reached the best accuracies of 0.79 and 0.83, whereas the NDVI and SWIR1 achieved 0.96 and −0.91, respectively. The agreement between the VH and VV indices with NDII is satisfactory, because they are all derived from sensors that record different vegetation properties when compared to NDVI and SWIR1, which were not only recorded by the same sensor, but also implicitly contain NDII bands. Nevertheless, the VH and VV bands both consistently detected clear-cut events, the VH record indicated stable low backscattering over clear-cut stands, while the VV band appeared more sensitive, but it is affected by moisture stress and led to weaker backscattering signals even over uncut forest areas. Further research is desirable in exploring this method over broader areas and for different forest types. The responsiveness of the VV to moisture variations could also be tested in drought affected stands. We hope that our results will encourage remote sensing-based forest researchers in order to exploit the freely available GEE resources further and develop methodologies, in support of sustainable forest management, using high-density Sentinel-1 data to characterize forest changes, particularly in regions with poor image availability due to persistent cloud cover.
Author Contributions: Conceptualization, S.X., N.M., K.P. and M.G.; methodology, formal analysis, resources, writing of original draft, S.X., N.M., K.P. and M.G., and review and editing of final manuscript: S.X., N.M., K.P. and M.G. All authors have read and agreed to the published version of the manuscript.