SuperDARN cross polar cap potential dependence on the solar wind conditions and comparisons with models

In this study SuperDARN Cross Polar Cap Potentials (CPCPs), collected over the year 2000, are investigated with a goal to statistically assess its relationship with various parameters of the solar wind and Interplanetary Magnetic Field (IMF). We show that SuperDARN CPCPs tend to cluster around discrete values, prescribed by the statistical model, unless the amount of points on each convection map is above ∼300. By selecting CPCP data obtained with radar coverage of >300 points, we investigate the CPCP relationship with IMF B_z and B_y, IMF clock angle, solar wind speed and dynamic pressure, Alfven velocity, Alfven–Mach number, and interplanetary electric field. Some reported tendencies, such as dependence upon IMF B_z, were found to be consistent with measurements by other instruments. We demonstrate that SuperDARN CPCPs show consistency with several theories/empirical models (predicting the CPCP) in terms of a linear trend but, on average, the slopes of the dependencies are at least two times smaller. We also determine the coupling function, out of those published in literature, best correlating with SuperDARN CPCPs.     

地球磁场对太阳风扰动响应的研究

分别对行星际激波、太阳风动压增大事件和减小事件的地球磁场响应进行了比较. 分析结果表明, 同步轨道磁场对太阳风扰动在向阳面产生较强的正响应, 在背阳面 响应较弱且有时会出现负响应, 地磁指数SYM-H对太阳风扰动的响应为正响应. 同时还得出, 向阳侧同步轨道磁场响应幅度d Bz与地磁指数响应幅度d SYM-H、上下游动压均方差均具有较好的相关性. 地磁指数响应幅度与同步轨道磁场响应幅度相关关系在激波和动压增大事件中具有一致性, 动压减小事件出 现明显差异, 这说明激波和动压增大事件在影响地球磁场方面具有某种共性.      

A solar wind driven model of geosynchronous plasma moments

We describe a tabular specification model of the density and temperature of ions and electrons at geosynchronous orbit as a function of magnetic local time and solar wind parameters. This model can be used to provide boundary conditions for numerical ring current models. Unlike previous specification models of geosynchronous plasma moments, this model is parameterized by upstream solar wind conditions. We find that solar wind parameters are a better predictor of geosynchronous ion density than magnetospheric indices, and as upstream parameters they are often more appropriate as model inputs since they causally precede the model outputs. Of the upstream parameters that were tested, the best predictors of geosynchronous conditions were the solar wind flow pressure and the magnitude and Z-component of the interplanetary magnetic field.     

Magnetospheric disturbances associated with the 13 December 2006 solar flare and their ionospheric effects over North-East Asia

We present an observational study of magnetospheric and ionospheric disturbances during the December 2006 intense magnetic storm associated with the 4В/Х3.4 class solar flare. To perform the study we utilize the ground data from North–East Asian ionospheric and magnetic observatories (60–72°N, 88–152°E) and in situ measurements from LANL, GOES, Geotail and ACE satellites. The comparative analysis of ionospheric, magnetospheric and heliospheric disturbances shows that the interaction of the magnetosphere with heavily compressed solar wind and interplanetary magnetic field caused the initial phase of the magnetic storm. It was accompanied by the intense sporadic E and F2 layers and the total black-out in the nocturnal subauroral ionosphere. During the storm main phase, LANL-97A, LANL 1994_084, LANL 1989-046 and GOES_11 satellites registered a compression of the dayside magnetosphere up to their orbits. In the morning–noon sector the compression was accompanied by an absence of reflections from ionosphere over subauroral ionospheric station Zhigansk (66.8°N, 123.3°E), and a drastic decrease in the F2 layer critical frequency (foF2) up to 54% of the quite one over subauroral Yakutsk station (62°N, 129.7°E). At the end of the main phase, these stations registered a sharp foF2 increase in the afternoon sector. At Yakutsk the peak foF2 was 1.9 time higher than the undisturbed one. The mentioned ionospheric disturbances occurred simultaneously with changes in the temperature, density and temperature anisotropy of particles at geosynchronous orbit, registered by the LANL-97A satellite nearby the meridian of ionospheric and magnetic measurements. The whole complex of disturbances may be caused by radial displacement of the main magnetospheric domains (magnetopause, cusp/cleft, plasma sheet) with respect to the observation points, caused by changes in the solar wind dynamic pressure, the field of magnetospheric convection, and rotation of the Earth.     

A working hypothesis for connections between electrically-induced changes in cloud microphysics and storm vorticity,with possible effects on circulation

This paper outlines, and explores the uncertainties in, hypothesized connections between a series of processes that could explain two long-standing puzzles; those of (1) the observed winter storm vorticity responses to atmospheric energy inputs that change the ionosphere–earth current density, J_z, that appear to involve storm invigoration, and (2) changes in anti-cyclonic blocking and circulation that include the observed colder winters in Great Britain and western Europe at solar minima, and especially at extended solar minima. A working hypothesis for the mechanism responsible for (1) is that the flow of J_z through conductivity gradients, as in stratified cloud layers and fog, especially with sea-salt aerosol haze over the high latitude winter oceans, deposits electric changes on droplets and aerosol particles; most importantly on cloud condensation nuclei (CCN). These electric charges modulate scavenging of the particles in clouds and haze layers, increasing the concentration of small CCN and decreasing the concentration of large CCN. When further cloud formation occurs there is increased concentration of small droplets and decreased concentration of large ones, reducing coalescence and the production of rain. Thus updrafts carry more liquid water above the freezing level, and there the increased production of ice releases more latent heat and invigorates the updraft (the Rosenfeld mechanism), leading to increased vorticity. Here we explore the major uncertainties for the reality of the above chain of physical processes. A consequence of cumulative cyclonic vorticity increases is increases in downstream anti-cyclonic blocking. A further working hypothesis for (2) is that the invigoration may be large enough to contribute to the observed increases in blocking in winters at solar minima (high J_z) in the North Atlantic, that result in colder winters in the UK and northern Europe.     

The role of volcanic aerosols and relativistic electrons in modulating winter storm vorticity

Small changes in the vorticity of winter storms, responding to solar wind variations, are found in winters from 1957 to 2011, and are greater for winters with higher levels of stratospheric volcanic aerosols. Using 1993–2011 data, the response of the vorticity area index (VAI) is shown to be of larger amplitude when the days of minima in the relativistic electron flux (REF) precipitating from the radiation belts are used, instead of heliospheric current sheet (HCS) crossings, as key days in superposed epoch analyses. The HCS crossings mostly occur within a few days of the REF minima. The VAI is an objective measure of the area of high cyclonic vorticity, and for the present work is derived from ERA-40 and ERA-Interim reanalyses of global meteorological data. The VAI dependencies on the stratospheric aerosol content (SAC) and the REF are consistent with a model in which the ionosphere-earth current density (J_z) affects cloud microphysics. One of the ways in which J_z is modulated is by changes in stratospheric column resistance (S), which is increased by stratospheric aerosols. Because S is in series with the tropospheric column resistance (T), J_z modulation by REF requires that S be not negligible with respect to T. So the J_z modulation and the VAI response appear when the SAC is very high, or the REF reductions (which also increase S) are very deep, and when the product of the SAC and the reciprocal of the REF exceeds a threshold value dependent on T.     

Comparison of the magnetic field before the subsolar magnetopause with the magnetic field in the solar wind before the bow shock

Crossings of the magnetopause near the subsolar point are analyzed using data of THEMIS mission. Variations of the magnetic field near magnetopause measured by one of THEMIS satellites are studied and compared with simultaneous measurements in the solar wind by another THEMIS satellite. The time delay of the solar wind arrival at the subsolar point of the magnetopause is taken into account. 30 and 90 s averaging of the magnetic field in the magnetosheath is produced. The results of averaging are compared with the results of measurements in the solar wind before the bow shock and foreshock. It is shown, that _Bx$B_x$ component of the magnetic field near magnetopause is near to zero, which supports the possibility to consider the magnetopause as the tangential discontinuity. Comparatively good correlation of _By$B_y$ component in the solar wind and near the magnetopause is observed. The correlation of _Bz$B_z$ component near the magnetopause and IMF is practically absent, the sign of the _Bz$B_z$ near the subsolar point does not coincide with the sign of IMF _Bz$B_z$ in ∼$\sim$30% cases.     

Contribution of geometry of interaction between interplanetary and terrestrial magnetic fields into global magnetospheric state and geomagnetic activity

The paper presents results of our study of dependence of geomagnetic activity from geoeffective parameters taking into account mutual orientation of the interplanetary magnetic field, electric field of the solar wind and geomagnetic moment. We attract a reconnection model elaborated by us made allowance for changes of geometry of the solar wind–magnetosphere interaction during annual and diurnal motions of the Earth. We take as our data base the interplanetary magnetic field and solar wind velocity measured at 1 a.u. at ecliptic plane for the period of 1963–2005 and K_p, D_st, a_m indices. Taken as a whole a geoeffective parameter suggested by us explains 95% of observed variations of the indices. Changes of the geometric factor determined by mutual orientation of the solar wind electric field and geomagnetic moment explain larger than 75% of observed statistical variations of D_st and a_m indices. Based on our results we suggest a new explanation of semi-annual and UT variation of geomagnetic activity.     

Moderate geomagnetic storms of January 22–25, 2012 and their influences on the wave components in ionosphere and upper stratosphere-mesosphere regions

Moderate geomagnetic storms occurred during January 22–25, 2012 period. The geomagnetic storms are characterized by different indices and parameters. The SYM-H value on January 22 increased abruptly to 67 nT at sudden storm commencement (SSC), followed by a sharp decrease to −87 nT. A second SSC on January 24 followed by a shock on January 25 was also observed. These SSCs before the main storms and the short recovery periods imply the geomagnetic storms are CME  -driven. The sudden jump of solar wind dynamic pressure and IMF _Bz$B_z$ are also consistent with occurrence of CMEs. This is also reflected in the change in total electron content (TEC) during the storm relative to quiet days globally. The response of the ionospheric to geomagnetic storms can also be detected from wave components that account for the majority of TEC variance during the period. The dominant coherent modes of TEC variability are diurnal and semidiurnal signals which account upto 83% and 30% of the total TEC variance over fairly exclusive ionospheric regions respectively. Comparison of TEC anomalies attributed to diurnal (DW1) and semidiurnal (SW2) tides, as well as stationary planetary waves (SPW1) at 12 UTC shows enhancement in the positive anomalies following the storm. Moreover, the impact of the geomagnetic storms are distinctly marked in the daily time series of amplitudes of DW1, SW2 and SPW1. The abrupt changes in amplitudes of DW1 (5 TECU) and SW2 (2 TECU) are observed within 20°S–20°N latitude band and along 20°N respectively while that of SPW1 is about 3 TECU. Coherent oscillation with a period of 2.4 days between interplanetary magnetic field and TEC was detected during the storm. This oscillation is also detected in the amplitudes of DW1 over EIA regions in both hemispheres. Eventhough upward coupling of quasi two day wave (QTDWs) of the same periodicity, known to have caused such oscillation, are detected in both ionosphere and upper stratosphere, this one can likely be attributed to the geomagnetic storm as it happens after the storm commencement. Moreover, further analysis has indicated that QTDWs in the ionosphere are strengthened as a result of coherent oscillation of interplanetary magnetic field with the same frequency as QTDWs. On the otherhand, occurrences of minor SSW and geomagnetic storms in quick succession complicated clear demarcation of attribution of the respective events to variability of QTDWs amplitudes over upper stratosphere.     

Helio-geomagnetic influence in cardiological cases

The effects of the energetic phenomena of the Sun, flares and coronal mass ejections (CMEs) on the Earth’s ionosphere–magnetosphere, through the solar wind, are the sources of the geomagnetic disturbances and storms collectively known as Space Weather. The research on the influence of Space Weather on biological and physiological systems is open. In this work we study the Space Weather   impact on Acute Coronary Syndromes (ACS) distinguishing between ST-segment elevation acute coronary syndromes (STE–ACS) and non-ST-segment elevation acute coronary syndromes (NSTE–ACS) cases. We compare detailed patient records from the 2nd Cardiologic Department of the General Hospital of Nicaea (Piraeus, Greece) with characteristics of geomagnetic storms (_DST${\mathrm{D}}_{\mathit{ST}}$), solar wind speed and statistics of flares and CMEs which cover the entire solar cycle 23 (1997–2007). Our results indicate a relationship of ACS to helio-geomagnetic activity as the maximum of the ACS cases follows closely the maximum of the solar cycle. Furthermore, within very active periods, the ratio NSTE–ACS to STE–ACS, which is almost constant during periods of low to medium activity, changes favouring the NSTE–ACS. Most of the ACS cases exhibit a high degree of association with the recovery phase of the geomagnetic storms; a smaller, yet significant, part was found associated with periods of fast solar wind without a storm.     

Probing the method of correcting Faraday rotation in vector magnetograms

By analyzing the vector magnetograms of Huairou Solar Observing Station (HSOS) taken at the line center (0.0 Å) and the line wing (−0.12 Å) of FeI λ5324.19 Å, we make an estimate of the measured errors in transversal azimuths (δ?) caused by Faraday rotation. Since many factors, such as the magnetic saturation and scattered light, can affect the measurement accuracy of the longitudinal magnetic field in the umbrae of sunspots, we limit our study in the region ∣B_z∣ < 800 G. The main mean azimuth rotations are about 4°, 6°, 7° and 9°, while ∣B_z∣ are in the ranges of 400–500 G, 500–600 G, 600–700 G and 700–800 G, respectively. Moreover, we find there is also an azimuth rotation of about 8° at the wavelength offset −0.12 Å of the line compared against a previous numerical simulation.     

27-day variations of the galactic cosmic ray intensity and anisotropy in the minimum of the 23rd solar activity cycle

We study the 27-day variations of the solar wind velocity, galactic cosmic ray (GCR) intensity and anisotropy in the last minimum epoch of solar activity (2007–2009, A < 0). The average amplitude of the 27-day variation of the galactic cosmic ray anisotropy (A27A) in the current minimum epoch of solar activity (2007–2009, A < 0) is lesser than in previous positive polarity period as it is expected from the drift theory. So, polarity dependence rule for the 27-day variation of the GCR anisotropy is fully kept. It is a universal principle for the amplitudes of the 27-day variation of the GCR anisotropy. At the same time, the average amplitude of the 27-day variation of the GCR intensity (A27I) remains at the same level as for previous minimum epoch 1995–1997 (A > 0) showing by the same token an violation of its polarity dependence rule established earlier. We assume that this phenomenon could be generally related with the well established 27-day variation of the solar wind velocity being in anti-correlation with the similar changes of the 27-day variation of the GCR intensity. Generally, a character of the heliolongitudinal asymmetry of spatial large-scale structure of the solar wind velocity (SWV) established in the recent minimum epoch, preferentially pronounces in the behavior of the 27-day variation of the GCR intensity than anisotropy. The formation of the 27-day variation of the GCR anisotropy preferentially takes place in a restricted disk like local vicinity in the helioequatorial region, whilst the 27-day variation of the GCR intensity is formed in the global three dimensional vicinity of the heliosphere.     

Prediction of geosynchronous electron fluxes using an artificial neural network driven by solar wind parameters

There are hundreds of satellites operating at the geosynchronous (GEO) orbit where relativistic electrons can cause severe damage. Thus, predicting relativistic electron fluxes is significant for spacecraft safety. In this study, using GOES satellite data during 2011–2020, we propose two neural network models with two hidden layers to predict geosynchronous relativistic electron fluxes at two energy channels (>0.8 MeV and > 2 MeV). The number of input neurons of the two channels (>0.8 MeV and > 2 MeV) are determined to be 36 and 44, respectively. The > 0.8 MeV model has 22 and 9 neurons in the hidden layers, while the > 2 MeV model has 25 and 15 neurons in the hidden layers. The input parameters include the north–south component of the interplanetary magnetic field, solar wind speed, solar wind dynamic pressure and solar wind proton density. Through the analysis of different time delays, we determine that the optimal time delays of two energy channels (>0.8 MeV and > 2 MeV) are 8 days and 10 days, respectively. The training set and validation set (Jan 2011-Dec 2018) are divided by the 10-fold cross-validation method, and the remaining data (Jan 2019-Feb 2020) is used to analyze the model performance as a test set. The prediction results of both energy channels show good agreement with satellite observations indicated by low RMSE (～0.3 cm^-2sr^-1s^?1), high PE (～0.8) and CC (～0.9). These results suggest that only using solar wind parameters is capable of obtaining reasonable predictions of geosynchronous relativistic electron fluxes.     

On the time lag between solar wind dynamic parameters and solar activity UV proxies

The solar activity displays variability and periodic behaviours over a wide range of timescales, with the presence of a most prominent cycle with a mean length of 11 years. Such variability is transported within the heliosphere by solar wind, radiation and other processes, affecting the properties of the interplanetary medium. The presence of solar activity–related periodicities is well visible in different solar wind and geomagnetic indices, although their time lags with respect to the solar cycle lead to hysteresis cycles. Here, we investigate the time lag behaviour between a physical proxy of the solar activity, the Ca II K index, and two solar wind parameters (speed and dynamic pressure), studying how their pairwise relative lags vary over almost five solar cycles. We find that the lag between Ca II K index and solar wind speed is not constant over the whole time interval investigated, with values ranging from 6 years to $～$1 year (average 3.2 years). A similar behaviour is found also for the solar wind dynamic pressure. Then, by using a Lomb-Scargle periodogram analysis we obtain a 10.21-year mean periodicity for the speed and 10.30-year for the dynamic pressure. We speculate that the different periodicities of the solar wind parameters with respect to the solar 11-year cycle may be related to the overall observed temporal evolution of the time lags. Finally, by accounting for them, we obtain empirical relations that link the amplitude of the Ca II K index to the two solar wind parameters.     

Dynamics of the convection in the inner magnetosphere by observations of the diffuse aurora,SAR arc and ionospheric drift

The first results of the comparison of subauroral luminosity dynamics in 557,7 and 630,0 nm emission with simultaneous measurements of the ionospheric drift in the F2 region with a digisonde DPS-4 at the Yakutsk meridian (CGMC: 55–60N, 200°E) at Kp = 2–6 are presented. It is shown from the analysis of individual events that during the magnetospheric convection intensification after the turn of the IMF B_z – component to the south the equatorward extension of diffuse aurora takes place. At the same time the westward ionospheric drift velocity increases both in the diffuse aurora region and much equatorward of it due to the occurrence of the northward polarization electric field. We suppose that the generation of polarization field can be associated with the development of the region 2 FAC during the intensification of magnetospheric convection. The comparison of ground-based observations with measurements of the plasma drift aboard the DMSP-F15 satellite has been carried out.     

On the possibility to check the magnetosphere’s model by CR: The strong geomagnetic storm in November 2003

We report on a study of cosmic ray cutoff rigidity variations during the strong geomagnetic storm of 18–24 November 2003. We employed the most recent Tsyganenko magnetospheric model to represent the very strong disturbed Magnetosphere. We used this magnetic field for the cosmic ray trajectory calculations to determine the geomagnetic cutoff rigidity throughout this period of severe geomagnetic disturbance. We determine the cutoff rigidity changes during this period by two methods, by trajectory calculations and by the spectrographic global survey method. The values of the change of cutoff rigidities obtained by two different methods are correlated with the D_st and interplanetary magnetic field and plasma parameters and result in correlation coefficients in the range 0.63–0.84 for the various cosmic ray stations. The result of this study indicates that the most significant contributions to the cutoff rigidity changes are due to D_st variation although the influence of solar wind density and B_z and B_y components of IMF variations is significant.     

Modeling and observations of the north–south ionospheric asymmetry at low latitudes at long deep solar minimum

Using the physics based model SUPIM and FORMOSAT-3/COSMIC electron density data measured at the long deep solar minimum (2008–2010) we investigate the longitude variations of the north–south asymmetry of the ionosphere at low latitudes (±30° magnetic). The data at around diurnal maximum (12:30–13:30 LT) for magnetically quiet (Ap ? 15) equinoctial conditions (March–April and September–October) are presented for three longitude sectors (a) 60°E–120°E, (b) 60°W–120°W and (c) 15°W–75°W. The sectors (a) and (b) have large displacements of the geomagnetic equator from geographic equator but in opposite hemispheres with small magnetic declination angles; and sector (c) has large declination angle with small displacement of the equators; vertical E × B drift velocities also have differences in the three longitude sectors. SUPIM investigates the importance of the displacement of the equators, magnetic declination angle, and E × B drift on the north–south asymmetry. The data and model qualitatively agree; and indicate that depending on longitudes both the displacement of the equators and declination angle are important in producing the north–south asymmetry though the displacement of the equators seems most effective. This seems to be because it is the displacement of the equators more than the declination angle that produces large north–south difference in the effective magnetic meridional neutral wind velocity, which is the main cause of the ionospheric asymmetry. For the strong control of the neutral wind, east–west electric field has only a small effect on the longitude variation of the ionospheric asymmetry. Though the study is for the long deep solar minimum the conclusions seem valid for all levels of solar activity since the displacement of the equators and declination angle are independent of solar activity.     

Dependence of the 27-day variation of cosmic rays on the global magnetic field of the Sun

We show that the higher range of the heliolongitudinal asymmetry of the solar wind speed in the positive polarity period (A > 0) than in the negative polarity period (A < 0) is one of the important reasons of the larger amplitudes of the 27-day variation of the galactic cosmic ray (GCR) intensity in the period of 1995–1997 (A > 0) than in 1985–1987 (A < 0). Subsequently, different ranges of the heliolongitudinal asymmetry of the solar wind speed jointly with equally important corresponding drift effect are general causes of the polarity dependence of the amplitudes of the 27-day variation of the GCR intensity. At the same time, we show that the polarity dependence is feeble for the last unusual minimum epoch of solar activity 2007–2009 (A < 0); the amplitude of the 27-day variation of the GCR intensity shows only a tendency of the polarity dependence. We present a three dimensional (3-D) model of the 27-day variation of GCR based on the Parker’s transport equation. In the 3-D model is implemented a longitudinal variation of the solar wind speed reproducing in situ measurements and corresponding divergence-free interplanetary magnetic field (IMF) derived from the Maxwell’s equations. We show that results of the proposed 3-D modeling of the 27-day variation of GCR intensity for different polarities of the solar magnetic cycle are in good agreement with the neutron monitors experimental data. To reach a compatibility of the theoretical modeling with observations for the last minimum epoch of solar activity 2007–2009 (A < 0) a parallel diffusion coefficient was increased by ∼40%.     

Plasma electron temperature variability in lunar surface potential and in electric field under average solar wind conditions

The Moon is immersed in plasma environment. The most interesting challenge of the lunar plasma– field environment is that it is alternatively dominated by the extended but variable outer atmosphere of the Earth – the magnetosphere – and by the extended but highly variable solar atmosphere – the solar wind. Understanding the plasma environment and its interaction with the lunar surface will be beneficial to both manned and robotic surface exploration activities and to scientific investigations. Presented is a preliminary map of variations of lunar surface electric potential over the day side and night side using probe equations and a discussion on dust dynamics in this E-field structure using the data from Electron Reflectometer in Lunar Prospector spacecraft during 1998–1999. On the day side, potential is around 5 V and on the night side it reaches up to −82 V. On the night side region, only highly energetic electrons can overcome this large negative potential. The variation at electron temperature (T_e) strongly reflects in the surface potential. The potential reaches to a value of −82 V for T_e = 58 eV. Surface charging causes the electrostatic transport of charged dust grains. Dust grain size of 0.1 μm shows a levitation height of 4.92 m on lunar day side, 748 m on terminator region and 3.7 km on the night side. The radius of maximum sized grain to be lofted, R_max, peaks at the terminator region (R_max = 0.83 μm). At the transition region dust levitation is almost absent. This region is most suited for exploration activities as the region is free from hazards caused by lunar dust.     

Variation of dayside chorus waves associated with solar wind dynamic pressure based on MMS observations

The whistler-mode chorus waves are one of the most important plasma waves in the Earth’s magnetosphere. Generally, the amplitude of whistler-mode chorus waves prefers to strengthen when the energetic fluxes of anisotropic electrons increase outside the plasmapause. This characteristic is commonly associated with the geomagnetic storms or substorms. However, the relationship between the solar wind dynamic pressure (P_sw) and the long-time variation of chorus waves during the quiet period of the geomagnetic activity still needs more detailed investigations. In this paper, based on MMS observations, we present a chorus event just observed in the inner side of magnetopause without obvious geomagnetic storms or substroms. Interestingly, during this time interval, some P_sw fluctuations were recorded. Both the amplitudes and frequencies of chorus waves changed as a response to the variation in P_sw. It proved that the enhancement of P_sw increases the energetic electrons fluxes, which provides free energies for the chorus amplification. Furthermore, the wave growth rates calculated using linear theory increases and the central frequency of the chorus waves shifts to a higher frequency when the P_sw enhancement is greater, which are also consistent well with the observations. The results provide a direct evidence that the P_sw play an important role in the long-time variation of whistler-mode chorus waves inside the magnetopause.