AbstractIn this work, we present a statistical study of substorms covering a five‐year period 2016–2020. Substorm phases were identified from time series of the SuperMAG AL (SML) index using a list of 5,077 previously identified substorm onsets, the SML peak value marking transition from expansion to recovery phase, and the recovery identified as return to activity less than −100 nT in the SML index. Magnetic field observations from THEMIS, RBSP, and MMS missions were used to study the magnetotail characteristics during the substorm evolution. A superposed epoch analysis indicates that the substorm onset occurs almost simultaneously with a few minutes of uncertainty throughout the magnetotail, ranging from geostationary orbit to 20 RE. The onset in the transition region precedes the ground onset by a few minutes. The peak SML time coincides with the peak of the outer transition region ΔBZ, which suggests that the field‐aligned currents driving the SML activity arise from the outer transition region. Analysis of 2D maps of the tail magnetic field shows that the magnetotail current changes are limited to the center of the tail within |Y| < 10RE. The substorm recovery is fastest in the inner transition region and lasts longer when moving further out. We did not find major asymmetries in the substorm signatures associated with IMF BY or BZ.