“…Then, SHA is performed on ΔH, ΔD and ΔZ components at each geomagnetic colatitude based on dominant harmonics rank m in the range 0 ≤ m ≤ 4 to calculate the spherical harmonic coefficients. Consequently, the Sq current functions over the two longitudinal sectors were obtained for each IQD at every hour of LT and geomagnetic colatitude, respectively based on Equation (e.g., Owolabi et al.,
2021; Yamazaki & Maute,
2017).
where
is assumed in the calculations because
,
is the associated Legendre function of geomagnetic colatitude,
and
are the degree and order at which the series of harmonics function is truncated,
…”