We perform various lattice numerical analyses with the energy-momentum tensor (EMT). The EMT operator defined on the basis of the gradient flow is used in these analyses. We explore the spatial distribution of the stress tensor in static quark-anti-quark systems and thermodynamic quantities at nonzero temperature, as well as the correlation functions of EMT. The stress tensor distribution is also studied in the Abelian-Higgs model, which is compared with the lattice result.