Temporal variability of soil CO2 emission after conventional and reduced tillage described by an exponential decay in time model

Variabilidade temporal da emissão de CO2 do solo após preparo convencional e reduzido descrito por modelo exponencial decrescente no tempo

Abstracts

To study Assessing the impact of tillage practices on soil carbon losses dependents it is necessary to describe the temporal variability of soil CO2 emission after tillage. It has been argued that large amounts of CO2 emitted after tillage may serve as an indicator for longer-term changes in soil carbon stocks. Here we present a two-step function model based on soil temperature and soil moisture including an exponential decay in time component that is efficient in fitting intermediate-term emission after disk plow followed by a leveling harrow (conventional), and chisel plow coupled with a roller for clod breaking (reduced) tillage. Emission after reduced tillage was described using a non-linear estimator with determination coefficient (R²) as high as 0.98. Results indicate that when emission after tillage is addressed it is important to consider an exponential decay in time in order to predict the impact of tillage in short-term emissions.

soil respiration; soil CO2 emission; soil tillage; greenhouse effect

A quantificação do impacto das práticas de preparo sobre as perdas de carbono do solo é dependente da habilidade de se descrever a variabilidade temporal da emissão de CO2 do solo após preparo. Tem sido sugerido que as grandes quantidades de CO2 emitido após o preparo do solo podem servir como um indicador das modificações nos estoques de carbono do solo em longo termo. Neste trabalho é apresentado um modelo de duas partes baseado na temperatura e na umidade do solo e que inclui um termo exponencial decrescente do tempo que é eficiente no ajuste das emissões intermediárias após preparo: arado de disco seguido de uma passagem com a grade niveladora (convencional) e escarificador de arrasto seguido da passagem com rolo destorroador (reduzido). As emissões após o preparo do solo são descritas utilizando-se estimativa não linear com um coeficiente de determinação (R²) tão alto quanto 0.98 após preparo reduzido. Os resultados indicam que nas previsões da emissão de CO2 após o preparo do solo é importante considerar um termo exponencial decrescente no tempo após preparo.

respiração do solo; emissão de CO2 do solo; preparo do solo; efeito estufa

SCIENTIFIC PAPERS

SOIL AND WATER ENGINEERING

Temporal variability of soil CO₂ emission after conventional and reduced tillage described by an exponential decay in time model

Variabilidade temporal da emissão de CO₂ do solo após preparo convencional e reduzido descrito por modelo exponencial decrescente no tempo

Luis G. Teixeira^I; Afonso Lopes^II; Newton La Scala Jr^III

^IAgrônomo, Pós Graduando Do PPG Em Produção Vegetal, FCAV/UNESP, [e-mail]

^IIProf. Adjunto, Departamento De Engenharia Rural, FCAV/UNESP, [e-mail]. Bolsista de Produtividade do CNPq

^IIIProf. Adjunto, Departamento De Ciências Exatas, FCAV/UNESP, [e-mail]. Bolsista de Produtividade do CNPq

ABSTRACT

To study Assessing the impact of tillage practices on soil carbon losses dependents it is necessary to describe the temporal variability of soil CO₂ emission after tillage. It has been argued that large amounts of CO₂ emitted after tillage may serve as an indicator for longer-term changes in soil carbon stocks. Here we present a two-step function model based on soil temperature and soil moisture including an exponential decay in time component that is efficient in fitting intermediate-term emission after disk plow followed by a leveling harrow (conventional), and chisel plow coupled with a roller for clod breaking (reduced) tillage. Emission after reduced tillage was described using a non-linear estimator with determination coefficient (R²) as high as 0.98. Results indicate that when emission after tillage is addressed it is important to consider an exponential decay in time in order to predict the impact of tillage in short-term emissions.

Keywords: soil respiration, soil CO₂ emission, soil tillage, greenhouse effect.

RESUMO

A quantificação do impacto das práticas de preparo sobre as perdas de carbono do solo é dependente da habilidade de se descrever a variabilidade temporal da emissão de CO₂ do solo após preparo. Tem sido sugerido que as grandes quantidades de CO₂ emitido após o preparo do solo podem servir como um indicador das modificações nos estoques de carbono do solo em longo termo. Neste trabalho é apresentado um modelo de duas partes baseado na temperatura e na umidade do solo e que inclui um termo exponencial decrescente do tempo que é eficiente no ajuste das emissões intermediárias após preparo: arado de disco seguido de uma passagem com a grade niveladora (convencional) e escarificador de arrasto seguido da passagem com rolo destorroador (reduzido). As emissões após o preparo do solo são descritas utilizando-se estimativa não linear com um coeficiente de determinação (R²) tão alto quanto 0.98 após preparo reduzido. Os resultados indicam que nas previsões da emissão de CO₂ após o preparo do solo é importante considerar um termo exponencial decrescente no tempo após preparo.

Palavras-chave: respiração do solo, emissão de CO₂ do solo, preparo do solo, efeito estufa.

INTRODUCTION

There has been a great deal of interest in studying soil CO₂ emission of agricultural areas, especially when the use and management of soils are considered. Due to the high spatial and temporal variability of soil CO₂ emission the uncertainty on this topic is significant; however, it is already accepted that reductions in tillage intensity and frequency could help achieving soil carbon recovery (BRONICK & LAL, 2005). Tillage causes higher rates of soil CO₂ emission in short and intermediate periods persisting up to some weeks after tillage and in longer term it has been observed soil carbon restore after conservation practices are applied (LAL, 2004). As a consequence, farmers have recently moved towards conservation tillage, and the number and frequency of conventional tillage events have declined over the last years (REICOSKY & ALLMARAS, 2003).

When tillage is applied, the protected organic matter becomes unprotected, in practice increasing the labile organic matter available to microbial activity (SIX et al., 1999). In addition, changes in soil carbon decay rate are expected with decreases in soil bulk density and increases in soil pore volume. Recently, based on these assumptions, conceptual models have been proposed to describe the tillage induced emission. LA SCALA et al. (2008) used no-till emission as a reference in order to model emission after tillage. REICOSKY & ARCHER (2007) suggested that the process is described by an initially rapid decline in CO₂ flux followed by a slower flux, requiring the use of a two-step function.A new two-step model is presented in order to fit emissions after tillage. The main hypothesis is that soil CO₂ emission after tillage is modulated by two components: the first directly related to soil temperature and soil moisture, and the second related to an expected rapid dynamics of the additional labile carbon induced by tillage. Thus, it is suggested the following function to describe the temporal variability of soil CO₂ emissions after tillage (FCO₂):

The first term of equation 1 ( F_I) incorporates soil temperature and soil moisture variability. The second term (FII) is related to the rapid changes due to the additional labile carbon (Clabile) made available to microbes by tillage. When this term is taken into account, it is considered that the labile carbon dynamics is governed by first-order decay equation:

In the equation above k (day^-1) is the so-called decay constant and CT is the labile carbon made available to microbial activity (in g m^-2) by tillage. The increase of labile carbon could be due to the crop residuals incorporated into the soil by tillage or the break of aggregates and exposition of protected fresh organic carbon to microbial activity. The labile carbon was modeled by the following exponential decay equation:

Where C0T is the labile carbon available at the initial time and induced by tillage (g CO₂ m^-2), k is the decay rate, and t is the time (hours). One reasonable assumption is that soil CO₂ losses is proportional to the soil carbon decayed after tillage. Mathematically, this can be written as:

Considering the efficiency of microbial activity to convert soil carbon in CO₂ as µ, equation 2 is converted to:

Making a₃=βC_0Tk e a₄= k, leads to:

Equation (4) is similar to the one suggested by ELLERT & JANZEN (1999) that quantified the tillage-induced emission equalizing tilled minus untilled plot fluxes to an exponential decay in time function. Most Recently, LA SCALA et al. (2009a, b) have also used this approach in order to model emission after application of several tillage systems and adjustments. Their main hypothesis was that soil contains a carbon pool which is susceptible to release by tillage according to a first order kinetics model.

The dependence of soil CO₂ emission on soil temperature and soil moisture has been studied by several authors. Linear, exponential, and logarithmic models have been proposed in vegetated and bare soils (JABRO et al., 2008). In this work a linear dependence of emission on soil moisture and soil temperature is considered in a semi-empirical model. Therefore, the first term of equation 1 can be written as:

And the two-term proposed model written as:

Two tillage systems are mainly used in Southern Brazil: disk plow followed by a leveling harrow (conventional) or chisel plow coupled with a roller for clod breaking (reduced). The objective of this study is to apply equations 4 and 5 separately to describe emission after tillage, and the two-step model (equation 6) to assess soil CO₂ emission in a two-week time frame.

MATERIAL AND METHODS

The experiment was conducted in the experimental station of the FCAV/UNESP, located at 21º14´S and 48º17´W, in the State of São Paulo, Brazil. The climate of the region is classified as Cwa, according to Köppen, having rainy summers and dry winters. Trials started on December 15, 2005, in two adjacent plots measuring 10 x 2 meters each. On a bare Latosol (pH around 5) previously cultivated with maize (Zea mays, L.) two tillage systems were applied: disk plow followed by a leveling harrow (conventional), and chisel plow coupled with a roller for clod breaking (reduced). The tillage depth in both plots was 30 cm, and the roller for clod breaking worked at a 6 cm depth. Soil was kept bare during the two week study as no seedling occurred. For those operations, it was used a Valtra tractor (140 cv) running at 4 km h^-1 and coupled to a harrow or a chisel.

Just after tillage, six PVC soil collars having 10 cm diameter each were installed in each plot, inserted 3 cm inside the soil. Measurements of soil CO₂ emission, soil temperature and soil moisture started few minutes after tillage during a period of nine days, from the 15^th to the 17^th , from the 19^th to the 22^nd , and at the 28^th and the 29^th of December 2005, always in the mornings. Emissions were measured using a LI-6400 system (Li-Cor, NE, USA) coupled with a soil temperature probe (LI-8100-102 model, thermistor based) that was inserted into the soil in all studied point at a 0-20 cm depth. This system uses infrared sensors to compute the time changes of CO₂ concentration inside a chamber, and calculates emission in each of the studied points. One single measurement takes around 1.5 minutes, based on raw measurements of CO₂ concentration changes inside the chamber every 3 seconds. Soil moisture (% volume) was monitored (0-20 cm depth) in each studied point by using TDR probes (Campbell Scientific, Utah, USA). The sofware Statistica (StatSoft, Inc.) was used for fitting soil CO₂ emission after conventional and reduced tillage using equations 4 to 6. Model parameters were estimated using non-linear least square estimation with the Gauss-Newton algorithm.

RESULTS AND DISCUSSION

Results of modeling soil CO₂ emission after conventional and reduced tillage by using FI and FII functions (equations 5 and 4) separately can be found in Table 1. Function FI is a simple linear dependence on soil temperature and soil moisture and explains 9 and 2% of emission variability in the reduced and conventional treatments. When FII only is applied to data, around 30% of emission variability is explained with just one non-significant (p>0.10) parameter (a₄). It can be noticed that a4 values were positive, indicating a predominant exponential decay trend after tillage.

When the two-step function is applied to the data set from both treatments, fitting parameters are significantly improved (Table 2). The coefficients of determination (R²) for the fittings were 0.93 e 0.98, for conventional and reduced treatments, respectively. Therefore, the exponential decay in time after tillage coupled with a model for linear dependence on soil temperature and soil moisture is a more adequate model to describe, in this case, emission after conventional and reduced tillage. Extracted parameters (a1-a4) of model proposed by Equation 6 were significant (p<0.10) with an exception made to a₂ parameter for conventional treatment. The a1 and a2 parameters were mostly positive as expected for the dependence of soil CO₂ emission on soil temperature and soil moisture (LIU et al., 2008; SAVAGE & DAVIDSON, 2003; SALOKHE & RAMALINGAM, 2001). In addition, values of a₄ were positive, therefore, a typical exponential decay is adequate for describing soil emission after tillage as has been reported in many studies conducted at different soil and climate conditions (LA SCALA et al., 2001, 2005, 2006; REICOSKY, 2002; ROCHETTE & ANGERS, 1999; FORTIN et al., 1996).

Despite all the efforts, only a few authors have proposed models that have considered soil moisture and temperature in modeling CO₂ soil emission after tillage using an exponential decay in time model. The coefficients of determination found here are superior to the ones found by ELLERT & JANZEN (1999) when simple decay in time function was applied to fitting observed data. Also, it is noticeable that total emissions predicted (by integrating the area bellow the curves of Figure 1) show good agreement with observed data. Deviations between predicted and observed values were around 2% in both treatments: 260.65/254.70 gCO₂ m^-2 and 107.19/105.04 (observed/predicted) gCO₂ m^-2 after conventional and reduced tillage, respectively (Table 2).

Figure 1 presents the observed (solid lines) and predicted (dashed lines) soil CO₂ emission in conventional (triangles) and reduced (squares) tillage treatments during the studied period. In the m^-2 h-¹ in first studied day, some hours after tillage, emissions were as high as 1.74 and 0.78 g CO₂ mthe conventional and reduced treatments, respectively. It is possible to observe an exponential decay-like aspect of emissions that, in our model, is modeled by the exponential in time function. Small fluctuations throughout experiment were presumably due to the changes in soil temperature and soil moisture. Rainfall occurred in this experiment in the days 1, 3 to 5, 10 and 12, with daily precipitations as high as 10 mm. Changes in soil moisture (% volume) and soil temperature (ºC) can be observed in Figures 2a and 2b, respectively. Soil temperature was as low as 22 ºC in the 5^th day and as high as 31 ºC in the 13^th day after tillage (Figure 2a). Changes in soil moisture were also large in time, as in reduced treatments they were as low as 14%, just after tillage, and as high as 29% (Figure 2b). Soil moisture was kept higher in conventional then in the reduced treatments during almost the whole studied period, reaching highest values in the 4^th day after tillage (31.2%). By comparing predicted (dashed lines) to observed values (solid lines) in Figure 1, it is possible to notice that our model was able to capture minor fluctuations of emission in both plots during the experiment.

Table 3 and Figure 3a and 3b present parameters and regression adjustments (p<0.001) between predicted and observed soil CO₂ emission after tillage. As one can notice, both adjustments are close to the 1:1 straight line (bisector). The linear and angular coefficients have values near to 0 and 1, respectively (Figures 3a, and 3b). In a longer term study, YOO et al. (2006) reported that soil temperature and soil water contents, which are the most important factors influencing soil organic carbon mineralization, do not successfully explain tillage impacts on mean soil organic carbon mineralization rates after conventional tillage. However, the model based on Equation 6, taking into account soil temperature and soil moisture variability, added to an exponential decay in time function, was able to explain satisfactorily the intermediate-term soil CO₂ emission after tillage (conventional and reduced). It can be sustained that time after tillage is an important variable defining emission after tillage events. This is probably associated to the additional labile carbon induced after tillage event and how it decays over time.

CONCLUSION

A model including an exponential decay in time part was able to better describe short-term emissions after tillage. The model that takes into account temperature, moisture coupled with an exponential in time decay function is better suitable to describe CO₂ emission over time than a simple model that considers only soil temperature and moisture.

Additional efforts should be placed to understand the complexity of emission after tillage, but it is sustained that using an exponential decay in time after tillage is an important aspect that should be taking into account.

Recebido pelo Conselho Editorial em: 19-8-2008

Aprovado pelo Conselho Editorial em: 30-12-2009

• BRONICK, C.J.; LAL, R. Soil structure and management: a review. Geoderma, Amsterdam, v.124, p.3-22, 2005.
• ELLERT, B.H.; JANZEN, H.H. Short-term influence of tillage on CO₂ fluxes from a semi-arid soil on the Canadian prairies. Soil & Tillage Research, Amsterdam, v.50, p.21-32, 1999.
• FORTIN, M.C.; ROCHETTE, P.; PATTEY, E. Soil carbon dioxide fluxes from conventional and no-tillage small-grain cropping systems. Soil Science Society of America Journal, Madison, v.60, p.1.541-1.547, 1996.
• JABRO, J.D.; SAINJU, U.; STEVENS, W.B.; EVANS, R.G. Carbon dioxide flux as affected by tillage and irrigation in soil converted from perennial forages to annual crops. Journal of Environmental Management, Amsterdam, v.88, n.4, p.1.478-1.484, 2008.
• LAL, R. Soil carbon sequestration to mitigate climate change. Geoderma, Amsterdam, v. 123, p.122, 2004.
• LA SCALA, N.; BOLONHEZI, D.; PEREIRA, G.T. Short-term soil CO₂ emission after conventional and reduced tillage of a no-till sugar cane area in southern Brazil. Soil & Tillage Research, Amsterdam, v.91, p.244-248, 2006.
• LA SCALA, N.; LOPES, A.; MARQUES JR., J.; PEREIRA, G.T. Carbon dioxide emissions after application of tillage systems for a dark red latosol in souther Brazil. Soil & Tillage Research, Amsterdam, v.62, p.163-166, 2001.
• LA SCALA, N.; LOPES, A.; PANOSSO, A.R.; CÂMARA, F.T.; PEREITA, G.T. Soil CO₂ efflux following rotary tillage of a tropical soil. Soil & Tillage Research, Amsterdam, v.84, p.222-225, 2005.
• LA SCALA, N.; LOPES, A.; SPOKAS, K.; ARCHER, D.; REICOSKY, D.C. First-order decay models to describe soil C-CO₂ loss after rotary tiller. Scientia Agricola, Piracicaba, v.66, p.650-657, 2009b.
• LA SCALA, N.; LOPES, A.; SPOKAS, K.; ARCHER, D.; REICOSKY, D.C. Short-term temporal changes of bare soil CO₂ fluxes after tillage described by first-order decay models. European Journal of Soil Science, Edinburgh, v.60, p.258-264, 2009a.
• LA SCALA, N.; LOPES, A.; SPOKAS, K.; BOLONHEZI, D.; ARCHER, D.; REICOSKY, D.C. Short-term temporal changes of soil carbon losses after tillage described by a first-order decay model. Soil & Tillage Research, Amsterdam, v.99, p.108-118, 2008.
• LIU, H.; ZHAO, P.; LU, P.; WANG, Y.S.; LIN, Y.B.; RAO, X.Q. Greenhouse gas fluxes from soils of different land-use types in a hilly area of south China. Agriculture, Ecosystems and Environment, Amsterdam, v.124, n.1-2, p.125-135, 2008.
• REICOSKY, D.C. Long-term effect of moldboard plowing on tillage-induced CO₂ loss. In: KIMBLE, J.M.; LAL, R. Agricultural practices and policies for carbon sequestration in soil Boca Raton: CRC Press, 2002. p.87-97.
• REICOSKY, D.C.; ALLMARAS, R.R. Advances in tillage research in North American cropping systems. In: SHRESTHA, A. Cropping systems: trends and advances. New York: Haworth Press, 2003. p.75-125.
• REICOSKY, D.C.; ARCHER, D.W. Moldboard plow tillage depth and short-term carbon dioxide release. Soil & Tillage Research, Amsterdam, v.94, p.109-121, 2007.
• ROCHETTE, P.; ANGERS, D.A. Soil surface carbon dioxide fluxes induced by spring, summer and fall moldboard plowing in a sandy loam. Soil Science Society of America Journal, Madison, v.63, p.621-628, 1999.
• SALOKHE, V.M.; RAMALINGAM, N. Effects of direction of rotation of a rotary tiller on properties of Bangkok clay soil. Soil Tillage Research, Amsterdam, v.63, p.65-74, 2001.
• SAVAGE, K.E.; DAVIDSON, E.A. A comparison of manual and automaged systems for soil CO₂ flux measurements: trade-offs between spatial and temporal resolution. Journal of Experimental Botany, Oxford, v.54, p.891-899, 2003.
• SIX, J.; ELLIOTT, E.T.; PAUSTIAN, K. Aggregate and soil organic matter dynamics under conventional and no-tillage systems. Soil Science Society of America Journal, Madison, v.63, p.1.350-1.358, 1999.
• YOO G.; NISSEN, T.M.; WANDER, M.M. Use of physical properties to predict the effects of tillage practices on organic matter dynamics in three Illinois soils. Journal of Environmental Quality, Madison, v.35, p.1.576-1.583, 2006.

Publication Dates

• Publication in this collection
  26 May 2010
• Date of issue
  Apr 2010

History

• Accepted
  30 Dec 2009
• Received
  19 Aug 2008