Remote Sens., Volume 12, Issue 18 (September-2 2020) – 245 articles

Digital and scalable technologies are increasingly important for rapid and large-scale assessment and monitoring of land cover change. Until recently, little research has existed on how these technologies can be specifically applied to the monitoring of Reducing Emissions from Deforestation and Forest Degradation (REDD+) activities. Using the Google Earth Engine (GEE) cloud computing platform, we applied the recently developed phenology-based threshold classification method (PBTC) for detecting and mapping forest cover and carbon stock changes in Siem Reap province, Cambodia, between 1990 and 2018. The obtained PBTC maps were validated using Google Earth high resolution historical imagery and reference land cover maps by creating 3771 systematic 5 × 5 km spatial accuracy points. The overall cumulative accuracy of this study was 92.1% and its cumulative Kappa was 0.9, which are sufficiently high to apply the PBTC method to detect forest land cover change. Accordingly, we estimated the carbon stock changes over a 28-year period in accordance with the Good Practice Guidelines of the Intergovernmental Panel on Climate Change. We found that 322,694 ha of forest cover was lost in Siem Reap, representing an annual deforestation rate of 1.3% between 1990 and 2018. This loss of forest cover was responsible for carbon emissions of 143,729,440 MgCO₂ over the same period. If REDD+ activities are implemented during the implementation period of the Paris Climate Agreement between 2020 and 2030, about 8,256,746 MgCO₂ of carbon emissions could be reduced, equivalent to about USD 6-115 million annually depending on chosen carbon prices. Our case study demonstrates that the GEE and PBTC method can be used to detect and monitor forest cover change and carbon stock changes in the tropics with high accuracy. Full article

Although vegetation phenology thresholds have been developed for a wide range of mapping applications, their use for assessing the distribution of natural bamboo and the related carbon stocks is still limited, especially in Southeast Asia. Here, we used Google Earth Engine (GEE) to collect time-series of Landsat 8 Operational Land Imager (OLI) and Sentinel-2 images and employed a phenology-based threshold classification method (PBTC) to map the natural bamboo distribution and estimate carbon stocks in Siem Reap Province, Cambodia. We processed 337 collections of Landsat 8 OLI for phenological assessment and generated 121 phenological profiles of the average vegetation index for three vegetation land cover categories from 2015 to 2018. After determining the minimum and maximum threshold values for bamboo during the leaf-shedding phenology stage, the PBTC method was applied to produce a seasonal composite enhanced vegetation index (EVI) for Landsat collections and assess the bamboo distributions in 2015 and 2018. Bamboo distributions in 2019 were then mapped by applying the EVI phenological threshold values for 10 m resolution Sentinel-2 satellite imagery by accessing 442 tiles. The overall Landsat 8 OLI bamboo maps for 2015 and 2018 had user’s accuracies (UAs) of 86.6% and 87.9% and producer’s accuracies (PAs) of 95.7% and 97.8%, respectively, and a UA of 86.5% and PA of 91.7% were obtained from Sentinel-2 imagery for 2019. Accordingly, carbon stocks of natural bamboo by district in Siem Reap at the province level were estimated. Emission reductions from the protection of natural bamboo can be used to offset 6% of the carbon emissions from tourists who visit this tourism-destination province. It is concluded that a combination of GEE and PBTC and the increasing availability of remote sensing data make it possible to map the natural distribution of bamboo and carbon stocks. Full article

This study investigates coastal erosion, revetment as a shoreline protection strategy, and human perceptions of revetments in the Lower Meghna estuary of the Bangladesh where new revetments were recently constructed. Questions addressed were: (1) How do rates of shoreline change vary over the period 2011–2019? (2) Did new revetments effectively halt erosion and what were the magnitudes of erosion change? (3) How have erosion rates changed for shorelines within 1 km of revetments, and (4) How do households perceive revetments? High-resolution Planet Lab imagery was used to quantify shoreline change rates. Analysis of household survey data assessed human perceptions of the revetment’s desirability and efficacy. Results revealed high rates of erosion for 2011–2019 with declining erosion after 2013. New revetments effectively halted erosion for protected shorelines. Significant spatial trends for erosion rates existed for shorelines adjacent to revetments. Survey respondents overwhelmingly had positive attitudes about a desire for revetment protection; however, upstream respondents expressed a strong majority perception that revetment acts to make erosion worse. Highlights of the research include integration of remote sensing with social science methods, the timing of the social survey shortly after revetment construction, and results showing significant erosion change upstream and downstream of new revetments. Full article

Synthetic aperture radar (SAR) altimeters represent a new method of microwave remote sensing for ocean wave observations. The adoption of SAR technology in the azimuthal direction has the advantage of a high resolution. The Sentinel-3 altimeter is the first radar altimeter to acquire global observations in SAR mode; hence, the data quality needs to be assessed before extensively applying these data. The European Space Agency (ESA) evaluates the Sentinel-3 accuracy on a global scale but has yet to perform a detailed analysis in terms of different offshore distances and different water depths. In this paper, Sentinel-3 and Jason-2 significant wave height (SWH) data are matched in both time and space with buoy data from the United States East and West Coasts and the Central Pacific Ocean. The Sentinel-3 SWH data quality is evaluated according to different offshore distances and water depths in comparison with Jason-2 SWH data. In areas more than 50 km offshore, the Sentinel-3 SWH accuracy is generally high and less affected by the water depth and sea conditions (root-mean-square error of 0.28 m and correlation coefficient of 0.98); in areas less than 50 km offshore, the SWH data accuracy is slightly affected by water depth and sea conditions (especially the former). Compared with Jason-2, the observation ability of the Sentinel-3 altimeter in nearshore areas with water depths of 0 m-500 m is greatly improved, but in some deep water areas with stable sea conditions, the Jason-2 SWH data accuracy is higher than that of Sentinel-3. This work provides a reference for the refined application of Sentinel-3 SWH data in offshore deep water areas and nearshore shallow water areas. Full article

The semantic segmentation of underwater imagery is an important step in the ecological analysis of coral habitats. To date, scientists produce fine-scale area annotations manually, an exceptionally time-consuming task that could be efficiently automatized by modern CNNs. This paper extends our previous work presented at the 3DUW’19 conference, outlining the workflow for the automated annotation of imagery from the first step of dataset preparation, to the last step of prediction reassembly. In particular, we propose an ecologically inspired strategy for an efficient dataset partition, an over-sampling methodology targeted on ortho-imagery, and a score fusion strategy. We also investigate the use of different loss functions in the optimization of a Deeplab V3+ model, to mitigate the class-imbalance problem and improve prediction accuracy on coral instance boundaries. The experimental results demonstrate the effectiveness of the ecologically inspired split in improving model performance, and quantify the advantages and limitations of the proposed over-sampling strategy. The extensive comparison of the loss functions gives numerous insights on the segmentation task; the Focal Tversky, typically used in the context of medical imaging (but not in remote sensing), results in the most convenient choice. By improving the accuracy of automated ortho image processing, the results presented here promise to meet the fundamental challenge of increasing the spatial and temporal scale of coral reef research, allowing researchers greater predictive ability to better manage coral reef resilience in the context of a changing environment. Full article

This study is a follow-up of a full methodology for the homogenisation and harmonisation of the two Ocean and Land Colour Instrument (OLCI) payloads based on the OLCI-A/OLCI-B tandem phase analysis. This analysis provided cross-calibration factors between the two instruments with a very high precision, providing a ‘truth’ from the direct comparison of simultaneous and collocated acquisitions. The long-term monitoring of such cross-calibration is a prerequisite for an operational application of sensors harmonisation along the mission lifetime, no other tandem phase between OLCI-A and OLCI-B being foreseen due to the cost of such operation. This article presents a novel approach for the monitoring of the OLCI radiometry based on statistics of Deep Convective Clouds (DCC) observations, especially dedicated to accurately monitor the full across-track dependency of the cross-calibration of OLCI-A and OLCI-B. Specifically, the inflexion point of DCC reflectance distributions is used as an indicator of the absolute calibration for each subdivision of the OLCI Field-of-View. This inflexion point is shown to provide better precision than the mode of the distributions which is commonly used in the community. Excess of saturation in OLCI-A high radiances is handled through the analysis of interband relationships between impacted channels and reference channels that are not impacted by saturation. Such analysis also provides efficient insights on the variability of the target’s response as well as on the evolution of the interband calibration of each payload. First, cross-calibration factors obtained over the tandem period allows to develop and validate the approach, notably for the handling of the saturated pixels, based on the comparison with the ‘truth’ obtained from the tandem analysis. Factors obtained out of (and far from) the tandem period then provides evidence that the cross-calibration reported over the tandem period (1–2% bias between the instruments) as well as inter-camera calibration residuals persist with very similar proportions, to the exception of the 400 nm channel and with slightly less precision for the 1020 nm channel. For all OLCI channels, relative differences between the cross-calibration factors obtained from the tandem analysis and the factors obtained over the other period are below 1% from a monthly analysis, even below 0.5% from a multi-monthly analysis). This opens the way not only to an accurate long-term monitoring of the OLCI radiometry but also, and precisely targeted for this study, to the monitoring of the cross-calibration of the two sensors over the mission lifetime. It also provides complementary information to the tandem analysis as the calibration indicators are traced individually for each sensor across-track, confirming and quantifying inter-camera radiometric biases, independently for both sensors. Assumptions used in this study are discussed and validated, also providing a framework for the adaptation of the presented methodology to other optical sensors. Full article

Chlorophyll is an essential pigment for photosynthesis in crops, and leaf chlorophyll content can be used as an indicator for crop growth status and help guide nitrogen fertilizer applications. Estimating crop chlorophyll content plays an important role in precision agriculture. In this study, a variable, rate of change in reflectance between wavelengths ‘a’ and ‘b’ (RCRW_a-b), derived from in situ hyperspectral remote sensing data combined with four advanced machine learning techniques, Gaussian process regression (GPR), random forest regression (RFR), support vector regression (SVR), and gradient boosting regression tree (GBRT), were used to estimate the chlorophyll content (measured by a portable soil–plant analysis development meter) of rice. The performances of the four machine learning models were assessed and compared using root mean square error (RMSE), mean absolute error (MAE), and coefficient of determination (R²). The results revealed that four features of RCRW_a-b, RCRW_{551.0–565.6}, RCRW_{739.5–743.5}, RCRW_{684.4–687.1} and RCRW_{667.9–672.0}, were effective in estimating the chlorophyll content of rice, and the RFR model generated the highest prediction accuracy (training set: RMSE = 1.54, MAE =1.23 and R² = 0.95; validation set: RMSE = 2.64, MAE = 1.99 and R² = 0.80). The GPR model was found to have the strongest generalization (training set: RMSE = 2.83, MAE = 2.16 and R² = 0.77; validation set: RMSE = 2.97, MAE = 2.30 and R² = 0.76). We conclude that RCRW_a-b is a useful variable to estimate chlorophyll content of rice, and RFR and GPR are powerful machine learning algorithms for estimating the chlorophyll content of rice. Full article

This study aimed to assess the ability of using visible and near-infrared reflectance (Vis–NIR) spectroscopy to quantify soil erodibility factor (K) rapidly in an ecologically restored watershed. To achieve this goal, we explored the performance and transferability of the developed spectral models in multiple land-use types: woodland, shrubland, terrace, and slope farmland (the first two types are natural land and the latter two are cultivated land). Subsequently, we developed an improved approach by combining spectral data with related topographic variables (i.e., elevation, watershed location, slope height, and normalized height) to estimate K. The results indicate that the calibrated spectral model using total samples could estimate K factor effectively (R²_CV = 0.71, RMSE_CV = 0.0030 Mg h Mj⁻¹ mm⁻¹, and RPD_CV = 1.84). When predicting K in the new samples, models performed well in natural land soils (R²_P = 0.74, RPD_P = 1.93) but failed in cultivated land soils (R²_P = 0.24, RPD_P = 0.99). Furthermore, the developed models showed low transferability between the natural and cultivated land datasets. The results also indicate that the combination of spectral data with topographic variables could slightly increase the accuracies of K estimation in total and natural land datasets but did not work for cultivated land samples. This study demonstrated that the Vis–NIR spectroscopy could be used as an effective method in predicting K. However, the predictability and transferability of the calibrated models were land-use type dependent. Our study also revealed that the coupling of spectrum and environmental variable is an effective improvement of K estimation in natural landscape region. Full article

In August 2018, a group of experts working with terrestrial/marine geophysics and remote sensing methods to explore archaeological sites in Denmark, Finland, Norway, Scotland and Sweden gathered together for the first time at the Workshop ‘Sensing Archaeology in The North’. The goal was to exchange experiences, discuss challenges, and consider future directions for further developing these methods and strategies for their use in archaeology. After the event, this special journal issue was arranged to publish papers that are based on the workshop presentations, but also to incorporate work that is produced by other researchers in the field. This paper closes the special issue and further aims to provide current state-of-the-art for the methods represented by the workshop. Here, we introduce the aspects that inspired the organisation of the meeting, a summary of the 12 presentations and eight paper contributions, as well as a discussion about the main outcomes of the workshop roundtables, including the production of two searchable databases (online resources and equipment). We conclude with the position that the ‘North’, together with its unique cultural heritage and thriving research community, is at the forefront of good practice in the application and development of sensing methods in archaeological research and management. However, further method development is required, so we claim the support of funding bodies to back research efforts based on testing/experimental studies to: explore unknown survey environments and identify optimal survey conditions, as well as to monitor the preservation of archaeological remains, especially those that are at risk. It is demonstrated that remote sensing and geophysics not only have an important role in the safeguarding of archaeological sites from development and within prehistorical-historical research, but the methods can be especially useful in recording and monitoring the increased impact of climate change on sites in the North. Full article

Snow surface spectral reflectance is very important in the Earth’s climate system. Traditional land surface models with parameterized schemes can simulate broadband snow surface albedo but cannot accurately simulate snow surface spectral reflectance with continuous and fine spectral wavebands, which constitute the major observations of current satellite sensors; consequently, there is an obvious gap between land surface model simulations and remote sensing observations. Here, we suggest a new integrated scheme that couples a radiative transfer model with a land surface model to simulate high spectral resolution snow surface reflectance information specifically targeting multisource satellite remote sensing observations. Our results indicate that the new integrated model can accurately simulate snow surface reflectance information over a large spatial scale and continuous time series. The integrated model extends the range of snow spectral reflectance simulation to the whole shortwave band and can predict snow spectral reflectance changes in the solar spectrum region based on meteorological element data. The kappa coefficients (K) of both the narrowband snow albedo targeting Moderate Resolution Imaging Spectroradiometer (MODIS) data simulated by the new integrated model and the retrieved snow albedo based on MODIS reflectance data are 0.5, and both exhibit good spatial consistency. Our proposed narrowband snow albedo simulation scheme targeting satellite remote sensing observations is consistent with remote sensing satellite observations in time series and can predict narrowband snow albedo even during periods of missing remote sensing observations. This new integrated model is a significant improvement over traditional land surface models for the direct spectral observations of satellite remote sensing. The proposed model could contribute to the effective combination of snow surface reflectance information from multisource remote sensing observations with land surface models. Full article

Synthetic aperture radar (SAR) tomography (TomoSAR) is able to separate multiple scatterers layovered inside the same resolution cell in high-resolution SAR images of urban scenarios, usually with a large number of orbits, making it an expensive and unfeasible task for many practical applications. Targeting at finding out the minimum number of images necessary for tomographic reconstruction, this paper innovatively applies minimum redundancy array (MRA) for tomographic baseline array optimization. Monte Carlo simulations are conducted by means of Two-step Iterative Shrinkage/Thresholding (TWIST) and Truncated Singular Value Decomposition (TSVD) to fully evaluate the tomographic performance of MRA orbits in terms of detection rates, Cramer Rao Lower Bounds, as well as resistance against sidelobes. Experiments on COSMO-SkyMed and TerraSAR-X/TanDEM-X data are also conducted in this paper. The results from simulations and experiments on real data have both demonstrated that introducing MRA for baseline optimization in SAR tomography can benefit from the dramatic reduction of necessary orbit numbers, if the recently proposed TWIST method is used for tomographic reconstruction. Although the simulation and experiments in this manuscript are carried out using spaceborne data, the outcome of this paper can also give examples for airborne TomoSAR when designing flight orbits using airborne sensors. Full article

Due to a limited number of monitoring stations in Western Africa, the impact of mineral dust on PM10 surface concentrations is still poorly known. We propose a new method to retrieve PM10 dust surface concentrations from sun photometer aerosol optical depth (AOD) and CALIPSO/CALIOP Level 2 aerosol layer products. The method is based on a multi linear regression model that is trained using co-located PM10, AERONET and CALIOP observations at 3 different locations in the Sahel. In addition to the sun photometer AOD, the regression model uses the CALIOP-derived base and top altitude of the lowermost dust layer, its AOD, the columnar total and columnar dust AOD. Due to the low revisit period of the CALIPSO satellite, the monthly mean annual cycles of the parameters are used as predictor variables rather than instantaneous observations. The regression model improves the correlation coefficient between monthly mean PM10 and AOD from 0.15 (AERONET AOD only) to 0.75 (AERONET AOD and CALIOP parameters). The respective high and low PM10 concentration during the winter dry season and summer season are well produced. Days with surface PM10 above 100 μg/m³ are better identified when using the CALIOP parameters in the multi linear regression model. The number of true positives (actual and predicted concentrations above the threshold) is increased and leads to an improvement in the classification sensitivity (recall) by a factor 1.8. Our methodology can be extrapolated to the whole Sahel area provided that satellite derived AOD maps are used in order to create a new dataset on population exposure to dust events in this area. Full article

L-band brightness temperature ($T_b$) is one of the key remotely-sensed variables that provides information regarding surface soil moisture conditions. In order to harness the information in $T_b$ observations, a radiative transfer model (RTM) is investigated for eventual inclusion into a data assimilation framework. In this study, $T_b$ estimates from the RTM implemented in the NASA Goddard Earth Observing System (GEOS) were evaluated against the nearly four-year record of daily $T_b$ observations collected by L-band radiometers onboard the Aquarius satellite. Statistics between the modeled and observed $T_b$ were computed over North America as a function of soil hydraulic properties and vegetation types. Overall, statistics showed good agreement between the modeled and observed $T_b$ with a relatively low, domain-average bias (0.79 K (ascending) and −2.79 K (descending)), root mean squared error (11.0 K (ascending) and 11.7 K (descending)), and unbiased root mean squared error (8.14 K (ascending) and 8.28 K (descending)). In terms of soil hydraulic parameters, large porosity and large wilting point both lead to high uncertainty in modeled $T_b$ due to the large variability in dielectric constant and surface roughness used by the RTM. The performance of the RTM as a function of vegetation type suggests better agreement in regions with broadleaf deciduous and needleleaf forests while grassland regions exhibited the worst accuracy amongst the five different vegetation types. Full article

The purpose of this work is to perform the joint design of a classification system including both a radar sensor and an image processing software. Experimental data were generated with a three-dimensional scanner. The criterion which rules the design is the false recognition rate, which should be as small as possible. The classifier involved is support vector machines, combined with an error correcting code. We apply the proposed method to optimize security check. For this purpose we retain eight relevant parameters which impact the recognition performances. To estimate the best parameters, we adapt our adaptive mixed grey wolf algorithm. This is a computational technique inspired by nature to minimize a criterion. Our adaptive mixed grey wolf algorithmwas found to outperform comparative methods in terms of computational load on simulations and with real-world data. Full article

One of the key sources of uncertainties in sea ice freeboard and thickness estimates derived from satellite radar altimetry results from changes in sea ice surface properties. In this study, we analyse this effect, comparing upward-looking sonar (ULS) measurements in the Beaufort Sea over the period 2003–2018 to sea ice draft derived from Envisat and Cryosat-2 data. We show that the sea ice draft growth underestimation observed for the most of winter seasons depends on the surface properties preconditioned by the melt intensity during the preceding summer. The comparison of sea ice draft time series in the Cryosat-2 era indicates that applying 50% retracker thresholds, used to produce the European Space Agency’s Climate Change Initiative (CCI) product, provide better agreement between satellite retrievals and ULS data than the 80% threshold that is closer to the expected physical waveform interpretation. Our results, therefore, indicate compensating error contributions in the full end-to-end sea-ice thickness processing chain, which prevents the quantification of individual factors with sea-ice thickness/draft validation data alone. Full article

Image segmentation is a cost-effective way to obtain information about the sizes and structural composition of agricultural parcels in an area. To accurately obtain such information, the parameters of the segmentation algorithm ought to be optimized using supervised or unsupervised methods. The difficulty in obtaining reference data makes unsupervised methods indispensable. In this study, we evaluated an existing unsupervised evaluation metric that minimizes a global score (GS), which is computed by summing up the intra-segment uniformity and inter-segment dissimilarity within a segmentation output. We modified this metric and proposed a new metric that uses absolute difference to compute the GS. We compared this proposed metric with the existing metric in two optimization approaches based on the Multiresolution Segmentation (MRS) algorithm to optimally delineate agricultural parcels from Sentinel-2 images in Lower Saxony, Germany. The first approach searches for optimal scale while keeping shape and compactness constant, while the second approach uses Bayesian optimization to optimize the three main parameters of the MRS algorithm. Based on a reference data of agricultural parcels, the optimal segmentation result of each optimization approach was evaluated by calculating the quality rate, over-segmentation, and under-segmentation. For both approaches, our proposed metric outperformed the existing metric in different agricultural landscapes. The proposed metric identified optimal segmentations that were less under-segmented compared to the existing metric. A comparison of the optimal segmentation results obtained in this study to existing benchmark results generated via supervised optimization showed that the unsupervised Bayesian optimization approach based on our proposed metric can potentially be used as an alternative to supervised optimization, particularly in geographic regions where reference data is unavailable or an automated evaluation system is sought. Full article

Understanding the spatiotemporal dynamics of bamboo forests is of critical importance as it characterizes the interaction between forest and agricultural ecosystems and provides essential information for sustainable ecosystem management and decision-making. Thus far, the specific dynamics of moso bamboo (Phyllostachys edulis) and other bamboo are still unknown. In this study, we used temporal information extracted from Landsat time series stacks with Google Earth Engine (GEE) to characterize the spatiotemporal dynamics of bamboo forests, including moso bamboo and other bamboo, in Lin’an County, China, from 2000 to 2019. The bamboo forests were mapped in four periods: the early 2000s (2000–2004), the late 2000s (2005–2009), the early 2010s (2010–2014), and the late 2010s (2015–2019). The overall accuracy of these maps ranged from 97% to 99%. We then analyzed the spatiotemporal dynamics of the bamboo forests at the county and subdistrict/township scales, and probed the bamboo forest gain and loss with respect to the terrain features. Our findings show that bamboo forests increased by 4% from 2000 to 2014, followed by a sharp decrease of 13% in the late 2010s. The decrease was mainly caused by the loss of other bamboo. Approximately 69% of the bamboo forest gain occurred in non-bamboo forest areas, and the rest occupied non-forest areas. Bamboo forest loss was mainly due to conversion into orchard (59%) and forest plantation (22%). Compared to bamboo forest gain, bamboo forest loss was typically observed in areas with lower elevations and steeper slopes. Our study offers spatially explicit and timely insight into bamboo forest changes at the regional scale. The derived maps can be applied to study the drivers, consequences, and future trends of bamboo forest dynamics, which will contribute to sustainable ecosystem management. Full article

A cavity ring-down spectroscopy (CRDS) G-2401m analyzer onboard a Beechcraft King Air 350, a new Korean Meteorological Administration (KMA) research aircraft measurement platform since 2018, has been used to measure in situ CO₂, CH₄, and CO. We analyzed the aircraft measurements obtained in two campaigns: a within-boundary layer survey over the western Republic of Korea (hereafter Korea) for analyzing the CO₂ and CH₄ emission characteristics for each season (the climate change monitoring (CM) CM mission), and a low altitude survey over the Yellow Sea for monitoring the pollutant plumes transported into Korea from China (the environment monitoring (EM) mission). This study analyzed CO₂, CH₄, and CO data from a total of 14 flights during 2019 season. To characterize the regional combustion sources signatures of CO₂ and CH₄, we calculated the short-term (1-min slope based on one second data) regression slope of CO to CO₂ and CH₄ to CO enhancements (subtracted with background level, present as ∆CO, ∆CO₂, and ∆CH₄); slope filtered with correlation coefficients (R²) (<0.4 were ignored). These short-term slope analyses seem to be sensitive to aircraft measurements in which the instrument samples short-time varying mixtures of different air masses. The EM missions all of which were affected by pollutants emitted in China, show the regression slope between ∆CO and ∆CO₂ with of 1.8–6% and 0.3–0.7 between ∆CH₄ and ∆CO. In particular, the regression slope between ∆CO and ∆CO₂ increased to >4% when air flows from east-central China such as Hebei, Shandong, and Jiangsu provinces, etc., sustained for 1–3 days, suggesting pollutants from these regions were most likely characterized by incomplete fossil fuel combustions at the industries. Over 80% of the observations in the Western Korea missions were attributed to Korean emission sources with regression slope between ∆CO and ∆CO₂ of 0.5–1.9%. The CO₂ emissions hotspots were mainly located in the north-Western Korea of high population density and industrial activities. The higher CH₄ were observed during summer season with the increasing concentration of approximately 6% over the background level, it seems to be attributed to biogenic sources such as rice paddies, landfill, livestock, and so on. It is also noted that occurrences of high pollution episodes in North-Western Korea are more closely related to the emissions in China than in Korea. Full article

Species identification in Quebec, Canada, is usually performed with photo-interpretation at the stand level, and often results in a lack of precision which affects forest management. Very high spatial resolution imagery, such as WorldView-3 and Light Detection and Ranging have the potential to overcome this issue. The main objective of this study is to map 11 tree species at the tree level using an object-based approach. For modeling, 240 variables were derived from WorldView-3 with pixel-based and arithmetic feature calculation techniques. A global approach (11 species) was compared to a hierarchical approach at two levels: (1) tree type (broadleaf/conifer) and (2) individual broadleaf (five) and conifer (six) species. Five different model techniques were compared: support vector machine, classification and regression tree, random forest (RF), k-nearest neighbors, and linear discriminant analysis. Each model was assessed using 16-band or first 8-band derived variables, with the results indicating higher precision for the RF technique. Higher accuracies were found using 16-band instead of 8-band derived variables for the global approach (overall accuracy (OA): 75% vs. 71%, Kappa index of agreement (KIA): 0.72 vs. 0.67) and tree type level (OA: 99% vs. 97%, KIA: 0.97 vs. 0.95). For broadleaf individual species, higher accuracy was found using first 8-band derived variables (OA: 70% vs. 68%, KIA: 0.63 vs. 0.60). No distinction was found for individual conifer species (OA: 94%, KIA: 0.93). This paper demonstrates that a hierarchical classification approach gives better results for conifer species and that using an 8-band WorldView-3 instead of a 16-band is sufficient. Full article

Percentile features derived from Landsat time-series data are widely adopted in land-cover classification. However, the temporal distribution of Landsat valid observations is highly uneven across different pixels due to the gaps resulting from clouds, cloud shadows, snow, and the scan line corrector (SLC)-off problem. In addition, when applying percentile features, land-cover change in time-series data is usually not considered. In this paper, an improved percentile called the time-series model (TSM)-adjusted percentile is proposed for land-cover classification based on Landsat data. The Landsat data were first modeled using three different time-series models, and the land-cover changes were continuously monitored using the continuous change detection (CCD) algorithm. The TSM-adjusted percentiles for stable pixels were then derived from the synthetic time-series data without gaps. Finally, the TSM-adjusted percentiles were used for generating supervised random forest classifications. The proposed methods were implemented on Landsat time-series data of three study areas. The classification results were compared with those obtained using the original percentiles derived from the original time-series data with gaps. The results show that the land-cover classifications obtained using the proposed TSM-adjusted percentiles have significantly higher overall accuracies than those obtained using the original percentiles. The proposed method was more effective for forest types with obvious phenological characteristics and with fewer valid observations. In addition, it was also robust to the training data sampling strategy. Overall, the methods proposed in this work can provide accurate characterization of land cover and improve the overall classification accuracy based on such metrics. The findings are promising for percentile-based land cover classification using Landsat time series data, especially in the areas with frequent cloud coverage. Full article

Based on high resolution satellite observations of sea surface temperature (SST), warm filaments near the Kuroshio around the Luzon Strait were systematically identified. These filaments extend an average length of about 200 km from the Kuroshio. The occurrence and features of the warm filaments are highly associated with both mesoscale eddies and the intensity of the SST gradient of the Kuroshio. Warm filaments are formed by heat advection from the warm Kuroshio into the colder interior Pacific Ocean by anticyclonic eddies (∼58%), cyclonic eddies (∼10%), and the dipole eddies (∼16%). The large temperature gradient near the Batanes Islands may also contribute to the high frequency of warm filaments in their vicinity. This study will help elucidate the role of zonal heat transport associated with the Kuroshio–eddy interaction during filament formation. Full article

Forest above-ground biomass (AGB) can be estimated based on light detection and ranging (LiDAR) point clouds. This paper introduces an accurate and detailed quantitative structure model (AdQSM), which can estimate the AGB of large tropical trees. AdQSM is based on the reconstruction of 3D tree models from terrestrial laser scanning (TLS) point clouds. It represents a tree as a set of closed and complete convex polyhedra. We use AdQSM to model 29 trees of various species (total 18 species) scanned by TLS from three study sites (the dense tropical forests of Peru, Indonesia, and Guyana). The destructively sampled tree geometry measurement data is used as reference values to evaluate the accuracy of diameter at breast height (DBH), tree height, tree volume, branch volume, and AGB estimated from AdQSM. After AdQSM reconstructs the structure and volume of each tree, AGB is derived by combining the wood density of the specific tree species from destructive sampling. The AGB estimation from AdQSM and the post-harvest reference measurement data show a satisfying agreement. The coefficient of variation of root mean square error (CV-RMSE) and the concordance correlation coefficient (CCC) are 20.37% and 0.97, respectively. AdQSM provides accurate tree volume estimation, regardless of the characteristics of the tree structure, without major systematic deviations. We compared the accuracy of AdQSM and TreeQSM in modeling the volume of 29 trees. The tree volume from AdQSM is compared with the reference value, and the determination coefficient (R²), relative bias (rBias), and CV-RMSE of tree volume are 0.96, 6.98%, and 22.62%, respectively. The tree volume from TreeQSM is compared with the reference value, and the R², relative Bias (rBias), and CV-RMSE of tree volume are 0.94, −9.69%, and 23.20%, respectively. The CCCs between the volume estimates based on AdQSM, TreeQSM, and the reference values are 0.97 and 0.96. AdQSM also models the branches in detail. The volume of branches from AdQSM is compared with the destructive measurement reference data. The R², rBias, and CV-RMSE of the branches volume are 0.97, 12.38%, and 36.86%, respectively. The DBH and height of the harvested trees were used as reference values to test the accuracy of AdQSM’s estimation of DBH and tree height. The R², rBias, and CV-RMSE of DBH are 0.94, −5.01%, and 9.06%, respectively. The R², rBias, and CV-RMSE of the tree height were 0.95, 1.88%, and 5.79%, respectively. This paper provides not only a new QSM method for estimating AGB based on TLS point clouds but also the potential for further development and testing of allometric equations. Full article

The spatial and temporal scale of rainfall datasets is crucial in modeling hydrological processes. Recently, open-access satellite precipitation products with improved resolution have evolved as a potential alternative to sparsely distributed ground-based observations, which sometimes fail to capture the spatial variability of rainfall. However, the reliability and accuracy of the satellite precipitation products in simulating streamflow need to be verified. In this context, the objective of the current study is to assess the performance of three rainfall datasets in the prediction of daily and monthly streamflow using Soil and Water Assessment Tool (SWAT). We used rainfall data from three different sources: Climate Hazards Group InfraRed Rainfall with Station data (CHIRPS), Climate Forecast System Reanalysis (CFSR) and observed rain gauge data. Daily and monthly rainfall measurements from CHIRPS and CFSR were validated using widely accepted statistical measures, namely, correlation coefficient (CC), root mean squared error (RMSE), probability of detection (POD), false alarm ratio (FAR), and critical success index (CSI). The results showed that CHIRPS was in better agreement with ground-based rainfall at daily and monthly scale, with high rainfall detection ability, in comparison with the CFSR product. Streamflow prediction across multiple watersheds was also evaluated using Kling-Gupta Efficiency (KGE), Nash-Sutcliffe Efficiency (NSE) and Percent BIAS (PBIAS). Irrespective of the climatic characteristics, the hydrologic simulations of CHIRPS showed better agreement with the observed at the monthly scale with the majority of the NSE values ranging between 0.40 and 0.78, and KGE values ranging between 0.62 and 0.82. Overall, CHIRPS outperformed the CFSR rainfall product in driving SWAT for streamflow simulations across the multiple watersheds selected for the study. The results from the current study demonstrate the potential of CHIRPS as an alternate open access rainfall input to the hydrologic model. Full article

Remote sensing (RS) soil moisture (SM) products have been widely used in various environmental studies. Understanding the error structure of data is necessary to properly apply RS SM products in trend and variation analysis and data fusion. However, a spatially continuous assessment of RS SM datasets is impeded by the limited spatial distribution of ground-based observations. As an alternative, the RS apparent thermal inertia (ATI) data related to the SM are transformed into SM values to expand the validation space. To obtain error components, the ATI-based SM along with the Soil Moisture Active Passive Mission (SMAP) and Advanced Microwave Scanning Radiometer 2 (AMSR2) SM are applied with the triple-collocation (TC) method to evaluate the RS SM data regarding random errors and amplitude variances at the regional scale. When the ATI-based SM is regarded as the reference data, the amplitude biases of the other two datasets are determined. The mean bias is also estimated by calculating the mean value difference between the ATI-based and validated RS SM. The results show that the ATI-based SM is a reliable source of reference data that, when combined with the TC method, can correctly estimate the error structure of RS SM datasets in wide space, promoting the reasonable application and calibration of RS SM datasets. Full article

Global rapid expansion of the coastal aquaculture industry has made great contributions to enhance food security, but has also caused a series of ecological and environmental issues. Sustainable management of coastal areas requires the explicit and efficient mapping of the spatial distribution of aquaculture ponds. In this study, a Google Earth Engine (GEE) application was developed for mapping coastal aquaculture ponds at a national scale with a novel classification scheme using Sentinel-1 time series data. Relevant indices used in the classification mainly include the water index, texture, and geometric metrics derived from radar backscatter, which were then used to segment and classify aquaculture ponds. Using this approach, we classified aquaculture ponds for the full extent of the coastal area in Vietnam with an overall accuracy of 90.16% (based on independent sample evaluation). The approach, enabling wall-to-wall mapping and area estimation, is essential to the efficient monitoring and management of aquaculture ponds. The classification results showed that aquaculture ponds are widely distributed in Vietnam’s coastal area and are concentrated in the Mekong River Delta and Red River delta (85.14% of the total area), which are facing the increasing collective risk of climate change (e.g., sea level rise and salinity intrusion). Further investigation of the classification results also provides significant insights into the stability and deliverability of the approach. The water index derived from annual median radar backscatter intensity was determined to be efficient at mapping water bodies, likely due to its strong response to water bodies regardless of weather. The geometric metrics considering the spatial variation of radar backscatter patterns were effective at distinguishing aquaculture ponds from other water bodies. The primary use of GEE in this approach makes it replicable and transferable by other users. Our approach lays a solid foundation for intelligent monitoring and management of coastal ecosystems. Full article

The detection of gas plumes from multibeam water column (MWC) data is the most direct way to discover gas hydrate reservoirs, but current methods often have low reliability, leading to inefficient detections. Therefore, this paper proposes an automatic method for gas plume detection and segmentation by analyzing the characteristics of gas plumes in MWC images. This method is based on the AdaBoost cascade classifier, combining the Haar-like feature and Local Binary Patterns (LBP) feature. After obtaining the detected result from the above algorithm, a target localization algorithm, based on a histogram similarity calculation, is given to exactly localize the detected target boxes, by considering the differences in gas plume and background noise in the backscatter strength. On this basis, a real-time segmentation method is put forward to get the size of the detected gas plumes, by integration of the image intersection and subtraction operation. Through the shallow-water and deep-water experiment verification, the detection accuracy of this method reaches 95.8%, the precision reaches 99.35% and the recall rate reaches 82.7%. Integrated with principles and experiments, the performance of the proposed method is analyzed and discussed, and finally some conclusions are drawn. Full article

Detection of road pavement cracks is important and needed at an early stage to repair the road and extend its lifetime for maintaining city roads. Cracks are hard to detect from images taken with visible spectrum cameras due to noise and ambiguity with background textures besides the lack of distinct features in cracks. Hyperspectral images are sensitive to surface material changes and their potential for road crack detection is explored here. The key observation is that road cracks reveal the interior material that is different from the worn surface material. A novel asphalt crack index is introduced here as an additional clue that is sensitive to the spectra in the range 450–550 nm. The crack index is computed and found to be strongly correlated with the appearance of fresh asphalt cracks. The new index is then used to differentiate cracks from road surfaces. Several experiments have been made, which confirmed that the proposed index is effective for crack detection. The recall-precision analysis showed an increase in the associated F1-score by an average of 21.37% compared to the VIS2 metric in the literature (a metric used to classify pavement condition from hyperspectral data). Full article

Stable and efficient ground moving target tracking and refocusing is a hard task in synthetic aperture radar (SAR) data processing. Since shadows in video-SAR indicate the actual positions of moving targets at different moments without any displacement, shadow-based methods provide a new approach for ground moving target processing. This paper constructs a novel framework to refocus ground moving targets by using shadows in video-SAR. To this end, an automatic-registered SAR video is first obtained using the video-SAR back-projection (v-BP) algorithm. The shadows of multiple moving targets are then tracked using a learning-based tracker, and the moving targets are ultimately refocused via a proposed moving target back-projection (m-BP) algorithm. With this framework, we can perform detecting, tracking, imaging for multiple moving targets integratedly, which significantly improves the ability of moving-target surveillance for SAR systems. Furthermore, a detailed explanation of the shadow of a moving target is presented herein. We find that the shadow of ground moving targets is affected by a target’s size, radar pitch angle, carrier frequency, synthetic aperture time, etc. With an elaborate system design, we can obtain a clear shadow of moving targets even in X or C band. By numerical experiments, we find that a deep network, such as SiamFc, can easily track shadows and precisely estimate the trajectories that meet the accuracy requirement of the trajectories for m-BP. Full article

Spectroscopy has demonstrated the ability to predict specific soil properties. Consequently, it is a promising avenue to complement the traditional methods that are costly and time-consuming. In the visible-near infrared (Vis-NIR) region, spectroscopy has been widely used for the rapid determination of organic components, especially soil organic carbon (SOC) using laboratory dry (lab-dry) measurement. However, steps such as collecting, grinding, sieving and soil drying at ambient (room) temperature and humidity for several days, which is a vital process, make the lab-dry preparation a bit slow compared to the field or laboratory wet (lab-wet) measurement. The use of soil spectra measured directly in the field or on a wet sample remains challenging due to uncontrolled soil moisture variations and other environmental conditions. However, for direct and timely prediction and mapping of soil properties, especially SOC, the field or lab-wet measurement could be an option in place of the lab-dry measurement. This study focuses on comparison of field and naturally acquired laboratory measurement of wet samples in Visible (VIS), Near-Infrared (NIR) and Vis-NIR range using several pretreatment approaches including orthogonal signal correction (OSC). The comparison was concluded with the development of validation models for SOC prediction based on partial least squares regression (PLSR) and support vector machine (SVMR). Nonetheless, for the OSC implementation, we use principal component regression (PCR) together with PLSR as SVMR is not appropriate under OSC. For SOC prediction, the field measurement was better in the VIS range with R²_CV = 0.47 and RMSEPcv = 0.24, while in Vis-NIR range the lab-wet measurement was better with R²_CV = 0.44 and RMSEPcv = 0.25, both using the SVMR algorithm. However, the prediction accuracy improves with the introduction of OSC on both samples. The highest prediction was obtained with the lab-wet dataset (using PLSR) in the NIR and Vis-NIR range with R²_CV = 0.54/0.55 and RMSEPcv = 0.24. This result indicates that the field and, in particular, lab-wet measurements, which are not commonly used, can also be useful for SOC prediction, just as the lab-dry method, with some adjustments. Full article

Spaceborne synthetic aperture radar (SAR) has been treated as a weather independent system for a long time. However, with the development of advanced SAR configurations, e.g., high resolution, bistatic, geosynchronous (GEO), the influence of tropospheric propagation error, which strongly depends on the weather, has begun to receive attention. In this paper, we focus on the effect of deterministic background tropospheric delay (BTD) during the image formation of GEO SAR. First, the decorrelation problems caused by the spatial variation and BTD are presented. Second, by combining with the SAR imaging geometry, the BTD error is decomposed as constant error, spatially variant error, and time variant error, the influences of which are analyzed under different circumstances. Third, an imaging method starting from the meteorological parameters and the GEO SAR systematic parameters is proposed to deal with the decorrelation problems. Finally, simulations with the dot-matrix targets are performed to validate the imaging method. Full article