<div>The basic linear inverse problem of atmospheric release can be formulated as <strong>y</strong> = M <strong>x</strong> + <strong>e</strong> , where <strong>y</strong> is the measurement vector which is typically in the form of gamma dose rates or concentrations, M is the source-receptor-sensitivity (SRS) matrix, <strong>x</strong> is the unknown source term to be estimated, and <strong>e</strong> is the model residue. The SRS matrix M is computed using an atmospheric transport model coupled with meteorological reanalyses. The inverse problem is typically ill-conditioned due to number of uncertainties, hence, the estimation of the source term is not straightforward and additional information, e.g. in the form of regularization or the prior source term, is often needed. Besides, traditional techniques rely on assumption that the SRS matrix is correct which is not realistic due to the number of approximations made during its computation. Therefore, we propose relaxation of the inverse model using introduction of the term &#916;<sub>M</sub> such as <strong>y</strong> = ( M+ &#916;<sub>M</sub> ) <strong>x</strong> + <strong>e</strong> leading to non-linear inverse problem formulation, where &#916;<sub>M</sub> can be, as an example, parametric perturbation of the SRS matrix M in the spatial or temporal domain. We estimate parameters of this perturbation together with solving the inverse problem using variational Bayes procedure. The method will be validated on synthetic dataset as well as demonstrated on real case scenario such as the controlled tracer experiment ETEX or episode of ruthenium-106 release over the Europe on the fall of 2017.</div>