Magnetic signature of different vegetation species in polluted environment

Detailed magnetic study on vegetation samples from several strongly polluted and clean sites in Bulgaria is carried out in order to evaluate suitability of different species as passive dust collectors in magnetometry. From each location, available species among lichens, mosses, poplar leaves, dandelion, needles have been sampled. Magnetic susceptibility calculated on mass-specific basis shows wide variability between diamagnetic signal up to 846 × 10⁻⁸ m³/kg. Lichens and mosses are found to be the species, showing magnetic signals with the strongest contrast between clean and polluted environment. The main magnetic phase is magnetite-like according to the results from thermomagnetic analysis of susceptibility on magnetic extracts. Scanning electron microscopy (SEM) microphotographs reveal the presence of abundant particulate matter on vegetation surface both with anthropogenic (spherules) and lithogenic origin. Magnetic grain size deduced by the ratio of saturation remanent magnetization (SIRM) and mass-specific magnetic susceptibility (χ) and coercivities (B_c and B_cr) suggest that different species accumulate preferentially small SD-like grains from pollution emissions. Contrasting relationship of the ratio of anhysteretic remanent magnetization (ARM) and χ for polluted vs clean sites deduced by needles and lichens may be related to transformation of the accumulated dust particles within lichens’ tissue. This finding indicates that the exact species used as biological dust collector is of importance when studying spatial grain size distribution of magnetic dust particles. Pilot study on polycyclic aromatic hydrocarbons (PAH) content and its relation to magnetic parameters shows good correspondence between high levels of PAHs and high SIRM values for locations affected by non-ferrous industrial production.     

Carriers of magnetism and vertical inhomogeneities in ultramafic-mafic rocks of the Kondyor Massif (Aldan shield)

Dunite samples from a borehole drilled in the platiniferous concentrically-zoned Kondyor Massif are studied by electron spin resonance (ESR). The spectrum profiles, relative intensities I, and volume magnetic susceptibilities κ are analyzed. These values experience strong irregular variations, sometimes by an order of magnitude, in the upper and medium parts of the column, at depths from 100 to 400 m; and the variations decay at greater depths. The magnetic properties of the samples are determined by iron (II) ions in the olivine lattice and by iron (III) ions in the magnetite and pyrrhotite microphases and in the products of breakdown of the solid solution: chromiferous magnetite, chromoferrite, etc. The I and κ values are directly related: κ_max = 27.8 × 10⁻³ SI units, κ_min=2.63 × 10⁻³ SI units, and κ_mean = 12.7 × 10⁻³ SI units. The maximum κ values are found in the zones with elevated contents of magnetite and pyrrhotite particles, and the minimum ones, in zones with few medium and small clusters with Fe³⁺ ions. The uneven distributions of solid solutions and magnetic phases over depths are suggested to be related to the disturbances in the conditions of crystallization.     

Low-field variation of magnetic susceptibility measured by the KLY-4S Kappabridge and KLF-4A magnetic susceptibility meter: Accuracy and interpretational programme

The KLY-4S Kappabridge and KLF-4A Magnetic Susceptibility Meter enable automated measurement of susceptibility variation with field in the ranges of 2–450 A/m and 5–300 A/m (in effective values), respectively. Unfortunately, the measurement accuracy decreases with decreasing field and it is not easy to decide whether the susceptibility variation at the lowest fields is natural phenomenon or results from measuring errors. To overcome this problem, the accuracies of both the above instruments were investigated experimentally using artificial specimens (mixture of pure magnetite and plaster of Paris) with variable susceptibilities ranging from 1 × 10⁻⁵ to 5 × 10⁻². The complete curve of the field variation of susceptibility of each specimen was measured 10 times and the relative error was calculated for each field. In the KLY-4S Kappabridge, in specimens with susceptibilities higher than 100 × 10⁻⁶, the relative errors are lower than 3% in all fields and lower than 1% in the fields stronger than 10 A/m. In the KLF-4A Magnetic Susceptibility Meter, in relatively strongly magnetic specimens with susceptibilities 5 × 10⁻⁴ to 5 × 10⁻², the relative error is less than 1.5% in the entire field range. While the former instrument is convenient for investigating almost all rock types, the latter instrument is convenient for measuring moderately and strongly magnetic specimens. To facilitate work with field variation of susceptibility curves, showing variable accuracies with field, the programme FieldVar was written. One of its options is plotting the measured data with corresponding field-variable error bars. In this way, a tool is offered for interpreting such susceptibility changes that are sound and reasonable from the point of view of measuring accuracy.     

Magnetic characteristics of dustfall in urban area of north China and its environmental significance

Urban particulate pollutants not only affect urban air quality but also directly threaten public health. Lanzhou has been one of the most seriously polluted cities in the world due to its special geophysical location and weather conditions, which makes it an ideal place for urban particulate pollution study. In this study a set of environmental magnetic parameters (χlf, χfd%, χARM, SIRM, HIRM, SOFT and back field IRM) were measured and analyzed for the dustfall samples collected monthly during 1997.4―2000.5 and 2004.11―2005.10 at Lanzhou University campus (KLD and CCS) and at a clear site of Gaolanxian (GLX). Results demonstrate that the main magnetic minerals in the dustfall are magnetite, with some maghaemite and haematite. Cross plots of Mrs/Ms vs. Bcr/Bc and χfd% vs. χARM/SIRM indicate that the main magnetic grain sizes are pseudo-single-domain (PSD). The high concentration of magnetic min-erals in the dustfall indicates that Lanzhou is seriously polluted. The main pollution sources are an-thropogenic activities and nature dust, with the former posing a great threat during the whole year, which is different from the general view that the natural dust is the main contributor in summer. With great efforts of the government and local authority in the past decades, the air quality in Lanzhou has been improved in winter, but it is not so obvious in summer. This study also reveals that as a simple, rapid and accurate technology, environmental magnetism could be a valuable tool for urban pollution study.     

A high-resolution magnetic susceptibility study of a loess/palaeosol couplet at Baoji,China

Summary Magnetic susceptibility, and its frequency dependence, are reported for 288 individual samples spanning the 8.3 m — thick S3 palaeosol/L4 loess couplet at an important site in the Chinese Loess Plateau. The resulting profile demonstrates that there is a very close link between magnetic properties and soil development: soil and loess sub-divisions recognised visually in the field are clearly reflected in both the bulk susceptibility data and in its frequency dependence. As found at all Chinese sites, the distribution of susceptibilities is bi-modal, one peak representing loess (median=0.74×10⁻⁶m³/kg), the other representing palaeosol (median=2.99×10⁻⁶m³/kg). This is the basis of the climatic proxy information. The Baoji section is the site of one of the most detailed grain-size analyses available anywhere (Ding et al., 1994), and we find a strong correlation between the resulting profile and the susceptibility data reported here. Frequency-dependence of susceptibility exhibits slight differences between the loess and palaeosol populations which we interpret as reflecting reduced transport distances during glacial intervals when the Siberian High causes stronger winds and expands southwards.     

Influence of the modified global ocean tide model with local tides of East and South China Seas on load gravity in China and its neighbor area

Introduction With the development of science and technology, the accuracy of gravity measurement is im-proved. The gravity observation with FG5 gravimeter has reached accuracy of μ magnitude. And the high accurate superconducting gravimeter can detect the tiny signal of 10?2 μ magnitude in frequency domain. With the high-accuracy gravity observation on Earth′s surface, the Earth′s tidal parameters can be determined precisely. And the observations can also be used to invert the struc-ture…     

Sorption of Eu<Superscript>3+</Superscript>onto nano-size silica-water interfaces

The sorption of Eu species onto nano-size silica-water interfaces is investigated at pH range of 1―8.5 and the initial Eu concentrations (C_Eu) of 2×10⁻⁵, 2×10⁻⁴ and 2×10⁻³ M using fluorescence spectroscopy. The sorption rate of Eu is initially low, but significantly increases at pH > 4. For the initial C_Eu of 2×10⁻⁵, 2×10⁻⁴ and 2×10⁻³ M, the dissolved Eu species are completely sorbed onto silica-water interfaces at pH = 4.75, −5.8 and 6.6, respectively, with the respective sorption densities of −1.58×10⁻⁸, 1.58×10⁻⁷ and 1.58×10⁻⁶ mol/m². The sorbed Eu species at pH < 6 is aquo Eu³⁺, which is sorbed onto silica-water interfaces as an outer-sphere complex at pH < 5, but may be sorbed as an inner-sphere bidentate complex at 5 < pH < 6, due to the decrease of the N_H2O to −6 at pH = 6. At pH = 6 – 8, Eu(OH)²⁺, Eu(CO₃)⁺and Eu(CO₃)₂ ⁻ form in the solutions, and Eu(CO₃)⁺is dominant at pH = −7.5. These ions may be sorbed onto silica-water interfaces as inner-sphere bidentate complexes or multi-nuclear pre-cipitates.     

Geochemistry and Sr-Nd isotope characteristics of tonalites in Zêtang, Tibet: New evidence for intra-Tethyan subduction

Tonalites from the island arc rock assemblage in the Zêtang segment of the Yarlung Zangbo suture zone were analyzed for major, trace elements (including REE) and Sr-Nd isotope. The experimental datademonstrate that the tonalites have the adakite-like characteristics, including high SiO2 (58%-63%),Al2O3 (18.4%-22.4%), Sr (810×10-6-940×10-6), Sr/Y (77-106), low HREE (Y=9×10-6-11×10-6, Yb=1×10-6-1.3×10-6), with LREE enrichment and faint Eu positive anomaly. Isr (0.70421-0.70487) is relatively low whereas 143Nd/144Nd (0.512896-0.512929) and εNd(t) values ( 6.7- 7.3) are high. These feainvolvement of a small amount of oceanic sediments. The identification of Z(e)tang adakites, derived from slab melting, presents new evidence for the intra-Tethyan subduction and the previous suggestion about the existence of intra-oceanic island arc within Tethys.     

Soil organic carbon storage and soil CO<Subscript>2</Subscript> flux in the alpine meadow ecosystem

High-resolution sampling, measurements of organic carbon contents and ¹⁴C signatures of selected four soil profiles in the Haibei Station situated on the northeast Tibetan Plateau, and application of ¹⁴C tracing technology were conducted in an attempt to investigate the turnover times of soil organic carbon and the soil-CO₂ flux in the alpine meadow ecosystem. The results show that the organic carbon stored in the soils varies from 22.12×10⁴ kg C hm⁻² to 30.75×10⁴ kg C hm⁻² in the alpine meadow ecosystems, with an average of 26.86×10⁴ kg C hm⁻². Turnover times of organic carbon pools increase with depth from 45 a to 73 a in the surface soil horizon to hundreds of years or millennia or even longer at the deep soil horizons in the alpine meadow ecosystems. The soil-CO₂ flux ranges from 103.24 g C m⁻² a⁻¹ to 254.93 gC m⁻² a⁻¹, with an average of 191.23 g C m⁻² a⁻¹. The CO₂ efflux produced from microbial decomposition of organic matter varies from 73.3 g C m⁻² a⁻¹ to 181 g C m⁻² a⁻¹. More than 30% of total soil organic carbon resides in the active carbon pool and 72.8%281.23% of total CO₂ emitted from organic matter decomposition results from the topsoil horizon (from 0 cm to 10 cm) for the Kobresia meadow. Responding to global warming, the storage, volume of flow and fate of the soil organic carbon in the alpine meadow ecosystem of the Tibetan Plateau will be changed, which needs further research. Supported by the National Natural Science Foundation of China (Grant Nos. 40231015, 40471120 and 40473002) and the Guangdong Provincial Natural Science Foundation of China (Grant No. 06300102)     

Long-term variations in the second sectorial Stokes harmonics on the basis of TOPEX/POSEIDON altimetry between 1993 and 2000

The long-term variations in the second degree sectorial Stokes parameters of the geopotential have been determined from TOPEX-POSEIDON (T/P) satellite altimeter data, covering the period of January 1, 1993 to January 3, 2001 (T/P cycles 11-305). It is the first attempt to determine the variations in the second sectorial harmonics in the Earth’s inertia tensor due to the ocean dynamics. The variations amount to about 1 × 10⁻¹⁰ (J ₂ ⁽²⁾ ≈ 1.6 × 10⁻⁶ and S ₂ ⁽²⁾ ≈ −0.9 × 10⁻⁶). The variations are about 5% of the tidal effect. This corresponds to variations in the directions of the equatorial axes of the Earth’s inertia ellipsoid of up to 10 arc-seconds. Consequently, the annual and semi-annual variations of the Earth’s equatorial flattening is about 10⁻⁹; i.e. it corresponds to a change of 8 units of its denominator of 91 030. (The equatorial flattening ≈ 1/91 030). Since the coverage of the Earth’s ocean surface is not worldwide, and the inclination of T/P is i = 66°, it is only 58.2% (min. depth of the ocean 2 000 m) of the Earth’s surface which is processed, the torque, resulting from the seasonal transfer of masses within a sea surface layer, is not zero. It amounts up to 10¹⁶ kg m²s⁻², which is comparable to the total indirect tidal torque due to the Moon and the Sun, ∼ 3.9 × 10¹⁶ kg m²s⁻². However, the above estimate strongly depends on the adopted thickness of the sea surface layer, ΔR = 50 m. For a larger thickness of ΔR = 100 m, the seasonal torque amounts to about ∼ 2.3 × 10¹⁶ kg m²s⁻².     

Comparison between different earthquake magnitudes determined by China Seismograph Network

By linear regression and orthogonal regression methods, comparisons are made between different magnitudes (lo-cal magnitude ML, surface wave magnitudes MS and MS7, long-period body wave magnitude mB and short-period body wave magnitude mb) determined by Institute of Geophysics, China Earthquake Administration, on the basis of observation data collected by China Seismograph Network between 1983 and 2004. Empirical relations between different magnitudes have been obtained. The result shows that: 1 As different magnitude scales reflect radiated energy by seismic waves within different periods, earthquake magnitudes can be described more objectively by using different scales for earthquakes of different magnitudes. When the epicentral distance is less than 1 000 km, local magnitude ML can be a preferable scale; In case M<4.5, there is little difference between the magnitude scales; In case 4.5MS, i.e., MS underestimates magnitudes of such events, therefore, mB can be a better choice; In case M>6.0, MS>mB>mb, both mB and mb underestimate the magnitudes, so MS is a preferable scale for deter-mining magnitudes of such events (6.08.5, a saturation phenomenon appears in MS, which cannot give an accurate reflection of the magnitudes of such large events; 2 In China, when the epicentral distance is less than 1 000 km, there is almost no difference between ML and MS, and thus there is no need to convert be-tween the two magnitudes in practice; 3 Although MS and MS7 are both surface wave magnitudes, MS is in general greater than MS7 by 0.2~0.3 magnitude, because different instruments and calculation formulae are used; 4 mB is almost equal to mb for earthquakes around mB4.0, but mB is larger than mb for those of mB≥4.5, because the periods of seismic waves used for measuring mB and mb are different though the calculation formulae are the same.     

Simulating bubble number density of rhyolitic pumices from Plinian eruptions: constraints from fast decompression experiments

Decompression experiments of a crystal-free rhyolitic liquid with ≈ 6.6 wt. % H₂O were carried out at a pressure range from 250 MPa to 30–75 MPa in order to characterize effects of magma ascent rate and temperature on bubble nucleation kinetics, especially on the bubble number density (BND, the number of bubbles produced per unit volume of liquid). A first series of experiments at 800°C and fast decompression rates (10–90 MPa/s) produced huge BNDs (≈ 2 × 10¹⁴ m⁻³ at 10 MPa/s ; ≈ 2 × 10¹⁵ m⁻³ at 90 MPa/s), comparable to those in natural silicic pumices from Plinian eruptions (10¹⁵–10¹⁶ m⁻³). A second series of experiments at 700°C and 1 MPa/s produced BNDs (≈ 9×10¹² m⁻³) close to those observed at 800°C and 1 MPa/s (≈ 6 × 10¹² m⁻³), showing that temperature has an insignificant effect on BNDs at a given decompression rate. Our study strengthens the theory that the BNDs are good markers of the decompression rate of magmas in volcanic conduits, irrespective of temperature. Huge number densities of small bubbles in natural silicic pumices from Plinian eruptions imply that a major nucleation event occurs just below the fragmentation level, at which the decompression rate of ascending magmas is a maximum (≥ 1 MPa/s).     

Rock magnetic properties of a loess-paleosol couple along an N-S transect in the Chinese Loess Plateau

We have conducted detailed rock magnetic experiments on samples from loess unit 8 (L8) and paleosol unit 8 (S8) in Jingbian, Yichuan and Duanjiapo loess sections along an N-S transect in the Chinese Loess Plateau. Major rock magnetic results are as follows: (i) An increase of high field susceptibility (χ_h) in the same level of loess or paleosol from north to south is observed, suggesting an enhancement of pedogenesis. (ii) The low field susceptibility(χL) in loess unit L8 is almost the same in three sections. In contrast, the χ_L of paleosol unit S8 in Yichuan is highest, and the χ_L of Duanjiapo is lower than that in Yichuan section, suggesting that there is not correlation between the χ_L and the degree of pedogenesis in loess-paleosols. (iii) With the increasing of χ_L, both the contents of the superparamagnetic (SP) and the ferrimagnetic grains in loess-paleosol increase, however, the enhancement of magnetic susceptibility is probably dependent more on the increase of the ferrimagnetic concentration than on a change in the grain size. (iv) The content of the maghemite in loess unit 8 increases from the northern to the southern part of the Chinese Loess Plateau, and is positively correlated with the pedogenesis of the loess.     

Specific features of development of two-dimensional convection in a rotating compressible conducting fluid

An axisymmetric model of convection in a rotating cylinder in an external uniform magnetic field has been considered. In the considered model, the meridional circulation is created by a nonuniform rotation of the lower boundary relative to the other boundaries. In the considered model, the time of formation of the stationary regime in the magnetic field considerably increases if the vertical density (compressibility) inhomogeneity is taken into account for Ekman numbers of E = E _M = 3 × 10⁻³. This example shows that the compressibility of a medium should be taken into account in the convection and dynamics of the magnetic field when the magnetohydrodynamics of the Earth is analyzed.     

Magnetic,paleomagnetic and palynologic studies of Paleolithic depositions of the Akhshtyrskaya cave (Russia)

The paper presents the results of experimental rock-magnetic, paleomagnetic and palynologic study of Paleolithic sediments sampled along two profiles in the Akhshtyrskaya cave, situated in the vicinity of Black Sea shore. In the upper part of profiles, some magnetite was observed; in the middle and lower parts, strongly oxidized non-stoichiometric magnetite and hematite prevail. Thin maghemite covers on the surface of fine magnetite grains are present in the majority of specimens. Natural remanence has one characteristic component (CHRM), mostly of chemical origin, although in few specimens containing magnetite it may be sedimentary. Directions of CHRM obtained by standard paleomagnetic methods revealed anomalous pattern only in layer 3/2, which is slightly older than the overlying layer 3/1 whose age was established as (35±2)×10³ years BP by the U-Th method. This suggests that this paleomagnetic anomaly (PMA) can be correlated with Kargapolovo excursion dated on about (45−39)×10³ years BP. In the remaining overlying and underlying layers, directions of CHRM are grouped around the present geomagnetic field. Depth distributions of scalar magnetic parameters generally coincide with the lithological division of the profiles. Palynologic study revealed the presence of 22 pollen zones. Five thermomers separated with colder periods were found in the middle and lower parts of profile. The non-magnetite composition of magnetic fraction of the majority of studied sediments — oxidized nonstoichiometric magnetite and hematite — resulted in the lack of correlations between paleoclimatic and scalar magnetic characteristics.     

CO<Subscript>2</Subscript> and He degassing at El Chichón volcano,Chiapas, Mexico: gas flux,origin and relationship with local and regional tectonics

During 2007–2008, three CO₂ flux surveys were performed on El Chichón volcanic lake, Chiapas, Mexico, with an additional survey in April 2008 covering the entire crater floor (including the lake). The mean CO₂ flux calculated by sequential Gaussian simulation from the lake was 1,190 (March 2007), 730 (December 2007) and 1,134 g m⁻² day⁻¹ (April 2008) with total emission rates of 164 ± 9.5 (March 2007), 59 ± 2.5 (December 2007) and 109 ± 6.6 t day⁻¹ (April 2008). The mean CO₂ flux estimated from the entire crater floor area was 1,102 g m⁻² day⁻¹ for April 2008 with a total emission rate of 144 ± 5.9 t day⁻¹. Significant change in CO₂ flux was not detected during the period of survey, and the mapping of the CO₂ flux highlighted lineaments reflecting the main local and regional tectonic patterns. The ³He/⁴He ratio (as high as 8.1 R _A) for gases in the El Chichón crater is generally higher than those observed at the neighbouring Transmexican Volcanic Belt and the Central American Volcanic Arc. The CO₂/³He ratios for the high ³He/⁴He gases tend to have the MORB-like values (1.41 × 10⁹), and the CO₂/³He ratios for the lower ³He/⁴He gases fall within the range for the arc-type gases. The high ³He/⁴He ratios, the MORB-like CO₂/³He ratios for the high ³He/⁴He gases and high proportion of MORB-CO₂ (M = 25 ±15%) at El Chichón indicate a greater depth for the generation of magma when compared to typical arc volcanoes.     

Net autotrophy in a fluvial lake: the relative role of phytoplankton and floating-leaved macrophytes

This study combined water- and sediment flux measurements with mass balances of dissolved gas and inorganic matter to determine the importance of pelagic and benthic processes for whole-system metabolism in a eutrophic fluvial lake. Mass balances of dissolved O₂, inorganic carbon (DIC), nitrogen (DIN), phosphorous (SRP), particulate N (PN) and P (PP) and Chl a were calculated at a nearly monthly frequency by means of repeated sampling at the lake inlet and outlet. Simultaneously, benthic fluxes of gas and nutrients, including denitrification rates, and the biomass of the dominant pleustophyte (Trapa natans) were measured, and fluxes of O₂ and CO₂ across the water–atmosphere interface were estimated from diel changes in outlet concentrations. On an annual scale, Middle Lake exhibited CO₂ supersaturation, averaging 313% (range 86–562%), but was autotrophic with a net O₂ production (6.35 ± 2.05 mol m⁻² y⁻¹), DIC consumption (−31.18 ± 18.77 mol m⁻² y⁻¹) and net export of Chl a downstream (8.38 ± 0.95 mol C m⁻² y⁻¹). Phytoplankton was the main driver of Middle Lake metabolism, with a net primary production estimated at 33.24 mol O₂ m⁻² y⁻¹, corresponding to a sequestration of 4.18 and 0.26 mol m⁻² y⁻¹ of N and P, respectively. At peak biomass, T. natans covered about 18% of Middle Lake’s surface and fixed 2.46, 0.17 and 0.02 mol m⁻² of C, N and P, respectively. Surficial sediments were a sink for O₂ (−14.47 ± 0.65 mol O₂ m⁻² y⁻¹) and a source of DIC and NH₄ ⁺ (18.84 ± 2.80 mol DIC m⁻² y⁻¹ and 0.83 ± 0.16 mol NH₄ ⁺ m⁻² y⁻¹), and dissipated nitrate via denitrification (1.44 ± 0.11 mol NO₃ ⁻ m⁻² y⁻¹). Overall, nutrient uptake by primary producers and regeneration from sediments were a minor fraction of external loads. This work suggests that the creation of fluvial lakes can produce net autotrophic systems, with elevated rates of phytoplanktonic primary production, largely sustained by allochtonous nutrient inputs. These hypereutrophic aquatic bodies are net C sinks, although they simultaneously release CO₂ to the atmosphere.     

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)