星载D-INSAR技术及初步应用——以西藏玛尼地震为例

近 10a来 ,干涉合成孔径雷达 (INSAR ,InterferometricSyntheticApertureRadar;简称 :干涉雷达 )测量技术取得了令世人瞩目的成绩 ,已成为极具有潜力的空间对地观测新技术。较详细地介绍了干涉合成孔径雷达、差分干涉合成孔径雷达 (D -INSAR ,DifferentialINSAR ;简称 :差分干涉雷达 )技术的基本原理 ,并以 1997年 11月 8日西藏玛尼地震为例 ,通过三通差分干涉处理 ,获取了玛尼地震前后的地表变形场。通过分析可知 ,变形梯度带与发震断层平行 ,均沿NEE -SWW(2 5 0°)分布 ,断层水平错距近 5m ,最大隆起斜距向位移量为 98cm ,最大沉降斜距向位移量为 95cm     

Gully and soil and water conservation structure densities in semi‐arid northern Ethiopia over the last 80 years

Gullies have been a common phenomenon in semi‐arid northern Ethiopia for the last centuries. However, soil and water conservation (SWC) structures have been implemented for a long time to curb soil erosion. Though, like most of the affected areas worldwide, density and distribution of gullies and SWC structures, their causes and interrelations are poorly understood. The aims of this study were to develop a technique for mapping these densities of gullies and SWC structures, to explain their spatial distribution and to analyze changes over the period 1935–2014. Aerial photographs from 1935 to 1936 and Google Earth images from 2014 of the 5142 km² Geba catchment were used. Transect lines were established to count gullies and SWC structures in order to calculate densities. On average, a gully density of 1.14 km km^?2 was measured in 1935–1936 of which the larger portion (75%) were vegetated, indicating they were not very active. Over 80 years, gully density has significantly increased to 1.59 km km^?2 with less vegetation growing in their channel, but 66% of these gullies were treated with check dams. There was c. 3 km km^?2 of indigenous SWC structures (daget or lynchets) in 1935–1936 whereas a high density (20 km km^?2) of introduced SWC structures (mainly stone bunds and terraces) were observed in 2014. The density of gullies is positively correlated with slope gradient and shrubland cover and negatively with cropland cover, whereas the density of SWC structures significantly increased with increasing cropland cover. Density maps of gullies and SWC structures indicate sensitive areas to gully formation and priority areas for the implementation of SWC structures in Geba catchment. The obtained results illustrate the feasibility of the methods applied to map the density of gullies and SWC structures in mountainous areas. Copyright © 2018 John Wiley & Sons, Ltd.     

干涉成像雷达(INSAR)技术及其应用现状

介绍了干涉成像雷达（ＩＮＳＡＲ）技术的发展历史、基本原理和关键技术,以及星载ＩＮＳＡＲ的处理流程,还介绍了主要的应用领域,特别是ＩＮＳＡＲ在确定兰德斯地震位移场中的应用。说明ＩＮＳＡＲ在三维信息获取及地壳形变监测方面具有良好的应用前景     

The INSAR technique: its principle and applications to mapping the deformation field of earthquakes

ｎｔｒｏｄｕｃｔｉｏｎＡｆｔｅｒｍｏｒｅｔｈａｎ３０ｙｅａｒｓ′ｄｅｖｅｌｏｐｍｅｎｔ,ｔｈｅｒｅｍｏｔｅｓｅｎｓｉｎｇｔｅｃｈｎｉｑｕｅｈａｓｍｕｔｕａｌｙｉｎｆｉｌｔｒａｔｅｄａｎｄｍｉｘｅｄｗｉｔｈｔｈｏｓｅｄｉｓｃｉｐｌｉｎｅｓｓｕｃｈａｓ...     

Northridge, California, earthquake of 1994: density of red-tagged buildings versus peak horizontal velocity and intensity of shaking

The 1994 Northridge earthquake occurred underneath a densely populated metropolitan area, and was recorded by over 200 strong motion stations in the metropolitan area and vicinity. This rare coincidence made it an ideal case to study, in statistical sense, the correlation of damage to structures with the level of strong shaking, in particular with respect to (1) instrumental characteristics of shaking and (2) the reported site intensity scale. In this paper, statistics for the incidence of red-tagged building in 1 × 1 km² blocks in San Fernando Valley and Los Angeles is presented and analyzed, as function of the observed peak ground velocity or the local intensity of shaking. The ‘observed’ peak velocity is estimated from contour maps based on the recorded strong motion. The intensity of shaking is estimated from the published intensity map and from our modification of this map to make it more consistent with observed high damage to buildings in some localized areas. Finally, empirical scaling equations are derived which predict the average density of red-tagged buildings (per km²) as a function of peak ground velocity or site intensity of shaking. These scaling equations are specific to the region studied, and apply to Wooden Frame Construction, typical of post World War II period, which is the prevailing building type in the sample studied. These can be used to predict the density of red-tagged buildings per km² in San Fernando Valley and in Los Angeles for a scenario earthquake or for an ensemble of earthquakes during specified exposure, within the framework of probabilistic seismic hazard analysis. Such predictions will be useful to government officials for emergency planning, to the insurance industry for realistic assessment of insured losses, and to structural engineers for assessment of the overall performance of this type of buildings.     

Empirically-based modeling of spatial sampling uncertainties associated with rainfall measurements by rain gauges

In the quantitative evaluation of radar-rainfall products (maps), rain gauge data are generally used as a good approximation of the true ground rainfall. However, rain gauges provide accurate measurements for a specific location, while radar estimates represent areal averages. Because these sampling discrepancies could introduce noise into the comparisons between these two sensors, they need to be accounted for. In this study, the spatial sampling error is defined as the ratio between the measurements by a single rain gauge and the true areal rainfall, defined as the value obtained by averaging the measurements by an adequate number of gauges within a pixel. Using a non-parametric scheme, the authors characterize its full statistical distribution for several spatial (4, 16 and 36 km²) and temporal (15 min and hourly) scales.     

福建省泉州地区断裂带地壳形变PS-InSAR监测

PS-InSAR技术能够有效降低常规差分干涉雷达受时间失相干、空间失相干和大气效应的影响,在常规D-InSAR不能形成干涉条纹的情况下,可利用时间序列的雷达影像和相位稳定的永久散射体目标点获取离散的PS像素点形变速率。以福建省泉州地区的断裂带为研究对象,对1996-1999年的22景ERS SAR数据进行PS-InSAR处理,得出研究区主要断裂的视线向位移速率为3~5mm/a,表明该区断裂仍有一定的活动性,具有潜在的地震危险。     

基于小基线DInSAR技术监测太原市2003～2009年地表形变场

基于高相干点目标反演缓慢地表形变已成为当前DInSAR技术的研究热点.本文通过融合PS方法和相干目标法优点,采用小基线DInSAR技术提取城市地表形变场,并重点分析了地表线性形变的反演.在此基础上,以太原市为研究区,利用23景ENVISAT ASAR影像,提取了该市2003～2009年的地表形变场.研究结果表明:(1)...     

Small scale spatial variability of snow density and depth over complex alpine terrain: Implications for estimating snow water equivalent

This study analyzes spatial variability of snow depth and density from measurements made in February and April of 2010 and 2011 in three 1–2 km² areas within a valley of the central Spanish Pyrenees. Snow density was correlated with snow depth and different terrain characteristics. Regression models were used to predict the spatial variability of snow density, and to assess how the error in computed densities might influence estimates of snow water equivalent (SWE).The variability in snow depth was much greater than that of snow density. The average snow density was much greater in April than in February. The correlations between snow depth and density were generally statistically significant but typically not very high, and their magnitudes and signs were highly variable among sites and surveys. The correlation with other topographic variables showed the same variability in magnitude and sign, and consequently the resulting regression models were very inconsistent, and in general explained little of the variance. Antecedent climatic and snow conditions prior to each survey help highlight the main causes of the contrasting relation shown between snow depth, density and terrain. As a consequence of the moderate spatial variability of snow density relative to snow depth, the absolute error in the SWE estimated from computed densities using the regression models was generally less than 15%. The error was similar to that obtained by relating snow density measurements directly to adjacent snow depths.     

An intercomparison of remotely sensed soil moisture products at various spatial scales over the Iberian Peninsula

Soil moisture (SM) can be retrieved from active microwave (AM), passive microwave (PM) and thermal infrared (TIR) observations, each having unique spatial and temporal coverages. A limitation of TIR‐based retrievals is a dependence on cloud‐free conditions, whereas microwave retrievals are almost all weather proof. A downside of SM retrievals from PM is the coarse spatial resolution. Although SM retrievals at coarse spatial resolution proved to be valuable for global‐scale and continental‐scale studies, their value for regional‐scale studies remains limited. To increase the use of SM retrievals from PM observations, an existing method to enhance their spatial resolution was applied. We present an intercomparison study over the Iberian Peninsula for three SM products on two different spatial sampling grids. The remotely sensed SM products were also compared with in situ observations from the Remedhus network. Variations between ground data and satellite‐based SM are observed; all three remotely sensed SM products show good agreement to the ground observations. The comparison shows that these ground observations and satellite data are consistent, based on the correlation coefficient (R) and root mean square error (RMSE). The remotely sensed products were intercompared after sampling at 25 × 25 km² and after applying the smoothing filter‐based intensity modulation (SFIM) downscaling technique at 10 × 10 km² grids. After the application of the SFIM technique, the SM retrievals from PM observations show better agreement with the other remotely sensed SM products for approximately 40% of the study area. For another 40% of the study area, we found a similar agreement between these product combinations, whereas in extreme environments, both arid and densely vegetated regions, the agreement decreases after the application of the SFIM technique. Agreement between retrievals of absolute SM content from PM and TIR observations is generally high (R = 0.77 for semi‐arid areas). This study enhances our understanding of the remotely sensed SM products for improvements of SM retrieval and merging strategies. Copyright © 2013 John Wiley & Sons, Ltd.     

A new approach for estimating seismic damage of buried water supply pipelines

Conventional damage prediction methods for lifeline structures are primarily based on peak ground motion measurements. However, line structures such as lifelines suffer damage that is mainly induced by the strain of the ground and therefore are likely to be vulnerable to sharp spatial changes in the ground motion. In this study, we propose a measure for evaluating the damage incurred by underground water supply pipelines based on the spatial gradient of the peak ground velocity (PGV), in an attempt to quantify the effects of the geospatial variabilities in the ground motion on pipeline damage. We investigated the spatial distribution of the damage caused to water pipelines during the Niigata‐ken Chuetsu earthquake on October 10, 2004 (Japan Meteorological Agency magnitude (M_JMA) of 6.8) and the Kobe earthquake on January 17, 1995 (M_JMA7.3) and compared the surveyed damage with the PGV distribution as well as with the gradients of the PGV calculated around the damage areas. For the Kobe earthquake, we used the PGV distribution obtained by the strong‐motion simulation performed by Matsushima and Kawase 1 . In case of the Chuetsu earthquake, we estimated the ground motion using a broadband‐frequency‐based strong‐ground‐motion simulation method based on a multiasperity source model. In both cases, we calculated the gradients of the PGV along the geographical coordinates, with the amplitude of the PGV gradient vector being employed as the damage estimator. Our results show that the distribution of damage to underground water supply pipelines exhibits a greater correlation with the gradients of the PGV than with the PGV itself. Thus, the gradient of the PGV is a useful index for preparing initial‐screening hazard maps of underground facilities. Copyright © 2017 John Wiley & Sons, Ltd.     

基于ShakeMap_CNST的2014年云南鲁甸MS6.5地震震动图研究

分析了2014年8月3日云南鲁甸6.5级地震的发生背景及震害特征,介绍了ShakeMap_CNST震动图系统的设计思路和实现原理及地震动参数的校正方法,并将其应用于鲁甸6.5级地震的震动图预测中,进一步分析了震区场地条件对局部震害的影响。结果表明,震动图预测的地震动特征与现场宏观调查的结果总体上是相适应的,最后对震动图系统的优缺点及应用前景进行了探讨。     

DInSAR measurements of ground deformation by sinkholes,mining subsidence,and landslides,Ebro River,Spain

Differential Interferometric Synthetic Aperture Radar (DInSAR) has been applied to detect and measure ground deformation in a stretch of the Ebro River valley (Spain) excavated in salt‐bearing evaporites. The capability of the Small Baseline Subset (SBAS) DInSAR technique to detect ground displacement is analyzed comparing the DInSAR results with the available geomorphological information. The deformation map derived from 27 European Remote Sensing (ERS) satellite images covering more than five years provides sub‐centimeter displacement measurements in zones coinciding with known active sinkholes and landslides. Moreover the map provides the first account of mining subsidence in the area. The measured deformation rates reach 1·68 cm/y for the sinkholes, 0·80 cm/y for the landslides and 1·45 cm/y for the area affected by mining subsidence. The SBAS DInSAR technique provided deformation measurements in a small proportion (5–10%) of the known active sinkholes and landslides. This limitation is mainly due to the lack of coherence in agricultural areas, the spatial resolution of the deformation map (pixel size of 90 m), and the parallelism between the ERS satellite line of sight and the linear escarpment on which most of the landslides occur. Despite this, the interferometric technique provides valuable data that complement traditional geomorphological studies including the quantification of the deformation phenomena, the identification of mining subsidence otherwise only recognizable by geodetic methods, and the detection of creep deformation which might correspond to premonitory indicators of catastrophic sinkholes and landslides capable of causing the loss of lives. Detailed DInSAR studies combined with field data would be required to improve the analysis of each deformation area. Copyright © 2009 John Wiley & Sons, Ltd.     

Coseismic Three-Dimensional Displacements Determined Using SAR Data: Theory and an Application Test

We present a methodology to construct three-dimensional deformation maps using different parameters that can be extracted using SAR data. We apply the methodology using ENVISAT SAR data before and after the December 26^th, 2003 M_w 6.6 Bam earthquake in Iran to determine spatial quasi-continuous three-dimensional coseismic deformation maps. Two near vertical deformation offset components are computed using Envisat ASAR differential interferometry (DInSAR), while two horizontal deformation offset components are measured using sub-pixel correlation techniques applied to ASAR amplitude images. Using the presented methodology, we combine four unevenly precise independent projections of surface coseismic deformation to obtain the full three-dimensional coseismic deformation field caused by an earthquake. The full 3-D coseismic displacement vector is modeled using elastic deformation models constraining details of the fault geometry and slip distribution at depth. Results from the inversion are interpreted in the framework of the Iranian present-day tectonism. Full exploitation of dense 3-D coseismic surface deformation using SAR data, even for moderate earthquakes, should facilitate distinguishing between different interpretations of the mechanical properties of seismically active areas and within the inherent ambiguity of the geophysical inversion solutions.     

InSAR技术在监测形变中的干涉条件分析

针对InSAR技术在形变监测领域中的应用,本文从分析D-InSAR的干涉机理对其约束条件出发,重点研究空间去相关和形变梯度这两个关键指标对形变信息提取的影响.研究针对不同系统参数的传感器来展开,并论证这两个因素对干涉测量的限制.本文的创新性工作主要集中在以下几个方面: (1)分析了各传感器参数下的临界坡度范围,进而确立了SAR传感器监测的无信息区域;(2)证明了在植被覆盖区,相干性分解必须考虑体散射分量,并依此对模型进行了修正;(3)明确了形变梯度与相位梯度的差异以及两者的函数关系,并将入射角引入形变梯度函数模型.(4)通过模拟数据和真实数据比较各传感器在不同条件下形变监测的优缺点,在统计分析的基础上进一步验证结论.本文的研究是对干涉测量条件的进一步补充,对于理解干涉成像机理、完善InSAR数值计算模型(形变梯度模型和相干性分解模型)以及帮助研究人员根据研究区域特征选择合适的InSAR数据具有一定的参考价值.     

Are spatial patterns of soil moisture at plot scales generalisable across catchments,climates, and other characteristics? A synthesis of synoptic soil moisture across the Mid-Atlantic

The spatial variation of soil moisture over very small areas (<100 m²) can have nonlinear impacts on cycling and flux rates resulting in bias if it is not considered, but measuring this variation is difficult over extensive temporal and spatial scales. Most studies examining spatial variation of soil moisture were conducted at hillslope (0.01 km²) to multi-catchment spatial scales (1000 km²). They found the greatest variation at mid wetness levels and the smallest variation at wet and dry wetness levels forming a concave down relationship. There is growing evidence that concave down relationships formed between spatial variation of soil moisture and average soil moisture are consistent across spatial scales spanning several orders of magnitude, but more research is needed at very small, plot scales (<100 m²). The goal of this study was to characterise spatial variation in shallow soil moisture at the plot scale by relating the mean of measurements collected in a plot to the standard deviation (SD). We combined data from a previous study with thousands of new soil moisture measurements from 212 plots in eight catchments distributed across the US Mid-Atlantic Region to (1) test for a generalisable mean–SD relationship at plot scales, (2) characterise how landcover, land use, season, and hillslope position contribute to differences in mean–SD relationships, and (3) use these generalised mean–SD relationships to quantify their impacts on catchment scale nitrification and denitrification potential. Our study found that 98% of all measurements formed a generalised mean–SD relationship like those observed at hillslope and catchment spatial scales. The remaining 2% of data comprised a mean–SD relationship with greater spatial variation that originated from two riparian plots reported in a previous study. Incorporating the generalised mean–SD relationship into estimates of nitrification and denitrification potential revealed strong bias that was even greater when incorporating mean–SD observations from the two riparian plots with significantly greater spatial variation.     

Use of Geophysical and Remote Sensing Data for Assessment of Aquifer Depletion and Related Land Deformation

An integrated approach [field, Interferometric Synthetic Aperture Radar (InSAR), hydrogeology, geodesy, and spatial analysis] was adopted to identify the nature, intensity, and spatial distribution of deformational features (sinkholes, fissures, differential settling) reported over fossil aquifers in arid lands, their controlling factors, and possible remedies. The Lower Mega Aquifer System (area 2 × 10⁶ km²) in central and northern Arabia was used as a test site. Findings suggest that excessive groundwater extraction from the fossil aquifer is the main cause of deformation: (1) deformational features correlated spatially and/or temporally with increased agricultural development and groundwater extraction, and with a decline in water levels and groundwater storage (? 3.7 ± 0.6 km³/year); (2) earthquake events (years 1985–2016; magnitude 1–5) are largely (65% of reported earthquakes) shallow (1–5 km) and increased from 1 event/year in the early 1980s (extraction 1 km³/year), up to 13 events/year in the 1990s (average annual extraction > 6.4 km³). Results indicate that faults played a role in localizing deformation given that deformational sites and InSAR-based high subsidence rates (? 4 to ? 15 mm/year) were largely found within, but not outside of, NW–SE-trending grabens bound by the Kahf fault system. Findings from the analysis of Gravity Recovery and Climate Experiment solutions indicate that sustainable extraction could be attained if groundwater extraction was reduced by 3.5–4 km³/year. This study provides replicable and cost-effective methodologies for optimum utilization of fossil aquifers and for minimizing deformation associated with their use.     

基于Sentinel-1 SAR数据的黑河上游冻土形变时序InSAR监测

多年冻土活动层变化导致冻土区大范围地面变形,严重破坏区域内基础设施和水文地质条件,亟需加强活动层季节冻融过程的观测研究.本文提出一种基于分布式目标的小基线集时序InSAR(DSs-SBAS)的冻土形变监测方法.该方法采用分布式目标提取和特征值分解算法,并结合基于地温-形变约束关系的参考点选取新策略,提高了冻土形变监测结果的时空分辨率和可靠性.以祁连山黑河西支源头的野牛沟为研究区域,通过对27景Sentinel-1SAR影像进行时序InSAR分析,获取了2014-2016年该区多年冻土的形变时间序列和年均形变速率,并利用Stefan模型联合地温数据估算其季节性形变幅度.实地踏勘和结果分析表明:(1)研究区大部分多年冻土处于稳定状态(-1.0^+1.0cm·a^-1),在地形陡峭的南坡边缘及含冰量丰富的野牛沟河上游两侧沟底部分区域存在较大形变;(2)区域内冻土形变时间序列呈现年周期变化,冻土冻融形变存在季节性周期形变和季节性波动下沉两种形变特征,形变幅度和速率最大可达6.0cm和-3.0cm·a^-1;(3)不同区域的活动层冻结/融化始日和冻土形变存在明显差异,主要和冻土地貌、土壤类型以及活动层厚度有关.本文提出的方法在青藏高原多年冻土区大范围冻融监测和活动层厚度反演研究方面具有很大的应用潜力.     

Spatial prediction of temporal soil moisture dynamics using HYDRUS‐1D

This article investigates the soil moisture dynamics within two catchments (Stanley and Krui) in the Goulburn River in NSW during a 3‐year period (2005–2007) using the HYDRUS‐1D soil water model. Sensitivity analyses indicated that soil type, and leaf area index were the key parameters affecting model performance. The model was satisfactorily calibrated on the Stanley microcatchment sites with a single point rainfall record from this microcatchment for both surface 30 cm and full‐profile soil moisture measurements. Good correlations were obtained between observed and simulated soil water storage when calibrations for one site were applied to the other sites. We extended the predictions of soil moisture to a larger spatial scale using the calibrated soil and vegetation parameters to the sites in the Krui catchment where soil moisture measurement sites were up to 30 km distant from Stanley. Similarly good results show that it is possible to use a calibrated soil moisture model with measurements at a single site to extrapolate the soil moisture to other sites for a catchment with an area of up to 1000 km² given similar soils and vegetation and local rainfall data. Site predictions were effectively improved by our simple data assimilation method using only a few sample data collected from the site. This article demonstrates the potential usefulness of continuous time, point‐scale soil moisture data (typical of that measured by permanently installed TDR probes) and simulations for predicting the soil wetness status over a catchment of significant size (up to 1000 km²). Copyright © 2012 John Wiley & Sons, Ltd.