Guest Editors help invite many high-quality papers for Remote Sensing[...] Full article

Peer review is an essential part of the publication process, ensuring that Remote Sensing maintains the high standard of its published papers[...] Full article

Triggered by earthquakes, rainfall, or anthropogenic activities, landslides represent widespread and problematic geohazards worldwide. In recent years, multiple remote sensing techniques, including synthetic aperture radar, optical, and light detection and ranging measurements from spaceborne, airborne, and ground-based platforms, have been widely applied for the analysis of landslide processes. Current techniques include landslide detection, inventory mapping, surface deformation monitoring, trigger factor analysis and mechanism inversion. In addition, landslide susceptibility modelling, hazard assessment, and risk evaluation can be further analyzed using a synergic fusion of multiple remote sensing data and other factors affecting landslides. We summarize the 19 articles collected in this special issue of Remote Sensing of Landslide, in the terms of data, methods and applications used in the papers. Full article

The most commonly used approach to estimate soil variables from remote-sensing data entails time-consuming and expensive data collection including chemical and physical laboratory analysis. Large spectral libraries could be exploited to decrease the effort of soil variable estimation and obtain more widely applicable models. We investigated the feasibility of a new approach, referred to as bottom-up, to provide soil organic carbon (SOC) maps of bare cropland fields over a large area without recourse to chemical analyses, employing both the pan-European topsoil database from the Land Use/Cover Area frame statistical Survey (LUCAS) and Airborne Prism Experiment (APEX) hyperspectral airborne data. This approach was tested in two areas having different soil characteristics: the loam belt in Belgium, and the Gutland–Oesling region in Luxembourg. Partial least square regression (PLSR) models were used in each study area to estimate SOC content, using both bottom-up and traditional approaches. The PLSR model’s accuracy was tested on an independent validation dataset. Both approaches provide SOC maps having a satisfactory level of accuracy (RMSE = 1.5–4.9 g·kg⁻¹; ratio of performance to deviation (RPD) = 1.4–1.7) and the inter-comparison did not show differences in terms of RMSE and RPD either in the loam belt or in Luxembourg. Thus, the bottom-up approach based on APEX data provided high-resolution SOC maps over two large areas showing the within- and between-field SOC variability. Full article

The issue for the validation of land surface remote sensing albedo products over rugged terrain is the scale effects between the reference albedo measurements and coarse scale albedo products, which is caused by the complex topography. This paper illustrates a multi-scale validation strategy specified for coarse scale albedo validation over rugged terrain. A Mountain-Radiation-Transfer-based (MRT-based) albedo upscaling model was proposed in the process of multi-scale validation strategy for aggregating fine scale albedo to coarse scale. The simulated data of both the reference coarse scale albedo and fine scale albedo were used to assess the performance and uncertainties of the MRT-based albedo upscaling model. The results showed that the MRT-based model could reflect the albedo scale effects over rugged terrain and provided a robust solution for albedo upscaling from fine scale to coarse scale with different mean slopes and different solar zenith angles. The upscaled coarse scale albedos had the great agreements with the simulated coarse scale albedo with a Root-Mean-Square-Error (RMSE) of 0.0029 and 0.0017 for black sky albedo (BSA) and white sky albedo (WSA), respectively. Then the MRT-based model was preliminarily applied for the assessment of daily MODerate Resolution Imaging Spectroradiometer (MODIS) Albedo Collection V006 products (MCD43A3 C6) over rugged terrain. Results showed that the MRT-based model was effective and suitable for conducting the validation of MODIS albedo products over rugged terrain. In this research area, it was shown that the MCD43A3 C6 products with full inversion algorithm, were generally in agreement with the aggregated coarse scale reference albedos over rugged terrain in the Heihe River Basin, with the BSA RMSE of 0.0305 and WSA RMSE of 0.0321, respectively, which were slightly higher than those over flat terrain. Full article

Suspended particulate matter (SPM) concentrations ([SPM]) in the Yangtze estuary, which has third-order bifurcations and four outlets, exhibit large spatial and temporal variations. Studying the characteristics of these variations in [SPM] is important for understanding sediment transport and pollutant diffusion in the estuary as well as for the construction of port and estuarine engineering structures. The 1-h revisit frequency of the Geostationary Ocean Color Imager (GOCI) sensor and the 30-m spatial resolution of the Landsat 8 Operational Land Imager (L8/OLI) provide a new opportunity to study the large spatial and temporal variations in the [SPM] in the Yangtze estuary. In this study, [SPM] images with a temporal resolution of 1 h and a spatial resolution of 30 m are generated through the product-level fusion of [SPM] data derived from L8/OLI and GOCI images using the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The results show that the details and accuracy of the spatial and temporal variations are maintained well in the [SPM] images that are predicted based on the fused images. Compared to the [SPM] observations at fixed field stations, the mean relative error (MRE) of the predicted SPM is 17.7%, which is lower than that of the GOCI-derived [SPM] (27.5%). In addition, thanks to the derived high-resolution [SPM] with high spatiotemporal dynamic changes, both natural phenomena (dynamic variation of the maximum turbid zone) and human engineering changes leading to the dynamic variability of SPM in the channel are observed. Full article

Cropland abandonment is globally widespread and has strong repercussions for regional food security and the environment. Statistics suggest that one of the hotspots of abandoned cropland is located in the drylands of the Aral Sea Basin (ASB), which covers parts of post-Soviet Central Asia, Afghanistan and Iran. To date, the exact spatial and temporal extents of abandoned cropland remain unclear, which hampers land-use planning. Abandoned land is a potentially valuable resource for alternative land uses. Here, we mapped the abandoned cropland in the drylands of the ASB with a time series of the Normalized Difference Vegetation Index (NDVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) from 2003–2016. To overcome the restricted ability of a single classifier to accurately map land-use classes across large areas and agro-environmental gradients, “stratum-specific” classifiers were calibrated and classification results were fused based on a locally weighted decision fusion approach. Next, the agro-ecological suitability of abandoned cropland areas was evaluated. The stratum-specific classification approach yielded an overall accuracy of 0.879, which was significantly more accurate ( $p$ < 0.05) than a “global” classification without stratification, which had an accuracy of 0.811. In 2016, the classification results showed that 13% (1.15 Mha) of the observed irrigated cropland in the ASB was idle (abandoned). Cropland abandonment occurred mostly in the Amudarya and Syrdarya downstream regions and was associated with degraded land and areas prone to water stress. Despite the almost twofold population growth and increasing food demand in the ASB area from 1990 to 2016, abandoned cropland was also located in areas with high suitability for farming. The map of abandoned cropland areas provides a novel basis for assessing the causes leading to abandoned cropland in the ASB. This contributes to assessing the suitability of abandoned cropland for food or bioenergy production, carbon storage, or assessing the environmental trade-offs and social constraints of recultivation. Full article

The constellation of Landsat-8 and Sentinel-2 optical satellites offers opportunities for a wide range of Earth Observation (EO) applications and scientific studies in Earth sciences mainly related to geohazards. The multi-temporal co-registration accuracy of images provided by both missions is, however, currently not fully satisfactory for change detection, time-series analysis and in particular Earth surface motion measurements. The objective of this work is the development, implementation and test of an automatic processing chain for correcting co-registration artefacts targeting accurate alignment of Sentinel-2 and Landsat-8 imagery for time series analysis. The method relies on dense sub-pixel offset measurements and robust statistics to correct for systematic offsets and striping artefacts. Experimental evaluation at sites with diverse environmental settings is conducted to evaluate the efficiency of the processing chain in comparison with previously proposed routines. The experimental evaluation suggests lower residual offsets than existing methods ranging between $RMS{E}_{xy}$ = 2.30 and 2.91 m remaining stable for longer time series. A first case study demonstrates the utility of the processor for the monitoring of continuously active landslides. A second case study demonstrates the use of the processor for measuring co-seismic surface displacements indicating an accuracy of 1/5 th of a pixel after corrections and 1/10th of a pixel after calibration with ground measurements. The implemented processing chain is available as an open source tool to support a better exploitation of the growing archives of Sentinel-2 and Landsat-8. Full article

The reliability of airborne light detection and ranging (LiDAR) for delineating individual trees and estimating aboveground biomass (AGB) has been proven in a diverse range of ecosystems, but can be difficult and costly to commission. Point clouds derived from structure from motion (SfM) matching techniques obtained from unmanned aerial systems (UAS) could be a feasible low-cost alternative to airborne LiDAR scanning for canopy parameter retrieval. This study assesses the extent to which SfM three-dimensional (3D) point clouds—obtained from a light-weight mini-UAS quadcopter with an inexpensive consumer action GoPro camera—can efficiently and effectively detect individual trees, measure tree heights, and provide AGB estimates in Australian tropical savannas. Two well-established canopy maxima and watershed segmentation tree detection algorithms were tested on canopy height models (CHM) derived from SfM imagery. The influence of CHM spatial resolution on tree detection accuracy was analysed, and the results were validated against existing high-resolution airborne LiDAR data. We found that the canopy maxima and watershed segmentation routines produced similar tree detection rates (~70%) for dominant and co-dominant trees, but yielded low detection rates (<35%) for suppressed and small trees due to poor representativeness in point clouds and overstory occlusion. Although airborne LiDAR provides higher tree detection rates and more accurate estimates of tree heights, we found SfM image matching to be an adequate low-cost alternative for the detection of dominant and co-dominant tree stands. Full article

The Geostationary Environment Monitoring Spectrometer (GEMS) is scheduled to be in orbit in 2019 onboard the GEO-KOMPSAT 2B satellite and will continuously monitor air quality over Asia. The GEMS will make measurements in the UV spectrum (300–500 nm) with 0.6 nm resolution. In this study, an algorithm is developed to retrieve aerosol optical properties from UV-visible measurements for the future satellite instrument and is tested using 3 years of existing OMI L1B data. This algorithm provides aerosol optical depth (AOD), single scattering albedo (SSA) and aerosol layer height (ALH) using an optimized estimation method. The retrieved AOD shows good correlation with Aerosol Robotic Network (AERONET) AOD with correlation coefficients of 0.83, 0.73 and 0.80 for heavy-absorbing fine (HAF) particles, dust and non-absorbing (NA) particles, respectively. However, regression tests indicate underestimation and overestimation of HAF and NA AOD, respectively. In comparison with AOD from the OMI/Aura Near-UV Aerosol Optical Depth and Single Scattering Albedo 1-orbit L2 Swath 13 km × 24 km V003 (OMAERUV) algorithm, the retrieved AOD has a correlation coefficient of 0.86 and linear regression equation, AOD_GEMS = 1.18AOD_OMAERUV + 0.09. An uncertainty test based on a reference method, which estimates retrieval error by applying the algorithm to simulated radiance data, revealed that assumptions in the spectral dependency of aerosol absorptivity in the UV cause significant errors in aerosol property retrieval, particularly the SSA retrieval. Consequently, retrieved SSAs did not show good correlation with AERONET values. The ALH results were qualitatively compared with the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) products and were found to be well correlated for highly absorbing aerosols. The difference between the attenuated-backscatter-weighted height from CALIOP and retrieved ALH were mostly closed to zero when the retrieved AOD is higher than 0.8 and SSA is lower than 0.93. Although retrieval accuracy was not significantly improved, the simultaneous consistent retrieval of AOD, SSA and ALH alone demonstrates the value of this stand-alone algorithm, given their nature for error using other methods. The use of these properties as input parameters for the air mass factor calculation is expected to improve the retrieval of other trace gases over Asia. Full article

The India–France SARAL/AltiKa mission is the first Ka-band altimetric mission dedicated primarily to oceanography. The mission objectives were firstly the observation of the oceanic mesoscales but also global and regional sea level monitoring, including the coastal zone, data assimilation, and operational oceanography. SARAL/AltiKa proved also to be a great opportunity for inland waters applications, for observing ice sheet or icebergs, as well as for geodetic investigations. The mission ended its nominal phase after three years in orbit and began a new phase (drifting orbit) in July 2016. The objective of this paper is to highlight some of the most remarkable achievements of the SARAL/AltiKa mission in terms of scientific applications. Compared to the standard Ku-band altimetry measurements, the Ka-band provides substantial improvements in terms of spatial resolution and data accuracy. We show here that this leads to remarkable advances in terms of observation of the mesoscale and coastal ocean, waves, river water levels, ice sheets, icebergs, fine scale bathymetry features as well as for the many related applications. Full article

Mountain ecosystems are among the most fragile environments on Earth. The availability of timely updated information on forest 3D structure would improve our understanding of the dynamic and impact of recent disturbance and regeneration events including fire, insect damage, and drought. Airborne lidar is a critical tool for monitoring forest change at high resolution but it has been little used for this purpose due to the scarcity of long-term time-series of measurements over a common region. Here, we investigate the reliability of on-going, multi-year lidar observations from the NASA-JPL Airborne Snow Observatory (ASO) to characterize forest 3D structure at a fine spatial scale. In this study, weekly ASO measurements collected at ~1 pt/m², primarily acquired to quantify snow volume and dynamics, are coherently merged to produce high-resolution point clouds ( $~$ 12 pt/m²) that better describe forest structure. The merging methodology addresses the spatial bias in multi-temporal data due to uncertainties in platform trajectory and motion by collecting tie objects from isolated tree crown apexes in the lidar data. The tie objects locations are assigned to the centroid of multi-temporal lidar points to fuse and optimize the location of multiple measurements without the need for ancillary data or GPS control points. We apply the methodology to ASO lidar acquisitions over the Tuolumne River Basin in the Sierra Nevada, California, during the 2014 snow monitoring campaign and provide assessment of the fidelity of the fused point clouds for forest mountain ecosystem studies. The availability of ASO measurements that currently span 2013–2017 enable annual forest monitoring of important vegetated ecosystems that currently face ecological threads of great significance such as the Sierra Nevada (California) and Olympic National Forest (Washington). Full article

Cape Trafalgar has been highlighted as a hotspot of high chlorophyll concentrations, as well as a source of biomass for the Alborán Sea. It is located in an unique geographical framework between the Gulf of Cádiz (GoC), which is dominated by long-term seasonal variability, and the Strait of Gibraltar, which is mainly governed by short-term tidal variability. Furthermore, here bathymetry plays an important role in the upwelling of nutrient-rich waters. In order to study the spatial and temporal variability of chlorophyll-a in this region, 10 years of ocean colour observations using the MEdium Resolution Imaging Spectrometer (MERIS) were analysed through different approaches. An empirical orthogonal function decomposition distinguished two coastal zones with opposing phases that were analysed by wavelet methods in order to identify their temporal variability. In addition, to better understand the physical–biological interaction in these zones, the co-variation between chlorophyll-a and different environmental variables (wind, river discharge, and tidal current) was analysed. Zone 1, located on the GoC continental shelf, was characterised by a seasonal variability weakened by the influence of other environmental variables. Meanwhile, Zone 2, which represented the dynamics in Cape Trafalgar but did not show any clear pattern of variability, was strongly correlated with tidal current whose variability was probably determined by other drivers. Full article

Stable Water Isotopologues (SWIs) are important diagnostic tracers for understanding processes in the atmosphere and the global hydrological cycle. Using eight years (2002–2009) of retrievals from Odin/SMR (Sub-Millimetre Radiometer), the global climatological features of three SWIs, H₂¹⁶O, HDO and H₂¹⁸O, the isotopic composition δD and δ¹⁸O in the stratosphere are analysed for the first time. Spatially, SWIs are found to increase with altitude due to stratospheric methane oxidation. In the tropics, highly depleted SWIs in the lower stratosphere indicate the effect of dehydration when the air comes through the cold tropopause, while, at higher latitudes, more enriched SWIs in the upper stratosphere during summer are produced and transported to the other hemisphere via the Brewer–Dobson circulation. Furthermore, we found that more H₂¹⁶O is produced over summer Northern Hemisphere and more HDO is produced over summer Southern Hemisphere. Temporally, a tape recorder in H₂¹⁶O is observed in the lower tropical stratosphere, in addition to a pronounced downward propagating seasonal signal in SWIs from the upper to the lower stratosphere over the polar regions. These observed features in SWIs are further compared to SWI-enabled model outputs. This helped to identify possible causes of model deficiencies in reproducing main stratospheric features. For instance, choosing a better advection scheme and including methane oxidation process in a specific model immediately capture the main features of stratospheric water vapor. The representation of other features, such as the observed inter-hemispheric difference of isotopic component, is also discussed. Full article

In this study, we analyze interactions in lake and lake catchment systems of a continuous permafrost area. We assessed colored dissolved organic matter (CDOM) absorption at 440 nm (a(440)_CDOM) and absorption slope (S_300–500) in lakes using field sampling and optical remote sensing data for an area of 350 km² in Central Yamal, Siberia. Applying a CDOM algorithm (ratio of green and red band reflectance) for two high spatial resolution multispectral GeoEye-1 and Worldview-2 satellite images, we were able to extrapolate the a(λ)_CDOM data from 18 lakes sampled in the field to 356 lakes in the study area (model R² = 0.79). Values of a(440)_CDOM in 356 lakes varied from 0.48 to 8.35 m⁻¹ with a median of 1.43 m⁻¹. This a(λ)_CDOM dataset was used to relate lake CDOM to 17 lake and lake catchment parameters derived from optical and radar remote sensing data and from digital elevation model analysis in order to establish the parameters controlling CDOM in lakes on the Yamal Peninsula. Regression tree model and boosted regression tree analysis showed that the activity of cryogenic processes (thermocirques) in the lake shores and lake water level were the two most important controls, explaining 48.4% and 28.4% of lake CDOM, respectively (R² = 0.61). Activation of thermocirques led to a large input of terrestrial organic matter and sediments from catchments and thawed permafrost to lakes (n = 15, mean a(440)_CDOM = 5.3 m⁻¹). Large lakes on the floodplain with a connection to Mordy-Yakha River received more CDOM (n = 7, mean a(440)_CDOM = 3.8 m⁻¹) compared to lakes located on higher terraces. Full article

Precipitation is a key aspect of the climate system. In this paper, the dependability of five satellite precipitation products (TRMM [Tropical Rainfall Measuring Mission] 3BV42, PERSIANN [Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks] CDR, GSMaP [Global Satellite Mapping of Precipitation] RENALYSIS, CMORPH [Climate Prediction Center’s morphing technique] BLD and CMORPH_RAW) were compared with in situ measurements over China for the period of 2005 to 2013. To completely evaluate these precipitation products, the annual, seasonal and monthly precipitation averages were calculated. Overall, the Huaihe River and Qinlin mountains are shown to have heavy precipitation to the southeast and lighter precipitation to the northwest. The comparison results indicate that Gauge correction (CMORPH_BLD) improves the quality of the original satellite products (CMORPH_RAW), resulting in the higher correlation coefficient (CC), the low relative bias (BIAS) and root mean square error (RMSE). Over China, the GSMaP_RENALYSIS outperforms other products and shows the highest CC (0.91) and lowest RMSE (0.85 mm/day) and all products except for PERSIANN_CDR exhibit underestimation. GSMaP_RENALYSIS gives the highest of probability of detection (81%), critical success index (63%) and lowest false alarm ratio (36%) while TRMM3BV42 gives the highest of frequency bias index (1.00). Over Tibetan Plateau, CMORPH_RAW demonstrates the poorest performance with the biggest BIAS (4.2 mm/month) and lowest CC (0.22) in December 2013. GSMaP_RENALYSIS displays quite consistent with in situ measurements in summer. However, GSMaP_RENALYSIS and CMORPH_RAW underestimate precipitation over South China. CMORPH_BLD and TRMM3BV42 show consistent with high CC (>0.8) but relatively large RMSE in summer. Full article

A high-accuracy Global Ionosphere Model (GIM) is significant for precise positioning and navigating with the Global Navigation Satellite System (GNSS), as well as space weather applications. To obtain a precise GIM, it is critical to take both the ionospheric observable and mathematical model into consideration. In this contribution, the undifferenced ambiguity-fixed carrier-phase ionospheric observable is first determined from a global distribution of permanent receivers. Accuracy assessment with a co-located station experiment shows that the observational errors affecting the ambiguity-fixed carrier-phase ionospheric observables range from 0.10 to 0.35 Total Electron Content Units (TECUs, where 1 TECU $={10}^{16}{\mathrm{e}}^{-}/{\mathrm{m}}^2$ and corresponds to 0.162 m on the Global Positioning System, GPS L1 frequency), indicating that the ambiguity-fixed carrier-phase ionospheric observable is over one order of magnitude more accurate than the carrier-phase leveled-code one (from 1.21 to 3.77 TECUs). Second, to better model the structure of the ionosphere, a two-layer GIM has been built based on the above carrier-phase observable. Preliminary global accuracy evaluation demonstrates that the accuracy of the two-layer GIM is below 1 TECU and about 2 TECUs during low and high solar activity periods. Third, the single-frequency point positioning experiment is adopted to test the ionosphere mitigation effects of the GIMs. Positioning results demonstrate that the single-frequency positioning accuracy can be improved by more than 30% using the undifferenced ambiguity-fixed ionospheric observable-derived two-layer GIM, compared with that using the carrier-phase leveled-code ionospheric observable-based single-layer GIM. Full article

Current methods for retrieving SWE (snow water equivalent) from space rely on passive microwave sensors. Observations are limited by poor spatial resolution, ambiguities related to separation of snow microstructural properties from the total snow mass, and signal saturation when snow is deep (~>80 cm). The use of SAR (Synthetic Aperture Radar) at suitable frequencies has been suggested as a potential observation method to overcome the coarse resolution of passive microwave sensors. Nevertheless, suitable sensors operating from space are, up to now, unavailable. Active microwave retrievals suffer, however, from the same difficulties as the passive case in separating impacts of scattering efficiency from those of snow mass. In this study, we explore the potential of applying active (radar) and passive (radiometer) microwave observations in tandem, by using a dataset of co-incident tower-based active and passive microwave observations and detailed in situ data from a test site in Northern Finland. The dataset spans four winter seasons with daily coverage. In order to quantify the temporal variability of snow microstructure, we derive an effective correlation length for the snowpack (treated as a single layer), which matches the simulated microwave response of a semi-empirical radiative transfer model to observations. This effective parameter is derived from radiometer and radar observations at different frequencies and frequency combinations (10.2, 13.3 and 16.7 GHz for radar; 10.65, 18.7 and 37 GHz for radiometer). Under dry snow conditions, correlations are found between the effective correlation length retrieved from active and passive measurements. Consequently, the derived effective correlation length from passive microwave observations is applied to parameterize the retrieval of SWE using radar, improving retrieval skill compared to a case with no prior knowledge of snow-scattering efficiency. The same concept can be applied to future radar satellite mission concepts focused on retrieving SWE, exploiting existing methods for retrieval of snow microstructural parameters, as employed within the ESA (European Space Agency) GlobSnow SWE product. Using radar alone, a seasonally optimized value of effective correlation length to parameterize retrievals of SWE was sufficient to provide an accuracy of <25 mm (unbiased) Root-Mean Square Error using certain frequency combinations. A temporally dynamic value, derived from e.g., physical snow models, is necessary to further improve retrieval skill, in particular for snow regimes with larger temporal variability in snow microstructure and a more pronounced layered structure. Full article

This study analyzes the changes in drought patterns in China’s major crop-growing areas over the past 60 years. The analysis was done using both weather station data and Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) rainfall data to calculate the Standardized Precipitation Index (SPI). The results showed that the occurrences of extreme drought were the most serious in recent years in the Southwest China and Sichuan crop-growing areas. The Yangtze River (MLRY) and South China crop-growing areas experienced extreme droughts during 1960–1980, whereas the Northeast China and Huang–Huai–Hai crop-growing areas experienced extreme droughts around 2003. The analysis showed that the SPIs calculated by TRMM data at time scales of one, three, and six months were reliable for monitoring drought in the study regions, but for 12 months, the SPIs calculated by gauge and TRMM data showed less consistency. The analysis of the spatial distribution of droughts over the past 15 years using TMI rainfall data revealed that more than 60% of the area experienced extreme drought in 2011 over the MLRY region and in 1998 over the Huang–Huai–Hai region. The frequency of different intensity droughts presented significant spatial heterogeneity in each crop-growing region. Full article

The main objective of this research is to investigate the potential combination of Sentinel-2A and ALOS-2 PALSAR-2 (Advanced Land Observing Satellite -2 Phased Array type L-band Synthetic Aperture Radar-2) imagery for improving the accuracy of the Aboveground Biomass (AGB) measurement. According to the current literature, this kind of investigation has rarely been conducted. The Hyrcanian forest area (Iran) is selected as the case study. For this purpose, a total of 149 sample plots for the study area were documented through fieldwork. Using the imagery, three datasets were generated including the Sentinel-2A dataset, the ALOS-2 PALSAR-2 dataset, and the combination of the Sentinel-2A dataset and the ALOS-2 PALSAR-2 dataset (Sentinel-ALOS). Because the accuracy of the AGB estimation is dependent on the method used, in this research, four machine learning techniques were selected and compared, namely Random Forests (RF), Support Vector Regression (SVR), Multi-Layer Perceptron Neural Networks (MPL Neural Nets), and Gaussian Processes (GP). The performance of these AGB models was assessed using the coefficient of determination (R²), the root-mean-square error (RMSE), and the mean absolute error (MAE). The results showed that the AGB models derived from the combination of the Sentinel-2A and the ALOS-2 PALSAR-2 data had the highest accuracy, followed by models using the Sentinel-2A dataset and the ALOS-2 PALSAR-2 dataset. Among the four machine learning models, the SVR model (R² = 0.73, RMSE = 38.68, and MAE = 32.28) had the highest prediction accuracy, followed by the GP model (R² = 0.69, RMSE = 40.11, and MAE = 33.69), the RF model (R² = 0.62, RMSE = 43.13, and MAE = 35.83), and the MPL Neural Nets model (R² = 0.44, RMSE = 64.33, and MAE = 53.74). Overall, the Sentinel-2A imagery provides a reasonable result while the ALOS-2 PALSAR-2 imagery provides a poor result of the forest AGB estimation. The combination of the Sentinel-2A imagery and the ALOS-2 PALSAR-2 imagery improved the estimation accuracy of AGB compared to that of the Sentinel-2A imagery only. Full article

In the present study, we investigate the effects of signal to noise (SNR), slit function (FWHM), and aerosol optical depth (AOD) on the accuracy of formaldehyde (HCHO) vertical column density (HCHO_VCD) using the ground-based direct-sun synthetic radiance based on differential optical absorption spectroscopy (DOAS). We found that the effect of SNR on HCHO retrieval accuracy is larger than those of FWHM and AOD. When SNR = 650 (1300), FWHM = 0.6, and AOD = 0.2, the absolute percentage difference (APD) between the true HCHO_VCD values and those retrieved ranges from 54 (30%) to 5% (1%) for the HCHO_VCD of 5.0 × 10¹⁵ and 1.1 × 10¹⁷ molecules cm⁻², respectively. Interestingly, the maximum AOD effect on the HCHO accuracy was found for the HCHO_VCD of 3.0 × 10¹⁶ molecules cm⁻². In addition, we carried out the first ground-based direct-sun measurements in the ultraviolet (UV) wavelength range to retrieve the HCHO_VCD using Pandora in Seoul. The HCHO_VCD was low at 12:00 p.m. local time (LT) in all seasons, whereas it was high in the morning (10:00 a.m. LT) and late afternoon (4:00 p.m. LT), except in winter. The maximum HCHO_VCD values were 2.68 × 10¹⁶, 3.19 × 10¹⁶, 2.00 × 10¹⁶, and 1.63 × 10¹⁶ molecules cm⁻² at 10:00 a.m. LT in spring, 10:00 a.m. LT in summer, 1:00 p.m. LT in autumn, and 9:00 a.m. LT in winter, respectively. The minimum values of Pandora HCHO_VCD were 1.63 × 10¹⁶, 2.23 × 10¹⁶, 1.26 × 10¹⁶, and 0.82 × 10¹⁶ molecules cm⁻² at around 1:45 p.m. LT in spring, summer, autumn, and winter, respectively. This seasonal pattern of high values in summer and low values in winter implies that photo-oxidation plays an important role in HCHO production. The correlation coefficient (R) between the monthly HCHO_VCD values from Pandora and those from the Ozone Monitoring Instrument (OMI) is 0.61, and the slope is 1.25. Full article

Remote and extreme regions such as in the Arctic remain a challenging ground for geological mapping and mineral exploration. Coastal cliffs are often the only major well-exposed outcrops, but are mostly not observable by air/spaceborne nadir remote sensing sensors. Current outcrop mapping efforts rely on the interpretation of Terrestrial Laser Scanning and oblique photogrammetry, which have inadequate spectral resolution to allow for detection of subtle lithological differences. This study aims to integrate 3D-photogrammetry with vessel-based hyperspectral imaging to complement geological outcrop models with quantitative information regarding mineral variations and thus enables the differentiation of barren rocks from potential economic ore deposits. We propose an innovative workflow based on: (1) the correction of hyperspectral images by eliminating the distortion effects originating from the periodic movements of the vessel; (2) lithological mapping based on spectral information; and (3) accurate 3D integration of spectral products with photogrammetric terrain data. The method is tested using experimental data acquired from near-vertical cliff sections in two parts of Greenland, in Karrat (Central West) and Søndre Strømfjord (South West). Root-Mean-Square Error of (6.7, 8.4) pixels for Karrat and (3.9, 4.5) pixels for Søndre Strømfjord in X and Y directions demonstrate the geometric accuracy of final 3D products and allow a precise mapping of the targets identified using the hyperspectral data contents. This study highlights the potential of using other operational mobile platforms (e.g., unmanned systems) for regional mineral mapping based on horizontal viewing geometry and multi-source and multi-scale data fusion approaches. Full article

Recently, ground-based hyperspectral imaging has come to the fore, supporting the arduous task of mapping near-vertical, difficult-to-access geological outcrops. The application of outcrop sensing within a range of one to several hundred metres, including geometric corrections and integration with accurate terrestrial laser scanning models, is already developing rapidly. However, there are few studies dealing with ground-based imaging of distant targets (i.e., in the range of several kilometres) such as mountain ridges, cliffs, and pit walls. In particular, the extreme influence of atmospheric effects and topography-induced illumination differences have remained an unmet challenge on the spectral data. These effects cannot be corrected by means of common correction tools for nadir satellite or airborne data. Thus, this article presents an adapted workflow to overcome the challenges of long-range outcrop sensing, including straightforward atmospheric and topographic corrections. Using two datasets with different characteristics, we demonstrate the application of the workflow and highlight the importance of the presented corrections for a reliable geological interpretation. The achieved spectral mapping products are integrated with 3D photogrammetric data to create large-scale now-called “hyperclouds”, i.e., geometrically correct representations of the hyperspectral datacube. The presented workflow opens up a new range of application possibilities of hyperspectral imagery by significantly enlarging the scale of ground-based measurements. Full article

Circular scanning synthetic aperture radar (SAR) is a novel imaging mode wherein the radar antenna rotates from 0 degrees to 360 degrees along the platform flight direction, providing us with a potentially effective technique to estimate the sea surface current velocity. In this paper, we propose a novel method to estimate the sea surface current velocity utilizing the Doppler centroid shifts of different scan angles over 360 degrees after the airborne platform motion compensation. In this method, the Doppler centroid shifts of the sea clutter at different scan angles are first extracted, and the corresponding compensation errors caused by the azimuth pointing and the incidence angle of the radar beam are considered. Finally, the least squares (LS) technique is applied to estimate the along-track velocity component and the cross-track velocity component of the sea surface current. The effectiveness of the proposed method is verified by the real data recorded by an airborne circular scanning SAR system. Full article

An accurate estimation of the leaf area index (LAI) by satellite remote sensing is essential for studying the spatial variation of ecosystem structure. The goal of this study was to estimate the spatial variation of LAI over a forested catchment in a mountainous landscape (ca. 60 km²) in central Japan. We used a simple model to estimate LAI using spectral reflectance by adapting the Monsi-Saeki light attenuation theory for satellite remote sensing. First, we applied the model to Landsat Operational Land Imager (OLI) imagery to estimate the spatial variation of LAI in spring and summer. Second, we validated the model’s performance with in situ LAI estimates at four study plots that included deciduous broadleaf, deciduous coniferous, and evergreen coniferous forest types. Pre-processing of the Landsat OLI imagery, including atmospheric correction by elevation-dependent dark object subtraction and Minnaert topographic correction, together with application of the simple model, enabled a satisfactory 30-m spatial resolution estimation of forest LAI with a maximum of 5.5 ± 0.2 for deciduous broadleaf and 5.3 ± 0.2―for evergreen coniferous forest areas. The LAI variation in May (spring) suggested an altitudinal gradient in the degree of leaf expansion, whereas the LAI variation in August (mid-summer) suggested an altitudinal gradient of yearly maximum forest foliage density. This study demonstrated the importance of an accurate estimation of fine-resolution spatial LAI variations for ecological studies in mountainous landscapes, which are characterized by complex terrain and high vegetative heterogeneity. Full article

In the European Citclops project, with a prime aim of developing new tools to involve citizens in the water quality monitoring of natural waters, colour was identified as a simple property that can be measured via a smartphone app and by dedicated low-cost instruments. In a recent paper, we demonstrated that colour, as expressed mainly by the hue angle (α), can also be derived accurately and consistently from the ocean colour satellite instruments that have observed the Earth since 1997. These instruments provide superior temporal coverage of natural waters, albeit at a reduced spatial resolution of 300 m at best. In this paper, the list of algorithms is extended to the very first ocean colour instrument, and the Moderate Resolution Imaging Spectroradiometer (MODIS) 500-m resolution product. In addition, we explore the potential of the hue angle derivation from multispectral imaging instruments with a higher spatial resolution but reduced spectral resolution: the European Space Agency (ESA) multispectral imager (MSI) on Sentinel-2 A,B, the Operational Land Imager (OLI) on the National Aeronautics and Space Administration (NASA) Landsat-8, and its precursor, the Enhanced Thematic Mapper Plus (ETM+) on Landsat-7. These medium-resolution imagers might play a role in an upscaling from point measurements to the typical 1-km pixel size from ocean colour instruments. As the parameter α (the colour hue angle) is fairly new to the community of water remote sensing scientists, we present examples of how colour can help in the image analysis in terms of water-quality products. Full article

Remote sensing and field spectroscopy of natural waters is typically performed under clear skies, low wind speeds and low solar zenith angles. Such measurements can also be made, in principle, under clouds and mixed skies using airborne or in-situ measurements; however, variable illumination conditions pose a challenge to data analysis. In the present case study, we evaluated the inversion of hyperspectral in-situ measurements for water constituent retrieval acquired under variable cloud cover. First, we studied the retrieval of Chlorophyll-a (Chl-a) concentration and colored dissolved organic matter (CDOM) absorption from in-water irradiance measurements. Then, we evaluated the errors in the retrievals of the concentration of total suspended matter (TSM), Chl-a and the absorption coefficient of CDOM from above-water reflectance measurements due to highly variable reflections at the water surface. In order to approximate cloud reflections, we extended a recent three-component surface reflectance model for cloudless atmospheres by a constant offset and compared different surface reflectance correction procedures. Our findings suggest that in-water irradiance measurements may be used for the analysis of absorbing compounds even under highly variable weather conditions. The extended surface reflectance model proved to contribute to the analysis of above-water reflectance measurements with respect to Chl-a and TSM. Results indicate the potential of this approach for all-weather monitoring. Full article

Earth observation satellites with various spatial, spectral and temporal resolutions provide an invaluable means for mapping and monitoring the Earth’s environments. With the increasing demand of satellite images for remote and harsh environments and nature disaster areas such as earthquake, flooding, bushfires and other emergent events, quickly geo-positioning those images without using ground control points (GCPs) is much preferable and desirable. Built on the previously developed Spatial Triangulated Network (STN) concept by the first author, this paper presents a Rational Function Model (RFM) based geo-positioning method utilizing some pre-orientated image(s) as reference, instead of using GCPs. The experimental results indicate that the RFM method is more sensitive to the base-height ratio in the vertical accuracy than the physical model based geo-positioning method which was also developed by the first author. Compared to the traditional RFM based block adjustment using GCPs, the proposed RFM based method without GCP (using orientated images instead) can achieve similar accuracies when more than one orientated image, which have reasonable strong geometric relationships with the new images, are introduced into the proposed RFM based method. The proposed method is applicable to the scenarios in which geo-positioning is required for those new satellite images that only have RFM and no GCPs available, but where there exists some orientated images covering the same region. Full article

Multispectral Airborne Laser Scanning (ALS) is a new technology and its output data have not been fully explored for tree species classification purposes. The objective of this study was to investigate what type of features from multispectral ALS data (wavelengths of 1550 nm, 1064 nm and 532 nm) are best suited for tree species classification. Remote sensing data were gathered over hemi-boreal forest in southern Sweden (58°27′18.35″N, 13°39′8.03″E) on 21 July 2016. The field data consisted of 179 solitary trees from nine genera and ten species. Two new methods for feature extraction were tested and compared to features of height and intensity distributions. The features that were most important for tree species classification were intensity distribution features. Features from the upper part of the upper and outer parts of the crown were better for classification purposes than others. The best classification model was created using distribution features of both intensity and height in multispectral data, with a leave-one-out cross-validated accuracy of 76.5%. As a comparison, only structural features resulted in an highest accuracy of 43.0%. Picea abies and Pinus sylvestris had high user’s and producer’s accuracies and were not confused with any deciduous species. Tilia cordata was the deciduous species with a large sample that was most frequently confused with many other deciduous species. The results, although based on a small and special data set, suggest that multispectral ALS is a technology with great potential for tree species classification. Full article

To reduce the floor space of receiving antenna arrays, the Radio Ocean Remote SEnsing (RORSE) laboratory of Wuhan University developed a circular receiving array for a multi-frequency high frequency (MHF) radar system in 2014, consisting of seven uniformly spaced antenna elements positioned on a circle with a diameter of 5 m. The new system, which is abbreviated MHF-C radar, adopts frequency modulated interrupted continuous wave (FMICW) chirps and is capable of simultaneously operating at a maximum of four frequencies in the band of 7.5–25 MHz, and providing current, wave and wind maps every ten minutes. The phase direction-finding method is utilized to estimate the directions of the current signals, and array phase uncertainties are also taken into consideration in the signal model. This paper introduces the system in detail and investigates the performance of current measurements using MHF-C radars installed at Shengshan and Zhujiajian along the coast of the East China Sea. Radial current measurements derived from 8.27 MHz and 19.20 MHz at the same range are compared. Observations and comparisons between MHF-C radars and acoustic Doppler current profilers (ADCPs) are also presented in this paper. The results preliminarily demonstrate that the MHF-C radar system is capable of maintaining the same performance for current measurements whenever it steers to any other azimuth in the coverage and has a good ability to measure currents. Full article

Downward shortwave radiation (DSR) is an essential parameter in the terrestrial radiation budget and a necessary input for models of land-surface processes. Although several radiation products using satellite observations have been released, coarse spatial resolution and low accuracy limited their application. It is important to develop robust and accurate retrieval methods with higher spatial resolution. Machine learning methods may be powerful candidates for estimating the DSR from remotely sensed data because of their ability to perform adaptive, nonlinear data fitting. In this study, the gradient boosting regression tree (GBRT) was employed to retrieve DSR measurements with the ground observation data in China collected from the China Meteorological Administration (CMA) Meteorological Information Center and the satellite observations from the Advanced Very High Resolution Radiometer (AVHRR) at a spatial resolution of 5 km. The validation results of the DSR estimates based on the GBRT method in China at a daily time scale for clear sky conditions show an R² value of 0.82 and a root mean square error (RMSE) value of 27.71 W·m⁻² (38.38%). These values are 0.64 and 42.97 W·m⁻² (34.57%), respectively, for cloudy sky conditions. The monthly DSR estimates were also evaluated using ground measurements. The monthly DSR estimates have an overall R² value of 0.92 and an RMSE of 15.40 W·m⁻² (12.93%). Comparison of the DSR estimates with the reanalyzed and retrieved DSR measurements from satellite observations showed that the estimated DSR is reasonably accurate but has a higher spatial resolution. Moreover, the proposed GBRT method has good scalability and is easy to apply to other parameter inversion problems by changing the parameters and training data. Full article

Unmanned aerial vehicles (UAVs) and digital photogrammetric techniques are two recent advances in remote sensing (RS) technology that are emerging as alternatives to high-cost airborne laser scanning (ALS) data sources. Despite the potential of UAVs in forestry applications, very few studies have included detailed analyses of UAV photogrammetric products at larger scales or over a range of forest types, including mixed conifer–broadleaf forests. In this study, we assessed the performance of fixed-wing UAV photogrammetric products of a mixed conifer–broadleaf forest with varying levels of canopy structural complexity. We demonstrate that fixed-wing UAVs are capable of efficiently collecting image data at local scales and that UAV imagery can be effectively utilized with digital photogrammetric techniques to provide detailed automated reconstruction of the three-dimensional (3D) canopy surface of mixed conifer–broadleaf forests. When combined with an accurate digital terrain model (DTM), UAV photogrammetric products are promising for producing reliable structural measurements of the forest canopy. However, the performance of UAV photogrammetric products is likely to be influenced by the structural complexity of the forest canopy. Furthermore, we highlight the potential of fixed-wing UAVs in operational forest management at the forest management compartment level, for acquiring high-resolution imagery at low cost. A future direction of this research would be to address the issue of how well the photogrammetric products can predict the actual structure of mixed conifer–broadleaf forests. Full article

The contribution to sea level rise from Patagonian icefields is one of the largest mass losses outside the large ice sheets of Antarctica and Greenland. However, only a few studies have provided large-scale assessments in a spatially detailed way to address the reaction of individual glaciers in Patagonia and hence to better understand and explain the underlying processes. In this work, we use repeat radar interferometric measurements of the German TerraSAR-X-Add-on for Digital Elevation Measurements (TanDEM-X) satellite constellation between 2011/12 and 2016 together with the digital elevation model from the Shuttle Radar Topography Mission (SRTM) in 2000 in order to derive surface elevation and mass changes of the Southern Patagonia Icefield (SPI). Our results reveal a mass loss rate of −11.84 ± 3.3 Gt·a⁻¹ (corresponding to 0.033 ± 0.009 mm·a⁻¹ sea level rise) for an area of 12573 km² in the period 2000–2015/16. This equals a specific glacier mass balance of −0.941 ± 0.19 m w.e.·a⁻¹ for the whole SPI. These values are comparable with previous estimates since the 1970s, but a magnitude larger than mass change rates reported since the Little Ice Age. The spatial pattern reveals that not all glaciers respond similarly to changes and that various factors need to be considered in order to explain the observed changes. Our multi-temporal coverage of the southern part of the SPI (south of 50.3° S) shows that the mean elevation change rates do not vary significantly over time below the equilibrium line. However, we see indications for more positive mass balances due to possible precipitation increase in 2014 and 2015. We conclude that bi-static radar interferometry is a suitable tool to accurately measure glacier volume and mass changes in frequently cloudy regions. We recommend regular repeat TanDEM-X acquisitions to be scheduled for the maximum summer melt extent in order to minimize the effects of radar signal penetration and to increase product quality. Full article

The unique feature of the Visible Infrared Imager Radiometer Suite (VIIRS) day/night band (DNB) is its ability to take quantitative measurements of low-light scenes at night. In order to monitor the stability of the high gain stage (HGS) of the DNB, nighttime observations over the Dome C site under moonlight are analyzed in this study. The Miller and Turner 2009 (MT2009) lunar irradiance model has been used to simulate lunar illumination over Dome C. However, the MT2009 model does not differentiate the waxing and waning lunar phases. In this paper, the MT-SWC (SeaWiFS Corrected) lunar irradiance model differentiating the waxing and waning lunar phases is derived by correcting the MT2009 model using lunar observations made by the Sea-Viewing Wide Field-of-View Sensor (SeaWiFS). In addition, a top of atmosphere (TOA) bi-directional reflectance distribution function (BRDF) model during nighttime over Dome C is developed to remove the angular dependence from the nighttime TOA reflectance. The long-term stability monitoring of the DNB high-gain stage (HGS) reveals a lower reflectance factor in 2012 in comparison to the following years, which can be traced back to the change in relative spectral response (RSR) of National Oceanic & Atmospheric Administration’s (NOAA’s) Interface Data Processing Segment (IDPS) VIIRS DNB in April 2013. It also shows the radiometric stability of DNB data, with long-term stability of less than 1.58% over the periods from 2013 to 2016. This method can be used to monitor the radiometric stability of other low-light observing sensors using vicarious calibration sites under moonlight illumination. Full article

In this paper, we present a novel method for ship classification in synthetic aperture radar (SAR) images. The proposed method consists of feature extraction and classifier training. Inspired by SAR-HOG feature in automatic target recognition, we first design a novel feature named MSHOG by improving SAR-HOG, adapting it to ship classification, and employing manifold learning to achieve dimensionality reduction. Then, we train the classifier and dictionary jointly in task-driven dictionary learning (TDDL) framework. To further improve the performance of TDDL, we enforce structured incoherent constraints on it and develop an efficient algorithm for solving corresponding optimization problem. Extensive experiments performed on two datasets with TerraSAR-X images demonstrate that the proposed method, MSHOG feature and TDDL with structured incoherent constraints, outperforms other existing methods and achieves state-of-art performance. Full article

Accurate data about the spatial distribution and planting area of maize is important for policy making, economic development, environmental protection and food security under climate change. This paper proposes a new identification method for spring maize based on spectral and phenological features derived from the moderate resolution imaging spectroradiometer (MODIS) land surface reflectance time-series data. The method focused on the spectral differences of different land cover types in the specific phenological phases of spring maize by testing the selections and combinations of classification metrics, feature extraction methods and classifiers. Taking Liaoning province, a representative planting region of spring maize in Northeast China, as the study area, the results indicated that the combined multiple metrics, including the red reflectance, near-infrared reflectance and normalized difference vegetation index (NDVI), were conducive to the maize identification and were better than any single metric. With regard to the feature extraction and selection, maize identification based on different phenological features selected with prior knowledge was more efficient than that based on statistical features derived from the principal component analysis. Compared with the maximum likelihood classification method, the decision tree classification based on expert knowledge was more suitable for phenological features selected from some prior knowledge. In summary, discriminant rules were defined with those phenological features from multiple metrics, and the decision tree classification was used to identify maize in the study area. The producer’s accuracy of maize identification was 98.57%, and the user’s accuracy was 81.18%. This method can be potentially applied to an operational identification of maize at large scales based on remote sensing time-series data. Full article

Accurate and detailed monitoring of population distribution and evolution is of great significance in formulating a population planning strategy in China. The Defense Meteorological Satellite Program’s Operational Linescan System (DMSP/OLS) nighttime lights time-series (NLT) image products offer a good opportunity for detecting the population distribution owing to its high correlation to human activities. However, their detection capability is greatly limited owing to a lack of in-flight calibration. At present, the synergistic use of systematically-corrected NLT products and population spatialization is rarely applied. This work proposed a methodology to improve the application precision and versatility of NLT products, explored a feasible approach to quantitatively spatialize the population to grid units of $1\mathrm{km}\times 1\mathrm{km}$ , and revealed the spatio-temporal characteristics of population distribution from 2000 to 2010. Results indicated that, (1) after inter-calibration, geometric, incompatibility and discontinuity corrections, and adjustment based on vegetation information, the incompatibility and discontinuity of NTL products were successfully solved. Accordingly, detailed actual residential areas and luminance differences between the urban core and the peripheral regions could be obtained. (2) The population spatialization method could effectively acquire population information at per ${\mathrm{km}}^2$ with high accuracy and exhibit more details in the evolution of population distribution. (3) Obvious differences in spatio-temporal characteristics existed in four economic regions, from the aspects of population distribution and dynamics, as well as population-weighted centroids. The eastern region was the most populous with the largest increased magnitude, followed by the central, northeastern, and western regions. The population-weighted centroids of the eastern, western, and northeastern regions moved along the southwest direction, while the population-weighted centroid of the central region moved along the southeast direction. (4) The population distribution and dynamics in nine-level population density types were significantly different. In the period of 2000–2010, the population in the basic no-man and high concentration types presented a net decrease. The population in seven other regions all increased with a net increase ranging from 25 ${\mathrm{km}}^2$ (the moderate concentration type) to 245,668 ${\mathrm{km}}^2$ (the general transition type). Except those in the core concentration and extremely sparse types, the population-weighted centroids in all other population density types moved along the southwest direction. Full article

Actual evapotranspiration is assessed via surface energy balance at an hourly rate. However, a robust estimation of daily evapotranspiration from hourly values is required. Outcomes of surface energy balance are frequently determined via measures of eddy covariance latent heat flux. Surface energy balance can be applied on images acquired at different times and spatial resolutions. In addition, hourly actual evapotranspiration needs to be integrated at a daily rate for operational uses. Questions arise whether the validation of surface energy balance models can benefit from complementary in situ measures of latent heat flux and sap flow transpiration. Here, validation was driven by image acquisition time, spatial resolution, and temporal integration. Thermal and optical images were collected with a proximity-sensing platform on an olive orchard at different acquisition times. Actual latent heat fluxes from canopy and sap flux at tree trunks were measured with a flux tower and heat dissipation probes. The latent heat fluxes were then further analyzed. A surface energy balance was applied over proximity sensing images re-sampled at different spatial resolutions with resulting latent heat fluxes compared to in situ ones. A time lag was observed and quantified between actual latent heat fluxes from canopy and sap flux at the tree trunk. Results also indicate that a pixel resolution comparable to the average canopy size was suitable for estimating the actual evapotranspiration via a single source surface energy balance model. Images should not be acquired at the beginning or the end of the diurnal period. Findings imply that sap flow transpiration can be used to measure surface energy balance at a daily rate or when images are found at an hourly rate near noon, and the existing time lag between the latent heat flux at the canopy and the sap flow at the trunk does not need to be taken into account. Full article

In this paper, to break the limit of the traditional linear models for synthetic aperture radar (SAR) image despeckling, we propose a novel deep learning approach by learning a non-linear end-to-end mapping between the noisy and clean SAR images with a dilated residual network (SAR-DRN). SAR-DRN is based on dilated convolutions, which can both enlarge the receptive field and maintain the filter size and layer depth with a lightweight structure. In addition, skip connections and a residual learning strategy are added to the despeckling model to maintain the image details and reduce the vanishing gradient problem. Compared with the traditional despeckling methods, the proposed method shows a superior performance over the state-of-the-art methods in both quantitative and visual assessments, especially for strong speckle noise. Full article

Satellite remote sensing has been widely used to retrieve aerosol optical depth (AOD), which is an indicator of air quality as well as radiative forcing. The dark target (DT) algorithm is applied to low reflectance areas, such as dense vegetation, and the deep blue (DB) algorithm is adopted for bright-reflecting regions. However, both DT and DB algorithms ignore the effect of surface bidirectional reflectance. This paper provides a method for AOD retrieval in arid or semiarid areas, in which the key points are the accurate estimation of surface reflectance and reasonable assumptions of the aerosol model. To reduce the uncertainty in surface reflectance, a minimum land surface reflectance database at the spatial resolution of 500 m for each month was constructed based on the moderate-resolution imaging spectroradiometer (MODIS) surface reflectance product. Furthermore, a bidirectional reflectance distribution function (BRDF) correction model was adopted to compensate for the effect of surface reflectance anisotropy. The aerosol parameters, including AOD, single scattering albedo, asymmetric factor, Ångström exponent and complex refractive index, are determined based on the observation of two sunphotometers installed in northern Xinjiang from July to August 2014. The AOD retrieved from the MODIS images was validated with ground-based measurements and the Terra-MODIS aerosol product (MOD04). The 500 m AOD retrieved from the MODIS showed high consistency with ground-based AOD measurements, with an average correlation coefficient of ~0.928, root mean square error (RMSE) of ~0.042, mean absolute error (MAE) of ~0.032, and the percentage falling within the expected error (EE) of the collocations is higher than that for the MOD04 DB product. The results demonstrate that the new AOD algorithm is more suitable to represent aerosol conditions over Xinjiang than the DB standard product. Full article

Wind speed and direction are important essential climate variables (ECVs). GNSS-R is an emerging remote sensing technique that can be potentially used to retrieve wind speed from space. However, few studies have addressed the wind direction retrieval from spaceborne GNSS-R observables, namely the Delay Doppler map (DDM). In this study, the feasibility of retrieving wind direction from the synthetic DDMs is analyzed. First, the simulation tool P²EPS is used to generate the DDMs under different geometry configurations, wind speed, and wind direction. Then, DDM changes caused by the wind direction are investigated, and two metrics are proposed to retrieve the wind direction from the DDM shape changes. The influence on wind direction retrieval of the wind speed, receiver’s elevation, and azimuth is further discussed. Finally, the sensitivity of DDM changes to noise is investigated, as well as the impact of noise on these two metrics. Full article

Spatial information on forest functional composition is needed to inform management and conservation efforts, yet this information is lacking, particularly in tropical regions. Canopy foliar traits underpin the functional biodiversity of forests, and have been shown to be remotely measurable using airborne 350–2510 nm imaging spectrometers. We used newly acquired imaging spectroscopy data constrained with concurrent light detection and ranging (LiDAR) measurements from the Carnegie Airborne Observatory (CAO), and field measurements, to test the performance of the Spectranomics approach for foliar trait retrieval. The method was previously developed in Neotropical forests, and was tested here in the humid tropical forests of Malaysian Borneo. Multiple foliar chemical traits, as well as leaf mass per area (LMA), were estimated with demonstrable precision and accuracy. The results were similar to those observed for Neotropical forests, suggesting a more general use of the Spectranomics approach for mapping canopy traits in tropical forests. Future mapping studies using this approach can advance scientific investigations and applications based on imaging spectroscopy. Full article

The atmospheric correction of optical remote sensing data requires the determination of aerosol and gas optical properties. A method is presented which allows the detection of the aerosol scattering effects from optical remote sensing data at spatial sampling intervals below 5 m in cloud-free situations from cast shadow pixels. The derived aerosol optical thickness distribution is used for improved atmospheric compensation. In a first step, a novel spectral cast shadow detection algorithm determines the shadow areas using spectral indices. Evaluation of the cast shadow masks shows an overall classification accuracy on an 80% level. Using the such derived shadow map, the ATCOR atmospheric compensation method is iteratively applied on the shadow areas in order to find the optimum aerosol amount. The aerosol optical thickness is found by analyzing the physical atmospheric correction of fully shaded pixels in comparison to directly illuminated areas. The shadow based aerosol optical thickness estimation method (SHAOT) is tested on airborne imaging spectroscopy data as well as on photogrammetric data. The accuracy of the reflectance values from atmospheric correction using the such derived aerosol optical thickness could be improved from 3–4% to a level of better than 2% in reflectance for the investigated test cases. Full article

To obtain the high-resolution seabed topography and overcome the limitations of existing topography reconstruction methods in requiring external bathymetric data and ignoring the effects of sediment variations and Side-Scan Sonar (SSS) image quality, this study proposes a method of reconstructing seabed topography from SSS images with a self-constraint condition. A reconstruction model is deduced by Lambert’s law and the seabed scattering model. A bottom tracking method is put forward to get the along-track SSS towfish heights and the initial seabed topography in the SSS measuring area is established by combining the along-track towfish heights, towfish depths and tidal levels obtained from Global Navigation Satellite System (GNSS). The complete process of reconstructing seabed topography is given by taking the initial topography as self-constraint and the high-resolution seabed topography is finally obtained. Experiments verified the proposed method by the data measured in Zhujiang River, China. The standard deviation of less than 15 cm is achieved and the resolution of the reconstructed topography is about 60 times higher than that of the Digital Elevation Model (DEM) established by bathymetric data. The effects of noise, suspended bodies, refraction of wave in water column, sediment variation, the determination of iteration termination condition as well as the performance of the proposed method under these effects are discussed. Finally, the conclusions are drawn out according to the experiments and discussions. The proposed method provides a simple and efficient way to obtain high-resolution seabed topography from SSS images and is a supplement but not substitution for the existing bathymetric methods. Full article

The detection of water stress in vineyards plays an integral role in the sustainability of high-quality grapes and prevention of devastating crop loses. Hyperspectral remote sensing technologies combined with machine learning provides a practical means for modelling vineyard water stress. In this study, we applied two ensemble learners, i.e., random forest (RF) and extreme gradient boosting (XGBoost), for discriminating stressed and non-stressed Shiraz vines using terrestrial hyperspectral imaging. Additionally, we evaluated the utility of a spectral subset of wavebands, derived using RF mean decrease accuracy (MDA) and XGBoost gain. Our results show that both ensemble learners can effectively analyse the hyperspectral data. When using all wavebands (p = 176), RF produced a test accuracy of 83.3% (KHAT (kappa analysis) = 0.67), and XGBoost a test accuracy of 80.0% (KHAT = 0.6). Using the subset of wavebands (p = 18) produced slight increases in accuracy ranging from 1.7% to 5.5% for both RF and XGBoost. We further investigated the effect of smoothing the spectral data using the Savitzky-Golay filter. The results indicated that the Savitzky-Golay filter reduced model accuracies (ranging from 0.7% to 3.3%). The results demonstrate the feasibility of terrestrial hyperspectral imagery and machine learning to create a semi-automated framework for vineyard water stress modelling. Full article

The area of interest is located on the eastern flank of the Andean Cordillera, San Juan province, Argentina. The 3600 km² area is characterized by Siluro-Devonian to Neogene sedimentary and igneous rocks and unconsolidated Quaternary sediments. Epithermal, porphyry-related, and magmatic-hydrothermal breccia-hosted ore deposits, common in this part of the Frontal Cordillera, are associated with various types of hydrothermal alteration assemblages. Kaolinite – alunite-rich argillic, quartz – illite-rich phyllic, epidote – chlorite – calcite-rich propylitic and silicic are the most common hydrothermal alteration assemblages in the study area. VNIR, SWIR and TIR ASTER data were used to characterize geological features on a portion of the Frontal Cordillera. Red-green-blue band combinations, band ratios, logical operations, mineral indices and principal component analysis were applied to successfully identify rock types and hydrothermal alteration zones in the study area. These techniques were used to enhance geological features to contrast different lithologies and zones with high concentrations of argillic, phyllic, propylitic alteration mineral assemblages and silicic altered rocks. Alteration minerals detected with portable short-wave infrared spectrometry in hand specimens confirmed the capability of ASTER to identify hydrothermal alteration assemblages. The results from field control areas confirmed the presence of those minerals in the areas classified by ASTER processing techniques and allowed mapping the same mineralogy where pixels had similar information. The current study proved ASTER processing techniques to be valuable mapping tools for geological reconnaissance of a large area of the Argentinean Frontal Cordillera, providing preliminary lithologic and hydrothermal alteration maps that are accurate as well as cost and time effective. Full article

The integration of Global Positioning System (GPS) real-time kinematics (RTK) and an inertial navigation system (INS) has been widely used in many applications, such as mobile mapping and autonomous vehicle control. Such applications require high-accuracy position information. However, continuous and reliable high-accuracy positioning is still challenging for GPS/INS integration in urban environments because of the limited satellite visibility, increasing multipath, and frequent signal blockages. Recently, with the rapid deployment of multi-constellation Global Navigation Satellite System (multi-GNSS) and the great advances in low-cost micro-electro-mechanical-system (MEMS) inertial measurement units (IMUs), it is expected that the positioning performance could be improved significantly. In this contribution, the tightly-coupled single-frequency multi-GNSS RTK/MEMS-IMU integration is developed to provide precise and continuous positioning solutions in urban environments. The innovation-based outlier-resistant ambiguity resolution (AR) and Kalman filtering strategy are proposed specifically for the integrated system to resist the measurement outliers or poor-quality observations. A field vehicular experiment was conducted in Wuhan City to evaluate the performance of the proposed algorithm. Results indicate that it is feasible for the proposed algorithm to obtain high-accuracy positioning solutions in the presence of measurement outliers. Moreover, the tightly-coupled single-frequency multi-GNSS RTK/MEMS-IMU integration even outperforms the dual-frequency multi-GNSS RTK in terms of AR and positioning performance for short baselines in urban environments. Full article

In this study, 27 polarimetric parameters were extracted from Radarsat-2 polarimetric synthetic aperture radar (SAR) at each growth stage of the rape crop. The sensitivity to growth parameters such as stem height, leaf area index (LAI), and biomass were investigated as a function of days after sowing. Based on the sensitivity analysis, five empirical regression models were compared to determine the best model for stem height, LAI, and biomass inversion. Of these five models, quadratic models had higher R² values than other models in most cases of growth parameter inversions, but when these results were related to physical scattering mechanisms, the inversion results produced overestimation in the performance of some parameters. By contrast, linear and logarithmic models, which had lower R² values than the quadratic models, had stable performance for growth parameter inversions, particularly in terms of their performance at each growth stage. The best biomass inversion performance was acquired by the volume component of a quadratic model, with an R² value of 0.854 and root mean square error (RMSE) of 109.93 g m⁻². The best LAI inversion was also acquired by a quadratic model, but used the radar vegetation index (Cloude), with an R² value of 0.8706 and RMSE of 0.56 m² m⁻². Stem height was acquired by scattering angle alpha ( $\mathsf{\alpha }$ ) using a logarithmic model, with an R² of 0.926 value and RMSE of 11.09 cm. The performances of these models were also analysed for biomass estimation at the second growth stage (P2), third growth stage (P3), and fourth growth stage (P4). The results showed that the models built at the P3 stage had better substitutability with the models built during all of the growth stages. From the mapping results, we conclude that a model built at the P3 stage can be used for rape biomass inversion, with 90% of estimation errors being less than 100 g m⁻². Full article

This study addresses the potential contribution of the first pair of Galileo FOC satellites sent into incorrect highly eccentric orbits for geodetic and surveying applications. We began with an analysis of the carrier to noise density ratio and the stochastic properties of GNSS measurements. The investigations revealed that the signal power of E14 & E18 satellites is higher than for regular Galileo satellites, what is related to their lower altitude over the experiment area. With regard to the noise of the observables, there are no significant differences between all Galileo satellites. Furthermore, the study confirmed that the precision of Galileo data is higher than that of GPS, especially in the case of code measurements. Next analysis considered selected domains of precise instantaneous medium-range positioning: ambiguity resolution and coordinate accuracy as well as observable residuals. On the basis of test solutions, with and without E14 & E18 data, we found that these satellites did not noticeably influence the ambiguity resolution process. The discrepancy in ambiguity success rate between test solutions did not exceed 2%. The differences between standard deviations of the fixed coordinates did not exceed 1 mm for horizontal components. The standard deviation of the L1/E1 phase residuals, corresponding to regular GPS and Galileo, and E14 & E18 satellite signals, was at a comparable level, in the range of 6.5–8.7 mm. The study revealed that the Galileo satellites with incorrect orbits were fully usable in most geodetic, surveying and many other post-processed applications and may be beneficial especially for positioning during obstructed visibility of satellites. This claim holds true when providing precise ephemeris of satellites. Full article

The recently available Landsat Analysis Ready Data (ARD) are provided as top of atmosphere (TOA) and atmospherically corrected (surface) reflectance tiled products and are designed to make the U.S. Landsat archive for the United States straightforward to use. In this study, the utility of ARD for 30 m percent tree cover mapping is demonstrated and the impact of different ARD processing levels on mapping accuracy examined. Five years of Landsat 5 and 7 ARD over 12 tiles encompassing Washington State are considered using an established bagged regression tree methodology and training data derived from Goddard LiDAR Hyperspectral & Thermal Imager (G-LiHT) data. Sensitivity to the amount of training data is examined with increasing mapping accuracy observed as more training data are used. Four processing levels of ARD are considered independently and the mapped results are compared: (i) TOA ARD; (ii) surface ARD; (iii) bidirectional reflectance distribution function (BRDF) adjusted atmospherically corrected ARD; and (iv) weekly composited BRDF adjusted atmospherically corrected ARD. The atmospherically corrected ARD provide marginally the highest mapping accuracies, although accuracy differences are negligible among the four (≤0.07% RMSE) when modest amounts of training data are used. The TOA ARD provide the most accurate maps compared to the other input data when only small amounts of training data are used, and the least accurate maps otherwise. The results are illustrated and the implications discussed. Full article

Soil moisture, especially surface soil moisture (SSM), plays an important role in the development of various natural hazards that result from extreme weather events such as drought, flooding, and landslides. There have been many remote sensing methods for soil moisture retrieval based on microwave or optical thermal infrared (TIR) measurements. TIR remote sensing has been popular for SSM retrieval due to its fine spatial and temporal resolutions. However, because of limitations in the penetration of optical TIR radiation and cloud cover, TIR methods can only be used under clear sky conditions. Microwave SSM retrieval is based on solid physical principles, and has advantages in cases of cloud cover, but it has low spatial resolution. For applications at the local scale, SSM data at high spatial and temporal resolutions are important, especially for agricultural management and decision support systems. Current remote sensing measurements usually have either a high spatial resolution or a high temporal resolution, but not both. This study aims to retrieve SSM at both high spatial and temporal resolutions through the fusion of Moderate Resolution Imaging Spectroradiometer (MODIS) and Land Remote Sensing Satellite (Landsat) data. Based on the universal triangle trapezoid, this study investigated the relationship between land surface temperature (LST) and the normalized difference vegetation index (NDVI) under different soil moisture conditions to construct an improved nonlinear model for SSM retrieval with LST and NDVI. A case study was conducted in Iowa, in the United States (USA) (Lat: 42.2°~42.7°, Lon: −93.6°~−93.2°), from 1 May 2016 to 31 August 2016. Daily SSM in an agricultural area during the crop-growing season was downscaled to 120-m spatial resolution by fusing Landsat 8 with MODIS, with an R² of 0.5766, and RMSE from 0.0302 to 0.1124 m³/m³. Full article

The extraction of a valuable set of features and the design of a discriminative classifier are crucial for target recognition in SAR image. Although various features and classifiers have been proposed over the years, target recognition under extended operating conditions (EOCs) is still a challenging problem, e.g., target with configuration variation, different capture orientations, and articulation. To address these problems, this paper presents a new strategy for target recognition. We first propose a low-dimensional representation model via incorporating multi-manifold regularization term into the low-rank matrix factorization framework. Two rules, pairwise similarity and local linearity, are employed for constructing multiple manifold regularization. By alternately optimizing the matrix factorization and manifold selection, the feature representation model can not only acquire the optimal low-rank approximation of original samples, but also capture the intrinsic manifold structure information. Then, to take full advantage of the local structure property of features and further improve the discriminative ability, local sparse representation is proposed for classification. Finally, extensive experiments on moving and stationary target acquisition and recognition (MSTAR) database demonstrate the effectiveness of the proposed strategy, including target recognition under EOCs, as well as the capability of small training size. Full article

Over several decades, improving the accuracy of Sea-Surface Temperatures (SSTs) derived from satellites has been a subject of intense research, and continues to be so. Knowledge of the accuracy of the SSTs is critical for weather and climate predictions, and many research and operational applications. In 2015, the operational Japanese MTSAT-2 geostationary satellite was replaced by the Himawari-8, which has a visible and infrared imager with higher spatial and temporal resolutions than its predecessor. In this study, data from both satellites during a three-month overlap period were compared with subsurface in situ temperature measurements from the Tropical Atmosphere Ocean (TAO) array and self-recording thermometers at the depths of corals of the Great Barrier Reef. Results show that in general the Himawari-8 provides more accurate SST measurements compared to those from MTSAT-2. At various locations, where in situ measurements were taken, the mean Himawari-8 SST error shows an improvement of ~0.15 K. Sources of the differences between the satellite-derived SST and the in situ temperatures were related to wind speed and diurnal heating. Full article

This study developed a retrieval algorithm for reflected shortwave radiation at the top of the atmosphere (RSR). This algorithm is based on Himawari-8/AHI (Advanced Himawari Imager) whose sensor characteristics and observation area are similar to the next-generation Geostationary Korea Multi-Purpose Satellite/Advanced Meteorological Imager (GK-2A/AMI). This algorithm converts the radiance into reflectance for six shortwave channels and retrieves the RSR with a regression coefficient look-up-table according to geometry of the solar-viewing (solar zenith angle, viewing zenith angle, and relative azimuth angle) and atmospheric conditions (surface type and absence/presence of clouds), and removed sun glint with high uncertainty. The regression coefficients were calculated using numerical experiments from the radiative transfer model (SBDART), and ridge regression for broadband albedo at the top of the atmosphere (TOA albedo) and narrowband reflectance considering anisotropy. The retrieved RSR were validated using Terra, Aqua, and S-NPP/CERES data on the 15th day of every month from July 2015 to February 2017. The coefficient of determination (R²) between AHI and CERES for scene analysis was higher than 0.867 and the Bias and root mean square error (RMSE) were −21.34–5.52 and 51.74–59.28 Wm⁻². The R², Bias, and RMSE for the all cases were 0.903, −2.34, and 52.12 Wm⁻², respectively. Full article

This paper provides a comparative study between microwave tomography and synthetic time-reversal imaging techniques as applied to ground penetrating radar (GPR) surveys. The comparison is carried out by processing experimental data collected at a controlled test site, with different types of buried targets at given subsurface depths and representative soil conditions. It is shown that the two techniques allow us to obtain complementary information about position, depth and size of the targets from a single GPR survey. Full article

Permafrost distribution in the Qinghai-Tibet Engineering Corridor (QTEC) is of growing interest due to the increase in infrastructure development in this remote area. Empirical models of mountain permafrost distribution have been established based on field sampled data, as a tool for regional-scale assessments of its distribution. This kind of model approach has never been applied for a large portion of this engineering corridor. In the present study, this methodology is applied to map permafrost distribution throughout the QTEC. After spatial modelling of the mean annual air temperature distribution from MODIS-LST and DEM, using high-resolution satellite image to interpret land surface type, a permafrost probability index was obtained. The evaluation results indicate that the model has an acceptable performance. Conditions highly favorable to permafrost presence (≥70%) are predicted for 60.3% of the study area, declaring a discontinuous permafrost distribution in the QTEC. This map is useful for the infrastructure development along the QTEC. In the future, local ground-truth observations will be required to confirm permafrost presence in favorable areas and to monitor permafrost evolution under the influence of climate change. Full article

Ensuring access to water and sanitation for all is Goal No. 6 of the 17 UN Sustainability Development Goals to transform our world. As one step towards this goal, we present an approach that leverages remote sensing data to plan optimal water supply networks for informal urban settlements. The concept focuses on slums within large urban areas, which are often characterized by a lack of an appropriate water supply. We apply methods of mathematical optimization aiming to find a network describing the optimal supply infrastructure. Hereby, we choose between different decentral and central approaches combining supply by motorized vehicles with supply by pipe systems. For the purposes of illustration, we apply the approach to two small slum clusters in Dhaka and Dar es Salaam. We show our optimization results, which represent the lowest cost water supply systems possible. Additionally, we compare the optimal solutions of the two clusters (also for varying input parameters, such as population densities and slum size development over time) and describe how the result of the optimization depends on the entered remote sensing data. Full article

Reliable information about the spatial distribution of surface waters is critically important in various scientific disciplines. Synthetic Aperture Radar (SAR) is an effective way to detect floods and monitor water bodies over large areas. Sentinel-1 is a new available SAR and its spatial resolution and short temporal baselines have the potential to facilitate the monitoring of surface water changes, which are dynamic in space and time. While several methods and tools for flood detection and surface water extraction already exist, they often comprise a significant manual user interaction and do not specifically target the exploitation of Sentinel-1 data. The existing methods commonly rely on thresholding at the level of individual pixels, ignoring the correlation among nearby pixels. Thus, in this paper, we propose a fully automatic processing chain for rapid flood and surface water mapping with smooth labeling based on Sentinel-1 amplitude data. The method is applied to three different sites submitted to recent flooding events. The quantitative evaluation shows relevant results with overall accuracies of more than 98% and F-measure values ranging from 0.64 to 0.92. These results are encouraging and the first step to proposing operational image chain processing to help end-users quickly map flooding events or surface waters. Full article

Most high-altitude lakes are more sensitive to global warming than the regional atmosphere. However, most existing climate models produce unrealistic surface temperatures on the Tibetan Plateau (TP) lakes, and few studies have focused on the influence of ice surface albedo on high-altitude lakes. Based on field albedo measurements, moderate resolution imaging spectrometer (MODIS) albedo products and numerical simulation, this study evaluates the ice albedo parameterization schemes in existing lake models and investigates the characteristics of the ice surface albedo in six typical TP lakes, as well as the influence of ice albedo error in the FLake model. Compared with observations, several ice albedo schemes all clearly overestimate the lake ice albedo by 0.26 to 0.66, while the average bias of MODIS albedo products is only 0.07. The MODIS-observed albedo of a snow-covered lake varies with the snow proportion, and the lake surface albedo in a snow-free state is approximately 0.15 during the frozen period. The MODIS-observed ice surface (snow-free) albedos are concentrated within the ranges of 0.14–0.16, 0.08–0.10 and 0.10–0.12 in Aksai Chin Lake, Nam Co Lake and Ngoring Lake, respectively. The simulated lake surface temperature is sensitive to variations in lake ice albedo especially in the spring and winter. Full article

A key step in the processing of satellite imagery is the radiometric correction of images to account for reflectance that water vapor, atmospheric dust, and other atmospheric elements add to the images, causing imprecisions in variables of interest estimated at the earth’s surface level. That issue is important when performing spatiotemporal analyses to determine ecosystems’ productivity. In this study, three correction methods were applied to satellite images for the period 2010–2014. These methods were Atmospheric Correction for Flat Terrain 2 (ATCOR2), Fast Line-of-Sight Atmospheric Analysis of Spectral Hypercubes (FLAASH), and Dark Object Substract 1 (DOS1). The images included 12 sub-scenes from the Landsat Thematic Mapper (TM) and the Operational Land Imager (OLI) sensors. The images corresponded to three Permanent Monitoring Sites (PMS) of grasslands, ‘Teseachi’, ‘Eden’, and ‘El Sitio’, located in the state of Chihuahua, Mexico. After the corrections were applied to the images, they were evaluated in terms of their precision for biomass estimation. For that, biomass production was measured during the study period at the three PMS to calibrate production models developed with simple and multiple linear regression (SLR and MLR) techniques. When the estimations were made with MLR, DOS1 obtained an R² of 0.97 (p < 0.05) for 2012 and values greater than 0.70 (p < 0.05) during 2013–2014. The rest of the algorithms did not show significant results and DOS1, which is the simplest algorithm, resulted in the best biomass estimator. Thus, in the multitemporal analysis of grassland based on spectral information, it is not necessary to apply complex correction procedures. The maps of biomass production, elaborated from images corrected with DOS1, can be used as a reference point for the assessment of the grassland condition, as well as to determine the grazing capacity and thus the potential animal production in such ecosystems. Full article

This research focuses on the effect of trampling on vegetation in high-mountain ecosystems through the electromagnetic spectrum’s interaction with plant pigments, cell structure, water content and other substances that have a direct impact on leaf properties. The aim of the study was to confirm with the use of fluorescence methods of variability in the state of high-mountain vegetation previously measured spectrometrically. The most heavily visited part of the High Tatras in Poland was divided into polygons and, after selecting the dominant species within alpine swards, a detailed analysis of trampled and reference patterns was performed. The Analytical Spectral Devices (ASD) FieldSpec 3/4 were used to acquire high-resolution spectral properties of plants, their fluorescence and the leaf chlorophyll content with the difference between the plant surface temperature (ts), and the air temperature (ta) as well as fraction of Absorbed Photosynthetically Active Radiation (fAPAR) used as reference data. The results show that, along tourist trails, vegetation adapts to trampling with the impact depending on the species. A lower chlorophyll value was confirmed by a decrease in fluorescence, and the cellular structures were degraded in trampled compared to reference species, with a lower leaf reflectance. In addition, at the extreme, trampling can eliminate certain species such as Luzula alpino-pilosa, for which significant changes were noted due to trampling. Full article

The correction of spatial scaling bias on the estimate of leaf area index (LAI) retrieved from remotely sensed data is an essential issue in quantitative remote sensing for vegetation monitoring. We analyzed three techniques, including Taylor’s theorem (TT), Wavelet-Fractal technique (WF), and Fractal theory (FT), for correcting the scaling bias of LAI with empirical models in different functions (i.e., power, exponential, logarithmic and polynomial) on both simulated data and a real dataset over a cropland site. The results demonstrated that the scaling bias became greater when the model non-linearity increased. The spatial heterogeneity, which was characterized by the class-specific proportion, the between-class spectral difference and the number of classes within each coarse pixel, was found to be the primary factor in the scaling effect. These factors influenced the scaling effect collectively and existed dependently. With the RMSE less than 0.3 × 10⁻⁶ m²/m², TT was suggested for the correction with a polynomial LAI-NDVI functions. WF was preferred for neighboring scales rather than continuous scales. FT was not recommended for correcting the scaling bias caused by the significant non-linearity in LAI estimation models. This study illustrates the main causes of the scaling effect and provides a reference of technique selection for scaling bias correction to improve the application of remotely sensed estimates. Full article

In object-based image analysis (OBIA), the appropriate parametrization of segmentation algorithms is crucial for obtaining satisfactory image classification results. One of the ways this can be done is by unsupervised segmentation parameter optimization (USPO). A popular USPO method does this through the optimization of a “global score” (GS), which minimizes intrasegment heterogeneity and maximizes intersegment heterogeneity. However, the calculated GS values are sensitive to the minimum and maximum ranges of the candidate segmentations. Previous research proposed the use of fixed minimum/maximum threshold values for the intrasegment/intersegment heterogeneity measures to deal with the sensitivity of user-defined ranges, but the performance of this approach has not been investigated in detail. In the context of a remote sensing very-high-resolution urban application, we show the limitations of the fixed threshold approach, both in a theoretical and applied manner, and instead propose a novel solution to identify the range of candidate segmentations using local regression trend analysis. We found that the proposed approach showed significant improvements over the use of fixed minimum/maximum values, is less subjective than user-defined threshold values and, thus, can be of merit for a fully automated procedure and big data applications. Full article

The objective was to validate the carbon monoxide (CO) total column product inferred from Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) full-mission (2003–2011) short-wave infrared (SWIR) nadir observations using the Beer InfraRed Retrieval Algorithm (BIRRA). Globally distributed Network for the Detection of Atmospheric Composition Change (NDACC) and Total Carbon Column Observing Network (TCCON) ground-based (g-b) measurements were used as a true reference. Weighted averages of SCIAMACHY CO observations within a circle around the g-b observing system were utilized to minimize effects due to spatial mismatch of space-based (s-b) and g-b observations, i.e., disagreements due to representation errors rather than instrument and/or algorithm deficiencies. In addition, temporal weighted averages were examined and then the unweighted (classical) approach was compared to the weighted (non-classical) method. The delivered distance-based filtered SCIAMACHY data were in better agreement with respect to CO averages as compared to square-shaped sampling areas throughout the year. Errors in individual SCIAMACHY retrievals have increased substantially since 2005. The global bias was determined to be in the order of $-10$ parts per billion in volume (ppbv) depending on the reference network and validation strategy used. The largest negative bias was found to occur in the northern mid-latitudes in Europe and North America, and was partly caused by insufficient a priori estimates of CO and cloud shielding. Furthermore, no significant trend was identified in the global bias throughout the mission. The global analysis of the CO columns retrieved by the BIRRA shows results that are largely consistent with similar investigations in previous works. Full article

Machine learning techniques are attractive tools to establish statistical models with a high degree of non linearity. They require a large amount of data to be trained and are therefore particularly suited to analysing remote sensing data. This work is an attempt at using advanced statistical methods of machine learning to predict the bias between Sea Surface Temperature (SST) derived from infrared remote sensing and ground “truth” from drifting buoy measurements. A large dataset of collocation between satellite SST and in situ SST is explored. Four regression models are used: Simple multi-linear regression, Least Square Shrinkage and Selection Operator (LASSO), Generalised Additive Model (GAM) and random forest. In the case of geostationary satellites for which a large number of collocations is available, results show that the random forest model is the best model to predict the systematic errors and it is computationally fast, making it a good candidate for operational processing. It is able to explain nearly 31% of the total variance of the bias (in comparison to about 24% for the multi-linear regression model). Full article

Airborne LiDAR technology is widely used in archaeology and over the past decade has emerged as an accurate tool to describe anthropomorphic landforms. Archaeological features are traditionally emphasised on a LiDAR-derived Digital Terrain Model (DTM) using multiple Visualisation Techniques (VTs), and occasionally aided by automated feature detection or classification techniques. Such an approach offers limited results when applied to heterogeneous structures (different sizes, morphologies), which is often the case for archaeological remains that have been altered throughout the ages. This study proposes to overcome these limitations by developing a multi-scale analysis of topographic position combined with supervised machine learning algorithms (Random Forest). Rather than highlighting individual topographic anomalies, the multi-scalar approach allows archaeological features to be examined not only as individual objects, but within their broader spatial context. This innovative and straightforward method provides two levels of results: a composite image of topographic surface structure and a probability map of the presence of archaeological structures. The method was developed to detect and characterise megalithic funeral structures in the region of Carnac, the Bay of Quiberon, and the Gulf of Morbihan (France), which is currently considered for inclusion on the UNESCO World Heritage List. As a result, known archaeological sites have successfully been geo-referenced with a greater accuracy than before (even when located under dense vegetation) and a ground-check confirmed the identification of a previously unknown Neolithic burial mound in the commune of Carnac. Full article

Phytophthora root rot (PRR) infects the roots of avocado trees, resulting in reduced uptake of water and nutrients, canopy decline, defoliation, and, eventually, tree mortality. Typically, the severity of PRR disease (proportion of canopy decline) is assessed by visually comparing the canopy health of infected trees to a standardised set of photographs and a corresponding disease rating. Although this visual method provides some indication of the spatial variability of PRR disease across orchards, the accuracy and repeatability of the ranking is influenced by the experience of the assessor, the visibility of tree canopies, and the timing of the assessment. This study evaluates two image analysis methods that may serve as surrogates to the visual assessment of canopy decline in large avocado orchards. A smartphone camera was used to collect red, green, and blue (RGB) colour images of individual trees with varying degrees of canopy decline, with the digital photographs then analysed to derive a canopy porosity percentage using a combination of ‘Canny edge detection’ and ‘Otsu’s’ methods. Coinciding with the on-ground measure of canopy porosity, the canopy reflectance characteristics of the sampled trees measured by high resolution Worldview-3 (WV-3) satellite imagery was also correlated against the observed disease severity rankings. Canopy porosity values (ranging from 20–70%) derived from RGB images were found to be significantly different for most disease rankings (p < 0.05) and correlated well (R² = 0.89) with the differentiation of three disease severity levels identified to be optimal. From the WV-3 imagery, a multivariate stepwise regression of 18 structural and pigment-based vegetation indices found the simplified ratio vegetation index (SRVI) to be strongly correlated (R² = 0.96) with the disease rankings of PRR disease severity, with the differentiation of four levels of severity found to be optimal. Full article

The high nutrient concentrations coming from non-point and point pollution have been linked to algae blooms, especially in hydroelectric plant reservoirs, due to higher residence time compared to rivers. The monitoring of algae is important to prevent risk of contamination by toxins in reservoirs used for drinking water supply. In this context, a physical model-based approach was adopted to retrieve chlorophyll-a (chl a) concentration, a photosynthetic pigment found in all phytoplankton species. We assumed that a semi-analytical algorithm parameterized to a eutrophic reservoir could also be applied to other eutrophic reservoirs, at least the specific inherent optical properties (SIOPs) are not similar. The parameterization was carried out based on Ocean and Land Color Instrument (OLCI) bands aboard Sentinel-3 spacecraft. In our study, the semi-analytical approach showed good performance in retrieving chl a content, with a normalized root mean square error (NRMSE) of 18.7%. The findings encourage the use of a unique semi-analytical algorithm in a reservoir cascade, where the impoundments present similar bio-optical status. The good performance of the algorithm indicates that this approach is rather useful in predicting trophic status in reservoirs. Full article

Monitoring changes in forest height, biomass and carbon stock is important for understanding the drivers of forest change, clarifying the geography and magnitude of the fluxes of the global carbon budget and for providing input data to REDD+. The objective of this study was to investigate the feasibility of covering these monitoring needs using InSAR DEM changes over time and associated estimates of forest biomass change and corresponding net CO₂ emissions. A wall-to-wall map of net forest change for Uganda with its tropical forests was derived from two Digital Elevation Model (DEM) datasets, namely the SRTM acquired in 2000 and TanDEM-X acquired around 2012 based on Interferometric SAR (InSAR) and based on the height of the phase center. Errors in the form of bias, as well as parallel lines and belts having a certain height shift in the SRTM DEM were removed, and the penetration difference between X- and C-band SAR into the forest canopy was corrected. On average, we estimated X-band InSAR height to decrease by 7 cm during the period 2000–2012, corresponding to an estimated annual CO₂ emission of 5 Mt for the entirety of Uganda. The uncertainty of this estimate given as a 95% confidence interval was 2.9–7.1 Mt. The presented method has a number of issues that require further research, including the particular SRTM biases and artifact errors; the penetration difference between the X- and C-band; the final height adjustment; and the validity of a linear conversion from InSAR height change to AGB change. However, the results corresponded well to other datasets on forest change and AGB stocks, concerning both their geographical variation and their aggregated values. Full article

The Optimal Estimation (OE) technique is developed within the European Space Agency Climate Change Initiative (ESA-CCI) to retrieve subskin Sea Surface Temperature (SST) from AQUA’s Advanced Microwave Scanning Radiometer—Earth Observing System (AMSR-E). A comprehensive matchup database with drifting buoy observations is used to develop and test the OE setup. It is shown that it is essential to update the first guess atmospheric and oceanic state variables and to perform several iterations to reach an optimal retrieval. The optimal number of iterations is typically three to four in the current setup. In addition, updating the forward model, using a multivariate regression model is shown to improve the capability of the forward model to reproduce the observations. The average sensitivity of the OE retrieval is 0.5 and shows a latitudinal dependency with smaller sensitivity for cold waters and larger sensitivity for warmer waters. The OE SSTs are evaluated against drifting buoy measurements during 2010. The results show an average difference of 0.02 K with a standard deviation of 0.47 K when considering the 64% matchups, where the simulated and observed brightness temperatures are most consistent. The corresponding mean uncertainty is estimated to 0.48 K including the in situ and sampling uncertainties. An independent validation against Argo observations from 2009 to 2011 shows an average difference of 0.01 K, a standard deviation of 0.50 K and a mean uncertainty of 0.47 K, when considering the best 62% of retrievals. The satellite versus in situ discrepancies are highest in the dynamic oceanic regions due to the large satellite footprint size and the associated sampling effects. Uncertainty estimates are available for all retrievals and have been validated to be accurate. They can thus be used to obtain very good retrieval results. In general, the results from the OE retrieval are very encouraging and demonstrate that passive microwave observations provide a valuable alternative to infrared satellite observations for retrieving SST. Full article

An ice-free Arctic summer would have pronounced impacts on global climate, coastal habitats, national security, and the shipping industry. Rapid and accelerated Arctic sea ice loss has placed the reality of an ice-free Arctic summer even closer to the present day. Accurate projection of the first Arctic ice-free summer year is extremely important for business planning and climate change mitigation, but the projection can be affected by many factors. Using an inter-calibrated satellite sea ice product, this article examines the sensitivity of decadal trends of Arctic sea ice extent and statistical projections of the first occurrence of an ice-free Arctic summer. The projection based on the linear trend of the last 20 years of data places the first Arctic ice-free summer year at 2036, 12 years earlier compared to that of the trend over the last 30 years. The results from a sensitivity analysis of six commonly used curve-fitting models show that the projected timings of the first Arctic ice-free summer year tend to be earlier for exponential, Gompertz, quadratic, and linear with lag fittings, and later for linear and log fittings. Projections of the first Arctic ice-free summer year by all six statistical models appear to converge to the 2037 ± 6 timeframe, with a spread of 17 years, and the earliest first ice-free Arctic summer year at 2031. Full article

The correlation coefficient (CC) was substantiated as a simple, yet robust statistical tool in the quality assessment of hyperspectral imaging (HSI) data. The sensitivity of the metric was also characterized with respect to artificially-induced errors. The CC was found to be sensitive to spectral shifts and single feature modifications in hyperspectral ground data despite the high, artificially-induced, signal-to-noise ratio (SNR) of 100:1. The study evaluated eight airborne hyperspectral images that varied in acquisition spectrometer, acquisition date and processing methodology. For each image, we identified a uniform ground target region of interest (ROI) that was comprised of a single asphalt road pixel from each column within the sensor field-of-view (FOV). A CC was calculated between the spectra from each of the pixels in the ROI and the data from the center pixel. Potential errors were located by reductions in the CCs below a designated threshold, which was derived from the results of the sensitivity tests. The spectral range associated with each error was established using a windowing technique where the CCs were recalculated after removing the spectral data within various windows. Errors were isolated in the spectral window that removed the previously-identified reductions in the CCs. Finer errors were detected by calculating the CCs across the ROI in the spectral range surrounding various atmospheric absorption features. Despite only observing deviations in the CCs from the 3rd–6th decimal places, non-trivial errors were detected in the imagery. An error was detected within a single band of the shortwave infrared imagery. Errors were also observed throughout the visible-near-infrared imagery, especially in the blue end. With this methodology, it was possible to immediately gauge the spectral consistency of the HSI data across the FOV. Consequently, the effectiveness of various processing methodologies and the spectral consistency of the imaging spectrometers themselves could be studied. Overall, the research highlights the utility of the CC as a simple, low monetary cost, analytical tool for the localization of errors in spectroscopic imaging data. Full article

We developed a new strategy for a synchronous generation of real-time (RT) and near real-time (NRT) tropospheric products. It exploits the precise point positioning method with Kalman filtering and backward smoothing, both supported by real-time orbit and clock products. The strategy can be optimized for the latency or the accuracy of NRT production. In terms of precision, it is comparable to the traditional NRT network solution using deterministic models in the least-square adjustment. Both RT and NRT solutions provide a consistent set of tropospheric parameters such as zenith total delays, horizontal tropospheric gradients and slant delays, all with a high resolution and optimally exploiting all observations from available GNSS multi-constellations. As the new strategy exploits RT processing, we assessed publicly precise RT products and results of RT troposphere monitoring. The backward smoothing applied for NRT solution, when using an optimal latency of 30 min, reached an improvement of 20% when compared to RT products. Additionally, multi-GNSS solutions provided more accurate (by 25%) tropospheric parameters, and the impact will further increase when constellations are complete and supported with precise models and products. The new strategy is ready to replace our NRT contribution to the EUMETNET EIG GNSS Water Vapour Programme (E-GVAP) and effectively support all modern multi-GNSS tropospheric products. Full article

Due to the rich physical meaning of aurora morphology, the classification of aurora images is an important task for polar scientific expeditions. However, the traditional classification methods do not make full use of the different features of aurora images, and the dimension of the description features is usually so high that it reduces the efficiency. In this paper, through combining multiple features extracted from aurora images, an aurora image classification method based on multi-feature latent Dirichlet allocation (AI-MFLDA) is proposed. Different types of features, whether local or global, discrete or continuous, can be integrated after being transformed to one-dimensional (1-D) histograms, and the dimension of the description features can be reduced due to using only a few topics to represent the aurora images. In the experiments, according to the classification system provided by the Polar Research Institute of China, a four-class aurora image dataset was tested and three types of features (MeanStd, scale-invariant feature transform (SIFT), and shape-based invariant texture index (SITI)) were utilized. The experimental results showed that, compared to the traditional methods, the proposed AI-MFLDA is able to achieve a better performance with 98.2% average classification accuracy while maintaining a low feature dimension. Full article

Aside from the well-known applications (positioning, navigation and timing) brought by Global Navigation Satellite System (GNSS), reconstruction of tropospheric atmosphere distribution information using tomography technique based on the multi-GNSS observations has been developed as a research point in the fields of GNSS Meteorology. In this paper, an optimal tropospheric tomography method using observations from multi-GNSS (Global Navigation Satellite System) is proposed, which considers the reasonable weightings of observation equations derived from multi-GNSS as well as the various constraints. Comparing to the equal weighting strategy of multi-GNSS observations for the previously multi-GNSS tomography studies, the proposed method in this paper has the ability to tune the weightings for a different type of equations. Experiments show that the proposed method can improve the internal/external accuracy of GNSS tomography modeling with the GNSS precise point positioning (PPP)-estimated slant wet delay as reference when compared to the conventional method. In addition, the data derived from radiosonde is used as an external testing, and the result also expresses the superiority of the proposed method when compared to the conventional method. Full article

The p-wave double-frequency (DF) microseisms generated by super typhoon Lupit (14–26 October 2009) over the western Pacific Ocean were detected by an on-land seismological array deployed in Northeastern China. We applied a frequency-domain beamforming method to investigate their source regions. Comparing with the best-track data and satellite observations, the located source regions of the p-wave DF microseisms, which corresponded to the strongest ocean wave–wave interactions, were found to be comparable to the typhoon centers in the microseismic frequency band of ~0.18–0.21 Hz. The p-wave DF microseisms were probably excited by the nonlinear interaction of ocean waves generated by the typhoon at different times, in good agreement with the Longuet–Higgins theory for the generation of DF microseisms. The localization deviation, which was ~120 km for typhoon Lupit in this study, might depend on the speed and direction of typhoon movement, the geometry of the seismological array, and the heterogeneity of the solid Earth structure. The p-wave DF microseisms generated in coastal source regions were also observed in the beamformer outputs, but with relatively lower dominant frequency band of ~0.14–0.16 Hz. These observations show that the p-wave DF microseisms generated near typhoon centers could be used as a seismic remote sensing proxy to locate and track typhoons over the oceans from under water in a near-real-time and continuous manner. Full article

Sensitivity to weather conditions, and specially to clouds, is a severe limiting factor to the use of optical remote sensing for Earth monitoring applications. A possible alternative is to benefit from weather-insensitive synthetic aperture radar (SAR) images. In many real-world applications, critical decisions are made based on some informative optical or radar features related to items such as water, vegetation or soil. Under cloudy conditions, however, optical-based features are not available, and they are commonly reconstructed through linear interpolation between data available at temporally-close time instants. In this work, we propose to estimate missing optical features through data fusion and deep-learning. Several sources of information are taken into account—optical sequences, SAR sequences, digital elevation model—so as to exploit both temporal and cross-sensor dependencies. Based on these data and a tiny cloud-free fraction of the target image, a compact convolutional neural network (CNN) is trained to perform the desired estimation. To validate the proposed approach, we focus on the estimation of the normalized difference vegetation index (NDVI), using coupled Sentinel-1 and Sentinel-2 time-series acquired over an agricultural region of Burkina Faso from May–November 2016. Several fusion schemes are considered, causal and non-causal, single-sensor or joint-sensor, corresponding to different operating conditions. Experimental results are very promising, showing a significant gain over baseline methods according to all performance indicators. Full article

Nowadays, satellite images are considered as one of the most relevant sources of information in the context of major disasters management. Their availability in extreme weather conditions and their ability to cover wide geographic areas make them an indispensable tool toward an effective disaster response. Among the various available sensors, Synthetic Aperture Radar (SAR) is distinguished in the context of flood management by its ability to penetrate cloud cover and its robustness to unfavourable weather conditions. This work aims at developing a new technique for flooded areas extraction from high resolution time-series SAR images. The proposed approach is mainly based on three steps: first, homogeneous regions characterizing water surfaces are extracted from each SAR image using a local texture descriptor. Then, mathematical morphology is applied to filter tiny artifacts and small homogeneous areas present in the image. And finally, spatial and radiometric information embedded in each pixel are extracted and are fused with the same pixel information but from another image to decide if the current pixel belongs to a flooded region. In order to assess the performance of the proposed algorithm, our methodology was applied to time-series images acquired before and during three different flooding events: (1) Richelieu River and lake Champlain floods, Quebec, Canada in 2011; (2) Evros River floods, Greece in 2014 and (3) Western and southwestern of Iran floods in 2016. Experiments show that our approach gives very promising results compared to existing techniques. Full article

Two remote sensing methods, satellite altimetry and 4D-Var assimilation of satellite imagery, are used to compute surface velocity fields in the Black Sea region. Surface currents derived from the two methods are compared for several cases with intense mesoscale and large-scale dynamics during low wind conditions. Comparison shows that the obtained results coincide well quantitatively and qualitatively. However, satellite imagery provides more reasonable results on the spatial variability of coastal dynamics than altimetry data. In particular, this is related to the reconstruction of eddy coastal dynamics, such as Black Sea near-shore anticyclones. Current streamlines in these eddies are not closed near the coast in altimetry data, which we relate to the extrapolation during mapping procedure in the absence of coastal along-track measurements. On the other hand, in offshore areas, imagery-derived currents can be underestimated due to the absence of thermal contrasts and smoothing during the procedure of the 4D-Var assimilation. Wind drift currents are another source of inconsistency, as their impact is directly observed in satellite imagery but absent in altimetry data. The advantage of the 4D-Var method for reconstructing coastal dynamics is used to compute surface currents in the Marmara Sea on the base of 250 m resolution Modis optical data. The results reveal the very complex dynamics of the basin, with a large number of mesoscale and sub-mesoscale eddies. 4D-Var assimilation of Modis imagery is used to obtain information about dynamic characteristics of these small eddies with radiuses of 4–10 km. Full article