  This paper provides a local and global Calder\'on-Zygmund type estimate of a
weak solution to the parabolic double-phase system. The proof of local estimate
is based on comparison estimates and the scaling invariant property of the
parabolic double-phase system in the intrinsic cylinders of the stopping time
argument setting. For the proof of the global estimate, we have applied the
reflection and approximation techniques.
