联合卫星测高、GRACE与温盐数据分析红海海平面变化季节性特征

本文利用卫星测高、GRACE与温盐数据监测2003-2014年红海海平面变化,并分析了蒸发降水以及亚丁湾-红海质量交换对红海质量变化的影响。红海地区单一的温盐数据存在覆盖不全或质量不佳的问题,综合CORA、SODA与ORAS4温盐数据估算结果得到平均比容海平面变化,以改善比容信号的精度。针对GRACE数据处理过程中截断与空间平滑滤波引起的泄漏误差,提出改进尺度因子纠正泄漏误差,利用卫星测高数据进行模拟实验验证了改进尺度因子的有效性。利用传统尺度因子和改进尺度因子反演的红海质量变化周年振幅分别为16.1±1.3 cm和20.5±1.7 cm,利用卫星测高和温盐数据估算的质量变化周年振幅为20.2±1.0 cm,表明改进尺度因子可有效减小泄漏误差的影响,改善GRACE模型反演红海质量变化的精度。卫星测高、GRACE卫星重力数据以及平均温盐数据具有较好的一致性,联合GRACE和温盐数据估算的红海综合海平面变化周年振幅为16.6±1.7 cm,与卫星测高估算的总海平面变化周年振幅(16.2±0.9 cm)基本一致,表明多源数据可构成完整的红海海平面监测手段。相比于降水-蒸发作用,红海质量变化受红海与亚丁湾的海水质量交换的影响更为显著,其主导了红海质量的季节性变化。     

Role of ocean variability and dynamic ocean topography in the recovery of the mean sea surface and gravity anomalies from satellite altimeter data

Procedures to calculate mean sea surface heights and gravity anomalies from altimeter-derived sea surface heights and along-track sea surface slopes using the least-squares collocation procedure are derived. The slope data is used when repeat track averaging is not possible to reduce ocean variability effects. Tests were carried out using Topex, Geosat, ERS-1 [35-day and 168-day (2 cycle)] data. Calculations of gravity anomalies in the Gulf Stream region were made using the sea surface height and slope data. Tests were also made correcting the sea surface heights for dynamic ocean topography calculated from a degree 360 expansion of data from the POCM-4B global ocean circulation model. Comparisons of the anomaly predictions were carried out with ship data using anomalies calculated for this paper as well as others. Received: 19 August 1996 / Accepted: 14 April 1997     

Forecasting sea level anomalies from TOPEX/Poseidon and Jason-1 satellite altimetry

This paper aims at the prediction of both global mean sea level anomalies (SLAs) and gridded SLA data in the east equatorial Pacific obtained from TOPEX/Poseidon and Jason-1 altimetric measurements. The first prediction technique (denoted as LS) is based on the extrapolation of a polynomial-harmonic deterministic least-squares model describing a linear trend, annual and semi-annual oscillations. The second prediction method (denoted as LS + AR) is a combination of the extrapolation of a polynomial-harmonic model with the autoregressive forecast of LS residuals. In the case of forecasting global mean SLA data, both techniques allow one to compute the predictions of comparable accuracy (root mean square error for 1-month in the future is of 0.5 cm). In the case of predicting gridded SLA data, the LS + AR prediction method gains significantly better prediction accuracy than the accuracy obtained by the LS technique during El Niño 1997/1998, La Niña 1998/1999 and during normal conditions.     

Seasonal sea level change from TOPEX/Poseidon observation and thermal contribution

Seasonal steric sea-level change due to temperature variation in the mixing layer is assessed using space-measured sea-surface temperature data and historical in situ temperature measurements. The results are compared with TOPEX/Poseidon satellite altimeter measurement at different large spatial scales. It is indicated that thermal effect accounts for much of the observed seasonal variability, especially when averaging over zonal regions. Some regional seasonal patterns of sea-level anomalies in the tropical oceans are well represented by the thermal model prediction. Systematic differences are shown between TOPEX/Poseidon observation and thermal contribution at a 1–2 cm level. The potential causes for these differences are discussed, including water mass exchanges among the atmosphere, land, and oceans, and error sources in the steric result and geophysical corrections applied in TOPEX/Poseidon data. Received: 25 September 1998 / Accepted: 13 July 1999     

A spectral optimally interpolated mean sea surface from satellite radar altimetry

 A technique is presented for the development of a high-precision and high-resolution mean sea surface model utilising radar altimetric sea surface heights extracted from the geodetic phase of the European Space Agency (ESA) ERS-1 mission. The methodology uses a cubic-spline fit of dual ERS-1 and TOPEX crossovers for the minimisation of radial orbit error. Fourier domain processing techniques are used for spectral optimal interpolation of the mean sea surface in order to reduce residual errors within the initial model. The EGM96 gravity field and sea surface topography models are used as reference fields as part of the determination of spectral components required for the optimal interpolation algorithm. A comparison between the final model and 10 cycles of TOPEX sea surface heights shows differences of between 12.3 and 13.8 cm root mean square (RMS). An un-optimally interpolated surface comparison with TOPEX data gave differences of between 15.7 and 16.2 cm RMS. The methodology results in an approximately 10-cm improvement in accuracy. Further improvement will be attained with the inclusion of stacked altimetry from both current and future missions. Received: 22 December 1999 / Accepted: 6 November 2000     

GNSS多径反射探测海平面变化初探

随着GNSS研究与应用的不断深入,GNSS多径反射技术用于地表环境监测已成为一种新兴的遥感手段。针对目前国内对GNSS多径反射技术开展海平面变化探测研究较少的现状,该文通过分析信噪比值的变化特性,详细给出了基于信噪比观测值的GNSS多径反射技术探测海平面变化的基本原理。利用布设在美国华盛顿州Friday Harbor岸边的GPS连续运行跟踪站SC02观测数据反演了海平面变化,并与该站相距359m的Friday Harbor验潮站的监测数据进行了对比分析,误差均值为0.091m,相关系数为0.99。实验结果进一步验证了基于大地测量型GPS接收机用于海平面变化监测的有效性。     

GRACE卫星的冰川消融与水储量变化探讨

针对青藏高原地区的冰川消融会影响发源于该地区的长江流域水储量变化这一情况,该文对二者的相关性进行深入研究。基于GRACE卫星2004—2010年重力数据,联合水文模型,利用水量平衡法,计算得到青藏高原地区冰川消融的平均速率为0.369cm/a,表明青藏高原冰川正在以较大速率消融;进而反演建立青藏高原和长江流域两个区域质量变迁的时间序列,计算得到两个区域的平均变化速率分别为-0.095cm/a和0.381cm/a,表明青藏高原地区的质量在以缓慢的速率减少,而长江流域水储量具有逐年增加的趋势且变化速率更大;最后通过分析计算结果,得出青藏高原地区冰川消融对长江流域水储量变化的贡献占其整个水量变化的20%。     

联合卫星重力、卫星测高和海洋资料研究全球海平面变化

利用GRACE、卫星测高和海洋实测温盐数据,探讨了2003~2012年间全球海平面、比容海平面和海水质量等的变化特征,并讨论了南极冰盖和格陵兰冰盖消融对全球海平面变化的影响。全球海平面整体呈上升趋势,上升速度为2.72±0.07 mm/a,且存在显著的空间分布特征。全球海平面、比容海平面和海水质量等的变化还具有显著的季节性特征,其中全球海平面变化的年周期振幅为4.6±0.3 mm。使用经验正交函数分析(EOF)得到全球海平面和比容海平面的季节性变化在南北半球存在显著的差异,但海水质量季节性变化不存在这种差异。南极冰盖和格陵兰冰盖的消融速率分别为-75.7±12.3 Gt/a和-124.1±2.9 Gt/a,对海平面的长期趋势项贡献分别为0.21±0.03 mm/a和0.34±0.01 mm/a,仅占全球海水质量增加速度1.80±0.10 mm/a的12%和19%,总计占31%,因此,两极冰盖质量消融并不是2003-2012年间海水质量增加的最主要因素。     

卫星测高推求的中国近海与全球海平面变化的趋势及小波相关分析

基于多种卫星测高资料融合海面高数据对黄海、东海、南海和全球海平面异常在1993年至2005年期间的变化趋势进行分析,结果表明:黄海、东海、南海线性趋势分别为5.05mm/a,4.8mm/a,4.27mm/a,均大于全球平均趋势项2.66mm/a,各海域的周期振幅和相位变化也存在较大差异。近海各海域与全球的海平面异常的小波交叉谱和相关谱分析表明,近海各海域与全球海平面异常变化在年际周期上具有较高的谱能,但相位差异较大,短周期相关性不明显。     

联合Argo浮标、卫星测高和GRACE数据研究海平面变化

卫星测高、GRACE、Argo等数据为监测海平面变化提供了丰富的观测数据,利用Argo数据计算的比容海平面变化,可以更加深入地理解卫星测高以及卫星重力获得的海平面变化。利用2004年1月至2010年12月间Argo浮标采集的温度和盐度数据,通过数值积分方法计算了65°S~65°N间的比容海平面异常,并通过最小二乘拟合得到比容海平面变化的长期趋势为0.63±0.45 mm/a,与Llovel得到的结果吻合较好。利用卫星测高数据得到该时间段内海平面变化趋势为2.52±0.71 mm/a,GRACE反演得到的海水质量变化引起的海平面趋势为1.84±0.13mm/a,结果表明海水质量变化成为引起海平面变化的主要因素。最后对联合卫星测高、GRACE得到比容海平面变化与相应Argo浮标数据计算结果的空间分布特征进行了比较。     

Sea level change in Hong Kong from tide gauge measurements of 1954–1999

 Tide gauge records of Hong Kong covering the past 45 years (1954.0–1999.0) are adopted to analyze the basic features of sea level changes in the region. Data sets of atmospheric pressure, southern oscillation index and sea surface temperature during the same time span are also used to determine the possible link between the sea level changes in Hong Kong and local and global geophysical processes. Results indicate that the sea level of Hong Kong has a rising trend of 1.9 ± 0.4 mm per year, and that there is an upward offset of about 15 cm in the pre-1957.0 tide gauge records. The effect of local atmospheric pressure variations on the amplitude of the annual sea level change is about 30% of the amplitude that is calculated after the effect is corrected. It is also found that the interannual variations in the sea level of Hong Kong are related to El Ni?o and La Ni?a events that happen frequently in the tropical Pacific. Received: 27 October 1999 / Accepted: 15 August 2000     

利用GRACE数据检测日本Mw9.0地震同震和震后重力变化

基于德国地学中心（GFZ）发布的GRACE RL05月重力场模型数据,考虑全球陆地同化系统陆地水储量的影响,采用300 km的扇形滤波,利用叠积法提取了日本Mw9.0地震的同震和震后重力时变信号,并利用最小二乘拟合的方法计算了两个同震重力变化极值点在地震前后85个月的重力年变率。结果表明：土壤水分和雪水引起的重力变化为-0.34~0.09 μgal;利用GRACE数据检测到的同震重力变化为-5.3~4.2 μgal,与基于PSGRN/PSCMP模型计算的结果在空间分布和量级上具有较好的一致性;震后5 a震中附近区域的重力整体上呈现增加趋势,断层上盘所在的日本海域与下盘所在的太平洋区域重力增加的最大值分别为2.6、4.5 μgal,下盘重力增加较大,可能与断层下盘所在地层的黏滞性相对较低有关。     

High precision vertical deflection over China marginal sea and global sea derived from multi-satellite altimeter

On the basis of gravity field model (EIGEN_CG01C), together with multi-altimeter data, the improved deflection of the vertical gridded in 2′×2′ in China marginal sea and gridded in 5′×5′ in the global sea was determined by using the weighted method of along-track least squares, and the accuracy is better than 1.2″ in China marginal sea. As for the quality of the deflection of the vertical, it meets the challenge for the gravity field of high resolution and accuracy. It shows that, compared with the shipboard gravimetry in the sea, the accuracy of the gravity anomalies computed with the marine deflection of the vertical by inverse Vening-Meinesz formula is 7.75 m·s?2.     

联合多种测高数据确定中国边缘海及全球海域的垂线偏差

On the basis of gravity field model （EIGEN_CG01C）, together with multi-altimeter data, the improved deflection of the vertical gridded in 2＇×2＇ in China marginal sea and gridded in 5＇×5＇ in the global sea was determined by using the weighted method of along-track least squares, and the accuracy is better than 1.2^# in China marginal sea. As for the quality of the deflection of the vertical, it meets the challenge for the gravity field of high resolution and accuracy, it shows that, compared with the shipboard gravimetry in the sea, the accuracy of the gravity anomalies computed with the marine deflection of the vertical by inverse Vening-Meinesz formula is 7.75 m.s ^-2.     

Combination of temporal gravity variations resulting from superconducting gravimeter (SG) recordings, GRACE satellite observations and global hydrology models

Gravity recovery and climate experiment (GRACE)-derived temporal gravity variations can be resolved within the μgal (10^?8 m/s ²) range, if we restrict the spatial resolution to a half-wavelength of about 1,500 km and the temporal resolution to 1 month. For independent validations, a comparison with ground gravity measurements is of fundamental interest. For this purpose, data from selected superconducting gravimeter (SG) stations forming the Global Geodynamics Project (GGP) network are used. For comparison, GRACE and SG data sets are reduced for the same known gravity effects due to Earth and ocean tides, pole tide and atmosphere. In contrast to GRACE, the SG also measures gravity changes due to load-induced height variations, whereas the satellite-derived models do not contain this effect. For a solid spherical harmonic decomposition of the gravity field, this load effect can be modelled using degree-dependent load Love numbers, and this effect is added to the satellite-derived models. After reduction of the known gravity effects from both data sets, the remaining part can mainly be assumed to represent mass changes in terrestrial water storage. Therefore, gravity variations derived from global hydrological models are applied to verify the SG and GRACE results. Conversely, the hydrology models can be checked by gravity variations determined from GRACE and SG observations. Such a comparison shows quite a good agreement between gravity variation derived from SG, GRACE and hydrology models, which lie within their estimated error limits for most of the studied SG locations. It is shown that the SG gravity variations (point measurements) are representative for a large area within the accuracy, if local gravity effects are removed. The individual discrepancies between SG, GRACE and hydrology models may give hints for further investigations of each data series.     

GRACE时变重力场位系数相关误差的滤波消除技术

利用UTCSR发布的2003-01~2013-07GRACE RL05月平均重力场模型,分析比较了高斯滤波、各向异性滤波、扇形滤波和维纳滤波,并结合去相关滤波在反演南极地区冰盖质量变化方面的差异。通过计算得到以下结论：①基于121组重力场模型阶方差分布得到维纳滤波与半径为300 km的高斯滤波效果最为接近,说明300 km滤波半径完全可以满足质量变化信号的提取;②在一定范围内,提高滤波半径能提高反演结果信噪比,建议南极区域的滤波半径为500 km;采用相同滤波半径,不同空间滤波算法计算的质量变化率基本一致,在南极区域可以选取任一滤波方法;③与其他算法相比,去相关滤波算法能在一定程度消除球谐系数中存在的系统误差,改善反演结果。     

根据人卫激光测距、GRACE和地球物理模型求解地球低阶重力场季节变化

本文利用2002年4月至2010年10月的Lageos1和Lageos2两颗激光卫星观测数据、GRACE以及地球物理模型三种独立的方法计算地球低阶重力场系数J 的变化,根据大气压强数据计算 J 时分别按反变气压计（IB）和非反变气压计(NIB)两种假设进行计算。通过分析 J 的季节特性表明,大气在NIB假设下得到的周年振幅比在IB假设下得到的振幅大3倍左右,相位相差47°;大气和陆地水的质量变化对 J 周年变化的贡献占主导地位,海洋的影响最小;大气、海洋和陆地水得到 J 半年振幅和相位值与SLR得到的振幅和相位值吻合较差,尤其是在IB假设下大气得到的结果与SLR结果相差最大; SLR、GRACE和地球物理模型三种独立方法得到的 J 周年项之间吻合相对较好,GRACE得到的周年振幅比SLR得到的周年振幅大50%左右, SLR观测得到的 J 周年振幅介于在NIB和IB两种假设下地球物理模型得到的结果之间;GRACE与SLR得到的 J 半年项的振幅相同,在IB假设下AOW得到的 J 半年振幅和相位与SLR结果差异最大。     

像素级和对象级高分辨率遥感图像变化检测分析

高分辨率遥感图像具有丰富的纹理信息,而像素级变化检测方法主要分析图像的光谱信息,导致将像素级变化检测方法用于高分辨率遥感图像具有一定的局限性.因此,本文提出了一种像素级与对象级相结合的高分辨率遥感图像变化检测方法,解决了像素级与对象级变化检测方法中存在的椒盐现象、误检等问题.首先,结合高分辨率遥感图像的多维特征,构建遥...     

Glacier and glacial lake classification for change detection studies using satellite data: a case study from Baspa basin,western Himalaya

Using high-resolution Google Earth^TM images in conjunction with Landsat images, the glaciers and lakes in the Baspa basin are classified to explore the recent changes. A total number of 109 glaciers (187 ± 3.7 km²) are mapped and subsequently classified as compound valley glaciers, simple valley glaciers, cirques, niches, glacieretes and ice aprons. The compound and simple valley glaciers contribute 67.1 ± 1.3% and 19.8 ± 0.3% to the total glacier cover of the basin. Similarly, a total number of 129 glacial lakes (0.360 ± 0.007 km²) are identified. From 1976 to 2011, the compound valley glaciers have lost a small area of 10.3 ± 0.03% at a rate of 0.41 ± 0.002 km² a^-1, whereas the niche glaciers have lost higher area of 40.1 ± 0.001% at a rate of 0.04 ± 0.0001 km² a^-1. Change detection of two benchmark glacial lakes revealed a progressive expansion during recent decades. The Baspa Bamak proglacial lake has expanded from 0.020 ± 0.0004 km² (2000) to 0.069 ± 0.001 km² (2011). Due to the complete loss of source ice, another glacial lake has expanded from 0.09 ± 0.001 km² (1994) to 0.10 ± 0.002 km² (2011). During the study period, the mean annual temperature that is T_avg, T_min and T_max have increased significantly at the 95% confidence level by 1.5 ^oC (0.070 °C a^-1), 1.8 ^oC (0.076 °C a^-1) and 1.6 ^oC (0.0071 °C a^-1) from 1985 to 2008. However, the precipitation has decreased significantly from 1976 and 1985 to 2008.