Diurnal, seasonal and annual variation in net ecosystem CO2 exchange of an alpine shrubland on Qinghai-Tibetan plateau

Thus far, grassland ecosystem research has mainly been focused on low-lying grassland areas, whereas research on high-altitude grassland areas, especially on the carbon budget of remote areas like the Qinghai-Tibetan plateau is insufficient. To address this issue, flux of CO2 were measured over an alpine shrubland ecosystem (37 degrees 36'N, 101 degrees 18'E; 325 above sea level [a. s. l.]) on the Qinghai-Tibetan Plateau, China, for 2 years (2003 and 2004) with the eddy covariance method. The vegetation is dominated by formation Potentilla fruticosa L. The soil is Mol-Cryic Cambisols. To interpret the biotic and abiotic factors that modulate CO2 flux over the course of a year we decomposed net ecosystem CO2 exchange (NEE) into its constituent components, and ecosystem respiration (R-eco). Results showed that seasonal trends of annual total biomass and NEE followed closely the change in leaf area index. Integrated NEE were -58.5 and -75.5 g C m(-2), respectively, for the 2003 and 2004 years. Carbon uptake was mainly attributed from June, July, August, and September of the growing season. In July, NEE reached seasonal peaks of similar magnitude (4-5 g C m(-2) day(-1)) each of the 2 years. Also, the integrated night-time NEE reached comparable peak values (1.5-2 g C m(-2) day(-1)) in the 2 years of study. Despite the large difference in time between carbon uptake and release (carbon uptake time < release time), the alpine shrubland was carbon sink. This is probably because the ecosystem respiration at our site was confined significantly by low temperature and small biomass and large day/night temperature difference and usually soil moisture was not limiting factor for carbon uptake. In general, R-eco was an exponential function of soil temperature, but with season-dependent values of Q(10). The temperature-dependent respiration model failed immediately after rain events, when large pulses of R-eco were observed. Thus, for this alpine shrubland in Qinghai-Tibetan plateau, the timing of rain events had more impact than the total amount of precipitation on ecosystem R-eco and NEE.

Depression of net ecosystem CO2 exchange in semi-arid Leymus chinensis steppe and alpine shrub

Uptake and release of carbon in grassland ecosystems is very critical to the global carbon balance and carbon storage. In this study, the dynamics of net ecosystem CO2 exchange (FNEE) of two grassland ecosystems were observed continuously using the eddy covariance technique during the growing season of 2003. One is the alpine shrub on the Tibet Plateau, and the other is the sem-arid Leymus chinensis steppe in Inner Mongolia of China. It was found that the FNEE of both ecosystems was significantly depressed under high solar radiation. Comprehensive analysis indicates that the depression of FNEE in the L. chinensis steppe was the results of decreased plant photosynthesis and increased ecosystem respiration (R-eco) under high temperature. Soil water stress in addition to the high atmospheric demand under the strong radiation was the primary factor limiting the stomatal conductance. In contrast, the depression of FNEE in the alpine shrub was closely related to the effects of temperature on both photosynthesis and ecosystem respiration, coupled with the reduction of plant photosynthesis due to partial stomatal closure under high temperature at mid-day. The R,c of the alpine shrub was sensitive to soil temperature during high turbulence (u* > 0.2 m s(-1)) but its FNEE decreased markedly when the temperature was higher than the optimal value of about 12 degrees C. Such low optimal temperature contrasted the optimal value (about 20 degrees C) for the steppe, and was likely due to the acclimation of most alpine plants to the long-term low temperature on the Tibet Plateau. We inferred that water stress was the primary factor causing depression of the FNEE in the semi-arid steppe ecosystem, while relative high temperature under strong solar radiation was the main reason for the decrease of FNEE in the alpine shrub. This study implies that different grassland ecosystems may respond differently to climate change in the future. (c) 2006 Elsevier B.V All rights reserved.

Temperature and biomass influences on interannual changes in CO2 exchange in an alpine meadow on the Qinghai-Tibetan Plateau

Three years of eddy covariance measurements were used to characterize the seasonal and interannual variability of the CO2 fluxes above an alpine meadow (3250 m a.s.l.) on the Qinghai-Tibetan Plateau, China. This alpine meadow was a weak sink for atmospheric CO2, with a net ecosystem production (NEP) of 78.5, 91.7, and 192.5 g C m(-2) yr(-1) in 2002, 2003, and 2004, respectively. The prominent, high NEP in 2004 resulted from the combination of high gross primary production (GPP) and low ecosystem respiration (R-e) during the growing season. The period of net absorption of CO2 in 2004, 179 days, was 10 days longer than that in 2002 and 5 days longer than that in 2003. Moreover, the date on which the mean air temperature first exceeded 5.0 degrees C was 10 days earlier in 2004 (DOY110) than in 2002 or 2003. This date agrees well with that on which the green aboveground biomass (Green AGB) started to increase. The relationship between light-use efficiency and Green AGB was similar among the three years. In 2002, however, earlier senescence possibly caused low autumn GPP, and thus the annual NEP, to be lower. The low summertime R-e in 2004 was apparently caused by lower soil temperatures and the relatively lower temperature dependence of R-e in comparison with the other years. These results suggest that (1) the Qinghai-Tibetan Plateau plays a potentially significant role in global carbon sequestration, because alpine meadow covers about one-third of this vast plateau, and (2) the annual NEP in the alpine meadow was comprehensively controlled by the temperature environment, including its effect on biomass growth.

A stable carbon isotopic approach for understanding the CO2 flux at the Haibei Alpine Meadow Ecosystem - A simple model

Although respiration of organisms and biomass as well as fossil fuel burning industrial production are identified as the major sources, the CO2 flux is still unclear due to the lack of proper measurements. A mass-balance approach that exploits differences in the carbon isotopic signature (delta(13)C) of CO2 Sources and sinks was introduced and may provide a means of reducing uncertainties in the atmospheric budget. delta(13)C measurements of atmospheric CO2 yielded an average of - 10.3 parts per thousand relative to the Peedee Belemnite standard; soil and plants had a narrow range from -25.09 parts per thousand to -26.51 parts per thousand and averaged at -25.80 parts per thousand. Based on the fact of steady fractionation and enrichment during respiration of mitochondria, we obtained the emission Of CO2 of 35.451 mol m(-2) a(-1) and CO2 flux of 0.2149 mu mol m(-2) s(-)1. The positive CO2 flux indicated the Haibei Alpine Meadow Ecosystem a source rather than a sink. The mass-balance model can be applied for other ecosystem even global carbon cycles because it neglects the complicated process of carbon metabolism, however just focuses on stable carbon isotopic compositions in any of compartments of carbon sources and sinks. (C) 2005 Elsevier B.V. All rights reserved.

Strong temperature dependence and no moss photosynthesis in winter CO2 flux for a Kobresia meadow on the Qinghai-Tibetan plateau

We examined the CO2 exchange of a Kobresia meadow ecosystem on the Qinghai-Tibetan plateau using a chamber system. CO2 efflux from the ecosystem was strongly dependence on soil surface temperature. The COZ efflux-temperature relationship was identical under both light and dark conditions, indicating that no photosynthesis could be detected under light conditions during the measurement period. The temperature sensitivity (Q(10)) of the COZ efflux showed a marked transition around -1.0 degrees C; Q(10) was 2.14 at soil surface temperatures above and equal to -1.0 degrees C but was 15.3 at temperatures below -1.0 degrees C. Our findings suggest that soil surface temperature was the major factor controlling winter COZ flux for the alpine meadow ecosystem and that freeze-thaw cycles at the soil surface layer play an important role in the temperature dependence of winter CO2 flux. (c) 2005 Elsevier Ltd. All rights reserved.

Characterizing CO2 fluxes for growing and non-growing seasons in a shrub ecosystem on the Qinghai-Tibet Plateau

To assess carbon budget for shrub ecosystems on the Qinghai-Tibet Plateau, CO2 flux was measured with an open-path eddy covariance system for an alpine shrub ecosystem during growing and non-growing seasons. CO2 flux dynamics was distinct between the two seasons. During the growing season from May to September, the ecosystem exhibited net CO2 uptake from 08:00 to 19:00 (Beijing Standard Time), but net CO2 emission from 19:00 to 08:00. Maximum CO2 uptake appeared around 12:00 with values of 0.71, 1,19, 1.46 and 0.67 g CO2 m(-2) h(-1) for June, July, August and September, respectively. Diurnal fluctuation Of CO2 flux showed higher correlation with photosynthetic photon flux density than temperature. The maximum net CO2 influx occurred in August with a value of 247 g CO2 m(-2). The total CO2 uptake by the ecosystem was up to 583 g CO2 m(-2) for the growing season. During the non-growing season from January to April and from October to December, CO2 flux showed small fluctuation with the largest net CO2 efflux of 0.30 g CO2 m(-2) h(-1) in April. The diurnal CO2 flux was close to zero during most time of the day, but showed a small net CO2 eff lux from 11:00 to 18:00. Diurnal CO2 flux, is significantly correlated to diurnal temperature in the non-growing season. The maximum monthly net CO2 eff lux appeared in April, with a value of 105 g CO2 m(-2). The total net CO2 eff lux for the whole non-growing season was 356 g CO2 m(-2).

Energy exchange between the atmosphere and a meadow ecosystem on the Qinghai-Tibetan Plateau

To reveal the potential contribution of grassland ecosystems to climate change, we examined the energy exchange over an alpine Kobresia meadow on the northeastern Qinghai-Tibetan Plateau. The annual pattern of energy exchange showed a clear distinction between periods of frozen soil with the daily mean soil temperature at 5 cm (T-s5 &LE; 0 &DEG; C) and non-frozen soil (T-s5 > 0 &DEG; C). More than 80% of net radiation was converted to sensible heat (H) during the frozen soil period, but H varied considerably with the change in vegetation during the non-frozen soil period. Three different sub-periods were further distinguished for the later period: (1) the pre-growth period with Bowen ratio (β) > 1 was characterized by a high β of 3.0 in average and the rapid increase of net radiation associated with the increases of H, latent heat (LE) and soil heat; (2) during the Growth period when β &LE; 1, the LE was high but H fluxes was low with β changing between 0.3 and 0.4; (3) the post-growth period with average β of 3.6 when H increased again and reached a second maximum around early October. The seasonal pattern suggests that the phenology of the vegetation and the soil water content were the major factors affecting the energy partitioning in the alpine meadow ecosystem. © 2005 Elsevier B.V. All rights reserved.

Grazing intensity alters soil respiration in an alpine meadow on the Tibetan plateau

Grazing intensity may alter the soil respiration rate in grassland ecosystems. The objectives of our study were to (1) determine the influence of grazing intensity on temporal variations in soil respiration of an alpine meadow on the northeastern Tibetan Plateau; and (2) characterise, the temperature response of soil respiration under different grazing intensities. Diurnal and seasonal soil respiration rates were measured for two alpine meadow sites with different grazing intensities. The light grazing (LG) meadow site had a grazing intensity of 2.55 sheep ha(-1), while the grazing intensity of the heavy grazing (HG) meadow site, 5.35 sheep ha(-1), was approximately twice that of the LG site. Soil respiration measurements - showed that CO2 efflux was almost twice as great at the LG site as at the HG site during the growing season, but the diurnal and seasonal patterns of soil respiration rate were similar for the two sites. Both exhibited the highest annual soil respiration rate in mid-August and the lowest in January. Soil respiration rate was highly dependent on soil temperature. The Q(10) value for annual soil respiration was lower for the HG site (2.75) than for the LG site (3.22). Estimates of net ecosystem CO2 exchange from monthly measurements of biomass and soil respiration revealed that during the period from May 1998 to April 1999, the LG site released 2040 g CO2 m(-2) y(-1) to the atmosphere, which was about one third more than the 1530g CO2 m(-2) y(-1) released at the HG site. The results suggest that (1) grazing intensity alters not only soil respiration rate, but also the temperature dependence of soil CO2 efflux; and (2) soil temperature is the major environmental factor controlling the temporal variation of soil respiration rate in the alpine meadow ecosystem. (C) 2003 Elsevier Ltd. All fights reserved.

Short-term variation of CO2 flux in relation to environmental controls in an alpine meadow on the Qinghai-Tibetan Plateau

[1] The alpine meadow ecosystem on the Qinghai-Tibetan Plateau may play a significant role in the regional carbon cycle. To assess the CO2 flux and its relationship to environmental controls in the ecosystem, eddy covariance of CO2, H2O, and energy fluxes was measured with an open-path system in an alpine meadow on the plateau at an elevation of 3,250 m. Net ecosystem CO2 influx (Fc) averaged 8.8 g m(-2) day(-1) during the period from August 9 to 31, 2001, with a maximum of 15.9 g m(-2) day(-1) and a minimum of 2.3 g m(-2) day(-1). Daytime Fc averaged 16.7 g m(-2) day(-1) and ranged from 10.4 g m(-2) day(-1) to 21.7 g m(-2) day(-1) during the study period. For the same photosynthetic photon flux density (PPFD), gross CO2 uptake (Gc) was significantly higher on cloudy days than on clear days. However, mean daily Gc was higher on clear days than on cloudy days. With high PPFD, Fc decreased as air temperature increased from 10degreesC to 23degreesC. The greater the difference between daytime and nighttime air temperatures, the more the sink was strengthened. Daytime average water use efficiency of the ecosystem (WUEe) was 8.7 mg (CO2)(g H2O)(-1); WUEe values ranged from 5.8 to 15.3 mg (CO2)(g H2O)(-1). WUEe increased with the decrease in vapor pressure deficit. Daily albedo averaged 0.20, ranging from 0.19 to 0.22 during the study period, and was negatively correlated with daily Fc. Our measurements provided some of the first evidence on CO2 exchange for a temperate alpine meadow ecosystem on the Qinghai-Tibetan Plateau, which is necessary for assessing the carbon budget and carbon cycle processes for temperate grassland ecosystems.

中国北方表沙无机碳储量及次生碳酸盐固碳潜力研究

This paper selected the Taklamakan Desert and the Badain Jaran Desert as the research areas, tested the carbonate content of surface-sand samples of dunes using Eijkelkamp carbonate goniophotometer, and analyzed the spatial-distribution characteristics of carbonate and estimated the carbonate-stock and secondary carbonate-stock in 1m depth of surface sand in the Taklamakan Desert and the Badain Jaran Desert. In addition, the paper test XRD, SEM, TDA, stable carbon isotope and radioactive strontium isotope of lacustrine deposits in the Taklamakan Desert and carbonates, such as kunkar, root canal, lacustrine deposits, sinter and calcrete, in the Badain Jaran Desert. Resting on the achievements by our predecessors, it analyzed the mineral-composition differences of the carbonates, calculated the contents of secondary carbonate and, furthermore, evaluated their potential of sequestration of CO2 in the atmosphere. The overall goal of this study was to increase our understanding of soil carbonate in the context of carbon sequestration in the arid region in China. That is, to advance our understanding about whether or not secondary carbonate in desert is a sink for atmospheric CO2. The following viewpoints were obtained: 1 Carbonate contents of surface-sand samples decend from the south to the north of the Taklamakan Desert. The minimum lies in the south and the maxmum in the mid. Carbonate content of surface-sand of megadunes in the Badain Jaran Desert has low value generally in the dune-crest and the base of slope, and large value in the mid. The average of Carbonate contents of all sorts of collected samples in the same area of the Taklamakan Desert has small diffetences. The average is about 9%. 2 Using carbonate contents as key parameters, calculate the carbon-stock of carbonates in 1m depth of surface sand in the Taklamakan Desert and the Badain Jaran Deser.They are 1.13Pg and 0.19 Pg respectively. There are 0.53Pg and 0.088Pg carbon-stock of secondary-carbonates in 1m depth of surface sand in the Taklamakan Desert and the Badain Jaran Desert. 3 Through testing data from XRD （X-ray diffraction）and TAD ( Thermal Analysis Data), the most significant conclusion derived from is that the main mineral ingredient is calcite in different carbonate substances in arid regions, From the SEM(Scanning electron microscopy ) images, can obtains the information about the micro environment of different carbonate forms in which they can grow. 4 Selected gas by termal cracking and traditional phosphoric acid method, their δ13C show that δ13C is a good parameter to indicate the micro environment in which different secondary carbonate forms. From the δ13C of the same type samples, if the redeposit degree is hard, theδ13C is light, the redeposit degree is weak, the δ13C is heave. and the δ13C of the different type samples, δ13C is mainly controlled by the micro environment in which secondary formed. if the procedure is characterized by redeposit and dissolve of marine facies carbonate, δ13C is heavy, it is characterized by CO2 which produced by plant respiration,δ13C is light. 5 From the δ13C of lacustrine deposit in the different grain size, there exsit certain differences in their micro environment and secondary degree among different grain size in the same grade. 6 The secondary carbonate content of lacustrine deposits in Taklimakan Desert is 47.26%. And those of root canal, sinter, calcrete, kunkar, lacustrine deposit and surface sand in Badain Jaran Desert are 91.74%, 78.46%, 76.26%, 87.87%, 85.37%and 46.49%, respectively. Of different grain size samples, the secondary carbonate contents of coarse fraction (20-63μm), sub-coarse fraction (5-20μm) and fine fraction (<5μm) are 80.10%, 47.2%and 50.07%, respectively. 7 There is no obvious relevance betweenδ13C of secondary carbonate and the content of secondary carbonate,theδ13C of secondary carbonate mainly reflects the parameters of secondary process, the content of secondary carbonate reflects difference of secondary degree.. 8 Silicates potentially supply 3.4 pencent calcium source during forming process of lacustrine deposits in Taklimakan Desert. If calcium source is mainly supplied by goundwater, it can be calculated that about 5.18 %, 6.13%, 5.68%, 5.64 % and 6.82% silicates supply calcium source respectively for root canal, kunkar, lacustrine deposit, calcrete and sinter, during the forming process of different kinds of carbonates in Badain Jaran Desert.

华南埃迪卡拉纪末期硅质岩和深海的氧化还原状态

A suit of cherts deposited in deep-ocean basin of South China during terminal Ediacaran and the beginning of Early Cambrian (about 550～540Ma). The origin of these cherts are controversial, and contrary standpoints exist for the redox state of the terminal Ediacaran deep-ocean because of poor study. In this paper, a detailed sedimentology, element and stable isotope geochemistry study were conducted for cherts of Liuchapo Formation in Anhua County, Hunan Province, Laobao Formation in Sanjiang County, Guangxi Province, and Piyuancun Formation in Xiuning County, Anhui Province. Some conclusions were drawn: (1) These cherts were mainly formed through chemical deposition of dissolved silica derived from chemical weathering in continent. Biological processes played a part in the cherts formation, however, the influence of hydrothermal fluids and detritus materials were trivial. (2) The terminal Ediacaran deep-ocean was anoxic and gradually oxidized. (3) Intense bacterial sulfate reduction decreased sulfate concentration in the ocean and the water column became euxinic during Ediacaran-Cambrian transition period. (4) Due to a high CO2 concentration in the terminal Ediacaran atmosphere, chemical weathering in continent dramatically increased and huge amounts of nutrimental material were transferred into ocean. The biota bloomed because of the nutrimental material. Degradation of huge amounts of organic matter maintained dissolved organic carbon reservoir in the ocean, and prolonged the deep-ocean anoxia.

松辽盆地南部大情字井地区CO2埋存体地质格架及储层建模

Increasing attentions have been paid to the subsurface geological storage for CO2 in view of the huge storage capacity of subsurface reservoirs. The basic requirement for subsurface CO2 storage is that the CO2 should be sequestrated as supercritical fluids (physical trapping), which may also interact with ambient reservoir rocks and formation waters, forming new minerals (chemical trapping). In order to the effective, durable and safe storage for CO2, enough storage space and stable sealing caprock with strong sealing capacity are necessitated, in an appropriate geological framework. Up till now, hydrocarbon reservoirs are to the most valid and appropriate CO2 storage container, which is well proven as the favorable compartment with huge storage capacity and sealing condition. The thesis focuses on two principal issues related to the storage and sealing capacity of storage compartment for the Qingshankou and Yaojia formations in the Daqingzijing block, Southern Songliao Basin, which was selected as the pilot well site for CO2-EOR storage. In the operation area, three facies, including deltaic plain, deltaic front and subdeep-deep lake facies associations, are recognized, in which 11 subfacies such as subaqueous distributary channel, river- mouth bar, interdistributary bay, sheet sandbody, crevasse splay and overflooding plain are further identified. These subfacies are the basic genetic units in the reservoir and sealing rocks. These facies further comprise the retrogradational and progradational depositional cycles, which were formed base- level rise and fall, respectively. During the regressive or lake lowstand stage, various sands including some turbidites and fans occurred mostly at the bottom of the hinged slope. During the progradation stage, these sands became smaller in size and episodically stepped backwards upon the slope, with greatly expanded and deeped lake. However, most of Cretaceous strata in the study area, localized in the basin centre under this stage, are mainly composed of grey or grizzly siltstones and grey or dark grey mudstones intercalated with minor fine sandstones and purple mudstones. On the base of borehole and core data, these siltstones are widespread, thin from 10 to 50 m thick, good grain sorting, and have relative mature sedimentary structures with graded bedding and cross- lamination or crossbeds such as ripples, which reflect strong hydrodynamic causes. Due to late diagenesis, pores are not widespread in the reservoirs, especially the first member of Qingshankou formation. There are two types of pores: primary pore and secondary cores. The primary pores include intergranular pores and micropores, and the secondary pores include emposieus and fracture pores. Throat channels related to pores is also small and the radius of throat in the first, second and third member of Qingshankou formation is only 0.757 μm, 0.802 μm and 0.631 μm respectively. In addition, based on analyzing the probability plot according to frequency of occurrence of porosity and permeability, they appear single- peaked distribution, which reflects strong hetero- geneity. All these facts indicate that the conditions of physical property of reservoirs are not better. One reason may be provided to interpret this question is that physical property of reservoirs in the study area is strong controlled by the depositional microfacies. From the statistics, the average porosity and permeability of microfacies such as subaqueous distributary channel, channel mouth bar, turbidites, is more than 9 percent and 1md respectively. On the contrary, the average porosity and permeability of microfacies including sand sheet, flagstone and crevasse splay are less than 9 percent and 0.2md respectively. Basically, different hydrodynamic environment under different microfacies can decide different physical property. According to the reservoir models of the first member of Qingshankou formation in the No. well Hei47 block, the character of sedimentary according to the facies models is accord to regional disposition evolution. Meantime, the parameter models of physical property of reservoir indicate that low porosity and low permeability reservoirs widespread widely in the study area, but the sand reservoirs located in the channels are better than other places and they are the main sand reservoirs. The distribution and sealing ability of fault- fractures and caprock are the key aspects to evaluate the stable conditions of compartments to store CO2 in the study area. Based on the core observation, the fractures widespread in the study area, especially around the wells, and most of them are located in the first and second member of Qingshankou formation, almost very few in the third member of Qingshankou formation and Yaojia formation instead. In addition, analyzing the sealing ability of eleven faults in the three-dimensional area in the study area demonstrates that most of faults have strong sealing ability, especially in the No. well Hei56 and Qing90-27. To some extent, the sealing ability of faults in the No. well Hei49, Qing4-6 and Qing84-29 are worse than others. Besides, the deposition environment of most of formations in the study area belongs to moderately deep and deep lake facies, which undoubtedly take advantage to caprocks composed of mudstones widespread and large scale under this deposition environment. In the study area, these mudstones distribute widely in the third member of Qingshankou formation, Yaojia and Nenjiang formation. The effective thickness of mudstone is nearly ~550m on an average with few or simple faults and fractures. In addition, there are many reservoir beds with widely- developed insulated interbeds consist of mudstones or silty mudstone, which can be the valid barrier to CO2 upper movement or leakage through diffusion, dispersion and convection. Above all, the closed thick mud caprock with underdeveloped fractures and reservoir beds can be taken regard as the favorable caprocks to provide stable conditions to avoid CO2 leakage.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

湘西地区上震旦统上部－下寒武统底部黑色岩系沉积学及地球化学

Widespread black chert-shales occur in the Ediacaran-Cambrian(E-C) boundary successions along the flank of Yangtze Platform, South China, remarkable changes in sedimentology, geochemistry and biology were recorded. Although extensive studies were carried out upon this boundary succession, the origin of black chert-shales still remain controversial. This paper focuses on the E-C black chert-shales in western Hunan, South China, upon which detailed depositional and geochemical changes are documented, accordingly a depositional model for black chert-shales is proposed. Stratigraphic anatomy across the depositional strike demonstrates that the shallow-water Dengying dolostone along the platform margin sharply pass basinward into the Liuchapo chert successions, which indicate syndepositional extensional faulting at depth could have occurred along the platform margin. The deep-water Niutitang phosphorite-rich black shales are either underlain by the Dengying dolostones on the platform margin toward platform interior or directly by the Liuchaopo chert successions farther basinwards. By detailed investigation, silica chimneys are firsly identified approximately in the chert along platform margin; two types of silica chimneys, including mounded and splayed/funnelized chert(generally brecciated) bodies are further sorted out. The mounded chert are exitbited by domed or hummocky surfaces on the top and irregular spongy to digitiform internal fabrics; within the silica mounds, abundant original vesicles/voids and/or channels were mostly plugged by initial chalcedony, quartze crystals with minor dolomite and bladed barite crystals. Splayed/funnelized brecciated chert “intrusion” cross-cut the uppermost dolostones capping to the horizon underneath, and are directly overlain by the Niutitang phosphorite-rich black shales. Their similarities to the silica chimneys reported from the oceanic spreading centres suggest a similar origin responsible for these unique silica bodies which is also supported by the microthermonmetric data and element geochemistry. High P, Ba, Fe contents and positive correlation between Fe and TOC concentrations in the Niutitang black shales indicate a high palaeo-productivity in the Early Cambrian ocean. The low Th/U and the high V/Cr, V/Sc, V/(V+Ni) ratios in the black shales suggest an anoxic water condition during this interval. Furthermore, Positive Eu anomalies and high Ba contents in the sediments also imply a hydrothermal influence on the formation of Niutitang black shales. To better constrain the placement of deep-water successions straddling the E-C boundary and the timing of hydrothermal silica chimneys, sensitive high-resoluton ion microprobe(SHRIMP) U-Pb dating of zircon grains from tuffs within the chert succession of Liuchapo Formation at Ganziping was conducted and yields a weighted-mean 206Pb/238Pb age at 536.6±5.5Ma, younger than E-C boundary age(542.0±0.3Ma). This age combined with carbon isotopic data is then proposed to correspond to the U-Pb age of zircons(538.2±1.5Ma) from the Zhongyicun member of Meishucun Formation at Meishucun in eastern Yunna, thus, the E-C boundary in Gazngziping was placed between the Dengying formations and Liuchapo formatioms. therefore, the silica chimneys took place at the beginning of the Cambrian period. The temporal coincidence of silica chimneys and negative excursions of δ13C and δ34Spy pairs suggest hydrothermal activities were likely responsible for the isotopic changes. Under such a circumstance, vast amounts of greenhouse gases(CO2, CH4, H2S), with highly 13C-depleted carbon and 34S-depleted sulfur would be released into the ocean and atmosphere. A positive shift in δ34Scas and Δ34S values from the late Ediacaran to the Early Cambrian could be a reflection of enhanced bacterial sulfate reduction(BSR), strengthened by the intensified oceanic anoxia stimulated by hydrothermal activities. Based on the analyses of sedimentology and geochemistry, a model- “oceanic anoxia induced by hydrothermal–volcanic activies” was proposed to responsible for the formation of black chert-shales during this E-C transition. Under this case, hydrothermal-volcanic activies could release large large amount of greenhouse into atmosphere and metal micronutrients into the ocean, which may lead to global warming, stratified ocean, thereby a high palaeoproductivity; on the other hand, the massive releasing of reduced hydrothermal fluids with abundant H2S, could have in turn enhanced the ocean anoxia. All of these were favourable the for preservation of organic matter, and subsequent extensive deposition of black silica-shales.

西江流域环境河水地球化学研究

River is a major component of the global surface water and CO2 cycles. The chemistry of river waters reveals the nature of weathering on a basin-wide scale and helps us understand the exogenic cycles of elements in the continent-river-ocean system. In particular, geochemical investigation of large river gives important information on the biogeochemical cycles of the elements, chemical weathering rates, physical erosion rates and CO2 consumption during the weathering of the rocks within the drainage basin. Its importance has led to a number of detailed geochemical studies on some of the world's large and medium-size river systems. Flowing in the south of China, the Xijiang River is the second largest river in the China with respect to its discharge, after the Yangtze River. Its headwaters drain the YunGui Plateau, where altitude is approximately 2000 meters. Geologically, the carbonate rocks are widely spread in the river drainage basin, which covers an area of about 0.17xl06 km2, i.e., 39% of the whole drainage basin. This study focuses on the chemistry of the Xijiang river system and constitutes the first geochemical investigation into major and trace elements concentrations for both suspended and dissolved loads of this river and its main tributaries, and Sr isotopic composition of the dissolved load is also investigated, in order to determine both chemical weathering and mechanical erosion rates. As compared with the other large rivers of the world, the Xijiang River is characterized by higher major element concentration. The dissolved major cations average 1.17, 0.33, 0.15, and 0.04 mmol I"1 for Ca, Mg, Na, and K, respectively. The total cation concentrations (TZ+) in these rivers vary between 2.2 and 4.4 meq I'1. The high concentration of Ca and Mg, high (Ca+Mg)/(Na+K) ratio (7.9), enormous alkalinity and low dissolved SiO2/HCO3 ratio (0.05) in river waters reveal the importance of carbonate weathering and relatively weak silicate weathering over the river drainage basin. The major elements in river water, such as the alkalis and alkaline-earths, are of different origins: from rain water, silicate weathering, carbonate and evaporite weathering. A mixing model based on mass budget equation is used in this study, which allows the proportions of each element derived from the different source to be calculated. The carbonate weathering is the main source of these elements in the Xijiang drainage basin. The contribution of rainwater, especially for Na, reaches to approximately 50% in some tributaries. Dissolved elemental concentration of the river waters are corrected for rain inputs (mainly oceanic salts), the elemental concentrations derived from the different rock weathering are calculated. As a consequence, silicate, carbonate and total rock weathering rates, together with the consumption rates of atmospheric CO2 by weathering of each of these lithologies have been estimated. They provide specific chemical erosion rates varying between 5.1~17.8 t/km2/yr for silicate, 95.5~157.2 t/km2/yr for carbonate, and 100.6-169.1 t/km2/yr for total rock, respectively. CO2 consumptions by silicate and carbonate weathering approach 13><109and 270.5x10 mol/yr. Mechanical denudation rates deduced from the multi-year average of suspended load concentrations range from 92-874 t/km2/yr. The high denudation rates are mainly attributable to high relief and heavy rainfall, and acid rain is very frequent in the drainage basin, may exceed 50% and the pH value of rainwater may be <4.0, result from SO2 pollution in the atmosphere, results in the dissolution of carbonates and aluminosilicates and hence accelerates the chemical erosion rate. The compositions of minerals and elements of suspended particulate matter are also investigated. The most soluble elements (e.g. Ca, Na, Sr, Mg) are strongly depleted in the suspended phase with respect to upper continent crust, which reflects the high intensity of rock weathering in the drainage basin. Some elements (e.g. Pb, Cu, Co, Cr) show positive anomalies, Pb/Th ratios in suspended matter approach 7 times (Liu Jiang) to 10 times (Nanpan Jiang) the crustal value. The enrichment of these elements in suspended matter reflects the intensity both of anthropogenic pollution and adsorption processes onto particles. The contents of the soluble fraction of rare earth elements (REE) in the river are low, and REE mainly reside in particulate phase. In dissolved phase, the PAAS-normalized distribution patterns show significant HREE enrichment with (La/Yb) SN=0.26~0.94 and Ce depletion with (Ce/Ce*) SN=0.31-0.98, and the most pronounced negative Ce anomalies occur in rivers of high pH. In the suspended phase, the rivers have LREE-enriched patterns relative to PAAS, with (La/Yb) SN=1 -00-1 .40. The results suggest that pH is a major factor controlling both the absolute abundances of REE in solution and the fractionation of REE of dissolved phase. Ce depletion in river waters with high pH values results probably from both preferential removal of Ce onto Fe-Mn oxide coating of particles and CeC^ sedimentation. This process is known to occur in the marine environment and may also occur in high pH rivers. Positive correlations are also observed between La/Yb ratio and DOC, HCO3", PO4", suggesting that colloids and (or) adsorption processes play an important role in the control of these elements.