In this paper we establish a higher integrability result up to the boundary of weak solutions to doubly nonlinear parabolic systems. We show that the spatial gradient of a weak solution with vanishing lateral boundary values is integrable to a larger power than the natural power p, where the statement holds for parameters in the subquadratic case $$ \max \lbrace \frac{2n}{n+2}, 1 \rbrace < p \le 2$$
max
{
2
n
n
+
2
,
1
}
<
p
≤
2
.