Global Open Data Remote Sensing Satellite Missions for Land Monitoring and Conservation: A Review

: The application of global open data remote sensing satellite missions in land monitoring and conservation studies is in the state of rapid growth, ensuring an observation with high spatial and spectral resolution over large areas. The purpose of this study was to provide a review of the most important global open data remote sensing satellite missions, current state-of-the-art processing methods and applications in land monitoring and conservation studies. Multispectral (Landsat, Sentinel-2, and MODIS), radar (Sentinel-1), and digital elevation model missions (SRTM, ASTER) were analyzed, as the most often used global open data satellite missions, according to the number of scientiﬁc research articles published in Web of Science database. Processing methods of these missions’ data consisting of image preprocessing, spectral indices, image classiﬁcation methods, and modelling of terrain topographic parameters were analyzed and demonstrated. Possibilities of their application in land cover, land suitability, vegetation monitoring, and natural disaster management were evaluated, having high potential in broad use worldwide. Availability of free and complementary satellite missions, as well as the open-source software, ensures the basis of e ﬀ ective and sustainable land use management, with the prerequisite of the more extensive knowledge and expertise gathering at a global scale.


Introduction
Land conservation strategies arose as a growing necessity in countries worldwide to prevent overexploitation of natural resources and to protect the ecosystem. To satisfy these demands, remote sensing satellite missions became the fundamental data source for land conservation and sustainability globally [1]. The application of remote sensing in land monitoring and conservation is based on the analysis of the reflectance or emittance of individual parts of the electromagnetic spectrum with an observed land cover. Each land cover and its properties result in a unique set of values in the observed parts of the spectrum, being represented by spectral signatures. By analyzing the spectral signatures of the same area observed over different time periods, changes in land cover and its properties are detected. Satellite missions have been one of the main data sources in land monitoring and conservation for many decades, characterized by a wide imaging swath and the ability of integration with various compatible types of imaging sensors [2]. The types of satellite missions differ primarily according to the wavelengths of the electromagnetic spectrum to which sensors mounted on satellites are sensitive. The integration of remote sensing data with the geographical information system (GIS) creates a database of spatial data of the environment, enabling the implementation of spatial analyses for land-use planning [3]. The development of such methods through of multitemporal satellite images and their analysis in the GIS environment for monitoring changes in land cover, especially for the monitoring of vegetation properties [15]. Satellite mission data in the same document are classified in the category of primary data sources, which form the basis of the analysis of environmental indicators and future reports. The aim of this paper was to analyze the most commonly used global open data satellite missions, data processing methods and their application in land monitoring and conservation. The primary segment of their importance is accessibility to the broad scientific community for the application in land conservation studies and strategies. Further development of these missions and familiarization with their possibilities provides the basis for the upgrade of land monitoring services and their widespread use globally.

Global Open Data Remote Sensing Satellite Missions
The global open data remote sensing satellite missions are analyzed according to the classification to multispectral missions, radar missions, and DEM acquiring. Major missions per classification were highlighted according to their popularity in articles for land monitoring and conservation using remote sensing indexed in WoSCC database in the last ten years (2010-2019). All representative figures in the paper were created by the authors, based on the freely available global open data satellite mission imagery.

Multispectral Satellite Missions
Remote sensing using multispectral sensors is performed by registering incoming reflected or emitted energy from objects on the earth's surface, which is dispersed and registered in sensors sensitive to certain spectral bands [16]. These spectral bands are represented as a narrow part of the electromagnetic spectrum, defined by the smallest and largest wavelength of the sensor sensitivity, which results in one raster image per spectral band. Spectral bands of four most commonly used global open data multispectral satellite missions are displayed in Figure 2, where black squares represent their spectral coverage in visible, red-edge, and infrared parts of the spectrum. The intensity of the incoming energy on the sensor is converted to digital numbers, which also depend on atmospheric conditions at the sensing time. As a result, the numerical value shows the relative relationship of reflection at the time of observation and does not represent a measurable physical unit [17]. For the quantitative interpretation of such data, it is necessary to perform the conversion of numerical values into reflectance, or the ratio of the amount of reflected energy per spectral band and the total incoming energy at the sensor. Detection and identification of land cover and its properties are based on differences of spectral signatures per spectral band. Each specific property of land cover results in unique spectral signatures, with a noticeable difference in its reflectance in one or more spectral bands compared to the rest of the observed area [18]. By highlighting these differences using the processing methods described in the next chapter, these characteristics are quantified or extracted from the rest of the observed area, which serves as a basis for interpretation and decision-making in land conservation. Sentinel-2 is a multispectral satellite mission with moderate spatial resolution within the ESA Copernicus program. The main purpose of the mission is the observation of the Earth's surface to provide forest and vegetation monitoring services, detection of land cover changes and management of natural disasters [19]. The imaging is performed using a Multispectral Imager (MSI) instrument mounted on two satellites, Sentinel-2A, and Sentinel-2B. The satellites are polar-orbiting on diametrically opposite sides of the same sun-synchronous orbit, allowing the observation of the Earth's surface between 84 • N and 56 • S latitudes.
The Landsat 8 satellite mission, owned by the National Aeronautics and Space Administration (NASA), uses two multispectral sensors, the Operational Land Imager (OLI) and the Thermal Infrared Sensor (TIRS). OLI uses nine spectral bands in the visible, near-infrared, and short-wave infrared part of the spectrum for moderate spatial resolution observations, while TIRS complements multispectral observation with two thermal bands with a lower spatial resolution [20]. Landsat 8 is an upgrade to the previous mission of the same series, Landsat 7, through better spectral resolution and a rectified error in observing the edge parts of the imaging area.
MODIS is a multispectral satellite mission owned by NASA and named after a multispectral sensor mounted on two satellites, Terra and Aqua, orbiting on diametrically opposite sides of the same polar orbit. MODIS is intended for remote sensing of land, oceans, and lower atmosphere, enabling high temporal resolution and a very wide imaging swath [21]. Due to the complexity and diversity of MODIS products, MODIS09GA data were selected as a representative product regarding average spatial and spectral resolution in relation to other MODIS products.
The basic characteristics of the most commonly used satellite multispectral missions applied in environmental studies are presented in Table 1. The supported research area is classified according to the classification by Herod [22], arranged from smaller to larger in a local, national, regional, and global scale. The greatest impact for this classification was the spatial resolutions of each mission, which are comparatively represented in Figure 3. Landsat missions contain a panchromatic band of 15 m spatial resolution, which serves as the basis for pan-sharpening of other images in the visible and near-infrared part of the spectrum to 15 m spatial resolution while preserving spectral values. The same procedure is possible by applying a 10 m spatial resolution visible light bands of Sentinel-2 for pan-sharpening of the red-edge and near-infrared images [23].  local, national national, regional national, regional regional, global

Radar Satellite Missions
Microwave sensors, commonly known as radars, belong to a group of active sensors. The method of registering the intensity of reflected energy by radars depends on the radar system's properties and terrain characteristics, primarily surface roughness and water content in soil and vegetation [24]. The specific characteristics of each radar system in remote sensing are the polarization of the electromagnetic waves, the angle of energy transmittance and its wavelength [25]. Polarization of the emitted or reflected electromagnetic wave represents the orientation of the electric field at the time of transmission or registration, which may be horizontal (H) or vertical (V). If the radar is sensitive to both polarizations, one or more images of different combinations of transmitted and registered polarizations in combinations of HH, HV, VH, or VV are created. Due to the differences in the registration of reflected polarization, these images contain different, most commonly complementary information of land cover and its properties [26].
Sentinel-1 is a radar satellite mission, the first of six planned missions within the ESA Copernicus program. It consists of two Sentinel-1 satellites carrying a synthetic aperture radar (SAR) and orbiting on diametrically opposite sides of the same polar orbit, resulting in a temporal resolution of six days [27]. This SAR is sensitive to wavelengths in the C-band of the electromagnetic spectrum, which covers wavelengths ranging from 3.8 to 7.5 cm. The advantage of remote sensing using radar is the ability of observation regardless of atmospheric conditions, as microwaves pass through fog, rain, and clouds [28]. The radar images at the pre-processing Level-1 available for download contain raw SAR observation data which were previously internally calibrated to the mean signal frequency using the Doppler centroid method. Radar imaging is performed by single or double polarization, which results in a spatial resolution in range of 10 and 40 m. The basic two products of Sentinel-1 are single look complex (SLC) and ground range detected (GRD) images [27]. SLC images consist of directed SAR data containing the complete phase and related complex data. These images are georeferenced using the data of the satellite orbits and positions at the sensing time, corrected for oblique observations of the terrain. GRD images contain generalized directed SAR data, projected in the WGS84/UTM projection coordinate system using an ellipsoid model. These images contain pixels with approximately square shape and reduced noise, but with a lower spatial resolution than SLC images. Phase data is therefore lost, but GRD images take up to five times less hard disk memory than SLC images.

Digital Elevation Models
The DEMs from global satellite missions contain absolute elevation as pixel values in raster data, indicating the average elevation in the observed area. Two basic methods of DEM calculation using satellite missions are based on the radar interferometry and image stereo pairs. The two most widely used global digital height models are the SRTM and the ASTER owned by NASA, according to the number of published land monitoring and conservation research articles in WoSCC in the last decade. SRTM contains elevation data with a spatial resolution of approximately 30 m, which is currently the highest spatial resolution of global open data DEMs.
SRTM is a global satellite mission for radar imaging of the Earth's surface used to calculate the global set of digital elevation data. The observations were performed in the range of 60 • N to 56 • S latitudes [29]. The basic principle of data collection and processing is radar interferometry, based on the comparison of two radar images of the same area sensed at different angles. ASTER is also a global satellite mission with the goal of DEM calculation, named after an instrument mounted on the Terra satellite. To create a DEM, remote sensing was performed in the near-infrared spectral band and processed by the image stereo pairs analysis method [30]. An additional property of the ASTER mission is the corrections for the elevation values of water surfaces, whereby the seas and oceans are characterized by uniform elevation, while rivers were modelled by proportionally lower elevations downstream.

Methods of Processing and Analyzing Remotely Sensed Data
The processing and analysis of data collected by one of the previously described satellite missions can generally be divided into three directions: visual, quantitative, and qualitative. A prerequisite for processing and analyzing with any approach is the image pre-processing [31]. It represents a set of corrections for remotely sensed data, which is specific for each individual satellite mission. The visual analysis method is the simplest and relatively limited approach for the interpretation of remotely sensed data. It is commonly applied in terrain reconnaissance and the identification of generalized land cover classes [32]. The visual interpretation is applicable for multispectral images using RGB composites, where it is possible to display terrain characteristics specific for the selected combination of spectral bands. The RGB composites based on Sentinel-2 images for specialized purposes in visual analysis are displayed in Figure 4. All used spectral bands were resampled to 20 m spatial resolution prior to the creation of RGB composites.
The quantitative approach is dominantly based on the application of multispectral images, by combining pixel values of several different spectral bands. The selected property of the observed area is quantified by the application of spectral indices, where the values of that property are proportional to the values of the calculated spectral index [33]. Mathematical models for modelling of the topographic terrain parameters based on the DEM also quantitatively evaluate the calculated properties. Qualitative evaluation of the terrain is based on the classification of spectral values of image pixels or objects in the predefined number of classes. This approach is flexible in terms of selecting spectral values as an input in the classification and does not necessarily require calibration of remotely observed values to physical units, but rather implements the relative classification of parts of the observed area. Open source GIS software, SAGA GIS v7.4.0 and QGIS v3.8.3, were used for data processing in case studies of global open data satellite missions in the next subsections.

Image Pre-Processing
The basic pre-processing procedures of any data obtained by remote sensing are the transformation into a selected coordinate reference system and grid clipping to an area of interest. Remote sensing data of the missions described above are commonly distributed in WGS84 or WGS84/UTM projection coordinate systems worldwide.
Multispectral images of Sentinel-2, Landsat and MODIS missions are mostly available to users at level 1C pre-processing and contain orthorectified data with geometrically calibrated reflectance values at the top of atmosphere (TOA). Data at earlier pre-processing levels is not openly distributed to users and is performed by organizations managing the satellite mission. Pre-processing corrections for Sentinel-2 satellite images performed by the ESA are displayed in Figure 5. Using the atmospheric corrections, users have the possibility of independent conversion of TOA reflectance to the bottom of the atmosphere (BOA) reflectance, which corresponds to the 2A pre-processing level. This eliminates the influence of atmospheric factors on reflectance at the sensing time, which makes is available for quantitative processing and analysis of the observed area. A comparative display of Sentinel-2 RGB natural composites based on TOA reflectance from Level-1C data and BOA reflectance from Level-2A data is shown in Figure 6. Some satellite missions, primarily Sentinel-2 and Landsat 8, have enabled free downloading of Level-2A images in recent years, with the conversion to the BOA reflection coefficient carried out by data distributors. Further Level-2 and Level-3 data designate more complex remote sensing products, such as land cover classification.   The pre-processing of radar images consists of four basic steps: radiometric correction, signal noise filtering, geometric correction, and georeferencing, conversion of registered signal intensity value to decibels (dB) (Figure 7) [34]. Radiometric calibration results in the conversion of the originally registered intensity to absolute normalized radar backscatter intensity using coefficients from the image metadata. Filtering or radar images reduces the amount of noise on the images, caused by coherent processing of return radar signals from multiple different surfaces on the ground within a single pixel. Frequently used filters are speckle filters, which smooth the image by calculating the mean pixel value from the selected pixel matrix [35]. Selecting a larger filtering pixel matrix proportionally reduces noise and spatial resolution of the image. For the pre-processing of the images displayed in Figure 7, the Lee Sigma filter with a matrix of 7 × 7 pixels was applied. Geometric correction eliminates the influence of relief and angle of radar observation on image values using orthorectification [36].
During the georeferencing process, the radar images are transformed into the correct position in the selected coordinate reference system. Doppler's orthorectification method of SAR images into the Croatian Terrestrial Reference System (HTRS96/TM) was also performed, using a global SRTM 1-arc second DEM for the orthorectification. The pre-processing specific for hydrological analysis of the terrain using a DEM is a mathematical filling of depressions, on which it is not possible to determine the direction of water flow. This correction is usually carried out by the Fill Sinks method, where it is necessary to define the value of the slope tolerance of the elevation filling [37].

Spectral Indices
Spectral indices are calculated using arithmetic functions of pixel values from two or more remotely sensed spectral bands (Table 2; Figure 8). Each spectral index aims to extract a certain type of land cover or its specific property, which is characterized by distinct values relative to the rest of the observed area. Spectral indices are often calculated using normalized differences of reflectance in two spectral bands, with the highest and lowest sensed relative reflectance of a particular land cover property. The advantage of this approach is the known possible range of spectral index values at interval (-1, 1), which makes it possible to achieve repeatability of the values of specific land cover property and therefore monitoring its changes in time. Spectral indices whose primary purpose is to enable monitoring of vegetation properties are known as vegetation indices. Their quantity and the range of applications in environmental monitoring are the widest of all spectral indices, as there are currently about 250 unique vegetation indices that can be calculated using Sentinel-2 satellite images [46]. Some of the most commonly used vegetation indices are the normalized difference vegetation index (NDVI), soil-adjusted vegetation index (SAVI), enhanced vegetation index (EVI), normalized difference red edge (NDRE), and normalized green-red vegetation index (NGRDI). The calculation of vegetation indices is based on the emphasis of the difference in the vegetation reflectance from very low in red, to very high in the near-infrared or red-edge spectral bands. With the application of low-cost RGB cameras, which offer the possibility of remote sensing only in the visible part of the spectrum, a green band is often used as an alternative to the red-edge or near-infrared band [47]. Almost all vegetation indices have the same primary purpose, but each is specific regarding highlighting certain vegetation properties or noise reduction due to reflectiveness of the bare soil. Although NDVI is the world's most widely used vegetation index for a wide range of applications in the monitoring of vegetation properties [48], it has a limitation in the form of low sensitivity to minor differences in case of high chlorophyll content and biomass, known as the saturation effect [49]. In this case, NDRE provides the alternative due to the greater sensitivity of the red-edge band to these vegetation properties [50]. Due to mutual compatibility, NDVI and NDRE are often used for the remote sensing of agricultural crops by applying sensors on agricultural machinery during precise fertilization and crop protection. The SAVI uses a correction factor to reduce the influence of bare soil reflectance, while a blue band for the calculation of EVI is used to reduce the influence of atmospheric effects during remote sensing [51]. The secondary purpose of vegetation indices is to extract areas under vegetation from the soil and artificial land covers due to their sensitivity to the chlorophyll content and biomass. For this reason, they are often used in combination with other processing methods of remote sensing data, resulting in high processing efficiency and high accuracy of land observation models [47]. Modified normalized difference water index (MNDWI) is used as a complement to vegetation indices when monitoring vegetation since it is sensitive to water content in vegetation and soil. It is a modified version of another commonly used spectral index, normalized difference water index (NDWI), which was upgraded by reducing noise of the detection of soil and vegetation water content by using green, instead of near-infrared spectral band [43]. MNDWI is also applied for land cover mapping, as it achieves maximum values for water surfaces. Normalized difference built-up index (NDBI) and normalized difference soil index (NDSI) are used primarily in monitoring of land cover changes for built-up objects and soil and their distinction, as these land cover classes often have similar spectral signatures.

Image Classification
The classification methods in remote sensing are based on the classification of the fundamental image parts into a predefined number of classes. Classification is often an indispensable procedure in the remote sensing for land monitoring and conservation, as it allows determination of land cover classes and its properties with high accuracy over a large area [52]. The necessity for manual field observation is consequentially decreased to discrete locations, used as a test dataset for the accuracy assessment of classification results. The fundamental image parts in classification are either pixels or objects, defined as two or more neighboring pixels of similar spectral and textural properties. These parts are classified in the n-dimensional space, where n represents the number of input grids in the classification. Input grids in the classification are the most commonly original images and spectral indices. Various spectral indices, like NDVI, NDBI, and MNDWI, were successfully integrated with the classification using Landsat 7 and Landsat 8 images by Estoque and Murayama [53]. The object-oriented classification is the superior classification approach, with the advantage of lower noise in classification results caused by atmospheric effects, as well as a more generalized result that accurately represents the actual situation in the field than pixel-based classification [47].
Classification in the number of classes defined by the operator can be performed either without the training dataset (unsupervised classification) or by determining the training dataset for individual classes of selected land cover properties (supervised classification). Training datasets contain representative samples of each defined class in form of groups of pixels or objects, created based on field observations or visual identification on the images. Unsupervised classification algorithms are based on the iterative classification of the fundamental image parts into a defined number of classes according to a defined non-parametric criterion [54]. The most commonly used unsupervised classification algorithms are k-means, which is based on the minimum Euclid distance of spectral values to the center of the class and the iterative uniform adjustment of the class center and ISODATA, which similarly adjusts the centers of classes based on the criterion of the least distance to the center of the class. Conventional supervised classification algorithms operate on a similar principle, but class centers are defined by spectral values of the training dataset and are not subject to iterative adjustment to input data. The main property of this approach is the subjectivity of the operator, which results in its greater influence on the classification process due to the possibility of arbitrary formation of class centers in the n-dimensional system of input values. This can result in multiple classes of close spectral values, which often correspond to the actual conditions in the field. For an independent accuracy assessment of the classification results by supervised classification, the creation of a ground-truth dataset containing representative class samples using field identification or vectorization based on visual identification of images is common. This unique dataset is then randomly divided into training and test datasets in ratios of 30:70 or 20:80, resulting in the reduced operator's subjectivity during the classification process.
Some of the conventional algorithms of supervised classification are maximum likelihood and minimum distance to mean. Over the past decade, machine learning algorithms have completely replaced conventional algorithms of supervised classification from application in land monitoring and conservation. The superior accuracy of classification results, reduced overfitting and computing efficiency are the main advantages of machine learning algorithms over conventional supervised classification algorithms [55]. The most commonly used machine learning algorithms for supervised classification in remote sensing currently are random forest (RF) [56], artificial neural networks (ANN) [57] and support vector machines (SVM) [58]. The performance of these algorithms is visible in Figure 9, which contains the classification results calculated using the same training dataset on Sentinel-2 images. Resulting classes were represented according to boxplots containing NDVI values, as the values of spectral indices commonly differ with land cover classes [59].
The current trend of machine learning algorithms development in remote sensing are deep machine learning algorithms [60]. Their advantage over the machine learning algorithms described above is the higher classification accuracy and processing capabilities of a larger quantity of training samples, but their characteristics and performance still require additional research. The implementation of these algorithms requires expensive computer hardware, primarily in terms of the graphics processing unit (GPU), which currently makes them less available for use in a broad scientific community.

Modelling of Terrain Topographic Parameters
Modelling of terrain topographic parameters is based on using a DEM, according to the neighborhood values of individual pixels, as well as the latitude of the specific location. Calculated models represent theoretical values and approximation of actual conditions in the field at the micro-level. The significance of such models is in the relative relationship of calculated values, which for spatial planning and suitability analyses allows adequate representation of the terrain topography. Figure 10 displays various terrain topographic models calculated using SRTM DEM: slope [61], wind exposition index, potential incoming solar irradiation [62], flow accumulation model [63], and topographic openness [64]. Terrain slope is the basic topographical parameter of the terrain and is widely used in spatial planning in civil engineering, transport, and environmental studies [65]. Potential solar irradiation was calculated for the year 2020, using a 70% generalized atmospheric permeability parameter, which includes the influences of clouds, fog, and precipitation on the amount of incoming solar radiation. The wind exposure index is calculated by analyzing the individual elevations of pixels and adjacent pixels incrementally every 15 • , covering the entire horizon around a pixel of 360 • . The baseline value for the flat area is 1, with lower values proportionally indicating more covered areas, while higher values indicate a more exposed area to wind. The flow accumulation model was calculated using the top-down method and the multiple flow direction technique, gradually analyzing the relative elevation difference of individual pixels, starting with those with the highest elevation value. For each pixel, eight neighboring pixels were observed, and the final values of the flow accumulation model represent the total surface from which precipitation or water potentially flows from higher elevations. This model represents a single component of the hydrographic analysis of the terrain, as precipitation drainage and water flow depend on several factors, such as land cover class or soil texture [66]. The topographic openness index indicates the horizon exposure for a single pixel, which is applicable for spatial planning in the construction and analysis of animal and plant habitats.

Applications of Land Monitoring and Conservation Using Satellite Missions
The application of global open data satellite missions covers a wide range of possibilities, most commonly for monitoring of land cover and vegetation changes, environmental and natural disasters, as well as the land suitability calculation using GIS-based multicriteria analysis. The development of new satellite missions is also expected to improve existing atmospheric and ocean monitoring systems.

Land Cover Change Monitoring
The observation of land cover changes over a period of time is one of the fundamental sources of information in land conservation. The classification of the observed area into several generalized land cover classes creates a universal basis for spatial planning and monitoring of changes in the environment. Some of the most important environmental factors, such as climate change, deforestation, and urbanization, are determined by comparing land cover databases created for two or more different time periods [67]. The same data are used to quantify biodiversity and detect endangered areas due to deforestation or various socio-economic reasons [68]. Global satellite missions, especially multispectral missions, are the most common source of data for the classification of land cover due to the high temporal resolution and the suitability for efficient computer processing [69]. For the monitoring of land cover changes, standardization of land cover classes over different time periods is necessary for their mutual comparison [70]. The most commonly used land cover classes within the EU are artificial (urban) areas, agricultural areas, forests, wetlands, and water bodies, according to the Corine Land Cover (CLC) classification within the ESA Copernicus program. The latest version of the land cover database in the EU, CLC 2018, was produced using the Sentinel-2 and Landsat 8 multispectral images with thematic accuracy above 85% [71]. For specialized studies and applications in environment monitoring, an additional classification of generalized classes in multiple subclasses is performed. For example, the agricultural areas class is divided into four subclasses in the CLC 2018: arable land, permanent plantations, pastures, and heterogeneous agricultural land, and a further 11 subclasses at the third level [71]. Algorithms of supervised machine learning classification are currently a fundamental method of processing satellite images for the land cover classification [72], and more frequent application of deep machine learning algorithms for the same purpose is expected in the future [73]. Trends in the development of classification algorithms of land cover from satellite images are also directed to their automation, thus reducing the processing time, and eliminating the possibility of operator subjective error in determining the training dataset for classification [74].

Land Suitability Determination Studies Using GIS-Based Multicriteria Analysis
The GIS-based multicriteria analysis is a fundamental approach to the creation of land suitability studies, characterized by a high level of flexibility and computational efficiency. The premise for the implementation of multicriteria analysis is the spatial component of the study objective and the criteria on which land suitability depends. The suitability calculation of environmental factors using GIS multicriteria analysis generally consists of six basic procedures [75,76]: 1.
Definition of the study aim 2.
Selection of spatially related criteria affecting the land suitability 3.
Standardization of criteria values 4.
Determination of criteria weights 5.
Calculation of land suitability by combining standardized values and criteria weights 6.
Accuracy assessment and interpretation of results The application of GIS multicriteria analysis in land conservation is based on the optimal utilization of existing natural resources and the preservation of their quality and quantity. The objectives of GIS multicriteria analysis in land conservation are often based on multidisciplinary sciences, covering sustainable waste management, agricultural production planning, groundwater protection, natural disaster risk mapping and spatial planning of energy plants. Two main objectives in sustainable waste management are the suitability analyses of the existing and planned landfills [77], as well as the suitability of waste incinerator locations [78]. The appropriate waste management method, especially in urban environments, reduces air, water and soil pollution caused by the discharge of ammonia, heavy metals, and nitrates from waste [79]. The application of GIS-based multicriteria analysis in agriculture focuses on the calculation of cropland suitability according to climatic, pedological, and topographical terrain conditions [65]. The selection of the optimal area for crop cultivation reduces the need for the application of fertilizer and pesticides, resulting in a lower release of heavy metal contaminants into the environment. Groundwater is the most important natural resource for the global water supply [80]. Using GIS-based multicriteria analysis, zoning of groundwater properties is indirectly performed by several indicators affecting groundwater status, such as climatic, topographic, and geological criteria, as well as certain soil properties, such as drainage [81]. Natural disasters risk mapping using GIS-based multicriteria analysis enables the development of an emergency plan for forest wildfires, floods, and other disasters. The multiannual climate and hydrological criteria in such analyses provide the basis for the insurance of crops due to damage caused by hail or drought [82]. Construction planning of costly and environmentally-friendly power plants is based on the existing renewable energy sources and natural habitats, spatially modelled using satellite mission data. The GIS-based multicriteria analysis is one of the fundamental procedures for the energy plants construction according to land conservation requirements, such as for solar power plants [83] and windmills [84]. Criteria modelled using satellite missions used in GIS-based multicriteria analyses for the objectives described above are presented in Table 3. Table 3. Criteria used in geographical information system (GIS)-based multicriteria analyses in environmental protection modelled from satellite mission data.

Groundwater protection
Slope, flow accumulation model, land cover [80,81] Natural disasters risk mapping DEM, slope, total biomass, dry vegetation mass, land cover [82,85] Spatial planning in energetics Land cover, slope, DEM, solar irradiation, wind exposition index [83,84] By the standardization process, input criteria values are transformed into a uniform number interval. Conventional methods, such as linear stretching and stepwise classification, are most commonly used for standardization, while fuzzy methods allow for more advanced standardization and greater expert subjective influence [75]. Criteria weight calculation is usually performed by pairwise comparison of relative impact on suitability within the analytic hierarchy process (AHP). Other known methods of weight determination of criteria are TOPSIS, ELECTRE, and PROMETHEE [86]. The standardized criteria values and their weights are commonly combined using the weighted linear combination, resulting in the land suitability values of the unconstrained study area. Satellite mission data, aside for the spatial modelling of criteria, are used for accuracy assessment and the interpretation of land suitability results. Vegetation indices calculated using multitemporal satellite images were successfully used for the suitability model validation for soybean [75] and wheat land suitability [87].

Monitoring of Vegetation Properties
The determination and monitoring of vegetation properties using satellite missions are performed by a combination of classification methods and vegetation indices, calculated based on multitemporal satellite images [88]. Red-edge and near-infrared spectral bands of multispectral satellite missions are primarily used for the determination of the leaf area index (LAI), as a fundamental vegetation property with high accuracy [89]. By analyzing the spatiotemporal LAI values of forests and agricultural areas, it is possible to accurately detect anomalies in the vegetation biomass and to determine the nature of their causes through field inspection. Other specific vegetation properties that can be accurately determined using multispectral satellite images are crop yield, plant nitrogen, and chlorophyll contents, as well as the vegetation stress caused by soil contamination of heavy metals [90,91]. Shortwave and thermal infrared spectral bands from multispectral satellite missions and radar satellite mission data are sensitive to vegetation water content and are used to model evapotranspiration and detect areas prone to water stress due to drought [92]. These data serve as a basis for the management of forests and agricultural crops, as well as for suitability analyses for the establishment of irrigation systems. Landsat 8 images also enable the monitoring of multiple factors affecting the urban vegetation, such as land surface temperature, urban green spaces and their proximity to water bodies and built-up area [93]. Moderate and high-resolution multispectral satellite missions are successfully used in weed detection, precise fertilization, and the determination of agricultural crop density in combination with observations using UAVs [47,94]. The application of these data in precise agriculture creates optimal conditions for crop growth with significant savings in fertilizer and pesticides, compared to the conventional farming approach. Reducing their application also reduces agricultural production costs, while also benefiting the environment and biodiversity conservation through the reduction of agricultural land contamination with heavy metals [95].

Management of Ecological and Natural Disasters
Timely response and management of damage caused by ecological and natural disasters of damage are important environmental factors, with satellite missions providing the basis for effective decision-making. High temporal resolution is the most important feature of satellite missions in such cases, making MODIS, Sentinel-1 and Sentinel-2 satellite missions a core data source. Open data satellite missions for specialized applications are in constant development, such as the visible infrared imaging radiometer suite (VIIRS) mission owned by NASA to automatically report new wildfires globally in near-real-time. Ecological and natural disasters monitored by global satellite missions consist of oil spills in the seas and oceans [96], glacier melting in the Arctic [97], forest wildfires [98], floods [32], and damage to agricultural crops affected by hailstorms [99]. Monitoring of the spatial coverage of ecological and natural disasters based on the mentioned studies represents the basis for remediating the resulting damage. It is accomplished through the procedures of rescue and crisis services, determination of evacuation area, informing the general public about the negative human impact on the environment, and the objective compensation payments to farmers whose crops have been damaged or destroyed. In addition to monitoring in near-real-time, satellite mission data is used for strategic emergency planning based on computer simulations of damage caused by wildfires, floods, or storms [100]. The basis for the simulation of the wildfire spread is a vegetation fuel model classified according to biomass and moisture content, while for simulation of floods and storms the basic source of data from satellite missions is a DEM.

Conclusions
A wide range of available global open data satellite missions provides an abundant data source for various land monitoring and conservation studies. The availability of completely free images and open-source GIS software for their processing ensures the basis of worldwide land monitoring and conservation research. The main limiting factor for that is currently deficient knowledge and expertise regarding remote sensing and GIS processing, as well as the specific objectives within the land conservation. Remote sensing speeds up the data collection process over large areas in land monitoring and conservation studies but field identification is still a necessity in the calibration of spectral indices and accuracy assessment of image classification results.
The most suitable satellite mission for vegetation monitoring currently is Sentinel-2, due to its relatively high spatial and temporal resolution. Remote sensing in the red-edge and near-infrared spectral bands allows the calculation of numerous specialized vegetation indices. Landsat 8 mission also enables effective monitoring of vegetation properties, especially for modeling of the evapotranspiration and water content in plants using thermal bands, which are complementary to Sentinel-2 images. The fusion of radar Sentinel-1 and multispectral Sentinel-2 images due to mutual compatibility provides great potential in mapping and managing damage caused by ecological and natural disasters. With the emergence of new Sentinel missions in the ESA Copernicus program and Landsat 9 mission by NASA, the application of global remote sensing satellite missions in land monitoring and conservation studies is expected to maintain a rapid increase in the upcoming years. Some of the most common obstacles in the use of multispectral satellite images, like cloud cover and temporal resolution, can be resolved with the integration of UAVs in land monitoring and conservation studies, which can be effectively integrated with satellite missions' data. Terrain topographic models are the basis of numerous applications in GIS-based multicriteria analyses and simulations of natural disasters in the GIS environment. There is currently a significant potential of upgrades for global DEMs in terms of spatial resolution and geometric accuracy, as some new regional DEM datasets occurred, like ESA Copernicus EU-DEM. The main reasons for that are new radar imaging technology and interferometry processing methods that allowed the creation of more accurate DEMs than SRTM and ASTER, also being based on more recent observations.
The historic research of spectral indices resulted in few hundreds of vegetation, water, soil, and built-up indices in the past decades, reaching their development peak for the current satellite missions in land monitoring and conservation studies. However, with the emergence of new multispectral satellite missions and potential observation in new spectral bands, the development of spectral indices should be further evaluated and researched. The classification methods based on the machine learning algorithms, like RF, ANN, and SVM, are currently the primary focus of the research and are highly applicable in land monitoring and conservation studies worldwide. The trend of further development of the classification methods of satellite mission data is the development of deep machine learning algorithms and automated algorithms for specialized environmental applications, allowing their widespread implementation in government environmental services. With the rapid emergence of new remote sensing data processing methods, global specifications of remote sensing data processing and analysis for land conservation remains a necessity on a global scale. This primarily refers to the GIS-based multicriteria analyses for the calculation of land suitability, as a significant part of present studies lacks all necessary processes of data standardization and weight determination. AHP remains a stable base of these analyses due to its superiority in computational efficiency compared to similar methods and is likely to remain that way in the upcoming years.