Development of a Code to Analyze the Solar White-Light Images from the Kodaikanal Observatory: Detection of Sunspots,Computation of Heliographic Coordinates and Area

Sunspots are the most conspicuous aspects of the Sun. They have a lower temperature, as compared to the surrounding photosphere; hence, sunspots appear as dark regions on a brighter background. Sunspots cyclically appear and disappear with a 11-year periodicity and are associated with a strong magnetic field ( ～10³ G) structure. Sunspots consist of a dark umbra, surrounded by a lighter penumbra. Study of umbra–penumbra area ratio can be used to give a rough idea as to how the convective energy of the Sun is transported from the interior, as the sunspot’s thermal structure is related to this convective medium.An algorithm to extract sunspots from the white-light solar images obtained from the Kodaikanal Observatory is proposed. This algorithm computes the radius and center of the solar disk uniquely and removes the limb darkening from the image. It also separates the umbra and computes the position as well as the area of the sunspots. The estimated results are compared with the Debrecen photoheliographic results. It is shown that both area and position measurements are in quite good agreement.     

The Temperature – Magnetic Field Relation in Observed and Simulated Sunspots

Observations of the relation between continuum intensity and magnetic field strength in sunspots have been made for nearly five decades. This work presents full-Stokes measurements of the full-split (\(g = 3\)) line Fe i 1564.85 nm with a spatial resolution of \(0.5^{\prime\prime}\) obtained with the GREGOR Infrared Spectrograph in three large sunspots. The continuum intensity is corrected for instrumental scattered light, and the brightness temperature is calculated. Magnetic field strength and inclination are derived directly from the line split and the ratio of Stokes components. The continuum intensity (temperature) relations to the field strength are studied separately in the umbra, light bridges, and penumbra. The results are consistent with previous studies, and it was found that the scatter of values in the relations increases with increasing spatial resolution thanks to resolved fine structures. The observed relations show trends common for the umbra, light bridges, and the inner penumbra, while the outer penumbra has a weaker magnetic field than the inner penumbra at equal continuum intensities. This fact can be interpreted in terms of the interlocking comb magnetic structure of the penumbra. A comparison with data obtained from numerical simulations was made. The simulated data generally have a stronger magnetic field and a weaker continuum intensity than the observations, which may be explained by stray light and limited spatial resolution of the observations, and also by photometric inaccuracies of the simulations.     

The Pattern of Moving Magnetic Inhomogeneities in and Around Sunspots

High-resolution MDI magnetograms are used to study the pattern of moving magnetic inhomogeneities in sunspots. We examine the inward and outward moving features in sunspots. The velocity of these features is small in the umbra while it is about 0.5 km s^–1 in the outer penumbra. The inward and outward moving features may be the possible origin for the long-term fluctuations of magnetic field strength in sunspots.     

Fractal properties of solar magnetic fields

We study the spatial properties of solar magnetic fields using data from the Solar Vector Magnetograph of the Marshall Space Flight Center (MSFC) (FeI 5250.2 Å) and SOHO/MDI longitudinal magnetic field measurements (Ni 6767.8 Å) (96-min full-disk maps). Our study is focused on two objects: the fractal properties of sunspots and the fractal properties of the spatial magnetic field distribution of active and quiet regions considered as global structures. To study the spatial structure of sunspots, we use a well-known method of determining the fractal dimension based on an analysis of the perimeter—area relation. To analyze the fractal properties of the spatial magnetic field distribution over the solar surface, we use a technique developed by Higuchi. We have revealed the existence of three families of self-similar contours corresponding to the sunspot umbra, penumbra, and adjacent photosphere. The fractal coefficient has maxima near the umbra—penumbra and penumbra—photosphere boundaries. The fractal dependences of the longitudinal and transverse magnetic field distributions are similar, but the fractal numbers themselves for the transverse fields are larger than those for the longitudinal fields approximately by a factor of 1.5. The fractal numbers decrease with increasing mean magnetic field strength, implying that the magnetic field distribution is more regular in active regions.     

The relation between the intensities of umbra and penumbra of large sunspots

Based on pinhole photometer observations in 4–6 wavelength regions we have searched for a connection between the intensities of the umbra and the penumbra of sunspots. For the 1.67 μm wavelength region it is apparent that spots with dark umbrae also have dark penumbrae. In the other wavelength regions similar relations are found. The darkness of the spot is probably connected with the degree of complexity of the spot.     

Width of emission cores of the line K Ca ii in sunspots

Emission core widths of K Ca ii line in the umbra and penumbra of 9 sunspots and in their vicinity are measured. All sunspots are located near the solar disc center. Data on variation of widths W _K along the mean sunspot radius are obtained. Values W _K in the umbra and penumbra centers are equal or somewhat less than on stars of the same luminosity.     

A Curious History of Sunspot Penumbrae: An Update

The ratio of penumbral to umbral area of sunspots is an important topic for solar and geophysical studies. Hathaway (Solar Phys.286, 347, 2013) found a curious behaviour in this parameter for small sunspot groups (areas smaller than 100 millionths of solar hemisphere, msh) using records from Royal Greenwich Observatory (RGO). Hathaway showed that the penumbra–umbra ratio decreased smoothly from more than 7 in 1905 to lower than 3 by 1930 and then increased to almost 8 in 1961. Thus, Hathaway proposed the existence of a secular variation in the penumbra–umbra area ratio. In order to confirm that secular variation, we employ data of the sunspot catalogue published by the Coimbra Astronomical Observatory (COI) for the period 1929?–?1941. Our results disagree with the penumbra–umbra ratio found by Hathaway for that period. However, the behaviour of this ratio for large (areas greater or equal than 100 msh) and small groups registered in COI during 1929?–?1941 is similar to data available from RGO for the periods 1874?–?1914 and 1950?–?1976. Nevertheless, while the average values and time evolution of the ratio in large groups are similar those for small groups according to the Coimbra data (1929?–?1941) it is not analogous for RGO data for the same period. We also found that the behaviour of the penumbra–umbra area ratio for smaller groups in both observatories is significantly different. The main difference between the area measurements made in Coimbra and RGO is associated with the umbra measurements. We would like to stress that the two observatories used different methods of observation and while in COI both methodology and instruments did not change during the study period, some changes were carried out in RGO that could have affected measurements of umbra and penumbra. These facts illustrate the importance of the careful recovery of past solar data.     

Broad-band circular polarimetry of sunspots, 0.4–1.7 microns: Spatial scans with a 3.4 arc sec diameter aperture

Spatial scans with a resolution of 3.4 arc sec of the broad-band circular polarization of several sunspots have been made in five filter bands over the wavelength range 0.4–1.7µ with a sensitivity of 1 × 10^–6 fractional polarization. The scans, across a spot through the penumbra and umbra center, revealed two important features: (1) The broad-band circularly polarized fluxV reverses in sign, or diminishes to near zero, at the center of the umbral region relative to the outer penumbra. This effect was wavelength dependent and was most clearly detected as a definite reversal in a band at 1.2µ, although a reversal was also detected in a very broad band extending from 0.8 to 1.6µ. (2) There is a marked asymmetry: in all cases the limbward penumbral region exhibited strongerV values than did the disk-center (inward) side of the spot, at all observed wavelengths. Such previously unreported structure in the magnetic circular polarization of sunspots provides new clues for understanding the anomalous large broad-band polarization at short wavelengths and at the same time imposes new constraints on sunspot models. For example, the polarization reversal in the umbra relative to the penumbra can be naively explained by return-flux sunspot models; but this is not the only interpretation. Alternatively, it can relate to reversals in mass-flow velocities and/or vertical velocity gradients, as between the umbra and penumbra.     

Observational Evidence for Rising Penumbral Flux Tubes?

On 13 May 2000 parts of a penumbra were observed in an active region NOAA 8990 with the La Palma Stokes Polarimeter attached to the Swedish Vacuum Solar Telescope. The stratification over the solar atmosphere of different physical parameters is retrieved from these data by using the Stokes inversion based on response functions. The results confirm the previous findings of the penumbral structure. In general, the magnetic field becomes weaker and more horizontal with increasing distance from the umbra and the line-of-sight velocities are increasing towards the outer boundary of the penumbra. The results also suggest the existence of the unresolved fine structure of the penumbra. The stratifications of the temperature and of the magnetic field strength indicate the presence of rising flux tubes, which were predicted theoretically by Schlichenmaier, Jahn and Schmidt (1998, Astron. Astrophys. 337, 897).     

Automated Detection, Characterization, and Tracking of Sunspots from SoHO/MDI Continuum Images

We describe a procedure for automated detection of sunspots from SoHO/MDI full-disk continuum images. The MDI Level-1.8 continuum images were first corrected for the limb darkening and stray light, and then were flat-fielded. Sunspots were extracted using a newly developed automated sunspot detection procedure, which is based on the level set, namely the selective binary and Gaussian function regularized level set (SBGFRLS) method (Zhang et al., Image Vis. Comput. 28, 668, 2010). In this method we initialize a two-dimensional level-set function and evolve it using a signed pressure force (SPF) function. For sunspot detection, the level-set function was defined twice, first for umbra and then for penumbra extraction. Using this procedure, along with the characterization of detected sunspots we have also generated tracking reports of all sunspots in a fully unsupervised manner.     

Motions and lifetimes of the penumbral bright grains in sunspots

It is confirmed that the penumbral bright grains are moving towards the sunspots umbra. We find different proper motions of 0.08 to 0.33 km s^–1 for different penumbrae and different reduction methods. The lifetimes of these bright grains are about 1.5 to 3 hr depending on the position in the penumbra.     

Bimodal Distribution of Magnetic Fields and Areas of Sunspots

We applied automatic identification of sunspot umbrae and penumbrae to daily observations from the Helioseismic Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) to study their magnetic flux density (B) and area (A). The results confirm an already known logarithmic relationship between the area of sunspots and their maximum flux density. In addition, we find that the relation between average magnetic flux density ( $B_{\rm avg}$ ) and sunspot area shows a bimodal distribution: for small sunspots and pores (A≤20 millionth of solar hemisphere, MSH), $B_{\rm avg} \approx 800~\mbox{G}$ (gauss), and for large sunspots (A≥100 MSH), $B_{\rm avg}$ is about 600 G. For intermediate sunspots, average flux density linearly decreases from about 800 G to 600 G. A similar bimodal distribution was found in several other integral parameters of sunspots. We show that this bimodality can be related to different stages of sunspot penumbra formation and can be explained by the difference in average inclination of magnetic fields at the periphery of small and large sunspots.     

Investigation of Umbral Dots with the <Emphasis Type="Italic">New Vacuum Solar Telescope</Emphasis>

Umbral dots (UDs) are small isolated brightenings observed in sunspot umbrae. They are convective phenomena existing inside umbrae. UDs are usually divided into central UDs (CUDs) and peripheral UDs (PUDs) according to their positions inside an umbra. Our purpose is to investigate UD properties and analyze their relationships, and further to find whether or not the properties depend on umbral magnetic field strengths. Thus, we selected high-resolution TiO images of four active regions (ARs) taken under the best seeing conditions with the New Vacuum Solar Telescope in the Fuxian Solar Observatory of the Yunnan Astronomical Observatory, China. The four ARs (NOAA 11598, 11801, 12158, and 12178) include six sunspots. A total of 1220 CUDs and 603 PUDs were identified. Meanwhile, the radial component of the vector magnetic field of the sunspots taken with the Helioseismic and Magnetic Imager on-board the Solar Dynamics Observatory was used to analyze relationships between UD properties and umbral magnetic field strengths. We find that diameters and lifetimes of UDs exhibit an increasing trend with the brightness, but velocities do not. Moreover, diameters, intensities, lifetimes and velocities depend on the surrounding magnetic field. A CUD diameter was found larger, the CUD brighter, its lifetime longer, and its motion slower in a weak umbral magnetic field environment than in a strong one.     

Temporal fluctuations of the magnetic field in sunspots

The magnetic field strength in sunspots was derived from time series of two-dimensional spectra taken with the Göttingen 2D-spectrometer at the Vacuum Tower Telescope on Tenerife in August 1997. For the present measurements the magnetically sensitive line Fe?i 684.3 nm was selected. The main spot of the investigated sunspot group has a maximum magnetic field strength of 2270 G. Enhanced power of the magnetic field variations was found at the boundary between umbra and penumbra for all frequency ranges. These fluctuations are not well correlated with those of intensity variations or line shifts. Other spatial power peaks occur in a dark patch inside the centreside penumbra and at the centres of some accompanying small spots. Since no clear peaks at certain frequencies are found, the variations are not harmonic oscillations. A possible relation to Hα flares is investigated. There are several cases of published observations of magnetic field variations where flares occurred soon after the measurements, but very little before. Therefore it is not very probable that flares act as exciters of magnetic field variations.

     

A comparison of the oscillations in sunspot umbrae,penumbrae, and the surrounding photosphere

Time series of the nonsplit Fei 7090 Å line have been observed in several sunspots with a 100 x 100 diode array corresponding to 48 arc sec times 1.39 Å. The spatial behaviour of Doppler motions along one fixed slit position has been studied as a function of time. Former results are confirmed, that the power in the five minute range decreases from the photosphere to the umbra, where, however, values still well above the noise level are measured. Regarding the penumbra, the power tends to exhibit a maximum at locations where the line-of-sight component of a radial horizontal field should be maximal. This indicates that the direction of the oscillatory velocities might be influenced by the magnetic field or the Evershed flow. No significant power is found in the 3 min range. An exception might be seen in a small patch at the limb of the umbra of one spot.     

Flows around sunspots and pores

We report on three sequences of high-resolution white-light and magnetogram observations obtained in the summer of 1989. The duration of sub-arcsecond seeing was three to four hours on each day. Study of the white-light and magnetogram data yields the following results:

1. For all but one of the sunspots we have observed, both dark fibrils and bright grains in the inner part of the penumbra of sunspots move toward the umbra with a speed of about 0.5 km s^-1. In the outer part of the penumbra, movement is away from the umbra. The one exception is a newly formed spot, which has inflow only in its penumbra.
2. Granular flows converge toward almost every pore, even before its formation. Pores are observed to form by the concentration of magnetic flux already existing in the photosphere. The pores (or small sunspots), in turn, then move and concentrate to form bigger sunspot.
3. We followed an emerging flux region (EFR) from 29 to 31 July, 1989 that was composed of a large number of bipoles with magnetic polarities mixed over a large area in the first day of its birth. As time went on, polarities sorted out: the leading polarity elements moved in one direction; the following, the opposite. During the process a large number of cancellations occurred, with some sub-flares and surges observed simultaneously. After about 24 hours, the positive and negative fluxes were essentially separated.
4. We find two kinds of photospheric dark alignments in the region of new flux emergence: (a) alignments connecting two poles of opposite magnetic polarity form the tops of rising flux tubes; (b) alignments corresponding to the magnetic flux of one polarity, which we call elongated pores.

     

A Multi-wavelength Analysis of Active Regions and Sunspots by Comparison of Automatic Detection Algorithms

Since the Solar Dynamics Observatory (SDO) began recording ≈?1 TB of data per day, there has been an increased need to automatically extract features and events for further analysis. Here we compare the overall detection performance, correlations between extracted properties, and usability for feature tracking of four solar feature-detection algorithms: the Solar Monitor Active Region Tracker (SMART) detects active regions in line-of-sight magnetograms; the Automated Solar Activity Prediction code (ASAP) detects sunspots and pores in white-light continuum images; the Sunspot Tracking And Recognition Algorithm (STARA) detects sunspots in white-light continuum images; the Spatial Possibilistic Clustering Algorithm (SPoCA) automatically segments solar EUV images into active regions (AR), coronal holes (CH), and quiet Sun (QS). One month of data from the Solar and Heliospheric Observatory (SOHO)/Michelson Doppler Imager (MDI) and SOHO/Extreme Ultraviolet Imaging Telescope (EIT) instruments during 12 May?–?23 June 2003 is analysed. The overall detection performance of each algorithm is benchmarked against National Oceanic and Atmospheric Administration (NOAA) and Solar Influences Data Analysis Center (SIDC) catalogues using various feature properties such as total sunspot area, which shows good agreement, and the number of features detected, which shows poor agreement. Principal Component Analysis indicates a clear distinction between photospheric properties, which are highly correlated to the first component and account for 52.86% of variability in the data set, and coronal properties, which are moderately correlated to both the first and second principal components. Finally, case studies of NOAA 10377 and 10365 are conducted to determine algorithm stability for tracking the evolution of individual features. We find that magnetic flux and total sunspot area are the best indicators of active-region emergence. Additionally, for NOAA 10365, it is shown that the onset of flaring occurs during both periods of magnetic-flux emergence and complexity development.     

Statistical Properties Of Sunspots In 1996–2004: I. Detection, North–South Asymmetry And Area Distribution

The first statistical results in sunspot distributions in 1996–2004 obtained from the Solar Feature Catalogues (SFC) are presented. A novel robust technique is developed for automated identification of sunspots on SOHO/MDI white-light (WL) full-disk solar images. The technique applies image standardization procedures for elimination of the limb darkening and non-circular image shape, uses edge-detection methods to find the sunspot candidates and their edges and morphological operations to smooth the features and fill in gaps. The detected sunspots are verified with the SOHO/MDI magnetograms by strong magnetic fields being present in sunspots. A number of physical and geometrical parameters of the detected sunspot features are extracted and stored in the relational SFC database including umbra/penumbra masks in the form of run-length data encoding of sunspot bounding rectangles. The detection results are verified by comparison with the manual daily detection results in Meudon and Locarno Observatories in 2002 and by correlation (about 96%) with the 4 year sunspot areas produced manually at NOAA. Using the SFC data, sunspot area distributions are presented in different phases of the solar cycle and hemispheres which reveals a periodicity of the north–south asymmetry with a period of about 7–8 years. The number of sunspots increases exponentially with the area decrease with the index slightly increasing from −1.15 (1997) to −1.34 (2001).     

NOAA6361活动区现象分析

本文分析了ＮＯＡＡ６３６１活动区中的一些现象，发现该活动区在衰亡阶段经历了两次同极性黑子的复合过程，复合后的黑子本影间均有光桥存在，观测结果倾向于支持Ｐａｒｋｅｒ１９７９年提出的黑子多磁流管模型。１４日复合后的黑子本影还顺时针方向旋转了约７０度角，从半影纤维的同样顺时针旋转可以认为：该黑子的半影磁场并非是普遍认为的简单的本影磁场的发散部分。我们还观测到另外两个比较有趣的现象：①δ黑子中的ｐ极性黑子     

The depth of sunspots

Good quality photographs of many regular sunspots were obtained in three spectral regions, with the use of narrow band filters. Isodensity contours were used for measurements of the umbra and penumbra size along different axes of sunspots. A parameter based on the morphology of the whole sunspot was defined for a better study of the Wilson effect. Some interesting results are the following: (a) an east-west asymmetry of the Wilson effect was clearly observed, (b) differences in the phenomenon in the three spectral regions were barely significant compared with observational errors and (c) measurements indicate that larger sunspots have greater depth.