Stochastic Models of Soil Denitrification

Abstract

Soil denitrification is a highly variable process that appears to be lognormally distributed. This variability is manifested by large sample coefficients of variation for replicate estimates of soil core denitrification rates. Deterministic models for soil denitrification have been proposed in the past, but none of these models predicts the approximate lognormality exhibited by natural denitrification rate estimates. In this study, probabilistic (stochastic) models were developed to understand how positively skewed distributions for field denitrification rate estimates result from the combined influences of variables known to affect denitrification. Three stochastic models were developed to describe the distribution of measured soil core denitrification rates. The driving variables used for all the models were denitrification enzyme activity and CO(2) production rates. The three models were distinguished by the functional relationships combining these driving variables. The functional relationships used were (i) a second-order model (model 1), (ii) a second-order model with a threshold (model 2), and (iii) a second-order saturation model (model 3). The parameters of the models were estimated by using 12 separate data sets (24 replicates per set), and their abilities to predict denitrification rate distributions were evaluated by using three additional independent data sets of 180 replicates each. Model 2 was the best because it produced distributions of denitrification rate which were not significantly different (P > 0.1) from distributions of measured denitrification rates. The generality of this model is unknown, but it accurately predicted the mean denitrification rates and accounted for the stochastic nature of this variable at the site studied. The approach used in this study may be applicable to other areas of ecological research in which accounting for the high spatial variability of microbiological processes is of interest.