New Insight to the Surface Temperature of the Sun

Information is given on thermal radiation from the Sun, considered in practical engineering calculations of heat exchange. It was found that although the surface temperature of the Sun is assumed to be about 5800 K, the solar spectrum data measured by Kondratyev lead to a value of at least 7134 K. Such a higher value can be obtained by interpreting the Planck formula for the black radiation spectrum for the Kondratyev data. In addition, using the Stefan-Boltzmann law, the energetic emissivity of the Sun’s surface was determined to be 0.431. Furthermore, based on Petela’s formulae for exergy of thermal radiation, the exergetic emissivity of the Sun’s surface was also calculated at the level of 0.426.

Location of energy source for coronal heating on the photosphere

It is reported that ultra-fine dynamic ejections along magnetic loops of an active region origi- nate from intergranular lanes and they are associated with subsequent heating in the corona. As continu- ing work, we analyze the same set of data but focus on a quiet region and the overlying EUV/UV emis- sion as observed by the Atmospheric Imaging Assembly （AIA） on board Solar Dynamics Observatory （SDO）. We find that there appear to be dark patches scattered across the quiet region and the dark patches always stay along intergranular lanes. Over the dark patches, the average UV/EUV emission at 131, 17 1, 304 and 1600 A （middle temperature） is more intense than that of other regions and EUV brightness is negatively correlated with 10830A intensity, though, such a trend does not exist for high temperature lines at 94, 193, 211 and 335 A. For the same quiet region, where both TiO 7057 A broad band images and 10830A filtergrams are available, contours for the darkest lane areas on TiO images and clark patches on 10830A filtergrams frequently differ in space. The results suggest that the dark patches do not simply reflect the areas with the darkest lanes but are associated with a kind of enhanced absorption （EA） at 10830A,. A strict definition for EA with narrow band 10830A filtergrams is found to be difficult. In this paper, we define enhanced absorption patches （EAPs） of a quiet region as the areas where emission is less than ,-90% of the mean intensity of the region. The value is equivalent to the average intensity along thin dark loops connecting two moss regions of the active region. A more strict definition for EAPs, say 88%, gives even more intense UV/EUV emission over those in the middle temperature range. The results provide further observational evidence that energy for heating the upper solar atmosphere comes from the intergranular lane area where the magnetic field is constantly brought in by convection motion in granules.

Numerical experiments on the evolution in coronal magnetic configurations including a filament in response to the change in the photosphere

We investigate equilibrium height of a flux rope, and its internal equilib- rium in a realistic plasma environment by carrying out numerical simulations of the evolution of systems including a current-carrying flux rope. We find that the equilib- rium height of a flux rope is approximately described by a power-law function of the relative strength of the background field. Our simulations indicate that the flux rope can escape more easily from a weaker background field. This further confirms that a catastrophe in the magnetic configuration of interest can be triggered by a decrease in strength of the background field. Our results show that it takes some time to reach internal equilibrium depending on the initial state of the flux rope. The plasma flow inside the flux rope due to the adjustment for the internal equilibrium of the flux rope remains small and does not last very long when the initial state of the flux rope com- mences from the stable branch of the theoretical equilibrium curve. This work also confirms the influence of the initial radius of the flux rope in its evolution; the results indicate that a flux rope with a larger initial radius erupts more easily. In addition, by using a realistic plasma environment and a much higher resolution in our simulations, we notice some different characteristics compared to previous studies in Forbes.

Magnetic non-potentiality on the quiet Sun and the filigree

From the observed vector magnetic fields by the Solar Optical Telescope/ Spectro-Polarimeter aboard the satellite Hinode, we have examined whether or not the quiet Sun magnetic fields are non-potential, and how the G-band filigrees and Ca II network bright points （NBPs） are associated with the magnetic non-potentiality. A sizable quiet region in the disk center is selected for this study. The new findings by the study are as follows. （1） The magnetic fields of the quiet region are obviously non-potential. The region-average shear angle is 40°, the average vertical current is 0.016A m^-2, and the average free magnetic energy density, 2.7× 10^2erg cm^-3. The magnitude of these non-potential quantities is comparable to that in solar active regions. （2） There are overall correlations among current helicity, free magnetic energy and longitudinal fields. The magnetic non-potentiality is mostly concentrated in the close vicinity of network elements which have stronger longitudinal fields. （3） The filigrees and NBPs are magnetically characterized by strong longitudinal fields, large electric helicity, and high free energy density. Because the selected region is away from any enhanced network, these new results can generally be applied to the quiet Sun. The findings imply that stronger network elements play a role in high magnetic non-potentiality in heating the solar atmosphere and in conducting the solar wind.

A study of the relation between intensity oscillations and magnetic field parameters in a sunspot: Hinode observations

We present properties of intensity oscillations of a sunspot in the photo- sphere and chromosphere using G band and Ca u H filtergrams from Hinode. Intensity power maps as function of magnetic field strength and frequency reveal reduction of power in the G band with an increase in photospheric magnetic field strength at all frequencies. In Ca II H, however, stronger fields exhibit more power at high frequen- cies, particularly in the 4.5-8.0 mHz band. Power distributions in different locations of the active region show that the oscillations in Ca II H exhibit more power compared to that of the G band. We also relate the power in intensity oscillations with differ- ent components of the photospheric vector magnetic field using near simultaneous spectro-polarimetric observations of the sunspot from the Hinode spectropolarime- ter. The photospheric umbral power is strongly anti-correlated with the magnetic field strength and its line-of-sight component but there is a good correlation with the trans- verse component. A reversal of this trend is observed in the chromosphere except at low frequencies （V≤ 1.5 mHz）. The power in sunspot penumbrae is anti-correlated with the magnetic field parameters at all frequencies （1.0 ≤ v ≤ 8.0 mHz） in both the photosphere and chromosphere, except that the chromospheric power shows a strong correlation in the frequency range 3-3.5 mHz.

Data-driven Simulations of Magnetic Field Evolution in Active Region 11429:Magneto-frictional Method Using PENCIL CODE

Coronal magnetic fields evolve quasi-statically over long timescales and dynamically over short timescales.As of now there exist no regular measurements of coronal magnetic fields,and therefore generating the coronal magnetic field evolution using observations of the magnetic field at the photosphere is a fundamental requirement to understanding the origin of transient phenomena from solar active regions(ARs).Using the magneto-friction(MF)approach,we aim to simulate the coronal field evolution in the solar AR 11429.The MF method is implemented in the open source PENCIL CODE along with a driver module to drive the initial field with different boundary conditions prescribed from observed vector magnetic fields at the photosphere.In order to work with vector potential and the observations,we prescribe three types of bottom boundary drivers with varying free-magnetic energy.The MF simulation reproduces the magnetic structure,which better matches the sigmoidal morphology exhibited by Atmospheric Imaging Assembly(AIA)images at the pre-eruptive time.We found that the already sheared field further driven by the sheared magnetic field will maintain and further build the highly sheared coronal magnetic configuration,as seen in AR 11429.Data-driven MF simulation is a viable tool to generate the coronal magnetic field evolution,capturing the formation of the twisted flux rope and its eruption.

Lossless Compression Method for the Magnetic and Helioseismic Imager(MHI)Payload

The Solar Polar-orbit Observatory(SPO),proposed by Chinese scientists,is designed to observe the solar polar regions in an unprecedented way with a spacecraft traveling in a large solar inclination angle and a small ellipticity.However,one of the most significant challenges lies in ultra-long-distance data transmission,particularly for the Magnetic and Helioseismic Imager(MHI),which is the most important payload and generates the largest volume of data in SPO.In this paper,we propose a tailored lossless data compression method based on the measurement mode and characteristics of MHI data.The background out of the solar disk is removed to decrease the pixel number of an image under compression.Multiple predictive coding methods are combined to eliminate the redundancy utilizing the correlation(space,spectrum,and polarization)in data set,improving the compression ratio.Experimental results demonstrate that our method achieves an average compression ratio of 3.67.The compression time is also less than the general observation period.The method exhibits strong feasibility and can be easily adapted to MHI.

Characterizing motion types of G-band bright points in the quiet Sun

We study the motion of G-band bright points （GBPs） in the quiet Sun to obtain the characteristics of different motion types. A high resolution image sequence taken with the Hinode/Solar Optical Telescope （SOT） is used, and GBPs are automat- ically tracked by segmenting 3D evolutional structures in a space-time cube. After putting the GBPs that do not move during their lifetimes aside, the non-stationary GBPs are categorized into three types based on an index of their motion type. Most GBPs that move in straight or nearly straight lines are categorized as a straight mo- tion type, a few moving in rotary paths into rotary motion, and the others fall into a motion type we called erratic. The mean horizontal velocities are 2.18±0.08 km s-1, 1.63±0.09km s^-1 and 1.33±0.07 km s^-1 for straight, erratic and rotary motion types, respectively. We find that a GBP drifts at a higher and constant velocity during its whole life if it moves in a straight line. However, it has a lower and variational velocity if it moves on a rotary path. The diffusive process is ballistic-, super- and sub-diffusion for straight, erratic and rotary motion types, respectively. The corresponding diffusion index （γ） and coefficients （K） are 2.13±0.09 and 850±37km^2 s^-1, 1.82±0.07 and 331 ±24 km^2 s^-1, and 0.73±0.19 and 13±9 km^2 s^-1. In terms of direction of motion, it is homogeneous and isotropic, and usually persists between neighboring frames, no matter what motion type a GBP is classified as.

Studies of grand minima in sunspot cycles by using a flux transport solar dynamo model

We propose that grand minima in solar activity are caused by simultane- ous fluctuations in the meridional circulation and the Babcock-Leighton mechanism for the poloidal field generation in the flux transport dynamo model. We present the following results： （a） fluctuations in the meridional circulation are more effective in producing grand minima; （b） both sudden and gradual initiations of grand minima are possible; （c） distributions of durations and waiting times between grand minima seem to be exponential; （d） the coherence time of the meridional circulation has an effect on the number and the average duration of grand minima, with a coherence time of about 30 yr being consistent with observational data. We also study the occurrence of grand maxima and find that the distributions of durations and waiting times between grand maxima are also exponential, like the grand minima. Finally we address the question of whether the Babcock-Leighton mechanism can be operative during grand minima when there are no sunspots. We show that an a-effect restricted to the upper portions of the convection zone can pull the dynamo out of the grand minima and can match various observational requirements if the amplitude of this a-effect is suitably fine-tuned.

What does the Sun tell and hint now?

Some historical records, which have held since the beginning of modern solar activity cycles, are being broken by the present Sun： cycle 23 records the longest cycle length and fall time; latitudes of high-latitude sunspots belonging to a new cycle around the minimum time of the cycle are statistically the lowest at present, compared with those of other cycles; there are only one or no sunspots in a month appearing at high latitudes for 58 months, which is the first time that such a long duration has been observed. The solar dynamo is believed to be slowing down due to： （1） the minimum smoothed monthly mean sunspot number is the smallest since cycle 16 onwards, and even probably among all modern solar cycles; and （2） once the time interval between the first observations of two neighboring sunspot groups is larger than 14 d, it should be approximately regarded as an observation of no sunspots on the visible solar disk, called a spotless event. Spotless events occur with the highest frequency around the minimum time of cycle 24, and the longest spotless event also appears around the minimum time for observations of the Sun since cycle 16. Cycle 24 is expected to have the lowest level of sunspot activity from cycle 16 onwards and even probably for all of the modern solar cycles.

A Statistical Study of Rapid Sunspot Structure Change Associated with Flares

We reported recently some rapid changes of sunspot structure in white-light （WL） associated with major flares. We extend the study to smaller events and present here results of a statistical study of this phenomenon. In total, we investigate 403 events from 1998 May 9 to 2004 July 17, including 40 X-class, 174 M-class, and 189 C-class flares. By monitoring the structure of the flaring active regions using the WL observations from the Transition Region and Coronal Explorer （TRACE）, we find that segments in the outer sunspot structure decayed rapidly right after many flares; and that, on the other hand, the central part of sunspots near the flare-associated magnetic neutral line became darkened. These rapid and permanent changes are evidenced in the time profiles of WL mean intensity and are not likely resulted from the flare emissions. Our study further shows that the outer sunspot structure decay as well as the central structure darkening are more likely to be detected in larger solar flares. For X-class flares, over 40% events show distinct sunspot structure change. For M- and C-class flares, this percentage drops to 17% and 10%, respectively. The results of this statistical study support our previously proposed reconnection picture, i.e., the flare-related magnetic fields evolve from a highly inclined to a more vertical configuration.

Solar flare forecasting based on sequential sunspot data

It is widely believed that the evolution of solar active regions leads to solar flares. However, information about the evolution of solar active regions is not employed in most existing solar flare forecasting models. In the current work, a short- term solar flare forecasting model is proposed, in which sequential sunspot data, in- cluding three days of information about evolution from active regions, are taken as one of the basic predictors. The sunspot area, the Mclntosh classification, the mag- netic classification and the radio flux are extracted and converted to a numerical for- mat that is suitable for the current forecasting model. Based on these parameters, the sliding-window method is used to form the sequential data by adding three days of information about evolution. Then, multi-layer perceptron and learning vector quanti- zation are employed to predict the flare level within 48 h. Experimental results indicate that the performance of the proposed flare forecasting model works better than previ- ous models.

Phase analysis of sunspot group numbers on both solar hemispheres

Cross-correlation analysis and wavelet transform methods are proposed to investigate the phase relationship between the monthly sunspot group numbers in the solar northern and southern hemispheres. It is found that （1） the monthly sunspot group numbers in the northern hemisphere begin two months earlier than those in the southern one, which should lead to phase asynchrony between them but with a slight effect; （2） the Schwabe cycle length for the monthly sunspot group numbers in the two hemispheres obviously differs from each other, and the mean Schwabe cycle length of the monthly sunspot group numbers in the northern hemisphere is slightly larger than that in the southern one; （3） the monthly sunspot group numbers in the northern hemisphere precede those in the southern hemisphere during the years of about 1874- 1927, after which, the southern hemisphere leads the northern hemisphere in the years 1928-1964, and then the northern hemisphere leads in time till the present.

The relationships of solar flares with both sunspot and geomagnetic activity

The relationships between solar flare parameters （total importance, time duration, flare index, and flux） and sunspot activity （R z ） as well as those between geomagnetic activity （aa index） and the flare parameters can be well described by an integral response model with the response time scales of about 8 and 13 months, respectively. Compared with linear relationships, the correlation coefficients of the flare parameters with R z , of aa with the flare parameters, and of aa with R z based on this model have increased about 6%, 17%, and 47% on average, respectively. The time delays between the flare parameters with respect to R z , aa to the flare parameters, and aa to R z at their peaks in a solar cycle can be predicted in part by this model （82%, 47%, and 78%, respectively）. These results may be further improved when using a cosine filter with a wider window. It implies that solar flares are related to the accumulation of solar magnetic energy in the past through a time decay factor. The above results may help us to understand the mechanism of solar flares and to improve the prediction of the solar flares.

Analysis of sudden variations in photospheric magnetic fields during a large flare and their influences in the solar atmosphere

The solar active region NOAA 11719 produced a large two-ribbon flare on 2013 April 11. We have investigated sudden variations in the photospheric magnetic fields in this active region during the flare by employing magnetograms obtained in the spectral line Fe I 6173 A acquired by the Helioseismic and Magnetic Imager （HMI） onboard the Solar Dynamics Observatory （SDO） spacecraft. The analysis of the line-of-sight magnetograms from HMI show sudden and persistent magnetic field changes at different locations of the active region before the onset of the flare and during the flare. The vector magnetic field observations available from HMI also show coincident variations in the total magnetic field strength and its inclination angle at these locations. Using the simultaneous Dopplergrams obtained from HMI, we observe perturbations in the photospheric Doppler signals following the sudden changes in the magnetic fields in the aforementioned locations. The power spectrum analysis of these velocity signals shows enhanced acoustic power in these affected locations during the flare as compared to the pre-flare condition. Accompanying these observations, we have also used nearly simultaneous chromospheric observations obtained in the spectral line Ha 6562.8 A by the Global Oscillation Network Group （GONG） to study the evolution of flare- ribbons and intensity oscillations in this active region. The Ha intensity oscillations also show enhanced oscillatory power during the flare in the aforementioned locations. These results indicate that the transient Lorentz force associated with sudden changes in the magnetic fields could drive localized photospheric and chromospheric oscillations, like the flare-induced oscillations in the solar atmosphere.

Acoustic streaming and Sun's meridional circulation

A vast number of physical processes involving oscillations of a bounded viscous fluid are relevantly influenced by acoustic streaming. When this happens a steady circulation of fluid develops in a thin boundary adjacent to the interface. Some examples are refracted sound waves, a fluid inside a spherical cavity undergoing torsional oscillations or a pulsating liquid droplet. Steady streaming around circular interfaces consists of a hemispherically symmetric recirculation of fluid from the equatorial plane to the polar axes closely resembling the meridional circulation pattern observed in the Sun＇s convection zone that determines the solar cycle. In this paper, it is argued that the acoustic pulsations exhibited by the Sun would lead to acoustic streaming in the boundary of the convection zone. A simple estimation using a typical dominant frequency of 3 mHz and the observed surface oscillation amplitude yields a steady streaming velocity Us - 10 m s^-1, which is on the order of the meridional circulation velocity observed in the Sun＇s convection zone.

Length Scale of Photospheric Granules in Solar Active Regions

Investigating the length scales of granules could help understand the dynamics of granules in the photosphere.In this work,we detected and identified granules in an active region near disk center observed at wavelength of TiO(7057?)by the 1.6 m Goode Solar Telescope(GST).By a detailed analysis of the size distribution and flatness of granules,we found a critical size that divides the granules in motions into two regimes:convection and turbulence.The length scales of granules with sizes larger than 600 km follow Gauss function and demonstrate“flat”in flatness,which reveal that these granules are dominated by convection.Those with sizes smaller than 600 km follow power-law function and behave power-law tendency in flatness,which indicate that the small granules are dominated by turbulence.Hence,for the granules in active regions,they are originally convective in large length scale,and directly become turbulent once their sizes turn to small,likely below the critical size of 600 km.Comparing with the granules in quiet regions,they evolve with the absence of the mixing motions of convection and turbulence.Such a difference is probably caused by the interaction between fluid motions and strong magnetic fields in active regions.The strong magnetic fields make high magnetic pressure which creates pressure walls and slows down the evolution of convective granules.Such walls cause convective granules extending to smaller sizes on one hand,and cause wide intergranular lanes on the other hand.The small granules isolated in such wide intergranular lanes are continually sheared,rotated by strong downflows in surroundings and hereby become turbulent.

Chinese Sunspot Drawings and Their Digitization—(Ⅶ) Sunspot Penumbra to Umbra Area Ratio Using the Hand-Drawing Records from Yunnan Observatories

The ratio of penumbral to umbra area of sunspots plays a crucial role in the solar physics fields, especially for understanding the origin and evolution of the solar activity cycle. By analyzing the recently digitized sunspot drawings observed from Yunnan Observatories(1957–2021), we investigate the long-term variation of the penumbral to umbra area ratio of sunspots. An automatic extraction method, based on the maximum between-class variance and the morphological discrimination, is used to accurately extract penumbra and umbra and to calculate the ratio over six solar cycles(cycle 19–24). The expected value of the ratio of penumbra to umbra area is found to be 6.63 ± 0.98, and it does not exhibit any systematic variation with sunspot latitudes and phases. The average ratio fluctuates from 5 to 7.5 per year and the overall trend has decreased after 1999 compared to the previous one.The ratio of sunspot penumbra to umbra area satisfies the log-normal distribution, implying that its variation is related to the evolution of the photospheric magnetic field. Our results are consistent with previous works.

Signature of high-order azimuthal MHD body modes in sunspot's low atmosphere

The five-minute oscillations inside sunspots appear to be the absorption of the solar p-mode. It is a potential tool to probe a sunspot＇s sub-structure. We studied the collective properties of five-minute oscillations in the power and phase distribution at the sunspot＇s umbra-penumbra boundary. The azimuthal distributions of the power and phase of five-minute oscillations enclosing a sunspot＇s umbra were obtained with images taken with the Solar Dynamics Observatory/Atmospheric Imaging Assembly （SDO/AIA）. The azimuthal modes were quantified with periodogram analysis and justified with significance tests. The azimuthal nodal structures in an approximately ax- ially symmetric sunspot AR 11131 （2010 Dec 08） were investigated. Mode numbers ra = 2, 3, 4, 7, 10 were obtained in both 1700 A and 1600A bandpasses. The 1600A channel also revealed an extra mode at m = 9. In the upper atmosphere （304 A）, fewer modes were detected at m = 3, 4, 7. The azimuthal modes in the sunspot＇s low atmo- sphere could be interpreted as high-order azimuthal MHD body modes. They were detected in the power and phase of the five-minute oscillations in sunspot AR 11131 with SDO/AIA data. Fewer modes were detected in the sunspot＇s upper atmosphere.

Magnetic cycles of Sun-like stars with different levels of coronal and chromospheric activity——comparison with the Sun

The atmospheric activity of the Sun and Sun-like stars is analyzed involving observations from the HK-project at the Mount Wilson Observatory, the California and Carnegie Planet Search Program at the Keck and Lick Observatories and the Magellan Planet Search Program at the Las Campanas Observatory. We show that for stars of E G and K spectral classes, the cyclic activity, similar to the 11-yr solar cycle, is different： it becomes more prominent in K-stars. Comparative study of Sun-like stars with different levels of chromospheric and coronal activity confirms that the Sun belongs to stars with a low level of chromospheric activity and stands apart among these stars by its minimum level of coronal radiation and minimum level of variations in photospheric flux.