Automatic Detection of Sunspots and Extractionof Their Feature Parameters

Sunspots are solar features located in active regions of the Sun, whose number is an indicator of the Sun's magnetic activity. With a substantial increase in the quantity of solar image data, the automated detection and verification of various solar features have become increasingly important for the accurate and timely forecasts of solar activity and space weather. In order to use the high time-cadence SDO/HMI data to extract the main sunspot features for forecasting solar activities, we have established an automatic detection method of sunspots based on mathematical morphology, and calculated the sunspot group area and sunspot number. By comparing our results with those obtained from the Solar Region Summary compiled by NOAA/SWPC, it is found that the sunspot group areas and sunspot numbers computed with our algorithm are in good agreement with the active region values released by SWPC, and the corresponding correlation coefficients for the sunspot group area and sunspot number are 0.77 and 0.79, respectively. By using the method of this paper, the high time-cadence feature parameters can be obtained from the HMI data to provide the timely and accurate inputs for the solar activity forecast.     

A new interpretation of Christian Horrebow''s sunspot observations from 1761 to 1777

Christian Horrebow and his colleagues of Copenhagen, Denmark, actively observed sunspots from 1761 to 1777. These observations were examined by Thiele in 1859 and by d'Arrest in 1873 with markedly different conclusions. Thiele reported nearly twice as many sunspot groups as d'Arrest. To resolve this discrepancy, we have reexamined Horrebow's original notebooks. We find slightly more sunspot groups then did d'Arrest. Thiele apparently called individual sunspots sunspot groups, so he would call a bipolar group two groups. d'Arrest seems to have missed counting some of the smaller sunspot groups. A correct interpretation of Horrebow's observations is required in efforts to reconstruct solar activity. Wolf gave a sunspot number for 1769 of 106.1. On the basis of our re-examination of Horrebow's drawings and other observers, we deduce a sunspot number of about 80.5 for 1769.     

Automatic detection of sunspots on solar continuum HMI images blending local–global threshold

A huge collection of solar images to visualize sunspot are acquired by various solar observatories spread across the globe. This necessitates efficient tools for detecting and analyzing the sunspots encompassing diverse solar features. One such contribution is delivered in this work by exploiting the intrinsic intensity variations of solar images associated with sunspots and their attributes. The presented mechanism initially, pre-processes the acquired solar images by correcting the intensity variations introduced while profiling from the sun center to the limb followed by smoothening using a localized window. The resultant is then differenced from the global threshold that is obtained as a result of the statistical analysis computed over the probability distribution function of the input image. This arrangement offers higher discerning variations concerned with the local contextual structures related to sunspot, umbra, and penumbra. Also, it captures the major gradient variation between these regions that adds to the pixel heterogeneity surrounding them to finally render an automatic sunspot detection mechanism distinguishing the diverse solar regions. Receiver Operating Characteristics (ROC) investigation on annual solar images in Flexible Image Transport System (FITS) format reveals the presented method’s efficacy. Also, Pearson correlation analysis of the evaluated sunspot numbers from the detected sunspots with the solar catalog reveals the scheme’s detection closeness. Moreover, the model’s simplicity analyzed along the time and space dimensions affirms its extension to real-time analysis     

A new look at wolf sunspot numbers in the late 1700's

Long-term homogeneous observations of solar activity or many solar cycles are essential for investigating many problems in solar physics and climatology. The one key parameter used in most long-term studies is the Wolf sunspot number, which is susceptible to observer bias, particularly because it is highly sensitive to the observer's ability to see the smallest sunspots. In this paper we show how the Wolf sunspot number can be derived from the number of sunspot groups alone. We utilize this approach to obtain a Group Wolf number. This technique has advantages over the classical method of determining the Wolf number because corrections for observer differences are reduced and long-term self-consistent time series can be developed. The level of activity can be calculated to an accuracy of ± 5% using this method. Applying the technique to Christian Horrebow's observations of solar cycles 1, 2, and 3 (1761–1777), we find that the standard Wolf numbers are nearly homogeneous with sunspot numbers measured from 1875 to 1976 except the peak of solar cycle 2 is too low by 30%. This result suggests that further analyses of early sunspot observations could lead to significant improvements in the uniformity of the measurements of solar activity. Such improvements could have important impacts upon our understanding of long-term variations in solar activity, such as the Gleissberg cycle, or secular variations in the Earth's climate.     

Sunspot groups at high latitude

Data of sunspot groups at high latitude (35°), from the year 1874 to the present (2000 January), are collected to show their evolutional behaviour and to investigate features of the yearly number of sunspot groups at high latitude. Subsequently, an evolutional pattern of sunspot group number at high latitude is given in this paper. Results obtained show that the number of sunspot groups of a solar cycle at high latitude rises to a maximum value about 1 yr earlier than the time of the maximum of sunspot relative numbers of the solar cycle, and then falls to zero more rapidly. The results also show that, at the moment, solar activity described by the sunspot relative numbers has not yet reached its minimum. In general, sunspot groups at high latitude have not appeared on the solar disc during the last 3 yr of a Wolf solar cycle. The asymmetry of the high latitude sunspot group number of a Wolf solar cycle can reflect the asymmetry of solar activity in the Wolf solar cycle, and it is suggested that one could further use the high latitude sunspot group number during the rising time of a Wolf solar cycle, maximum year included, to judge the asymmetry of solar activity over the whole solar cycle.     

Solar Activity in the Past: From Different Proxies to Combined Reconstruction

Using proxy series that includes ancient observations of sunspots and auroras, concentrations of cosmogenic isotopes ¹⁴C and ¹⁰Be, we reconstruct sunspot activity level since 850 AD to the present. As a main reference index of solar activity we use the Wolf sunspot numbers, which, as we demonstrate, reflect true levels of the activity in the 18th and 19th centuries better than the Group sunspot numbers. We construct a set of linear and nonlinear inductive models, which are in a good agreement with each other and reproduce the known global solar activity extrema. According to our results, amplitudes of the global maxima are in intermediate positions in respect to estimations of other authors. It follows from our reconstructions that the global maximum of the 20th century is not much higher, if at all, than others.     

Survey and Merging of Sunspot Catalogs

In view of the construction of new sunspot-based activity indices and proxies, we conducted a comprehensive survey of all existing catalogs providing detailed parameters of photospheric features over long time intervals. Although there are a fair number of such catalogs, a global evaluation showed that they suffer from multiple limitations: finite or fragmented time coverage, limited temporal overlap between catalogs, and, more importantly, a mismatch in contents and conventions. Starting from the existing material, we demonstrate how the information from parallel catalogs can be merged to form a much more comprehensive record of sunspots and sunspot groups. To do this, we use the uniquely detailed Debrecen Photoheliographic Data (DPD), which is already a composite of several ground-based observatories and of SOHO data, and the USAF/Mount Wilson catalog from the Solar Observing Optical Network (SOON). We also outline our cross-identification method, which was needed to match the non-overlapping solar active-region nomenclature. This proved to be the most critical and subtle step when working with multiple catalogs. This effort, focused here first on the last two solar cycles, should lead to a better central database that collects all available sunspot group parameters to address future solar-cycle studies beyond the traditional sunspot-index time series [R _i].     

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).     

What do the solar activity indices represent?

Sunspot number, sunspot area, and radio flux at 10.7 cm are the indices which are most frequently used to describe the long‐term solar activity. The data of the daily solar full‐disk magnetograms measured at Mount Wilson Observatory from 19 January 1970 to 31 December 2012 are utilized together with the daily observations of the three indices to probe the relationship of the full‐disk magnetic activity respectively with the indices. Cross correlation analyses of the daily magnetic field measurements at Mount Wilson observatory are taken with the daily observations of the three indices, and the statistical significance of the difference of the obtained correlation coefficients is investigated. The following results are obtained: (1) The sunspot number should be preferred to represent/reflect the full‐disk magnetic activity of the Sun to which the weak magnetic fields (outside of sunspots) mainly contribute, the sunspot area should be recommended to represent the strong magnetic activity of the Sun (in sunspots), and the 10.7 cm radio flux should be preferred to represent the full‐disk magnetic activity of the Sun (both the weak and strong magnetic fields) to which the weak magnetic fields mainly contribute. (2) On the other hand, the most recommendable index that could be used to represent/reflect the weak magnetic activity is the 10.7 cm radio flux, the most recommendable index that could be used to represent the strong magnetic activity is the sunspot area, and the most recommendable index that could be used to represent the full‐disk magnetic activity of the Sun is the 10.7cm radio flux. Additionally, the cycle characteristics of the magnetic field strengths on the solar disk are given. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Group Sunspot Numbers: A New Solar Activity Reconstruction

In this paper, we construct a time series known as the Group Sunspot Number. The Group Sunspot Number is designed to be more internally self-consistent (i.e., less dependent upon seeing the tiniest spots) and less noisy than the Wolf Sunspot Number. It uses the number of sunspot groups observed, rather than groups and individual sunspots. Daily, monthly, and yearly means are derived from 1610 to the present. The Group Sunspot Numbers use 65941 observations from 117 observers active before 1874 that were not used by Wolf in constructing his time series. Hence, we have calculated daily values of solar activity on 111358 days for 1610–1995, compared to 66168 days for the Wolf Sunspot Numbers. The Group Sunspot Numbers also have estimates of their random and systematic errors tabulated. The generation and preliminary analysis of the Group Sunspot Numbers allow us to make several conclusions: (1) Solar activity before 1882 is lower than generally assumed and consequently solar activity in the last few decades is higher than it has been for several centuries. (2) There was a solar activity peak in 1801 and not 1805 so there is no long anomalous cycle of 17 years as reported in the Wolf Sunspot Numbers. The longest cycle now lasts no more than 15 years. (3) The Wolf Sunspot Numbers have many inhomogeneities in them arising from observer noise and this noise affects the daily, monthly, and yearly means. The Group Sunspot Numbers also have observer noise, but it is considerably less than the noise in the Wolf Sunspot Numbers. The Group Sunspot Number is designed to be similar to the Wolf Sunspot Number, but, even if both indices had perfect inputs, some differences are expected, primarily in the daily values.     

Overlooked sunspot observations by Hevelius in the early Maunder Minimum, 1653–1684

In the bookMachina Coelestis (1679), Johannes Hevelius lists his daily solar observations from 1653 to 1679. He mentions 19 sunspot groups during this interval, of which 14 are unique to Hevelius and five are confirmed by other observers. There are an additional 9 sunspot groups during this interval that were not observed by Hevelius. In five cases he was not observing, but in the other four cases he did observe but failed to comment upon sunspots. The spots he missed or failed to observe tend to occur near the end of his career. This suggests Hevelius occasionally missed sunspots but usually was a reliable observer. These observations are important because they provide us the only known daily listing of solar observations during the early years of the Maunder Minimum. They are also important because they were overlooked by Wolf, Spoerer, Maunder, Eddy, and others in their study of solar activity in the seventeenth century. They provide us the best record of the sunspot maximum of 1660 when one sunspot lasted at least 86 days as it traversed the solar disk four times. The same region was active for seven solar rotations.     

Group Sunspot Numbers: A New Solar Activity Reconstruction

In this paper, we construct a time series known as the Group Sunspot Number. The Group Sunspot Number is designed to be more internally self-consistent (i.e., less dependent upon seeing the tiniest spots) and less noisy than the Wolf Sunspot Number. It uses the number of sunspot groups observed, rather than groups and individual sunspots. Daily, monthly, and yearly means are derived from 1610 to the present. The Group Sunspot Numbers use 65941 observations from 117 observers active before 1874 that were not used by Wolf in constructing his time series. Hence, we have calculated daily values of solar activity on 111358 days for 1610–1995, compared to 66168 days for the Wolf Sunspot Numbers. The Group Sunspot Numbers also have estimates of their random and systematic errors tabulated. The generation and preliminary analysis of the Group Sunspot Numbers allow us to make several conclusions: (1) Solar activity before 1882 is lower than generally assumed and consequently solar activity in the last few decades is higher than it has been for several centuries. (2) There was a solar activity peak in 1801 and not 1805 so there is no long anomalous cycle of 17 years as reported in the Wolf Sunspot Numbers. The longest cycle now lasts no more than 15 years. (3) The Wolf Sunspot Numbers have many inhomogeneities in them arising from observer noise and this noise affects the daily, monthly, and yearly means. The Group Sunspot Numbers also have observer noise, but it is considerably less than the noise in the Wolf Sunspot Numbers. The Group Sunspot Number is designed to be similar to the Wolf Sunspot Number, but, even if both indices had perfect inputs, some differences are expected, primarily in the daily values.     

Relative sunspot numbers in the first half of eighteenth century

We revised relative sunspot numbers in the time interval 1700–1748 for which Wolf derived their annual means. The frequency of daily observations, counting simultaneously the number of sunspots and the number of sunspot groups necessary for determinating Wolf's relative sunspot numbers, is in this time interval very low and covers, on average, 4.8% of the number of all days only. There also exist incomplete observations not convenient to determine relative sunspot numbers. To enlarge the number of daily relative sunspot numbers we used the nonlinear, two-step interpolation method derived earlier by Letfus (1996, 1999). After interpolation, the mean value increased to 13.8%. Waldmeier (1968) found that the scaling factor k can be derived directly from the observed number of spots f and from the number of sunspot groups g. From the observations made at Zürich (Wolf and his assistants, Wolfer), at Peckeloh, and at Moncalieri during the years 1861–1928, we derived a new, more correct empirical relation. The resulting annual relative sunspot numbers are given in Table II. However, only for 26 years (53.0%) from the total number of 49 years was it possible to derive annual relative sunspot numbers. The observations were missing for the other years. This corresponds with results of Wolf, which gives the annual relative sunspot numbers for all 49 years. For the years when the data were missing, he marked these values as interpolated or very uncertain ones. Most of the observations originate from two data series (Kirch, Plantade), for which Wolf derived a higher scaling factor (k=2.0) than followed from the newly derived relation (k=1.40). The investigated time interval covers four solar cycles. After our results, the height of the first cycle (No. –4), given by Wolf, should be lowered by about two-thirds, the following two cycles (Nos. –3 and –2) lowered by one-third, as given by Wolf, and only the height of the fourth one (No. –1) should be unchanged. The activity levels of the cycles, as represented by group sunspot numbers, are lower by about one-fourth and, in the case of the first one (No. –4) even by two-thirds of the levels derived by us. The group sunspot numbers, derived from a much greater number of observations, have also greater credibility than other estimates. The shapes of the cycles, as given by Wolf, can be considered only as their more or less idealized form.     

Comparison of the Variations of CMEs and ICMEs with those of other Solar and Interplanetary Parameters During Solar Cycle 23

This paper examines the variations of coronal mass ejections (CMEs) and interplanetary CMEs (ICMEs) during solar cycle 23 and compares these with those of several other indices. During cycle 23, solar and interplanetary parameters had an increase from 1996 (sunspot minimum) to ∼2000, but the interval 1998–2002 had short-term fluctuations. Sunspot numbers had peaks in 1998, 1999, 2000 (largest), 2001 (second largest), and 2002. Other solar indices had matching peaks, but the peak in 2000 was larger than the peak in 2001 only for a few indices, and smaller or equal for other solar indices. The solar open magnetic flux had very different characteristics for different solar latitudes. The high solar latitudes (45^∘–90^∘) in both N and S hemispheres had flux evolutions anti-parallel to sunspot activity. Fluxes in low solar latitudes (0^∘–45^∘) evolved roughly parallel to sunspot activity, but the finer structures (peaks etc. during sunspot maximum years) did not match with sunspot peaks. Also, the low latitude fluxes had considerable N–S asymmetry. For CMEs and ICMEs, there were increases similar to sunspots during 1996–2000, and during 2000–2002, there was good matching of peaks. But the peaks in 2000 and 2001 for CMEs and ICMEs had similar sizes, in contrast to the 2000 peak being greater than the 2001 peak for sunspots. Whereas ICMEs started decreasing from 2001 onwards, CMEs continued to remain high in 2002, probably due to extra contribution from high-latitude prominences, which had no equivalent interplanetary ICMEs or shocks. Cosmic ray intensity had features matching with those of sunspots during 2000–2001, with the 2000 peak (on a reverse scale, actually a cosmic ray decrease or trough) larger than the 2001 peak. However, cosmic ray decreases started with a delay and ended with a delay with respect to sunspot activity.     

The He I 10830 Å Line as an Indicator of Solar Activity

The behavior of the He I 10830 Å infrared triplet parameters in active and quiet solar regions was traced from 1976 until 2000. We analyze the correlation between the central depth of the main He I line component and other solar activity indices: the Wolf number, the radiation flux at a frequency of 2800 MHz, the mean number of flares in sunspot groups, and the mean solar magnetic field. We show that the strong correlation between the He I 10830 Å line depth and the phase of the 11-year solar cycle allows this depth to be effectively used as a new solar activity index both on long time scales (years) and on times scales of the order of a month or even days. The suggested new activity index is shown to have advantages over the universally accepted indices. The depth of the He I 10830 Å line in quiet regions was found to increase from the phase of minimum solar activity to the phase of maximum by a factor of about 2. In active regions, this increase is less than 30%. The differences between the cyclic variations of the chromospheric He I 10830 Å line radiation in active and quiet structures on the solar disk are indicative of the probable differences in the nature of cyclicity and its manifestations in magnetic fields of different spatial scales. The background magnetic fields appear to vary during the solar cycle more strongly than do the local fields associated with sunspots, faculae, and activity complexes. We suggest using regular observations in the He I 10830 Å line to predict solar activity.     

New Evidence for Long-Term Persistence in the Sun's Activity

Possible persistence of sunspot activity was studied using rescaled range and detrended fluctuation analyses. In addition to actual Wolf numbers (1700–2000 A.D.), two solar proxies were used in this research, viz., an annual sunspot proxy obtained for 1090–1700 A.D. and sunspot numbers reconstructed from the decadal radiocarbon series (8005 B.C. – 1895 A.D). The reconstruction was made using a five-box carbon exchange model. Analyses showed that in all cases the scaling exponent is significantly higher than 0.5 in the range of scales from 25 yr up to 3000 yr. This indicates the existence of a long-term memory in solar activity, in agreement with results obtained for other solar indices.     

The time-delayed solar cycle luminosity modulation by sub-surface magnetic flux tubes

In order to explore the mechanism of the solar cycle luminosity change observed by the Active Cavity Radiometer Irradiance Monitor (ACRIM) I experiment on board of the spacecraft Solar Maximum Mission, we examined running mean time profiles of the daily ACRIM data from the declining phase of solar cycle 21 to the rising phase of solar cycle 22. By comparing them with those of the daily sunspot number, integrated surface magnetic field flux, integrated He I 10830 Å line equivalent width data, and two kinds of data sets of the daily integrated Ca II K line index as indices of the surface magnetic activities, we found (i) that the running mean time profiles of the six independent data sets have several peaks and valleys in common in one solar cycle with time intervals on the order of a few hundreds of days, and (ii) that the peaks and valleys of the ACRIM data profiles followed the peaks and valleys of all the other five indices of the surface activities by 40 to 60 days. This time delay phenomenon suggests (i) that the luminosity modulation was not directly caused by dark and bright features of the surface magnetic activities that the other five indices represent, and (ii) that the missing sunspot radiative flux which was blocked by sub-surface magnetic flux tubes of sunspots and sunspot groups should be re-radiated 40 to 60 days after the surface emergence of the magnetic flux tubes. The concept of the time delay resolves the enigma of the missing sunspot radiative flux and the enigma of the ACRIM experiment that the luminosity dropped when a sunspot or a sunspot group appeared on the surface while the yearly mean of the luminosity decreased and increased along with the decrease and increase of the yearly sunspot number of the 11-year solar cycle. A model of the mechanism to understand these phenomena is presented and its application to other stars is suggested.     

Applying an Automatic Image-Processing Method to Synoptic Observations

We used an automatic image-processing method to detect solar-activity features observed in white light at the Kislovodsk Solar Station. This technique was applied to automatically or semi-automatically detect sunspots and active regions. The results of this automated recognition were verified with statistical data available from other observatories and revealed a high detection accuracy. We also provide parameters of sunspot areas, of the umbra, and of faculae as observed in Solar Cycle 23 as well as the magnetic flux of these active elements, calculated at the Kislovodsk Solar Station, together with white-light images and magnetograms from the Michaelson Doppler Imager onboard the Solar and Heliospheric Observatory (SOHO/MDI). The ratio of umbral and total sunspot areas during Solar Cycle 23 is ≈?0.19. The area of sunspots of the leading polarity was approximately 2.5 times the area of sunspots of the trailing polarity.     

太阳黑子面积数与黑子数指标在大气预报中的应用对比分析

分别应用太阳黑子视面积数和太阳黑子相对数代表太阳活动水平与天津夏季降水总量进行相关分析,结果表明黑子面积指标明显优于黑子数.     

Solar cycles during the Maunder minimum

We have obtained new consistent versions of the 400-yr time series of the Wolf sunspot number W, the sunspot group number G, and the total sunspot area S (or the total sunspot magnetic flux Φ). We show that the 11-yr cycle did not cease during the Maunder minimum of solar activity. The characteristics of the extrema of individual 11-yr cycles in 1600–2005 have been determined in terms of the total sunspot area index. We provide arguments for using alternating (“magnetic”) time series of indices in investigating the solar cyclicity.