Arid and semi-arid regions have different spectral characteristics from other climatic regions. Therefore, appropriate remotely sensed indicators of land use and land cover types need to be defined for arid and semi-arid lands, as indices developed for other climatic regions may not give plausible results in arid and semi-arid regions. For instance, the normalized difference built-up index (NDBI) and normalized difference bareness index (NDBaI) are unable to distinguish between built-up areas and bare and dry soil that surrounds many cities in dry climates. This paper proposes the application of two newly developed indices, the dry built-up index (DBI) and dry bare-soil index (DBSI) to map built-up and bare areas in a dry climate from Landsat 8. The developed DBI and DBSI were applied to map urban areas and bare soil in the city of Erbil, Iraq. The results show an overall classification accuracy of 93% (κ = 0.86) and 92% (κ = 0.84) for DBI and DBSI, respectively. The results indicate the suitability of the proposed indices to discriminate between urban areas and bare soil in arid and semi-arid climates.

Land use and land cover (LULC) mapping can be of great help in changing land use decisions, but accurate mapping of LULC categories is challenging, especially in semi-arid areas with extensive farming systems and seasonal vegetation phenology. Machine learning algorithms are now widely used for LULC mapping because they provide analytical capabilities for LULC classification. However, the use of machine learning algorithms to improve classification performance is still being explored. The objective of this study is to investigate how to improve the performance of LULC models to reduce prediction errors. To address this question, the study applied a Random Forest (RF) based feature selection approach using Sentinel-1, -2, and Shuttle Radar Topographic Mission (SRTM) data. Results from RF show that the Sentinel-2 data only achieved an out-of-bag overall accuracy of 84.2%, while the Sentinel-1 and SRTM data achieved 83% and 76.44%, respectively. Classification accuracy improved to 89.1% when Sentinel-2, Sentinel-1 backscatter, and SRTM data were combined. This represents a 4.9% improvement in overall accuracy compared to Sentinel-2 alone and a 6.1% and 12.66% improvement compared to Sentinel-1 and SRTM data, respectively. Further independent validation, based on equally sized stratified random samples, consistently found a 5.3% difference between the Sentinel-2 and the combined datasets. This study demonstrates the importance of the synergy between optical, radar, and elevation data in improving the accuracy of LULC maps. In principle, the LULC maps produced in this study could help decision-makers in a wide range of spatial planning applications.

Although the way in which vegetation phenology mediates the feedback of vegetation to climate systems is now well understood, the magnitude of these changes is still unknown. A thorough understanding of how the recent shift in phenology may impact on, for example, land surface temperature (LST) is important. To address this knowledge gap, it is important to quantify these impacts and identify patterns from the global to the regional scale. This study examines the trend and linear regression modeling of the leaf area index (LAI) and LST derived from the moderate resolution imaging spectroradiometer (MODIS) data, specifically to assess their spatial distribution and changing trends at the continental and regional scales. The change detection analysis of interannual variability in the global LAI and LST between two periods (2003–2010 and 2011–2018) demonstrates more positive LAI trends than negative, while for LST most changes were not significant. The relationships between LAI and LST were assessed across the continents to ascertain the response of vegetation to changes in LST. The regression between LAI and LST was negative in Australia (R² = 0.487 ***), positive but minimal in Africa (R² = 0.001), positive in North America (R² = 0.641 ***), negative in Central America (R² = 0.119), positive in South America (R² = 0.253 *) and positive in Europe (R² = 0.740 ***). Medium temperatures enhance photosynthesis and lengthen the growing season in Europe. We also found a significant greening trend in China (trendp = 0.16 ***) and India (trendp = 0.13 ***). The relationships between LAI and LST in these most prominent greening countries of the world are R² = 0.06 and R² = 0.25 for China and India, respectively. Our deductions here are twofold—(1) In China, an insignificant association appeared between greening trend and temperature. (2) In India, the significant greening trend may be a factor in lowering temperatures. Therefore, temperature may stabilize if the greening trend continues. We attribute the trends in both countries to the different land use management and climate mitigation policies adopted by these countries.

Phenology is an important component in the climate system, and play a key role in controlling many feedbacks of vegetation to the climate systems. Differences in phenology of plant functional types (PFTs) due to the variation in seasonal cycles (e.g. changes in weather variability), the impact from land-use activities (e.g. fire) and their mechanisms for adaptation (e.g. climate change, regrowth/post-fire regeneration) have major implications in conservation planning and monitoring actions. Detecting phenological variations in PFTs is therefore of great importance to quantitative remote sensing applications, especially in a biome as diverse and complex as savannah. In this study, we implement Savitzky–Golay filtering and the Breaks For Additive Seasonal and Trend (BFAST) algorithms to detect changes in PFTs based on their phenological events using Moderate Resolution Imaging Spectroradiometer (MODIS), normalized difference vegetation index (NDVI) data from 2001 to 2018. In this region, PFTs present distinct seasonal, annual and interannual variability. Woody phenology presents early green-up dates with a prolonged senescence period and invariably observed the longest growing season length (GSL). The relationship between the start of season (SOS) or end of season (EOS) was assessed for each PFT to find out the extent to which they determine GSL for different PFTs. GSL is mostly determined by the SOS in woody savannah (coefficient of determination (R²) =0.41, p <0.01), open shrubs and (R² =0.79, p <0.001), grassland (R² =0.35, p <0.01) while EOS determined the GSL for both dryland crops (R² =0.75, p <0.001) and paddy rice (R² =0.69, p <0.001). We compared the interannual variability of woody savannah and other PFTs by measuring the differences of their phenological indicators using Welch’s t-test. All PFTs show statistically significant difference with the GSL of woody savannah except open shrubs (p = 0.23). The abrupt vegetation changes estimated with BFAST varied by PFT. Some PFTs are more resilient to harsh environmental conditions than others. While woody species present a few abrupt changes, grass phenology is more vulnerable to disturbance, seasonally as well as in the trend components (large number of abrupt changes) with browning trend following an abrupt negative change. In all PFTs, breakpoint (disturbance) assessed using BFAST negatively correlate with precipitation data which means the magnitude of disturbance decreases with increasing precipitation. Woody species had an r (correlation coefficient) value of -0.5 while grassland had r = -0.7 which is a further indication that grass phenology responds more strongly to annual precipitation than the woody species. These results show that MODIS NDVI time-series data can be used to distinguish the phenological events of different PFTs in West African savannah dominated landscape.

Assessments of tree/grass fractional cover in savannahs using remote sensing are challenging due to the heterogeneous mixture of the two plant functional types. Time-series decomposition models can be used to characterize vegetation phenology from satellite data, but have rarely been used for attributing phenological signal components to different plant functional types. Here, tree/grass dynamics are assessed in savannah ecosystems using time-series decomposition of 14 years of Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index data acquired from 2002 to 2015. The decomposition method uses harmonic analysis and tests the individual harmonic terms for statistical significance. Field data of fractional cover of trees and grasses were collected for 28 plots in Kruger National Park, South Africa. Matching MODIS pixels were analysed for their tree/grass phenological signals. Tree/grass annual and interannual variability were then assessed based on the harmonic models. In most harmonic cycles, grass-dominated sites had higher amplitudes than tree-dominated sites, while the tree green-up started earlier than grasses, before the start of the wet season. While changes in tree phenology are gradual, grasses present higher variability over time. Tree cover showed a significant correlation with the amplitude (r (correlation coefficient) = −0.59, p = 0.001) and phase of the first harmonic term (r = −0.73, p = 0.0001) and the number of cycles of the second harmonic term (r = 0. 56, p = 0.002). Grass cover was also significantly correlated with the amplitude (r = 0. 51, p = 0.005) and phase of the first harmonic term (r = 0.55, p = 0.002) and the number of cycles of the second harmonic term (r = −0.52, p = 0.005). The positive correlation of grass cover with phase and negative correlation with number of cycles is indicating a late greening period and higher variability, respectively. Tree cover estimated from the phase of the strongest harmonic term showed a positive correlation with field-measured tree cover (R² (coefficient of determination) = 0.55, p < 0.01, slope = 0.93, root mean square error = 13.26%). The estimated tree cover also had a strong correlation with the woody cover map (r = 0.78, p < 0.01) produced by Bucini. The results show that MODIS time-series data can be used to estimate the fractional tree cover in heterogeneous savannahs from the phase of the plant functional type’s phenological behaviour. This study shows that harmonic analysis is able to discriminate between fractional cover by trees and grasses in savannahs. The quantitative analysis of tree/grass phenology from satellite time-series data enables a better understanding of the dynamics of the tree/grass competition and coexistence.

Monitoring the ecological and socioeconomic impacts of wildfires using traditional approaches requires significant financial resources, time, and sampling expertise. However, not only are resources scarce, but the spatial and temporal extent of forest fires can also make it impractical to assess large areas over time. Thus, fire monitoring initiatives are often not realized. This has inevitably made the remote sensing approach an interesting option for fire protection managers and decision-makers due to its ability to measure large areas and its temporal capabilities. In this study, burn spectral indices derived from Landsat 8 (difference normalized vegetation index (dNDVI) and difference normalized burn ratio (dNBR)) were used to assess the ecological and socioeconomic impacts of forest fires based on an existing land use/land cover dataset. The relationships between estimated fire severity/area and environmental and anthropogenic factors were also evaluated. The results show that more than 700 hectares of forest and other land use categories were burned. Fires adversely affect high forests, thickets, degraded forests, and most cultivated and rural areas. The study also revealed a moderate positive relationship between burn severity and pre-fire vegetation (R ² = 0.48 and R ² = 0.49 for the dNDVI and dNBR, respectively). This result suggested that the fuel amount is the main driver of burn severity during the fire season in this particular ecosystem. Topography has been shown to affect fire behavior in the study area, where fires occur primarily at elevations averaging 400-800 meters above mean sea level. In contrast, there is a weak positive relationship between population density and burnt area. This phenomenon is commonly observed in specific regions, where the incidence of fire is directly proportional to the density of the population. However, the severity decreases when burning exceeds a threshold. This study has shown that Landsat 8 data-derived burn spectral indices (dNDVI and dNBR) have high potential for the spatial analysis of wildfires.