Total Electron Content (TEC) is central to characterizing ionospheric response to solar and geomagnetic activity. Variations in TEC structures over time provide insight into underlying physical processes and inform monitoring of space weather events, which pose a risk to navigation and communication systems. JPL processed GNSS observations over 20 years provide a series of 15-minute Global Ionospheric Maps (GIMs) of spatial resolution 1x1 degree longitude/latitude. We translate these into geomagnetic coordinates centered about the sub-solar point and we isolate the top 1% of TEC values in each map to define High Density Regions (HDRs) of TEC. Image processing tools are used to develop an algorithm that detects and tracks these to compile a set of contiguous, uniquely labelled space-time TEC HDRs. We find that HDRs naturally divide into two populations by peak area, separated by a size of 8.0x106 km2, which is around the continental scale. These populations are studied for different storm conditions - quiet (Kp<4), moderate (4 ≤ Kp < 7) and extreme (Kp ≥ 7): small HDRs form primarily around four magnetic latitude bands and move roughly parallel to lines of constant magnetic latitude towards later MLT. Large HDRs form around the same latitude bands but follow more complex paths. The statistical nature of these results could be used in predictive ionospheric models and identify reproducible trends on these spatial/temporal scales.