Phenology is one of the ubiquitous fingerprints of climate change on our ecosystems. Monitoring the spatiotemporal patterns of vegetation phenology is thus critical. A wide range of sensors have been used to monitor vegetation phenology. Sensor point of view and resolution can potentially impact estimates of phenology. We compared three different sensors from three different remote sensing platforms—a UAV mounted RGB camera, an under canopy, upward facing hemispherical camera with R, G and NIR capabilities, and a tower mounted RGB PhenoCam—to estimate spring phenological transition in a mixed-species temperate forest in central Virginia, USA. Our study had two objectives: 1) to compare the above- and below- canopy inference of canopy greenness (green chromatic coordinate and normalized difference vegetation index) and canopy structural attributes (leaf area and gap fraction) by matching under-canopy hemispherical photos with high spatial resolution (0.03 m) drone imagery to find the appropriate spatial coverage and resolution for comparison; 2) to compare how each sensor performed in estimating the temporality of the spring phenological transition. We find that a spatial buffer of 20 m radius for UAV imagery is most closely comparable to under-canopy imagery in this system. Sensors and platforms agree within +/- 5 days of when canopy greenness stabilizes from the spring phenophase into the growing season. This work has implications for paring UAV imagery with both tower-based observation platforms, as well as plot-based studies (e.g. long-term monitoring, existing research networks, permanent plots).

A neglected dimension of the fire regime concept is fire patchiness. Habitat mosaics that emerge from the grain of burned and unburned patches (pyrodiversity) are critical for the persistence of a diverse range of plant and animal species. This issue is of particular importance in frequently burned tropical Eucalyptus savannas, where coarse fire mosaics have been hypothesized to have caused the recent drastic population declines of small mammals. Satellites routinely used for fire mapping in these systems are unable to accurately map fine-grained fire mosaics, frustrating our ability to determine whether declines in biodiversity are associated with local pyrodiversity. To advance this problem, we have developed a novel method (we call ‘double-differenced dNBR’) that combines the infrequent (c. bi-monthly) detailed spatial resolution Landsat with daily coarse scale coverage of MODIS and VIIRS to map pyrodiversity in the savannas of Kakadu National Park. We used seasonal Landsat mosaics and differenced Normalized Burn Ratio (dNBR) to define burned areas, with a modification to dNBR that subtracts long-term average dNBR to increase contrast. Our results show this approach is effective in mapping fine-scale fire mosaics in the homogenous lowland savannas, although inappropriate for nearby heterogenous landscapes. Comparison of this methods to other fire metrics (e.g., area burned, seasonality) based on Landsat and MODIS imagery suggest this method is likely accurate and better at quantifying fine-scale patchiness of fire, albeit it demands detailed field validation.

This paper presents ensemble learning of multi-source satellite sensors dataset to obtain better predictive performance of the forest biomass. Spectral, spectral-indices, and spectral-textural features were generated from two optical satellite sensors, Landsat 8 Operational Land Imager (OLI) and Sentinel-2 Multispectral Instrument (MSI). In addition, two radar satellite sensors, Sentinel-1 C-band Synthetic Aperture Radar (CSAR), and Advanced Land Observing Satellite (ALOS-2) Phased Array type L-band Synthetic Aperture Radar (PALSAR-2) were utilized to generate backscattering and backscattering-textural features. The plot-wise above ground biomass data available from five forests in New England region were utilized. Ensemble learning of multi-source satellite sensors dataset was carried out by employing four machine learning regressors namely, Support Vector Machines (SVM), Random Forests (RF), Gradient Boosting (GB), and Multilayer Perceptron (MLP). A five-fold cross-validation method was used to evaluate predictive performance of the multi-source satellite sensors. The integration of multi-source satellite features, comprising of spectral, spectral-indices, backscattering, spectral-textural, and backscattering-textural information, through ensemble learning and cross-validation approach implemented in the research showed promising results (R² = 0.81, RMSE = 46.2 Mg/ha) for the estimation of plots-level forest biomass in New England region.

This study explores the changes in the landmass bounded by the coast of India during 1975-2005 by using on-screen visual interpretation technique (with 100m resolution and 1:50,000 scale) from NASA Landsat Imagery in three different time periods viz. 1975, 1990, and 2005. The result indicated an overall expansion of 130 sq. km area of the landmass that surrounded by the Indian coast during 1975-2005 (74 sq. km during 1975-1990 and 56 sq. km during 1991-2005). These estimations are based on the preliminary analysis and may be estimated more accurately by reducing the scale and using further higher resolution images.

Relationship between caveman prehistoric life in terms of heat induced food processing and its geological ecosystems have received many attentions. Previous studies have investigated the sources of heat included using Fourier transform infrared spectroscopy and biomarker approaches. Here this study proposes the use of remote sensing to identify the relationship of 9500 year old (9.5 ka) prehistoric mongoloid occupancy with hydrothermal manifestations at Pawon cave of West Java. The hydrothermal manifestations around Pawon cave were identified using Landsat 8 band combinations, land surface temperature, and sedimentary lithology. The results showed the hydrothermal manifestations surrounding Pawon cave were within a distance of 0.5-2 km. The results also showed bones representing 12 animal taxon groups with high abundance of rodents. To conclude this study sheds the light of proximity and preferences of mongoloid prehistoric occupancy towards hydrothermal landscape due to its advantage as heat sources for food processing purposes.

Data that have been processed to allow analysis with a minimum of additional user effort are often referred to as Analysis Ready Data (ARD). The ability to perform large scale Landsat analysis relies on the ability to access observations that are geometrically and radiometrically consistent, and have had non-target features (clouds) and poor quality observations flagged so that they can be excluded. The United States Geological Survey (USGS) has processed all of the Landsat 4 and 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper Plus (ETM+), Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) archive over the conterminous United States (CONUS), Alaska, and Hawaii, into Landsat ARD. The ARD are available to significantly reduce the burden of pre-processing on users of Landsat data. Provision of pre-prepared ARD is intended to make it easier for users to produce Landsat-based maps of land cover and land-cover change and other derived geophysical and biophysical products. The ARD are provided as tiled, georegistered, top of atmosphere and atmospherically corrected products defined in a common equal area projection, accompanied by spatially explicit quality assessment information, and

Suitable measures of grazing impacts on ground cover, that enable separation of the effects of climatic variations, are needed to inform land managers and policy makers across the arid rangelands of the Northern Territory of Australia. This work developed and tested a time-series, change-point detection method for application to time series of vegetation fractional cover derived from Landsat data to identify irregular and episodic ground-cover growth cycles. These cycles were classified to distinguish grazing impacts from that of climate variability. A measure of grazing impact was developed using a multivariate technique to quantify the rate and degree of ground cover change. The method was successful in detecting both long term (> 3 years) and short term (< 3 years) growth cycles. Growth cycle detection was assessed against rainfall surplus measures indicating a relationship with high rainfall periods. Ground cover change associated with grazing impacts was also assessed against field measurements of ground cover indicating a relationship between both field and remotely sensed ground cover. Cause and effects between grazing practices and ground cover resilience can now be explored in isolation to climatic drivers. This is important to the long term balance between ground cover utilisation and overall landscape function and resilience.

Agriculture in Morocco has been extensive until the middle of the 20th century due to the distribution of rainfall and the availability of water. In the middle of the last century hydraulic works were built that allowed the transition to intensive agriculture by the increase of irrigated areas, allowing that in the territories where there is water for irrigation and the climate allows it, the crops adapt to the demands of the market. The objective of the study is to assess by satellite images the land cover between 1985 and 2020, analyzing the changes in cultivation areas, as well as the changes in desert, sub-desert and forest areas of the Oum Er Rbia hydrological basin in Morocco. Landsat satellite images have been used since 1984 by the US government (Aerospace and Geological Agencies). A series of vegetation indices (NDVI, RVI, TNDVI and EVI) have been used; among which TNDVI (Transformed Normalized Vegetation Index) stands out for its better accuracy, which has allowed us to distinguish vegetation in cultivated and forest areas, as well as arid zones. In addition, the study has compared the use of two methodologies to calculate changes in the coverage of the Earth’s surface, has used local image processing from the Sentinel Application Platform tool and has also used the Google Earth Engine tool. The latter being the most optimal, although at the moment it has great limitations. In both methodologies and in the different indices it has been possible to observe during these 35 years as the cultivated area has increased (related to the availability of water by the construction of reservoirs and canals), how plant cover has improved in forest areas, and a range of variations in arid areas.

In this study, our goal was to research land-use change by combining spatio–temporal land use/land cover monitoring (LULC (1989–2019) and urban growth modeling (1999–2039) in Benslimane, Morocco, to determine the effect of urban growth on different groups based on cellular automata (CA) and geospatial methods. A further goal was to test the reliability of the AC algorithm for urban expansion modeling. To do this, four years of satellite data were used at the same time as population density, downtown distance, slope, and ground road distance. The LULC satellite reported a rise of 3.8 km2 (318% variation) during 1989–2019. Spatial transformation analysis reveals a good classification similarity ranging from 89% to 91% with the main component analysis (PCA) technique. The statistical accuracy between the satellite scale and the replicated built region of 2019 gave 97.23 %t of the confusion matrix overall accuracy, and the region under the receiver operational characteristics (ROC) curve to 0.94, suggesting the model's high accuracy. Although the constructed area remains low relative to the total area of the municipality's territory, the LULC project shows that the urban area will extend to 5,044 km2 in 2019, principally in the western and southwestern sections. In 2019–2039, urban development is expected to lead to a transformation of the other class (loss of 1,364 km2), followed by vegetation cover (loss of 0.345 km2). In spatial modeling and statistical calculations, the GDAL and NumPy Python 3.8 libraries were successful.

Due to increase demand of food grain in the world, assessment of yield before actual production is important in making policies and decisions in agricultural production system. For a large area, forecast models developed from vegetation indices derived from remote sensing satellite data possesses the potential to give quantitative and timely information on crops over large areas. Different vegetation indices are being made used for this purpose, however, their efficiency in estimating crop yield is needed to be certainly tested. In this study, wheat yield forecast was derived by regressing ground truthing yield data against time series of spatial vegetation indices for the 2013 to 2019 growing seasons. These spatial vegetation indices derived from Landsat 8 image data: Normalized Difference Vegetation Index (NDVI) and Soil Adjusted Vegetation Index (SAVI) were compared to evaluate the most appropriate index that performs better in forecasting wheat production at Karcag, Kunhegyes and Ecsegfalva settlements in Jász-Nagykun-Szolnok county, in the Northern Great Plain region of central Hungary. The best time for making wheat yield prediction with Landsat 8- SAVI and NDVI was found to be the beginning of ripening period (160th day of the year) with higher correlation between the vegetation indices and the wheat yield. The validation results revealed that the model from SAVI provides more consistent and accurate forecasts yield compared to NDVI. The SAVI model forecast yield for the validation years, 2018 and 2019 were within 6.00% and 4.41% of the final reported values while that of NDVI model were within 8.31% and 6.27%. Nash-Sutcliffe efficiency index is positive with E1= 0.99 for the model from SAVI and for NDVI, E1=0.57, which connote that the forecasting method developed and evaluated performs acceptable forecast efficiency.

Oil spills from offshore drilling and coastal refineries often cause degradation of coastal wetlands that can take a long time to recover. Early oil detection may prevent losses and speed up recovery if monitoring of the initial oil extent, oil impact, and recovery are in place. Satellite imagery data can provide a cost-effective alternative to expensive airborne imagery or labor intensive field campaigns for monitoring effects of oil spills on wetlands. However, these satellite data may be restricted in their ability to detect and map ecosystem recovery post-spill given their spectral measurement properties and temporal frequency. In this study, we assessed whether spatial and spectral resolution, and other sensor characteristics influence the ability to detect and map vegetation stress and die-off due to oil. We compared how well three satellite multispectral sensors: WorldView2, RapidEye and Landsat EMT+, match the ability of the airborne hyperspectral AVIRIS sensor to map oil-induced vegetation stress, recovery, and die-off after the DeepWater Horizon oil spill in the Gulf of Mexico in 2010. We found that finer spatial resolution (3.5m) provided better delineation of the oil-impacted wetlands and better detection of vegetation stress along oiled shorelines in saltmarsh wetland ecosystems. As spatial resolution become coarser (3.5m to 30m) the ability to accurately detect and map stressed vegetation decreased. Spectral resolution did improve the detection and mapping of oil-impacted wetlands but less strongly than spatial resolution, suggesting that broad-band data may be sufficient to detect and map oil-impacted wetlands. AVIRIS narrow-band data performs better detecting vegetation stress, followed by WorldView2, RapidEye and then Landsat 15m (pan sharpened) data. Higher quality sensor optics and higher signal-to-noise ratio (SNR) may also improve detection and mapping of oil-impacted wetlands; we found that resampled coarser resolution AVIRIS data with higher SNR performed better than either of the three satellite sensors. The ability to acquire imagery during certain times (midday, low tide, etc.) or a certain date is also important in these tidal wetlands; WorldView2 imagery captured at high-tide detected a narrower band of shoreline affected by oil likely because some of the impacted wetland was below the tideline. These results suggest that while multispectral data may be sufficient for detecting the extent of oil-impacted wetlands, high spectral and spatial resolution, high-quality sensor characteristics, and the ability to control time of image acquisition may improve assessment and monitoring of vegetation stress and recovery post oil spills.

Unlike the contiguous windthrows, the diffuse windthrows occurred as a result of wind gusts of lower speed (100-140 km/h) than in the first case (>140 km/h) are much more difficult to detect due to their much lower areas and due to their very large number, of several hundreds in the wooded mountain massifs. The objective of this research is to present a rapid procedure for the detection of the diffuse windthrows based on low cost, Landsat type images, knowing that certain sensors cannot be accessed without significant investments. Our application is based on the study of effects caused by the Xynthia storm in the Vosges Mountains in the North-East of France, on 28 February 2010. Thus, based on two sets of Landsat satellite images, we used the “dark object” approach and the Disturbance Index, as well as a classification of the images before and after the storm, resulting in a change map. Following the detection process, 257 scattered polygons were detected, totalling 229 ha. For validation purposes, high-resolution images and orthophotoplans taken before and after storm were used. The error matrix was calculated, achieving an overall accuracy of 86%, which confirms the quality of our analysis and supports this procedure for detecting diffuse windthrow based on low cost resources.

Glaciers in northern Pakistan are a prime source of freshwater, providing headwater in the Indus river system and serving as a lifeline to millions of people in the region. These glaciers undergo continuous changes by melting due to global warming or accumulation due to snowfall/precipitation at higher altitudes. In this study, we used remote sensing data to quantify glacier changes in spatiotemporal variability in the past three decades. Five glaciers in the Gilgit region (near the junction of the Hindukush and Karakoram Mountains) with an extent of less than 5 square kilometers were selected, namely Phakor glacier, Karamber glacier, East Gammu glacier, Bhort glacier, and Bad-e-Swat glacier. The fluctuations in these glaciers were monitored using a digital elevation model (DEM) and a cloud-free continuous series of Landsat satellite pictures from the minimal snow cover season. The annual climatic trends were studied through spatially interpolated gridded climate data WοrldClim version-1 climate database for 1970 – 2000. We used it to study the variations of minimum and maximum temperature, solar radiation, and precipitation through the preparation of sub-sets from the original global grids. Because of its exact delineation in the Gilgit sub-basin, the characterized watersheds were visually compared to optical Landsat 8 OLI data for mountainous ridge matching, revealing that SRTM 30m (radar-based) demonstrated greater accuracy than other DEMs. The temporal assessment of Bhort, Bad-e-Sawat, East Gammu, Karamber, and other rivers was also carried out. It is observed that the glaciers in the Gilgit watershed are rather stable. The little variability of glaciers is due to their geographic condition, altitude, topography, and orientation. Validation of the mapped glacier classes has been performed to check the accuracy assessment through an error matrix method. The kappa coefficient from the error matrix has been calculated to be 84 %. The study makes a critical input to a greater understanding of watershed controlling and hydrological processes in the upper Indus catchment's Gilgit watershed.

Sun glint, i.e., direct solar radiation reflected from a water surface, negatively affects the accuracy of ocean color retrieval schemes if entering the field-of-view of the observing instrument. Herein, a simple and robust method to quantify the sun glint contribution to top-of-atmosphere (TOA) reflectances in the visible (VIS) and near-infrared (NIR) is proposed, exploiting concomitant observations of the sun glint’s morphology in the shortwave infrared (SWIR) characterized by reflectance contrasts typically higher than those resulting from other in-water or atmospheric processes. The proposed method, termed Glint Removal through Contrast Minimization (GRCM), requires high spatial resolution (ca. 10–50 m) imagery to resolve the sun glint’s characteristic morphology, meeting additional criteria on radiometric resolution and temporal delay between the individual band’s acquisitions. It has been applied with good success to a selection of Landsat 8 (L8) Operational Land Imager (OLI) scenes encompassing a wide range of environmental conditions in terms of observation geometry and glint intensity, as well as aerosol and Rayleigh optical depth. The method proposed herein is entirely image based and does not require ancillary information on the sea surface roughness or related parameters (e.g., surface wind), neither the presence of clear water areas in the image under consideration. Limitations of the proposed method are discussed, and its potential for sensors other than OLI and applications beyond glint removal is sketched.

Soil salinity is a severe soil degradation problem mainly faced in arid and semi-arid regions. About 11 million ha of land in the arid, semi-arid, and desert parts of Ethiopia is salt-affected, especially in the Awash River basin, including Afambo irrigated area. Remote sensing approaches are significant tools for accurately predicting and modeling accurately predicting and modeling soil salinity in various world regions. This study aims to analyze and model soil salinity status in the case of Afambo irrigated areas using Landsat-8 and sentinel-2A, Afar region, Ethiopia, by applying remote sensing with field measurements. Thirty-two soil samples were collected from the topsoil (0-30 cm); out of these, 25 soil samples with various EC ranges were selected for modeling, and the remaining 7 samples were utilized to validate the model. Landsat-8 and Sentinel-2A images acquired in the same month were used to extract soil salinity indices. Linear regression analyses correlated the EC data with corresponding soil salinity spectral index values derived from satellite images. The best-performing model was selected for salinity mapping. The soil salinity indices extracted from both Landsat-8 and Sentinel-2A bands estimated soil salinity with high acceptable accuracy of R² values of SI, 0.78 and 0.81, respectively. The model results in three salinity classes with varying degree of salinity, namely, highly saline, moderately saline, and slightly saline, which covers 15.1%, 39.8% and 45.1% of the total area for Landsat-8, respectively and 26.1%, 32%, and 41.9% for sentinel 2A, respectively. Generally, the results revealed that the expansion rate of salt-affected soils has been increasing. From this study, it is possible to infer that if the present irrigation practice continues, it is expected that total the cultivated lands will become sterile within a short period. Thus, it needs to be monitored regularly to secure up-to-date knowledge of their extent to improve management practices and take appropriate actions.

Forests are critical in regulating the world’s climate and they maintain overall Earth’s energy balance. The variability in forest canopy structure, topography and underneath vegetation background condition creates uncertainty in the estimation and modelling of Earth’s surface radiation particularly for boreal regions in high latitude. We studied seasonal variation in surface reflectance with respect to land cover classes, canopy structures, and topography in a boreal region of Alaska by fusing together Landsat 8 surface reflectance and LiDAR-derived canopy matrices. Our study shows that canopy structure and topography interplay and influence surface reflectance in a complex way particularly during the snow season. Topographic aspect and elevation control vegetation growth, type and structure. The southern slope is featured with more deciduous and taller trees having greater rugosity than the northern slope. Higher elevation is associated with taller trees for both vegetation types, particularly in the southern slope. In general, surface reflectance shows similar relationships with canopy cover, height and rugosity, mainly due to close relationships between these parameters. Surface reflectance decreases with canopy cover, tree height, and rugosity especially for evergreen forest. Deciduous forest shows larger variability of surface reflectance, particularly in March, mainly due to the mixing effect of snow and vegetation. The relationship between vegetation structure and surface reflectance is greatly impacted by topography. The negative relationship between elevation and surface reflectance may be due to taller and denser vegetation distribution in higher elevation. Surface reflectance in the southern slope is slightly larger than the northern slope for both deciduous and evergreen forest. The shadow effect from topography and tree crowns on surface reflectance play a different role for deciduous and evergreen forests. For deciduous forest, topographic shadow effect on surface reflectance is stronger than from tree shadowing in all seasons. For evergreen forest, shadow effects from topography and tree crowns on surface reflectance are both equally dominant, however tree shadow effect is more significant in March than in May and August. The generalized additive models (GAM) based on non-linear relationships between response (surface reflectance) and predictor (canopy structures and topography) variables confirms such observations. Our study not only provides accurate quantification of surface radiation budget but also helps in parametrization of climate change models.

In general, remote sensing has proven to be a powerful tool in the overall understanding of natural and anthropogenic phenomena. Satellites have become useful tools for tasks such as characterization, monitoring, and the continuous prospecting of natural resources. This research aims to analyze spatial dynamic and destructive on coral reefs area and correlation between live coral reduction and population on small islands. Landsat MSS, TM, ETM, and OLI-TIRS are used to spatial analyze of coral reef dynamics from 1972 to 2016. The image processing includes gap-filling, atmospheric correction, geometric correction, image composite (true color), water column correction, unsupervised classification, reclassification, accuracy assessment. The statistical analysis identifies the relationship between dynamic population data with a reduction of live coral, namely Principal Component Analysis (PCA) and Multiple Regression Analysis. The effect of the population shows a positive correlation with the reduction in the area of live coral, although it is significant. The fact is the practice of coral destruction on an island; it is usually not only caused or carried out by residents who live on the island but also carried out by other residents of different islands.

In the oil palm industry, stands age is an important parameter to monitor the sustainability of cultivation, to develop the growth yield model, to identify the disease or stressed area, and to estimate the carbon storage capacity. This research is focused to estimate and distinguish oil palm stands age based on crown/ canopy density obtained using Forest Canopy Density (FCD) model derived from four indices as follows; Advanced Vegetation Index, Bare Soil Index, Shadow Index, and Thermal Index. FCD model employs multi temporal image analysis resulting four classes of oil palm stands age categorized as seed with FCD value of 29–56% (0 years), young with FCD value of 56–63% (1–9 years), teen with FCD value of 63–80% (10–15 years), and mature with FCD value of >80% (>15 years). Minimum canopy density value is 29% even in the zero years old indicates incomplete land clearance or the type of seed planted in the land.

Many remote sensing analytical data products are most useful when they are in an appropriate regional or national projection, rather than globally based projections like Universal Transverse Mercator (UTM) or geographic coordinates, i.e., latitude and longitude. Furthermore, leaving data in the global systems can create problems, either due to misprojection of imagery because of UTM zone boundaries, or because said projections are not optimised for local use. We developed the open-source Irish Earth Observation (IEO) Python module to maintain a local remote sensing data library for Ireland. This pure Python module, in conjunction with the IEOtools Python scripts, utilises the Geospatial Data Abstraction Library (GDAL) for its geoprocessing functionality. At present, the module supports only Landsat TM/ETM+/OLI/TIRS data that have been corrected to surface reflectance using the USGS/ESPA LEDAPS/ LaSRC Collection 1 architecture. This module and the IEOtools catalogue available Landsat data from the USGS/EROS archive, and includes functions for the importation of imagery into a defined local projection and calculation of cloud-free vegetation indices. While this module is distributed with default values and data for Ireland, it can be adapted for other regions with simple modifications to the configuration files and geospatial data sets.

The Operational Remote Sensing (ORS) program leverages Landsat and MODIS data to detect forest disturbances across the conterminous United States (CONUS). The ORS program was initiated in 2014 as a collaboration between the US Department of Agriculture Forest Service Geospatial Technology and Applications Center (GTAC) and the Forest Health Assessment and Applied Sciences Team (FHAAST). The goal of the ORS program is to supplement the Insect and Disease Survey (IDS) and MODIS Real-Time Forest Disturbance (RTFD) programs with imagery-derived forest disturbance data that can be used to augment traditional IDS data. We developed three algorithms and produced ORS forest change products using both Landsat and MODIS data. These were assessed over Southern New England and the Rio Grande National Forest. Reference data were acquired using TimeSync to conduct an independent accuracy assessment of IDS, RTFD, and ORS products. Overall accuracy for all products ranged from 77.64% to 93.51% (kappa 0.09–0.59) in the Southern New England study area and 59.57% to 79.57% (kappa 0.09–0.45) in the Rio Grande National Forest study area. In general, ORS products met or exceeded the overall accuracy and kappa of IDS and RTFD products. This demonstrates the current implementation of ORS is sufficient to provide data to augment IDS data.

Soil salinity, whether natural or human induced, is a major geo-hazard in arid and semi-arid landscapes. In agricultural lands, it negatively affects plant growth, crop yields, whereas in semi-arid and arid non-agricultural areas it affects urban structures due to subsidence, corrosion and ground water quality, leading to further soil erosion and land degradation Accurately mapping soil salinity through remote sensing techniques has been an active area of research in the past few decades particularly for agricultural lands. Most of this research has focused on the utilization and development of salinity indices for properly mapping and identifying saline agricultural soils. This research study develops a soil salinity index and model using Landsat 8 OLI image data from the near infra-red and shortwave infra-red spectral information with emphasis on soil salinity mapping and assessment in non-agricultural desert arid and semi-arid surfaces. The developed index when integrated into a semi-empirical model outperformed in its soil salinity mapping overall accuracy (60%) in comparison to other salinity indices (~50%). The newly developed index further outperformed other indices in its accuracy in mapping and identifying high saline soils (67%) and excessively high saline soils (90%).

The Brazilian livestock is predominantly extensive, with approximately 90% of the production being sustained on pasture, which occupies around 20% of the territory. In the current climate change scenario and where cropland is becoming a limited resource, there is a growing need for a more efficient land use and occupation. It is estimated that more than half of the Brazilian pastures have some level of degradation; however there is still no mapping of the quality of pastures on a national scale. In this study, we mapped and evaluated the spatio-temporal dynamics of pasture quality in Brazil, between 2010 and 2018, considering three classes of degradation: Absent (D0), Intermediate (D1), and Severe (D2). There was no variation in the total area occupied by pastures in the evaluated period, in spite of the accentuated spatial dynamics, with a retraction in the center-south and expansion to the north, over areas of ​​native vegetation. The percentage of non-degraded pastures increased ~12%, due to the recovery of degraded areas and the emergence of new pasture areas as a result of the prevailing spatial dynamics. However, about 44 Mha of the pasture area is currently severely degraded. The dynamics in pasture quality were not homogeneous in property size classes. We observed that in the approximately 2.68 million properties with livestock activity, the proportion with quality gains was twice as low in small properties compared to large ones, and the proportion with losses was three times greater, showing an increase in inequality between properties with more and less resources (large and small, respectively). The areas occupied by pastures in Brazil present an unique opportunity to increase livestock production and make available areas for agriculture, without the need for new deforestation in the coming decades.

Proper satellite-based crop monitoring applications at the farm-level often require near-daily imagery at medium to high spatial resolution. The synthesizing of ongoing satellite missions by ESA (Sentinel 2) and NASA (Landsat7/8) provides this unprecedented opportunity at a global scale; nonetheless, this is rarely implemented because these procedures are data demanding and computationally intensive. This study developed a complete stream processing in the Google Earth Engine cloud platform to generate harmonized surface reflectance images of Landsat7,8 and Sentinel 2 missions. The harmonized images were generated for two agriculture schemes in Bekaa (Lebanon) and Ninh Thuan (Vietnam) during the period 2018-2019. We evaluated the performance of several pre-processing steps needed for the harmonization including image co-registration, brdf correction, topographic correction, and band adjustment. This study found that the miss-registration between Landsat 8 and Sentinel 2 images, varied from 10 meters in Ninh Thuan, Vietnam to 32 meters in Bekaa, Lebanon, and if not treated, posed a great impact on the quality of the harmonized dataset. Analysis of a pair overlapped L8-S2 images over the Bekaa region showed that after the harmonization, all band-to-band spatial correlations were greatly improved from (0.57, 0.64, 0.67, 0.75, 0.76, 0.75, 0.79) to (0.87, 0.91, 0.92, 0.94, 0.97, 0.97, 0.96) in bands (blue, green, red, nir,swir1,swir2, ndvi) respectively. We demonstrated that dense observation of the harmonized dataset can be very helpful for characterizing cropland in highly dynamic areas. We detected unimodal, bimodal and trimodal shapes in the temporal NDVI patterns (likely cycles of paddy rice) in Ninh Thuan province only during the year 2018. We fitted the temporal signatures of the NDVI time series using harmonic (Fourier) analysis. Derived phase (angle from the starting point to the cycle's peak) and amplitude (the cycle's height) were combined with max-NDVI to generate an R-G-B image. This image highlighted croplands as colored pixels (high phase and amplitude) and other types of land as grey/dark pixels (low phase/amplitude). Generated harmonized datasets that contain surface reflectance images (bands blue, green, red, nir, swir1, swir2, and ndvi at 30 meters) over the two studied sites are provided for public usage and testing.

While deforestation rates decline globally they are rising in the Western Amazon. Artisanal-scale gold mining (ASGM) is a large cause of this deforestation and brings with it extensive environmental, social, governance, and public health impacts, including large carbon emissions and mercury pollution. Underlying ASGM is a broad network of factors that influence its growth, distribution, and practices such as poverty, flows of legal and illegal capital, conflicting governance, and global economic trends. Despite its central role in land use and land cover change in the Western Amazon and the severity of its social and environmental impacts, it is relatively poorly studied. While ASGM in Southeastern Peru has been quantified previously, doing so is difficult due to the heterogeneous nature of the resulting landscape. Using a novel approach to classify mining that relies on a fusion of CLASlite and the Global Forest Change dataset, two Landsat-based deforestation detection tools, we sought to quantify ASGM-caused deforestation in the period 1984–2017 in the southern Peruvian Amazon and examine trends in the geography, methods, and impacts of ASGM across that time. We identify nearly 100,000 ha of deforestation due to ASGM in the 34-year study period, an increase of 21% compared to previous estimates. Further, we find that 10% of that deforestation occurred in 2017, the highest annual amount of deforestation in the study period, with 53% occurring since 2011. Finally, we demonstrate that not all mining is created equal by examining key patterns and changes in ASGM activity and techniques through time and space. We discuss their connections with, and impacts on, socio-economic factors, such as land tenure, infrastructure, international markets, governance efforts, and social and environmental impacts.

Remote sensing datasets are increasingly being used to provide spatially explicit large scale evapotranspiration (ET) estimates. The focus of this study was to estimate and thematically map on a pixel-by-pixel basis, the actual evapotranspiration (ET_a) of the Wonji Shoa Sugarcane Estate using the Surface Energy Balance Algorithm for Land (SEBAL), Simplified Surface Energy Balance (SSEB) and Operational Simplified Surface Energy Balance (SSEBop) algorithms. The results obtained revealed that the ranges of the daily ETa estimated on January 25, February 26, September 06 and October 08, 2002 using SEBAL were 0.0 - 6.85, 0.0 – 9.36, 0.0 – 3.61, 0.0 – 6.83 mm/day; using SSEB 0.0 - 6.78, 0.0 – 7.81, 0.0 – 3.65, 0.0 – 6.46 mm/day, and SSEBop were 0.05 - 8.25, 0.0 – 8.82, 0.2 – 4.0, 0.0 – 7.40 mm/day, respectively. The Root Mean Square Error (RMSE) values between SSEB and SEBAL, SSEBop and SEBAL, and SSEB and SSEBop were 0.548, 0.548, and 0.99 for January 25, 2002; 0.739, 0.753, and 0.994 for February 26, 2002;0.847, 0.846, and 0.999 for September 06, 2002; 0.573, 0.573, and 1.00 for October 08, 2002, respectively. The standard deviation of ET_a over the sugarcane estate showed high spatio-temporal variability perhaps due to soil moisture variability and surface cover. The three algorithm results showed that well watered sugarcane fields in the mid-season growing stage of the crop had higher ET_a values compared with the other dry agricultural fields confirming that they consumptively use more water. Generally during the dry season, ET_a is limited to water surplus areas only and in wet season, ET_a was high throughout the entire sugarcane estate. The evaporation fraction (ETrF) results also followed the same pattern as the daily ET_a over the sugarcane estate. The total crop and irrigation water requirement and effective rainfall estimated using the Cropwat model were 2468.8, 2061.6 and 423.8 mm/yr for January 2001 planted and 2281.9, 1851.0 and 437.8 mm/yr for March 2001 planted sugarcanes, respectively. The mean annual ET_a estimated for the whole estate were 107 Mm³, 140 Mm³, and 178 Mm³ using SEBAL, SSEB, and SSEBop, respectively. Even though the algorithms should be validated through field observation, they have potential to be used for effective estimation of ET in the sugarcane estate.

VNREDSat-1 is the first Vietnamese satellite allowing the survey of environmental parameters such as vegetation and water coverages, or surface water quality at medium spatial resolution (from 2.5 to 10 meters depending on the considered channel). The NAOMI sensor on board VNREDSat-1 has the required spectral bands to assess the suspended particulate matter concentration, SPM. Because recent studies have shown that the remote sensing reflectance, Rrs(), at the blue (450 – 520 nm), green (530 – 600 nm), and red (620 – 690 nm) spectral bands can be assessed from NAOMI with a good accuracy, the present study is dedicated to the development and validation of an algorithm (hereafter referred to as V1SPM) to assess SPM from Rrs() over inland and coastal waters of Vietnam. For that purpose, an in situ data set of hyper-spectral Rrs() and SPM (from 0.47 to 240.14 g.m-3) measurements collected at 205 coastal and inland stations has been gathered. Among the different approaches, including 4 historical algorithms, the polynomial algorithms involving the red-to-green reflectance ratio presents the best performance on the validation data set (MAPD of 18,7%). Compared to the use of a single spectral band, the band ratio allows to reduce the scatter around the polynomial fit as well as the impact of imperfect atmospheric corrections. Due to the lack of matchup data points with VNREDSat-1, the full VNREDSat-1 processing chain (RED-NIR and V1SPM) aiming at estimate SPM from the top-of-atmosphere signal has been applied to the Landsat-8/OLI match-up data points with relatively low to moderate SPM concentration (3.33-15.25 g.m-3) showing a MAPD of 15,8%. An illustration of the use of this VNREDSat-1 processing chain during a flooding event occurring in Vietnam is provided.

A remote sensing-based evapotranspiration (ET) study was conducted over the Central Arizona Irrigation and Drainage District (CAIDD), an Arizona agricultural region. ET was assessed means for 137 wheat plots, 183 cotton plots, and 225 alfalfa plots. The remote sensing ET models were the Satellite-Based Energy Balance for Mapping Evapotranspiration with Internalized Calibration (METRIC), the Two Source Energy Balance (TSEB), and Vegetation Index ET for the US Southwest (VISW). Remote sensing data were principally Landsat 5, supplemented by Landsat 7, MODIS Terra, MODIS Aqua, and ASTER. The models produced similar daily ET for wheat, with 6–8 mm/d mid-season. For cotton and alfalfa daily ET showed greater differences, where TSEB produced largest daily ET, METRIC the least, and VISW in the midrange. Modeled cotton ET at mid-season ranged from 9.5 mm/d (TSEB), to 8 mm/d (VISW), and 6 mm/d (METRIC). For alfalfa ET, values at peak cover ranged from 8 mm/d (TSEB), 6 mm/d (VISW), and 5 mm/d (METRIC). Model bias ranged −10% to +18%. Relative to potential ET, FAO-56 ET, and USDA-SW gravimetric-ET, model variability ranged from negligible to 35% of annual crop water use. Model averaging was found a useful way to consider and reconcile all ET estimates.

Reliable evapotranspiration (ET) estimation is a key factor for water resources planning, attaining sustainable water resources use, irrigation water management, and water regulation. During the past few decades, researchers have developed a variety of remote sensing techniques to estimate ET. The Earth Engine Evapotranspiration Flux (EEFlux) application uses Landsat imagery archives on the Google Earth Engine platform to calculate the daily evapotranspiration at the local field scale (30 m). Automatically calibrated for each Landsat image, the EEFlux application design is based on the widely vetted Mapping Evapotranspiration at high Resolution with Internalized Calibration (METRIC) model and produces ET estimation maps for any Landsat 5, 7 or 8 scene in a matter of seconds. In this research we evaluate the consistency and accuracy of EEFlux products that are produced when standard US and global assets are used. Processed METRIC products for 58 scenes distributed around the western and central United States were used as the baseline for comparison. The goal of this paper is to compare the results from EEFlux with the standard METRIC applications to illustrate the utility of the EEFlux products as they currently stand. Given that EEFlux is derived from METRIC, differences are expected to occur due to differing calibration methods (automatic versus manual) and differing input datasets. The products compared include the fraction of reference ET (ETrF), actual ET (ETa), and surface energy balance components net radiation (Rn), ground heat flux (G), and sensible heat flux (H), as well as Ts, albedo and NDVI. The product comparisons show that the intermediate products of Ts, Albedo, and NDVI, and also Rn have similar values and behavior for both EEFlux and METRIC. Larger differences were found for H and G. Despite the more significant differences in H and G, results show that EEFlux is able to calculate ETrF and ETa values comparable to the values from trained expert METRIC users for agricultural areas. For non-agricultural areas such as semi-arid rangeland and forests, the automated EEFlux calibration algorithm needs to be improved in order to be able to reproduce ETrF and ETa that is similar to the manually calibrated METRIC products.

Relationships among the population abundance of four predator groups for rice insect pests, namely: carabid beetles, staphylinid beetles, green mirid bugs, and spiders in three landscape categories were evaluated. Both rice plots and the associated bund margins of these rice plots found among three Bangladesh landscape categories were sampled by sweep net. The results revealed that the abundance significantly varied across landscapes. The rice landscape of one location harbored higher numbers of a specific predator than other location in other regions of Bangladesh. The results also showed a dependency on the width of the rice bund margins of the rice plots, where spiders populations increased with increased bund widths, but the population abundance of these predators did not depend on the diversity of the number of weed species found on the rice bund margins. The relative abundance of predator populations also significantly differed among the three landscapes, with the green mirid bug having the highest number among the four predators. This study indicates that predators of rice insect pests are highly landscape specific. In order to design integrated pest management systems for different Bangladeshi rice production locales, considerations unique to the characteristics of each locale are necessary. Preliminary efforts to apply variography analyses to the RED spectral band of LANDSAT 8 imagery from December 2016 are presented as first step toward learning a suite of methods which describe useful local characteristics affecting rice pest predators.

The monitoring of lake waterbody area in the Qinghai Tibet Plateau (QTP) is of great significance to deal with global climate change. As the latest generation of Landsat series satellites, Landsat-9 data not only have higher radiometric resolution, but also cooperate with other Landsat satellites to greatly improve the temporal resolution. It has great application potential in lake waterbody area monitoring. In order to explore the performance of different algorithms for extracting waterbody and lake waterbody area in Landsat-9 data under large-scale QTP regions, this study relies on Google Earth Engine (GEE) platform and selects 10 waterbody extraction algorithms as the basis to realize the quantitative evaluation of QTP lake waterbody area extraction results. The results show that the Random Forest (RF) algorithm performs best in all models. The overall accuracy of waterbody extraction is 95.84%, and the average error of lake waterbody area extraction is 1.505%. Among the traditional threshold segmentation waterbody extraction algorithms, the overall accuracy of the NDWI waterbody extraction method is 89.89%, and the average error of lake waterbody area extraction is 3.501%, which is the highest performance model in this kind of algorithms. This study proves that Landsat-9 data can effectively classify QTP waterbodies. With the development of cloud computing technologies such as Gee, more complex models such as RF can be selected to improve the extraction accuracy of water body and Lake area in large-scale research.

Urban cities are the major drivers of economic growth and development. Economic growth and development however results in considerable land cover land use dynamics. This study assessed the dynamics in land cover land use that have occurred in New Braunfels, Texas in the last 7 years (2013 - 2020) to observe areas in the city that had experienced considerable shifts in land cover and land use. A 30-meter resolution Landsat images were used to examine possible changes in land cover land use. New Braunfels was observed to have experienced significant changes in land use especially in developed areas. This change can be attributed to the influx of people into the city, contributing to the need for increased urban development. Analysis of this study shows that about 16% (about 553 hectares) of forest land cover class and 28% (about 1,139 hectares) of grassland class in time 1 (August 31, 2013) changed to built-up land use class in time 2 (November 5, 2020). A limitation to this study was that of the spatial resolution of images used. Higher spatial resolution images could impact the producers, users, and overall accuracy assessment. Results from this study can aid in supporting better decision-making for sustainable urban development and climate change mitigation.

LULC changes are major environmental challenges in many parts of the world which are adversely affecting ecosystem services. This study was aimed to analyze LULC changes in the ecological landscape of Ethiopia CRV areas from 1985 to 2015. Satellite images were accessed and pre-processing and classification is done. Major LULC types were detected and change analysis was executed. Nine LULC changes were successfully evaluated. The classification result revealed that in 1985, 44.34% of the land was covered with small scale farming followed by mixed cultivated/acacia (21.89%), open woodland (11.96%), and water bodies (9.77%). Whereas for the same study year open grazing land, forest, degraded savannah and settlements accounted the smallest proportion. Though the area varied among land use classes, the trend of share occupied by the LULC types in the study area remained the same in 1995 and 2015. Increase in small and large scale farming, settlements and mixed cultivation/acacia while a decrease in water bodies, forest, and open woodlands is noted. About 86.11% of the land showed major changes in land use/cover. Lastly, DPSIR framework analysis was done and integrated land use and development planning and policy reform are suggested for sustainable land use planning and management.

Land Use Land Cover (LULC) change and urban growth have a significant influence on local climate of cities. From 1985 to 2021 the population of Baghdad increased by 103%. Therefore, the risen question is how this expansion influences the temperature of the city. The study aims to identify urban growth of Baghdad, investigate its influence on variation of Land Surface Temperature (LST) and identify the main factors that control the surface temperature of the city. Three Landsat images from 1985 to 2021, in addition to sixteen potential factors, were used in the study. Our findings suggest that during the study period, vegetated areas declined by 39% while built-up class increased by 139%. Bare soil recorded the highest surface temperature. The study found that surface temperature has a strong inverse relationship with vegetation (Normalized Difference Vegetation Index (NDVI): r = -0.62, p < 0.001) and moisture (Normalized Difference Moisture Index (NDMI): r = -0.65, p < 0.001). Therefore, increasing vegetation and water body lead to decrease temperature of the city. Our findings help policymakers to deal with climatic issues rising from urban growth of the city.

Impervious surface area (ISA) is a crucial indicator for quantitative urban studies. It is also important for land use land cover classification, groundwater recharge, sustainable development, urban heat island effects, and more. Spectral ISA mapping suffers from mixed pixel problems, especially with bare soil. This study aims to develop an ISA index for spatiotemporal urban mapping from common multispectral bands by reducing soil signature better than in previous studies. This study proposed a global impervious surface area index (GISAI) enhancing ISA mapping accuracy using a temporal parameter of the remote sensing (RS) dataset. Bare soil spectral reflectance shows more fluctuation than urban ISA. Therefore, the study uses minimum composites of earlier urban indices to compile minimum soil signature. It is later improved by removing water, increasing the contrast between bare soil and urban ISA and reducing bright bare soil area. This study maps the ISA of all 12 megacities using the annual RS image collection from 2021. GISAI reduced the bare soil signature and achieved an overall accuracy of 87.29%, F1-score of 0.84, and Kappa coefficient of 0.75. However, it has some limitations with grey bare soil and barren hilly areas. By limiting bare soil signature, GISAI broadens the scope of inter-urban studies globally and lengthens potential urban time-series analysis using common multispectral bands.

Land use land cover (LULC) changes are highly pronounced in African countries, as they are characterized by an agriculture-based economy and a rapidly growing population. Understanding how land use/cover change (LULCC) influence watershed hydrology will enable local governments and policymakers to formulate and implement effective and appropriate response strategies to minimize the undesirable effects of future land use/cover change or modification and sustain the local socio-economic situation. The hydrological response of the Ethiopia Fincha’a watershed to LULCC happened during the last 30 years was investigated comparing the situation in three reference years: 1994, 2004 and 2018. The information was derived from Landsat sensors, respectively Landsat 5 TM, Landsat 7 ETM and Landsat 8 OLI/TIRS. The various LULC classes were derived via ArcGIS using a supervised classification system, and the accuracy assessment was done using confusion matrixes. For all the years investigated the overall accuracies and the kappa coefficients were higher than 80%, with 2018 as the more accurate year. The analysis of LULCC revealed that forest decreased by 19.99% between the years 1994-2004, and it decreased by 11.85% in the following period 2004-2018. Such decline in areas covered by forest is correlated to an expansion of cultivated land by 16.4% and 10.81%, respectively. After having evaluated the LULCC at the basin scale, the watershed was divided into 18 sub-watersheds, which contained 176 Hydrologic Response Units (HRUs), having a specific LULC. Accounting for such a detailed subdivision of the Fincha’a watershed, the SWAT model was firstly calibrated and validated on past data, and then applied to infer information on the hydrological response of each HRU on LULCC. The modelling results pointed out a general increase of average water flow, both during dry and wet periods, as a consequence of a shift of land coverage from forest and grass towards settlements and build-up areas. The present analysis pointed out the need of accounting for past and future LULCC in modelling the hydrological responses of rivers at the watershed scale.

Spatio-temporal Land-Use and Land-Cover (LULC) changes have been affecting geo-environmental and climate change globally. This study aims to analyze LULC changes in Bahir Dar city and its surrounds. Landsat 5 TM (1987), Landsat 7 ETM+ (2002) and Landsat 8 OLI (2017) and SPOT images, and aerial photographs, master plan map and Google Earth Landsat images were used to analyze changes. In Bahir Dar city and its surrounds, LULC has been changing in space and time. During 1987-2017, more than 50% of the study area was covered with cropland. Settlement areas have increased from 3.3% in 1987 to 9.13% in 2017. However, wetland vegetation, shrubland, grassland, forest, and waterbodies have degraded. These changes are mainly attributed to population growth and its effect on the environment. Land-use and land-cover is a serious problem and it causes land and environmental degradation, climate change and loss of the biological environment.

If given the correct remotely sensed information, machine learning can accurately describe soil moisture conditions in a heterogeneous region at the large scale based on soil moisture readings at the small scale through rule transference across scale. This paper reviews an approach to increase soil moisture resolution over a sample region over Australia using the Soil Moisture Active Passive (SMAP) sensor and Landsat 8 only and a validation experiment using Sentinal-2 and the Advanced Microwave Scanning Radiometer (AMSR-E) over Nevada. This approach uses an inductive localized approach, replacing the need to obtain a deterministic model in favor of a learning model. This model is adaptable to heterogeneous conditions within a single scene unlike traditional polynomial fitting models and has fixed variables unlike most learning models. For the purposes of this analysis, the SMAP 36 km soil moisture product is considered fully valid and accurate. Landsat bands coinciding in collection date with a SMAP capture are down sampled to match the resolution of the SMAP product. A series of indices describing the Soil-Vegetation-Atmosphere Triangle (SVAT) relationship are then produced, including two novel variables, using the down sampled Landsat bands. These indices are then related to the local coincident SMAP values to identify a series of rules or trees to identify the local rules defining the relationship between soil moisture and the indices. The defined rules are then applied to the Landsat image in the native Landsat resolution to determine local soil moisture. Ground truth comparison is done via a series of grids using point soil moisture samples and air-borne L-band Multibeam Radiometer (PLMR) observations done under the SMAPEx-5 campaign. This paper uses a random forest due to its highly accurate learning against local ground truth data yet easily understandable rules. The predictive power of the inferred learning soil moisture algorithm did well with a mean absolute error of 0.054 over an airborne L-band retrieved surface over the same region.

There is a growing need for an area-wide knowledge of SOC contents in agricultural soils at field scale for food security, monitoring long-term changes related to soil health and climate change. In Germany, large-scale SOC maps are mostly available with a spatial resolution of 250 m to 1 km2. The nationwide availability of both digital elevation models at various spatial resolutions and multi-temporal satellite imagery enables the derivation of multi-scale terrain attributes and Landsat-based multi-temporal soil reflectance composites (SRC) as explanatory variables. On the example of an Bavarian test of about 8000 km2, the scale-specific dependencies between the representativeness of 220 soil samples and different aggregation levels of the explanatory variables were analyzed for their scale-specific predictive power. The aggregation levels were generated by applying a region-growing segmentation procedure, the SOC content prediction was realized by the Random Forest algorithm. In doing so, established approaches of (geographic) object-based image analysis (GEOBIA) and machine learning were combined. The modeling results revealed scale-specific differences. Compared to terrain attributes, the use of SRC parameters lead to a significant model improvement at large field-related scale levels. The joint use of both terrain attributes and SRC parameters resulted in further model improvements. The best modeling variant is characterized by an accuracy of R2=0.84 and RMSE=1.99.

Accurate inundation maps for flooded wetlands and rivers are a critical resource for their management and conservation. In this paper we automate a method (thresholding of the short-wave infrared band) for classifying inundation, using Landsat imagery and Google Earth Engine. We demonstrate the method in the Okavango Delta, northern Botswana, a complex case study due to the spectral overlap between inundated areas covered with aquatic vegetation and dryland vegetation classes on satellite imagery. Inundation classifications in the Okavango Delta have predominately been implemented on broad spatial resolution images. We present the longest time series to date (1990-2019) of inundation maps at high spatial resolution (30m) for the Okavango Delta. We validated the maps using image-based and in situ data accuracy assessments, with accuracy ranging from 91.5 - 98.1%. Use of Landsat imagery resulted in consistently lower estimates of inundation extent than previous studies, likely due to the increased number of mixed pixels that occur when using broad spatial resolution imagery, which can lead to overestimations of the size of inundated areas. We provide the inundation maps and Google Earth Engine code for public use.

To address three important issues related to extraction of water features from Landsat imagery, i.e., selection of water indexes and classification algorithms for image classification, collection of ground truth data for accuracy assessment, this study applied four sets (ultra-blue, blue, green, and red light based) of water indexes (NWDI, MNDWI, MNDWI2, AWEI_ns, and AWEI_s) combined with three types of image classification methods (zero-water index threshold, Otsu, and kNN) to 24 selected lakes across the globe to extract water features from Landsat-8 OLI imagery. 1440 (4x5x3x24) image classification results were compared with the extracted water features from high resolution Google Earth images with the same (or ±1 day) acquisition dates through computing the Kappa coefficients. Results show the kNN method is better than the Otsu method, and the Otsu method is better than the zero-water index threshold method. If the computational cost is not an issue, the kNN method combined with the ultra-blue light based AWEI_ns is the best method for extracting water features from Landsat imagery because it produced the highest Kappa coefficients. If the computational cost is taken into account, the Otsu method is a good choice. AWEI_ns and AWEI_s are better than NDWI, MNDWI and MNDWI2. AWEI_ns works better than AWEI_s under the Otsu method, and the average rank of the image classification accuracy from high to low is the ultra-blue, blue, green, and red light-based AWEI_ns.

An oceanic radar altimeter such as TOPEX/Poseidon (T/P) is typically for observing elevation changes over the open oceans or large inland lakes/rivers, with limited applications over solid earth due to its large footprint and susceptibility to waveform contamination and slope effect. Here we demonstrate that it is possible to construct a long-term time series of glacier elevation change from T/P-series radar altimeters over two flat surfaces near a glacier terminus and an icefield (Sites A and B, with slopes of 2° and 0.8°) in Mt. Tanggula, Tibet, at elevations over 5400 m. We retracked radar waveforms using the subwaveform threshold algorithm, selected quality altimeter data (1/4 of the original) with nearly the same slope and adjusted the original elevations by fitting with a time-varying, 2nd order surface. The glacier elevation changes at the two sites from T/P (1993–2002) show seasonal elevation oscillations with linear rates at about −3 m/year and abnormal seasonal changes around the 1997–98 El Niño. Site A is over a deep valley in southern Tanggula. Its elevation dropped about 30 m over 1993–2002 (from T/P) and the glacier almost disappeared by 2016 (from altimeters and satellite images). Despite the sporadic Jason-2 and Jason-3 altimeter data, we also derived long-term rates of glacier elevation change over 1993–2017. Landsat-derived glacier area and elevation changes near the two sites confirm the rapid glacier thinning from the altimeters. The glacier meltwater near Site A supplied increasing source water to Chibuzhang Co west of Mt. Tanggula, contributing partially to its accelerated rising lake level. The glacier thinning at Site B (icefield) was correlated with the increased discharge of the Tuotuo River in eastern Mt. Tanggula, a source region of the Yangtze River. The successful detection of glacier thinning at the two sites shows that T/P-series altimeters can serve as a virtual station at a flat glacier spot to monitor long-term glacier elevation changes in connection to climate change. This virtual station concept is particularly useful for inaccessible glaciers, but its implementation faces two challenging issues: increasing the volume of quality altimeter data and improving the ranging accuracy over a targeted mountain glacier spot.

The rapid urbanization that takes place in developing regions such as Sub-Saharan Africa is associated with a large range of environmental and social issues. In this context, remote sensing is essential to provide accurate and up-to-date spatial information to support risk assessment and decision making. However, mapping urban areas remains a challenge because of their heterogeneity, especially in developing regions where the highest rates of misclassification are observed. Nevertheless, urban areas located in arid climates --- which are among the most vulnerables to anthropogenic impacts, suffer from the spectral confusion occurring between built-up and bare soil areas when using optical imagery. Today, the increasing availability of satellite imagery from multiple sensors allow to tackle the aforementioned issues by combining optical data with Synthetic Aperture Radar (SAR). In this paper, we assess the complementarity of the Landsat 8 and Sentinel-1 sensors to map built-up areas in twelve Sub-Saharan African urban areas, using a pixel-level supervised classification based on the Random Forest classifier. We make use of textural information extracted from SAR backscattering data in order to reduce the speckle noise and to introduce contextual information at the pixel level. Results suggest that combining both optical and SAR features consistently improves classification performances, mainly by enhancing the differentiation between built-up and bare lands. However, the fusion was less beneficial in mountainous case studies, suggesting that including features derived from a Digital Elevation Model (DEM) could improve the reliability of the proposed approach. As suggested by previous studies, combining features computed from both VV and VH polarizations consistently led to better classification performances. On the contrary, introducing textures computed from different spatial scales did not improve the classification performances.

The Global Human Settlement Layer (GHSL) produces new global spatial information, evidence-based analytics and knowledge describing the human presence on the planet based mainly on two quantitative factors: i) the spatial distribution (density) of built-up structures and ii) the spatial distribution (density) of resident people. Both factors are observed in the long-term temporal domain and per uniform surface units in order to support trends and indicators for monitoring the implementation of international framework agreements. The GHSL uses various input data including global, multi-temporal archives of fine-scale satellite imagery, census data, and volunteered geographic information. In this paper, we present the characteristics of GHSL information to demonstrate how original frameworks of data and tools rooted on Earth Observation could support Sustainable Development Goals monitoring. In particular, we demonstrate the reach of gridded, open and free, local yet globally consistent, multi-temporal data in filling the data gap for the Sustainable Development Goal 11. Our experiments produce a global estimate for the Land Use Efficiency indicator (SDG 11.3.1) for 10,000 urban centers, calculating the ratio of land consumption to population growth rate that took place between 1990 and 2015. The results of our research demonstrate that there is a potential to lift SDG 11.3.1 from a tier 2 as GHSL provides a global baseline for the essential variables called by the SDG 11.3.1 metadata.

The key anthropogenic effects on climate include the changes in land use and emission of greenhouse gases into the atmosphere. Depletion of vegetation poses serious threat that speeds the process of climate change and reduces carbon sequestration by the environment. Thus, the preservation of natural environment in urban areas is an essential component of the garden city model, proposed by Sir Ebenezer Howard in 1898, to ensure ecological balance. Recent Landsat images showed that Kumasi does not have the required percentage of green vegetation as was stipulated in the garden city model on which the city was built. It was observed that most parts of Kumasi's green vegetation have been lost to built environments. This study was conducted to assess the impact of urbanization on the garden city status and its effect on the micro-climate of the city. Significant changes in the vegetation cover of the city was evaluated from Landsat-TM imagery and analysis of a long term climatic data of Kumasi carried out over a 55-year period (1960 to 2015). It was observed that, climatic conditions have slightly changed, as mean surface temperature of has increased by 1.2 °C/ 55 years, due to the significant landuse changes from development of non-transpiring, reduced evaporative urban surfaces. However, the impact is not greatly felt due to the geographical location of the city on the globe despite the evidence of a considerable temperature change. Green vegetation conservation for the city is recommended as a top priority in future for city authorities and planners.

Effectively mapping and monitoring rubber plantation is still changing. Previous studies have explored the potential of phenology features for rubber plantation mapping through a pixel-based approach (pixel-based phenology approach). However, in fragmented mountainous Xishuangbanna, it could lead to noises and low accuracy of resultant maps. In this study, we investigated the capability of an integrated approach by combining phenology information with an object-based approach (object-based phenology approach) to map rubber plantations in Xishuangbanna. Moderate Resolution Imaging Spectroradiometer (MODIS) data were firstly used to acquire the temporal profile and phenological features of rubber plantations and natural forests, which delineates the time windows of defoliation and foliation phases. Landsat images were then used to extract a phenology algorithm comparing three different approaches: pixel-based phenology, object-based phenology, and extended object-based phenology to separate rubber plantations and natural forests. The results showed that the two object-based approaches achieved higher accuracy than the pixel-based approach, having overall accuracies of 96.4%, 97.4%, and 95.5%, respectively. This study proved the reliability of a phenology-based rubber mapping in fragmented landscapes with a distinct dry/cool season using Landsat images. This study indicated that the object-based phenology approaches can effectively improve the accuracy of the resultant maps in fragmented landscapes.

Wheat and barley are among the primary food resources of the world population; therefore, their growth and observation are essential in farms to enhance food security worldwide. On top of that, careful observation of the product is essential to find solutions for the issues faced during their production and to reduce the impacts of weather changes. With the advancement of Remote Sensing technology, the observation and estimation process has increased. In this study, numbers of spectral vegetation indices was used along with canopy biophysical properties ( LAI ) and biochemical properties (chlorophyll), there calculated from (Landsat 8 and Sentinel-2) satellite data. The wheat and barley samples were collected before were be ready for harvest, and a relation with the vegetarian indices was established using the Multi-Linear Regression module, in which the equations used in predicting the harvest were developed and used to create a graph for expected harvest. The result indicated that there is a strong relationship between the vegetation indices of Sentinel-2 and Landsat images and the actual grain yield with R² of 0.77 and 0.71, respectively. The results show that the strongest correlation is observed between the LAI data obtained from Sentinel data and cereal yield data, with an R² 0.68, and the highest correlation for the indices of Landsat images is observed in the NDWI with R² 0.59 and the lowest degree of error was in the root mean square error (RMSE) for the Sentinel-2 and Landsat 8 with 0.57 and 1.54. In addition, this study also showed that the least relationship for grain yield prediction was observed between the NDRI for Sentinel-2 (R² 0.1) and SAVI for Landsat image (R² 0.47).

Detailed built-up area information is valuable for mapping complex urban environments. Although a large number of classification algorithms about built-up areas have been developed, they are rarely tested from the perspective of feature engineering and feature learning. Therefore we launched a unique investigation to provide a full test of the OLI imagery for 15-m resolution built-up area classification in 2015, in Beijing, China. Training a classifier requires many sample points, and we propose a method based on the ESA's 38-meter global built-up area data of 2014, Open Street Map and MOD13Q1-NDVI to achieve rapid and automatic generation of a large number of sample points. Our aim is to examine the influence of a single pixel and image patch under traditional feature engineering and modern feature learning strategies. In feature engineering, we consider spectra, shape and texture as the input features, and SVM, random forest (RF) and AdaBoost as the classification algorithms. In feature learning, the convolution neural network (CNN) is used as the classification algorithm. In total, 26 built-up land cover maps were produced. Experimental results show that: (1) the approaches based on feature learning are generally better than those based on feature engineering in terms of classification accuracy, and the performance of ensemble classifiers e.g., RF, is comparable to that of CNN. Two dimensional CNN and the 7 neighborhood RF have the highest classification accuracy of nearly 91%. (2) Overall, the classification effect and accuracy based on image patches are better than those based on single pixels. The features that can highlight the information of the target category (for example, PanTex and EMBI) can help improve classification accuracy.

We aimed at analyzing Chlorophyll-a and CDOM dynamics from field measurements and at assessing the potential of multispectral satellite data for retrieving water-quality parameters in three small surface reservoirs in the Brazilian semiarid region. More specifically, this work comprises i) analysis of Chl-a and trophic dynamics; ii) characterization of CDOM; iii) estimation of Chl-a and CDOM from OLI/Landsat-8 and RapidEye imagery. The monitoring lasted 20 months within a multi-year drought, which contributed to water-quality deterioration. Chl-a and trophic state analysis showed a highly eutrophic status for the perennial reservoir during the entire study period, while the non-perennial reservoirs ranged from oligotrophic to eutrophic, with changes associated with the first events of the rainy season. CDOM characterization suggests that the perennial reservoir is mostly influenced by autochthonous sources, while allochthonous sources dominate the non-perennial ones. Spectral-group classification assigned the perennial as CDOM-moderate and highly eutrophic reservoir, whereas the non-perennial ones were assigned as CDOM-rich and oligotrophic-dystrophic reservoirs. The remote sensing initiative was partially successful: the Chl-a was best modelled using RapidEye for the perennial; whereas CDOM performed best with Landsat-8 for non-perennial reservoirs. This investigation showed high potential for retrieving water quality parameters in dry areas with small reservoirs.

Water quality is the measure of chemical, physical and biological suitability of water in relation to natural effects and intended purpose which may affect human health and aquatic life. Assessment of water quality is very essential for the management of water resources and human health. Traditionally, in-situ measurements have been used to obtain the water quality parameters of the water bodies. However, with the availability of satellite images, researchers have shown that satellite images are a reliable tool that can be used to estimate water quality. Satellite image-derived water quality parameters provide extensive spatial extent and large temporal variations when compared to traditional in situ sample collection and laboratory measurements. The present work estimated several parameters for quality of water in the Kamuzu reservoir of Lilongwe River for the 2013-2020 period using Sentinel-2 and Landsat-8 satellite images. The band ratio algorithms were used to retrieve Chlorophyll a (Chl-a), Turbidity, Total Suspended Matter (TSM), Secchi depth, Coloured Dissolved Organic Matter (CDOM), and Cyanobacteria from the reservoir. Turbidity and TSM were compared with the in-situ data collected over the same period. The comparison indicated R² of 0.9 and 0.69 for TSM and Turbidity respectively from Sentinel-2 images whereas R² of 0.56 and 0.61 was obtained using Landsat 8 images which are quite encouraging. The other set of results included the spatial distribution maps of water quality parameters using Landsat-8 and Sentinel-2 satellite data. It was observed that the spatial distribution of water quality parameters, except for CDOM and Cyanobacteria, showed very good distribution and matches with the theoretical results. However, for CDOM and Cyanobacteria, the distribution was almost similar for the entire study area and the band ratio algorithm may not be able to estimate them quite reasonably. This research reiterates the need for the use of remote sensing in estimating the water quality parameters and may be a substitute to the in-situ data, in terms of spread and frequency, which is very common to most of the water bodies, across the globe.

Satellite-based remote sensing technologies are utilized extensively to investigate urban thermal environments under rapid urban expansion. Current MODIS data is, however, unable to adequately represent the spatially detailed information because of its relatively coarser spatial resolution, while Landsat data can’t explore temporally the refined analysis due to the low temporal resolution. In order to resolve this situation, we used MODIS and Landsat data to generate “Landsat-like” data by using the flexible spatiotemporal data fusion method (FSDAF), and then studied spatiotemporal variation of land surface temperature (LST) and its driving factors. The results showed that 1) The estimated “Landsat-like” data have high precision; 2) By comparing 2013 and 2016 datasets, LST increases ranging from 1.8°C to 4°C were measurable in areas where the impervious surface area (ISA) increased, while LST decreases ranging from -3.52°C to -0.70°C were detected in areas where ISA decreased; 3) LST has a strongly negative relationship with the Normalized Difference Vegetation Index (NDVI), and a strongly positive relationship with Normalized Difference Built Index (NDBI) in summer; and 4) LST is well correlated with Building density (BD), in a complex conic mode, and LST may increase by 0.460°C to 0.786°C when BD increases by 0.1. Our findings can provide information useful for mitigating undesirable thermal conditions and for long-term urban thermal environmental management.

Based on 541 Landsat images between 1988 and 2016, the normalized difference vegetation indices (NDVIs) of the wetland vegetation at Xitugou (XTG) and Wowachi (WWC) inside the Dunhuang Yangguan National Nature Reserve (YNNR) in northwest China were calculated for assessing impacts of climate change on wetland vegetation in the YNNR. It was found that the wetland vegetation at the XTG and WWC both had shown a significant increasing trend in the past 30 years, and the increase in both annual mean temperature and peak snow depth over the Altun Mountains led to the increase of wetland vegetation. The influence of local precipitation on the XTG wetland vegetation was greater than on the WWC wetland vegetation, which demonstrates that in extremely arid regions, the major constrain to the wetland vegetation is water availability in soils which is greatly related to the surface water detention and discharge of groundwater. At both XTG and WWC, snowmelt from the Altun Mountains is the main contributor to the groundwater discharge, while local precipitation plays a less role in influencing the wetland vegetation at the WWC than at the XTG, because the wetland vegetation grows on a relatively flat terrain at the WWC, while in a stream channel at the XTG.

Both climate change and anthropogenic activities contribute to the deterioration of terrestrial water resources and ecosystems worldwide. Central Asian endorheic basins are among the most affected regions through both climate and human impacts. Here, we used a digital elevation model, digitized bathymetry maps and Landsat images to estimate the areal water cover extent and volumetric storage changes in small terminal lakes in Burabay National Nature Park (BNNP), located in Northern Central Asia (CA), for the period of 1986 to 2016. Based on the analysis of long-term climatic data from meteorological stations, short-term hydrometeorological network observations, gridded climate datasets (CRU) and global atmospheric reanalysis (ERA Interim), we have evaluated the impacts of historical climatic conditions on the water balance of BNNP lake catchments. We also discuss the future based on regional climate model projections. We attribute the overall decline of BNNP lakes to long-term deficit of water balance with lake evaporation loss exceeding precipitation inputs. Direct anthropogenic water abstraction has a minor importance in water balance. However, the changes in watersheds caused by the expansion of human settlements and roads disrupting water drainage may play a more significant role in lake water storage decline. More precise water resources assessment at the local scale will be facilitated by further development of freely available higher spatial resolution remote sensing products. In addition, the results of this work can be used for the development of lake/reservoir evaporation models driven by remote sensing and atmospheric reanalysis data without the direct use of ground observations.

Satellite derived bathymetry methods over coastal areas were born to deliver basic and useful information like bathymetry. However, the process is not straightforward, the main limitation being the need of field data. The Self-calibrated Spectral Supervised Shallow-water Modeler (4SM) method was tested to obtain coastal bathymetry without the use of any field data. Using LANDSAT-8 multispectral images from 2013 to 2016, a bathymetric time series was produced. Groundtruthed depths and an alternative method, Stumpf’s Band Ratio Algorithm, were used to verify the results. Retrieved (4SM) vs groundtruthed depths scored an average r2 (0.90), and a low error (RMSE = 1.47 m). Also 4SM showed, over the whole time series, the same average accuracy of the control method (40%). Advantages, limitations and operability under complex atmosphere and water column conditions, and high and low-albedo bottom processing capabilities of 4SM are discussed. In conclusion, the findings suggest that 4SM is equally accurate as the commonly used Stumpf’s method, the only difference being the independence of 4SM to previous field data, and the potential to deliver bottom spectral characteristics for further modelling. 4SM represents a significative advance in coastal remote sensing potential to obtain bathymetry and optical properties of the marine bottom.

This study is aimed at demonstrating application of vegetation spectral techniques for detection and monitoring of impact of oil spills on vegetation. Vegetation spectral reflectance from Landsat 8 data were used in the calculation of five vegetation indices (normalized difference vegetation index (NDVI), soil adjusted vegetation index (SAVI), adjusted resistant vegetation index 2 (ARVI2), green-infrared index (G/NIR) and green-shortwave infrared (G/SWIR) from the spill sites (SS) and non-spill (NSS) sites in 2013 (pre-oil spill), 2014 (oil spill date) and 2015 (post-oil spill) for statistical comparison. The result shows that NDVI, SAVI, ARVI2, G/NIR and G/SWIR indicated certain level difference between vegetation condition at the SS and the NSS were significant with p-value <0.5 in December 2013. In December 2014 vegetation conditions indicated higher level of significant difference between the vegetation at the SS and NSS as follows where NDVI, SAVI and ARVI2 with p-value 0.005, G/NIR - p-value 0.01 and GSWIR p-value 0.05. Similarly, in January 2015 a very significant difference with p-value <0.005. Three indices NDVI, ARVI2 and G/NIR indicated highly significant difference in vegetation conditions with p-value <0.005 between December 2013 and December 2014 at the same sites. Post—spill analysis show that NDVI and ARVI2 indicated low level of significance difference p-value <0.05 suggesting subtle change in vegetation conditions between December 2014 and January 2015. This technique is essential for real time detection, response and monitoring of oil spills from pipelines for mitigation of pollution at the affected sites in the mangrove forest.

1) Landsat operational land imager (OLI) data and consequent laboratory measurements were used to predict Chlorophyll-a (Chl-a) concentration and the trophic states for an inland lake within the East Kolkata Wetland, India; 2) The most suitable band ratio was identified by performing Pearson correlation analysis between Chl-a concentrations and possible OLI band and band ratios from the study points; 3) The results showed highest correlation coefficient from the band ratio OLI5/OLI4 with an R value of 0.85. The prediction model was then developed by applying regression analysis between the band ratio OLI5/OLI4 and Chl-a concentration of the study points. The reflectance ratios of the validation points were given as input on the prediction model and the model output was considered as predicted Chl-a values of the validation points to check the efficiency of the prediction model. The regression model between laboratory-derived Chl-a value and model-fitted Chl-a value of the validation points revealed a high correlation with an R² value of 0.78. Trophic State Index (TSI) of the lake was also calculated from laboratory-derived Chl-a value and model-fitted Chl-a value of the validation points. The study presented a high correlation of TSI determined from predicted data with TSI from laboratory reference data (R = 0.88). The TSI values of the lake ranged from 65 to 75 which indicate that the lake is appeared to be eutrophic to hypereutrophic conditions. 4) This empirical study showed that Landsat 8 OLI imagery can be effectively applied to estimate Chl-a levels and trophic states for inland lakes.

The changing of land use and land cover (LULC) are both affected by climate and human activity and affect climate, biological diversity, and human well-being. Accurate and timely information about the LULC pattern and change is crucial for land management decision-making, ecosystem monitoring, and urban planning, especially in developing economies undergoing industrialization, urbanization, and globalization. Biodiversity degradation and urban expansion in eastern China are research hot-spots. However, the influence of LULC changes on the region remains largely unexplored. Here, an object-based and multi-temporal image analysis approach was developed to detect how LULC changes during 1985-2015 in the Tiaoxi watershed (Zhejiang province, eastern China) using Landsat TM and OLI data. The main objective of this study is to improve the accuracy of unsupervised change detection from object-based and multi-temporal images. To this end, a total of seven LULC maps are generated with multi-temporal images. A random stratified sample design was used for assessing change detection accuracy. The proposed method achieved an overall accuracy of 91.86%, 92.14%, 92.00%, and 93.86% for 2000, 2005, 2010, and 2015, respectively. Nevertheless, the proposed method, in conjunction with object-oriented and multi-temporal satellite images, offers a robust and flexible approach to LULC changes mapping that helps with emergency response and government management. Urbanization and agriculture efficiency are the main reasons for LULC changes in the region. We anticipate that this freely available data will improve the modeling for surface forcing, provide evidence of changes in LULC, and inform water-management decision-making.

Remote sensing data and techniques utilized for various purposes including natural disasters such as earthquake as well as flood. The research aims to consume liberates Landsat 8 images for investigating crashed airplanes such as MH370. Overall approximately 300 Landsat images with less than 10% clouds utilized in addition processed through Google Engine Platform. Due to the materials as well as the color of airplane body different from the area which is a plane crashed there, moreover, it should be the characteristics of the plane shapefile different in terms of albedo, temperature as well as vegetation index value. The research observed Landsat 8 data as well as methods utilized in this research, especially, NDVI, albedo in addition to band 4, capable to distinguish between the plane and its surrounding green area. Therefore, our result confirms during the research period, there was no plane on the location as well as MH370 not crashed in this site.

Surface water distribution extracted from remote sensing data has been used in water resource assessment, coastal management, and environmental change studies. Traditional manual methods for extracting water bodies cannot satisfy the requirements for mass processing of remote sensing data; therefore, accurate automated extraction of such water bodies has remained a challenge. The histogram bimodal method (HBM) is a frequently used objective tool for threshold selection in image segmentation. The threshold is determined by seeking twin peaks, and the valley values between them; however, automatically calculating the threshold is difficult because complex surfaces and image noise which lead to not perfect twin peaks (single or multiple peaks). We developed an operational automated water extraction method, the modified histogram bimodal method (MHBM). The MHBM defines the threshold range of water extraction through mass static data; therefore, it does not require the identification of twin histogram peaks. It then seeks the minimum values in the threshold range to achieve automated threshold. We calibrated the MHBM for many lakes in China using Landsat 8 Operational Land Imager (OLI) images, for which the relative error (RE) and squared correlation coefficient (R2) for threshold accuracy were found to be 2.1% and 0.96, respectively. The RE and root-mean-square error (RMSE) for the area accuracy of MHBM were 0.59% and 7.4 km2. The results show that the MHBM could easily be applied to mass time-series remote sensing data to calculate water thresholds within water index images and successfully extract the spatial distribution of large water bodies automatically.

Abstract: This paper presents a comprehensive analysis of links between satellite-measured vegetation vigor and climate variables in Armenian mountain grassland ecosystems in years 1984–2018. NDVI is derived from MODIS and Landsat data, temperature and precipitation data are from meteorological stations. Two study sites were selected, representing arid and semi-arid grassland vegetation types, respectively. Various trend estimators including Mann-Kendall (MK) and derivatives were combined for vegetation change analysis at different time scales. Results suggest that temperature and precipitation had negative and positive impacts on vegetation growth, respectively, in both areas. NDVI-to-precipitation correlation was significant but with an apparent time-lag effect that was further investigated. No significant general changes were observed in vegetation along the observed period. Further comparisons between results from corrected and uncorrected data led us to conclude that MODIS and Landsat data with BRDF, topographic and atmospheric corrections applied are best suited for analyzing relationships between NDVI and climatic factors for the 2000-2018 period in grassland at a very local scale, but in the absence of correction tools and information, uncorrected data can still provide meaningful results. Future refinements will include removal of anthropogenic impact, and deeper investigation of time-lag effects of climatic factors on vegetation dynamics.

Land Surface Temperature (LST) estimation has been studied for several purposes, while the optimal method of estimating the LST has not been criticized yet. This research explores the optimum method in Land Surface Temperature (LST) estimation using LANDSAT-8 imagery data. Four different LST retrieval approaches, the Radiative Transfer Equation-based method (RTE), the Improved Mono-Window method (IMW), the Generalized Single-Channel method (GSC), and the Split-Window algorithm (SW), were calculated to present the LSTs over Buriram Town Municipality, Thailand. The calculated LSTs from these four methods were compared with the ground-based temperature data, taken on the same date and time of the employed LANDSAT-8 images. For this reason, the optimum method of the LST calculation was justified by considering the lowest normalized root means square error (NRMSE) values. As a result, the SW algorithm presents an optimum method in LST estimation. Regarding the SW, this algorithm requires not only the atmospheric profiles during satellite acquisition but also the retrieval of several coefficients. Besides, the LST retrieval method based on the SW algorithm is sensitive to water vapor content and coefficients. Although the SW algorithm is an optimum method explored in this study, it is emphasized that the adjustable values of coefficient response to the atmospheric state may be recommended. With these conditions, the SW algorithm can generate the land-surface temperature over the mixed land-use and land cover on the LANDSAT-8 images.

Farm activities continued sand winning operations and the allocation of plots of land to prospective developers in Ghana pose a serious threat to the forest covers and lifespan of the Forest and game reserves. With all the positive add ups to the country from forests, Ghana has lost more than 33.7%(equivalent to 2,500,000 hectares) of its forest, since the early 1990s between 2005 and 2010, the rate of deforestation in Ghana was estimated at 2.19% per annum; the sixth highest deforestation rate globally for that period. This shows how important forest monitoring can be to the forestry commission in Ghana. Despite the frameworks which have been developed to help Ghana to protect and restore its forest resources, inadequate monitoring systems remain a barrier to effective implementation. In this study, Google earth engine was used to map and analyze the structural changes of forest cover using JavaScript to query and compute Landsat, MODIS and NOAA AVHRR satellite imageries of the study area (Ghana) with spatial resolutions 30m, 250m and 7km respectively. A supervised classification was performed on three multi-temporal satellite imageries and a total of six major land use and land cover classes were identified and mapped. By using random Forest-classification technique, from 1985 to 2018 recorded by NOAA AVHRR, forest cover has decreased by 66% and 2000 to 2018 recorded by Landsat and MODIS 61% and 47% respectively. A decrease in the forest has been as a result of anthropogenic activities in Ghana. A change detection analysis was performed on these images and it was noted that Ghana is losing forest reserves in every 5years. Overlay of the reserved forest of the 2000 and the classified map of 2018 shows vegetation changed during 2000-2018 remarkably. Therefore, forest-related institutions like the Forestry Commission can employ and use this monitoring system on Google Earth Engine for processing satellite images particularly Landsat, MODIS and NOAA AVHRR for forest cover monitoring and analysis for fast, efficient and reliable results.

Salinity intrusion is one of the most serious consequences of climate change coupled with rising sea level that significantly affects agricultural activities in many parts of the world. This phenomenon has increasingly become more serious and frequently occurred in the Mekong Delta of Vietnam. As a result, Vietnam has been ranked among top five countries where have been devastatingly impacted by climate change, in particular, its Tra Vinh Province characterized by coastal plain and alluvial deposit. In addition, this area is of the tropical monsoon zone of long rainy season with source of salt brought from the sea by the tides and sea level rise. Regions that are contaminated by salt are located in lowland and often suffer from floods linking to tidal effects with salty water from river systems and channels. Soil salinity evaluation is critical for coastal protection, restoration, and agricultural planning since it can be considered as an agricultural indicator to evaluate quality of soil. Here, we attempt to estimate the soil salinity in Tra Vinh Province, in the Mekong Delta of Vietnam. Landsat 8 OLI images are utilized to derive indices for soil salinity evaluation including single bands, Vegetation Soil Salinity Index (VSSI), Soil Adjusted Vegetation Index (SAVI), Normalized Difference Vegetation Index (NDVI), and Normalized Difference Salinity Index (NDSI). Subsequently, satistical analysis between soil salinity, electrical conductivity (EC, dS/m), and environmental indices derived from Landsat 8 OLI image is performed. Results indicate that spectral value of Near Infrared (NIR) band and VSSI are highly correlated with EC (R2 = 0.7779 and R2 = 0.6957, respectively) in comparison with the other indices. Comparative results show that soil salinity derived from Landsat 8 is consistent with in situ data. Findings of this study demonstrate that Landsat 8 OLI images reveal a high potential for spatiotemporally monitoring the magnitude of soil salinity at the top soil layer. Outcomes of this study are useful for agricultural activities, planners, and farmers by providing the base map of soil salinity contamination for better selection of accomodating crop types to reduce economical lost in the context of climate change. Our proposed method that estimates soil salinity using satellite-derived variables can be applied in the other regions.

Many papers in recent years have been devoted to estimating the per pixel proportions of three broad classes of materials (e.g. photosynthetic vegetation, non-photosynthetic vegetation and bare soil) using data from multispectral sensors. Many of these papers use estimation methods based on the linear mixture model. Very few of these papers assess the accuracy of their estimators. I show how to produce confidence intervals (CIs) and joint confidence regions (JCRs) for the proportions associated with various linear mixture models. There are two main models, both of which assume that the coefficients in the model are non-negative. The first model assumes that the coefficients sum to 1. The second does not, but uses rescaling of the estimated coefficients to produce estimated proportions. Three variants of these two models are also analysed. JCRs are shown to be particularly informative, because they are typically better at localising the information than CIs are. The methodology is illustrated using examples from Landsat Thematic Mapper data at 1169 locations across Australia, each of which has associated field observations. There is also discussion about the extent to which the methodology can be extended to hyperspectral data.

The work consisted in identifying possible effects from heavy metals (HMs) pollution due to waste disposal activities in three potentially polluted sites located in Basilicata (Italy), where a release of pollutants with values over the thresholds allowed by the Italian legislation was detected. The potential variations in the physiological efficiency of vegetation have been analyzed through the multitemporal processing of satellite images. In detail, Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) images were used to calculate the Normalized Difference Vegetation Index (NDVI) trend over the years. Then, the multitemporal trends were analyzed using the median of Theil-Sen, a non-parametric estimator particularly suitable for the treatment of remote sensing data, being able to minimize the outlier effects due to exogenous factors. Finally, the subsequent application of the Mann-Kendall test on the trends identified by Theil-Sen slope allowed the evaluation of trends significance and, therefore, the areas characterized by the effects of contamination on vegetation. The application of the procedure to the three survey sites led to the exclusion of the presence of significant effects of HMs contamination on the vegetation surrounding the sites during the years of waste disposal activities.

In the remote sensing domain, it is crucial to complete semantic segmentation on the raster images, e.g., river, building, forest, etc, on raster images. A deep convolutional encoder--decoder (DCED) network is the state-of-the-art semantic segmentation method for remotely sensed images. However, the accuracy is still limited, since the network is not designed for remotely sensed images and the training data in this domain is deficient. In this paper, we aim to propose a novel CNN for semantic segmentation particularly for remote sensing corpora with three main contributions. First, we propose applying a recent CNN called a global convolutional network (GCN), since it can capture different resolutions by extracting multi-scale features from different stages of the network. Additionally, we further enhance the network by improving its backbone using larger numbers of layers, which is suitable for medium resolution remotely sensed images. Second, "channel attention'' is presented in our network in order to select the most discriminative filters (features). Third, "domain-specific transfer learning'' is introduced to alleviate the scarcity issue by utilizing other remotely sensed corpora with different resolutions as pre-trained data. The experiment was then conducted on two given datasets: (i) medium resolution data collected from Landsat-8 satellite and (ii) very high resolution data called the ISPRS Vaihingen Challenge Dataset. The results show that our networks outperformed DCED in terms of $F1$ for 17.48% and 2.49% on medium and very high resolution corpora, respectively.

Fire regimes in mountain landscapes of southern Europe have been shifting from their baselines due to the accumulation of fuel fostered by long-standing rural abandonment and fire exclusion policies. Understanding the role of fire on biodiversity is paramount to implement adequate management to mitigate the impacts of altered fire regimes and land abandonment on biodiversity. Here, we explored to what extent the spatiotemporal variation in burn severity has affected bird abundance of a mountain abandoned landscape located in the Atlantic-Mediterranean transition (NW Iberia). We took advantage of: (1) satellite images of Sentinel 2 and Landsat missions to compute burn severity indicators from 2010 to 2020, and (2) standardized bird surveys carried out over 206 point-counts along the breeding season of 2021. Bird abundance models were built from burn severity metrics together with well-known fire regime attributes (% of burnt area and time since fire). Our results showed that the spatiotemporal variation of burn severity significantly correlated with the abundance of the 39% of the modeled species, supporting the role of pyro-diversity in driving bird populations in our region. The burnt area also explained abundance patterns for 28% of species. Time since fire only correlated with the abundance of 3 species. Our findings confirm the importance of incorporating burn severity indicators into the toolkit of decision makers to anticipate the response of birds to fire management.

Forest degradation is known to be widespread in the tropics, but is currently very poorly mapped, in part because there is little quantitative data on which satellite sensor characteristics and analysis methods are best at detecting it. To improve this, we used data from the Tropical Forest Degradation Experiment (FODEX) plots in the southern Peruvian Amazon, where different numbers of trees had been removed from four 1 ha forest plots, carefully inventoried by hand and Terrestrial Laser Scanning before and after the logging to give a range of biomass change (ΔAGB) values. We conducted a comparative study of six multispectral optical satellite sensors (WorldView-3, SkySat, SPOT-7, PlanetScope, Sentinel-2 and Landsat 8) at 0.3 – 30 m spatial resolution, to find the best combination of sensor and remote sensing indicator for change detection. Spectral reflectance, the Normalized Difference Vegetation Index (NDVI) and texture parameters were extracted after radiometric calibration and image preprocessing. The strength of the relationships between the change in these values and field-measured ΔAGB (computed in % ha−1) was analysed. The results demonstrate that: (a) texture measures correlates more with ΔAGB than simple spectral parameters; (b) the strongest correlations are achieved for those sensors with spatial resolutions in the intermediate range (1.5 - 10 m), with finer or coarser resolutions producing worse results, and (c) when texture is computed using a moving square window ranging between 9 - 14 m in length. Maps predicting ΔAGB showed very promising results using a NIR-derived texture parameter for 3 m resolution PlanetScope (R2 = 0.97 and RMSE = 1.80 % ha−1), followed by 1.5 m SPOT-7 (R2 = 0.74 and RMSE = 5.25 % ha−1) and 10 m Sentinel-2 (R2 = 0.71 and RMSE = 5.55 % ha−1). Texture models derived from 0.3 m WorldView-3 improved with increasing window size, with highest R2 of 0.62 and RMSE = 6.35 % ha−1 for a window of 14 m in length. The degradation in our field plots is invisible to the 30 m resolution Landsat data. Our findings imply that, at least for lowland Peru, low-medium intensity disturbance can be detected best in optical wavelengths using a texture measure derived from 3 m PlanetScope data. That such data are being collected daily, and currently released free as monthly mosaics over tropical forests as part of the Norway’s International Climate and Forest Initiative (NICFI), is excellent news for monitoring such degradation.

A simple approach to enable water-management agencies employing free data to achieve the goal of using a single set of predictive equations for water-quality retrievals with satisfactory accuracy is proposed. Multiple regression-derived equations based on surface reflectance, band ratios, and environmental factors as predictor variables for concentrations of Total Suspended Solids (TSS), Total Nitrogen (TN), and Total Phosphorus (TP) were derived using a hybrid forward-selection method that considers Variance Inflation Factor (VIF) in the forward-selection process. Landsat TM, ETM+, and OLI/TIRS images were jointly utilized with environmental factors, such as wind speed and water surface temperature, to derive the single set of equations. The coefficients of determination of the best-fitting resultant equations varied from 0.62 to 0.79. Among all chosen predictor variables, ratio of reflectance of visible red (Band 3 for Landsat TM and ETM+, or Band 4 for Landsat OLI/TIRS) to visible blue (Band 1 for Landsat TM and ETM+, or Band 2 for Landsat OLI/TIRS) has a strong influence on the predictive power for TSS retrieval. Environmental factors including wind speed, remote sensing-derived water surface temperature, solar altitude, and time difference (in days) between the image acquisition and water sampling were found important in water-quality parameter estimation.

Densely urbanized areas, with a low percentage of green vegetation, are highly exposed to Heat Waves (HW) which nowadays are increasing in terms of frequency and intensity also in the middle-latitude regions, due to ongoing Climate Change (CC). Their negative effects may combine with those of the UHI (Urban Heat Island), a local phenomenon where air temperatures in the compact built up cores of towns increase more than those in the surrounding rural areas, with significant impact on the quality of urban environment, on citizens health and energy consumption and transport, as it has occurred in the summer of 2003 on France and Italian central-northern areas. In this context this work aims at designing and developing a methodology based on aero-spatial remote sensing (EO) at medium-high resolution and most recent GIS techniques, for the extensive characterization of the urban fabric response to these climatic impacts related to the temperature within the general framework of supporting local and national strategies and policies of adaptation to CC. Due to its extension and variety of built-up typologies, the municipality of Rome was selected as test area for the methodology development and validation. First of all, we started by operating through photointerpretation of cartography at detailed scale (CTR 1: 5000) on a reference area consisting of a transect of about 5x20 km, extending from the downtown to the suburbs and including all the built-up classes of interest. The reference built-up vulnerability classes found inside the transect were then exploited as training areas to classify the entire territory of Rome municipality. To this end, the satellite EO HR (High Resolution) multispectral data, provided by the Landsat sensors were used within a on purpose developed "supervised" classification procedure, based on data mining and “object-classification” techniques. The classification results were then exploited for implementing a calibration method, based on a typical UHI temperature distribution, derived from MODIS satellite sensor LST (Land Surface Temperature) data of the summer 2003, to obtain an analytical expression of the vulnerability model, previously introduced on a semi-empirical basis.