基于SAR干涉点目标分析技术的城市地表形变监测

通过深入研究干涉点目标的相位模型,提出基于空间搜索的邻近点目标干涉相位差解缠方法,用以计算点目标的地形残差和线性形变,以及分离点目标大气延迟相位和非线性形变相位的时空域滤波方法,解决干涉点目标分析中的关键问题.最后,以苏州地区地表沉降监测为应用试验,利用形成的SAR干涉点目标形变信息提取技术,获取苏州市区1992-2002年间的地表沉降信息.研究结果与已有文献记录保持比较好的一致性,证明SAR干涉点目标技术完全可以发展成为应用于城市地表形变监测的实用化技术.

Abstract：

Interferometric point target analysis (IPTA) is one of the latest developments in radar interferometric processing. It is achieved by analysis of the interferometric phases of some individual point targets, which are discrete and present temporarily stable backscattering characteristics, in long temporal series of interferometric SAR images. This paper analyzes the interferometric phase model of point targets, and then addresses two key issues within IPTA process. Firstly, a spatial searching method is proposed to unwrap the interferometric phase difference between two neighboring point targets. The height residual error and linear deformation rate of each point target can then be calculated, when a global reference point with known height correction and deformation history is chosen. Secondly, a spatial-temporal filtering scheme is proposed to further separate the atmosphere phase and nonlinear deformation phase from the residual interferometric phase. Finally, an experiment of the developed IPTA methodology is conducted over Suzhou urban area. Totally 38 ERS-1/2 SAR scenes are analyzed, and the deformation information over 3 546 point targets in the time span of 1992-2002 are generated. The IPTA-derivecl deformation shows very good agreement with the published result, which demonstrates that the IPTA technique can be developed into an operational tool to map the ground subsidence over urban area.     

Variation in tree mortality and regeneration affect forest carbon recovery following fuel treatments and wildfire in the Lake Tahoe Basin, California, USA

Background

Forest fuel treatments have been proposed as tools to stabilize carbon stocks in fire-prone forests in the Western U.S.A. Although fuel treatments such as thinning and burning are known to immediately reduce forest carbon stocks, there are suggestions that these losses may be paid back over the long-term if treatments sufficiently reduce future wildfire severity, or prevent deforestation. Although fire severity and post-fire tree regeneration have been indicated as important influences on long-term carbon dynamics, it remains unclear how natural variability in these processes might affect the ability of fuel treatments to protect forest carbon resources. We surveyed a wildfire where fuel treatments were put in place before fire and estimated the short-term impact of treatment and wildfire on aboveground carbon stocks at our study site. We then used a common vegetation growth simulator in conjunction with sensitivity analysis techniques to assess how predicted timescales of carbon recovery after fire are sensitive to variation in rates of fire-related tree mortality, and post-fire tree regeneration.

Results

We found that fuel reduction treatments were successful at ameliorating fire severity at our study site by removing an estimated 36% of aboveground biomass. Treated and untreated stands stored similar amounts of carbon three years after wildfire, but differences in fire severity were such that untreated stands maintained only 7% of aboveground carbon as live trees, versus 51% in treated stands. Over the long-term, our simulations suggest that treated stands in our study area will recover baseline carbon storage 10?C35?years more quickly than untreated stands. Our sensitivity analysis found that rates of fire-related tree mortality strongly influence estimates of post-fire carbon recovery. Rates of regeneration were less influential on recovery timing, except when fire severity was high.

Conclusions

Our ability to predict the response of forest carbon resources to anthropogenic and natural disturbances requires models that incorporate uncertainty in processes important to long-term forest carbon dynamics. To the extent that fuel treatments are able to ameliorate tree mortality rates or prevent deforestation resulting from wildfire, our results suggest that treatments may be a viable strategy to stabilize existing forest carbon stocks.     

Use of the L2C signal for inversions of GPS radio occultation data in the neutral atmosphere

Results from processing FORMOSAT-3/COSMIC radio occultations (RO) with the new GPS L2C signal acquired both in phase locked loop (PLL) and open loop (OL) modes are presented. Analysis of L2P, L2C, and L1CA signals acquired in PLL mode shows that in the presence of strong ionospheric scintillation not only L2P tracking, but also L1CA tracking often fails, while L2C tracking is most stable. The use of L2C improves current RO processing in the neutral atmosphere mainly by increasing the number of processed occultations (due to significant reduction in the number of L2 tracking failures) and marginally by a reduction in noise in statistics. The latter is due to the combination of reduced L2C noise (compared to L2P) and increased L1CA noise in those occultations where L2P would have failed. This result suggests application of OL tracking for L1CA and L2C signals throughout an entire occultation to optimally acquire RO data. Two methods of concurrent processing of L1CA and L2C RO signals are considered. Based on testing of individual occultations, these methods allow: (1) reduction in uncertainty of bending angles retrieved by wave optics in the lower troposphere and (2) reduction in small-scale residual errors of the ionospheric correction in the stratosphere.     

EXPERIENCE IN COMPILING LARGE-SCALE MAPS OF EROSION-PRONE AREAS

Specific criteria used in the compilation of large-scale maps of erosion-prone areas are defined and evaluated for tracts of land representing a number of different agricultural cover types. Formulas are derived which facilitate the mapping of such areas based on such readily observable morphometric characteristics as slope length and steepness. Translated by Edward Torrey, Alexandria, VA 22308 from: Vestnik Moskovskogo Universiteta, geografiya, 1988, No. 5, pp. 25-31.     

Spatial-Correlation Based Persistent Scatterer Interferometric Study for Ground Deformation

Interferometric Synthetic Aperture Radar (InSAR), nowadays, is a precise technique for monitoring and detecting ground deformation at a millimetric level over large areas using multi-temporal SAR images. Persistent Scatterer Interferometric SAR (PSInSAR), an advanced version of InSAR, is an effective tool for measuring ground deformation using temporally stable reference points or persistent scatterers. We have applied both PSInSAR and Small Baseline Subset (SBAS) methods, based on the spatial correlation of interferometric phase, to estimate the ground deformation and time-series analysis. In this study, we select Las Vegas, Nevada, USA as our test area to detect the ground deformation along satellite line-of-sight (LOS) during November 1992–September 2000 using 44 C-band SAR images of the European Remote Sensing (ERS-1 and ERS-2) satellites. We observe the ground displacement rate of Las Vegas is in the range of ?19 to 8 mm/year in the same period. We also cross-compare PSInSAR and SBAS using mean LOS velocity and time-series. The comparison shows a correlation coefficient of 0.9467 in the case of mean LOS velocity. Along this study, we validate the ground deformation results from the satellite with the ground water depth of Las Vegas using time-series analysis, and the InSAR measurements show similar patterns with ground water data.     

Ionospheric and receiver DCB-constrained multi-GNSS single-frequency PPP integrated with MEMS inertial measurements

Single-frequency precise point positioning (SF-PPP) is a potential precise positioning technique due to the advantages of the high accuracy in positioning after convergence and the low cost in operation. However, there are still challenges limiting its applications at present, such as the long convergence time, the low reliability, and the poor satellite availability and continuity in kinematic applications. In recent years, the achievements in the dual-frequency PPP have confirmed that its performance can be significantly enhanced by employing the slant ionospheric delay and receiver differential code bias (DCB) constraint model, and the multi-constellation Global Navigation Satellite Systems (GNSS) data. Accordingly, we introduce the slant ionospheric delay and receiver DCB constraint model, and the multi-GNSS data in SF-PPP modular together. In order to further overcome the drawbacks of SF-PPP in terms of reliability, continuity, and accuracy in the signal easily blocking environments, the inertial measurements are also adopted in this paper. Finally, we form a new approach to tightly integrate the multi-GNSS single-frequency observations and inertial measurements together to ameliorate the performance of the ionospheric delay and receiver DCB-constrained SF-PPP. In such model, the inter-system bias between each two GNSS systems, the inter-frequency bias between each two GLONASS frequencies, the hardware errors of the inertial sensors, the slant ionospheric delays of each user-satellite pair, and the receiver DCB are estimated together with other parameters in a unique Kalman filter. To demonstrate its performance, the multi-GNSS and low-cost inertial data from a land-borne experiment are analyzed. The results indicate that visible positioning improvements in terms of accuracy, continuity, and reliability can be achieved in both open-sky and complex conditions while using the proposed model in this study compared to the conventional GPS SF-PPP.     

A two-step ionospheric modeling algorithm considering the impact of GLONASS pseudo-range inter-channel biases

The Global Navigation Satellite System presents a plausible and cost-effective way of computing the total electron content (TEC). But TEC estimated value could be seriously affected by the differential code biases (DCB) of frequency-dependent satellites and receivers. Unlike GPS and other satellite systems, GLONASS adopts a frequency-division multiplexing access mode to distinguish different satellites. This strategy leads to different wavelengths and inter-frequency biases (IFBs) for both pseudo-range and carrier phase observations, whose impacts are rarely considered in ionospheric modeling. We obtained observations from four groups of co-stations to analyze the characteristics of the GLONASS receiver P1P2 pseudo-range IFB with a double-difference method. The results showed that the GLONASS P1P2 pseudo-range IFB remained stable for a period of time and could catch up to several meters, which cannot be absorbed by the receiver DCB during ionospheric modeling. Given the characteristics of the GLONASS P1P2 pseudo-range IFB, we proposed a two-step ionosphere modeling method with the priori IFB information. The experimental analysis showed that the new algorithm can effectively eliminate the adverse effects on ionospheric model and hardware delay parameters estimation in different space environments. During high solar activity period, compared to the traditional GPS + GLONASS modeling algorithm, the absolute average deviation of TEC decreased from 2.17 to 2.07 TECu (TEC unit); simultaneously, the average RMS of GPS satellite DCB decreased from 0.225 to 0.219 ns, and the average deviation of GLONASS satellite DCB decreased from 0.253 to 0.113 ns with a great improvement in over 55%.     

Tropospheric refractivity and zenith path delays from least-squares collocation of meteorological and GNSS data

Precise positioning requires an accurate a priori troposphere model to enhance the solution quality. Several empirical models are available, but they may not properly characterize the state of troposphere, especially in severe weather conditions. Another possible solution is to use regional troposphere models based on real-time or near-real time measurements. In this study, we present the total refractivity and zenith total delay (ZTD) models based on a numerical weather prediction (NWP) model, Global Navigation Satellite System (GNSS) data and ground-based meteorological observations. We reconstruct the total refractivity profiles over the western part of Switzerland and the total refractivity profiles as well as ZTDs over Poland using the least-squares collocation software COMEDIE (Collocation of Meteorological Data for Interpretation and Estimation of Tropospheric Pathdelays) developed at ETH Zürich. In these two case studies, profiles of the total refractivity and ZTDs are calculated from different data sets. For Switzerland, the data set with the best agreement with the reference radiosonde (RS) measurements is the combination of ground-based meteorological observations and GNSS ZTDs. Introducing the horizontal gradients does not improve the vertical interpolation, and results in slightly larger biases and standard deviations. For Poland, the data set based on meteorological parameters from the NWP Weather Research and Forecasting (WRF) model and from a combination of the NWP model and GNSS ZTDs shows the best agreement with the reference RS data. In terms of ZTD, the combined NWP-GNSS observations and GNSS-only data set exhibit the best accuracy with an average bias (from all stations) of 3.7 mm and average standard deviations of 17.0 mm w.r.t. the reference GNSS stations.     

Group delay variations of GPS transmitting and receiving antennas

GPS code pseudorange measurements exhibit group delay variations at the transmitting and the receiving antenna. We calibrated C1 and P2 delay variations with respect to dual-frequency carrier phase observations and obtained nadir-dependent corrections for 32 satellites of the GPS constellation in early 2015 as well as elevation-dependent corrections for 13 receiving antenna models. The combined delay variations reach up to 1.0 m (3.3 ns) in the ionosphere-free linear combination for specific pairs of satellite and receiving antennas. Applying these corrections to the code measurements improves code/carrier single-frequency precise point positioning, ambiguity fixing based on the Melbourne–Wübbena linear combination, and determination of ionospheric total electron content. It also affects fractional cycle biases and differential code biases.