In this chapter we describe a set of finite element formulations employed to approximate the solution to an extension to partially saturated porous media of Pride's equations, which is one of the most widespread theoretical frames that model the coupled electromagnetics and acoustics of saturated porous media. We also show how these numerical algorithms can be implemented on multi-core computers, and analyse their performance, observing that the parallel efficiency of the algorithms decays slower when the resources of each computing core are used as close to their maximum as possible. Further, we perform a study of two dimensional PSVTM electroseismic conversions in a time-lapse carbon dioxide geological deposition monitoring scenario, setting the electromagnetic source in depth, and measuring the seismic responses both in wells and at the surface. We observe that, contrary to what is known from standard field settings, the interface responses play no important role in elucidating the CO 2 plume behaviour, and that the in-depth converted signals convey information about the carbon dioxide saturation in the plume.