ALBERT

Description: East Asia is the most complex region in the world for aerosol studies, as it encounters a lot of varieties of aerosols, and aerosol classification can be a challenge in this region. In the present study, we focused on the relationship between aerosol types and aerosol optical properties. We analyzed the long-term (2005–2012) data of vertical profiles of aerosol extinction coefficients, lidar ratio (Sp), and other aerosol optical properties obtained from a NASA Micro-Pulse Lidar Network and Aerosol Robotic Network site in northern Taiwan, which frequently receives Asian continental outflows. Based on aerosol extinction vertical profiles, the profiles were classified into two types: type 1 (single-layer structure) and type 2 (two-layer structure). Fall season (October–November) was the prevailing season for the Type 1, whereas type 2 mainly happened in spring (March–April). In type 1, air masses normally originated from three regional sectors, i.e., Asia continental (AC), Pacific Ocean (PO), and Southeast Asia (SA). The mean Sp values were 39 ± 17 sr, 30 ± 12 sr, and 38 ± 18 sr for the AC, PO, and SA sectors, respectively. The Sp results suggested that aerosols from the AC sector contained dust and anthropogenic particles, and aerosols from the PO sector were most likely sea salts. We further combined the EPA dust event database and backward trajectory analysis for type 2. Results showed that Sp was 41 ± 14 sr and 53 ± 21 sr for dust storm and biomass-burning events, respectively. The Sp for biomass-burning events in type 2 showed two peaks patterns. The first peak occurred within range of 30–50 sr corresponding to urban pollutant, and the second peak occurred within range of 60–80 sr in relation to biomass burning. Finally, our study summarized the Sp values for four major aerosol types over northern Taiwan, viz., urban (42 ± 18 sr), dust (34 ± 6 sr), biomass-burning (69 ± 12 sr), and oceanic (30 ± 12 sr). Our findings provide useful references for aerosol classification and air pollution identification over the western North Pacific.

Description: Urbanization is a complex process closely involving the economy, society, and population. While monitoring urban development and exploring the industry-driving force in a real-time and effective way are the prerequisites for optimizing industry structure, narrowing the urban development gap, and achieving sustainable development. Nighttime light is an effective tool to monitor urban development from a macro perspective. However, the systematic research of nighttime light spatiotemporal variation modes and the industry-driving force of urban nighttime light are still unknown. Considering these issues, this paper analyzes the spatiotemporal variation modes of the average light index (ALI) and investigates the industry-driving force of ALI in 100 major prefecture-level cities across China mainland based on National Polar-Orbiting Partnership Satellite Visible Infrared Imaging Radiometer Suite (NPP VIIRS). The conclusions are as following three aspects. First, ALI is observed a funnel pattern among four regions in spatial dimension, with low in center and high in the surrounding, and it shows 5 variation modes (“W,” “√,” “Exponent,” “Logarithm,” and “N”) in temporal dimension, of which the “√” mode accounts for the highest proportion (60%). Second, the industry structure is closely related to ALI. Besides, the factor analysis result illustrates that the secondary and tertiary industry are the driving industries of ALI. Third, the classification result based on the industry contribution rate indicates that cities driven by different industries show significant spatial distribution differences. The three major industry-driving cities are mainly distributed in central and western regions, the secondary and tertiary industry-driving cities are evenly distributed, and the tertiary industry-driving cities are mainly distributed in provincial capitals. From 2013 to 2018, the fluctuation of city distribution driven by different industries changes obviously. The number of tertiary industry-driving cities increases steadily and the three major industry-driving cities are distributed wider spatially. Additionally, the impacts of location and raw coal on ALI are discussed. In general, these findings are essential to further research urban development mode and can be considered as the reference to narrow urban development gap.

Description: A low-mass and low-volume dual-polarization L-band radiometer is introduced that has applications for ground-based remote sensing or unmanned aerial vehicle (UAV)-based mapping. With prominent use aboard the ESA Soil Moisture and Ocean Salinity (SMOS) and NASA Soil Moisture Active Passive (SMAP) satellites, L-band radiometry can be used to retrieve environmental parameters, including soil moisture, sea surface salinity, snow liquid water content, snow density, vegetation optical depth, etc. The design and testing of the air-gapped patch array antenna is introduced and is shown to provide a 3-dB full power beamwidth of 37°. We present the radio-frequency (RF) front end design, which uses direct detection architecture and a square-law power detector. Calibration is performed using two internal references, including a matched resistive source (RS) at ambient temperature and an active cold source (ACS). The radio-frequency (RF) front end does not require temperature stabilization, due to characterization of the ACS noise temperature by sky measurements. The ACS characterization procedure is presented. The noise equivalent delta (Δ) temperature (NEΔT) of the radiometer is ~0.14 K at 1 s integration time. The total antenna temperature uncertainty ranges from 0.6 to 1.5 K.

Description: Avalanche disasters are extremely destructive and catastrophic, often causing serious casualties, economic losses and surface erosion. However, far too little attention has been paid to utilizing remote sensing mapping avalanches quickly and automatically to mitigate calamity. Such endeavors are limited by formidable natural conditions, human subjective judgement and insufficient understanding of avalanches, so they have been incomplete and inaccurate. This paper presents an objective and widely serviceable method for regional auto-detection using the scattering and interference characteristics of avalanches extracted from Sentinel-1 SLC images. Six indices are established to distinguish avalanches from surrounding undisturbed snow. The active avalanche belts in Kizilkeya and Aktep of the Western TianShan Mountains in China lend urgency to this research. Implementation found that smaller avalanches can be consistently identified more accurately in descending images. Specifically, 281 and 311 avalanches were detected in the ascending and descending of Kizilkeya, respectively. The corresponding numbers on Aktep are 104 and 114, respectively. The resolution area of single avalanche detection can reach 0.09 km2. The performance of the model was excellent in all cases (areas under the curve are 0.831 and 0.940 in descending and ascending of Kizilkeya, respectively; and 0.807 and 0.938 of Aktep, respectively). Overall, the evaluation of statistical indices are POD 〉 0.75, FAR 〈 0.34, FOM 〈 0.13 and TSS 〉 0.75. The results indicate that the performance of the innovation proposed in this paper, which employs multivariate comprehensive descriptions of avalanche characteristics to actualize regional automatic detection, can be more objective, accurate, applicable and robust to a certain extent. The latest and more complete avalanche inventory generated by this design can effectively assist in addressing the increasingly severe avalanche disasters and improving public awareness of avalanches in alpine areas.

Description: An extremely heavy rainfall event hit Guangdong province, China, from 27 August to 1 September 2018. There were two different extreme rain regions, respectively, at the Pearl River estuary and eastern Guangdong, and a record-breaking daily precipitation of 1056.7 mm was observed at Gaotan station on 30 August. This paper utilizes a suite of observations from soundings, a gauge network, disdrometers, and polarimetric radars to gain insights to the two rainfall centers. The large-scale meteorological forcing, rainfall patterns, and microphysical processes, as well as radar-based precipitation signatures are investigated. It is concluded that a west-moving monsoon depression played a critical role in sustaining the moisture supply to the two extreme rain regions, and the combined orographic enhancement further contributed to the torrential rainfall over Gaotan station. The raindrop size distributions (DSD) observed at Zhuhai and Huidong stations, as well as the observed polarimetric radar signatures indicate that the rainfall at Doumen region was characterized by larger raindrops but a lower number concentration compared with that at Gaotan region. In addition, the dual-polarization radars are used to quantify precipitation intensity during this extreme event, providing timely information for flood warning and emergency management decision-making.

Description: Earth remote sensing optical satellite systems are often divided into two categories—geosynchronous and sun-synchronous. Geosynchronous systems essentially rotate with the Earth and continuously observe the same region of the Earth. Sun-synchronous systems are generally in a polar orbit and view differing regions of the Earth at the same local time. Although similar in instrument design, there are enough differences in these two types of missions that often the calibration of the instruments can be substantially different. Thus, respective calibration teams develop independent methods and do not interact regularly or often. Yet, there are numerous areas of overlap and much to learn from one another. To address this issue, a panel of experts from both types of systems was convened to discover common areas of concern, areas where improvements can be made, and recommendations for the future. As a result of the panelist’s efforts, a set of eight recommendations were developed. Those that are related to improvements of current technologies include maintaining sun-synchronous orbits (not allowing orbital decay), standardization of spectral bandpasses, and expanded use of well-developed calibration techniques such as deep convective clouds, pseudo invariant calibration sites, and lunar methodologies. New techniques for expanded calibration capability include using geosynchronous instruments as transfer radiometers, continued development of ground-based prelaunch calibration technologies, expansion of RadCalNet, and development of space-based calibration radiometer systems.

Description: As satellite observation technology improves, the number of remote sensing images significantly and rapidly increases. Therefore, a growing number of studies are focusing on remote sensing image retrieval. However, having a large number of remote sensing images considerably slows the retrieval time and takes up a great deal of memory space. The hash method is being increasingly used for rapid image retrieval because of its remarkably fast performance. At the same time, selecting samples that contain more information and greater stability to train the network has gradually become the key to improving retrieval performance. Given the above considerations, we propose a deep hash remote sensing image retrieval method, called the hard probability sampling hash retrieval method (HPSH), which combines hash code learning with hard probability sampling in a deep network. Specifically, we used a probability sampling method to select training samples, and we designed one novel hash loss function to better train the network parameters and reduce the hashing accuracy loss due to quantization. Our experimental results demonstrate that HPSH could yield an excellent representation compared with other state-of-the-art hash approaches. For the university of California, merced (UCMD) dataset, HPSH+S resulted in a mean average precision (mAP) of up to 90.9% on 16 hash bits, 92.2% on 24 hash bits, and 92.8% on 32 hash bits. For the aerial image dataset (AID), HPSH+S achieved a mAP of up to 89.8% on 16 hash bits, 93.6% on 24 hash bits, and 95.5% on 32 hash bits. For the UCMD dataset, with the use of data augmentation, our proposed approach achieved a mAP of up to 99.6% on 32 hash bits and 99.7% on 64 hash bits.

Description: Upwelling and downwelling processes play a critical role in the connectivity between offshore waters and coastal ecosystems, having relevant implications in terms of intense biogeochemical activity and global fisheries production. A variety of in situ and remote-sensing networks were used in concert with the Iberia–Biscay–Ireland (IBI) circulation forecast system, in order to investigate two persistent upwelling and downwelling events that occurred in the Northwestern (NW) Iberian coastal system during summer 2014. Special emphasis was placed on quality-controlled surface currents provided by a high-frequency radar (HFR), since this land-based technology can effectively monitor the upper layer flow over broad coastal areas in near-real time. The low-frequency spatiotemporal response of the ocean was explored in terms of wind-induced currents’ structures and immediacy of reaction. Mean kinetic energy, divergence and vorticity maps were also calculated for upwelling and downwelling favorable events, in order to verify HFR and IBI capabilities, to accurately resolve the prevailing surface circulation features, such as the locus of a persistent upwelling maximum in the vicinity of Cape Finisterre. This integrated approach proved to be well-founded to efficiently portray the three-dimensional characteristics of the NW Iberian coastal upwelling system regardless of few shortcomings detected in IBI performance, such as the misrepresentation of the most energetic surface dynamics or the overestimation of the cooling and warming associated with upwelling and downwelling conditions, respectively. Finally, the variability of the NW Iberian upwelling system was characterized by means of the development of a novel ocean-based coastal upwelling index (UI), constructed from HFR-derived hourly surface current observations (UIHFR). The proposed UIHFR was validated against two traditional UIs for 2014, to assess its credibility. Results suggest that UIHFR was able to adequately categorize and characterize a wealth of summer upwelling and downwelling events of diverse length and strength, paving the way for future investigations of the subsequent biophysical implications.

Description: Shore-based phased-array HF radars have been widely used for remotely sensing ocean surface current, wave, and wind around the world. However, phase uncertainties, especially phase distortions, in receiving elements significantly degrade the performance of beam forming and direction-of-arrival (DOA) estimation for phased-array HF radar. To address this problem, the conventional array signal model is modified by adding a direction-based phase error matrix. Subsequently, an array phase manifold calibration method using antenna responses of incoming ship echoes is proposed. Later, an assessment on the proposed array calibration method is made based on the DOA estimations and current measurements that are obtained from the datasets that were collected with a multi-frequency HF (MHF) radar. MHF radar-estimated DOAs using three calibration strategies are compared with the ship directions that are provided by an Automatic Identification System (AIS). Additionally, comparisons between the MHF radar-derived currents while using three calibration strategies and Acoustic Doppler Current Profilers (ADCP)-measured currents are made. The results indicate that the proposed array calibration method is effective in DOA estimation and current measurement for phased-array HF radars, especially in the phase distortion situation.

Description: In this study we (1) mapped the areal extent of current dust sources over Northern Africa between 8°W–31°E and 22°N - Mediterranean coast; and (2) identified and characterized the geomorphic units and soil types that emit dust from these areas. We used the full resolution (3 km) data from the MSG-SEVIRI to map dust sources over a 2-year period between 2005–2006, and examined these regions with remotely sensed images and geomorphic and soil maps. A total of 〉2600 individual dust emission events were mapped; with frequency up to 34 events in the 2-year study period. The areal extent of dust emission sources exhibited a lognormal distribution with most sources ranging from 20 to 130 km2. Most dust events were singular and related to a variety of specific geomorphic units. Dust events that created hotspots were mostly located over playas and fluvial landforms, and to a lesser extent over sand dunes and anthropogenic affected regions. About 20% of dust hotspots were offset a few kilometers from clear geomorphic units. Quantitative analysis of emissions revealed that dust sourced from various geomorphic units, among them playas (12%) and fluvial systems (10%). The importance of sand dunes as dust-emission sources greatly differs between examined datasets (7% vs. 30%). Our study emphasizes the importance of scattered dust emission events that are not considered as hotspots, as these sources are usually neglected in dust emission modeling.

Description: Surface all-wave net radiation (Rn) is a crucial variable driving many terrestrial latent heat (LE) models that estimate global LE. However, the differences between different Rn products and their impact on global LE estimates still remain unclear. In this study, we evaluated two Rn products, Global LAnd Surface Satellite (GLASS) beta version Rn and Modern-Era Retrospective Analysis for Research and Applications-version 2 (MERRA-2) Rn, from 2007–2017 using ground-measured data from 240 globally distributed in-situ radiation measurements provided by FLUXNET projects. The GLASS Rn product had higher accuracy (R2 increased by 0.04–0.26, and RMSE decreased by 2–13.3 W/m2) than the MERRA-2 Rn product for all land cover types on a daily scale, and the two Rn products differed greatly in spatial distribution and variations. We then determined the resulting discrepancies in simulated annual global LE using a simple averaging model by merging five diagnostic LE models: RS-PM model, SW model, PT-JPL model, MS-PT model, and SIM model. The validation results showed that the estimated LE from the GLASS Rn had higher accuracy (R2 increased by 0.04–0.14, and RMSE decreased by 3–8.4 W/m2) than that from the MERRA-2 Rn for different land cover types at daily scale. Importantly, the mean annual global terrestrial LE from GLASS Rn was 2.1% lower than that from the MERRA-2 Rn. Our study showed that large differences in satellite and reanalysis Rn products could lead to substantial uncertainties in estimating global terrestrial LE.

Description: Detection of terrain features (ridges, spurs, cliffs, and peaks) is a basic research topic in digital elevation model (DEM) analysis and is essential for learning about factors that influence terrain surfaces, such as geologic structures and geomorphologic processes. Detection of terrain features based on general geomorphometry is challenging and has a high degree of uncertainty, mostly due to a variety of controlling factors on surface evolution in different regions. Currently, there are different computational techniques for obtaining detailed information about terrain features using DEM analysis. One of the most common techniques is numerically identifying or classifying terrain elements where regional topologies of the land surface are constructed by using DEMs or by combining derivatives of DEM. The main drawbacks of these techniques are that they cannot differentiate between ridges, spurs, and cliffs, or result in a high degree of false positives when detecting spur lines. In this paper, we propose a new method for automatically detecting terrain features such as ridges, spurs, cliffs, and peaks, using shaded relief by controlling altitude and azimuth of illumination sources on both smooth and rough surfaces. In our proposed method, we use edge detection filters based on azimuth angle on shaded relief to identify specific terrain features. Results show that the proposed method performs similar to or in some cases better (when detecting spurs than current terrain features detection methods, such as geomorphon, curvature, and probabilistic methods.

Description: The shipborne high-frequency surface wave radar (HFSWR) platform produces six degrees of freedom (DOF) motion at sea, which affects the performance of radar target detection and remote sensing of ocean surface dynamics parameters. Motion compensation can mitigate the effect of six-DOF motion, but motion parameters (including amplitude and angular frequency) need to be known. Motion parameters obtained by using high precision sensors are affected by the precision error and time delay, thus affecting the effect of motion compensation. To obtain the motion parameters accurately and in real time, a method of identifying the motion parameters by using an artificially transmitted reference radio frequency (RF) signal generated at the shore is proposed. Based on the results of the parameter identification, the reference RF signal and the first-order radar cross-sections (RCSs) modulated by six-DOF motion of the shipborne HFSWR platform can be compensated. The identification of angular frequency is divided into two steps: (1) Preliminary identification results are obtained by using the reference RF signal; (2) the pattern search method is used to further improve the identification accuracy of angular frequency. The amplitude of translation (including surge and sway) can be identified accurately through the reference RF signal. Due to the small amplitude of rotation (including roll, pitch, and yaw), it needs to be identified by the reference RF signal and pattern search method. After identifying the motion parameters, division in the time domain is used for motion compensation. Through the simulation results, both translation and rotation have good motion compensation effects. In addition, the method of using high precision sensors to obtain motion parameters and compensation is compared with the method in this paper, the simulation results of motion compensation show that the latter is better.

Description: The information of building types is highly needed for urban planning and management, especially in high resolution building modeling in which buildings are the basic spatial unit. However, in many parts of the world, this information is still missing. In this paper, we proposed a framework to derive the information of building type using geospatial data, including point-of-interest (POI) data, building footprints, land use polygons, and roads, from Gaode and Baidu Maps. First, we used natural language processing (NLP)-based approaches (i.e., text similarity measurement and topic modeling) to automatically reclassify POI categories into which can be used to directly infer building types. Second, based on the relationship between building footprints and POIs, we identified building types using two indicators of type ratio and area ratio. The proposed framework was tested using over 440,000 building footprints in Beijing, China. Our NLP-based approaches and building type identification methods show overall accuracies of 89.0% and 78.2%, and kappa coefficient of 0.83 and 0.71, respectively. The proposed framework is transferrable to other China cities for deriving the information of building types from web mapping platforms. The data products generated from this study are of great use for quantitative urban studies at the building level.

Description: Despite the potential implications of a cropland canopy water content (CCWC) thematic product, no global remotely sensed CCWC product is currently generated. The successful launch of the Landsat-8 Operational Land Imager (OLI) in 2012, Sentinel-2A Multispectral Instrument (MSI) in 2015, followed by Sentinel-2B in 2017, make possible the opportunity for CCWC estimation at a spatial and temporal scale that can meet the demands of potential operational users. In this study, we designed and tested a novel radiative transfer model (RTM) inversion technique to combine multiple sources of a priori data in a look-up table (LUT) for inverting the NASA Harmonized Landsat Sentinel-2 (HLS) product for CCWC estimation. This study directly builds on previous research for testing the constraint of the leaf parameter (Ns) in PROSPECT, by applying those constraints in PRO4SAIL in an agricultural setting where the variability of canopy parameters are relatively minimal. In total, 225 independent leaf measurements were used to train the LUTs, and 102 field data points were collected over the 2015–2017 growing seasons for validating the inversions. The results confirm increasing a priori information and regularization yielded the best performance for CCWC estimation. Despite the relatively low variable canopy conditions, the inclusion of Ns constraints did not improve the LUT inversion. Finally, the inversion of Sentinel-2 data outperformed the inversion of Landsat-8 in the HLS product. The method demonstrated ability for HLS inversion for CCWC estimation, resulting in the first HLS-based CCWC product generated through RTM inversion.

Description: The technological growth and accessibility of Unoccupied Aerial Systems (UAS) have revolutionized the way geographic data are collected. Digital Surface Models (DSMs) are an integral component of geospatial analyses and are now easily produced at a high resolution from UAS images and photogrammetric software. Systematic testing is required to understand the strengths and weaknesses of DSMs produced from various UAS. Thus, in this study, we used photogrammetry to create DSMs using four UAS (DJI Inspire 1, DJI Phantom 4 Pro, DJI Mavic Pro, and DJI Matrice 210) to test the overall accuracy of DSM outputs across a mixed land cover study area. The accuracy and spatial variability of these DSMs were determined by comparing them to (1) 12 high-precision GPS targets (checkpoints) in the field, and (2) a DSM created from Light Detection and Ranging (LiDAR) (Velodyne VLP-16 Puck Lite) on a fifth UAS, a DJI Matrice 600 Pro. Data were collected on July 20, 2018 over a site with mixed land cover near Middleton, NS, Canada. The study site comprised an area of eight hectares (~20 acres) with land cover types including forest, vines, dirt road, bare soil, long grass, and mowed grass. The LiDAR point cloud was used to create a 0.10 m DSM which had an overall Root Mean Square Error (RMSE) accuracy of ±0.04 m compared to 12 checkpoints spread throughout the study area. UAS were flown three times each and DSMs were created with the use of Ground Control Points (GCPs), also at 0.10 m resolution. The overall RMSE values of UAS DSMs ranged from ±0.03 to ±0.06 m compared to 12 checkpoints. Next, DSMs of Difference (DoDs) compared UAS DSMs to the LiDAR DSM, with results ranging from ±1.97 m to ±2.09 m overall. Upon further investigation over respective land covers, high discrepancies occurred over vegetated terrain and in areas outside the extent of GCPs. This indicated LiDAR’s superiority in mapping complex vegetation surfaces and stressed the importance of a complete GCP network spanning the entirety of the study area. While UAS DSMs and LiDAR DSM were of comparable high quality when evaluated based on checkpoints, further examination of the DoDs exposed critical discrepancies across the study site, namely in vegetated areas. Each of the four test UAS performed consistently well, with P4P as the clear front runner in overall ranking.

Description: Radar images suffer from the impact of sidelobes. Several sidelobe-suppressing methods including the convolutional neural network (CNN)-based one has been proposed. However, the point spread function (PSF) in the radar images is sometimes spatially variant and affects the performance of the CNN. We propose the spatial-variant convolutional neural network (SV-CNN) aimed at this problem. It will also perform well in other conditions when there are spatially variant features. The convolutional kernels of the CNN can detect motifs with some distinctive features and are invariant to the local position of the motifs. This makes the convolutional neural networks widely used in image processing fields such as image recognition, handwriting recognition, image super-resolution, and semantic segmentation. They also perform well in radar image enhancement. However, the local position invariant character might not be good for radar image enhancement, when features of motifs (also known as the point spread function in the radar imaging field) vary with the positions. In this paper, we proposed an SV-CNN with spatial-variant convolution kernels (SV-CK). Its function is illustrated through a special application of enhancing the radar images. After being trained using radar images with position-codings as the samples, the SV-CNN can enhance the radar images. Because the SV-CNN reads information of the local position contained in the position-coding, it performs better than the conventional CNN. The advance of the proposed SV-CNN is tested using both simulated and real radar images.

Description: In this paper we document the design, development, results, performance and field applications of a compact directive transmit antenna for the long-range High Frequency ocean RADAR (HFR) systems operating in the International Telecommunication Union (ITU) designated 4MHz and 5MHz radiodetermination bands. The antenna design is based on the combination of the concepts of an electrically small loop with that of travelling wave antenna. This has the effect of inducing a radiated wave predominantly in a direction opposed to that of energy flow on the antenna structures. We demonstrate here that travelling wave design allows for a more compact antenna than other directive options, it has straightforward feed-point matching arrangements, and a flat frequency and phase response over an entire radiodetermination band. In situ measurements of the antenna radiation pattern, obtained with the aid of a drone, correlate well with those obtained from simulations, and show between 8dB and 30dB front-to-back suppression, with a 3dB beam width in the forward lobe of 100∘ or more. The broad-beam radiation pattern ensures proper illumination over the ocean and the significant front-to-back suppression guarantees reduced interference to terrestrial services. The proposed antenna design is compact and straight forward and can be easily deployed by minimal modifications of an existing transmission antenna. The design may be readily adapted to different environments due to the relative insensitivity of its radiation pattern and frequency response to geometric detail. The only downside to these antennas is their relatively low radiation efficiency which, however, may easily be compensated for by the available power output of a typical HFR transmitter. Antennas based on this design are currently deployed at the SeaSonde HFR sites in New South Wales, Australia, with operational ranges up to 200 km offshore despite their low radiating efficiency and the extremely low output power in use at these installations. Due to their directional pattern, it is also planned to test these antennas in phased-array Wellen RADAR (WERA) systems in both the standard receive arrays: where in-band radio frequency noise of terrestrial origin is impacting on data quality, and in the transmit array: to possibly simplify splitting, phasing and tuning requirements.

Description: Supraglacial liquid water at the margins of ice sheets has an important impact on the surface energy balance and can also influence the ice flow when supraglacial lakes drain to the bed. Optical imagery is able to monitor supraglacial lakes during the summer season. Here we developed an alternative method using polarimetric SAR from Sentinel-1 during 2017–2020 to distinguish between liquid water and other surface types at the margin of the Northeast Greenland Ice Stream. This allows the supraglacial hydrology to be monitored during the winter months too. We found that the majority of supraglacial lakes persist over winter. When comparing our results to optical data, we found significantly more water. Even during summer, many lakes are partly or fully covered by a lid of ice and snow. We used our classification results to automatically map the outlines of supraglacial lakes, create time series of water area for each lake, and hence detect drainage events. We even found several winter time drainages, which might have an important effect on ice flow. Our method has problems during the peak of the melt season, but for the rest of the year it provides crucial information for better understanding the component of supraglacial hydrology in the glaciological system.

Description: To investigate the effect of wave-induced mixing on the upper ocean structure, especially under typhoon conditions, an ocean-wave coupled model is used in this study. Two physical processes, wave-induced turbulence mixing and wave transport flux residue, are introduced. We select tropical cyclone (TC) Nepartak in the Northwest Pacific ocean as a TC example. The results show that during the TC period, the wave-induced turbulence mixing effectively increases the cooling area and cooling amplitude of the sea surface temperature (SST). The wave transport flux residue plays a positive role in reproducing the distribution of the SST cooling area. From the intercomparisons among experiments, it is also found that the wave-induced turbulence mixing has an important effect on the formation of mixed layer depth (MLD). The simulated maximum MLD is increased to 54 m and is only 1 m less than the observed value. The wave transport flux residue shows a dominant role in the mixed layer temperature (MLT) changing. The mean error of the MLT is reduced by 0.19 °C compared with the control experiment without wave mixing effects. The study shows that the effect of wave mixing should be included in the upper ocean structure modeling.

Description: The current techniques used for monitoring the blasting process in open pit mines are manual, intermittent and inefficient and can expose technical manpower to hazardous conditions. This study presents the application of unmanned aerial vehicle (UAV) systems for monitoring and improving the blasting process in open pit mines. Field experiments were conducted in different open pit mines to assess rock fragmentation, blast-induced damage on final pit walls, blast dynamics and the accuracy of blastholes including production and pre-split holes. The UAV-based monitoring was done in three different stages, including pre-blasting, blasting and post-blasting. In the pre-blasting stage, pit walls were mapped to collect structural data to predict in situ block size distribution and to develop as-built pit wall digital elevation models (DEM) to assess blast-induced damage. This was followed by mapping the production blasthole patterns implemented in the mine to investigate drillhole alignment. To monitor the blasting process, a high-speed camera was mounted on the UAV to investigate blast initiation, sequencing, misfired holes and stemming ejection. In the post-blast stage, the blasted rock pile (muck pile) was monitored to estimate fragmentation and assess muck pile configuration, heave and throw. The collected aerial data provide detailed information and high spatial and temporal resolution on the quality of the blasting process and significant opportunities for process improvement. The current challenges with regards to the application of UAVs for blasting process monitoring are discussed, and recommendations for obtaining the most value out of an UAV application are provided.

Description: A method for estimation of the turbulent energy dissipation rate from measurements by a conically scanning pulsed coherent Doppler lidar (PCDL), with allowance for the wind transport of turbulent velocity fluctuations, has been developed. The method has been tested in comparative atmospheric experiments with a Stream Line PCDL (Halo Photonics, Brockamin, Worcester, United Kingdom) and a sonic anemometer. It has been demonstrated that the method provides unbiased estimates of the dissipation rate at arbitrarily large ratios of the mean wind velocity to the linear scanning speed.

Description: Pansharpening is a typical image fusion problem, which aims to produce a high resolution multispectral (HRMS) image by integrating a high spatial resolution panchromatic (PAN) image with a low spatial resolution multispectral (MS) image. Prior arts have used either component substitution (CS)-based methods or multiresolution analysis (MRA)-based methods for this propose. Although they are simple and easy to implement, they usually suffer from spatial or spectral distortions and could not fully exploit the spatial and/or spectral information existed in PAN and MS images. By considering their complementary performances and with the goal of combining their advantages, we propose a pansharpening weight network (PWNet) to adaptively average the fusion results obtained by different methods. The proposed PWNet works by learning adaptive weight maps for different CS-based and MRA-based methods through an end-to-end trainable neural network (NN). As a result, the proposed PWN inherits the data adaptability or flexibility of NN, while maintaining the advantages of traditional methods. Extensive experiments on data sets acquired by three different kinds of satellites demonstrate the superiority of the proposed PWNet and its competitiveness with the state-of-the-art methods.

Description: Digital Elevation Models (DEMs) are widely used as a proxy for bathymetric data and several studies have attempted to improve DEM accuracy for hydrodynamic (HD) modeling. Most of these studies attempted to quantitatively improve estimates of channel conveyance (assuming a non-braided morphology) rather than accounting for the actual channel planform. Accurate representation of river conveyance and planform in a DEM is critical to HD modeling and can be achieved with a combination of remote sensing (e.g., satellite image) and field data, such as water surface elevation (WSE). Therefore, the objectives of this study are (i) to develop an algorithm for predicting channel conveyance and characterizing planform via satellite images and in situ WSE and (ii) to estimate discharge using the predicted conveyance via an HD model. The algorithm is named River Bathymetry via Satellite Image Compilation (RiBaSIC) and uses Landsat satellite imagery, Shuttle Radar Topography Mission (SRTM) DEM, Multi-Error-Removed Improved-Terrain (MERIT) DEM, and observed WSE. The algorithm is tested on four study areas along the Willamette River, Kushiyara River, Jamuna River, and Solimoes River. Channel slope and predicted hydraulic radius are subsequently estimated for approximating Manning’s roughness factor. Two-dimensional HD models using DEMs modified by the RiBaSIC algorithm and corresponding Manning’s roughness factors are employed for discharge estimation. The proposed algorithm can represent river planform and conveyance in single-channeled, meandering, wandering, and braided river reaches. Additionally, the HD models estimated discharge within 14–19% relative root mean squared error (RRMSE) in simulation of five years period.

Description: Quantifying river discharge is a critical component for hydrological studies, floodplain ecological conservation research, and water resources management. In recent years, a series of remote sensing-based discharge estimation methods have been developed. An example is the use of the near infrared (NIR) band of optical satellite images, with the principle of calculating the ratio between a stable land pixel for calibration (C) and a pixel within the river for measurement (M), applying a linear regression between C/M series and observed discharge series. This study trialed the C/M method, utilizing the Harmonized Landsat and Sentinel-2 (HLS) surface reflectance product on relatively small rivers with 30~100 m widths. Two study sites with different river characteristics and geographic settings in the Murray-Darling Basin (MDB) of Australia were selected as case studies. Two independent sets of HLS data and gauged discharge data for the 2017 and 2018 water years were acquired for modeling and validation, respectively. Results reveal high consistency between the HLS-derived discharge and gauged discharge at both sites. The Relative Root Mean Square Errors are 53% and 19%, and the Nash-Sutcliffe Efficiency coefficients are 0.24 and 0.69 for the two sites. This study supports the effectiveness of applying the fine-resolution HLS for modeling discharge on small rivers based on the C/M methodology, which also provides evidence of using multisource synthesized datasets as the input for discharge estimation.

Description: On 2 February 2018, the China Seismo-Electromagnetic Satellite (CSES) ZhangHeng 01 (ZH-01) was successfully launched, carrying on board, in addition to a suite of plasma and particle physics instruments, a high precision magnetometer package (HPM), able to observe the ultra-low frequency (ULF) waves. In this paper, a night time Pi2 pulsation observed by CSES is reported for the first time. This Pi2 event occurred on 3 September 2018, and began at 14:30 UT (02:37 magnetic local time), when the satellite was in the southern hemisphere between −49 and −13 magnetic latitude (MLAT). Kakioka (KAK) ground station in Japan detected the same Pi2 between 14:30–14:42 UT (23:30–23:42 local time). The Pi2 oscillations in the compressional, toroidal, and poloidal components at the CSES satellite and the H-component at the KAK station are investigated by estimating coherence, amplitude, and cross-phase. We noticed a high degree of similarity between the Pi2 event in the horizontal component at KAK and the ionospheric fluctuations in the compressional component at CSES. This high correlation indicated the magnetospheric source of the Pi2 event. In addition, Pi2 is exhibited clearly in the δBy component at CSES, which is highly correlated with the ground H component, so the Pi2 event could be explained by the Substorm Current Wedge (SCW). This interpretation is further confirmed by checking the compressional component of Van Allen Probe (VAP) B satellite inside the plasmasphere, which, for the first time, gives observational support for an earlier model. This ULF wave observation shows the consistency and reliability of the high precision magnetometer (HPM) equipped by two fluxgate magnetometers (FGM1 and FGM2) onboard CSES.

Description: High-lying vegetated marshes and low-lying bare mudflats have been suggested to be two stable states in intertidal ecosystems. Being able to identify the conditions enabling the shifts between these two stable states is of great importance for ecosystem management in general and the restoration of tidal marsh ecosystems in particular. However, the number of studies investigating the conditions for state shifts from bare mudflats to vegetated marshes remains relatively low. We developed a GIS approach to identify the locations of expected shifts from bare intertidal flats to vegetated marshes along a large estuary (Western Scheldt estuary, SW Netherlands), by analyzing the interactions between spatial patterns of vegetation biomass, elevation, tidal currents, and wind waves. We analyzed false-color aerial images for locating marshes, LIDAR-based digital elevation models, and spatial model simulations of tidal currents and wind waves at the whole estuary scale (~326 km²). Our results demonstrate that: (1) Bimodality in vegetation biomass and intertidal elevation co-occur; (2) the tidal currents and wind waves change abruptly at the transitions between the low-elevation bare state and high-elevation vegetated state. These findings suggest that biogeomorphic feedback between vegetation growth, currents, waves, and sediment dynamics causes the state shifts from bare mudflats to vegetated marshes. Our findings are translated into a GIS approach (logistic regression) to identify the locations of shifts from bare to vegetated states during the studied period based on spatial patterns of elevation, current, and wave orbital velocities. This GIS approach can provide a scientific basis for the management and restoration of tidal marshes.

Description: High frequency and spatially explicit irrigated land maps are important for understanding the patterns and impacts of consumptive water use by agriculture. We built annual, 30 m resolution irrigation maps using Google Earth Engine for the years 1986–2018 for 11 western states within the conterminous U.S. Our map classifies lands into four classes: irrigated agriculture, dryland agriculture, uncultivated land, and wetlands. We built an extensive geospatial database of land cover from each class, including over 50,000 human-verified irrigated fields, 38,000 dryland fields, and over 500,000 km 2 of uncultivated lands. We used 60,000 point samples from 28 years to extract Landsat satellite imagery, as well as climate, meteorology, and terrain data to train a Random Forest classifier. Using a spatially independent validation dataset of 40,000 points, we found our classifier has an overall binary classification (irrigated vs. unirrigated) accuracy of 97.8%, and a four-class overall accuracy of 90.8%. We compared our results to Census of Agriculture irrigation estimates over the seven years of available data and found good overall agreement between the 2832 county-level estimates (r 2 = 0.90), and high agreement when estimates are aggregated to the state level (r 2 = 0.94). We analyzed trends over the 33-year study period, finding an increase of 15% (15,000 km 2 ) in irrigated area in our study region. We found notable decreases in irrigated area in developing urban areas and in the southern Central Valley of California and increases in the plains of eastern Colorado, the Columbia River Basin, the Snake River Plain, and northern California.

Description: Due to the influence of equipment instability and surveying environment, scattering echoes and other factors, it is sometimes difficult to obtain high-quality sub-bottom profile (SBP) images by traditional denoising methods. In this paper, a novel SBP image denoising method is developed for obtaining underlying clean images based on a non-local low-rank framework. Firstly, to take advantage of the inherent layering structures of the SBP image, a direction image is obtained and used as a guidance image. Secondly, the robust guidance weight for accurately selecting the similar patches is given. A novel denoising method combining the weight and a non-local low-rank filtering framework is proposed. Thirdly, after discussing the filtering parameter settings, the proposed method is tested in actual measurements of sub-bottom, both in deep water and shallow water. Experimental results validate the excellent performance of the proposed method. Finally, the proposed method is verified and compared with other methods quantificationally based on the synthetic images and has achieved the total average peak signal-to-noise ratio (PSNR) of 21.77 and structural similarity index (SSIM) of 0.573, which is far better than other methods.

Description: River discharge is a fundamental hydrologic quantity that summarizes how a watershed transforms the input of precipitation into output as channelized streamflow [...]

Description: Dryland ecosystems are frequently struck by droughts. Yet, woody vegetation is often able to recover from mortality events once precipitation returns to pre-drought conditions. Climate change, however, may impact woody vegetation resilience due to more extreme and frequent droughts. Thus, better understanding how woody vegetation responds to drought events is essential. We used a phenology-based remote sensing approach coupled with field data to estimate the severity and recovery rates of a large scale die-off event that occurred in 2014–2015 in Senegal. Novel low (L-band) and high-frequency (Ku-band) passive microwave vegetation optical depth (VOD), and optical MODIS data, were used to estimate woody vegetation dynamics. The relative importance of soil, human-pressure, and before-drought vegetation dynamics influencing the woody vegetation response to the drought were assessed. The die-off in 2014–2015 represented the highest dry season VOD drop for the studied period (1989–2017), even though the 2014 drought was not as severe as the droughts in the 1980s and 1990s. The spatially explicit Die-off Severity Index derived in this study, at 500 m resolution, highlights woody plants mortality in the study area. Soil physical characteristics highly affected die-off severity and post-disturbance recovery, but pre-drought biomass accumulation (i.e., in areas that benefited from above-normal rainfall conditions before the 2014 drought) was the most important variable in explaining die-off severity. This study provides new evidence supporting a better understanding of the “greening Sahel”, suggesting that a sudden increase in woody vegetation biomass does not necessarily imply a stable ecosystem recovery from the droughts in the 1980s. Instead, prolonged above-normal rainfall conditions prior to a drought may result in the accumulation of woody biomass, creating the basis for potentially large-scale woody vegetation die-off events due to even moderate dry spells.

Description: Augmented reality can improve construction and facility management by visualizing an as-planned model on its corresponding surface for fast, easy, and correct information retrieval. This requires the localization registration of an as-built model in an as-planned model. However, the localization and registration of indoor environments fail, owing to self-similarity in an indoor environment, relatively large as-planned models, and the presence of additional unplanned objects. Therefore, this paper proposes a computer vision-based method to (1) homogenize indoor as-planned and as-built models, (2) reduce the search space of model matching, and (3) localize the structure (e.g., room) for registration of the scanned area in its as-planned model. This method extracts a representative horizontal cross section from the as-built and as-planned point clouds to make these models similar, restricts unnecessary transformation to reduce the search space, and corresponds the line features for the estimation of the registration transformation matrix. The performance of this method, in terms of registration accuracy, is evaluated on as-built point clouds of rooms and a hallway on a building floor. A rotational error of 0.005 rad and a translational error of 0.088 m are observed in the experiments. Hence, the geometric feature described on a representative cross section with transformation restrictions can be a computationally cost-effective solution for indoor localization and registration.

Description: Salt marshes provide important services to coastal ecosystems in the southeastern United States. In many locations, salt marsh habitats are threatened by coastal development and erosion, necessitating large-scale monitoring. Assessing vegetation height across the extent of a marsh can provide a comprehensive analysis of its health, as vegetation height is associated with Above Ground Biomass (AGB) and can be used to track degradation or growth over time. Traditional methods to do this, however, rely on manual measurements of stem heights that can cause harm to the marsh ecosystem. Moreover, manual measurements are limited in scale and are often time and labor intensive. Unoccupied Aircraft Systems (UAS) can provide an alternative to manual measurements and generate continuous results across a large spatial extent in a short period of time. In this study, a multirotor UAS equipped with optical Red Green Blue (RGB) and multispectral sensors was used to survey five salt marshes in Beaufort, North Carolina. Structure-from-Motion (SfM) photogrammetry of the resultant imagery allowed for continuous modeling of the entire marsh ecosystem in a three-dimensional space. From these models, vegetation height was extracted and compared to ground-based manual measurements. Vegetation heights generated from UAS data consistently under-predicted true vegetation height proportionally and a transformation was developed to predict true vegetation height. Vegetation height may be used as a proxy for Above Ground Biomass (AGB) and contribute to blue carbon estimates, which describe the carbon sequestered in marine ecosystems. Employing this transformation, our results indicate that UAS and SfM are capable of producing accurate assessments of salt marsh health via consistent and accurate vegetation height measurements.

Description: Accounting for endmember variability is a challenging issue when unmixing hyperspectral data. This paper models the variability that is associated with each endmember as a conical hull defined by extremal pixels from the data set. These extremal pixels are considered as so-called prototypal endmember spectra that have meaningful physical interpretation. Capitalizing on this data-driven modeling, the pixels of the hyperspectral image are then described as combinations of these prototypal endmember spectra weighted by bundling coefficients and spatial abundances. The proposed unmixing model not only extracts and clusters the prototypal endmember spectra, but also estimates the abundances of each endmember. The performance of the approach is illustrated thanks to experiments conducted on simulated and real hyperspectral data and it outperforms state-of-the-art methods.

Description: Assessments of long-term changes of air quality and global radiative forcing at a large scale heavily rely on satellite aerosol optical depth (AOD) datasets, particularly their temporal binning products. Although some attempts focusing on the validation of long-term satellite AOD have been conducted, there is still a lack of comprehensive quantification and understanding of the representativeness of satellite AOD at different temporal binning scales. Here, we evaluated the performances of the Moderate Resolution Imaging Spectroradiometer (MODIS) AOD products at various temporal scales by comparing the MODIS AOD datasets from both the Terra and Aqua satellites with the entire global AErosol RObotic NETwork (AERONET) observation archive between 2000 and 2017. The uncertainty levels of the MODIS hourly and daily AOD products were similarly high, indicating that MODIS AOD retrievals could be used to represent daily aerosol conditions. The MODIS data showed the reduced quality when integrated from the daily to monthly scale, where the relative mean bias (RMB) changed from 1.09 to 1.21 for MODIS Terra and from 1.04 to 1.17 for MODIS Aqua, respectively. The limitation of valid data availability within a month appeared to be the primary reason for the increased uncertainties in the monthly binning products, and the monthly data associated uncertainties could be reduced when the number of valid AOD retrievals reached 15 times in one month. At all three temporal scales, the uncertainty levels of satellite AOD products decreased with increasing AOD values. The results of this study could provide crucial information for satellite AOD users to better understand the reliability of different temporal AOD binning products and associated uncertainties in their derived long-term trends.

Description: Human activities are mainly responsible for the Aral Sea crisis, and excessive farmland expansion and unreasonable irrigation regimes are the main manifestations. The conflicting needs of agricultural water consumption and ecological water demand of the Aral Sea are increasingly prominent. However, the quantitative relationship among the water balance elements in the oasis located in the lower reaches of the Amu Darya River Basin and their impact on the retreat of the Aral Sea remain unclear. Therefore, this study focused on the water consumption of the Nukus irrigation area in the delta of the Amu Darya River and analyzed the water balance variations and their impacts on the Aral Sea. The surface energy balance algorithm for land (SEBAL) was employed to retrieve daily and seasonal evapotranspiration (ET) levels from 1992 to 2018, and a water balance equation was established based on the results of a remote sensing evapotranspiration inversion. The results indicated that the actual evapotranspiration (ETa) simulated by the SEBAL model matched the crop evapotranspiration (ETc) calculated by the Penman–Monteith method well, and the correlation coefficients between the two ETa sources were greater than 0.8. The total ETa levels in the growing seasons decreased from 1992 to 2005 and increased from 2005 to 2015, which is consistent with the changes in the cultivated land area and inflows from the Amu Darya River. In 2000, 2005 and 2010, the groundwater recharge volumes into the Aral Sea during the growing season were 6.74×109 m3, 1.56×109 m3 and 8.40×109 m3; respectively; in the dry year of 2012, regional ET exceeded the river inflow, and 2.36×109 m3 of groundwater was extracted to supplement the shortage of irrigation water. There is a significant two-year lag correlation between the groundwater level and the area of the southern Aral Sea. This study can provide useful information for water resources management in the Aral Sea region.

Description: In the literature of pan-sharpening based on neural networks, high resolution multispectral images as ground-truth labels generally are unavailable. To tackle the issue, a common method is to degrade original images into a lower resolution space for supervised training under the Wald’s protocol. In this paper, we propose an unsupervised pan-sharpening framework, referred to as “perceptual pan-sharpening”. This novel method is based on auto-encoder and perceptual loss, and it does not need the degradation step for training. For performance boosting, we also suggest a novel training paradigm, called “first supervised pre-training and then unsupervised fine-tuning”, to train the unsupervised framework. Experiments on the QuickBird dataset show that the framework with different generator architectures could get comparable results with the traditional supervised counterpart, and the novel training paradigm performs better than random initialization. When generalizing to the IKONOS dataset, the unsupervised framework could still get competitive results over the supervised ones.

Description: This study investigates the robustness of the physically-based hydrological model WaSiM (water balance and flow simulation model) for simulating hydrological processes in two data sparse small-scale inland valley catchments (Bankandi-Loffing and Mebar) in Burkina Faso. An intensive instrumentation with two weather stations, three rain recorders, 43 piezometers, and one soil moisture station was part of the general effort to reduce the scarcity of hydrological data in West Africa. The data allowed us to successfully parameterize, calibrate (2014–2015), and validate (2016) WaSiM for the Bankandi-Loffing catchment. Good model performance concerning discharge in the calibration period (R2 = 0.91, NSE = 0.88, and KGE = 0.82) and validation period (R2 = 0.82, NSE = 0.77, and KGE = 0.57) was obtained. The soil moisture (R2 = 0.7, NSE = 0.7, and KGE = 0.8) and the groundwater table (R2 = 0.3, NSE = 0.2, and KGE = 0.5) were well simulated, although not explicitly calibrated. The spatial transposability of the model parameters from the Bankandi-Loffing model was investigated by applying the best parameter-set to the Mebar catchment without any recalibration. This resulted in good model performance in 2014–2015 (R2 = 0.93, NSE = 0.92, and KGE = 0.84) and in 2016 (R2 = 0.65, NSE = 0.64, and KGE = 0.59). This suggests that the parameter-set achieved in this study can be useful for modeling ungauged inland valley catchments in the region. The water balance shows that evaporation is more important than transpiration (76% and 24%, respectively, of evapotranspiration losses) and the surface flow is very sensitive to the observed high interannual variability of rainfall. Interflow dominates the uplands, but base flow is the major component of stream flow in inland valleys. This study provides useful information for the better management of soil and scarce water resources for smallholder farming in the area.

Description: The frequency and intensity of flood quantiles and its attendant damage in agricultural establishments have generated a lot of issues in Ethiopia. Moreover, precise estimates of flood quantiles are needed for efficient design of hydraulic structures; however, quantification of these quantiles in data-scarce regions has been a continuing challenge in hydrologic design. Flood frequency analysis is thus essential to reduce possible flood damage by investigating the most suitable flood prediction model. The annual maximum discharges from six representative stations in the Upper Blue Nile River Basin were fitted to the commonly used nine statistical distributions. This study also assessed the performance evolution of the probability distributions with varying spatial scales, such that three different spatial scales of small-, medium-, and large-scale basins in the Blue Nile River Basin were considered. The performances of the candidate probability distributions were assessed using three goodness-of-fit test statistics, root mean square error, and graphical interpretation approaches to investigate the robust probability distribution for flood frequency analysis over different basin spatial scales. Based on the overall analyses, the generalized extreme value distribution was proven to be a robust model for flood frequency analysis in the study region. The generalized extreme value distribution significantly improved the performance of the flood prediction over different spatial scales. The generalized extreme value flood prediction performance improvement measured in root mean square error varied between 5.84 and 67.91% over other commonly used probability distribution models. Thus, the flood frequency analysis using the generalized extreme value distribution could be essential for the efficient planning and design of hydraulic structures in the Blue Nile River Basin. Furthermore, this study suggests that, in the future, significant efforts should be put to conduct similar flood frequency analyses over the other major river basins of Ethiopia.

Description: Many developing nations are facing severe food insecurity partly because of their dependence on rainfed agriculture. Climate variability threatens agriculture-based community livelihoods. With booming population growth, agricultural land expands, and natural resource extraction increases, leading to changes in land use and land cover characterized by persistent vegetation greening and browning. This can modify local climate variability due to changing land–atmosphere interactions. Yet, for landscapes with significant interannual variability, such as the Mount Elgon ecosystem in Kenya and Uganda, characterizing these changes is a difficult task and more robust methods have been recommended. The current study combined trend (Mann–Kendall and Sen’s slope) and breakpoint (bfast) analysis methods to comprehensively examine recent vegetation greening and browning in Mount Elgon at multiple time scales. The study used both Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) and Climate Hazards group Infrared Precipitation with Stations (CHIRPS) data and attempted to disentangle nature- versus human-driven vegetation greening and browning. Inferences from a 2019 field study were valuable in explaining some of the observed patterns. The results indicate that Mount Elgon vegetation is highly variable with both greening and browning observable at all time scales. Mann–Kendall and Sen’s slope revealed major changes (including deforestation and reforestation), while bfast detected most of the subtle vegetation changes (such as vegetation degradation), especially in the savanna and grasslands in the northeastern parts of Mount Elgon. Precipitation in the area had significantly changed (increased) in the post-2000 era than before, particularly in 2006–2010, thus influencing greening and browning during this period. The greenness–precipitation relationship was weak in other periods. The integration of Mann–Kendall and bfast proved useful in comprehensively characterizing vegetation greenness. Such a comprehensive description of Mount Elgon vegetation dynamics is an important first step to instigate policy changes for simultaneously conserving the environment and improving livelihoods that are dependent on it.

Description: Airborne laser scanning (ALS) systems tuned to the near-infrared (NIR; 1064 nm) wavelength have become the best available data source for characterizing vegetation structure. Proliferation of multi-spectral ALS (M-ALS) data with lasers tuned at two additional wavelengths (commonly 532 nm; green, and 1550 nm; short-wave infrared (SWIR)) has promoted interest in the benefit of additional wavelengths for forest inventory modelling. In this study, structural and intensity based M-ALS metrics were derived from wavelengths independently and combined to assess their value for modelling forest inventory attributes (Lorey’s height (HL), gross volume (V), and basal area (BA)) and overstorey species diversity (Shannon index (H), Simpson index (D), and species richness (R)) in a diverse mixed-wood forest in Ontario, Canada. The area-based approach (ABA) to forest attribute modelling was used, where structural- and intensity-based metrics were calculated and used as inputs for random forest models. Structural metrics from the SWIR channel (SWIRstruc) were found to be the most accurate for H and R (%RMSE = 14.3 and 14.9), and NIRstruc were most accurate for V (%RMSE = 20.4). The addition of intensity metrics marginally increased the accuracy of HL models for SWIR and combined channels (%RMSE = 7.5). Additionally, a multi-resolution (0.5, 1, 2 m) voxel analysis was performed, where intensity data were used to calculate a suite of spectral indices. Plot-level summaries of spectral indices from each voxel resolution alone, as well as combined with structural metrics from the NIR wavelength, were used as random forest predictors. The addition of structural metrics from the NIR band reduced %RMSE for all models with HL, BA, and V realizing the largest improvements. Intensity metrics were found to be important variables in the 1 m and 2 m voxel models for D and H. Overall, results indicated that structural metrics were the most appropriate. However, the inclusion of intensity metrics, and continued testing of their potential for modelling diversity indices is warranted, given minor improvements when included. Continued analyses using M-ALS intensity metrics and voxel-based indices would help to better understand the value of these data, and their future role in forest management.

Description: Transformation to Continuous Cover Forestry (CCF) is a long and difficult process in which frequent management interventions rapidly alter forest structure and dynamics with long lasting impacts. Therefore, a critical component of transformation is the acquisition of up-to-date forest inventory data to direct future management decisions. Recently, the use of single tree detection methods derived from unmanned aerial vehicle (UAV) has been identified as being a cost effective method for inventorying forests. However, the rapidly changing structure of forest stands in transformation amplifies the difficultly in transferability of current individual tree detection (ITD) methods. This study presents a novel ITD Bayesian parameter optimisation approach that uses quantile regression and external biophysical tree data sets to provide a transferable and low cost ITD approach to monitoring stands in transformation. We applied this novel method to 5 stands in a variety of transformation stages in the UK and to a independent test study site in California, USA, to assess the accuracy and transferability of this method. Requiring small amounts of training data (15 reference trees) this approach had a mean test accuracy (F-score = 0.88) and provided mean tree diameter estimates (RMSE = 5.6 cm) with differences that were not significance to the ground data (p 〈 0.05). We conclude that this method can be used to monitor forests stands in transformation and thus can also be applied to a wide range of forest structures with limited manual parameterisation between sites.

Description: Glacier surface temperature (GST) is influenced by both the energy flux from the atmosphere above and the thermal dynamics at the ice–water–debris interfaces. However, previous studies on GST are inadequate in time series research and mountain glacier surface temperature retrieval. We evaluate the GST variability at Hailuogou glacier, a temperate glacier located in Southeastern Tibetan Plateau, from 1990 to 2018. We utilized a modified mono-window algorithm to calculate the GST using the Landsat 8 thermal infrared sensor (TIRS) band 10 data and Landsat 5 thematic mapper (TM) band 6 data. Three essential parameters, including the emissivity of ice and snow, atmospheric transmittance, and effective mean atmospheric temperature, were employed in the GST algorithm. The remotely-sensed temperatures were compared with two other single-channel algorithms to validate GST algorithm’s accuracy. Results from different algorithms showed a good agreement, with a mean difference of about 0.6 ℃. Our results showed that the GST of the Hailuogou glacier, both in the upper debris-free part and the lower debris-covered tongue, has experienced a slightly increasing trend at a rate of 0.054 ℃ a−1 during the past decades. Atmospheric warming, expanding debris cover in the lower part, and a darkening debris-free accumulation area are the main causes of the warming of the glacier surface.

Description: Tourism is a primary socio-economic factor on many coastal islands. Tourism contributes to the livelihoods of the residents, but also influences natural resources and energy consumption and can become a significant driver of land conversion and environmental change. Understanding the influence of tourist-related activities is vital for sustainable tourism development. We chose Hainan Island in South China as a research area to study the influence of tourist-driven activities on environmental variables (as Land Surface Temperatures (LST) and related ecosystem variables) during the period of 2000 to 2019. In Hainan, the local economy relies heavily on tourism, with an ever-growing influx of tourists each year. We categorised location-based points of interest (POIs) into two classes, non-tourism sites and tourism-related sites, and utilised satellite data from the cloud-based platform Google Earth Engine (GEE) to extract LST and Normalized Difference Vegetation Index (NDVI) data. We analysed the LST variations, NDVI changes and the land use/land cover (LULC) changes and compared the relative difference in LST and NDVI between the tourism-related sites and non-tourism-related sites. The main findings of this study were: (1) The median LST in the tourism-related sites was relatively higher (1.3) than the LST in the non-tourism-related sites for the 20 years. Moreover, every annual mean LST of tourism-related sites was higher than the LST values in non-tourism-related sites, with an average difference of 1.2 °C for the 20 years and a maximum difference of 1.7 °C. We found higher annual LST anomalies for tourist-related sites compared to non-tourism sites after 2010, which indicated the likely positive differences in LST above the average LST during 20 years for tourism-related sites when compared against the non-tourism related sites, thus highlighting the potential influence of tourism activities on LST. (2) The annual mean NDVI value for tourism-related sites was significantly lower than for non-tourism places every year, with an average NDVI difference of 0.26 between the two sites. (3) The land cover changed significantly: croplands and forests reduced by 3.5% and 2.8% respectively, while the areas covered by orchards and urban areas increased by 2% and 72.3% respectively. These results indicate the influence of the tourism-driven activities includes the relatively high LST, vegetation degradation and land-use conversion particular to urban cover type. The outcome of this work provides a method that combines cloud-based satellite-derived data with location-based POIs data for quantifying the long-term influence of tourism-related activities on sensitive coastal ecosystems. It contributes to designing evidence-driven management plans and policies for the sustainable tourism development in coastal areas.

Description: The variation of soil moisture (SM) is a complex and synthetic process, which is impacted by numerous factors. The effects of these factors on soil moisture are dynamic. As a result, the relationship between soil moisture and explanatory variables varies with time and season. This kind of change should be considered in obtaining fine spatial resolution soil moisture products. We chose a study area with four distinct seasons in the temperate monsoon region. In this research, we established seasonal downscaling models to avoid the influence of seasonal differences. Precipitation, land surface temperature, evapotranspiration, vegetation index, land cover, elevation, slope, aspect and soil texture were taken as explanatory variables to produce fine spatial resolution SM. SM products derived from Advanced Microwave Scanning Radiometer–Earth Observing System (AMSR-E) and Advanced Microwave Scanning Radiometer 2 (AMSR2) were downscaled with the help of machine learning algorithms. We compared three machine learning algorithms of random forest (RF), support vector machine (SVM), and K-nearest neighbors (KNN) to determine the most suitable algorithm for this study. The results show that season-based downscaling is even better than continuous time series. In the analysis of seasonal differences, precipitation plays a dominant role, but its contribution rate is different in each season. Moreover, the influence of vegetation is more prominent in winter, while the influence of terrain is more important in the other three seasons. It could be noted that the accuracy of the RF model is the best among three machine learning algorithms, and the RF-downscaled products have superior matching performance to both AMSR (AMSR-E and AMSR2) SM products and in-situ measurements. The analysis indicates considering seasonal difference and the application of machine learning has high potential for spatial downscaling in remote sensing applications.

Description: Satellite remote sensing and model data play an important role in research and applications of tropical meteorology and climatology over vast, data-sparse oceans and remote continents. Since the first weather satellite was launched by NASA in 1960, a large collection of NASA’s Earth science data is freely available to the research and application communities around the world, significantly improving our overall understanding of the Earth system and environment. Established in the mid-1980s, the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC), located in Maryland, USA, is a data archive center for multidisciplinary, satellite and model assimilation data products. As one of the 12 NASA data centers in Earth sciences, GES DISC hosts several important NASA satellite missions for tropical meteorology and climatology such as the Tropical Rainfall Measuring Mission (TRMM), the Global Precipitation Measurement (GPM) Mission and the Modern-Era Retrospective analysis for Research and Applications (MERRA). Over the years, GES DISC has developed data services to facilitate data discovery, access, distribution, analysis and visualization, including Giovanni, an online analysis and visualization tool without the need to download data and software. Despite many efforts for improving data access, a significant number of challenges remain, such as finding datasets and services for a specific research topic or project, especially for inexperienced users or users outside the remote sensing community. In this article, we list and describe major NASA satellite remote sensing and model datasets and services for tropical meteorology and climatology along with examples of using the data and services, so this may help users better utilize the information in their research and applications.

Description: Data acquisition and an efficient processing method for hydrological model initialization, such as soil moisture and parameter value identification are critical for a physics-based distributed watershed modelling of flood and flood related disasters such as sediment and debris flow. Site measurements can provide accurate estimates of soil moisture, but such techniques are limited due to the number of physical sensors required to cover a large area effectively. Available satellite-based digital soil moisture data ranges from 9 km to 20 km in resolution which obscures the soil moisture details of a hill slope scale. This resolution limitation of available satellite-based distributed soil moisture data has impacted critical analysis of soil moisture resolution variance on physics-based distributed simulation results. Moreover, available satellite-based digital soil moisture data represents only a few centimeters of the top soil column and that would inform little about the effective root-zone wetness. A recently developed soil moisture estimation method called SERVES (Soil moisture Estimation of Root zone through Vegetation index-based Evapotranspiration fraction and Soil properties) overcomes this limitation of satellite-based soil moisture data by estimating distributed effective root zone soil moisture at 30 m resolution. In this study, a distributed watershed hydrological model of a sub-catchment of Reynolds Creek Experimental Watershed was developed with the GSSHA (Gridded Surface Sub-surface Hydrological Analysis) Model. SERVES soil moisture estimated at 30 m resolution was deployed in the watershed hydrological parameter value calibration and identification process. The 30 m resolution SERVES soil moisture data was resampled to 4500 m and 9000 m resolutions and was separately employed in the calibrated hydrological model to determine the soil moisture resolution effect on the model simulated outputs and the model parameter values. It was found that the simulated discharge is underestimated, infiltration rate/volume is overestimated and higher soil moisture state distribution is filtered out as the initial soil moisture resolution was coarsened. To compensate for this disparity in the simulated results, the soil saturated hydraulic conductivity value decreased with respect to the decreased resolutions.

Description: The impact of aerosol spatio-temporal variability on the Arctic radiative budget is not fully constrained. This case study focuses on the intra-Arctic modification of long-range transported aerosol and its direct aerosol radiative effect (ARE). Different types of air-borne and ground-based remote sensing observations (from Lidar and sun-photometer) revealed a high tropospheric aerosol transport episode over two parts of the European Arctic in April 2018. By incorporating the derived aerosol optical and microphysical properties into a radiative transfer model, we assessed the ARE over the two locations. Our study displayed that even in neighboring Arctic upper tropospheric levels, aged aerosol was transformed due to the interplay of removal processes (nucleation scavenging and dry deposition) and alteration of the aerosol source regions (northeast Asia and north Europe). Along the intra-Arctic transport, the coarse aerosol mode was depleted and the visible wavelength Lidar ratio (LR) increased significantly (from 15 to 64–82 sr). However, the aerosol modifications were not reflected on the ARE. More specifically, the short-wave (SW) atmospheric column ARE amounted to +4.4 - +4.9 W m−2 over the ice-covered Fram Strait and +4.5 W m−2 over the snow-covered Ny-Ålesund. Over both locations, top-of-atmosphere (TOA) warming was accompanied by surface cooling. These similarities can be attributed to the predominant accumulation mode, which drives the SW radiative budget, as well as to the similar layer altitude, solar geometry, and surface albedo conditions over both locations. However, in the context of retreating sea ice, the ARE may change even along individual transport episodes due to the ice albedo feedback.

Description: This study aimed to understand which vegetation indices (VIs) are an ideal proxy for describing phenology and interannual variability of Gross Primary Productivity (GPP) in short-rotation coppice (SRC) plantations. Canopy structure- and chlorophyll-sensitive VIs derived from Sentinel-2 images were used to estimate the start and end of the growing season (SOS and EOS, respectively) during the period 2016–2018, for an SRC poplar (Populus spp.) plantation in Lochristi (Belgium). Three different filtering methods (Savitzky–Golay (SavGol), polynomial (Polyfit) and Harmonic Analysis of Time Series (HANTS)) and five SOS- and EOS threshold methods (first derivative function, 10% and 20% percentages and 10% and 20% percentiles) were applied to identify the optimal methods for the determination of phenophases. Our results showed that the MEdium Resolution Imaging Spectrometer (MERIS) Terrestrial Chlorophyll Index (MTCI) had the best fit with GPP phenology, as derived from eddy covariance measurements, in identifying SOS- and EOS-dates. For SOS, the performance was only slightly better than for several other indices, whereas for EOS, MTCI performed markedly better. The relationship between SOS/EOS derived from GPP and VIs varied interannually. MTCI described best the seasonal pattern of the SRC plantation’s GPP (R2 = 0.52 when combining all three years). However, during the extreme dry year 2018, the Chlorophyll Red Edge Index performed slightly better in reproducing growing season GPP variability than MTCI (R2 = 0.59; R2 = 0.49, respectively). Regarding smoothing functions, Polyfit and HANTS methods showed the best (and very similar) performances. We further found that defining SOS as the date at which the 10% or 20% percentile occurred, yielded the best agreement between the VIs and the GPP; while for EOS the dates of the 10% percentile threshold came out as the best.

Description: High-resolution real-time satellite-based precipitation estimation datasets can play a more essential role in flood forecasting and risk analysis of infrastructures. This is particularly true for extended deserts or mountainous areas with sparse rain gauges like Iran. However, there are discrepancies between these satellite-based estimations and ground measurements, and it is necessary to apply adjustment methods to reduce systematic bias in these products. In this study, we apply a quantile mapping method with gauge information to reduce the systematic error of the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Cloud Classification System (PERSIANN-CCS). Due to the availability and quality of the ground-based measurements, we divide Iran into seven climate regions to increase the sample size for generating cumulative probability distributions within each region. The cumulative distribution functions (CDFs) are then employed with a quantile mapping 0.6° × 0.6° filter to adjust the values of PERSIANN-CCS. We use eight years (2009–2016) of historical data to calibrate our method, generating nonparametric cumulative distribution functions of ground-based measurements and satellite estimations for each climate region, as well as two years (2017–2018) of additional data to validate our approach. The results show that the bias correction approach improves PERSIANN-CCS data at aggregated to monthly, seasonal and annual scales for both the calibration and validation periods. The areal average of the annual bias and annual root mean square errors are reduced by 98% and 56% during the calibration and validation periods, respectively. Furthermore, the averages of the bias and root mean square error of the monthly time series decrease by 96% and 26% during the calibration and validation periods, respectively. There are some limitations in bias correction in the Southern region of the Caspian Sea because of shortcomings of the satellite-based products in recognizing orographic clouds.

Description: Collecting in situ observations from remote, high mountain rivers presents major challenges, yet real-time, high temporal resolution (e.g., daily) discharge data are critical for flood hazard mitigation and river management. In this study, we propose a method for estimating daily river discharge (RD) based on free, operational remote sensing precipitation data (Tropical Rainfall Measuring Mission (TRMM), since 2001). In this method, an exponential filter was implemented to produce a new precipitation time series from daily basin-averaged precipitation data to model the time lag of precipitation in supplying RD, and a linear-regression relationship was constructed between the filtered precipitation time series and observed discharge records. Because of different time lags in the wet season (rainfall-dominant) and dry season (snowfall-dominant), the precipitation data were processed in a segmented way (from June to October and from November to May). The method was evaluated at two hydrological gauging stations in the Upper Brahmaputra (UB) river basin, where Nash–Sutcliffe Efficiency (NSE) coefficients for Nuxia (〉0.85) and Yangcun (〉0.80) indicate good performance. By using the degree-day method to estimate the snowmelt and acquire the time series of new active precipitation (rainfall plus snowmelt) in the target basins, the discharge estimations were improved (NSE 〉 0.9 for Nuxia) compared to the original data. This makes the method applicable for most rivers on the Tibetan Plateau, which are fed mainly by precipitation (including snowfall) and are subject to limited human interference. The method also performs well for reanalysis precipitation data (Chinese Meteorological Forcing Dataset (CMFD), 1980–2000). The real-time or historical discharges can be derived from satellite precipitation data (or reanalysis data for earlier historical years) by using our method.

Description: The European Space Agency (ESA) Soil Moisture and Ocean Salinity (SMOS) mission with the MIRAS (Microwave Imaging Radiometer using Aperture Synthesis) L-band radiometer provides global soil moisture (SM) data. SM data and products from remote sensing are relatively new, but they are providing significant observations for weather forecasting, water resources management, agriculture, land surface, and climate models assessment, etc. However, the accuracy of satellite measurements is still subject to error from the retrieval algorithms and vegetation cover. Therefore, the validation of satellite measurements is crucial to understand the quality of retrieval products. The objectives of this study, precisely framed within this mission, are (i) validation of the SMOS Level 1C Brightness Temperature (TBSMOS) products in comparison with simulated products from the L-MEB model (TBL-MEB) and (ii) validation of the SMOS Level 2 SM (SMSMOS) products against ground-based measurements at 10 significant Iranian agrometeorological stations. The validations were performed for the period of January 2012 to May 2015 over the Southwest and West of Iran. The results of the validation analysis showed an RMSE ranging between 9 to 13 K and a strong correlation (R = 0.61–0.84) between TBSMOS and TBL-MEB at all stations. The bias values (0.1 to 7.5 K) showed a slight overestimation for TBSMOS at most of the stations. The results of SMSMOS validation indicated a high agreement (RMSE = 0.046–0.079 m3 m−3 and R = 0.65–0.84) between the satellite SM and in situ measurements over all the stations. The findings of this research indicated that SMSMOS shows high accuracy and agreement with in situ measurements which validate its potential. Due to the limitation of SM measurements in Iran, the SMOS products can be used in different scientific and practical applications at different Iranian study areas.

Description: European satellite missions Sentinel-1 (S1) and Sentinel-2 (S2) provide at high spatial resolution and high revisit time, respectively, radar and optical images that support a wide range of Earth surface monitoring tasks, such as Land Use/Land Cover mapping. A long-standing challenge in the remote sensing community is about how to efficiently exploit multiple sources of information and leverage their complementarity, in order to obtain the most out of radar and optical data. In this work, we propose to deal with land cover mapping in an object-based image analysis (OBIA) setting via a deep learning framework designed to leverage the multi-source complementarity provided by radar and optical satellite image time series (SITS). The proposed architecture is based on an extension of Recurrent Neural Network (RNN) enriched via a modified attention mechanism capable to fit the specificity of SITS data. Our framework also integrates a pretraining strategy that allows to exploit specific domain knowledge, shaped as hierarchy over the set of land cover classes, to guide the model training. Thorough experimental evaluations, involving several competitive approaches were conducted on two study sites, namely the Reunion island and a part of the Senegalese groundnut basin. Classification results, 79% of global accuracy on the Reunion island and 90% on the Senegalese site, respectively, have demonstrated the suitability of the proposal.

Description: Spaceborne lidar (light detection and ranging) is a very promising tool for the optical properties of global atmosphere and ocean detection. Although some studies have shown spaceborne lidar’s potential in ocean application, there is no spaceborne lidar specifically designed for ocean studies at present. In order to investigate the detection mechanism of the spaceborne lidar and analyze its detection performance, a spaceborne oceanic lidar simulator is established based on the semianalytic Monte Carlo (MC) method. The basic principle, the main framework, and the preliminary results of the simulator are presented. The whole process of the laser emitting, transmitting, and receiving is executed by the simulator with specific atmosphere–ocean optical properties and lidar system parameters. It is the first spaceborne oceanic lidar simulator for both atmosphere and ocean. The abilities of this simulator to characterize the effect of multiple scattering on the lidar signals of different aerosols, clouds, and seawaters with different scattering phase functions are presented. Some of the results of this simulator are verified by the lidar equation. It is confirmed that the simulator is beneficial to study the principle of spaceborne oceanic lidar and it can help develop a high-precision retrieval algorithm for the inherent optical properties (IOPs) of seawater.

Description: We evaluated six empirical and semianalytical models of the diffuse attenuation coefficient at 490 nm (Kd(490)) using an in situ dataset collected in the Pearl River estuary (PRE). A combined model with the most accurate performance (correlation coefficient, R2 = 0.92) was selected and applied for long-term estimation from 2003 to 2017. Physical and biological processes in the PRE over the 14-year period were investigated by applying satellite observations (MODIS/Aqua data) and season-reliant empirical orthogonal function analysis (S-EOF). In winter, the average Kd(490) was significantly higher than in the other three seasons. A slight increasing trend was observed in spring and summer, whereas a decreasing trend was observed in winter. In summer, a tongue with a relatively high Kd(490) was found in southeastern Lingdingyang Bay. In Eastern Guangdong province (GDP), the relatively higher Kd(490) value was found in autumn and winter. Based on the second mode of S-EOF, we found that the higher values in the eastern GDP extended westward and formed a distinguishable tongue in winter. The grey relational analysis revealed that chlorophyll-a concentration (Cchla) and total suspended sediment concentration (Ctsm) were two dominant contributors determining the magnitude of Kd(490) values. The Ctsm-dominated waters were generally located in coastal and estuarine turbid waters; the Cchla-dominated waters were observed in open clear ocean. The distribution of constituents-dominated area was different in the four seasons, which was affected by physical forces, including wind field, river runoff, and sea surface temperature.

Description: The authors wish to make the following corrections to this paper [...]

Description: Ocean surface whitecaps manifest surface wave breaking. Most of the whitecap data reported in the literature are based on optical observations through photographic or video recording. The air in whitecaps modifies the dielectric properties of microwave emissions and scattering. Therefore, whitecap information is intrinsic to microwave signals. This paper discusses a method to retrieve the ocean surface whitecap coverage from microwave radiometer signals.

Description: Urban expansions to adjoining greenfield sites, particularly in metropolitan regions, have become a global occurrence. Such urbanization practice results in a significant loss in ecosystem services and triggers climate change—where these changes in land cover and emissions of certain pollutants are the fundamental drivers of climate change. Despite its crucial importance, little is known on how to quantify the impact of local drivers on anthropogenic climate change. This study aims to address the question of how the impacts of local drivers on anthropogenic climate change can be measured. The study utilizes a remote sensing approach to investigate the impacts of a period of over 30 years (1989–2019) in Brisbane, Australia and its adjoining local government areas. The methodological steps of the study are two-fold. First, we measure the greenfield development and corresponding ecosystem services losses and, then, we quantify the risk of such losses attributable to direct and indirect anthropogenic climate change. The findings of the study reveal the followings: (a) the utilized remote sensing method is a useful technique in quantifying the impacts of climate change; (b) over the last 30-year period, Brisbane and its adjoining areas encountered a total loss of about USD 4.5 billion in ecosystem services, due to direct and indirect anthropogenic climate change; (c) peri-urban areas encountered the biggest losses in ecosystem service values; (d) peri-urban areas experienced the highest greenhouse gas emission production levels, and; (e) ecosystem services should be backed up by robust urban management policies—this is critical for mitigating climate change.

Description: Reliable methods for estimating wheat grain yield before harvest could help improve farm management and, if applied on a regional level, also help identify spatial factors that influence yield. Regional grain yield can be estimated using conventional methods, but the typical process is complex and labor-intensive. Here we describe the development of a streamlined approach using publicly accessible agricultural data, field-level yield, and remote sensing data from Sentinel-2 satellite to estimate regional wheat grain yield. We validated our method on wheat croplands in Navarre in northern Spain, which features heterogeneous topography and rainfall. First, this study developed stepwise multilinear equations to estimate grain yield based on various vegetation indices, which were measured at various phenological stages in order to determine the optimal timings. Second, the most suitable model was used to estimate grain yield in wheat parcels mapped from Sentinel-2 satellite images. We used a supervised pixel-based random forest classification and the estimates were compared to government-published post-harvest yield statistics. When tested, the model achieved an R2 of 0.83 in predicting grain yield at field level. The wheat parcels were mapped with an accuracy close to 86% for both overall accuracy and compared to official statistics. Third, the validated model was used to explore potential relationships of the mapped per-parcel grain yield estimation with topographic features and rainfall by using geographically weighted regressions. Topographic features and rainfall together accounted for an average for 11 to 20% of the observed spatial variation in grain yield in Navarre. These results highlight the ability of our method for estimating wheat grain yield before harvest and determining spatial factors that influence yield at the regional scale.

Description: Pine Wilt Disease is one of the most destructive pests affecting coniferous forests. After being infected by the harmful Bursaphelenchus xylophilus nematode, most trees die within one year. The complex spreading pattern of the disease and the tedious hard labor process of diagnosis involving field wood sampling followed by laboratory analysis call for alternative methods to detect and manage the infected areas. Remote sensing comes naturally into play owing to the possibility of covering relatively large areas and the ability to discriminate healthy from sick trees based on spectral characteristics. This paper presents the development of machine learning classification algorithms for the detection of Pine Wilt Disease in Pinus pinaster, performed in the framework of the European Commission’s Horizon 2020 project “Operational Forest Monitoring using Copernicus and UAV Hyperspectral Data” (FOCUS) in two provinces of central Portugal. Five flight campaigns have been carried out in two consecutive years in order to capture a multitemporal variation of disease distribution. Classification algorithms based on a Random Forest approach were separately designed for the acquired very-high-resolution multispectral and hyperspectral data, respectively. Both algorithms achieved overall accuracies higher than 0.91 in test data. Furthermore, our study shows that the early detection of decaying trees is feasible, even before symptoms are visible in the field.

Description: This study identifies the evolution of tropical vertical cloud regimes (CRs) and their associated heating structures on the intraseasonal time scales. Using the cloud classification retrievals of CloudSat during boreal winter between 2006 and 2017, the CR index is defined as the leading pair of the combined multivariate empirical orthogonal functions of the daily mean frequency of deep, high, and low clouds over the tropical Indian Ocean, Maritime Continents, and the Western Pacific. The principal components of the CR index exhibit robust temporal variance in the 30 to 80 day intraseasonal band. Based on the propagation stages of the CRs, the coherent vertical structures of cloud composition and large-scale moisture and vertical motion exhibit a westward-tilted structure. The associated Q1-QR diabatic heating and cloud radiative forcing are consistent with the key characteristics of the Madden Julian Oscillation (MJO) documented in the previous studies. Lastly, an MJO case study showcases that the presented approach characteristically captures the propagation of moisture, cloud vertical structure, and precipitation activity across spatial and temporal scales. The current results suggest that the CR index can potentially serve as an evaluation metric to cloud-associated processes in the simulated tropical intraseasonal variability in global climate models.

Description: Unmanned aerial vehicles (UAVs) are quickly emerging as a popular platform for 3D reconstruction/modeling in various applications such as precision agriculture, coastal monitoring, and emergency management. For such applications, LiDAR and frame cameras are the two most commonly used sensors for 3D mapping of the object space. For example, point clouds for the area of interest can be directly derived from LiDAR sensors onboard UAVs equipped with integrated global navigation satellite systems and inertial navigation systems (GNSS/INS). Imagery-based mapping, on the other hand, is considered to be a cost-effective and practical option and is often conducted by generating point clouds and orthophotos using structure from motion (SfM) techniques. Mapping with photogrammetric approaches requires accurate camera interior orientation parameters (IOPs), especially when direct georeferencing is utilized. Most state-of-the-art approaches for determining/refining camera IOPs depend on ground control points (GCPs). However, establishing GCPs is expensive and labor-intensive, and more importantly, the distribution and number of GCPs are usually less than optimal to provide adequate control for determining and/or refining camera IOPs. Moreover, consumer-grade cameras with unstable IOPs have been widely used for mapping applications. Therefore, in such scenarios, where frequent camera calibration or IOP refinement is required, GCP-based approaches are impractical. To eliminate the need for GCPs, this study uses LiDAR data as a reference surface to perform in situ refinement of camera IOPs. The proposed refinement strategy is conducted in three main steps. An image-based sparse point cloud is first generated via a GNSS/INS-assisted SfM strategy. Then, LiDAR points corresponding to the resultant image-based sparse point cloud are identified through an iterative plane fitting approach and are referred to as LiDAR control points (LCPs). Finally, IOPs of the utilized camera are refined through a GNSS/INS-assisted bundle adjustment procedure using LCPs. Seven datasets over two study sites with a variety of geomorphic features are used to evaluate the performance of the developed strategy. The results illustrate the ability of the proposed approach to achieve an object space absolute accuracy of 3–5 cm (i.e., 5–10 times the ground sampling distance) at a 41 m flying height.

Description: With the increasing utilization of satellite-based soil moisture products, a primary challenge is knowing their accuracy and robustness. This study presents a comprehensive assessment over China of three widely used global satellite soil moisture products, i.e., Soil Moisture Active Passive (SMAP), European Space Agency (ESA) Climate Change Initiative (CCI) Soil Moisture, Soil Moisture and Ocean Salinity (SMOS). In situ soil moisture from 1682 stations and Variable Infiltration Capacity (VIC) model are used to evaluate the performance of SMAP_L3, ESA_CCI_SM_COMBINED, SMOS_CATDS_L3 from 31 March 2015 to 3 June 2018. The Triple Collocation (TC) approach is used to minimize the uncertainty (e.g., scale issue) during the validation process. The TC analysis is conducted using three triplets, i.e., [SMAP-Insitu-VIC], [CCI-Insitu-VIC], [SMOS-Insitu-VIC]. In general, SMAP is the most reliable product, reflecting the main spatiotemporal characteristics of soil moisture, while SMOS has the lowest accuracy. The results demonstrate that the overall root mean square error of SMAP, CCI, SMOS is 0.040, 0.028, 0.107 m3m−3, respectively. The overall temporal correlation coefficient of SMAP, CCI, SMOS is 0.68, 0.65, 0.38, respectively. The overall fractional root mean square error of SMAP, CCI, SMOS is 0.707, 0.750, 0.897, respectively. In irrigated areas, the accuracy of CCI is reduced due to the land surface model (which does not consider irrigation) used for the rescaling of the CCI_COMBINED soil moisture product during the merging process, while SMAP and SMOS preserve the irrigation signal. The quality of SMOS is most strongly impacted by land surface temperature, vegetation, and soil texture, while the quality of CCI is the least affected by these factors. With the increase of Radio Frequency Interference, the accuracy of SMOS decreases dramatically, followed by SMAP and CCI. Higher representativeness error of in situ stations is noted in regions with higher topographic complexity. This study helps to provide a guideline for the application of satellite soil moisture products in scientific research and gives some references (e.g., modify data algorithm according to the main error sources) for improving the data quality.

Description: Non-contact vital signs monitoring using microwave Doppler radar has shown great promise in healthcare applications. Recently, this unobtrusive form of physiological sensing has also been gaining attention for its potential for continuous identity authentication, which can reduce the vulnerability of traditional one-pass validation authentication systems. Physiological Doppler radar is an attractive approach for continuous identity authentication as it requires neither contact nor line-of-sight and does not give rise to privacy concerns associated with video imaging. This paper presents a review of recent advances in radar-based identity authentication systems. It includes an evaluation of the applicability of different research efforts in authentication using respiratory patterns and heart-based dynamics. It also identifies aspects of future research required to address remaining challenges in applying unobtrusive respiration-based or heart-based identity authentication to practical systems. With the advancement of machine learning and artificial intelligence, radar-based continuous authentication can grow to serve a wide range of valuable functions in society.

Description: In marine habitat mapping, a demand exists for high-resolution maps of the seafloor both for marine spatial planning and research. One topic of interest is the detection of boulders in side scan sonar backscatter mosaics of continental shelf seas. Boulders are oftentimes numerous, but encompass few pixels in backscatter mosaics. Therefore, both their automatic and manual detection is difficult. In this study, located in the German Baltic Sea, the use of super resolution by deep learning to improve the manual and automatic detection of boulders in backscatter mosaics is explored. It is found that upscaling of mosaics by a factor of 2 to 0.25 m or 0.125 m resolution increases the performance of small boulder detection and boulder density grids. Upscaling mosaics with 1.0 m pixel resolution by a factor of 4 improved performance, but the results are not sufficient for practical application. It is suggested that mosaics of 0.5 m resolution can be used to create boulder density grids in the Baltic Sea in line with current standards following upscaling.

Description: Highly fragmented land property hinders the planning and management of single species tree plantations. In such situations, acquiring information about the available resources is challenging. This study aims to propose a method to locate and characterize tree plantations in these cases. Galicia (Northwest of Spain) is an area where property is extremely divided into small parcels. European chestnut (Castanea sativa) plantations are an important source of income there; however, it is often difficult to obtain information about them due to their small size and scattered distribution. Therefore, we selected a Galician region with a high presence of chestnut plantations as a case study area in order to locate and characterize small plantations using open-access data. First, we detected the location of chestnut plantations applying a supervised classification for a combination of: Sentinel-2 images and the open-access low-density Light Detection and Ranging (LiDAR) point clouds, obtained from the untapped open-access LiDAR Spanish national database. Three classification algorithms were used: Random Forest (RF), Support Vector Machine (SVM), and XGBoost. We later characterized the plots at the tree-level using the LiDAR point-cloud. We detected individual trees and obtained their height applying a local maxima algorithm to a point-cloud-derived Canopy Height Model (CHM). We also calculated the crown surface of each tree by applying a method based on two-dimensional (2D) tree shape reconstruction and canopy segmentation to a projection of the LiDAR point cloud. Chestnut plantations were detected with an overall accuracy of 81.5%. Individual trees were identified with a detection rate of 96%. The coefficient of determination R2 value for tree height estimation was 0.83, while for the crown surface calculation it was 0.74. The accuracy achieved with these open-access databases makes the proposed procedure suitable for acquiring knowledge about the location and state of chestnut plantations as well as for monitoring their evolution.

Description: We introduce a soft computing approach for automatically selecting and combining indices from remote sensing multispectral images that can be used for classification tasks. The proposed approach is based on a Genetic-Programming (GP) framework, a technique successfully used in a wide variety of optimization problems. Through GP, it is possible to learn indices that maximize the separability of samples from two different classes. Once the indices specialized for all the pairs of classes are obtained, they are used in pixelwise classification tasks. We used the GP-based solution to evaluate complex classification problems, such as those that are related to the discrimination of vegetation types within and between tropical biomes. Using time series defined in terms of the learned spectral indices, we show that the GP framework leads to superior results than other indices that are used to discriminate and classify tropical biomes.

Description: The authors are sorry to report that some of the validation data used in their recently published paper [...]

Description: High-resolution terrain models of open-pit mine highwalls and benches are essential in developing new automated slope monitoring systems for operational optimization. This paper presents several contributions to the field of remote sensing in surface mines providing a practical framework for generating high-resolution images using low-trim Unmanned Aerial Vehicles (UAVs). First, a novel mobile application was developed for autonomous drone flights to follow mine terrain and capture high-resolution images of the mine surface. In this article, case study is presented showcasing the ability of developed software to import area terrain, plan the flight accordingly, and finally execute the area mapping mission autonomously. Next, to model the drone’s battery performance, empirical studies were conducted considering various flight scenarios. A multivariate linear regression model for drone power consumption was derived from experimental data. The model has also been validated using data from a test flight. Finally, a genetic algorithm for solving the problem of flight planning and optimization has been employed. The developed power consumption model was used as the fitness function in the genetic algorithm. The designed algorithm was then validated using simulation studies. It is shown that the offered path optimization can reduce the time and energy of high-resolution imagery missions by over 50%. The current work provides a practical framework for stability monitoring of open-pit highwalls while achieving required energy optimization and imagery performance.

Description: We present technical advances and methods to measure effective broadband physical albedo in snowy mountain headwaters using a prototype dual-sensor pyranometer mounted on an Autonomous Aerial Vehicle (an AAV). Our test flights over snowy meadows and forested areas performed well during both clear sky and snowy/windy conditions at an elevation of ~2650 m above mean sea level (MSL). Our AAV-pyranometer platform provided high spatial (m) and temporal resolution (sec) measurements of effective broadband (310–2700 nm) surface albedo. The AAV-based measurements reveal spatially explicit changes in landscape albedo that are not present in concurrent satellite measurements from Landsat and MODIS due to a higher spatial resolution. This AAV capability is needed for validation of satellite snow albedo products, especially over variable montane landscapes at spatial scales of critical importance to hydrological applications. Effectively measuring albedo is important, as annually the seasonal accumulation and melt of mountain snowpack represent a dramatic transformation of Earth’s albedo, which directly affects headwaters’ water and energy cycles.

Description: Current cardiopulmonary activity monitoring is based on contact devices which cannot be used in extreme cases such as premature infants, burnt victims or rescue operations. In order to overcome these limitations, the use of radar technologies emerges as an alternative. This paper aims to enhance the comprehension that non-contact technologies, in particular radar techniques, offer as a monitoring tool. For this purpose, a modified low cost commercial 122 GHz frequency-modulated continuous-wave (FMCW) radar is used to better fit the current application domain. The radar signals obtained are processed using a classic linear filtering algorithm aiming to separate the breathing from the heartbeat component while preserving signals integrity. In a standoff configuration and with different subject orientations, results show that the signal obtained with the radar can be used to extract not only the respiratory and heartbeat rates, but also the heart rate variability (HRV) sequence. Moreover, results evidence the coupling between breathing and heartbeat, also showing that the HRV sequence obtained can identify the respiratory sinus arrhythmia (RSA) effect. Finally, the radar is tested in a simultaneous multi-target scenario, demonstrating its monitoring capabilities in more complex situations. Nevertheless, there are some challenges left to use the system in a real-life monitoring environments, such as the removal of random body movements.

Description: This work presents a regridding procedure applied to the nitrogen dioxide (NO2) tropospheric column data, derived from the Copernicus Sentinel 5 Precursor Tropospheric Monitoring Instrument (S5P/TROPOMI). The regridding has been performed to provide a better comparison with punctual surface observations. It will be demonstrated that TROPOMI NO2 tropospheric column data show improved consistency with in situ surface measurements once the satellite retrievals are scaled to 1 km spatial sampling. A geostatistical technique, i.e., the ordinary kriging, has been applied to improve the spatial distribution of Level 2 TROPOMI NO2 data, which is originally sparse and uneven because of gaps introduced by clouds, to a final spatial, regular, sampling of 1 km × 1 km. The analysis has been performed for two study areas, one in the North and the other in the South of Italy, and for May 2018-April 2020, which also covers the period January 2020-April 2020 of COVID-19 diffusion over the Po Valley. The higher spatial sampling NO2 dataset indicated as Level 3 data, allowed us to explore spatial and seasonal data variability, obtaining better information on NO2 sources. In this respect, it will be shown that NO2 concentrations in March 2020 have likely decreased as a consequence of the lockdown because of COVID-19, although the far warmest winter season ever recorded over Europe in 2020 has favored a general NO2 decrease in comparison to the 2019 winter. Moreover, the comparison between NO2 concentrations related to weekdays and weekend days allowed us to show the strong correlation of NO2 emissions with traffic and industrial activities. To assess the quality and capability of TROPOMI NO2 observations, we have studied their relationship and correlation with in situ NO2 concentrations measured at air quality monitoring stations. We have found that the correlation increases when we pass from Level 2 to Level 3 data, showing the importance of regridding the satellite data. In particular, correlation coefficients of Level 3 data, which range between 0.50–0.90 have been found with higher correlation applying to urban, polluted locations and/or cities.

Description: Convolutional Neural Networks (CNNs) consistently proved state-of-the-art results in image Super-resolution (SR), representing an exceptional opportunity for the remote sensing field to extract further information and knowledge from captured data. However, most of the works published in the literature focused on the Single-image Super-resolution problem so far. At present, satellite-based remote sensing platforms offer huge data availability with high temporal resolution and low spatial resolution. In this context, the presented research proposes a novel residual attention model (RAMS) that efficiently tackles the Multi-image Super-resolution task, simultaneously exploiting spatial and temporal correlations to combine multiple images. We introduce the mechanism of visual feature attention with 3D convolutions in order to obtain an aware data fusion and information extraction of the multiple low-resolution images, transcending limitations of the local region of convolutional operations. Moreover, having multiple inputs with the same scene, our representation learning network makes extensive use of nestled residual connections to let flow redundant low-frequency signals and focus the computation on more important high-frequency components. Extensive experimentation and evaluations against other available solutions, either for Single or Multi-image Super-resolution, demonstrated that the proposed deep learning-based solution can be considered state-of-the-art for Multi-image Super-resolution for remote sensing applications.

Description: The purpose of this study was to develop an optimal estimation model for rainfall rate retrievals using radar reflectivity, thereby gaining an effective grasp of rainfall information for disaster prevention uses. A process was designed for evaluating the optimal retrieval models using various dataset combinations with radar reflectivity and ground meteorological attributes. Various ground meteorological attributes (such as relative humidity, wind speed, precipitation, etc.) were obtained using the land-based weather stations affiliated with Taiwan’s Central Weather Bureau (CWB). This study used nine radar reflectivity provided by the Hualien weather surveillance radar station’s Volume Cover Pattern 21 system. The developed models are built using multiple machine learning algorithms, including linear regression (REG), support vector regression (SVR), and extreme gradient boosting (XGBoost), in addition to the Marshall–Palmer formula (MP). The study examined 14 typhoons that occurred from 2008 to 2017 at Chenggong station in southeast Taiwan, and Lanyu station in the outlying islands, and the top four major rainfall events were designated as test typhoons—Nanmadol (2011), Tembin (2012), Matmo (2014), and Nepartak (2016). The results indicated that for rainfall retrievals, radar reflectivity at a scanning (elevation) angle of 6.0° combined with ground meteorological attributes were the optimal input variables for the Chenggong station, whereas radar reflectivity at an elevation angle of 4.3° combined with ground meteorological attributes were optimal for the Lanyu station. In terms of model performance, XGBoost models had the lowest error index at Chenggong and Lanyu stations compared with MP, REG, and SVR models. XGBoost models at Lanyu station had the highest efficiency coefficient (0.903), and those at Chenggong station had the second highest (0.885). As a result, pairing the combination of optimal radar reflectivity and ground meteorological attributes, as verified by the evaluation process, with a high-efficiency algorithm (XGBoost) can effectively increase the accuracy of rainfall retrieval during typhoons.

Description: The increasing demand for high-resolution hyperspectral images from nano and microsatellites conflicts with the strict bandwidth constraints for downlink transmission. A possible approach to mitigate this problem consists in reducing the amount of data to transmit to ground through on-board processing of hyperspectral images. In this paper, we propose a custom Convolutional Neural Network (CNN) deployed for a nanosatellite payload to select images eligible for transmission to ground, called CloudScout. The latter is installed on the Hyperscout-2, in the frame of the Phisat-1 ESA mission, which exploits a hyperspectral camera to classify cloud-covered images and clear ones. The images transmitted to ground are those that present less than 70% of cloudiness in a frame. We train and test the network against an extracted dataset from the Sentinel-2 mission, which was appropriately pre-processed to emulate the Hyperscout-2 hyperspectral sensor. On the test set we achieve 92% of accuracy with 1% of False Positives (FP). The Phisat-1 mission will start in 2020 and will operate for about 6 months. It represents the first in-orbit demonstration of Deep Neural Network (DNN) for data processing on the edge. The innovation aspect of our work concerns not only cloud detection but in general low power, low latency, and embedded applications. Our work should enable a new era of edge applications and enhance remote sensing applications directly on-board satellite.

Description: Geomagnetic storms are extreme space weather events, which have considerable impacts on the ionosphere and power transmission systems. In this paper, the ionospheric responses to the geomagnetic storm on 22 June 2015, are analyzed from ground-based and satellite-based Global Navigation Satellite System (GNSS) observations as well as observational data of digital ionosondes, and the main physical mechanisms of the ionospheric disturbances observed during the geomagnetic storm are discussed. Salient positive and negative storms are observed from vertical total electron content (VTEC) based on ground-based GNSS observations at different stages of the storm. Combining topside observations of Low-Earth-Orbit (LEO) satellites (GRACE and MetOp satellites) with different orbital altitudes and corresponding ground-based observations, the ionospheric responses above and below the orbits are studied during the storm. To obtain VTEC from the slant TEC between Global Positioning System (GPS) and LEO satellites, we employ a multi-layer mapping function, which can effectively reduce the overall error caused by the single-layer geometric assumption where the horizontal gradient of the ionosphere is not considered. The results show that the topside observations of the GRACE satellite with a lower orbit can intuitively detect the impact caused by the fluctuation of the F2 peak height (hmF2). At the same time, the latitude range corresponding to the peak value of the up-looking VTEC on the event day becomes wider, which is the precursor of the Equatorial Ionization Anomaly (EIA). However, no obvious response is observed in the up-looking VTEC from MetOp satellites with higher orbits, which indicates that the VTEC responses to the geomagnetic storm mainly take place below the orbit of MetOp satellites.

Description: Rapid response mapping of floodwater extents in urbanized areas, while essential for early damage assessment and rescue operations, also presents significant image interpretation challenges. Images from visible band (red–green–blue (RGB)) remote sensors are the most common and cost-effective for real-time applications. Based on an understanding of the differing characteristics of turbid floodwater and urban land surface classes, a robust method was developed and automatized to extract visible floodwater using RGB band digital numbers. The methodology was applied to delineate visible floodwater distribution from very high-resolution aerial image data acquired during the 2013 Calgary flood event. The methodology development involved segment- and pixel-based feature analysis, rule development, automated feature extraction, and result validation processing. The accuracies for the visible floodwater class were above 0.8394% and the overall accuracies were above 0.9668% at both pixel and segment levels for three test sites with diverse urban landscapes.

Description: Ecosystem services refer to the environmental conditions and utilities provided and maintained by ecosystems, which are the basis for the survival and development of human society. The studies on ecosystem services in quantitative assessments, driving mechanisms, and correlation with human well-being, based on remote sensing, have increased in recent years. Various applications of remote sensing in ecosystem services are reported in six papers published in this Special Issue. The major research topics covered by this Special Issue include the multi-method analysis (e.g., linear regression, geographical detector, and geographically weighted regression methodology) of the normalized difference vegetation index (NDVI) to reflect ecosystem structure, the dynamic changing process of ecosystem services, and the determinants, which include a new image-analysis method based on a time series of a biophysical variable and the application of fractional vegetation cover (FVC) to analyze the spatiotemporal relationship between ecosystem structure and function and the comprehensive study on ecosystem function and service based on multi-source remote sensing data. The application of remote sensing data to ecosystem services research has the advantage of monitoring ecological structure and functions at multi-scales. Furthermore, the quantitative calculation of ecosystem services, based on remote sensing, can provide a scientific basis for enhancing land use optimization and sustainable development.

Description: The extensive artificially accelerated erosion of spoil heaps on newly engineered landforms is a key ecological management point requiring better understanding. Soil surface roughness is a crucial factor influencing erosion processes; however, study on spoil heap erosion with a view of surface roughness is lacking. This study investigated the erosion processes and the spatiotemporal variation of surface roughness on spoil heaps, and then, analyzed how the roughness affected the hydrological and sediment yield characteristics. Sequences of four artificial rainstorms with constant rainfall intensity (90 mm/h) were applied to cone-shaped spoil heaps (ground radius 3.5 m, height 2.3 m) of a loess soil containing 30 mass percent rock fragments. The surface elevation was sampled by a laser scanner. For the surface roughness indicators, the root mean square height (rmsh) and the correlation length (cl) increased sharply during the first rainfall event, and in the last three rainfall events, rmsh increased slightly and cl showed a relative decrease. The initial rmsh/cl of the whole slope surface ranged from 0.063 to 0.135, and increased with the rainfall sequence, thus, indicating that the spoil heap surface became rougher. Increasing soil roughness in the rainfall sequence delayed the initial runoff time and increased the runoff yield. The average runoff coefficient of the spoil heaps was 0.658. The average erosion rate of each rainfall event can be simulated by a regression equation of the corresponding average runoff rate and median cl (R-square of 0.816). Soil slumping with an average volume of 0.014 m3 occurred in the first two rainfall events, thus, significantly changing the roughness and peak instant erosion rate. Together, the results revealed the effects of surface roughness on the erosion of spoil heaps and would provide a useful reference for soil loss prediction and control.

Description: Understanding the factors that drive green space composition and richness in heterogeneous urban landscapes is critical for maintaining important ecosystem services and biodiversity. Few studies have been conducted on urban greening and plant diversity at the urban functional unit (UFU) level, although a handful of studies have explored the drivers of greening percentage and its relationships with plant richness in tropical cities. In this study, we conducted field surveys, compiled census and remote sensing data, and performed spatial analyses to investigate the interrelationship between greening percentages, plant diversity, and the socioeconomic variables of different primary and secondary UFUs in the cities of Beijing, Zhanjiang, and Haikou in China. We found that these relationships did not differ significantly between primary and secondary UFUs, and that Parks represented the largest areas of forested land, grassland, and water bodies across all three cities. Moreover, the greening percentages of all UFUs across these three cities were positively correlated with both socioeconomic variables and plant species richness. These relationships can be utilized to guide future green space planning within urban ecosystems.

Description: Eight years after the first record in Italy, Kiwifruit Decline (KD), a destructive disease causing root rot, has already affected more than 25% of the area under kiwifruit cultivation in Italy. Diseased plants are characterised by severe decay of the fine roots and sudden wilting of the canopy, which is only visible after the season’s first period of heat (July–August). The swiftness of symptom appearance prevents correct timing and positioning for sampling of the disease, and is therefore a barrier to aetiological studies. The aim of this study is to test the feasibility of thermal and multispectral imaging for the detection of KD using an unsupervised classifier. Thus, RGB, multispectral and thermal data from a kiwifruit orchard, with healthy and diseased plants, were acquired simultaneously during two consecutive growing seasons (2017–2018) using an Unmanned Aerial Vehicle (UAV) platform. Data reduction was applied to the clipped areas of the multispectral and thermal data from the 2017 survey. Reduced data were then classified with two unsupervised algorithms, a K-means and a hierarchical method. The plant vigour (canopy size and presence/absence of wilted leaves) and the health shifts exhibited by asymptomatic plants between 2017 and 2018 were evaluated from RGB data via expert assessment and used as the ground truth for cluster interpretation. Multispectral data showed a high correlation with plant vigour, while temperature data demonstrated a good potential use in predicting health shifts, especially in highly vigorous plants that were asymptomatic in 2017 and became symptomatic in 2018. The accuracy of plant vigour assessment was above 73% when using multispectral data, while clustering of the temperature data allowed the prediction of disease outbreak one year in advance, with an accuracy of 71%. Based on our results, the unsupervised clustering of remote sensing data could be a reliable tool for the identification of sampling areas, and can greatly improve aetiological studies of this new disease in kiwifruit.

Description: Estimating extreme precipitation events over complex terrain is challenging but crucial for evaluating the performance of climate models for the present climate and expected changes of the climate in the future. New satellites operating in the microwave wavelengths have started to open new opportunities for performing such estimation at adequate temporal and spatial scales and within sensible error limits. This paper illustrates the feasibility and limits of estimating precipitation extremes from satellite data for climatological applications. Using a high-resolution gauge database as ground truth, it was found that global precipitation measurement (GPM) constellation data can provide valuable estimates of extreme precipitation over the southern slopes of the Pyrenees, a region comprising several climates and a very diverse terrain (a challenge for satellite precipitation algorithms). Validation using an object-based quality measure showed reasonable performance, suggesting that GPM estimates can be advantageous reference data for climate model evaluation.

Description: Unsupervised machine learning algorithms (clustering, genetic, and principal component analysis) automate Unmanned Aerial Vehicle (UAV) missions as well as the creation and refinement of iterative 3D photogrammetric models with a next best view (NBV) approach. The novel approach uses Structure-from-Motion (SfM) to achieve convergence to a specified orthomosaic resolution by identifying edges in the point cloud and planning cameras that “view” the holes identified by edges without requiring an initial model. This iterative UAV photogrammetric method successfully runs in various Microsoft AirSim environments. Simulated ground sampling distance (GSD) of models reaches as low as 3.4 cm per pixel, and generally, successive iterations improve resolution. Besides analogous application in simulated environments, a field study of a retired municipal water tank illustrates the practical application and advantages of automated UAV iterative inspection of infrastructure using 63 % fewer photographs than a comparable manual flight with analogous density point clouds obtaining a GSD of less than 3 cm per pixel. Each iteration qualitatively increases resolution according to a logarithmic regression, reduces holes in models, and adds details to model edges.

Description: The estuarine turbidity maximum (ETM) zone occurs in river estuaries due to the effects of tidal dynamics, density-driven residual circulation and deposition/erosion of fine sediments. Even though tropical river estuaries contribute proportionally more to the sediment supply of coastal areas, the ETM in them has been hardly studied. In this study, surface suspended particulate matter (SPM) determined from OLI (Operational Land Imager)-Landsat 8images was used to gain a better understanding of the spatio-temporal dynamics of the ETM of the tropical Maroni estuary (located on the Guianas coast, South America). A method to estimate the remotely-sensed ETM location and its spatiotemporal evolution between 2013 and 2019 was developed. Each ETM was defined from an envelope of normalized SPM values 〉 0.6 calculated from images of the estuary. The results show the influence of the well-marked seasonal river discharge and of tides, especially during the dry season. The ETM is located in the middle estuary during low river-flow conditions, whereas it shifts towards the mouth during high river flow. Neap–spring tidal cycles result in a push of the ETM closer to the mouth under spring-tide conditions or even outside the mouth during the rainy season. An increase in SPM, especially since 2017, coincident with an extension of the ETM, is shown to reflect the periodic influence of mud banks originating from the mouth of the Amazon and migrating along the coast towards the Orinoco (Venezuela). These results demonstrate the advantages of ocean color data in an exploratory study of the spatio-temporal dynamics of the ETM of a tropical estuary, such as that of the Maroni.

Description: As a key component of terrestrial water cycle, evapotranspiration (ET), specifically over the Amazon River basin, is of high scientific significance. However, due to the sparse observation network and relatively short observational period of eddy covariance data, large uncertainties remain in the spatial-temporal characteristics of ET over the Amazon. Recently, a great number of long-term global remotely sensed ET products have been developed to fill the observation gap. However, the reliabilities of these global ET products over the Amazon are unknown. In this study, we assessed the consistency of the magnitude, trend and spatial pattern of Amazon ET among five global remotely sensed ET reconstructions. The magnitudes of these products are similar but the long-term trends from 1982 to 2011 are completely divergent. Validation from the eddy covariance data and water balance method proves a better performance of a product grounded on local measurements, highlighting the importance of local measurements in the ET reconstruction. We also examined four hypotheses dealing with the response of ET to brightening, warming, greening and deforestation, which shows that in general, these ET products respond better to warming and greening than to brightening and deforestation. This large uncertainty highlights the need for future studies focusing on ET issues over the Amazon.

Description: In this study we present the retrieval of the column-averaged dry air mole fraction of carbon dioxide (XCO2) from the TanSat observations using the ACOS (Atmospheric CO2 Observations from Space) algorithm. The XCO2 product has been validated with collocated ground-based measurements from the Total Carbon Column Observing Network (TCCON) for 2 years of TanSat data from 2017 to 2018. Based on the correlation of the XCO2 error over land with goodness of fit in three spectral bands at 0.76, 1.61 and 2.06 μm, we applied an a posteriori bias correction to TanSat retrievals. For overpass averaged results, XCO2 retrievals show a standard deviation (SD) of ~2.45 ppm and a positive bias of ~0.27 ppm compared to collocated TCCON sites. The validation also shows a relatively higher positive bias and variance against TCCON over high-latitude regions. Three cases to evaluate TanSat target mode retrievals are investigated, including one field campaign at Dunhuang with measurements by a greenhouse gas analyzer deployed on an unmanned aerial vehicle and two cases with measurements by a ground-based Fourier-transform spectrometer in Beijing. The results show the retrievals of all footprints, except footprint-6, have relatively low bias (within ~2 ppm). In addition, the orbital XCO2 distributions over Australia and Northeast China between TanSat and the second Orbiting Carbon Observatory (OCO-2) on 20 April 2017 are compared. It shows that the mean XCO2 from TanSat is slightly lower than that of OCO-2 with an average difference of ~0.85 ppm. A reasonable agreement in XCO2 distribution is found over Australia and Northeast China between TanSat and OCO-2.

Description: Sugarcane (complex hybrids of Saccharum spp., C4 plant) croplands provide cane stalk feedstock for sugar and biofuel (ethanol) production. It is critical for us to analyze the phenology and gross primary production (GPP) of sugarcane croplands, which would help us to better understand and monitor the sugarcane growing condition and the carbon cycle. In this study, we combined the data from two sugarcane EC flux tower sites in Brazil and the USA, images from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor, and data-driven models to study the phenology and GPP of sugarcane croplands. The seasonal dynamics of climate, vegetation indices from MODIS images, and GPP from two sugarcane flux tower sites (GPPEC) reveal the temporal consistency in sugarcane phenology (crop calendar: green-up dates and harvesting dates) as estimated by the vegetation indices and GPPEC data. The Land Surface Water Index (LSWI) is shown to be useful to delineate the phenology of sugarcane croplands. The relationship between the sugarcane GPPEC and the Enhanced Vegetation Index (EVI) is stronger than the relationship between the GPPEC and the Normalized Difference Vegetation Index (NDVI). We ran the Vegetation Photosynthesis Model (VPM), which uses the light use efficiency (LUE) concept and is driven by climate data and MODIS images, to estimate the daily GPP at the two sugarcane sites (GPPVPM). The seasonal dynamics of the GPPVPM and GPPEC at the two sites agreed reasonably well with each other, which indicates that VPM is a powerful tool for estimating the GPP of sugarcane croplands in Brazil and the USA. This study clearly highlights the potential of combining eddy covariance technology, satellite-based remote sensing technology, and data-driven models for better understanding and monitoring the phenology and GPP of sugarcane croplands under different climate and management practices.

Description: Detailed analysis of five gas emission craters (GEC) found in the north of West Siberia is presented. Remote sensing data used in the study is verified by field surveys. Previous studies show that all of the GECs were preceded by mounds 2 to 6 m high and 20 to 55 m in diameter. GECs initially were 20–25 m in diameter, which increased in the first years of their existence. GECs are found in various environmental (shrublands or moss-grass tundra) and geomorphic (river valley, terrace, slopes) conditions. The objective of the paper is to identify common and differing geomorphologic and environmental characteristics of all the five GEC, and their mound-predecessors. The study is based on a compilation of DSMs before and after the GEC formation using very high-resolution satellite imagery stereo pairs compared to ArcticDEM project data. Diversity of terrain and environmental settings along with rather a narrow range of GEC and mound-predecessor morphometric parameters allows concluding that the mechanism of GEC formation is most likely similar for all the GEC and is controlled rather by internal geologic and cryolithologic structure than by any surface properties.

Description: The availability of large amounts of Sentinel-2 data has been a trigger for its increasing exploitation in various types of applications. It is, therefore, of importance to understand the limits above which these data still guarantee a meaningful outcome. This paper proposes a new method to quantify and specify restrictions of the Sentinel-2 imagery in the context of checks by monitoring, a newly introduced control approach within the European Common Agriculture Policy framework. The method consists of a comparison of normalized difference vegetation index (NDVI) time series constructed from data of different spatial resolution to estimate the performance and limits of the coarser one. Using similarity assessment of Sentinel-2 (10 m pixel size) and PlanetScope (3 m pixel size) NDVI time series, it was estimated that for 10% out of 867 fields less than 0.5 ha in size, Sentinel-2 data did not provide reliable evidence of the activity or state of the agriculture field over a given timeframe. Statistical analysis revealed that the number of clean or full pixels and the proportion of pixels lost after an application of a 5-m (1/2 pixel) negative buffer are the geospatial parameters of the field that have the highest influence on the ability of the Sentinel-2 data to qualify the field’s state in time. We specified the following limiting criteria: at least 8 full pixels inside a border and less than 60% of pixels lost. It was concluded that compliance with the criteria still assures a high level of extracted information reliability. Our research proved the promising potential, which was higher than anticipated, of Sentinel-2 data for the continuous state assessment of small fields. The method could be applied to other sensors and indicators.

Description: Shrub-dominated ecosystems support biodiversity and play an important storage role in the global carbon cycle. However, it is challenging to characterize biophysical properties of low-stature vegetation like shrubs from conventional ground-based or remotely sensed data. We used spectral and structural variables derived from high-resolution unmanned aerial system (UAS) imagery to estimate the aboveground biomass of shrubs in the Betula and Salix genera in a montane meadow in Banff National Park, Canada using an area-based approach. In single-variable linear regression models, visible light (RGB) indices outperformed multispectral or structural data. A linear model based on the red ratio vegetation index (VI) accumulated over shrub area could model biomass (calibration R2 = 0.888; validation R2 = 0.774) nearly as well as the top multivariate linear regression models (calibration R2 = 0.896; validation R2 〉 0.750), which combined an accumulated RGB VI with a multispectral metric. The excellent performance of accumulated RGB VIs represents a novel approach to fine-scale vegetation biomass estimation, fusing spectral and spatial information into a single parsimonious metric that rivals the performance of more complex multivariate models. Methods developed in this study will be relevant to researchers interested in estimating fine-scale shrub aboveground biomass within a range of ecosystems.

Description: Despite providing vital ecosystem services, wetlands are increasingly threatened across the globe by both anthropogenic activities and natural processes. Synthetic aperture radar (SAR) has emerged as a promising tool for rapid and accurate monitoring of wetland extent and type. By acquiring information on the roughness and moisture content of the surface, SAR offers unique potential for wetland monitoring. However, there are still challenges in applying SAR for mapping complex wetland environments. The backscattering similarity of different wetland classes is one of the challenges. Choosing the appropriate SAR specifications (incidence angle, frequency and polarization), based on the wetland type, is also a subject of debate and should be investigated more thoroughly. The geometric distortion of SAR imagery and loss of coherency are other remaining challenges in applying SAR and its processing techniques for wetland studies. Hence, this study provides a systematic meta-analysis based on compilation and analysis of indexed research studies that used SAR for wetland monitoring. This meta-analysis reviewed 172 papers and documented an upward trend in usage of SAR data, increasing usage of multi-sensor data, increasing integration of C- and L- bands over other configurations and higher classification accuracy with multi-frequency and multi-polarized SAR data. The highest number of wetland research studies using SAR data came from the USA, Canada and China. This meta-analysis highlighted the current challenges and solutions for wetland monitoring using SAR sensors.

Description: The world’s most extensive tropical peatlands occur in the Cuvette Centrale depression in the Congo Basin, which stores 30.6 petagrams of carbon (95% CI, 6.3–46.8). Improving our understanding of the genesis, development and functioning of these under-studied peatlands requires knowledge of their topography and, in particular, whether the peat surface is domed, as this implies a rain-fed system. Here we use a laser altimeter mounted on an unmanned airborne vehicle (UAV) to measure peat surface elevation along two transects at the edges of a peatland, in the northern Republic of Congo, to centimetre accuracy and compare the results with an analysis of nearby satellite LiDAR data (ICESat and ICESat-2). The LiDAR elevations on both transects show an upward slope from the peatland edge, suggesting a surface elevation peak of around 1.8 m over ~20 km. While modest, this domed shape is consistent with the peatland being rainfed. In-situ peat depth measurements and our LiDAR results indicate that this peatland likely formed at least 10,000 years BP in a large shallow basin ~40 km wide and ~3 m deep.

Description: The Direct Broadcast Network (DBNet) provides near-real-time delivery of low-earth-orbiting (LEO) meteorological satellites to operational numerical weather prediction (NWP) systems that need short data cut-off times to allow for the assimilation of the most recent satellite measurements. The NWP model requires timely delivery of observations including atmospheric temperature, humidity, and surface wind vectors. The World Meteorological Organization (WMO) Space Program (WSP) recommends the data latency of no more than 20 min for the satellite measurements. Currently, not all DBNet stations are delivering satellite data within the 20-min time frame. In this study, the forecast impact of the observations of LEO satellite sounders with data latency of 20 min or less was evaluated using the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS). Reducing the data latency up to 5 min increases the number of LEO infrared (IR) and microwave (MW) sounder observations delivered to the NCEP GFS data assimilation system by more than 20%. Overall, this study demonstrates a positive impact on the global weather forecasts when the IR and MW sounder data are delivered by 20 min anywhere in the world. Additional forecast benefits are not obvious for shorter data latency. Results from this study support the WSP recommendation of 20–minute data latency.

Description: Drought is a natural phenomenon that originates from the absence of precipitation over a certain period and is capable of causing damage to societal development. With the advent of orbital remote sensing, rainfall estimates from satellites have appeared as viable alternatives to monitor natural hazards in ungauged basins and complex areas of the world; however, the accuracies of these orbital products still need to be verified. Thus, this work aims to evaluate the performance of Tropical Rainfall Measuring Mission (TRMM) satellite rainfall estimates in monitoring the spatiotemporal behavior of droughts at multiple temporal scales over Paraíba State based on the standardized precipitation index (SPI) over 20 years (1998–2017). For this purpose, rainfall data from 78 rain gauges and 187 equally spaced TRMM cell grids throughout the region are used, and accuracy analyses are performed at the single-gauge level and in four mesoregions at eight different time scales based on 11 statistical metrics calculations divided into three different categories. The results show that in the mesoregions close to the coast, the satellite-based product is less accurate in capturing the drought behavior regardless of the evaluated statistical metrics. At the temporal scale, the TRMM is more accurate in identifying the pattern of medium-term droughts; however, there is considerable spatial variation in the accuracy of the product depending on the performance index. Therefore, it is concluded that rainfall estimates from the TRMM satellite are a valuable source of data to identify drought behavior in a large part of Paraíba State at different time scales, and further multidisciplinary studies should be conducted to monitor these phenomena more accurately based on satellite data.

Description: In this paper, we applied the re-analysis data cobe-SST (cobe-sea surface temperature) and Global Land Data Assimilation System (GLDAS) surface soil moisture (SM) data from 1961 to 2011 by using regional correlation analysis and time series causality analysis to trace annual variations in and identify the abnormal relationship of sea surface temperature (SST) in the eastern China Sea and SM in eastern China (EC). We also used satellite Moderate Resolution Imaging Spectroradiometer (MODIS) SST and AMSR-E SM data to examine the correlation of SST and SM in EC from 2004–2009. The results show that the SST in the eastern China Sea has experienced a warming trend since 1987, whereas the SM in EC has shown a drying trend since 1978. Before 1967 and after 1997, SST and SM changed during opposite phases, whereas from 1967 to 1997 they changed during the same phase. The differences between them may result from the abnormal summer precipitation causing abnormal SM. According to the regional correlation analysis, SST of the East China Sea is significantly related to SM in the southeast coastal area, and temporal sequence causality analysis shows that SST is correlated with and has higher influence on SM than vice versa. SM during spring and autumn shows a similar correlation with SST during the four seasons, so that SM in spring and autumn is positively correlated with SST in autumn and negatively correlated with SST in other seasons. SM in summer and winter correlated with SST in the four seasons, contradicting the foregoing conclusions. All these findings indicate that the thermodynamic state of the eastern China Sea has affected SM in EC.

Description: Archaeologists engaging with Airborne Laser Scanning (ALS) data rely heavily on manual inspection of various derived visualizations. However, manual inspection of ALS data is extremely time-consuming and as such presents a major bottleneck in the data analysis workflow. We have therefore set out to learn and test a deep neural network model for classifying from previously manually annotated ancient Maya structures of the Chactún archaeological site in Campeche, Mexico. We considered several variations of the VGG-19 Convolutional Neural Network (CNN) to solve the task of classifying visualized example structures from previously manually annotated ALS images of man-made aguadas, buildings and platforms, as well as images of surrounding terrain (four classes and over 12,000 anthropogenic structures). We investigated how various parameters impact model performance, using: (a) six different visualization blends, (b) two different edge buffer sizes, (c) additional data augmentation and (d) architectures with different numbers of untrainable, frozen layers at the beginning of the network. Many of the models learned under the different scenarios exceeded the overall classification accuracy of 95%. Using overall accuracy, terrain precision and recall (detection rate) per class of anthropogenic structure as criteria, we selected visualization with slope, sky-view factor and positive openness in separate bands; image samples with a two-pixels edge buffer; Keras data augmentation; and five frozen layers as the optimal combination of building blocks for learning our CNN model.