固体潮位移InSAR相位模拟及对广域地表形变监测的影响初探

固体地球潮汐(Solid earth tide,SET)在地表产生的径向位移可达40～50 cm,形变梯度可达2 cm/100 km,是各种精密大地测量和地球物理观测必须考虑的因素之一.随着星载合成孔径雷达干涉测量(Interferometric synthetic aperture radar,InSAR)地表形变监测范围的不断增大以及对精度要求不断提高,固体潮位移对InSAR观测的影响不容忽视.本文利用固体潮位移理论模型,根据InSAR测量基本原理和Sentinel-1卫星成像参数,模拟了固体潮位移InSAR相位,定量分析了其时空分布特征,并以我国江汉平原和北美大平原的Sentinel-1数据为例,探讨了固体潮位移InSAR相位对广域地表形变监测的影响.结果表明:(1)固体潮位移对InSAR广域地表形变监测存在较大影响,在250 km×250 km范围中,以C波段为例,其相位变化可达12 rad(对应52 mm视线向形变);(2)固体潮位移相位在中低纬度(60°S—60°N)地区变化较大,两极地区较小,且在时间上具有明显的周期性;(3)在Sentinel-1 InSAR观测中,通过固体潮位移相位改正去除了干涉图中的部分低频相位偏差,相比原始干涉图,改正后的解缠相位标准差减小了约29%.本研究对于认识固体潮位移InSAR相位的时空分布特征以及提高星载InSAR广域地表形变监测的准确度与可靠性均具有重要意义.     

Toward a High-Resolution Monitoring of Continental Surface Water Extent and Dynamics,at Global Scale: from GIEMS (Global Inundation Extent from Multi-Satellites) to SWOT (Surface Water Ocean Topography)

Up to now, high-resolution mapping of surface water extent from satellites has only been available for a few regions, over limited time periods. The extension of the temporal and spatial coverage was difficult, due to the limitation of the remote sensing technique [e.g., the interaction of the radiation with vegetation or cloud for visible observations or the temporal sampling with the synthetic aperture radar (SAR)]. The advantages and the limitations of the various satellite techniques are reviewed. The need to have a global and consistent estimate of the water surfaces over long time periods triggered the development of a multi-satellite methodology to obtain consistent surface water all over the globe, regardless of the environments. The Global Inundation Extent from Multi-satellites (GIEMS) combines the complementary strengths of satellite observations from the visible to the microwave, to produce a low-resolution monthly dataset (\(0.25^\circ \,\times \,0.25^\circ\)) of surface water extent and dynamics. Downscaling algorithms are now developed and applied to GIEMS, using high-spatial-resolution information from visible, near-infrared, and synthetic aperture radar (SAR) satellite images, or from digital elevation models. Preliminary products are available down to 500-m spatial resolution. This work bridges the gaps and prepares for the future NASA/CNES Surface Water Ocean Topography (SWOT) mission to be launched in 2020. SWOT will delineate surface water extent estimates and their water storage with an unprecedented spatial resolution and accuracy, thanks to a SAR in an interferometry mode. When available, the SWOT data will be adopted to downscale GIEMS, to produce a long time series of water surfaces at global scale, consistent with the SWOT observations.     

Timely Low Resolution SAR Imagery To Support Floodplain Modelling: a Case Study Review

It is widely recognised that remote sensing can support flood monitoring, modelling and management. In particular, satellites carrying Synthetic Aperture Radar (SAR) sensors are valuable as radar wavelengths can penetrate cloud cover and are insensitive to daylight. However, given the strong inverse relationship between spatial resolution and revisit time, monitoring floods from space in near real time is currently only possible through low resolution (about 100 m pixel size) SAR imagery. For instance, ENVISAT-ASAR (Advanced Synthetic Aperture Radar) in WSM (wide swath mode) revisit times are of the order of 3 days and the data can be obtained within 24 h at no (or low) cost. Hence, this type of space-borne data can be used for monitoring major floods on medium-to-large rivers. This paper aims to discuss the potential for, and uncertainties of, coarse resolution SAR imagery to monitor floods and support hydraulic modelling. The paper first describes the potential of globally and freely available space-borne data to support flood inundation modelling in near real time. Then, the uncertainty of SAR-derived flood extent maps is discussed and the need to move from deterministic binary maps (wet/dry) of flood extent to uncertain flood inundation maps is highlighted.     

Oil spill detection on the ocean surface using hybrid polarimetric SAR imagery

Hybrid-polarimetric SAR (synthetic aperture radar) is a new SAR mode, with relatively simple architecture, low cost, and wide swath, which will be carried by several Earth-observing systems from now to the near future. Here, we show how the second Stokes parameter of hybrid-polarimetric SAR can be employed to detect oil on the ocean surface using the classic well-known Otsu threshold methodology, in relation to contributions from different polarizations and dampening effects on backscatter intensity, neglecting the specific scattering mechanisms and oil types for an oil-covered surface. The detection methodology is demonstrated to be reliable in three example cases: oil-on-water experiments conducted by the Norwegian Clean Seas Association, natural oil seeps from the Gulf of Mexico, and observations from the Deep Water Horizon oil spill disaster in 2010.     

Statistical interpolation of ozone measurements from satellite data (TOMS, SBUV and SAGE II) using the kriging method

This study demonstrates that ordinary kriging in spherical coordinates using experimental semi-variograms provides highly usable results, especially near the pole in winter and/or where there could be data missing over large areas. In addition, kriging allows display of the spatial variability of daily ozone measurements at different pressure levels. Three satellite data sets were used: Total Ozone Mapping Spectrometer (TOMS) data, Solar Backscattered UltraViolet (SBUV), and the Stratospheric Aerosol and Gas Experiment (SAGE II) ozone profiles. Since SBUV is a nadir-viewing instrument, measurements are only taken along the sun-synchronous polar orbits of the satellite. SAGE II is a limb-viewing solar occultation instrument, and measurements have high vertical resolution but poor daily coverage. TOMS has wider coverage with equidistant distribution of data (resolution 1° × 1.25°) but provides no vertical information. Comparisons of the resulting SBUV-interpolated (column-integrated) ozone field with TOMS data are strongly in agreement, with a global correlation of close to 98%. Comparisons of SBUV-interpolated ozone profiles with daily SAGE II profiles are relatively good, and comparable to those found in the literature. The interpolated ozone layers at different pressure levels are shown.     

基于星载被动微波遥感的青藏高原湖冰物候监测方法

湖冰物候影响着区域及全球气候,是全球变化的敏感因子,青藏高原湖泊众多,冻融现场监测数据缺乏,而微波具有对冰水相变敏感、时间分辨率高、历史存档数据长等特点,这对于长时间序列湖冰物候研究具有重要意义.然而,被动微波遥感空间分辨率低、湖泊亮温的精准定位难.论文通过获取AMSR-E/Aqua和AMSR-2/Gcom-W1的亮温数据,构建了基于轨道亮温数据的阈值判别法,通过对青藏高原不同区域和不同大小的青海湖、色林错、哈拉湖以及阿其克库勒湖进行测试研究:与青海湖现场观测对比,湖泊完全冻结日期与开始融化日期最大误差小于3天;与无云光学遥感判别结果相比,4个湖泊的冻融参数误差为2~4天.结果表明,被动微波轨道亮温数据可实现青藏高原地区亚像元级中大型湖泊冻融信息的获取,历史卫星资料可为湖冰物候的监测提供重要的支撑.     

The Representation of Tropospheric Water Vapor Over Low-Latitude Oceans in (Re-)analysis: Errors,Impacts, and the Ability to Exploit Current and Prospective Observations

This paper addresses the representation of lower tropospheric water vapor in the meteorological analyses—fully detailed estimates of atmospheric state—providing the wide temporal and spatial coverage used in many process studies. Analyses are produced in a cycle combining short forecasts from initial conditions with data assimilation that optimally estimates the state of the atmosphere from the previous forecasts and new observations, providing initial conditions for the next set of forecasts. Estimates of water vapor are among the less certain aspects of the state because the quantity poses special challenges for data assimilation while being particularly sensitive to the details of model parameterizations. Over remote tropical oceans observations of water vapor come from two sources: passive observations at microwave or infrared wavelengths that provide relatively strong constraints over large areas on column-integrated moisture but relatively coarse vertical resolution, and occultations of Global Positioning System provide much higher accuracy and vertical resolution but are relatively spatially coarse. Over low-latitude oceans, experiences with two systems suggest that current analyses reproduce much of the large-scale variability in integrated water vapor but have systematic errors in the representation of the boundary layer with compensating errors in the free troposphere; these errors introduce errors of order 10% in radiative heating rates through the free troposphere. New observations, such as might be obtained by future observing systems, improve the estimates of water vapor but this improvement is lost relatively quickly, suggesting that exploiting better observations will require targeted improvements to global forecast models.     

Accuracy assessment of the MODIS snow products

A suite of Moderate‐Resolution Imaging Spectroradiometer (MODIS) snow products at various spatial and temporal resolutions from the Terra satellite has been available since February 2000. Standard products include daily and 8‐day composite 500 m resolution swath and tile products (which include fractional snow cover (FSC) and snow albedo), and 0·05° resolution products on a climate‐modelling grid (CMG) (which also include FSC). These snow products (from Collection 4 (C4) reprocessing) are mature and most have been validated to varying degrees and are available to order through the National Snow and Ice Data Center. The overall absolute accuracy of the well‐studied 500 m resolution swath (MOD10_L2) and daily tile (MOD10A1) products is ～93%, but varies by land‐cover type and snow condition. The most frequent errors are due to snow/cloud discrimination problems, however, improvements in the MODIS cloud mask, an input product, have occurred in ‘Collection 5’ reprocessing. Detection of very thin snow (<1 cm thick) can also be problematic. Validation of MOD10_L2 and MOD10A1 applies to all higher‐level products because all the higher‐level products are all created from these products. The composited products may have larger errors due, in part, to errors propagated from daily products. Recently, new products have been developed. A fractional snow cover algorithm for the 500 m resolution products was developed, and is part of the C5 daily swath and tile products; a monthly CMG snow product at 0·05° resolution and a daily 0·25° resolution CMG snow product are also now available. Similar, but not identical products are also produced from the MODIS on the Aqua satellite, launched in May 2002, but the accuracy of those products has not yet been assessed in detail. Published in 2007 by John Wiley & Sons, Ltd.     

Mantle structure from inter-station Rayleigh wave dispersion and its tectonic implication in western China and neighboring regions

We processed more than 3000 inter-station great circle paths to determine the phase velocity for the fundamental mode of Rayleigh wave, and finally arrived at 110 paths of high quality dispersion data, which show good spatial coverage in western China and neighboring regions. Rayleigh wave phase velocity dispersion model WChina1D was obtained and compared with previous global and regional models. Phase velocity maps from 15 to 120 s were inverted and the maps of 20, 40, 80, and 120 s are presented in this paper. Checkerboard tests show the average lateral resolution in our area of interest is about 7°. Our tomographic results corroborate a prominent low-velocity anomaly lying mainly in the lower crust and uppermost mantle in the Chang Thang terrane. The apparent low-velocity anomaly also appears in the wide area of northeastern Tibet in the crust and upper mantle. The low-velocity area around southeastern Tibet may be created by the southeastern migration of the low-velocity mass of the Tibetan plateau. The eastern Tarim shows structure with higher velocities relative to that of central Tarim. A large-scale low-velocity anomaly is clearly seen in central and western Mongolia. Our high quality measurements were also used to evaluate the CUB global shear velocity model [Shapiro, N., Ritzwoller, M., 2002. Monte-Carlo inversion for a global shear-velocity model of the crust and upper mantle. Geophys. J. Int. 151, 88-105] of the crust and upper mantle. The 40 s Rayleigh phase velocity map predicted from CUB model shows an apparent discrepancy with our measurements in western China and western Mongolia, which implies a higher estimated (about +1-2%) phase velocity model in these regions, probably due to the Gaussian smoothing condition in their tomography inversion.     

降雨对大柏舍地电台地电阻率观测的影响

选取大柏舍地电台2011—2019年地电阻率观测数据和降雨量资料,分析降雨对地电阻率观测的影响,利用影响系数理论,对结果进行分析探讨。结果表明,由于大柏舍地电台表层介质影响系数较小且为负值,小幅降雨对该台地电阻率短期影响不明显,仅自然电位差受到干扰;大范围强降雨时,地电阻率和自然电位差均受到影响;因ZD8B仪地电阻率观测系统架空线路漏电,导致EW向分量呈阶降变化,且比用于背景场观测的ZD8M仪变化幅度大。     

青藏高原及其东南前陆地壳运动的GPS监测

作为世界上最高最大最活动的青藏高原,以其独特的自然地理条件吸引了广大学者和科研人员的关注,在此地区开展的各项研究工作.本文主要介绍了利用GPS监测青藏高原及邻区所得到的速度矢量场.观测结果表明：高原南北向收缩、高原物质东流、高原东部地区地壳物质相对喜马拉雅汇聚区作大规模的顺时针旋转等高原动力学特征.并就其产生原因进行了初步分析.     

Requirements of a global near-surface soil moisture satellite mission: accuracy,repeat time,and spatial resolution

Soil moisture satellite mission accuracy, repeat time and spatial resolution requirements are addressed through a numerical twin data assimilation study. Simulated soil moisture profile retrievals were made by assimilating near-surface soil moisture observations with various accuracy (0, 1, 2, 3, 4, 5 and 10%v/v standard deviation) repeat time (1, 2, 3, 5, 10, 15, 20 and 30 days), and spatial resolution (0.5, 6, 12 18, 30, 60 and 120 arc-min). This study found that near-surface soil moisture observation error must be less than the model forecast error required for a specific application when used as data assimilation input, else slight model forecast degradation may result. It also found that near-surface soil moisture observations must have an accuracy better than 5%v/v to positively impact soil moisture forecasts, and that daily near-surface soil moisture observations achieved the best soil moisture and evapotranspiration forecasts for the repeat times assessed, with 1–5 day repeat times having the greatest impact. Near-surface soil moisture observations with a spatial resolution finer than the land surface model resolution (∼30 arc-min) produced the best results, with spatial resolutions coarser than the model resolution yielding only a slight degradation. Observations at half the land surface model spatial resolution were found to be appropriate for our application. Moreover, it was found that satisfying the spatial resolution and accuracy requirements was much more important than repeat time.     

The spatial structure of local surficial sediment characteristics on Georges Bank,USA

Georges Bank is one of the world’s most productive marine ecosystems, but the lack of accurate broad-scale sediment maps presently limits habitat assessments and spatial fisheries management. From 1999 to 2009 we surveyed 36,669 km² of Georges Bank using 2.8 and 0.6 m² quadrats viewed with live underwater video (video quadrats). The sediment types observed in 61,604 quadrats were used to map and evaluate spatial structure of local surficial sediment coarseness, dominance, heterogeneity, and maximum size characteristics at a 1 km² spatial resolution. Sand dominated sediment covered 62% of the study area, and there was a logarithmic decline in coverage by larger, coarser and more heterogeneous sediments. Gravel dominated sediments covered 38% of the study area and were more than twice as abundant as previously estimated. A 12,890 km² swath of gravel dominated seabed stretched from Cape Cod to northeastern Georges Bank consistent with estimates of prehistoric glacial extent. Within the swath there were 14 large gravel outcrops (15–2743 km²) . This work increases the spatial resolution of sediment information available for habitat assessments and spatial fisheries management on Georges Bank by two orders of magnitude. The four sediment characteristics we evaluated support further detailed investigations of the Bank’s benthos, including the influences of surficial sediment characteristics on species and community distributions, and more spatially accurate estimates of seabed roughness. Finally, this work demonstrates the use of video quadrats as an alternative to traditional grab sampling and modern acoustic sampling for continental shelf-scale mapping.     

中国大陆及邻区Rayleigh面波相速度分布特征

本文根据102个数字化台站记录的长周期垂直向面波资料,利用双台互相关方法测量了538条独立路径的基阶Rayleigh面波相速度频散资料,反演获得了中国大陆及边邻地区（70°E～140°E,18°N～55°N）20～120 s（周期间隔为5 s） 共21个周期的Rayleigh波相速度空间分布图像. 检测板测试结果显示中国大陆中东部地区横向分辨率可达3°,而西部及边邻地区大约5°. 研究表明,中国大陆地区的Rayleigh波相速度分布横向差异显著,大致以104°E为界,可分成具有不同速度结构特征的东、西两部分. 一般较短周期（20～35 s）的相速度分布受地形和地壳厚度的影响较大,总体表现为东部速度高,西部速度低;塔里木盆地、青藏地块及其东缘的松潘-甘孜地块形成整个研究区内最为突出的低速异常体,蒙古西部低速特征也较清晰;东部的四川盆地、扬子地块、华南地块、松辽盆地、日本海及蒙古东部高速特征明显. 随着周期的增大,青藏地块中部的低速异常体横向尺度逐渐缩小,而喜马拉雅冲断带、塔里木盆地相速度不断升高,意味着青藏低速区受到南、西北、东三个方向的高速区夹击,可能导致高原中部软弱的低速物质向东南方向迁移;同时,东部地区由高速逐渐转变为大面积的低速分布,反映东部地区岩石圈较薄而软流圈发育. 随着青藏地块低速特征的减弱,印支地块北部及相邻海域、东海、东北吉林深震区、日本海、中-朝地块至蒙古东部成为120 s周期上突出的低速异常体,而上扬子地块包括四川盆地高速特征依然明显,显示出稳定的古板块特征. 南北地震带始终呈现出相对较低的速度特征,并成为划分中国大陆具有不同岩石圈相速度特征的东部与西部的天然分界.     

海潮负荷对沿海地区宽幅InSAR形变监测的影响

海岸带地区是全球自然生态环境最为复杂和脆弱的地域之一,合成孔径雷达干涉测量（InSAR）技术可以为全球人类活动、气候变暖和俯冲带剧烈构造运动等背景下的大范围海岸带地理环境变化研究提供重要观测资料.海洋潮汐导致固体地球长周期形变,波长尺度为10²~10³ km的海潮负荷引入mm级至cm级的形变梯度,此类非构造信号对海岸带InSAR精密形变分析（如：大范围、微小、缓慢且非稳态构造过程等）造成显著影响.本文以宽幅模式SAR数据为例,基于多种海潮模型研究了全球典型海岸带地区（福建、智利和阿拉斯加湾）海潮负荷效应对宽幅InSAR形变监测的影响,给出了宽幅InSAR海潮负荷三维分量估计与差分相位提取方法,并进一步讨论了基于不同海潮模型估计海潮负荷位移的差异.海潮负荷影响不仅与研究范围大小有关,其形变梯度变化与研究区域地形特征存在强相关,对于长波长形变分析而言,传统平面或者曲面拟合方法难以有效分离海潮负荷位移.     

基于全球气象再分析资料的InSAR对流层延迟改正研究

气压、温度和水汽含量等大气物理参数的时空变化导致的对流层延迟是制约合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar,InSAR)高精度应用的重要因素之一.最新研究显示气象再分析资料在补偿对流层延迟影响方面具有巨大的应用潜力,这促使我们对其有效性和鲁棒性做进一步的研究和探索.本文首先推导了利用气象再分析资料对InSAR进行对流层延迟校正的算法;然后以美国南加州地区的ENVISAT ASAR数据为例,分析了基于两种气象再分析资料(ERA-Interim和North American Regional Reanalysis,NARR)校正InSAR对流层延迟改正的效果;通过与MERIS水汽延迟改正结果比较,验证了该方法的有效性.实验结果表明:(1)不能简单忽略干延迟,可通过气象再分析资料进行有效估计;(2)通过与MERIS水汽产品获得的对流层延迟比较发现,气象再分析资料能够取得接近于MERIS的改善效果;(3)对ERA-Interim和NARR两种气象再分析资料而言,虽然后者具有更高的时间和空间分辨率,但在改正InSAR对流层延迟方面并没有表现出比前者更明显的优势;(4)气象再分析资料可以很好地估计与地形强相关的垂直分层延迟,但对于小尺度的湍流混合延迟的捕捉能力有限.综合分析认为,气象再分析资料的优势在于其数据可随时获得、免费和全球覆盖,它可以显著减弱大尺度的垂直分层延迟对干涉图相位的影响,从而有助于InSAR获取更真实可靠的地形高程和地表形变信息.     

基于BFGS法融合InSAR和GPS技术监测地表三维形变

虽然InSAR技术具有高精度、大范围和高空间分辨率等优点,但只能监测雷达视线方向上的一维地表形变;而GPS技术虽可以监测地表的三维形变,但其空间分辨率很低.本文针对融合InSAR和GPS技术监测地表高空间分辨率三维形变展开研究.首先证明了简单的局部最优化迭代算法就能求得综合InSAR和GPS监测地表形变速率的能量函数模型的全局最优估值.随后提出了利用BFGS局部最优算法反演最优的地表三维形变速率.该方法既能避免全局最优化算法计算复杂且难以收敛的问题,又能克服传统的解析法中数值计算不稳定的缺点.最后,通过模拟实验和美国南加州真实数据实验表明,该方法能够得到高精度的地表三维形变速率场.而且当观测或插值误差导致解析法误差较大时,BFGS方法仍能得到高精度、稳定的全局最优解.     

基于激光雷达手段的海南地区重力波与其波谱的季节分布特性研究

利用子午工程海南激光雷达对我国海南地区上空进行持续观测,通过3年的累积观测数据对我国低纬度地区重力波活动的季节分布特性进行研究,依据重力波线性理论对海南地区上空的大气密度扰动规律、空间功率谱及时间频率谱进行分析,并通过选择波长在1km至8km范围内具有特定波长以及具有波动周期为60 min至25min的特定频率的重力波辅助研究大气密度扰动的季节变化规律,总结得出海南地区重力波活动具有夏季大、春秋季小、而冬季依然频繁的季节性分布规律.结合海南地区特殊的地理位置与当地季节性气候特征分析得出海南地区上空重力波活动季节性变化的可能原因为青藏高原地形及我国南海地区存在的热带强对流与赤道潜流共同作用的结果.     

The European Upper Mantle as Seen by Surface Waves

We derive a global, three-dimensional tomographic model of horizontally and vertically polarized shear velocities in the upper mantle. The model is based on a recently updated global database of Love- and Rayleigh-wave fundamental-mode phase-anomaly observations, with a good global coverage and a particularly dense coverage over Europe and the Mediterranean basin (broadband stations from the Swiss and German seismic networks). The model parameterization is accordingly finer within this region than over the rest of the globe. The large-scale, global structure of our model is very well correlated with that of earlier shear-velocity tomography models, based both on body- and surface-wave observations. At the regional scale, within the region of interest, correlation is complicated by the different resolution limits associated to different databases (surface waves, compressional waves, shear waves), and, accordingly, to different models; while a certain agreement appears to exist for what concerns the grand tectonic features in the area, heterogeneities of smaller scale are less robustly determined. Our new model is only one step towards the identification of a consensus model of European/Mediterranean upper-mantle structure: on the basis of the findings discussed here, we expect that important improvements will soon result from the combination, in new tomographic inversions, of fundamental-mode phase-anomaly data like ours with observations of surface-wave overtones, of body-wave travel times, of ambient “noise”, and by accounting for an a-priori model of crustal structure more highly resolved than the one employed here.     

A comprehensive physical pattern of land-air dynamic and thermal structure on the Qinghai-Xizang Plateau

According to the boundary layer observations of three stations (Garze, Damxung and Qamdu) and relevant earth satellite, radiosonde and surface observations during the intensive observational period (IOP) of the second Tibetan (Qinghai-Xizang) Plateau Experiment of atmospheric science (TIPEX), the land-air physical process and dynamic model on the Tibetan Plateau were comprehensively analyzed in this study. The dynamic characteristics of boundary layer and the rules of turbulent motion on the plateau were illustrated. The characteristics of distributions of wind speed and direction with mutiple-layer structure and deep convective mixed layer on the plateau, the strong buoyancy effect in turbulent motion on the plateau on which the air density is obviously smaller than on the plain, and the Ekman spiral and its dynamic pump effect of the plateau deep boundary layer have been found. The local static distribution of water vapor and the horizontal advection of water vapor in the plateau boundary layer were studied. The abnomal thermodynamic structure on the plateau surface and boundary layer, including the plateau strong radiation phenomenon and strong heating source characteristics of the middle plateau, was also analyzed. The authors synthesized the above dynamic and thermodynamic structures of both surface and boundary layers on the plateau and posed the comprehensive physical model of the turbulence and convective mixture mechanism on the plateau boundary layer. The characteristics of formation, development and movement for convective cloud cluster over the plateau influencing floods in the Yangtze River area of China were studied. The conceptual model of dynamic and thermodynamic structures of turbulent motion and convective plume related to the frequent occurrence of "pop-corn-like" cloud system is given as well.