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.     

大黑子群对太阳总辐照度的影响

本文用云南天文台在第２２周太阳活动峰年期间拍摄到的大太阳黑子群照相资料，太阳黑子目视描述资料，以及Ｎｉｍｂｕｓ—７卫星上辐射计测量的太阳总辐照度，分别计算了太阳总辐射照度与大黑子群的本影视面积，大黑子群全群视面积和日面上全部黑子的总视面积的相关系数。结果表明，太阳总辐射照度与这三种视面积均存在强的负相关。其中与大黑子群本影视面积的相关最强，其次是与全群视面积的相关，最后是与日面上全部黑子的总视面积的相关。并对以上结果和其它有关结果进行了分析和讨论。     

Changes in measured vector magnetic fields when transformed into heliographic coordinates

In this paper we investigate the changes that occur in measured magnetic fields when they are transformed into a heliographic coordinate system. To carry out this investigation we took measurements of the vector magnetic field of an active region that was observed at 1/3 the solar radius from disk center and transformed the observed field into heliographic coordinates. We also examined differences in the calculated potential field that occur when the heliographic normal component of the field is used as the boundary condition rather than the observed line-of-sight component. The results of this analysis show (1) that the observed fields of sunspots more closely resemble the generally accepted picture of the distribution of umbral fields if they are displayed in heliographic coordinates, (2) that the differences in the potential calculations are less than 200 G in field strength and 20° in field azimuth outside sunspots, and (3) that differences in the two potential calculations in the sunspot areas are no more than 400 G in field strength but range from 60 to 80° in field azimuth in localized umbral areas.     

The Greenwich Photo-heliographic Results (1874 – 1976): Initial Corrections to the Printed Publications

A new sunspot and faculae digital dataset for the interval 1874?–?1955 has been prepared under the auspices of the NOAA National Geophysical Data Center (NGDC). This digital dataset contains measurements of the positions and areas of both sunspots and faculae published initially by the Royal Observatory, Greenwich, and subsequently by the Royal Greenwich Observatory (RGO), under the title Greenwich Photo-heliographic Results (GPR), 1874?–?1976. Quality control (QC) procedures based on logical consistency have been used to identify the more obvious errors in the RGO publications. Typical examples of identifiable errors are North versus South errors in specifying heliographic latitude, errors in specifying heliographic (Carrington) longitude, errors in the dates and times, errors in sunspot group numbers, arithmetic errors in the summation process, and the occasional omission of solar ephemerides. Although the number of errors in the RGO publications is remarkably small, an initial table of necessary corrections is provided for the interval 1874?–?1917. Moreover, as noted in the preceding companion papers, the existence of two independently prepared digital datasets, which both contain information on sunspot positions and areas, makes it possible to outline a preliminary strategy for the development of an even more accurate digital dataset. Further work is in progress to generate an extremely reliable sunspot digital dataset, based on the long programme of solar observations supported first by the Royal Observatory, Greenwich, and then by the Royal Greenwich Observatory.     

Total Solar Irradiance Variation During Rapid Sunspot Growth

Large sunspot areas correspond to dips in the total solar irradiance (TSI), a phenomenon associated with the local suppression of convective energy transport in the spot region. This results in a strong correlation between sunspot area and TSI. During the growth phase of a sunspot other physics may affect this correlation; if the physical growth of the sunspot resulted in surface flows affecting the temperature, for example, we might expect to see an anomalous variation in TSI. In this paper we study NOAA active region 8179, in which large sunspots suddenly appeared near disk center, at a time (March 1998) when few competing sunspots or plage regions were present on the visible hemisphere. We find that the area/TSI correlation does not significantly differ from the expected pattern of correlation, a result consistent with a large thermal conductivity in solar convection zone. In addition we have searched for a smaller-scale effect by analyzing white-light images from MDI (the Michelson Doppler Imager) on SOHO. A representative upper-limit energy consistent with the images is on the order of 3×10³¹ ergs, assuming the time scale of the actual spot area growth. This is of the same order of magnitude as the buoyant energy of the spot emergence even if it is shallow. We suggest that detailed image analyses of sunspot growth may therefore show `transient bright rings' at a detectable level.     

Subsurface Flows in and Around Active Regions with Rotating and Non-rotating Sunspots

The temporal variation of the horizontal velocity in sub-surface layers beneath three different types of active region is studied using the technique of ring diagrams. In this study, we select active regions (ARs) 10923, 10930, 10935 from three consecutive Carrington rotations: AR 10930 contains a fast-rotating sunspot in a strong emerging active region while other two have non-rotating sunspots with emerging flux in AR 10923 and decaying flux in AR 10935. The depth range covered is from the surface to about 12 Mm. In order to minimize the influence of systematic effects, the selection of active and quiet regions is made so that these were observed at the same heliographic locations on the solar disk. We find a significant variation in both components of the horizontal velocity in active regions as compared to quiet regions. The magnitude is higher in emerging-flux regions than in the decaying-flux region, in agreement with earlier findings. Further, we clearly see a significant temporal variation in depth profiles of both zonal and meridional flow components in AR 10930, with the variation in the zonal component being more pronounced. We also notice a significant influence of the plasma motion in areas closest to the rotating sunspot in AR 10930, while areas surrounding the non-rotating sunspots in all three cases are least affected by the presence of the active region in their neighborhood.     

Identification of Sunspots on Full-Disk Solar Images Using Wavelet Analysis

A new method to detect sunspots on full-disk solar images, recorded in an optical wide spectral band and in the line of Ca ii K₁, is presented. It is based on the compact wavelet transform and a process that automatically identifies the noise in the images of the Sun. This method estimates also its center and the limb positions with a precision of a few milliarcseconds. Precisions of 1 and 8 milliarcseconds are, respectively, found for the disk center position and the radius of simulated solar images. The results obtained on sunspot identification on solar images recorded during April 2002 at Meudon Observatory (France) reveal that we are able to detect 3.35 % additional sunspots compared to a manual method. Our method preserves the sunspot areas better than 95 %.     

Digitization of sunspot drawings by Spörer made in 1861–1894

Most of our knowledge about the Sun's activity cycle arises from sunspot observations over the last centuries since telescopes have been used for astronomy. The German astronomer Gustav Spörer observed almost daily the Sun from 1861 until the beginning of 1894 and assembled a 33‐year collection of sunspot data covering a total of 445 solar rotation periods. These sunspot drawings were carefully placed on an equidistant grid of heliographic longitude and latitude for each rotation period, which were then copied to copper plates for a lithographic reproduction of the drawings in astronomical journals. In this article, we describe in detail the process of capturing these data as digital images, correcting for various effects of the aging print materials, and preparing the data for contemporary scientific analysis based on advanced image processing techniques. With the processed data we create a butterfly diagram aggregating sunspot areas, and we present methods to measure the size of sunspots (umbra and penumbra) and to determine tilt angles of active regions. A probability density function of the sunspot area is computed, which conforms to contemporary data after rescaling. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

A Curious History of Sunspot Penumbrae

Daily records of sunspot group areas compiled by the Royal Observatory, Greenwich, from May of 1874 through 1976 indicate a curious history for the penumbral areas of the smaller sunspot groups. On average, the ratio of penumbral area to umbral area in a sunspot group increases from 5 to 6 as the total sunspot group area increases from 100 to 2000 μHem (a μHem is 10^?6 the area of a solar hemisphere). This relationship does not vary substantially with sunspot group latitude or with the phase of the sunspot cycle. However, for the sunspot groups with total areas <?100 μHem, this ratio changes dramatically and systematically through this historical record. The ratio for these smallest sunspots is near 5.5 from 1874 to 1900. After a rapid rise to more than 7 in 1905, it drops smoothly to less than 3 by 1930 and then rises smoothly back to more than 7 in 1961. It then returns to near 5.5 from 1965 to 1976. The smooth variation from 1905 to 1961 shows no indication of any step-like changes that might be attributed to changes in equipment or personnel. The overall level of solar activity was increasing monotonically during this time period when the penumbra-to-umbra area ratio dropped to less than half its peak value and then returned. If this history can be confirmed by other observations (e.g. Mt. Wilson or Kodaikanal), it may impact our understanding of penumbra formation, our dynamo models, and our estimates of historical changes in the solar irradiance.     

Sunspot Positions and Areas from Observations by Galileo Galilei

Sunspot records in the seventeenth century provide important information on the solar activity before the Maunder minimum, yielding reliable sunspot indices and the solar butterfly diagram. Galilei’s letters to Cardinal Francesco Barberini and Marcus Welser contain daily solar observations on 3?–?11 May, 2 June?–?8 July, and 19?–?21 August 1612. These historical archives do not provide the time of observation, which results in uncertainty in the sunspot coordinates. To obtain them, we present a method that minimizes the discrepancy between the sunspot latitudes. We provide areas and heliographic coordinates of 82 sunspot groups. In contrast to Sheiner’s butterfly diagram, we found only one sunspot group near the Equator. This provides a higher reliability of Galilei’s drawings. Large sunspot groups are found to emerge at the same longitude in the northern hemisphere from 3 May to 21 August, which indicates an active longitude.     

Automatic Sunspots Detection on Full-Disk Solar Images using Mathematical Morphology

Sunspots are solar features located in active regions of the Sun, whose number is an indicator of the Sun’s magnetic activity. Therefore accurate detection and classification of sunspots are fundamental for the elaboration of solar activity indices such as the Wolf number. However, irregularities in the shape of the sunspots and their variable intensity and contrast with the surroundings, make their automated detection from digital images difficult. Here, we present a morphological tool that has allowed us to construct a simple and automatic procedure to treat digital photographs obtained from a solar telescope, and to extract the main features of sunspots. Comparing the solar indices computed with our algorithm against those obtained with the previous method exhibit an obvious improvement. A favorable comparison of the Wolf sunspot number time series obtained with our methodology and from other reference observatories is also presented. Finally, we compare our sunspot and group detection to that of other observatories.     

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)     

On the Reliability of the de la rue Sunspot Area Measurements

In the 19th century De la Rue, Stewart, and Loewy carried out a compilation of drawings and photographs of the solar sunspots in the period 1832–1868. From these drawings and photographs, they determined fortnightly values of the sunspot areas. In this work, monthly values of the sunspot areas for the period 1832–1868 are calculated and the reliability of these data in terms of the solar activity indices is discussed.     

Sunspot Group Decay

We examine daily records of sunspot group areas (measured in millionths of a solar hemisphere or μHem) for the last 130 years to determine the rate of decay of sunspot group areas. We exclude observations of groups when they are more than 60° in longitude from the central meridian and only include data when at least three days of observations are available following the date of maximum area for a group’s disk passage. This leaves data for over 18 000 measurements of sunspot group decay. We find that the decay rate increases linearly from 28 μHem day⁻¹ to about 140 μHem day⁻¹ for groups with areas increasing from 35 μHem to 1000 μHem. The decay rate tends to level off for groups with areas larger than 1000 μHem. This behavior is very similar to the increase in the number of sunspots per group as the area of the group increases. Calculating the decay rate per individual sunspot gives a decay rate of about 3.65 μHem day⁻¹ with little dependence upon the area of the group. This suggests that sunspots decay by a Fickian diffusion process with a diffusion coefficient of about 10 km² s⁻¹. Although the 18 000 decay rate measurements are lognormally distributed, this can be attributed to the lognormal distribution of sunspot group areas and the linear relationship between area and decay rate for the vast majority of groups. We find weak evidence for variations in decay rates from one solar cycle to another and for different phases of each sunspot cycle. However, the strongest evidence for variations is with latitude and the variations with cycle and phase of each cycle can be attributed to this variation. High latitude spots tend to decay faster than low latitude spots.     

On umbral flashes in different sunspot groups

The results of a statistical investigation of the occurrence of umbral flashes for 40 sunspot groups are reported for the period 1966–1983. The following characteristics were chosen for the analysis: (a) position on the solar disk; (b) group area; (c) sunspot area; (d) maximum magnetic field strength of a sunspot; (e) modified Zürich class; (f) sunspot age; (g) magnetic structure; and (h) flare activity of a group. The dependence of umbral flashes on magnetic structure of a sunspot is the most essential feature. The absence of umbral flashes in the umbrae of main sunspots perhaps may be used as one of the predictors of flare activity.     

The effects of sunspots on solar irradiance

Sunspots have an obvious direct effect upon the visible radiant energy falling upon the Earth. We show how to estimate this effect and compare it quantitatively with recent observations of the solar total irradiance (Willson et al., 1981). The sunspots explain about half of the total observed variance of one-day averages. Since the sunspot effect on irradiance produces an asymmetry of the solar radiation, rather than (necessarily) a variation of the total luminosity, we have also estimated the sunspot population on the invisible hemisphere. This extrapolation allows us to estimate the true luminosity deficit produced by sunspots, in a manner that tends toward the correct long-term average value. We find no evidence for instantaneous global re-emission to compensate for the sunspot flux deficit.     

Statistics of the largest sunspot and facular areas per solar cycle

The statistics of extreme values is used to investigate the statistical properties of the largest areas of sunspots and photospheric faculae per solar cycle. The largest values of the synodic-solar-rotation mean areas of umbrae, whole spots and faculae, which have been recorded for nine solar cycles, are each shown to comply with the general form of the extreme value probability function. Empirical expressions are derived for the three extreme value populations from which the characteristic statistical parameters, namely the mode, median, mean and standard deviation, can be calculated for each population. These three extreme value populations are also used to find the expected ranges of the extreme areas in a group of solar cycles as a function of the number of cycles in the group. The extreme areas of umbrae and whole spots have a dispersion comparable to that found by Siscoe for the extreme values of sunspot number, whereas the extreme areas of faculae have a smaller dispersion which is comparable to that found by Siscoe for the largest geomagnetic storm per solar cycle. The expected range of the largest sunspot area per solar cycle for a group of one hundred cycles appears to be inconsistent with the existence of the prolonged periods of sunspot minima that have been inferred from the historical information on solar variability. This inconsistency supports the contention that there are temporal changes of solar-cycle statistics during protracted periods of sunspot minima (or maxima). Indeed, without such temporal changes, photospheric faculae should have been continually observable throughout the lifetime of the Sun.     

10.7-cm Solar radio flux and the magnetic complexity of active regions

During sunspot cycles 20 and 21, the maximum in smoothed 10.7-cm solar radio flux occurred about 1.5 yr after the maximum smoothed sunspot number, whereas during cycles 18 and 19 no lag was observed. Thus, although 10.7-cm radio flux and Zürich suspot number are highly correlated, they are not interchangeable, especially near solar maximum. The 10.7-cm flux more closely follows the number of sunspots visible on the solar disk, while the Zürich sunspot number more closely follows the number of sunspot groups. The number of sunspots in an active region is one measure of the complexity of the magnetic structure of the region, and the coincidence in the maxima of radio flux and number of sunspots apparently reflects higher radio emission from active regions of greater magnetic complexity. The presence of a lag between sunspot-number maximum and radio-flux maximum in some cycles but not in others argues that some aspect of the average magnetic complexity near solar maximum must vary from cycle to cycle. A speculative possibility is that the radio-flux lag discriminates between long-period and short-period cycles, being another indicator that the solar cycle switches between long-period and short-period modes.Operated by the Association of Universities for Research in Astronomy, Inc. under contract with the National Science Foundation.     

Sunspot Numbers and Areas from the Madrid Astronomical Observatory (1876 – 1986)

The solar program of the Astronomical Observatory of Madrid started in 1876. Observations were made in this institution to determine sunspot numbers and areas for ten solar cycles. The program was completed in 1986 and the resulting data have been published in various Spanish scientific publications. Four periods of this program (with different observers and instruments) were identified with the aid of the interesting metadata that has been made available. In the present work, the published data were retrieved and digitized. Their subsequent analysis showed that most of these data could be considered reliable given their very high correlation with reference indices (international sunspot number, group sunspot number, and sunspot area). An abrupt change emerged in the sunspots/groups ratio in 1946, which lasted until 1972.