Soil H<Subscript>2</Subscript> and CO<Subscript>2</Subscript> Surveys at Several Active Faults in Japan

Soil H₂ and CO₂ surveys were carried out along seven active faults and around the aftershock region of the 2000 Tottori-ken Seibu earthquake in Japan. Diffuse CO₂ effluxes were also measured along one fault and around the 2000 aftershock region. The results show highly variable H₂ concentration in space and time and it seems that the maximum H₂ concentration at each active fault correlates with fault activity as exemplified by the time of the latest big earthquakes. Even though observed H₂ concentrations in four faults were markedly lower than those collected previously in the latter half of the 1970s, it is evident that the higher H₂ concentrations in this study are due to the addition of the fault gases. Comparing the chemical composition of trapped gases (H₂: 5–20% and CO₂/H₂: 0.5–12) in fractured rocks of drill cores bored at the Nojima fault, a soil gas sample with the highest H₂ concentration showed large amounts of the trapped fault gas, diluted with atmospheric component. The profile experiment across a fracture zone at the Yamasaki fault showed higher H₂ concentrations and lower CO₂/H₂ ratios as was observed in soil gas from the fracture zone. A few days after the 2000 Tottori-kei Seibu earthquake, no CO₂ effluxes related to the occurrence of earthquakes were observed at the aftershock region. However, only above the epicenter zone, relatively high H₂ concentrations in soil gases were observed.     

Diffuse CO<Subscript>2</Subscript> degassing at Vesuvio,Italy

At Vesuvio, a significant fraction of the rising hydrothermal–volcanic fluids is subjected to a condensation and separation process producing a CO₂–rich gas phase, mainly expulsed through soil diffuse degassing from well defined areas called diffuse degassing structures (DDS), and a liquid phase that flows towards the outer part of the volcanic cone. A large amount of thermal energy is associated with the steam condensation process and subsequent cooling of the liquid phase. The total amount of volcanic–hydrothermal CO₂ discharged through diffuse degassing has been computed through a sequential Gaussian simulation (sGs) approach based on several hundred accumulation chamber measurements and, at the time of the survey, amounted to 151 t d^–1. The steam associated with the CO₂ output, computed assuming that the original H₂O/CO₂ ratio of hydrothermal fluids is preserved in fumarolic effluents, is 553 t d^–1, and the energy produced by the steam condensation and cooling of the liquid phase is 1.47×10¹² J d^–1 (17 MW). The location of the CO₂ and temperature anomalies show that most of the gas is discharged from the inner part of the crater and suggests that crater morphology and local stratigraphy exert strong control on CO₂ degassing and subsurface steam condensation. The amounts of gas and energy released by Vesuvio are comparable to those released by other volcanic degassing areas of the world and their estimates, through periodic surveys of soil CO₂ flux, can constitute a useful tool to monitor volcanic activity.Editorial responsibility: H. Shinohara     

Variations of near-surface atmospheric CO2 and H2O concentrations during summer on Muztagata

Expeditions during the summers of 2002 and 2003 implemented continuous monitoring of near-surface (2 m height) atmospheric CO₂ and H₂O concentrations at the 4500 m elevation on Muztagata. The resultant data sets reveal a slight decrease of CO₂ concentrations (of about 5 μmol·mol^-1) and changes in the diurnal variations from the end of June to the middle August. The daily maximum CO₂ concentrations occur between 02:30-05:30 AM (local time) and the minimum levels occur between 12:00-15:30 PM. The atmospheric CO₂ concentrations in the summer of 2002 were around 5 μmol·mol^-1 lower than those during the same period of 2003, whereas the diurnal amplitude was higher. In contrast, we found that the daily mean atmospheric H₂O content in 2003 was much lower than that in 2002 and there exists a striking negative correlation between CO₂ and H₂O concentrations. We therefore suggest that the near-surface atmospheric CO₂ concentration is affected not only by photosynthesis and respiration, but also by the air H₂O content in the glaciated region around Muztagata.     

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.     

Seasonal variations and mechanism for environmental control of NEE of CO<Subscript>2</Subscript> concerning the <Emphasis Type="Italic">Potentilla fruticosa</Emphasis> in alpine shrub meadow of Qinghai-Tibet Plateau

The study by the eddy covariance technique in the alpine shrub meadow of the Qinghai-Tibet Plateau in 2003 and 2004 showed that the net ecosystem carbon dioxide exchange (NEE) exhibited noticeable diurnal and annual variations, with more distinct daily changes during the warmer seasons. The CO₂ emission of the shrub ecosystem culminated in April and September while the CO₂ absorption capacity reached a maximum in July and August. The absorbed carbon dioxide during the two consecutive years was 231.4 and 274.8 g CO₂·m^?2 respectively, yielding an average of 253.1 gCO₂·m^?2 per year: that accounts for a large proportion of absorbed CO₂ in the region. Obviously, the diurnal carbon flux was negatively related to temperature, radiation and other atmospheric factors. Still, minute discrepancies in kurtosis and duration of carbon emission/absorption were detected between 2003 and 2004. It was found that the CO₂ flux in the daytime was similarly affected by photosynthetic photon flux density in both years. Temperature appears to be the most important determinant of CO₂ flux: specifically, the high temperature during the plant growing season inhibits the carbon absorption capacity. One potential explanation is that soil respiration is enhanced under such condition. Analysis of biomass revealed that the annual net carbon fixed capacity of aboveground and belowground biomass was 544.0 in 2003 and 559.4 g C·m^?2 in 2004, which coincided with the NEE absorption capacity (63.1 g C·m^?2 in 2003 and 74.9 g C·m^?2 in 2004) in the corresponding plant growing season.     

Physics and Seismic Modeling for Monitoring CO<Subscript>2</Subscript> Storage

We present a new petro-elastical and numerical-simulation methodology to compute synthetic seismograms for reservoirs subject to CO₂ sequestration. The petro-elastical equations model the seismic properties of reservoir rocks saturated with CO₂, methane, oil and brine. The gas properties are obtained from the van der Waals equation and we take into account the absorption of gas by oil and brine, as a function of the in situ pore pressure and temperature. The dry-rock bulk and shear moduli can be obtained either by calibration from real data or by using rock-physics models based on the Hertz-Mindlin and Hashin-Shtrikman theories. Mesoscopic attenuation due to fluids effects is quantified by using White's model of patchy saturation, and the wet-rock velocities are calculated with Gassmann equations by using an effective fluid modulus to describe the velocities predicted by White's model. The simulations are performed with a poro-viscoelastic modeling code based on Biot's theory, where viscoelasticity is described by generalizing the solid/fluid coupling modulus to a relaxation function. Using the pseudo-spectral method, which allows general material variability, a complete and accurate characterization of the reservoir can be obtained. A simulation, that considers the Utsira sand of the North Sea, illustrates the methodology.     

The characteristics of atmospheric CO<Subscript>2</Subscript> concentration variation of four national background stations in China

Atmospheric carbon dioxide is an important kind of greenhouse gas which influences global temperature. Its concentration variation could indicate the distribution of human and natural activities in various regions. Through the non-dispersive infrared method, flask sampling of atmospheric CO₂ concentration was measured weekly at four national background stations including Waliguan, Shangdianzi, Lin’an, and Longfengshan. Based on the data collected from September 2006 to August 2007, along with the Waliguan station’s experience on in situ observational data processing, the selection methods for sampling data through the atmospheric background CO₂ concentration analysis were preliminarily discussed. On the basis of this result, the variation features of the four typical regions’ atmospheric background CO₂ concentration was analyzed for the first time. The results show that the atmospheric CO₂ concentration at Waliguan, Shangdianzi, Lin’an, and Longfengshan is 383.5, 385.9, 387.8, and 384.3 ppm, respectively. During the research period, CO₂ concentration at the Waliguan station changed slightly. However, the CO₂ concentration changed sharply at the Shangdianzi and the Lin’an stations due to the great influence of human activities in the Jingjinji and the Changjiang Delta economic zones, and changed regularly with seasons at Longfengshan station under dual influences of human activities and plant photosynthesis. The results from this study can lay the foundation for more profound studies on atmospheric CO₂ concentration level of different areas in China, and could be used to improve the understanding of carbon source and sink distribution.     

A comparison of atmospheric CO2 concentration GOSAT-based observations and model simulations

Satellite observations of atmospheric CO2 are able to truly capture the variation of global and regional CO2 concentration.The model simulations based on atmospheric transport models can also assess variations of atmospheric CO2 concentrations in a continuous space and time,which is one of approaches for qualitatively and quantitatively studying the atmospheric transport mechanism and spatio-temporal variation of atmospheric CO2 in a global scale.Satellite observations and model simulations of CO2 offer us two different approaches to understand the atmospheric CO2.However,the difference between them has not been comprehensively compared and assessed for revealing the global and regional features of atmospheric CO2.In this study,we compared and assessed the spatio-temporal variation of atmospheric CO2 using two datasets of the column-averaged dry air mole fractions of atmospheric CO2(XCO2)in a year from June 2009 to May 2010,respectively from GOSAT retrievals(V02.xx)and from Goddard Earth Observing System-Chemistry(GEOS-Chem),which is a global 3-D chemistry transport model.In addition to the global comparison,we further compared and analyzed the difference of CO2 between the China land region and the United States(US)land region from two datasets,and demonstrated the reasonability and uncertainty of satellite observations and model simulations.The results show that the XCO2 retrieved from GOSAT is globally lower than GEOS-Chem model simulation by 2 ppm on average,which is close to the validation conclusion for GOSAT by ground measures.This difference of XCO2 between the two datasets,however,changes with the different regions.In China land region,the difference is large,from 0.6 to 5.6 ppm,whereas it is 1.6 to 3.7 ppm in the global land region and 1.4 to 2.7 ppm in the US land region.The goodness of fit test between the two datasets is 0.81 in the US land region,which is higher than that in the global land region(0.67)and China land region(0.68).The analysis results further indicate that the inconsistency of CO2concentration between satellite observations and model simulations in China is larger than that in the US and the globe.This inconsistency is related to the GOSAT retrieval error of CO2 caused by the interference among input parameters of satellite retrieval algorithm,and the uncertainty of driving parameters in GEOS-Chem model.     

Continuous gas monitoring in the West Bohemian earthquake area,Czech Republic: First results

Two stations monitoring concentrations of carbon dioxide and radon in soil gas (Oldřišská and Novy Kostel) and one station monitoring flow of carbon dioxide at a mofette (Soos) have been operated in the area of the West Bohemian earthquake swarms. We present preliminary results obtained on the base of four-year observations. We found that data are not influenced considerably by barometric pressure. Although the CO₂ concentration varies greatly, the long-term trends at stations Oldřišská and Novy Kostel are similar, which indicates that the CO₂ flow is controlled by common geogenic processes. Also temporal trends of CO₂ and Rn concentrations in soil gas at individual stations are analogous. We found diurnal variations of both CO₂ concentration in soil gas and the CO₂ flow at mofettes due to the earth tides. A response to tides of semi-diurnal period is insignificant in CO₂ concentration and only weak in the CO₂ flow. We also examined possible pre-seismic, co-seismic and post-seismic effects of the intensive 2008 earthquake swarm on the CO₂ concentration at Oldřišská and Novy Kostel, and on the CO₂ flow at Soos. However, all potential indications were insignificant and there has not been proven any influence of the swarm on the CO₂ concentration as well as on the CO₂ flow. Nevertheless, a gradual decrease of amplitudes of diurnal variations before the swarm and the lowest amplitudes during the swarm is a noteworthy phenomenon, which might indicate the strain changes of the rock associated with earthquake swarm.     

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.     

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.     

Structure and CO<Subscript>2</Subscript> budget of Merapi volcano during inter-eruptive periods

Soil temperature and gas (CO₂ concentration and flux) have been investigated at Merapi volcano (Indonesia) during two inter-eruptive periods (2002 and 2007). Precise imaging of the summit crater and the spatial pattern of diffuse degassing along a gas traverse on the southern slope are interpreted in terms of summit structure and major caldera organization. The summit area is characterized by decreasing CO₂ concentrations with distance from the 1932 crater rim, down to atmospheric levels at the base of the terminal cone. Similar patterns are measured on any transect down the slopes of the cone. The spatial distribution of soil gas anomalies suggests that soil degassing is controlled by structures identified as concentric historical caldera rims (1932, 1872, and 1768), which have undergone severe hydrothermal self-sealing processes that dramatically lower the permeability and porosity of soils. Temperature and CO₂ flux measurements in soils near the dome display heterogeneous distributions which are consistent with a fracture network identified by previous geophysical studies. These data support the idea that the summit is made of isolated and mobile blocks, whose boundaries are either sealed by depositional processes or used as pathways for significant soil degassing. Within this context, self-sealing both prevents long-distance soil degassing and controls heat and volatile transfers near the dome. A rough estimate of the CO₂ output through soils near the dome is 200–230 t day⁻¹, i.e. 50% of the estimated total gas output from the volcano summit during these quiescent periods. On Merapi’s southern slope, a 2,500 m long CO₂ traverse shows low-amplitude anomalies that fit well with a recently observed electromagnetic anomaly, consistent with a faulted structure related to an ancient avalanche caldera rim. Sub-surface soil permeability is the key parameter that controls the transfer of heat and volatiles within the volcano, allowing its major tectonic architecture to be revealed by soil gas and soil temperature surveys.     

Time-dependent CO<Subscript>2</Subscript> variations in Lake Albano associated with seismic activity

Lake Albano (Alban Hills volcanic complex, Central Italy) is located in a densely populated area near Rome. The deep lake waters have significant dissolved CO₂ concentrations, probably related to sub-lacustrine fluid discharges fed by a pressurized CO₂-rich reservoir. The analytical results of geochemical surveys carried out in 1989–2010 highlight the episodes of CO₂ removal from the lake. The total mass of dissolved CO₂ decreased from ∼5.8 × 10⁷ kg in 1989 to ∼0.5 × 10⁷ kg in 2010, following an exponential decreasing trend. Calculated values of both dissolved inorganic carbon and CO₂ concentrations along the vertical profile of the lake indicate that this decrease is caused by CO₂ release from the epilimnion, at depth <9 m, combined with (1) water circulation at depth <95 m and (2) CO₂ diffusion from the deeper lake layers. According to this model, Lake Albano was affected by a large CO₂ input that coincided with the last important seismic swarm at Alban Hills in 1989, suggesting an intimate relationship between the addition of deep-originated CO₂ to the lake and seismic activity. In the case of a CO₂ degassing event of an order of magnitude larger than the one that occurred in 1989, the deepest part of Lake Albano would become CO₂-saturated, resulting in conditions compatible with the occurrence of a gas outburst. These results reinforce the idea that a sudden CO₂ input into the lake may cause the release of a dense gas cloud, presently representing the major volcanic threat for this densely populated area.     

Influence of small-scale heterogeneity on upward CO2 plume migration in storage aquifers

Recent advancements in experimental techniques allow sub-core-scale heterogeneities to be quantified in a high resolution. Based on the observations of heterogeneity distributions in natural core samples, we perform simulations to study the influence of small-scale heterogeneities on large-scale CO₂ migration during geological storage. We observe that even the heterogeneities at millimeter scale (the scale of a Representative Elementary Volume for sandstones) can affect large-scale buoyancy-driven upward CO₂ migration. For the representative examples we study, ignoring small-scale heterogeneities can lead to an overestimation of the migration speed by a factor of two.To analyze the cause of such overestimation, we introduce a dimensionless heterogeneity factor to characterize different levels of heterogeneity. The influence on CO₂ migration is quantified with respect to a variety of heterogeneity factors, correlation lengths, and fluid viscosity ratios for isotropic and anisotropic media. Our findings suggest that small-scale heterogeneities should not be ignored in core analysis, even for cores that appear relatively homogeneous and do not have distinguishable heterogeneous patterns. In addition, relative-permeability curves measured from core-flood experiments under high-flow-velocity conditions (a common practice to eliminate capillary end-effects) should not be directly used in modeling low-velocity CO₂ migration if small-scale heterogeneities are present.     

Links between global CO<Subscript>2</Subscript> variability and climate anomalies of biomes

The global rate of fossil fuel combustion continues to rise, but the amount of CO₂ accumulating in the atmosphere has not increased accordingly. The causes for this discrepancy are widely debated. Particularly, the location and drivers for the interannual variability of atmospheric CO₂ are highly uncertain. Here we examine links between global atmospheric CO₂ growth rate (CGR) and the climate anomalies of biomes based on (1986–1995) global climate data of ten years and accompanying satellite data sets. Our results show that four biomes, the tropical rainforest, tropical savanna, C₄ grassland and boreal forest, and their responses to climate anomalies, are the major climate-sensitive CO₂ sinks/sources that control the CGR. The nature and magnitude by which these biomes respond to climate anomalies are generally not the same. However, one common influence did emerge from our analysis; the extremely high CGR observed for the one extreme El Niño year was caused by the response of the tropical biomes (rainforest, savanna and C₄ grassland) to temperature.     

Seasonal variation of the <Emphasis Type="Italic">M</Emphasis><Subscript>2</Subscript> tide

The seasonal cycle of the main lunar tidal constituent M ₂ is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M ₂ tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M ₂ tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M ₂ internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M ₂ surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.     

Quantitative assessment of the role of doubled CO<Subscript>2</Subscript> and associated climate change in the vegetation dynamics and hydrological cycle in the Sino-Mongolia arid and semi-arid region

The effects of elevated CO₂ on vegetation dynamics and the hydrological cycle have been widely studied at the site level. However, quantitative assessments of these effects on a regional scale remain a challenge. We conducted numerical simulations to predict the possible responses of vegetation and the hydrological cycle in the Sino-Mongolia arid and semi-arid region (SMASR) to doubled CO₂ and its associated climate change using the Community Earth System Model in tandem with a dynamic global vegetation model. The results showed that the doubled CO₂ had a positive effect on the leaf area index of the SMASR, but its associated climate change exerted a negative effect in most parts of the SMASR. Although climate change had a weak negative effect on ground runoff at the regional scale, a 4.74 mm increase was predicted under the combined effect of doubled CO₂ and climate change, largely due to the positive effect of doubled CO₂. Spatially, the evident increase in ground runoff, which primarily occurred in the southeastern part of the SMASR, resulted from decreased ground evaporation and canopy transpiration under the doubled CO₂ condition. A negative effect was predicted in the central west as a result of increased temperature and a changed precipitation under doubled CO₂. These findings implied that the condition of water resources would be improved slightly under a doubled CO₂ condition, whereas there would be a larger spatial heterogeneity in relation to different sensitivities of vegetation and hydrological variables to doubled CO₂ and associated climate change.     

Geophysical monitoring technology for CO<Subscript>2</Subscript> sequestration

Geophysical techniques play key roles in the measuring, monitoring, and verifying the safety of CO₂ sequestration and in identifying the efficiency of CO₂-enhanced oil recovery. Although geophysical monitoring techniques for CO₂ sequestration have grown out of conventional oil and gas geophysical exploration techniques, it takes a long time to conduct geophysical monitoring, and there are many barriers and challenges. In this paper, with the initial objective of performing CO₂ sequestration, we studied the geophysical tasks associated with evaluating geological storage sites and monitoring CO₂ sequestration. Based on our review of the scope of geophysical monitoring techniques and our experience in domestic and international carbon capture and sequestration projects, we analyzed the inherent difficulties and our experiences in geophysical monitoring techniques, especially, with respect to 4D seismic acquisition, processing, and interpretation.     

Seasonal dynamics of soil CO<Subscript>2</Subscript> effluxes with responses to environmental factors in lower subtropical forests of China

Seasonal metrics and environmental responses to forestry soil surface CO₂ emission effluxes among three types of lower subtropical forests were consistently monitored over two years with static chamber-gas chromatograph techniques among three types of lower subtropical forests. Results showed that annual CO₂ effluxes (S+L) reached 3942.20, 3422.36 and 2163.02 CO₂ g·m^?2·a^?1, respectively in the monsoon evergreen broadleaf forest, mixed broadleaf-coniferous forest and coniferous forest. All the three types of forests revealed the same characteristics of seasonal changes with the CO₂ effluxes peaking throughout June to August. During this peaking period, the effluxes were 35.9%, 38.1% and 40.2% of the total annual effluxes, respectively. The CO₂ emission process responding to the environmental factors displayed significantly different patterns in forestry soils of the three types of forests. The coniferous forest (CF) was more sensitive to temperature than the other two types. The Q ₁₀ values were higher, along with greater seasonal variations of the CO₂ efflux, indicating that the structurally unique forestry ecosystem has disadvantage against interferences. All the three types of forestry CO₂ effluxes showed significant correlation with the soil temperature (T _s), soil water content (M _s) and air pressure (P _a). However, stepwise regression analysis indicated no significant correlation between air pressure and the soil CO₂ efflux. With an empirical model to measure soil temperature and water content in 5 cm beneath the soil surface, the CO₂ effluxes accounting for 75.7%, 77.8% and 86.5% of the efflux variability respectively in soils of BF, MF and PF were calculated. This model can be better used to evaluate the CO₂ emission of soils under water stress and arid or semi-arid conditions.     

Cluster analysis of soil CO<Subscript>2</Subscript> data from Mt. Etna (Italy) reveals volcanic influences on temporal and spatial patterns of degassing

Soil CO₂ concentration data were collected periodically from July 2001 to June 2005 from sampling site grids in two areas located on the lower flanks of Mt. Etna volcano (Paternò and Zafferana Etnea–Santa Venerina). Cluster analysis was performed on the acquired data in order to identify possible groups of sites where soil degassing could be fed by different sources. In both areas three clusters were recognised, whose average CO₂ concentration values throughout the whole study period remained significantly different from one another. The clusters with the lowest CO₂ concentrations showed time-averaged values ranging from 980 to 1,170 ppm vol, whereas those with intermediate CO₂ concentrations showed time-averaged values ranging from 1,400 to 2,320 ppm vol, and those with the highest concentrations showed time-averaged values between 1,960 and 55,430 ppm vol. We attribute the lowest CO₂ concentrations largely to a biogenic source of CO₂. Conversely, the highest CO₂ concentrations are attributed to a magmatic source, whereas the intermediate values are due to a variable mixing of the two sources described above. The spatial distribution of the CO₂ values related to the magmatic source define a clear direction of anomalous degassing in the Zafferana Etnea–Santa Venerina area, which we attribute to the presence of a hidden fault, whereas in the Paternò area no such oriented anomalies were observed, probably because of the lower permeability of local soil. Time-series analysis shows that most of the variations observed in the soil CO₂ data from both areas were related to changes in the volcanic activity of Mt. Etna. Seasonal influences were only observed in the time patterns of the clusters characterised by low CO₂ concentrations, and no significant interdependence was found between soil CO₂ concentrations and meteorological parameters. The largest observed temporal anomalies are interpreted as release of CO₂ from magma batches that migrated from deeper to shallower portions of Etna’s feeder system. The pattern of occurrence of such episodes of anomalous gas release during the observation period was quite different between the two studied areas. This pattern highlighted an evident change in the mechanism of magma transport and storage within the volcano’s feeder system after June 2003, interpreted as magma accumulation into a shallow (<8 km depth) reservoir.