The effects of ground deformation pose a significant geo-hazard to the environment and infrastructure in Wuhan, the most populous city in Central China, in the eastern Jianghan Plain at the intersection of the Yangtze and Han rivers. Prior to this study, however, rates and patterns of region-wide ground deformation in Wuhan were little known. Here we employ multi-temporal SAR interferometry to detect and characterize spatiotemporal variations of ground deformation in major metropolitan areas in Wuhan. A total of twelve TerraSAR-X images acquired during 2009–2010 are used in the InSAR time series analysis. InSAR-derived results are validated by levelling survey measurements and reveal a distinct subsidence pattern within six zones in major commercial and industrial areas, with a maximum subsidence rate up to −67.3 mm/year. A comparison analysis between subsiding patterns and urban developments as well as geological conditions suggests that land subsidence in Wuhan is mainly attributed to anthropogenic activities, natural compaction of soft soil, and karst dissolution of subsurface carbonate rocks. However, anthropogenic activities related to intensive municipal construction and industrial production have more significant impacts on the measured subsidence than natural factors. Moreover, remarkable signals of secular land uplift are found along both banks of the Yangtze River, especially along the southern bank, with deformation rates ranging mostly from +5 mm/year to +17.5 mm/year. A strong temporal correlation is highlighted between the detected displacement evolutions and the water level records of the Yangtze River, inferring that this previously unknown deformation phenomenon is likely related to seasonal fluctuations in water levels of the Yangtze River.
A hybrid method is proposed to calculate complete synthetic seismograms based on a spherically symmetric and self-gravitating Earth with a multilayered structure of atmosphere, ocean, mantle, liquid core and solid core. For large wavelengths, a numerical scheme is used to solve the geodynamic boundary-value problem without any approximation on the deformation and gravity coupling. With decreasing wavelength, the gravity effect on the deformation becomes negligible and the analytical propagator scheme can be used. Many useful approaches are used to overcome the numerical problems that may arise in both analytical and numerical schemes. Some of these approaches have been established in the seismological community and the others are developed for the first time. Based on the stable and efficient hybrid algorithm, an all-in-one code QSSP is implemented to cover the complete spectrum of seismological interests. The performance of the code is demonstrated by various tests including the curvature effect on teleseismic body and surface waves, the appearance of multiple reflected, teleseismic core phases, the gravity effect on long period surface waves and free oscillations, the simulation of near-field displacement seismograms with the static offset, the coupling of tsunami and infrasound waves, and free oscillations of the solid Earth, the atmosphere and the ocean. QSSP is open source software that can be used as a stand-alone FORTRAN code or may be applied in combination with a Python toolbox to calculate and handle Green's function databases for efficient coding of source inversion problems.
ABSTRACT Measurements of short-interval variations in glacier surface velocity, which contribute to our understanding of ice motion mechanisms, remain scarce on the Tibetan Plateau. Here we present sub-hourly measurements of glacier surface motion variations at the terminus region of Laohugou No. 12 Glacier. Field observations were collected over 4 d in July 2015 from terrestrial radar interferometry. The observed glacier displacement time series are generally in agreement with the results measured by differential GPS and highlight that glacier surface velocity is characterized by clear diurnal fluctuations in the study period. During day-time hours, glacier flow speeds were higher than 3.0 mm h −1 , whereas they were below 1.0 mm h −1 during night-time hours. The large diurnal fluctuations of glacier surface velocity indicate that variations in basal slip are the dominant motion mechanism. Moreover, a positive correlation ( R = 0.82, P < 0.001) between air temperature and glacier surface velocity suggests that glacier motion variations are probably affected by changes in air temperature during the ablation season.
Radar penetration correction is essential for accurately estimating glacier mass balance when using the geodetic methods based on the radar-derived digital elevation model (DEM). Due to heterogeneous snow distribution and snowpack properties, the radar penetration depth varies by region and is basically dependent on the altitudes. However, a constant value is usually used to correct the radar penetration for regional glacier mass balance estimation in high-mountain Asia (HMA), because little is known about the detailed distribution of radar penetration depth of the glacierized area in HMA. In this study, the penetration depth difference (hereafter referred to as PDD) of the C/X-band radar signals on glaciers over the whole HMA region was first estimated by comparing the Space Shuttle Radar Topographic Mission (SRTM) C/X-band digital elevation models (DEMs), and the spatial characteristics of the altitude-dependent PDD were analyzed on a 1° × 1° grid. The C/X-band radar PDD (the X minus C band penetration) was calculated in each 50-m elevation bin, and a linear model was used to study the correlation between the altitude and the penetration difference. The results show that the difference of the SRTM C/X-band penetration depth in the HMA region is between -1.67 ± 0.65 m and 7.63 ± 0.99 m, with an average value of 2.40±0.04 m. A significant linear positive correlation between elevation and the PDD was found in most of 1° × 1° grid regions in HMA. These results can be used as important reference data for estimating glacier mass balance in HMA based on SRTM DEMs.
Advanced techniques of multi-temporal InSAR (MT-InSAR) represent a valuable tool in ground subsidence studies allowing remote investigation of the behavior of mass movements in long time intervals by using large datasets of SAR images covering the same area and acquired at different epochs. Beijing is susceptible to subsidence, producing undesirable environmental impacts and affecting dense population. Excessive groundwater withdrawal is thought to be the primary cause of land subsidence, and rapid urbanization and economic development, mass construction of skyscrapers, highways and underground engineering facilities (e.g., subway) are also contributing factors. In this paper, a spatial–temporal analysis of the land subsidence in Beijing was performed using one of the MT-InSAR techniques, referred to as Small Baseline Subset (SBAS). This technique allows monitoring the temporal evolution of a deformation phenomenon, via the generation of mean deformation velocity maps and displacement time series from a data set of acquired SAR images. 52 C-band ENVISAT ASAR images acquired from June 2003 to August 2010 were used to produce a linear deformation rate map and to derive time series of ground deformation. The results show that there are three large subsidence funnels within this study area, which separately located in Balizhuang-Dajiaoting in Chaoyang district, Wangjing-Laiguangying Chaoyang district, Gaoliying Shunyi district. The maximum settlement center is Wangsiying-Tongzhou along the Beijing express; the subsidence velocity exceeds 110 mm/y in the LOS direction. In particular, we compared the achieved results with leveling measurements that are assumed as reference. The estimated long-term subsidence results obtained by SBAS approach agree well with the development of the over-exploitation of ground water, indicating that SBAS techniques is adequate for the retrieval of land subsidence in Beijing from multi-temporal SAR data.
The high spatial-temporal variability of soil moisture necessitates monitoring at a high resolution in order to improve our understanding of Earth system processes. Current large-scale soil moistures inferred from the microwave satellites have limited spatial resolution, typically in the range of tens of kilometers. Recent studies have revealed that synthetic aperture radar (SAR) backscatter exhibits qualitative relationships with soil moisture, suggesting the potential for large-scale high-resolution mapping of soil moisture. Here, we proposed a method for directly estimating soil moisture content based on the Advanced Integral Equation Model (AIEM) and Mironov dielectric model. The approach involves establishing a series of semi-empirical models, independent of preceding surface roughness determination, using two Envisat ASAR alternating polarization (AP) model precision products. We generate a time series of high-resolution soil moisture using Envisat ASAR AP data acquired from 2004 to 2011, with an uncertainty of approximately 0.05 m 3 /m 3 . Our soil moisture retrievals demonstrate very good agreement with European Space Agency (ESA) Climate Change Initiative (CCI) soil moisture products and the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 re-analysis hourly products, even in the absence of synchronous ground measurements. Furthermore, our study reveals good temporal coherence between drought and heavy rainfall events, and SAR-derived soil moisture, which suggests a potential to capture heavy rainfall and drought events. We conclude that SAR-derived soil moisture is a more direct and efficient method in quantifying soil moisture at a high spatial resolution, making it more suitable for watershed scale hydrological studies.
SUMMARY Analysis of data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission allows us to identify regions of long-term mass changes such as the areas of Glacial Isostatic Adjustment (GIA) in North America and Fennoscandia. As there are now more than 7 yr of data available, the determined trends are robust enough for the inference of viscosity structure of the Earth’s mantle. In this study, we focus on the Fennoscandian rebound area as there are abundant high-quality terrestrial data to use as ground-truth. In the first step, GRACE data are taken to fix the optimal radial (1-D) viscosity profile and the lithospheric thickness combination, which are needed as background parameters in 3-D earth modelling. The results are in basic agreement with results based upon relative sea level and GPS data, showing a lithospheric thickness in Fennoscandia between 90 and 160 km and an upper mantle viscosity of about [2–4] × 10 20 Pa s. The lower mantle viscosity is poorly resolved, however. In the second step, GRACE data are used to constrain the 3-D viscosity using spherical finite element modelling. In this case, the results also agree with past investigations, but GRACE data alone cannot discriminate between lateral heterogeneities in the mantle that are thermal in origin from those due to changes in chemical composition. More notably, we treat in detail GRACE-related questions such as implementation of an adequate Level-2 filter technique and identification of the best reduction method for hydrological mass change signals. It turns out that the Gaussian filter technique is the best for this type of investigation. Even the best global hydrology models used in GRACE investigations still fail to improve the mismatches—thus one should be careful not to blindly use them for ‘improving’ GIA models in North America or other centres of rebound. In conclusion, our study clearly shows that GRACE data greatly complement the study of GIA. As there are new GRACE releases in progress, and also in light of a new generation of ice history models, GRACE is likely to sharpen our insights concerning earth structure and rheology within the next few years.
Accurate measurements of glacier surface topography and their changes play an essential role in various glaciological studies related to glacier dynamics and mass balance. The focus of this study is on mapping glacier digital elevation model (DEM) and elevation changes in the western Qilian Mountains, northern Tibetan Plateau, by synergistically using the TanDEM-X (TDX) bistatic Interferometric Synthetic Aperture Radar (InSAR) data in 2013 and Shuttle Radar Topography Mission (SRTM) DEM in 2000. The first high-resolution and high-precision glacier DEM is derived in this region by a TDX InSAR procedure with a nonlocal (NL) filter. Validated against the Ice, Cloud, and land Elevation Satellite (ICESat) height references, the absolute height error of the TanDEM-X DEM derived with the NL filter and the Goldstein filter with the parameters investigated is, respectively, 1.493 ± 0.747 and 1.857 ± 1.709 m. Further, four combinations of differential phase method (DiffPha) and DEM differencing method (DiffDem) with Goldstein filter and NL filter are applied to estimate glacier elevation changes between 2000 and 2013. The synergistic use of the DiffPha method and the NL filter is superior to other three combinations in terms of uncertainty and noise reduction. Generally, a clear surface thinning can be found in most glacier tongue regions, the maximum value of elevation lowering up to approximately -40 m, whereas a slight thickening is detected in accumulation areas, which are in agreement with the height difference results between GPS measurements and SRTM DEM over Laohugou Glacier No.12. This study demonstrates the potential of the TanDEM-X bistatic InSAR in mapping surface topography and elevation changes of valley glaciers in the Tibetan Plateau.
'/%3e%3cuse xlink:href='%23a' transform='rotate(23.036 .093 25.536)'/%3e%3cuse xlink:href='%23a' transform='rotate(45.87 1.273 16.18)'/%3e%3cuse xlink:href='%23a' transform='rotate(69.945 .996 12.078)'/%3e%3cuse xlink:href='%23a' transform='rotate(20.66 -19.689 31.932)'/%3e%3c/svg%3e)
