At present, accessing and processing Earth Observation (EO) data on different cloud platforms requires users to exercise distinct communication strategies as each backend platform is designed differently. The openEO API (Application Programming Interface) standardises EO-related contracts between local clients (R, Python, and JavaScript) and cloud service providers regarding data access and processing, simplifying their direct comparability. Independent of the providers’ data storage system, the API mimics the functionalities of a virtual EO raster data cube. This article introduces the communication strategy and aspects of the data cube model applied by the openEO API. Two test cases show the potential and current limitations of processing similar workflows on different cloud platforms and a comparison of the result of a locally running workflow and its openEO-dependent cloud equivalent. The outcomes demonstrate the flexibility of the openEO API in enabling complex scientific analysis of EO data collections on cloud platforms in a homogenised way.

The Sentinel-1 satellites provide the formerly unprecedented combination of high spatial and temporal resolution of dual polarization synthetic aperture radar data. The availability of dense time series enables the derivation and analysis of temporally filtered annual backscatter signals. The study concentrates on the use of Sentinel-1 seasonal backscatter signatures for forest area estimation and forest type classification. A classification method based on time series similarity measures is introduced and tested in three test areas covered by various forest types including broadleaf temperate, boreal and montane forests. The results are compared with two European-wide Copernicus high resolution layers, namely forest type and tree cover density (TCD). The correspondence of forest/non-forest maps and TCD is high in all test areas, with overall accuracies for forest/non-forest classification between 86% and 91% and Pearson correlation coefficients for TCD between 0.68 and 0.74. The forest type classification (non-forest, coniferous and broadleaf forest classes) provides best results in temperate forests with an overall accuracy of 85%; in boreal forest, the accuracy decreases to only 65%. Generally, the method provides reliable results for forest area estimation, including regions where methods based on static parameters are often problematic (mountainous areas), and it enables forest type classification in temperate forests.

As the Earth system is exposed to large anthropogenic interferences, it becomes ever more important to assess the resilience of natural systems, i.e., their ability to recover from natural and human-induced perturbations. Several, often related, measures of resilience have been proposed and applied to modeled and observed data, often by different scientific communities. Focusing on terrestrial ecosystems as a key component of the Earth system, we review methods that can detect large perturbations (temporary excursions from a reference state as well as abrupt shifts to a new reference state) in spatio-temporal datasets, estimate the recovery rate after such perturbations, or assess resilience changes indirectly from stationary time series via indicators of critical slowing down. We present here a sequence of ideal methodological steps in the field of resilience science, and argue how to obtain a consistent and multi-faceted view on ecosystem or climate resilience from Earth observation (EO) data. While EO data offers unique potential to study ecosystem resilience globally at high spatial and temporal scale, we emphasize some important limitations, which are associated with the theoretical assumptions behind diagnostic methods and with the measurement process and pre-processing steps of EO data. The latter class of limitations include gaps in time series, the disparity of scales, and issues arising from aggregating time series from multiple sensors. Based on this assessment, we formulate specific recommendations to the EO community in order to improve the observational basis for ecosystem resilience research.

Terrestrial laser scanning can provide high-resolution, two-dimensional sampling of soil surface roughness. While previous studies demonstrated the usefulness of these roughness measurements in geophysical applications, questions about the number of required scans and their resolution were not investigated thoroughly. Here, we suggest a method to generate digital elevation models, while preserving the surface’s stochastic properties at high frequencies and additionally providing an estimate of their spatial resolution. We also study the impact of the number and positions of scans on roughness indices’ estimates. An experiment over a smooth and isotropic soil plot accompanies the analysis, where scanning results are compared to results from active triangulation. The roughness measurement conditions for ideal sampling are revisited and updated for diffraction-limited sampling valid for close-range laser scanning over smooth and isotropic soil roughness. Our results show that terrestrial laser scanning can be readily used for roughness assessment on scales larger than 5 cm, while for smaller scales, special processing is required to mitigate the effect of the laser beam footprint. Interestingly, classical roughness parametrization (correlation length, root mean square height (RMSh)) was not sensitive to these effects. Furthermore, comparing the classical roughness parametrization between one- and four-scan setups shows that the one-scan data can replace the four-scan setup with a relative loss of accuracy below 1% for ranges up to 3 m and incidence angles no larger than 50°, while two opposite scans can replace it over the whole plot. The incidence angle limit for the spectral slope is even stronger and is 40°. These findings are valid for scanning over smooth and isotropic soil roughness.