Direct current (DC) resistivity has been widely investigated for non‐invasive mapping of dense non‐aqueous phase liquids (DNAPLs); however, due to its difficulty in distinguishing DNAPLs from adjacent soils, the DC method is limited for static detection of DNAPLs and more often employed for monitoring DNAPL changes over time. Time‐domain induced polarization (IP) can provide complementary information to better discriminate between DNAPL, water and surrounding soils. Since highly resistive DNAPLs tend to laterally spread and pool on polarizable (chargeable) clay lenses, combined DC and IP (DCIP) has the potential to enhance static detection, and also monitoring of DNAPL source zones (SZs). The objective of this study is to assess DCIP for characterizing and monitoring DNAPL SZs at the field‐scale. A new DNAPL‐DCIP numerical model was first developed that couples a 3D multiphase flow model, which simulates DNAPL release and remediation scenarios, with a 3D DCIP model, which calculates the corresponding resistivity and chargeability response. The sensitivity of the DCIP response to key DNAPL and soil properties was then analyzed at a single subsurface location, closely matching previous theoretical and experimental observations. Finally, a field‐scale simulation of DNAPL release and remediation was conducted, with simultaneous mapping by DCIP surveys. Results demonstrate that chargeability can provide enhanced understanding of the lithological distribution that controls the variability in the DNAPL SZ, with time‐lapse resistivity being used to monitor DNAPL mass changes during SZ remediation. This numerical study suggests that combined DCIP can be valuable for site characterization and time‐lapse monitoring performance at DNAPL‐impacted sites.