Bayesian area-to-point kriging using expert knowledge as informative priors

Area-to-point (ATP) kriging is a common geostatistical framework to address the problem of spatial disaggregation or downscaling from block support observations (BSO) to point support (PoS) predictions for continuous variables. This approach requires that the PoS variogram is known. Without PoS observations, the parameters of the PoS variogram cannot be deterministically estimated from BSO, and as a result, the PoS variogram parameters are uncertain. In this research, we used Bayesian ATP conditional simulation to estimate the PoS variogram parameters from expert knowledge and BSO, and quantify uncertainty of the PoS variogram parameters and disaggregation outcomes. We first clarified that the nugget parameter of the PoS variogram cannot be estimated from only BSO. Next, we used statistical expert elicitation techniques to elicit the PoS variogram parameters from expert knowledge. These were used as informative priors in a Bayesian inference of the PoS variogram from BSO and implemented using a Markov chain Monte Carlo algorithm. ATP conditional simulation was done to obtain stochastic simulations at point support. MODIS (Moderate Resolution Imaging Spectroradiometer) atmospheric temperature profile data were used in an illustrative example. The outcomes from the Bayesian ATP inference for the Matérn variogram model parameters confirmed that the posterior distribution of the nugget parameter was effectively the same as its prior distribution; for the other parameters, the uncertainty was substantially decreased when BSO were introduced to the Bayesian ATP estimator. This confirmed that expert knowledge brought new information to infer the nugget effect at PoS while BSO only brought new information to infer the other parameters. Bayesian ATP conditional simulations provided a satisfactory way to quantify parameters and model uncertainty propagation through spatial disaggregation.     

非线性半参数模型最小二乘核估计法方程系数阵的非奇异性质

本文基于非线性半参数模型最小二乘核估计的迭代解法,证明了非线性半参数模型最小二乘类估计法方程系数阵在一定条件下的非奇异性。这种性质可类推到非线性半参数模型的其他最小二乘类估计方法和其他非线性解算方法,这是对非线性半参数模型估计理论的初步推导。     

A theoretical study on terrestrial gravimetric data refinement by Earth gravity models

The idea of this paper is to present estimators for combining terrestrial gravity data with Earth gravity models and produce a high‐quality source of the Earth's gravity field data through all wavelengths. To do so, integral and point‐wise estimators are mathematically developed, based on the spectral combination theory, in such a way that they combine terrestrial data with one and/or two Earth gravity models. The integral estimators are developed so that they become biased or unbiased to a priori information. For testing the quality of the estimators, their global mean square errors are generated using an Earth gravity model08 model and one of the recent products of the gravity field and steady‐state ocean circulation explorer mission. Numerical results show that the integral estimators have smaller global root mean square errors than the point‐wise ones but they are not efficient practically. The integral estimator of the biased type is the most suited due to its smallest global root mean square error comparing to the rest of the estimators. Due largely to the omission errors of Earth gravity models the point‐wise estimators are not sensitive to the Earth gravity model commission error; therefore, the use of high‐degree Earth gravity models is very influential for reduction of their root mean square errors. Also it is shown that the use of the ocean circulation explorer Earth gravity model does not significantly reduce the root mean square errors of the presented estimators in the presence of Earth gravity model08. All estimators are applied in the region of Fennoscandia and a cap size of 2° for numerical integration and a maximum degree of 2500 for generation of band‐limited kernels are found suitable for the integral estimators.     

采用变系数回归模型提取GPS坐标序列季节性信号

提出采用变系数回归模型提取连续GPS坐标序列中包含的振幅时变季节性信号。对模拟数据及实际GPS坐标序列两组数据的分析结果表明,变系数回归模型在提取GPS坐标序列季节信号方面比传统模型更有效。经处理,GPS坐标序列能获得更合理的速度及噪声估计。     

测流不确定度的插值方差估计法

目前,测流不确定度通过误差试验或通过经验数值来确定,但这些方式存在着工作量大或不确定估计不足等局限性。为解决此问题,对基于实测数据和统计理论的插值方差估计法在不同测流条件下进行了验证,选取白河、襄阳和沙洋3个流量站进行了实测数据的不确定度分析,同时对白河站进行了Monte Carlo试验,比较插值方差估计法得到的不确定度与真实误差的差异。结果表明,插值方差估计法能较好地反映水位变化的影响,插值方差估计法所得到的不确定度与真实测流误差的相关系数达0.64,与断面水位变化的Spearman相关系数达0.79,高、中水位情况下插值方差估计法的不确定度估计结果较为合理,低水位情况下偏高。     

Estimation of rainfall elasticity of streamflow in Australia

Abstract

Estimates of rainfall elasticity of streamflow in 219 catchments across Australia are presented. The rainfall elasticity of streamflow is defined here as the proportional change in mean annual streamflow divided by the proportional change in mean annual rainfall. The elasticity is therefore a simple estimate of the sensitivity of long-term streamflow to changes in long-term rainfall, and is particularly useful as an initial estimate of climate change impact in land and water resources projects. The rainfall elasticity of streamflow is estimated here using a hydrological modelling approach and a nonparametric estimator. The results indicate that the rainfall elasticity of streamflow (? _P ) in Australia is about 2.0–3.5 (observed in about 70% of the catchments), that is, a 1% change in mean annual rainfall results in a 2.0–3.5% change in mean annual streamflow. The rainfall elasticity of streamflow is strongly correlated to runoff coefficient and mean annual rainfall and streamflow, where streamflow is more sensitive to rainfall in drier catchments, and those with low runoff coefficients. There is a clear relation-ship between the ? _P values estimated using the hydrological modelling approach and those estimated using the nonparametric estimator for the 219 catchments, although the values estimated by the hydrological modelling approach are, on average, slightly higher. The modelling approach is useful where a detailed study is required and where there are sufficient data to reliably develop and calibrate a hydrological model. The nonparametric estimator is useful where consistent estimates of the sensitivity of long-term streamflow to climate are required, because it is simple to use and estimates the elasticity directly from the historical data. The nonparametric method, being model independent, can also be easily applied in comparative studies to data sets from many catchments across large regions.     

平衡损失下回归系数stein压缩估计的性质

在平衡损失风险函数准则下,研究线性模型中回归系数的stein估计优于最小二乘估计(LS)的充分必要条件,然后在pitman closeness(PC)准则下比较了stein估计相对于最小二乘估计的优良性。     

神经网络辅助的GPS/INS组合导航故障检测算法

针对GPS/INS松组合导航系统观测信息无冗余,而且观测信息可能存在故障的情形,提出一种神经网络辅助的组合导航故障检测算法。该算法克服了基于模型的故障检测算法受模型误差影响的局限性;能够自动地对观测信息进行故障的检测、定位和剔除;能够基于故障检测后可靠的观测信息进一步调整动力学模型信息对导航解的贡献;能够在GPS失锁时,较好地进行导航预报。最后利用车载实测数据进行验证,结果表明该算法能够很好地从模型误差中分离出观测信息含有的故障信息,降低了故障检测算法存在的虚警率,避免故障信息对导航解的影响;且GPS失锁时,神经网络的预报输出在一定程度上能够进一步提高导航解的精度。     

公开DEM辅助无地面控制点国产卫星影像定位方法

在全球测绘的背景下,为实现无地面控制点的卫星摄影测量,提出了利用公开DEM辅助国产卫星影像进行无控定位的方法。为充分利用公开DEM在大范围内具有一致稳定的高精度特性,将其作为基准与从影像中提取的待定位DEM进行表面匹配,匹配时采用引入截尾最小二乘估计的最小高差(LZD)法,并根据对应点高差的分布自适应探测及剔除DEM之间的差异;匹配确定的变换参数用于对直接定位结果进行物方改正。设计了多组针对天绘一号和资源三号国产卫星影像的对比试验,结果表明本方法切实可行,能充分利用基准DEM的优势,具有较好的稳健性;定位精度很大程度上取决于但不限于基准DEM的精度,基本不受其分辨率的影响;当待定位DEM分辨率较高时,单景影像利用SRTM DEM即可得到能较好满足1∶5万比例尺地形图测制要求的无控定位精度。本方法还为卫星影像定位精度的检核提供了一种新的有效手段。