Let $J_{s,k}(P)$ be the number of integer solutions of
\begin{align*}
x_1^j+\cdots+x_s^j=y_1^j+\cdots+y_s^j \qquad (1\le j\le k)
\end{align*}
with $1\le x_i,y_i\le P$. For every $s,k\in\mathbb N$ and every $\varepsilon>0$,
\begin{align*}
J_{s,k}(P)\lesssim_{s,k,\varepsilon} P^\varepsilon\left(P^s+P^{2s-k(k+1)/2}\right).
\end{align*}