Solar Wind Conditions Research Articles

Global MHD Simulations of the Earth's Bow Shock Shape and Motion Under Variable Solar Wind Conditions

AbstractEmpirical models of the Earth's bow shock are often used to place in situ measurements in context and to understand the global behavior of the foreshock/bow shock system. They are derived statistically from spacecraft bow shock crossings and typically treat the shock surface as a conic section parameterized according to a uniform solar wind ram pressure, although more complex models exist. Here a global magnetohydrodynamic simulation is used to analyze the variability of the Earth's bow shock under real solar wind conditions. The shape and location of the bow shock is found as a function of time, and this is used to calculate the shock velocity over the shock surface. The results are compared to existing empirical models. Good agreement is found in the variability of the subsolar shock location. However, empirical models fail to reproduce the two‐dimensional shape of the shock in the simulation. This is because significant solar wind variability occurs on timescales less than the transit time of a single solar wind phase front over the curved shock surface. Empirical models must therefore be used with care when interpreting spacecraft data, especially when observations are made far from the Sun‐Earth line. Further analysis reveals a bias to higher shock speeds when measured by virtual spacecraft. This is attributed to the fact that the spacecraft only observes the shock when it is in motion. This must be accounted for when studying bow shock motion and variability with spacecraft data.

An operational solar wind prediction system transitioning fundamental science to operations

We present in this paper an operational solar wind prediction system. The system is an outcome of the collaborative efforts between scientists in research communities and forecasters at Space Environment Prediction Center (SEPC) in China. This system is mainly composed of three modules: (1) a photospheric magnetic field extrapolation module, along with the Wang-Sheeley-Arge (WSA) empirical method, to obtain the background solar wind speed and the magnetic field strength on the source surface; (2) a modified Hakamada-Akasofu-Fry (HAF) kinematic module for simulating the propagation of solar wind structures in the interplanetary space; and (3) a coronal mass ejection (CME) detection module, which derives CME parameters using the ice-cream cone model based on coronagraph images. By bridging the gap between fundamental science and operational requirements, our system is finally capable of predicting solar wind conditions near Earth, especially the arrival times of the co-rotating interaction regions (CIRs) and CMEs. Our test against historical solar wind data from 2007 to 2016 shows that the hit rate (HR) of the high-speed enhancements (HSEs) is 0.60 and the false alarm rate (FAR) is 0.30. The mean error (ME) and the mean absolute error (MAE) of the maximum speed for the same period are −73.9 km s−1 and 101.2 km s−1, respectively. Meanwhile, the ME and MAE of the arrival time of the maximum speed are 0.15 days and 1.27 days, respectively. There are 25 CMEs simulated and the MAE of the arrival time is 18.0 h.

Ionospheric Irregularities at Mars Probed by MARSIS Topside Sounding

AbstractThe upper ionosphere of Mars contains a variety of perturbations driven by solar wind forcing from above and upward propagating atmospheric waves from below. Here we explore the global distribution and variability of ionospheric irregularities around the exobase at Mars by analyzing topside sounding data from the Mars Advanced Radar for Subsurface and Ionosphere Sounding (MARSIS) instrument on board Mars Express. As irregular structure gives rise to off‐vertical echoes with excess propagation time, the diffuseness of ionospheric echo traces can be used as a diagnostic tool for perturbed reflection surfaces. The observed properties of diffuse echoes above unmagnetized regions suggest that ionospheric irregularities with horizontal wavelengths of tens to hundreds of kilometers are particularly enhanced in the winter hemisphere and at high solar zenith angles. Given the known inverse dependence of neutral gravity wave amplitudes on the background atmospheric temperature, the ionospheric irregularities probed by MARSIS are most likely associated with plasma perturbations driven by atmospheric gravity waves. Though extreme events with unusually diffuse echoes are more frequently observed for high solar wind dynamic pressures during some time intervals, the vast majority of the diffuse echo events are unaffected by varying solar wind conditions, implying limited influence of solar wind forcing on the generation of ionospheric irregularities. Combination of remote and in situ measurements of ionospheric irregularities would offer the opportunity for a better understanding of the ionospheric dynamics at Mars.

On fast and slow Earth’s magnetospheric dynamics during geomagnetic storms: a stochastic Langevin approach

The Earth’s magnetosphere responds to the external changes of interplanetary magnetic field and solar wind conditions showing a multiscale dynamics, manifesting in the occurrence of fluctuations over a very wide range of timescales. Here, using an approach based on a Langevin/Fokker-Planck description we investigate the nature of the fast (short-) and slow (long-timescale) fluctuations of SYM-H index during geomagnetic storms. The results point towards a different origin of the fast (τ < 200 min) and slow (τ > 200 min) fluctuations, which are characterized by state functions of different nature. In detail, the state function associated with the slow dynamics shows the evidence of the occurrence of first-order-like topological phase transition during the different phases of a geomagnetic storm, while the fast dynamics seems to be characterized by a quasi-invariant quadratic state function. A modeling in terms of stochastic Langevin equation is discussed and the relevance of our results in the framework of Space Weather studies is outlined.

A new short-term forecasting model for the total electron content storm time disturbances

This paper aims to introduce a new model for the short-term forecast of the vertical Total Electron Content (vTEC). The basic idea of the proposed model lies on the concept of the Solar Wind driven autoregressive model for Ionospheric short-term Forecast (SWIF). In its original version, the model is operationally implemented in the DIAS system (http://dias.space.noa.gr) and provides alerts and warnings for upcoming ionospheric disturbances, as well as single site and regional forecasts of the foF2 critical frequency over Europe up to 24 h in advance. The forecasts are driven by the real time assessment of the solar wind conditions at ACE location. The comparative analysis of the variations in foF2 and vTEC during eleven geomagnetic storm events that occurred in the present solar cycle 24 reveals similarities but also differences in the storm-time response of the two characteristics with respect to the local time and the latitude of the observation point. Since the aforementioned dependences drive the storm-time forecasts of the SWIF model, the results obtained here support the upgrade of the SWIF’s modeling technique in forecasting the storm-time vTEC variation from its onset to full development and recovery. According to the proposed approach, the vTEC storm-time response can be forecasted from 1 to 12–13 h before its onset, depending on the local time of the observation point at storm onset at L1. Preliminary results on the assessment of the performance of the proposed model and further considerations on its potential implementation in operational mode are also discussed.

Dawn‐Dusk Auroral Oval Oscillations Associated With High‐Speed Solar Wind

AbstractWe report evidence of global‐scale auroral oval oscillations in the millihertz range, using global auroral images acquired from the Ultraviolet Imager on board the decommissioned Polar satellite and concurrent solar wind measurements. On the basis of two events (15 January 1999 and 6 January 2000) studied, it is found that (1) quasi‐periodic auroral oval oscillations (~3 mHz) can occur when solar wind speeds are high at northward or southward interplanetary magnetic field turning, (2) the oscillation amplitudes range from a few to more than 10° in latitudes, (3) the oscillation frequency is the same for each event irrespective of local time and without any azimuthal phase shift (i.e., propagation), (4) the auroral oscillations occur in phase within both the dawn and dusk sectors but 180° out of phase between the dawn and dusk sectors, and (5) no micropulsations on the ground match the auroral oscillation periods. While solar wind conditions favor the growth of the Kelvin‐Helmholtz (K‐H) instability on the magnetopause as often suggested, the observed wave characteristics are not consistent with predictions for K‐H waves. The in‐phase and out‐of‐phase features found in the dawn‐dusk auroral oval oscillations suggest that wiggling motions of the magnetotail associated with fast solar winds might be the direct cause of the global‐scale millihertz auroral oval oscillations.

A comparative study of sawtooth events and substorm onsets triggered by interplanetary shocks

comparison of solar wind conditions, geomagnetic response, and auroral boundary movement, using a similar number of sawtooth events and shock induced substorm triggers is carried out. 81% of the sawtooth onsets were triggered at low latitudes compared to 33% for onsets during shocks. Results of superposed epoch analysis indicated that the mean interplanetary magnetic field (IMF) B z remained strongly southward during sawtooth events while it was southward during the loading period and turned northward 20 minutes prior to the onset, during the shock triggered events. During both the event types, the mean energy available in the solar wind was found to be above the substorm threshold level. Relatively high levels of solar wind density indicated to the magnetosphere being compressed to a larger degree during substorms initiated by shocks. The high latitude indices were elevated during sawtooth events with 2-3 hr fluctuations was observed in the AL and PCN indices but AU index after the onset remained at similar levels. SYM-H remained below storm time values during shock related onsets but was much stronger and remained above the storm threshold during sawtooth events, that they occurred embedded in a geomagnetic storm. Comparatively higher mean values of ASYM-H indicated to a much stronger asymmetric ring current during sawtooth events. Mean boundary locations during sawtooth events were located a few degrees equatorward of those during the shock related onsets. Oval in the dawn, dusk, and midnight sectors was much thicker during sawtooth events with clear widening around the onset during both event groups but the thickness in noon sector was similar and remained steady during the epoch window during both the event groups. Elevated driving during sawtooth events was also indicative in a greater stretching and relaxation of the magnetic field lines at geosynchronous orbit in the midnight sector.

The Science of Sungrazers, Sunskirters, and Other Near-Sun Comets

This review addresses our current understanding of comets that venture close to the Sun, and are hence exposed to much more extreme conditions than comets that are typically studied from Earth. The extreme solar heating and plasma environments that these objects encounter change many aspects of their behaviour, thus yielding valuable information on both the comets themselves that complements other data we have on primitive solar system bodies, as well as on the near-solar environment which they traverse. We propose clear definitions for these comets: We use the term near-Sun comets to encompass all objects that pass sunward of the perihelion distance of planet Mercury (0.307 AU). Sunskirters are defined as objects that pass within 33 solar radii of the Sun’s centre, equal to half of Mercury’s perihelion distance, and the commonly-used phrase sungrazers to be objects that reach perihelion within 3.45 solar radii, i.e. the fluid Roche limit. Finally, comets with orbits that intersect the solar photosphere are termed sundivers. We summarize past studies of these objects, as well as the instruments and facilities used to study them, including space-based platforms that have led to a recent revolution in the quantity and quality of relevant observations. Relevant comet populations are described, including the Kreutz, Marsden, Kracht, and Meyer groups, near-Sun asteroids, and a brief discussion of their origins. The importance of light curves and the clues they provide on cometary composition are emphasized, together with what information has been gleaned about nucleus parameters, including the sizes and masses of objects and their families, and their tensile strengths. The physical processes occurring at these objects are considered in some detail, including the disruption of nuclei, sublimation, and ionisation, and we consider the mass, momentum, and energy loss of comets in the corona and those that venture to lower altitudes. The different components of comae and tails are described, including dust, neutral and ionised gases, their chemical reactions, and their contributions to the near-Sun environment. Comet-solar wind interactions are discussed, including the use of comets as probes of solar wind and coronal conditions in their vicinities. We address the relevance of work on comets near the Sun to similar objects orbiting other stars, and conclude with a discussion of future directions for the field and the planned ground- and space-based facilities that will allow us to address those science topics.

STUDY OF THE ENERGY POTENTIAL OF RENEWABLE ENERGY SOURCES IN THE REGION OF GABROVO FOR USE IN HYBRID WIND-PHOTOVOLTAIC STAND-ALONE SYSTEM

This paper presents the measurements conducted for two periods (winter and summer) in the region of Gabrovo, Bulgaria for solar radiation and wind conditions using the digital measuring equipment. The aim is to evaluate and predict the possible production of electricity from already installed low power hybrid system. The conclusion is that a combination of a wind generator and a photovoltaic system in the hybrid system for independent power supply has a good economic effect and allows fuller use of energy from renewable sources.

Long‐Term Variability of Jupiter's Magnetodisk and Implications for the Aurora

AbstractObservations of Jupiter's UV auroral emissions collected over several years show that the ionospheric positions of the main emission and the Ganymede footprint can vary by as much as 3° in latitude. One explanation for this shift is a change of Jupiter's current sheet current density, which would alter the amount of field line stretching and displace the ionospheric mapping of field lines from a given radial distance in the magnetosphere. In this study we measure the long‐term variability of Jupiter's magnetodisk using Galileo magnetometer data collected from 1996 to 2003. Using the Connerney et al. (1981) current sheet model, we calculate the current sheet density parameter that gives the best fit to the data from each orbit and find that the current density parameter varies by about 15% of its average value during the Galileo era. We investigate possible relationships between the observed current sheet variability and quantities such as Io's plasma torus production rate inferred from volcanic activity and external solar wind conditions extrapolated from data at 1 AU but find only a weak correlation. Finally, we trace Khurana (1997) model field lines to show that the observed changes in Jupiter's current sheet are sufficient to shift the ionospheric footprint of Ganymede and main auroral emission by a few degrees of latitude, consistent with the magnitude of auroral variability observed by Hubble Space Telescope (HST). However, we find that the measured auroral shifts in HST images are not consistent with concurrent changes in the current density parameter measured by Galileo.

Mars Under Primordial Solar Wind Conditions: Mars Express Observations of the Strongest CME Detected at Mars Under Solar Cycle #24 and its Impact on Atmospheric Ion Escape

AbstractAn extremely strong Coronal Mass Ejection (CME) impacted Mars on 12 July 2011, while the Mars Express spacecraft was present inside the nightside ionosphere. Estimated solar wind density and speed during the event are 39 particles cm−3 and 730 km/s, corresponding to nominal solar wind flux at Mars when the solar system was ∼1.1 Ga old. Comparing with expected average atmospheric heavy ion fluxes under similar XUV conditions, the CME impact is found to have no significant effect on the escape rate 3.3 × 1024s−1, with an upper limit at 1025s−1 if the observed tail contraction is not taken into account. On the subsequent orbit, 7 h later after magnetosphere response, fluxes were only 2.4% of average. As such, even under primordial solar wind conditions we are unable to find support for a strong solar wind‐driven ion escape, rather the main effect appears to be acceleration of the escaping ions by ×10–×20 typical characteristic energy.

Investigating the Development of Localized Neutral Density Increases During the 24 August 2005 Geomagnetic Storm

AbstractWe investigate localized thermospheric neutral density increases (or density spikes) developed during the 24 August 2005 geomagnetic storm and utilize a multi‐instrument database for monitoring the underlying thermospheric, ionospheric, magnetospheric, and solar wind conditions. We study five scenarios occurring under eastward/westward BY dominations and demonstrating the development of density spikes during various events. These scenarios show (1) a northern density spike associated with enhanced antisunward polar flows driven by dayside merging along old‐open field lines during the initial phase, (2) some southern density spikes associated with enhanced sunward auroral flows caused by auroral brightening events during the storm initial phase and main phase, (3) two pairs of southern dayside‐nightside density spike configuration associated with enhanced antisunward polar flows caused by dayside merging and nightside magnetotail reconnections along old‐open field lines during the underlying substorm activity, and (4) a series of enhanced polar sunward flows developed in the polar cap driven by different processes during the recovery phase. We provide a detailed analysis of these events including the specification of underlying field aligned current systems, flow channel types, and auroral forms. Observational results demonstrate the various density spikes and their associated auroral forms and/or flow channels and thus provide a better understanding of their development during this intense geomagnetic storm.

Probabilistic Solar Wind Forecasting Using Large Ensembles of Near-Sun Conditions With a Simple One-Dimensional "Upwind" Scheme.

Long lead‐time space‐weather forecasting requires accurate prediction of the near‐Earth solar wind. The current state of the art uses a coronal model to extrapolate the observed photospheric magnetic field to the upper corona, where it is related to solar wind speed through empirical relations. These near‐Sun solar wind and magnetic field conditions provide the inner boundary condition to three‐dimensional numerical magnetohydrodynamic (MHD) models of the heliosphere out to 1 AU. This physics‐based approach can capture dynamic processes within the solar wind, which affect the resulting conditions in near‐Earth space. However, this deterministic approach lacks a quantification of forecast uncertainty. Here we describe a complementary method to exploit the near‐Sun solar wind information produced by coronal models and provide a quantitative estimate of forecast uncertainty. By sampling the near‐Sun solar wind speed at a range of latitudes about the sub‐Earth point, we produce a large ensemble (N = 576) of time series at the base of the Sun‐Earth line. Propagating these conditions to Earth by a three‐dimensional MHD model would be computationally prohibitive; thus, a computationally efficient one‐dimensional “upwind” scheme is used. The variance in the resulting near‐Earth solar wind speed ensemble is shown to provide an accurate measure of the forecast uncertainty. Applying this technique over 1996–2016, the upwind ensemble is found to provide a more “actionable” forecast than a single deterministic forecast; potential economic value is increased for all operational scenarios, but particularly when false alarms are important (i.e., where the cost of taking mitigating action is relatively large).

Weakly turbulent plasma processes in the presence of inverse power-law velocity tail population

Observations show that plasma particles in the solar wind frequently display power-law velocity distributions, which can be isotropic or anisotropic. Particularly, the velocity distribution functions of solar wind electrons are frequently modeled as a combination of a background Maxwellian distribution and a non-thermal distribution which is known as the “halo” distribution. For fast solar wind conditions, highly anisotropic field-aligned electrons, denominated as the “strahl” distribution, are also present. Motivated by these observations, the present paper considers a tenuous plasma with Maxwellian ions, and electrons described by a summation of an isotropic Maxwellian distribution and an isotropic Kappa distribution. The formalism of weak turbulence theory is utilized in order to discuss the spectra of electrostatic waves that must be present in such a plasma, satisfying the conditions of quasi-equilibrium between the processes of spontaneous fluctuations and of induced emission. The kappa index and relative density of the Kappa electron distribution are varied. By taking into account the effects due to electromagnetic waves into the weak turbulence formalism, we investigate the electromagnetic spectra that satisfy the conditions of “turbulent equilibrium,” and also the time evolution of the wave spectra and of the electron distribution, which occurs in the case of the presence of an electron beam in the electron distribution.

Flows, Fields, and Forces in the Mars‐Solar Wind Interaction

AbstractWe utilize suprathermal ion and magnetic field measurements from the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, organized by the upstream magnetic field, to investigate the morphology and variability of flows, fields, and forces in the Mars‐solar wind interaction. We employ a combination of case studies and statistical investigations to characterize the interaction in both quasi‐parallel and quasi‐perpendicular regions and under high and low solar wind Mach number conditions. For the first time, we include a detailed investigation of suprathermal ion temperature and anisotropy. We find that the observed magnetic fields and suprathermal ion moments in the magnetosheath, bow shock, and upstream regions have observable asymmetries controlled by the interplanetary magnetic field, with particularly large asymmetries found in the ion parallel temperature and anisotropy. The greatest temperature anisotropies occur in quasi‐perpendicular regions of the magnetosheath and under low Mach number conditions. These results have implications for the growth and evolution of wave‐particle instabilities and their role in energy transport and dissipation. We utilize the measured parameters to estimate the average ion pressure gradient, J × B, and v × B macroscopic force terms. The pressure gradient force maintains nearly cylindrical symmetry, while the J × B force has larger asymmetries and varies in magnitude in comparison to the pressure gradient force. The v × B force felt by newly produced planetary ions exceeds the other forces in magnitude in the magnetosheath and upstream regions for all solar wind conditions.

Global Structure and Sodium Ion Dynamics in Mercury's Magnetosphere With the Offset Dipole

AbstractWe conducted global magnetohydrodynamics (MHD) simulation of Mercury's magnetosphere with the dipole offset, which was revealed by MESSENGER (Mercury Surface, Space Environment, Geochemistry, and Ranging) observations, in order to investigate its global structure under northward interplanetary magnetic field conditions. Sodium ion dynamics originating from the Mercury's exosphere is also investigated based on statistical trajectory tracing in the electric and magnetic fields obtained from the MHD simulations. The results reveal a north‐south asymmetry characterized by open field lines around the southern polar region and northward deflection of the plasma sheet in the far tail. The asymmetry of magnetic field structure near the planet drastically affects trajectories of sodium ion and thus their pressure distributions and precipitation pattern onto the planet. Weaker magnetic field strength in the southern hemisphere than in the north increases ion loss by precipitation onto the planetary surface in the southern hemisphere. The “sodium ring,” which is formed by high‐energy sodium ions drifting around the planet, is also found in the vicinity of the planet. The sodium ring is almost circular under nominal solar wind conditions. The ring becomes partial under high solar wind density, because dayside magnetosphere is so compressed that there is no space for the sodium ions to drift around. In both cases, the sodium ring is formed by sodium ions that are picked up, accelerated in the magnetosheath just outside the magnetopause, and reentered into the magnetosphere due to combined effects of finite Larmor radius and convection electric field in the dawnside magnetosphere.

Model‐Based Assessments of Magnetic Reconnection and Kelvin‐Helmholtz Instability at Jupiter's Magnetopause

AbstractThe interaction between the solar wind and Jupiter's magnetic field confines the planetary field to the largest magnetosphere in the Solar System. However, the full picture of when and where key processes operate at the magnetopause boundary of the system remains unclear. This is essential for testing understanding with observations and for determining the relative importance of different drivers of Jovian magnetospheric dynamics. Here we present a global analytical model of Jovian magnetopause conditions under steady state, which forms the basis of boundary process assessments. Sites of magnetic reconnection at Jupiter's magnetopause are expected to be in regions of sufficiently high magnetic shear across the boundary, controlled by the orientation of the interplanetary magnetic field. Reconnection rates are also most sensitive to changes in the highly variable IMF, followed by changes in the solar wind plasma mass density. The largest plasma flow shear across the boundary is in the equatorial dawn region, producing a region that is typically unstable to growth of the Kelvin‐Helmholtz (K‐H) instability. Compared to magnetopause reconnection site locations, this K‐H‐unstable region at dawn is less sensitive to changing conditions. Motion of K‐H boundary perturbations typically includes dawn‐to‐dusk motion across the subsolar region. Model‐predicted reconnection voltages are typically hundreds of kV but rely on steady solar wind conditions on a time scale that is longer than typical at Jupiter's orbit. How the reconnection voltage compares to the voltage applied due to the “viscous‐like” interaction involving K‐H instability remains unclear.

MMS Observations and Hybrid Simulations of Surface Ripples at a Marginally Quasi‐Parallel Shock

AbstractSimulations and observations of collisionless shocks have shown that deviations of the nominal local shock normal orientation, that is, surface waves or ripples, are expected to propagate in the ramp and overshoot of quasi‐perpendicular shocks. Here we identify signatures of a surface ripple propagating during a crossing of Earth's marginally quasi‐parallel (θBn∼45∘) or quasi‐parallel bow shock on 27 November 2015 06:01:44 UTC by the Magnetospheric Multiscale (MMS) mission and determine the ripple's properties using multispacecraft methods. Using two‐dimensional hybrid simulations, we confirm that surface ripples are a feature of marginally quasi‐parallel and quasi‐parallel shocks under the observed solar wind conditions. In addition, since these marginally quasi‐parallel and quasi‐parallel shocks are expected to undergo a cyclic reformation of the shock front, we discuss the impact of multiple sources of nonstationarity on shock structure. Importantly, ripples are shown to be transient phenomena, developing faster than an ion gyroperiod and only during the period of the reformation cycle when a newly developed shock ramp is unaffected by turbulence in the foot. We conclude that the change in properties of the ripple observed by MMS is consistent with the reformation of the shock front over a time scale of an ion gyroperiod.

High‐Latitude Neutral Mass Density Maxima

AbstractRecent studies have reported that thermospheric effects due to solar wind driving can be observed poleward of auroral latitudes. In these papers, the measured neutral mass density perturbations appear as narrow, localized maxima in the cusp and polar cap. They conclude that Joule heating below the spacecraft is the cause of the mass density increases, which are sometimes associated with local field‐aligned current structures, but not always. In this paper we investigate neutral mass densities measured by accelerometers on the CHAllenging Minisatellite Payload (CHAMP) and Gravity Recovery and Climate Experiment (GRACE) spacecraft from launch until years 2010 (CHAMP) and 2012 (GRACE), approximately 10 years of observations from each satellite. We extract local maxima in neutral mass densities over the background using a smoothing window with size of one quarter of the orbit. The maxima have been analyzed for each year and also for the duration of each set of satellite observations. We show where they occur, under what solar wind conditions, and their relation to magnetic activity. The region with the highest frequency of occurrence coincides approximately with the cusp and mantle, with little direct evidence of an auroral zone source. Our conclusions agree with the “hot polar cap” observations that have been reported and studied in the past.

Climatology of the Auroral Electrojets Derived From the Along‐Track Gradient of Magnetic Field Intensity Measured by POGO, Magsat, CHAMP, and Swarm

AbstractThe auroral electrojets (AEJs) are complex and dynamic horizontal ionospheric electric currents which form ovals around Earth's poles, being controlled by the morphology of the main magnetic field and the energy input from the solar wind interaction with the magnetosphere. The strength and location of the AEJ varies with solar wind conditions and the solar cycle but should also be controlled on decadal timescales by main field secular variation. To determine the AEJ climatology, we use data from four polar Low Earth Orbit magnetic satellite missions: POGO, Magsat, CHAMP, and Swarm. A simple estimation of the AEJ strength and latitude is made from each pass of the satellites, from peaks in the along‐track gradient of the magnetic field intensity after subtracting a core and crustal magnetic field model. This measure of the AEJ activity is used to study the response in different sectors of magnetic local time (MLT) during different seasons and directions of the interplanetary magnetic field (IMF). We find a season‐dependent hemispherical asymmetry in the AEJ response to IMF By, with a tendency toward stronger (weaker) AEJ currents in the north than the south during By>0 (By<0) around local winter. This effect disappears during local summer when we find a tendency toward stronger currents in the south than the north. The solar cycle modulation of the AEJ and the long‐term shifting of its position and strength due to the core field variation are presented as challenges to internal field modeling.