Mass renormalization actually has classical underpinning, going back to the 19th century, and what you see in quantum field theory is a quantized form of it. It is not something introduced by quantum theory, but inherited by quantum theory from classical physics.
As for field and charge renormalization: the question you should be asking is now "why" but "what" is being renormalized? The "bare" versus "dressed" distinction for charges and fields are those taken with respect to them referred to in Gaussian units. This unfortunate choice of units, which today survives only in the theoretical literature as one of its few remaining hold-outs, leads to confusion, false attribution: this very issue being a case in point.
This is a follow-up on the comments I made in my reply here
Gaussian Unit Of Charge And Force
For charges: a Gaussian charge, which we'll denote $q$ - as in the reply - is related to ordinary charge, which we'll denote $e$, as measured (say) in Coulombs by $q = e/\sqrt{4πε_0}$, where $ε_0$ is the permittivity of the vacuum.
When we associate a renormalization coefficient $\sqrt{Z_3}$ with $q$, what we're actually doing is renormalizing $ε_0$ by $1/Z_3$. It's not the charge that has the "bare" versus "dressed" distinction and is not the charge that's being renormalized, but the permittivity and the vacuum, itself!
For the electromagnetic field, denoting the Gaussian version of the potential by $a$ and the Si version by $A$, the relation between the two is $a = \sqrt{4πε_0}c A$. So, when we associate a renormalization coefficient - the same one (in fact) - $\sqrt{Z_3}$ with $a$, what we're actually doing is associating $1/Z_3$ with $ε_0$ and renormalizing it.
In both cases, it's the permittivity that's being renormalized, not the fields or charges. The use of Gaussian units obscures this and confuses the issue.
An analogous situation occurs with other gauge fields. To see what the analogue is, write out the Maxwell-Lorentz Lagrangian density in ordinary, SI, units:
$$_{ML} = -¼ ε_0 c \sqrt{|g|} g^{μρ} g^{νσ} F_{μν} F_{ρσ}.$$
By comparison, the Yang-Mills Lagrangian density has the form:
$$_{YM} = -¼ k_{ab} \sqrt{|g|} g^{μρ} g^{νσ} F^a_{μν} F^b_{ρσ},$$
where the indices $a,b,...$ are taken with respect to a basis of the Lie algebra associated with the gauge field (e.g. color $SU(3)$ for the QCD field), and $k_{ab}$ is the adjoint-invariant metric associated with the gauge group and its Lie algebra. It is usually (but not always) suppressed in the theoretical literature which - correspondingly - plays fast and loose with its careless handling of the indices failing to distinguish covariant from contravariant, or even recognize that there is a difference! For simple gauge groups, like $SU(3)$, an adjoint invariant metric is determined uniquely, up to scale, as a multiple of the Killing metric, $κ_{ab} = \sum_{0≤c,d<n} f^c_{ad} f^d_{bc}$, where the $f^c_{ab}$ are the structure coefficients associated with the Lie algebra basis $\left(Y_a: 0 ≤ a < n\right)$, i.e. the coefficients in the Lie bracket relations $\left[Y_a, Y_b\right] = \sum_{0≤c<n} f^c_{ab} Y_c$. (For $SU(3)$, $n = 8$.) The extra coefficient is the "coupling coefficient" $g$, itself, with $k = κ/g^2$ up to signs and factors (I forget whether there's a $4π$ or not) and is frequently, in the theoretical physics literature, migrated over into the fields as part of their definition by rescaling them in much the same way that the Gaussian $a$ is rescaled relative to SI $A$.
The role analogous to the vacuum permittivity $ε_0$ in $ε_0 c$ is played by $k_{ab}$ and, by extension, by its coefficient $g$. So, we could consider $k_{ab}/c$ as the gauge version of "vacuum permittivity". The role played by $g$ is thus analogous to $1/\sqrt{4πε_0}$.
So, when we're talking about the "running of the couplings", once again, what we're actually referring to is the renormalization of the vacuum and the running of its coefficients: here, the metric $k_{ab}$ and its coefficient $g$.
For semi-simple gauge groups, and products of them with Abelian gauge groups, the chief example in mind being $SU(3)⊗SU(2)⊗U(1)$, the metric is a direct sum of the metrics associated with each factor. So, you have several metric independent coefficients in place of just one: here $g_s$ for color $SU(3)$, $g$ for isospin $SU(2)$ and $g'$ for hypercharge $U(1)$.