ALBERT

Description: In order to provide a good theoretical guidance for the development and utilization of weathered phosphorite resources, we investigated the geochemical and mineralogical characteristics of primary and weathered phosphorites. The analysis of trace elements showed that the primary ore has hydrothermal sedimentation effect in the later stage, the weathered ore has obvious residual enrichment and the phosphate ore belongs to clastic lithologic phosphate rock. In addition, through leaching test method, it was shown that rare earth elements are present in fluorapatite in the form of isomorphic substitution, and the proportion of rare earth elements adsorbed on clay and other minerals was likely to be between 2% and 3%. The light rare earth elements are relatively enriched in both primary and weathered phosphorite, and Ce and Eu have obvious negative anomalies. The primary phosphorite is a dolomitic phosphorite containing rare earth elements, which are naturally enriched by weathering, and its weathered ore has obvious residual enrichment, while the deposit was characterized by normal marine sedimentation and hydrothermal action.

Description: Haiyang-2A (HY-2A) has been working in-flight for over seven years, and the accuracy of HY-2A calibration microwave radiometer (CMR) data is extremely important for the wet troposphere delay correction (WTC) in sea surface height (SSH) determination. We present a comprehensive evaluation of the HY-2A CMR observation using the numerical weather model (NWM) for all the data available period from October 2011 to February 2018, including the WTC and the precipitable water vapor (PWV). The ERA(ECMWF Re-Analysis)-Interim products from European Centre for Medium-Range Weather Forecasts (ECMWF) are used for the validation of HY-2A WTC and PWV products. In general, a global agreement of root-mean-square (RMS) of 2.3 cm in WTC and 3.6 mm in PWV are demonstrated between HY-2A observation and ERA-Interim products. Systematic biases are revealed where before 2014 there was a positive WTC/PWV bias and after that, a negative one. Spatially, HY-2A CMR products show a larger bias in polar regions compared with mid-latitude regions and tropical regions and agree better in the Antarctic than in the Arctic with NWM. Moreover, HY-2A CMR products have larger biases in the coastal area, which are all caused by the brightness temperature (TB) contamination from land or sea ice. Temporally, the WTC/PWV biases increase from October 2011 to March 2014 with a systematic bias over 1 cm in WTC and 2 mm in PWV, and the maximum RMS values of 4.62 cm in WTC and 7.61 mm in PWV occur in August 2013, which is because of the unsuitable retrieval coefficients and systematic TB measurements biases from 37 GHz band. After April 2014, the TB bias is corrected, HY-2A CMR products agree very well with NWM from April 2014 to May 2017 with the average RMS of 1.68 cm in WTC and 2.65 mm in PWV. However, since June 2017, TB measurements from the 18.7 GHz band become unstable, which led to the huge differences between HY-2A CMR products and the NWM with an average RMS of 2.62 cm in WTC and 4.33 mm in PWV. HY-2A CMR shows high accuracy when three bands work normally and further calibration for HY-2A CMR is in urgent need. Furtherly, 137 global coastal radiosonde stations were used to validate HY-2A CMR. The validation based on radiosonde data shows the same variation trend in time of HY-2A CMR compared to the results from ECMWF, which verifies the results from ECMWF.

Description: Crop yield estimation at a regional scale over a long period of time is of great significance to food security. In past decades, the integration of remote sensing observations and crop growth models has been recognized as a promising approach for crop growth monitoring and yield estimation. Optical remote sensing data are susceptible to cloud and rain, while synthetic aperture radar (SAR) can penetrate through clouds and has all-weather capabilities. This allows for more reliable and consistent crop monitoring and yield estimation in terms of radar sensor data. The aim of this study is to improve the accuracy for winter wheat yield estimation by assimilating time series soil moisture images, which are retrieved by a water cloud model using SAR and optical data as input, into the crop model. In this study, SAR images were acquired by C-band SAR sensors boarded on Sentinel-1 satellites and optical images were obtained from a Sentinel-2 multi-spectral instrument (MSI) for Hengshui city of Hebei province in China. Remote sensing data and ground data were all collected during the main growing season of winter wheat. Both the normalized difference vegetation index (NDVI), derived from Sentinel-2, and backscattering coefficients and polarimetric indicators, computed from Sentinel-1, were used in the water cloud model to derive time series soil moisture (SM) images. To improve the prediction of crop yields at the field scale, we incorporated remotely sensed soil moisture into the World Food Studies (WOFOST) model using the Ensemble Kalman Filter (EnKF) algorithm. In general, the trend of soil moisture inversion was consistent with the ground measurements, with the coefficient of determination (R2) equal to 0.45, 0.53, and 0.49, respectively, and RMSE was 9.16%, 7.43%, and 8.53%, respectively, for three observation dates. The winter wheat yield estimation results showed that the assimilation of remotely sensed soil moisture improved the correlation of observed and simulated yields (R2 = 0.35; RMSE =934 kg/ha) compared to the situation without data assimilation (R2 = 0.21; RMSE = 1330 kg/ha). Consequently, the results of this study demonstrated the potential and usefulness of assimilating SM retrieved from both Sentinel-1 C-band SAR and Sentinel-2 MSI optical remote sensing data into WOFOST model for winter wheat yield estimation and could also provide a reference for crop yield estimation with data assimilation for other crop types.

Description: Urban areas globally are vulnerable to warming climate trends exacerbated by their growing populations and heat island effects. The Local Climate Zone (LCZ) typology has become a popular framework for characterizing urban microclimates in different regions using various classification methods, including a widely adopted pixel-based protocol by the World Urban Database and Access Portal Tools (WUDAPT) Project. However, few studies to date have explored the potential of object-based image analysis (OBIA) to facilitate classification of LCZs given their inherent complexity, and few studies have further used the LCZ framework to analyze land cover changes in urban areas over time. This study classified LCZs in the Salt Lake Metro Region, Utah, USA for 1993 and 2017 using a supervised object-based analysis of Landsat satellite imagery and assessed their change during this time frame. The overall accuracy, measured for the most recent classification period (2017), was equal to 64% across 12 LCZs, with most of the error resulting from similarities among highly developed LCZs and non-developed classes with sparse or low-stature vegetation. The observed 1993–2017 changes in LCZs indicated a regional tendency towards primarily suburban, open low-rise development, and large low-rise and paved classes. However, despite the potential for local cooling with landscape transitions likely to increase vegetation cover and irrigation compared to pre-development conditions, summer averages of Landsat-derived top-of-atmosphere brightness temperatures showed a pronounced warming between 1992–1994 and 2016–2018 across the study region, with a 0.1–2.9 °C increase among individual LCZs. Our results indicate that future applications of LCZs towards urban change analyses should develop a stronger understanding of LCZ microclimate sensitivity to changes in size and configuration of urban neighborhoods and regions. Furthermore, while OBIA is promising for capturing the heterogeneous and multi-scale nature of LCZs, its applications could be strengthened by adopting more generalizable approaches for LCZ-relevant segmentation and validation, and by incorporating active remote sensing data to account for the 3D complexity of urban areas.

Description: An important application of airborne- and satellite-based hyperspectral imaging is the mapping of the spatial distribution of vegetation biophysical and biochemical parameters in an environment. Statistical models, such as Gaussian processes, have been very successful for modeling vegetation parameters from captured spectra, however their performance is highly dependent on the amount of available ground truth. This is a problem because it is generally expensive to obtain ground truth information due to difficulties and costs associated with sample collection and analysis. In this paper, we present two Gaussian processes based approaches for improving the accuracy of vegetation parameter retrieval when ground truth is limited. The first is the adoption of covariance functions based on well-established metrics, such as, spectral angle and spectral correlation, which are known to be better measures of similarity for spectral data owing to their resilience to spectral variabilities. The second is the joint modeling of related vegetation parameters by multitask Gaussian processes so that the prediction accuracy of the vegetation parameter of interest can be improved with the aid of related vegetation parameters for which a larger set of ground truth is available. We experimentally demonstrate the efficacy of the proposed methods against existing approaches on three real-world hyperspectral datasets and one synthetic dataset.

Description: Spatial features retrieved from satellite data play an important role for improving crop classification. In this study, we proposed a deep-learning-based time-series analysis method to extract and organize spatial features to improve parcel-based crop classification using high-resolution optical images and multi-temporal synthetic aperture radar (SAR) data. Central to this method is the use of multiple deep convolutional networks (DCNs) to extract spatial features and to use the long short-term memory (LSTM) network to organize spatial features. First, a precise farmland parcel map was delineated from optical images. Second, hundreds of spatial features were retrieved using multiple DCNs from preprocessed SAR images and overlaid onto the parcel map to construct multivariate time-series of crop growth for parcels. Third, LSTM-based network structures for organizing these time-series features were constructed to produce a final parcel-based classification map. The method was applied to a dataset of high-resolution ZY-3 optical images and multi-temporal Sentinel-1A SAR data to classify crop types in the Hunan Province of China. The classification results, showing an improvement of greater than 5.0% in overall accuracy relative to methods without spatial features, demonstrated the effectiveness of the proposed method in extracting and organizing spatial features for improving parcel-based crop classification.

Description: The classification of very-high-resolution (VHR) remote sensing images is essential in many applications. However, high intraclass and low interclass variations in these kinds of images pose serious challenges. Fully convolutional network (FCN) models, which benefit from a powerful feature learning ability, have shown impressive performance and great potential. Nevertheless, only classification results with coarse resolution can be obtained from the original FCN method. Deep feature fusion is often employed to improve the resolution of outputs. Existing strategies for such fusion are not capable of properly utilizing the low-level features and considering the importance of features at different scales. This paper proposes a novel, end-to-end, fully convolutional network to integrate a multiconnection ResNet model and a class-specific attention model into a unified framework to overcome these problems. The former fuses multilevel deep features without introducing any redundant information from low-level features. The latter can learn the contributions from different features of each geo-object at each scale. Extensive experiments on two open datasets indicate that the proposed method can achieve class-specific scale-adaptive classification results and it outperforms other state-of-the-art methods. The results were submitted to the International Society for Photogrammetry and Remote Sensing (ISPRS) online contest for comparison with more than 50 other methods. The results indicate that the proposed method (ID: SWJ_2) ranks #1 in terms of overall accuracy, even though no additional digital surface model (DSM) data that were offered by ISPRS were used and no postprocessing was applied.

Description: Low Earth orbit (LEO) satellites play an important role in human space activities, and market demands for commercial uses of LEO satellites have been increasing rapidly in recent years. LEO satellites mainly consist of Earth observation satellites (EOSs), the major commercial applications of which are various sorts of Earth observations, such as map making, crop growth assessment, and disaster surveillance. However, the success rates of observation tasks are influenced considerably by uncertainties in local weather conditions, inadequate sunlight, observation dip angle, and other practical factors. The available time windows (ATWs) suitable for observing given types of targets and for transmitting data back to ground receiver stations are relatively narrow. In order to utilize limited satellite resources efficiently and maximize their commercial benefits, it is necessary to evaluate the overall effectiveness of satellites and planned tasks considering various factors. In this paper, we propose a method for determining the ATWs considering the influence of sunlight angle, elevation angle, and the type of sensor equipped on the satellite. After that, we develop a satellite effectiveness evaluation (SEE) model for satellite observation and data-downlink scheduling (SODS) based on the Availability–Capacity–Profitability (ACP) framework, which is designed to evaluate the overall performance of satellites from the perspective of time resource utilization, the success rate of tasks, and profit return. The effects of weather uncertainties on the tasks’ success are considered in the SEE model, and the model can be applied to support the decision-makers on optimizing and improving task arrangements for EOSs. Finally, a case study is presented to demonstrate the effectiveness of the proposed method and verify the ACP-based SEE model. The obtained ATWs by the proposed method are compared with those by the Systems Tool Kit (STK), and the correctness of the method is thus validated.

Description: In real-time onboard hyperspectral-image(HSI) anomalous targets detection, processing speed and accuracy are equivalently desirable which is hard to satisfy at the same time. To improve detection accuracy, deep learning based HSI anomaly detectors (ADs) are widely studied. However, their large scale network results in a massive computational burden. In this paper, to improve the detection throughput without sacrificing the accuracy, a pruning–quantization–anomaly–detector (P-Q-AD) is proposed by building an underlying constraint formulation to make a trade-off between accuracy and throughput. To solve this formulation, multi-objective optimization with nondominated sorting genetic algorithm II (NSGA-II) is employed to shrink the network. As a result, the redundant neurons are removed. A mixed precision network is implemented with a delicate customized fixed-point data expression to further improve the efficiency. In the experiments, the proposed P-Q-AD is implemented on two real HSI data sets and compared with three types of detectors. The results show that the performance of the proposed approach is no worse than those comparison detectors in terms of the receiver operating characteristic curve (ROC) and area under curve (AUC) value. For the onboard mission, the proposed P-Q-AD reaches over 4.5 × speedup with less than 0.5 % AUC loss compared with the floating-based detector. The pruning and the quantization approach in this paper can be referenced for designing the anomalous targets detectors for high efficiency.

Description: To date, several algorithms for the retrieval of cyanobacterial phycocyanin (PC) from ocean colour sensors have been presented for inland waters, all of which claim to be robust models. To address this, we conducted a comprehensive comparison to identify the optimal algorithm for retrieval of PC concentrations in the highly optically complex waters of Lake Balaton (Hungary). MEdium Resolution Imaging Spectrometer (MERIS) top-of-atmosphere radiances were first atmospherically corrected using the Self-Contained Atmospheric Parameters Estimation for MERIS data v.B2 (SCAPE-M_B2). Overall, the Simis05 semi-analytical algorithm outperformed more complex inversion algorithms, providing accurate estimates of PC up to ±7 days from the time of satellite overpass during summer cyanobacteria blooms (RMSElog 〈 0.33). Same-day retrieval of PC also showed good agreement with cyanobacteria biomass (R2 〉 0.66, p 〈 0.001). In-depth analysis of the Simis05 algorithm using in situ measurements of inherent optical properties (IOPs) revealed that the Simis05 model overestimated the phytoplankton absorption coefficient [aph(λ)] by a factor of ~2. However, these errors were compensated for by underestimation of the mass-specific chlorophyll absorption coefficient [a*chla(λ)]. This study reinforces the need for further validation of algorithms over a range of optical water types in the context of the recently launched Ocean Land Colour Instrument (OLCI) onboard Sentinel-3.

Description: In recent years, with more open data platforms and tools available to store and process satellite imagery, Earth Observation data have become widely accessible and usable especially for countries previously not in the possession of tasking rights to satellites and the needed processing capacity. Due to its ideal scanning and acquisition conditions for low cloud coverage imagery, Namibia aims to make use of this new development and integrate Earth Observation data into its national monitoring system of sustainable development goals (SDG). The purpose of this study is to assess the potential of open source tools and global datasets to estimate the national SDG indicators on Change of water-related ecosystems (6.6.1), Rural population with access to roads (9.1.1), Forest coverage (15.1.1) and Land degradation (15.3.1). The results are set into perspective of existing information in each particular sector. The study shows that, in the absence of in-situ measurements or data collected through surveys, the Earth Observation-based results represent a high potential to supplement the national statistics for Namibia or to serve as primary data sources once validated through ground-truthing. Furthermore, examples are given for the limitations of the assessed Earth Observation solutions in the context of Namibia. Hence, the study also serves as valuable input for discussions on a consensus on national definitions and standards by all stakeholders responsible for releasing official statistics.

Description: Despite many efforts of the radar community, quantitative precipitation estimation (QPE) from weather radar data remains a challenging topic. The high resolution of X-band radar imagery in space and time comes with an intricate correction process of reflectivity. The steep and high mountain topography of the Andes enhances its complexity. This study aims to optimize the rainfall derivation of the highest X-band radar in the world (4450 m a.s.l.) by using a random forest (RF) model and single Plan Position Indicator (PPI) scans. The performance of the RF model was evaluated in comparison with the traditional step-wise approach by using both, the Marshall-Palmer and a site-specific Z–R relationship. Since rain gauge networks are frequently unevenly distributed and hardly available at real time in mountain regions, bias adjustment was neglected. Results showed an improvement in the step-wise approach by using the site-specific (instead of the Marshall-Palmer) Z–R relationship. However, both models highly underestimate the rainfall rate (correlation coefficient 〈 0.69; slope up to 12). Contrary, the RF model greatly outperformed the step-wise approach in all testing locations and on different rainfall events (correlation coefficient up to 0.83; slope = 1.04). The results are promising and unveil a different approach to overcome the high attenuation issues inherent to X-band radars.

Description: Microbes can mediate the precipitation of primary dolomite under surface conditions. Meanwhile, primary dolomite mediated by microbes often contains more Fe2+ than standard dolomite in modern microbial culture experiments. Ferroan dolomite and ankerite have been regarded as secondary products. This paper reviews the process and possible mechanisms of microbial mediated precipitation of primary ferroan dolomite and/or ankerite. In the microbial geochemical Fe cycle, many dissimilatory iron-reducing bacteria (DIRB), sulfate-reducing bacteria (SRB), and methanogens can reduce Fe3+ to Fe2+, while SRB and methanogens can also promote the precipitation of primary dolomite. There are an oxygen respiration zone (ORZ), an iron reduction zone (IRZ), a sulfate reduction zone (SRZ), and a methanogenesis zone (MZ) from top to bottom in the muddy sediment diagenesis zone. DIRB in IRZ provide the lower section with Fe2+, which composes many enzymes and proteins to participate in metabolic processes of SRB and methanogens. Lastly, heterogeneous nucleation of ferroan dolomite on extracellular polymeric substances (EPS) and cell surfaces is mediated by SRB and methanogens. Exploring the origin of microbial ferroan dolomite may help to solve the “dolomite problem”.

Description: Pyrometallurgical processing of ore from the Zeehan mineral field was performed intermittently between 1896 and 1948, primarily recovering Pb, Ag and Cu. While Zn recovery was attempted at the time, it was unsuccessful using the available technology. Consequently, Zn reported to the slag during the smelting process. Today, the former smelter site consists of two large slag piles (North and South). Using a range of techniques (including X-ray diffractometry, scanning electron microscopy, laser ablation inductively coupled plasma mass spectrometry, and static testing) the geometallurgical and geo-environmental properties of these slag materials (n = 280) were determined. The South and North piles contain on average 15% and 11% Zn, respectively. A range of complex mineral phases were identified, and are dominated by glass, silicates (i.e., monticellite–kirschsteinite and hardystonite), oxides (gahnite and hercynite) and minor sulfides (sphalerite and wurtzite). Microtextural examinations defined nine mineral phases (Glass A, Silicates A to D, Oxides A and B, Sulfides A and B). Zn was concentrated in Sulfide A (26%), Glass A (24%) and the Silicates (43%), while Pb was concentrated in Oxide B (76%), with Sulfide B host to the highest Ag (45%) and Cu (65%). Considering this, recovery of Zn using conventional hydrometallurgical processes (i.e., sulfuric acid leaching) is suitable, however the application of unconventional biohydrometallurgical techniques could be explored, as well re-smelting. These slag materials are classified geo-environmentally as potentially acid forming, with leachate concentrations of Zn, Pb consistently above ANZECC (2000) aquatic ecosystem 80% protection guideline values, and, for the majority of samples, exceedances of Cu, Ni and Cd were also measured. Considering these findings, reprocessing of these historic slags for Zn extraction may provide an economically feasible management option for rehabilitating this historical site.

Description: New findings of silicate-melt inclusions in two alluvial diamonds (from the Kholomolokh placer, northeastern Siberian Platform) are reported. Both diamonds exhibit a high degree of N aggregation state (60–70% B) suggesting their long residence in the mantle. Raman spectral analysis revealed that the composite inclusions consist of clinopyroxene and silicate glass. Hopper crystals of clinopyroxene were observed using scanning electron microscopy and energy-dispersive spectroscopic analyses; these are different in composition from the omphacite inclusions that co-exist in the same diamonds. The glasses in these inclusions contain relatively high SiO2, Al2O3, Na2O and, K2O. These composite inclusions are primary melt that partially crystallised at the cooling stage. Hopper crystals of clinopyroxene imply rapid cooling rates, likely related to the uplift of crystals in the kimberlite melt. The reconstructed composition of such primary melts suggests that they were formed as the product of metasomatised mantle. One of the most likely source of melts/fluids metasomatising the mantle could be a subducted slab.

Description: The nature of upper mantle is important to understand the evolution of the South China Sea (SCS); thus, we need better constrains on its mantle heterogeneity. Magma water concentration is a good indicator, but few data have been reported. However, the rarity of glass and melt inclusions and the special genesis for phenocrysts in SCS basalts present challenges to analyzing magmatic water content. Therefore, it is possible to estimate the water variations through the characteristics of partial melting and magma crystallization. We evaluated variations in Fe depletion, degree of melt fractions, and mantle source composition along the fossil spreading ridge (FSR) using SCS basalt data from published papers. We found that lava from the FSR 116.2° E, FSR 117.7° E, and non-FSR regions can be considered normal lava with normal water content; in contrast, lava from the FSR 117° E-carbonatite and 114.9–115.0° E basalts have higher water content and show evidence of strong Fe depletion during the fractional crystallization after elimination of the effects of plagioclase oversaturation. The enriched water in the 117° E-carbonatite basalts is contained in carbonated silicate melts, and that in the 114.9–115.0° E basalts results from mantle contamination with the lower continental crust. The lava from the 117° E-normal basalt has much lower water content because of the lesser influence of the Hainan plume. Therefore, there must be a mantle source compositional transition area between the southwestern and eastern sub-basins of the SCS, which have different mantle evolution histories. The mantle in the west is more affected by contamination with continental materials, while that in the east is more affected by the Hainan mantle plume.

Description: The large-scale Maoping W–Sn deposit in the Gannan metallogenic belt of the eastern Nanling Range, South China, spatially associated with the Maoping granite pluton, hosts total ore reserves of 103,000 t WO3 and 50,000 t Sn. Two different types of mineralization developed in this deposit: Upper quartz vein-type mineralization, mostly within the Cambrian metamorphosed sandstone and slate, and underneath greisen-type mineralization within the Maoping granite. Cassiterites from both types of mineralization coexist with wolframite. Here we report for the first time in situ U–Pb data on cassiterite and zircon of the Maoping deposit obtained by LA-ICP-MS. Cassiterite from quartz vein and greisen yielded weighted average 206Pb/238U ages of 156.8 ± 1.5 Ma and 156.9 ± 1.4 Ma, respectively, which indicates that the two types of mineralization formed roughly at the same time. In addition, the two mineralization ages are consistent with the emplacement age of the Maoping granite (159.0 ± 1.5 Ma) within error, suggesting a close temporal and genetic link between W–Sn mineralization and granitic magmatism. The two types of mineralization formed at the same magmatic-hydrothermal event. Cassiterite from both types of mineralization shows high Fe, Ta, and Zr contents with a low Zr/Hf ratio, suggesting that the ore-forming fluid should be derived from the highly differentiated Maoping granite pluton. Cassiterite in greisen has higher contents of Nb and Ta but a lower concentration of Ti compared with that in quartz vein, indicating that the formation temperature of greisen-type mineralization is little higher than that of quartz-vein-type mineralization.

Description: The new Chinese Ka-band solid-state transmitter cloud radar (CR) uses four operational modes with different pulse widths and coherent integration and non-coherent integration numbers to meet long-term cloud measurement requirements. The CR and an instrument-equipped aircraft were used to observe clouds and precipitation on the east side of Taihang Mountain in Hebei Province in 2018. To resolve the data quality problems caused by attenuation in the precipitation area; we focused on developing an algorithm for attenuation correction based on rain drop size distribution (DSD) retrieved from the merged Doppler spectral density data of the four operational modes following data quality control (QC). After dealiasing Doppler velocity and removal of range sidelobe artifacts; we merged the four types of Doppler spectral density data. Vertical air speed and DSD are retrieved from the merged Doppler spectral density data. Finally, we conducted attenuation correction of Doppler spectral density data and recalculated Doppler moments such as reflectivity; radial velocity; and spectral width. We evaluated the consistencies of reflectivity spectra from the four operational modes and DSD retrieval performance using airborne in situ observation. We drew three conclusions: First, the four operational modes observed similar reflectivity and velocity for clouds and low-velocity solid hydrometeors; however; three times of coherent integration underestimated Doppler reflectivity spectra for velocities greater than 2 m s−1. Reflectivity spectra were also underestimated for low signal-to-noise ratios in the low-sensitivity operational mode. Second, QC successfully dealiased Doppler velocity and removed range sidelobe artifacts; and merging of the reflectivity spectra mitigated the effects of coherent integration and pulse compression on radar data. Lastly, the CR observed similar DSD and liquid water content vertical profiles to airborne in situ observations. Comparing CR and aircraft data yielded uncertainty due to differences in observation space and temporal and spatial resolutions of the data.

Description: The Postmasburg Manganese Field (PMF), Northern Cape Province, South Africa, once represented one of the largest sources of manganese ore worldwide. Two belts of manganese ore deposits have been distinguished in the PMF, namely the Western Belt of ferruginous manganese ores and the Eastern Belt of siliceous manganese ores. Prevailing models of ore formation in these two belts invoke karstification of manganese-rich dolomites and residual accumulation of manganese wad which later underwent diagenetic and low-grade metamorphic processes. For the most part, the role of hydrothermal processes and metasomatic alteration towards ore formation has not been adequately discussed. Here we report an abundance of common and some rare Al-, Na-, K- and Ba-bearing minerals, particularly aegirine, albite, microcline, banalsite, sérandite-pectolite, paragonite and natrolite in Mn ores of the PMF, indicative of hydrothermal influence. Enrichments in Na, K and/or Ba in the ores are generally on a percentage level for most samples analysed through bulk-rock techniques. The presence of As-rich tokyoite also suggests the presence of As and V in the hydrothermal fluid. The fluid was likely oxidized and alkaline in nature, akin to a mature basinal brine. Various replacement textures, particularly of Na- and K- rich minerals by Ba-bearing phases, suggest sequential deposition of gangue as well as ore-minerals from the hydrothermal fluid, with Ba phases being deposited at a later stage. The stratigraphic variability of the studied ores and their deviation from the strict classification of ferruginous and siliceous ores in the literature, suggests that a re-evaluation of genetic models is warranted. New Ar-Ar ages for K-feldspars suggest a late Neoproterozoic timing for hydrothermal activity. This corroborates previous geochronological evidence for regional hydrothermal activity that affected Mn ores at the PMF but also, possibly, the high-grade Mn ores of the Kalahari Manganese Field to the north. A revised, all-encompassing model for the development of the manganese deposits of the PMF is then proposed, whereby the source of metals is attributed to underlying carbonate rocks beyond the Reivilo Formation of the Campbellrand Subgroup. The main process by which metals are primarily accumulated is attributed to karstification of the dolomitic substrate. The overlying Asbestos Hills Subgroup banded iron formation (BIF) is suggested as a potential source of alkali metals, which also provides a mechanism for leaching of these BIFs to form high-grade residual iron ore deposits.

Description: Mixed cyanobacteria-dominated biofilms, enriched from a tributary of the Mérantaise (France) were used to conduct laboratory experiments in order to understand the relationship between the morphology of carbonate precipitates and the biological activity (e.g., cyanobacterial exopolymeric substances (EPS) production, photosynthetic pH increases). DNA sequencing data showed that the enriched biofilm was composed predominantly of two types of filamentous cyanobacteria that belonged to the Oscillatoriaceae and Phormidiaceae families, respectively. Microscopic analysis also indicated the presence of some coccoid cyanobacteria resembling Gloeocapsa. Analysis of carbonate precipitates in experimental biofilms showed three main morphologies: micro-peloids with different shapes of mesocrystals associated with Oscillatoriaceae filaments and theirs EPS, lamellae of carbonate formed directly on Phormidiaceae filaments, and rhombic sparite crystals wrapped in EPS. All crystals were identified by FT-IR spectroscopy as calcite. Similar structures as those that formed in laboratory conditions were observed in the microbial-tufa deposits collected in the stream. Microscopic and spectroscopic analysis of laboratory and natural samples indicated a close proximity of the cyanobacterial EPS and precipitated carbonates in both. Based on the laboratory experiments, we conclude that the microbial tufa in the stream is in an early stage of formation.

Description: We introduce the fully automatic design of a numerically optimized decision-tree algorithm and demonstrate its application to sea ice classification from SAR data. In the decision tree, an initial multi-class classification problem is split up into a sequence of binary problems. Each branch of the tree separates one single class from all other remaining classes, using a class-specific selected feature set. We optimize the order of classification steps and the feature sets by combining classification accuracy and sequential search algorithms, looping over all remaining features in each branch. The proposed strategy can be adapted to different types of classifiers and measures for the class separability. In this study, we use a Bayesian classifier with non-parametric kernel density estimation of the probability density functions. We test our algorithm on simulated data as well as airborne and spaceborne SAR data over sea ice. For the simulated cases, average per-class classification accuracy is improved between 0.5% and 4% compared to traditional all-at-once classification. Classification accuracy for the airborne and spaceborne SAR datasets was improved by 2.5% and 1%, respectively. In all cases, individual classes can show larger improvements up to 8%. Furthermore, the selection of individual feature sets for each single class can provide additional insights into physical interpretation of different features. The improvement in classification results comes at the cost of longer computation time, in particular during the design and training stage. The final choice of the optimal algorithm therefore depends on time constraints and application purpose.

Description: As a state-of-the-art method, the digital image correlation (DIC) technique is used to capture the fracture properties of wood along the longitudinal direction, such as the crack propagation, the strain field, and the fracture process zone (FPZ). Single-edge notched (SEN) specimens made of Douglas fir (Pseudotsuga menziesii) from Canada with different notch-to-depth ratios are tested by three-point-bending (3-p-b) experiment. The crack mouth opening displacements (CMOD) measured by the clip gauge and DIC technique agree well with each other, verifying the applicability of the DIC technique. Then, the quasi-brittle fracture process of wood is analyzed by combing the load-CMOD curve and the strain field in front of the preformed crack. Additionally, the equivalent elastic crack length is calculated using the linear superposition hypothesis. The comparison between the FPZ evolution and the equivalent elastic crack shows that specimens with higher notch-to-depth ratios have better cohesive effect and higher cracking resistance.

Description: The simultaneous availability of observations from space by remote sensing platforms operating at multiple frequencies in the microwave domain suggests investigating their complementarity in thematic mapping and retrieval of biophysical parameters. In particular, there is an interest to understand whether the wealth of short wavelength Synthetic Aperture Radar (SAR) backscatter observations at X-, C-, and L-band from currently operating spaceborne missions can improve the retrieval of forest stem volume, i.e., above-ground biomass, in the boreal zone with respect to a single frequency band. To this scope, repeated observations from TerraSAR-X, Sentinel-1 and ALOS-2 PALSAR-2 from the test sites of Remningstorp and Krycklan, Sweden, have been analyzed and used to estimate stem volume with a retrieval framework based on the Water Cloud Model. Individual estimates of stem volume were then combined linearly to form single-frequency and multi-frequency estimates. The retrieval was assessed at large 0.5 ha forest inventory plots (Remningstorp) and small 0.03 ha forest inventory plots (Krycklan). The relationship between SAR backscatter and stem volume differed depending on forest structure and environmental conditions, in particular at X- and C-band. The highest retrieval accuracy was obtained at both test sites at L-band. The combination of stem volume estimates from data acquired at two or three frequencies achieved an accuracy that was superior to values obtained at a single frequency. When combining estimates from X-, C-, and L-band data, the relative RMSE for the 0.5 ha inventory plots at Remningstorp was 31.3%. For the 0.03 ha inventory plots at Krycklan, the relative RMSE was above 50%. In a retrieval scenario involving short wavelength SAR backscatter data, these results suggest combining multiple frequencies to ensure the highest possible retrieval accuracy achievable. Retrievals should be undertaken to target spatial scales well above the size of a pixel.

Description: Mixed pixels in medium spatial resolution imagery create major challenges in acquiring accurate pixel-based land use and land cover information. Deep belief network (DBN), which can provide joint probabilities in land use and land cover classification, may serve as an alternative tool to address this mixed pixel issue. Since DBN performs well in pixel-based classification and object-based identification, examining its performance in subpixel unmixing with medium spatial resolution multispectral image in urban environments would be of value. In this study, (1) we examined DBN’s ability in subpixel unmixing with Landsat imagery, (2) explored the best-fit parameter setting for the DBN model and (3) evaluated its performance by comparing DBN with random forest (RF), support vector machine (SVM) and multiple endmember spectral mixture analysis (MESMA). The results illustrated that (1) DBN performs well in subpixel unmixing with a mean absolute error (MAE) of 0.06 and a root mean square error (RMSE) of 0.0077. (2) A larger sample size (e.g., greater than 3000) can provide stable and high accuracy while two-RBM-layer and 50 batch sizes are the best parameters for DBN in this study. Epoch size and learning rate should be decided by specific applications since there is not a consistent pattern in our experiments. Finally, (3) DBN can provide comparable results compared to RF, SVM and MESMA. We concluded that DBN can be viewed as an alternative method for subpixel unmixing with Landsat imagery and this study provides references for other scholars to use DBN in subpixel unmixing in urban environments.

Description: As the world urbanizes and builds more infrastructure, the extraction of built-up areas using remote sensing is crucial for monitoring land cover changes and understanding urban environments. Previous studies have proposed a variety of methods for mapping regional and global built-up areas. However, most of these methods rely on manual selection of training samples and classification thresholds, leading to low extraction efficiency. Furthermore, thematic accuracy is limited by interference from other land cover types like bare land, which hinder accurate and timely extraction and monitoring of dynamic changes in built-up areas. This study proposes a new method to map built-up areas by combining VIIRS (Visible Infrared Imaging Radiometer Suite) nighttime lights (NTL) data and Landsat-8 multispectral imagery. First, an adaptive NTL threshold was established, vegetation and water masks were superimposed, and built-up training samples were automatically acquired. Second, the training samples were employed to perform supervised classification of Landsat-8 data before deriving the preliminary built-up areas. Third, VIIRS NTL data were used to obtain the built-up target areas, which were superimposed onto the built-up preliminary classification results to obtain the built-up area fine classification results. Four major metropolitan areas in Eurasia formed the study areas, and the high spatial resolution (20 m) built-up area product High Resolution Layer Imperviousness Degree (HRL IMD) 2015 served as the reference data. The results indicate that our method can accurately and automatically acquire built-up training samples and adaptive thresholds, allowing for accurate estimates of the spatial distribution of built-up areas. With an overall accuracy exceeding 94.7%, our method exceeded accuracy levels of the FROM-GLC and GUL built-up area products and the PII built-up index. The accuracy and efficiency of our proposed method have significant potential for global built-up area mapping and dynamic change monitoring.

Description: We present a joint 2D inversion approach for magnetotelluric (MT) and gravity data with elastic-net regularization and cross-gradient constraints. We describe the main features of the approach and verify the inversion results against a synthetic model. The results indicate that the best fit solution using the L2 is overly smooth, while the best fit solution for the L1 norm is too sparse. However, the elastic-net regularization method, a convex combination term of L2 norm and L1 norm, can not only enforce the stability to preserve local smoothness, but can also enforce the sparsity to preserve sharp boundaries. Cross-gradient constraints lead to models with close structural resemblance and improve the estimates of the resistivity and density of the synthetic dataset. We apply the novel approach to field datasets from a copper mining area in the northeast of China. Our results show that the method can generate much more detail and a sharper boundary as well as better depth resolution. Relative to the existing solution, the large area divergence phenomenon under the anomalous bodies is eliminated, and the fine anomalous bodies boundary appeared in the smooth region. This method can provide important technical support for detecting deep concealed deposits.

Description: Reliable image classification results are crucial for the application of remote sensing images, but the reliability of image classification has received less attention. In particular, the inherent uncertainty of remote sensing images has been disregarded. The uncertainty of a remote sensing image accumulates and propagates continuously in the classification process and ultimately affects the reliability of the classification results. Therefore, quantitative description and investigation of the inherent uncertainty of remote sensing images are crucial in achieving reliable remote sensing image classification. In this study, we analyze the sources of uncertainty of remote sensing images in detail and propose a quantitative descriptor for measuring image uncertainty comprehensively and effectively. In addition, we also design two verification schemes to verify the validity of the proposed uncertainty descriptor. Finally, the validity of the proposed uncertainty descriptor is confirmed by experimental results on three real remote sensing images. Our study on the uncertainty of remote sensing images may help the development of uncertainty control methods and reliable classification schemes of remote sensing images.

Description: Current PM2.5 retrieval maps have many missing values, which seriously hinders their performance in real applications. This paper presents a framework to map full-coverage daily average PM2.5 concentrations from MODIS C6 aerosol optical depth (AOD) products and fill missing pixels in both the AOD and PM2.5 maps. First, a two-stage inversed variance weights (IVW) algorithm was adopted to fuse the MODIS C6 Terra and Aqua AOD products, which fills missing data in MODIS standard AOD data and obtains a high coverage daily average. After that, using the fused MODIS daily average AOD and ground-level PM2.5 in all grid cells, a two-stage generalized additive model (GAM) was implemented to obtain the full-coverage PM2.5 concentrations. Experiments on the Yangtze River Delta (YRD) in 2013–2016 were carefully designed to validate the performance of our proposed framework. The results show that the two-stage IVW could not only improve the spatial coverage of MODIS AOD against the original standard product by 230%, but could also keep its data accuracy. When compared with the ground-level measurements, the two-stage GAM can obtain accurate PM2.5 concentration estimates (R2 = 0.78, RMSE = 19.177 μg/m3, and RPE = 28.9%). Moreover, our method performs better than the inverse distance weighted method and kriging methods in mapping full-coverage daily PM2.5 concentrations. Therefore, the proposed framework provides a good methodology for retrieving full-coverage daily average PM2.5 concentrations from MODIS standard AOD products.

Description: The Chinese BeiDou Navigation Satellite System has provided a global-coverage service since 27 December 2018. Eighteen BD3 MEO satellites have been launched into space during 2017 and 2018. The signal constitution has been redesigned and four open service signals are used for transmission, including B1I, B1C, B2a and B3I. This paper focuses on the signal performance, Precise Orbit Determination (POD) and the atomic clock’s frequency stability issues of the BD3 satellites. The satellite-induced code bias issue found in BD2 satellites multipath combination has been proven to be eliminated in BD3 satellites. However, the pseudorange code of B1C is much noisier than that of other three frequencies, which may be related to the signal constitution and power distribution, as the minimum received power levels on the ground of B1C is 3 dB lower than that of the B2a signal. Similar results were achieved by the Ionosphere-Free combination residuals in POD using four signals, B1I-B3I, B1I-B2a, B1C-B3I and B1C-B2a, and the phase residual of B1C-B2a combination performed best. Considering the noise amplitude and compatibility with other GNSS (Global Navigation Satellite System), the B1C-B2a combination is recommended in priority for precise GNSS data processing. GFIFP combinations were also implemented to evaluate the inter-frequency phase bias of the four signals. The experimental results showed that the systematic signal with an amplitude of about 2 cm could be found in the GFIFP series. In addition, multi-GNSS POD was performed and analyzed as well, using about a hundred global-distributed IGS and iGMAS stations. Furthermore, the atomic clock’s frequency stability was estimated using the parameters of clock bias calculated in POD and the Overlap Allan Deviations showed that the frequency stability of BD3 reached approximately 2.43 × 10−14 at intervals of 10,000 s and 2.51 × 10−15 at intervals of 86,400 s, which was better than that of the GPS BLOCK IIF satellites but worse than that of Galileo satellites.

Description: Removal of calcium and magnesium ions through biomineralization induced by bacteria has been proven to be an effective and environmentally friendly method to improve water quality, but the process and mechanism are far from fully understood. In this study, a newly isolated probiotic Bacillus licheniformis SRB2 (GenBank: KM884945.1) was used to induce the bio-precipitation of calcium and magnesium at various Mg/Ca molar ratios (0, 6, 8, 10, and 12) in medium with 30 g L−1 sodium chloride. Due to the increasing pH and HCO3− and CO32− concentrations caused by NH3 and carbonic anhydrase, about 98% Ca2+ and 50% Mg2+ were precipitated in 12 days. The pathways of bio-precipitation include extracellular and intracellular processes. Biominerals with more negative δ13C values (−16‰ to −18‰) were formed including calcite, vaterite, monohydrocalcite, and nesquehonite with preferred orientation. The nucleation on extracellular polymeric substances was controlled by the negatively charged amino acids and organic functional groups. The intracellular amorphous inclusions containing calcium and magnesium also contributed to the bio-precipitation. This study reveals the process and mechanism of microbial desalination for the removal of calcium and magnesium, and provides some references to explain the formation of the nesquehonite and other carbonate minerals in a natural and ancient earth surface environment.

Description: Bangladesh is one of the most flood-affected countries in the world. In the last few decades, flood frequency, intensity, duration, and devastation have increased in Bangladesh. Identifying flood-damaged areas is highly essential for an effective flood response. This study aimed at developing an operational methodology for rapid flood inundation and potential flood damaged area mapping to support a quick and effective event response. Sentinel-1 images from March, April, June, and August 2017 were used to generate inundation extents of the corresponding months. The 2017 pre-flood land cover maps were prepared using Landsat-8 images to identify major land cover on the ground before flooding. The overall accuracy of flood inundation mapping was 96.44% and the accuracy of the land cover map was 87.51%. The total flood inundated area corresponded to 2.01%, 4.53%, and 7.01% for the months April, June, and August 2017, respectively. Based on the Landsat-8 derived land cover information, the study determined that cropland damaged by floods was 1.51% in April, 3.46% in June, 5.30% in August, located mostly in the Sylhet and Rangpur divisions. Finally, flood inundation maps were distributed to the broader user community to aid in hazard response. The data and methodology of the study can be replicated for every year to map flooding in Bangladesh.

Description: Despite recent research on the potential of dual- (DP) and full-polarimetry (FP) Synthetic Aperture Radar (SAR) data for crop mapping, the capability of compact polarimetry (CP) SAR data has not yet been thoroughly investigated. This is of particular concern, given the availability of such data from RADARSAT Constellation Mission (RCM) shortly. Previous studies have illustrated potential for accurate crop mapping using DP and FP SAR features, yet what contribution each feature makes to the model accuracy is not well investigated. Accordingly, this study examined the potential of the early- to mid-season (i.e., May to July) RADARSAT-2 SAR images for crop mapping in an agricultural region in Manitoba, Canada. Various classification scenarios were defined based on the extracted features from FP SAR data, as well as simulated DP and CP SAR data at two different noise floors. Both overall and individual class accuracies were compared for multi-temporal, multi-polarization SAR data using the pixel- and object-based random forest (RF) classification schemes. The late July C-band SAR observation was the most useful data for crop mapping, but the accuracy of single-date image classification was insufficient. Polarimetric decomposition features extracted from CP and FP SAR data produced relatively equal or slightly better classification accuracies compared to the SAR backscattering intensity features. The RF variable importance analysis revealed features that were sensitive to depolarization due to the volume scattering are the most important FP and CP SAR data. Synergistic use of all features resulted in a marginal improvement in overall classification accuracies, given that several extracted features were highly correlated. A reduction of highly correlated features based on integrating the Spearman correlation coefficient and the RF variable importance analyses boosted the accuracy of crop classification. In particular, overall accuracies of 88.23%, 82.12%, and 77.35% were achieved using the optimized features of FP, CP, and DP SAR data, respectively, using the object-based RF algorithm.

Description: Coffee grounds are the most significant production waste in the coffee industry and contain about 15% coffee oil. Coffee oil is rich in fatty acids and polyphenols, which have great application potential in the flotation of oxidized minerals. In this study, coffee oil as a green flotation collector for ilmenite was investigated by micro-flotation, zeta potential measurement, and foam stability analysis. The results of zeta potential reveal that both coffee oil and MOH can be adsorbed on the ilmenite surface at pH 6.7, and the chemical adsorption mode is dominant. However, when the pH is 2.8, the adsorption capacity of coffee oil on the ilmenite surface is much larger than that of MOH. The pH value of the pulp has little effect on the foam properties in the coffee oil solution and has a great influence on the foaming performance and foam stability of the MOH solution. When coffee oil is used as a collector, the grade of TiO2 in ilmenite concentrate is increased from 21.68% to 46.83%, and the recovery is 90.22%, indicating that the potential of coffee oil in the application of ilmenite flotation is large.

Description: The controlled crystallisation of struvite (MgNH4PO4∙6H2O) is a viable means for the recovery and recycling of phosphorus (P) from municipal and industrial wastewaters. However, an efficient implementation of this recovery method in water treatment systems requires a fundamental understanding of struvite crystallisation mechanisms, including the behavior and effect of metal contaminants during struvite precipitation. Here, we studied the crystallisation pathways of struvite from aqueous solutions using a combination of ex situ and in situ time-resolved synthesis and characterization techniques, including synchrotron-based small- and wide-angle X-ray scattering (SAXS/WAXS) and cryogenic transmission electron microscopy (cryo-TEM). Struvite syntheses were performed both in the pure Mg-NH4-PO4 system as well as in the presence of cobalt (Co), which, among other metals, is typically present in waste streams targeted for P-recovery. Our results show that in the pure system and at Co concentrations 〈 0.5 mM, struvite crystals nucleate and grow directly from solution, much in accordance with the classical notion of crystal formation. In contrast, at Co concentrations ≥ 1 mM, crystallisation was preceded by the transient formation of an amorphous nanoparticulate phosphate phase. Depending on the aqueous Co/P ratio, this amorphous precursor was found to transform into either (i) Co-bearing struvite (at Co/P 〈 0.3) or (ii) cobalt phosphate octahydrate (at Co/P 〉 0.3). These amorphous-to-crystalline transformations were accompanied by a marked colour change from blue to pink, indicating a change in Co2+ coordination in the formed solid from tetrahedral to octahedral. Our findings have implications for the recovery of nutrients and metals during struvite crystallisation and contribute to the ongoing general discussion about the mechanisms of crystal formation.

Description: This work aims to provide daily estimates of the evolution of popcorn dry masses at the field scale using an agro-meteorological model, named the simple algorithm for yield model combined with a water balance model (SAFY-WB), controlled by the Green Area Index (GAI), derived from satellite images acquired in the microwave and optical domains. Synthetic aperture radar (SAR) satellite information (σ°VH/VV) was provided by the Sentinel-1A (S1-A) mission through two orbits (30 and 132), with a repetitiveness of six days. The optical data were obtained from the Landsat-8 mission. SAR and optical data were acquired over one complete agricultural season, in 2016, over a test site located in the southwest of France. The results show that the total dry masses of corn can be estimated accurately (R² = 0.92) at daily time steps due to a combination of satellite and model data. The SAR data are more suitable for characterizing the first period of crop development (until the end of flowering), whereas the optical data can be used throughout the crop cycle. Moreover, the model offers good performances in plant (R2 = 0.90) and ear (R2 = 0.93) mass retrieval, irrespective of the phenological stage. The results also reveal that four phenological stages (four to five leaves, flowering, ripening, and harvest) can be accurately predicted by the proposed approach (R2 = 0.98; root-mean-square error (RMSE) = seven days). Nevertheless, some important points must be taken into account before assimilation, namely the SAR signal must be corrected with respect to thermal noise before being assimilated, and the relationship estimated between the GAI and SAR signal must be performed over fields cultivated without intercrops. These results are unique in the literature and provide a new way to better monitor corn production over time.

Description: Surface and deep potential geophysical signals respond to the spatial redistribution of global mass variations, which may be monitored by geodetic observations. In this study, we analyze dense Global Positioning System (GPS) time series in the Eastern Tibetan Plateau using principal component analysis (PCA) and wavelet time-frequency spectra. The oscillations of interannual and residual signals are clearly identified in the common mode component (CMC) decomposed from the dense GPS time series from 2000 to 2018. The newly developed spherical harmonic coefficients of the Gravity Recovery and Climate Experiment Release-06 (GRACE RL06) are adopted to estimate the seasonal and interannual patterns in this region, revealing hydrologic and atmospheric/nontidal ocean loads. We stack the averaged elastic GRACE-derived loading displacements to identify the potential physical significance of the CMC in the GPS time series. Interannual nonlinear signals with a period of ~3 to ~4 years in the CMC (the scaled principal components from PC1 to PC3) are found to be predominantly related to hydrologic loading displacements, which respond to signals (El Niño/La Niña) of global climate change. We find an obvious signal with a period of ~6 yr on the vertical component that could be caused by mantle-inner core gravity coupling. Moreover, we evaluate the CMC’s effect on the GPS-derived velocities and confirm that removing the CMC can improve the recognition of nontectonic crustal deformation, especially on the vertical component. Furthermore, the effects of the CMC on the three-dimensional velocity and uncertainty are presented to reveal the significant crustal deformation and dynamic processes of the Eastern Tibetan Plateau.

Description: The accurate prediction of surface solar irradiance is of great significance for the generation of photovoltaic power. Surface solar irradiance is affected by many random mutation factors, which means that there are great challenges faced in short-term prediction. In Northwest China, there are abundant solar energy resources and large desert areas, which have broad prospects for the development of photovoltaic (PV) systems. For the desert areas in Northwest China, where meteorological stations are scarce, satellite remote sensing data are extremely precious exploration data. In this paper, we present a model using FY-4A satellite images to forecast (up to 15–180 min ahead) global horizontal solar irradiance (GHI), at a 15 min temporal resolution in desert areas under different sky conditions, and compare it with the persistence model (SP). The spatial resolution of the FY-4A satellite images we used was 1 km × 1 km. Particle image velocimetry (PIV) was used to derive the cloud motion vector (CMV) field from the satellite cloud images. The accuracy of the forecast model was evaluated by the ground observed GHI data. The results showed that the normalized root mean square error (nRMSE) ranged from 18.9% to 21.6% and the normalized mean bias error (nMBE) ranged from 3.2% to 4.9% for time horizons from 15 to 180 min under all sky conditions. Compared with the SP model, the nRMSE value was reduced by about 6%, 8%, and 14% with the time horizons of 60, 120, and 180 min, respectively.

Description: Traditional ways of validating satellite-derived sea surface temperature (SST) and sea surface salinity (SSS) products by comparing with buoy measurements, do not allow for evaluating the impact of mesoscale-to-submesoscale variability. We present the validation of remotely sensed SST and SSS data against the unmanned surface vehicle (USV)—called Saildrone—measurements from the 60 day 2018 Baja California campaign. More specifically, biases and root mean square differences (RMSDs) were calculated between USV-derived SST and SSS values, and six satellite-derived SST (MUR, OSTIA, CMC, K10, REMSS, and DMI) and three SSS (JPLSMAP, RSS40, RSS70) products. Biases between the USV SST and OSTIA/CMC/DMI were approximately zero, while MUR showed a bias of 0.3 °C. The OSTIA showed the smallest RMSD of 0.39 °C, while DMI had the largest RMSD of 0.5 °C. An RMSD of 0.4 °C between Saildrone SST and the satellite-derived products could be explained by the diurnal and sub-daily variability in USV SST, which currently cannot be resolved by remote sensing measurements. SSS showed fresh biases of 0.1 PSU for JPLSMAP and 0.2 PSU and 0.3 PSU for RMSS40 and RSS70 respectively. SST and SSS showed peaks in coherence at 100 km, most likely associated with the variability of the California Current System.

Description: The paper describes a new tool called JupyTEP integrated development environment (IDE), which is an online integrated development environment for earth observation data processing available in the cloud. This work is a result of the project entitled “JupyTEP IDE—Jupyter-based IDE as an interactive and collaborative environment for the development of notebook style EO algorithms on network of exploitation platforms infrastructure” carried out in cooperation with European Space Agency. The main goal of this project was to provide a universal earth observation data processing tool to the community. JupyTEP IDE is an extension of Jupyter software ecosystem with customization of existing components for the needs of earth observation scientists and other professional and non-professional users. The approach is based on configuration, customization, adaptation, and extension of Jupyter, Jupyter Hub, and Docker components on earth observation data cloud infrastructure in the most flexible way; integration with accessible libraries and earth observation data tools (sentinel application platform (SNAP), geospatial data abstraction library (GDAL), etc.); adaptation of existing web processing service (WPS)-oriented earth observation services. The user-oriented product is based on a web-related user interface in the form of extended and modified Jupyter user interface (frontend) with customized layout, earth observation data processing extension, and a set of predefined notebooks, widgets, and tools. The final IDE is addressed to the remote sensing experts and other users who intend to develop Jupyter notebooks with the reuse of embedded tools, common WPS interfaces, and existing notebooks. The paper describes the background of the system, its architecture, and possible use cases.

Description: Hyperspectral and light detection and ranging (LiDAR) data fusion and classification has been an active research topic, and intensive studies have been made based on mathematical morphology. However, matrix-based concatenation of morphological features may not be so distinctive, compact, and optimal for classification. In this work, we propose a novel Coupled Higher-Order Tensor Factorization (CHOTF) model for hyperspectral and LiDAR data classification. The innovative contributions of our work are that we model different features as multiple third-order tensors, and we formulate a CHOTF model to jointly factorize those tensors. Firstly, third-order tensors are built based on spectral-spatial features extracted via attribute profiles (APs). Secondly, the CHOTF model is defined to jointly factorize the multiple higher-order tensors. Then, the latent features are generated by mode-n tensor-matrix product based on the shared and unshared factors. Lastly, classification is conducted by using sparse multinomial logistic regression (SMLR). Experimental results, conducted with two popular hyperspectral and LiDAR data sets collected over the University of Houston and the city of Trento, respectively, indicate that the proposed framework outperforms the other methods, i.e., different dimensionality-reduction-based methods, independent third-order tensor factorization based methods, and some recently proposed hyperspectral and LiDAR data fusion and classification methods.

Description: Floods frequently occur in Nigeria. The catastrophic 2012 flood in Nigeria claimed 363 lives and affected about seven million people. A total loss of about 2.29 trillion Naira (7.2 billion US Dollars) was estimated. The effect of flooding in the country has been devastating because of sparse to no flood monitoring, and a lack of an effective early flood warning system in the country. Here, we evaluated the efficacy of using the Gravity Recovery and Climate Experiment (GRACE) terrestrial water storage anomaly (TWSA) to evaluate the hydrological conditions of the Lower Niger River Basin (LNRB) in Nigeria in terms of precipitation and antecedent terrestrial water storage prior to the 2012 flood event. Furthermore, we accessed the potential of the GRACE-based flood potential index (FPI) at correctly predicting previous floods, especially the devastating 2012 flood event. For validation, we compared the GRACE terrestrial water storage capacity (TWSC) quantitatively and qualitatively to the water budget of TWSC and Dartmouth Flood Observatory (DFO) respectively. Furthermore, we derived a water budget-based FPI using Reager’s methodology and compared it to the GRACE-derived FPI quantitatively. Generally, the GRACE TWSC estimates showed seasonal consistency with the water budget TWSC estimates with a correlation coefficient of 0.8. The comparison between the GRACE-derived FPI and water budget-derived FPI gave a correlation coefficient of 0.9 and also agreed well with the flood reported by the DFO. Also, the FPI showed a marked increase with precipitation which implies that rainfall is the main cause of flooding in the study area. Additionally, the computed GRACE-based storage deficit revealed that there was a decrease in water storage prior to the flooding month while the FPI increased. Hence, the GRACE-based FPI and storage deficit when supplemented with water budget-based FPI could suggest a potential for flood prediction and water storage monitoring respectively.

Description: Hyperspectral imagery contains abundant spectral information. Each band contains some specific characteristics closely related to target objects. Therefore, using these characteristics, hyperspectral imagery can be used for anomaly detection. Recently, with the development of compressed sensing, low-rank-representation-based methods have been applied to hyperspectral anomaly detection. In this study, novel low-rank representation methods were developed for anomaly detection from hyperspectral images based on the assumption that hyperspectral pixels can be effectively decomposed into a low-rank component (for background) and a sparse component (for anomalies). In order to improve detection performance, we imposed a spatial constraint on the low-rank representation coefficients, and single or multiple local window strategies was applied to smooth the coefficients. Experiments on both simulated and real hyperspectral datasets demonstrated that the proposed approaches can effectively improve hyperspectral anomaly detection performance.

Description: Bauxites in southern France (Provence and Languedoc) have been exploited since the beginning of the last century. Though most of the deposits are now subeconomic or mined-out, these bauxites represent model analogs for other economic bauxites of the world. These Cretaceous karst-type deposits lie directly on Jurassic carbonates, and have been formed through a combination of different processes: in-situ alteration of siliciclastic sediments deposited on carbonate platforms, and reworking of early bauxites in the karst network. In this study, we present preliminary bulk rock geochemical and in-situ laser ablation (LA) -ICP-MS analyses on Al- and Fe-oxy-hydroxides of Provence (Les Baux-de-Provence) and Languedoc (Villeveyrac, Loupian) bauxites, with the aim of evaluating the concentrations of rare earth elements (REEs) and their deportment in these minerals. REEs have total average concentrations of 700 mg/kg in the analyzed samples, which are mostly composed of boehmite, γ-AlO(OH), and Fe-oxy-hydroxides (hematite and goethite). Maximum REEs concentrations are commonly associated with positive Ce anomalies in chondrite-normalized patterns. In contrast with other examples from the literature, it has been observed that high REE concentrations also occur in samples apparently devoid or poor of REE-minerals. In these samples, the total amount of REEs is positively correlated with that of Ga (commonly contained in boehmite). LA-ICP-MS trace element analyses on boehmite and Fe-oxy-hydroxides have shown that while the Al-hydroxide contains the suite of REEs, goethite and hematite are preferentially enriched only in Ce. Considering that Al-hydroxides are digested during the Bayer process, an interesting issue to develop in the future is whether (and how) REEs released during Al-hydroxide digestion could be recovered together with Al from the pregnant leach liquor, as routinely done for Ga.

Description: River deltas and estuaries belong to the most significant coastal landforms on our planet and are usually very densely populated. Nearly 600 million people live in river deltas, benefiting from the large variety of locational advantages and rich resources. Deltas are highly dynamic and vulnerable environments that are exposed to a wide range of natural and man-made threats. Sustainable management of river deltas therefore requires a holistic assessment of historic and recent ongoing changes and the dynamics in settlement sprawl, land cover and land use change, ecosystem development, as well as river and coastline geomorphology, all of which is difficult to achieve solely with traditional land-based surveying techniques. This review paper presents the potential of Earth Observation for analyses and quantification of land surface dynamics in the large river deltas globally, emphasizing the different geo-information products that can be derived from medium resolution, high resolution and highest resolution optical, multispectral, thermal and SAR data. Over 200 journal papers on remote sensing related studies for large river deltas and estuaries have been analyzed and categorized into thematic fields such as river course morphology, coastline changes, erosion and accretion processes, flood and inundation dynamics, regional land cover and land use dynamics, as well as the monitoring of compliance with respect to anthropogenic activity such as industry expansion-related habitat destruction. Additionally, our own exemplary analyses are interwoven into the review to visualize related delta work.

Description: Migration is a valuable life history strategy for many species because it enables individuals to exploit spatially and temporally variable resources. Globally, the prevalence of species’ migratory behavior is decreasing as individuals forgo migration to remain resident year-round, an effect hypothesized to result from anthropogenic changes to landscape dynamics. Efforts to conserve and restore migrations require an understanding of the ecological characteristics driving the behavioral tradeoff between migration and residence. We identified migratory and resident behaviors of 42 mule deer (Odocoileus hemionus) based on GPS locations and correlated their locations to remotely sensed indicators of forage quality, land cover, snow cover, and human land use. The model classified mule deer seasonal migratory and resident niches with an overall accuracy of 97.8% and cross-validated accuracy of 81.2%. The distance to development was the most important variable in discriminating in which environments these behaviors occur, with resident niche space most often closer to developed areas than migratory niches. Additionally, snow cover in December was important for discriminating summer migratory niches. This approach demonstrates the utility of niche analysis based on remotely sensed environmental datasets and provides empirical evidence of human land use impacts on large-scale wildlife migrations.

Description: Timely and accurate information about spatial distribution of tree species in urban areas provides crucial data for sustainable urban development, management and planning. Very high spatial resolution data collected by sensors onboard Unmanned Aerial Vehicles (UAV) systems provide rich data sources for mapping tree species. This paper proposes a method of tree species mapping from UAV images over urban areas using similarity in tree-crown object histograms and a simple thresholding method. Tree-crown objects are first extracted and used as processing units in subsequent steps. Tree-crown object histograms of multiple features, i.e., spectral and height related features, are generated to quantify within-object variability. A specific tree species is extracted by comparing similarity in histogram between a target tree-crown object and reference objects. The proposed method is evaluated in mapping four different tree species using UAV multispectral ortho-images and derived Digital Surface Model (DSM) data collected in Shanghai urban area, by comparing with an existing method. The results demonstrate that the proposed method outperforms the comparative method for all four tree species, with improvements of 0.61–5.81% in overall accuracy. The proposed method provides a simple and effective way of mapping tree species over urban area.

Description: Storms can cause significant damage to forest areas, affecting biodiversity and infrastructure and leading to economic loss. Thus, rapid detection and mapping of windthrows are crucially important for forest management. Recent advances in computer vision have led to highly-accurate image classification algorithms such as Convolutional Neural Network (CNN) architectures. In this study, we tested and implemented an algorithm based on CNNs in an ArcGIS environment for automatic detection and mapping of damaged areas. The algorithm was trained and tested on a forest area in Bavaria, Germany. . It is a based on a modified U-Net architecture that was optimized for the pixelwise classification of multispectral aerial remote sensing data. The neural network was trained on labeled damaged areas from after-storm aerial orthophotos of a ca. 109 k m 2 forest area with RGB and NIR bands and 0.2-m spatial resolution. Around 10 7 pixels of labeled data were used in the process. Once the network is trained, predictions on further datasets can be computed within seconds, depending on the size of the input raster and the computational power used. The overall accuracy on our test dataset was 92 % . During visual validation, labeling errors were found in the reference data that somewhat biased the results because the algorithm in some instance performed better than the human labeling procedure, while missing areas affected by shadows. Our results are very good in terms of precision, and the methods introduced in this paper have several additional advantages compared to traditional methods: CNNs automatically detect high- and low-level features in the data, leading to high classification accuracies, while only one after-storm image is needed in comparison to two images for approaches based on change detection. Furthermore, flight parameters do not affect the results in the same way as for approaches that require DSMs and DTMs as the classification is only based on the image data themselves, and errors occurring in the computation of DSMs and DTMs do not affect the results with respect to the z component. The integration into the ArcGIS Platform allows a streamlined workflow for forest management, as the results can be accessed by mobile devices in the field to allow for high-accuracy ground-truthing and additional mapping that can be synchronized back into the database. Our results and the provided automatic workflow highlight the potential of deep learning on high-resolution imagery and GIS for fast and efficient post-disaster damage assessment as a first step of disaster management.

Description: Different types of methods have been developed to retrieve vegetation attributes from remote sensing data, including conventional empirical regressions (i.e., linear regression (LR)), advanced empirical regressions (e.g., multivariable linear regression (MLR), partial least square regression (PLSR)), machine learning (e.g., random forest regression (RFR), decision tree regression (DTR)), and radiative transfer modelling (RTM, e.g., PROSAIL). Given that each algorithm has its own strengths and weaknesses, it is essential to compare them and evaluate their effectiveness. Previous studies have mainly used single-date multispectral imagery or ground-based hyperspectral reflectance data for evaluating the models, while multi-seasonal hyperspectral images have been rarely used. Extensive spectral and spatial information in hyperspectral images, as well as temporal variations of landscapes, potentially influence the model performance. In this research, LR, PLSR, RFR, and PROSAIL, representing different types of methods, were evaluated for estimating vegetation chlorophyll content from bi-seasonal hyperspectral images (i.e., a middle- and a late-growing season image, respectively). Results show that the PLSR and RFR generally performed better than LR and PROSAIL. RFR achieved the highest accuracy for both images. This research provides insights on the effectiveness of different models for estimating vegetation chlorophyll content using hyperspectral images, aiming to support future vegetation monitoring research.

Description: Globally, the number of dams increased dramatically during the 20th century. As a result, monitoring water levels and storage volume of dam-reservoirs has become essential in order to understand water resource availability amid changing climate and drought patterns. Recent advancements in remote sensing data show great potential for studies pertaining to long-term monitoring of reservoir water volume variations. In this study, we used freely available remote sensing products to assess volume variations for Lake Mead, Lake Powell and reservoirs in California between 1984 and 2015. Additionally, we provided insights on reservoir water volume fluctuations and hydrological drought patterns in the region. We based our volumetric estimations on the area–elevation hypsometry relationship, by combining water areas from the Global Surface Water (GSW) monthly water history (MWH) product with corresponding water surface median elevation values from three different digital elevation models (DEM) into a regression analysis. Using Lake Mead and Lake Powell as our validation reservoirs, we calculated a volumetric time series for the GSWMWH–DEMmedian elevation combinations that showed a strong linear ‘area (WA) – elevation (WH)’ (R2 〉 0.75) hypsometry. Based on ‘WA-WH’ linearity and correlation analysis between the estimated and in situ volumetric time series, the methodology was expanded to reservoirs in California. Our volumetric results detected four distinct periods of water volume declines: 1987–1992, 2000–2004, 2007–2009 and 2012–2015 for Lake Mead, Lake Powell and in 40 reservoirs in California. We also used multiscalar Standardized Precipitation Evapotranspiration Index (SPEI) for San Joaquin drainage in California to assess regional links between the drought indicators and reservoir volume fluctuations. We found highest correlations between reservoir volume variations and the SPEI at medium time scales (12–18–24–36 months). Our work demonstrates the potential of processed, open source remote sensing products for reservoir water volume variations and provides insights on usability of these variations in hydrological drought monitoring. Furthermore, the spatial coverage and long-term temporal availability of our data presents an opportunity to transfer these methods for volumetric analyses on a global scale.

Description: High-resolution optical remote sensing data can be utilized to investigate the human behavior and the activities of artificial targets, for example ship detection on the sea. Recently, the deep convolutional neural network (DCNN) in the field of deep learning is widely used in image processing, especially in target detection tasks. Therefore, a complete processing system called the broad area target search (BATS) is proposed based on DCNN in this paper, which contains data import, processing and storage steps. In this system, aiming at the problem of onshore false alarms, a method named as Mask-Faster R-CNN is proposed to differentiate the target and non-target areas by introducing a semantic segmentation sub network into the Faster R-CNN. In addition, we propose a DCNN framework named as Saliency-Faster R-CNN to deal with the problem of multi-scale ships detection, which solves the problem of missing detection caused by the inconsistency between large-scale targets and training samples. Based on these DCNN-based methods, the BATS system is tested to verify that our system can integrate different ship detection methods to effectively solve the problems that existed in the ship detection task. Furthermore, our system provides an interface for users, as a data-driven learning, to optimize the DCNN-based methods.

Description: Power corridor classification using LiDAR (light detection and ranging) point clouds is an important means for power line inspection. Many supervised classification methods have been used for classifying power corridor scenes, such as using random forest (RF) and JointBoost. However, these studies did not systematically analyze all the relevant factors that affect the classification, including the class distribution, feature selection, classifier type and neighborhood radius for classification feature extraction. In this study, we examine these factors using point clouds collected by an airborne laser scanning system (ALS). Random forest shows strong robustness to various pylon types. When classifying complex scenes, the gradient boosting decision tree (GBDT) shows good generalization. Synthetically, considering performance and efficiency, RF is very suitable for power corridor classification. This study shows that balanced learning leads to poor classification performance in the current scene. Data resampling for the original unbalanced dataset may not be necessary. The sensitivity analysis shows that the optimal neighborhood radius for feature extraction of different objects may be different. Scale invariance and automatic scale selection methods should be further studied. Finally, it is suggested that RF, original unbalanced class distribution, and complete feature set should be considered for power corridor classification in most cases.

Description: Scan planning of buildings under construction is a key issue for an efficient assessment of work progress. This work presents an automatic method aimed to determinate the optimal scan positions and the optimal route based on the use of Building Information Models (BIM) and considering data completeness as stopping criteria. The method is considered for a Terrestrial Laser Scanner mounted on a mobile robot following a stop & go procedure. The method starts by extracting floor plans from the BIM model according to the planned construction status, and including geometry and semantics of the building elements considered for construction control. The navigable space is defined from a binary map considering a security distance to building elements. After a grid-based and a triangulation-based distribution are implemented for generating scan position candidates, a visibility analysis is carried out to determine the optimal number and position of scans. The optimal route to visit all scan positions is addressed by using a probabilistic ant colony optimization algorithm. The method has been tested in simulated and real buildings under very dissimilar conditions and structural construction elements. The two approaches for generating scan position candidates are evaluated and results show the triangulation-based distribution as the more efficient approach in terms of processing and acquisition time, especially for large-scale buildings.

Description: The Qarhan Salt Lake has attracted increasing attention due to its significant national economic status and increased human activity, especially mining. Therefore, a sediment core collected from the confluence of the Golmud River to the Qarhan Salt Lake was chosen to investigate the concentrations, pollution levels, and ecological assessment of nine targeted elements (Al, As, Cd, Cr, Cu, Ni, P, Pb, and Zn). The excess 210Pb activities were calculated and a sedimentation rate of approximately 0.041 cm/y was estimated. Elements sources were identified, and the results show that Al, As, Cu, Ni, Pb, and Zn were mainly from natural sources, Cd and P were mainly from human input, and Cr appeared to have both sources. For Cd and P there was an increasing trend from 1987 and 1975, respectively, coinciding with the Chinese economic reform, Qarhan Salt Lake development and utilization, and also with the gross domestic product of Haixi State, Qinghai Province. Though the pollution and ecological assessment showed that there was nil to very low contamination and ecological risk, which is different from previous assumptions, the obviously increasing trend of Cd and P in the surface is still a concern. More attention should be paid to Cd and P in the further development of the Qarhan Salt Lake and the Golmud City.

Description: Accurate built-up area extraction is one of the most critical issues in land-cover classification. In previous studies, various techniques have been developed for built-up area extraction using Landsat images. However, the efficiency of these techniques under different technical and geographical conditions, especially for bare and sandy areas, is not optimal. One of the main challenges of built-up area extraction techniques is to determine an optimum and stable threshold with the highest possible accuracy. In many of these techniques, the optimum threshold value fluctuates substantially in different parts of the image scene. The purpose of this study is to provide a new index to improve built-up area extraction with a stable optimum threshold for different environments. In this study, the developed Automated Built-up Extraction Index (ABEI) is presented to improve the classification accuracy in areas containing bare and sandy surfaces. To develop and evaluate the accuracy of the new method for built-up area extraction with Landsat 8 OLI reflective bands, five test sites located in the Iranian cities (Babol, Naqadeh, Kashmar, Bam and Masjed Soleyman), eleven European cities (Athens, Brussels, Bucharest, Budapest, Ciechanow, Hamburg, Lyon, Madrid, Riga, Rome and Porto) and high resolution layer imperviousness (HRLI) data were used. Each site has varying environmental and complex surface coverage conditions. To determine the optimal weights for each of the Landsat 8 OLI reflective bands, the pure pixel sets for different classes and the improved gravitational search algorithm (IGSA) optimization were used. The Kappa coefficient and overall error were calculated to evaluate the accuracy of the built-up extraction map. Additionally, the ABEI performance was compared with the urban index (UI) and normalized difference built-up index (NDBI) performances. In each of the five test sites and eleven cities, the extraction accuracy of the built-up areas using the ABEI was higher than that using the UI, and NDBI (P-value of 0.01). The relative standard deviations of the optimal threshold values for the ABEI and UI were 27 and 155% (at five test sites) and were 16 and 37% (at eleven European cities), respectively, which indicates the stability of the ABEI threshold value when the location and environmental conditions change. The results of this study demonstrated that the ABEI can be used to extract built-up areas from other land covers. This index is effective even in bare soil and sandy areas, where other indices experience major challenges.

Description: Many existing satellite evapotranspiration (ET)-based drought indices have characterized regional drought condition successfully, but the relatively short time span of ET products limits their use in long-term climatological drought assessment. In this study, we assess Evaporative Drought Index (EDI) as a drought monitoring indicator over Northeast China through a retrospective comparison with drought-related indicators. After verifying its utility for detecting documented regional drought events and impacts of drought on crop production, we apply it to improve our understanding of the variation in dryness over Northeast China from 1982 to 2015. Our results illustrate that EDI is generally effective for characterizing terrestrial moisture condition and its standardized formula, namely, Standardized Evaporative Drought Index (sEDI) corresponds well with historical drought events and inter-annual grain crop yields over Northeast China. Although the calculation of sEDI does not directly incorporate precipitation and soil moisture, statistical analyses indicate sEDI can detect drought in accordance with the Standardized Precipitation Index (SPI) and Palmer Drought Severity Index (PDSI), with the highest correlations found in the west part of Northeast China (R 〈 −0.7). Further analysis illustrates sEDI is more related to commonly-used drought metrics over areas with short canopy vegetation (R 〈 −0.5) than woodland (R 〈 −0.2), which suggests precipitation may not be a good representative of drought condition over areas with deep-rooted vegetation. Then, we find 56.5% of Northeast China shows an upward dry trend from 1982 to 2015, which mainly concentrates in the west part of the study area. Conversely, 14.4% of Northeast China shows a significant wetted trend and most of them locate at cropland areas, due to the improved water management. This study suggests that EDI is a feasible method to monitor spatially distributed drought condition and can provide unique drought information not reflected by rainfall deficits, which also can be used to evaluate traditional precipitation-based indicators.

Description: Previous research has recognized the importance of edges to crime. Various scholars have explored how one specific type of edges such as physical edges or social edges affect crime, but rarely investigated the importance of the composite edge effect. To address this gap, this study introduces nightlight data from the Visible Infrared Imaging Radiometer Suite sensor on the Suomi National Polar-orbiting Partnership Satellite (NPP-VIIRS) to measure composite edges. This study defines edges as nightlight gradients—the maximum change of nightlight from a pixel to its neighbors. Using nightlight gradients and other control variables at the tract level, this study applies negative binomial regression models to investigate the effects of edges on the street robbery rate and the burglary rate in Cincinnati. The Akaike Information Criterion (AIC) of models show that nightlight gradients improve the fitness of models of street robbery and burglary. Also, nightlight gradients make a positive impact on the street robbery rate whilst a negative impact on the burglary rate, both of which are statistically significant under the alpha level of 0.05. The different impacts on these two types of crimes may be explained by the nature of crimes and the in-situ characteristics, including nightlight.

Description: Thermal infrared (TIR) target tracking is a challenging task as it entails learning an effective model to identify the target in the situation of poor target visibility and clutter background. The sparse representation, as a typical appearance modeling approach, has been successfully exploited in the TIR target tracking. However, the discriminative information of the target and its surrounding background is usually neglected in the sparse coding process. To address this issue, we propose a mask sparse representation (MaskSR) model, which combines sparse coding together with high-level semantic features for TIR target tracking. We first obtain the pixel-wise labeling results of the target and its surrounding background in the last frame, and then use such results to train target-specific deep networks using a supervised manner. According to the output features of the deep networks, the high-level pixel-wise discriminative map of the target area is obtained. We introduce the binarized discriminative map as a mask template to the sparse representation and develop a novel algorithm to collaboratively represent the reliable target part and unreliable target part partitioned with the mask template, which explicitly indicates different discriminant capabilities by label 1 and 0. The proposed MaskSR model controls the superiority of the reliable target part in the reconstruction process via a weighted scheme. We solve this multi-parameter constrained problem by a customized alternating direction method of multipliers (ADMM) method. This model is applied to achieve TIR target tracking in the particle filter framework. To improve the sampling effectiveness and decrease the computation cost at the same time, a discriminative particle selection strategy based on kernelized correlation filter is proposed to replace the previous random sampling for searching useful candidates. Our proposed tracking method was tested on the VOT-TIR2016 benchmark. The experiment results show that the proposed method has a significant superiority compared with various state-of-the-art methods in TIR target tracking.

Description: Observation of the spatial distribution of cloud optical thickness (COT) is useful for the prediction and diagnosis of photovoltaic power generation. However, there is not a one-to-one relationship between transmitted radiance and COT (so-called COT ambiguity), and it is difficult to estimate COT because of three-dimensional (3D) radiative transfer effects. We propose a method to train a convolutional neural network (CNN) based on a 3D radiative transfer model, which enables the quick estimation of the slant-column COT (SCOT) distribution from the image of a ground-mounted radiometrically calibrated digital camera. The CNN retrieves the SCOT spatial distribution using spectral features and spatial contexts. An evaluation of the method using synthetic data shows a high accuracy with a mean absolute percentage error of 18% in the SCOT range of 1–100, greatly reducing the influence of the 3D radiative effect. As an initial analysis result, COT is estimated from a sky image taken by a digital camera, and a high correlation is shown with the effective COT estimated using a pyranometer. The discrepancy between the two is reasonable, considering the difference in the size of the field of view, the space–time averaging method, and the 3D radiative effect.

Description: Monitoring soil moisture at the Earth’surface is of great importance for drought early warnings. Spaceborne remote sensing is a keystone in monitoring at continental scale, as satellites can make observations of locations which are scarcely monitored by ground-based techniques. In recent years, several soil moisture products for continental scale monitoring became available from the main space agencies around the world. Making use of sensors aboard polar satellites sampling in the microwave spectrum, soil moisture can be measured and mapped globally every few days at a spatial resolution as fine as 25 km. However, complementarity of satellite observations is a crucial issue to improve the quality of the estimations provided. In this context, measurements within the visible and infrared from geostationary satellites provide information on the surface from a totally different perspective. In this study, we design a new retrieval algorithm for daily soil moisture monitoring based only on the land surface temperature observations derived from the METEOSAT second generation geostationary satellites. Soil moisture has been retrieved from the retrieval algorithm for an eight years period over Europe and Africa at the SEVIRI sensor spatial resolution (3 km at the sub-satellite point). The results, only available for clear sky and partly cloudy conditions, are for the first time extensively evaluated against in-situ observations provided by the International Soil Moisture Network and FLUXNET at sites across Europe and Africa. The soil moisture retrievals have approximately the same accuracy as the soil moisture products derived from microwave sensors, with the most accurate estimations for semi-arid regions of Europe and Africa, and a progressive degradation of the accuracy towards northern latitudes of Europe. Although some possible improvements can be expected by a better use of other products derived from SEVIRI, the new approach developped and assessed here is a valuable alternative to microwave sensors to monitor daily soil moisture at the resolution of few kilometers over entire continents and could reveal a good complementarity to an improved monitoring system, as the algorithm can produce surface soil moisture with less than 1 day delay over clear sky and non-steady cloudy conditions (over 10% of the time).

Description: In-situ observation, climate reanalyses, and satellite remote sensing are used to study the annual cycle of turbulent latent heat flux (LHF) in the Agulhas Current system. We assess if the datasets do represent the intense exchange of moisture that occurs above the Agulhas Current and the Retroflection region, especially the new reanalyses as the former, the National Centers for Environmental Prediction Reanalysis 2 (NCEP2) and the European Centre for Medium-Range Weather Forecast (ECMWF) reanalysis second-generation reanalysis (ERA-40) have lower sea and less distinct surface temperature (SST) in the Agulhas Current system due to their low spatial resolution thus do not adequately represent the Agulhas Current LHF. We use monthly fields of LHF, SST, surface wind speed, saturated specific humidity at the sea surface (Qss), and specific humidity at 10 m (Qa). The Climate Forecast System Reanalysis (CFSR), the European Centre for Medium-Range Weather Forecast fifth generation (ERA-5), and the Modern-Era Retrospective analysis for Research and Applications version-2 (MERRA-2) are similar to the air–sea turbulent fluxes (SEAFLUX) and do represent the signature of the Agulhas Current. ERA-Interim underestimates the LHF due to lower surface wind speeds than other datasets. The observation-based National Oceanography Center Southampton (NOCS) dataset is different from all other datasets. The highest LHF of 250 W/m2 is found in the Retroflection in winter. The lowest LHF (~100 W/m2) is off Port Elizabeth in summer. East of the Agulhas Current, Qss-Qa is the main driver of the amplitude of the annual cycle of LHF, while it is the wind speed in the Retroflection and both Qss-Qa and wind speed in between. The difference in LHF between product are due to differences in Qss-Qa wind speed and resolution of datasets.

Description: Spectral-spatial classification of hyperspectral images (HSIs) has recently attracted great attention in the research domain of remote sensing. It is well-known that, in remote sensing applications, spectral features are the fundamental information and spatial patterns provide the complementary information. With both spectral features and spatial patterns, hyperspectral image (HSI) applications can be fully explored and the classification performance can be greatly improved. In reality, spatial patterns can be extracted to represent a line, a clustering of points or image texture, which denote the local or global spatial characteristic of HSIs. In this paper, we propose a spectral-spatial HSI classification model based on superpixel pattern (SP) and kernel based extreme learning machine (KELM), called SP-KELM, to identify the land covers of pixels in HSIs. In the proposed SP-KELM model, superpixel pattern features are extracted by an advanced principal component analysis (PCA), which is based on superpixel segmentation in HSIs and used to denote spatial information. The KELM method is then employed to be a classifier in the proposed spectral-spatial model with both the original spectral features and the extracted spatial pattern features. Experimental results on three publicly available HSI datasets verify the effectiveness of the proposed SP-KELM model, with the performance improvement of 10% over the spectral approaches.

Description: Migratory insect identification has been concerning entomology and pest managers for a long time. Their nocturnal behavior, as well as very small radar cross-section (RCS), makes individual detection challenging for any radar network. Typical entomological radars work at the X-band (9.4 GHz) with a vertical pencil beam. The measured RCS can be used to estimate insect mass and wingbeat frequency, and then migratory insects can be categorized into broad taxon classes using the estimated parameters. However, current entomological radars cannot achieve species identification with any higher precision or confidence. The limited frequency range of current insect radars have precluded the acquisition of more information useful for the identification of individual insects. In this paper, we report an improved measurement method of insect mass and body length using a radar with many more measurement frequencies than current entomological radars. The insect mass and body length can be extracted from the multi-frequency RCSs with uncertainties of 16.31% and 10.74%, respectively. The estimation of the thorax width and aspect ratio can also be achieved with uncertainties of 13.37% and 7.99%, respectively. Furthermore, by analyzing the statistical data of 5532 insects representing 23 species in East China, we found that the correct identification probabilities exceed 0.5 for all of the 23 species and are higher than 0.8 for 15 of the 23 species under the achievable measurement precision of the proposed technique. These findings provide promising improvements of individual parameter measurement for entomological radars and imply a possibility of species identification with higher precision.

Description: Atmospheric motion vectors (AMVs) derived from geostationary meteorological satellites have long stood as an important observational contributor to analyses of global-scale tropospheric wind patterns. This paradigm is evolving as numerical weather prediction (NWP) models and associated data assimilation systems are at the point of trying to better resolve finer scales. Understanding the physical processes that govern convectively-driven weather systems is usually hindered by a lack of observations on the scales necessary to adequately describe these events. Fortunately, satellite sensors and associated scanning strategies have improved and are now able to resolve convective-scale flow fields. Coupled with the increased availability of computing capacity and more sophisticated algorithms to track cloud motions, we are now poised to investigate the development and application of AMVs to convective-scale weather events. Our study explores this frontier using new-generation GOES-R Series imagery with a focus on hurricane applications. A proposed procedure for processing enhanced AMV datasets derived from multispectral geostationary satellite imagery for hurricane-scale analyses is described. We focus on the use of the recently available GOES-16 mesoscale domain sector rapid-scan (1-min) imagery, and emerging methods to optimally extract wind estimates (atmospheric motion vectors (AMVs)) from close-in-time sequences. It is shown that AMV datasets can be generated on spatiotemporal scales not only useful for global applications, but for mesoscale applications such as hurricanes as well.

Description: Nighttime light observations from remote sensing provide us with a timely and spatially explicit measure of human activities, and therefore enable a host of applications such as tracking urbanization and socioeconomic dynamics, evaluating armed conflicts and disasters, investigating fisheries, assessing greenhouse gas emissions and energy use, and analyzing light pollution and health effects. The new and improved sensors, algorithms, and products for nighttime lights, in association with other Earth observations and ancillary data (e.g., geo-located big data), together offer great potential for a deep understanding of human activities and related environmental consequences in a changing world. This paper reviews the advances of nighttime light sensors and products and examines the contributions of nighttime light remote sensing to perceiving the changing world from two aspects (i.e., human activities and environmental changes). Based on the historical review of the advances in nighttime light remote sensing, we summarize the challenges in current nighttime light remote sensing research and propose four strategic directions, including: Improving nighttime light data; developing a long time series of consistent nighttime light data; integrating nighttime light observations with other data and knowledge; and promoting multidisciplinary and interdisciplinary analyses of nighttime light observations.

Description: With the pressure of the increasing density of urban areas, some public infrastructures are moving to the underground to free up space above, such as utility lines, rail lines and roads. In the big data era, the three-dimensional (3D) data can be beneficial to understand the complex urban area. Comparing to spatial data and information of the above ground, we lack the precise and detailed information about underground infrastructures, such as the spatial information of underground infrastructure, the ownership of underground objects and the interdependence of infrastructures in the above and below ground. How can we map reliable 3D underground utility networks and use them in the land administration? First, to explain the importance of this work and find a possible solution, this paper observes the current issues of the existing underground utility database in Singapore. A framework for utility data governance is proposed to manage the work process from the underground utility data capture to data usage. This is the backbone to support the coordination of different roles in the utility data governance and usage. Then, an initial design of the 3D underground utility data model is introduced to describe the 3D geometric and spatial information about underground utility data and connect it to the cadastral parcel for land administration. In the case study, the newly collected data from mobile Ground Penetrating Radar is integrated with the existing utility data for 3D modelling. It is expected to explore the integration of new collected 3D data, the existing 2D data and cadastral information for land administration of underground utilities.

Description: 3D Cadastre models capture both the complex interrelations between physical objects and their corresponding legal rights, restrictions, and responsibilities. Most of the ongoing research on 3D Cadastre worldwide is focused on interrelations at the level of buildings and infrastructures. So far, the analysis of such interrelations in terms of indoor spaces, considering the time aspect, has not been explored yet. In The Netherlands, there are many examples of changes in the functionality of buildings over time. Tracking these changes is challenging, especially when the geometry of the spaces changes as well; for example, a change in functionality, from administrative to residential use of the space or a change in the geometry when merging two spaces in a building without modifying the functionality. To record the changes, a common practice is to use 2D plans for subdivisions and assign new rights, restrictions, and responsibilities to the changed spaces in a building. In the meantime, with the advances of 3D data collection techniques, the benefits of 3D models in various forms are increasingly being researched. This work explores the opportunities for using 3D point clouds to establish a platform for 3D Cadastre studies in indoor environments. We investigate the changes in time of the geometry of the building that can be automatically detected from point clouds, and how they can be linked with a Land Administration Model (LADM) and included in a 3D spatial database, to update the 3D indoor Cadastre. The results we have obtained are promising. The permanent changes (e.g., walls, rooms) are automatically distinguished from dynamic changes (e.g., human, furniture) and are linked to the space subdivisions.

Description: Differential SAR Interferometry (DInSAR) has proven its unprecedented ability and merits of monitoring ground deformation on a large scale with centimeter to millimeter accuracy. However, atmospheric artifacts due to spatial and temporal variations of the atmospheric state often affect the reliability and accuracy of its results. The commonly-known Atmospheric Phase Screen (APS) appears in the interferograms as ghost fringes not related to either topography or deformation. Atmospheric artifact mitigation remains one of the biggest challenges to be addressed within the DInSAR community. State-of-the-art research works have revealed that atmospheric artifacts can be partially compensated with empirical models, point-wise GPS zenith path delay, and numerical weather prediction models. In this study, we implement an accurate and realistic computing strategy using atmospheric reanalysis ERA5 data to estimate atmospheric artifacts. With this approach, the Line-of-Sight (LOS) path along the satellite trajectory and the monitored points is considered, rather than estimating it from the zenith path delay. Compared with the zenith delay-based method, the key advantage is that it can avoid errors caused by any anisotropic atmospheric phenomena. The accurate method is validated with Sentinel-1 data in three different test sites: Tenerife island (Spain), Almería (Spain), and Crete island (Greece). The effectiveness and performance of the method to remove APS from interferograms is evaluated in the three test sites showing a great improvement with respect to the zenith-based approach.

Description: After a period of mild eruptive activity, Stromboli showed between 2017 and 2018 a reawakening phase, with an increase in the eruptive activity starting in May 2017. The alert level of the volcano was raised from “green” (base) to “yellow” (attention) on 7 December 2017, and a small lava overflowed the crater rim on 15 December 2017. Between July 2017 and August 2018 the monitoring networks recorded nine major explosions, which are a serious hazard for Stromboli because they affect the summit area, crowded by tourists. We studied the 2017–2018 eruptive phase through the analysis of multidisciplinary data comprising thermal video-camera images, seismic, geodetic and geochemical data. We focused on the major explosion mechanism analyzing the well-recorded 1 December 2017 major explosion as a case study. We found that the 2017–2018 eruptive phase is consistent with a greater gas-rich magma supply in the shallow system. Furthermore, through the analysis of the case study major explosion, we identified precursory phases in the strainmeter and seismic data occurring 77 and 38 s before the explosive jet reached the eruptive vent, respectively. On the basis of these short-term precursors, we propose an automatic timely alarm system for major explosions at Stromboli volcano.

Description: Vegetation mapping, identifying the type and distribution of plant species, is important for analysing vegetation dynamics, quantifying spatial patterns of vegetation evolution, analysing the effects of environmental changes and predicting spatial patterns of species diversity. Such analysis can contribute to the development of targeted land management actions that maintain biodiversity and ecological functions. This paper presents a methodology for 3D vegetation mapping of a coastal dune complex using a multispectral camera mounted on an unmanned aerial system with particular reference to the Buckroney dune complex in Co. Wicklow, Ireland. Unmanned aerial systems (UAS), also known as unmanned aerial vehicles (UAV) or drones, have enabled high-resolution and high-accuracy ground-based data to be gathered quickly and easily on-site. The Sequoia multispectral sensor used in this study has green, red, red edge and near-infrared wavebands, and a regular camer with red, green and blue wavebands (RGB camera), to capture both visible and near-infrared (NIR) imagery of the land surface. The workflow of 3D vegetation mapping of the study site included establishing coordinated ground control points, planning the flight mission and camera parameters, acquiring the imagery, processing the image data and performing features classification. The data processing outcomes included an orthomosaic model, a 3D surface model and multispectral imagery of the study site, in the Irish Transverse Mercator (ITM) coordinate system. The planimetric resolution of the RGB sensor-based outcomes was 0.024 m while multispectral sensor-based outcomes had a planimetric resolution of 0.096 m. High-resolution vegetation mapping was successfully generated from these data processing outcomes. There were 235 sample areas (1 m × 1 m) used for the accuracy assessment of the classification of the vegetation mapping. Feature classification was conducted using nine different classification strategies to examine the efficiency of multispectral sensor data for vegetation and contiguous land cover mapping. The nine classification strategies included combinations of spectral bands and vegetation indices. Results show classification accuracies, based on the nine different classification strategies, ranging from 52% to 75%.

Description: Satellite-based precipitation products, especially those with high temporal and spatial resolution, constitute a potential alternative to sparse rain gauge networks for multidisciplinary research and applications. In this study, the validation of the 30-year Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record (PERSIANN-CDR) daily precipitation dataset was conducted over the Huai River Basin (HRB) of China. Based on daily precipitation data from 182 rain gauges, several continuous and categorical validation statistics combined with bias and error decomposition techniques were employed to quantitatively dissect the PERSIANN-CDR performance on daily, monthly, and annual scales. With and without consideration of non-rainfall data, this product reproduces adequate climatologic precipitation characteristics in the HRB, such as intra-annual cycles and spatial distributions. Bias analyses show that PERSIANN-CDR overestimates daily, monthly, and annual precipitation with a regional mean percent total bias of 11%. This is related closely to the larger positive false bias on the daily scale, while the negative non-false bias comes from a large underestimation of high percentile data despite overestimating lower percentile data. The systematic sub-component (error from high precipitation), which is independent of timescale, mainly leads to the PERSIANN-CDR total Mean-Square-Error (TMSE). Moreover, the daily TMSE is attributed to non-false error. The correlation coefficient (R) and Kling–Gupta Efficiency (KGE) respectively suggest that this product can well capture the temporal variability of precipitation and has a moderate-to-high overall performance skill in reproducing precipitation. The corresponding capabilities increase from the daily to annual scale, but decrease with the specified precipitation thresholds. Overall, the PERSIANN-CDR product has good (poor) performance in detecting daily low (high) rainfall events on the basis of Probability of Detection, and it has a False Alarm Ratio of above 50% for each precipitation threshold. The Equitable Threat Score and Heidke Skill Score both suggest that PERSIANN-CDR has a certain ability to detect precipitation between the second and eighth percentiles. According to the Hanssen–Kuipers Discriminant, this product can generally discriminate rainfall events between two thresholds. The Frequency Bias Index indicates an overestimation (underestimation) of precipitation totals in thresholds below (above) the seventh percentile. Also, continuous and categorical statistics for each month show evident intra-annual fluctuations. In brief, the comprehensive dissection of PERSIANN-CDR performance reported herein facilitates a valuable reference for decision-makers seeking to mitigate the adverse impacts of water deficit in the HRB and algorithm improvements in this product.

Description: Persistent scatterers interferometric Synthetic Aperture Radar (PS-InSAR) is capable of precise topography measurement up to sub-meter scale and monitoring subtle deformation up to mm/year scale for all the radar image pixels with stable radiometric characteristics. As a representative PS-InSAR method, the Stanford Method for Persistent Scatterers (StaMPS) is widely used due to its high density of PS points for both rural and urban areas. However, when it comes to layover regions, which usually happen in urban areas, the StaMPS is limited locally. Moreover, the measurement points are greatly reduced due to the removal of adjacent PS pixels. In this paper, an improved StaMPS method, called IStaMPS, is proposed. The PS pixels are selected with high density by the improved PS selection strategy. Moreover, the topography information not provided in StaMPS can be accurately measured in IStaMPS. Based on the data acquired by TerraSAR-X/TanDEM-X over the Terminal 3 E (T3 E) site of Beijing Capital International Airport and the Chaobai River of Beijing Shunyi District, a comparison between StaMPS-retrieved results and IStaMPS-retrieved ones was performed, which demonstrated that the density of PS points detected by IStaMPS is increased by about 1.8 and 1.6 times for these two areas respectively. Through comparisons of local statistical results of topography estimation and mean deformation rate, the improvement granted by the proposed IStaMPS was demonstrated for both urban areas with complex buildings or man-made targets and non-urban areas with natural targets. In terms of the spatiotemporal deformation variation, the northwest region of T3 E experienced an exceptional uplift during the period from June 2012 to August 2015, and the maximum uplift rate is approximately 4.2 mm per year.

Description: Unprecedented human-induced land cover changes happened in China after the Reform and Opening-up in 1978, matching with the era of Landsat satellite series. However, it is still unknown whether Landsat data can effectively support retrospective analysis of land cover changes in China over the past four decades. Here, for the first time, we conduct a systematic investigation on the availability of Landsat data in China, targeting its application for retrospective and continuous monitoring of land cover changes. The latter is significant to assess impact of land cover changes, and consequences of past land policy and management interventions. The total and valid observations (excluding clouds, cloud shadows, and terrain shadows) from Landsat 5/7/8 from 1984 to 2017 were quantified at pixel scale, based on the cloud computing platform Google Earth Engine (GEE). The results show higher intensity of Landsat observation in the northern part of China as compared to the southern part. The study provides an overall picture of Landsat observations suitable for satellite-based annual land cover monitoring over the entire country. We uncover that two sub-regions of China (i.e., Northeast China-Inner Mongolia-Northwest China, and North China Plain) have sufficient valid observations for retrospective analysis of land cover over 30 years (1987–2017) at an annual interval; whereas the Middle-Lower Yangtze Plain (MLYP) and Xinjiang (XJ) have sufficient observations for annual analyses for the periods 1989–2017 and 2004–2017, respectively. Retrospective analysis of land cover is possible only at a two-year time interval in South China (SC) for the years 1988–2017, Xinjiang (XJ) for the period 1992–2003, and the Tibetan Plateau (TP) during 2004–2017. For the latter geographic regions, land cover dynamics can be analyzed only at a three-year interval prior to 2004. Our retrospective analysis suggest that Landsat-based analysis of land cover dynamics at an annual interval for the whole country is not feasible; instead, national monitoring at two- or three-year intervals could be achievable. This study provides a preliminary assessment of data availability, targeting future continuous land cover monitoring in China; and the code is released to the public to facilitate similar data inventory in other regions of the world.

Description: Millimeter wave cloud radar (MMCR) is one of the primary instruments employed to observe cloud–precipitation. With appropriate data processing, measurements of the Doppler spectra, spectral moments, and retrievals can be used to study the physical processes of cloud–precipitation. This study mainly analyzed the vertical structures and microphysical characteristics of different kinds of convective cloud–precipitation in South China during the pre-flood season using a vertical pointing Ka-band MMCR. Four kinds of convection, namely, multi-cell, isolated-cell, convective–stratiform mixed, and warm-cell convection, are discussed herein. The results show that the multi-cell and convective–stratiform mixed convections had similar vertical structures, and experienced nearly the same microphysical processes in terms of particle phase change, particle size distribution, hydrometeor growth, and breaking. A forward pattern was proposed to specifically characterize the vertical structure and provide radar spectra models reflecting the different microphysical and dynamic features and variations in different parts of the cloud body. Vertical air motion played key roles in the microphysical processes of the isolated- and warm-cell convections, and deeply affected the ground rainfall properties. Stronger, thicker, and slanted updrafts caused heavier showers with stronger rain rates and groups of larger raindrops. The microphysical parameters for the warm-cell cloud–precipitation were retrieved from the radar data and further compared with the ground-measured results from a disdrometer. The comparisons indicated that the radar retrievals were basically reliable; however, the radar signal weakening caused biases to some extent, especially for the particle number concentration. Note that the differences in sensitivity and detectable height of the two instruments also contributed to the compared deviation.

Description: The non-uniformity of the relationships between urban temperature and landscape has attracted board attention. The non-uniformity in urban areas is reflected in the spatial landscape’s heterogeneity and the difference of socio-economic functions. The former is shown as the spatial differentiation of land-cover, land-use, landscape composition, and configuration, while the latter leads to the difference of the intensity of human activities and population density, which are closely related with anthropogenic heat emission. Therefore, this study introduces urban functional zones (UFZs) to express urban spatial heterogeneity. This study also attempts to comprehend urban heat island (UHI) effects and discloses the variability of urban surface temperature (LST)–landscape relationships in different kinds of UFZs. There are two main technical difficulties—how to characterize the spatial heterogeneity of UFZs and how to quantify non-uniform LST effects. A three-level variable system is established from their attributes, inner structures, and interrelationships to characterize UFZs and their LST effects hierarchically. Considering the multi-collinearity among high-dimensional variables, the Elastic Net regression method is selected for quantitative analysis. The experimental results reveal the deficiency of uniform LST analysis for heterogeneous urban areas and verify the variable relationships of LST-landscaped with different kinds of UFZs.

Description: This study improved on previous efforts to map longleaf pine (Pinus palustris) over large areas in the southeastern United States of America by developing new methods that integrate forest inventory data, aerial photography and Landsat 8 imagery to model forest characteristics. Spatial, statistical and machine learning algorithms were used to relate United States Forest Service Forest Inventory and Analysis (FIA) field plot data to relatively normalized Landsat 8 imagery based texture. Modeling algorithms employed include softmax neural networks and multiple hurdle models that combine softmax neural network predictions with linear regression models to estimate key forest characteristics across 2.3 million ha in Georgia, USA. Forest metrics include forest type, basal area and stand density. Results show strong relationships between Landsat 8 imagery based texture and field data (map accuracy 〉 0.80; square root basal area per ha residual standard errors 〈 1; natural log transformed trees per ha 〈 1.081). Model estimates depicting spatially explicit, fine resolution raster surfaces of forest characteristics for multiple coniferous and deciduous species across the study area were created and made available to the public in an online raster database. These products can be integrated with existing tabular, vector and raster databases already being used to guide longleaf pine conservation and restoration in the region.

Description: In this paper, a method for extracting Fenglin and Fengcong landform units based on karst topographic feature points is proposed. First, the variable analysis window method is used to extract peaks, nadirs, and saddle points in the karst area based on digital elevation model (DEM) data. Thiessen polygons that cover the karst surface area are constructed according to the locations of the peaks and nadirs, and the attributes of the saddles are assigned to corresponding polygons. The polygons are automatically classified via grouping analysis according to the corresponding spatial combinations of peaks, saddles, and nadirs in the Fenglin and Fengcong landform units. Then, a detailed division of the surface morphology of the karst area is achieved by distinguishing various types of Fenglin or Fengcong landform units. Experiments in the Guilin research area show that the proposed method successfully distinguishes the Fenglin and Fengcong terrain areas and extracts Fengcong landform units, individual Fenglin units, and Fenglin chains. The Fengcong area covers approximately two-thirds of the whole area, the individual Fenglin area covers approximately one-fourth, and the Fenglin chain area covers approximately one-tenth. The development of Fenglin has different stages in the Guilin area. This study provides data support for the detailed morphological study of karst terrain, and proposes a new research idea for the division and extraction of karst landform units.

Description: A scheme to semi-analytically derive waters’ Secchi depth (Zsd) from remote sensing reflectance (Rrs) considering the effects of the residual errors in satellite Rrs was developed for the China Eastern Coastal Zone (CECZ). This approach was evaluated and compared against three existing algorithms using field measurements. As it was challenging to provide the accurately inherent optical properties data for running the three existing algorithms in the extremely turbid waters, the new developed algorithm worked more effective than the latter. Moreover, with both synthetic and match-up data, the results indicated that the proposed algorithm was able to minimize some residual errors in Rrs, and thus could generate inter-mission consistent Zsd results from two ocean color missions. Finally, after application of new model to satellite images, we presented the spatial and temporal variations of Secchi depth and trophic state in the CECZ during 2002–2014. The study led to several findings: Firstly, the Zsd-based trophic state index (TSI) in the East China Sea first increased since 2002, and then gradually dropped during 2008–2014. Secondly, more and more waters within 30–35 m and 20–25 m isobaths were deteriorating from oligotrophic to mesotrophic type and from mesotrophic to eutrophic water, respectively, during 2002–2014. Lastly, the TSI increased on average 0.091 and 0.286 m per year respectively in Bohai Sea and Yellow Sea since 2002, and it might only take 14 and 67 years for Bohai Sea and Yellow Sea to deteriorate from mesotrophic to eutrophic water, following their current yearly deterioration rate and trophic trend. These results highlighted the importance to make some strict regulations for protecting the aquatic environment in the CECZ.

Description: Leaf dry matter content (LDMC), the ratio of leaf dry mass to its fresh mass, is a key plant trait, which is an indicator for many critical aspects of plant growth and survival. Accurate and fast detection of the spatiotemporal dynamics of LDMC would help understanding plants’ carbon assimilation and relative growth rate, and may then be used as an input for vegetation process models to monitor ecosystems. Satellite remote sensing is an effective tool for predicting such plant traits non-destructively. However, studies on the applicability of remote sensing for LDMC retrieval are scarce. Only a few studies have looked into the practicality of using remotely sensed data for the prediction of LDMC in a forest ecosystem. In this study, we assessed the performance of partial least squares regression (PLSR) plus 11 widely used vegetation indices (VIs), calculated based on different combinations of Sentinel-2 bands, in predicting LDMC in a coastal wetland. The accuracy of the selected methods was validated using LDMC, destructively measured in 50 randomly distributed sample plots at the study site in Schiermonnikoog, the Netherlands. The PLSR applied to canopy reflectance of Sentinel-2 bands resulted in accurate prediction of LDMC (coefficient of determination (R2) = 0.71, RMSE = 0.033). PLSR applied to the studied VIs provided an R2 of 0.70 and RMSE of 0.033. Four vegetation indices (enhanced vegetation index(EVI), specific leaf area vegetation index (SLAVI), simple ratio vegetation index (SRVI), and visible atmospherically resistant index (VARI)) computed using band 3 (green) and band 11 of the Sentinel-2 performed equally well and achieved a good measure of accuracy (R2 = 0.67, RMSE = 0.034). Our findings demonstrate the feasibility of using Sentinel-2 surface reflectance data to map LDMC in a coastal wetland.

Description: In order to improve the accuracy of subsurface target classification with ground penetrating radar (GPR) systems, it is desired to transmit and receive ultra-wide band pulses with varying combinations of polarization (a technique referred to as polarimetry). The sinuous antenna exhibits such desirable properties as ultra-wide bandwidth, polarization diversity, and low-profile form factor, making it an excellent candidate for the radiating element of such systems. However, sinuous antennas are dispersive since the active region moves with frequency along the structure, resulting in the distortion of radiated pulses. This distortion may be compensated in signal processing with accurately simulated or measured antenna phase information. However, in a practical GPR, the antenna performance may deviate from that simulated, accurate measurements may be impractical, and/or the dielectric loading of the environment may cause deviations. In such cases, it may be desirable to employ a simple dispersion model based on antenna design parameters which may be optimized in situ. This paper explores the dispersive properties of the sinuous antenna and presents a simple, adjustable, model that may be used to correct dispersed pulses. The dispersion model is successfully applied to both simulated and measured scenarios, thereby enabling the use of sinuous antennas in polarimetric GPR applications.

Description: The spectral reflectance of crop canopy is a spectral mixture, which includes soil background as one of the components. However, as soil is characterized by substantial spatial variability and temporal dynamics, its contribution to the spectral reflectance of crops will also vary. The aim of the research was to determine the impact of soil background on spectral reflectance of crop canopy in visible and near-infrared parts of the spectrum at different stages of crop development and how the soil type factor and the dynamics of soil surface affect vegetation indices calculated for crop assessment. The study was conducted on three test plots with winter wheat located in the Tula region of Russia and occupied by three contrasting types of soil. During field trips, information was collected on the spectral reflectance of winter wheat crop canopy, winter wheat leaves, weeds and open soil surface for three phenological phases (tillering, shooting stage, milky ripeness). The assessment of the soil contribution to the spectral reflectance of winter wheat crop canopy was based on a linear spectral mixture model constructed from field data. This showed that the soil background effect is most pronounced in the regions of 350–500 nm and 620–690 nm. In the shooting stage, the contribution of the soil prevails in the 620–690 nm range of the spectrum and the phase of milky ripeness in the region of 350–500 nm. The minimum contribution at all stages of winter wheat development was observed at wavelengths longer than 750 nm. The degree of soil influence varies with soil type. Analysis of variance showed that normalized difference vegetation index (NDVI) was least affected by soil type factor, the influence of which was about 30%–50%, depending on the stage of winter wheat development. The influence of soil type on soil-adjusted vegetation index (SAVI) and enhanced vegetation index (EVI2) was approximately equal and varied from 60% (shooting phase) to 80% (tillering phase). According to the discriminant analysis, the ability of vegetation indices calculated for winter wheat crop canopy to distinguish between winter wheat crops growing on different soil types changed from the classification accuracy of 94.1% (EVI2) in the tillering stage to 75% (EVI2 and SAVI) in the shooting stage to 82.6% in the milky ripeness stage (EVI2, SAVI, NDVI). The range of the sensitivity of the vegetation indices to the soil background depended on soil type. The indices showed the greatest sensitivity on gray forest soil when the wheat was in the phase of milky ripeness, and on leached chernozem when the wheat was in the tillering phase. The observed patterns can be used to develop vegetation indices, invariant to second-type soil variations caused by soil type factor, which can be applied for the remote assessment of the state of winter wheat crops.

Description: Our study aims to provide a comparison of the P- and L-band TomoSAR profiles, Land Vegetation and Ice Sensor (LVIS), and discrete return LiDAR to assess the ability for TomoSAR to monitor and estimate the tropical forest structure parameters for enhanced forest management and to support biomass missions. The comparison relies on the unique UAVSAR Jet propulsion Laboratory (JPL)/NASA L-band data, P-band data acquired by ONERA airborne system (SETHI), Small Footprint LiDAR (SFL), and NASA Land, Vegetation and Ice Sensor (LVIS) LiDAR datasets acquired in 2015 and 2016 in the frame of the AfriSAR campaign. Prior to multi-baseline data processing, a phase residual correction methodology based on phase calibration via phase center double localization has been implemented to improve the phase measurements and compensate for the phase perturbations, and disturbances originated from uncertainties in allocating flight trajectories. First, the vertical structure was estimated from L- and P-band corrected Tomography SAR data measurements, then compared with the canopy height model from SFL data. After that, the SAR and LiDAR three-dimensional (3D) datasets are compared and discussed at a qualitative basis at the region of interest. The L- and P-band’s performance for canopy penetration was assessed to determine the underlying ground locations. Additionally, the 3D records for each configuration were compared with their ability to derive forest vertical structure. Finally, the vertical structure extracted from the 3D radar reflectivity from L- and P-band are compared with SFL data, resulting in a root mean square error of 3.02 m and 3.68 m, where the coefficient of determination shows a value of 0.95 and 0.93 for P- and L-band, respectively. The results demonstrate that TomoSAR holds promise for a scientific basis in forest management activities.

Description: Although subsidence has been observed at the San Emidio geothermal field in Nevada using interferometric synthetic aperture radar since the early 1990s, the spatial extent and temporal evolution of the subsidence have not heretofore been quantified. Furthermore, the weather conditions and geographic location of San Emidio negatively affect interferometric image quality, causing low correlation amongst pairs. To address this, we introduce a new method for selecting pairs in areas of low correlation and small deformation signal using a minimum spanning tree method with a measure of image quality as the weighting criterion. We validate our pair selection approach by comparing our data products to SqueeSAR TM data products from a previous study at San Emidio. We also develop a deformation model which characterizes the spatial extent of subsidence at San Emidio in terms of volume change of the reservoir. After applying this deformation model to our data set of interferometric pairs, we examine the temporal relationship of the observed deformation with production and injection operations associated with geothermal power production.

Description: As an important tropospheric trace gas and precursor of photochemical smog, the accumulation of NO2 will cause serious air pollution. China, as the largest developing country in the world, has experienced a large amount of NO2 emissions in recent decades due to the rapid economic growth. Compared with the traditional air pollution monitoring technology, the rapid development of the remote sensing monitoring method of atmospheric satellite has gradually become the critical technical means of global atmospheric environmental monitoring. To reveal the NO2 pollution situation in China, based on the latest NO2 products from Sentinel-5P TROPOMI, the spatial–temporal characteristics and impact factors of troposphere NO2 column concentration of mainland China in the past year (February 2018 to January 2019) were analyzed on two administrative levels for the first time. Results show that the monthly fluctuation of tropospheric NO2 column concentration has obvious characteristics of “high in winter and low in summer”, while the spatial distribution forms a “high in East and low in west” pattern, bounded by Hu Line. The comparison of Coefficient of Variation (CV) and spatial autocorrelation models at two kinds of administrative scales indicates that although the spatial heterogeneity of NO2 column concentration is less affected by the observed scale, there is a “delayed effect” of about one month in the process of NO2 column concentration fluctuation. Besides, the impact factors analysis based on Spatial Lag Model (SLM) and Geographic Weighted Regression (GWR) reveals that there is a positive correlation between nighttime light intensity, the secondary and tertiary industries proportion and NO2 column concentration. Furthermore, for regions with serious NO2 pollution in North China Plain, the whole society electricity consumption and vehicle ownership also play a positive role in increasing the NO2 column concentration. This study will enlighten the government and policy makers to formulate policies tailored to local conditions, to more effectively implement NO2 emission reduction and air pollution prevention.

Description: Polarimetric synthetic aperture radar (PolSAR) image classification is a recent technology with great practical value in the field of remote sensing. However, due to the time-consuming and labor-intensive data collection, there are few labeled datasets available. Furthermore, most available state-of-the-art classification methods heavily suffer from the speckle noise. To solve these problems, in this paper, a novel semi-supervised algorithm based on self-training and superpixels is proposed. First, the Pauli-RGB image is over-segmented into superpixels to obtain a large number of homogeneous areas. Then, features that can mitigate the effects of the speckle noise are obtained using spatial weighting in the same superpixel. Next, the training set is expanded iteratively utilizing a semi-supervised unlabeled sample selection strategy that elaborately makes use of spatial relations provided by superpixels. In addition, a stacked sparse auto-encoder is self-trained using the expanded training set to obtain classification results. Experiments on two typical PolSAR datasets verified its capability of suppressing the speckle noise and showed excellent classification performance with limited labeled data.

Description: Terrestrial hydrocarbon spills have the potential to cause significant soil degradation across large areas. Identification and remedial measures taken at an early stage are therefore important. Reflectance spectroscopy is a rapid remote sensing method that has proven capable of characterizing hydrocarbon-contaminated soils. In this paper, we develop a deep learning approach to estimate the amount of Hydrocarbon (HC) mixed with different soil samples using a three-term backpropagation algorithm with dropout. The dropout was used to avoid overfitting and reduce computational complexity. A Hyspex SWIR 384 m camera measured the reflectance of the samples obtained by mixing and homogenizing four different soil types with four different HC substances, respectively. The datasets were fed into the proposed deep learning neural network to quantify the amount of HCs in each dataset. Individual validation of all the dataset shows excellent prediction estimation of the HC content with an average mean square error of ~2.2 × 10−4. The results with remote sensed data captured by an airborne system validate the approach. This demonstrates that a deep learning approach coupled with hyperspectral imaging techniques can be used for rapid identification and estimation of HCs in soils, which could be useful in estimating the quantity of HC spills at an early stage.

Description: Dritsite, ideally Li2Al4(OH)12Cl2·3H2O, is a new hydrotalcite supergroup mineral formed as a result of diagenesis in the halite−carnallite rock of the Verkhnekamskoe salt deposit, Perm Krai, Russia. Dritsite forms single lamellar or tabular hexagonal crystals up to 0.25 mm across. The mineral is transparent and colourless, with perfect cleavage on {001}. The chemical composition of dritsite (wt. %; by combination of electron microprobe and ICP−MS; H2O calculated by structure refinement) is: Li2O 6.6, Al2O3 45.42, SiO2 0.11, Cl 14.33, SO3 0.21, H2Ocalc. 34.86, O = Cl − 3.24, total 98.29. The empirical formula based on Li + Al + Si = 6 apfu (atom per formula unit) is Li1.99Al4.00Si0.01[(OH)12.19Cl1.82(SO4)0.01]Σ14.02·2.60(H2O). The Raman spectroscopic data indicate the presence of O–H bonding in the mineral, whereas CO32– groups are absent. The crystal structure has been refined in the space group P63/mcm, a = 5.0960(3), c = 15.3578(13) Å, and V = 345.4(5) Å3, to R1 = 0.088 using single-crystal data. The strongest lines of the powder X-ray diffraction pattern (d, Å (I, %) (hkl)) are: 7.68 (100) (002), 4.422 (61) (010), 3.832 (99) (004, 012), 2.561 (30) (006), 2.283 (25) (113), and 1.445 (26) (032). Dritsite was found as 2H polytype, which is isotypic with synthetic material and shows strong similarity to chlormagalumite-2H. The mineral is named in honour of the Russian crystallographer and mineralogist Prof. Victor Anatol`evich Drits.

Description: Previous literature has compared the performance of existing ground point classification (GPC) techniques on airborne LiDAR (ALS) data (LiDAR—light detection and ranging); however, their performance when applied to terrestrial LiDAR (TLS) data has not yet been addressed. This research tested the classification accuracy of five openly-available GPC algorithms on seven TLS datasets: Zhang et al.’s inverted cloth simulation (CSF), Kraus and Pfeiffer’s hierarchical weighted robust interpolation classifier (HWRI), Axelsson’s progressive TIN densification filter (TIN), Evans and Hudak’s multiscale curvature classification (MCC), and Vosselman’s modified slope-based filter (MSBF). Classification performance was analyzed using the kappa index of agreement (KIA) and rasterized spatial distribution of classification accuracy datasets generated through comparisons with manually classified reference datasets. The results identified a decrease in classification accuracy for the CSF and HWRI classification of low vegetation, for the HWRI and MCC classifications of variably sloped terrain, for the HWRI and TIN classifications of low outlier points, and for the TIN and MSBF classifications of off-terrain (OT) points without any ground points beneath. Additionally, the results show that while no single algorithm was suitable for use on all datasets containing varying terrain characteristics and OT object types, in general, a mathematical-morphology/slope-based method outperformed other methods, reporting a kappa score of 0.902.

Description: Aerial surveys in coastal areas using Unmanned Aerial Vehicles (UAVs) present many limitations. However, the need for detailed and accurate information in a marine environment has made UAVs very popular. The aim of this paper is to present a protocol which summarizes the parameters that affect the reliability of the data acquisition process over the marine environment using Unmanned Aerial Systems (UAS). The proposed UAS Data Acquisition Protocol consists of three main categories: (i) Morphology of the study area, (ii) Environmental conditions, (iii) Flight parameters. These categories include the parameters prevailing in the study area during a UAV mission and affect the quality of marine data. Furthermore, a UAS toolbox, which combines forecast weather data values with predefined thresholds and calculates the optimal flight window times in a day, was developed. The UAS toolbox was tested in two case studies with data acquisition over a coastal study area. The first UAS survey was operated under optimal conditions while the second was realized under non-optimal conditions. The acquired images and the produced orthophoto maps from both surveys present significant differences in quality. Moreover, a comparison between the classified maps of the case studies showed the underestimation of some habitats in the area at the non-optimal survey day. The UAS toolbox is expected to contribute to proper flight planning in marine applications. The UAS protocol can provide valuable information for mapping, monitoring, and management of the coastal and marine environment, which can be used globally in research and a variety of marine applications.

Description: Atmospheric scattering and surface polarization affect radiance measurements of polarization-sensitive instruments on orbit. Neglecting the polarization effects may lead to an inaccurate radiance/reflectance determination and underestimated radiance/reflectance uncertainty. Of the two instruments, CERES and VIIRS, slated to be intercalibrated by the CLARREO Pathfinder (CPF), the latter is known to be sensitive to polarization. The Pathfinder mission is tasked with accurately determining the uncertainty contribution of polarization and will provide the benchmark for the determination of the polarization correction factor for polarization-sensitive instruments. In this article, we show the formalism necessary to correct the reflectance for sensitivity to polarization after the CLARREO Pathfinder/VIIRS intercalibration, as well as the associated polarization uncertainty contribution to the overall intercalibrated reflectance error. To illustrate its usage, the formalism is applied to three dominant scene types.

Description: One of the unavoidable bottlenecks in the public application of passive signal (e.g., received signal strength, magnetic) fingerprinting-based indoor localization technologies is the extensive human effort that is required to construct and update database for indoor positioning. In this paper, we propose an accurate visual-inertial integrated geo-tagging method that can be used to collect fingerprints and construct the radio map by exploiting the crowdsourced trajectory of smartphone users. By integrating multisource information from the smartphone sensors (e.g., camera, accelerometer, and gyroscope), this system can accurately reconstruct the geometry of trajectories. An algorithm is proposed to estimate the spatial location of trajectories in the reference coordinate system and construct the radio map and geo-tagged image database for indoor positioning. With the help of several initial reference points, this algorithm can be implemented in an unknown indoor environment without any prior knowledge of the floorplan or the initial location of crowdsourced trajectories. The experimental results show that the average calibration error of the fingerprints is 0.67 m. A weighted k-nearest neighbor method (without any optimization) and the image matching method are used to evaluate the performance of constructed multisource database. The average localization error of received signal strength (RSS) based indoor positioning and image based positioning are 3.2 m and 1.2 m, respectively, showing that the quality of the constructed indoor radio map is at the same level as those that were constructed by site surveying. Compared with the traditional site survey based positioning cost, this system can greatly reduce the human labor cost, with the least external information.

Description: Unmanned aircraft systems (UAS) allow us to collect aerial data at high spatial and temporal resolution. Raw images are taken along a predetermined flight path and processed into a single raster file covering the entire study area. Radiometric calibration using empirical or manufacturer methods is required to convert raw digital numbers into reflectance and to ensure data accuracy. The performance of five radiometric calibration methods commonly used was investigated in this study. Multispectral imagery was collected using a Parrot Sequoia camera. No method maximized data accuracy in all bands. Data accuracy was higher when the empirical calibration was applied to the processed raster rather than the raw images. Data accuracy achieved with the manufacturer-recommended method was comparable to the one achieved with the best empirical method. Radiometric error in each band varied linearly with pixel radiometric values. Smallest radiometric errors were obtained in the red-edge and near-infrared (NIR) bands. Accuracy of the composite indices was higher for the pixels representing a dense vegetative cover in comparison to a lighter cover or bare soil. Results provided a better understanding of the advantages and limitations of existing radiometric calibration methods as well as the impact of the radiometric error on data quality. The authors recommend that researchers evaluate the performance of their radiometric calibration before analyzing UAS imagery and interpreting the results.

Description: Amphiboles are an important family of rock forming minerals, whose identification is crucial in provenance studies as well as in many other fields of geology, archaeology and environmental sciences. This study is aimed to find a quick way to characterize Ca-amphiboles in the tremolite (Ca2Mg5Si8O22(OH)2)–ferro–actinolite (Ca2Fe5Si8O22(OH)2) series. Raman spectroscopy is established as technique to perform non-destructive and quick analysis, with micrometric resolution, able to give the composition in terms of Mg/(Mg + Fe2+) ratio. To exploit the method, a preliminary characterization is performed by Scanning Electron Microscopy coupled with Energy-dispersed X-ray Spectroscopy (SEM-EDS). Two independent methods to evaluate the composition from the Raman data (aiming to an accuracy of about 5%), using the low-wavenumbers part of the spectrum and the OH stretching bands, are developed. The application of the proposed method to micro-Raman mappings and the possible use of handheld Raman spectroscopy to have compositional information on Ca-amphiboles are discussed.

Description: Pyrochlore group minerals are the main raw phases in granitic rocks of the Katugin complex-ore deposit that stores Nb, Ta, Y, REE, U, Th, Zr, and cryolite. There are three main types: Primary magmatic, early postmagmatic (secondary-I), and late hydrothermal (secondary-II) pyrochlores. The primary magmatic phase is fluornatropyrochlore, which has high concentrations of Na2O (to 10.5 wt.%), F (to 5.4 wt.%), and REE2O3 (to 17.3 wt.%) but also low CaO (0.6–4.3 wt.%), UO2 (to 2.6 wt.%), ThO2 (to 1.8 wt.%), and PbO (to 1.4 wt.%). Pyrochlore of this type is very rare in nature and is limited to a few occurrences: Rare-metal deposits of Nechalacho in syenite and nepheline syenite (Canada) and Mariupol in nepheline syenite (Ukraine). It may have crystallized synchronously with or slightly later than melanocratic minerals (aegirine, biotite, and arfvedsonite) at the late magmatic stage when Fe from the melt became bound, which hindered the crystallization of columbite. Secondary-I pyrochlore follows cracks or replaces primary pyrochlore in grain rims and is compositionally similar to the early phase, except for lower Na2O concentrations (2.8 wt.%), relatively low F (4 wt.%), and less complete A- and Y-sites occupancy. Secondary-II pyrochlore is a product of late hydrothermal alteration, which postdated the formation of the Katugin deposit. It differs in large ranges of elements and contains minor K, Ba, Pb, Fe, and significant Si concentrations but also low Na and F. Its composition mostly falls within the field of hydro- and keno-pyrochlore.

Description: Accurately estimating aboveground biomass (AGB) is important in many applications, including monitoring carbon stocks, investigating deforestation and forest degradation, and designing sustainable forest management strategies. Although lidar provides critical three-dimensional forest structure information for estimating AGB, acquiring comprehensive lidar coverage is often cost prohibitive. This research focused on developing a lidar sampling framework to support AGB estimation from Landsat images. Two sampling strategies, systematic and classification-based, were tested and compared. The proposed strategies were implemented over a temperate forest study site in northern New York State and the processes were then validated at a similar site located in central New York State. Our results demonstrated that while the inclusion of lidar data using systematic or classification-based sampling supports AGB estimation, the systematic sampling selection method was highly dependent on site conditions and had higher accuracy variability. Of the 12 systematic sampling plans, R2 values ranged from 0.14 to 0.41 and plot root mean square error (RMSE) ranged from 84.2 to 93.9 Mg ha−1. The classification-based sampling outperformed 75% of the systematic sampling strategies at the primary site with R2 of 0.26 and RMSE of 70.1 Mg ha−1. The classification-based lidar sampling strategy was relatively easy to apply and was readily transferable to a new study site. Adopting this method at the validation site, the classification-based sampling also worked effectively, with an R2 of 0.40 and an RMSE of 108.2 Mg ha−1 compared to the full lidar coverage model with an R2 of 0.58 and an RMSE of 96.0 Mg ha−1. This study evaluated different lidar sample selection methods to identify an efficient and effective approach to reduce the volume and cost of lidar acquisitions. The forest type classification-based sampling method described in this study could facilitate cost-effective lidar data collection in future studies.

Description: Effective monitoring of changes in the geographic distribution of cryospheric vegetation requires high-resolution and accurate baseline maps. The rationale of the present study is to compare multiple feature extraction approaches to remotely mapping vegetation in Antarctica, assessing which give the greatest accuracy and reproducibility relative to those currently available. This study provides precise, high-resolution, and refined baseline information on vegetation distribution as is required to enable future spatiotemporal change analyses of the vegetation in Antarctica. We designed and implemented a semiautomated customized normalized difference vegetation index (NDVI) approach for extracting cryospheric vegetation by incorporating very high resolution (VHR) 8-band WorldView-2 (WV-2) satellite data. The viability of state-of-the-art target detection, spectral processing/matching, and pixel-wise supervised classification feature extraction techniques are compared with the customized NDVI approach devised in this study. An extensive quantitative and comparative assessment was made by evaluating four semiautomatic feature extraction approaches consisting of 16 feature extraction standalone methods (four customized NDVI plus 12 existing methods) for mapping vegetation on Fisher Island and Stornes Peninsula in the Larsemann Hills, situated on continental east Antarctica. The results indicated that the customized NDVI approach achieved superior performance (average bias error ranged from ~6.44 ± 1.34% to ~11.55 ± 1.34%) and highest statistical stability in terms of performance when compared with existing feature extraction approaches. Overall, the accuracy analysis of the vegetation mapping relative to manually digitized reference data (supplemented by validation with ground truthing) indicated that the 16 semi-automatic mapping methods representing four general feature extraction approaches extracted vegetated area from Fisher Island and Stornes Peninsula totalling between 2.38 and 3.72 km2 (2.85 ± 0.10 km2 on average) with bias values ranging from 3.49 to 31.39% (average 12.81 ± 1.88%) and average root mean square error (RMSE) of 0.41 km2 (14.73 ± 1.88%). Further, the robustness of the analyses and results were endorsed by a cross-validation experiment conducted to map vegetation from the Schirmacher Oasis, East Antarctica. Based on the robust comparative analysis of these 16 methods, vegetation maps of the Larsemann Hills and Schirmacher Oasis were derived by ensemble merging of the five top-performing methods (Mixture Tuned Matched Filtering, Matched Filtering, Matched Filtering/Spectral Angle Mapper Ratio, NDVI-2, and NDVI-4). This study is the first of its kind to detect and map sparse and isolated vegetated patches (with smallest area of 0.25 m2) in East Antarctica using VHR data and to use ensemble merging of feature extraction methods, and provides access to an important indicator for environmental change.

Description: Under climate change and increasing water demands, groundwater depletion has become regional and global threats for water security, which is an indispensable target to achieving sustainable developments of human society and ecosystems, especially in arid and semiarid regions where groundwater is a major water source. In this study, groundwater depletion of 2003–2016 over Xinjiang in China, a typical arid region of Central Asia, is assessed using the gravity recovery and climate experiment (GRACE) satellite and the global land data assimilation system (GLDAS) datasets. In the transition of a warm-dry to a warm-wet climate in Xinjiang, increases in precipitation, soil moisture and snow water equivalent are detected, while GRACE-based groundwater storage anomalies (GWSA) exhibit significant decreasing trends with rates between-3.61 ± 0.85 mm/a of CSR-GWSA and −3.10 ± 0.91 mm/a of JPL-GWSA. Groundwater depletion is more severe in autumn and winter. The decreases in GRACE-based GWSA are in a good agreement with the groundwater statistics collected from local authorities. However, at the same time, groundwater abstraction in Xinjiang doubled, and the water supplies get more dependent on groundwater. The magnitude of groundwater depletion is about that of annual groundwater abstraction, suggesting that scientific exploitation of groundwater is the key to ensure the sustainability of freshwater withdrawals and supplies. Furthermore, GWSA changes can be well estimated by the partial least square regression (PLSR) method based on inputs of climate data. Therefore, GRACE observations provide a feasible approach for local policy makers to monitor and forecast groundwater changes to control groundwater depletion.

Description: Indigenous forests cover 24% of New Zealand and provide valuable ecosystem services. However, a national map of forest types, that is, physiognomic types, which would benefit conservation management, does not currently exist at an appropriate level of detail. While traditional forest classification approaches from remote sensing data are based on spectral information alone, the joint use of space-based optical imagery and structural information from synthetic aperture radar (SAR) and canopy metrics from air-borne Light Detection and Ranging (LiDAR) facilitates more detailed and accurate classifications of forest structure. We present a support vector machine (SVM) classification using data from the European Space Agency (ESA) Sentinel-1 and 2 missions, Advanced Land Orbiting Satellite (ALOS) PALSAR, and airborne LiDAR to produce a regional map of physiognomic types of indigenous forest. A five-fold cross-validation (repeated 100 times) of ground data showed that the highest classification accuracy of 80.5% is achieved for bands 2, 3, 4, 8, 11, and 12 from Sentinel-2, the ratio of bands VH (vertical transmit and horizontal receive) and VV (vertical transmit and vertical receive) from Sentinel-1, and mean canopy height and 97th percentile canopy height from LiDAR. The classification based on optical bands alone was 72.7% accurate and the addition of structural metrics from SAR and LiDAR increased accuracy by 7.4%. The classification accuracy is sufficient for many management applications for indigenous forest, including biodiversity management, carbon inventory, pest control, ungulate management, and disease management.

Description: How to efficiently utilize vast amounts of easily accessed aerial imageries is a critical challenge for researchers with the proliferation of high-resolution remote sensing sensors and platforms. Recently, the rapid development of deep neural networks (DNN) has been a focus in remote sensing, and the networks have achieved remarkable progress in image classification and segmentation tasks. However, the current DNN models inevitably lose the local cues during the downsampling operation. Additionally, even with skip connections, the upsampling methods cannot properly recover the structural information, such as the edge intersections, parallelism, and symmetry. In this paper, we propose the Web-Net, which is a nested network architecture with hierarchical dense connections, to handle these issues. We design the Ultra-Hierarchical Sampling (UHS) block to absorb and fuse the inter-level feature maps to propagate the feature maps among different levels. The position-wise downsampling/upsampling methods in the UHS iteratively change the shape of the inputs while preserving the number of their parameters, so that the low-level local cues and high-level semantic cues are properly preserved. We verify the effectiveness of the proposed Web-Net in the Inria Aerial Dataset and WHU Dataset. The results of the proposed Web-Net achieve an overall accuracy of 96.97% and an IoU (Intersection over Union) of 80.10% on the Inria Aerial Dataset, which surpasses the state-of-the-art SegNet 1.8% and 9.96%, respectively; the results on the WHU Dataset also support the effectiveness of the proposed Web-Net. Additionally, benefitting from the nested network architecture and the UHS block, the extracted buildings on the prediction maps are obviously sharper and more accurately identified, and even the building areas that are covered by shadows can also be correctly extracted. The verified results indicate that the proposed Web-Net is both effective and efficient for building extraction from high-resolution remote sensing images.

Description: Impervious surfaces are commonly acknowledged as major components of human settlements. The expansion of impervious surfaces could lead to a series of human−dominated environmental and ecological issues. Tracing impervious surface dynamics at a finer temporal−spatial scale is a critical way to better understand the increasingly human-dominated system of Earth. In this study, we put forward a new scheme to conduct long-term monitoring of impervious−relevant land disturbances using high frequency Landsat archives and the Google Earth Engine (GEE). First, the developed region was identified using a classification-based approach. Then, the GEE-version LandTrendr (Landsat-based detection of Trends in Disturbance and Recovery) was used to detect land disturbances, characterizing the conversion from vegetation to impervious surfaces. Finally, the actual disturbance areas within the developed regions were derived and quantitatively evaluated. A case study was conducted to detect impervious surface dynamics in Nanjing, China, from 1988 to 2018. Results show that our scheme can efficiently monitor impervious surface dynamics at yearly intervals with good accuracy. The overall accuracy (OA) of the classification results for 1988 and 2018 are 95.86% and 94.14%. Based on temporal−spatial accuracy assessments of the final detection result, the temporal accuracy is 90.75%, and the average detection time deviation is −1.28 a. The OA, precision, and recall of the sampling inspection, respectively, are 84.34%, 85.43%, and 96.37%. This scheme provides new insights into capturing the expansion of impervious−relevant land disturbances with high frequency Landsat archives in an efficient way.

Description: Precise modeling of tropospheric delay and weighted mean temperature (Tm) is critical for Global Navigation Satellite System (GNSS) positioning and meteorology. However, the model data in previous models cover a limited time span, which limits the accuracy of these models. Besides, the vertical variations of tropospheric delay and Tm are not perfectly modeled in previous studies, which affects the performance of height corrections. In this study, we used the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis from 1979 to 2017 to build a new empirical model. We first carefully modeled the lapse rates of tropospheric delay and Tm. Then we considered the temporal variations by linear trends, annual, and semi-annual variations and the spatial variations by grids. This new model can provide zenith hydrostatic delay (ZHD), zenith wet delay (ZWD), and Tm worldwide with a spatial resolution of 1° × 1°. We used the ECMWF ERA-Interim data and the radiosonde data in 2018 to validate this new model in comparison with the canonical GPT2w model. The results show that the new model has higher accuracies than the GPT2w model in all parameters. Particularly, this new model largely improves the accuracy in estimating ZHD and Tm at high-altitude (relative to the grid point height) regions.