Inverse problems that are governed by partial differential equations arise in many applications in computational science and engineering. Solving these large scale problems is a real challenge to the existing numerical methods, as they are generally highly ill-posed and non-convex. These difficulties are usually handled by introducing statistical Bayesian estimation methods that promote a-priori k ...