Monitoring forest–agriculture mosaics is crucial for understanding landscape heterogeneity and managing biodiversity. Mapping these mosaics from remotely sensed imagery remains challenging, since ecological gradients from forested to agricultural areas make characterizing vegetation more difficult. The recent synthetic aperture radar (SAR) Sentinel-1 (S-1) and optical Sentinel-2 (S-2) time series provide a great opportunity to monitor forest–agriculture mosaics due to their high spatial and temporal resolutions. However, while a few studies have used the temporal resolution of S-2 time series alone to map land cover and land use in cropland and/or forested areas, S-1 time series have not yet been investigated alone for this purpose. The combined use of S-1 & S-2 time series has been assessed for only one or a few land cover classes. In this study, we assessed the potential of S-1 data alone, S-2 data alone, and their combined use for mapping forest–agriculture mosaics over two study areas: a temperate mountainous landscape in the Cantabrian Range (Spain) and a tropical forested landscape in Paragominas (Brazil). Satellite images were classified using an incremental procedure based on an importance rank of the input features. The classifications obtained with S-2 data alone (mean kappa index = 0.59–0.83) were more accurate than those obtained with S-1 data alone (mean kappa index = 0.28–0.72). Accuracy increased when combining S-1 and 2 data (mean kappa index = 0.55–0.85). The method enables defining the number and type of features that discriminate land cover classes in an optimal manner according to the type of landscape considered. The best configuration for the Spanish and Brazilian study areas included 5 and 10 features, respectively, for S-2 data alone and 10 and 20 features, respectively, for S-1 data alone. Short-wave infrared and VV and VH polarizations were key features of S-2 and S-1 data, respectively. In addition, the method enables defining key periods that discriminate land cover classes according to the type of images used. For example, in the Cantabrian Range, winter and summer were key for S-2 time series, while spring and winter were key for S-1 time series.

There is an urgent need to quantify anthropogenic influence on forest carbon stocks. Using satellite-based radar imagery for such purposes has been challenged by the apparent loss of signal sensitivity to changes in forest aboveground volume (AGV) above a certain ‘saturation’ point. The causes of saturation are debated and often inadequately addressed, posing a major limitation to mapping AGV with the latest radar satellites. Using ground- and lidar-measurements across La Rioja province (Spain) and Denmark, we investigate how various properties of forest structure (average stem height, size and number density; proportion of canopy and understory cover) simultaneously influence radar backscatter. It is found that increases in backscatter due to changes in some properties (e.g. increasing stem sizes) are often compensated by equal magnitude decreases caused by other properties (e.g. decreasing stem numbers and increasing heights), contributing to the apparent saturation of the AGV-backscatter trend. Thus, knowledge of the impact of management practices and disturbances on forest structure may allow the use of radar imagery for forest biomass estimates beyond commonly reported saturation points.