Reliable crop type maps from satellite data are an essential prerequisite for quantifying crop growth, health, and yields. However, such maps do not exist for most parts of Africa, where smallholder farming is the dominant system. Prevalent cloud cover, small farm sizes, and mixed cropping systems pose substantial challenges when creating crop type maps for sub-Saharan Africa. In this study, we provide a mapping scheme based on freely available Sentinel-2A/B (S2) time series and very high-resolution SkySat data to map the main crops—maize and potato—and intercropping systems including these two crops on the Jos Plateau, Nigeria. We analyzed the spectral-temporal behavior of mixed crop classes to improve our understanding of inter-class spectral mixing. Building on the Framework for Operational Radiometric Correction for Environmental monitoring (FORCE), we preprocessed S2 time series and derived spectral-temporal metrics from S2 spectral bands for the main temporal cropping windows. These STMs were used as input features in a hierarchical random forest classification. Our results provide the first wall-to-wall crop type map for this key agricultural region of Nigeria. Our cropland identification had an overall accuracy of 84%, while the crop type map achieved an average accuracy of 72% for the five relevant crop classes. Our crop type map shows distinctive regional variations in the distribution of crop types. Maize is the dominant crop, followed by mixed cropping systems, including maize–cereals and potato–maize cropping; potato was found to be the least prevalent class. Plot analyses based on a sample of 1166 fields revealed largely homogeneous mapping patterns, demonstrating the effectiveness of our classification system also for intercropped classes, which are temporally and spatially highly heterogeneous. Moreover, we found that small field sizes were dominant in all crop types, regardless of whether or not intercropping was used. Maize–legume and maize exhibited the largest plots, with an area of up to 3 ha and slightly more than 10 ha, respectively; potato was mainly cultivated on fields smaller than 0.5 ha and only a few plots were larger than 1 ha. Besides providing the first spatially explicit map of cropping practices in the core production area of the Jos Plateau, Nigeria, the study also offers guidance for the creation of crop type maps for smallholder-dominated systems with intercropping. Critical temporal windows for crop type differentiation will enable the creation of mapping approaches in support of future smart agricultural practices for aspects such as food security, early warning systems, policies, and extension services.

Air temperatures in the Arctic have increased substantially over the last decades, which has extensively altered the properties of the land surface. Capturing the state and dynamics of Land Surface Temperatures (LSTs) at high spatial detail is of high interest as LST is dependent on a variety of surficial properties and characterizes the land–atmosphere exchange of energy. Accordingly, this study analyses the influence of different physical surface properties on the long-term mean of the summer LST in the Arctic Mackenzie Delta Region (MDR) using Landsat 30 m-resolution imagery between 1985 and 2018 by taking advantage of the cloud computing capabilities of the Google Earth Engine. Multispectral indices, including the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI) and Tasseled Cap greenness (TCG), brightness (TCB), and wetness (TCW) as well as topographic features derived from the TanDEM-X digital elevation model are used in correlation and multiple linear regression analyses to reveal their influence on the LST. Furthermore, surface alteration trends of the LST, NDVI, and NDWI are revealed using the Theil-Sen (T-S) regression method. The results indicate that the mean summer LST appears to be mostly influenced by the topographic exposition as well as the prevalent moisture regime where higher evapotranspiration rates increase the latent heat flux and cause a cooling of the surface, as the variance is best explained by the TCW and northness of the terrain. However, fairly diverse model outcomes for different regions of the MDR (R² from 0.31 to 0.74 and RMSE from 0.51 °C to 1.73 °C) highlight the heterogeneity of the landscape in terms of influential factors and suggests accounting for a broad spectrum of different factors when modeling mean LSTs. The T-S analysis revealed large-scale wetting and greening trends with a mean decadal increase of the NDVI/NDWI of approximately +0.03 between 1985 and 2018, which was mostly accompanied by a cooling of the land surface given the inverse relationship between mean LSTs and vegetation and moisture conditions. Disturbance through wildfires intensifies the surface alterations locally and lead to significantly cooler LSTs in the long-term compared to the undisturbed surroundings.