Anthropogenically-driven climate change, land-use changes, and related biodiversity losses are threatening the capability of forests to provide a variety of valuable ecosystem services. The magnitude and diversity of these services are governed by tree species richness and structural complexity as essential regulators of forest biodiversity. Sound conservation and sustainable management strategies rely on information from biodiversity indicators that is conventionally derived by field-based, periodical inventory campaigns. However, these data are usually site-specific and not spatially explicit, hampering their use for large-scale monitoring applications. Therefore, the main objective of our study was to build a robust method for spatially explicit modeling of biodiversity variables across temperate forest types using open-access satellite data and deep learning models. Field data were obtained from the Biodiversity Exploratories, a research infrastructure platform that supports ecological research in Germany. A total of 150 forest plots were sampled between 2014 and 2018, covering a broad range of environmental and forest management gradients across Germany. From field data, we derived key indicators of tree species diversity (Shannon Wiener Index) and structural heterogeneity (standard deviation of tree diameter) as proxies of forest biodiversity. Deep neural networks were used to predict the selected biodiversity variables based on Sentinel-1 and Sentinel-2 images from 2017. Predictions of tree diameter variation achieved good accuracy (r² = 0.51) using Sentinel-1 winter-based backscatter data. The best models of species diversity used a set of Sentinel-1 and Sentinel-2 features but achieved lower accuracies (r² = 0.25). Our results demonstrate the potential of deep learning and satellite remote sensing to predict forest parameters across a broad range of environmental and management gradients at the landscape scale, in contrast to most studies that focus on very homogeneous settings. These highly generalizable and spatially continuous models can be used for monitoring ecosystem status and functions, contributing to sustainable management practices, and answering complex ecological questions. Full article

Timely crop yield forecasts at a national level are substantial to support food policies, to assess agricultural production, and to subsidize regions affected by food shortage. This study presents an operational crop yield forecasting system for Poland that employs freely available satellite and agro-meteorological products provided by the Copernicus programme. The crop yield predictors consist of: (1) Vegetation condition indicators provided daily by Sentinel-3 OLCI (optical) and SLSTR (thermal) imagery, (2) a backward extension of Sentinel-3 data (before 2018) derived from cross-calibrated MODIS data, and (3) air temperature, total precipitation, surface radiation, and soil moisture derived from ERA-5 climate reanalysis generated by the European Centre for Medium-Range Weather Forecasts. The crop yield forecasting algorithm is based on thermal time (growing degree days derived from ERA-5 data) to better follow the crop development stage. The recursive feature elimination is used to derive an optimal set of predictors for each administrative unit, which are ultimately employed by the Extreme Gradient Boosting regressor to forecast yields using official yield statistics as a reference. According to intensive leave-one-year-out cross validation for the 2000–2019 period, the relative RMSE for voivodships (NUTS-2) are: 8% for winter wheat, and 13% for winter rapeseed and maize. Respectively, for municipalities (LAU) it equals 14% for winter wheat, 19% for winter rapeseed, and 27% for maize. The system is designed to be easily applicable in other regions and to be easily adaptable to cloud computing environments such as Data and Information Access Services (DIAS) or Amazon AWS, where data sets from the Copernicus programme are directly accessible. Full article

Indicator 11.3. 1 of the 2030 sustainable development goals (SDG) 11, i.e., the ratio of the land use to the population growth rate, is currently classified by the United Nations as a Tier II indicator, as there is a globally-accepted methodology for its calculation, but the data are not available, nor are not regularly updated. Recently, the increased availability of remotely sensed data and products allows not only for the calculation of the SDG 11.3. 1, but also for its monitoring at different levels of detail. That is why this study aims to address the interrelationships between population development and land use changes in Poland and Lithuania, two neighboring countries in Central and Eastern Europe, using the publicly available remotely sensed products, CORINE land cover and GHS-POP. The paper introduces a map modelling process that starts with data transformation through GIS analyses and results in the geovisualisation of the LCRPGR (land use efficiency), the PGR (population growth rate), and the LCR (land use rate). We investigated the spatial patterns of the index values by utilizing hotspot analyses, autocorrelations, and outlier analyses. The results show how the indicators’ values were concentrated in both countries; the average value of SDG 11.3. 1, from 2000 to 2018 in Poland amounted to 0.115 and, in Lithuania, to −0.054. The average population growth ratio (PGR) in Poland equaled 0.0132, and in Lithuania, it was −0.0067, while the average land consumption ratios (LCRs) were 0.0462 and 0.0067, respectively. Areas with an increase in built-up areas were concentrated mainly on the outskirts of large cities, whereas outliers of the LCRPGR index were mainly caused by the uncertainty of the source data and the way the indicator is interpreted. Full article

Over the past few decades, the global population and the built environment’s vulnerability to natural hazards have risen dramatically. As a result, decisive actions, such as the SENDAI framework, have emerged to foster a global culture of successful disaster risk reduction policies, including actions to mitigate the social and economic impact of geohazards. The effective study of natural disasters requires meticulous and precise monitoring of their triggering factors, with ground- and space-based techniques. The integration of GNSS and SAR observations through the establishment of permanent infrastructures, i.e., Continuously Operating Reference Stations (CORS) networks and arrays of Corner Reflectors (CRs), may form a seamless ground displacement monitoring system. The current research literature provides fragmented guidelines, regarding the co-location of SAR and GNSS permanent infrastructures. Furthermore, there exist no guidelines for the determination of the most suitable locations using a holistic approach, in terms of criteria and required data. The purpose of this paper is to present a semi-automatic multicriteria site suitability analysis and evaluation of candidate sites for the installation of a permanent CORS and two CRs; one for each pass, taking into account various parameters and criteria. The first results demonstrate that the collocation of SAR and GNSS permanent infrastructures, utilizing a holistic criteria-based approach, is successful and complies with all the literature’s requirements. Full article

Monitoring of land cover plays an important role in effective environmental management, assessment of natural resources, environmental protection, urban planning and sustainable development. Increasing demand for accurate and repeatable information on land cover and land cover changes causes rapid development of the advanced, machine learning algorithms dedicated to land cover mapping using satellite images. Free and open access to Sentinel-2 data, characterized with high spatial and temporal resolution, increased the potential to map and to monitor land surface with high accuracy and frequency. Despite a considerable number of approaches towards land cover classification based on satellite data, there is still a challenge to clearly separate complex land cover classes, for example grasslands, arable land and wetlands. The aim of this study is to examine, whether a hierarchal classification of Sentinel-2 data can improve the accuracy of land cover mapping and delineation of complex land cover classes. The study is conducted in the Lodz Province, in central Poland. The pixel-based land cover classification is carried out using the machine learning Random Forest (RF) algorithm, based on a time series of Sentinel-2 imagery acquired in 2020. The following nine land cover classes are mapped: sealed surfaces, woodland broadleaved, woodland coniferous, shrubs, permanent herbaceous (grassy cover), periodically herbaceous (i.e., arable land), mosses, non-vegetated (bare soil) and water bodies. The land cover classification is conducted following two approaches: (1) flat, where all land cover classes are classified together, and (2) hierarchical, where the stratification is applied to first separate the most stable land cover classes and then classifying the most problematic once. The national databases served as the source of the reference sampling plots for the classification process. The process of selection and verification of the reference sampling plots is performed automatically. To assess the stability of the classification models the classification processes are performed iteratively. The results of this study confirmed that the hierarchical approach gave more accurate results compared to the commonly used flat approach. The median of the overall accuracy (OA) of the hierarchical classification was higher by 3–9 percentage points compared to the flat one. Of interest, the OA of the hierarchical classification reached 0.93–0.99, whereas the flat approach reached 0.90. Individual classes are also better classified in the hierarchical approach. Full article

Riparian zones fulfill diverse ecological and economic functions. Sustainable management requires detailed spatial information about vegetation and hydromorphological properties. In this study, we propose a machine learning classification workflow to map classes of the thematic levels Basic surface types (BA), Vegetation units (VE), Dominant stands (DO) and Substrate types (SU) based on multispectral imagery from an unmanned aerial system (UAS). A case study was carried out in Emmericher Ward on the river Rhine, Germany. The results showed that: (I) In terms of overall accuracy, classification results decreased with increasing detail of classes from BA (88.9%) and VE (88.4%) to DO (74.8%) or SU (62%), respectively. (II) The use of Support Vector Machines and Extreme Gradient Boost algorithms did not increase classification performance in comparison to Random Forest. (III) Based on probability maps, classification performance was lower in areas of shaded vegetation and in the transition zones. (IV) In order to cover larger areas, a gyrocopter can be used applying the same workflow and achieving comparable results as by UAS for thematic levels BA, VE and homogeneous classes covering larger areas. The generated classification maps are a valuable tool for ecologically integrated water management. Full article

Earth-observation-based mapping plays a critical role in humanitarian responses by providing timely and accurate information in inaccessible areas, or in situations where frequent updates and monitoring are required, such as in internally displaced population (IDP)/refugee settlements. Manual information extraction pipelines are slow and resource inefficient. Advances in deep learning, especially convolutional neural networks (CNNs), are providing state-of-the-art possibilities for automation in information extraction. This study investigates a deep convolutional neural network-based Mask R-CNN model for dwelling extractions in IDP/refugee settlements. The study uses a time series of very high-resolution satellite images from WorldView-2 and WorldView-3. The model was trained with transfer learning through domain adaptation from nonremote sensing tasks. The capability of a model trained on historical images to detect dwelling features on completely unseen newly obtained images through temporal transfer was investigated. The results show that transfer learning provides better performance than training the model from scratch, with an MIoU range of 4.5 to 15.3%, and a range of 18.6 to 25.6% for the overall quality of the extracted dwellings, which varied on the bases of the source of the pretrained weight and the input image. Once it was trained on historical images, the model achieved 62.9, 89.3, and 77% for the object-based mean intersection over union (MIoU), completeness, and quality metrics, respectively, on completely unseen images. Full article

The fraction of absorbed photosynthetic active radiation (FAPAR) is an essential climate variable for assessing the productivity of ecosystems. Satellite remote sensing provides spatially distributed FAPAR products, but their accurate and efficient validation is challenging in forest environments. As the FAPAR is linked to the canopy structure, it may be approximated by the fractional vegetation cover (FCOVER) under the assumption that incoming radiation is either absorbed or passed through gaps in the canopy. With FCOVER being easier to retrieve, FAPAR validation activities could benefit from a priori information on FCOVER. Spatially distributed FCOVER is available from satellite remote sensing or can be retrieved from imagery of Unmanned Aerial Vehicles (UAVs) at a centimetric resolution. We investigated remote sensing-derived FCOVER as a proxy for in situ FAPAR in a dense mixed-coniferous forest, considering both absolute values and spatiotemporal variability. Therefore, direct FAPAR measurements, acquired with a Wireless Sensor Network, were related to FCOVER derived from UAV and Sentinel-2 (S2) imagery at different seasons. The results indicated that spatially aggregated UAV-derived FCOVER was close (RMSE = 0.02) to in situ FAPAR during the peak vegetation period when the canopy was almost closed. The S2 FCOVER product underestimated both the in situ FAPAR and UAV-derived FCOVER (RMSE > 0.3), which we attributed to the generic nature of the retrieval algorithm and the coarser resolution of the product. We concluded that UAV-derived FCOVER may be used as a proxy for direct FAPAR measurements in dense canopies. As another key finding, the spatial variability of the FCOVER consistently surpassed that of the in situ FAPAR, which was also well-reflected in the S2 FAPAR and FCOVER products. We recommend integrating this experimental finding as consistency criteria in the context of ECV quality assessments. To facilitate the FAPAR sampling activities, we further suggest assessing the spatial variability of UAV-derived FCOVER to benchmark sampling sizes for in situ FAPAR measurements. Finally, our study contributes to refining the FAPAR sampling protocols needed for the validation and improvement of FAPAR estimates in forest environments. Full article

This contribution focuses on the utilization of very-high-resolution (VHR) images to identify construction areas and their temporal changes aiming to estimate the investment in construction as a basis for economic forecasts. Triggered by the need to improve macroeconomic forecasts and reduce their time intervals, the idea arose to use frequently available information derived from satellite imagery. For the improvement of macroeconomic forecasts, the period to detect changes between two points in time needs to be rather short because early identification of such investments is beneficial. Therefore, in this study, it is of interest to identify and quantify new construction areas, which will turn into build-up areas later. A multiresolution segmentation followed by a kNN classification is applied to WorldView images from an area around the southern part of Berlin, Germany. Specific material compositions of construction areas result in typical classification patterns different from other land cover classes. A GIS-based analysis follows to extract specific temporal “patterns of life” in construction areas. With the early identification of such patterns of life, it is possible to predict construction areas that will turn into real estate later. This information serves as an input for macroeconomic forecasts to support quicker forecasts in future. Full article

The Global Human Settlement Population Grid (GHS-POP) the latest released global gridded population dataset based on remotely sensed data and developed by the EU Joint Research Centre, depicts the distribution and density of the total population as the number of people per grid cell. This study aims to assess the GHS-POP data accuracy based on root mean square error (RMSE), mean absolute error (MAE), and mean absolute percentage error (MAPE) and the correlation coefficient. The study was conducted for Poland and Portugal, countries characterized by different population distribution as well as two spatial resolutions of 250 m and 1 km on the GHS-POP. The main findings show that as the size of administrative zones decreases (from NUTS (Nomenclature of Territorial Units for Statistics) to LAU (local administrative unit)) and the size of the GHS-POP increases, the difference between the population counts reported by the European Statistical Office and estimated by the GHS-POP algorithm becomes larger. At the national level, MAPE ranges from 1.8% to 4.5% for the 250 m and 1 km resolutions of GHS-POP data in Portugal and 1.5% to 1.6%, respectively in Poland. At the local level, however, the error rates range from 4.5% to 5.8% in Poland, for 250 m and 1 km, and 5.7% to 11.6% in Portugal, respectively. Moreover, the results show that for densely populated regions the GHS-POP underestimates the population number, while for thinly populated regions it overestimates. The conclusions of this study are expected to serve as a quality reference for potential users and producers of population density datasets. Full article

Due to the pathogenic nature of asbestos, a statutory ban on asbestos-containing products has been in place in Poland since 1997. In order to protect human health and the environment, it is crucial to estimate the quantity of asbestos–cement products in use. It has been evaluated that about 90% of them are roof coverings. Different methods are used to estimate the amount of asbestos–cement products, such as the use of indicators, field inventory, remote sensing data, and multi- and hyperspectral images; the latter are used for relatively small areas. Other methods are sought for the reliable estimation of the quantity of asbestos-containing products, as well as their spatial distribution. The objective of this paper is to present the use of convolutional neural networks for the identification of asbestos–cement roofing on aerial photographs in natural color (RGB) and color infrared (CIR) compositions. The study was conducted for the Chęciny commune. Aerial photographs, each with the spatial resolution of 25 cm in RGB and CIR compositions, were used, and field studies were conducted to verify data and to develop a database for Convolutional Neural Networks (CNNs) training. Network training was carried out using the TensorFlow and R-Keras libraries in the R programming environment. The classification was carried out using a convolutional neural network consisting of two convolutional blocks, a spatial dropout layer, and two blocks of fully connected perceptrons. Asbestos–cement roofing products were classified with the producer’s accuracy of 89% and overall accuracy of 87% and 89%, depending on the image composition used. Attempts have been made at the identification of asbestos–cement roofing. They focus primarily on the use of hyperspectral data and multispectral imagery. The following classification algorithms were usually employed: Spectral Angle Mapper, Support Vector Machine, object classification, Spectral Feature Fitting, and decision trees. Previous studies undertaken by other researchers showed that low spectral resolution only allowed for a rough classification of roofing materials. The use of one coherent method would allow data comparison between regions. Determining the amount of asbestos–cement products in use is important for assessing environmental exposure to asbestos fibres, determining patterns of disease, and ultimately modelling potential solutions to counteract threats. Full article

Accurately identifying Natura 2000 habitat areas with the support of remote sensing techniques is becoming increasingly feasible. Various data types and methods are used for this purpose, and the fusion of data from various sensors and temporal periods (terms) within the phenological cycle allows natural habitats to be precisely identified. This research was aimed at selecting optimal datasets to classify three grassland Natura 2000 habitats (codes 6210, 6410 and 6510) in the Ostoja Nidziańska Natura 2000 site in Poland based on hyperspectral imagery and botanical on-ground reference data acquired in three terms during one vegetative period in 2017 (May, July and September), as well as a digital terrain model (DTM) obtained by airborne laser scanning (ALS). The classifications were carried out using a random forest (RF) algorithm on minimum noise fraction (MNF) transform output bands obtained for single terms, as well as data fusion combining the topographic indices (TOPO) calculated from the DTM, multitemporal hyperspectral data, or a combination of the two. The classification accuracy statistics were analysed in various combinations based on the datasets and their terms of acquisition. Topographic indices improved the classification accuracy of habitats 6210 and 6410, with the greatest impact noted in increased classification accuracy of xerothermic grasslands. The best terms for identifying specific habitats were autumn for 6510 and summer for 6210 and 6410, while the best results overall were obtained by combining data from all terms. The highest obtained values of the F1 coefficient were 84.5% for habitat 6210, 83.2% for habitat 6410, and 69.9% for habitat 6510. Comparing the data fusion results for habitats 6210 and 6410, greater accuracy was obtained by adding topographic indices to multitemporal hyperspectral data, while for habitat 6510, greater accuracy was obtained by fusing only multitemporal hyperspectral data. Full article

As a result of tightened waste regulation across Europe, reports of waste crime have been on the rise. Significant stockpiles of tyres and plastic materials have been identified as a threat to both human and environmental health, leading to water and livestock contamination, providing substantial fuel for fires, and cultivating a variety of disease vectors. Traditional methods of identifying illegal stockpiles usually involve laborious field surveys, which are unsuitable for national scale management. Remotely-sensed investigations to tackle waste have been less explored due to the spectrally variable and complex nature of tyres and plastics, as well as their similarity to other land covers such as water and shadow. Therefore, the overall objective of this study was to develop an accurate classification method for both tyre and plastic waste to provide a viable platform for repeatable, cost-effective, and large-scale monitoring. An augmented land cover classification is presented that combines Copernicus Sentinel-2 optical imagery with thematic indices and Copernicus Sentinel-1 microwave data, and two random forests land cover classification algorithms were trained for the detection of tyres and plastics across Scotland. Testing of the method identified 211 confirmed tyre and plastic stockpiles, with overall classification accuracies calculated above 90%. Full article