When do we need a trend model in kriging?

Under usual estimation practice with local search windows for data and for interpolation situations, universal kriging and ordinary kriging yield the same estimates, using a data set with apparent trend, for both the unknown attribute and its trend component. Modeling the trend matters only in extrapolation situations. Because conditions of the case study presented arise most frequently in practice, the simpler ordinary kriging is the preferred option.     

High time resolution monitoring of tropospheric temperature with a Radio Acoustic Sounding System (RASS)

We have observed the time-height variation of the temperature field in the upper troposphere using a Radio Acoustic Sounding System (RASS) which consists of the MU radar and a high-power acoustic transmitter. The fast beam steerability of the MU radar has made it possible to measure temperature profiles in a fairly wide height range in the upper troposphere (5–11 km), even under intense wind conditions. Observations were continued for about 32 hr on 24–26 December, 1986 with a time-height resolution of 30 min and 150 m. During the observation period, the tropospheric jet was so intense that the acoustic wavefronts were severely distorted. Using wind velocity profiles observed by the MU radar we have numerically estimated the propagation of acoustic wavefronts, and further determined favorable pointing directions for the MU radar to receive significant backscattering from refractive index fluctuations produced by the acoustic waves. Conventional radiosonde soundings were carried out every 6 hr, which showed a temperature decrease of 4 K/day in the upper troposphere during the observation period. Temperature profiles taken by RASS agree well with the radiosonde results.     

Theoretical background for the inversion of seismic waveforms including elasticity and attenuation

To account for elastic and attenuating effects in the elastic wave equation, the stress-strain relationship can be defined through a general, anisotropic, causal relaxation function ^ijkl (x, ). Then, the wave equation operator is not necessarily symmetric (self-adjoint), but the reciprocity property is still satisfied. The representation theorem contains a term proportional to the history of strain. The dual problem consists of solving the wave equation withfinal time conditions and an anti-causal relaxation function. The problem of interpretation of seismic waveforms can be set as the nonlinear inverse problem of estimating the matter density (x) and all the functions ^ijkl (x, ). This inverse problem can be solved using iterative gradient methods, each iteration consisting of the propagation of the actual source in the current medium, with causal attenuation, the propagation of the residuals—acting as if they were sources—backwards in time, with anti-causal attenuation, and the correlation of the two wavefields thus obtained.     

行业空间配置的线性对偶模式

本文研究诸如钢铁、有色金属、基本化工等一类行业的空间配置。由于韦勃、廖施的理论的一些前提条件在这里不存在，使得他们的研究方法显得十分不足。本文作者建立了具有对偶的线性规划模式来解决这类行业的空间配置。这一对模式的解给出了具有最小行业总成本的配置方案和与之相应的产品和原料在各个区位的影子价格。尤其是作者依据空间价格均衡找出了原模式与对偶模式的解之间的关系，这些关系是行业空间配置优化的效率条件。     

探地雷达图像的正演合成与偏移处理

本文讨论了射线追踪法实现二维地电断面的探地雷达图像的正演合成问题,其中包括二维地电断面的数字拟合、多层连续界面介质中高频电磁波的反射与透射统一的运动学边界方程;根据Huygen's-Fresnel原理导出了均匀介质中二维管状体模型的衍射雷达波的计算公式;给出了用射线追踪法合成的二维地电模型的探地雷达图像;用衍射波公式合成了二维管状模型的雷达图像,并与物理模拟作了比较.在此基础上,用Kirchhoff积分法对数字模拟与物理模拟记录的雷达图像作了偏移处理,并给出了实例.     

用数字测高仪漂移测量研究电离层声重波扰动

数字测高仪Digisonde中的漂移测量,常用来研究小尺度电离层扰动,如电离层小不均匀体的漂移。文中提出一种新的分析方法,利用漂移测量数据中多普勒频移和到达角参量的最大熵动态功率谱,估算声重波一类大尺度电离层扰动的水平传播速度和传播方向。作为实例,研究了Millstone Hill测高仪站的漂移观测资料,并对处理结果进行了初步分析。分析结果表明,从数字测高仪漂移测量数据中,可有效地提取声重波一类大尺度电离层扰动的传播参量,在电离层动力过程的研究中很有意义。     

裂隙岩石介质中核素迁移的模型及其算子分裂、迎风、均衡格式

建立了放射性核素在裂隙岩石介质中迁移的双重介质模型,对模型的求解提出了一种新的数值方法—Galerkin有限元法与算子分裂、迎风、均衡格式相结合的新方法,给出了水质模型算子分裂、迎风、均衡格式的稳定性条件,且所得到的计算格式是非负的。最后通过对核素90Sr 100年、99Tc 1000年的预测计算,验证了本文所提方法的有效性和稳定性,并得出了一些有重要意义的结论。     

Rainfall estimation using raingages and radar — A Bayesian approach: 1. Derivation of estimators

Procedures for estimating rainfall from radar and raingage observations are constructed in a Bayesian framework. Given that the number of raingage measurements is typically very small, mean and variance of gage rainfall are treated as uncertain parameters. Under the assumption that log gage rainfall and log radar rainfall are jointly multivariate normal, the estimation problem is equivalent to lognormal co-kriging with uncertain mean and variance of the gage rainfall field.The posterior distribution is obtained under the assumption that the prior for the mean and inverse of the variance of log gage rainfall is normal-gamma 2. Estimate and estimation variance do not have closed-form expressions, but can be easily evaluated by numerically integrating two single integrals. To reduce computational burden associated with evaluating sufficient statistics for the likelihood function, an approximate form of parameter updating is given. Also, as a further approximation, the parameters are updated using raingage measurements only, yielding closed-form expressions for estimate and estimation variance in the Gaussian domain.With a reduction in the number of radar rainfall data in constructing covariance matrices, computational requirements for the estimation procedures are not significantly greater than those for simple co-kriging. Given their generality, the estimation procedures constructed in this work are considered to be applicable in various estimation problems involving an undersampled main variable and a densely sampled auxiliary variable.     

SAR图像纠正的数学模型

本文对Ｇ．Ｋｏｎｅｃｎｙ等提出的平距投影的合成孔径雷达（ＳＡＲ）图像纠正的数学模型给子证明，并指出其近似性。随后推导了斜距投影的雷达图像纠正的数学模型，且以雷达图像的模拟数据，验证了该数学模型与Ｇ．Ｋｏｎｅｃｎｙ的数学模型具有相同的数学精度。     

多卜勒雷达经向速度资料的预处理方法

利用变分法原理提出了对雷达速度资料进行预处理的一种方法，数值试验表明该法能有效地消除雷达速度资料中的小尺度运动径向速度分量及“噪声”。