Search

For InSAR topographic mapping, multi-baseline InSAR height estimation is known to be an effective way to facilitate phase unwrapping by significantly increasing the ambiguity intervals and maintaining good height measurement sensitivity, especially in mountainous areas. In this paper, an efficient multi-baseline SAR interferometry approach based on maximum likelihood estimation is developed for refined topographic mapping in mountainous areas. In the algorithm, maximum likelihood (ML) height estimation is used to measure the topographic details and avoid the complicated phase unwrapping process. In order to be well-adapted to the mountainous terrain conditions, the prior height probability is re-defined to take the local terrain conditions and neighboring height constraint into consideration in the algorithm. In addition, three strategies are used to optimize the maximum likelihood height estimation process to obtain higher computational efficiency, so that this method is more suitable for spaceborne InSAR data. The strategies include substituting a rational function model into the complicated conversion process from candidate height to interferometric phase, discretizing the continuous height likelihood probability, and searching for the maximum likelihood height with a flexible step length. The experiment with simulated data is designed to verify the improvement of the ML height estimation accuracy with the re-defined prior height distribution. Then the optimized processing procedure is tested with the multi-baseline L-band ALOS/PALSAR data covering the Mount Tai area in China. The height accuracy of the generated multi-baseline InSAR DEM can meet both standards of American DTED-2 and Chinese national 1:50,000 DEM (mountain) Level 2. [ABSTRACT FROM AUTHOR]

Surface soil moisture (SSM) is a key variable for many environmental studies, including hydrology and agriculture. Synthetic aperture radar (SAR) data in the C-band are widely used nowadays to estimate SSM since the Sentinel-1 provides free-of-charge C-band SAR images at high spatial resolution with high revisit time, whereas the use of L-band is limited due to the low data availability. In this context, the main objective of this paper is to develop an operational approach for SSM estimation that mainly uses data in the C-band (Sentinel-1) with L-bands (ALOS/PALSAR) as additional data to improve SSM estimation accuracy. The approach is based on the use of the artificial neural networks (NNs) technique to firstly derive the soil roughness (Hrms) from the L-band (HH polarization) to then consider the L-band-derived Hrms and C-band SAR data (VV and VH polarizations) in the input vectors of NNs for SSM estimation. Thus, the Hrms estimated from the L-band at a given date is assumed to be constant for a given times series of C-band images. The NNs were trained and validated using synthetic and real databases. The results showed that the use of the L-band-derived Hrms in the input vector of NN in addition to C-band SAR data improved SSM estimation by decreasing the error (bias and RMSE), mainly for SSM values lower than 15 vol.% and regardless of Hrms values. Based on the synthetic database, the NNs that neglect the Hrms underestimate and overestimate the SSM (bias ranges between −8.0 and 4.0 vol.%) for Hrms values lower and higher than 1.5 cm, respectively. For Hrms <1.5 cm and most SSM values higher than 10 vol.%, the use of Hrms as an input in the NNs decreases the underestimation of the SSM (bias ranges from −4.5 to 0 vol.%) and provides a more accurate estimation of the SSM with a decrease in the RMSE by approximately 2 vol.%. Moreover, for Hrms values between 1.5 and 2.0 cm, the overestimation of SSM slightly decreases (bias decreased by around 1.0 vol.%) without a significant improvement of the RMSE. In addition, for Hrms >2.0 cm and SSM between 8 to 22 vol.%, the accuracy on the SSM estimation improved and the overestimation decreased by 2.2 vol.% (from 4.5 to 2.3 vol.%). From the real database, the use of Hrms estimated from the L-band brought a significant improvement of the SSM estimation accuracy. For in situ SSM less than 15 vol.%, the RMSE decreased by 1.5 and 2.2 vol.% and the bias by 1.2 and 2.6 vol.%, for Hrms values lower and higher than 1.5 cm, respectively. [ABSTRACT FROM AUTHOR]

Chiayi County is located in the largest alluvial plain of Taiwan with extensive aquaculture and rice farming sustained by water extracted from groundwater wells. Chiayi is a typical aquaculture area affected by land subsidence, yet such lands worldwide combine to provide nearly 90% of global aquaculture products, greatly reducing oceanic overfishing problems. This study uses precision leveling, multi-layer compaction monitoring well (MLCW) and spaceborne SAR interferometry (InSAR) to examine the cause and effect of land subsidence in Chiayi associated with groundwater extractions and changes. Heights at benchmarks in a leveling network are measured annually and soil compactions at 24-26 layers up to 300-m depths at 7 MLCWs are collected at one-month intervals. Over 2007-2011, 15 ALOS/PALSAR images are processed by the method of TCPInSAR to produce subsidence rates. All sensors show that land subsidence occur in most parts of Chiayi, with rates reaching 4.5 cm/year around its coast, a result of groundwater pumping from shallow to deep aquifers. MLCWs detect mm-accuracy seasonal soil compactions coinciding with groundwater level fluctuations and causing dynamic compactions. Compactions near Taiwan High Speed Rail may reduce the strength of the rail's supporting columns to degrade its safety. The SAR images yield subsidence rates consistent with those from leveling and compaction wells after corrections for systematic errors by the leveling result. Subsidence in Chiayi's coastal area leads to relative sea level rises at rates up to 15 times larger than the global eustatic sea level rising rate, a risk typical for world's aquaculture-rich regions. At the fish pond-covered Budai Township, InSAR identifies subsidence spots not detected by leveling, providing crucial geo-information for a sustainable land management for aquaculture industry. [ABSTRACT FROM AUTHOR]

Spaceborne Synthetic Aperture Radar (SAR) interferometry has been utilised to acquire high-resolution Digital Elevation Models (DEMs) with wide coverage, particularly for persistently cloud-covered regions where stereophotogrammetry is hard to apply. Since the discovery of sand buried drainage systems by the Shuttle Imaging Radar-A (SIR-A) L-band mission in 1982, radar images have been exploited to map subsurface features beneath a sandy cover of extremely low loss and low bulk humidity in some hyper-arid regions such as from the Japanese Earth Resources Satellite 1 (JERS-1) and Advanced Land Observing Satellite/Phased Array type L-band Synthetic Aperture Radar (ALOS/PALSAR). Therefore, we hypothesise that a Digital Elevation Model (DEM) derived by InSAR in hyper-arid regions is likely to represent a subsurface elevation model, especially for lower frequency radar systems, such as the L-band system (1.25 GHz). In this paper, we compare the surface appearance of radar images (L-band and C-band) with that of optical images to demonstrate their different abilities to show subsurface features. Moreover, we present an application of L-band InSAR to measure penetration depths in the eastern Sahara Desert. We demonstrate how the retrieved L-band InSAR DEM appears to be of a consistently 1-2 m lower elevation than the C-band Shuttle Radar Topography Mission (SRTM) DEM over sandy covered areas, which indicates the occurrence of penetration and confirms previous studies. [ABSTRACT FROM AUTHOR]

Remote sensing data proved to be a valuable resource in a variety of earth science applications. Using high-dimensional data with advanced methods such as machine learning algorithms (MLAs), a sub-domain of artificial intelligence, enhances lithological mapping by spectral classification. Support vector machines (SVM) are one of the most popular MLAs with the ability to define non-linear decision boundaries in high-dimensional feature space by solving a quadratic optimization problem. This paper describes a supervised classification method considering SVM for lithological mapping in the region of Souk Arbaa Sahel belonging to the Sidi Ifni inlier, located in southern Morocco (Western Anti-Atlas). The aims of this study were (1) to refine the existing lithological map of this region, and (2) to evaluate and study the performance of the SVM approach by using combined spectral features of Landsat 8 OLI with digital elevation model (DEM) geomorphometric attributes of ALOS/PALSAR data. We performed an SVM classification method to allow the joint use of geomorphometric features and multispectral data of Landsat 8 OLI. The results indicated an overall classification accuracy of 85%. From the results obtained, we can conclude that the classification approach produced an image containing lithological units which easily identified formations such as silt, alluvium, limestone, dolomite, conglomerate, sandstone, rhyolite, andesite, granodiorite, quartzite, lutite, and ignimbrite, coinciding with those already existing on the published geological map. This result confirms the ability of SVM as a supervised learning algorithm for lithological mapping purposes. [ABSTRACT FROM AUTHOR]

Forest environment classification in mountain regions based on single-sensor remote sensing approaches is hindered by forest complexity and topographic effects. Temperate broadleaf forests in western Asia such as the Hyrcanian forest in northern Iran have already suffered from intense anthropogenic activities. In those regions, forests mainly extend in rough terrain and comprise different stand structures, which are difficult to discriminate. This paper explores the joint analysis of Landsat7/ETM+, L-band SAR and their derived parameters and the effect of terrain corrections to overcome the challenges of discriminating forest stand age classes in mountain regions. We also verified the performances of three machine learning methods which have recently shown promising results using multisource data; support vector machines (SVM), neural networks (NN), random forest (RF) and one traditional classifier (i.e., maximum likelihood classification (MLC)) as a benchmark. The non-topographically corrected ETM+ data failed to differentiate among different forest stand age classes (average classification accuracy (OA) = 65%). This confirms the need to reduce relief effects prior data classification in mountain regions. SAR backscattering alone cannot properly differentiate among different forest stand age classes (OA = 62%). However, textures and PolSAR features are very efficient for the separation of forest classes (OA = 82%). The highest classification accuracy was achieved by the joint usage of SAR and ETM+ (OA = 86%). However, this shows a slight improvement compared to the ETM+ classification (OA = 84%). The machine learning classifiers proved t o be more robust and accurate compared to MLC. SVM and RF statistically produced better classification results than NN in the exploitation of the considered multi-source data. [ABSTRACT FROM AUTHOR]

This study focuses on evaluating the potential of ALOS/PALSAR time-series data to analyze the activation of deep-seated landslides in the foothill zone of the high mountain Alai range in the southern Tien Shan (Kyrgyzstan). Most previous field-based landslide investigations have revealed that many landslides have indicators for ongoing slow movements in the form of migrating and newly developing cracks. L-band ALOS/PALSAR data for the period between 2007 and 2010 are available for the 484 km2 area in this study. We analyzed these data using the Small Baseline Subset (SBAS) time-series technique to assess the surface deformation related to the activation of landslides. We observed up to ±17 mm/year of LOS velocity deformation rates, which were projected along the local steepest slope and resulted in velocity rates of up to -63 mm/year. The obtained rates indicate very slow movement of the deep-seated landslides during the observation time. We also compared these movements with precipitation and earthquake records. The results suggest that the deformation peaks correlate with rainfall in the 3 preceding months and with an earthquake event. Overall, the results of this study indicated the great potential of L-band InSAR time series analysis for efficient spatiotemporal identification and monitoring of slope activations in this region of high landslide activity in Southern Kyrgyzstan. [ABSTRACT FROM AUTHOR]

The Global Digital Elevation Model produced from stereo images of Advanced Spaceborne Thermal Emission and Reflection Radiometer data (ASTER GDEM) covers land surfaces between latitudes of 83°N and 83°S. The Phased Array type L-band Synthetic Aperture Radar (PALSAR) onboard Advanced Land Observing Satellite (ALOS) collected many SAR images since it was launched on 24 January 2006. The combination of ALOS/PALSAR interferometric data and ASTER GDEM should provide the penetration depth of SAR data assuming ASTER GDEM was the elevation of vegetation canopy top. It would be correlated with forest biomass because penetration depth could be affected by forest density and forest canopy height. Their combination held great promises for the forest biomass mapping over large area. The feasibility of forest biomass mapping through the data synthesis of ALOS/PALSAR InSAR data and ASTER GDEM was investigated in this study. A procedure for the extraction of penetration depth was firstly proposed. Then three models were built for biomass estimation: (I) model only using backscattering coefficients of ALOS/PALSAR data; (II) model only using penetration depth; (III) model using both of them. The biomass estimated from Lidar data was taken as reference data to evaluate the three different models. The results showed that the combination of backscattering coefficients and penetration depth gave the best accuracy. The forest disturbance has to be considered in forest biomass estimation because of the long time span of ASTER data for generating ASTER GDEM. The spatial homogeneity could be used to improve estimation accuracy. [ABSTRACT FROM AUTHOR]

A decomposition scheme was applied to ALOS/PALSAR data obtained from a fast-growing tree plantation in Sumatra, Indonesia to extract tree stem information and then estimate the forest stand volume. The scattering power decomposition of the polarimetric SAR data was performed both with and without a rotation matrix and compared to the following field-measured forest biometric parameters: tree diameter, tree height and stand volume. The analytical results involving the rotation matrix correlated better than those without the rotation matrix even for natural scattering surfaces within the forests. Our primary finding was that all of the decomposition powers from the rotated matrix correlated significantly to the forest biometric parameters when divided by the total power. The surface scattering ratio of the total power markedly decreased with the forest growth, whereas the canopy and double-bounce scattering ratios increased. The observations of the decomposition powers were consistent with the tree growth characteristics. Consequently, we found a significant logarithmic relationship between the decomposition powers and the forest biometric parameters that can potentially be used to estimate the forest stand volume. [ABSTRACT FROM AUTHOR]

At present, the Arabian Peninsula is one of the driest regions on Earth; however, this area experienced heavy rainfall in the past thousand years. During this period, catchments received substantial amounts of surface water and sustained vast networks of streams and paleolakes, which are currently inactive. The Advanced Land Observing Satellite (ALOS) Phased Array Type L-band Synthetic Aperture Radar (PALSAR) data reveal paleohydrologic features buried under shallow aeolian deposits in many areas of the ad-Dawasir, Sahba, Rimah/Batin, and as-Sirhan wadis. Optical remote-sensing data support that the middle of the trans-peninsula Wadi Rimah/Batin, which extends for ~1200 km from the Arabian Shield to Kuwait and covers ~200,000 km2, is dammed by linear sand dunes formed by changes in climate conditions. Integrating Landsat 8 Operational Land Imager (OLI), Geo-Eye, Shuttle Radar Topography Mission (SRTM) digital elevation model, and ALOS/PALSAR data allowed for the characterization of paleodrainage reversals and diversions shaped by structural and volcanic activity. Evidence of streams abruptly shifting from one catchment to another is preserved in Wadi ad-Dawasir along the fault trace. Volcanic activity in the past few thousand years in northern Saudi Arabia has also changed the slope of the land and reversed drainage systems. Relics of earlier drainage directions are well maintained as paleoslopes and wide upstream patterns. This study found that paleohydrologic activity in Saudi Arabia is impacted by changes in climate and by structural and volcanic activity, resulting in changes to stream direction and activity. Overall, the integration of radar and optical remote-sensing data is significant for deciphering past hydrologic activity and for predicting potential water resource areas. [ABSTRACT FROM AUTHOR]

In this study, we compared the accuracies of above-ground biomass (AGB) estimated by integrating ALOS (Advanced Land Observing Satellite) PALSAR (Phased-Array-Type L-Band Synthetic Aperture Radar) data and TanDEM-X-derived forest heights (TDX heights) at four scales from 1/4 to 25 ha in a hemi-boreal forest in Japan. The TDX heights developed in this study included nine canopy height models (CHMs) and three model-based forest heights (ModelHs); the nine CHMs were derived from the three digital surface models (DSMs) of (I) TDX 12 m DEM (digital elevation model) product, (II) TDX 90 m DEM product and (III) TDX 5 m DSM, which we developed from two TDX–TSX (TerraSAR-X) image pairs for reference, and the three digital terrain models (DTMs) of (i) an airborne Light Detection and Ranging (LiDAR)-based DTM (LiDAR DTM), (ii) a topography-based DTM and (iii) the Shuttle Radar Topography Mission (SRTM) DEM; the three ModelHs were developed from the two TDX-TSX image pairs used in (III) and the three DTMs (i to iii) with the Sinc inversion model. In total, 12 AGB estimation models were developed for comparison. In this study, we included the C-band SRTM DEM as one of the DTMs. According to Walker et al. (2007), the SRTM DEM serves as a DTM for most of the Earth's surface, except for the areas with extensive tree and/or shrub coverage, e.g., the boreal and Amazon regions. As our test site is located in a hemi-boreal zone with medium forest cover, we tested the ability of the SRTM DEM to serve as a DTM in our test site. This study especially aimed to analyze the capability of the two TDX DEM products (I and II) to estimate AGB in practice in the hemi-boreal region, and to examine how the different forest height creation methods (the simple DSM and DTM subtraction for the nine CHMs and the Sinc inversion model-based approach for the three ModelHs) and the different spatial resolutions of the three DSMs and three DTMs affected the AGB estimation results. We also conducted the slope-class analysis to see how the varying slopes influenced the AGB estimation accuracies. The results show that the combined use of the PALSAR data and the CHM derived from (I) TDX 12 m DEM and (i) LiDAR DTM achieved the highest AGB estimation accuracies across the scales (R2 ranged from 0.82 to 0.97), but the CHMs derived from (I) TDX 12 m DEM and another two DTMs, (ii) and (iii), showed low R2 values at any scales. In contrast, the two CHMs derived from (II) TDX 90 m DEM and both (i) LiDAR DTM and (iii) SRTM DEM showed high R2 values > 0.87 and 0.78, respectively, at the scales > 9.0 ha, but they yielded much lower R2 values at smaller scales. The three ModelHs gave the lowest R2 values across the scales (R2 ranged from 0.39 to 0.60). Analyzed by slope class at the 1.0 ha scale, however, all the 12 AGB estimation models yielded high R2 values > 0.66 at the lowest slope class (0° to 9.9°), including the three ModelHs (R2 ranged between 0.68 to 0.69). The two CHMs derived from (II) TDX 90 m DEM and both (i) LiDAR DTM and (iii) SRTM DEM showed R2 values of 0.80 and 0.71, respectively, at the lowest slope class, while the CHM derived from (I) TDX 12 m DEM and (i) LiDAR DTM showed high R2 values across the slope classes (R2 > 0.82). The results show that (I) TDX 12 m DEM had a high capability to estimate AGB, with a high accuracy across the scales and the slope classes in the form of CHM, but the use of (i) LiDAR DTM was required. On the other hand, (II) TDX 90 m DEM was able to achieve high AGB estimation accuracies not only with (i) LiDAR DTM, but also with (iii) SRTM DEM in the form of CHM, but it was limited to large scales > 9.0 ha; however, all the models developed in this study have the possibility to achieve higher AGB estimation accuracies at the 1.0 ha scale in flat terrains with slope < 10°. The analysis showed the strengths and limitations of each model, and it also indicates that the data creation methods, the spatial resolutions of datasets and topographic features affects the effective spatial scales for AGB mapping, and the optimal combinations of these features should be chosen to obtain high AGB estimation accuracies. [ABSTRACT FROM AUTHOR]

Robust remote monitoring of land cover changes is essential for a range of studies such as climate modeling, ecosystems, and environmental protection. However, since each satellite data has its own effective features, it is difficult to obtain high accuracy land cover products derived from a single satellite's data, perhaps because of cloud cover, suboptimal acquisition schedules, and the restriction of data accessibility. In this study, we integrated Landsat 5, 7, and 8, Sentinel-2, Advanced Land Observing Satellite Advanced Visual, and Near Infrared Radiometer type 2 (ALOS/AVNIR-2), ALOS Phased Array L-band Synthetic Aperture Radar (PALSAR) Mosaic, ALOS-2/PALSAR-2 Mosaic, Shuttle Radar Topography Mission (SRTM), and ancillary data, using kernel density estimation to map and analyze land use/cover change (LUCC) over Central Vietnam from 2007 to 2017. The region was classified into nine categories, i.e., water, urban, rice paddy, upland crops, grassland, orchard, forest, mangrove, and bare land by an automatic model which was trained and tested by 98,000 reference data collected from field surveys and visual interpretations. Results were the 2007 and 2017 classified maps with the same spatial resolutions of 10 m and the overall accuracies of 90.5% and 90.6%, respectively. They indicated that Central Vietnam experienced an extensive change in land cover (33 ± 18% of the total area) during the study period. Gross gains in forests (2680 km2) and water bodies (570 km2) were primarily from conversion of orchards, paddy fields, and crops. Total losses in bare land (495 km2) and paddy (485 km2) were largely to due transformation to croplands and urban & other infrastructure lands. In addition, the results demonstrated that using global land cover products for specific applications is impaired because of uncertainties and inconsistencies. These findings are essential for the development of resource management strategy and environmental studies. [ABSTRACT FROM AUTHOR]