Author:
Glagolev M. V.,Suvorov G. G.,Il’yasov D. V.,Sabrekov A. F.,Terentieva I. E.
Abstract
The spread of published values of the rate of methane uptake by soils makes up several orders of magnitude from 0.0001 to 1 mgm-2h-1, which is comparable in magnitude to the spread of estimates of the release of CH4 out of waterlogged soils. The high values of CH4 emissions out of waterlogged soils are well explained, since with high methane production, it can be removed from the soil at almost any speed through a convective (most often bubble) transport mechanism. But when being absorbed by the soil, methane can penetrate in it only due to an apparently slow diffusion mechanism. Thus, the question arises of the maximum theoretically justified assessment of methane consumption by the soil. The aim of our work was to try to quantify the maximum possible amount of CH4 consumption by the soil relying on a strict basis of soil biokinetics and physics.
To estimate the maximum specific absorption flux of CH4 by the soil, we used the "mass conservation equation" [Walter et al., 1996; Zhuang et al., 2004; Глаголев, 2006, p. 316; 2010, p. 35-36]:
C/t = -F/z + Qebull + Qplant + Rprod + Roxid,
where C (mg/m3) is the concentration of methane at time t at depth z; F (mgm-2h-1) is the specific flux of methane due to diffusion; Qebull and Qplant (mgm3h-1) are the rates of change in methane concentration at time t at depth z due to the formation of bubbles and drainage through the roots of plants, respectively; Rprod and Roxid (mgm-3 h-1) are the rates of formation and consumption of methane, respectively.
Since we going to estimate the flux of CH4 only at its maximum possible consumption, the equation is simplified, as far as its terms accounted for the formation and transport of methane (Rprod, Qebull, Qplant) will be equal to 0. Finally, we will consider the system in a steady state, i.e. C/t=0. Thus:F(t,z)/z = Roxid(t,z).
Using Fick's first law to calculate the diffusion flux (used with a modified sign compared to its traditional form):
F(t,z) = D(z)C/z,
where D(z) is the diffusion coefficient [Zhuangetal.,2004]; and the modified Michaelis-Menten equation for calculating methane oxidation is:Roxid(t,z) = -Vmax(C-CTh)/(KM + C-CTh), where CTh (mgm-3) is the threshold concentration [Panikov, 1995, p. 151]; Vmax (mgm-3h-1) is the maximum specific consumption rate; KM (mgm-3) is the halfsaturation constant, and also under assumptions, (i) the concentration of CH4 is approximately equal to atmospheric (CA=1.29mg/m3) at the upper boundary (soil/atmosphere); (ii) the flux of CH4 can be assumed to be zero at an infinitely great depth [Born et al., 1990]; (iii) D, Vmax and KM (C- CTh) do not change with depth. Therefore, the absolute value of the specific flux from the atmosphere to the soil is:
|F(0)|= (CA-CTh)(VmaxD/KM).
The maximum value of the diffusion coefficient can be estimated by the Penman equation: D=D oPa0.66, where Do is the diffusion coefficient in air; Pa is the porosity of aeration [Смагин,2005,p.165]. Since we are going to estimate the maximum value of diffusion, we will take the limit value of porosity, which is 1, but as far as the proportion of pores of stable aeration accounts for half of the total pore volume [Растворова,1983,p.52], then for further calculations we will take Pa = 0.5, hence D = D o0.33. According to [Arah andStephen, 1998], for CH4
Do = 1.910-5∙(T/273)1.82 m2/s = 6.810-2∙(T/273)1.82 m2/h,
where T is temperature (K). When solving our diffusion problem, we assumed that the temperature is the same throughout the soil profile, and is 293 K. then D=6.810-2∙(293/273)1.820.33 = 2.5510-2m2/h.
The maximum rate of CH4 oxidation by soil was experimentally estimated in [Bender and Conrad, 1992] and was 57.3 mg/(hm3), which is in good agreement with the value of Vmax=47mg/(hm3) obtained at T=32 C according to the temperature dependence for automorphic soils of boreal forests Vmax=1.5(T ‑5.4)/10mmol/(hL), given in the work of Zhuang et al. [2004].
The halfsaturation constant is the concentration of the substrate, at which the specific growth rate of microorganisms takes a value equal to a half of the maximum. Summaries of the values KM have been repeatedly published (see, for example, [King,1992,Tab.II; Segers,1998,Tab.4; Глаголев,2006,pp.324-325]). For our purposes, we should take the KM obtained directly in the experiments with substrate concentrations (CH4) closest to those found in natural conditions. The minimum value (310-8mol/L) is given in [Bender and Conrad, 1992]. This value corresponds to the methane concentration in the air of about 20 ppm (14.3 mg/m3). This КМ value will be taken for further calculations.
The threshold concentration of CH4 for methanotrophs in the upper soil layer, given in the scientific literature, varies from 0.1 to 3.5 ppm [Crill, 1991; Bender and Conrad, 1992; Kravchenko et al., 2010]. Since we are interested in the minimum value of this indicator, we will bring it to the minimum temperature (273 K or 0 C): CTh=0.0714 mg/m3.
Now, having all the necessary numerical values, we can estimate the maximum intensity of methane consumption by natural soils:
|F(0)|= 1.2186(57.32.5510-2/14.3)0.39mg/(m2h).
Thus, for a certain "ideal" soil (evenly warmed throughout the profile, perfectly aerated, and at the same time containing enough moisture to create optimal living conditions for methanotrophs, which, by the way, are extremely numerous in the soil, and their methane halfsaturation constant is very low, etc.) we obtained an absorption intensity of CH4 of about 0.39 mg/(m2h). Since the combination of optimal values of all factors affecting methane consumption is very unlikely (or, rather, even improbable) in real soils, the resulting value can be considered extremely possible. And in view of this, the empirical generalization made in [Crill, 1991] becomes clear: "From the Amazon floodplain to the Arctic, the most rapid rates rarely exceed 6 mgCH4m-2d-1" i.e. 0.25 mg/(m2h).
Conclusion. So, we considered the absorption of methane as a biochemical process (following the Michaelis-Menten law with certain kinetic parameters), limited by diffusion in porous medium (soil). Based on this theoretical analysis, we came to the conclusion that the extremely large values of the specific absorption flux of CH4 (about 0.4 mgm-2h-1 and more), which are sometimes found in the literature, are unrealistic, if we are talking about the soils, which are always under methane concentrations no greater than atmospheric 1.8 ppmv. This applies to the vast majority of soils almost all, except for wetlands and soils covering landfills, underground gas storage facilities or other powerful sources of methane.
Reference63 articles.
1. Бронштейн И.Н., Семендяев К.А. 1986. Справочник по математике для инженеров и учащихся втузов. М.: Наука. 544 с. [Bronshtein I.N., Semendyaev K.A. 1986. Spravochnik po matematike dlya inzhenerov i uchashchikhsya vtuzov. M.: Nauka. 544 p. (In Russian)]
2. Вавилин В.А., Васильев В.Б. 1979. Математическое моделирование процессов биологической очистки сточных вод активным илом. М.: Наука. 119 с. [Vavilin V.A., Vasil'ev V.B. 1979. Matematicheskoe modelirovanie protsessov biologicheskoi ochistki stochnykh vod aktivnym ilom. M.: Nauka. 119 p. (In Russian)]
3. Вадюнина А.Ф., Корчагина З.А. 1973. Методы исследования физических свойств почв и грунтов. М.: Высшая школа. 399 с. [Vadyunina A.F., Korchagina Z.A. 1973. Metody issledovaniya fizicheskikh svoistv pochv i gruntov. M.: Vysshaya shkola. 399 p. (In Russian)]
4. Варфоломеев С.Д., Калюжный С.В. 1990. Биотехнология: Кинетические основы микробиологических процессов. М.: Высш. шк. 296 с. [Varfolomeev S.D., Kalyuzhnyi S.V. 1990. Biotekhnologiya: Kineticheskie osnovy mikrobiologicheskikh protsessov. M.: Vyssh. shk. 296 p. (In Russian)]
5. Воронин А.Д. 1986. Основы физики почв. М.: Изд-во МГУ. 244 с. [Voronin A.D. 1986. Osnovy fiziki pochv. M.: Publishing house MSU. 244 p. (In Russian)]