  Subsystem symmetry protected topological (SSPT) order is a type of quantum
order that is protected by symmetries acting on lower-dimensional subsystems of
the entire system. In this paper, we show how SSPT order can be characterized
and detected by a constant correction to the entanglement area law, similar to
the topological entanglement entropy. Focusing on the paradigmatic
two-dimensional cluster phase as an example, we use tensor network methods to
give an analytic argument that almost all states in the phase exhibit the same
correction to the area law, such that this correction may be used to reliably
detect the SSPT order of the cluster phase. Based on this idea, we formulate a
numerical method that uses tensor networks to extract this correction from
ground-state wave functions. We use this method to study the fate of the SSPT
order of the cluster state under various external fields and interactions, and
find that the correction persists unless a phase transition is crossed, or the
subsystem symmetry is explicitly broken. Surprisingly, these results uncover
that the SSPT order of the cluster state persists beyond the cluster phase,
thanks to a new type of subsystem time-reversal symmetry. Finally, we discuss
the correction to the area law found in three-dimensional cluster states on
different lattices, indicating rich behavior for general subsystem symmetries
