The equivalence of predictions from universal kriging and intrinsic random-function kriging

A proof is provided that the predictions obtained from kriging based on intrinsic random functions of orderk are identical to those obtained from anappropriate universal kriging model. This is a theoretical result based on known variability measures. It does not imply that people performing traditional universal kriging will get the same predictions as those using intrinsic random functions, because traditionally these methods differ in how variability is modeled. For intrinsic random functions, the same proof shows that predictions do not depend on the specific choice of the generalized covariance function. It is argued that the choice between these methods is really one of modeling and estimating the variability in the data.     

Estimation of a multidimensional covariance function in case of anisotropy

A model of a multivariate covariance function with an ellipsoidal directional correlation scale has been developed. The axes of the ellipsoidal scale are related to the eigenvalues and eigenvectors of a matrix B which characterizes the ellipsoid of the range of influence. The matrix B is found to be related to a matrix T which can be estimated directly from sparse sampling data and can be used to determine estimates of the matrix B. The method has been applied to both two-dimensional and three-dimensional cases. The numerical results show that the satisfactory accuracy is obtained with sparse sampling data from an anisotropic random function.     

大地测量相关观测抗差估计理论

相关观测异常诊断、质量控制是测量数据处理领域亟待解决的难题之一。分别从方差膨胀模型和相关权元素压缩模型入手研究了相关观测的质量控制理论和方法；给出了误差影响函数；构造了方差膨胀函数和权因子收缩函数；利用观测量的等价协方差阵和等价权矩阵讨论了相关观测质量控制的计算方法。该等价协方差矩阵和等价权矩阵不仅保持了原有协方差矩阵和权矩阵的对称性，而且保持了原有协方差矩阵的相关性不变。计算结果表明异常观测的方差膨胀法和等价权法能有效地控制异常观测对参数估值的影响。     

相关稳健估计及其在测量数据

在对现有的基于等价权的稳健估计的弊端进行分析的基础上,直接从方差膨胀模型出发,建立起一套基于调整观测量方差-协方差阵的相关稳健估计方法,并通过算例证明其估计思想、模型的正确性和可行性.     

海洋风暴形成的一种动力学机制

文中从观测统计学、瞬变涡动能量学和 MM5中尺度数值模拟角度 ,研究了海洋风暴 (爆发性气旋 )形成的气候特征及其可能的动力学机制 ,揭示了一幅爆发性发展的物理图像。结果表明 ,在冷季大气特别是日本以东洋面上大气特有的热力气候背景下 ,通过同海洋风暴过程相联系的涡动热通量 vθ的向极地输送 (- vθ· θm>0 ) ,将季节尺度的时间平均有效位能向瞬变涡旋时间尺度的涡动有效位能转换 ,是海洋风暴形成的主要动力机制。在该过程中转换来的具有最大贡献的涡动有效位能 ,连同具有次大贡献的积云加热制造的涡动有效位能(q3 )一起 ,通过暖异常区 (α >0 )暖湿空气上升运动 (-ω >0 )的斜压转换 (-ωα) ,促使涡动动能增长。同时 ,补充的涡动有效位能又加强了暖异常区的暖湿空气上升运动 ,进而产生积云对流活动及其潜热释放的正反馈过程 ,最终导致涡动动能急剧增长和海洋风暴的形成。海-气潜热输送的作用是在风暴形成初期提供后来积云尺度对流活动及潜热释放的水汽潜力。研究还表明 ,海洋风暴主要发生在冷季月份 1 3 0°E以东的中高纬洋面上 ,这种对特定季节和特定海域的依赖性是大气和海洋气候背景的动力 /热力共同作用的结果     

变化涡动沾性系数的Ekman动量近似模式风场结构

在边界层动力学中，涡动粘性系数是影响边界层风场结构的一个重要参数。本文利用边界层动力学中的Ekman动量近似理论，给出了涡动粘性系数随高度缓变条件下的Ekman动量近似边界层模式解，着重讨论了边界层的风场结构、水平散度、垂直涡度以及边界层顶部的垂直速度。结果分析表明：与常值涡动粘性系数情况相比，在边界层低层随高度增加的涡动粘性系数可以导致低层边界层风速随高度迅速增加，即风速垂直切变增加，同时风速矢与地转风之间的夹角减小。惯性项作用可以导致上述作用在气旋性区域减小、而在反气旋性区域增大。随高度增加的涡动粘性系数导致水平散度绝对值、垂直涡度绝对值以及边界层顶部的垂直速度绝对值在气旋性区域减小，而在反气性旋区域增大。涡动粘性系数与惯性之间的非线性相互作用是边界层动力学中重要过程。     

A SIMPLIFIED APPROACH TO MODELING 3D SEDIMENT-LADEN TURBULENT FLOWS

ＡＳＩＭＰＬＩＦＩＥＤＡＰＰＲＯＡＣＨＴＯＭＯＤＥＬＩＮＧ３ＤＳＥＤＩＭＥＮＴ－ＬＡＤＥＮＴＵＲＢＵＬＥＮＴＦＬＯＷＳＤｏｎｇｈｕｏＺＨＯＵａｎｄＳａｍＳ．Ｙ．ＷＡＮＧＡｂｓｔｒａｃｔ：Ａ３－ｄｎｕｍｅｒｉｃａｌｍｏｄｅｌｆｏｒｓｉｍｕｌａｔｉｎｇｓ...     

花岗岩的成岩作用和形成环境——存在的问题和解决的途径

本文针对当前花岗岩研究中存在的一些问题,结合笔者等近年来的工作,概述了进一步研究花岗岩成岩作用和形成环境的新思路。主要包括:(1)利用成分协变图,结合岩石包体成因类型的认识,探讨花岗岩化学演化进程;(2)在全面考察花岗岩的构造样式、岩石组合、岩浆演化和区域地壳不均一性基础上,分析花岗岩形成的构造背景。     

Generalized cross-validation for covariance model selection

A weighted cross-validation technique known in the spline literature as generalized cross-validation (GCV), is proposed for covariance model selection and parameter estimation. Weights for prediction errors are selected to give more importance to a cluster of points than isolated points. Clustered points are estimated better by their neighbors and are more sensitive to model parameters. This rational weighting scheme also provides a simplifying significantly the computation of the cross-validation mean square error of prediction. With small- to medium-size datasets, GCV is performed in a global neighborhood. Optimization of usual isotropic models requires only a small number of matrix inversions. A small dataset and a simulation are used to compare performances of GCV to ordinary cross-validation (OCV) and least-squares filling (LS).     

Study on the processing method of nighttime CO_2 eddy covariance flux data in ChinaFLUX

At present, using Eddy Covariance (EC) method to estimate the "true value" of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO2 flux data) appearing in post-processing CO2 flux data. The ratio of effective and reliable nighttime EC CO2 flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO2 flux data and its effect analysis on estimating CO2 flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO2 flux using EC method, and introduce the general theory and method for processing nighttime CO2 flux data. By analyzing the relationship between nighttime CO2 flux and air fraction velocity u., we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u.c) for screening the effective nighttime CO2 flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO2 flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.