Geophysical inversion by iterative modeling involves fitting observations by adjusting model parameters. Both seismic and potential‐field model responses can be influenced by the adjustment of the parameters of the rock properties. The objective of this “cooperative inversion” is to obtain a model which is consistent with all available surface and borehole geophysical data. Although inversion of geophysical data is generally non‐unique and ambiguous, we can lessen the ambiguities by inverting all available surface and borehole data. This paper illustrates this concept with a case history in which surface seismic data, sonic logs, surface gravity data, and borehole gravity meter (BHGM) data are adequately modeled by using least‐squares inversion and a series of forward modeling steps.