Study on the processing method of nighttime CO<Subscript>2</Subscript> eddy covariance flux data in ChinaFLUX

At present, using Eddy Covariance (EC) method to estimate the “true value” of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO₂ flux data) appearing in post-processing CO₂ flux data. The ratio of effective and reliable nighttime EC CO₂ flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO₂ flux data and its effect analysis on estimating CO₂ flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO₂ flux using EC method, and introduce the general theory and method for processing nighttime CO₂ flux data. By analyzing the relationship between nighttime CO₂ flux and air fraction velocity u*, we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u*_c) for screening the effective nighttime CO₂ flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO₂ flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.     

Study on the processing method of nighttime CO_2 eddy covariance flux data in ChinaFLUX

At present, using Eddy Covariance (EC) method to estimate the "true value" of carbon sequestration in terrestrial ecosystem arrests more attention. However, one issue is how to solve the uncertainty of observations (especially the nighttime CO2 flux data) appearing in post-processing CO2 flux data. The ratio of effective and reliable nighttime EC CO2 flux data to all nighttime data is relatively low (commonly, less than 50%) for all the long-term and continuous observation stations in the world. Thus, the processing method of nighttime CO2 flux data and its effect analysis on estimating CO2 flux annual sums are very important. In this paper, the authors analyze and discuss the reasons for underestimating nighttime CO2 flux using EC method, and introduce the general theory and method for processing nighttime CO2 flux data. By analyzing the relationship between nighttime CO2 flux and air fraction velocity u., we present an alternate method, Average Values Test (AVT), to determine the thresholds of fraction velocity (u.c) for screening the effective nighttime CO2 flux data. Meanwhile, taking the data observed in Yucheng and Changbai Mountains stations for an example, we analyze and discuss the effects of different methods or parameters on nighttime CO2 flux estimations. Finally, based on the data of part ChinaFLUX stations and related literatures, empirical models of nighttime respiration at different sites in ChinaFLUX are summarized.     

CO2 flux evaluation over the evergreen coniferous and broad-leaved mixed forest in Dinghushan,China

The Dinghushan flux observation site, as one of the four forest sites of ChinaFLUX, aims to acquire long-term measurements of CO₂ flux over a typical southern subtropical evergreen coniferous and broad-leaved mixed forest ecosystem using the open path eddy covariance method. Based on two years of data from 2003 to 2004, the characteristics of temporal variation in CO₂ flux and its response to environmental factors in the forest ecosystem are analyzed. Provided two-dimensional coordinate rotation, WPL correction and quality control, poor energy-balance and underestimation of ecosystem respiration during nighttime implied that there could be a CO₂ leak during the nighttime at the site. Using daytime (PAR > 1.0 μmol⁻¹·m⁻²·s⁻¹) flux data during windy conditions (u* > 0.2 m·s⁻¹), monthly ecosystem respiration (Reco) was derived through the Michaelis-Menten equation modeling the relationship between net ecosystem C0₂ exchange (NEE) and photosynthetically active radiation (PAR). Exponential function was employed to describe the relationship between Reco and soil temperature at 5 cm depth (Ts05), then Reco of both daytime and nighttime was calculated respectively by the function. The major results are: (i) Derived from the Michaelis-Menten equation, the apparent quantum yield (α) was 0.0027±0.0011 mgCO₂·μmol⁻¹ photons, and the maximum photosynthetic assimilation rate (Amax) was 1.102±0.288 mgCO₂·m⁻²·s⁻¹. Indistinctive seasonal variation of α or Amax was consistent with weak seasonal dynamics of leaf area index (LAf) in such a lower subtropical evergreen mixed forest, (ii) Monthly accumulated Reco was estimated as 95.3±21.1 gC·m⁻²mon⁻¹, accounting for about 68% of the gross primary product (GPP). Monthly accumulated WEE was estimated as −43.2±29.6 gC·m⁻²·mon⁻¹. The forest ecosystem acted as carbon sink all year round without any seasonal carbon efflux period. Annual NEE of 2003 and 2004 was estimated as −563.0 and −441.2 gC·m⁻²·a⁻¹ respectively, accounting for about 32% of GPP.

     

Simulating CO2 flux of three different ecosystems in ChinaFLUX based on artificial neural networks

The nonlinearity of the relationship between CO₂ flux and other micrometeorological variables flux parameters limits the applicability of carbon flux models to accurately estimate the flux dynamics. However, the need for carbon dioxide (CO₂) estimations covering larger areas and the limitations of the point eddy covariance technique to address this requirement necessitates the modeling of CO₂ flux from other micrometeorological variables. Artificial neural networks (ANN) are used because of their power to fit highly nonlinear relations between input and output variables without explaining the nature of the phenomena. This paper applied a multilayer perception ANN technique with error back propagation algorithm to simulate CO₂ flux on three different ecosystems (forest, grassland and cropland) in ChinaFLUX. Energy flux (net radiation, latent heat, sensible heat and soil heat flux) and temperature (air and soil) and soil moisture were used to train the ANN and predict the CO₂ flux. Diurnal half-hourly fluxes data of observations from June to August in 2003 were divided into training, validating and testing. Results of the CO₂ flux simulation show that the technique can successfully predict the observed values with R ² value between 0.75 and 0.866. It is also found that the soil moisture could not improve the simulative accuracy without water stress. The analysis of the contribution of input variables in ANN shows that the ANN is not a black box model, it can tell us about the controlling parameters of NEE in different ecosystems and micrometeorological environment. The results indicate the ANN is not only a reliable, efficient technique to estimate regional or global CO₂ flux from point measurements and understand the spatiotemporal budget of the CO₂ fluxes, but also can identify the relations between the CO₂ flux and micrometeorological variables.

     

CO<Subscript>2</Subscript> flux evaluation over the evergreen coniferous and broad-leaved mixed forest in Dinghushan,China

The Dinghushan flux observation site, as one of the four forest sites of ChinaFLUX, aims to acquire long-term measurements of CO₂ flux over a typical southern subtropical evergreen coniferous and broad-leaved mixed forest ecosystem using the open path eddy covariance method. Based on two years of data from 2003 to 2004, the characteristics of temporal variation in CO₂ flux and its response to environmental factors in the forest ecosystem are analyzed. Provided two-dimensional coordinate rotation, WPL correction and quality control, poor energy-balance and underestimation of ecosystem respiration during nighttime implied that there could be a CO₂ leak during the nighttime at the site. Using daytime (PAR > 1.0 μmol^?1·m^?2·s^?1) flux data during windy conditions (u* > 0.2 m·s^?1), monthly ecosystem respiration (Reco) was derived through the Michaelis-Menten equation modeling the relationship between net ecosystem C0₂ exchange (NEE) and photosynthetically active radiation (PAR). Exponential function was employed to describe the relationship between Reco and soil temperature at 5 cm depth (Ts05), then Reco of both daytime and nighttime was calculated respectively by the function. The major results are: (i) Derived from the Michaelis-Menten equation, the apparent quantum yield (α) was 0.0027±0.0011 mgCO₂·μmol^?1 photons, and the maximum photosynthetic assimilation rate (Amax) was 1.102±0.288 mgCO₂·m^?2·s^?1. Indistinctive seasonal variation of α or Amax was consistent with weak seasonal dynamics of leaf area index (LAf) in such a lower subtropical evergreen mixed forest, (ii) Monthly accumulated Reco was estimated as 95.3±21.1 gC·m^?2mon^?1, accounting for about 68% of the gross primary product (GPP). Monthly accumulated WEE was estimated as ?43.2±29.6 gC·m^?2·mon^?1. The forest ecosystem acted as carbon sink all year round without any seasonal carbon efflux period. Annual NEE of 2003 and 2004 was estimated as ?563.0 and ?441.2 gC·m^?2·a^?1 respectively, accounting for about 32% of GPP.     

Eddy covariance imaging of diffuse volcanic CO<Subscript>2</Subscript> emissions at Mammoth Mountain,CA, USA

Use of eddy covariance (EC) techniques to map the spatial distribution of diffuse volcanic CO₂ fluxes and quantify CO₂ emission rate was tested at the Horseshoe Lake tree-kill area on Mammoth Mountain, California, USA. EC measurements of CO₂ flux were made during September–October 2010 and ranged from 85 to 1,766 g m⁻² day⁻¹. Comparative maps of soil CO₂ flux were simulated and CO₂ emission rates estimated from three accumulation chamber (AC) CO₂ flux surveys. Least-squares inversion of measured eddy covariance CO₂ fluxes and corresponding modeled source weight functions recovered 58–77% of the CO₂ emission rates estimated based on simulated AC soil CO₂ fluxes. Spatial distributions of modeled surface CO₂ fluxes based on EC and AC observations showed moderate to good correspondence (R ² = 0.36 to 0.70). Results provide a framework for automated monitoring of volcanic CO₂ emissions over relatively large areas.     

Diurnal variation of ozone flux over corn field in Northwestern Shandong Plain of China

High concentration ground-level ozone(O3)has adverse effects on plant growth and photosynthesis.Compared to the O3concentration-based index,the O3flux-based(especially stomatal O3uptake)index has been considered the better criterion for assessing the impact of ozone on vegetation and ecosystems.This paper reports on a study of O3flux using the eddy covariance technique over a corn field in the Northwestern Shandong Plain of China.Diurnal variation of atmospheric O3concentration,deposition velocity and flux,and their relationships to environmental factors are analyzed.The results show that:(1)During the observation period(9 August–28 September,2011),there was a strong diurnal variation of O3concentration,with low(16.5 nL L?1)and high(60.1 nL L?1)O3mean concentrations observed around 6:30 and 16:00,respectively.Mean O3concentrations during daytime(6:00–18:00)and nighttime(18:00–6:00)were 39.8±23.1 and 20.7±14.1 nL L?1(mean±std),respectively.The maximum observed concentration was 97.5 nL L?1.The concentration was mainly affected by solar radiation and air temperature.(2)Whether daytime or nighttime,ground-level O3flux is always downward.The diurnal course of mean deposition velocity was divided into 4 phases:a low and stable process during nighttime,fast increasing in early morning,relatively large and steady changes around noon,and quickly decreasing in later afternoon.Daytime and nighttime mean deposition velocities were 0.29 and 0.09 cm s?1,respectively.The maximum deposition velocity was 0.81 cm s?1.The magnitude of deposition velocity was influenced by the corn growth period,and its diurnal variation was significantly correlated with global radiation and relative humidity.(3)O3flux was affected by variations of both O3concentration and deposition velocity,with mean O3fluxes-317.7 and?70.2 ng m?2s?1during daytime and nighttime,respectively.There was strong correlation between O3flux and CO2flux or latent heat flux.By comparing the deposition velocities of daytime and nighttime,we infer that stomatal uptake was probably the main sink of ground-level O3.     

Estimates of Horizontal Groundwater Flow Velocities in Boreholes

Currently, monitoring tools can be deployed in observation boreholes to better assess groundwater flow, flux of dissolved contaminants and their mass discharge in an aquifer. The relationship between horizontal water velocity in observation boreholes and Darcy fluxes in the surrounding aquifer has been studied for natural flow conditions (i.e., no pumping). Interpretation of measurements taken with dilution tests, the colloidal borescope, the Heat Pulse Flowmeter, and other techniques require the conversion of observed borehole velocity u to aquifer Darcy flux q_∞ . This conversion is typically done through a proportionality factor α = u/q_∞ . In experimental studies as well as in theoretical developments, reported values of α vary almost three orders of magnitude (from 0.5 to 10). This large variability in reported values of α could be explained by: (1) unclear distinction between Darcy flux and water seepage velocity, (2) unclear definition of water velocity in the borehole, (3) effects of well screen and the presence of the measurement device itself on the observable velocities, and (4) hydraulic conditions in the borehole annulus. We address (1), (2) from a conceptual/theoretical perspective, and (3) by means of numerical simulations. We show that issue (1) in low porosity aquifers can yield to order-of-magnitude discrepancies in estimates of q_∞ ; (2) may result in discrepancies of up to 50%, and (3) can cause differences up to 20% of water velocity in the borehole void space compared to the theoretical case of an open borehole.     

Recent progress and future directions of ChinaFLUX

The eddy covariance technique has emerged as an important tool to directly measure carbon dioxide, water vapor and heat fluxes between the terrestrial ecosystem and the atmosphere after a long history of fundamental research and technological developments. With the realization of regional networks of flux measurements in North American, European, Asia, Brazil, Australia and Africa, a global-scale network of micrometeorological flux measurement (FLUXNET) was established in 1998. FLUXNET has made great progresses in investigating the environmental mechanisms controlling carbon and water cycles, quantifying spatial-temporal patterns of carbon budget and seeking the “missing carbon sink” in global terrestrial ecosystems in the past ten years. The global-scale flux measurement also built a platform for international communication in the fields of resource, ecology and environment sciences. With the continuous development of flux research, FLUXNET will introduce and explore new techniques to extend the application fields of flux measurement and to answer questions in the fields of bio-geography, eco-hydrology, meteorology, climate change, remote sensing and modeling with eddy covariance flux data. As an important part of FLUXNET, ChinaFLUX has made significant progresses in the past three years on the methodology and technique of eddy covariance flux measurement, on the responses of CO₂ and H₂O exchange between the terrestrial ecosystem and the atmosphere to environmental change, and on flux modeling development. Results showed that the major forests on the North-South Transect of Eastern China (NSTEC) were all carbon sinks during 2003 to 2005, and the alpine meadows on the Tibet Plateau were also small carbon sinks. However, the reserved natural grassland, Leymus chinensis steppe in Inner Mongolia, was a carbon source. On a regional scale, temperature and precipitation are the primary climatic factors that determined the carbon balance in major terrestrial ecosystems in China. Finally, the current research emphasis and future directions of ChinaFLUX were presented. By combining flux network and terrestrial transect, ChinaFLUX will develop integrated research with multi-scale, multi-process, multi-subject observations, placing emphasis on the mechanism and coupling relationships between water, carbon and nitrogen cycles in terrestrial ecosystems.     

A preliminary study for spatial representiveness of flux observation at ChinaFLUX sites

The results of eddy covariance observation system could represent the physical process at certain area of the surface. Thus point-to-area representativeness was of primary interest in flux observation. This research presents a preliminary study for flux observation at ChinaFLUX sites by the use of observation data and Flux Source Area Model (FSAM). Results show that the footprint expands and is further away from flux tower when atmosphere becomes more stable, the observation height increases, or the surfaces become smoother. This suggests that the area represented by the flux observation becomes larger. The distances from the reference point to the maximum point S _max and the minimum point x ₁ of source weight function (D _max and D _min, respectively) can be influenced by atmosphere stability which becomes longer when atmosphere is more stable. For more rough surfaces and lower observation point D _max and D _min become shorter. This research gives the footprint at level P=90% at ChinaFLUX sites at different atmosphere stability. The preliminary results of spatial representiveness at ChinaFLUX sites were given based on the dominant wind direction and footprint response to various factors. The study also provides some theoretical basis for data quality control and evaluating data uncertainty.     

Morphological analysis of the winter nighttime ionosphere at mid-latitudes of the Asian region

We present a study of peculiarities of the winter nighttime maximum in the critical frequencies f ₀ F2 at mid-latitudes of the Asian region. The data of stations located at different longitudes and close latitudes have been used in the analysis: Novosibirsk (54.8°N, 83.2°E), Irkutsk (52.5°N, 104.0°E), and Khabarovsk (48.5°N, 135.1°E). It has been found that the nighttime maximum in f ₀ F2 is observed after midnight (∼0200–0400 LT) and is a stable feature of the quiet ionosphere from the middle of October to the middle of March at low solar activity (SA) at all analyzed stations. This interval decreases with increasing SA. The difference between the maximal and minimal f ₀ F2 values in nighttime hours is the largest in December–January, and its amplitude is almost independent of SA. Variations in the critical frequency of the h _m F2 layer are inversely related to those in the height of the maximum. We have studied periods when the difference between the daytime and nighttime values of f ₀ F2 is less than 2 MHz. The intervals of observations of such events at different longitudes do not coincide. No dependence of the winter nighttime maximum amplitude on magnetic activity has been found.     

Environmental influences on soil CO2 degassing at Furnas and Fogo volcanoes (São Miguel Island,Azores archipelago)

Since October 2001, four soil CO₂ flux stations were installed in the island of São Miguel (Azores archipelago), at Fogo and Furnas quiescent central volcanoes. These stations perform measurements by the accumulation chamber method and, as the gas flux may be influenced by external variables, the stations are equipped with several meteorological sensors. Multivariate regression analysis applied to the large datasets obtained allowed observing that the meteorological variables may influence the soil CO₂ flux oscillations from 18% to 50.5% at the different monitoring sites. Additionally, it was observed that meteorological variables (mainly soil water content, barometric pressure, wind speed and rainfall) play a different role in the control of the gas flux, depending on the selected monitoring site and may cause significant short-term (spike-like) fluctuations. These divergences may be potentially explained by the porosity and hydraulic conductivity of the soils, topographic effects, drainage area and different exposure of the monitoring sites to the weather conditions. Seasonal effects are responsible for long-term oscillations on the gas flux.     

Carbon dioxide diffuse emission from the soil: ten years of observations at Vesuvio and Campi Flegrei (Pozzuoli), and linkages with volcanic activity

Carbon dioxide flux from the soil is regularly monitored in selected areas of Vesuvio and Solfatara (Campi Flegrei, Pozzuoli) with the twofold aim of i) monitoring spatial and temporal variations of the degassing process and ii) investigating if the surface phenomena could provide information about the processes occurring at depth. At present, the surveyed areas include 15 fixed points around the rim of Vesuvio and 71 fixed points in the floor of Solfatara crater. Soil CO₂ flux has been measured since 1998, at least once a month, in both areas. In addition, two automatic permanent stations, located at Vesuvio and Solfatara, measure the CO₂ flux and some environmental parameters that can potentially influence the CO₂ diffuse degassing. Series acquired by continuous stations are characterized by an annual periodicity that is related to the typical periodicities of some meteorological parameters. Conversely, series of CO₂ flux data arising from periodic measurements over the arrays of Vesuvio and Solfatara are less dependent on external factors such as meteorological parameters, local soil properties (porosity, hydraulic conductivity) and topographic effects (high or low ground). Therefore we argue that the long-term trend of this signal contains the “best” possible representation of the endogenous signal related to the upflow of deep hydrothermal fluids.     

Near-bed mass flux profiles in aeolian sand transport: high-resolution measurements in a wind tunnel

Vertical profiles of the streamwise mass flux of blown sand in the near-bed (< 17 mm) region are analysed from high-resolution measurements made using an optical sensor in a wind tunnel. This analysis is complemented by detailed measurements of mass flux and mean velocity profiles throughout the boundary layer depth (0·17 m) using passive, chambered sand traps of small dimensions and armoured thermal anemometers, respectively. The data permit a preliminary analysis of the relations between the observed forms of the profiles of near-bed fluid stress and horizontal mass flux within a carefully conditioned boundary layer. Profiles of mass flux density are found to be characterized by three regions of differing gradient with transitions at about 2 mm and 19 mm above the bed. The exponential decay of mass flux with height is confirmed for elevations above 19 mm, and when plotted as a function of u_*²/g (a parameter of mean vertical trajectory height in saltation), the gradient of mass flux in this region scales with the wake-corrected friction velocity (u), where u > 0·30 m s⁻¹. A separate near-bed region of more intense transport below 19 mm is identified which carries 80 per cent of the total mass flux. This region is evident in some previous field and wind tunnel data but not in profiles simulated by numerical models. Ventilated passive sand traps underestimate mass flux in this region by 37 per cent. At slow or moderate wind speeds a third significant region below 2 mm is observed. These regions are likely to be related to grain populations in successive saltation, low-energy ejections and intermittent bed contact, respectively. Optical measurements reveal locally high grain concentrations at some elevations below 5 mm; these heights scale with transport rate, mass flux gradient and wind speed. Copyright © 1999 John Wiley & Sons, Ltd.     

The response of CO2 flux to rain pulses at a saline desert

As the substantial component of the ecosystem respiration, soil CO₂ flux is strongly influenced by infrequent and unpredictable precipitation in arid region. In the current study, we investigated the response of soil CO₂ flux to rain pulses at a saline desert in western China. Soil CO₂ flux was measured continuously during the whole growing season of 2009 at six sites. We found that there were remarkable changes in amplitude or diurnal patterns of soil CO₂ flux induced by rainfall events: from bimodal before rain to a single peak after that. Further analysis indicated that there is a significant linear relationship (P < 0.001) between soil CO₂ flux and soil temperature (T_soil). However, a hysteresis between the waveform of diurnal course of CO₂ flux and T_soil was observed: with soil CO₂ flux always peaked earlier than T_soil. Furthermore, a double exponential decay function was fitted to the soil CO₂ flux after rainfall, and total carbon (C) releases were estimated by numerical integration for rainfall events. The relative enhancement and total C release, in association with the rain pulses, was linearly related to the amount of precipitation. According to the size and frequency of rainfall events, the total amount of C release induced by rain pulses was computed as much as 7.88 g C·m^–2 in 2009, equivalent to 10.25% of gross primary production. These results indicated that rain pulses played a significant role in the carbon budget of this saline desert ecosystem, and the size of them was a good indicator of rain‐induced flux enhancement. Copyright © 2012 John Wiley & Sons, Ltd.     

CO2 absorption of sandy soil induced by rainfall pulses in a desert ecosystem

Soil CO₂ flux is strongly influenced by precipitation in many ecosystem types, yet knowledge of the effects of precipitation on soil CO₂ flux in semi‐arid desert ecosystems remains insufficient, particularly for sandy soils. To address this, we investigated the response of sandy soil CO₂ flux to rainfall pulses in a desert ecosystem in northern China during August–September 2011. Significant changes (P < 0.05) were found in diel patterns of soil CO₂ flux induced by small (2.1 mm), moderate (12.4 mm) and large (19.7 mm) precipitation events. Further analysis indicated that rainfall pulses modified the response of soil CO₂ flux to soil temperature, including hysteresis between soil CO₂ flux and soil temperature, with F_s higher when T_s was increasing than when T_s was decreasing, and the linear relationship between them. Moreover, our results showed that rainfall could result in absorption of atmospheric CO₂ by soil, possibly owing to mass flow of CO₂ induced by a gradient of gas pressure between atmosphere and soil. After each precipitation event, soil CO₂ flux recovered exponentially to pre‐rainfall levels with time, with the recovery times exhibiting a positive correlation with precipitation amount. On the basis of the amounts of precipitation that occurred at our site during the measurement period (August–September), the accumulated rain‐induced carbon absorption evaluated for rainy days was 1.068 g C m^?2; this corresponds approximately to 0.5–2.1% of the net primary production of a typical desert ecosystem. Thus, our results suggest that rainfall pulses can strongly influence carbon fluxes in desert ecosystems. Copyright © 2014 John Wiley & Sons, Ltd.     

Seasonal change in CO2 and H2O exchange between grassland and atmosphere

The seasonal change in CO₂ flux over an artificial grassland was analyzed from the ecological and meteorological point of view. This grassland contains C₃ and C₄ plants; the three dominant species belonging to the Gramineae; Festuca elatior (C₃) dominated in early spring, and Imperata cylindrica (C₄) and Andropogon virginicus (C₄) grew during early summer and became dominant in mid-summer. CO₂ flux was measured by the gradient method, and the routinely observed data for the surface-heat budget were used to analyze the CO₂ and H₂O exchange between the grassland and atmosphere. From August to October in 1993, CO₂ flux was reduced to around half under the same solar-radiation conditions, while H₂O flux decreased 20% during the same period. The monthly values of water use efficiency, i.e., ratio of CO₂ flux to H₂O flux decreased from 5.8 to 3.3 mg CO₂/g H₂O from August to October, the Bowen ratio increased from 0.20 to 0.30, and the ratio of the bulk latent heat transfer coefficient C_E to the sensible heat transfer coefficient C_H was maintained around 0.40-0.50. The increase in the Bowen ratio was explained by the decrease in air temperature from 22.3 °C in August to 16.6 °C in October without considering biological effects such as stomatal closure on the individual leaves. The nearly constant C_E/C_H ratios suggested that the contribution ratio of canopy resistance to aerodynamic resistance did not change markedly, although the meteorological conditions changed seasonally. The decrease in the water use efficiency, however, suggested that the photosynthetic rate decreased for individual leaves from August to October under the same radiation conditions. Diurnal variations of CO₂ exchange were simulated by the multi-layer canopy model taking into account the differences in the stomatal conductance and photosynthetic pathway between C₃ and C₄ plants. The results suggested that C₄ plants played a major role in the CO₂ exchange in August, the contribution of C₄ plants decreased in September, and daily variations of CO₂ exchange were mainly due to C₃ plants in October. The results also suggested that the decrease in the net canopy CO₂ exchange from August to October was induced partly by the decrease of net photosynthesis on the individual leaves in both C₄ and C₃ plants, which could be due to aging of the leaves.     

Continuous soil CO<Superscript>2</Superscript> and discrete plume SO<Subscript>2</Subscript> measurements at Mt. Etna (Italy) during 1997–2000: a contribution to volcano monitoring

Continuous monitoring of soil CO₂ dynamic concentration (which is proportional to the CO₂ flux through the soil) was carried out at a peripheral site of Mt. Etna during the period November 1997–September 2000 using an automated station. The acquired data were compared with SO₂ flux from the summit craters measured two to three times a week during the same period. The high frequency of data acquisition with both methods allowed us to analyze in detail the time variations of both parameters. Anomalous high values of soil CO₂ dynamic concentration always preceded periods of increased flux of plume SO₂, and these in turn were followed by periods of summit eruptions. The variations were modeled in terms of gas efflux increase due to magma ascent to shallow depth and its consequent depressurization and degassing. This model is supported by data from other geophysical and volcanological parameters. The rates of increase both of soil CO₂ dynamic concentration and of plume SO₂ flux are interpreted to be positively correlated both to the velocity of magma ascent within the volcano and to lava effusion rate once magma is erupted at the surface. Low rates of the increase were recorded before the nine-month-long 1999 subterminal eruption. Higher rates of increase were observed before the violent summit eruption of September-November 1999, and the highest rates were observed during shorter and very frequent spike-like anomalies that preceded the sequence of short-lived but very violent summit eruptions that started in late January 2000 and continued until late June of the same year. Furthermore, the time interval between the peaks of CO₂ and SO₂ in a single sequence of gas anomalies is likely to be controlled by magma ascent velocity.Editorial responsibility: H. Shinohara     

Shear wave velocity prediction during CO<Subscript>2</Subscript>-EOR and sequestration in the Gao89 well block of the Shengli Oilfield

Shear-wave velocity is a key parameter for calibrating monitoring time-lapse 4D seismic data during CO₂-EOR (Enhanced Oil Recovery) and CO₂ sequestration. However, actual S-wave velocity data are lacking, especially in 4D data for CO₂ sequestration because wells are closed after the CO₂ injection and seismic monitoring is continued but no well log data are acquired. When CO₂ is injected into a reservoir, the pressure and saturation of the reservoirs change as well as the elastic parameters of the reservoir rocks. We propose a method to predict the S-wave velocity in reservoirs at different pressures and porosities based on the Hertz–Mindlin and Gassmann equations. Because the coordination number is unknown in the Hertz–Mindlin equation, we propose a new method to predict it. Thus, we use data at different CO₂ injection stages in the Gao89 well block, Shengli Oilfield. First, the sand and mud beds are separated based on the structural characteristics of the thin sand beds and then the S-wave velocity as a function of reservoir pressure and porosity is calculated. Finally, synthetic seismic seismograms are generated based on the predicted P- and S-wave velocities at different stages of CO₂ injection.     

Dynamics of carbon dioxide emissions,crystallization, and magma ascent: hypotheses,theory, and applications to volcano monitoring at Mount St. Helens

 Measurements of CO₂ fluxes from open-vent volcanos are rare, yet may offer special capabilities for monitoring volcanos and forecasting activity. The measured fluxes of CO₂ and SO₂ from Mount St. Helens decreased from July through November 1980, but the record includes variations of CO₂/SO₂ in the emitted gas and episodes of greatly increased fluxes of CO₂. We propose that the CO₂ flux variations reflect two gas components: (a) a component whose flux decreased in proportion to 1/ √t with a CO₂/SO₂ mass ratio of 1.7, and (b) a residual flux of CO₂ consisting of short-lived, large peaks with a CO₂/SO₂ mass ratio of 15. We propose two hypotheses: (a) the 1/ √t dependence was generated by crystallization in a deep magma body at rates governed by diffusion-limited heat transfer, and (b) the gas component with the higher CO₂/SO₂ was released from ascending magma, which replenished the same magma body. The separation of the total CO₂ flux into contributions from known processes permits quantitative inferences about the replenishment and crystallization rates of open-system magma bodies beneath volcanos. The flux separations obtained by using two gas sources with distinct CO₂/SO₂ ratios and a peak minus background approach to obtain the CO₂ contributions from an intermittent source and a continuously emitting source are similar. The flux separation results support the hypothesis that the second component was generated by episodic magma ascent and replenishment of the magma body. The diffusion-limited crystallization hypothesis is supported by the decay of minimum CO₂ and SO₂ fluxes with 1/ √t after 1 July 1980. We infer that the magma body at Mount St. Helens was replenished at an average rate (2.8×10⁶m³d^–1) which varied by less than 5% during July, August, and September 1980. The magma body volume (2.4–3.0 km³) in early 1982 was estimated by integrating a crystallization rate function inferred from CO₂ fluxes to maximum times (20±4 years) estimated from the increase of sample crystallinity with time. These new volcanic gas flux separation methods and the existence of relations among the CO₂ flux, crystallization rates, and magma body replenishment rates yield new information about the dynamics of an open-vent, replenished magma body. Received: 15 February 1995 / Accepted: 30 March 1996