First-order decay models to describe soil C-CO2 Loss after rotary tillage

To further understand the impact of tillage on CO2 emission, the applicability of two conceptual models was tested, which describe the CO2 emission after tillage as a function of the non-tilled emission plus a correction due to the tillage disturbance. Models assume that C in readily decomposable organic matter follows a first-order reaction kinetics equation as: dCsoil (t) / dt = -k Csoil (t), and that soil C-CO2 emission is proportional to the C decay rate in soil, where Csoil(t) is the available labile soil C (g m-2) at any time (t) and k is the decay constant (time-1). Two possible assumptions were tested to determine the tilled (F T) fluxes: the decay constants (k) of labile soil C before and after tillage are different (Model 1) or not (Model 2). Accordingly, C flux relationships between non-tilled (F NT) and tilled (F T) conditions are given by: F T = F NT + a1 e-a2t (model 1) and F T = a3 F NT e-a4t (model 2), where t is time after tillage. Predicted and observed CO2 fluxes presented good agreement based on the coefficient of determination (R² = 0.91). Model comparison revealed a slightly improved statistical fit of model 2, where all C pools are assigned with the same k constant. Rotary speed was related to increases in the amount of labile C available and to changes of the mean resident labile C pool available after tillage. This approach allows describing the temporal variability of tillage-induced emissions by a simple analytical function, including non-tilled emission plus an exponential term modulated by tillage and environmentally dependent parameters.Para entendimento do impacto do preparo do solo sobre as emissões de CO2 desenvolvemos e aplicamos dois modelos conceituais que são capazes de prever a emissão de CO2 do solo após seu preparo em função da emissão da parcela sem distúrbio, acrescida de uma correção devido ao preparo. Os modelos assumem que o carbono presente na matéria orgânica lábil segue uma cinética de decaimento de primeira ordem, dada pela seguinte equação: dCsoil (t) / dt = -k Csoil (t), e que a emissão de C-CO2 é proporcional a taxa de decaimento do C no solo, onde Csolo(t) é a quantidade de carbono lábil disponível no tempo (t) e k é a constante de decaimento (tempo-1). Duas suposições foram testadas para determinação das emissões após o preparo do solo (Fp): a constante de decaimento do carbono lábil do solo (k) antes e após o preparo é igual (Modelo 1) ou desigual (Modelo 2). Conseqüentemente, a relação entre os fluxos de C das parcelas sem distúrbio (F SD) e onde o preparo do solo foi conduzido (F P) são dadas por: F P = F SD + a1 e-a2t (modelo 1) e F P = a3 F SD e-a4t (modelo 2), onde t é o tempo após o preparo. Fluxos de CO2 previstos e observados relevam um bom ajuste dos resultados com coeficiente de determinação (R²) tão alto quanto 0,91. O modelo 2 produz um ajuste ligeiramente superior quando comparado com o outro modelo. A velocidade das pás da enxada rotativa foi relacionada a um aumento na quantidade de carbono lábil e nas modificações do tempo de residência médio do carbono lábil do solo após preparo. A vantagem desta metodologia é que a variabilidade temporal das emissões induzidas pelo preparo do solo pode ser descrita a partir de uma função analítica simples, que inclui a emissão da parcela sem distúrbio e um termo exponencial modulado por parâmetros dependentes do preparo e de condições ambientais onde o experimento foi conduzido