Polar Cap Boundary Research Articles

Do the Throat Auroras Create Polar Cap Patches?

AbstractThroat auroras and polar cap patches are common phenomena in the polar ionosphere resulting from magnetosphere‐ionosphere coupling. A campaign was organized, with all‐sky imagers at Yellow River Station, the European Incoherent Scatter Svalbard Radar, and coordinated low‐altitude spacecraft observations. During periods of radial interplanetary magnetic field (IMF), observations showed that, as poleward moving throat auroras faded around the polar cap boundary, they linked to poleward moving ionization patches. The throat auroras were produced by soft‐electron precipitation associated with dayside magnetic reconnection. The red line emission intensity of throat auroras was correlated with dayside reconnection events. Dense plasma from lower latitudes was transported poleward via enhanced convection in the throat auroras to form patches. This is a potentially new formation mechanism for patches associated with throat auroras and magnetic reconnection for radial IMF. Moreover, the patches move anti‐sunward due to the E × B drift.

Estimating the Evolution of Sparks in the Partially Screened Gap of Pulsars from Subpulse Drifting

A novel scheme has been developed to show that the observed phase behavior associated with subpulse drifting from two pulsars, J1034−3224 and J1720−2933, can be used to obtain the magnetic field configuration in the partially screened gap (PSG). The outflowing plasma along the open magnetic field line region of pulsars is generated as a result of spark discharges in an inner acceleration region (IAR) above the polar cap. The IAR has been modeled as a PSG with a steady supply of positively charged ions emitted from the heated polar cap surface dominated by strong non-dipolar magnetic fields. In a PSG the sparks are tightly packed and constrained to be present along the polar cap boundary. The sparks lag behind the rotation of the star during their lifetimes. As a result, the sparking pattern evolves along two different directions in a clockwise and counterclockwise manner around a stationary central spark and can be associated with the observed phenomenon of subpulse drifting. PSR J1034−3224 has four prominent components and exhibits bi-drifting where alternate components show the opposite sense of drifting, while PSR J1720−2933 has a single component profile and shows systematically coherent drift bands. We show that the differences in their drifting behavior can be directly linked to different natures of the non-dipolar surface magnetic field configurations.

Storm time polar cap expansion: interplanetary magnetic field clock angle dependence

Abstract. It is well known that the polar cap, delineated by the open–closed field line boundary (OCB), responds to changes in the interplanetary magnetic field (IMF). In general, the boundary moves equatorward when the IMF turns southward and contracts poleward when the IMF turns northward. However, observations of the OCB are spotty and limited in local time, making more detailed studies of its IMF dependence difficult. Here, we simulate five solar storm periods with the coupled model consisting of the Open Geospace General Circulation Model (OpenGGCM) coupled with the Coupled Thermosphere Ionosphere Model (CTIM) and the Rice Convection Model (RCM), i.e., the OpenGGCM-CTIM-RCM, to estimate the location and dynamics of the OCB. For these events, polar cap boundary location observations are also obtained from Defense Meteorological Satellite Program (DMSP) precipitation spectrograms and compared with the model output. There is a large scatter in the DMSP observations and in the model output. Although the model does not predict the OCB with high fidelity for every observation, it does reproduce the general trend as a function of IMF clock angle. On average, the model overestimates the latitude of the open–closed field line boundary by 1.61∘. Additional analysis of the simulated polar cap boundary dynamics across all local times shows that the MLT of the largest polar cap expansion closely correlates with the IMF clock angle, that the strongest correlation occurs when the IMF is southward, that during strong southward IMF the polar cap shifts sunward, and that the polar cap rapidly contracts at all local times when the IMF turns northward.

On the morphology of electrostatic solitary waves in the Earth’s aurora

Electrostatic solitary waves (ESWs) have been detected in abundance in Space plasma observations, both by satellites in near-Earth plasma environments as well as by planetary missions, e.g. Cassini in Saturn or MAVEN in Mars. In their usual form, these are manifested as a bipolar electric field corresponding to a bell-shaped pulse in the electrostatic potential. Recent studies have suggested the existence of alternative forms of ESWs, including flat-top solitary waves (FTSWs) and supersolitary waves (SSWs), both of which are often encountered in Space observations such as in polar cap boundary layer, the auroral acceleration region and elsewhere. This article focuses on the existence and characterization of different types of electrostatic solitary waves in multicomponent Space plasmas. Relying on a multi-fluid plasma model, comprising two types of ions and two different electron populations, we have identified the conditions for existence of flat-top solitary waves and supersolitons, in contrast to “standard" solitary waves. Both ion species are models as cold fluids, for simplicity. Our analysis reveals that the coexistence of the two electron populations is pivotal for the formation of such non-standard electrostatic structures, and that their characteristic parameters (temperature, density ratio) plays a decisive role in their generation and structural characteristics. Nonetheless, while supersolitary waves may exist in a wide range of parameter values (as confirmed by earlier theoretical studies), it appears that flat-top solitary waves will occur in a narrow window in the parameter region, which may explain their scarce (but non-negligible) frequency of observation. Our theoretical findings confirm and validate the existence of alternative (non-conventional) ESW waveforms in auroral plasma (in addition to the ubiquitous bipolar electric field form), where such an electron coexistence is typically observed.

Cellular patterns and dry convection in textured dust storms at the edge of Mars North Polar Cap

We present a study of textured local dust storms that develop at the northern polar cap boundary on Mars springtime. We have used images obtained with VMC and HRSC cameras onboard Mars Express and MARCI on MRO to analyze dust storms captured from March to July 2019 (Ls = 350° in MY 34–Ls = 54° in MY 35). The textured storms grow in the longitude sector 150°E-210°E centered at latitude ~60°N and exhibit spiral, filamentary and compact shapes that change and evolve rapidly in a daily basis. The storms translate by prevailing east and southeast winds with speeds 15–45 ms−1. In some areas of their interiors they show organized clusters of cells formed typically by 100 elements with sizes ~5–30 km with a length/width ratio ~ 1.2–3 in the wind direction. The cells have elongated downwind tails with lengths 4–8 times the cell size. The cells top altitudes are ~6–11 km above their surroundings. We propose that the spirals grow as baroclinic vortices within a vertically sheared eastward jet present at this epoch in Mars due to the intense meridional temperature gradient at the polar cap edge. We show using a simple one-dimensional model that the cells can be produced by shallow dry convection with dust acting as the heating source to generate the updrafts. These patterns resemble those seen in laboratory experiments and on clouds in Earth's atmosphere and can serve to comparatively elucidate and discern the different mechanisms at work in each case.

Polar Cap Boundary Identification Using Redline Optical Data and DMSP Satellite Particle Data

AbstractThe location of the polar cap boundary is typically determined by using low‐orbit satellite measurements in which the boundary is identified by its unique signature of a sharp decrease in energy and particle flux poleward of the auroral oval. In principle, this decrease in precipitating particles should appear as a concomitant sharp change in auroral luminosity. Based on a few events, Blanchard et al. (1995, https://doi.org/10.1029/94JA02631) suggested that a dramatic gradient in redline aurora may also be an indicator of the polar cap boundary. In recent years, advances in capabilities and the deployment of ground‐based all‐sky imagers have ushered in a new era of auroral measurements. Auroral imaging has moved well beyond the capabilities of the instrumentation in the previous study in terms of both spatial and temporal resolution. We now have access to decades of optical data from arrays spanning a huge spatial range, enabling a fresh examination of the relationship between redline aurora, particle precipitation, and the polar cap open closed boundary. In this study, we use data from the Defense Meteorological Satellite Program satellites in conjunction with the University of Calgary's Redline Emission Geospace Observatory (630.0 nm) data to assess the viability of automated detection of the 2‐dimensional polar cap boundary. Our results exhibit good agreement between the optical and particle polar cap boundary and suggest that a luminosity in redline emission could serve as a reasonable proxy for the location of the electron poleward boundary, while providing both high temporal and spatial resolution maps of the open‐closed boundary.

Quantifying the Lobe Reconnection Rate During Dominant IMF By Periods and Different Dipole Tilt Orientations

AbstractLobe reconnection is usually thought to play an important role in geospace dynamics only when the Interplanetary Magnetic Field (IMF) is mainly northward. This is because the most common and unambiguous signature of lobe reconnection is the strong sunward convection in the polar cap ionosphere observed during these conditions. During more typical conditions, when the IMF is mainly oriented in a dawn‐dusk direction, plasma flows initiated by dayside and lobe reconnection both map to high‐latitude ionospheric locations in close proximity to each other on the dayside. This makes the distinction of the source of the observed dayside polar cap convection ambiguous, as the flow magnitude and direction are similar from the two topologically different source regions. We here overcome this challenge by normalizing the ionospheric convection observed by the Super Dual Aurora Radar Network (SuperDARN) to the polar cap boundary, inferred from simultaneous observations from the Active Magnetosphere and Planetary Electrodynamics Response Experiment (AMPERE). This new method enable us to separate and quantify the relative contribution of both lobe reconnection and dayside/nightside (Dungey cycle) reconnection during periods of dominating IMF By. Our main findings are twofold. First, the lobe reconnection rate can typically account for 20% of the Dungey cycle flux transport during local summer when IMF By is dominating and IMF Bz ≥ 0. Second, the dayside convection relative to the open/closed boundary is vastly different in local summer versus local winter, as defined by the dipole tilt angle.

Changes in the Magnetic Field Topology and the Dayside/Nightside Reconnection Rates in Response to a Solar Wind Dynamic Pressure Front: A Case Study

AbstractOne of the most significant observations associated with a sharp enhancement in solar wind dynamic pressure, , is the poleward expansion of the auroral oval and the closing of the polar cap. The polar cap shrinking over a wide range of magnetic local times (MLTs), in connection with an observed increase in ionospheric convection and the transpolar potential, led to the conclusion that the nightside reconnection rate is significantly enhanced after a pressure front impact. However, this enhanced tail reconnection has never been directly measured. We demonstrate the effect of a solar wind dynamic pressure front on the polar cap closure, and for the first time, measure the enhanced reconnection rate in the magnetotail, for a case occurring during southward background Interplanetary Magnetic Field (IMF) conditions. We use Polar Ultra‐Violet Imager (UVI) measurements to detect the location of the open‐closed field line boundary, and combine them with Assimilative Mapping of Ionospheric Electrodynamics (AMIE) potentials to calculate the ionospheric electric field along the polar cap boundary, and thus evaluate the variation of the dayside/nightside reconnection rates. We find a strong response of the polar cap boundary at all available MLTs, exhibiting a significant reduction of the open flux content. We also observe an immediate response of the dayside reconnection rate, plus a phased response, delayed by ∼15–20 min, of the nightside reconnection rate. Finally, we provide comparison of the observations with the results of the Open Geospace General Circulation Model (OpenGGCM), elucidating significant agreements and disagreements.

Dynamics of field-aligned currents in two hemispheres during a magnetospheric storm from magnetogram inversion technique data

We continue to study the physical processes occurring during the August 17, 2001 magnetospheric storm by analyzing the dynamics of the intensity of field-aligned currents (FACs) in Iijima—Potemra Region 1 in the polar ionospheres of the two hemispheres, using the modernized magnetogram inversion technique. The results obtained on the dynamics of two types of FAC asymmetry (dawn-dusk and interhemispheric), as well as the previously obtained regularities in the behavior of Hall currents and the polar cap boundaries depending on the large azimuthal component of the interplanetary magnetic field (IMF), observed during the storm, and the seasonal behavior of the conductivity are consistent with the open magnetosphere model and with satellite observations of auroras in two hemispheres. We have shown that the weakening of the asymmetry of two types in the FAC distribution during substorms in the storm under study occurs almost completely in the winter hemisphere and is much weaker in the summer one. We associate this phenomenon with the predominance of the effect of long-term exposure to the azimuthal IMF component in the sunlit polar ionosphere of the summer hemisphere over the substorm symmetrization effect of the night magnetosphere. A symmetrization effect of the polar cap and FACs, created by the solar wind pressure pulse at the end of the storm, is observed. We propose a qualitative explanation of this effect.

Динамика продольных токов в двух полушариях Земли в ходе магнитосферной бури по данным техники инверсии магнитограмм

We continue to study the physical processes occurring during the August 17, 2001 magnetospheric storm by analyzing the dynamics of the intensity of field-aligned currents (FACs) in Iijima—Potemra Region 1 in the polar ionospheres of the two hemispheres, using the modernized magnetogram inversion technique. The results obtained on the dynamics of two types of FAC asymmetry (dawn-dusk and interhemispheric), as well as the previously obtained regularities in the behavior of Hall currents and the polar cap boundaries depending on the large azimuthal component of the interplanetary magnetic field (IMF), observed during the storm, and the seasonal behavior of the conductivity are consistent with the open magnetosphere model and with satellite observations of auroras in two hemispheres. We have shown that the weakening of the asymmetry of two types in the FAC distribution during substorms in the storm under study occurs almost completely in the winter hemisphere and is much weaker in the summer one. We associate this phenomenon with the predominance of the effect of long-term exposure to the azimuthal IMF component in the sunlit polar ionosphere of the summer hemisphere over the substorm symmetrization effect of the night magnetosphere. A symmetrization effect of the polar cap and FACs, created by the solar wind pressure pulse at the end of the storm, is observed. We propose a qualitative explanation of this effect.

Dynamics of the field-aligned currents distribution asymmetry during substorms in the equinox season

We continue to study the physical processes occurring during the August 17, 2001 magnetospheric storm by analyzing the dynamics of the intensity of field-aligned currents (FACs) in Iijima—Potemra Region 1 in the polar ionospheres of two hemispheres, using the modernized magnetogram inversion technique. The results obtained on the dynamics of the FAC asymmetry of two types (dawn–dusk and interhemispheric), as well as the previously obtained regularities in the behavior of Hall currents and polar cap boundaries depending on the large azimuthal component of the interplanetary magnetic field (IMF), observed during the storm, and the seasonal behavior of the conductivity are consistent with the open magnetosphere model and with satellite observations of auroras in two hemispheres. We have shown that the weakening of the asymmetry of two types in the FAC distribution during substorms in the storm under study occurs almost completely in the winter hemisphere and is much weaker in the summer one. We associate this phenomenon with the predominance of the effect of long-term exposure to the azimuthal IMF component in the sunlit polar ionosphere of the summer hemisphere over the substorm symmetrization effect of the night magnetosphere. A symmetrization effect of the polar cap and FACs, created by the solar wind pressure pulse at the end of the storm, is observed. We propose a qualitative explanation of this effect.

Динамика асимметрии распределения продольных токов во время суббурь в сезон равноденствия

We continue to study the physical processes occurring during the August 17, 2001 magnetospheric storm by analyzing the dynamics of the intensity of field-aligned currents (FACs) in Iijima—Potemra Region 1 in the polar ionospheres of two hemispheres, using the modernized magnetogram inversion technique. The results obtained on the dynamics of the FAC asymmetry of two types (dawn–dusk and interhemispheric), as well as the previously obtained regularities in the behavior of Hall currents and polar cap boundaries depending on the large azimuthal component of the interplanetary magnetic field (IMF), observed during the storm, and the seasonal behavior of the conductivity are consistent with the open magnetosphere model and with satellite observations of auroras in two hemispheres. We have shown that the weakening of the asymmetry of two types in the FAC distribution during substorms in the storm under study occurs almost completely in the winter hemisphere and is much weaker in the summer one. We associate this phenomenon with the predominance of the effect of long-term exposure to the azimuthal IMF component in the sunlit polar ionosphere of the summer hemisphere over the substorm symmetrization effect of the night magnetosphere. A symmetrization effect of the polar cap and FACs, created by the solar wind pressure pulse at the end of the storm, is observed. We propose a qualitative explanation of this effect.

Model‐Based Properties of the Dayside Open/Closed Boundary: Is There a UT‐Dependent Variation?

AbstractThe open‐closed boundary (OCB) defines a region of significant transformation in Earth's protective magnetic shield. Principle among these changes is the transition of magnetic field lines from having two foot points, one in each hemisphere, to one foot point at Earth, the other mapping to the solar wind. Charged particles in the solar wind are able to follow these open field lines into Earth's upper atmosphere. The OCB also defines the polar cap boundary. Being able to identify and track the OCB allows study of several components of the geomagnetic system. Among them are the electrodynamics of the geomagnetic field and the reconnection balance between the dayside and nightside of the geomagnetic field. Furthermore, the OCB can provide insights into the precipitation of energetic protons into the ionosphere. Using the Tsyganenko model of the geomagnetic field (T96), we demonstrate a diurnal fluctuation which we call the Universal Time (UT) effect of the OCB. This UT effect is independent of all other inputs. We anticipate this UT effect to have important consequences in modeling the OCB and other polar cap‐associated structures, especially polar cap absorption events that adversely affect high‐frequency radio wave propagation in polar regions.

Диагностика границ аврорального овала на основе техники инверсии магнитограмм

It is shown that the convection reversal boundary is a fundamental parameter of the magnetosphere-ionosphere coupling, which determines a strong analogy between the electrostatic potential of the ionosphere and the equivalent current function in the dipole geomagnetic field approximation and the uniform ionospheric conductance. We have developed a new ground-based method for automatically diagnosing boundaries of the auroral oval using output data obtained with the magnetogram inversion technique (MIT). Using maps of the current function and field-aligned currents, calculated at the first stage of MIT with uniform ionospheric conductance, we determine the convection reversal boundary, polar cap boundary, equatorial boundary of the auroral oval, and line of maximum density of auroral electrojets. These parameters have previously been determined by a visual-manual method: analyzing maps of field-aligned and equivalent currents on the monitor screen and carrying out predetermined boundaries with the mouse — this took a very long time (weeks and months). The comparison between manually and automatically obtained boundaries has shown that the correlation coefficient between the two boundaries is, on average, 0.85, and the root-mean-square deviation does not exceed 2° in latitude. By providing an adequate accuracy for the boundary determination, the automatic method reduces the time for map processing by a factor of 2–3 (to minutes and hours), releasing a researcher from laborious visual work. The new method is implemented as one of the important modules in the updated MIT software.

Diagnostics of auroral oval boundaries on the basis of the magnetogram inversion technique

It is shown that the convection reversal boundary is a fundamental parameter of the magnetosphere-ionosphere coupling, which determines a strong analogy between the electrostatic potential of the ionosphere and the equivalent current function in the dipole geomagnetic field approximation and the uniform ionospheric conductance. We have developed a new ground-based method for automatically diagnosing boundaries of the auroral oval using output data obtained with the magnetogram inversion technique (MIT). Using maps of the current function and field-aligned currents, calculated at the first stage of MIT with uniform ionospheric conductance, we determine the convection reversal boundary, polar cap boundary, equatorial boundary of the auroral oval, and line of maximum density of auroral electrojets. These parameters have previously been determined by a visual-manual method: analyzing maps of field-aligned and equivalent currents on the monitor screen and carrying out predetermined boundaries with the mouse — this took a very long time (weeks and months). The comparison between manually and automatically obtained boundaries has shown that the correlation coefficient between the two boundaries is, on average, 0.85, and the root-mean-square deviation does not exceed 2° in latitude. By providing an adequate accuracy for the boundary determination, the automatic method reduces the time for map processing by a factor of 2–3 (to minutes and hours), releasing a researcher from laborious visual work. The new method is implemented as one of the important modules in the updated MIT software.

Magnetosphere response to the IMF turning from north to south

In this paper, the Space Weather Modeling Framework (SWMF) is used to simulate the real-time response of the magnetosphere to a solar wind event on June 5, 1998, in which the interplanetary magnetic field shifted its direction from north to south. Since most current models do not take into account convective effects of the inner magnetosphere, we first study the importance of Rice Convection Model (RCM) in the global model. We then focus on the following four aspects of the magnetosphere’s response: the magnetosphere’s density distribution, the structure of its magnetic field lines, the area of the polar cap boundary, and the corresponding ionospheric current change. We find that (1) when the IMF changes from north to south in this event, high magnetosheath density is observed to flow downstream along the magnetopause with the solar wind; low-latitude reconnection at dayside occurs under the southward IMF, while the magnetic field lines in the tail lobe caudal, caused by the nightside high latitude reconnection, extend into the interplanetary space. Open magnetic field lines exist simultaneously at both high and low latitudes at the magnetopause; (2) the area of the polar cap is obviously increased if the IMF turns from the north to the south; this observation is highly consistent with empirical observations; (3) the ionospheric field align current in the northern hemisphere is stronger than in the southern hemisphere and also increases as the IMF changes from north to south. SWMF with the Rice Convection effect provides reliable modeling of the magnetospheric and ionospheric response to this solar wind variation.

The Ion/Electron Temperature Characteristics of Polar Cap Classical and Hot Patches and Their Influence on Ion Upflow

AbstractThe term of “polar cap hot patch” is a newly identified high‐density plasma irregularity at high latitudes, which is associated with high electron temperature and particle precipitation, while a classical polar cap patch has lower electron temperature. To investigate characteristics of hot patches versus classical patches, five years of in situ database of plasma observations from the DMSP satellites was analyzed. For the first time, we show how the ion/electron temperature ratio (or temperature difference) can be used to distinguish between classical and hot patches. For classical patches (Ti/Te > 0.8 or Te < Ti + 600 K), the vertical ion flux is generally downward. For hot patches (Ti/Te < 0.8 or Te > Ti + 600 K), the vertical ion flux is generally upward. The highest upflow occurrence was found near the polar cap boundary, associated with hot patches, particle precipitation, strong convection speed, and localized field‐aligned currents. This result shows that the polar cap hot patches may play a very important role in solar wind‐magnetosphere‐ionosphere coupling processes.

Poynting Flux in the Dayside Polar Cap Boundary Regions From DMSP F15 Satellite Measurements

AbstractPoynting flux, which describes electromagnetic energy flux, is an important energy source for the high‐latitude upper atmosphere. After the launch of Defense Meteorological Satellite Program (DMSP) F15 spacecraft with a boom‐mounted magnetometer on board, there was a new opportunity to calculate Earth‐directed Poynting flux at satellite altitudes (~850 km) in the upper atmosphere. A persistent enhancement of thermospheric density in the dayside polar cap boundary regions has been reported in the CHAMP satellite observations. To understand the significance of different physical mechanisms including Poynting flux and particle precipitation, and the correlation between them, a statistical study of Poytning flux and particle energy flux in the dayside cusp and low‐latitude boundary layer (LLBL) regions has been conducted based on DMSP F15 measurements. DMSP satellite observations showed a dominate downward Poynting flux for most cases in the cusp region. Our analysis of DMSP F15 data for five years (2000–2004) reveals that approximately 53% of 660 cusp crossings at 800–850 km showed strong downward Poynting flux (S > 10 mW/m2), 32% of the crossings had noticeable downward Poynting flux (S > 3 mW/m2), and 7% of the crossings did not show clear Poynting flux (S < 1 mW/m2). Only 13 out of 660 cusp crossings (~2%) showed noticeable upward Poynting flux. In the LLBL region, 35% of 11,641 LLBL crossings showed significant downward Poynting flux, 34% of the crossings had noticeable downward Poynting flux, and only 13% of the crossings did not show clear Poynting flux. On average, Poynting flux in LLBL is smaller than that in the cusp. The results show a slightly negative correlation between Poynting flux and particle precipitation energy flux in the dayside polar cap boundary regions. Statistically, Poynting flux in the cusp is enhanced during interplanetary magnetic field By positive conditions.

Swarm Satellite and EISCAT Radar Observations of a Plasma Flow Channel in the Auroral Oval Near Magnetic Midnight

AbstractWe present Swarm satellite and EISCAT radar observations of electrodynamical parameters in the midnight sector at high latitudes. The most striking feature is a plasma flow channel located equatorward of the polar cap boundary within the dawn convection cell. The flow channel is 1.5° wide in latitude and contains southward electric field of 150 mV/m, corresponding to eastward plasma velocities of 3,300 m/s in the F‐region ionosphere. The theoretically computed ion temperature enhancement produced by the observed ion velocity is in accordance with the measured one by the EISCAT radar. The total width of the auroral oval is about 10° in latitude. While the poleward part is electric field dominant with low conductivity and the flow channel, the equatorward part is conductivity dominant with at least five auroral arcs. The main part of the westward electrojet flows in the conductivity dominant part, but it extends to the electric field dominant part. According to Kamide and Kokubun (1996), the whole midnight sector westward electrojet is expected to be conductivity dominant, so the studied event challenges the traditional view. The flow channel is observed after substorm onset. We suggest that the observed flow channel, which is associated with a 13‐kV horizontal potential difference, accommodates increased nightside plasma flows during the substorm expansion phase as a result of reconnection in the near‐Earth magnetotail.

Determination of Polar Cap Boundary for the Substorm Event of 8 March 2008

The polar cap boundary (PCB) is a fundamental indicator of magnetospheric activities especially during a substorm cycle. Taking a period on 8 March 2008 as an example, we investigate the location of PCB and its dynamics during a substorm event. The PCB location is determined from the Piecewise Parabolic Method with a Lagrangian Remap (PPMLR) -Magnetohydrodynamic (MHD) simulation data and Defense Meteorological Satellite Program (DMSP) observations, respectively. Model-observation comparison indicates that the PPMLR-MHD model gives a reliable estimate of PCB location during a complex substorm sequence. We further analyze the evolution of PCB in that period. The polar cap expands under southward interplanetary magnetic field (IMF), since the low-latitude dayside reconnection produces new open magnetic flux. Meanwhile, more solar wind energy enters and stores in the magnetosphere with the decreasing SML (SuperMAG Auroral Lower) index. After the substorm expansion onset, the polar cap contracts for a while due to the explosive increase of nightside reconnection. When the IMF direction turns northward, the polar cap contracts continuously, since the dayside reconnection ceases and no more open magnetic flux are supplied, and the storage energy in the magnetosphere releases with the increasing SML index. The model results are in good accord with the features from observations.