复杂地质地貌区地球化学异常识别非线性研究

应用现代数学非线性理论与方法对西南地区三江北段恒星错测区1：50000水系沉积物测量数据进行了异常识别的非线性研究。率先尝试用人工免疫算法对地球化学数据降噪，并采用分形含量面积法对降噪数据寻求异常下限，从而弥补传统方法确定异常下限出现的缺陷。在此基础上构造了测区地球化学元素分布的分形插值曲面，做出异常等值线图和异常分形插值图。处理结果令人满意。     

A method for assessing local variation among scattered measurements

Automatic triangulation of scattered locations permits analysis of local variation in a dependent variable through calculation of a roughness index. This is approached by treating triangles of the triangulation (including the dependent variable) as vectorial structures, and accumulating at each data point the vector sum of the cluster of triangles surrounding it. The roughness index is defined as the complement of the ratio of the area of a triangle cluster to the area of component triangles as projected onto a gradient plane defined by their vector sum. The roughness index provides a measure of consistency of data values relative to surrounding observations and can be interpreted as a local index of reliability of interpolation.     

Quality control and homogeneity of Turkish precipitation data

Outlier trimming and homogeneity checking/correction were performed on the monthly precipitation time series of various lengths from 267 stations in Turkey. Outlier values are usually found during dry summer months, and are concentrated mostly over the southern parts of the country, where the dry period is most pronounced, implying natural extremes rather than wrong measurements. Homogeneity analysis was done using the Standard Normal Homogeneity Test, on an individual monthly basis, which led to many non‐testable series due to lack of reference stations, especially during summer months. Yet, remaining testable months were usually helpful for the assessment of homogenity, revealing a well distributed set of stations that proved to be homogeneous. There were still a number of stations which either could not be tested efficiently, or were classified as inhomogeneous. Lack of metadata is argued to be largely responsible for inefficient homogeneity testing. Copyright © 2008 John Wiley & Sons, Ltd.     

FPGA在时间测量中的应用

利用现场可编程门阵列(FPGA)的内部资源，在FPGA中实现了高精度时间间隔测量系统的设计。     

多波束海底地形数据空值区域的插值填补方法分析

由于边界误差剔除、测线布设过于稀疏等因素影响,利用多波束数据生成高分辨率DEM时,常常会存在一些水深空值单元或区域。为了对这些数据空值区域进行可靠填补,在考虑地形特征及空值单元大小的情况下,利用IDW、Spline、Kriging、TIN四种插值方法对某区域多波束海底DEM空值区域进行直接内插填补试验,结果表明插值算法的选取不仅受地形类型限制,而且与空值区域的大小也密切相关。     

基于改进的双线性旅行时插值的三维射线追踪

旅行时射线追踪的精度和效率一直是影响地震层析成像质量和效果的关键因素。这里在三维近地表速度结构层析成像应用中,结合旅行时插值算法与最短路径算法,提出了一种三维初至波射线追踪算法。通过对双线性旅行时插值算法进行改进,并运用判定条件与简化插值计算公式进行快速计算,减少了插值次数,降低了运算量,同时也保证了较高的射线精度,有效地解决了三维射线追踪算法的计算低效率,射线精度不高的问题。     

Propagation of vertical and horizontal source data errors into a TIN with linear interpolation

Digital elevation models (DEMs) have been widely used for a range of applications and form the basis of many GIS-related tasks. An essential aspect of a DEM is its accuracy, which depends on a variety of factors, such as source data quality, interpolation methods, data sampling density and the surface topographical characteristics. In recent years, point measurements acquired directly from land surveying such as differential global positioning system and light detection and ranging have become increasingly popular. These topographical data points can be used as the source data for the creation of DEMs at a local or regional scale. The errors in point measurements can be estimated in some cases. The focus of this article is on how the errors in the source data propagate into DEMs. The interpolation method considered is a triangulated irregular network (TIN) with linear interpolation. Both horizontal and vertical errors in source data points are considered in this study. An analytical method is derived for the error propagation into any particular point of interest within a TIN model. The solution is validated using Monte Carlo simulations and survey data obtained from a terrestrial laser scanner.     

斜压罗斯贝变形半径优化的误差相关尺度及其对最优插值效果的改进

利用最优插值方法进行三维温盐场重构时,如何正确估算背景场误差协方差至关重要。针对背景场误差协方差主要取决于误差相关尺度的问题,本文提出了用中国东部海域及其附近海域的最新的高分辨率气候态数据进行斜压罗斯贝变形半径优化其误差相关尺度计算的研究方案和技术途径;对比分析了均一化相关尺度方案和法国ISAS系统尺度方案,讨论了变形半径对最优插值的影响。结果表明:均一化相关尺度方案的均方根误差小于ISAS方案,但温度场过于平滑,难以刻画一些重要的物理现象;相比而言,本文提出的基于变形半径的相关尺度方案在取2倍变形半径时,不仅均方根误差在各水平层较小,且温度场能够更好地刻画四国海盆海域涡旋及黑潮影响的温度场三维结构。由于实际海洋中各层物理过程的尺度存在差异,实际应用时各层的最优尺度设置也应有所不同。     

基于ROMS模式利用集合最优插值同化沿轨海表面高度异常数据的研究

The ensemble optimal interpolation （EnOI） is applied to the regional ocean modeling system （ROMS） with the ability to assimilate the along-track sea level anomaly （TSLA）. This system is tested with an eddy-resolving system of the South China Sea （SCS）. Background errors are derived from a running seasonal ensemble to account for the seasonal variability within the SCS. A fifth-order localization function with a 250 km localization radius is chosen to reduce the negative effects of sampling errors. The data assimilation system is tested from January 2004 to December 2006. The results show that the root mean square deviation （RMSD） of the sea level anomaly decreased from 10.57 to 6.70 cm, which represents a 36.6% reduction of error. The data assimilation reduces error for temperature within the upper 800 m and for salinity within the upper 200 m, although error degrades slightly at deeper depths. Surface currents are in better agreement with trajectories of surface drifters after data assimilation. The variance of sea level improves significantly in terms of both the amplitude and position of the strong and weak variance regions after assimilating TSLA. Results with AGE error （AGE） perform better than no AGE error （NoAGE） when considering the improvements of the temperature and the salinity. Furthermore, reasons for the extremely strong variability in the northern SCS in high resolution models are investigated. The results demonstrate that the strong variability of sea level in the high resolution model is caused by an extremely strong Kuroshio intrusion. Therefore, it is demonstrated that it is necessary to assimilate the TSLA in order to better simulate the SCS with high resolution models.     

Improvement of the semi-Lagrangian advection scheme in the GRAPES model: Theoretical analysis and idealized tests

ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion （referred to as the SL_CL scheme）. The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich＇s semi-Lagrangian advection scheme （referred to as the R2007 scheme）, the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that： （1） the damping effects were remarkably reduced with the Re_R2007 scheme; and （2） the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme＇s advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided.