Solar particle precipitation into the polar atmosphere and their dependence on hemisphere and local time

The precipitation of solar energetic particles, protons as well as electrons, at high latitudes is commonly assumed to be homogeneous across both polar caps. Using Low-Earth Orbit POES (Polar Orbiting Environmental Satellites) we determine particle penetration ratios into the polar atmosphere for protons ranging from about 0.1 MeV to 500 MeV and for electrons spanning about one order of magnitude in energy with a maximum of 0.3 MeV. Based on power law fits for the POES spectrum we show, that for energies interesting for middle and lower atmosphere chemistry, particle flux over the poles is comparable in magnitude to flux at the geostationary orbit or at L1 in interplanetary space. The time period under study are the solar energetic particle (SEP) event series of October/November 2003 and January 2005.     

Variation of energetic particle precipitation with local magnetic time

The detailed study of the precipitation of magnetospheric particles into the atmosphere is complicated by the rather complex spatial configuration of the precipitation region and its variability with geomagnetic activity. In this paper we will introduce polar oval coordinates and apply them to POES observations of 30 keV to 2.5 MeV electrons and comparable protons to illustrate the dependence of particle precipitation on local time and geomagnetic activity. These coordinates also allow an easy separation of the spatial precipitation patterns of solar and magnetospheric particles. The results indicate that (a) the spatial precipitation pattern of energetic magnetospheric electrons basically follows the pattern of the field parallel Birkeland currents up to MeV energies and (b) at least in the mesosphere the influence of magnetospheric electrons is comparable to the one of solar electrons. Implications for modeling of atmospheric chemistry will be sketched.     

Some peculiarities of longitudinal distribution of proton fluxes at high latitudes

Dynamical features of proton fluxes at high and middle latitudes were studied based on data measured by Sun-synchronous low-altitude (800 km height) polar-orbiting three NOAA series satellites: POES 15, 16, and 17 during the geomagnetic storm on January, 21–22, 2005. Data from three satellites that passed the Northern hemisphere along different MLTs allow reconstructing the longitudinal distribution of the proton fluxes. Measurements of protons with energies of 30–80 keV and 80–240 keV (the ring current energy range) by 0- and 90-detectors were used to evaluate and compare the longitudinal asymmetry of proton flux distribution measured in the regions equatorward and poleward of the isotropic boundary. It was found that during all the phases of the geomagnetic storm distribution of the maximum flux of precipitating protons (0-detector data) is sufficiently asymmetric. The maximal flux position along MLT is moving from pre-midnight sector in quiet time to post-midnight one before and during SSC and moving back during recovery phase. The longitudinal distribution of precipitation maxima demonstrates the local increase in afternoon sector (approximately at 13:30 MLT) and decrease in the dusk one during SSC. These features are evident consequence of the magnetosphere compression. To identify the origin of the particles, the locations of maximum fluxes have been projected to the magnetosphere. It was determined that during geomagnetic storm main and recovery phases maximum fluxes were measured at latitudes poleward of the isotropic boundary. To evaluate the trapped particle flux asymmetry, the particle fluences (90-detector data) were calculated along the satellite orbit from L = 2 to the isotropic boundary. The total fluences of trapped particles calculated along the satellite orbit show regular asymmetry between dusk and dawn during main and recovery phases. The maximal intensity of proton fluxes of both investigated populations located poleward and equatorward of the isotropic boundary is achieved during SSC. The total flux measured during crossing the anisotropic region can be considered as a proxy for ring current injection rate.     

Enhancement of solar wind low-energy energetic particles as precursor of geomagnetic disturbance in operational geomagnetic forecast

A study of the relationship between solar wind low-energy energetic particles using data from the Electron, Proton, and Alpha Monitor (EPAM) onboard the Advanced Compositional Explorer spacecraft (ACE) and geomagnetic activity using data from Canadian magnetic observatories in Canada’s polar cap, auroral zone, and subauroral zone was carried out for a period spanning 1997–2005. Full halo coronal mass ejections (CMEs) were used to gauge the initial particle enhancements and the subsequent geomagnetic activity. It was found that maximum geomagnetic activity is related to maximum particle enhancements in a non-linear fashion. Quadratic fit of the data results in expressions that can be easily used in an operational space weather setting to forecast geomagnetic disturbance quantitatively. A superposed epoch analysis shows increase in particle flux level starts hours before geomagnetic activity attains its peak, affirming the precursory nature of EPAM particles for the impending geomagnetic impact of CME. This can supplement the decision process in formulating geomagnetic warning after the launch of CME from the Sun but before the arrival of shock at Earth. The empirical relationships between solar wind low-energy energetic particles and geomagnetic activity revealed in this statistical study can be easily codified, and thus utilized in operational space weather forecast to appraise the geoeffectiveness of the CME and to provide a quantitative forecast for maximum geomagnetic activity in Canada’s polar cap, auroral zone, and subauroral zone after the occurrence of a CME.     

Statistical survey on sawtooth events,SMCs and isolated substorms

Solar wind driving can cause a variety of different responses in the magnetosphere. Strong and steady driving during geomagnetic storms may result in sawtooth events. Strong to moderate driving may be followed by either sawtooth events or steady magnetospheric convection (SMC) events. Lower solar wind energy input typically leads to the formation of isolated non-storm substorms. This study uses superposed epoch analysis to reveal the typical properties of these three event groups as well as their similarities and differences. We use IMF and solar wind parameters, as well as ground-based indices (AL, SYM-H, ASY-H, PCN) to examine the level of solar wind driving and its response in the magnetosphere. Our results show that sawtooth events are associated with the strongest ionospheric activity. The subgroups of events during constant solar wind _EY$E_Y$ show that the key difference between the events is the average solar wind speed. Particularly, the high activity during sawtooth events is driven by high solar wind speed, while the lowest average speed during the SMCs may explain the lack of substorm activity during the steady convection periods.     

Changes in sea-level pressure over South Korea associated with high-speed solar wind events

We explore a possibility that the daily sea-level pressure (SLP) over South Korea responds to the high-speed solar wind event. This is of interest in two aspects: first, if there is a statistical association this can be another piece of evidence showing that various meteorological observables indeed respond to variations in the interplanetary environment. Second, this can be a very crucial observational constraint since most models proposed so far are expected to preferentially work in higher latitude regions than the low latitude region studied here. We have examined daily solar wind speed V, daily SLP difference ΔSLP, and daily log(BV²) using the superposed epoch analysis in which the key date is set such that the daily solar wind speed exceeds 800 km s⁻¹. We find that the daily ΔSLP averaged out of 12 events reaches its peak at day +1 and gradually decreases back to its normal level. The amount of positive deviation of ΔSLP is +2.5 hPa. The duration of deviation is a few days. We also find that ΔSLP is well correlated with both the speed of solar wind and log(BV²). The obtained linear correlation coefficients and chance probabilities with one-day lag for two cases are r ? 0.81 with P > 99.9%, and r ? 0.84 with P > 99.9%, respectively. We conclude by briefly discussing future direction to pursue.     

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.     

Long-term observations of the trapped high-energy proton population (L < 4) by the NOAA Polar Orbiting Environmental Satellites (POES)

The Space Environment Monitor (SEM) onboard the NOAA POES satellites has been measuring the near-Earth charged particle environment since 1978, providing an extensive database that can be used for studying the long-term behavior of this population of trapped particles. POES stands for Polar Orbiting Environmental Satellite. These satellites orbit at ∼840 km altitude and at an inclination of 98°. The SEM-1 instrument was flown on the POES satellites beginning in 1978 with TIROS-N and NOAA-6 in 1979 and continuing to NOAA-14 launched in 1995 with the exception of NOAA-9 and NOAA-11 (NOAA-13 failed shortly after launch). Its replacement, SEM-2, has flown on the POES NOAA-15, -16, and -17 satellites (from 1998). Here we present the results of a study on the statistical variations of the high-energy trapped proton environment. Among the detectors in SEM, the four SEM-2 omni-directional proton detectors for energies >16 MeV, >36 MeV, >70 MeV, and >140 MeV provide the data most relevant to this study.     

Simultaneous observations of Pc 1 micropulsation activity and stratospheric electrodynamic perturbations on 27 January 2003

The 2nd Polar Patrol Balloon campaign (2nd-PPB) was carried out at Syowa Station in Antarctica during 2002–2003. Identical stratospheric balloon payloads were launched as close together in time as allowed by weather conditions to constitute a cluster of balloons during their flights. A very pronounced negative ion conductivity enhancement was observed at 32 km in the stratosphere below the auroral zone on 27 January 2003 from 1500 to 2200 UT. During this event, the conductivity doubled for an interval of about 7 h. This perturbation was associated with an extensive Pc 1 or Pi 1 wave event that was observed by several Antarctic ground stations, balloon PPB 10, and the Polar spacecraft. No appreciable X-ray precipitation was observed in association with this event, which would point to >60 Mev proton precipitation as a possible magnetosphere–stratosphere coupling mechanism responsible for the conductivity enhancement. Such precipitation is consistent with the wave data. During the latter half of the event, E_z was briefly positive. There was a tropospheric Southern Ocean storm system underneath the balloon during this interval. If the event was associated with this storm system and not energetic proton precipitation, the observations imply an electrified Southern Ocean storm and major perturbations in stratospheric conductivity driven by a tropospheric disturbance. This event represents a poorly understood source for global circuit current. Precipitating energetic proton data from Akebono and NOAA POES spacecraft show significant >16 MeV precipitation was occurring at the location of PPB 8 but not PPB 10, suggesting that proton precipitation was, in fact, the responsible coupling mechanism.     

Different space weather effects in anomalies of the high and low orbital satellites

Preliminary results of the EU INTAS Project 00810, which aims to improve the methods of safeguarding satellites in the Earth’s magnetosphere from the negative effects of the space environment, are presented. Anomaly data from the “Kosmos” series satellites in the period 1971–1999 are combined in one database, together with similar information on other spacecraft. This database contains, beyond the anomaly information, various characteristics of the space weather: geomagnetic activity indices (Ap, AE and Dst), fluxes and fluences of electrons and protons at different energies, high energy cosmic ray variations and other solar, interplanetary and solar wind data. A comparative analysis of the distribution of each of these parameters relative to satellite anomalies was carried out for the total number of anomalies (about 6000 events), and separately for high (5000 events) and low (about 800 events) altitude orbit satellites. No relation was found between low and high altitude satellite anomalies. Daily numbers of satellite anomalies, averaged by a superposed epoch method around sudden storm commencements and proton event onsets for high (>1500 km) and low (<1500 km) altitude orbits revealed a big difference in a behavior. Satellites were divided on several groups according to the orbital characteristics (altitude and inclination). The relation of satellite anomalies to the environmental parameters was found to be different for various orbits that should be taken into account under developing of the anomaly frequency models.     

Ion production and ionization effect in the atmosphere during the Bastille day GLE 59 due to high energy SEPs

The influence of high energy particles, specifically cosmic rays, on atmospheric physics and chemistry is highly discussed. In most of the proposed models the role of ionization in the atmosphere due to cosmic rays is not negligible. Moreover, effect(s) on minor constituents and aerosols are recently observed, specifically over the polar regions during strong solar particle events. According to the recent findings for such effects it is necessary an essential increase of ion production, specifically during the winter period. The galactic cosmic rays are the main source of ionization in the Earth’s stratosphere and troposphere. Occasionally, the atmospheric ionization is significantly enhanced during strong solar energetic particles events, specifically over the polar caps. During the solar cycle 23 several strong ground level enhancements were observed. One of the strongest was the Bastille day event occurred on 14 July 2000. Using a full Monte Carlo 3-D model, we compute the atmospheric ionization, considering explicitly the contribution of cosmic rays with galactic and solar origin, focusing on high energy particles. The model is based on atmospheric cascade simulation with the PLANETOCOSMICS code. The ion production rate is computed as a function of the altitude above the sea level. The ion production rate is computed on a step ranging from 10 to 30?min throughout the event, considering explicitly the spectral and angular characteristics of the high energy part of solar protons as well as their time evolution. The corresponding event averaged ionization effect relative to the average due to galactic cosmic rays is computed in lower stratosphere and upper troposphere at various altitudes, namely 20?km, 15?km, 12?km and 8?km above the sea level in a sub-polar and polar regions. The 24^h and the weekly ionization effects are also computed in the troposphere and low stratosphere. Several applications are discussed.     

Evolution of solar wind structures from 0.72 to 1 AU

Understanding the evolution of solar wind structures in the inner heliosphere as they approach the Earth is important to space weather prediction. From the in situ solar wind plasma and magnetic field measurements of Pioneer Venus Orbiter (PVO) at 0.72 AU (1979–1988), and of Wind/Advanced Composition Explorer (ACE) missions at 1 AU (1995–2004), we identify and characterize two major solar wind structures, stream interaction regions (SIRs) and interplanetary coronal mass ejections (ICMEs). The average percentage of SIRs occurring with shocks increases significantly from 3% to 24% as they evolve from 0.72 to 1 AU. The average occurrence rate, radial extent, and bulk velocity variation of SIRs do not change from 0.72 to 1 AU, while peak pressure and magnetic field strength both decrease with the radial evolution of SIRs. Within the 0.28 AU distance from the orbit of Venus to that of Earth, the average fraction of ICMEs with shocks increases from 49% to 66%, and the typical radial extent of ICMEs expands by about a fraction of 1.4, with peak pressure and magnetic field strength decreasing significantly. The mean occurrence rate and expansion velocity of ICMEs do not change from 0.72 to 1 AU.     

太阳风-磁层-电离层耦合过程中的能量收支

太阳风向磁层电离层(Magnetosphere and Ionosphere, MI)系统输入能量, 而输入的能量随后在MI系统中消耗. 本文从能量守恒原理出发, 讨论太阳风-磁层-电离层 (SMI)耦合过程中的能流路径和能量收支的定量关系. 主要讨论9个问题: (1) 太阳风向MI系统的能量输入, (2) MI系统对能量输入的响应, (3) 环电流的能量消耗, (4) 极区电离层焦耳加热的能量消耗, (5) 极光粒子沉降的能量消耗, (6) 磁尾能量的消耗、储存以及返回下游太阳风, (7)平静期间的能量积累与释放, (8)能量在不同能汇中的分配, (9)评价能量函数的准则和方法.      

天气尺度到气候尺度太阳风变速对中高纬大气环流的影响

通过对冬季太阳风短时(天气尺度)降速与北大西洋涛动和北极涛动等北半球中高纬度环流指数的时序重叠分析,结合对1963年以来48个冬季太阳风平均速度与北极涛动等指数的相关分析发现,从短时太阳风降速到向亚极光带沉降的辐射带高能电子通量显著下降,北极涛动也有迅速的响应,这预示着从太阳风到大气环流存在天气尺度的短时关系链,在这一时间尺度现有理论中仅有“空间粒子-大气电-云微物理”联系机制能较好地解释;太阳风速度与北极涛动的正相关信号在气候尺度上也有显著体现,太阳风可能通过高能电子沉降与北半球冬季中高纬度环流相联系,这表明太阳风通过大气电-云微物理过程驱动的过程是太阳活动影响气候变化的不可忽视的途径;开展太阳风起源、空间环境与大气(环流、电场)和地磁系统的联合观测及数值模拟是揭示日地天气与气候联系的重要研究内容之一.      

Semi-empirical model of middle atmosphere wind from the ground to the lower thermosphere

During recent years, special attention has been paid to understanding the background circulation of the middle atmosphere. Particularly in the mesosphere/lower thermosphere (MLT) region, this has involved including data from a range of new radar measurements. It has also involved the comparison of existing empirical middle atmosphere wind models, such as CIRA-86 and HWM-93 to the new data. This has led to the construction of empirical models of MLT winds such as the Global Empirical Wind Model (GEWM). Further investigations are aimed at the construction of new empirical and semi-empirical wind models of the entire middle atmosphere including these new experimental results. The results of a new wind climatology (0–100 km) are presented here, based upon the GEWM, a reanalysis of stratospheric data, and a numerical model which is used to fill the gap between data from the stratospheric and MLT regions.     

Statistical properties of the IMF clock angle in the solar wind with northward and southward interplanetary magnetic field based on ACE observation from 1998 to 2009: Dependence on the temporal scale of the solar wind

Utilizing ACE satellite observations from 1998 to 2009, we performed the elaborate study on the properties of the clock angle θ_CA (arctan(By/Bz) (?90° to 90°) of the interplanetary magnetic field (IMF) in the solar wind at 1?AU. The solar wind with northward IMF (NW-IMF) and southward IMF (SW-IMF) are analyzed, independently. Statistical analysis shows that the solar wind with SW-IMF and NW-IMF has similar properties in general, including their durations, the IMF Bz and By components, and the IMF θ_CA. Then, the solar wind with NW-IMF (SW-IMF) is classified into five different temporal scales according to the duration of the NW-IMF (SW-IMF), i.e., very-short wind of 10–30?min, short-scale wind of 0.5–1?h, moderate-scale wind of 1–3?h, long-scale wind of 3–5?h, and super-long wind >5?h. Our analysis reveals that the IMF θ_CA has a distinct decrease with increase of the temporal scale of the solar wind. Next, the solar wind is classified into two groups, i.e., the high-speed solar wind (>450?km/s) and the low-speed solar wind (<450?km/s). Our analysis indicates that the IMF θ_CA depends highly on the solar wind speed. Statistically, high-speed solar wind tends to have larger IMF θ_CA than low-speed solar wind. The evolutions of the solar wind and IMF with the solar activity are further studied, revealing no clear solar variation of the IMF θ_CA. Finally, we analyze the monthly variation of the IMF θ_CA. Superposed epoch result strongly suggests the seasonal variation of the IMF θ_CA.     

Tidal associated temperature disturbances observed at the middle atmosphere (30–65 km) by a Rayleigh lidar at 23°S

A lidar has been operated in São José dos Campos, Brazil (23.2°S, 45.8°W) since 1972, mainly dedicated to the study of mesospheric sodium at the 589 nm resonant line. The molecular Rayleigh scattering can also be used provided we limit the height to ∼75 km where the sodium scattering begins. Nevertheless, the weak signal obtained only permits the determination of density and temperature profiles by accumulating a large number of shots giving only nocturnal average profiles. Temporal variations in density and temperature on the scale of hours can however, be obtained by performing a superposed epoch analysis for a given time interval and covering a period of several days. In this way we obtained hourly mean profiles grouped by months, seasons and overall, with data acquired from 1993 to 2004. The difference between the hourly temperatures and the nocturnal means shows for some months, with enough data coverage, downward propagating structures that apparently have tidal origin. The seasonal averages show a recurrent feature with high temperatures before and low temperatures after midnight above 50 km. Some similarity is found with the GSWM model, but the observed temperature amplitudes are twice of that for the model.     

Gravity waves and wind-shear in the MLT at 23°S

We have used the technique suggested by Hocking [Hocking, W. A new approach to momentum flux determinations using SKiYMET meteor radars. Ann. Geophys. 23, 2005.] to derive short period wind variances in the 80–100 km region from meteor radar data. We find that these fluctuating winds, assumed to correspond to gravity waves and turbulence, are closely correlated with the vertical shear of the horizontal tidal winds. This close correlation suggests that in situ wind shear may be a major source of gravity waves and turbulence in the MLT. If this is the case, gravity waves generated in the troposphere and propagating up to the MLT region, generally assumed to constitute an important influence on the climatology of the region, may be a less important source of energy and momentum in the 80–100 km region than has been hitherto believed.     

Neutral winds in the polar thermosphere as measured from dynamics explorer

Remote sensing measurements of the meridional thermospheric neutral wind using the Fabry-Perot Interferometer on Dynamics Explorer have been combined within-situ measurements of the zonal component using the Wind and Temperature Spectrometer on the same spacecraft. The two data sets with appropriate spatial phasing and averaging determine the vector wind along the track of the polar orbiting spacecraft. A study of fifty-eight passes over the Southern (sunlit) pole has enabled the average Universal Time dependence of the wind field to be determined for essentially a single solar local time cut. The results show the presence of a “back-ground” wind field driven by solar EUV heating upon which is superposed a circulating wind field driven by high latitude momentum and energy sources.     

Wave perturbations in the MLT at high northern latitudes in winter,observed by two SATI instruments

To investigate the Mesosphere and Lower Thermosphere (MLT) region, several ground-based instruments called SATI (Spectral Airglow Temperature Imager) were designed and built to measure airglow emission and temperature in the upper mesosphere. One SATI instrument was installed at Resolute Bay (74.7°N, 94.9°W) and has monitored the polar MLT region since November, 2001. In October 2007 another SATI instrument was installed at Eureka (80.0°N, 86.3°W) at the Polar Environment Atmospheric Research Laboratory (PEARL) as part of the Canadian Network for the Detection of Atmospheric Change (CANDAC) project. SATI is a spatial scanning Fabry–Perot spectrometer measuring column emission rates for several rotational lines of OH and O₂ airglow at 87 and 94 km height. The rotational temperatures are inferred from the ratios of these lines. The measurements are divided into 12 sectors with an annular field of view. The phase differences between the sectors yield information on the horizontal atmospheric wave direction and wavelength. Horizontal perturbations of 2–8 h period have correlatively been observed and investigated at both locations. Short-periodic oscillations identified as gravity waves with periods between 2 and 8 h propagate in southward and eastward directions, but in opposite directions in some cases. The wave propagation characteristics are often different at the two locations; the relationship with the lower mean wind is considered.