### Abstract

Markov chain Monte Carlo methods such as the Gibbs sampler and the Metropolis algorithm can be used to sample the solutions to non-linear inverse problems. In principle these methods allow incorporation of arbitrarily complex a priori information, but current methods allow only relatively simple priors to be used. We demonstrate how sequential simulation can be seen as an application of the Gibbs sampler, and how such a Gibbs sampler assisted by sequential simulation can be used to perform a random walk generating realizations of a relatively complex random function. We propose to combine this algorithm with the Metropolis algorithm to obtain an efficient method for sampling posterior probability densities for nonlinear inverse problems.

Original language | English |
---|---|

Title of host publication | IAMG Extended Abstracts |

Number of pages | 8 |

Publication date | 2010 |

Publication status | Published - 2010 |

Event | 14th Annual Conference of the International Association for Mathematical Geosciences - Budapest, Hungary Duration: 29 Aug 2010 → 2 Sep 2010 Conference number: 14 http://www.iamg2010-budapest.hu/ |

### Conference

Conference | 14th Annual Conference of the International Association for Mathematical Geosciences |
---|---|

Number | 14 |

Country | Hungary |

City | Budapest |

Period | 29/08/2010 → 02/09/2010 |

Internet address |

### Keywords

- Monte Carlo
- Prior
- Sequential simulation
- Inverse problem

## Cite this

Hansen, T. M., Mosegaard, K., & Cordua, K. S. (2010). Sampling informative/complex a priori probability distributions using Gibbs sampling assisted by sequential simulation. In

*IAMG Extended Abstracts*