Models of subsurface structures are important for successful deposit exploration, but are challenged by the need to integrate data from different geophysical methods. In the present study, we evaluated a method of joint inversion in which audio-magneto telluric (AMT) and dual frequency induced polarization (DFIP) data sets are inverted simultaneously to produce a consistent 2D resistivity model to show a clear image of subsurface structures. To achieve the objectives, we conducted AMT and DFIP surveys along the same survey line within the Dongjun lead–zinc deposit in inner Mongolia by measuring 31 AMT survey sites with a station spacing of 40 m on a 1440 m survey track and operated in fifty-three frequencies in the range of 1–10,400 Hz to record the resistivity distribution of subsurface to depths exceeding 800 m. The same survey setup up was applied to the DFIP method using a pole–dipole array configuration and operating frequencies of 4 Hz and 4/13 Hz. The two-dimensional (2D) model obtained from AMT data revealed distinct low-resistivity anomalies in the middle of the 2D inversion model. In contrast, the DFIP inversion model showed a high resistive body in the same region with relatively high percent frequency effect (PFE) indicating high chargeability. In response to the discrepancies observed in the separate 2D inversion models, we implemented a joint inversion for both the AMT and DFIP data sets. The joint inversion resistivity model shows surficial conducting bodies and a high conductive body along the profile with relatively high PFE, indicating high chargeability. The final joint inversion resistivity model clearly images the large silica alteration zone and the Pb-Zn mineralization. This study demonstrates the feasibility of a joint inversion methodology and highlights the value of integrating geophysical methods through joint inversion for enhanced characterization and exploration of lead–zinc ores.