Extinction pattern and process of siliceous sponge spicules in deep-water during the latest Permian in South China

Diverse and abundant siliceous sponge spicules were found in the latest Permian beds, Dongpan and Ma'anying sections, South China, including 52 types and 85 forms. Further investigation on these spicules allows us to understand extinction patterns and processes of deep-water sponges. These sponge spicules rapidly decreased below the Permian/Triassic boundary (PTB), and the extinction rates reach up to 88%-90% for types and 88%-92% for forms. Their extinction pattern is a gradual one that consists of two stages: the first is characterized by a gentle and slow extinction speed and low extinction rate, and the second by sharp and fast extinction speed and high extinction rate. The morphological extinction process is involved in the disappearance first of the triaxons and tetraxons, then of the polyaxons and demas, and last of monaxons. In exterior structure extinction, the complex spicules with branches and spines became extinct more easily than did smooth spicules. After the end-Permian mass extinction, only five common and smooth forms survived: Oxeas A, Oxeas B, Strongles B, Oxy-orthpentactines and Oxy-orthohexactines A.     

Devonian Frasnian-Famennian transitional event deposits of Guangxi, South China and their possible tsunami origin

A suite of event deposits, isochronous and interrelated in origin, occurs widely near the Frasnian- Famennian boundary in Guangxi, South China. It is mainly distributed in the facies areas of the platform-margin slope and inter-platform rift-trough. The rudstone or calcirudite occur in the Liujing section of Hengxian and Luoxiu and Baqi sections of Xiangzhou in the facies area of platform-margin slope. The turbidites are observed in the sections of Yangdi and Baisha of Guilin, Du'an of Debao, Nandong and Sanli of Wuxuan, Xiangtian, Ma'anshan of Xiangzhou, Nayi of Chongzuo, Yunpan of Shanglin in the facies area of inter-platform rift-trough. The massive homogenites occur in sections of Mangchang, Luofu and Road from Nandan to Tian'e and Ma'anshan etc in the facies area of inter-platform rift-trough. Herein event deposits can be correlated in stratohorizon to the turbidite in the bottom of the Lower triangularis Zone in Hony railroad cut of Belgium, Devils Gate of Nevada, USA, Atrous of Morocco, South Urals and Fore-Kolyma of northeastern Siberia of Russia, and erosional discontinuities and brecciation in the Frasnian-Famennian boundary in South Polish-Moravian shelf. The event deposits could be caused by a violent tsunami related to bolide impacts into ocean.     

Two new species of <Emphasis Type="Italic">Vestrogothia</Emphasis> (Phosphatocopina,Crustacea) of Orsten-type preservation from the Upper Cambrian in western Hunan,South China

Although fossils of Orsten-type preservation represented by Skaracarida and Phosphatocopina were first reported from the Middle and Upper Cambrian in western Hunan, South China in 2005, diversified phosphatocopine species have never been appropriately described and elucidated in terms of their evolutionary relationships. Here, we described two new species of Phosphatocopina, Vestrogothia anterispinata sp. nov. and V. bispinata sp. nov. The evolutionary relationship among all the valid species of Phosphatocopina is discussed and the phylogeny of Phosphatocopina is reestablished using cladistic analysis. Accordingly, phosphatocopines primitively evolved along two lineages: one evolved towards the emergence of lobes; the other evolved towards the emergence of dorsal spines. Respectively, lobes originated independently three times, and dorsal spines originated only once. With regard to the dorsal rims, interdorsum originated first, whereas dorsal furrow originated independently twice on the basis of interdorsum. Probably Hesslandona may represent a polyphyletic group, whereas Vestrogothia a monophyletic group. Supported by National Natural Science Foundation of China (Grant Nos. 40772008, 40572003), Research Fund for Doctoral Program of High Education (Grant No. 20060001059), State Key Laboratory of Palaeobiology and Stratigraphy, Nanjing Institute of Geology and Palaeontology, Chinese Academy of Sciences (Grant No. 083101)     

Observation of baroclinic eddies southeast of Okinawa Island

In the region southeast of Okinawa, during May to July 2001, a cyclonic and an anticyclonic eddy were observed from combined measurements of hydrocasts, an upward-looking moored acoustic Doppler current profiler (MADCP), pressure-recording inverted echo sounders (PIESs), satellite altimetry, and a coastal tide gauge. The hydrographic data showed that the lowest/highest temperature (T) and salinity (S) anomalies from a 13-year mean for the same season were respectively -3.0/ 2.5℃ and -0.20/ 0.15 psu at 380/500 dbar for the cyclonic/anticyclonic eddies. From the PIES data, using a gravest empirical mode method, we estimated time-varying surface dynamic height (D) anomaly referred to 2000 dbar changing from -20 to 30 cm, and time-varying T and S anomalies at 500 dbar ranging through about ±2 ℃ and ±0.2 psu, respectively. The passage of the eddies caused variations of both satellite-measured sea surface height anomaly (SSHA) and tide-gauge-measured sea level anomaly to change from about –20 to 30 cm, consistent with the D anomaly from the PIESs. Bottom pressure sensors measured no variation related to these eddy activities, which indicated that the two eddies were dominated by baro-clinicity. Time series of SSHA map confirmed that the two eddies, originating from the North Pacific Subtropical Countercurrent region near 20°―30°N and 150°―160°E, traveled about 3000 km for about 18 months with mean westward propagation speed of about 6 cm/s, before arriving at the region southeast of Okinawa Island.     

Discovery of Cu-Ni-Zn-Sn-Fe intermetallic compounds and S-bearing alloys in the Zhaishang gold deposit, southern Gansu Province and its geological significance

Examination of ores by optical microscope and EPMA from the Zhaishang gold deposit, southern Gansu Province, has revealed an abundance of rare minerals. These include native metals, Cu-Ni-Zn-Sn-Fe polymetallic compounds and S-bearing alloys of Ni, Fe, Zn, Cu and Sn, occurring as native nickel, Zn-Cu alloy, Ni-Zn-Cu alloy, Sn-Zn-Ni-Cu alloy, Zn-Cu-Ni alloy, Zn-Fe-Cu-Sn-Ni alloy, Fe-Ni-S alloy, Sn-Fe-Ni-S alloy, Fe-Zn-Cu-Ni-S alloy, Zn-Ni-Cu-Fe-S alloy and others. Compared with the Zn-Cu alloy minerals discovered previously, these Zn-Cu minerals fall in the α or α β portion in Zn-Cu alloy phase diagram, and the α portion has higher Cu content. Cu-Ni-Zn-Sn-Fe intermetallic compounds and S-bearing alloy minerals have not been previously reported in the literature. These rare alloys formed in a strongly reducing environment with absent oxygen and low sulfur activities.     

Retrospective forecasting test of a statistical physics model for earthquakes in Sichuan-Yunnan region

Pattern informatics (PI) model is one of the recently developed predictive models of earthquake phys- ics based on the statistical mechanics of complex systems. In this paper, retrospective forecast test of the PI model was conducted for the earthquakes in Sichuan-Yunnan region since 1988, exploring the possibility to apply this model to the estimation of time-dependent seismic hazard in continental China. Regional earthquake catalogue down to ML3.0 from 1970 to 2007 was used. The ‘target magnitude’ for the forecast test was MS5.5. Fifteen-year long ‘sliding time window’ was used in the PI calculation, with ‘anomaly training time window’ being 5 years and ‘forecast time window’ being 5 years, respectively. Receiver operating characteristic (ROC) test was conducted for the evaluation of the forecast result, showing that the PI forecast outperforms not only random guess but also the simple number counting approach based on the clustering hypothesis of earthquakes (the RI forecast). If the ‘forecast time window’ was shortened to 3 years and 1 year, respectively, the forecast capability of the PI model de- creased significantly, albeit outperformed random forecast. For the one year ‘forecast time window’, the PI result was almost comparable to the RI result, indicating that clustering properties play a more important role at this time scale.     

New discovery of Palaeontinid fossils from the Middle Jurassic in Daohugou, Inner Mongolia (Homoptera, Palaeontinidae)

Three new species of fossil Palaeontinidae are described from Daohugou Village, Inner Mongolia, China: Daohugoucossus shii sp. nov., D. parallelivenius sp. nov., D. lii sp. nov. The diagnosis of Dao- hugoucossus Wang, Zhang and Fang, 2006 is revised. These materials are the first complete fossil palaeontinids in the Middle Jurassic of the world. Based on observation of these new specimens, Sc of this new genus has the following combined characters: joined with R at base, diverged near M separa- tion, coalesced with R before the nodal line, extended beyond the nodal line and fused with R1 before distal end. Due to the Sc characteristics mentioned above, the genus Liaocossus Ren, Yin and Dou, 1998 can be well distinguished from “Ilerdocossus-complex”.     

A primitive confuciusornithid bird from China and its implications for early avian flight

Confuciusornithids, lived from 120―125 million years ago, form a basal bird group and include the oldest birds with horny beaks. Here we describe Eoconfuciusornis zhengi, gen. et sp. nov. from the Early Cretaceous Dabeigou Formation (131 Ma) in Fengning, Hebei Province, northern China. It repre- sents a new and, more primitive than other known, member of this group and extends the lifespan of this family to 11 Ma, the longest of any known Early Cretaceous avian lineages. Furthermore, Eoconfuciusornis and its relatives present many osteological transformations, such as the size increase of the deltopectoral crest of the humerus and the keel of the sternum, apparently an adaptation toward improved flight in the evolution of the Confuciusornithidae.     

Temporal variations of <Emphasis Type="Italic">δ</Emphasis><Superscript>18</Superscript>O of atmospheric water vapor at Delingha

Oxygen stable isotope of atmospheric water vapor is widely used to study the modern process of climate. Atmospheric water vapor samples were collected at Dlingha, northeast of Tibetan Plateau during the period from July 2005 to February 2006. The variation of δ ¹⁸O and the relationships between δ ¹⁸O and both the temperature and specific humidity are analyzed in this paper. Results show that the seasonal variation of δ ¹⁸O of atmospheric water vapor at Delingha is remarkable with higher δ ¹⁸O in summer and lower δ ¹⁸O in winter. The temporal variation of vapor δ ¹⁸O shows obvious fluctuations, with magnitude of over 37‰ The daily variation of the δ ¹⁸O is highly correlated with air temperature. The relationship between δ ¹⁸O and atmospheric water vapor content is complex. Study shows that δ ¹⁸O of atmospheric water vapor is positively correlated with specific humidity in winter in seasonal scale and inversely correlated with specific humidity in summer rainy period. The δ ¹⁸O values of atmospheric water vapor are lower than those of precipitation at Delingha, and the average difference is 10.7‰ Variations of δ ¹⁸O of atmospheric water vapor is also found to be affected by precipitation events, The model results show that the precipitation effect could have caused the vapor δ ¹⁸O in the raining season to lower by 7% in average in July and August. Supported by the National Natural Science Foundation of China (Grant Nos. 40671043, 40571039 and 40771048) and National Basic Research Program of China (Grant No. 2005CB422002)     

Vp and Vp/Vs structures in the crust and upper mantle of the Taiwan region, China

A tomographic study of the V _p and V _p/V _s structures in the crust and upper mantle beneath the Taiwan region of China is conducted by simultaneous inversion of P and S arrival times. Compared with the previous tomographic results, the spherical finite difference technique is suitable for the strong heterogeneous velocity structure, and may improve the accuracy in the travel time and three-dimensional ray tracing calculations. The V _p and V _p/V _s structures derived from joint inversion and the relocated earthquakes can provide better constraints for analyzing the lateral heterogeneity and deep tectonic characters in the crust and upper mantle. Our tomographic results reveal significant relations between the seismic wavespeed structure and the tectonic characters. In the shallow depth, sedimentary basins and orogen show distinct wavespeed anomalies, with low V _p, high V _p/V _s in basins and high V _p, low V _p/V _s in orogen. As the suture zone of Eurasian Plate and Philippine Sea Plate, Longitudinal Valley is characterized by a significant high V _p/V _s anomaly extending to the middle-lower crust and upper mantle, which reflects the impact of rock cracking, partial melting, and the presence of fluids. In the northeast Taiwan, the V _p, V _p/V _s anomalies and relocated earthquakes depict the subducting Philippine Sea Plate under the Eurasian Plate. The high V _p of oceanic plate and the low V _p, high V _p/V _s atop the subducted oceanic plate extend to 80 km depth. Along the east-west profiles, the thickness of crust reaches 60 km at the east of Central Range with eastward dipping trend, which reveals the eastward subduction of the thickened and deformed crust of the Eurasian continental plate. Supported by Knowledge Innovation Program of the Chinese Academy of Sciences (Grant No. KZCX3-SW-234-2), National Basic Research Program of China (Grant No. 2007CB411701), National High Technology Research and Development Program of China (Grant No. 2006AA09A101-0201) and National Natural Science Foundation of China (Grant Nos. 40804016, 40704013)     

Geochemical characteristics and genesis of Dachadaban ophiolite in North Qilian area

There are three groups of pillow volcanics in the Dachadaban ophiolite. Group 1 is typical boninite, enriched in Si, Mg and depleted in Ti, HREE and HFSE; group 2 is the evolved boninite, slightly higher abundance of Ti, HREE, HFSE and large variation of Mg’ due to fractional crystallization; and group 3 is tholeiite with MORB character. The existence of MORB and boninite indicates that the Dachadaban ophiolite was probably formed in islandare and back-arc environments. Project supported by both the Department of Foundation Research and High Technology, State Science and Technology Commission of PRC, and the National Natural Science Foundation of China (Grant No. 49472101).     

Trace elements in tree rings and their environmental effects: A case study in Xi’an City

Using wet digested method and ICP mass spectrometer, we analyzed the concentration of five trace elements (Cd, Mn, P, Zn and Pb) for the tree rings from both urban and suburbs of Xi’an. At the urban sampling site, one Chinese mahogany (Toona sinensis) disc and one phoenix tree (Firmiana simplex) disc were sampled from a steelworks in Xi’an City. At the suburb site, a Chinese mahogany disc was collected from a village in the south of the City. In addition, some soils near the roots of the sampled trees were ...     

Microbial and molecular fossils from the Permian Zoophycos in South China

The trace fossil Zoophycos was enriched in the lower part of the Guadalupian Maokou Formation of Permian in Laibin, Guangxi of South China, and characterized by observable lamellae within its spreiten. Associated with these lamellae are the prolific coccus- and spirillum-shaped microbe fossils. These microbes morphologically identified are 0.4―3.0 μm in diameter or length. Molecular fossils, including normal alkanes (dominated by C18, without an odd-over-even predominance), acyclic isoprenoids (such as pristane (Pr) and phytane (Ph)), extended tricyclic terpanes, pentacyclic triterpanes, steranes, al-kylcyclohexanes, dibenzothiophenes, benzonaphthothiophenes, benzobisbenzothiophenes, were in-strumentally identified in the wackestone characterized by the occurrence of abundant Zoophycos (composite ichnofabric indices are 4) and crowded microbial fossils. The value of the Pr/Ph ratio is less than 1, indicative of a dysoxic condition. The identification of abundant sulfur compounds (the thio-phene series) related to the contribution of reductive sulfur favors the occurrence of sulphate-reducing bacteria. The association of abundant microbial colonies with microbial molecular fossils within the spreiten suggests that trace fossils Zoophycos would be a multifunctional garden carefully constructed by the Zoophycos-producer, where different microbial colonies were orderly and carefully planted and cultured in different minor lamellae within spreiten. Hence, it is proposed that the Zoophycos-producer symbiosed with microbial colonies on the mutual basis of food supply and the redox conditions.     

Response of needle dark respiration of Pinus koraiensis and Pinus sylvestriformis to elevated CO2 concentrations for four growing seasons’ exposure

The long-term effect of elevated CO2 concentrations on needle dark respiration of two coniferous species-Pinus koraiensis and Pinus sylvestriformis on the Changbai Mountain was investigated using open-top chambers. P. Koraiensis and P. Sylvestriformis were exposed to 700,500μmol·mol-1 CO2 and ambient CO2(approx.350 μmol·mol-1)for four growing seasons. Needle dark respiration was measurd during the second, third and fourth growing seasons' exposure to elevated CO2.The results showed that needle dark respiration rate increased for P. Koraiensis and P. Sylvestriformis grown at elevated CO2 concentrations during the second growing season, could be attributed to the change of carbohydrate and/or nitrogen content of needles. Needle dark respiration of P. Koraiensis was stimulated and that of P. Sylvestriformis was inhibited by elevated CO2 concentrations during the third growing season. Different response of the two tree species to elevated CO2 mainly resulted from the difference in the growth rate. Elevated CO2 concentrations inhibited needle dark respiration of both P. Koraiensis and P. Sylvestriformis during the fourth growing season. There was consistent trend between the short-term effect and the long-term effect of elevated CO2 on needle dark respiration in P. Sylvestriformis during the third growing season by changing measurement CO2 concentrations. However, the short-term effect was different from the long-term effect for P. Koraiensis. Response of dark respiration of P. Koraiensis and P. Sylvestriformis to elevated CO2 concentrations was related to the treatment time of CO2 and the stage of growth and development of plant. The change of dark respiration for the two tree species was determined by the direct effect of CO2 and long-term acclimation. The prediction of the long-term response of needle dark respiration to elevated CO2 concentration based on the short-term response is in dispute.     

Modeling gross primary production of a temperate grassland ecosystem in Inner Mongolia,China, using MODIS imagery and climate data

Carbon fluxes in temperate grassland ecosystems are characterized by large inter-annual variations due to fluctuations in precipitation and land water availability. Since an eddy flux tower has been in operation in the Xilin Gol grassland, which belongs to typical temperate grassland in North China, in this study, observed eddy covariance flux data were used to critically evaluate the biophysical per- formance of different remote sensing vegetation indices in relation to carbon fluxes. Furthermore, vegetation photosynthesis model (VPM) was introduced to estimate gross primary production (GPP) of the grassland ecosystem for assessing its dependability. As defined by the input variables of VPM, Moderate Resolution Imaging Spectroradimeter (MODIS) and standard data product MOD09A1 were downloaded for calculating enhanced vegetation index (EVI) and land surface water index (LSWI). Measured air temperature (Ta) and photosynthetically active radiation (PAR) data were also included for model simulating. Field CO2 flux data, during the period from May, 2003 to September, 2005, were used to estimate the "observed" GPP (GPPobs) for validation. The seasonal dynamics of GPP predicted from VPM (GPPVPM) was compared quite well (R2=0.903, N=111, p<0.0001) with the observed GPP. The ag- gregate GPPVPM for the study period was 641.5 g C·m?2, representing a ~6% over-estimation, compared with GPPobs. Additionally, GPP predicted from other two typical production efficiency model (PEM) represents either higher overestimation or lower underestimation to GPPobs. Results of this study demonstrate that VPM has potential for estimating site-level or regional grassland GPP, and might be an effective tool for scaling-up carbon fluxes.     

Development and deployment of an instrument for measurement of atmospheric peroxy radical by chemical amplification

An instrument for measurement of peroxy radical by chemical amplification (PERCA) has been devel- oped at Peking University (PKU) and deployed in several field campaigns. PKU PERCA measures per- oxy radical via amplification of NO2 by peroxy radical in the presence of NO and CO through a chain reaction. The amount of amplified NO2 is detected by a NO2-luminal chemiluminescence detector. The chain length (CL) of 75±20 (1σ ) was determined routinely during field campaigns using a HO2 source from photolysis of...     

The Pengguan tectonic dome of Longmen Mountains, Sichuan Province: Mesozoic denudation of a Neoproterozoic magmatic arc-basin system

Neoproterozoic igneous and metamorphic complexes occur as tectonic domes in the Longmen Mountains of the western margin of the Yangtze Block, and are important in reconstructing the Rodinian supercontinent and constraining the timing and mechanism of tectonic denudational processes. The Pengguan dome consists of granitic intrusions and metamorphic rocks of the Huangshuihe Group and is tectonically overlain by ductilly deformed Sinian to Paleozoic strata. The plutonic intrusions consist of granites with abundant amphibolite enclaves. New LA-ICP-MS zircon U-Pb dating yielded an emplacement age of 809±3 Ma and a protolith age of 844±6 Ma for the granite. The granitic rocks have geochemical signatures typical of A-type granites, indicating their formation under an extensional environment, by melting of newly formed tonalite-trondhjemite-granodiorite (TTG) rocks. A detachment fault, characterized by variable ductile shear deformation of S-C fabric and ESE-ward kinematics, separates the Pengguan dome from the Sinian-Paleozoic cover. 40Ar/39Ar dating of muscovite from the mylonite in the detachment fault of the dome demonstrates that ductile deformation occurred at ~160 Ma. This study indicates the existence of a Neoproterozoic magmatic arc-basin system, which was denudated by a Jurassic middle crustal ductile channel flow along the Longmenshan thrust belt.     

Evolving patterns of the fluids within the TAG hydrothermal field

The mixing of seawater/hydrothermal fluid within the large seafloor hydrothermal sulfide deposits plays a key role in the formation processes of the sulfide deposits. Some issues attract considerable attentions in the study of seafloor hydrothermal system in recent years, such as the relationships among different types of vent fluids, the characteristics of chemical compositions and mineral assemblages of the hydrothermal deposits and their governing factors. Combined with the measured data of hydrothermal fluid in the TAG field, the thermodynamic model of mixing processes of the heated seawater at different temperatures and the hydrothermal fluid is calculated to understand the precipitation mechanism of anhydrite and the genetic relationships between the black and white smoker fluids within the TAG mound. The results indicate that the heating of seawater and the mixing of hydrothermal fluid/seawater are largely responsible for anhydrite precipitation and the temperature of the heated seawater is not higher than 150°C and the temperature of the end-member hydrothermal fluid is not lower than 400°C. Based on the simulated results, the evolving patterns of fluids within the TAG deposit are discussed. The mixed fluid of the end-member hydrothermal fluid and the seawater heated by wall rock undergoes conductive cooling during upflowing within the deposit and forms “White Smoker” eventually. In addition, the end-member hydrothermal fluid without mixed with seawater, but undergoing conductive cooling, vents out of the deposit and forms “Black Smoker”. Supported by China Ocean Mineral Resources Research and Development Association Program (Grant No. DY115-02-1-01) and National Basic Research Program of China (Grant No. G2000078503)     

Relations between oxygen stable isotopic ratios in precipitation and relevant meteorological factors in southwest China

The correlations of isotopic ratios in precipitation with temperature, air pressure and humidity at dif- ferent altitudes, in southwest China, are analyzed. There appear marked negative correlations for the δ 18O in precipitation with precipitation amount, vapor pressure and atmospheric precipitable water (PW) at Mengzi, Simao and Tengchong stations on synoptic timescale; the marked negative correlations between the δ 18O in precipitation and the diurnal mean temperature at 400 hPa, 500 hPa, 700 hPa and 850 hPa are different from the temperature effect in middle-high-latitude inland. Moreover, the notable positive correlation between the δ 18O in precipitation and the dew-point deficit △Td at different altitudes is found at the three stations. On annual timescale, the annual precipitation amount weighted mean δ 18O display the negative correlations not only with annual precipitation but also with annual mean temperature at 500 hPa. It can be deduced that, in the years with abnormally strong summer monsoon, more warm and wet air from low-latitude oceans is transported northward along the vapor channel located in southwest China and generates abnormally strong rainfall on the way. Meanwhile, the ab- normally strong condensation process will release more condensed latent heat into atmosphere, and lead to the rise of atmospheric temperature during rainfall, but decline of the δ 18O in precipitation. On the contrary, in the years with abnormally weak summer monsoon, the abnormally weak condensation process will release less condensed latent heat into atmosphere, and lead to the decline of atmos- pheric temperature during rainfall, but increase of the δ 18O in precipitation.     

Uncertainty analysis of CO<Subscript>2</Subscript> flux components in subtropical evergreen coniferous plantation

We present an uncertainty analysis of ecological process parameters and CO₂ flux components (R _eco, NEE and gross ecosystem exchange (GEE)) derived from 3 years’ continuous eddy covariance measurements of CO₂ fluxes at subtropical evergreen coniferous plantation, Qianyanzhou of ChinaFlux. Daily-differencing approach was used to analyze the random error of CO₂ fluxes measurements and bootstrapping method was used to quantify the uncertainties of three CO₂ flux components. In addition, we evaluated different models and optimization methods in influencing estimation of key parameters and CO₂ flux components. The results show that: (1) Random flux error more closely follows a double-exponential (Laplace), rather than a normal (Gaussian) distribution. (2) Different optimization methods result in different estimates of model parameters. Uncertainties of parameters estimated by the maximum likelihood estimation (MLE) are lower than those derived from ordinary least square method (OLS). (3) The differences between simulated R_eco, NEE and GEE derived from MLE and those derived from OLS are 12.18% (176 g C·m⁻²·a⁻¹), 34.33% (79 g C·m⁻²·a⁻¹) and 5.4% (92 g C·m⁻²·a⁻¹). However, for a given parameter optimization method, a temperature-dependent model (T_model) and the models derived from a temperature and water-dependent model (TW_model) are 1.31% (17.8 g C·m⁻²·a⁻¹), 2.1% (5.7 g C·m⁻²·a⁻¹), and 0.26% (4.3 g C·m⁻²·a⁻¹), respectively, which suggested that the optimization methods are more important than the ecological models in influencing uncertainty in estimated carbon fluxes. (4) The relative uncertainty of CO₂ flux derived from OLS is higher than that from MLE, and the uncertainty is related to timescale, that is, the larger the timescale, the smaller the uncertainty. The relative uncertainties of R_eco, NEE and GEE are 4%−8%, 7%−22% and 2%−4% respectively at annual timescale. Supported by the National Natural Science Foundation of China (Grant No. 30570347), Innovative Research International Partnership Project of the Chinese Academy of Sciences (Grant No. CXTD-Z2005-1) and National Basic Research Program of China (Grant No. 2002CB412502)