Subsurface simulator PFLOTRAN is internally coupled with geophysical simulator E4D. E4D runs as a module to PFLOTRAN, using a dedicated set of parallel compute nodes. The coupled code enables the simulation of changes in bulk electrical conductivity. The coupled code enables the simulation of time-lapse electrical resistivity surveys. A simulation of geophysical monitoring of groundwater/river water interaction is demonstrated.