In this poster we describe a statistical analysis of extreme temporal changes in the horizontal component of the geomagnetic field (dB/dt) – an important indicator of geomagnetically induced currents. Extreme value theory is applied to data from 125 magnetometers in the global SuperMAG archive – with an average of 28 years of measurements per site – to determine return levels (RL) of |dB/dt| expected over periods of 100 years or more. This is achieved by fitting generalized Pareto (GP) distributions to declustered measurements of |dB/dt| above a 99.97-percentile threshold. Since large fluctuations are driven by diverse magneto-ionospheric driving processes (substorm expansions, sudden commencements, Pc5 ULF waves, etc.), the occurrence rate and high percentiles of |dB/dt| vary with geomagnetic latitude, magnetic local time (MLT), season, and with the compass direction of the fluctuation, dB. The interplanetary magnetic field orientation also exerts a strong influence on these patterns of occurrence, which we present alongside examples of fitted smoothing functions (splines and limited-order spherical harmonics and polynomials). By adapting statistical methods developed for the directional analysis of extreme ocean wave heights, we show how |dB/dt| GP distributions and RLs are generated for data sets sectored by compass direction (or equivalently by season, or MLT sector) and then combined and compared with GP distributions (and RLs) that ignore directionality or seasonal/diurnal variation.