基于方差递推法的Kalman滤波在钟差预报中的应用

卫星导航系统中星载原子钟的钟差预报对于导航、定位及授时具有重要的作用。为了提高卫星钟差预报的精度,设计了一种两步确定卫星钟噪声协方差矩阵的Kalman滤波钟差预报模型。该方法首先基于Hadamard总方差确定卫星钟噪声协方差矩阵的初值,然后,使用方差递推法得到滤波过程中卫星钟的噪声协方差矩阵。使用GPS系统的星载铷钟数据进行短期预报,并与常用的二次多项式模型、灰色模型进行对比,结果表明:本文中提出的方法可以实现高精度的卫星钟差预报且预报效果优于两种常用模型,同时,该方法能够在一定程度上弥补预报误差随预报时间增加而不断变大的不足。     

局部区域模型重力异常快速算法

地球重力场位系数模型可以用于计算局部重力扰动场元。然而随着地球重力场模型阶次的提高、局域重力场计算范围的增大,其计算速度往往不能满足工程需求。针对这一问题,在对位系数模型泰勒级数展开的基础上提出了采用向量运算、混合编程的方法,同时对连带勒让德函数Belikov递推方法中与经纬度无关的量进行了预先计算,有效提高了计算速度。提出的方法对于利用超高阶次重力场模型快速解算大范围、高分辨率重力场元数据以及累加求和计算具有一定的参考与借鉴意义。     

利用Abel-Poisson径向基函数模型化局部重力场

多种类型高分辨率重力场数据的不断增加,使得在局部范围内精化重力场模型成为了可能。本文采用Abel-Poisson核将重力场量表示成有限个径向基函数线性求和的形式,对局部区域的多种重力场数据进行联合建模。为了提高运算速度,运用了基于自适应精化格网算法的最小均方根误差准则(RMS)来求解径向基函数平均带宽。以南海核心地区为例,联合两种不同类型、不同分辨率的重力场资料(大地水准面起伏6'×6'、重力异常2'×2'),构建了局部区域高分辨率的重力场模型。所建模型表示的重力场参量达到了2'×2'的分辨率,对原始的重力异常数据(2'×2')拟合的符合程度达到±0.8×10-5m/s2。结果表明,利用径向基函数方法进行局部重力场建模,避免了球谐函数建模收敛慢的问题,有效提高了模型表示重力场的分辨率。     

顾及周期项误差和起点偏差修正的北斗卫星钟差预报

根据北斗卫星导航系统星载原子钟自身的物理特性,采用武汉大学IGS数据中心发布的2016年1月1日至2016年11月1日共306天的事后钟差产品进行谱分析。分析结果表明:北斗卫星导航系统的3类卫星钟都存在一定的周期特性;其中GEO和IGSO卫星钟的主周期相对较为明显;GEO卫星钟的主周期依次为12、24、8和6h;IGSO的主周期为24、12、8和6h;而MEO的3个主周期为12.9、6.4和24h。依据各类原子钟的周期特性,同时对各天之间钟差的起始点偏差进行修正,并利用修正模型对北斗卫星钟差进行预报和精度分析。试验结果表明,改进的预报模型能显著提升北斗卫星的钟差预报精度,北斗卫星钟差24h、12h、6h的平均预报精度分别能达到6.55ns、3.17ns和1.76ns。     

Trends in Extreme Mean Sea Level Quantiles from Satellite Altimetry

Satellite altimetry allows the study of sea-level long-term variability on a global and spatially uniform basis. Here quantile regression is applied to derive robust median regression trends of mean sea level as well as trends in extreme quantiles from radar altimetry time series. In contrast with ordinary least squares regression, which only provides an estimate on the rate of change of the mean of data distribution, quantile regression allows the estimation of trends at different quantiles of the data distribution, yielding a more complete picture of long-term variability. Trends derived from basin-wide averaged regional mean sea level time series are robust and similar for all quantiles, indicating that all parts of the data distribution are changing at the same rate. In contrast, trends are not robust and diverge across quantiles in the case of local time series. Trends are under- (over-)estimated in the western (eastern) equatorial Pacific. Furthermore, trends in the lowermost quantile (0.05) are larger than the median trend in the western Pacific, while trends in the uppermost quantile (0.95) are lower than the median trend in the eastern Pacific. These differences in trends in extreme mean sea level quantiles are explained by the exceptional effect of the strong 1997–1998 El Niño–Southern Oscillation (ENSO) event.     

强潮流作用下桥墩不对称“双肾型”冲刷地貌特征与机理

本文在海图地形资料分析桥轴线附近的海床自然冲刷的基础上,利用多波束测深技术研究大桥主墩附近局部冲刷地形。结果表明,该大桥桥位附近地形冲刷较显著,且大桥主墩位置有持续冲刷的趋势;主墩上、下游群桩最大冲刷深度呈上游最深、中部淤积、下游渐深的不对称形态,最大局部冲刷深度为4 m;桥墩整体冲刷坑形态呈南北“双肾型”;潮流流向与桥墩迎流面存在偏南的入射角,使得各桥墩南侧的最大冲深和冲刷范围均大于北侧。     

Quality control methods for KOOS operational sea surface temperature products

Sea surface temperature SST obtained from the initial version of the Korea Operational Oceanographic System(KOOS) SST satellite have low accuracy during summer and daytime. This is attributed to the diurnal warming effect. Error estimation of SST data must be carried out to use the real-time forecasting numerical model of the KOOS. This study suggests two quality control methods for the KOOS SST system. To minimize the diurnal warming effect, SSTs of areas where wind speed is higher than 5 m/s were used. Depending on the wind threshold value, KOOS SST data for August 2014 were reduced by 0.15°C. Errors in SST data are considered to be a combination of random, sampling, and bias errors. To estimate bias error, the standard deviation of bias between KOOS SSTs and climatology SSTs were used. KOOS SST data yielded an analysis error standard deviation value similar to OSTIA and NOAA NCDC(OISST) data. The KOOS SST shows lower random and sampling errors with increasing number of observations using six satellite datasets. In further studies, the proposed quality control methods for the KOOS SST system will be applied through more long-term case studies and comparisons with other SST systems.     

Bias estimation and correction for triangle-based surface area calculations

The calculation of surface area is meaningful for a variety of space-filling phenomena, e.g., the packing of plants or animals within an area of land. With Digital Elevation Model (DEM) data we can calculate the surface area by using a continuous surface model, such as by the Triangulated Irregular Network (TIN). However, just as the triangle-based surface area discussed in this paper, the surface area is generally biased because it is a nonlinear mapping about the DEM data which contain measurement errors. To reduce the bias in the surface area, we propose a second-order bias correction by applying nonlinear error propagation to the triangle-based surface area. This process reveals that the random errors in the DEM data result in a bias in the triangle-based surface area while the systematic errors in the DEM data can be reduced by using the height differences. The bias is theoretically given by a probability integral which can be approximated by numerical approaches including the numerical integral and the Monte Carlo method; but these approaches need a theoretical distribution assumption about the DEM measurement errors, and have a very high computational cost. In most cases, we only have variance information on the measurement errors; thus, a bias estimation based on nonlinear error propagation is proposed. Based on the second-order bias estimation proposed, the variance of the surface area can be improved immediately by removing the bias from the original variance estimation. The main results are verified by the Monte Carlo method and by the numerical integral. They show that an unbiased surface area can be obtained by removing the proposed bias estimation from the triangle-based surface area originally calculated from the DEM data.     

中国地区集合预报产品自适应递减平均偏差订正法的改进研究

基于中国地区T213集合预报产品2 m温度预报数据,采用卡尔曼滤波类型的自适应递减平均法进行偏差订正处理,原方案在剧烈降温天气订正效果表现不理想。通过对递减平均参数w的重新构建得到改进的订正方案w(i,p)(i为站点信息,p为天气过程信息),在此基础上进一步优化对历史信息的有效提取,得到改进的方案w(i,p)相似法和w(i,p)统计法,并进行效果检验。结果表明:改进为包含空间和天气过程信息的函数w(i,p)后方案的订正效果得到不同程度的提高,其中24 h剧烈降温预报各成员预报均方根误差平均减小了0. 15℃;而进一步改进的w(i,p)统计法在当前几种剧烈降温预报中订正效果最优,其集合平均偏差与w(i,p)方案相比减小2. 54℃。     

Latitudinal and Scan-dependent Biases of Microwave Humidity Sounder Measurements and Their Dependences on Cloud Ice Water Path

The relationship between differences in microwave humidity sounder(MHS)–channel biases which represent measured brightness temperatures and model-simulated brightness temperatures, and cloud ice water path(IWP) as well as the influence of the cloud liquid water path(LWP) on the relationship is examined. Seven years(2011–17) of NOAA-18 MHS-derived measured brightness temperatures and IWP/LWP data generated by the NOAA Comprehensive Large Array-data Stewardship System Microwave Surface and Precipitation Products System are used. The Community Radiative Transfer Model, version2.2.4, is used to simulate model-simulated brightness temperatures using European Center for Medium-Range Weather Forecasts reanalysis data as background fields. Scan-angle deviations of the MHS window channel biases range from-1.7 K to1.0 K. The relationships between channels 2, 4, and 5 biases and scan angle are symmetrical about the nadir. The latitudedependent deviations of MHS window channel biases are positive and range from 0–7 K. For MHS non-window channels,the latitudinal deviations between measured brightness temperatures and model-simulated brightness temperatures are larger when the detection height is higher. No systematic warm or cold deviations are found in the global spatial distribution of difference between measured brightness temperatures and model-simulated brightness temperatures over oceans after removing scan-angle and latitudinal deviations. The corrected biases of five different MHS channels decrease differently with respect to the increase in IWP. This decrease is stronger when LWP values are higher.