一种从经纬度网格资料获得正方形网格资料的插值方案

本文提出了一种把经纬度网格点上的气象要素值转化到正方形网格点上的多项式插值方案。对插值方案作了稳定性讨论和定量及定性分析。结果表明,这种插值方案方便可行,而且有广泛适用性。     

拉格朗日插值法在渔业资源评估上的应用

应用拉格朗日插值法拟合鱼类体长与体重之间的函数关系 ,使用计算机VB语言进行编写程序 ,与常用的线性回归法、Ricker法进行比较 ,得出体长与体重的关系式为Pn(x) =∑nk=0yklk(x)=∑nk=0(∏nj=0j≠kx-xjxk-xj)yk 。     

The state of land subsidence and prediction approaches due to groundwater withdrawal in China

This article gives a general introduction to land subsidence with the prediction approaches due to withdrawal of groundwater in three subsided/subsiding regions in China: the deltaic plain of Yangtse River (YRDP), North China Plain (NCP), and Fenwei Plain (FP). On YRDP, Shanghai is the typical subsided/subsiding city; on NCP Tianjin is the typical subsided/subsiding city, and on FP Taiyuan is the typical subsided/subsiding city. The subsided area with subsidence over 200 mm on YRDP is about 10,000 km² and the maximum subsided value reached 2.9 m at Shanghai; on NCP the subsided area reached 60,000 km² with the maximum subsidence of 3.9 m at Tianjing; on FP the subsided area is relatively smaller than that on the other two plains and is about 1,135 km² with maximum subsidence of 3.7 m at Taiyuan city. In order to protect the civil and industrial facilities, it is necessary to predict the future development of land subsidence based on present state. Many researchers proposed several approaches to predict the land subsidence due to groundwater withdrawal according to different geological conditions and groundwater withdrawal practice. This article classifies these approaches into five categories: (i) statistical methods; (ii) 1D numerical method; (iii) quasi-3D seepage model; (iv) 3D seepage model; (v) fully coupled 3D model. In China, the former four categories are presently employed in the prediction practice and their merits and demerits are discussed. According to the prediction practice, 3D seepage model is the best method presently.     

Climate change mitigation scenarios and policies and measures: the case of Kazakhstan

This article illustrates the main difficulties encountered in the preparation of GHG emission projections and climate change mitigation policies and measures (P&M) for Kazakhstan. Difficulties in representing the system with an economic model have been overcome by representing the energy system with a technical-economic growth model (MARKAL-TIMES) based on the stock of existing plants, transformation processes, and end-use devices. GHG emission scenarios depend mainly on the pace of transition in Kazakhstan from a planned economy to a market economy. Three scenarios are portrayed: an incomplete transition, a fast and successful one, and even more advanced participation in global climate change mitigation, including participation in some emission trading schemes. If the transition to a market economy is completed by 2020, P&M already adopted may reduce emissions of CO₂ from combustion by about 85 MtCO₂ by 2030 – 17% of the emissions in the baseline (WOM) scenario. One-third of these reductions are likely to be obtained from the demand sectors, and two-thirds from the supply sectors. If every tonne of CO₂ not emitted is valued up to US$10 in 2020 and $20 in 2030, additional P&M may further reduce emissions by 110 MtCO₂ by 2030.     

GPS精密星历的轨道内插方法比较

在高精度的GPS测量数据处理过程中,获取高精度卫星轨道是重要的环节,内插是获取任意历元的精密轨道信息的重要手段,对拉格朗日插值法、牛顿插值法、Neville逐次线性内插法和切比雪夫插值法的效果做了细致的比较,并在此基础上进一步探讨了短时间外推的效果。     

基于EMD方法的地心运动时间序列分析

针对时间序列混有的高频信息会影响地心运动规律分析的问题,采用网平移法对IGS提供的GNSS周解进行解算,得到2007-2017年地心运动时间序列,对其进行分解重构,剔除高频项,并利用重构时序对地心运动规律作进一步分析探讨。结果表明：本文解算的地心运动在T_x、T_y和T_z方向的精度均为毫米级。EMD方法重构的时序保留了原序列的基本信息,且抑制了高频项的影响,提高了周期贡献率,3个方向的贡献率分别提高了12.3%、16.7%及6.3%。通过分析重构后的时序发现,周年项振幅为各周期对应振幅的最大值,分别为2.32、1.89和2.07 mm;T_y和T_z方向长期变化趋势较T_x更为明显,分别为0.13和-0.27 mm/a;半年项较小,且在T_x和T_y方向上具有时变性。此外,还发现了一些其他较小的年际变化。     

非滑动式与滑动式拉格朗日插值法在BDS精密星历内插中的比较分析

在高精度卫星定位数据处理中,对精密星历内插处理是一项基础工作。利用拉格朗日多项式插值法,从三方面对北斗导航系统精密星历进行了内插精度分析,分别是非滑动式拉格朗日插值法在不同插值阶数下的插值精度,滑动式拉格朗日插值法在不同插值阶数下的插值精度,以及非滑动式和滑动式拉格朗日插值法在同一插值阶数下的插值精度。算例结果表明,如果将插值阶数作为横坐标,且向右为正,插值精度为纵坐标,向上为正,那么非滑动式拉格朗日插值法的插值精度与插值阶数呈现一个开口向下的抛物线关系;滑动式拉格朗日插值法的插值阶数高于7后,插值精度处于稳定且高的状态,插值误差最大值为1. 1 mm左右,均方差在0. 4 mm左右;在同一插值阶数下,滑动式拉格朗日插值法的精度高于非滑动式拉格朗日插值法的精度。     

高度计波高数据同化对印度洋海域海浪模式预报影响研究

为提高海浪模式预报的精度,改善初始场是途径之一。研制了基于最优插值(OI)方法的海浪数据同化并行程序模块,并将其植入第三代海浪模式WAVEWATCH IIITM,建立了印度洋海域海浪同化预报方法,使用卫星高度计波高数据进行了同化预报试验。OI模块的并行设计使得植入同化模块的海浪模式仍能以并行方式运行。文中5°S以北印度洋海域为目标区域,嵌套在WAVE-WATCH IIITM的全球网格中,使得目标区域开边界条件得到较好解决。同化数据使用Jason-2高度计测量有效波高(SWH)沿轨数据。海浪同化预报模式由大气模式WRF(Weather Research andForecasting)输出的1小时一次的海面10 m风场驱动。将同化的模式结果(SWH)、无同化的模式结果(SWH)分别与高度计沿轨数据(SWH)进行比较,表明同化改善模式预报初始场的效果是明显的。以同化初始场出发进行海浪预报试验,结果表明,高度计波高数据同化在一定程度上可改进海浪短期预报的精度。     

位场曲化平的插值-迭代法

将起伏曲面B上的位场向下延拓至曲面最低点的平面A的插值-迭代法步骤是：1）将曲面B上的场值放置在水平面A上具有相同水平坐标的点上,作为A上的初值;2）用若干水平面切割B,从A的初值,用快速傅里叶变换法(FFT)向上延拓出这些平面的场值,用插值的方法从这些平面的场值计算曲面B的场值;3）根据B上的实测值与计算值的差值,对A上的值进行加权改正;4）重复步骤2）和3）,直到B上的差值小到可以忽略.这种插值-迭代法具有高的计算速度,比通常的FFT法延拓得更深,可以超过10倍点距.文中给出计算实例.     

Conditional Probability Approach of the Assessment of Tsunami Potential: Application in Three Tsunamigenic Regions of the Pacific Ocean

We develop stochastic approaches to determine the potential for tsunami generation from earthquakes by combining two interrelated time series, one for the earthquake events, and another for the tsunami events. Conditional probabilities for the occurrence of tsunamis as a function of time are calculated by assuming that the inter-arrival times of the past events are lognormally distributed and by taking into account the time of occurrence of the last event in the time series. An alternative approach is based on the total probabilitiy theorem. Then, the probability for the tsunami occurrence equals the product of the ratio, r (= tsunami generating earthquakes/total number of earthquakes) by the conditional probability for the occurrence of the next earthquake in the zone. The probabilities obtained by the total probability theorem are bounded upwards by the ratio r and, therefore, they are not comparable with the conditional probabilities. The two methods were successfully tested in three characteristic seismic zones of the Pacific Ocean: South America, Kuril-Kamchatka and Japan. For time intervals of about 20 years and over the probabilities exceed 0.50 in the three zones. It has been found that the results depend on the approach applied. In fact, the conditional probabilities of tsunami occurrence in Japan are slightly higher than in the South America region and in Kuril-Kamchatka they are clearly lower than in South America. Probabilities calculated by the total probability theorem are systematically higher in South America than in Japan while in Kuril-Kamchatka they are significantly lower than in Japan. The stochastic techniques tested in this paper are promising for the tsunami potential assessment in other tsunamigenic regions of the world.