Geospace Environment Research Articles

Citizen science: Development of a low-cost magnetometer system for a coordinated space weather monitoring

As part of Ham Radio Science Citizen Investigation (HamSCI) Personal Space Weather Station (PSWS) project, a low-cost, commercial off-the-shelf magnetometer has been developed to provide quantitative and qualitative measurements of the geospace environment from the ground for both scientific and operational purposes at a cost that will allow for crowd-sourced data contributions. The PSWS magnetometers employ a magneto-inductive sensor technology to record three-axis magnetic field variations with a field resolution of ∼3 nT at a 1 Hz sample rate. The measurement range of the sensor is ±1.1×106 nT) and is valid over a temperature range of −40 °C to +85 °C. Data from the PSWS network will combine these magnetometer measurements with high frequency (HF, 3–30 MHz) radio observations to monitor large-scale current systems and ionospheric disturbances due to drivers from both space and the atmosphere. A densely-spaced magnetometer array, once established, will demonstrate their space weather monitoring capability to an unprecedented spatial extent. Magnetic field data obtained by the magnetometers installed at various locations in the US are presented and compared with the existing magnetometers nearby, demonstrating that the performance is very adequate for scientific investigations.

Volumetric Reconstruction of Ionospheric Electric Currents From Tri‐Static Incoherent Scatter Radar Measurements

AbstractWe present a new technique for the upcoming tri‐static incoherent scatter radar system EISCAT 3D (E3D) to perform a volumetric reconstruction of the 3D ionospheric electric current density vector field, focusing on the feasibility of the E3D system. The input to our volumetric reconstruction technique are estimates of the 3D current density perpendicular to the main magnetic field, j⊥, and its covariance, to be obtained from E3D observations based on two main assumptions: (a) Ions fully magnetized above the E region, set to 200 km here. (b) Electrons fully magnetized above the base of our domain, set to 90 km. In this way, j⊥ estimates are obtained without assumptions about the neutral wind field, allowing it to be subsequently determined. The volumetric reconstruction of the full 3D current density is implemented as vertically coupled horizontal layers represented by Spherical Elementary Current Systems with a built‐in current continuity constraint. We demonstrate that our technique is able to retrieve the three dimensional nature of the currents in our idealized setup, taken from a simulation of an active auroral ionosphere using the Geospace Environment Model of Ion‐Neutral Interactions (GEMINI). The vertical current is typically less constrained than the horizontal, but we outline strategies for improvement by utilizing additional data sources in the inversion. The ability to reconstruct the neutral wind field perpendicular to the magnetic field in the E region is demonstrated to mostly be within ±50 m/s in a limited region above the radar system in our setup.

Infrared Radiation in the Thermosphere From 2002 to 2023

AbstractTwenty‐two years (2002–2023) of infrared radiative cooling rate data derived from the SABER instrument on the NASA TIMED satellite are presented. Global daily and global annual infrared power (Watts, W) emitted by nitric oxide (NO) and carbon dioxide (CO2) illustrate the variability of the geospace environment on timescales from days to decades. The 11‐year solar cycle (SC) is evident in the global power data and in vertical profiles of infrared cooling rates (nW/m3). The global annual power radiated by NO and CO2 are larger in 2023 than at any time since 2003 and 2002, respectively. The to‐date peak in NO infrared power in SC 25 is larger than in SC 24, is comparable to SC 20, but is less than in SCs 18–19 and 21–23. Two geomagnetic storms in 2023 radiated more than 1 TW and are in the top 10 strongest storms observed by SABER.

Validation of Simulated Statistical Characteristics of Magnetosphere‐Ionosphere Coupling in Global Geospace Simulations Over an Entire Carrington Rotation

AbstractWe study the statistical features of magnetosphere‐ionosphere (M‐I) coupling using a two‐way M‐I model, the GT configuration of the Multiscale Atmosphere Geospace Environment (MAGE) model. The M‐I coupling characteristics, such as field‐aligned current, polar cap potential, ionospheric Joule heating, and downward Alfvénic Poynting flux, are binned according to the interplanetary magnetic field clock angles over an entire Carrington Rotation event between 20 March and 16 April 2008. The MAGE model simulates similar distributions of field‐aligned currents compared to empirical Weimer/AMPS models and Iridium observations and reproduces the Region 0 current system. The simulated convection potential agrees well with the Weimer empirical model and displays consistent two‐cell patterns with SuperDARN observations, which benefit from more extensive data sets. The Joule heating structure in MAGE is generally consistent with both empirical Cosgrove and Weimer models. Moreover, our model reproduces Joule heating enhancements in the cusp region, as presented in the Cosgrove model and observations. The distribution of the simulated Alfvénic Poynting flux is consistent with that observed by the FAST satellite in the dispersive Alfvén wave regime. These M‐I coupling characteristics are also binned by the Kp indices, indicating that the Kp dependence of these patterns in the M‐I model is more effective than the empirical models within the Carrington Rotation. Furthermore, the MAGE simulation exhibits an improved M‐I current‐voltage relation that closely resembles the Weimer model, suggesting that the updated global model is significantly improved in terms of M‐I coupling.

Short Wave Infrared Imaging for Auroral Physics and Aeronomy Studies

Advances in solar-terrestrial physics are generally linked to the development of innovative new sensor technologies, affording us ever better sensitivity, higher resolution, and broader spectral response. Recent advances in low-noise InGaAs sensor technology have enabled the realization of low-light-level scientific imaging within the short-wave infrared (SWIR) region of the electromagnetic spectrum. This paper describes a new and highly sensitive ultra-wide angle imager that offers an expansion of auroral and airglow imaging capabilities into the SWIR spectral range of 900–1,700 nm. The imager has already proven successful in large-area remote sensing of mesospheric temperatures and in providing intensity maps showing the propagation and dissipation of atmospheric gravity waves and ripples. The addition of an automated filter wheel expands the range of applications of an already versatile SWIR detector. Several potential applications are proposed herein, with an emphasis on auroral science. The combined data from this type of instrument and other existing instrumentation holds a strong potential to further enhance our understanding of the geospace environment.

Accuracy of Global Geospace Simulations: Influence of Solar Wind Monitor Location and Solar Wind Driving

AbstractSome space weather models, such as the Space Weather Modeling Framework (SWMF) used in this study, use solar wind propagated from the first Lagrange point (L1) to the bow shock nose (BSN) to forecast geomagnetic storms. The SWMF is a highly coupled framework of space weather models that include multiple facets of the Geospace environment, such as the magnetosphere and ionosphere. The propagated solar wind measurements are used as a boundary condition for SWMF. The solar wind propagation method is a timeshift based on the calculated phase front normal (PFN) which leads to some uncertainties. For example, the propagated solar wind could have evolved during this timeshift. We use a data set of 123 geomagnetic storms between 2010 and 2019 run by the SWMF Geospace configuration to analyze the impact solar wind propagation and solar wind driving has on the geomagnetic indices. We look at the probability distributions of errors in SYM‐H, cross polar cap potential (CPCP), and auroral electrojet indices AL and AU. Through studying the median errors (MdE), standard deviations and standardized regression coefficients, we find that the errors depend on the propagation parameters. Among the results, we show that the accuracy of the simulated SYM‐H depends on the spacecraft distance from the Sun‐Earth line. We also quantify the dependence of the standard deviation in SYM‐H errors on the PFN and solar wind pressure. These statistics provide an insight into how the propagation method affects the final product of the simulation, which are the geomagnetic indices.

Studying the Impact of the Geospace Environment on Solar Lithosphere Coupling and Earthquake Activity

In solar–terrestrial physics, there is an open question: does a geomagnetic storm affect earthquakes? We expand research in this direction, analyzing the seismic situation after geomagnetic storms (GMs) accompanied by the precipitation of relativistic electrons from the outer radiation belt to form an additional radiation belt (RB) around lower geomagnetic lines. We consider four widely discussed cases in the literature for long-lived (weeks, months) RBs due to GMs and revealed that the 1/GMs 24 March 1991 with a new RB at L~2.6 was followed by an M7.0 earthquake in Alaska, 30 May 1991, near footprint L = 2.69; the 2/GMs 29 October 2003 (Ap = 204) with new RB first in the slot region at L = 2–2.5 cases followed by an M7.8 earthquake on 17 November 2003 at the Aleutian Islands near footprint L = 2.1, and after forming an RB at L~1.5 which lasted for ~26 months, two mega quakes, M9.1 in 2004 and M8.6 in 2005, occurred at the globe; the 3/GMs 3 September 2012 with a new RB at L= 3.0–3.5 was followed by an M7.8 earthquake in Canada near footprint L = 3.2; and the 4/GMs 21 June 2015 with a new RB at L = 1.5–1.8 was followed by an M6.3 earthquake on 7 September 2015 in New Zealand, near footprint L = 1.58. The obtained results suggest that (1) major earthquakes occur near the footprints of geomagnetic lines filled with relativistic electrons precipitating from the outer radiation belt due to geomagnetic storms, and (2) the time delay between geomagnetic storm onset and earthquake occurrence may vary from several weeks to several months. The results may expand the framework for developing mathematical magnetosphere–ionosphere coupling models.

The Contribution of Plasma Sheet Bubbles to Stormtime Ring Current Buildup and Evolution of Its Energy Composition

AbstractThe formation of the stormtime ring current is a result of the inward transport and energization of plasma sheet ions. Previous studies have demonstrated that a significant fraction of the total inward plasma sheet transport takes place in the form of bursty bulk flows, known theoretically as flux tube entropy‐depleted “bubbles.” However, it remains an open question to what extent bubbles contribute to the buildup of the stormtime ring current. Using the Multiscale Atmosphere Geospace Environment Model, we present a case study of the 17 March 2013 storm, including a quantitative analysis of the contribution of plasma transported by bubbles to the ring current. We show that bubbles are responsible for at least 50% of the plasma energy enhancement within 6 RE during this strong geomagnetic storm. The bubbles that penetrate within 6 RE transport energy primarily in the form of enthalpy flux, followed by Poynting flux and relatively little as bulk kinetic flux. Return flows can transport outwards a significant fraction of the plasma energy being transported by inward flows, and therefore must be considered when quantifying the net contribution of bubbles to the energy buildup. Data‐model comparison with proton intensities observed by the Van Allen Probes show that the model accurately reproduces both the bulk and spectral properties of the stormtime ring current. The evolution of the ring current energy spectra throughout the modeled storm is driven by both inward transport of an evolving plasma sheet population and by charge exchange with Earth's geocorona.

Investigating a Solar Wind Stream Interaction Region using Interplanetary Spacecraft Radio Signals: A Magnetohydrodynamic Simulation Study

Stream interaction regions (SIRs) are spiral heliospheric structures that arise at the interface between fast and preceding slow solar wind regions. SIR enhancements of density and magnetic field intensity, often with magnetic polarity inversion, are potentially geoeffective and therefore important in the analysis of space weather. We studied an MHD heliospheric simulation containing a well-defined SIR using a new instrument concept based on trans-heliospheric radio sensing: Faraday Effect Tracker of Coronal and Heliospheric structures (FETCH). FETCH uses line-of-sight radio propagation techniques to measure Faraday rotation and electron column density. Analysis of the simulated FETCH observations clearly demonstrated density and magnetic field enhancements, and magnetic polarity reversal, all of which were confirmed in Wind spacecraft measurements at 1 au. FETCH provided 4.5–5.7 days lead times for predicting the arrival of SIR features at Earth. The SIR radial speed was estimated to be 350–390 km s−1. These initial results hold promise that FETCH will be valuable in detecting and characterizing the inner heliosphere SIR properties well ahead of their presentation in the local geospace environment.

Effects of Cold Plasma on the Excitation of Internally Driven ULF Waves by Ring Current Ions Based On the Magnetosphere‐Ionosphere Coupled Model

AbstractInternally driven Pc4‐5 waves are excited in the plasmaspheric drainage plume (PDP) and near the plasmapause. The excitation of ultralow frequency (ULF) waves was investigated by using the magnetosphere‐ionosphere coupled model between Geospace Environment Modeling System for Integrated Studies‐Ring Current (GEMSIS‐RC) and GEMSIS‐POTential solver (GEMSIS‐POT). In order to investigate the effects of cold plasma on the wave excitation, the simulation code to describe the dynamics of cold plasma was included in the model. The model can reproduce the shrink of the plasmasphere on the nightside and the formation of the PDP on the dayside. First, fundamental Pc5 waves are excited through the drift resonance on the dayside. The waves are caused by positive energy gradient of ion phase space density (PSD) at 50–130 keV. Second harmonic waves (drift‐bounce resonance) are generated outside the plasmapause. These two types of ULF waves are also seen in the case of constant density. Unlike the case of constant density, localized eastward propagating Pc4 waves (drift resonance) are seen on the dawnside associated with the azimuthal density gradient. The azimuthal wave number of Pc4 waves is about 70 and anti‐earthward gradient of PSD about 10 keV contributes to wave growth. We also detect fundamental Pc4‐5 waves near the plasmapause on the nightside in the drift resonance with 100–150 keV ions. Simulation results suggest that the plasmapause has an effect to sustain the excitation of Pc4‐5 ULF wave through the drift resonance.

A Time-efficient, Data-driven Modeling Approach for Predicting the Geomagnetic Impact of Coronal Mass Ejections

To understand the global-scale physical processes behind coronal mass ejection (CME)–driven geomagnetic storms and predict their intensity as a space weather forecasting measure, we develop an interplanetary CME flux rope–magnetosphere interaction module using 3D magnetohydrodynamics. The simulations adequately describe CME-forced dynamics of the magnetosphere including the imposed magnetotail torsion. These interactions also result in induced currents, which are used to calculate the geomagnetic perturbation. Through a suitable calibration, we estimate a proxy of geoeffectiveness—the Storm Intensity index (STORMI)—that compares well with the Dst/SYM-H index. Simulated impacts of two contrasting CMEs quantified by the STORMI index exhibit a high linear correlation with the corresponding Dst and SYM-H indices. Our approach is relatively simple, has fewer parameters to be fine-tuned, and is time efficient compared to complex fluid-kinetic methods. Furthermore, we demonstrate that flux rope erosion does not significantly affect our results. Thus our method has the potential to significantly extend the time window for predictability—an outstanding challenge in geospace environment forecasting—if early predictions of near-Earth CME flux rope structures based on near-Sun observations are available as inputs. This study paves the way for early warnings based on operational predictions of CME-driven geomagnetic storms.

Causality in heliophysics: Magnetic fields as a bridge between the Sun’s interior and the Earth’s space environment

Our host star, the Sun, is a middle-aged main sequence G type star whose activity varies. These variations are primarily governed by solar magnetic fields which are produced in the Sun’s interior via a magnetohydrodynamic dynamo mechanism. Solar activity manifests across different timescales, spanning transient phenomena such as flares, energetic particle events and coronal mass ejections to short to long-term modulation of solar irradiance, plasma winds, open flux and cosmic ray flux in the heliosphere. Collectively, these phenomena define space weather and space climate, which impact the state of the near-Earth space environment, the Earth’s magnetosphere, atmosphere and our space-reliant technologies. Understanding physical processes that are at the heart of solar variability and which causally connect the Sun–Earth system is therefore of immense importance to humanity. Such understanding leads to predictions of the impact of solar activity on our planet and provides a window to explore the plasma universe and other star–planet systems, including assessing the habitability of (exo)planets. In this review, based on our research on the solar–terrestrial system and extant scientific literature, we illuminate processes related to the genesis of solar magnetic fields in the Sun’s interior, their emergence and evolution, their manifestation as solar eruptive events, and their eventual impact on the geospace environment mediated via solar winds and storms. We focus on few phenomena that establish causal connections and demonstrate how our current understanding can lead to development of predictive capabilities encompassing the domain of heliophysics.

Mesoscale phenomena and their contribution to the global response: a focus on the magnetotail transition region and magnetosphere-ionosphere coupling

An important question that is being increasingly studied across subdisciplines of Heliophysics is “how do mesoscale phenomena contribute to the global response of the system?” This review paper focuses on this question within two specific but interlinked regions in Near-Earth space: the magnetotail’s transition region to the inner magnetosphere and the ionosphere. There is a concerted effort within the Geospace Environment Modeling (GEM) community to understand the degree to which mesoscale transport in the magnetotail contributes to the global dynamics of magnetic flux transport and dipolarization, particle transport and injections contributing to the storm-time ring current development, and the substorm current wedge. Because the magnetosphere-ionosphere is a tightly coupled system, it is also important to understand how mesoscale transport in the magnetotail impacts auroral precipitation and the global ionospheric system response. Groups within the Coupling, Energetics and Dynamics of Atmospheric Regions Program (CEDAR) community have also been studying how the ionosphere-thermosphere responds to these mesoscale drivers. These specific open questions are part of a larger need to better characterize and quantify mesoscale “messengers” or “conduits” of information—magnetic flux, particle flux, current, and energy—which are key to understanding the global system. After reviewing recent progress and open questions, we suggest datasets that, if developed in the future, will help answer these questions.

Investigation of the Southern Hemisphere Mid‐High Latitude Thermospheric ∑O/N2 Responses to the Space‐X Storm

AbstractThe geomagnetic storm on February 3, 2022 caused the loss of 38 Starlink satellites of Space‐X. The Global‐scale Observations of the Limb and Disk (GOLD) observations and Multi‐Scale Atmosphere Geospace Environment (MAGE) model simulations are utilized to investigate the thermospheric composition responses to the Space‐X storm. The percentage difference of the GOLD observed thermospheric O and N2 column density ratio (∑O/N2) between the storm time (February 3, Day‐of‐Year [DOY] 34) and quiet time (DOY 32) shows a depletion region in the local noon sector mid‐high latitudes in the southern hemisphere, which corresponds to the east side of GOLD field‐of‐view (FOV). This is different from the classic theory of thermospheric composition disturbance during geomagnetic storms, under which the ∑O/N2 depletion is usually generated at local midnight and high latitudes, and thus, appear on the west side of GOLD FOV. MAGE simulations reproduce the observations qualitatively and indicate that the ∑O/N2 depletion is formed due to strong upwelling in the local morning caused by strong Joule heating. Interestingly, enhanced equatorward winds appear near local midnight, but also in the local morning sector, which transports ∑O/N2 depletion equatorward. The depletion corotates toward the local afternoon and is observed in the GOLD FOV. The equatorward winds in the local morning are due to the ion‐neutral coupling under the conditions of a dominant positive interplanetary magnetic field east‐west component (By) during the storm.

An Ionospheric Conductance Gradient Driver for Subauroral Picket Fence Visible Signatures Near STEVE Events

AbstractReported observations of picket fence signatures associated with subauroral Strong Thermal Emission Velocity Enhancement (STEVE) emission events look strikingly similar to rayed auroral curtains. Rayed auroral curtains are often the visible signatures of tearing‐mode‐unstable current sheets with precipitating auroral electron current carriers. Picket fence signatures are not located where auroral precipitation explanations apply, and are closely collocated with STEVE emission events. A similar tearing‐mode‐instability explanation can be invoked with a different source for the originating field‐aligned current (FAC) sheet. In this explanation, the FAC sheet is sourced by ionospheric conductance gradients adjacent to the localized flows of the STEVE event. Geospace Environment Model of Ion‐Neutral Interactions (GEMINI) models of the 3D ionosphere near counterstreaming STEVE‐associated flow structures show the development of sufficiently strong current sheets for tearing mode instabilities to take hold. These instabilities can locally accelerate ambient ionospheric thermal electrons to the few eV needed for the reported observed green picket fence signatures.

Geospace: The naturally occurring plasma laboratory

While the societal relevance of space physics drives most of the research in the field, this paper argues for the merit in addressing problems that are core to plasma physics, regardless of discipline or application to society. The geospace environment is a unique, naturally occurring laboratory that can be used for investigating basic plasma processes. Ground and satellite-based observations provide important data for understanding basic plasma processes such as transport, charge-neutral interactions, collisions, and turbulence. These are important topics in space physics, and we show how these topics are relevant to the broader plasma physics community. This paper recommends improved funding lines to examine fundamental plasma physics in the geospace environment. This can be accomplished by annual solicitation of the NASA ROSES B.3 Heliophysics Theory, Modeling, and Simulations (HTMS) program; consideration of plasma physics as relevant to all NASA proposals; and creating a pathway for high-risk research into fundamental science.

Python tools for ESA’s Swarm mission: VirES for Swarm and surrounding ecosystem

ESA’s Swarm mission is a constellation probing both Earth’s interior and geospace, delivering magnetic and plasma measurements which are used to generate many derived data products. From empirical magnetic field models of the core, crust, ionosphere, and magnetosphere, to multi-point estimates of ionospheric currents and in-situ plasma properties, these are challenging to navigate, process, and visualize. The VirES for Swarm platform (https://vires.services) has been built to tackle this problem, providing tools to increase usability of Swarm data products. The VirES (Virtual environments for Earth Scientists) platform provides both a graphical web interface and an API to access and visualise Swarm data and models. This is extended with a cloud-hosted development environment powered by JupyterHub (the “Virtual Research Environment/VRE”). VirES provides two API’s: the full VirES API for which a dedicated Python client is provided, viresclient, and the more interoperable Heliophysics API (HAPI). The VRE is furnished with a bespoke Python environment containing thematic libraries supporting science with Swarm. This service aims to ease the pathway for scientists writing computer code to analyze Swarm data products, increase opportunities for collaboration, and leverage cloud technologies. Beyond simply providing data and model access to Python users, it is extremely helpful to provide higher-level analysis and visualization tools, and ready-to-use code recipes that people can explore and extend. Critically for space physics, this involves crossover with many other datasets and so it is highly valuable to embed such tools within the wider data and software ecosystems. Through Swarm DISC (Data, Innovation, and Science Cluster), we are tackling this through cookbooks and Python libraries. Cookbooks are built and presented using Jupyter technologies, and tested to work within the VRE. A new library we are building is SwarmPAL, which includes tools for time-frequency analysis and inversion of magnetic field measurements for electric current systems, among others, while relying on the VirES server to provide data portability and other utilities. This paper reviews the current state of these tools and services for Swarm, particularly in the context of the Python in Heliophysics Community, and the wider heliophysics and geospace data environment.

Thermospheric Wind Observation and Simulation during the Nov 4, 2021 Geomagnetic Storm Event

Thermospheric wind observations from high to mid latitudes are compared with the newly developed Multiscale Atmosphere Geospace Environment (MAGE) model for the Nov 3–4 geomagnetic storm. The observation and simulation comparison shows a very good agreement and is better at high latitudes in general. We were able to identify a thermospheric poleward wind reduction possibly linked to a northward turning of the Interplanetary Magnetic Field (IMF) at ~22 UT on Nov 3 and an enhancement of the poleward wind to a southward turning near 10 UT on Nov 4 at high latitudes. An IMF southward turning may have led to an enhancement of equatorward winds at Boulder, Colorado near midnight. Simultaneous occurrence of aurora may be associated with an IMF By turning negative. The MAGE model wind simulations are consistent with observations in these cases. The results show the model can be a very useful tool to further study the magnetosphere and ionosphere coupling on short time scales.

Spatiotemporal Limitations of Data‐Driven Modeling: An ISINGLASS Case Study

AbstractIonospheric outflow remains difficult to predict because of myriad processes acting over a large range of altitudes, physical regimes, and time scales. This study uses the ionospheric model Geospace Environment Model of Ion‐Neutral Interactions with Transverse Ion Acceleration (GEMINI‐TIA) and sounding rocket campaign data to assess the impact of this spatially and temporally variable energy inputs on ionospheric upflow. The Ionospheric Structuring: In situ and Groundbased Low Altitude Studies B sounding rocket was launched on 03/02/17, into a substorm expansion, and included nearby ground‐based measurements which, along with in situ data, are used to construct electron precipitation and electric field inputs to GEMINI‐TIA for purposes of modeling complicated upflow patterns in two important limits: (a) time variable but spatially constant inputs or (b) spatially variable inputs that are constant in time. Realistic ionospheric responses will likely fall between these two limiting cases and may also include additional energy sources not constrained by the rocket. To investigate mixed spatial and temporal forcing we conduct a couple of simulations driven with ground‐based data that resolves both space and time variations in the forcing but at a much lower spatial and temporal resolution. Results demonstrate that the time history of the ionosphere also plays a significant role when accurately describing the upflow response to auroral heating. It takes up to 20 min for simulations initialized with different conditions to converge, suggesting that at the local scales and altitude range studied here it is necessary to include at least that much hysteresis in data analysis and modeling efforts.

Excitation of Two Types of Storm‐Time Pc5 ULF Waves by Ring Current Ions Based on the Magnetosphere‐Ionosphere Coupled Model

AbstractTwo types of storm‐time Pc5 waves are excited in the magnetosphere–ionosphere coupled model. We conduct magnetosphere–ionosphere (M–I) coupling between geospace environment modeling system for integrated studies‐ring current (GEMSIS‐RC) and GEMSIS‐POTential solver (GEMSIS‐POT) models in order to investigate the excitation of storm‐time Pc5 waves under more realistic conditions. The coupled model enables us to simulate ion transport from the plasma sheet by the convection electric field. First, we find the excitation of fundamental mode Pc5 waves on the dayside. These waves are generated by the drift resonance and driven by positive energy gradient of ion phase space density (PSD). The wave excitation is similar with the case without M–I coupling. Second, we detect second harmonic mode Pc5 waves with high azimuthal wave number in the dusk and premidnight region. The PSD oscillates associated with the waves, and the drift‐bounce resonance occurs at 50–80 keV ions. We find that these ions drift around the Earth and reach the nightside again at t ∼ 4,500 s, whose time scale is comparable to the drift period of these ions. At that time, inward gradient of the PSD is created, which contributes to the wave growth. This effect does not happen in the case without dawn‐to‐dusk convection electric field driven by Region 1 field aligned current (FAC). We find that the shielding by Region 2 FAC suppresses the flux of newly injected ions at the nightside, inward gradient of the PSD, and power spectra of second harmonic mode waves.