广西沿岸海平面变化的关联维数研究

采用相空间重构和关联维数方法，研究了广西沿岸月平均海平面变化的分形特征，结果表明，相空间重构的延迟时间τ为3△t，当关联维数趋于饱和时，石头埠，北海和龙尾的相空间维数为7，而涠洲岛为9，它们的关联维数平均值平均值分别为6.6507,6.5414,6.7059.     

龙口平均海平面变化原因的研究

本文通过对实测海平面高度的对比，发现龙口平均海平面从１９９１年起有明显的升高趋势。作者猜测此乃龙口在码头扩建后有系统下沉所致。     

千里岩的气温变化特征及其与西北太平洋副热带高压的关系

本文根据１９６０～１９９４年千里岩和青岛月平均气温资料，分析了千里岩与青岛３５年来气温的年、季变化特征及差别，发现３５年来，千里岩和青岛的年、季平均气温均呈变暖趋势，但千里岩增温幅度较青岛小。从分析中发现，与气温变暖趋势相对应，西北太平洋副热带高压（以下简称副高）面积、强度指数的年、季变化也呈明显的增强、增大趋势。且年平均气温与冬季副高面积指数相关最显著。最后依此关系建立回归方程，并预测了１９９６年千里岩年平均气温的变化趋势。     

A New Methodology for Incorporating Tide Gauge Data in Sea Surface Topography Models

As part of the Vertical Offshore Reference Frames (VORF) project sponsored by the U. K. Hydrographic Office, a new model for Sea Surface Topography (SST) around the British Isles has been developed. For offshore areas (greater than 30 km from the coast), this model is largely derived from satellite altimetry. However, its accuracy and level of detail have been enhanced in coastal areas by the inclusion of not only the 60 PSMSL tide gauges with long-term records around the coasts of the United Kingdom and Ireland but also some 385 gauges established at different epochs and for different observation spans by the U. K. Admiralty. All tide gauge data were brought into a common reference frame by a combination of datum models and direct GPS observations, but a more significant challenge was to bring all short-term sea level observations to an unbiased value at a common epoch. This was achieved through developing a spatial-temporal correlation model for the variations in mean sea level around the British Isles, which in turn meant that gauges with long-term observation spans could be used as control points to improve the accuracy of Admiralty gauges. It is demonstrated that the latter can contribute point observations of mean sea level (MSL) with a precision of 0.078 m. A combination of least squares collocation and interpolation was developed to merge the coastal point and offshore gridded data sets, with particular algorithms having to be developed for different configurations of coastal topology. The resulting model of sea surface topography is shown to present a smooth transition from inshore coastal areas to offshore zones. Further benefits of the techniques developed include an enhanced methodology for detecting datum discontinuities at permanent tide gauges.     

Effect of Sea Level Variation on Tidal Characteristic Values for the East China Sea

1 .IntroductionTheglobalairtemperatureroseabout 0 .5～ 0 .6°Coverthepast 2 0thcentury ,andtheglobalmeansealevelincreasedbyabout2 0cmduringtheperiod .Theregionalmeansealevelriseswiththerisingglobalmeansealevel.Zuoetal.( 1 997)indicatedthatthemeanrisingrateofabsolutemeansealevelalongtheChinacoastontheassumptionofunifiedisostaticdatumis 2mm a .Woodworth( 1 999)analyzedsealevelspanning 1 76 8tothepresentinLiverpool,andobtainedaseculartrendforheperiodupto 1 880of0 .39± 0 .1 7mm a ,andatrendfort…     

成都地区大气平均温度建模及其在GPS/PWV计算中的应用研究

地基GPS气象学的关键技术是在于通过垂直方向上GPS信号的湿分量延迟值来确定大气可降水量-PWV,而这两个物理量间进行转换时需要用到一个关键性的参数——大气加权平均温度。本文首先讨论了估算加权平均温度的几种方法,然后利用成都地区2005年全年的大气探空数据,采用回归分析方法,建立了适合成都地区的大气加权平均温度模型,精度为±2.21K。最后将该公式应用到成都地区PWV的计算中,得到了非常理想的效果。     

利用Jason-1卫星确定全球平均海面高

介绍了卫星测高数据处理的基本方法,分析研究了卫星测高的主要误差,提出了Jason-1的数据编辑准则,并用其测高数据确定了全球平均海面高,与CLS01模型进行了比较分析.     

利用DInSAR和GPS综合方法估计地表3维形变速率

在分析基于贝叶斯统计和马尔科夫随机场理论的解析优化法的基础上,提出以GPS提供的高精度的水平形变速率作为约束,利用直接分解法,将DInSAR的斜距向形变速率分解成垂直形变速率。试验结果表明,直接分解法有效地综合了GPS的高水平精度与DInSAR的高垂直精度优点,且计算简单,克服了解析优化法可能出现的数值不稳定问题。     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.