Solar EUV flux variation and flare brightenings

On 2010 February 8, the Extreme ultraviolet (EUV) flux variation in 195 Å and flare brightening has been examined in different sizes of active regions by using SOHO/EIT, MDI and Hα$\alpha$ observational data. These three active regions represent a large active region with a sunspot group, a moderate active region without a sunspot and a small region with weak plage in Hα$\alpha$ band respectively. Our study shows that the main full disk EUV flux comes from active regions, especially from large active regions. The sudden increases of EUV flux are corresponding to the EUV flare brightenings. For the large active region, the local EUV 195 Å flux peaks are well correlated to that of the GOES X-ray flux. The EUV 195 Å flux peaking time of M-class flares delay GOES X-ray flux a few minutes. For the moderate active region, the local EUV 195 Å flux is not well correlated to GOES X-ray flux. The EUV 195 Å flare brightenings in the moderate active region appeared in the duration of sudden increase of its own local EUV flux. For the small active region, the local EUV 195  Å flux varied almost independently of the GOES X-ray flux. Our study suggests that for an active region its local EUV 195 Å flux is more closely correlated to the EUV flare brightening than the full disk GOES X-ray flux.     

Solar magnetic feature detection and tracking for space weather monitoring

We present an automated system for detecting, tracking, and cataloging emerging active regions throughout their evolution and decay using SOHO Michelson Doppler Interferometer (MDI) magnetograms. The SolarMonitor Active Region Tracking (SMART) algorithm relies on consecutive image differencing to remove both quiet-Sun and transient magnetic features, and region-growing techniques to group flux concentrations into classifiable features. We determine magnetic properties such as region size, total flux, flux imbalance, flux emergence rate, Schrijver’s R-value, R^∗ (a modified version of R), and Falconer’s measurement of non-potentiality. A persistence algorithm is used to associate developed active regions with emerging flux regions in previous measurements, and to track regions beyond the limb through multiple solar rotations. We find that the total number and area of magnetic regions on disk vary with the sunspot cycle. While sunspot numbers are a proxy to the solar magnetic field, SMART offers a direct diagnostic of the surface magnetic field and its variation over timescale of hours to years. SMART will form the basis of the active region extraction and tracking algorithm for the Heliophysics Integrated Observatory (HELIO).     

Flux emergence,flux imbalance,magnetic free energy and solar flares

Emergence of complex magnetic flux in the solar active regions lead to several observational effects such as a change in sunspot area and flux embalance in photospheric magnetograms. The flux emergence also results in twisted magnetic field lines that add to free energy content. The magnetic field configuration of these active regions relax to near potential-field configuration after energy release through solar flares and coronal mass ejections. In this paper, we study the relation of flare productivity of active regions with their evolution of magnetic flux emergence, flux imbalance and free energy content. We use the sunspot area and number for flux emergence study as they contain most of the concentrated magnetic flux in the active region. The magnetic flux imbalance and the free energy are estimated using the HMI/SDO magnetograms and Virial theorem method. We find that the active regions that undergo large changes in sunspot area are most flare productive. The active regions become flary when the free energy content exceeds 50% of the total energy. Although, the flary active regions show magnetic flux imbalance, it is hard to predict flare activity based on this parameter alone.     

Variation in the ionospheric propagating factor M(3000)F2 at Ouagadougou,Burkina Faso

We use hourly monthly median values of propagation factor M(3000)F2 data observed at Ouagadougou Ionospheric Observatory (geographic12.4°N, 1.5°W; 5.9^o dip), Burkina Faso (West Africa) during the years Januar1987–December1988 (average F10.7 < 130 × 10⁻²² W/m²/Hz, representative of low solar flux conditions) and for January 1989–December1990 (average F10.7 ? 130 × 10⁻²² W/m²/Hz, representative of high solar epoch) for magnetically quiet conditions to describe local time, seasonal and solar cycle variations of equatorial ionospheric propagation factor M(3000)F2 in the African region. We show that that seasonal trend between solar maximum and solar minimum curves display simple patterns for all seasons and exhibits reasonable disparity with root mean square error (RMSE) of about 0.31, 0.29 and 0.26 for December solstice, June solstice and equinox, respectively. Variability Σ defined by the percentage ratio of the absolute standard deviation to the mean indicates significant dissimilarity for the two solar flux levels. Solar maximum day (10–14 LT) and night (22–02 LT) values show considerable variations than the solar minimum day and night values. We compare our observations with those of the IRI 2007 to validate the prediction capacity of the empirical model. We find that the IRI model tends to underestimate and overestimate the observed values of M(3000)F2, in particular, during June solstice season. There are large discrepancies, mainly during high solar flux equinox and December solstice between dawn and local midnight. On the other hand, IRI provides a slightly better predictions for M(3000)F2 between 0900 and 1500 LT during equinox low and high solar activity and equinox high sunspot number. Our data are of great importance in the area of short-wave telecommunication and ionospheric modeling.     

Reconstruction of solar irradiance using the Group sunspot number

We present a reconstruction of total solar irradiance since 1610 to the present based on variations of the surface distribution of the solar magnetic field. The latter is calculated from the historical record of the Group sunspot number using a simple but consistent physical model. Our model successfully reproduces three independent data sets: total solar irradiance measurements available since 1978, total photospheric magnetic flux from 1974 and the open magnetic flux since 1868 (as empirically reconstructed from the geomagnetic aa-index). The model predicts an increase in the total solar irradiance since the Maunder Minimum of about 1.3 Wm⁻².     

Photospheric flows in the active regions (asymmetric and localized Doppler velocities)

We observed 10 active regions through their disk passage during June 25–August 25, 1988, with the Tower Vector Magnetograph (TVM) of Marshall Space Flight Center. The TVM was used in scanning mode to measure the photospheric Doppler velocities with the Line-Center-Magnetogram (LCM) technique in the spectral line of FeI 5250.2 Å. In this paper we present the result of a subset of observations obtained while the active regions were situated away from the solar limb. A wide range of magnetic complexity and associated chromospheric activity characterized these active regions. It was found that the value of zero-crossing wavelength of the integrated Stokes-V profile of two opposite magnetic polarities were different, corresponding to Doppler velocities ranging from ∼100 m s⁻¹ to ∼1475 m s⁻¹. The measurements of relative velocities between different locations, connected by magnetic flux tubes as inferred from YOHKOH soft X-ray and TRACE 171 Å Fe IX images, showed widely different values of dominant localized flows. The region of parasite polarity, which showed recurrent chromospheric activity, was blue shifted with respect to the main “magnetic element” of the same polarity. Some of them were also the sites of sheared magnetic field configuration. The magnitude of the relative velocity between the leading and following polarity is more for the active regions of higher “field asymmetry”.     

基于PCA方法对太阳活动区主要参量的分析

太阳耀斑与太阳质子事件的发生通常与太阳活动区存在非常密切的关系, 对这种关系的深入分析有助于太阳耀斑和太阳质子事件预报模型的建立. 本文利用主成分分析(Principal Component Analysis, PCA)方法对1997-2010年太阳质子事件所在活动区的主要参量进行分析, 选取的参量包括黑子磁分类、 McIntosh分类、太阳黑子群面积、10.7 cm射电流量、耀斑指数、质子耀斑位置和软X射线耀斑强度. 结果得到81个太阳活动主成分得分值排序(得分值代表每个事件的强弱), 与太阳质子事件峰值流量、太阳黑子年均值以及10.7 cm射电流量年均值的对比显示相似度非常高, 表明主成分得分值一定程度上可以反映太阳活动的强弱规律.      

Tilts and polarity separation in the sunspot groups and active regions at the ascending and descending phases of the cycle 23

The Solar Feature Catalogues for sunspots and active regions measured with SOHO/MDI instrument and Ca II K3 spectroheliograph of the Paris-Meudon Observatory are analyzed with the automated classification technique for sunspot groups and active region polarities. We report the first classification results for daily variations of tilt angles (normal and trigonometric ones) in sunspot groups (SG) and active (AR) regions in the cycle 23. The average normal tilts are presented for every year at the ascending and descending phases of the cycle 23 which are similar to those deduced by other authors for the cycles 19–22. The normal tilts of both the sunspot groups and active regions are shown to increase in the ascending phase and a decrease in the descending phase. Similar to SG and AR areas, the trigonometric tilts are shown to have the noticeable North–South asymmetry with the Southern hemisphere dominant in the selected ascending and descending periods. The normal tilt variations with latitude follow Joy’s law revealing a periodicity along the meridian of about 10° and reaching the maximum of 14° at the latitude of about 32° corresponding to the top of the ‘royal zone’ where the sunspots appear. The variations of polarity separation with a latitude are in an anti-phase with those of the tilts reaching a maximum at the latitude of 35° and showing a small positive separation for the groups/active regions in a vicinity of the average tilts ±40°. The ratio R of the polarity separation to the trigonometric tilt fits the linear function of a latitude φ as R = −0.0213φ − 0.1245 confirming positive separation for the polarities of active regions with the average tilts, or the dominance of activity in the Southern hemisphere activity, for the selected period of observations.     

1986年2月太阳的高活动I活动区4711的演化和特征

本文使用太阳黑子、磁场、Hα色球、10.7cm射电及软X射线流量等观测资料,对太阳活动谷期的高活动区4711(SESC编号)从光球、色球和日冕三个方面做了综述.指出该活动区演化过程的特征是:(1)黑子群在主要发展阶段呈一个紧密的结构复杂的强磁区;(2)两次大的太阳爆发均发生在黑子群面积衰减阶段的初期;(3)黑子群的转动可能是活动区日冕加热和耀斑活动的主要供能机制;(4)色球暗条的频繁活动是爆发的先兆;(5) 10.7cm射电辐射和软X射线辐射的逐日流量有彼此不重合的双峰.      

Actors of the main activity in large complex centres during the 23 solar cycle maximum

During the maximum of Solar Cycle 23, large active regions had a long life, spanning several solar rotations, and produced large numbers of X-class flares and CMEs, some of them associated to magnetic clouds (MCs). This is the case for the Halloween active regions in 2003. The most geoeffective MC of the cycle (Dst = −457) had its source during the disk passage of one of these active regions (NOAA 10501) on 18 November 2003. Such an activity was presumably due to continuous emerging magnetic flux that was observed during this passage. Moreover, the region exhibited a complex topology with multiple domains of different magnetic helicities. The complexity was observed to reach such unprecedented levels that a detailed multi-wavelength analysis is necessary to precisely identify the solar sources of CMEs and MCs. Magnetic clouds are identified using in situ measurements and interplanetary scintillation (IPS) data. Results from these two different sets of data are also compared.     

Results from the study of an M7.6 flare and its associated CME

An M7.6 flare was well observed on October 24, 2003 in active region 10486 by a few instruments and satellites, including GOES, TRACE, SOHO, RHESSI and NoRH. Multi-wavelength study shows that this flare underwent two episodes. During the first episode, only a looptop source of <40 keV was observed in reconstructed RHESSI images, which showed shrinkage with a velocity of 12–14 km s⁻¹ in a period of about 12 min. During the second process, in addition to the looptop source, two footpoint sources were observed in energy channel of as high as ∼200 keV. One of them showed fast propagation along one of the two TRACE 1600 Å flare ribbons and the 195 Å loop footpoints, which could be explained by successive magnetic reconnection. The associated CME showed a mass pickup process with decreasing center-of-mass velocity. The decrease of the CME kinetic energy and the increase of its potential energy lead to an almost constant total energy during the CME propagation. Our results reveal that the flare and its associated CME have comparable energy content, and the flare is of non-thermal property.     

Mid-latitude thermospheric zonal winds during the equinoxes

The diurnal variation of the mid-latitude upper thermosphere zonal winds during equinoxes has been studied using data recently generated from CHAMP measurements from 2002 to 2004 using an iterative algorithm. The wind data was separated into two geomagnetic activity levels, representing high geomagnetic activity level (Ap > 8) and low geomagnetic activity level (Ap ? 8). The data were further separated into two solar flux levels; with F10.7 > 140 for high and F10.7 ? 140 for low. Geomagnetic activity is a correlator just as significant as solar activity. The response of mid-latitude thermospheric zonal winds to increases in geomagnetic disturbances and solar flux is evident. With increase in geomagnetic activity, midday to midnight winds are generally less eastward and generally more westward after the about midnight transitions. The results show that east west transitions generally occurred about midnight hours for all the situations analyzed. The west to east transition occurs from 1400–1500 MLT. Enhanced westward averaged zonal wind speeds going above 150 ms⁻¹ are observed in the north hemisphere mid-latitude about sunrise hours (∼0700–1100 MLT). Nighttime winds in the north hemisphere are in good agreement with previous single station ground observations over Millstone Hill. Improved ground observations and multi satellite observations from space will greatly improve temporal coverage of the Earth’s thermosphere.     

Cluster observations of earthward propagating plasmoid and flux ropes in the near-tail during the course of a substorm event

We analyze the magnetic structures in the near-tail at Xgsm = −17.5 Re on September 19, 2003 by Cluster. During the course of a substorm event, the earthward propagating plasmoid and flux ropes in the near-tail are observed. The earthward propagating plasmoid is associated with the bipolar Bz and By signatures. The two flux ropes are embedded within the earthward plasma flows, which might be referred to the population as ‘‘BBF-type’’ flux ropes. The first flux rope diameter is about 0.7 Re and duration based upon the Bz signature is ∼20 s, while the second one diameter is about 1.4 Re and duration is ∼30 s. The earthward propagating plasmoid and flux ropes could have influence upon the dipolarization and injection in inner magnetosphere. The Cluster observations of earthward propagating plasmoid and flux ropes can be interpreted as strong evidence for multiple X-lines. Our observations are consistent with that multiple plasmoids or flux ropes are formed repeatedly and ejected tailward in the course of geomagnetically active time.     

Analytical-HZETRN model for rapid assessment of active magnetic radiation shielding

The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA’s radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than ∼15 Tm and passive/structural shielding thicknesses greater than 40 g/cm² have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.     

Solar cycle effect on temporal and spatial variation of topside ion density measured by SROSS C2 and ROCSAT 1 over the Indian longitude sector

We investigated the diurnal, seasonal and latitudinal variations of ion density Ni over the Indian low and equatorial topside ionosphere within 17.5°S to 17.5°N magnetic latitudes by combining the data from SROSS C2 and ROCSAT 1 for the 9 year period from 1995 to 2003 during solar cycle 23. The diurnal maximum density is found in the local noon or in the afternoon hours and the minimum occurs in the pre sunrise hours. The density is higher during the equinoxes as compared to that in the June and December solstice. The local time spread of the daytime maximum ion density increases with increase in solar activity. A north south asymmetry with higher ion density over northern hemisphere in the June solstice and over southern hemisphere in December solstice has been observed in moderate and high solar activity years. The crest to crest distance increases with increase in solar flux. Ion density bears a nonlinear relationship with F_10.7 cm solar flux and EUV flux in general. The density increases linearly with solar flux up to ∼150 sfu (1 sfu = 10⁻²²Wm⁻²Hz⁻¹) and EUV flux up to ∼50 units (10⁹ photons cm⁻² s⁻¹). But beyond this the density saturates. Inverse saturation and linear relationship have been observed in some season or latitude also. Inter-comparison of the three solar activity indices F_10.7 cm flux, EUV flux and F_10.7P (= (F_10.7 + F_10.7A)/2, where F_10.7A is the 81 day running average value of F_10.7) shows that the ion density correlates better with F_10.7P and F_10.7 cm fluxes. The annual average daytime total ion density from 1995 to 2003 follows a hysteresis loop as the solar cycle reverses. The ion density at 500 km over the Indian longitude sector as obtained by the international reference ionosphere is in general lower than the measured densities during moderate and high solar activity years. In low solar activity years the model densities are equal or higher than measured densities. The IRI EIA peaks are symmetric (±10°) in equinox while densities are higher at 10°N in June solstice and at 10°S in the December solstice. The model density follows F_10.7 linearly up to about F_10.7 > ∼150 sfu and then saturates.     

Coordinated in situ measurements of plasma irregularities and ground based scintillation observations at the crest of equatorial anomaly

First comparison of in situ density fluctuations measured by the DEMETER satellite with ground based GPS receiver measurements at the equatorial anomaly station Bhopal (geographic coordinates (23.2°N, 77.6°E); geomagnetic coordinates (14.29°N, 151.12°E)) for the low solar activity year 2005, are presented in this paper. Calculation of the diurnal maximum of the strength of the equatorial electrojet, which can serve as precursor to ionospheric scintillations in the anomaly region is also done. The Langmuir Probe experiment and Plasma Analyzer onboard DEMETER measure the electron and ion densities respectively. Irregularities in electron density distribution cause scintillations on transionospheric links and there exists a close relationship between an irregularity and scintillation. In 40% of the cases, DEMETER detects the irregularity structures (dNe/Ne ? 5% and dNi/Ni (O⁺) ? 5%) and GPS L band scintillations (S4 ? 0.2) are also observed around the same time, for the low solar activity period. It is found that maximum irregularity intensity is obtained in the geomagnetic latitude range of 10–20° for both electron density and ion density. As the GPS signals pass through this irregularity structure, scintillations are recorded by the GPS receiver installed at the equatorial anomaly station, Bhopal it is interesting to note that in situ density fluctuations observed on magnetic flux tubes that pass over Bhopal can be used as indicator of ionospheric scintillations at that site. Many cases of density fluctuations and associated scintillations have been observed during the descending low solar activity period. The percentage occurrence of density irregularities and scintillations shows good correspondence with diurnal maximum of the strength of electrojet, however this varies with different seasons with maximum correspondence in summer (up to 66%) followed by equinox (up to 50%) and winter (up to 46%). Also, there is a threshold value of EEJ strength to produce density irregularities ((dNe/Ne)_max ? 5%) and for moderate to strong scintillations (S4 ? 0.3) to occur. For winter this value is found to be ∼40 nT whereas for equinox and summer it is around 50 nT.     

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.     

Latitudinal variation of F2-region response to geomagnetic disturbance

The effect of geomagnetic storms on the F2 region was studied by calculating the deviation, ΔfoF2, of foF2 during 40 magnetic storms, ranging from moderate (Dst < −50 nT) to very intense (Dst < −200 nT) of the 21st solar cycle. In order to study the variation of storm-time foF2 with latitude, season and storm strength, ionosonde data were obtained from eight stations spanning a latitudinal range of +60–−60°. The stations chosen lay in a narrow longitudinal range of 140–151°, so that local time difference between the stations is practically negligible. The features exhibited by positive and negative phases were essentially different. The storm time ΔfoF2 clearly exhibited a latitudinal variation and this variation were found to be coupled with the seasonal variation. As for the variation with storm intensity, though ΔfoF2 was found to vary even between two storms of almost equal intensity, the amplitude of a positive or negative phase, |ΔfoF2max| showed a distinct upper limit for each intensity category of storms.     

Study of the 27-day variations of the galactic cosmic ray intensity and anisotropy

We demonstrate that the general features of the radial and azimuthal components of the anisotropy of galactic cosmic rays can be studied by the harmonic analysis method using data from an individual neutron monitor with cut off rigidity <5 GV. In particular, we study the characteristics of the 27-day (solar rotation period) variations of the galactic cosmic ray intensity and anisotropy, solar wind velocity, interplanetary magnetic field strength and sunspot number. The amplitudes of the 27-day variations of the galactic cosmic ray anisotropy are greater, and the phases more clearly established, in A > 0 polarity periods than in A < 0 polarity periods at times of minimum solar activity. The phases of the 27-day variations of the galactic cosmic rays intensity and anisotropy are opposite with respect to the similar changes of the solar wind velocity in A > 0 polarity periods. No significant dependence of the amplitude of the 27-day variation of the galactic cosmic ray anisotropy on the tilt angle of the heliospheric neutral sheet is found. Daily epicyclegrams obtained by Chree’s method show that the 27-day variations of the galactic cosmic ray anisotropy during A > 0 polarity periods follow elliptical paths with the major axes oriented approximately along the interplanetary magnetic field. The paths are more irregular during A < 0 polarity periods.