Anthropogenic stress and disturbance of forest ecosystems (FES) has been increasing at all scales from local to global. In rapidly changing environments, in-situ terrestrial FES monitoring approaches ...have made tremendous progress but they are intensive and often integrate subjective indicators for forest health (FH). Remote sensing (RS) bridges the gaps of these limitations, by monitoring indicators of FH on different spatio-temporal scales, and in a cost-effective, rapid, repetitive and objective manner. In this paper, we provide an overview of the definitions of FH, discussing the drivers, processes, stress and adaptation mechanisms of forest plants, and how we can observe FH with RS. We introduce the concept of spectral traits (ST) and spectral trait variations (STV) in the context of FH monitoring and discuss the prospects, limitations and constraints. Stress, disturbances and resource limitations can cause changes in FES taxonomic, structural and functional diversity; we provide examples how the ST/STV approach can be used for monitoring these FES characteristics. We show that RS based assessments of FH indicators using the ST/STV approach is a competent, affordable, repetitive and objective technique for monitoring. Even though the possibilities for observing the taxonomic diversity of animal species is limited with RS, the taxonomy of forest tree species can be recorded with RS, even though its accuracy is subject to certain constraints. RS has proved successful for monitoring the impacts from stress on structural and functional diversity. In particular, it has proven to be very suitable for recording the short-term dynamics of stress on FH, which cannot be cost-effectively recorded using in-situ methods. This paper gives an overview of the ST/STV approach, whereas the second paper of this series concentrates on discussing in-situ terrestrial monitoring, in-situ RS approaches and RS sensors and techniques for measuring ST/STV for FH.
Agricultural vegetation development and harvest date monitoring over large areas requires frequent remote sensing observations. In regions with persistent cloud coverage during the vegetation season ...this is only feasible with active systems, such as SAR, and is limited for optical data. To date, optical remote sensing vegetation indices are more frequently used to monitor agricultural vegetation status because they are easily processed, and the characteristics are widely known. This study evaluated the correlations of three Sentinel-2 optical indices with Sentinel-1 SAR indices over agricultural areas to gain knowledge about their relationship. We compared Sentinel-2 Normalized Difference Vegetation Index, Normalized Difference Water Index, and Plant Senescence Radiation Index with Sentinel-1 SAR VV and VH backscatter, VH/VV ratio, and Sentinel-1 Radar Vegetation Index. The study was conducted on 22 test sites covering approximately 35,000 ha of four different main European agricultural land use types, namely grassland, maize, spring barley, and winter wheat, in Lower Saxony, Germany, in 2018. We investigated the relationship between Sentinel-1 and Sentinel-2 indices for each land use type considering three phenophases (growing, green, senescence). The strength of the correlations of optical and SAR indices differed among land use type and phenophase. There was no generic correlation between optical and SAR indices in our study. However, when the data were split by land use types and phenophases, the correlations increased remarkably. Overall, the highest correlations were found for the Radar Vegetation Index and VH backscatter. Correlations for grassland were lower than for the other land use types. Adding auxiliary data to a multiple linear regression analysis revealed that, in addition to land use type and phenophase information, the lower quartile and median SAR values per field, and a spatial variable, improved the models. Other auxiliary data retrieved from a digital elevation model, Sentinel-1 orbit direction, soil type information, and other SAR values had minor impacts on the model performance. In conclusion, despite the different nature of the signal generation, there were distinct relationships between optical and SAR indices which were independent of environmental variables but could be stratified by land use type and phenophase. These relationships showed similar patterns across different test sites. However, a regional clustering of landscapes would significantly improve the relationships.
Stress in forest ecosystems (FES) occurs as a result of land-use intensification, disturbances, resource limitations or unsustainable management, causing changes in forest health (FH) at various ...scales from the local to the global scale. Reactions to such stress depend on the phylogeny of forest species or communities and the characteristics of their impacting drivers and processes. There are many approaches to monitor indicators of FH using in-situ forest inventory and experimental studies, but they are generally limited to sample points or small areas, as well as being time- and labour-intensive. Long-term monitoring based on forest inventories provides valuable information about changes and trends of FH. However, abrupt short-term changes cannot sufficiently be assessed through in-situ forest inventories as they usually have repetition periods of multiple years. Furthermore, numerous FH indicators monitored in in-situ surveys are based on expert judgement. Remote sensing (RS) technologies offer means to monitor FH indicators in an effective, repetitive and comparative way. This paper reviews techniques that are currently used for monitoring, including close-range RS, airborne and satellite approaches. The implementation of optical, RADAR and LiDAR RS-techniques to assess spectral traits/spectral trait variations (ST/STV) is described in detail. We found that ST/STV can be used to record indicators of FH based on RS. Therefore, the ST/STV approach provides a framework to develop a standardized monitoring concept for FH indicators using RS techniques that is applicable to future monitoring programs. It is only through linking in-situ and RS approaches that we will be able to improve our understanding of the relationship between stressors, and the associated spectral responses in order to develop robust FH indicators.
Local and landscape-scale agricultural intensification is a major driver of global biodiversity loss. Controversially discussed solutions include wildlife-friendly farming or combining high-intensity ...farming with land-sparing for nature. Here, we integrate biodiversity and crop productivity data for smallholder cacao in Indonesia to exemplify for tropical agroforests that there is little relationship between yield and biodiversity under current management, opening substantial opportunities for wildlife-friendly management. Species richness of trees, fungi, invertebrates, and vertebrates did not decrease with yield. Moderate shade, adequate labor, and input level can be combined with a complex habitat structure to provide high biodiversity as well as high yields. Although livelihood impacts are held up as a major obstacle for wildlife-friendly farming in the tropics, our results suggest that in some situations, agroforests can be designed to optimize both biodiversity and crop production benefits without adding pressure to convert natural habitat to farmland.
Vegetation canopy height is a relevant proxy for aboveground biomass, carbon stock, and biodiversity. Wall‐to‐wall information of canopy height with high spatial resolution and accuracy is not yet ...available on large scales. For the globally consistent TanDEM‐X data, simplifications are necessary to estimate canopy height with semi‐empirical models based on polarimetric synthetic aperture radar interferometry (PolInSAR).
We trained the semi‐empirical models with sampled GEDI data, because the assumptions behind the application of such simplifications are not always valid for TanDEM‐X. General linear as well as sinc models and empirical parameterizations of these models were applied to estimate the canopy height in tropical landscapes of Sumatra, Indonesia. Airborne laser scanning (ALS) data were consistently used as an independent reference. The general simplified models were compared with the trained empirical versions to assess the potential improvement of the empirical parameterization of the models. The residuals of the different canopy height models were further evaluated in relation to land use and structural information of the vegetation.
Our results indicated that the empirical parameters substantially improved the estimation from a root‐mean‐square‐error (RMSE) of 10.3 m (55.8%) to 8.8 m (47.7%), when using the linear model. In contrast, the improvement of the sinc model with empirical parameters was not substantial compared to the general sinc model (7.4 m 40.4% vs. 6.9 m 37.5%). A consistent improvement was observed in the linear model, whereas the improvement of the sinc model was dependent on the land‐use type. Structural attributes like the canopy height itself and vegetation cover had a significant effect on the accuracies, with higher and denser vegetation generally resulting in higher residuals.
We demonstrate the potential of the combined exploitation of the TanDEM‐X and GEDI missions for a wall‐to‐wall canopy height estimation in a tropical region. This study provides relevant findings for a consistent mapping of vegetation canopy height in tropical landscapes and on large scales with spaceborne laser and SAR data.
The potential of time series metrics based on optical sensors for monitoring of crops and their phenological stages has long been recognized. SAR (synthetic aperture radar) sensors such as Sentinel-1 ...can provide dense time series, but their potential for retrieving time series metrics is to date underexplored. The backscatter coefficient in VH (vertical-horizontal) and VV (vertical-vertical) polarization as well as the interferometric coherence of Sentinel-1 was retrieved in a German agricultural area in order to obtain a time series with a temporal resolution of 6 days. The backscatter coefficient ratio of VH and VV polarization was also calculated. Local extrema and breakpoints were detected in the locally smoothed time series data in wheat fields for the years 2017, 2018 and 2019. Breakpoints in the VH/VV ratio time series have high potential for detecting the date of shooting and harvesting of the studied wheat fields. The root mean squared differences between detected date and in situ observed date were 4.6 days in 2017, 5.3 days in 2018 and 9.5 days in 2019 for the shooting and 5.1 days in 2017, 8.2 days in 2018 and 10.4 days in 2019 for the harvest. Ripeness stages could not be detected with a high degree of accuracy. The results suggest that the VH/VV ratio has a high potential for monitoring the phenological stages of wheat fields, which can be used for large-scale crop productivity monitoring and risk assessment in agricultural cropping systems.
•We analyzed time series data from Sentinel-1 to monitor the phenology of wheat.•Breakpoints and local extrema of the time series were used as metrics.•Breakpoints in the VH/VV backscatter ratio time series achieved best results.•Dates of shooting and harvest were detectable with high accuracy.
•Crop types are mapped at high accuracy on a regional scale for 12 consecutive years.•Temporal encoding (TE) of annual time series into 365 equidistant features is introduced.•TE with two data ...augmentation techniques improves temporal transfer of deep learning models.•A single-year model trained on Landsat and Sentinel-2 data from 2018 is transferable to 2010 – 2021.•The single-year model performs consistently well with varying sensors constellations.
Detailed maps on the spatial and temporal distribution of crops are key for a better understanding of agricultural practices and for food security management. Multi-temporal remote sensing data and deep learning (DL) have been extensively studied for deriving accurate crop maps. However, strategies to solve the problem of transferring crop classification models over time, e.g., training the model with data for a recent year and mapping back to the past, have not been fully explored. This is due to the lack of a generalized method for aggregating optical data with regard to the irregularity in annual clear sky observations and the scarcity of multi-annual crop reference data to support a more generalized DL model. In this study, we tackled these challenges by introducing a method namely Temporal Encoding (TE) to capture the irregular phenological information. Subsequently, we adapted and integrated two methods, i.e., Random Observations Selection (ROS) and Random Day Shifting (RDS) to simulate the variability of temporal sparsity as well as the shifts of crop phenology over different years. We tested this approach with a 1-dimensional Convolutional Neural Network (1D-CNN) and a Transformer Network models. Our results for both classifiers showed that models trained with crop reference data from 2018 and a dense time series of Landsat 7/8 and Sentinel-2 A/B data can be transferred with little decreases in accuracy to map 12 consecutive years from 2010 to 2021. The Transformer Network was slightly more accurate, while the 1D-CNN was much three times faster. Furthermore, the proposed models could achieve similar performances in the same years with and without fully available satellite information. The TE with ROS and RDS appears well suited for improving temporal transferability to support long term historic crop mapping.
Many Indonesian forests have been cleared and replaced by fast-growing cash crops (e.g., oil palm and rubber plantations), altering the vegetation structure of entire regions. Complex vegetation ...structure provides habitat niches to a large number of native species. Airborne laser scanning (ALS) can provide detailed three-dimensional information on vegetation structure. Here, we investigate the potential of ALS metrics to highlight differences across a gradient of land-use management intensities in Sumatra, Indonesia. We focused on tropical rainforests, jungle rubber, rubber plantations, oil palm plantations and transitional lands. Twenty-two ALS metrics were extracted from 183 plots. Analysis included a principal component analysis (PCA), analysis of variance (ANOVAs) and random forest (RF) characterization of the land use/land cover (LULC). Results from the PCA indicated that a greater number of canopy gaps are associated with oil palm plantations, while a taller stand height and higher vegetation structural metrics were linked with rainforest and jungle rubber. A clear separation in metrics performance between forest (including rainforest and jungle rubber) and oil palm was evident from the metrics pairwise comparison, with rubber plantations and transitional land behaving similar to forests (rainforest and jungle rubber) and oil palm plantations, according to different metrics. Lastly, two RF models were carried out: one using all five land uses (5LU), and one using four, merging jungle rubber with rainforest (4LU). The 5LU model resulted in a lower overall accuracy (51.1%) due to mismatches between jungle rubber and forest, while the 4LU model resulted in a higher accuracy (72.2%). Our results show the potential of ALS metrics to characterize different LULCs, which can be used to track changes in land use and their effect on ecosystem functioning, biodiversity and climate.
Losses of biodiversity and ecosystem functioning due to rainforest destruction and agricultural intensification are prime concerns for science and society alike. Potentially, ecosystems show ...nonlinear responses to land-use intensification that would open management options with limited ecological losses but satisfying economic gains. However, multidisciplinary studies to quantify ecological losses and socioeconomic tradeoffs under different management options are rare. Here, we evaluate opposing land use strategies in cacao agroforestry in Sulawesi, Indonesia, by using data on species richness of nine plant and animal taxa, six related ecosystem functions, and on socioeconomic drivers of agroforestry expansion. Expansion of cacao cultivation by 230% in the last two decades was triggered not only by economic market mechanisms, but also by rarely considered cultural factors. Transformation from near-primary forest to agroforestry had little effect on overall species richness, but reduced plant biomass and carbon storage by almost equal to75% and species richness of forest-using species by almost equal to60%. In contrast, increased land use intensity in cacao agroforestry, coupled with a reduction in shade tree cover from 80% to 40%, caused only minor quantitative changes in biodiversity and maintained high levels of ecosystem functioning while doubling farmers' net income. However, unshaded systems further increased income by almost equal to40%, implying that current economic incentives and cultural preferences for new intensification practices put shaded systems at risk. We conclude that low-shade agroforestry provides the best available compromise between economic forces and ecological needs. Certification schemes for shade-grown crops may provide a market-based mechanism to slow down current intensification trends.
Image segmentation is a cost-effective way to obtain information about the sizes and structural composition of agricultural parcels in an area. To accurately obtain such information, the parameters ...of the segmentation algorithm ought to be optimized using supervised or unsupervised methods. The difficulty in obtaining reference data makes unsupervised methods indispensable. In this study, we evaluated an existing unsupervised evaluation metric that minimizes a global score (GS), which is computed by summing up the intra-segment uniformity and inter-segment dissimilarity within a segmentation output. We modified this metric and proposed a new metric that uses absolute difference to compute the GS. We compared this proposed metric with the existing metric in two optimization approaches based on the Multiresolution Segmentation (MRS) algorithm to optimally delineate agricultural parcels from Sentinel-2 images in Lower Saxony, Germany. The first approach searches for optimal scale while keeping shape and compactness constant, while the second approach uses Bayesian optimization to optimize the three main parameters of the MRS algorithm. Based on a reference data of agricultural parcels, the optimal segmentation result of each optimization approach was evaluated by calculating the quality rate, over-segmentation, and under-segmentation. For both approaches, our proposed metric outperformed the existing metric in different agricultural landscapes. The proposed metric identified optimal segmentations that were less under-segmented compared to the existing metric. A comparison of the optimal segmentation results obtained in this study to existing benchmark results generated via supervised optimization showed that the unsupervised Bayesian optimization approach based on our proposed metric can potentially be used as an alternative to supervised optimization, particularly in geographic regions where reference data is unavailable or an automated evaluation system is sought.