IMD is what THD looks like with multiple tones. I don't know why they are treated as being intrinsically different, when all that's changed is the input signal.
Non-linearity of the system transfer function is the underlying, common cause. THD and IMD is what it looks like when passing one or multiple tones, respectively, through the same non-linearity.