汶川M_S8.0地震河谷地形对汉源县城高烈度异常的影响

汉源县城位于汶川M_S8.0特大地震六度区内唯一的八度异常区,为典型且罕见的远震高烈度异常区.汉源县城处于流沙河的河流阶地之上,河谷地形对地震震害有显著的影响.为定量分析河谷地形对汶川大地震汉源县高烈度异常的影响,基于汉源县城震害科学考察和场地勘察获取的资料,根据震害分布特征和流沙河谷地形的特点,选取1条实测得到的横切汉源县城并垂直流沙河河谷方向的典型剖面作为计算模型,以脉冲作为基底输入,采用有限差分方法研究了该剖面的场地放大效应,分析了地形对高烈度异常的影响.计算结果表明:汉源县城场地对地震动放大效应的显著频段是1.0~10Hz,且这一频段老县城场地的放大效应比新县城场地显著;汉源老县城场地对汶川M_S8.0地震主震的地震动有显著的放大效应,地表峰值加速度大大超过了抗震设计规范的规定值;汉源场地地形放大效应显著频段与汶川M_S8.0地震的能量集中频段基本吻合,汉源老县城建筑物的自振频率恰恰位于该频段,产生共振效应,从而造成更显著的放大效应,这也是汉源震害异常的主要原因之一.由此可见,河谷地形对地表地震动有重要的影响,在工程选址和抗震设计时应考虑其影响.     

汶川大地震汉源县城建筑物震害调查

汶川地震中,距离震中约200公里的汉源县出现了明显烈度异常现象,汉源县城建筑物震害表现尤为明显.为了解释这一现象,在汶川大地震第一阶段科学考察中,相关调查小组对汉源县城建筑物震害开展了相应的调查.调查结果表明,该处烈度异常现象确实存在,烈度异常与该处局部场地条件有一定的联系.     

“5·12”汶川8.0级地震汉源烈度异常机理的初步探讨

本文通过总结分析汉源烈度异常主要特征,初步分析了导致烈度异常的可能原因。初步研究结果表明,“5·12”汶川8．0级地震震源过程及其能量释放空间分布和场地条件是导致汉源烈度异常的主要因素,同时莫霍界面反射波S11断层破裂带的阻震及激震效应等,也可能是汉源烈度的影响因素。探索汉源烈度异常对总结特大地震烈度分布、震害特征和成灾机理具有重要意义。为此作者提出,应在震害特征总结提炼的基础上,考虑震源机制、地壳介质、场地非线性等因素,对汉源烈度异常区进行全面的数值模拟,进一步确定汉源地震烈度异常的机理。     

汶川特大地震汉源烈度异常原因的初步分析

远离震中的汉源县老县城在汶川大地震中遭到了严重震害,是Ⅵ度区内的唯一的Ⅷ度烈度异常区.汉源的地震烈度异常现象引起学术界的广泛关注.在汉源科学考察基础上,介绍了汉源震害情况,收集和整理了汉源县城及附近有关工程场地的强震和土体力学参数资料.通过背后山滑坡体有限元分析结果与汉源附近强震台记录的比较,从传递函数角度分析了背后山滑坡复活对烈度异常的影响,研究了汉源特殊土层结构对地震动的影响.结果表明,背后山滑坡的复活和特殊土层结构的放大作用等场地条件的影响是引起烈度异常的主要原因.所得结论对进一步分析汉源烈度异常的原因具有一定的参考价值.     

土层结构对汉源烈度异常的影响

汶川MS8.0大地震在远离震中近200 km的汉源县县城产生了高烈度异常, 其原因比较复杂．为了分析土层结构对汶川大地震中汉源县老县城高烈度异常的影响, 在汉源县城震害科学考察基础上, 依据地震烈度异常的分布情况, 在背后山滑坡前缘地带布设5个工程地质勘察钻孔, 获得了汉源县老县城场地土层结构资料．在现场测试和室内试验的基础上, 给出了汉源县老县城场地各层土体动力学参数．本文利用汶川地震九襄强震台的强震记录, 结合其台站的场地资料反演给出了汉源县老县城的基岩地震动时程, 作为地震反应分析的基岩输入地震动．在此基础上利用土层地震反应一维等效线性方法对汉源县老县城场地进行了地震反应分析, 并将计算所得结果与Ⅵ度区其它强震台站获得的强震记录进行比较．研究结果表明, 汉源县老县城土层结构对地震动的放大作用导致地表地震动异常, 是汉源县老县城高烈度异常的主要原因之一．      

Spectral element analysis on the characteristics of seismic wave propagation triggered by Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 earthquake

The 2008 Wenchuan earthquake occurred in an active earthquake zone, i.e., Longmenshan tectonic zone. Seismic waves triggered by this earthquake can be used to explore the characteristics of the fault rupture process and the hierarchical structure of the Earth’s interior. 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. The AK135 model is employed as a prototype of our numerical global Earth model. The Earth’s ellipticity, Earth’s medium attenuation, and topography data are taken into consideration. These wave propagation processes are simulated by solving three-dimensional elastic wave governing equations. Three-dimensional visualization of our numerical results displays the profile of the seismic wave propagation. The three-point source, which is proposed from the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture process than the one-point source. We take comparison of synthetic seismograms with observational data recorded at 16 observatory stations. Primary results show that the synthetic seismograms calculated from three-point source agree well with the observations. This can further reveal that the source rupture process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. Supported by National Basic Research Program of China (Grant No. 2004CB418406), National Natural Science Foundation of China (Grant Nos. 40774049 and 40474038), and Computer Network Information Center, Chinese Academy of Sciences (Grant No. INF105-SCE-02-12)     

A new technique for moment tensor inversion with applications to the 2008 Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>8.0 earthquake sequence

A M_S8.0 earthquake occurred in Wenchuan County, Sichuan Province, China, on May 12, 2008, and subsequently, numerous aftershocks followed. We obtained the moment tensor solutions and source time functions (STFs) for the Wenchuan earthquake and its seven larger aftershocks (M_S5.0～6.0) by a new technique of moment tensor inversion using the broadband and long-period seismic waveform data from the Global Seismic Network (GSN). Firstly, the theoretical background and technical flow of the new technique was briefly introduced, and an aftershock of the Wenchuan earthquake sequence was employed to illustrate the real procedure for inverting the moment tensor; secondly, the moment tensor solutions and STFs of the eight events, including the main shock, were presented, and finally, the interpretation of the results was made. The agreement of our results with the GCMT results indicates the new approach is efficient and feasible. By using this approach, not only the moment tensor solution can be obtained but also the STF can be retrieved; the inverted STFs indicate that the source rupture process may be complicated even for the moderate earthquakes. The inverted focal mechanisms of the Wenchuan earthquake sequence show that the most of the aftershocks occurred in the main faults of the Longmenshan fault zone with predominantly thrustingwith minor right-lateral strike-slip component, but some of them may have occurred in the subfaults with strike-slip faulting in the vicinity of the main faults.     

Stochastic finite-fault simulation of ground motion from the August 11, 2012, <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 6.4 Ahar earthquake,northwestern Iran

In this study, the 11 August 2012 M _w 6.4 Ahar earthquake is investigated using the ground motion simulation based on the stochastic finite-fault model. The earthquake occurred in northwestern Iran and causing extensive damage in the city of Ahar and surrounding areas. A network consisting of 58 acceleration stations recorded the earthquake within 8–217 km of the epicenter. Strong ground motion records from six significant well-recorded stations close to the epicenter have been simulated. These stations are installed in areas which experienced significant structural damage and humanity loss during the earthquake. The simulation is carried out using the dynamic corner frequency model of rupture propagation by extended fault simulation program (EXSIM). For this purpose, the propagation features of shear-wave including \( {Q}_s \) value, kappa value \( {k}_0 \), and soil amplification coefficients at each site are required. The kappa values are obtained from the slope of smoothed amplitude of Fourier spectra of acceleration at higher frequencies. The determined kappa values for vertical and horizontal components are 0.02 and 0.05 s, respectively. Furthermore, an anelastic attenuation parameter is derived from energy decay of a seismic wave by using continuous wavelet transform (CWT) for each station. The average frequency-dependent relation estimated for the region is \( Q=\left(122\pm 38\right){f}^{\left(1.40\pm 0.16\right)}. \) Moreover, the horizontal to vertical spectral ratio \( H/V \) is applied to estimate the site effects at stations. Spectral analysis of the data indicates that the best match between the observed and simulated spectra occurs for an average stress drop of 70 bars. Finally, the simulated and observed results are compared with pseudo acceleration spectra and peak ground motions. The comparison of time series spectra shows good agreement between the observed and the simulated waveforms at frequencies of engineering interest.     

Detecting changes in long-period site responses after the <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 7.6 Chi-Chi earthquake,Taiwan, using strong motion records

Temporal changes in site effects are obtained using the HVSR(horizontal-to-vertical spectral ratio) method and strong motion records after the M w 7.6 Chi-Chi earthquake, Taiwan. Seismic data recorded between 1995 and 2010 are used, comprising 3,708 data from 15 stations adjacent to the Chelungpu fault. Temporal fl uctuations are determined by analyzing the site effect variation using a time–frequency variation(TFV) diagram based on these seismic data. Stations adjacent to the fault show signifi cant disturbances in the resonance frequency at 16–26 Hz. Station TCU129 shows a 40% drop in fundamental frequency after the main shock, and a gradual return to the original state over nine years. For stations located farther from the fault zone, sudden changes in tectonic stress play a dominant role in temporal changes to the HVSR. An impact analysis of the directional factor confi rms our fi nding that the proximity of the fault to seismic stations has the most infl uence on data.     

Study on relationship between features of strain anomaly obtained by borehole strain meter at Wushi station in Xinjiang and strong earthquake: Taking Jiashi <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>=6.8 earthquake as an example

The borehole strain meter at Wushi station recorded obvious anomaly before Jiashi M _s=6.8 earthquake occurred on February 24, 2003. Its features are as follows. (1) Anomaly types are complete. The trend anomaly, short-term anomaly, short-imminent anomaly and exponential anomaly appeared 19 months, 56 days, 4 days, and more than 1 month before the event, respectively; (2) Anomaly is large in magnitude. The maximal magnitude of strain anomaly is 1.7×10^?5, which is rare in the past 20-year observation records at Wushi station; (3) Strain rate fluctuates sharply with obvious alternation of tension and compression. According to the magnitude of strain anomaly, time of exponential anomaly appearance and regional features of recorded anomaly, we could predict the magnitude, occurrence time and potential region to a certain degree.     

Joint inversion of GNSS and teleseismic data for the rupture process of the 2017 <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript>6.5 Jiuzhaigou,China, earthquake

The spatio-temporal slip distribution of the earthquake that occurred on 8 August 2017 in Jiuzhaigou, China, was estimated from the teleseismic body wave and near-field Global Navigation Satellite System (GNSS) data (coseismic displacements and high-rate GPS data) based on a finite fault model. Compared with the inversion results from the teleseismic body waves, the near-field GNSS data can better restrain the rupture area, the maximum slip, the source time function, and the surface rupture. The results show that the maximum slip of the earthquake approaches 1.4 m, the scalar seismic moment is ~ 8.0 × 10¹⁸ N·m (M_w?≈?6.5), and the centroid depth is ~ 15 km. The slip is mainly driven by the left-lateral strike-slip and it is initially inferred that the seismogenic fault occurs in the south branch of the Tazang fault or an undetectable fault, a NW-trending left-lateral strike-slip fault, and belongs to one of the tail structures at the easternmost end of the eastern Kunlun fault zone. The earthquake rupture is mainly concentrated at depths of 5–15 km, which results in the complete rupture of the seismic gap left by the previous four earthquakes with magnitudes >?6.0 in 1973 and 1976. Therefore, the possibility of a strong aftershock on the Huya fault is low. The source duration is ~ 30 s and there are two major ruptures. The main rupture occurs in the first 10 s, 4 s after the earthquake; the second rupture peak arrives in ~ 17 s. In addition, the Coulomb stress study shows that the epicenter of the earthquake is located in the area where the static Coulomb stress change increased because of the 12 May 2017 M_w7.9 Wenchuan, China, earthquake. Therefore, the Wenchuan earthquake promoted the occurrence of the 8 August 2017 Jiuzhaigou earthquake.     

Microtremor Measurements for the Microzonation of Dinar

v--vThe geotechnical site conditions in Dinar town located in western Turkey were investigated after the 1995 Dinar earthquake based on borings, in situ penetration tests, seismic wave velocity measurements, and microtremor records. The variation of damage distribution within the town was evaluated with respect to 23 district damage ratios calculated, based on the detailed damage survey conducted by the General Directorate of Disaster Affairs. Site amplifications were estimated from microtremor spectral ratios and microzonation was performed using a GIS methodology. The results of in situ penetration tests and seismic wave velocity measurements as well as the damage distribution were compared with the amplification zonation obtained from microtremor records. The results indicate the applicability of microtremor spectral ratios for assessing the local site conditions and site amplifications.     

Distribution and characteristics of gravelly soil liquefaction in the Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript> 8.0 earthquake

In this paper,a distribution map of gravelly soil liquefaction that was caused by the Wenchuan M_s 8.0 earthquake in China is proposed based on a detailed field investigation and an analysis of geological soil profiles. The geological background of the earthquake disaster region is summarized by compiling geological cross sections and borehole logs. Meanwhile,four typical liquefied sites were selected to conduct sample drillings,dynamic penetration tests (DPT),and shear wave velocity tests,to understand the features of liquefied gravelly soil. One hundred and eighteen (118) liquefied sites were investigated shortly after the earthquake. The field investigation showed:(1) sandboils and waterspouts occurred extensively,involving thousands of miles of farmland,120 villages,eight schools and five factories,which caused damage to some rural houses,schools,manufacturing facilities and wells,etc.; (2) the Chengdu plain is covered by a gravelly soil layer with a thickness of 0 m to 541 m according to the geological cross sections; (3) there were 80 gravelly soil liquefied sites in the Chengdu plain,shaped as five belt areas that varied from 20 km to 40 km in length,and about ten gravelly soil liquefied sites distributed within Mianyang area; and (4) the grain sizes of the sampled soil were relative larger than the ejected soil on the ground,thus the type of liquefied soil cannot be determined by the ejected soil. The gravelly soil liquefied sites are helpful in enriching the global database of gravelly soil liquefaction and developing a corresponding evaluation method in further research efforts.     

Damage and Seismic Intensity in the Meizoseismal Region of the M_S8.0 Wenchuan Earthquake,Sichuan,China

This paper presents the damage in the meizoseismal region of the Ms8.0 Wenchuan earthquake, Sichuan, China, and the seismic intensifies determined according to "the Chinese Seismic Intensity Scale", and discusses briefly the types of earthquake-generating faults and some features of seismic damage.     

Source process of the 14 November 2001 western Kunlun Mountain <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 earthquake

Based on digital teleseismic P-wave seismograms recorded by 28 long-period seismograph stations of the global seismic network, source process of the November 14, 2001 western Kunlun Mountain M _S=8.1 (M _W=7.8) earthquake is estimated by a new inversion method. The result shows that the earthquake is a very complex rupture event. The source rupture initiated at the hypocenter (35.95°N, 90.54°E, focal depth 10 km, by USGS NEIC), and propagated to the west at first. Then, in several minutes to a hundred minutes and over a large spatial range, several rupture growth points emerged in succession at the eastern end and in the central part of the finite fault. And then the source rupture propagated from these rupture growth points successively and, finally, stopped in the area within 50 km to the east of the centroid position (35.80°N, 92.91°E, focal depth 15 km, by Harvard CMT). The entire rupture lasted for 142 s, and the source process could be roughly separated into three stages: The first stage started at the 0 s and ended at the 52 s, lasting for 52 s and releasing approximately 24.4% of the total moment; The second stage started at the 55 s and ended at the 113 s, lasting for 58 s and releasing approximately 56.5% of the total moment; The third stage started at the 122 s and ended at the 142 s, lasting for 20 s and releasing approximately 19.1% of the total moment. The length of the ruptured fault plane is about 490 km. The maximum width of the ruptured fault plane is about 45 km. The rupture mainly occurred within 30 km in depth under the surface of the Earth. The average static slip in the underground rocky crust is about 1.2 m with the maximum static slip 3.6 m. The average static stress drop is about 5 MPa with the maximum static stress drop 18 MPa. The maximum static slip and the maximum stress drop occurred in an area within 50 km to the east of the centroid position.     

The seismotectonic and macroseismic features of the Wenchuan (Sichuan) earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">S</Emphasis></Subscript> = 8.0) of May 12, 2008

A disastrous earthquake with a magnitude M _S = 8.0 (M _W = 7.9), in China called “the 5.12 Wenchuan earthquake,” occurred on May 12, 2008, in Sichuan province on the border between the Sino-Tibetan Mountains and the Sichuan depression. The instrumental epicenter was registered in the southeastern part of Wenchuan county, and the hypocenter depth was 14 km. As the strongest and most destructive earthquake within mainland China, it caused numerous human losses and destruction of buildings and infrastructure. The seismic effect from the main shock and aftershocks was felt in many counties, towns, and villages, though Sichuan province suffered the most. The maximum intensity of the shocks was estimated at 11 degrees, according to the Chinese macroseismic scale. In the process of source opening, from the southern part of Wenchuan county to the vicinities of Quingchuan, a seismic fault system with a total length up to 240 km out-cropped on the earth’s surface, confined to the Longmenshan fault belt. The seismic fault system disturbed the original ground, resulting in the collapse or damage to various constructions, such as buildings, homes, bridges, roads, etc. Fault offsets had a dextral strike-slip and thrust kinematic combination. The earthquake generated several tens of thousands of landslides, rockfalls, and debris flows. Many dammed ponds appeared in the epicentral zone due to the activation of landslides. Thus, the geological effects turned out to be the most destructive factor in this case. At the same time, the seismic intensity of surface shaking was abnormally low even in direct proximity to the seismic fault system. Usually it was no more than 7–8 degrees. This macroseismic phenomenon may turn out to be rather typical for many major earthquakes.     

The response law of far-field seismic ground motion of the Wenchuan earthquake and its damaging mechanism in the Loess Plateau

A series of housing collapses and other serious damage was caused by the 2008 Wenchuan M_S 8.0 earthquake in the seismic intensity Ⅵ areas of the Loess Plateau, which is hundreds of kilometers away from the epicenter, and which showed a remarkable seismic intensity anomaly. The seismic disasters are closely related to the seismic response characteristics of the site, therefore, the systematic study of the far-field seismic response law of the Wenchuan earthquake in the Loess Plateau is of great significance to prevent the far-field disaster of great earthquake. In this paper, the seismic acceleration records of several bedrock stations and loess stations from the seismogenic fault of the Wenchuan earthquake to the Loess Plateau were collected, and the attenuation law of ground motion along the propagation path and the characteristics of seismic response on the loess site are studied, and the mechanism of amplification effect of ground motion is analyzed based on the dynamic feature parameters of the loess site obtained through the HVSR method. Taking a typical loess site of thick deposit as the prototype, a series of shaking table tests of dynamic response of loess site models with different thicknesses were carried out. Amplification effect, spectral characteristics of acceleration in model sites were analyzed under the action of a far-field seismic wave of the Wenchuan earthquake. The results show that seismic attenuation on the propagation path along the NE strike of the seismogenic fault to the Loess Plateau is slower than that in other directions, and the predominant period range of ground motion on bedrock site of the Loess Plateau presents broadband characteristics. Because the natural periods of loess sites with thick deposits are within the predominant period range of bedrock input wave, loess sites appear significant amplification effect of ground motion, the horizontal acceleration of ground motion exceeds 0.1 ?g, the seismic intensity reaches 7°. The thicker the loess deposit is, the more significant the change of spectral characteristics of ground motion on loess sites, and the narrower the predominant period range of ground motion becomes, and the closer it is to the natural period of loess sites. Therefore, for some old houses on thick loess sites, the poor seismic performance and strong seismic response eventually led to their collapses and damages because their natural periods are very close to the predominant period of ground motion of the Wenchuan earthquake on thick loess sites; For these damaged high-rise buildings, the resonance effect might be the main reason for their damages because their natural periods are included in the predominant period range of ground motion of the Wenchuan earthquake on thick loess sites.These research results would provide a basis for seismic disasters prediction and evaluation and seismic design of construction engineering in the Loess Plateau.     

3D <Emphasis Type="Italic">v</Emphasis><Subscript>P</Subscript> and <Emphasis Type="Italic">v</Emphasis><Subscript>S</Subscript> models of southeastern margin of the Tibetan plateau from joint inversion of body-wave arrival times and surface-wave dispersion data

A new 3D velocity model of the crust and upper mantle in the southeastern (SE) margin of the Tibetan plateau was obtained by joint inversion of body- and surface-wave data. For the body-wave data, we used 7190 events recorded by 102 stations in the SE margin of the Tibetan plateau. The surface-wave data consist of Rayleigh wave phase velocity dispersion curves obtained from ambient noise cross-correlation analysis recorded by a dense array in the SE margin of the Tibetan plateau. The joint inversion clearly improves the v _S model because it is constrained by both data types. The results show that at around 10 km depth there are two low-velocity anomalies embedded within three high-velocity bodies along the Longmenshan fault system. These high-velocity bodies correspond well with the Precambrian massifs, and the two located to the northeast of 2013 M _S 7.0 Lushan earthquake are associated with high fault slip areas during the 2008 Wenchuan earthquake. The aftershock gap between 2013 Lushan earthquake and 2008 Wenchuan earthquake is associated with low-velocity anomalies, which also acts as a barrier zone for ruptures of two earthquakes. Generally large earthquakes (M ≥ 5) in the region occurring from 2008 to 2015 are located around the high-velocity zones, indicating that they may act as asperities for these large earthquakes. Joint inversion results also clearly show that there exist low-velocity or weak zones in the mid-lower crust, which are not evenly distributed beneath the SE margin of Tibetan plateau.     

汶川8．0级地震前电磁辐射异常变化特征初步研究

汶川8．0级地震前,距离震中只有35km的金河电磁波台记录到了宝贵的电磁扰动信息。本文分析了汶川5．12大地震前后该台的电磁波异常变化特征,表明：①强震前确实有电磁辐射异常存在,震前存在丰富的电磁扰动信息;②电磁扰动前兆异常在时间进程上有明显的阶段性变化特征,这可能是地震孕育不同阶段特征的反映。     

Numerical simulation of strong ground motion for the <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 Wenchuan earthquake of 12 May 2008

The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morphological factors affecting the seismic hazard, we simulated the strong ground motion caused by the earthquake, incorporating three-dimensional (3D) earth structure, finite-fault rupture, and realistic surface topography. The simulated ground motions reveal that the fault rupture and basin structure control the overall pattern of the peak ground shaking. Large peak ground velocity (PGV) is distributed in two narrow areas: one with the largest PGV values is above the hanging wall of the fault and attributed to the locations of fault asperities and rupture directivity; the other is along the northwestern margin of the Sichuan Basin and caused by both the directivity of fault rupture and the amplification in the thick sediment basin. Rough topography above the rupture fault causes wave scattering, resulting in significantly larger peak ground motion on the apex of topographic relief than in the valley. Topography and scattering also reduce the wave energy in the forward direction of fault rupture but increase the PGV in other parts of the basin. These results suggest the need for a localized hazard assessment in places of rough topography that takes the topographic effects into account. Finally, had the earthquake started at the northeast end of the fault zone and ruptured to the southwest, Chengdu would have suffered a much stronger shaking than it experienced on 12 May, 2008. Supported by the U.S. National Science Foundation (Grant Nos. EAR 0738779 and OCE 0727919), the National Basic Research Program of China (Grant No. 2004CB418404), and partially by the National Nature Science Foundation of China (Grant No. 40521002)