Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data

Magnetotelluric (MT)data are inverted for smooth 2-D models using an extension of the existing 1-D algorithm, Occam's inversion. Since an MT data set consists of a finite number of imprecise data, an infinity of solutions to the inverse problem exists. Fitting field or synthetic electromagnetic data as closely as possible results in theoretical models with a maximum amount of roughness, or structure. However, by relaxing the misfit criterion only a small amount, models which are maximally smooth may be generated. The models are composed of a large number of rectangular prisms, each having a constant conductivity. A prior information, in the form of boundary locations only or both boundary locations and conductivity, may be included, providing a powerful tool for improving the resolving power of the data. Joint inversion of TE and TM synthetic data generated from known models allows comparison of smooth models with the true structure.