Satellite-station time synchronization information based real-time orbit error monitoring and correction of navigation satellite in Beidou System

Satellite-station two-way time comparison is a typical design in Beidou System(BDS)which is significantly different from other satellite navigation systems.As a type of two-way time comparison method,BDS time synchronization is hardly influenced by satellite orbit error,atmosphere delay,tracking station coordinate error and measurement model error.Meanwhile,single-way time comparison can be realized through the method of Multi-satellite Precision Orbit Determination(MPOD)with pseudo-range and carrier phase of monitor receiver.It is proved in the constellation of 3GEO/2IGSO that the radial orbit error can be reflected in the difference between two-way time comparison and single-way time comparison,and that may lead to a substitute for orbit evaluation by SLR.In this article,the relation between orbit error and difference of two-way and single-way time comparison is illustrated based on the whole constellation of BDS.Considering the all-weather and real-time operation mode of two-way time comparison,the orbit error could be quantifiably monitored in a real-time mode through comparing two-way and single-way time synchronization.In addition,the orbit error can be predicted and corrected in a short time based on its periodic characteristic.It is described in the experiments of GEO and IGSO that the prediction accuracy of space signal can be obviously improved when the prediction orbit error is sent to the users through navigation message,and then the UERE including terminal error can be reduced from 0.1 m to 0.4 m while the average accuracy can be improved more than 27%.Though it is still hard to make accuracy improvement for Precision Orbit Determination(POD)and orbit prediction because of the confined tracking net and the difficulties in dynamic model optimization,in this paper,a practical method for orbit accuracy improvement is proposed based on two-way time comparison which can result in the reflection of orbit error.     

The method and experiment analysis of two-way common-view satellite time transfer for compass system

Time synchronization between ground and satellites is a key technology for satellite navigation system. With dual-channel satellite, a method called Two-Way Common-View(TWCV) satellite time transfer for Compass system is proposed, which combines both characteristics of satellite common-view and two-way satellite-ground time transfer. By satellite-ground two-way pseudo-range differencing and two stations common-view differencing, this TWCV method can completely eliminate the influence of common errors, such as satellite clock offset, ephemeris errors, troposphere delay and station coordinates errors. At the same time, ionosphere delay related to signal frequency is also weakened significantly. So the precision of time transfer is improved much more greatly than before. In this paper, the basic principle is introduced in detail, the effect of major errors is analyzed and the practical calculation model in the Earth-fixed coordinate system for this new method is provided. Finally, experiment analysis is conducted with actual Compass observing data. The results show that the deviation and the stability of the satellite dual channel can be better than 0.1 ns, and the accuracy of the two-way common-view satellite time transfer can achieve 0.4 ns. All these results have verified the correctness of this TWCV method and model. In addition, we compare this TWCV satellite time transfer with the independent C-band TWSTFT(Two-Way Satellite Time and Frequency Transfer). It shows that the result of the TWCV satellite time transfer is in accordance with the C-band TWSTFT result, which further suggests that the TWCV method is a remote high precision time transfer technique. The research results in this paper are very important references for the development and application of Compass satellite navigation system.     

风云三号C星全球导航卫星掩星探测仪首次实现北斗掩星探测

全球导航卫星掩星探测仪(GNOS)是国际首台北斗系统(BDS)和全球定位系统(GPS)双系统兼容掩星探测仪, 2013年9月23日随风云三号卫星C星(FY3 C)发射入轨, 目前已测得大量掩星数据. 介绍了FY3 C-GNOS的组成; 统计了2013年10月12日全天的FY3 C-GNOS掩星事件及其全球分布情况; 通过与GPS精密定轨结果进行对比分析, 测试了BDS在轨实时定位精度, 检验了BDS掩星产品的可靠性和一致性. 初步分析结果表明: 14颗BDS卫星在轨运营的条件下, BDS和GPS兼容掩星探测仪可将掩星事件数提高约33.3%; BDS实时定位平均偏差优于6 m, 标准偏差优于7 m; 5-25 km高度范围内, BDS与GPS内符合精度大气折射率优于2%, 温度优于2 K, 湿度优于1.5 g/kg, 压强优于2%, 电离层峰值密度优于15.6%. GNOS的在轨正常运行及BDS与GPS掩星定位精度及反演产品的一致性为北斗掩星探测的业务化运行奠定了基础.     

Orbit determination and prediction for Beidou GEO satellites at the time of the spring/autumn equinox

Geostationary(GEO) satellites form an indispensable component of the constellation of Beidou navigation system(BDS). The ephemerides, or predicted orbits of these GEO satellites(GEOs), are broadcast to positioning, navigation, and timing users. User equivalent ranging error(UERE) based on broadcast message is better than 1.5 m(root formal errors: RMS) for GEO satellites. However, monitoring of UERE indicates that the orbital prediction precision is significantly degraded when the Sun is close to the Earth's equatorial plane(or near spring or autumn Equinox). Error source analysis shows that the complicated solar radiation pressure on satellite buses and the simple box-wing model maybe the major contributor to the deterioration of orbital precision. With the aid of BDS' two-way frequency and time transfer between the GEOs and Beidou time(BDT, that is maintained at the master control station), we propose a new orbit determination strategy, namely three-step approach of the multi-satellite precise orbit determination(MPOD). Pseudo-range(carrier phase) data are transformed to geometric range(biased geometric range) data without clock offsets; and reasonable empirical acceleration parameters are estimated along with orbital elements to account for the error in solar radiation pressure modeling. Experiments with Beidou data show that using the proposed approach, the GEOs' UERE when near the autumn Equinox of 2012 can be improved to 1.3 m from 2.5 m(RMS), and the probability of user equivalent range error(UERE)2.0 m can be improved from 50% to above 85%.     

基于STK的BDS星座仿真和性能分析

利用卫星仿真工具包(STK)设计和建立了北斗卫星导航系统(BDS)的全球星座。在此基础上,选取南北半球不同经纬度的11个城市作为观测点,分别获得了BDS和全球定位系统(GPS)的可见星数和几何精度衰减因子(GDOP)的值,并按照5°×5°的空间分辨率,仿真了全球范围内BDS和GPS的可见星数和GDOP值。仿真结果表明,全球星座的BDS具有与GPS同样优良的导航性能,其在亚太地区的性能甚至更优于GPS,完全可以为用户提供高性能的满足所需导航性能(RNP)的导航服务。     

基于GPS/北斗共视技术的铷钟驯服方法

在分析了现有铷钟驯服方法的基础上,提出了一种基于GPS/北斗共视技术的铷钟驯服方法,将铷钟的输出频率驯服到时间频率的最高国家基准UTC(NIM)上。试验结果表明,本文提出的方法可以有效实现对铷钟的驯服,经驯服的铷钟频率信号指标良好,能够提供铯钟量级的高精度频率标准。     

GNSS clock corrections densification at SHAO: from 5 min to 30 s

High frequency multi-GNSS zero-difference applications like Precise Orbit Determination(POD)for Low Earth Orbiters(LEO)and high frequency kinematic positioning require corresponding high-rate GNSS clock corrections.The determination of the GNSS clocks in the orbit determination process is time consuming,especially in the combined GPS/GLONASS processing.At present,a large number of IGS Analysis Centers(AC)provide clock corrections in 5-min sampling and only a few ACs provide clocks in 30-s sampling for both GPS and GLONASS.In this paper,an efficient epoch-difference GNSS clock determination algorithm is adopted based on the algorithm used by the Center for Orbit Determination in Europe(CODE).The clock determination procedure of the GNSS Analysis Center at Shanghai Astronomical Observatory(SHAO)and the algorithm is described in detail.It is shown that the approach greatly speeds up the processing,and the densified 30-s clocks have the same quality as the 5-min clocks estimated based on a zero-difference solution.Comparing the densified 30-s GNSS clocks provided by SHAO with that of IGS and its ACs,results show that our 30-s GNSS clocks are of the same quality as that of IGS.Allan deviation also gives the same conclusion.Further validation of the SHAO 30-s clock product is performed in kinematic PPP and LEO POD.Results indicate that the positions have the same accuracy when using SHAO 30-s GNSS clocks or IGS(and its AC)finals.The robustness of the algorithm and processing approach ensure its extension to provide clocks in 5-s or even higher frequencies.The implementation of the new approach is simple and it could be delivered as a black-box to the current scientific software packages.     

Orbit determination for geostationary satellites with the combination of transfer ranging and pseudorange data

Geostationary satellites(GEOs) play a significant role in the regional satellite navigation system.Simulation experiments show that the clock corrections could be mitigated through a single strategy or double differencing strategies for a navigation constellation,but for the mode of individual GEO orbit determination,high precision orbit and clock correction could not be obtained in the orbit determination based on the pseudorange data.A new GEO combined precise orbit determination(POD) strategy is studied in this paper,which combines pseudorange data and C-band transfer ranging data.This strategy overcomes the deficiency of C-band transfer ranging caused by limited stations and tracking time available.With the combination of transfer ranging and pseudorange data,clock corrections between the GEO and the stations can be estimated simultaneously along with orbital parameters,maintaining self-consistency between the satellite ephemeris and clock correction parameters.The error covariance analysis is conducted to illuminate the contributions from the transfer ranging data and the psudoranging data.Using data collected for a Chinese GEO satellite with 3 C-band transfer ranging stations and 4 L-band pseudorange tracking stations,POD experiments indicate that a meter-level accuracy is achievable.The root-mean-square(RMS) of the post-fit C-band ranging data is about 0.203 m,and the RMS of the post-fit pseudorange is 0.408 m.Radial component errors of the POD experiments are independently evaluated with the satellite laser ranging(SLR) data from a station in Beijing,with the residual RMS of 0.076 m.The SLR evaluation also suggests that for 2-h orbital predication,the predicted radial error is about 0.404 m,and the clock correction error is about 1.38 ns.Even for the combination of one C-band transfer ranging station and 4 pseudorange stations,POD is able to achieve a reasonable accuracy with the radial error of 0.280 m and the 2-h predicted radial error of 0.888 m.Clock synchronization between the GEO and tracking stations is achieved with an estimated accuracy of about 1.55 ns,meeting the navigation service requirements.     

Orbit determination and time synchronization for a GEO/IGSO satellite navigation constellation with regional tracking network

Aiming at regional services,the space segment of COMPASS (Phase I) satellite navigation system is a constellation of Geostationary Earth Orbit (GEO),Inclined Geostationary Earth Orbit (IGSO) and Medium Earth Orbit (MEO) satellites.Precise orbit determination (POD) for the satellites is limited by the geographic distribution of regional tracking stations.Independent time synchronization (TS) system is developed to supplement the regional tracking network,and satellite clock errors and orbit data may be obtai...     

Analysis of RDSS positioning accuracy based on RNSS wide area differential technique

The BeiDou Navigation Satellite System(BDS) provides Radio Navigation Service System(RNSS) as well as Radio Determination Service System(RDSS).RDSS users can obtain positioning by responding the Master Control Center(MCC) inquiries to signal transmitted via GEO satellite transponder.The positioning result can be calculated with elevation constraint by MCC.The primary error sources affecting the RDSS positioning accuracy are the RDSS signal transceiver delay,atmospheric trans-mission delay and GEO satellite position error.During GEO orbit maneuver,poor orbit forecast accuracy significantly impacts RDSS services.A real-time 3-D orbital correction method based on wide-area differential technique is raised to correct the orbital error.Results from the observation shows that the method can successfully improve positioning precision during orbital maneuver,independent from the RDSS reference station.This improvement can reach 50% in maximum.Accurate calibration of the RDSS signal transceiver delay precision and digital elevation map may have a critical role in high precise RDSS positioning services.     

Satellite virtual atomic clock with pseudorange difference function

Satellite atomic clocks are the basis of GPS for the control of time and frequency of navigation signals. In the Chinese Area Positioning System (CAPS), a satellite navigation system without the satellite atomic clocks onboard is successfully developed. Thus, the method of time synchronization based on satellite atomic clocks in GPS is not suitable. Satellite virtual atomic clocks are used to implement satellite navigation. With the satellite virtual atomic clocks, the time at which the signals are transmitted from the ground can be delayed into the time that the signals are transmitted from the satellites and the pseudorange measuring can be fulfilled as in GPS. Satellite virtual atomic clocks can implement the navigation, make a pseudorange difference, remove the ephemeris error, and improve the accuracy of navigation positioning. They not only provide a navigation system without satellite clocks, but also a navigation system with pseudorange difference. Supported by the National Basic Research Program of China (Grant No. 2007CB815502) and the National High Technology Research and Development Program of China (Grant No. 2007AA12Z300)     

A new method for determination of satellite orbits by transfer

The original idea of a new method for determination of satellite orbits by transfer is from Two-Way Satellite Time and Frequency Transfer (TWSTFT). The original method is called “determination of satellite orbit by transfer”. The method is not only for determination of satellite orbit but also for the time transfer with high accuracy and precision. The advantage is that the accuracy and the precision for determination of satellite orbit are very high and the new method is favorable for various applications. The combination of various signals disseminated and received forms various modes of satellite orbit determinations. If receivers at stations receive the own station-disseminated signals via a satellite transponder, it forms an orbit determination mode called “receiving the own station-disseminated signals mode”. If receivers at all stations receive the signals disseminated from the master station via satellite transponders, it forms an orbit determination mode called “receiving the master station-disseminated signals mode”. If all of receivers at stations receive all stations-disseminated signals via satellite transponders, it forms an orbit determination mode called “receiving all stations-disseminated signals mode”. Also there are other combinations of signals for satellite orbit determination. For different orbit determination modes with different signal combinations, their rigorous formulae of processing are hereby presented in this paper. The accurate and the precise satellite orbit determination for both of the modes, “receiving the own station-disseminated signals mode” and “receiving the master station-disseminated signals mode” is attempted. It shows that the accuracy and precision for both of modes are nearly the same, the ranging accuracy is better than 1 cm, and the observation residuals of satellite orbit determination are better than 9 cm in the observation duration of 1 day. Supported by the National Basic Research and Development Program of China (Grant No. 2007CB815503100453001)     

SIS accuracy and service performance of the BDS-3 basic system

On December 27,2018,the basic system of the third-generation BeiDou navigation satellite system(BDS-3)completed the deployment of its constellation of 18 MEO networking satellites as well as the construction of the operation control system(OCS)and began to provide basic navigation services to users worldwide.Compared with BDS-2,BDS-3 aims to offer users better navigation signals and higher precision with a series of new technologies.For example,the spaceborne atomic clock of BDS-3 is upgraded for higher performance,the Ka-band inter-satellite link is adopted for inter-satellite ranging and communication,and new B1C and B2a signals are broadcast in addition to B1I and B3I signals(compatible with BDS-2).In addition,a 9-parameter model based on a spherical harmonic function is employed for ionospheric delay corrections.Using the observation data from 18 satellites of the basic system,this paper conducts a comprehensive evaluation of the pseudorange measurement characteristics,signal-in-space(SIS)accuracy of navigation messages and global service capability of BDS-3.The results indicate that the pseudorange measurement multipath effect and observation noise of BDS-3 satellites are better than those of BDS-2;additionally,with the support of inter-satellite links,the user range error(URE)of the BDS-3 satellite broadcast ephemeris is better than 10 cm,the precision of the broadcast clock parameter is better than 1.5 ns,and the SIS accuracy is better than 0.6 m overall.Different from the traditional Klobuchar model,the BeiDou global broadcast ionospheric delay correction model(BDGIM)can provide ionospheric delay corrections better than 70%for worldwide single-frequency users.The service capability evaluation of the basic system consists mainly of the accuracy improvement of the B1I and B3I signals according to BDS-2 as well as the global positioning accuracy of the new signals.These results prove that the BDS-3 basic system has achieved the design goal;that is,both the horizontal and the vertical global positioning accuracies are better than 10 m(95%).In the future,6 MEO satellites as well as 3 GEO satellites and 3 IGSO satellites for regional enhancement purposes will be deployed for full operation;consequently,BDS-3 will definitely provide a higher SIS accuracy and better service capability.     

COMPASS time synchronization and dissemination—Toward centimetre positioning accuracy

In this paper we investigate methods to achieve highly accurate time synchronization among the satellites of the COMPASS global navigation satellite system (GNSS). Owing to the special design of COMPASS which implements several geo-stationary satellites (GEO), time synchronization can be highly accurate via microwave links between ground stations to the GEO satellites. Serving as space-borne relay stations, the GEO satellites can further disseminate time and frequency signals to other satellites such as the inclined geo-synchronous (IGSO) and mid-earth orbit (MEO) satellites within the system. It is shown that, because of the accuracy in clock synchronization, the theoretical accuracy of COMPASS positioning and navigation will surpass that of the GPS. In addition, the COMPASS system can function with its entire positioning, navigation, and time-dissemination services even without the ground link, thus making it much more robust and secure. We further show that time dissemination using the COMPASS-GEO satellites to earth-fixed stations can achieve very high accuracy, to reach 100 ps in time dissemination and 3 cm in positioning accuracy, respectively. In this paper, we also analyze two feasible synchronization plans. All special and general relativistic effects related to COMPASS clocks frequency and time shifts are given. We conclude that COMPASS can reach centimeter-level positioning accuracy and discuss potential applications.     

Orbit determination and prediction accuracy analysis for a regional tracking network

China's COMPASS satellite navigation system relies on a regional tracking network to provide navigation services. Limited by its geographic border,the regional network is able to cover only 30% of the medium-earth-orbits(MEO). Accuracy of determined and predicted orbits is not able to satisfy system requirements if the tracking data processing strategy for global tracking network processing is used for the regional network. Two major error sources for orbital prediction are accuracy of initial orbital elements and dynamical modeling. To achieve better prediction accuracy,we propose a two-step orbit determination and prediction strategy. For step 1,only solar radiation pressure(SRP) parameters are estimated along with the orbital elements and other parameters; for step 2,all parameters are estimated but the SRP parameters are tightly constrained to their step 1 estimates. Experimenting with data from a regional GPS network,we conclude for orbital prediction using the proposed two-step strategy,the average user range error(URE) for 24-h prediction arcs is better than 0.6 m.     

The principle of the positioning system based on communication satellites

It is a long dream to realize the communication and navigation functionality in a satellite system in the world. This paper introduces how to establish the system, a positioning system based on communication satellites called Chinese Area Positioning System (CAPS). Instead of the typical navigation satellites, the communication satellites are configured firstly to transfer navigation signals from ground stations, and can be used to obtain service of the positioning, velocity and time, and to achieve the function of navigation and positioning. Some key technique issues should be first solved; they include the accuracy position determination and orbit prediction of the communication satellites, the measuring and calculation of transfer time of the signals, the carrier frequency drift in communication satellite signal transfer, how to improve the geometrical configuration of the constellation in the system, and the integration of navigation & communication. Several innovative methods are developed to make the new system have full functions of navigation and communication. Based on the development of crucial techniques and methods, the CAPS demonstration system has been designed and developed. Four communication satellites in the geosynchronous orbit (GEO) located at 87.5°E, 110.5°E, 134°E, 142°E and barometric altimetry are used in the CAPS system. The GEO satellites located at 134°E and 142°E are decommissioned GEO (DGEO) satellites. C-band is used as the navigation band. Dual frequency at C1=4143.15 MHz and C2=3826.02 MHz as well as dual codes with standard code (CA code and precision code (P code)) are adopted. The ground segment consists of five ground stations; the master station is in Lintong, Xi’an. The ground stations take a lot of responsibilities, including monitor and management of the operation of all system components, determination of the satellite position and prediction of the satellite orbit, accomplishment of the virtual atomic clock measurement, transmission and receiving navigation signals to and from each satellite. In the north, the south, the east, the west and the center of Chinese main land, the function of CAPS demonstration system is checked and measured. In cars and on board the system is also checked and measured. The results are as follow: CA-code, horizontal positioning accuracy, 15–25 m (1 σ), vertical, 1–3 m; P-code, horizontal positioning accuracy, 8–10 m (1 σ), vertical, 1–3 m; velocity accuracy, CA-code, 0.13–0.30 m/s, P-code, 0.15–0.17 m/s; time accuracy, CA-code, 160 ns, P-code, 13 ns; determination accuracy of orbit ≤2 m. About 20 million US $ and two years are spent for the development of demonstration. A complete CAPS system is now being established. Supported by the National Natural Science Foundation of China (Grant No. 10453001), the National Basic Research Program of China (Grant No. 2007CB815500), the National High Technology Research and Development Program of China (Grant No. 2004AA105030), and the Funds of the Chinese Academy of Sciences for Key Topics in Innovation Engineering (Grant No. KGCXI-21)     

Methods of rapid orbit forecasting after maneuvers for geostationary satellites

A geostationary (GEO) satellite may serve as a navigation satellite, but there is a problem that maneuvers frequently occur and the forces are difficult to model. Based on the technique of determining satellite orbits by transfer, a predicted orbit with high accuracy may be achieved by the method of statistical orbit determination in case of no maneuver force. The predicted orbit will soon be invalid after the maneuver starts, and it takes a long time to get a valid orbit after the maneuver ends. In order to improve ephemeris usability, the method of rapid orbit forecasting after maneuvers is studied. First, GEO satellite movement is analyzed in case of maneuvers based on the observation from the orbit measurement system by transfer. Then when a GEO satellite is in the free status just after maneuvers, the short arc observation is used to forecast the orbit. It is assumed that the common system bias and biases of each station are constant, which can be obtained from orbit determination with long arc observations. In this way, only 6 orbit elements would be solved by the method of statistical orbit determination, and the ephemeris with high accuracy may be soon obtained. Actual orbit forecasting with short arc observation for SINOSAT-1 satellite shows that, with the tracking network available, the precision of the predicted orbit (RMS of O-C) can reach about 5 m with 15 min arc observation, and about 3 m with 30 min arc observation. Supported by the National High Technology Research and Development Program of China (Grant No. 2006AA12Z322), the National Basic Research Program of China (Grant No. 2007CB815503), and the West Light Program of Chinese Academy of Sciences (Grant No. 2007LH01)     

空间弥散X射线对脉冲星导航精度的影响

基于RXTE卫星天基数据,建立了时空坐标系转换、时间修正及历元折叠方法,构建了用于提取导航信息的Crab脉冲星轮廓,剖析了该天基载荷结构及特性,对卫星运行空间背景辐射进行了模拟计算。结果表明,在设定导航条件下,空间弥散X射线对航天器单星定轨及多星定位影响在km量级以上。同时阐述了实用化脉冲星导航探测中,改进导航定位精度急需注意的技术问题。     

A new method of ionospheric-free hybrid differential positioning based on a double-antenna CAPS receiver

Chinese Area Positioning System (CAPS) is a transmitted satellite navigation system moved by the Chinese Academy of Sciences. Three basic modes of navigation and positioning with CAPS are given, and then a comparative analysis is made in this paper. In terms of the principle that the ionospheric delay is at an inverse ratio to the frequency square, a new ionospheric-free positioning method based on a double-antenna CAPS receiver is put forward. Then the hybrid differential observations and the solving equations and algorithms for one epoch and multi epochs are deduced according to the basic principle of the method. The method may remove the global errors in signal emission, propagation, transmission and receiving (e.g., ionospheric delay, hardware delay, and clock error). So it is very convenient for the single-epoch solution and multi-epoch navigation and positioning, and may efficiently improve the precision of real time CAPS navigation. Furthermore, the method can be used not only for the geometric orbit determination of CAPS GEO and IGSO satellites and the navigation and positioning, but also for the estimation of the tropospheric zenith delay, which is useful for the study of water vapor changes in the atmosphere. Polynomials are used in this method to express the tropospheric zenith delay and CAPS satellite orbits within the limited time interval, which reduces the number of unknown parameters and thus speeds the computation. Supported by the Knowledge Innovation Project of the Chinese Academy of Sciences (Grant No. KGCX1-21), the National Basic Research Program of China (Grant No. 2007CB815500), the National High Technology Research and Development Program of China (Grant No. 2006AA12z303), the National Natural Science Foundation of China (Grant No. 40774009), and the Special Project of Taishan Scholars of Shandong Province of China (Grant No. TSXZ0502)     

An analysis of the wide area differential method of geostationary orbit satellites

This work aims to obtain a wide area differential method for geostationary orbit (GEO) constellation. A comparison between the dilution of precision (DOP) of four-dimensional (4D) calculation including satellite clock errors and ephemeris errors and that of three-dimensional (3D) calculation only including ephemeris errors with the inverse positioning theory of GPS shows the conclusion that all the 3D PDOPs are greatly reduced. Based on this, a basic idea of correcting satellite clock errors and ephemeris errors apart is put forward, and moreover, a specific method of separation is proposed. Satellite clock errors are separated in a master station with time synchronization, and all the remaining pseudo-range errors after the satellite clock errors have been deducted are used to work out ephemeris corrections of all GEO satellites. By a comparative analysis of user positioning accuracy before and after differential, the wide area differential method is verified to be quite valid for GEO constellation. Supported by the National Natural Science Foundation of China (Grant No. 10778715), the National Key Basic Research Development Program of China (Grant No. 2007CB815502), and the Scientific Research Fund of Hunan Provincial Education Department (Grant No. 08B039)