ALBERT

Description: Comprehensive research on glacier changes in the Tian Shan is available for the current decade; however, there is limited information about glacier investigations of previous decades and especially before the mid 1970s. The earliest stereo images from the Corona missions were acquired in the 1960s but existing studies dealing with these images focus on single glaciers or small areas only. We developed a workflow to generate digital terrain models (DTMs) and orthophotos from 1964 Corona KH-4 for an entire mountain range (Ak-Shirak) located in the Central Tian Shan. From these DTMs and orthoimages, we calculated geodetic mass balances and length changes in comparison to 1973 and 1980 Hexagon KH-9 data. We found mass budgets between −0.4 ± 0.1 m·w.e.a−1 (1964–1980) and −0.9 ± 0.4 m·w.e.a−1 (1973–1980) for the whole region and individual glaciers. The length changes, on the other hand, vary heterogeneously between +624 ± 18 m (+39.0 ± 1.1 m·a−1) and −923 ± 18 m (−57.7 ± 1.1 m·a−1) for 1964–1980. An automation of the processing line can successively lead to region-wide Corona data processing allowing the analysis and interpretation of glacier changes on a larger scale and supporting a refinement of glacier modelling.

Description: Geographic Object-Based Image Analysis (GEOBIA) mostly uses proprietary software,but the interest in Free and Open-Source Software (FOSS) for GEOBIA is growing. This interest stems not only from cost savings, but also from benefits concerning reproducibility and collaboration. Technical challenges hamper practical reproducibility, especially when multiple software packages are required to conduct an analysis. In this study, we use containerization to package a GEOBIA workflow in a well-defined FOSS environment. We explore the approach using two software stacks to perform an exemplary analysis detecting destruction of buildings in bi-temporal images of a conflict area. The analysis combines feature extraction techniques with segmentation and object-based analysis to detect changes using automatically-defined local reference values and to distinguish disappeared buildings from non-target structures. The resulting workflow is published as FOSS comprising both the model and data in a ready to use Docker image and a user interface for interaction with the containerized workflow. The presented solution advances GEOBIA in the following aspects: higher transparency of methodology; easier reuse and adaption of workflows; better transferability between operating systems; complete description of the software environment; and easy application of workflows by image analysis experts and non-experts. As a result, it promotes not only the reproducibility of GEOBIA, but also its practical adoption.

Description: Digital surface models (DSMs) derived from spaceborne and airborne sensors enable the monitoring of the vertical structures for forests in large areas. Nevertheless, due to the lack of an objective performance assessment for this task, it is difficult to select the most appropriate data source for DSM generation. In order to fill this gap, this paper performs change detection analysis including forest decrease and tree growth. The accuracy of the DSMs is evaluated by comparison with measured tree heights from inventory plots (field data). In addition, the DSMs are compared with LiDAR data to perform a pixel-wise quality assessment. DSMs from four different satellite stereo sensors (ALOS/PRISM, Cartosat-1, RapidEye and WorldView-2), one satellite InSAR sensor (TanDEM-X), two aerial stereo camera systems (HRSC and UltraCam) and two airborne laser scanning datasets with different point densities are adopted for the comparison. The case study is a complex central European temperate forest close to Traunstein in Bavaria, Germany. As a major experimental result, the quality of the DSM is found to be robust to variations in image resolution, especially when the forest density is high. The forest decrease results confirm that besides aerial photogrammetry data, very high resolution satellite data, such as WorldView-2, can deliver results with comparable quality as the ones derived from LiDAR, followed by TanDEM-X and Cartosat DSMs. The quality of the DSMs derived from ALOS and Rapid-Eye data is lower, but the main changes are still correctly highlighted. Moreover, the vertical tree growth and their relationship with tree height are analyzed. The major tree height in the study site is between 15 and 30 m and the periodic annual increments (PAIs) are in the range of 0.30–0.50 m.

Description: This paper presents an automated and effective framework for classifying airborne laser scanning (ALS) point clouds. The framework is composed of four stages: (i) step-wise point cloud segmentation, (ii) feature extraction, (iii) Random Forests (RF) based feature selection and classification, and (iv) post-processing. First, a step-wise point cloud segmentation method is proposed to extract three kinds of segments, including planar, smooth and rough surfaces. Second, a segment, rather than an individual point, is taken as the basic processing unit to extract features. Third, RF is employed to select features and classify these segments. Finally, semantic rules are employed to optimize the classification result. Three datasets provided by Open Topography are utilized to test the proposed method. Experiments show that our method achieves a superior classification result with an overall classification accuracy larger than 91.17%, and kappa coefficient larger than 83.79%.

Description: Nitrogen is an essential nutrient in many terrestrial ecosystems because it affects vegetation’s primary production. Due to the variety of nitrogen-containing substances and the differences in their composition across species, statistical approaches are now dominant in remote sensing retrieval of leaf nitrogen content. Many studies remove spectral regions characterized by strong water absorptions before retrieving nitrogen content, because water is believed to mask the absorption features of nitrogen. The objectives of this study are to discuss the necessity of this practice and to explore how water absorption affects leaf nitrogen estimation. Spectral measurements and chemical analyses for Maize, Sawtooth Oak, and Sweetgum leaves were carried out in 2014. The leaf optical properties model PROSPECT5 was used to eliminate the influences of water on the measured reflectance spectra. The inversion accuracy of PROPECT5 for chlorophyll, carotenoid, water, and dry matter of Maize was also discussed. Measured, simulated, and water-removed spectra were used to: (1) find the optimal nitrogen-related spectral index; and (2) regress with the area-based leaf nitrogen concentration (LNC) using the partial least square regression technique (PLSR). Two types of spectral indices were selected in this study: Normalized Difference Spectral Index (NDSI) and Ratio Spectral Index (RSI). Additionally, first-order derivative forms of measured, simulated, and water-removed spectra were devised to search for the optimal spectral indices. Finally, species-specific optimal indices and cross-species optimal indices, as well as their root mean square errors (RMSE) and coefficients of determination (R2), were obtained. The Ending Top Percentile (ETP), an indicator of the performance of cross-species optimal indices, was also calculated. PLSR was combined with leave-one-out cross validation (LOOCV) for each species. The predicted root mean square errors (RMSEP) and predicted R2 were finally calculated. The results showed that chlorophyll, carotenoid, and water contents could be estimated with R2 of 0.75, 0.59, and 0.69, respectively, which were acceptable for fresh leaves. The dry matter was retrieved with a relatively lower accuracy because of the fixed absorption coefficients adopted by PROSPECT5. The performances of species-specific optimal indices using water-free spectra were comparable to or worse than the corresponding indices derived with measured or simulated spectra. Compared with measured spectra, ETP did not change much after the effects of water were removed, and the R2 between cross-species optimal spectral indices and area-based LNC for Sawtooth Oak and Sweetgum decreased while it remained almost the same for Maize, suggesting that the water-removed cross-species optimal indices were inferior to the corresponding optimal indices found without water removal. ETP was larger than 30% for all spectra, demonstrating the non-existence of common optimal NDSI or RSI for the three species. After water removal, the accuracy of PLSR for Sawtooth Oak and Sweetgum decreased and increased negligibly for Maize. The results suggest that water absorption has limited effects on reducing the accuracy of leaf nitrogen estimation. On the contrary, the accuracy may decrease due to the loss of spectral information caused by the removal of water-sensitive spectral regions.

Description: Accurate and high resolution bathymetric data is a necessity for a wide range of coastal oceanographic research topics. Active sensing methods, such as ship-based soundings and Light Detection and Ranging (LiDAR), are expensive and time consuming solutions. Therefore, the significance of Satellite-Derived Bathymetry (SDB) has increased in the last ten years due to the availability of multi-constellation, multi-temporal, and multi-resolution remote sensing data as Open Data. Effective SDB algorithms have been proposed by many authors, but there is no ready-to-use software module available in the Geographical Information System (GIS) environment as yet. Hence, this study implements a Geographically Weighted Regression (GWR) based SDB workflow as a Geographic Resources Analysis Support System (GRASS) GIS module (i.image.bathymetry). Several case studies were carried out to examine the performance of the module in multi-constellation and multi-resolution satellite imageries for different study areas. The results indicate a strong correlation between SDB and reference depth. For instance, case study 1 (Puerto Rico, Northeastern Caribbean Sea) has shown an coefficient of determination (R2) of 0.98 and an Root Mean Square Error (RMSE) of 0.61 m, case study 2 (Iwate, Japan) has shown an R2 of 0.94 and an RMSE of 1.50 m, and case study 3 (Miyagi, Japan) has shown an R2 of 0.93 and an RMSE of 1.65 m. The reference depths were acquired by using LiDAR for case study 1 and an echo-sounder for case studies 2 and 3. Further, the estimated SDB has been used as one of the inputs for the Australian National University and Geoscience Australia (ANUGA) tsunami simulation model. The tsunami simulation results also show close agreement with post-tsunami survey data. The i.mage.bathymetry module developed as a part of this study is made available as an extension for the Open Source GRASS GIS to facilitate wide use and future improvements.

Description: The Soil Moisture Active Passive (SMAP) satellite makes coincident global measurements of soil moisture using an L-band radar instrument and an L-band radiometer. It is crucial to evaluate the errors in the newest L-band SMAP satellite-derived soil moisture products, before they are routinely used in scientific research and applications. This study represents the first evaluation of the SMAP radiometer soil moisture product over China. In this paper, a preliminary evaluation was performed using sparse in situ measurements from 655 China Meteorological Administration (CMA) monitoring stations between 1 April 2015 and 31 August 2016. The SMAP radiometer-derived soil moisture product was evaluated against two schemes of original soil moisture and the soil moisture anomaly in different geographical zones and land cover types. Four performance metrics, i.e., bias, root mean square error (RMSE), unbiased root mean square error (ubRMSE), and the correlation coefficient (R), were used in the accuracy evaluation. The results indicated that the SMAP radiometer-derived soil moisture product agreed relatively well with the in situ measurements, with ubRMSE values of 0.058 cm3·cm−3 and 0.039 cm3·cm−3 based on original data and anomaly data, respectively. The values of the SMAP radiometer-based soil moisture product were overestimated in wet areas, especially in the Southwest China, South China, Southeast China, East China, and Central China zones. The accuracies over croplands and in Northeast China were the worst. Soil moisture, surface roughness, and vegetation are crucial factors contributing to the error in the soil moisture product. Moreover, radio frequency interference contributes to the overestimation over the northern portion of the East China zone. This study provides guidelines for the application of the SMAP-derived soil moisture product in China and acts as a reference for improving the retrieval algorithm.

Description: Vegetation is often represented by the leaf area index (LAI) in many ecological, hydrological and meteorological land surface models. However, the spatio-temporal dynamics of the vegetation are important to represent in these models. While the widely applied methods, such as the Canopy Structure Dynamic Model (CSDM) and the Double Logistic Model (DLM) are solely based on cumulative daily mean temperature data as input, a new spatio-temporal LAI prediction model referred to as the Temperature Precipitation Vegetation Model (TPVM) is developed that also considers cumulative precipitation data as input into the modelling process. TPVM as well as CDSM and DLM model performances are compared and evaluated against filtered LAI data from the Moderate Resolution Imaging Spectroradiometer (MODIS). The calibration/validation results of a cross-validation performed in the meso-scale Attert catchment in Luxembourg indicated that the DLM and TPVM generally provided more realistic and accurate LAI data. The TPVM performed superiorly for the agricultural land cover types compared to the other two models, which only used the temperature data. The Pearson's correlation coefficient (CC) between TPVM and the field measurement is 0.78, compared to 0.73 and 0.69 for the DLM and CSDM, respectively. The phenological metrics were derived from the TPVM model to investigate the interaction between the climate variables and the LAI variations. These interactions illustrated the dominant control of temperature on the LAI dynamics for deciduous forest cover, and a combined influence of temperature with precipitation for the agricultural land use areas.

Description: Water vapour (H2O) is the dominant species in volcanic gas plumes. Therefore, measurements of H2O fluxes could provide valuable constraints on subsurface degassing and magmatic processes. However, due to the large and variable concentration of this species in the background atmosphere, little attention has been devoted to monitoring the emission rates of this species from volcanoes. Instead, the focus has been placed on remote measurements of SO2, which is present in far lower abundances in plumes, and therefore provides poorer single flux proxies for overall degassing conditions. Here, we present a new technique for the measurement of H2O emissions at degassing volcanoes at high temporal resolution (≈1 Hz), via remote sensing with low cost digital cameras. This approach is analogous to the use of dual band ultraviolet (UV) cameras for measurements of volcanic SO2 release, but is focused on near infrared absorption by H2O. We report on the field deployment of these devices on La Fossa crater, Vulcano Island, and the North East Crater of Mt. Etna, during which in-plume calibration was performed using a humidity sensor, resulting in estimated mean H2O fluxes of ≈15 kg·s−1 and ≈34 kg·s−1, respectively, in accordance with previously reported literature values. By combining the Etna data with parallel UV camera and Multi-GAS observations, we also derived, for the first time, a combined record of 1 Hz gas fluxes for the three most abundant volcanic gas species: H2O, CO2, and SO2. Spectral analysis of the Etna data revealed oscillations in the passive emissions of all three species, with periods spanning ≈40–175 s, and a strong degree of correlation between the periodicity manifested in the SO2 and H2O data, potentially related to the similar exsolution depths of these two gases. In contrast, there was a poorer linkage between oscillations in these species and those of CO2, possibly due to the deeper exsolution of carbon dioxide, giving rise to distinct periodic degassing behaviour.

Description: Identifying individual trees and delineating their canopy structures from the forest point clouddataacquiredbyanairborneLiDAR(LightDetectionAndRanging)hassignificantimplications in forestry inventory. Once accurately identified, tree structural attributes such as tree height, crown diameter, canopy based height and diameter at breast height can be derived. This paper focuses on a novel computationally efficient method to adaptively calibrate the kernel bandwidth of a computational scheme based on mean shift—a non-parametric probability density-based clustering technique—to segment the 3D (three-dimensional) forest point clouds and identify individual tree crowns. The basic concept of this method is to partition the 3D space over each test plot into small vertical units (irregular columns containing 3D spatial features from one or more trees) first, by using a fixed bandwidth mean shift procedure and a small square grouping technique, and then rough estimation of crown sizes for distinct trees within a unit, based on an original 2D (two-dimensional) incremental grid projection technique, is applied to provide a basis for dynamical calibration of the kernel bandwidth for an adaptive mean shift procedure performed in each partition. The adaptive mean shift-based scheme, which incorporates our proposed bandwidth calibration method, is validated on 10 test plots of a dense, multi-layered evergreen broad-leaved forest located in South China. Experimental results reveal that this approach can work effectively and when compared to the conventional point-based approaches (e.g., region growing, k-means clustering, fixed bandwidth or multi-scale mean shift), its accuracies are relatively high: it detects 86 percent of the trees (“recall”) and 92 percent of the identified trees are correct (“precision”), showing good potential for use in the area of forest inventory.

Description: The Advanced Along-Track Scanning Radiometer (AATSR) land surface temperature (LST) product has a long-term time series of data from 20 May 2002 to 8 April 2012 and is a crucial dataset for global change studies. Accuracy and uncertainty assessment of satellite derived LST is important for its use in studying land–surface–atmosphere interactions. However, the validation of AATSR-derived LST products is scarce in China, especially in arid and semi-arid areas. In this study, we evaluated the accuracy of the AATSR LST product using ground-based measurements from 2007 to 2011 in the Heihe River Basin (HRB), China. The AATSR-derived LST results over Yingke site are closer to ground measurements than those over A’rou site for both daytime and nighttime temperatures. For nighttime, the averaged bias, STD, RMSE and R2 over both sites are 0.67 K, 3.03 K, 3.13 K and 0.93 K, respectively. Based on the accuracy assessment, we analyzed the AATSR-derived annual LST variations both in the HRB region and the two validation sites for the period of 2003 to 2011. The results at the A’rou site show an obvious increasing trend for daytime from 2003 to 2011. For the whole HRB region, the warming trend is clearly shown in the downstream of HRB.

Description: Recently, we reported on the development of low-cost ultraviolet (UV) cameras, based on the modification of sensors designed for the smartphone market. These units are built around modified Raspberry Pi cameras (PiCams; ≈USD 25), and usable system sensitivity was demonstrated in the UVA and UVB spectral regions, of relevance to a number of application areas. Here, we report on the first deployment of PiCam devices in one such field: UV remote sensing of sulphur dioxide emissions from volcanoes; such data provide important insights into magmatic processes and are applied in hazard assessments. In particular, we report on field trials on Mt. Etna, where the utility of these devices in quantifying volcanic sulphur dioxide (SO2) emissions was validated. We furthermore performed side-by-side trials of these units against scientific grade cameras, which are currently used in this application, finding that the two systems gave virtually identical flux time series outputs, and that signal-to-noise characteristics of the PiCam units appeared to be more than adequate for volcanological applications. Given the low cost of these sensors, allowing two-filter SO2 camera systems to be assembled for ≈USD 500, they could be suitable for widespread dissemination in volcanic SO2 monitoring internationally.

Description: The assessment of the soil redistribution and real long-term soil degradation due to erosion on agriculture land is still insufficient in spite of being essential for soil conservation policy. Imaging spectroscopy has been recognized as a suitable tool for soil erosion assessment in recent years. In our study, we bring an approach for assessment of soil degradation by erosion by means of determining soil erosion classes representing soils differently influenced by erosion impact. The adopted methods include extensive field sampling, laboratory analysis, predictive modelling of selected soil surface properties using aerial hyperspectral data and the digital elevation model and fuzzy classification. Different multivariate regression techniques (Partial Least Square, Support Vector Machine, Random forest and Artificial neural network) were applied in the predictive modelling of soil properties. The properties with satisfying performance (R2 > 0.5) were used as input data in erosion classes determination by fuzzy C-means classification method. The study was performed at four study sites about 1 km2 large representing the most extensive soil units of the agricultural land in the Czech Republic (Chernozems and Luvisols on loess and Cambisols and Stagnosols on crystalline rocks). The influence of site-specific conditions on prediction of soil properties and classification of erosion classes was assessed. The prediction accuracy (R2) of the best performing models predicting the soil properties varies in range 0.8–0.91 for soil organic carbon content, 0.21–0.67 for sand content, 0.4–0.92 for silt content, 0.38–0.89 for clay content, 0.73–089 for Feox, 0.59–0.78 for Fed and 0.82 for CaCO3. The performance and suitability of different properties for erosion classes’ classification are highly variable at the study sites. Soil organic carbon was the most frequently used as the erosion classes’ predictor, while the textural classes showed lower applicability. The presented approach was successfully applied in Chernozem and Luvisol loess regions where the erosion classes were assessed with a good overall accuracy (82% and 67%, respectively). The model performance in two Cambisol/Stagnosol regions was rather poor (51%–52%). The results showed that the presented method can be directly and with a good performance applied in pedologically and geologically homogeneous areas. The sites with heterogeneous structure of the soil cover and parent material will require more precise local-fitted models and use of further auxiliary information such as terrain or geological data. The future application of presented approach at a regional scale promises to produce valuable data on actual soil degradation by erosion usable for soil conservation policy purposes.

Description: The methods used for image contrast enhancement in the wavelet domain have been previously documented. The essence of these methods lies in the manipulation of the image during the reconstruction process, by changing the relationship between the components that require transformation. This paper proposes a new variant based on using undecimated wavelet transform and adapting the Gaussian function for scaling the coefficients of detail wavelet components, so that the role of low coefficients in the reconstructed image is greater. The enhanced image is then created by combining the new components. Applying the Haar wavelet minimises the effects of the relationship disturbance between components, and creates a small buffer around the edge. The proposed method was tested using six images at different scales, collected with handheld photo cameras, and aerial and satellite optical sensors. The results of the tests indicate that the method can achieve comparable, or even better enhancement effects for weak edges, than the well-known unsharp masking and Retinex methods. The proposed method can be applied in order to improve the visual interpretation of remote sensing images taken by various sensors at different scales.

Description: Remote Sensing, Vol. 9, Pages 820: Soil Moisture Data Assimilation in a Hydrological Model: A Case Study in Belgium Using Large-Scale Satellite Data Remote Sensing doi: 10.3390/rs9080820 Authors: Pierre Baguis Emmanuel Roulin In the present study, we focus on the assimilation of satellite observations for Surface Soil Moisture (SSM) in a hydrological model. The satellite data are produced in the framework of the EUMETSAT project H-SAF and are based on measurements with the Advanced radar Scatterometer (ASCAT), embarked on the Meteorological Operational satellites (MetOp). The product generated with these measurements has a horizontal resolution of 25 km and represents the upper few centimeters of soil. Our approach is based on the Ensemble Kalman Filter technique (EnKF), where observation and model uncertainties are taken into account, implemented in a conceptual hydrological model. The analysis is carried out in the Demer catchment of the Scheldt River Basin in Belgium, for the period from June 2013–May 2016. In this context, two methodological advances are being proposed. First, the generation of stochastic terms, necessary for the EnKF, of bounded variables like SSM is addressed with the aid of specially-designed probability distributions, so that the bounds are never exceeded. Second, bias due to the assimilation procedure itself is removed using a post-processing technique. Subsequently, the impact of SSM assimilation on the simulated streamflow is estimated using a series of statistical measures based on the ensemble average. The differences from the control simulation are then assessed using a two-dimensional bootstrap sampling on the ensemble generated by the assimilation procedure. Our analysis shows that data assimilation combined with bias correction can improve the streamflow estimations or, at a minimum, produce results statistically indistinguishable from the control run of the hydrological model.

Description: Remote Sensing, Vol. 9, Pages 823: Low-Altitude Aerial Methane Concentration Mapping Remote Sensing doi: 10.3390/rs9080823 Authors: Bara Emran Dwayne Tannant Homayoun Najjaran Detection of leaks of fugitive greenhouse gases (GHGs) from landfills and natural gas infrastructure is critical for not only their safe operation but also for protecting the environment. Current inspection practices involve moving a methane detector within the target area by a person or vehicle. This procedure is dangerous, time consuming, labor intensive and above all unavailable when access to the desired area is limited. Remote sensing by an unmanned aerial vehicle (UAV) equipped with a methane detector is a cost-effective and fast method for methane detection and monitoring, especially for vast and remote areas. This paper describes the integration of an off-the-shelf laser-based methane detector into a multi-rotor UAV and demonstrates its efficacy in generating an aerial methane concentration map of a landfill. The UAV flies a preset flight path measuring methane concentrations in a vertical air column between the UAV and the ground surface. Measurements were taken at 10 Hz giving a typical distance between measurements of 0.2 m when flying at 2 m/s. The UAV was set to fly at 25 to 30 m above the ground. We conclude that besides its utility in landfill monitoring, the proposed method is ready for other environmental applications as well as the inspection of natural gas infrastructure that can release methane with much higher concentrations.

Description: Remote Sensing, Vol. 9, Pages 815: Validation of Automatically Generated Global and Regional Cropland Data Sets: The Case of Tanzania Remote Sensing doi: 10.3390/rs9080815 Authors: Juan Laso Bayas Linda See Christoph Perger Christina Justice Catherine Nakalembe Jan Dempewolf Steffen Fritz There is a need to validate existing global cropland maps since they are used for different purposes including agricultural monitoring and assessment. In this paper we validate three recent global products (ESA-CCI, GlobeLand30, FROM-GC) and one regional product (Tanzania Land Cover 2010 Scheme II) using a validation data set that was collected by students through the Geo-Wiki tool. The ultimate aim was to understand the usefulness of these products for agricultural monitoring. Data were collected wall-to-wall for Kilosa district and for a sample across Tanzania. The results show that the amount of and spatial extent of cropland in the different products differs considerably from 8% to 42% for Tanzania, with similar values for Kilosa district. The agreement of the validation data with the four different products varied between 36% and 54% and highlighted that cropland is overestimated by the ESA-CCI and underestimated by FROM-GC. The validation data were also analyzed for consistency between the student interpreters and also compared with a sample interpreted by five experts for quality assurance. Regarding consistency between the students, there was more than 80% agreement if one difference in cropland category was considered (e.g., between low and medium cropland) while most of the confusion with the experts was also within one category difference. In addition to the validation of current cropland products, the data set collected by the students also has potential value as a training set for improving future cropland products.

Description: Remote Sensing, Vol. 9, Pages 816: Mapping Aboveground Carbon in Oil Palm Plantations Using LiDAR: A Comparison of Tree-Centric versus Area-Based Approaches Remote Sensing doi: 10.3390/rs9080816 Authors: Matheus Nunes Robert Ewers Edgar Turner David Coomes Southeast Asia is the epicentre of world palm oil production. Plantations in Malaysia have increased 150% in area within the last decade, mostly at the expense of tropical forests. Maps of the aboveground carbon density (ACD) of vegetation generated by remote sensing technologies, such as airborne LiDAR, are vital for quantifying the effects of land use change for greenhouse gas emissions, and many papers have developed methods for mapping forests. However, nobody has yet mapped oil palm ACD from LiDAR. The development of carbon prediction models would open doors to remote monitoring of plantations as part of efforts to make the industry more environmentally sustainable. This paper compares the performance of tree-centric and area-based approaches to mapping ACD in oil palm plantations. We find that an area-based approach gave more accurate estimates of carbon density than tree-centric methods and that the most accurate estimation model includes LiDAR measurements of top-of-canopy height and canopy cover. We show that tree crown segmentation is sensitive to crown density, resulting in less accurate tree density and ACD predictions, but argue that tree-centric approach can nevertheless be useful for monitoring purposes, providing a method to detect, extract and count oil palm trees automatically from images.

Description: Remote Sensing, Vol. 9, Pages 818: Water Optics and Water Colour Remote Sensing Remote Sensing doi: 10.3390/rs9080818 Authors: Yunlin Zhang Claudia Giardino Linhai Li The editorial paper aims to highlight the main topics investigated in the Special Issue (SI) “Water Optics and Water Colour Remote Sensing”. The outcomes of the 21 papers published in the SI are presented, along with a bibliometric analysis in the same field, namely, water optics and water colour remote sensing. This editorial summarises how the research articles of the SI approach the study of bio-optical properties of aquatic systems, the development of remote sensing algorithms, and the application of time-series satellite data for assessing long-term and temporal-spatial dynamics in inland, coastal, and oceanic waters. The SI shows the progress with a focus on: (1) bio-optical properties (three papers); (2) atmospheric correction and data uncertainties (five papers); (3) remote sensing estimation of chlorophyll-a (Chl-a) (eight papers); (4) remote sensing estimation of suspended matter and chromophoric dissolved organic matter (CDOM) (four papers); and (5) water quality and water ecology remote sensing (four papers). Overall, the SI presents a variety of applications at the global scale (with case studies in Europe, Asia, South and North America, and the Antarctic), achieved with different remote sensing instruments, such as hyperspectral field and airborne sensors, ocean colour radiometry, geostationary platforms, and the multispectral Landsat and Sentinel-2 satellites. The bibliometric analysis, carried out to include research articles published from 1900 to 2016, indicates that “chlorophyll-a”, “ocean colour”, “phytoplankton”, “SeaWiFS” (Sea-Viewing Wide Field-of-View Sensor), and “chromophoric dissolved organic matter” were the five most frequently used keywords in the field. The SI contents, along with the bibliometric analysis, clearly suggest that remote sensing of Chl-a is one of the topmost investigated subjects in the field.

Description: Remote Sensing, Vol. 9, Pages 811: Radiometric Cross-Calibration of GF-4 PMS Sensor Based on Assimilation of Landsat-8 OLI Images Remote Sensing doi: 10.3390/rs9080811 Authors: Yepei Chen Kaimin Sun Deren Li Ting Bai Chengquan Huang Earth observation data obtained from remote sensors must undergo radiometric calibration before use in quantitative applications. However, the large view angles of the panchromatic multispectral sensor (PMS) aboard the GF-4 satellite pose challenges for cross-calibration due to the effects of atmospheric radiation transfer and the bidirectional reflectance distribution function (BRDF). To address this problem, this paper introduces a novel cross-calibration method based on data assimilation considering cross-calibration as an optimal approximation problem. The GF-4 PMS was cross-calibrated with the well-calibrated Landsat-8 Operational Land Imager (OLI) as the reference sensor. In order to correct unequal bidirectional reflection effects, an adjustment factor for the BRDF was established, making complex models unnecessary. The proposed method employed the Shuffled Complex Evolution-University of Arizona (SCE-UA) algorithm to find the optimal calibration coefficients and BRDF adjustment factor through an iterative process. The validation results revealed a surface reflectance error of <5% for the new cross-calibration coefficients. The accuracy of calibration coefficients were significantly improved when compared to the officially published coefficients as well as those derived using conventional methods. The uncertainty produced by the proposed method was less than 7%, meeting the demands for future quantitative applications and research. This method is also applicable to other sensors with large view angles.

Description: Remote Sensing, Vol. 9, Pages 817: Estimation of Satellite-Based SO42− and NH4+ Composition of Ambient Fine Particulate Matter over China Using Chemical Transport Model Remote Sensing doi: 10.3390/rs9080817 Authors: Yidan Si Shenshen Li Liangfu Chen Chao Yu Wende Zhu Epidemiologic and health impact studies have examined the chemical composition of ambient PM2.5 in China but have been constrained by the paucity of long-term ground measurements. Using the GEOS-Chem chemical transport model and satellite-derived PM2.5 data, sulfate and ammonium levels were estimated over China from 2004 to 2014. A comparison of the satellite-estimated dataset with model simulations based on ground measurements obtained from the literature indicated our results are more accurate. Using satellite-derived PM2.5 data with a spatial resolution of 0.1 × 0.1°, we further presented finer satellite-estimated sulfate and ammonium concentrations in anthropogenic polluted regions, including the NCP (the North China Plain), the SCB (the Sichuan Basin) and the PRD (the Pearl River Delta). Linear regression results obtained on a national scale yielded an r value of 0.62, NMB of −35.9%, NME of 48.2%, ARB_50% of 53.68% for sulfate and an r value of 0.63, slope of 0.67, and intercept of 5.14 for ammonium. In typical regions, the satellite-derived dataset was significantly robust. Based on the satellite-derived dataset, the spatial-temporal variation of 11-year annual average satellite-derived SO42− and NH4+ concentrations and time series of monthly average concentrations were also investigated. On a national scale, both exhibited a downward trend each year between 2004 and 2014 (SO42−: −0.61%; NH4+: −0.21%), large values were mainly concentrated in the NCP and SCB. For regions captured at a finer resolution, the inter-annual variation trends presented a positive trend over the periods 2004–2007 and 2008–2011, followed by a negative trend over the period 2012–2014, and sulfate concentrations varied appreciably. Moreover, the seasonal distributions of the 11-year satellite-derived dataset over China were presented. The distribution of both sulfate and ammonium concentrations exhibited seasonal characteristics, with the seasonal concentrations ranking as follows: winter > summer > autumn > spring. High concentrations of these species were concentrated in the NCP and SCB, originating from coal-fired power plants and agricultural activities, respectively. Efforts to reduce sulfur dioxide (SO2) emissions have yielded remarkable results since the government has adopted stricter control measures in recent years. Moreover, ammonia emissions should be controlled while reducing the concentration of sulfur, nitrogen and particulate matter. This study provides an assessment of the population’s exposure to certain chemical components.

Description: Remote Sensing, Vol. 9, Pages 837: National BDS Augmentation Service System (NBASS) of China: Progress and Assessment Remote Sensing doi: 10.3390/rs9080837 Authors: Chuang Shi Fu Zheng Yidong Lou Shengfeng Gu Weixing Zhang Xiaolei Dai Xianjie Li Hailin Guo Xiaopeng Gong Abstract: In this contribution, the processing strategies of real-time BeiDou System (BDS) precise orbits, clocks, and ionospheric corrections in the National BDS Augmentation Service System (NBASS) are briefly introduced. The Root Mean Square (RMS) of BDS predicted orbits are better than 10 cm in radial and cross-track components, and the accuracy of the BDS real-time clock is better than 0.5 ns for Inclined Geosynchronous Orbit (IGSO) and Mid Earth Orbit (MEO) satellites. The accuracy of BDS Geostationary Earth Orbit (GEO) orbits and clocks are worse than the IGSO and MEO satellites due to its poor geometry conditions. The real-time ionospheric correction is evaluated by cross-validation, and the average accuracy in the vertical direction is about 4 TECU. With these real-time corrections, the overall single and dual-frequency kinematic precise point positioning (PPP) performance in China are evaluated in terms of positioning accuracy at the 95% confidence level and convergence time. The BDS PPP positioning accuracy shows significant regional characteristics due to the geometry distribution of BDS satellites and the accuracy of ionospheric model in different regions. The BDS dual-frequency PPP positioning accuracy in high-latitude and western fringe region is about 0.5 m and 1.0 m in the horizontal and vertical component, respectively, while the horizontal accuracy is better than 0.2 m and the vertical accuracy is better than 0.3 m in the midlands. The convergence time of the BDS PPP is much longer than the GPS PPP and it needs more than 60 min to achieve the accuracy better than 10 cm in both horizontal and vertical directions for dual-frequency PPP. Similar with dual-frequency PPP, the positioning accuracy of the BDS single-frequency PPP in the fringe region is worse than other regions. The positioning in the midlands can achieve 0.5 m in horizontal component and 1.0 m in the vertical component. In addition, when GPS and BDS are combined, the positioning performance of both single-frequency and dual-frequency PPP can be greatly improved.

Description: Remote Sensing, Vol. 9, Pages 833: A Study of Landfast Ice with Sentinel-1 Repeat-Pass Interferometry over the Baltic Sea Remote Sensing doi: 10.3390/rs9080833 Authors: Marjan Marbouti Jaan Praks Oleg Antropov Eero Rinne Matti Leppäranta Mapping of fast ice displacement and investigating sea ice rheological behavior is a major open topic in coastal ice engineering and sea ice modeling. This study presents first results on Sentinel-1 repeat-pass space borne synthetic aperture radar interferometry (InSAR) in the Gulf of Bothnia over the fast ice areas. An InSAR pair acquired in February 2015 with a temporal baseline of 12 days has been studied here in detail. According to our results, the surface of landfast ice in the study area was stable enough to preserve coherence over the 12-day baseline, while previous InSAR studies over the fast ice used much shorter temporal baselines. The advantage of longer temporal baseline is in separating the fast ice from drift ice and detecting long term trends in deformation maps. The interferogram showed displacement of fast ice on the order of 40 cm in the study area. Parts of the displacements were attributed to forces caused by sea level tilt, currents, and thermal expansion, but the main factor of the displacement seemed to be due to compression of the drift ice driven by southwest winds with high speed. Further interferometric phase and the coherence measurements over the fast ice are needed in the future for understanding sea ice mechanism and establishing sustainability of the presented InSAR approach for monitoring dynamics of the landfast ice with Sentinel-1 data.

Description: IJGI, Vol. 6, Pages 249: Evaluating the Impact of Meteorological Factors on Water Demand in the Las Vegas Valley Using Time-Series Analysis: 1990–2014 ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6080249 Authors: Patcha Huntra Tim Keener Many factors impact a city’s water consumption, including population distribution, average household income, water prices, water conservation programs, and climate. Of these, however, meteorological effects are considered to be the primary determinants of water consumption. In this study, the effects of climate on residential water consumption in Las Vegas, Nevada, were examined during the period from 1990 to 2014. The investigations found that climatic variables, including maximum temperature, minimum temperature, average temperature, precipitation, diurnal temperature, dew point depression, wind speed, wind direction, and percent of calm wind influenced water use. The multivariate autoregressive integrated moving average (ARIMAX) model found that the historical data of water consumption and dew point depression explain the highest percentage of variance (98.88%) in water use when dew point depression is used as an explanatory variable. Our results indicate that the ARIMAX model with dew point depression input, and average temperature, play a significant role in predicting long-term water consumption rates in Las Vegas. The sensitivity analysis results also show that the changes in average temperature impacted water demand three times more than dew point depression. The accuracy performance, specifically the mean average percentage error (MAPE), of the model’s forecasting is found to be about 2–3% from five years out. This study can be adapted and utilized for the long-term forecasting of water demand in other regions. By using one significant climate factor and historical water demand for the forecasting, the ARIMAX model gives a forecast with high accuracy and provides an effective technique for monitoring the effects of climate change on water demand in the area.

Description: Remote Sensing, Vol. 9, Pages 845: Assimilation of Sentinel-1 Derived Sea Surface Winds for Typhoon Forecasting Remote Sensing doi: 10.3390/rs9080845 Authors: Yi Yu Xiaofeng Yang Weimin Zhang Boheng Duan Xiaoqun Cao Hongze Leng High-resolution synthetic aperture radar (SAR) wind observations provide fine structural information for tropical cycles and could be assimilated into numerical weather prediction (NWP) models. However, in the conventional method assimilating the u and v components for SAR wind observations (SAR_uv), the wind direction is not a state vector and its observational error is not considered during the assimilation calculation. In this paper, an improved method for wind observation directly assimilates the SAR wind observations in the form of speed and direction (SAR_sd). This method was implemented to assimilate the sea surface wind retrieved from Sentinel-1 synthetic aperture radar (SAR) in the basic three-dimensional variational system for the Weather Research and Forecasting Model (WRF 3DVAR). Furthermore, a new quality control scheme for wind observations is also presented. Typhoon Lionrock in August 2016 is chosen as a case study to investigate and compare both assimilation methods. The experimental results show that the SAR wind observations can increase the number of the effective observations in the area of a typhoon and have a positive impact on the assimilation analysis. The numerical forecast results for this case show better results for the SAR_sd method than for the SAR_uv method. The SAR_sd method looks very promising for winds assimilation under typhoon conditions, but more cases need to be considered to draw final conclusions.

Description: Remote Sensing, Vol. 9, Pages 840: Topic Modelling for Object-Based Unsupervised Classification of VHR Panchromatic Satellite Images Based on Multiscale Image Segmentation Remote Sensing doi: 10.3390/rs9080840 Authors: Li Shen Linmei Wu Yanshuai Dai Wenfan Qiao Ying Wang Image segmentation is a key prerequisite for object-based classification. However, it is often difficult, or even impossible, to determine a unique optimal segmentation scale due to the fact that various geo-objects, and even an identical geo-object, present at multiple scales in very high resolution (VHR) satellite images. To address this problem, this paper presents a novel unsupervised object-based classification for VHR panchromatic satellite images using multiple segmentations via the latent Dirichlet allocation (LDA) model. Firstly, multiple segmentation maps of the original satellite image are produced by means of a common multiscale segmentation technique. Then, the LDA model is utilized to learn the grayscale histogram distribution for each geo-object and the mixture distribution of geo-objects within each segment. Thirdly, the histogram distribution of each segment is compared with that of each geo-object using the Kullback-Leibler (KL) divergence measure, which is weighted with a constraint specified by the mixture distribution of geo-objects. Each segment is allocated a geo-object category label with the minimum KL divergence. Finally, the final classification map is achieved by integrating the multiple classification results at different scales. Extensive experimental evaluations are designed to compare the performance of our method with those of some state-of-the-art methods for three different types of images. The experimental results over three different types of VHR panchromatic satellite images demonstrate the proposed method is able to achieve scale-adaptive classification results, and improve the ability to differentiate the geo-objects with spectral overlap, such as water and grass, and water and shadow, in terms of both spatial consistency and semantic consistency.

Description: IJGI, Vol. 6, Pages 252: A Novel Approach for Publishing Linked Open Geodata from National Registries with the Use of Semantically Annotated Context Dependent Web Pages ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6080252 Authors: Adam Iwaniak Marta Leszczuk Marek Strzelecki Francis Harvey Iwona Kaczmarek Many of the standards used to build spatial data infrastructure (SDI), such as Web Map Service (WMS) or Web Feature Service (WFS), have become outdated. They do not follow current web technology development and do not fully exploit its capabilities. Spatial data often remains available only through application programming interfaces (APIs), reflecting the persistence of organizational silos. The potential of the web for discovering knowledge hidden in data and discoverable through integration and fusion remains very difficult. This article presents a strategy to take advantage of these newer semantic web technologies for SDI. We describe the implementation of a public registry in the age of Web 3.0. Our goal is to convert existing geographic information systems (GIS) data into explicit knowledge that can be easily used for a variety of purposes. This turns SDI into a framework to utilize the many advantages of the web. In this paper we present the working prototype system developed for the province of Mazowieckie in Poland and describes the underlying concepts. Further development of this approach comes from using linked data (LD) with expert systems to support analysis functions and tasks.

Description: Remote Sensing, Vol. 9, Pages 851: Reflectance Intensity Assisted Automatic and Accurate Extrinsic Calibration of 3D LiDAR and Panoramic Camera Using a Printed Chessboard Remote Sensing doi: 10.3390/rs9080851 Authors: Weimin Wang Ken Sakurada Nobuo Kawaguchi This paper presents a novel method for fully automatic and convenient extrinsic calibration of a 3D LiDAR and a panoramic camera with a normally printed chessboard. The proposed method is based on the 3D corner estimation of the chessboard from the sparse point cloud generated by one frame scan of the LiDAR. To estimate the corners, we formulate a full-scale model of the chessboard and fit it to the segmented 3D points of the chessboard. The model is fitted by optimizing the cost function under constraints of correlation between the reflectance intensity of laser and the color of the chessboard’s patterns. Powell’s method is introduced for resolving the discontinuity problem in optimization. The corners of the fitted model are considered as the 3D corners of the chessboard. Once the corners of the chessboard in the 3D point cloud are estimated, the extrinsic calibration of the two sensors is converted to a 3D-2D matching problem. The corresponding 3D-2D points are used to calculate the absolute pose of the two sensors with Unified Perspective-n-Point (UPnP). Further, the calculated parameters are regarded as initial values and are refined using the Levenberg-Marquardt method. The performance of the proposed corner detection method from the 3D point cloud is evaluated using simulations. The results of experiments, conducted on a Velodyne HDL-32e LiDAR and a Ladybug3 camera under the proposed re-projection error metric, qualitatively and quantitatively demonstrate the accuracy and stability of the final extrinsic calibration parameters.

Description: Remote Sensing, Vol. 9, Pages 850: What is the Direction of Land Change? A New Approach to Land-Change Analysis Remote Sensing doi: 10.3390/rs9080850 Authors: Mingde You Anthony Filippi İnci Güneralp Burak Güneralp Accurate characterization of the direction of land change is a neglected aspect of land dynamics. Knowledge on direction of historical land change can be useful information when understanding relative influence of different land-change drivers is of interest. In this study, we present a novel perspective on land-change analysis by focusing on directionality of change. To this end, we employed Maximum Cross-Correlation (MCC) approach to estimate the directional change in land cover in a dynamic river floodplain environment using Landsat 5 Thematic Mapper (TM) images. This approach has previously been used for detecting and measuring fluid and ice motions but not to study directional changes in land cover. We applied the MCC approach on land-cover class membership layers derived from fuzzy remote-sensing image classification. We tested the sensitivity of the resulting displacement vectors to three user-defined parameters—template size, search window size, and a threshold parameter to determine valid (non-noisy) displacement vectors—that directly affect the generation of change, or displacement, vectors; this has not previously been thoroughly investigated in any application domain. The results demonstrate that it is possible to quantitatively measure the rate of directional change in land cover in this floodplain environment using this particular approach. Sensitivity analyses indicate that template size and MCC threshold parameter are more influential on the displacement vectors than search window size. The results vary by land-cover class, suggesting that spatial configuration of land-cover classes should be taken into consideration in the implementation of the method.

Description: IJGI, Vol. 6, Pages 253: Estimation of Travel Time Distributions in Urban Road Networks Using Low-Frequency Floating Car Data ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6080253 Authors: Chaoyang Shi Bi Chen Qingquan Li Travel times in urban road networks are highly stochastic. However, most existing travel time estimation methods only estimate the mean travel times, while ignoring travel time variances. To this end, this paper proposes a robust travel time distribution estimation method to estimate both the mean and variance of travel times by using emerging low-frequency floating car data. Different from the existing studies, the path travel time distribution in this study is formulated as the sum of the deterministic link travel times and stochastic turning delays at intersections. Using this formulation, distinct travel time delays for different turning movements at the same intersection can be well captured. In this study, a speed estimation algorithm is developed to estimate the deterministic link travel times, and a distribution estimation algorithm is proposed to estimate the stochastic turning delays. Considering the low sampling rate of the floating car data, a weighted moving average algorithm is further developed for a robust estimation of the path travel time distribution. A real-world case study in Wuhan, China is carried out to validate the applicability of the proposed method. The results of the case study show that the proposed method can obtain a reliable and accurate estimation of path travel time distribution in congested urban road networks.

Description: IJGI, Vol. 6, Pages 258: Spatial-Spectral Graph Regularized Kernel Sparse Representation for Hyperspectral Image Classification ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6080258 Authors: Jianjun Liu Zhiyong Xiao Yufeng Chen Jinlong Yang This paper presents a spatial-spectral method for hyperspectral image classification in the regularization framework of kernel sparse representation. First, two spatial-spectral constraint terms are appended to the sparse recovery model of kernel sparse representation. The first one is a graph-based spatially-smooth constraint which is utilized to describe the contextual information of hyperspectral images. The second one is a spatial location constraint, which is exploited to incorporate the prior knowledge of the location information of training pixels. Then, an efficient alternating direction method of multipliers is developed to solve the corresponding minimization problem. At last, the recovered sparse coefficient vectors are used to determine the labels of test pixels. Experimental results carried out on three real hyperspectral images point out the effectiveness of the proposed method.

Description: Remote Sensing, Vol. 9, Pages 867: The Effects of Aerosol on the Retrieval Accuracy of NO2 Slant Column Density Remote Sensing doi: 10.3390/rs9080867 Authors: Hyunkee Hong Jhoon Kim Ukkyo Jeong Kyung-soo Han Hanlim Lee We investigate the effects of aerosol optical depth (AOD), single scattering albedo (SSA), aerosol peak height (APH), measurement geometry (solar zenith angle (SZA) and viewing zenith angle (VZA)), relative azimuth angle, and surface reflectance on the accuracy of NO2 slant column density using synthetic radiance. High AOD and APH are found to decrease NO2 SCD retrieval accuracy. In moderately polluted (5 × 1015 molecules cm−2 < NO2 vertical column density (VCD) < 2 × 1016 molecules cm−2) and clean regions (NO2 VCD < 5 × 1015 molecules cm−2), the correlation coefficient (R) between true NO2 SCDs and those retrieved is 0.88 and 0.79, respectively, and AOD and APH are about 0.1 and is 0 km, respectively. However, when AOD and APH are about 1.0 and 4 km, respectively, the R decreases to 0.84 and 0.53 in moderately polluted and clean regions, respectively. On the other hand, in heavily polluted regions (NO2 VCD > 2 × 1016 molecules cm−2), even high AOD and APH values are found to have a negligible effect on NO2 SCD precision. In high AOD and APH conditions in clean NO2 regions, the R between true NO2 SCDs and those retrieved increases from 0.53 to 0.58 via co-adding four pixels spatially, showing the improvement in accuracy of NO2 SCD retrieval. In addition, the high SZA and VZA are also found to decrease the accuracy of the NO2 SCD retrieval.

Description: IJGI, Vol. 6, Pages 259: Identifying and Analyzing the Prevalent Regions of a Co-Location Pattern Using Polygons Clustering Approach ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090259 Authors: Wenhao Yu Given a co-location pattern consisting of spatial features, the prevalent region mining process identifies local areas in which these features are co-located with a high probability. Many approaches have been proposed for co-location mining due to its key role in public safety, social-economic development and environmental management. However, traditionally, most of the solutions focus on itemsets mining and results outputting in a textual format, which fail to adequately treat all the spatial nature of the underlying entities and processes. In this paper, we propose a new co-location analysis approach to find the prevalent regions of a pattern. The approach combines kernel density estimation and polygons clustering techniques to specifically consider the correlation, heterogeneity and contextual information existing within complex spatial interactions. A kernel density estimation surface is created for each feature and subsequently the generated multiple surfaces are combined into a final surface with cell attribute representing the pattern prevalence measure value. Polygons consisting of cells are then extracted according to the predefined threshold. Through adding appended environmental data to the polygons, an outcome of similar groups is achieved using polygons clustering approach. The effectiveness of our approach is evaluated using Points-of-Interest datasets in Shenzhen, China.

Description: Remote Sensing, Vol. 9, Pages 880: A Microtopographic Feature Analysis-Based LiDAR Data Processing Approach for the Identification of Chu Tombs Remote Sensing doi: 10.3390/rs9090880 Authors: Shaohua Wang Qingwu Hu Fengzhu Wang Mingyao Ai Ruofei Zhong Most of the cultural sites hidden under dense vegetation in the mountains of China have been destroyed. In this paper, we present a microtopographic feature analysis (MFA)-based Light Detection and Ranging (LiDAR) data processing approach and an archaeological pattern-oriented point cloud segmentation (APoPCS) algorithm that we developed for the classification of archaeological objects and terrain points and the detection of archaeological remains. The archaeological features and patterns are interpreted and extracted from LiDAR point cloud data to construct an archaeological object pattern database. A microtopographic factor is calculated based on the archaeological object patterns, and this factor converts the massive point cloud data into a raster feature image. A fuzzy clustering algorithm based on the archaeological object patterns is presented for raster feature image segmentation and the detection of archaeological remains. Using the proposed approach, we investigated four typical areas with different types of Chu tombs in Central China, which had dense vegetation and high population densities. Our research results show that the proposed LiDAR data processing approach can identify archaeological remains from large-volume and massive LiDAR data, as well as in areas with dense vegetation and trees. The studies of different archaeological object patterns are important for improving the robustness of the proposed APoPCS algorithm for the extraction of archaeological remains.

Description: IJGI, Vol. 6, Pages 272: An Adaptive Sweep-Circle Spatial Clustering Algorithm Based on Gestalt ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090272 Authors: Qingming Zhan Shuguang Deng Zhihua Zheng An adaptive spatial clustering (ASC) algorithm is proposed in this present study, which employs sweep-circle techniques and a dynamic threshold setting based on the Gestalt theory to detect spatial clusters. The proposed algorithm can automatically discover clusters in one pass, rather than through the modification of the initial model (for example, a minimal spanning tree, Delaunay triangulation, or Voronoi diagram). It can quickly identify arbitrarily-shaped clusters while adapting efficiently to non-homogeneous density characteristics of spatial data, without the need for prior knowledge or parameters. The proposed algorithm is also ideal for use in data streaming technology with dynamic characteristics flowing in the form of spatial clustering in large data sets.

Description: IJGI, Vol. 6, Pages 271: Contextual Building Selection Based on a Genetic Algorithm in Map Generalization ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090271 Authors: Lin Wang Qingsheng Guo Yuangang Liu Yageng Sun Zhiwei Wei In map generalization, scale reduction and feature symbolization inevitably generate problems of overlapping objects or map congestion. To solve the legibility problem with respect to the generalization of dispersed rural buildings, selection of buildings is necessary and can be transformed into an optimization problem. In this paper, an improved genetic algorithm for building selection is designed to be able to incorporate cartographic constraints related to the building selection problem. Part of the local constraints for building selection is used to constrain the encoding and genetic operation. To satisfy other local constraints, a preparation phase is necessary before building selection, which includes building enlargement, local displacement, conflict detection, and attribute enrichment. The contextual constraints are used to ascertain a fitness function. The experimental results indicate that the algorithm proposed in this article can obtain good results for building selection whilst preserving the spatial distribution characteristics of buildings.

Description: Remote Sensing, Vol. 9, Pages 898: Characteristics of Aerosol Types in Beijing and the Associations with Air Pollution from 2004 to 2015 Remote Sensing doi: 10.3390/rs9090898 Authors: Yang Ou Wenhui Zhao Junqian Wang Wenji Zhao Bo Zhang With the fast development of the economy and expansion, a large number of people have concentrated in Beijing over the past few decades, leading to the result that Beijing has become home to one of the most complex mixtures of aerosol types in the world. The various aerosol types play different roles in the determination of global climate change, visibility, and human health. However, to the best of our knowledge, research has rarely analyzed the correlation between aerosol types and air quality index (AQI) in Beijing (urban and suburban) over a long-term series of observations. Therefore, in this study, we aim to identify and discuss the different aerosol types and AQI in Beijing from 2004 to 2015. The aerosol types are classified into six categories: dust, mixed, highly-absorbing, moderately-absorbing, slightly-absorbing, and scattering by a multiple clustering method with the fine mode fraction (FMF) and single scattering albedo (SSA) data of retrievals from the global Aerosol Robotic Network (AERONET) sun photometer sites. The AQI levels: are good (0–50); moderate (51–100); unhealthy for sensitive groups (101–150); unhealthy (151–200); very unhealthy (201–300); and hazardous (>300). The results show that a significant FMF variability occurred among different seasons in Beijing, with maximum values present in spring and minimum values in winter. The SSA values exhibit variation, with small fluctuations from season to season. In the case of BJ station, the scattering aerosols are more frequent in summer (39%) and less in winter (1%), while the coarse particles (dust) are more frequent in spring (18%) and less in autumn (6%). In contrast, the absorbing aerosols (especially slightly-absorbing) are more frequent in summer (35%) and winter (15%). However, the mixed aerosol types are more frequent in spring (38%) and less in summer (8%). There is a similar seasonal variation in XH. In the past 12 years, the slightly-absorbing aerosol type in Beijing has increased by approximately 14%, which is believed to be due to the rapid development of industrial cities. In addition, comparing the urban and suburban regions, the slightly-absorbing aerosol type is the dominant aerosol type in both areas. Furthermore, to identify the dominant aerosol types which lead to air pollution, a related analysis was carried out by analyzing different aerosol types and the relationship between aerosol types and AQI. The results indicate that the air pollution was strongly correlated to slightly-absorbing aerosols, in which the percentage of slightly-absorbing aerosols was about 49% during the hazardous days in 2013–2015, and the correlation between AQI and aerosol types is also strong (R2 = 0.76 and 0.97, in Beijing and Xianghe).

Description: Remote Sensing, Vol. 9, Pages 897: Elevation Extraction and Deformation Monitoring by Multitemporal InSAR of Lupu Bridge in Shanghai Remote Sensing doi: 10.3390/rs9090897 Authors: Jingwen Zhao Jicang Wu Xiaoli Ding Mingzhou Wang Monitoring, assessing, and understanding the structural health of large infrastructures, such as buildings, bridges, dams, tunnels, and highways, is important for urban development and management, as the gradual deterioration of such structures may result in catastrophic structural failure leading to high personal and economic losses. With a higher spatial resolution and a shorter revisit period, interferometric synthetic aperture radar (InSAR) plays an increasing role in the deformation monitoring and height extraction of structures. As a focal point of the InSAR data processing chain, phase unwrapping has a direct impact on the accuracy of the results. In complex urban areas, large elevation differences between the top and bottom parts of a large structure combined with a long interferometric baseline can result in a serious phase-wrapping problem. Here, with no accurate digital surface model (DSM) available, we handle the large phase gradients of arcs in multitemporal InSAR processing using a long–short baseline iteration method. Specifically, groups of interferometric pairs with short baselines are processed to obtain the rough initial elevation estimations of the persistent scatterers (PSs). The baseline threshold is then loosened in subsequent iterations to improve the accuracy of the elevation estimates step by step. The LLL lattice reduction algorithm (by Lenstra, Lenstra, and Lovász) is applied in the InSAR phase unwrapping process to rapidly reduce the search radius, compress the search space, and improve the success rate in resolving the phase ambiguities. Once the elevations of the selected PSs are determined, they are used in the following two-dimensional phase regression involving both elevations and deformations. A case study of Lupu Bridge in Shanghai is carried out for the algorithm’s verification. The estimated PS elevations agree well (within 1 m) with the official Lupu Bridge model data, while the PS deformation time series confirms that the bridge exhibits some symmetric progressive deformation, at 4–7 mm per year on both arches and 4–9 mm per year on the bridge deck during the SAR image acquisition period.

Description: Remote Sensing, Vol. 9, Pages 907: Transfer Learning with Deep Convolutional Neural Network for SAR Target Classification with Limited Labeled Data Remote Sensing doi: 10.3390/rs9090907 Authors: Zhongling Huang Zongxu Pan Bin Lei Tremendous progress has been made in object recognition with deep convolutional neural networks (CNNs), thanks to the availability of large-scale annotated dataset. With the ability of learning highly hierarchical image feature extractors, deep CNNs are also expected to solve the Synthetic Aperture Radar (SAR) target classification problems. However, the limited labeled SAR target data becomes a handicap to train a deep CNN. To solve this problem, we propose a transfer learning based method, making knowledge learned from sufficient unlabeled SAR scene images transferrable to labeled SAR target data. We design an assembled CNN architecture consisting of a classification pathway and a reconstruction pathway, together with a feedback bypass additionally. Instead of training a deep network with limited dataset from scratch, a large number of unlabeled SAR scene images are used to train the reconstruction pathway with stacked convolutional auto-encoders (SCAE) at first. Then, these pre-trained convolutional layers are reused to transfer knowledge to SAR target classification tasks, with feedback bypass introducing the reconstruction loss simultaneously. The experimental results demonstrate that transfer learning leads to a better performance in the case of scarce labeled training data and the additional feedback bypass with reconstruction loss helps to boost the capability of classification pathway.

Description: Remote Sensing, Vol. 9, Pages 900: Irrigation-Induced Environmental Changes around the Aral Sea: An Integrated View from Multiple Satellite Observations Remote Sensing doi: 10.3390/rs9090900 Authors: Qinjian Jin Jiangfeng Wei Zong-Liang Yang Peirong Lin The Aral Sea basin (ASB) is one of the most environmentally vulnerable regions to climate change and human activities. During the past 60 years, irrigation has greatly changed the water distribution and caused severe environmental issues in the ASB. Using remote sensing data, this study investigated the environmental changes induced by irrigation activities in this region. The results show that, in the past decade, land water storage has significantly increased in the irrigated upstream regions (13 km3 year−1) but decreased in the downstream regions (−27 km3 year−1) of the Amu Darya River basin, causing a water storage decrease in the whole basin (−20 km3 year−1). As a result, the water surface area of the Aral Sea has decreased from 32,000 in 2000 to 10,000 km2 in 2015. The shrinking Aral Sea exposed a large portion of the lake bottom to the air, increasing (decreasing) the daytime (nighttime) temperatures by about 1 °C year−1 (0.5 °C year−1). Moreover, there were other potential environmental changes, including drier soil, less vegetation, decreasing cloud and precipitation, and more severe and frequent dust storms. Possible biases in the remote sensing data due to the neglect of the shrinking water surface area of the Aral Sea were identified. These findings highlight the severe environmental threats caused by irrigation in Central Asia and call attention to sustainable water use in such dryland regions.

Description: Remote Sensing, Vol. 9, Pages 906: Evaluation of Sentinel-2A Satellite Imagery for Mapping Cotton Root Rot Remote Sensing doi: 10.3390/rs9090906 Authors: Xiaoyu Song Chenghai Yang Mingquan Wu Chunjiang Zhao Guijun Yang Wesley Hoffmann Wenjiang Huang Cotton (Gossypium hirsutum L.) is an economically important crop that is highly susceptible to cotton root rot. Remote sensing technology provides a useful and effective means for detecting and mapping cotton root rot infestations in cotton fields. This research assessed the potential of 10-m Sentinel-2A satellite imagery for cotton root rot detection and compared it with airborne multispectral imagery using unsupervised classification at both field and regional levels. Accuracy assessment showed that the classification maps from the Sentinel-2A imagery had an overall accuracy of 94.1% for field subset images and 91.2% for the whole image, compared with the airborne image classification results. However, some small cotton root rot areas were undetectable and some non-infested areas within large root rot areas were incorrectly classified as infested due to the images’ coarse spatial resolution. Classification maps based on field subset Sentinel-2A images missed 16.6% of the infested areas and the classification map based on the whole Sentinel-2A image for the study area omitted 19.7% of the infested areas. These results demonstrate that freely-available Sentinel-2 imagery can be used as an alternative data source for identifying cotton root rot and creating prescription maps for site-specific management of the disease.

Description: IJGI, Vol. 6, Pages 274: Using Eye Tracking to Explore the Guidance and Constancy of Visual Variables in 3D Visualization ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090274 Authors: Bing Liu Weihua Dong Liqiu Meng An understanding of guidance, which means guiding attention, and constancy, meaning that an area can be perceived for what it is despite environmental changes, of the visual variables related to three-dimensional (3D) symbols is essential to ensure rapid and consistent human perception in 3D visualization. Previous studies have focused on the guidance and constancy of visual variables related to two-dimensional (2D) symbols, but these aspects are not well documented for 3D symbols. In this study, we used eye tracking to analyze the visual guidance from shapes, hues and sizes, and the visual constancy that is related to the shape, color saturation and size of 3D symbols in different locations. Thirty-six subjects (24 females and 12 males) participated in the study. The results indicate that hue and shape provide a high level of visual guidance, whereas guidance from size, a variable that predominantly guides attention in 2D visualization, is much more limited in 3D visualization. Additionally, constancy of shape and saturation are perceived with relatively high accuracy, whereas constancy of size is perceived with only low accuracy. These first empirical studies are intended to pave the way for a more comprehensive user understanding of 3D visualization design.

Description: IJGI, Vol. 6, Pages 277: Closing Data Gaps with Citizen Science? Findings from the Danube Region ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090277 Authors: Josip Lisjak Sven Schade Alexander Kotsev Although data is increasingly shared online and accessible for re-use, we still witness heterogeneous coverage of thematic areas and geographic regions. This especially becomes an issue when data is needed for large territories and including different nations, as, for example, required to support macro-regional development policies. Once identified, data gaps might be closed using different approaches. Existing—but so far non accessible—data might be made available; new public sector information could be gathered; or data might be acquired from the private sector. Our work explores a fourth option: closing data gaps with direct contributions from citizen (Citizen Science). This work summarizes a particular case study that was conducted in 2016 in the Danube Region. We provide a gap analysis over an existing macro-regional data infrastructure, and examine potential Citizen Science approaches that might help to close these gaps. We highlight already existing Citizen Science projects that could address a large part of the identified gaps, and suggest one particular new application in order to indicate how a—so far uncovered—gap might be approached. This new application addresses bioenergy as a particular field of the circular economy. On this basis we discuss the emerging opportunities and challenges for this particular way of public participation in regional development policy. We close by highlighting areas for future research.

Description: Remote Sensing, Vol. 9, Pages 916: “Kill Two Birds with One Stone”: Urban Tree Species Classification Using Bi-Temporal Pléiades Images to Study Nesting Preferences of an Invasive Bird Remote Sensing doi: 10.3390/rs9090916 Authors: Marine Le Louarn Philippe Clergeau Elodie Briche Magali Deschamps-Cottin This study presents the results of object-based classifications assessing the potential of bi-temporal Pléiades images for mapping broadleaf and coniferous tree species potentially used by the ring-necked parakeet Psittacula krameri for nesting in the urban area of Marseille, France. The first classification was performed based solely on a summer Pléiades image (acquired on 28 July 2015) and the second classification based on bi-temporal Pléiades images (a spring image acquired on 24 March 2016 and the summer image). An ensemble of spectral and textural features was extracted from both images and two machine-learning classifiers were used, Random Forest (RF) and Support Vector Machine (SVM). Regardless of the classifiers, model results suggest that classification based on bi-temporal Pléiades images produces more satisfying results, with an overall accuracy 11.5–13.9% higher than classification using the single-date image. Textural and spectral features extracted from the blue and the NIR bands were consistently ranked among the most important features. Regardless of the classification scheme, RF slightly outperforms SVM. RF classification using bi-temporal Pléiades images allows identifying 98.5% of the tree species used by the ring-necked parakeet for nesting, highlighting the promising value of remote sensing techniques to assess the ecological requirements of fauna in urban areas.

Description: Remote Sensing, Vol. 9, Pages 919: Diversification of Land Surface Temperature Change under Urban Landscape Renewal: A Case Study in the Main City of Shenzhen, China Remote Sensing doi: 10.3390/rs9090919 Authors: Yanxu Liu Jian Peng Yanglin Wang Unprecedented rapid urbanization in China during the past several decades has been accompanied by extensive urban landscape renewal, which has increased the urban thermal environmental risk. However, landscape change is a sufficient but not necessary condition for land surface temperature (LST) variation. Many studies have merely highlighted the correlation between landscape pattern and LST, while neglecting to comprehensively present the spatiotemporal diversification of LST change under urban landscape renewal. Taking the main city of Shenzhen as a case study area, this study tracked the landscape renewal and LST variation for the period 1987–2015 using 49 Landsat images. A decision tree algorithm suitable for fast landscape type interpretation was developed to map the landscape renewal. Analytical tools that identified hot-cold spots, the gravity center, and transect of LST movement were adopted to identify LST changes. The results showed that the spatial variation of LST was not completely consistent with landscape change. The transformation from Green landscape to Grey landscape usually increased the LST within a median of 0.2 °C, while the reverse transformation did not obviously decrease the LST (the median was nearly 0 °C). The median of LST change from Blue landscape to Grey landscape was 1.0 °C, corresponding to 0.5 °C in the reverse transformation. The imbalance of LST change between the loss and gain of Green or Blue landscape indicates the importance of protecting natural space, where the benefits in terms of temperature mitigation cannot be completely substituted by reverse transformation.

Description: Remote Sensing, Vol. 9, Pages 920: Understanding How Low-Level Clouds and Fog Modify the Diurnal Cycle of Orographic Precipitation Using In Situ and Satellite Observations Remote Sensing doi: 10.3390/rs9090920 Authors: Yajuan Duan Ana Barros Satellite orographic precipitation estimates exhibit large errors with space-time structure tied to landform. Observations in the Southern Appalachian Mountains (SAM) suggest that low-level clouds and fog (LLCF) amplify mid-day rainfall via seeder-feeder interactions (SFI) at both high and low elevations. Here, a rainfall microphysics model constrained by fog observations was used first to reveal that fast SFI (2–5 min time-scales) modify the rain drop size distributions by increasing coalescence efficiency among small drops (<0.7 mm diameter), whereas competition between coalescence and filament-only breakup dominates for larger drops (3–5 mm diameter). The net result is a large increase in the number concentrations of intermediate size raindrops in the 0.7–3 mm range and up to a ten-fold increase in rainfall intensity. Next, a 10-year climatology of satellite observations was developed to map LLCF. Combined estimates from CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) and CloudSat products reveal persistent shallower cloud base heights at high elevations enveloping the terrain. The regional cloud top height climatology derived from the MODIS (Moderate Resolution Imaging Spectroradiometer) shows high-frequency daytime LLCF over mountain ridges in the warm season shifting to river valleys at nighttime. In fall and winter, LLCF patterns define a cloud-shadow region east of the continental divide, consistent with downwind rain-shadow effects. Optical and microphysical properties from collocated MODIS and ground ceilometers indicate small values of vertically integrated cloud water path (CWP < 100 g/m2), optical thickness (COT < 15), and particle effective radius (CER) < 15 μm near cloud top whereas surface observed CER ~25 μm changes to ~150 μm and higher prior to the mid-day rainfall. The vertical stratification of LLCF microphysics and SFI at low levels pose a significant challenge to satellite-based remote sensing in complex topography.

Description: Remote Sensing, Vol. 9, Pages 918: Scale- and Region-Dependence in Landscape-PM2.5 Correlation: Implications for Urban Planning Remote Sensing doi: 10.3390/rs9090918 Authors: Huihui Feng Bin Zou Yumeng Tang Under rapid urbanization, many cities in China suffer from serious fine particulate matter (PM2.5) pollution. As the emission sources or adsorption sinks, land use and the corresponding landscape pattern unavoidably affect the concentration. However, the correlation varies with different regions and scales, leaving a significant gap for urban planning. This study clarifies the correlation with the aid of in situ and satellite-based spatial datasets over six urban agglomerations in China. Two coverage and four landscape indices are adopted to represent land use and landscape pattern. Specifically, the coverage indices include the area ratios of forest (F_PLAND) and built-up areas (C_PLAND). The landscape indices refer to the perimeter-area fractal dimension index (PAFRAC), interspersion and juxtaposition index (IJI), aggregation index (AI), Shannon’s diversity index (SHDI). Then, the correlation between PM2.5 concentration with the selected indices are evaluated from supporting the potential urban planning. Results show that the correlations are weak with the in situ PM2.5 concentration, which are significant with the regional value. It means that land use coverage and landscape pattern affect PM2.5 at a relatively large scale. Furthermore, regional PM2.5 concentration negatively correlate to F_PLAND and positively to C_PLAND (significance at p < 0.05), indicating that forest helps to improve air quality, while built-up areas worsen the pollution. Finally, the heterogeneous landscape presents positive correlation to the regional PM2.5 concentration in most regions, except for the urban agglomeration with highly-developed urban (i.e., the Jing-Jin-Ji and Chengdu-Chongqing urban agglomerations). It suggests that centralized urbanization would be helpful for PM2.5 pollution controlling by reducing the emission sources in most regions. Based on the results, the potential urban planning is proposed for controlling PM2.5 pollution for each urban agglomeration.

Description: IJGI, Vol. 6, Pages 275: Mapping Forest Species in the Central Middle Atlas of Morocco (Azrou Forest) through Remote Sensing Techniques ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090275 Authors: Meriame Mohajane Ali Essahlaoui Fatiha Oudija Mohammed El Hafyani Ana Cláudia Teodoro The studies of forest ecosystems from remotely-sensed data are of great interest to researchers because of ecosystem services provided by this ecosystem, including protection of soils and vegetation, climate stabilization, and regulation of the hydrological cycle. In this context, our study focused on the use of a spectral angle mapper (SAM) classification method for mapping species in the Azrou Forest, Central Middle Atlas, Morocco. A Sentinel-2A image combined with ground reference data were used in this research. Four classes (holm oak, cedar forest, bare soil, and others-unclassified) were selected; they represent, respectively, 27, 11, 24, and 38% of the study area. The overall accuracy of classification was estimated to be around 99.72%. This work explored the potential of the SAM classification combined with Sentinel-2A data for mapping land cover in the Azrou Forest ecosystem.

Description: Remote Sensing, Vol. 9, Pages 923: Assessing Lodging Severity over an Experimental Maize (Zea mays L.) Field Using UAS Images Remote Sensing doi: 10.3390/rs9090923 Authors: Tianxing Chu Michael Starek Michael Brewer Seth Murray Luke Pruter Lodging has been recognized as one of the major destructive factors for crop quality and yield, resulting in an increasing need to develop cost-efficient and accurate methods for detecting crop lodging in a routine manner. Using structure-from-motion (SfM) and novel geospatial computing algorithms, this study investigated the potential of high resolution imaging with unmanned aircraft system (UAS) technology for detecting and assessing lodging severity over an experimental maize field at the Texas A&M AgriLife Research and Extension Center in Corpus Christi, Texas, during the 2016 growing season. The method was proposed to not only detect the occurrence of lodging at the field scale, but also to quantitatively estimate the number of lodged plants and the lodging rate within individual rows. Nadir-view images of the field trial were taken by multiple UAS platforms equipped with consumer grade red, green, and blue (RGB), and near-infrared (NIR) cameras on a routine basis, enabling a timely observation of the plant growth until harvesting. Models of canopy structure were reconstructed via an SfM photogrammetric workflow. The UAS-estimated maize height was characterized by polygons developed and expanded from individual row centerlines, and produced reliable accuracy when compared against field measures of height obtained from multiple dates. The proposed method then segmented the individual maize rows into multiple grid cells and determined the lodging severity based on the height percentiles against preset thresholds within individual grid cells. From the analysis derived from this method, the UAS-based lodging results were generally comparable in accuracy to those measured by a human data collector on the ground, measuring the number of lodging plants (R2 = 0.48) and the lodging rate (R2 = 0.50) on a per-row basis. The results also displayed a negative relationship of ground-measured yield with UAS-estimated and ground-measured lodging rate.

Description: Remote Sensing, Vol. 9, Pages 931: Mapping Smallholder Yield Heterogeneity at Multiple Scales in Eastern Africa Remote Sensing doi: 10.3390/rs9090931 Authors: Zhenong Jin George Azzari Marshall Burke Stephen Aston David Lobell Accurate measurements of crop production in smallholder farming systems are critical to the understanding of yield constraints and, thus, setting the appropriate agronomic investments and policies for improving food security and reducing poverty. Nevertheless, mapping the yields of smallholder farms is challenging because of factors such as small field sizes and heterogeneous landscapes. Recent advances in fine-resolution satellite sensors offer promise for monitoring and characterizing the production of smallholder farms. In this study, we investigated the utility of different sensors, including the commercial Skysat and RapidEye satellites and the publicly accessible Sentinel-2, for tracking smallholder maize yield variation throughout a ~40,000 km2 western Kenya region. We tested the potential of two types of multiple regression models for predicting yield: (i) a “calibrated model”, which required ground-measured yield and weather data for calibration, and (ii) an “uncalibrated model”, which used a process-based crop model to generate daily vegetation index and end-of-season biomass and/or yield as pseudo training samples. Model performance was evaluated at the field, division, and district scales using a combination of farmer surveys and crop cuts across thousands of smallholder plots in western Kenya. Results show that the “calibrated” approach captured a significant fraction (R2 between 0.3 and 0.6) of yield variations at aggregated administrative units (e.g., districts and divisions), while the “uncalibrated” approach performed only slightly worse. For both approaches, we found that predictions using the MERIS Terrestrial Chlorophyll Index (MTCI), which included the red edge band available in RapidEye and Sentinel-2, were superior to those made using other commonly used vegetation indices. We also found that multiple refinements to the crop simulation procedures led to improvements in the “uncalibrated” approach. We identified the prevalence of small field sizes, intercropping management, and cloudy satellite images as major challenges to improve the model performance. Overall, this study suggested that high-resolution satellite imagery can be used to map yields of smallholder farming systems, and the methodology presented in this study could serve as a good foundation for other smallholder farming systems in the world.

Description: Remote Sensing, Vol. 9, Pages 932: Hindcasting and Forecasting of Surface Flow Fields through Assimilating High Frequency Remotely Sensing Radar Data Remote Sensing doi: 10.3390/rs9090932 Authors: Lei Ren Michael Hartnett In order to improve the forecasting ability of numerical models, a sequential data assimilation scheme, nudging, was applied to blend remotely sensing high-frequency (HF) radar surface currents with results from a three-dimensional numerical, EFDC (Environmental Fluid Dynamics Code) model. For the first time, this research presents the most appropriate nudging parameters, which were determined from sensitivity experiments. To examine the influence of data assimilation cycle lengths on forecasts and to extend forecasting improvements, the duration of data assimilation cycles was studied through assimilating linearly interpolated temporal radar data. Data assimilation nudging parameters have not been previously analyzed. Assimilation of HF radar measurements at each model computational timestep outperformed those assimilation models using longer data assimilation cycle lengths; root-mean-square error (RMSE) values of both surface velocity components during a 12 h model forecasting period indicated that surface flow fields were significantly improved when implementing nudging assimilation at each model computational timestep. The Data Assimilation Skill Score (DASS) technique was used to quantitatively evaluate forecast improvements. The averaged values of DASS over the data assimilation domain were 26% and 33% for east–west and north–south velocity components, respectively, over the half-day forecasting period. Correlation of Averaged Kinetic Energy (AKE) was improved by more than 10% in the best data assimilation model. Time series of velocity components and surface flow fields were presented to illustrate the improvement resulting from data assimilation application over time.

Description: Remote Sensing, Vol. 9, Pages 929: Detecting Drought-Induced Tree Mortality in Sierra Nevada Forests with Time Series of Satellite Data Remote Sensing doi: 10.3390/rs9090929 Authors: Sarah Byer Yufang Jin A five-year drought in California led to a significant increase in tree mortality in the Sierra Nevada forests from 2012 to 2016. Landscape level monitoring of forest health and tree dieback is critical for vegetation and disaster management strategies. We examined the capability of multispectral imagery from the Moderate Resolution Imaging Spectroradiometer (MODIS) in detecting and explaining the impacts of the recent severe drought in Sierra Nevada forests. Remote sensing metrics were developed to represent baseline forest health conditions and drought stress using time series of MODIS vegetation indices (VIs) and a water index. We used Random Forest algorithms, trained with forest aerial detection surveys data, to detect tree mortality based on the remote sensing metrics and topographical variables. Map estimates of tree mortality demonstrated that our two-stage Random Forest models were capable of detecting the spatial patterns and severity of tree mortality, with an overall producer’s accuracy of 96.3% for the classification Random Forest (CRF) and a RMSE of 7.19 dead trees per acre for the regression Random Forest (RRF). The overall omission errors of the CRF ranged from 19% for the severe mortality class to 27% for the low mortality class. Interpretations of the models revealed that forests with higher productivity preceding the onset of drought were more vulnerable to drought stress and, consequently, more likely to experience tree mortality. This method highlights the importance of incorporating baseline forest health data and measurements of drought stress in understanding forest response to severe drought.

Description: Remote Sensing, Vol. 9, Pages 930: A Feasibility Study of Sea Ice Motion and Deformation Measurements Using Multi-Sensor High-Resolution Optical Satellite Images Remote Sensing doi: 10.3390/rs9090930 Authors: Chang-Uk Hyun Hyun-cheol Kim Sea ice motion and deformation have generally been measured using low-resolution passive microwave or mid-resolution radar remote sensing datasets of daily (or few days) intervals to monitor long-term trends over a wide polar area. This feasibility study presents an application of high-resolution optical images from operational satellites, which have become more available in polar regions, for sea ice motion and deformation measurements. The sea ice motion, i.e., Lagrangian vector, is measured by using a maximum cross-correlation (MCC) technique and multi-temporal high-resolution images acquired on 14–15 August 2014 from multiple spaceborne sensors on board Korea Multi-Purpose Satellites (KOMPSATs) with short acquisition time intervals. The sea ice motion extracted from the six image pairs of the spatial resolutions were resampled to 4 m and 15 m yields with vector length measurements of 57.7 m root mean square error (RMSE) and −11.4 m bias and 60.7 m RMSE and −13.5 m bias, respectively, compared with buoy location records. The errors from both resolutions indicate more accurate measurements than from conventional sea ice motion datasets from passive microwave and radar data in ice and water mixed surface conditions. In the results of sea ice deformation caused by interaction of individual ice floes, while free drift patterns of ice floes were delineated from the 4 m spatial resolution images, the deformation was less revealing in the 15 m spatial resolution image pairs due to emphasized discretization uncertainty from coarser pixel sizes. The results demonstrate that using multi-temporal high-resolution optical satellite images enabled precise image block matching in the melting season, thus this approach could be used for expanding sea ice motion and deformation dataset, with an advantage of frequent image acquisition capability in multiple areas by means of many operational satellites.

Description: Remote Sensing, Vol. 9, Pages 937: Landslide Displacement Monitoring with Split-Bandwidth Interferometry: A Case Study of the Shuping Landslide in the Three Gorges Area Remote Sensing doi: 10.3390/rs9090937 Authors: Xuguo Shi Houjun Jiang Lu Zhang Mingsheng Liao Landslides constitute a major threat to people’s lives and property in mountainous regions such, as in the Three Gorges area in China. Synthetic Aperture Radar Interferometry (InSAR) with its wide coverage and unprecedented displacement measuring capabilities has been widely used in landslide monitoring. However, it is difficult to apply traditional InSAR techniques to investigate landslides having large deformation gradients or moving primarily in north-south direction. In this study, we propose a time series split-bandwidth interferometry (SBI) procedure to measure two dimensional (azimuth and range) displacements of the Shuping landslide in the Three Gorges area with 36 TerraSAR-X high resolution spotlight (HS) images acquired from February 2009 to April 2010. Since the phase based SBI procedure is sensitive to noise, we focused on extracting displacements of corner reflectors (CRs) installed on or surrounding the Shuping landslide. Our results agreed well with measurements obtained by the point-like targets offset tracking (PTOT) technique and in-situ GPS stations. Centimeter level accuracy could be achieved with SBI on CRs which shows great potential in futures studies on fast moving geohazards.

Description: Remote Sensing, Vol. 9, Pages 938: Integration of Information Theory, K-Means Cluster Analysis and the Logistic Regression Model for Landslide Susceptibility Mapping in the Three Gorges Area, China Remote Sensing doi: 10.3390/rs9090938 Authors: Qian Wang Yi Wang Ruiqing Niu Ling Peng In this work, an effective framework for landslide susceptibility mapping (LSM) is presented by integrating information theory, K-means cluster analysis and statistical models. In general, landslides are triggered by many causative factors at a local scale, and the impact of these factors is closely related to geographic locations and spatial neighborhoods. Based on these facts, the main idea of this research is to group a study area into several clusters to ensure that landslides in each cluster are affected by the same set of selected causative factors. Based on this idea, the proposed predictive method is constructed for accurate LSM at a regional scale by applying a statistical model to each cluster of the study area. Specifically, each causative factor is first classified by the natural breaks method with the optimal number of classes, which is determined by adopting Shannon’s entropy index. Then, a certainty factor (CF) for each class of factors is estimated. The selection of the causative factors for each cluster is determined based on the CF values of each factor. Furthermore, the logistic regression model is used as an example of statistical models in each cluster using the selected causative factors for landslide prediction. Finally, a global landslide susceptibility map is obtained by combining the regional maps. Experimental results based on both qualitative and quantitative analysis indicated that the proposed framework can achieve more accurate landslide susceptibility maps when compared to some existing methods, e.g., the proposed framework can achieve an overall prediction accuracy of 91.76%, which is 7.63–11.5% higher than those existing methods. Therefore, the local scale LSM technique is very promising for further improvement of landslide prediction.

Description: Remote Sensing, Vol. 9, Pages 944: An Alternative Approach to Using LiDAR Remote Sensing Data to Predict Stem Diameter Distributions across a Temperate Forest Landscape Remote Sensing doi: 10.3390/rs9090944 Authors: Rebecca Spriggs David Coomes Trevor Jones John Caspersen Mark Vanderwel We apply a spatially-implicit, allometry-based modelling approach to predict stem diameter distributions (SDDs) from low density airborne LiDAR data in a heterogeneous, temperate forest in Ontario, Canada. Using a recently published algorithm that relates the density, size, and species of individual trees to the height distribution of first returns, we estimated parameters that succinctly describe SDDs that are most consistent with each 0.25-ha LiDAR tile across a 30,000 ha forest landscape. Tests with independent validation plots showed that the diameter distribution of stems was predicted with reasonable accuracy in most cases (half of validation plots had R2 ≥ 0.75, and another 23% had 0.5 ≤ R2 < 0.75). The predicted frequency of larger stems was much better than that of small stems (8 ≤ x < 11 cm diameter), particularly small conifers. We used the predicted SDDs to calculate aboveground carbon density (ACD; RMSE = 21.4 Mg C/ha), quadratic mean diameter (RMSE = 3.64 cm), basal area (RMSE = 6.99 m2/ha) and stem number (RMSE = 272 stems/ha). The accuracy of our predictions compared favorably with previous studies that have generally been undertaken in simpler conifer-dominated forest types. We demonstrate the utility of our results to spatial forest management planning by mapping SDDs, the proportion of broadleaves, and ACD at a 0.25 ha resolution.

Description: Remote Sensing, Vol. 9, Pages 946: Automatic Detection and Parameter Estimation of Trees for Forest Inventory Applications Using 3D Terrestrial LiDAR Remote Sensing doi: 10.3390/rs9090946 Authors: Ahmad Aijazi Paul Checchin Laurent Malaterre Laurent Trassoudaine Forest inventory plays an important role in the management and planning of forests. In this study, we present a method for automatic detection and estimation of trees, especially in forest environments using 3D terrestrial LiDAR data. The proposed method does not rely on any predefined tree shape or model. It uses the vertical distribution of the 3D points partitioned in a gridded Digital Elevation Model (DEM) to extract out ground points. The cells of the DEM are then clustered together to form super-clusters representing potential tree objects. The 3D points contained in each of these super-clusters are then classified into trunk and vegetation classes using a super-voxel based segmentation method. Different attributes (such as diameter at breast height, basal area, height and volume) are then estimated at individual tree levels which are then aggregated to generate metrics for forest inventory applications. The method is validated and evaluated on three different data sets obtained from three different types of terrestrial sensors (vehicle-borne, handheld and static) to demonstrate its applicability and feasibility for a wide range of applications. The results are evaluated by comparing the estimated parameters with real field observations/measurements to demonstrate the efficacy of the proposed method. Overall segmentation and classification accuracies greater than 84 % while average parameter estimation error ranging from 1 . 6 to 9 % were observed.

Description: Remote Sensing, Vol. 9, Pages 949: Spatial Disaggregation of Latent Heat Flux Using Contextual Models over India Remote Sensing doi: 10.3390/rs9090949 Authors: Rajasekaran Eswar Muddu Sekhar Bimal Bhattacharya Soumya Bandyopadhyay Estimation of latent heat flux at the agricultural field scale is required for proper water management. The current generation thermal sensors except Landsat-8 provide data on the order of 1000 m. The aim of this study is to test three approaches based on contextual models using only remote sensing datasets for the disaggregation of latent heat flux over India. The first two approaches are, respectively, based on the estimation of the evaporative fraction (EF) and solar radiation ratio at coarser resolution and disaggregating them to yield the latent heat flux at a finer resolution. The third approach is based on disaggregation of the thermal data and estimating a finer resolution latent heat flux. The three approaches were tested using MODIS datasets and the validation was done using the Bowen Ratio energy balance observations at five sites across India. From the validation, it was observed that the first two approaches performed similarly and better than the third approach at all five sites. The third approach, based on the disaggregation of the thermal data, yielded larger errors. In addition to better performance, the second approach based on the disaggregation of solar radiation ratio was simpler and required lesser data processing than the other approaches. In addition, the first two approaches captured the spatial pattern of latent heat flux without introducing any artefacts in the final output.

Description: Remote Sensing, Vol. 9, Pages 951: Estimating Rice Leaf Nitrogen Concentration: Influence of Regression Algorithms Based on Passive and Active Leaf Reflectance Remote Sensing doi: 10.3390/rs9090951 Authors: Jia Sun Jian Yang Shuo Shi Biwu Chen Lin Du Wei Gong Shalei Song Nitrogen (N) is important for the growth of crops. Estimating leaf nitrogen concentration (LNC) accurately and nondestructively is important for precision agriculture, reduces environmental pollution, and helps model global carbon and N cycles. Leaf reflectance, especially in the visible and near-infrared regions, has been identified as a useful indicator of LNC. Except reflectance passively acquired by spectrometers, the newly developed multispectral LiDAR and hyperspectral LiDAR provide possibilities for measuring leaf spectra actively. The regression relationship between leaf reflectance spectra and rice (Oryza sativa) LNC relies greatly on the algorithm adopted. It would be preferable to find one algorithm that performs well with respect to passive and active leaf spectra. Thus, this study assesses the influence of six popular linear and nonlinear methods on rice LNC retrieval, namely, partial least-square regression, least squares boosting, bagging, random forest, back-propagation neural network (BPNN), and support vector regression of different types/kernels/parameter values. The R2, root mean square error and relative error in rice LNC estimation using these different methods were compared through the passive and active spectral measurements of rice leaves of different varieties at different locations and time (Yongyou 4949, Suizhou, 2014, Yangliangyou 6, Wuhan, 2015). Results demonstrate that BPNN provided generally satisfactory performance in estimating rice LNC using the three kinds of passive and active reflectance spectra.

Description: Remote Sensing, Vol. 9, Pages 948: Inter-System Differencing between GPS and BDS for Medium-Baseline RTK Positioning Remote Sensing doi: 10.3390/rs9090948 Authors: Wang Gao Chengfa Gao Shuguo Pan Xiaolin Meng Yan Xia An inter-system differencing model between two Global Navigation Satellite Systems (GNSS) enables only one reference satellite for all observations. If the associated differential inter-system biases (DISBs) are priori known, double-differenced (DD) ambiguities between overlapping frequencies from different GNSS constellations can also be fixed to integers. This can provide more redundancies for the observation model, and thus will be beneficial to ambiguity resolution (AR) and real-time kinematic (RTK) positioning. However, for Global Positioning System (GPS) and the regional BeiDou Navigation Satellite System (BDS-2), there are no overlapping frequencies. Tight combination of GPS and BDS needs to process not only the DISBs but also the single-difference ambiguity of the reference satellite, which is caused by the influence of different frequencies. In this paper, we propose a tightly combined dual-frequency GPS and BDS RTK positioning model for medium baselines with real-time estimation of DISBs. The stability of the pseudorange and phase DISBs is analyzed firstly using several baselines with the same or different receiver types. The dual-frequency ionosphere-free model with parameterization of GPS-BDS DISBs is proposed, where the single-difference ambiguity is estimated jointly with the phase DISB parameter from epoch to epoch. The performance of combined GPS and BDS RTK positioning for medium baselines is evaluated with simulated obstructed environments. Experimental results show that with the inter-system differencing model, the accuracy and reliability of RTK positioning can be effectively improved, especially for the obstructed environments with a small number of satellites available.

Description: Remote Sensing, Vol. 9, Pages 960: An Enhanced IT2FCM* Algorithm Integrating Spectral Indices and Spatial Information for Multi-Spectral Remote Sensing Image Clustering Remote Sensing doi: 10.3390/rs9090960 Authors: Jifa Guo Hongyuan Huo Interval type-2 fuzzy c-means (IT2FCM) clustering methods for remote-sensing data classification are based on interval type-2 fuzzy sets and can effectively handle uncertainty of membership grade. However, most of these methods neglect the spatial information when they are used in image clustering. The spatial information and spectral indices are useful in remote-sensing data classification. Thus, determining how to integrate them into IT2FCM to improve the quality and accuracy of the classification is very important. This paper proposes an enhanced IT2FCM* (EnIT2FCM*) algorithm by combining spatial information and spectral indices for remote-sensing data classification. First, the new comprehensive spatial information is defined as the combination of the local spatial distance and attribute distance or membership-grade distance. Then, a novel distance metric is proposed by combining this new spatial information and the selected spectral indices; these selected spectral indices are treated as another dataset in this distance metric. To test the effectiveness of the EnIT2FCM* algorithm, four typical validity indices along with the confusion matrix and kappa coefficient are used. The experimental results show that the spatial information definition proposed here is effective, and some spectral indices and their combinations improve the performance of the EnIT2FCM*. Thus, the selection of suitable spectral indices is crucial, and the combination of soil adjusted vegetation index (SAVI) and the Automated Water Extraction Index (AWEIsh) is the best choice of spectral indices for this method.

Description: Remote Sensing, Vol. 9, Pages 959: Estimating High Resolution Daily Air Temperature Based on Remote Sensing Products and Climate Reanalysis Datasets over Glacierized Basins: A Case Study in the Langtang Valley, Nepal Remote Sensing doi: 10.3390/rs9090959 Authors: Wang Zhou Bin Peng Jiancheng Shi Tianxing Wang Yam Dhital Ruzhen Yao Yuechi Yu Zhongteng Lei Rui Zhao Near surface air temperature (Ta) is one of the key input parameters in land surface models and hydrological models as it affects most biogeophysical and biogeochemical processes of the earth surface system. For distributed hydrological modeling over glacierized basins, obtaining high resolution Ta forcing is one of the major challenges. In this study, we proposed a new high resolution daily Ta estimation scheme under both clear and cloudy sky conditions through integrating the moderate-resolution imaging spectroradiometer (MODIS) land surface temperature (LST) and China Meteorological Administration (CMA) land data assimilation system (CLDAS) reanalyzed daily Ta. Spatio-temporal continuous MODIS LST was reconstructed through the data interpolating empirical orthogonal functions (DINEOF) method. Multi-variable regression models were developed at CLDAS scale and then used to estimate Ta at MODIS scale. The new Ta estimation scheme was tested over the Langtang Valley, Nepal as a demonstrating case study. Observations from two automatic weather stations at Kyanging and Yala located in the Langtang Valley from 2012 to 2014 were used to validate the accuracy of Ta estimation. The RMSEs are 2.05, 1.88, and 3.63 K, and the biases are 0.42, −0.68 and −2.86 K for daily maximum, mean and minimum Ta, respectively, at the Kyanging station. At the Yala station, the RMSE values are 4.53, 2.68 and 2.36 K, and biases are 4.03, 1.96 and −0.35 K for the estimated daily maximum, mean and minimum Ta, respectively. Moreover, the proposed scheme can produce reasonable spatial distribution pattern of Ta at the Langtang Valley. Our results show the proposed Ta estimation scheme is promising for integration with distributed hydrological model for glacier melting simulation over glacierized basins.

Description: Remote Sensing, Vol. 9, Pages 829: An Improved Vegetation Adjusted Nighttime Light Urban Index and Its Application in Quantifying Spatiotemporal Dynamics of Carbon Emissions in China Remote Sensing doi: 10.3390/rs9080829 Authors: Xing Meng Ji Han Cheng Huang Nighttime Light (NTL) has been widely used as a proxy of many socio-environmental issues. However, the limited range of sensor radiance of NTL prevents its further application and estimation accuracy. To improve the performance, we developed an improved Vegetation Adjusted Nighttime light Urban Index (VANUI) through fusing multi-year NTL with population density, the Normalized Difference Vegetation Index and water body data and applied it to fine-scaled carbon emission analysis in China. The results proved that our proposed index could reflect more spatial variation of human activities. It is also prominent in reducing the carbon modeling error at the inter-city level and distinguishing the emission heterogeneity at the intra-city level. Between 1995 and 2013, CO2 emissions increased significantly in China, but were distributed unevenly in space with high density emissions mainly located in metropolitan areas and provincial capitals. In addition to Beijing-Tianjin-Hebei, the Yangzi River Delta and the Pearl River Delta, the Shandong Peninsula has become a new emission hotspot that needs special attention in carbon mitigation. The improved VANUI and its application to the carbon emission issue not only broadened our understanding of the spatiotemporal dynamics of fine-scaled CO2 emission, but also provided implications for low-carbon and sustainable development plans.

Description: Remote Sensing, Vol. 9, Pages 825: Semi-Analytical Retrieval of the Diffuse Attenuation Coefficient in Large and Shallow Lakes from GOCI, a High Temporal-Resolution Satellite Remote Sensing doi: 10.3390/rs9080825 Authors: Changchun Huang Ling Yao Monitoring the dynamic characteristics of the diffuse attenuation coefficient (Kd(490)) on the basis of the high temporal-resolution satellite data is critical for regulating the ecological environment of lake. By measuring the in-situ Kd(490) and the remote-sensing reflectance, a semi-analytical algorithm for Kd(490) was developed to determine the short-term variation of Kd(490). From 2006 to 2014, the data about 412 samples (among which 60 were used as match-up points, 282 for calibrating dataset and the remaining 70 for validating dataset) were gathered from nine expeditions to calibrate and validate the aforesaid semi-analytical algorithm. The root mean square percentage error (RMSP) and the mean absolute relative error (MAPE) of validation datasets were respectively 27.44% and 22.60 ± 15.57%, while that of the match-up datasets were respectively 34.29% and 27.57 ± 20.56%. These percentages indicate that the semi-analytical algorithm and Geostationary Ocean Color Imager (GOCI) data are applicable to obtain the short-term variation of Kd(490) in the turbid shallow inland waters. The short-term GOCI-observed Kd(490) shows a significant seasonal and spatial variation and a similar distribution to the matching Moderate Resolution Imaging Spectroradiometer (MODIS) which derived Kd(490). A comparative analysis on wind (observed by buoys) and GOCI-derived Kd(490) suggests that wind is a primary driving factor of Kd(490) variation, but the lacustrine morphometry affects the wind force that is contributing to Kd(490) variation.

Description: Remote Sensing, Vol. 9, Pages 828: Adaptive Estimation of Crop Water Stress in Nectarine and Peach Orchards Using High-Resolution Imagery from an Unmanned Aerial Vehicle (UAV) Remote Sensing doi: 10.3390/rs9080828 Authors: Suyoung Park Dongryeol Ryu Sigfredo Fuentes Hoam Chung Esther Hernández-Montes Mark O’Connell The capability to monitor water status from crops on a regular basis can enhance productivity and water use efficiency. In this paper, high-resolution thermal imagery acquired by an unmanned aerial vehicle (UAV) was used to map plant water stress and its spatial variability, including sectors under full irrigation and deficit irrigation over nectarine and peach orchards at 6.12 cm ground sample distance. The study site was classified into sub-regions based on crop properties, such as cultivars and tree training systems. In order to enhance the accuracy of the mapping, edge extraction and filtering were conducted prior to the probability modelling employed to obtain crop-property-specific (‘adaptive’ hereafter) lower and higher temperature references (Twet and Tdry respectively). Direct measurements of stem water potential (SWP, ψstem) and stomatal conductance (gs) were collected concurrently with UAV remote sensing and used to validate the thermal index as crop biophysical parameters. The adaptive crop water stress index (CWSI) presented a better agreement with both ψstem and gs with determination coefficients (R2) of 0.72 and 0.82, respectively, while the conventional CWSI applied by a single set of hot and cold references resulted in biased estimates with R2 of 0.27 and 0.34, respectively. Using a small number of ground-based measurements of SWP, CWSI was converted to a high-resolution SWP map to visualize spatial distribution of the water status at field scale. The results have important implications for the optimal management of irrigation for crops.

Description: Remote Sensing, Vol. 9, Pages 826: LiDAR Validation of a Video-Derived Beachface Topography on a Tidal Flat Remote Sensing doi: 10.3390/rs9080826 Authors: David Didier Pascal Bernatchez Emmanuel Augereau Charles Caulet Dany Dumont Eliott Bismuth Louis Cormier France Floc’h Christophe Delacourt Increasingly used shore-based video stations enable a high spatiotemporal frequency analysis of shoreline migration. Shoreline detection techniques combined with hydrodynamic conditions enable the creation of digital elevation models (DEMs). However, shoreline elevations are often estimated based on nearshore process empirical equations leading to uncertainties in video-based topography. To achieve high DEM correspondence between both techniques, we assessed video-derived DEMs against LiDAR surveys during low energy conditions. A newly installed video system on a tidal flat in the St. Lawrence Estuary, Atlantic Canada, served as a test case. Shorelines were automatically detected from time-averaged (TIMEX) images using color ratios in low energy conditions synchronously with mobile terrestrial LiDAR during two different surveys. Hydrodynamic (waves and tides) data were recorded in-situ, and established two different cases of water elevation models as a basis for shoreline elevations. DEMs were created and tested against LiDAR. Statistical analysis of shoreline elevations and migrations were made, and morphological variability was assessed between both surveys. Results indicate that the best shoreline elevation model includes both the significant wave height and the mean water level. Low energy conditions and in-situ hydrodynamic measurements made it possible to produce video-derived DEMs virtually as accurate as a LiDAR product, and therefore make an effective tool for coastal managers.

Description: Remote Sensing, Vol. 9, Pages 841: A Probabilistic Weighted Archetypal Analysis Method with Earth Mover’s Distance for Endmember Extraction from Hyperspectral Imagery Remote Sensing doi: 10.3390/rs9080841 Authors: Weiwei Sun Dianfa Zhang Yan Xu Long Tian Gang Yang Weiyue Li A Probabilistic Weighted Archetypal Analysis method with Earth Mover’s Distance (PWAA-EMD) is proposed to extract endmembers from hyperspectral imagery (HSI). The PWAA-EMD first utilizes the EMD dissimilarity matrix to weight the coefficient matrix in the regular Archetypal Analysis (AA). The EMD metric considers manifold structures of spectral signatures in the HSI data and could better quantify the dissimilarity features among pairwise pixels. Second, the PWAA-EMD adopts the Bayesian framework and formulates the improved AA into a probabilistic inference problem by maximizing a joint posterior density. Third, the optimization problem is solved by the iterative multiplicative update scheme, with a careful initialization from the two-stage algorithm and the proper endmembers are finally obtained. The synthetic and real Cuprite Hyperspectral datasets are utilized to verify the performance of PWAA-EMD and five popular methods are implemented to make comparisons. The results show that PWAA-EMD surpasses all the five methods in the average results of spectral angle distance (SAD) and root-mean-square-error (RMSE). Especially, the PWAA-EMD obtains more accurate estimation than AA in almost all the classes of endmembers including two similar ones. Therefore, the PWAA-EMD could be an alternative choice for endmember extraction on the hyperspectral data.

Description: Remote Sensing, Vol. 9, Pages 842: Long-Term Water Storage Changes of Lake Volta from GRACE and Satellite Altimetry and Connections with Regional Climate Remote Sensing doi: 10.3390/rs9080842 Authors: Shengnan Ni Jianli Chen Clark Wilson Xiaogong Hu Satellite gravity data from the Gravity Recovery and Climate Experiment (GRACE) provides a quantitative measure of terrestrial water storage (TWS) change at different temporal and spatial scales. In this study, we investigate the ability of GRACE to quantitatively monitor long-term hydrological characteristics over the Lake Volta region. Principal component analysis (PCA) is employed to study temporal and spatial variability of long-term TWS changes. Long-term Lake Volta water storage change appears to be the dominant long-term TWS change signal in the Volta basin. GRACE-derived TWS changes and precipitation variations compiled by the Global Precipitation Climatology Centre (GPCC) are related both temporally and spatially, but spatial leakage attenuates the magnitude of GRACE estimates, especially at small regional scales. Using constrained forward modeling, we successfully remove leakage error in GRACE estimates. After this leakage correction, GRACE-derived Lake Volta water storage changes agree remarkably well with independent estimates from satellite altimetry at interannual and longer time scales. This demonstrates the value of GRACE estimates to monitor and quantify water storage changes in lakes, especially in relatively small regions with complicated topography.

Description: Remote Sensing, Vol. 9, Pages 839: Multiscale Remote Sensing to Map the Spatial Distribution and Extent of Cropland in the Sudanian Savanna of West Africa Remote Sensing doi: 10.3390/rs9080839 Authors: Gerald Forkuor Christopher Conrad Michael Thiel Benewinde Zoungrana Jérôme Tondoh Food security is the topmost priority on the global agenda. Accurate agricultural statistics (i.e., cropped area) are essential for decision making and developing appropriate programs to achieve food security. However, derivation of these essential agricultural statistics, especially in developing countries, is fraught with many challenges including financial, logistical and human capacity limitations. This study investigated the use of fractional cover approaches in mapping cropland area in the heterogeneous landscape of West Africa. Discrete cropland areas identified from multi-temporal Landsat data were upscaled to MODIS resolution using random forest regression. Producer’s accuracy and user’s accuracy of the cropland class in the Landsat scale analysis averaged 95% and 94%, respectively, indicating good separability between crop and non-crop land. Validation of the fractional cropland cover map at MODIS resolution (MODIS_FCM) revealed an overall mean absolute error of 19%. Comparison of MODIS_FCM with the MODIS land cover product (e.g., MODIS_LCP) demonstrate the suitability of the proposed approach to cropped area estimation in smallholder dominant heterogeneous landscapes over existing global solutions. Comparison with official government statistics (i.e., cropped area) revealed variable levels of agreement and partly enormous disagreements, which clearly indicate the need to integrate remote sensing approaches and ground based surveys conducted by agricultural ministries in improving cropped area estimation. The recent availability of a wide range of open access remote sensing data is expected to expedite this integration and contribute missing information urgently required for regional assessments of food security in West Africa and beyond.

Description: Remote Sensing, Vol. 9, Pages 844: Extension of a Fast GLRT Algorithm to 5D SAR Tomography of Urban Areas Remote Sensing doi: 10.3390/rs9080844 Authors: Alessandra Budillon Angel Caroline Johnsy Gilda Schirinzi This paper analyzes a method for Synthetic Aperture Radar (SAR) Tomographic (TomoSAR) imaging, allowing the detection of multiple scatterers that can exhibit time deformation and thermal dilation by using a CFAR (Constant False Alarm Rate) approach. In the last decade, several methods for TomoSAR have been proposed. The objective of this paper is to present the results obtained on high resolution tomographic SAR data of urban areas, by using a statistical test for detecting multiple scatterers that takes into account phase variations due to possible deformations and/or thermal dilation. The test can be evaluated in terms of probability of detection (PD) and probability of false alarm (PFA), and is based on an approximation of a Generalized Likelihood Ratio Test (GLRT), denoted as Fast-Sup-GLRT. It was already applied and validated by the authors in the 3D case, while here it is extended and experimented in the 5D case. Numerical experiments on simulated and on StripMap TerraSAR-X (TSX) data have been carried out. The presented results show that the adopted method allows the detection of a large number of scatterers and the estimation of their position with a good accuracy, and that the consideration of the thermal dilation and surface deformation helps in recovering more single and double scatterers, with respect to the case in which these contributions are not taken into account. Moreover, the capability of method to provide reliable estimates of the deformations in urban structure suggests its use in structure stress monitoring.

Description: Remote Sensing, Vol. 9, Pages 848: Pre-Trained AlexNet Architecture with Pyramid Pooling and Supervision for High Spatial Resolution Remote Sensing Image Scene Classification Remote Sensing doi: 10.3390/rs9080848 Authors: Xiaobing Han Yanfei Zhong Liqin Cao Liangpei Zhang The rapid development of high spatial resolution (HSR) remote sensing imagery techniques not only provide a considerable amount of datasets for scene classification tasks but also request an appropriate scene classification choice when facing with finite labeled samples. AlexNet, as a relatively simple convolutional neural network (CNN) architecture, has obtained great success in scene classification tasks and has been proven to be an excellent foundational hierarchical and automatic scene classification technique. However, current HSR remote sensing imagery scene classification datasets always have the characteristics of small quantities and simple categories, where the limited annotated labeling samples easily cause non-convergence. For HSR remote sensing imagery, multi-scale information of the same scenes can represent the scene semantics to a certain extent but lacks an efficient fusion expression manner. Meanwhile, the current pre-trained AlexNet architecture lacks a kind of appropriate supervision for enhancing the performance of this model, which easily causes overfitting. In this paper, an improved pre-trained AlexNet architecture named pre-trained AlexNet-SPP-SS has been proposed, which incorporates the scale pooling—spatial pyramid pooling (SPP) and side supervision (SS) to improve the above two situations. Extensive experimental results conducted on the UC Merced dataset and the Google Image dataset of SIRI-WHU have demonstrated that the proposed pre-trained AlexNet-SPP-SS model is superior to the original AlexNet architecture as well as the traditional scene classification methods.

Description: Remote Sensing, Vol. 9, Pages 849: Correction: Yao, P. et al. Rebuilding Long Time Series Global Soil Moisture Products Using the Neural Network Adopted the Microwave Vegetation Index. Remote Sens. 2017, 9, 35 Remote Sensing doi: 10.3390/rs9080849 Authors: Panpan Yao Jiancheng Shi Tianjie Zhao Hui Lu Amen Al-Yaari After publication of the research paper [1], the authors wish to make the following correction to this paper. In the fourth line from the bottom in abstract, due to a typing error, “RMSE = 0.84 m3/m3” should be replaced with “RMSE = 0.084 m3/m3”.[...]

Description: Remote Sensing, Vol. 9, Pages 846: A Novel Method of Change Detection in Bi-Temporal PolSAR Data Using a Joint-Classification Classifier Based on a Similarity Measure Remote Sensing doi: 10.3390/rs9080846 Authors: Jinqi Zhao Jie Yang Zhong Lu Pingxiang Li Wensong Liu Le Yang Accurate and timely change detection of the Earth’s surface features is extremely important for understanding the relationships and interactions between people and natural phenomena. Owing to the all-weather response capability, polarimetric synthetic aperture radar (PolSAR) has become a key tool for change detection. Change detection includes both unsupervised and supervised methods. Unsupervised change detection is simple and effective, but cannot detect the type of land cover change. Supervised change detection can detect the type of land cover change, but is easily affected and depended by the human interventions. To solve these problems, a novel method of change detection using a joint-classification classifier (JCC) based on a similarity measure is introduced. The similarity measure is obtained by a test statistic and the Kittler and Illingworth (TSKI) minimum-error thresholding algorithm, which is used to automatically control the JCC. The efficiency of the proposed method is demonstrated by the use of bi-temporal PolSAR images acquired by RADARSAT-2 over Wuhan, China. The experimental results show that the proposed method can identify the different types of land cover change and can reduce both the false detection rate and false alarm rate in the change detection.

Description: Remote Sensing, Vol. 9, Pages 852: Influence of Droughts on Mid-Tropospheric CO2 Remote Sensing doi: 10.3390/rs9080852 Authors: Xun Jiang Angela Kao Abigail Corbett Edward Olsen Thomas Pagano Albert Zhai Sally Newman Liming Li Yuk Yung Using CO2 data from the Atmospheric Infrared Sounder (AIRS), it is found for the first time that the mid-tropospheric CO2 concentration is ~1 part per million by volume higher during dry years than wet years over the southwestern USA from June to September. The mid-tropospheric CO2 differences between dry and wet years are related to circulation and CO2 surface fluxes. During drought conditions, vertical pressure velocity from NCEP2 suggests that there is more rising air over most regions, which can help bring high surface concentrations of CO2 to the mid-troposphere. In addition to the circulation, there is more CO2 emitted from the biosphere to the atmosphere during droughts in some regions, which can contribute to higher concentrations of CO2 in the atmosphere. Results obtained from this study demonstrate the significant impact of droughts on atmospheric CO2 and therefore on a feedback cycle contributing to greenhouse gas warming. It can also help us better understand atmospheric CO2, which plays a critical role in our climate system.

Description: Remote Sensing, Vol. 9, Pages 853: Specular Reflection Effects Elimination in Terrestrial Laser Scanning Intensity Data Using Phong Model Remote Sensing doi: 10.3390/rs9080853 Authors: Kai Tan Xiaojun Cheng The intensity value recorded by terrestrial laser scanning (TLS) systems is significantly influenced by the incidence angle. The incidence angle effect is an object property, which is mainly related to target scattering properties, surface structures, and even some instrumental effects. Most existing models focus on diffuse reflections of rough surfaces and ignore specular reflections, despite that both reflections simultaneously exist in all natural surfaces. Due to the coincidence of the emitter and receiver in TLS, specular reflections can be ignored at large incidence angles. On the contrary, at small incidence angles, TLS detectors can receive a portion of specular reflections. The received specular reflections can trigger highlight phenomenon (hot-spot effects) in the intensity data of the scanned targets, particularly those with a relatively smooth or highly-reflective surface. In this study, a new method that takes diffuse and specular reflections, as well as the instrumental effects into consideration, is proposed to eliminate the specular reflection effects in TLS intensity data. Diffuse reflections and instrumental effects are modeled by a polynomial based on Lambertian reference targets, whereas specular reflections are modeled by the Phong model. The proposed method is tested and validated on different targets scanned by the Faro Focus3D 120 terrestrial scanner. Results imply that the coefficient of variation of the intensity data from a homogeneous surface is reduced by approximately 38% when specular reflections are considered. Compared with existing methods, the proposed method exhibits good feasibility and high accuracy in eliminating the specular reflection effects for intensity image interpretation and 3D point cloud representation by intensity.

Description: IJGI, Vol. 6, Pages 256: Spatiotemporal Assessment of Littoral Waterbirds for Establishing Ecological Indicators of Mediterranean Coastal Lagoons ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6080256 Authors: Pablo Farinós-Celdrán Francisco Robledano-Aymerich María Francisca Carreño Javier Martínez-López Waterbirds are vital indicators of anthropogenic influence on the ecological status of Mediterranean coastal lagoons. Our study relates temporal waterbird data to key environmental gradients at catchment scale that have a structural or functional influence on littoral waterbird assemblages at different scales. During two full-year cycles and two additional wintering seasons, the nearshore waterbird assemblages of the Mar Menor coastal lagoon (Murcia Region, SE Spain) were monitored monthly. Several biological indicator variables were related to the anthropogenic environmental gradient in the catchment area. Results showed that there was a strong dependence of waterbird assemblages on the distance to shore, emphasizing the importance of the first 100-m band, in which many species relevant to conservation converge on food resources. Well-preserved shoreline tracts therefore had a clear positive effect on community richness and diversity values, and were correlated with the occurrence of some species. These results clearly support the need for effective protection and restoration measures of such littoral habitats. Specific responses to local disturbing processes were nested within habitat and landscape preferences, supporting the value of aquatic birds as integrative ecological signals in semi-enclosed coastal systems. Moreover, waterbird-based indicators responded positively to environmental improvements both qualitatively and quantitatively.

Description: Remote Sensing, Vol. 9, Pages 857: A Robust Algorithm for Estimating Surface Fractional Vegetation Cover from Landsat Data Remote Sensing doi: 10.3390/rs9080857 Authors: Linqing Yang Kun Jia Shunlin Liang Xiangqin Wei Yunjun Yao Xiaotong Zhang Fractional vegetation cover (FVC) is an essential land surface parameter for Earth surface process simulations and global change studies. The currently existing FVC products are mostly obtained from low or medium resolution remotely sensed data, while many applications require the fine spatial resolution FVC product. The availability of well-calibrated coverage of Landsat imagery over large areas offers an opportunity for the production of FVC at fine spatial resolution. Therefore, the objective of this study is to develop a general and reliable land surface FVC estimation algorithm for Landsat surface reflectance data under various land surface conditions. Two machine learning methods multivariate adaptive regression splines (MARS) model and back-propagation neural networks (BPNNs) were trained using samples from PROSPECT leaf optical properties model and the scattering by arbitrarily inclined leaves (SAIL) model simulations, which included Landsat reflectance and corresponding FVC values, and evaluated to choose the method which had better performance. Thereafter, the MARS model, which had better performance in the independent validation, was evaluated using ground FVC measurements from two case study areas. The direct validation of the FVC estimated using the proposed algorithm (Heihe: R2 = 0.8825, RMSE = 0.097; Chengde using Landsat 7 ETM+: R2 = 0.8571, RMSE = 0.078) (Chengde using Landsat 8 OLI: R2 = 0.8598, RMSE = 0.078) showed the proposed method had good performance. Spatial-temporal assessment of the estimated FVC from Landsat 7 ETM+ and Landsat 8 OLI data confirmed the robustness and consistency of the proposed method. All these results indicated that the proposed algorithm could obtain satisfactory accuracy and had the potential for the production of high-quality FVC estimates from Landsat surface reflectance data.

Description: Remote Sensing, Vol. 9, Pages 858: Deriving Hourly PM2.5 Concentrations from Himawari-8 AODs over Beijing–Tianjin–Hebei in China Remote Sensing doi: 10.3390/rs9080858 Authors: Wei Wang Feiyue Mao Lin Du Zengxin Pan Wei Gong Shenghui Fang Monitoring fine particulate matter with diameters of less than 2.5 μm (PM2.5) is a critical endeavor in the Beijing–Tianjin–Hebei (BTH) region, which is one of the most polluted areas in China. Polar orbit satellites are limited by observation frequency, which is insufficient for understanding PM2.5 evolution. As a geostationary satellite, Himawari-8 can obtain hourly optical depths (AODs) and overcome the estimated PM2.5 concentrations with low time resolution. In this study, the evaluation of Himawari-8 AODs by comparing with Aerosol Robotic Network (AERONET) measurements showed Himawari-8 retrievals (Level 3) with a mild underestimate of about −0.06 and approximately 57% of AODs falling within the expected error established by the Moderate-resolution Imaging Spectroradiometer (MODIS) (±(0.05 + 0.15AOD)). Furthermore, the improved linear mixed-effect model was proposed to derive the surface hourly PM2.5 from Himawari-8 AODs from July 2015 to March 2017. The estimated hourly PM2.5 concentrations agreed well with the surface PM2.5 measurements with high R2 (0.86) and low RMSE (24.5 μg/m3). The average estimated PM2.5 in the BTH region during the study time range was about 55 μg/m3. The estimated hourly PM2.5 concentrations ranged extensively from 35.2 ± 26.9 μg/m3 (1600 local time) to 65.5 ± 54.6 μg/m3 (1100 local time) at different hours.

Description: Remote Sensing, Vol. 9, Pages 869: Detection of Asian Dust Storm Using MODIS Measurements Remote Sensing doi: 10.3390/rs9080869 Authors: Yong Xie Wenhao Zhang John Qu Every year, a large number of aerosols are released from dust storms into the atmosphere, which may have potential impacts on the climate, environment, and air quality. Detecting dust aerosols and monitoring their movements and evolutions in a timely manner is a very significant task. Satellite remote sensing has been demonstrated as an effective means for observing dust aerosols. In this paper, an algorithm based on the multi-spectral technique for detecting dust aerosols was developed by combining measurements of moderate resolution imaging spectroradiometer (MODIS) reflective solar bands and thermal emissive bands. Data from dust events that occurred during the past several years were collected as training data for spectral and statistical analyses. According to the spectral curves of various scene types, a series of spectral bands was selected individually or jointly, and corresponding thresholds were defined for step-by-step scene classification. The multi-spectral algorithm was applied mainly to detect dust storms in Asia. The detection results were validated not only visually with MODIS true color images, but also quantitatively with products of Ozone Monitoring Instrument (OMI) and Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP). The validations showed that this multi-spectral detection algorithm was suitable to monitor dust aerosols in the selected study areas.

Description: Remote Sensing, Vol. 9, Pages 866: Azimuth Ambiguities Removal in Littoral Zones Based on Multi-Temporal SAR Images Remote Sensing doi: 10.3390/rs9080866 Authors: Xiangguang Leng Kefeng Ji Shilin Zhou Huanxin Zou Synthetic aperture radar (SAR) is one of the most important techniques for ocean monitoring. Azimuth ambiguities are a real problem in SAR images today, which can cause performance degradation in SAR ocean applications. In particular, littoral zones can be strongly affected by land-based sources, whereas they are usually regions of interest (ROI). Given the presence of complexity and diversity in littoral zones, azimuth ambiguities removal is a tough problem. As SAR sensors can have a repeat cycle, multi-temporal SAR images provide new insight into this problem. A method for azimuth ambiguities removal in littoral zones based on multi-temporal SAR images is proposed in this paper. The proposed processing chain includes co-registration, local correlation, binarization, masking, and restoration steps. It is designed to remove azimuth ambiguities caused by fixed land-based sources. The idea underlying the proposed method is that sea surface is dynamic, whereas azimuth ambiguities caused by land-based sources are constant. Thus, the temporal consistence of azimuth ambiguities is higher than sea clutter. It opens up the possibilities to use multi-temporal SAR data to remove azimuth ambiguities. The design of the method and the experimental procedure are based on images from the Sentinel data hub of Europe Space Agency (ESA). Both Interferometric Wide Swath (IW) and Stripmap (SM) mode images are taken into account to validate the proposed method. This paper also presents two RGB composition methods for better azimuth ambiguities visualization. Experimental results show that the proposed method can remove azimuth ambiguities in littoral zones effectively.

Description: Remote Sensing, Vol. 9, Pages 876: Mapping Vulnerable Urban Areas Affected by Slow-Moving Landslides Using Sentinel-1 InSAR Data Remote Sensing doi: 10.3390/rs9090876 Authors: Marta Béjar-Pizarro Davide Notti Rosa M. Mateos Pablo Ezquerro Giuseppe Centolanza Gerardo Herrera Guadalupe Bru Margarita Sanabria Lorenzo Solari Javier Duro José Fernández Landslides are widespread natural hazards that generate considerable damage and economic losses worldwide. Detecting terrain movements caused by these phenomena and characterizing affected urban areas is critical to reduce their impact. Here we present a fast and simple methodology to create maps of vulnerable buildings affected by slow-moving landslides, based on two parameters: (1) the deformation rate associated to each building, measured from Sentinel-1 SAR data, and (2) the building damage generated by the landslide movement and recorded during a field campaign. We apply this method to Arcos de la Frontera, a monumental town in South Spain affected by a slow-moving landslide that has caused severe damage to buildings, forcing the evacuation of some of them. Our results show that maximum deformation rates of 4 cm/year in the line-of-sight (LOS) of the satellite, affects La Verbena, a newly-developed area, and displacements are mostly horizontal, as expected for a planar-landslide. Our building damage assessment reveals that most of the building blocks in La Verbena present moderate to severe damages. According to our vulnerability scale, 93% of the building blocks analysed present high vulnerability and, thus, should be the focus of more in-depth local studies to evaluate the serviceability of buildings, prior to adopting the necessary mitigation measures to reduce or cope with the negative consequences of this landslide. This methodology can be applied to slow-moving landslides worldwide thanks to the global availability of Sentinel-1 SAR data.

Description: IJGI, Vol. 6, Pages 268: Public Transit Route Mapping for Large-Scale Multimodal Networks ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090268 Authors: Flavio Poletti Patrick Bösch Francesco Ciari Kay Axhausen For the simulation of public transport, next to a schedule, knowledge of the public transport routes is required. While the schedules are becoming available, the precise network routes often remain unknown and must be reconstructed. For large-scale networks, however, a manual reconstruction becomes unfeasible. This paper presents a route reconstruction algorithm, which requires only the sequence and positions of the public transport stops and the street network. It uses an abstract graph to calculate the least-cost path from a route’s first to its last stop, with the constraint that the path must contain a so-called link candidate for every stop of the route’s stop sequence. The proposed algorithm is implemented explicitly for large-scale, real life networks. The algorithm is able to handle multiple lines and modes, to combine them at the same stop location (e.g., train and bus lines coming together at a train station), to automatically reconstruct missing links in the network, and to provide intelligent and efficient feedback if apparent errors occur. GPS or OSM tracks of the lines can be used to improve results, if available. The open-source algorithm has been tested for Zurich for mapping accuracy. In summary, the new algorithm and its MATSim-based implementation is a powerful, tested tool to reconstruct public transport network routes for large-scale systems.

Description: Remote Sensing, Vol. 9, Pages 887: Sea Wind Measurement by Doppler Navigation System with X-Configured Beams in Rectilinear Flight Remote Sensing doi: 10.3390/rs9090887 Authors: Alexey Nekrasov Alena Khachaturian Mária Gamcová Pavol Kurdel Viktor Obukhovets Vladimir Veremyev Mikhail Bogachev We suggest a conceptual approach to the measurement of the near-surface wind vector over water using a Doppler navigation system, in addition to its standard navigation capabilities. We consider a Doppler navigation system with a track-stabilized antenna and x-configuration of its beams. For the measurement of the sea-surface wind, the system operates in the multi-beam scatterometer mode in rectilinear flight. The proposed conceptual design has been validated, and its accuracy for the wind vector measurement has been estimated using Monte Carlo computational simulations.

Description: Remote Sensing, Vol. 9, Pages 888: The 2015 Surge of Hispar Glacier in the Karakoram Remote Sensing doi: 10.3390/rs9090888 Authors: Frank Paul Tazio Strozzi Thomas Schellenberger Andreas Kääb The Karakoram mountain range is well known for its numerous surge-type glaciers of which several have recently surged or are still doing so. Analysis of multi-temporal satellite images and digital elevation models have revealed impressive details about the related changes (e.g., in glacier length, surface elevation and flow velocities) and considerably expanded the database of known surge-type glaciers. One glacier that has so far only been reported as impacted by surging tributaries, rather than surging itself, is the 50 km long main trunk of Hispar Glacier in the Hunza catchment. We here present the evolution of flow velocities and surface features from its 2015/16 surge as revealed from a dense time series of SAR and optical images along with an analysis of historic satellite images. We observed maximum flow velocities of up to 14 m d−1 (5 km a−1) in spring 2015, sudden drops in summer velocities, a second increase in winter 2015/16 and a total advance of the surge front of about 6 km. During a few months the surge front velocity was much higher (about 90 m d−1) than the maximum flow velocity. We assume that one of its northern tributary glaciers, Yutmaru, initiated the surge at the end of summer 2014 and that the variability in flow velocities was driven by changes in the basal hydrologic regime (Alaska-type surge). We further provide evidence that Hispar Glacier has surged before (around 1960) over a distance of about 10 km so that it can also be regarded as a surge-type glacier.

Description: Remote Sensing, Vol. 9, Pages 883: Quantifying Snow Albedo Radiative Forcing and Its Feedback during 2003–2016 Remote Sensing doi: 10.3390/rs9090883 Authors: Lin Xiao Tao Che Linling Chen Hongjie Xie Liyun Dai Snow albedo feedback is one of the most crucial feedback processes that control equilibrium climate sensitivity, which is a central parameter for better prediction of future climate change. However, persistent large discrepancies and uncertainties are found in snow albedo feedback estimations. Remotely sensed snow cover products, atmospheric reanalysis data and radiative kernel data are used in this study to quantify snow albedo radiative forcing and its feedback on both hemispheric and global scales during 2003–2016. The strongest snow albedo radiative forcing is located north of 30°N, apart from Antarctica. In general, it has large monthly variation and peaks in spring. Snow albedo feedback is estimated to be 0.18 ± 0.08 W∙m−2∙°C−1 and 0.04 ± 0.02 W∙m−2∙°C−1 on hemispheric and global scales, respectively. Compared to previous studies, this paper focuses specifically on quantifying snow albedo feedback and demonstrates three improvements: (1) used high spatial and temporal resolution satellite-based snow cover data to determine the areas of snow albedo radiative forcing and feedback; (2) provided detailed information for model parameterization by using the results from (1), together with accurate description of snow cover change and constrained snow albedo and snow-free albedo data; and (3) effectively reduced the uncertainty of snow albedo feedback and increased its confidence level through the block bootstrap test. Our results of snow albedo feedback agreed well with other partially observation-based studies and indicate that the 25 Coupled Model Intercomparison Project Phase 5 (CMIP5) models might have overestimated the snow albedo feedback, largely due to the overestimation of surface albedo change between snow-covered and snow-free surface in these models.

Description: Remote Sensing, Vol. 9, Pages 882: Robust Feature Matching Method for SAR and Optical Images by Using Gaussian-Gamma-Shaped Bi-Windows-Based Descriptor and Geometric Constraint Remote Sensing doi: 10.3390/rs9090882 Authors: Min Chen Ayman Habib Haiqing He Qing Zhu Wei Zhang Improving the matching reliability of multi-sensor imagery is one of the most challenging issues in recent years, particularly for synthetic aperture radar (SAR) and optical images. It is difficult to deal with the noise influence, geometric distortions, and nonlinear radiometric difference between SAR and optical images. In this paper, a method for SAR and optical images matching is proposed. First, interest points that are robust to speckle noise in SAR images are detected by improving the original phase-congruency-based detector. Second, feature descriptors are constructed for all interest points by combining a new Gaussian-Gamma-shaped bi-windows-based gradient operator and the histogram of oriented gradient pattern. Third, descriptor similarity and geometrical relationship are combined to constrain the matching processing. Finally, an approach based on global and local constraints is proposed to eliminate outliers. In the experiments, SAR images including COSMO-Skymed, RADARSAT-2, TerraSAR-X and HJ-1C images, and optical images including ZY-3 and Google Earth images are used to evaluate the performance of the proposed method. The experimental results demonstrate that the proposed method provides significant improvements in the number of correct matches and matching precision compared with the state-of-the-art SIFT-like methods. Near 1 pixel registration accuracy is obtained based on the matching results of the proposed method.

Description: Remote Sensing, Vol. 9, Pages 884: Reducing the Effect of the Endmembers’ Spectral Variability by Selecting the Optimal Spectral Bands Remote Sensing doi: 10.3390/rs9090884 Authors: Omid Ghaffari Mohammad Javad Valadan Zoej Mehdi Mokhtarzade Variable environmental conditions cause different spectral responses of scene endmembers. Ignoring these variations affects the accuracy of fractional abundances obtained from linear spectral unmixing. On the other hand, the correlation between the bands of hyperspectral data is not considered by conventional methods developed for dealing with spectral variability. In this paper, a novel approach is proposed to simultaneously mitigate spectral variability and reduce correlation among different endmembers in hyperspectral datasets. The idea of the proposed method is to utilize the angular discrepancy of bands in the Prototype Space (PS), which is constructed using the endmembers of the image. Using the concepts of PS, in which each band is treated as a space point, we proposed a method to identify independent bands according to their angles. The proposed method comprised two main steps. In the first step, which aims to alleviate the spectral variability issue, image bands are prioritized based on their standard deviations computed over some sets of endmembers. Independent bands are then recognized in the prototype space, employing the angles between the prioritized bands. Finally, the unmixing process is done using the selected bands. In addition, the paper presents a technique to form a spectral library of endmembers’ variability (sets of endmembers). The proposed method extracts endmembers sets directly from the image data via a modified version of unsupervised spatial–spectral preprocessing. The performance of the proposed method was evaluated by five simulated images and three real hyperspectral datasets. The experiments show that the proposed method—using both groups of spectral variability reduction methods and independent band selection methods—produces better results compared to the conventional methods of each group. The improvement in the performance of the proposed method is observed in terms of more appropriate bands being selected and more accurate fractional abundance values being estimated.

Description: Remote Sensing, Vol. 9, Pages 881: Use of Proper Orthogonal Decomposition for Extraction of Ocean Surface Wave Fields from X-Band Radar Measurements of the Sea Surface Remote Sensing doi: 10.3390/rs9090881 Authors: Andrew Kammerer Erin Hackett Radar remote sensing of the sea surface for the extraction of ocean surface wave fields requires separating wave and non-wave contributions to the sea surface measurement. Conventional methods of extracting wave information from radar measurements of the sea surface rely on filtering the wavenumber-frequency spectrum using the linear dispersion relationship for ocean surface waves. However, this technique has limitations, e.g., it isn’t suited for the inclusion of non-linear wave features. This study evaluates an alternative method called proper orthogonal decomposition (POD) for the extraction of ocean surface wave fields remotely sensed by marine radar. POD is an empirical and optimal linear method for representing non-linear processes. The method was applied to Doppler velocity data of the sea surface collected using two different radar systems during two different experiments that spanned a variety of environmental conditions. During both experiments, GPS mini-buoys simultaneously collected wave data. The POD method was used to generate phase-resolved wave orbital velocity maps that are statistically evaluated by comparing wave statistics computed from the buoy data to those obtained from these maps. The results show that leading POD modes contain energy associated with the peak wavelength(s) of the measured wave field, and consequently, wave contributions to the radar measurement of the sea surface can be separated based on modes. Wave statistics calculated from optimized POD reconstructions are comparable to those calculated from GPS wave buoys. The accuracy of the wave statistics generated from POD-reconstructed orbital velocity maps was not sensitive to the radar configuration or environmental conditions examined. Further research is needed to determine a rigorous method for selecting modes a priori.

Description: Remote Sensing, Vol. 9, Pages 885: Comparing Fuzzy Sets and Random Sets to Model the Uncertainty of Fuzzy Shorelines Remote Sensing doi: 10.3390/rs9090885 Authors: Ratna Sari Dewi Wietske Bijker Alfred Stein This paper addresses uncertainty modelling of shorelines by comparing fuzzy sets and random sets. Both methods quantify extensional uncertainty of shorelines extracted from remote sensing images. Two datasets were tested: pan-sharpened Pleiades with four bands (Pleiades) and pan-sharpened Pleiades stacked with elevation data as the fifth band (Pleiades + DTM). Both fuzzy sets and random sets model the spatial extent of shoreline including its uncertainty. Fuzzy sets represent shorelines as a margin determined by upper and lower thresholds and their uncertainty as confusion indices. They do not consider randomness. Random sets fit the mixed Gaussian model to the image histogram. It represents shorelines as a transition zone between water and non-water. Their extensional uncertainty is assessed by the covering function. The results show that fuzzy sets and random sets resulted in shorelines that were closely similar. Kappa (κ) values were slightly different and McNemar’s test showed high p-values indicating a similar accuracy. Inclusion of the DTM (digital terrain model) improved the classification results, especially for roofs, inundated houses and inundated land. The shoreline model using Pleiades + DTM performed better than that of using Pleiades only, when using either fuzzy sets or random sets. It achieved κ values above 80%.

Description: Remote Sensing, Vol. 9, Pages 886: An Improved Tomography Approach Based on Adaptive Smoothing and Ground Meteorological Observations Remote Sensing doi: 10.3390/rs9090886 Authors: Bao Zhang Qingbiao Fan Yibin Yao Caijun Xu Xingxing Li Using the Global Navigation Satellite System (GNSS) to sense three-dimensional water vapor (WV) has been intensively investigated. However, this technique still heavily relies on the a priori information. In this study, we propose an improved tomography approach based on adaptive Laplacian smoothing (ALS) and ground meteorological observations. By using the proposed approach, the troposphere tomography is less dependent on a priori information and the ALS constraints match better with the actual situation than the constant constraints. Tomography experiments in Hong Kong during a heavy rainy period and a rainless period show that the ALS method gets superior results compared with the constant Laplacian smoothing (CLS) method. By validation with radiosonde and European Centre for Medium-Range Weather Forecasts (ECMWF) data, we found that the introduction of ground meteorological observations into tomography can solve the perennial problem of resolving the wet refractivity in the lower troposphere and thus significantly improve the tomography results. However, bad data quality and incompatibility of the ground meteorological observations may introduce errors into tomography results.

Description: Remote Sensing, Vol. 9, Pages 891: Linear and Non-Linear Trends for Seasonal NO2 and SO2 Concentrations in the Southern Hemisphere (2004−2016) Remote Sensing doi: 10.3390/rs9090891 Authors: Adrián Yuchechen Susan Lakkis Pablo Canziani In order to address the behaviour of nitrogen dioxide (NO2) and sulphur dioxide (SO2) in the context of a changing climate, linear and non-linear trends for the concentrations of these two trace gases were estimated over their seasonal standardised variables in the Southern Hemisphere—between the Equator and 60° S—using data retrieved by the Ozone Monitoring Instrument, for the period 2004–2016. A rescaling was applied to the calculated linear trends so that they are expressed in Dobson units (DU) per decade. Separately, the existence of monotonic—not necessarily linear—trends was addressed by means of the Mann-Kendall test. Results indicate that the SO2 exhibits significant linear trends in the planetary boundary layer only; they are present in all the analysed seasons but just in a small number of grid cells that are generally located over the landmasses or close to them. The SO2 concentrations in the quarterly time series exhibit, on average, a linear trend that is just below 0.08 DU decade−1 when significant and not significant values are considered altogether, but this figure increases to 0.80 DU decade−1 when only the significant trends are included. On the other hand, an important number of pixels in the lower troposphere, the middle troposphere, and the lower stratosphere have significant monotonic upward or downward trends. As for the NO2, no significant linear trends were found either in the troposphere or in the stratosphere, yet monotonic upward and downward trends were observed in the former and latter layers, respectively. Unlike the linear trends, semi-linear and non-linear trends were seen over the continents and in remote regions over the oceans. This suggests that pollutants are transported away from their sources by large-scale circulation and redistributed hemispherically. The combination of regional meteorological phenomena with atmospheric chemistry was raised as a possible explanation for the observed trends. If extrapolated, these trends are in an overall contradiction with the projected emissions of both gases for the current century.

Description: Remote Sensing, Vol. 9, Pages 894: Spatial Variability of L-Band Brightness Temperature during Freeze/Thaw Events over a Prairie Environment Remote Sensing doi: 10.3390/rs9090894 Authors: Alexandre Roy Peter Toose Chris Derksen Tracy Rowlandson Aaron Berg Juha Lemmetyinen Alain Royer Erica Tetlock Warren Helgason Oliver Sonnentag Passive microwave measurements from space are known to be sensitive to the freeze/thaw (F/T) state of the land surface. These measurements are at a coarse spatial resolution (~15–50 km) and the spatial variability of the microwave emissions within a pixel can have important effects on the interpretation of the signal. An L-band ground-based microwave radiometer campaign was conducted in the Canadian Prairies during winter 2014–2015 to examine the spatial variability of surface emissions during frozen and thawed periods. Seven different sites within the Kenaston soil monitoring network were sampled five times between October 2014 and April 2015 with a mobile ground-based L-band radiometer system at approximately monthly intervals. The radiometer measurements showed that in a seemingly homogenous prairie landscape, the spatial variability of brightness temperature (TB) is non-negligible during both frozen and unfrozen soil conditions. Under frozen soil conditions, TB was negatively correlated with soil permittivity (εG). This correlation was related to soil moisture conditions before the main freezing event, showing that the soil ice volumetric content at least partly affects TB. However, because of the effect of snow on L-Band emission, the correlation between TB and εG decreased with snow accumulation. When compared to satellite measurements, the average TB of the seven plots were well correlated with the Soil Moisture Ocean Salinity (SMOS) TB with a root mean square difference of 8.1 K and consistent representation of the strong F/T signal (i.e., TB increases and decreases when soil freezing and thawing, respectively). This study allows better quantitative understanding of the spatial variability in L-Band emissions related to landscape F/T, and will help the calibration and validation of satellite-based F/T retrieval algorithms.

Description: Remote Sensing, Vol. 9, Pages 915: Diurnal Air Temperature Modeling Based on the Land Surface Temperature Remote Sensing doi: 10.3390/rs9090915 Authors: Mehdi Gholamnia Seyed Kazem Alavipanah Ali Darvishi Boloorani Saeid Hamzeh Majid Kiavarz The air temperature is an essential variable in many applications related to Earth science. Sporadic spatial distribution of weather stations causes a low spatial resolution of measured air temperatures. This study focused on modeling the air diurnal temperature cycle (DTC) based on the land surface temperature (LST) DTC. The air DTC model parameters were estimated from LST DTC model parameters by a regression analysis. Here, the LST obtained from the INSAT-3D geostationary satellite and the air temperature extracted from weather stations were used within the time frame of 4 March 2015 to 22 May 2017 across Iran. Constant parameters of the air DTC model for each weather station were estimated based on an experimental approach over the time period. Results showed these parameters decrease as elevation increases. The mean absolute error (MAE) and the root mean square error (RMSE) for three hours sampling were calculated. The MAE and RMSE ranges were between [0.1, 4] °C and [0.1, 3.3] °C, respectively. Additionally, 95% of MAEs and RMSEs were less than 2.9 °C and 2.4 °C values, correspondingly. The range of the mean values of MAEs and RMSEs for a three-hour sampling time were [−0.29, 0.6] °C and [2, 2.11] °C. The DTC model results showed a meaningful statistical fitting in both air DTCs modeled from LST and weather station-based DTCs. The variability of mean error and RMSE in different land covers and elevation classes were also investigated. In spite of the complex behavior of the environmental variables in the study area, the model error bar did not show significantly biased estimations for various classes. Therefore, the developed model was less sensitive to variations of land covers and elevation changes. It can be conclude that the coefficients of regression between LST and air DTC could model properly the environmental factors.

Description: IJGI, Vol. 6, Pages 279: Making Spatial Decisions Using ArcGIS Pro: A Workbook. By Kathryn Keranen and Robert Kolvoord, Esri Press, 2017; 376 Pages. Price $69.99, ISBN 9781589484849 ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090279 Authors: IJGI Editorial Office n/a

Description: IJGI, Vol. 6, Pages 281: Mapping Parallels between Outdoor Urban Environments and Indoor Manufacturing Environments ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090281 Authors: Stefan Schabus Johannes Scholz Thomas Lampoltshammer The concepts of “Smart Cities” and “Smart Manufacturing” are different data-driven domains, although both rely on intelligent information technology and data analysis. With the application of linked data and affordance-based approaches, both domains converge, paving the way for new and innovative viewpoints regarding the comparison of urban tasks with indoor manufacturing tasks. The present study builds on the work, who state that cities are scaled versions of each other, by extending this thesis towards indoor manufacturing environments. Based on their structure and complexity, these environments are considered to form ecosystems of their own, comparable to “small cities”. This conceptual idea is demonstrated by examining the process of human problem-solving in transportation situations from both perspectives (i.e., city-level and manufacturing-level). In particular, the authors model tasks of human operators that are used to support transportation processes in indoor manufacturing environments based on affordances and spatial-temporal data. This paper introduces the fundamentals of the transformation process of outdoor tasks and process planning activities to indoor environments, particularly to semiconductor manufacturing environments. The idea is to examine the mapping of outdoor tasks and applications to indoor environments, and vice-versa, based on an example focusing on the autonomous transportation of production assets in a manufacturing environment. The approach is based on a spatial graph database, populated with an indoor navigation ontology and instances of indoor and outdoor objects. The results indicate that human problem-solving strategies can be applied to indoor manufacturing environments to support decision-making in autonomous transportation tasks.

Description: Remote Sensing, Vol. 9, Pages 939: Feature Selection Solution with High Dimensionality and Low-Sample Size for Land Cover Classification in Object-Based Image Analysis Remote Sensing doi: 10.3390/rs9090939 Authors: Yaohuan Huang Chuanpeng Zhao Haijun Yang Xiaoyang Song Jie Chen Zhonghua Li Land cover information extraction through object-based image analysis (OBIA) has become an important trend in remote sensing, thanks to the increasing availability of high-resolution imagery. Segmented objects have a large number of features that cause high-dimension and low-sample size problems in the classification process. In this study, on the basis of a partial least squares generalized linear regression (PLSGLR), we propose a group corrected PLSGLR, known as G-PLSGLR, that aims to reduce the redundancy of object features for land cover identifications. Using Gaofen-2 images, the area of interest was segmented and sampled to generate small sample-size training datasets with 51 object features. The features selected by G-PLSGLR were compared against a guided regularized random forest (GRRF) in metrics of reduction rate, feature redundancy, and accuracy assessment of classification. Three indicators of overall accuracy (OA), user’s accuracy (UA), and producer’s accuracy (PA) were applied for accuracy assessment in this paper. The result shows that the G-PLSGLR achieved a reduction rate of 9.27 with a feature redundancy of 0.29, and a value of OA 90.63%. The GRRF achieved a reduction rate of 1.61 with a feature redundancy of 0.42, and a value of OA 85.56%. The PA of each land cover category was more than 95% using features selected by G-PLSGLR, while the PA ranged from 77 to 96% using features selected by GRRF. The UA of G-PLSGLR-selected features ranged from 70 to 80% except for grass land and bare land, which achieved 10% higher UA than GRRF-selected features. The G-PLSGLR method we proposed has the advantages of a large reduction rate, low feature redundancy, and high classification performance, which can be applied in OBIA-based land cover classification.

Description: Remote Sensing, Vol. 9, Pages 947: Circum-Arctic Changes in the Flow of Glaciers and Ice Caps from Satellite SAR Data between the 1990s and 2017 Remote Sensing doi: 10.3390/rs9090947 Authors: Tazio Strozzi Frank Paul Andreas Wiesmann Thomas Schellenberger Andreas Kääb We computed circum-Arctic surface velocity maps of glaciers and ice caps over the Canadian Arctic, Svalbard and the Russian Arctic for at least two times between the 1990s and 2017 using satellite SAR data. Our analyses are mainly performed with offset-tracking of ALOS-1 PALSAR-1 (2007–2011) and Sentinel-1 (2015–2017) data. In certain cases JERS-1 SAR (1994–1998), TerraSAR-X (2008–2012), Radarsat-2 (2009–2016) and ALOS-2 PALSAR-2 (2015–2016) data were used to fill-in spatial or temporal gaps. Validation of the latest Sentinel-1 results was accomplished by means of SAR data at higher spatial resolution (Radarsat-2 Wide Ultra Fine) and ground-based measurements. In general, we observe a deceleration of flow velocities for the major tidewater glaciers in the Canadian Arctic and an increase in frontal velocity along with a retreat of frontal positions over Svalbard and the Russian Arctic. However, all regions have strong accelerations for selected glaciers. The latter developments can be well traced based on the very high temporal sampling of Sentinel-1 acquisitions since 2015, revealing new insights in glacier dynamics. For example, surges on Spitsbergen (e.g., Negribreen, Nathorsbreen, Penckbreen and Strongbreen) have a different characteristic and timing than those over Eastern Austfonna and Edgeoya (e.g., Basin 3, Basin 2 and Stonebreen). Events similar to those ongoing on Eastern Austofonna were also observed over the Vavilov Ice Cap on Severnaya Zemlya and possibly Simony Glacier on Franz-Josef Land. Collectively, there seems to be a recently increasing number of glaciers with frontal destabilization over Eastern Svalbard and the Russian Arctic compared to the 1990s.

Description: IJGI, Vol. 6, Pages 290: Selecting the Best Band Ratio to Estimate Chlorophyll-a Concentration in a Tropical Freshwater Lake Using Sentinel 2A Images from a Case Study of Lake Ba Be (Northern Vietnam) ISPRS International Journal of Geo-Information doi: 10.3390/ijgi6090290 Authors: Nguyen Thi Thu Ha Nguyen Thien Phuong Thao Katsuaki Koike Mai Trong Nhuan This study aims to develop a method to estimate chlorophyll-a concentration (Chla) in tropical freshwater lake waters using in situ data of Chla, water reflectance, and concurrent Sentinel 2A MSI imagery (S2A) over Lake Ba Be, a Ramsar site and the largest natural freshwater lake in Vietnam. Data from 30 surveyed sampling sites over the lake water in June 2016 and May 2017 demonstrated the appropriateness of S2A green-red band ratio (band 3 versus band 4) for estimating Chla. This was shown through a strong correlation of corresponded field measured reflectance ratio with Chla by an exponential curve (r2 = 0.68; the mean standard error of the estimates corresponding to 5% of the mean value of in situ Chla). The small error between in situ Chla, and estimated Chla from S2A acquired concurrently, confirmed the S2A green-red band ratio as the most suitable option for monitoring Chla in Lake Ba Be water. Resultant Chla distribution maps over time described a partially-seasonal pattern and also displayed the spatial dynamic of Chla in the lake. This allows a better understanding of the lake’s limnological processes to be developed and provides an insight into the factors that affect lake water quality. The results also confirmed the potential of S2A to be used as a free tool for lake monitoring and research due to high spatial resolution data (10 m pixel size).

Description: Remote Sensing, Vol. 9, Pages 830: An Adaptive Offset Tracking Method with SAR Images for Landslide Displacement Monitoring Remote Sensing doi: 10.3390/rs9080830 Authors: Jiehua Cai Changcheng Wang Xiaokang Mao Qijie Wang With the development of high-resolution Synthetic Aperture Radar (SAR) systems, researchers are increasingly paying attention to the application of SAR offset tracking methods in ground deformation estimation. The traditional normalized cross correlation (NCC) tracking method is based on regular matching windows. For areas with different moving characteristics, especially the landslide boundary areas, the NCC method will produce incorrect results. This is because in landslide boundary areas, the pixels of the regular matching window include two or more types of moving characteristics: some pixels with large displacement, and others with small or no displacement. These two kinds of pixels are uncorrelated, which result in inaccurate estimations. This paper proposes a new offset tracking method with SAR images based on the adaptive matching window to improve the accuracy of landslide displacement estimation. The proposed method generates an adaptive matching window that only contains pixels with similar moving characteristics. Three SAR images acquired by the Jet Propulsion Laboratory’s Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) system are selected to estimate the surface deformation of the Slumgullion landslide located in the southwestern Colorado, USA. The results show that the proposed method has higher accuracy than the traditional NCC method, especially in landslide boundary areas. Furthermore, it can obtain more detailed displacement information in landslide boundary areas.

Description: Remote Sensing, Vol. 9, Pages 831: Identifying Droughts Affecting Agriculture in Africa Based on Remote Sensing Time Series between 2000–2016: Rainfall Anomalies and Vegetation Condition in the Context of ENSO Remote Sensing doi: 10.3390/rs9080831 Authors: Karina Winkler Ursula Gessner Volker Hochschild Droughts are amongst the most destructive natural disasters in the world. In large regions of Africa, where water is a limiting factor and people strongly rely on rain-fed agriculture, droughts have frequently led to crop failure, food shortages and even humanitarian crises. In eastern and southern Africa, major drought episodes have been linked to El Niño-Southern Oscillation (ENSO) events. In this context and with limited in-situ data available, remote sensing provides valuable opportunities for continent-wide assessment of droughts with high spatial and temporal resolutions. This study aimed to monitor agriculturally relevant droughts over Africa between 2000–2016 with a specific focus on growing seasons using remote sensing-based drought indices. Special attention was paid to the observation of drought dynamics during major ENSO episodes to illuminate the connection between ENSO and droughts in eastern and southern Africa. We utilized Tropical Rainfall Measuring Mission (TRMM)-based Standardized Precipitation Index (SPI) with 0 . 25 ∘ resolution and Moderate-resolution Imaging Spectroradiometer (MODIS)-derived Vegetation Condition Index (VCI) with 500 m resolution as indices for analysing the spatio-temporal patterns of droughts. We combined the drought indices with information on the timing of site-specific growing seasons derived from MODIS-based multi-annual average of Normalized Difference Vegetation Index (NDVI). We proved the applicability of SPI-3 and VCI as indices for a comprehensive continental-scale monitoring of agriculturally relevant droughts. The years 2009 and 2011 could be revealed as major drought years in eastern Africa, whereas southern Africa was affected by severe droughts in 2003 and 2015/2016. Drought episodes occurred over large parts of southern Africa during strong El Niño events. We observed a mixed drought pattern in eastern Africa, where areas with two growing seasons were frequently affected by droughts during La Niña and zones of unimodal rainfall regimes showed droughts during the onset of El Niño. During La Niña 2010/2011, large parts of cropland areas in Somalia (88%), Sudan (64%) and South Sudan (51%) were affected by severe to extreme droughts during the growing seasons. However, no universal El Niño- or La Niña-related response pattern of droughts could be deduced for the observation period of 16 years. In this regard, we discussed multi-year atmospheric fluctuations and characteristics of ENSO variants as further influences on the interconnection between ENSO and droughts. By utilizing remote sensing-based drought indices focussed on agricultural zones and periods, this study attempts to contribute to a better understanding of spatio-temporal patterns of droughts affecting agriculture in Africa, which can be essential for implementing strategies of drought hazard mitigation.