On the capability of Monte Carlo and adjoint inversion techniques to derive posterior parameter uncertainties in terrestrial ecosystem models
(2012) In Global Biogeochemical Cycles 26. p.3025-3025- Abstract
- Terrestrial ecosystem models (TEMs) contain the coupling of many biogeochemical processes with a large number of parameters involved. In many cases those parameters are highly uncertain. In order to reduce those uncertainties, parameter estimation methods can be applied, which allow the model to be constrained against observations. We compare the performance and results of two such parameter estimation techniques - the Metropolis algorithm (MA) which is a Markov Chain Monte Carlo (MCMC) method and the adjoint approach as it is used in the Carbon Cycle Data Assimilation System (CCDAS). Both techniques are applied here to derive the posterior probability density function (PDF) for 19 parameters of the Biosphere Energy Transfer and Hydrology... (More)
- Terrestrial ecosystem models (TEMs) contain the coupling of many biogeochemical processes with a large number of parameters involved. In many cases those parameters are highly uncertain. In order to reduce those uncertainties, parameter estimation methods can be applied, which allow the model to be constrained against observations. We compare the performance and results of two such parameter estimation techniques - the Metropolis algorithm (MA) which is a Markov Chain Monte Carlo (MCMC) method and the adjoint approach as it is used in the Carbon Cycle Data Assimilation System (CCDAS). Both techniques are applied here to derive the posterior probability density function (PDF) for 19 parameters of the Biosphere Energy Transfer and Hydrology (BETHY) scheme. We also use the MA to sample the posterior parameter distribution from the adjoint inversion. This allows us to assess if the commonly made assumption in variational data assimilation, that everything is normally distributed, holds. The comparison of the posterior parameter PDF derived by both methods shows that in most cases an approximation of the PDF by a normal distribution as used by the adjoint approach is a valid assumption. The results also indicate that the global minimum has been identified by both methods for the given set up. However, the adjoint approach outperforms the MA by several orders of magnitude in terms of computational time. Both methods show good agreement in the PDF of estimated net carbon fluxes for the decades of the 1980s and 1990s. (Less)
Please use this url to cite or link to this publication:
https://lup.lub.lu.se/record/4637510
- author
- Ziehn, T. ; Scholze, Marko LU and Knorr, W. LU
- organization
- publishing date
- 2012
- type
- Contribution to journal
- publication status
- published
- subject
- in
- Global Biogeochemical Cycles
- volume
- 26
- pages
- 3025 - 3025
- publisher
- American Geophysical Union (AGU)
- external identifiers
-
- wos:000308893500001
- ISSN
- 0886-6236
- DOI
- 10.1029/2011GB004185
- language
- English
- LU publication?
- yes
- id
- 7e1e2309-0577-4dff-aaaf-3b71ef5ae8d6 (old id 4637510)
- date added to LUP
- 2016-04-01 13:24:28
- date last changed
- 2024-02-14 13:33:07
@article{7e1e2309-0577-4dff-aaaf-3b71ef5ae8d6, abstract = {{Terrestrial ecosystem models (TEMs) contain the coupling of many biogeochemical processes with a large number of parameters involved. In many cases those parameters are highly uncertain. In order to reduce those uncertainties, parameter estimation methods can be applied, which allow the model to be constrained against observations. We compare the performance and results of two such parameter estimation techniques - the Metropolis algorithm (MA) which is a Markov Chain Monte Carlo (MCMC) method and the adjoint approach as it is used in the Carbon Cycle Data Assimilation System (CCDAS). Both techniques are applied here to derive the posterior probability density function (PDF) for 19 parameters of the Biosphere Energy Transfer and Hydrology (BETHY) scheme. We also use the MA to sample the posterior parameter distribution from the adjoint inversion. This allows us to assess if the commonly made assumption in variational data assimilation, that everything is normally distributed, holds. The comparison of the posterior parameter PDF derived by both methods shows that in most cases an approximation of the PDF by a normal distribution as used by the adjoint approach is a valid assumption. The results also indicate that the global minimum has been identified by both methods for the given set up. However, the adjoint approach outperforms the MA by several orders of magnitude in terms of computational time. Both methods show good agreement in the PDF of estimated net carbon fluxes for the decades of the 1980s and 1990s.}}, author = {{Ziehn, T. and Scholze, Marko and Knorr, W.}}, issn = {{0886-6236}}, language = {{eng}}, pages = {{3025--3025}}, publisher = {{American Geophysical Union (AGU)}}, series = {{Global Biogeochemical Cycles}}, title = {{On the capability of Monte Carlo and adjoint inversion techniques to derive posterior parameter uncertainties in terrestrial ecosystem models}}, url = {{http://dx.doi.org/10.1029/2011GB004185}}, doi = {{10.1029/2011GB004185}}, volume = {{26}}, year = {{2012}}, }