Inverse method for estimating respiration rates from decay time series
Name
Forney-2012-Inverse method for e.pdf
Size
601.77 KB
Format
Adobe PDF
Checksum (MD5)
5944f444277e547c64628e4ec6c66b13
Author(s)
Forney, David C.
Rothman, Daniel H.
Date Issued
September 2012
Journal
Biogeosciences
Publisher
Copernicus GmbH
Citation
Forney, D. C., and D. H. Rothman. “Inverse Method for Estimating Respiration Rates from Decay Time Series.” Biogeosciences 9.9 (2012): 3601–3612.
Version
Final published version
Abstract
Long-term organic matter decomposition experiments typically measure the mass lost from decaying organic matter as a function of time. These experiments can provide information about the dynamics of carbon dioxide input to the atmosphere and controls on natural respiration processes. Decay slows down with time, suggesting that organic matter is composed of components (pools) with varied lability. Yet it is unclear how the appropriate rates, sizes, and number of pools vary with organic matter type, climate, and ecosystem. To better understand these relations, it is necessary to properly extract the decay rates from decomposition data. Here we present a regularized inverse method to identify an optimally-fitting distribution of decay rates associated with a decay time series. We motivate our study by first evaluating a standard, direct inversion of the data. The direct inversion identifies a discrete distribution of decay rates, where mass is concentrated in just a small number of discrete pools. It is consistent with identifying the best fitting "multi-pool" model, without prior assumption of the number of pools. However we find these multi-pool solutions are not robust to noise and are over-parametrized. We therefore introduce a method of regularized inversion, which identifies the solution which best fits the data but not the noise. This method shows that the data are described by a continuous distribution of rates, which we find is well approximated by a lognormal distribution, and consistent with the idea that decomposition results from a continuum of processes at different rates. The ubiquity of the lognormal distribution suggest that decay may be simply described by just two parameters: a mean and a variance of log rates. We conclude by describing a procedure that estimates these two lognormal parameters from decay data. Matlab codes for all numerical methods and procedures are provided.
MIT Department
Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences
Massachusetts Institute of Technology. Department of Mechanical Engineering
Terms of Use
Creative Commons Attribution 3.0
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.5194/bg-9-3601-2012
