EOF analysis of a time series with application to tsunami detection

Fragments of deep-ocean tidal records up to 3 days long belong to the same functional sub-space, regardless of the record’s origin. The tidal sub-space basis can be derived via Empirical Orthogonal Function (EOF) analysis of a tidal record of a single buoy. Decomposition of a tsunami buoy record in a functional space of tidal EOFs presents an efficient tool for a short-term tidal forecast, as well as for an accurate tidal removal needed for early tsunami detection and quantification [Tolkova, E., 2009. Principal component analysis of tsunami buoy record: tide prediction and removal. Dyn. Atmos. Oceans 46 (1–4), 62–82] EOF analysis of a time series, however, assumes that the time series represents a stationary (in the weak sense) process. In the present work, a modification of one-dimensional EOF formalism not restricted to stationary processes is introduced. With this modification, the EOF-based de-tiding/forecasting technique can be interpreted in terms of a signal passage through a filter bank, which is unique for the sub-space spanned by the EOFs. This interpretation helps to identify a harmonic content of a continuous process whose fragments are decomposed by given EOFs. In particular, seven EOFs and a constant function are proved to decompose 1-day-long tidal fragments at any location. Filtering by projection into a reduced sub-space of the above EOFs is capable of isolating a tsunami wave within a few millimeter accuracy from the first minutes of the tsunami appearance on a tsunami buoy record, and is reliable in the presence of data gaps. EOFs with ∼$\sim$3-day duration (a reciprocal of either tidal band width) allow short-term (24.75 h in advance) tidal predictions using the inherent structure of a tidal signal. The predictions do not require any a priori knowledge of tidal processes at a particular location, except for recent 49.5 h long recordings at the location.     

纬向平均的风向改变指数与季节转换

通过定义一个能客观定量描述大气环流四维时空变化的风向改变指数WI,用以研究大气环流的时空演变规律和季节转换。在对纬向平均的WI做了经验正交函数(EOF)分析后,得到了其前四个模态。第一模态揭示的是全年平均的WI空间分布的基态。第二模态反映的是WI场的偏差中呈现准调和变化的部分,热带、副热带、温寒带季风区在该模态中均有明显体现。第三模态反映的是WI场的偏差中呈现非调和变化的部分,该部分揭示了因南北半球海陆分布的差异和因大气、海洋流体的非线性效应,其所造成的从春到秋与从秋到春的季节变化的不对称性以及平流层的二月突变现象。第四模态反映的是全球各层盛行风反向区域从春到夏的北进和从夏到秋南撤的现象。     

EXPERIMENTS TO TRACK STORMS USING MODERN OPTIMIZATION ALGORITHMS

Storm identification and tracking based on weather radar data are essential to nowcasting and severe weather warning. A new two-dimensional storm identification method simultaneously seeking in two directions is proposed, and identification results are used to discuss storm tracking algorithms. Three modern optimization algorithms (simulated annealing algorithm, genetic algorithm and ant colony algorithm) are tested to match storms in successive time intervals. Preliminary results indicate that the simulated annealing algorithm and ant colony algorithm are effective and have intuitionally adjustable parameters, whereas the genetic algorithm is unsatisfactorily constrained by the mode of genetic operations. Experiments provide not only the feasibility and characteristics of storm tracking with modern optimization algorithms, but also references for studies and applications in relevant fields.     

基于5分钟雷达资料降水估测模型的比较

For the Z-R relationship in radar-based rainfall estimation, the distribution of corresponding R values for a given Z value （or the corresponding Z value for a given R value） may be highly skewed. However, the traditional power-law model is physically deduced and fitted under the normal-distribution presumption of radar wave echoes associated with a rain rate value, and it may not be very appropriate. Considering this problem, the authors devised several generalized linear models with different forms and distribution presumptions to represent the Z-R relationship. Radar-reflectivity scans observed by a CINRAD/SC Doppler radar and 5-minute rainfall accumulation recorded by 10 ground gauges were used to fit these models. All data used in this study were collected during some large rainfalls of the period from 2005 to 2007. The radar and all gauges were installed in the catchment of the Yishu River, a branch of the Huaihe River in China. Three models based on normal distribution and a dBZ presumption of gamma distribution were fitted using maximum-likelihood techniques, which were resolved by genetic algorithms. Comparisons of estimated maximized likelihoods based on assumptions of gamma and normal distribution showed that all generalized linear models （GLMs） of presumed gamma distribution were better fitted than GLMs based on normal distribution. In a comparison of maximum-likelihood, the differences between these three models were small. Three error statistics were used to assess the agreement between radar estimated rainfall and gauge rainfall： relative bias （B）, root mean square error （RMSE）, and correlation coefficient （r）. The results showed that no one model was excellent in all criteria. On the whole, the GLM-based models gave smaller relative bias than the traditional power-law model. It is suggested that validations conducted in many previous works should have been made against a specific criterion but overlooked others.     

Integration schemes with substepping algorithms for creep analysis in geomaterials

This paper describes a robust and efficient methodology for predicting displacements, deformations, and stresses in geomaterials that are susceptible to creep. The methodology is based on two integration schemes, which consider substepping algorithms. The first scheme is used for integrating space-time relations in a global sense, whereas the second scheme is used for integrating stress-stain relations in a local sense. Different from previous studies, both integration schemes are easy to implement and general in the sense that they can be applied to any type of creep law. Through an in-house finite element simulator, several numerical tests are performed. They include triaxial and wellbore closure analyses considering soft soils and salt rocks. The results show that the combination of both schemes leads to stable and accurate solutions with reduced computational time.     

顶点Pk覆盖问题的研究综述

顶点覆盖问题是经典的NP完全问题,在排序、计算机网络等现实生活中有许多的应用.近几年来,许多研究者开始探究它的推广形式——顶点P_k覆盖（VCP_k）问题,即寻找一个顶点子集,从拓扑结构图中删除后使得剩下的顶点导出的子图不包含P_k路,其中P_k是指包含k个顶点的路.本文简单介绍了VCP_k问题的应用背景,归纳了它在近似算法、精确算法、参数化算法3个方面的主要研究进展,并分析了一些主要的方法和技巧.在此基础上,对VCP_k问题及其相关问题的研究前景进行了展望.     

中国近117年年平均气温变化的区域特征研究

用ＲＥＯＦ法以１６０站近４７年的年平均气温变化为依据将我国分为８个年平均气温变化区．以国内及临近国家长资料序列为基础,建立了反映我国各区及全国近１１７年的年平均气温标准化序列．从近１１７年气温变化来看,Ｉ～ＶＩＩ区都有升温趋势,全国平均气温的升温趋势为０．７６℃／１００ａ,北半球陆地则为０．６４℃／１００ａ．最大的为东北,华北、新疆,最小的为四川、贵州．从年代际尺度看我国与北半球陆地的平均气温的变化是相当一致的,相关系数高达０．９３,主要差别在于年际变化．     

关于-遗传算法收敛性的注记

遗传算法是一种受到广泛注意的全局优化算法,已经在包括地震工程的很多领域中获得应用.本文将结合这一算法的实际操作步骤,对简单遗传算法的不收敛性和包含最优个体保护策略的遗传算法的收敛性给出一个简要的证明.     

Dislocation parameters of Gonghe earthquake jointly inferred by using genetic algorithms and least squares method

IntroductionOnApril26,1990,anMs=7.0strongewthquakeoccurrednearTanggannofarminnerGonghebasin,QinghaiProvince.Aftertheevent,theSecondCrustalDeformationMonitoringCented,ChinaSeismologicalBureau,immediatelyconductedarelevelingsurveyalongthetwolinesacrosstheearthquakeregioninMayandJuneof1990,thereweretotally27levelingpoilltsreoccupied.Alongwiththelastpre-earthquakelevelingsurveyin1978,abigupliftwiththemaximum345mmatDH20aswellasasmallsubsidencewiththemaximum45.1mmatDH16havebeendetected(Figur…     

基于匹配追踪和遗传算法的大地电磁噪声压制

针对匹配追踪计算量大、大地电磁数据处理效率低的问题,提出基于匹配追踪和遗传算法的大地电磁噪声压制方法.首先,利用Gabor原子构建过完备原子库,并对过完备原子库集合进行划分.然后,借助遗传算法的自适应性,快速搜寻最优匹配原子及所在位置.最后,运用最优匹配原子对待处理信号进行稀疏分解,重构有用信号.通过对计算机模拟的典型强干扰和矿集区实测大地电磁数据进行分析处理,实验结果表明,相对于匹配追踪和正交匹配追踪,文中所提方法能从过完备原子库中快速、自适应地选取最优匹配原子与不同噪声干扰类型高精度的匹配,极大地提升了计算效率;大地电磁时间域序列中的大尺度强干扰被有效剔除,视电阻率曲线更为光滑、连续,低频段的数据质量得到明显改善.