I am searching for a reference to use to confirm that mediation with multiple IV through Hayes' PROCESS macro on SPSS is (or is not) based on multiple regression.
Absolutely. PROCESS is based on OLS regression/General Linear Modeling. You can use Hayes' book on PROCESS as reference (see link), if a reference is really needed. Andy Field's Discovering Statistics Using SPSS 4th ed. also discusses this from what I remember.
It's actually more of a "system" of regressions.
For example, for a simple mediation of X->M->Y, you have (excuse the non-math ability of the textbox):
M ~ b0M + b1M*X (so a simple regression with X as a predictor)
&
Y ~ b0Y + b1Y* X + b2Y*M (so a multiple regression with X and M as predictors)
Then the direct effect is b1Y , and the indirect is (b1M * b2Y).
I would use his book as reference (see link), though honestly, most (if not all) path analysis (mediation, moderation, etc.) with (assumed) normally distributed outcomes use OLS (e.g., the R package Lavaan), as the only option, or at least as default.
I am using both Hayes' and Field's books. I just couldn't find the answer to my particular question, as Hayes is really not understandable for me and Field do not cover mediation.
Oh yes, I think I've read through million of Hayes' publications, and he always says linear regression, but not really if it is simple or multiple (with multiple IVs).
It is multiple. I think he felt he was conveying that in the formulas, but it is never made explicit. That does mean the regression weights are contextual, not zero order.