Abstract
Kodaikanal Solar Observatory (KoSO) possesses one of world’s longest and homogeneous records of sunspot observations that span more than a century (1904–2017). Interestingly, these observations (originally recorded in photographic plates/films) were taken with the same setup over this entire time period which makes this data unique and best suitable for long-term solar variability studies. A large part of this data, between 1921–2011, were digitized earlier and a catalog containing the detected sunspot parameters (e.g., area and location) was published in Mandal et al. (2017). In this article, we extend the earlier catalog by including new sets of data between 1904–1921 and 2011–2017. To this end, we digitize and calibrate these new datasets which include resolving the issue of random image orientation. We fix this by comparing the KoSO images with co-temporal data from Royal Greenwich Observatory. Following that, a semi-automated sunspot detection and automated umbra detection algorithm are implemented onto these calibrated images to detect sunspots and umbra. Additionally, during this catalog updation, we also filled data gaps in the existing KoSO sunspot catalog (1921–2011) by virtue of re-calibrating the “rouge” plates. This updated sunspot area series covering nearly 115 years (1904–2017) are being made available to the community and will be a unique source to study the long term variability of the Sun.
1 Introduction
Sunspots have always been a central part of our understanding of the Sun and its long-term variability (). The systematic and methodical observation of these spots has revealed that their appearance is periodic, with a periodicity of around 11 years, known as the solar cycle or solar activity cycle (; ). Furthermore, the magnetic nature of these spots (; ) makes them an ideal proxy for understanding solar magnetism and its complex variability (; ). Today, it is well established that the solar activity cycle is governed by the solar dynamo process operating in the convection zone of the Sun (; ; ). In the dynamo process, the solar activity cycle is the manifestation of the periodic nature of the large-scale solar magnetic field (poloidal ⇌ toroidal). As a consequence, the number of sunspots and the corresponding area covered by these spots on the solar surface is dictated by the strength of the toroidal field generated in the dynamo process. Hence, the historical observation of these spots carries a vital information about the nature of toroidal field in the past and it will be crucial for the reconstruction of historical global solar magnetic field (, ). Apart from that, sunspots and solar activity are also intrinsically linked to solar transient events such as solar flares and coronal mass ejections (CMEs) and their frequency of occurrence. Since these transients are the primary drivers for the space weather condition, the historical observation will have crucial role for the understanding of space weather condition in past.
Over the last 150 years, many observatories around the world have begun regular observations of the Sun. The Royal Observatory of Greenwich (RGO; ) had been the leader in such a campaign having the record of white-light observations from 1874 to 1976, which was later continued using the Solar Optical Observing Network (SOON) by the US Air Force (USAF). After a few decades, in 1904, Kodaikanal Solar Observatory (KoSO; ) also joined this campaign and started regular observation of the Sun on photographic plates/films in multi-wavelengths (white-light since 1904, Ca-K since 1904 and H-α since 1912), independently in India. In particular, the white-light observation taken at KoSO provides one of the most homogeneous data series for over 100 years. Since these observations are taken from the same location and using the same telescope (since 1918) for such an extended period, KoSO provides a unique data series ideal for the long-term study of the Sun. Although the white-light data had been digitized () since 1904; owing to calibration issues, it has not been utilized for the period 1904–1920. Therefore, the area series (), the study of penumbra to umbra area ratio (1923–2011; ; ) and solar differential rotation (1923–2011; ) were limited to the period 1921–2011.
In-spite of the fact that RGO is the only observatory in the world with white-light data for 1904–1920 (Cycle-14 and Cycle-15), these images are available either in low resolution or in the form of drawings1. In this regard, the high-resolution white-light data available from KoSO serves as the only data set for this time frame. It is almost impossible for the ground based observatories to have a uniform and homogeneous data because of the varying atmospheric conditions and bad weather. Hence, in the last 30 years, there has been considerable effort to make the homogeneous sunspot area series by cross-correlating the sunspot area from various observatories, e.g., , , , and recently . In all these studies, RGO data has been used primarily as a reference to cross-correlate other data because of its outstanding data coverage. Now, the availability of KoSO data will complement the existing data series and will be helpful in cross-correlating the RGO and KoSO data in these initial overlapping periods. In this article, we present the extension of the KoSO sunspot area series reported in for 1904–2017 along with the umbral area after resolving the calibration issues for the initial 17 years of data. We have also included data for the period of 2012–2017, which was previously not reported due to ongoing digitization. In this article we present the updated data statistics, the issues with the calibration and the resolution of these issues in Section 2. In Section 3, we will discuss the updated area series and its comparison with the earlier ones, and finally, in Section 4 we will summarize our findings.
2 Data
The white-light observation at Kodaikanal started in 1904 using a 10 cm objective lens telescope which was later replaced in 1912 by a better quality lens while maintaining same size. A few years later, a 15 cm achromatic lens was installed, and the same setup has been used to take observations since 13 June 1918 (), and is still in use. Since 2017, the unavailability of photographic films has interrupted regular observation, but it is still taken whenever the films are available. These observations at KoSO, initially taken on the photographic plates/films, have been digitized using 4 k×4 k CCD at Kodaikanal and made available for the community (; ). Here, we use the white-light digitized data for the period of 1904–1920 and 2012–2018, and will extend the sunspot-area series reported earlier (1921–2011; ) for the period of 1904–2017 (114 years, covering ≈11 solar cycles). An example of a digitized white-light image from the very initial period, is shown in Figure 1A. Furthermore, we discovered that there are a few observations during 1921–2011, that were missed in earlier area series (), so we included them in the updated and revised series.
FIGURE 1
2.1 Updated data statistics
In Figure 2A, we present the number of observations per year from 1904 to 2017, with grey histograms representing the observations used in area-series as reported by and blue ones representing the additional observations included in this new series. Since, Figure 2A includes the multiple observation from the same day hence, to get an idea about the data coverage, we counted the number of observing days in each year and plotted it against the years in Figure 2B. From Figure 2B, we infer that the KoSO provides continuous white-light observations in the last 114 years with coverage of ≈ 73%. However, Kodaikanal experiences a rainy season each year during July to November and this has an impact on the data coverage (see ).
FIGURE 2
2.1.1 Where does the Kodaikanal Solar Observatory stand?
In Figure 2C, we compare the extent of data from the various observatories, showing that the KoSO has one of the most extended white-light data series (114 years) in the world. After KoSO, RGO provides the longest white-light sunspot data series (103 years). Although KoSO offers the most extended series, if we look at the percentage data coverage, RGO stands far ahead of the KoSO with data coverage of compared to 73.5% of KoSO. Despite the fact that RGO has better coverage of the data but it has compiled the data from various observatories around the globe including KoSO (
2.2 Data calibration
We follow the same calibration steps, which include flat fielding, disk detection and disk centring–to bring the disk centre to the image centre, as described in
Now, the next step is to correctly orient these images to get the solar North at the top of the image. In the earlier work by
2.2.1 Orientation correction
So far, we have discussed the issue of image orientation, and now we will look at how we can get the correct orientation of the images for the aforementioned period. Here, we cannot use EW or NS reference line for the first 5 years of data (1904–1908) since we do not know which of them represents EW or NS. In principle, only the EW line in the rest of the data can be used as a reference to get proper orientation. However, these images were flipped in the EW or NS direction during the digitization process, so we have to look for an alternate method. For this we used the already available sunspot location information from the RGO digitized full-disk images. In this method, firstly we chose the closest observation from the RGO data series, noted the time difference between observations is less than 12 h. Then using the available sunspot heliographic location information, we create a dummy mask with the same image size as the KoSO sunspot detected binary mask. After that, the sunspot binary mask obtained from the KoSO data is overlapped with the dummy mask created from RGO data for five different possibilities 1) no change, i.e., correct orientation, 2) North-South flip, 3) East-West flip, 4) 90°clockwise and, 5) 90°anti-clockwise. We also looked for the overlap and mark the observations with the appropriate flag. There are few cases where no overlap is seen in any of the five mentioned cases; hence we flag these observations as “others.” An example of the steps mentioned above is represented in Figures 3A–H.
FIGURE 3

A representative example shows the steps we use to identify the correct orientation. (A) Shows the image we have taken from KoSO; here, red contours denote the identified sunspot regions. (B) Shows the dummy mask we create from the near-simultaneous RGO sunspot data. Here, small dark regions represent the location of sunspots as per RGO observation. (C–G) represents the binary mask of KoSO overlapped with RGO (dummy) mask for all five possible orientations. (H) Shows the corrected digitized image; in this case, it is N-S flipped.
Now, following the aforementioned steps, we go through 3,565 observations during 1904–1920 and flag them accordingly, as represented in Table 1. We see in Table 1 that the majority (, predominantly in 1904–1912) of the observation are NS flipped (for yearly distribution of different orientations, see Supplementary Figure S3). In Figures 4A,B, we show two such cases where the images are EW and NS flipped, respectively. Apart from that, we also notice that there are 73 (2%) observations for which we cannot get the correct orientation. When we carefully looked at them, we found that the most probable reason is the incorrect time of observation of these images, which leads to fallacious overlapping pair. The unavailability of the correct time of observation of these images makes it difficult to get an accurate orientation, and therefore, we do not include these observations in our analysis.
TABLE 1
| Orientation | Number of observations | Percentage |
|---|---|---|
| Total Number of Observations | 3,565 | |
| Correct | 1,576 | 44.2% |
| North-South Flip | 1892 | 53.1% |
| East-West Flip | 19 | 0.53% |
| 90°(Clock-wise) | 3 | % |
| −90°(Anti-clockwise) | 2 | % |
| Others | 73 | 2.0% |
Number of observations for different flags representing the orientation of the images during 1904–1920.
FIGURE 4

Two representative examples from 1904 to 1909 show the orientation correction (A) for East-West flip and (B) North-South flip images.
3 Results
To detect the sunspots, we use a modified version of the sunspot Tracking And Recognition Algorithm (STARA,
First and foremost, we compare the sunspot area obtained from KoSO digitized white-light data with the existing sunspot area series. Then we discuss the sunspot umbra area extracted from them.
3.1 Sunspot area series
3.1.1 Comparison with existing series
In Figures 5A,B we compare the daily sunspot area obtained from KoSO with the recently cross calibrated composite sunspot area series from
FIGURE 5

The scatter plot showing the comparison of the daily sunspot area, panels (A–C); and monthly averaged sunspot area, panels (D–F) with M2020 for 1904–1920, 2012–2017 and for the extended period of 1904–2017. In each panel, the line with unit slop and the linear fit with zero intercepts are represented by dashed red and solid blue lines, respectively.
3.1.2 Extended sunspot area series
Now, in Figure 6A, we show the variation of yearly averaged sunspot area from KoSO for the extended period of 1904–2017 (covering 11 solar cycles, Cycle 14 to Cycle 24) along with M2020. It is evident from Figure 6A that the new addition in the series shows a good match with M2020 during 1904–1920 and 2012–2017 whereas in 1921–2011, the yearly averaged value is lower than M2020, and the possible reasons have been already been reported and discussed in
FIGURE 6

In panel (A) yearly averaged sunspot area obtained from KoSO white-light digitized data (solid blue) is plotted as a function of time along with the yearly averaged M2020 (dashed red). Panel (B) show the latitude time plot for the extended period obtained from KoSO sunspot data.
3.2 Umbra area series
Following the automatic umbra detection method as explained in
FIGURE 7

In panels (A) and (B), we compare the KoSO monthly averaged umbra area with RGO and DPD umbra, respectively, in their corresponding periods. Panel (C) shows the solar cycle variation of yearly averaged umbra area compared to RGO and DPD umbra area. Panel (D) and (E) show the variation of q for smaller and larger sunspots for an extended period, respectively. The shaded regions represent the data used in
3.2.1 Penumbra to umbra area ratio
In this subsection, we extend the work presented in
In Figure 7D, we do not notice any variation in q for smaller spots even in the extended period, which is consistent with the findings of
4 Summary and conclusion
In this article, we present the sunspot area series from KoSO, for the period of 1904–1920, and 2012–2017, which were not included due to calibration issues and incomplete digitization in the earlier published series (
We have also compared the yearly averaged umbra area series with RGO (1904–1976) and DPD (1977–2017) in their corresponding periods, which show a good agreement between the umbral areas (cc) except in three cycles (Cycle 15, Cycle 16, and Cycle 17) where we notice KoSO umbral area significantly lower than RGO umbral area. We also calculated the penumbra to umbra area ratio (q) and compared it with RGO data. We noted that KoSO data do not show any long-term trend in the ratio for smaller (Area hem) as well as for larger (Area hem), which further supports the findings of
The availability of high-resolution white-light digitized data for Cycle 14 and Cycle 15 is a key asset for the long-term studies of the Sun, and it provides an excellent opportunity to cross-correlate the sunspot area data for observatories. In addition, the accessibility of a homogeneous and uniform white-light digitized data for such a long period observed from the same location and setup will benefit the community for the long-term studies of the Sun and its global magnetic field variability. Moreover, KoSO also provides Sun charts, which combine the multi-wavelength observation in a single drawing. These Sun charts are getting digitized and will help fill the gaps and make the series even more homogeneous. In future, we will be looking for machine learning (ML), artificial intelligence (AI) and deep learning based methods to further improve the sunspot detection method by making it completely automatic. Furthermore, ML and AI based method can be expended to the study of historical global solar magnetic field and space weather conditions.
The digitized white-light data and the area series for the extended period will be available at https://kso.iiap.res.in/new/data. The sunspot area series for the extended period presented here is also available at https://github.com/bibhuraushan/KoSoDigitalArchive.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://kso.iiap.res.in/new/data or https://github.com/bibhuraushan/KoSoDigitalArchive.
Author contributions
BJ has done the data calibration (orientation correction), umbra detection using an automatic algorithm, all the analysis presented and manuscript drafting. MH has done 80% of sunspot detection, and AP has done the rest 20% of sunspot detection and contributed to the calibration step of orientation correction. BR has performed the data’s basic calibration (flat fielding and disk detection) step. SM and DB have helped in the structuring of the draft and helped in the language correction.
Funding
This project is partially funded by the project grant DST/ICPS/CLUSTER/DataScience/2018/General/Sl. No. 18.
Acknowledgments
Kodaikanal Solar Observatory is a facility of the Indian Institute of Astrophysics, Bangalore, India. This data is now available for public use at http://kso.iiap.res.in through a service developed at IUCAA under the Data Driven Initiatives project funded by the National Knowledge Network. BJ would like to thank the entire digitization team for their tireless effort to digitize the historical data and make data publically available. A special mention to Satish, who has digitized all the plates/films after 2011.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2022.1019751/full#supplementary-material
References
1
BalmacedaL. A.SolankiS. K.KrivovaN. A.FosterS. (2009). A homogeneous database of sunspot areas covering more than 130 years. J. Geophys. Res.114, A07104. 10.1029/2009JA014299
2
BappuM. K. V. (1967). Solar Physics at kodaikanal. Sol. Phys.1, 151–156. 10.1007/BF00150312
3
BaranyiT.GyoriL.LudmányA.CoffeyH. E. (2001). Comparison of sunspot area data bases. Mon. Not. R. Astron. Soc.323, 223–230. 10.1046/j.1365-8711.2001.04195.x
4
BaranyiT.KirályS.CoffeyH. E. (2013). Indirect comparison of Debrecen and Greenwich daily sums of sunspot areas. Mon. Not. R. Astron. Soc.434, 1713–1720. 10.1093/mnras/stt1134
5
CarrascoV. M. S.VaqueroJ. M.TrigoR. M.GallegoM. C. (2018). A curious history of sunspot penumbrae: An update. Sol. Phys.293, 104. 10.1007/s11207-018-1328-z
6
CharbonneauP. (2010). Dynamo models of the solar cycle. Living Rev. Sol. Phys.7, 3. 10.12942/lrsp-2010-3
7
FliggeM.SolankiS. K. (1997). Inter-cycle variations of solar irradiance: Sunspot areas as a pointer. Sol. Phys.173, 427–439. 10.1023/A:1004971807172
8
HaleG. E. (1908). On the probable existence of a magnetic field in sun-spots. Astrophys. J.28, 315. 10.1086/141602
9
HasanS. S.MallikD. C. V.BagareS. P.RajaguruS. P. (2010). “Solar Physics at the kodaikanal observatory: A historical perspective,” in Magnetic coupling between the interior and atmosphere of the sun. Vol. 19 of Astrophysics and space science proceedings (Springer), 12–36. 10.1007/978-3-642-02859-5_3
10
HathawayD. H. (2013). A curious history of sunspot penumbrae. Sol. Phys.286, 347–356. 10.1007/s11207-013-0291-y
11
HathawayD. H. (2015). The solar cycle. Living Rev. Sol. Phys.12, 4. 10.1007/lrsp-2015-4
12
HouJ.-W.ZengS.-G.ZhengS.LuoX.-Y.DengL.-H.LiY.-Y.et al (2022). Chinese sunspot drawings and their digitization-(VII) sunspot penumbra to Umbra area ratio using the hand-drawing records from yunnan observatories. Res. Astron. Astrophys.22, 095012. 10.1088/1674-4527/ac7f87
13
JhaB. K.MandalS.BanerjeeD. (2018). Long-term variation of sunspot penumbra to umbra area ratio: A study using kodaikanal white-light digitized data. Proc. Int. Astron. Union13, 185–186. 10.1017/S1743921318001989
14
JhaB. K.MandalS.BanerjeeD. (2019). Study of sunspot penumbra to Umbra area ratio using kodaikanal white-light digitised data. Sol. Phys.294, 72. 10.1007/s11207-019-1462-2
15
JhaB. K.PriyadarshiA.MandalS.ChatterjeeS.BanerjeeD. (2021). Measurements of solar differential rotation using the century long kodaikanal sunspot data. Sol. Phys.296, 25. 10.1007/s11207-021-01767-8
16
JiangJ.CameronR. H.SchmittD.SchüsslerM. (2011). The solar magnetic field since 1700. II. Physical reconstruction of total, polar and open flux. Astron. Astrophys.528, A83. 10.1051/0004-6361/201016168
17
JiangJ.HathawayD. H.CameronR. H.SolankiS. K.GizonL.UptonL. (2014). Magnetic flux transport at the solar surface. Space Sci. Rev.186, 491–523. 10.1007/s11214-014-0083-1
18
MandalS.HegdeM.SamantaT.HazraG.BanerjeeD.RavindraB. (2017). Kodaikanal digitized white-light data archive (1921-2011): Analysis of various solar cycle features. Astron. Astrophys.601, A106. 10.1051/0004-6361/201628651
19
MandalS.KrivovaN. A.SolankiS. K.SinhaN.BanerjeeD. (2020). Sunspot area catalog revisited: Daily cross-calibrated areas since 1874. Astron. Astrophys.640, A78. 10.1051/0004-6361/202037547
20
NagovitsynY. A.PevtsovA. A.OsipovaA. A. (2017). Long-term variations in sunspot magnetic field-area relation. Astron. Nachr.338, 26–34. 10.1002/asna.201613035
21
ParkerE. N. (1955a). Hydromagnetic dynamo models. Astrophys. J.122, 293. 10.1086/146087
22
ParkerE. N. (1955b). The formation of sunspots from the solar toroidal field. Astrophys. J.121, 491. 10.1086/146010
23
ParkerE. N. (1975). The generation of magnetic fields in astrophysical bodies. X - magnetic buoyancy and the solar dynamo. Astrophys. J.198, 205–209. 10.1086/153593
24
RavindraB.PriyaT. G.AmareswariK.PriyalM.NaziaA. A.BanerjeeD. (2013). Digitized archive of the Kodaikanal images: Representative results of solar cycle variation from sunspot area determination. Astron. Astrophys.550, A19. 10.1051/0004-6361/201220416
25
SchwabeH. (1844). Sonnenbeobachtungen im Jahre 1843. Von Herrn Hofrath Schwabe in Dessau. Astr. Nachr.21, 234–235. 10.1002/asna.18440211505
26
SivaramanK. R.GuptaS. S.HowardR. F. (1993). Measurement of Kodaikanal white-light images. Sol. Phys.146, 27–47. 10.1007/BF00662168
27
SolankiS. K. (2003). Sunspots: An overview. Astronomy Astrophysics Rev.11, 153–286. 10.1007/s00159-003-0018-4
28
TlatovA. G.PevtsovA. A. (2014). Bimodal distribution of magnetic fields and areas of sunspots. Sol. Phys.289, 1143–1152. 10.1007/s11207-013-0382-9
29
WatsonF. T.FletcherL.MarshallS. (2011). Evolution of sunspot properties during solar cycle 23. Astron. Astrophys.533, A14. 10.1051/0004-6361/201116655
30
WillisD. M.CoffeyH. E.HenwoodR.ErwinE. H.HoytD. V.WildM. N.et al (2013). The Greenwich photo-heliographic results (1874 - 1976): Summary of the observations, applications, datasets, definitions and errors. Sol. Phys.288, 117–139. 10.1007/s11207-013-0311-y
Summary
Keywords
Sun, sunspots, solar cycle, Kodaikanal Solar Observatory, white-light, sunspot area, umbra area
Citation
Jha BK, Hegde M, Priyadarshi A, Mandal S, Ravindra B and Banerjee D (2022) Extending the sunspot area series from Kodaikanal Solar Observatory. Front. Astron. Space Sci. 9:1019751. doi: 10.3389/fspas.2022.1019751
Received
15 August 2022
Accepted
13 October 2022
Published
25 October 2022
Volume
9 - 2022
Edited by
Maria Ann Weber, Delta State University, United States
Reviewed by
Laure Lefevre, Royal Observatory of Belgium, Belgium
Sergei Zharkov, University of Hull, United Kingdom
Updates

Check for updates
Copyright
© 2022 Jha, Hegde, Priyadarshi, Mandal, Ravindra and Banerjee.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bibhuti Kumar Jha, bibhuti@aries.res.in; Dipankar Banerjee, dipu@aries.res.in
This article was submitted to Stellar and Solar Physics, a section of the journal Frontiers in Astronomy and Space Sciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.