Variation in joint mode of the tropical Indian-Pacific Ocean thermodynamic anomaly

Previous research has defined the index of the Indian-Pacific thermodynamic anomaly joint mode (IPTAJM) and suggested that the winter IPTAJM has an important impact on summer rainfall over China. However, the possible causes for the interannual and decadal variability of the IPTAJM are still unclear. Therefore, this work investigates zonal displacements of both the western Pacific warm pool (WPWP) and the eastern Indian Ocean warm pool (EIOWP). The relationships between the WPWP and the EIOWP and the IPTAJM are each examined, and then the impacts of the zonal wind anomalies over the equatorial Pacific and Indian Oceans on the IPTAJM are studied. The WPWP eastern edge anomaly displays significant interannual and decadal variability and experienced a regime shift in about 1976 and 1998, whereas the EIOWP western edge exhibits only distinct interannual variability. The decadal variability of the IPTAJM may be mainly caused by both the zonal migration of the WPWP and the 850 hPa zonal wind anomaly over the central equatorial Pacific. On the other hand, the zonal migrations of both the WPWP and the EIOWP and the zonal wind anomalies over the central equatorial Pacific and the eastern equatorial Indian Ocean may be all responsible for the interannual variability of the IPTAJM.

Mixed-layer water oscillations in tropical Pacific for ENSO cycle

The main modes of interannal variabilities of thermocline and sea surface wind stress in the tropical Pacific and their interactions are investigated, which show the following results. (1) The thermocline anomalies in the tropical Pacific have a zonal dipole pattern with 160 W as its axis and a meridional seesaw pattern with 6-8 degrees N as its transverse axis. The meridional oscillation has a phase lag of about 90 to the zonal oscillation, both oscillations get together to form the El Nino/La Nina cycle, which behaves as a mixed layer water oscillates anticlockwise within the tropical Pacific basin between equator and 12 degrees N. (2) There are two main patterns of wind stress anomalies in the tropical Pacific, of which the first component caused by trade wind anomaly is characterized by the zonal wind stress anomalies and its corresponding divergences field in the equatorial Pacific, and the abnormal cross- equatorial flow wind stress and its corresponding divergence field, which has a sign opposite to that of the equatorial region, in the off-equator of the tropical North Pacific, and the second component represents the wind stress anomalies and corresponding divergences caused by the ITCZ anomaly. (3) The trade winds anomaly plays a decisive role in the strength and phase transition of the ENSO cycle, which results in the sea level tilting, provides an initial potential energy to the mixed layer water oscillation, and causes the opposite thermocline displacement between the west side and east side of the equator and also between the equator and 12 degrees N of the North Pacific basin, therefore determines the amplitude and route for ENSO cycle. The ITCZ anomaly has some effects on the phase transition. (4) The thermal anomaly of the tropical western Pacific causes the wind stress anomaly and extends eastward along the equator accompanied with the mixed layer water oscillation in the equatorial Pacific, which causes the trade winds anomaly and produces the anomalous wind stress and the corresponding divergence in favor to conduce the oscillation, which in turn intensifies the oscillation. The coupled system of ocean-atmosphere interactions and the inertia gravity of the mixed layer water oscillation provide together a phase-switching mechanism and interannual memory for the ENSO cycle. In conclusion, the ENSO cycle essentially is an inertial oscillation of the mixed layer water induced by both the trade winds anomaly and the coupled ocean-atmosphere interaction in the tropical Pacific basin between the equator and 12 degrees N. When the force produced by the coupled ocean-atmosphere interaction is larger than or equal to the resistance caused by the mixed layer water oscillation, the oscillation will be stronger or maintain as it is, while when the force is less than the resistance, the oscillation will be weaker, even break.

大气季节内振荡对东亚夏季风活动的影响研究

利用ERA40逐日再分析资料、NCEP/NCAR2逐日再分析资料、中国740个测站日降水资料、上海台风研究所提供的西太平洋热带气旋资料、Kaplan等重建的月平均SSTA资料、NOAA逐日长波辐射（OLR）等资料，应用离散功率谱分析、带通滤波、EOF分析等统计方法，研究了东亚夏季风(EASM）的移动特征、东亚地区季节内振荡(ISO)的基本特征、季节内振荡对东亚夏季风活动的影响、季节内振荡对东亚夏季风异常活动的影响机理。主要结论如下： （1）综合动力和热力因素定义了可动态描述东亚夏季风移动和强度的指数，并利用该指数研究了东亚夏季风的爆发和移动的季节内变化及其年际和年代际变化特征。研究发现，气候平均东亚夏季风前沿分别在28候、33候、36候、38候、40候、44候出现了明显的跳跃。东亚夏季风活动具有显著的年际变率，主要由于季风前沿在某些区域异常停滞和突然跨越北跳或南撤引起，造成中国东部旱涝灾害频繁发生。东亚夏季风的活动具有明显的年代际变化，在1965年、1980年、1994年发生了突变，造成中国东部降水由“南旱北涝”向“南涝北旱”的转变。 （2）东亚季风区季节内变化具有10~25d和30~60d两个波段的季节内振荡周期，以30-60d为主。存在三个主要低频模态，第一模态主要表征了EASM在长江中下游和华北地区活动期间的低频形势；第二模态印度洋-菲律宾由低频气旋式环流控制，主要表现了ISO在EASM爆发期间的低频形势；第三模态主要出现在EASM在华南和淮河活动期间的低频形势。第一模态和第三模态是代表东亚夏季风活动异常的主要低频形势。 （3）热带和副热带地区ISO总是沿垂直切变风的垂直方向传播。因此，在南海-菲律宾东北风垂直切变和副热带西太平洋北风垂直切变下，大气热源激发菲律宾附近交替出现的低频气旋和低频反气旋不断向西北传播，副热带西太平洋ISO以向西传播为主。中高纬度地区，乌拉尔山附近ISO以向东、向南移动或局地振荡为主；北太平洋中部ISO在某些情况下向南、向西传播。 （4）季风爆发期，伴随着热带东印度洋到菲律宾一系列低频气旋和低频反气旋， 冷空气向南输送，10~25天和30~60天季节内振荡低频气旋同时传入南海加快了南海夏季风的爆发。在气候态下，ISO活动表现的欧亚- 太平洋（EAP）以及太平洋-北美（PNA）低频波列分布特征(本文提出的EAP和PNA低频波列与传统意义上的二维定点相关得到的波列不同)。这种低频分布形式使得欧亚和太平洋中高纬度的槽、脊及太平洋副热带高压稳定、加强，东亚地区的低频波列则成为热带和中高纬度ISO相互作用影响东亚夏季风活动的纽带。不同的阶段表现不同的低频模态，30~60d低频模态的转变加快了EASM推进过程中跳跃性；30-60d低频模态的维持使得EASM前沿相对停滞。 （5）30-60d滤波场，菲律宾海域交替出现的低频气旋和低频反气旋不断向西北传播到南海-西太平洋一带。当南海-西太平洋地区低频气旋活跃时，季风槽加强、东伸，季风槽内热带气旋（TC）频数增加；当南海-西太平洋低频反气旋活跃时，季风槽减弱、西退，TC处于间歇期，生成位置不集中。 （6）在El Nino态下，大气季节内振荡偏弱，北传特征不明显，但ISO由中高纬度北太平洋中部向南和副热带西太平洋向西的传播特征显著，东亚地区ISO活动以第三模态为主，EASM集中停滞在华南和淮河流域，常伴随着持续性区域暴雨的出现，易造成华南和江淮流域洪涝灾害，长江和华北持续干旱。在La Nina态下，大气季节内振荡活跃，且具有明显的向北传播特征，PNA低频波列显著，东亚地区ISO活动以第一模态单峰为主；EASM主要停滞在长江中下游和华北地区，这些地区出现异常持续强降水，华南和淮河流域多干旱；在El Nino态向La Nina态转换期，ISO活动以第一模态双峰为主，长江中下游常常出现二度梅。

中国近海和西北太平洋温跃层时空变化分析、模拟及预报

温度跃层是反映海洋温度场的重要物理特性指标，对水下通讯、潜艇活动及渔业养殖、捕捞等有重要影响。本文利用中国科学院海洋研究所“中国海洋科学数据库”在中国近海及西北太平洋（110ºE-140ºE,10ºN-40ºN）的多年历史资料（1930-2002年，510143站次），基于一种改进的温跃层判定方法，分析了该海域温跃层特征量的时空分布状况。同时利用Princeton Ocean Model(POM)，对中国近海，特别是东南沿海的水文结构进行了模拟，研究了海洋水文环境对逆温跃层的影响。最后根据历史海温观测资料，利用EOF分解统计技术，提出了一种适于我国近海及毗邻海域，基于现场有限层实测海温数据，快速重构海洋水温垂直结构的统计预报方法，以达到对现场温跃层的快速估计。 历史资料分析结果表明，受太阳辐射和风应力的影响，20°N以北研究海域，温跃层季节变化明显，夏季温跃层最浅、最强，冬季相反，温跃层厚度的相位明显滞后于其他变量，其在春季最薄、秋季最厚。12月份到翌年3月份，渤、黄及东海西岸，呈无跃层结构，西北太平洋部分海域从1月到3月份，也基本无跃层结构。在黄海西和东岸以及台湾海峡附近的浅滩海域，由于风力搅拌和潮混合作用，温跃层出现概率常年较低。夏季，海水层化现象在近海陆架海域得到了加强，陆架海域温跃层强度季节性变化幅度(0.31°C/m)明显大于深水区(约0.05°C/m)，而前者温跃层深度和厚度的季节性变化幅度小于后者。20°N以南研究海域，温跃层季节变化不明显。逆温跃层主要出现在冬、春季节（10月-翌年5月）。受长江冲淡水和台湾暖流的影响，东南沿海区域逆温跃层持续时间最长，出现概率最大，而在山东半岛北及东沿岸、朝鲜半岛西及北岸，逆温跃层消长过程似乎和黄海暖流有关。多温跃层结构常年出现于北赤道流及对马暖流区。在黑潮入侵黄、东、南海的区域，多温跃层呈现明显不同的季节变化。在黄海中部，春季多温跃层发生概率高于夏季和秋季，在东海西部，多跃层主要出现在夏季，在南海北部，冬季和春季多温跃层发生概率大于夏季和秋季。这些变化可能主要受海表面温度变化和风力驱动的表层流的影响。 利用Princeton Ocean Model(POM)，对中国东南沿海逆温跃层结构进行了模拟，模拟结果显示，长江冲淡水的季节性变化以及夏季转向与实际结果符合较好，基本再现了渤、黄、东海海域主要的环流、温盐场以及逆温跃层的分布特征和季节变化。通过数值实验发现，若无长江、黄河淡水输入，则在整个研究海域基本无逆温跃层出现，因此陆源淡水可能是河口附近逆温跃层出现的基本因素之一。长江以及暖流（黑潮和台湾暖流）流量的增加，均可在不同程度上使逆温跃层出现概率及强度、深度和厚度增加，且暖流的影响更加明显。长江对东南沿海逆温跃层的出现，特别是秋季到冬季初期，有明显的影响，使长江口海域逆温跃层位置偏向东南。暖流对于中国东南沿海的逆温跃层结构，特别是初春时期，有较大影响，使长江口海域的逆温跃层位置向东北偏移。 通过对温跃层长期变化分析得出，黄海冷水团区域，夏季温跃层强度存在3.8年左右的年际变化及18.9年左右的年代际变化，此变化可能主要表现为对当年夏季和前冬东亚地区大气气温的热力响应。东海冷涡区域，夏季温跃层强度存在3.7年的年际变化，在El Nino年为正的强度异常，其可能主要受局地气旋式大气环流变异所影响。谱分析同时表明，该海域夏季温跃层强度还存在33.2年的年代际变化，上世纪70年代中期，温跃层强度由弱转强，而此变化可能与黑潮流量的年代际变化有关。 海洋水温垂直结构的统计预报结果显示，EOF分解的前四个主分量即能够解释原空间点温度距平总方差的95%以上，以海洋表层附近观测资料求解的特征系数推断温度垂直结构分布的结果最稳定。利用东海陆架区、南海深水区和台湾周边海域三个不同区域的实测CTD样本廓线资料，对重构模型的检验结果表明，重构与实测廓线的相关程度超过95%的置信水平。三个区重构与实测温度廓线值的平均误差分别为0.69℃，0.52℃，1.18℃，平均重构廓线误差小于平均气候偏差，统计模式可以很好的估算温度廓线垂直结构。东海陆架海区温度垂直重构廓线与CTD观测廓线获得的温跃层结果对比表明，重构温跃层上界、下界深度和强度的平均绝对误差分别为1.51m、1.36m和0.17℃/m，它们的平均相对误差分别为24.7%、8.9%和22.6%，虽然温跃层深度和强度的平均相对误差较大，但其绝对误差量值较小。而在南海海区，模型重构温跃层上界、下界和强度的平均绝对预报误差分别为4.1m、27.7m和0.007℃/m，它们的平均相对误差分别为16.1%、16.8%和9.5%，重构温跃层各特征值的平均相对误差都在20%以内。虽然南海区温跃层下界深度平均绝对预报误差较大，但相对于温跃层下界深度的空间尺度变化而言（平均温跃层下界深度为168m），平均相对误差仅为16.8%。因此说模型重构的温度廓线可以达到对我国陆架海域、深水区温跃层的较好估算。 基于对历史水文温度廓线观测资料的分析及自主温跃层统计预报模型，研制了实时可利用微机简单、快捷地进行温跃层估算及查询的可视化系统，这是迄今进行大范围海域温跃层统计与实时预报研究的较系统成果。

Typhoon effects on DOC dynamics in a phosphate-limited reservoir

To explore typhoon effects on dissolved organic carbon (DOC) dynamics, field investigations (tributary and dam site) and laboratory experiments (bioassay and DOC consumption) were conducted in a subtropical reservoir. A tributary survey indicated that after typhoon disruption, upstream areas were the sources of phosphate (P) but not DOC for the dam site located downstream. Bioassay experiments verified P-limitation on bacteria and phytoplankton during summer stratification, and bacteria showed a faster response than algae to added P. Experiments indicated that DOC consumption was determined by the availability of P. The 4 yr typhoon period (June-September) data of the dam site denoted that DOC concentration (27 to 270 mu M C) and its rate of change (-13 to 24 mu M C d(-1)) varied more dramatically in the weak (2006 and 2007) than in the strong (2004 and 2005) typhoon years. The negative correlation of DOC with the ratio of bacterial production (BP) to primary production (PP) in the euphotic zone (0 to 10 m) signified the interactive effects of auto- and heterotrophic processes on DOC variation. In the aphotic zone, the variation of DOC could be ascribed to the change of BP, which showed a positive correlation with P concentrations. This study documents that DOC concentration in the studied system varied at multiple time scales. Such variation can be explained by the decoupling between BP and PP, which is believed to be a function of the limiting nutrient's availability. More importantly, this study suggests that the P supply introduced by strong typhoons might have substantiated a tighter coupling between BP and PP, so that the amplitude of DOC oscillation during the summer period was effectively reduced.

Genetic variation and phylogeographic history of Picea likiangensis revealed by RAPD markers

Repeated cycles of retreat and recolonization during the Quaternary ice ages are thought to have greatly influenced current species distributions and their genetic diversity. It remains unclear how this climatic oscillation has affected the distribution of genetic diversity between populations of wind-pollinated conifers in the Qinghai-Tibetan region. In this study, we investigated the within-species genetic diversity and phylogenetic relationships of Picea likiangensis, a dominant forest species in this region using polymorphic DNA (RAPD) markers. Our results suggest that this species has high overall genetic diversity, with 85.42% of loci being polymorphic and an average expected heterozygosity (H (E)) of 0.239. However, there were relatively low levels of polymorphism at population levels and the differences between populations were not significant, with percentages of polymorphic bands (PPB) ranging from 46.88 to 69.76%, Nei's gene diversity (H (E)) from 0.179 to 0.289 and Shannon's indices (Hpop) from 0.267 to 0.421. In accordance with our proposed hypothesis, a high level of genetic differentiation among populations was detected based on Nei's genetic diversity (G (ST) = 0.256) and AMOVA analysis (Phi (st) = 0.236). Gene flow between populations was found to be limited (Nm = 1.4532) and far lower than reported for other conifer species with wide distribution ranges from other regions. No clusters corresponding to three morphological varieties found in the south, north and west, respectively, were detected in either UPGMA or PCO analyses. Our results suggest that this species may have had different refugia during the glacial stages in the southern region and that the northern variety may have multiple origins from these different refugia.

Stability of alpine meadow ecosystem on the Qinghai-Tibetan Plateau

The meadow ecosystem on the Qinghai-Tibetan Plateau is considered to be sensitive to climate change. An understanding of the alpine meadow ecosystem is therefore important for predicting the response of ecosystems to climate change. In this study, we use the coefficients of variation (Cv) and stability (E) obtained from the Haibei Alpine Meadow Ecosystem Research Station to characterize the ecosystem stability. The results suggest that the net primary production of the alpine meadow ecosystem was more stable (Cv = 13.18%) than annual precipitation (Cv = 16.55%) and annual mean air temperature (Cv= 28.82%). The net primary production was insensitive to either the precipitation (E = 0.0782) or air temperature (E = 0.1113). In summary, the alpine meadow ecosystem on the Qinghai-Tibetan Plateau is much stable. Comparison of alpine meadow ecosystem stability with other five natural grassland ecosystems in Israel and southern African indicates that the alpine meadow ecosystem on the Qinghai-Tibetan Plateau is the most stable ecosystem. The alpine meadow ecosystem with relatively simple structure has high stability, which indicates that community stability is not only correlated with biodiversity and community complicity but also with environmental stability. An average oscillation cycles of 3-4 years existed in annual precipitation, annual mean air temperature, net primary production and the population size of consumers at the Haibei natural ecosystem. The high stability of the alpine meadow ecosystem may be resulting also from the adaptation of the ecosystem to the alpine environment.

Origin of the Qinghai-Tibetan Plateau endemic Milula (Liliaceae): further insights from karyological comparisons with Allium

Milula, a monotypic genus endemic to the Qinghai-Tibetan Plateau, was found to be nested deeply within Allium by the molecular phylogeny despite the aberrant morphology. It remains unknown what had contributed to the rapid evolution of morphology and origin of this exceptional species. In contrast to a previous report of its karyotypes with 2n = 16 = 8M+8SM (2SAT), similar to most species of Allium, a rather different karyotype, 2n = 20 = 4M +10SM+6T (2SAT), was found in examined 31 individuals from 6 populations of M. spicata distributed in the central Tibet. Karyotypes of 7 Allium species occurring in the Qinghai-Tibetan Plateau were further reported. The basic number x = 8 was confirmed for all of them and their karyotypes consist mainly of metacentric and submetacentric chromosomes with rare subterminal and terminal chromosomes. The karyotype of M. spicata is distinctly different from that of most Allium species occurring in the plateau through a complete comparison of all available species in this region and adjacent areas. However, the same chromosome number and similar karyotypic structure were found in A. fasciculatum of Sect. Bromatorrhiza, indicating a possible close relationship between them. But this similarity is contradictory to the preliminary molecular phylogenetic analysis that Milula was closely related to A. cyathophorum of Sect. Bromatorrhiza with x=8, but the other species with x=10 and 11 in this section were clearly placed in the other clade. We therefore suggested that the paralleling evolution from x=8 to x=9, 10 and 11 with increasing asymmetry of karyotype possibly due to the chromosomal Robertsonian translocation might occur separately in the two recognized phylogenetic lineages of Allium. In addition to aneuploidy and following change of the chromosomal structures, the habitat isolation due to the recent uplift of the Qinghai-Tibetan Plateau and the Quaternary climatic oscillation, plays a greater role in origin of Milula and other endemic species (genera) with aberrant morphology from their progenitors.

一种用于异步电机效率优化的混合搜索方法

针对异步电机效率优化问题,提出了一种混合搜索方法。该方法起始于模糊自适应搜索,然后切换至黄金分割法以获取确定收敛速度。这样的搜索步骤能够降低转矩波动,避免在最优点附近发生振荡。利用一个包含铁损和机械损耗的异步电机模型,对该方法进行了矢量控制下的性能验证。仿真结果验证了该方法的有效性。

Smith预估器在纳米操作系统中的应用

针对纳米操作系统中存在的滞后性质,提出了带Smith预估器的PID控制方法。从理论上解决了纳米操作系统中由于纯滞后性质而引起的系统超调或振荡。从而保证了纳米操作中观测的精度和准确性。

转子系统油膜振荡的小波包分解与频带能量比例特征分析

在介绍小波包分解原理的基础上 ,对试验测得的单盘单跨转子系统的油膜振荡非平稳信号用小波包分解方法进行了研究。采用db44小波基函数进行 4层小波包分解。给出了各频带内分解信号的特点及频带能量比例 ,其中第 3频带是该转子系统在 960 0r/min时产生油膜振荡的特征频带。得到的试验数据及其分析结果对转子系统油膜振荡研究和旋转机械状态监测等具有重要意义。

液压反铲挖掘机工作装置有限元动态分析

针对挖掘机工作装置的结构特点 ,建立了其系统 3维有限元模型 ,通过分析得出了系统的固有频率和振型 ,以及结构动态响应各时刻的变形和应力分布规律。比较了工作装置在开挖难度较小和较大 2种情况下 ,挖掘阻力对系统动力响应的影响

一种具有波浪补偿和防晃功能的船用起重机

船用特种起重机是指在海浪环境中进行起吊作业的专用起重机器人，其需要解决的主要问题是波浪运动的补偿和摆动的抑制。本文针对此采用并串联结构机器人的理论和技术，提出了能实现海浪补偿与防晃一体化的起重机方案，对方案的机构原理、测量方法、控制方法和非对称液压缸的建模等都进行了具体的分析与论证。最后得出结论：现有技术水平能够实现海浪补偿与防晃一体化，研制特种起重机是可行的。本项研究的科学意义在于探索在海浪环境

jλ变换与离散振荡指标Md在海洋机器人控制系统动态综合中的应用

本文给出由连续系统的S域,到离散领域?的转换的通式。并给出了离散传递函数W(?)的基本属性、离散振荡指标Md的基本概念、?变换的稳定判据、二阶无静差数字随动系统的动态综合、Md综合在海洋机器人航向控制中的应用。

特大地震激发地球自由振荡的数值模拟研究

Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.