Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.     

Effet du vieillissement des céramiques poreuses sur leur capacité à évaluer la concentration de pesticide en solution

The soil solution sampling by ceramic cups allows pesticide transfer monitoring in the soil during long times. The ageing of material involves a bias in the sampling results. In laboratory, the comparison of two types of ceramic suction cups, new and installed in situ during four years, shows a modification of the hydrodynamic properties and a possible evolution of the adsorption capacity of the matrix. The passage rate, as well qualitative as quantitative, is better for the old material. Recommendations about site management are finally exposed. To cite this article: N. Domange et al., C. R. Geoscience 336 (2004).     

利用现有重力场模型求定CHAMP卫星加速度计修正参数

CHAMP卫星加速度计数据的标定是通过确定其尺度因子和偏差参数来完成的.本文基于能量守恒方程,给出利用现有重力场模型标定CHAMP卫星加速度数据的基本原理和数学模型;提出相邻历元间差分算法,大大简化了观测方程,同时避免积分常量的计算.该算法既能同时解算尺度因子和偏差参数,也可任意求解其中之一.基于实测的CHAMP卫星加速度数据,利用EGM96模型和最新公布的EIGEN-2模型进行计算与比较,验证该方法的有效性.     

Using raw regional climate model outputs for quantifying climate change impacts on hydrology

General circulation model outputs are rarely used directly for quantifying climate change impacts on hydrology, due to their coarse resolution and inherent bias. Bias correction methods are usually applied to correct the statistical deviations of climate model outputs from the observed data. However, the use of bias correction methods for impact studies is often disputable, due to the lack of physical basis and the bias nonstationarity of climate model outputs. With the improvement in model resolution and reliability, it is now possible to investigate the direct use of regional climate model (RCM) outputs for impact studies. This study proposes an approach to use RCM simulations directly for quantifying the hydrological impacts of climate change over North America. With this method, a hydrological model (HSAMI) is specifically calibrated using the RCM simulations at the recent past period. The change in hydrological regimes for a future period (2041–2065) over the reference (1971–1995), simulated using bias‐corrected and nonbias‐corrected simulations, is compared using mean flow, spring high flow, and summer–autumn low flow as indicators. Three RCMs driven by three different general circulation models are used to investigate the uncertainty of hydrological simulations associated with the choice of a bias‐corrected or nonbias‐corrected RCM simulation. The results indicate that the uncertainty envelope is generally watershed and indicator dependent. It is difficult to draw a firm conclusion about whether one method is better than the other. In other words, the bias correction method could bring further uncertainty to future hydrological simulations, in addition to uncertainty related to the choice of a bias correction method. This implies that the nonbias‐corrected results should be provided to end users along with the bias‐corrected ones, along with a detailed explanation of the bias correction procedure. This information would be especially helpful to assist end users in making the most informed decisions.     

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

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

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

根据北斗卫星导航系统星载原子钟自身的物理特性,采用武汉大学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。     

基于稳定同位素的东太平洋大青鲨(Prionace glauca)营养级研究

碳、氮稳定同位素为分析淡水和海洋生态系统的生物营养关系提供了有力手段,但在大洋鱼类中的应用还很少。本研究以热带东太平洋大青鲨(Prionace glauca)整节脊椎骨作为样品(叉长范围153—242cm),运用稳定同位素技术,分析大青鲨的δ~(13)C值和δ~(15)N值及其变化,比较不同基线生物的选择对营养级计算的影响。结果表明,大青鲨δ~(13)C范围为–15.76‰—–13.41‰,最大差值2.35‰;δ~(15)N范围为10.62‰—17.72‰,最大差值7.1‰;δ~(13)C值和δ~(15)N值随个体长度的变化不明显(可能原因是样品鱼均为较大个体),性别间的差异不显著,但不同基线生物的选择对营养级计算值的影响较大。本文研究表明,通过整节脊椎骨的稳定同位素测定来获得大青鲨的摄食特征信息是可行的。由于基线生物对营养级计算的影响问题难以在短期内解决,今后可将相对营养级或营养级的变化作为重点,研究人类活动(如捕捞)对大洋鱼类摄食和营养级的影响。     

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.