Asymmetric readout error
While its hard to make strong claims about the noise characteristics of any given quantum device, one explanation for what you're observing is readout error. For superconducting qubits readout error tends to be asymmetric: the probability $p(0|1)$ of observing a "0" after performing measurement on a computational basis state $|1\rangle$ tends to be higher than $p(1|0)$, the probability of observing "1" after measuring $|0\rangle$. You can verify this for yourself by looking at the Prob meas0 prep1
and Prob meas1 prep0
for the different systems offered by IBM.
Physically, this makes sense since usually $|0\rangle$ is defined as the ground state while $|1\rangle$ is an excited state. During the measurement process (for superconducting qubits this is usually on the order of $\sim 1\, \mu s$, which is quite long compared to gate durations) there is a tendency for the excited state $|1\rangle$ to relax to its ground state during the measurement process. However, as long as the the environment (i.e. dilution refrigerator) is cool enough the probability of exciting the ground state $|0\rangle$ tends to be much lower. Considering just thermalization effects, a back-of-the-envelope approach to this is to compute the Boltzmann ratio, which describes the equilibrium distribution for a two level system in terms of environmental temperature:
\begin{equation}
\frac{p_0}{p_1} = \exp\left(hf_{01} / kT\right) \tag{1}
\end{equation}
Plugging in some very rough values for superconducting qubits, $T\approx 10 \, \text{mK}$, $f_{01} \gtrsim 1 \, \text{GHz}$, Equation (1) gives something like $p_0 \gtrsim 100 p_1$, showing a strong tendency towards relaxation due to thermal effects (on a long enough timescale!). Of course this relaxation tendency will compete with other decoherence effects that may excite the qubit (e.g. $T_1$ decay will tend to depolarize your qubit state towards the maximally mixed state, which won't actually change the statistics you're trying to observe).
Readout error mitigation
It is useful to consider a model for readout error to check whether this could be the source of your issue. Let $\mathbf{p}' = [p_0, p_1]^T$ be the vector of observed probabilities for the bits "0" and "1". A rudimentary model for readout error is to consider the measured bits for each qubit as resulting from a stochastic process described by the equation
\begin{equation}
\mathbf{p}' = R \mathbf{p} \tag{2}
\end{equation}
where $\mathbf{p}$ is the vector of "true" bit probabilities that we would have observed with perfect measurement, and $R$ is the matrix
\begin{equation}
R = \begin{pmatrix}
p(0 |0) & p(0 |1) \\
p(1 |0) & p(1 |1)
\end{pmatrix}\tag{3}
\end{equation}
From this expression we can think about how readout error would affect the observed distribution when the true distribution for the state $|+\rangle$ given as $\mathbf{p}_{|+\rangle} = [0.5, 0.5]^T$:
\begin{equation}
\mathbf{p}_{|+\rangle}' = R
\begin{pmatrix}
0.5 \\
0.5\end{pmatrix}
=\frac{1}{2} \begin{pmatrix}
1-p(1|0) + p(0|1) \\
1 - p(0|1) + p(1|0)
\end{pmatrix}\tag{4}
\end{equation}
where I've simplified using $p(0|0) + p(1|0) = 1$. From this we can see that $p(0|1) \geq p(1|0)$ will result in $p_0' \geq p'_1$ - the observed distributions will contain too many zeros. We can further check how well readout error explains the discrepancies in the measurements by computing
$$
\mathbf{p} = R^{-1} \mathbf{p}'\tag{5}
$$
This depends a lot on the accuracy of device diagnostics, when you performed you experiments, etc. but computing Equation (4) using $p(0|1)$ and $p(1|0)$ diagnoistics taken from IBM's page just now (averaged over each device) and the probabilities $\mathbf{p}'$ you provided for Santiago and Lima I found $p_0$ for the corrected distributions to be 0.49 and 0.51 respectively, which suggests that readout error could explain much of the discrepancy in the observed distributions.
For better accuracy this calculation should be run using recent readout error diagnostics for $p(0|1)$ and $p(1|0)$ collected at the same time you run your experiment. I believe IBM also offers some readout error mitigation tools in their software package but you should make sure it uses fresh diagnostic data since these values tend to drift on very short timescales.