Next Article in Journal
A Novel Roll and Pitch Estimation Approach for a Ground Vehicle Stability Improvement Using a Low Cost IMU
Previous Article in Journal
A Novel Piezoresistive MEMS Pressure Sensors Based on Temporary Bonding Technology
Previous Article in Special Issue
A Decade of Ground Deformation in Kunming (China) Revealed by Multi-Temporal Synthetic Aperture Radar Interferometry (InSAR) Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of the Trend of Deformation around the Kanto Region Estimated Using the Time Series of PALSAR-2 Data †

1
College of Industrial Technology, Nihon University, Chiba 2758575, Japan
2
Harris Geospatial Co., Tokyo 113-0033, Japan
*
Author to whom correspondence should be addressed.
This paper is an extension version of conference paper Nonaka, T.; Asaka, T.; Iwashita, K.; Ogushi, F. ANALYSIS OF THE TREND OF THE DEFORMATION AROUND KANTO REGION ESTIMATEDBY TIME SERIES OF PALSAR-2 DATA. In Proceedings of the International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences—GeoInformation for Disaster Management, Prague, Czech Republic, 3–6 September 2019.
Submission received: 8 October 2019 / Revised: 24 December 2019 / Accepted: 27 December 2019 / Published: 7 January 2020
(This article belongs to the Special Issue SAR and Optical Data for Crustal Deformation Monitoring)

Abstract

:
In the Kanto region of Japan, a large quantity of natural gas is dissolved in brine. The large-scale production of gas and iodine in the region has caused large-scale land subsidence in the past. Therefore, continuous and accurate monitoring for subsidence using satellite remote sensing is essential to prevent extreme subsidence and ensure the safety of residences. This study focused on the small baseline subset (SBAS) method to assess ground deformation trends around the Kanto region. Data for the SBAS method was acquired by the Advanced Land Observing Satellite (ALOS)-2 Phased Array type L-band Synthetic Aperture Radar (PALSAR)-2 from 2015 to 2019. A comparison of our results with reference levelling data shows that the SBAS method underestimates displacement. We corrected our results using linear regression and determined the maximum displacement around the Kujyukuri area to be approximately 20 mm/year; the mean displacement rate for 2015–2019 was −7.9 ± 2.9 mm/year. These values exceed those obtained using past PALSAR observations owing to the horizontal displacement after the Great East Japan Earthquake of 2011. Moreover, fewer points were acquired, and the root mean-squared error of each time-series displacement value was larger in our results. Further analysis is needed to address these bias errors.

1. Introduction

The south Kanto gas field contains a large quantity of natural gas in the water, which distributed widely across the Kanto region of Japan [1]. The natural gas produced in the South Kanto gas field accounts for more than 90% of the total Japanese production of natural gas dissolved in water. Since the 1930s, gas-production activities have been undertaken by private companies [2]. Large-scale production of gas and iodine began in the mid-1950s in Mobara-shi (the Kujyukuri area) of the Kanto region. This production involved the withdrawal of groundwater from depths of 500 to 2400 m [1]. Consequently, large-scale land subsidence began to occur in the 1960s, and serious land-subsidence surveys were started in the late 1960s [3]. From 1970 to 2010, the maximum accumulated subsidence was 1.0 m, and the subsided area spanned 800 km2 [3]. It was reported that the trend of land subsidence has reduced in recent years owing to the efforts of administrations and companies [4]. Land subsidence causes the collapse or tilting of buildings, discontinuities in pipelines, and intrusion of seawater into freshwater aquifers. That is, in addition to the destruction of residential and industrial facilities, it also causes environmental problems and geological disasters. Therefore, understanding land deformation trends caused by human action is significant for sustainable development.
One of the methods for monitoring large-scale ongoing deformation effectively at a low cost is remote sensing. In contrast, conventional levelling can provide precise measurements. The use of the differential interferometry synthetic aperture radar (DInSAR) is a well-known technique for measuring surface deformation [5,6]. However, spatial and temporal decorrelation severely restrict the use of 2-, 3-, and 4-pass DInSAR [7]. The small baseline subset (SBAS) technique is an extension of conventional DInSAR methods that reduces the effects of the decorrelation [7,8]. This technique is intended to be used for the analysis of distributed targets. Although the result of SBAS resembles that generated by conventional DInSAR, the key difference is that SBAS enables the analysis of large time series data. Another deformation monitoring method, especially in urban areas, is the persistent scatterers (PS) technique [9,10]. SBAS is less sensitive to the number of acquisitions than PS because it exploits spatially distributed coherence instead of exclusively considering single points. In addition, SBAS is not limited to linear displacement.
A radar sensor operates in the microwave domain. Different satellite SAR systems operate at different frequencies, which range from 3 MHz to 300 GHz. The longer the wavelength, the greater the ability to penetrate the materials. In previous studies, the SBAS technique has been applied to Phased Array type L-band Synthetic Aperture Radar (PALSAR) data to evaluate the surface deformation around the Chiba Prefecture [11,12]. For example, 26 PALSAR scenes were utilized to monitor deformation around the Kujyukuri area from 2006 to 2010 [11]. According to the analysis of the data, a displacement of more than 15 mm/year was observed around the cities of Yachimata and Oamishirasato. These displacement values were compared with GNSS Earth Observation Network System (GEONET) data at 12 validation points; the results showed that the root mean square error (RMSE) was approximately 10 mm.
PALSAR-2 on the Advanced Land Observing Satellite (ALOS)-2 has been operated since 2014 and is expected to be used for many applications, especially disaster management [13,14]. The orbit can be controlled with high accuracy, enabling improved coherence in PALSAR-2 in comparison with PALSAR [15]. In previous studies, displacement was evaluated using DInSAR, and large displacement errors occurred depending on the pairs of images [16]. However, knowledge regarding the applicability of the SBAS technique and the errors in the mean displacement velocity, as well as the difference between the results acquired from PALSAR-2 and PALSAR, were not gathered fully in previous studies. Therefore, the purpose of this study was to evaluate the displacement accuracy using PALSAR-2 data and reveal the displacement in the Kujyukuri area by this method. Additionally, we discussed the applicability of the SBAS technique to the PALSAR-2 data by comparing the acquired results with those of PALSAR.

2. Study Area and Used Data

2.1. Study Area

Figure 1 shows the study area (demarcated by a yellow box in the figure), which encompasses most of the South Kanto gas field. It includes a large part of the Chiba Prefecture and small parts of Tokyo, Kanagawa, Saitama, and Ibaraki. A large part of the prefecture is made up of the Boso Peninsula, which shields Tokyo Bay from the Pacific Ocean. The northern and western parts of the Boso Peninsula are highly urbanized, which has caused overpopulation in this area (the total population of the Chiba Prefecture was more than 6.2 million in October 2018 [17]). The elevation of this area ranges only from 0 m to approximately 300 m, and the topography is steep. Even today, most of the natural-gas development in Chiba is concentrated in the Kujyukuri area, which is on the east coast of the Boso peninsula.

2.2. Used Data

The ALOS-2, developed by Japan Aerospace Exploration Agency (JAXA), was launched on 24 May 2014 as the successor to ALOS-1 (launched on 24 January 2006 and decommissioned in 12 May 2011). ALOS-2 was designed for many applications, such as disaster monitoring, land monitoring, agricultural monitoring, natural resource exploitation, and global forest monitoring. ALOS-2 is equipped with PALSAR-2, which provides imaging with a high spatial resolution of 1–3 m via the spotlight and Stripmap mode. ALOS-2 has a 14 day revisit cycle, which makes it three times faster than ALOS (with a 46 day cycle). Furthermore, the high-accuracy orbit control, which is maintained within a 500 m radius tube for the reference orbit, enables high coherence for InSAR processing [16]. Table 1 lists the acquisition parameters used in this study. We collected 15 PALSAR-2 scenes from 15 January 2015 to 7 March 2019 (Table 2) for the InSAR stacking process. All the scenes were observed in the ultrafine single mode (with a 3 m spatial resolution), the descending direction, and right-looking mode.
The Geospatial Information Authority of Japan (GSI) has been operating a network of continuous Global Navigation Satellite Systems (GNSS) stations called GEONET [18]. We exploited GEONET and the levelling survey measurement data [19] to validate the InSAR derived deformation. Four points (93018: Shiroi, 93020: Tomisato, 93030: Ichihara, and 93033: Chosei) were selected for validation from the daily GPS data provided by the GSI (Figure 1). Levelling survey measurements were conducted in the Chiba Prefecture annually to understand the status of land deformation in the Chiba Prefecture. Figure 2 shows the mean velocity map derived from the levelling data for the study site from 2015 to 2018 and from 2006 to 2010. It can be seen that the spatial trend almost agrees between these durations, although the mean velocity decreased for the former period. Furthermore, subsidences of more than 10 mm/year appear to have occurred at several points around the Kujyukuri area [20].

3. Method

Figure 3 shows the flow of the study. We used the SARScape module of the ENVI 5.5 software, provided by Harris Geospatial Solutions, to perform the SBAS analysis. The SBAS workflow included the following seven steps [21,22]: the creation of a connection graph (generating all differential interferograms from the input image stack that fulfil the criteria for temporal and geometric baselines within a given interval in time or for a normal baseline with respect to the critical one), differential interferogram generation (spectral shift and adaptive filtering), phase unwrapping (either with a conventional 2D approach or a combined 3D approach), orbit refinement and reflattening, first estimation of the average displacement, atmospheric phase screen, and final estimation of the average displacement. In this study, we achieved a ground resolution of 8 m by using a multi-looking factor of three in the range and four in the azimuth. The topographic phase was removed using a 30 m digital elevation model (DEM) developed by the Shuttle Radar Topography Mission (SRTM) [23], which has been used for processing previous PALSAR data, as well. We applied the Goldstein filter [24] to smooth the differential phase and applied the minimum cost flow (MCF) network [25] along with the Delaunay method to unwrap the differential interferograms. In this stage, the unwrapping coherence threshold was 0.30. SBAS estimated displacement relative to the ground control points (GCPs) used during phase-displacement conversion. A previous study stated that the selection of the reference points is not real issue [7]; we selected approximately one hundred GCPs evenly across the entire image; these GCPs were also used for the PALSAR SBAS processing, refinement, and re-flattening. Subsequently, the average displacement rate and residual height-correction factors were estimated by inverting a linear system, which included all measures (one for each interferogram) together with their proportionality coefficients, depending on the temporal and geometric baseline of each pair [8]. Then, low-pass and high-pass spatial filters were used for the time-series images, to screen and remove the atmospheric phase components. All of the final displacement measurements were conducted in the satellite line of sight (LOS) direction and were geocoded in the WGS84 reference ellipsoid. To validate the displacement rates, we used the GEONET (Section 4.2.1) and levelling (Section 4.2.2) data as references. We established the approximate equation to correct the values based on the relationships between SBAS and levelling. Finally, the corrected displacement value was acquired for each point (Section 5).
Figure 4 shows the temporal and spatial baseline distributions of the SAR interferogram from the ALOS-2 PALSAR-2 dataset, in which each acquisition is represented by a diamond associated with an ID number; the green diamonds represent valid acquisitions, while the yellow diamonds represent the super master image (this image is the reference image for the entire process, and all the processed slant range pairs are co-registered on this reference geometry). Using the PALSAR-2 data, 45 interferograms were generated for the SBAS processing with respect to the multi-master images. The threshold criterion for obtaining the absolute maximum normal baselines was half that of the critical baseline (≒19.3 km), and that for the absolute maximum temporal baselines was 365 days. The image acquired on 9 March 2017, was chosen as the super master image. The maximum baseline among all the pairs was approximately 500 m, which satisfied the threshold criterion of the spatial baseline.

4. Results

4.1. Displacement Rates

Figure 5 shows the estimated mean displacement velocity from 2015 to 2019. Positive values indicate movement toward the satellite. Displacement was rarely observed around the coast of Tokyo Bay, while obvious displacement was observed around the Kujyukuri area. The mean displacement rate in the study area was −1.4 ± 3.2 mm/year. The results show that 75% of the points indicate displacement rates between −3 mm/year to 3 mm/year [26]. From 2006 to 2010, the displacement rates for this area derived from PALSAR data were −3.3 ± 5.8 mm/year in previous studies [11].
Some studies have obtained displacement values even for vegetated areas using longer L-band wavelengths [27,28]. However, using PALSAR-2, vegetated areas could not be observed due to lower coherence. Further, PALSAR-2 operates by acquiring data on fewer points than that acquired by PALSAR. A reason for this could be the larger RMSE of each time-series displacement value in the PALSAR-2 dataset (Section 4.2). Further analysis is required to understand the pertinent conditions, such as the number of images, threshold coherence, and the threshold for the spatial and temporal baseline.

4.2. Validation

4.2.1. Comparison with GEONET Data

We conducted a quantitative comparison of the time-series displacement values with the GEONET data for four validation points (Figure 1) to assess the accuracy of the InSAR-derived results. We did not use the values of several GEONET points inside the study area because negligible coherence was observed. The values derived from the GEONET data points were converted into values in the line of sight (LOS) direction for comparison (Figure 6). The overall displacement for these four stations indicates that the SBAS method (shown by green squares) underestimated the overall displacement. The displacement in Tomisato and Chosei was continuous, and the total displacement was approximately 50 mm from 2015 to 2019.
The estimated RMSE of the displacement value, which was 20 mm from 2015 to 2019, is twice that of the displacement from 2007 to 2010, and was derived from PALSAR data. Several factors caused errors in the displacement values: atmospheric disturbance [29], interferogram noise [30], and geometric errors in the sensor (specification: 5 m [31]). In the case of InSAR, the interferogram noise (α) is estimated by the ratio of DEM error (DEM ERR) to the height of ambiguity (dh/), which is defined by Equations (1) and (2):
D E M E R R = α d h d ϕ + β ,
d h d ϕ = λ R s i n θ 2 B p
where λ is the wavelength, R is the slant-range length, θ is the angle of incidence, and Bp is the perpendicular baseline. Previous studies have revealed that the interferometric phase noise of PALSAR-2 is 2.0%, half that of PALSAR [15,32,33]. The phase noise displacement error was estimated to be 3 mm by using a noise ratio of 2.0% in the case of PALSAR-2, which is much smaller than the RMSE displacement error of 20 mm. The geometric errors were estimated using a spatial deviation of 3 × 3 pixels around the validation points and were found to be less than 0.6 mm.

4.2.2. Comparison with Levelling Data

The mean velocity was quantitatively compared with the levelling data for 310 points. First, we took the leveling data for the same periods of study and converted the data into values in the LOS direction using the horizontal displacement from the nearest GEONET data (Figure 7). It should be noted that the reference levelling data were temporally too sparse (one data per year) to enable acquisition of the mean velocity. Figure 8 shows the comparison between SBAS-derived linear subsidence rates and the levelling data. The results show that the SBAS method underestimated the velocity for the whole range from −25 to 0 mm/year, and the RMSE is 6.8 mm/year. This trend can be expressed by the following linear regression, and the correlation coefficient is 0.82:
V L E V E L = 0.89 V S B A S 6.69
where VLEVEL is the mean velocity per level, and VSBAS is the velocity calculated by SBAS.

5. Discussion

In this article, we present a method for correcting the offset in SBAS-derived mean displacement velocity values using levelling data. According to the analysis presented in the previous section, we acquired the relationships (Equation (3)) between the SBAS and levelling data. Using Equation (3), we calculated the corrected mean displacement velocity values, as shown in Figure 9. The average and standard deviation of the corrected mean velocities were −7.9 mm/year and 2.9 mm/year, respectively. The spatial distribution of the mean velocity agreed with that of the levelling data (Figure 7). This implies that the SBAS-derived mean displacement velocity can supplement the levelling data subsequent to correction, even when the number of images is not sufficient. The time evolutions of SBAS were also corrected using the coefficient values derived in Equation (3). The amount of displacement can be obtained by calculating the difference in the mean velocity between the corrected and original SBAS values. The difference in the mean velocity can then be converted to offset the displacement by multiplying time, and the corrected SBAS displacement, DispCORRECTED SBAS, can be calculated, as described in Equation (4):
D i s p C O R R E C T E D   S B A S = D i s p S B A S + { ( 0.89 V S B A S 6.69 ) V S B A S } T
where DispSBAS is the original SBAS displacement, VSBAS is the mean velocity of the original SBAS displacement, and T is the length of time from the start of the observation. The corrected SBAS displacements are shown in red in Figure 6, and the time evolutions were better fitted to the GEONET data points than the original SBAS displacements in Figure 6a,b,d. However, the corrected SBAS displacement in Figure 6c was more overestimated than the GEONET displacement. According to the GEONET displacements in Figure 6, there are some types of displacements (linear and non-linear) in this study area, and these types are locally changed within a scene. This could be a potential reason why the estimated displacements by SBAS have some offsets with respect to the ground truth data.
Table 3 displays the differences between the validated mean displacement velocities of the values obtained from the corrected PALSAR-2, the original PALSAR-2, and the PALSAR data points. The RMSE of the PALSAR-2 data was twice that of the PALSAR data because of the bias in the displacement derived from the PALSAR-2 data. However, the high correlation coefficient for the linear regressions indicates that the deviation from the regression is relatively small. Therefore, after correcting the mean displacement velocity using linear regression, the RMSE of the PALSAR-2 data was reduced to 2.0 mm/year, almost half that of the PALSAR data. Based on this result, we assumed that the temporal displacement errors of PALSAR-2 were also reduced to half those of PALSAR. If this assumption is correct, our results are in agreement with the results of a previous study in which the phase noise of PALSAR-2 was found to be half that of PALSAR [15]. These results indicate the feasibility of using both PALSAR-2 and PALSAR data for displacement estimation using the SBAS technique. However, the accuracy of the obtained results is not good compared to those obtained from previous studies conducted on C-band satellites due to the following reasons: (1) The sensitivity of the L-band interferograms to surface deformation is much lower than that of the C-band [7]. (2) The observation frequency of PALSAR-2 is low (approximately four observations per year) compared with other C-band satellites, such as ERS-1, ERS-2, and Envisat [34,35].
Finally, it is worth determining whether the deformation patterns of PALSAR-2 and PALSAR are consistent with each other, although the acquisition period does not overlap. Therefore, the displacement values during the two time periods were compared. The acquisition conditions for the PALSAR data were very similar to those of the PALSAR-2 data: a 34.3° incidence angle and descending path. Figure 10 shows a comparison between the deformation rates in the LOS direction obtained from the PALSAR-2 and PALSAR data. A positive correlation can clearly be seen between them. Most of the subsidence values ranged from −5 to 5 mm/year from 2006 to 2010 and shifted to −10 to −5 mm/year from 2015 to 2019. The reason for the increase in the deformation rates despite a decrease in the subsidence measured by levelling during this period (Figure 2) is that the horizontal displacement increased after the extensive Great East Japan Earthquake in 2011.

6. Conclusions

To monitor the temporal changes in land subsidence caused by the extraction of natural gas, this study estimated the land displacement within the Chiba Prefecture from 2015 to 2019 using time-series satellite SAR data. We adopted the SBAS technique, and fifteen PALSAR-2 data points were used. The estimated displacement had a similar spatial trend to that derived from PALSAR. In contrast to previous studies using PALSAR data, the number of observed points in the present study was low. This can be explained by the observation that a substantial number of points of PALSAR-2 depicted a non-linear trend for each time series displacement value. This, in turn, was attributed to the insufficient number of images. This study quantitatively evaluated the accuracy of the annual displacement rate determined using PALSAR-2 data by comparing this data with in-situ levelling survey data. The estimation accuracy was found to be 6.8 mm/year, and PALSAR-2 underestimated the displacement values. After correction of the bias via linear regression, the RMSE of the displacement velocity derived using the SBAS method was reduced to 2.0 mm/year, approximately half that of PALSAR. The maximum velocity of displacement derived from the map subsequent to the bias correction was found to be approximately 20 mm/year in the Kujyukuri area for the period 2015–2019; however, this value is higher than that obtained for the 2006–2010 period. This can be explained by an increase in the horizontal displacement after the Great East Japan Earthquake that occurred in 2011. Through this study, we confirmed the applicability of the SBAS-derived displacement rates to PALSAR-2, as well as their reliability despite an insufficient number of PALSAR-2 images. In addition, this study quantitatively revealed the displacement errors using an extensive network of levelling data. Additional acquisition results in better product quality because it allows for better estimation and a reduction of the atmospheric component. However, to determine the causes of the bias errors obtained, further analyses are required. We intend to continue monitoring the temporal changes in displacement in the study area through SBAS analysis using PALSAR-2 data in future studies.

Author Contributions

Methodology, T.N.; software, T.N., T.A., and F.O.; formal analysis, T.N. and F.O.; writing—original draft preparation, T.N.; writing—review and editing, T.A. and K.I.; visualization, K.I.; project administration, T.N. All authors have read and agreed to the published version of the manuscript.

Funding

This study was a part of the project “Temporal analysis of ALOS and ALOS-2 data for monitoring the urban environment”, funded by the Japan Aerospace Exploration Agency (JAXA). The authors are grateful to JAXA for providing ALOS-2 PALSAR-2 data. This work was supported by JSPS KAKENHI grant number 17K06468.

Acknowledgments

The authors are grateful to JAXA for providing ALOS-2 PALSAR-2 data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Horiguchi, Y. Outline of water-dissolved natural gas field in Chiba Prefecture. J. Jpn. Assoc. Pet. Technol. 1998, 63, 475–484. (In Japanese) [Google Scholar] [CrossRef]
  2. Kanto Natural Gas Development Co., Ltd. History of Lodine. Available online: http://www.gasukai.co.jp/english/iodine/index2.html (accessed on 25 June 2019).
  3. Chiba Prefecture. Land Subsidence in Chiba Prefecture. Available online: https://www.pref.chiba.lg.jp/suiho/jibanchinka/torikumi/genkyou.html (accessed on 25 June 2019). (In Japanese).
  4. Furuno, K.; Kagawa, A.; Kazaoka, O.; Kusuda, T.; Nirei, K. Ground water management based on monitoring of land subsidence and groundwater levels in the Kanto groundwater basin, central Japan. Proc. IAHS 2015, 372, 53–57. [Google Scholar] [CrossRef] [Green Version]
  5. Gabriel, A.; Goldstein, R.; Zebker, H. Mapping small elevation changes over large areas: Differential radar interferometry. J. Geophys. Res. 1989, 94, 9183–9191. [Google Scholar] [CrossRef]
  6. Massonnet, D.; Rossi, M.; Carmona, C.; Adragna, F.; Peltzer, G.; Feigl, K.; Rabaute, T. The displacement field of the Landers earthquake mapped by radar interferometry. Nature 1993, 364, 138–142. [Google Scholar] [CrossRef]
  7. Ferretti, A. Satellite InSAR Data—Reservoir Monitoring from Space; EAGE Publications: Houten, The Netherlands, 2014; p. 159. ISBN 978-90-73834-71-2. [Google Scholar]
  8. Berardino, P.; Fornaro, G.; Lanari, R.; Sansosti, E. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms. IEEE Trans. Geosci. Remote Sens. 2002, 40, 23275–23283. [Google Scholar] [CrossRef] [Green Version]
  9. Ferretti, A.; Prati, C.; Rocca, F. Permanent scatterers in SAR interferometry. IEEE Trans. Geosci. Remote Sens. 2001, 39, 8–20. [Google Scholar] [CrossRef]
  10. Hooper, A.; Zebker, H.; Segall, P.; Kampes, B. A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers. Geophys. Res. Lett. 2004, 31, 1–5. [Google Scholar] [CrossRef]
  11. Ogushi, F. Monitoring of the Kanto plane in Japan with PALSAR-1 and PALSAR-2 with SBAS method. In Proceedings of the JAXA PI Meeting of Global Environment Observation Mission, Tokyo, Japan, 22–26 January 2018. [Google Scholar]
  12. Asaka, T.; Nonaka, T.; Hashiba, H.; Iwashita, K.; Sugimura, T. ALOS/PALSAR SBAS analysis of surface deformation in the Pacific side of Chiba prefecture, Japan. Trans. Jpn. Soc. Aeronaut. Space Sci. 2018, 16, 593–598. [Google Scholar] [CrossRef]
  13. Rosenqvist, A.; Shimada, M.; Suzuki, S.; Ogushi, F.; Tadono, T.; Watanabe, M.; Tsuzuku, K.; Watanabe, T.; Kamijyo, S.; Aoki, E. Operational performance of the ALOS global systematic acquisition strategy and observation plans for ALOS-2 PALSAR-2. Remote Sens. Environ. 2014, 155, 3–12. [Google Scholar] [CrossRef]
  14. Japan Aerospace Exploration Agency. ALOS-2 Project. Available online: https://www.eorc.jaxa.jp/ALOS/en/top/about_top.htm (accessed on 22 July 2019).
  15. Nonaka, T.; Asaka, T.; Iwashita, K.; Ogushi, F. Estimation of phase noise by time series analysis of ALOS-2 PALSAR-2 data. In Proceedings of the IGARSS2017 International Geoscience and Remote Sensing Symposium, Fort Worth, TX, USA, 23–28 July 2017. [Google Scholar]
  16. Asaka, T.; Nonaka, T.; Uchida, Y.; Sugimura, T. Evaluation of reinforced slope dynamics using ALOS-2/PALSAR-2. In Proceedings of the IGARSS2019 International Geoscience and Remote Sensing Symposium, Yokohama, Japan, 29 July–2 August 2019. [Google Scholar]
  17. Chiba Prefecture. Monthly Residents Survey Statics. Available online: https://www.pref.chiba.lg.jp/toukei/toukeidata/joujuu/geppou/2019/index.html (accessed on 22 July 2019). (In Japanese).
  18. Geospatial Information Authority of Japan, Service to Provide GEONET Data. Available online: https://terras.gsi.go.jp/ (accessed on 25 June 2019).
  19. Chiba Prefecture. Chiba Prefecture Levelling Survey Results. Available online: https://www.pref.chiba.lg.jp/suiho/jibanchinka/torikumi/seikaomote.html (accessed on 25 June 2019). (In Japanese).
  20. Tanno, Y.; Maruyama, M.; Kannno, R. Location of gas-chemical and metal machinery industries and their business change in Kujukuri area Chiba Prefecture. Ann. Hum. Reg. Geogr. 2006, 28, 101–126. (In Japanese) [Google Scholar]
  21. SARMap SBAS Tutorial. Available online: http://sarmap.ch/tutorials/sbas_tutorial_V_2_0.pdf (accessed on 25 June 2019).
  22. Pasquali, P.; Cantone, A.; Riccardi, P.; Defilippi, M.; Ogushi, F.; Gagliano, S.; Tamura, M. Mapping of Ground Deformations with Interferometric Stacking Techniques. In Land Application of Radar Remote Sensing; Holecz, F., Pasquali, P., Milisavljevic, N., Closson, D., Eds.; IntechOpen: London, UK, 2014; Available online: https://www.intechopen.com/books/land-applications-of-radar-remote-sensing/mapping-of-ground-deformations-with-interferometric-stacking-techniques (accessed on 25 June 2019). [CrossRef] [Green Version]
  23. Jarvis, A.; Reuter, H.; Nelson, A.; Guevara, E. Hole-filled SRTM for the Globe Version 4, Available from the CGIAR-CSI SRTM 90m Database. Available online: http://srtm.csi.cgiar.org (accessed on 25 June 2019).
  24. Goldstein, R.; Werner, C. Radar interferogram filtering for geophysical applications. Geophys. Res. Lett. 1998, 25, 4035–4038. [Google Scholar] [CrossRef] [Green Version]
  25. Costantini, M. A novel phase unwrapping method based on network programming. IEEE Trans. Geosci. Remote Sens. 1998, 36, 813–821. [Google Scholar] [CrossRef]
  26. Nonaka, T.; Asaka, T.; Iwashita, K.; Ogushi, F. Analysis of the trend of the deformation around Kanto region estimated by the time series of PALSAR-2 data. In Proceedings of the Geoinformation for Disaster Management 2019, Prague, Czech Republic, 3–6 September 2019. [Google Scholar]
  27. Aimaiti, Y.; Yamazaki, F.; Liu, W. Multi-sensor InSAR analysis of progressive land subsidence over coastal city of Urayasu, Japan. Remote Sens. 2018, 10, 1304. [Google Scholar] [CrossRef] [Green Version]
  28. Lu, Y.; Ke, C.; Zhou, X.; Wang, M.; Lin, H.; Chen, D.; Jiang, H. Monitoring land deformation in Changzhou city (China) with multi-band InSAR data sets from 2006 to 2012. Int. J. Remote Sens. 2018, 39, 1151–1174. [Google Scholar] [CrossRef]
  29. Hansen, R. Radar Interferometry: Data Interpretation and Error Analysis; Springer Netherlands: Dordrecht, The Netherlands, 2001; p. 308. ISBN 978-90-481-5696-2. [Google Scholar]
  30. Zebker, H.; Villasenor, J. Decorrelation in interferometric radar echoes. IEEE Trans. Geosci. Remote Sens. 1992, 30, 950–959. [Google Scholar] [CrossRef] [Green Version]
  31. Motooka, T.; Isoguchi, O.; Sakashita, M.; Shimada, M. Results of ALOS-2 PALSAR-2 calibration and validation after 3 years of operation. In Proceedings of the IGARSS2018 International Geoscience and Remote Sensing Symposium, Valencia, Spain, 22–27 July 2018. [Google Scholar]
  32. Nonaka, T.; Asaka, T.; Iwashita, K. Evaluation of atmospheric effects on interferograms using DEM errors of fixed ground points. Sensors 2018, 18, 2336. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Nonaka, T.; Asaka, T.; Iwashita, K. Feature of DEM produced by Interferometric technique by ALOS PALSAR data. In Proceedings of the 60th (2016 Spring) Meeting of the Remote Sensing Society of Japan, Chiba, Japan, 12–13 May 2016. (In Japanese). [Google Scholar]
  34. Casu, F.; Manzo, M.; Lanari, R. A quantitative assessment of SBAS algorithm performance for surface deformation retrieval from DInSAR data. Remote Sens. Environ. 2006, 102, 195–210. [Google Scholar] [CrossRef]
  35. Pepe, A.; Bonano, M.; Zhao, Q.; Yang, T.; Wang, H. The use of C-/X-band time-gapped SAR data and geotechnical models for the study of Shanghai’s ocean-reclaimed lands through the SBAS-DInSAR technique. Remote Sens. 2016, 8, 911. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of the study site (yellow box). The red stars indicate GeoNET stations.
Figure 1. Location of the study site (yellow box). The red stars indicate GeoNET stations.
Sensors 20 00339 g001
Figure 2. Mean subsidence velocity acquired from the levelling data for Chiba Prefecture from 2015 to 2019 (left) and from 2006 to 2010 (right).
Figure 2. Mean subsidence velocity acquired from the levelling data for Chiba Prefecture from 2015 to 2019 (left) and from 2006 to 2010 (right).
Sensors 20 00339 g002
Figure 3. Flowchart of the short baseline subset method used in this study. DEM, digital elevation model; GCP, ground control point.
Figure 3. Flowchart of the short baseline subset method used in this study. DEM, digital elevation model; GCP, ground control point.
Sensors 20 00339 g003
Figure 4. Temporal and spatial baseline distributions of the SAR interferogram, acquired from 15 ALOS-2 PALSAR-2 data points (the normal baseline restriction used was up to 50% of the critical baseline, and the time interval between the interferometric pair was up to 365 days). Each number shows the acquisition date (Table 2).
Figure 4. Temporal and spatial baseline distributions of the SAR interferogram, acquired from 15 ALOS-2 PALSAR-2 data points (the normal baseline restriction used was up to 50% of the critical baseline, and the time interval between the interferometric pair was up to 365 days). Each number shows the acquisition date (Table 2).
Sensors 20 00339 g004
Figure 5. Mean displacement velocity from 2015 to 2019, derived from the PALSAR-2 data [26].
Figure 5. Mean displacement velocity from 2015 to 2019, derived from the PALSAR-2 data [26].
Sensors 20 00339 g005
Figure 6. Comparison between the displacement calculated by the small baseline subset (SBAS) method and that obtained from the GEONET data for (a) Shiroi, (b) Tomisato, (c) Ichihara, and (d) Chosei. The green and red squares and black crosses represent the SBAS, corrected SBAS, and GEONET data, respectively.
Figure 6. Comparison between the displacement calculated by the small baseline subset (SBAS) method and that obtained from the GEONET data for (a) Shiroi, (b) Tomisato, (c) Ichihara, and (d) Chosei. The green and red squares and black crosses represent the SBAS, corrected SBAS, and GEONET data, respectively.
Sensors 20 00339 g006aSensors 20 00339 g006b
Figure 7. The reference mean velocities in the land observing satellite (LOS) direction calculated by the levelling data and GEONET data.
Figure 7. The reference mean velocities in the land observing satellite (LOS) direction calculated by the levelling data and GEONET data.
Sensors 20 00339 g007
Figure 8. Comparison between the mean velocities obtained using the SBAS method and the levelling data.
Figure 8. Comparison between the mean velocities obtained using the SBAS method and the levelling data.
Sensors 20 00339 g008
Figure 9. Corrected mean displacement velocity from 2015 to 2019, derived from the PALSAR-2 data.
Figure 9. Corrected mean displacement velocity from 2015 to 2019, derived from the PALSAR-2 data.
Sensors 20 00339 g009
Figure 10. Comparison between the mean displacement velocity from 2015 to 2019, derived from PALSAR-2 data, and from 2006 to 2010, derived from PALSAR data.
Figure 10. Comparison between the mean displacement velocity from 2015 to 2019, derived from PALSAR-2 data, and from 2006 to 2010, derived from PALSAR data.
Sensors 20 00339 g010
Table 1. Acquisition parameters of the Advanced Land Observing Satellite (ALOS-2) Phased Array Table 2. dataset.
Table 1. Acquisition parameters of the Advanced Land Observing Satellite (ALOS-2) Phased Array Table 2. dataset.
Band (wavelength)L-band (24 cm)
Operation modeUltrafine
PolarizationHH
Resolution3 m
Swathe50 km
Off-nadir angle32.8°
Orbit directionDescending
Revisit cycle14 days
Table 2. List of used PALSAR-2 data (path number: 18).
Table 2. List of used PALSAR-2 data (path number: 18).
NoAcquisition Date
115 January 2015
224 September 2015
310 March 2016
416 June 2016
511 August 2016
617 November 2016
79 March 2017
815 June 2017
924 August 2017
1016 November 2017
118 March 2018
1214 June 2018
1323 August 2018
1415 November 2018
157 March 2019
Table 3. Comparison between the validated mean displacement velocities of the values obtained Figure 2. original PALSAR-2, and PALSAR data.
Table 3. Comparison between the validated mean displacement velocities of the values obtained Figure 2. original PALSAR-2, and PALSAR data.
Corrected PALSAR-2Original PALSAR-2PALSAR
Time2015–20192006–2010
Number of points310519
RMSE of the mean velocity2.0 mm/year6.8 mm/year3.5 mm/year
Linear regression slope1.000.890.96
Linear regression intercept−0.03−6.69−1.44
Correlation coefficients0.820.70

Share and Cite

MDPI and ACS Style

Nonaka, T.; Asaka, T.; Iwashita, K.; Ogushi, F. Evaluation of the Trend of Deformation around the Kanto Region Estimated Using the Time Series of PALSAR-2 Data. Sensors 2020, 20, 339. https://0-doi-org.brum.beds.ac.uk/10.3390/s20020339

AMA Style

Nonaka T, Asaka T, Iwashita K, Ogushi F. Evaluation of the Trend of Deformation around the Kanto Region Estimated Using the Time Series of PALSAR-2 Data. Sensors. 2020; 20(2):339. https://0-doi-org.brum.beds.ac.uk/10.3390/s20020339

Chicago/Turabian Style

Nonaka, Takashi, Tomohito Asaka, Keishi Iwashita, and Fumitaka Ogushi. 2020. "Evaluation of the Trend of Deformation around the Kanto Region Estimated Using the Time Series of PALSAR-2 Data" Sensors 20, no. 2: 339. https://0-doi-org.brum.beds.ac.uk/10.3390/s20020339

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop