“…It produces an image of M voxels of electrical resistivities (ρ i with i = 1,2, …, M ) given a set of N four electrode resistances ( R j [Ω] with j = 1,2, …, N ) by minimizing the objective function Ψ: where d is the data vector, m is the model given by the parameters of the inversion m j = log(ρ j ) (Ω m), f ( m ), is the forward model for parameters m , m 0 is a homogeneous starting model vector, W ε is a smoothing operator, λ is the regularization parameter that determines the amount of smoothing imposed on m , and W s is an error weighting matrix. We acknowledge the importance of choosing appropriate regularization parameters for accurate inversion results (Rao, Lesparre, Orozco, Wagner, & Javaux, 2020). After a series of synthetic model simulations, we concluded that 100 is an appropriate value for λ.…”