-
Differential Rotation of the Solar Chromosphere: A Century-long Perspective from Kodaikanal Solar Observatory Ca II K Data
Authors:
Dibya Kirti Mishra,
Srinjana Routh,
Bibhuti Kumar Jha,
Theodosios Chatzistergos,
Judhajeet Basu,
Subhamoy Chatterjee,
Dipankar Banerjee,
Ilaria Ermolli
Abstract:
Chromospheric differential rotation is a key component in comprehending the atmospheric coupling between the chromosphere and the photosphere at different phases of the solar cycle. In this study, we therefore utilize the newly calibrated multidecadal Ca II K spectroheliograms (1907-2007) from the Kodaikanal Solar Observatory (KoSO) to investigate the differential rotation of the solar chromospher…
▽ More
Chromospheric differential rotation is a key component in comprehending the atmospheric coupling between the chromosphere and the photosphere at different phases of the solar cycle. In this study, we therefore utilize the newly calibrated multidecadal Ca II K spectroheliograms (1907-2007) from the Kodaikanal Solar Observatory (KoSO) to investigate the differential rotation of the solar chromosphere using the technique of image cross-correlation. Our analysis yields the chromospheric differential rotation rate $Ω(θ) = (14.61\pm 0.04 - 2.18\pm 0.37\sin^2θ - 1.10 \pm 0.61\sin^4θ)^\circ{\rm /day}$. These results suggest the chromospheric plages exhibit an equatorial rotation rate 1.59% faster than the photosphere when compared with the differential rotation rate measured using sunspots and also a smaller latitudinal gradient compared to the same. To compare our results to those from other observatories, we have applied our method on a small sample of Ca II K data from Rome, Meudon, and Mt. Wilson observatories, which support our findings from KoSO data. Additionally, we have not found any significant north-south asymmetry or any systematic variation in chromospheric differential rotation over the last century.
△ Less
Submitted 30 November, 2023;
originally announced November 2023.
-
Analysis of full-disc H alpha observations: Carrington maps and filament properties in 1909-2022
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Dipankar Banerjee,
Teresa Barata,
Ioannis Chouinavas,
Mariachiara Falco,
Ricardo Gafeira,
Fabrizio Giorgi,
Yoichiro Hanaoka,
Natalie A. Krivova,
Viktor V. Korokhin,
Ana Lourenço,
Gennady P. Marchenko,
Jean-Marie Malherbe,
Nuno Peixinho,
Paolo Romano,
Takashi Sakurai
Abstract:
Full disc observations of the Sun in the H$α$ line provide information about the solar chromosphere and in particular about the filaments, which are dark and elongated features that lie along magnetic field polarity inversion lines. This makes them important for studies of solar magnetism. Since full disc H$α$ observations have been performed at various sites since 1800s, with regular photographic…
▽ More
Full disc observations of the Sun in the H$α$ line provide information about the solar chromosphere and in particular about the filaments, which are dark and elongated features that lie along magnetic field polarity inversion lines. This makes them important for studies of solar magnetism. Since full disc H$α$ observations have been performed at various sites since 1800s, with regular photographic data having started in the beginning of the 20th century, they are an invaluable source of information on past solar magnetism. In this work we aimed at deriving accurate information about filaments from historical and modern full disc H$α$ observations. We have consistently processed observations from 15 H$α$ archives spanning 1909-2022. Our data processing includes photometric calibration of the data stored on photographic plates. We have constructed also Carrington maps from the calibrated H$α$ images. We find that filament areas are affected by the bandwidth of the observation. Thus, cross-calibration of the filament areas derived from different archives is needed. We have produced a composite of filament areas from individual archives by scaling all of them to the Meudon series. Our composite butterfly diagram shows very distinctly the common features of filament evolution, that is the poleward migration as well as a decrease in the mean latitude of filaments as the cycle progresses. We also find that during activity maxima, filaments on average cover about 1% of the solar surface. We see only a weak change in the amplitude of cycles in filament areas, in contrast to sunspot and plage areas. Analysis of H$α$ data for archives with contemporaneous Ca II K observations allowed us to identify and verify archive inconsistencies, which will also have implications for reconstructions of past solar magnetism and irradiance from Ca II K data.
△ Less
Submitted 18 October, 2023; v1 submitted 18 September, 2023;
originally announced September 2023.
-
Investigating the effect of solar ambient and data characteristics on Ca II K observations and line profile measurements
Authors:
M. Murabito,
I. Ermolli,
T. Chatzistergos,
S. Jafarzadeh,
F. Giorgi,
L. Rouppe van der Voort
Abstract:
We analysed state-of-the-art observations of the solar atmosphere to investigate the dependence of the \ca brightness of several solar features on spectral bandwidth and spatial resolution of the data. Specifically, we study data obtained at the Swedish Solar Telescope with the CRiSP and CHROMIS instruments. The analyzed data, which are characterized by spectral bandwidth of 0.12 Å\ and spatial re…
▽ More
We analysed state-of-the-art observations of the solar atmosphere to investigate the dependence of the \ca brightness of several solar features on spectral bandwidth and spatial resolution of the data. Specifically, we study data obtained at the Swedish Solar Telescope with the CRiSP and CHROMIS instruments. The analyzed data, which are characterized by spectral bandwidth of 0.12 Å\ and spatial resolution of 0.078\arcsec, were acquired close to disc center by targeting a quiet Sun area and an active region. We convolved the original observations with Gaussian kernels to degrade their spectral bandwidth and spatial resolution to the instrumental characteristics of the most prominent series of \ca observations available to date. We then studied the effect of data degradation on the observed regions and on parameters derived from \ca line measurements that are largely employed as diagnostics of the solar and stellar chromospheres. We find that the effect of degrading the spectral resolution of \ca observations and line profiles depends on both the employed bandwidth and observed solar region. Besides, we found that the spatial degradation impacts the data characterized by a broad bandwidth to a larger extent compared to those acquired with a narrow band. However, the appearance of the observed solar regions is only slightly affected by the spatial resolution of data with bandwidths up to 1 Å\ and in the range [3,10] Å. Finally, we derived relationships that can be used to intercalibrate results from observations taken with different instruments in diverse regions of the solar atmosphere.
△ Less
Submitted 30 March, 2023;
originally announced March 2023.
-
Long-term changes in solar activity and irradiance
Authors:
Theodosios Chatzistergos,
Natalie A. Krivova,
Kok Leng Yeo
Abstract:
The Sun is the main energy source to Earth, and understanding its variability is of direct relevance to climate studies. Measurements of total solar irradiance exist since 1978, but this is too short compared to climate-relevant time scales. Coming from a number of different instruments, these measurements require a cross-calibration, which is not straightforward, and thus several composite record…
▽ More
The Sun is the main energy source to Earth, and understanding its variability is of direct relevance to climate studies. Measurements of total solar irradiance exist since 1978, but this is too short compared to climate-relevant time scales. Coming from a number of different instruments, these measurements require a cross-calibration, which is not straightforward, and thus several composite records have been created. All of them suggest a marginally decreasing trend since 1996. Most composites also feature a weak decrease over the entire period of observations, which is also seen in observations of the solar surface magnetic field and is further supported by Ca II K data. Some inconsistencies, however, remain and overall the magnitude and even the presence of the long-term trend remain uncertain. Different models have been developed, which are used to understand the irradiance variability over the satellite period and to extend the records of solar irradiance back in time. Differing in their methodologies, all models require proxies of solar magnetic activity as input. The most widely used proxies are sunspot records and cosmogenic isotope data on centennial and millennial time scale, respectively. None of this, however, offers a sufficiently good, independent description of the long-term evolution of faculae and network responsible for solar brightening. This leads to uncertainty in the amplitude of the long-term changes in solar irradiance. Here we review recent efforts to improve irradiance reconstructions on time scales longer than the solar cycle and to reduce the existing uncertainty in the magnitude of the long-term variability. In particular, we highlight the potential of using 3D magnetohydrodynamical simulations of the solar atmosphere as input to more physical irradiance models and of historical full-disc Ca II K observations encrypting direct facular information back to 1892.
△ Less
Submitted 16 October, 2023; v1 submitted 6 March, 2023;
originally announced March 2023.
-
Re-calibration of the Sunspot Number: Status Report
Authors:
F. Clette,
L. Lefèvre,
T. Chatzistergos,
H. Hayakawa,
V. M. Carrasco,
R. Arlt,
E. W. Cliver,
T. Dudok de Wit,
T. Friedli,
N. Karachik,
G. Kopp,
M. Lockwood,
S. Mathieu,
A. Muñoz-Jaramillo,
M. Owens,
D. Pesnell,
A. Pevtsov,
L. Svalgaard,
I. G. Usoskin,
L. van Driel-Gesztelyi,
J. M. Vaquero
Abstract:
We report progress on the ongoing recalibration of the Wolf sunspot number (SN) and Group sunspot number (GN) following the release of version 2.0 of SN in 2015. This report constitutes both an update of the efforts reported in the 2016 Topical Issue of Solar Physics and a summary of work by the International Space Science Institute (ISSI) International Team formed in 2017 to develop optimal SN an…
▽ More
We report progress on the ongoing recalibration of the Wolf sunspot number (SN) and Group sunspot number (GN) following the release of version 2.0 of SN in 2015. This report constitutes both an update of the efforts reported in the 2016 Topical Issue of Solar Physics and a summary of work by the International Space Science Institute (ISSI) International Team formed in 2017 to develop optimal SN and GN re-construction methods while continuing to expand the historical sunspot number database. Significant progress has been made on the database side while more work is needed to bring the various proposed SN and (primarily) GN reconstruction methods closer to maturity, after which the new reconstructions (or combinations thereof) can be compared with (a) ``benchmark'' expectations for any normalization scheme (e.g., a general increase in observer normalization factors going back in time), and (b) independent proxy data series such as F10.7 and the daily range of variations of Earth's undisturbed magnetic field. New versions of the underlying databases for SN and GN will shortly become available for years through 2022 and we anticipate the release of next versions of these two time series in 2024.
△ Less
Submitted 5 January, 2023;
originally announced January 2023.
-
Rome Precision Solar Photometric Telescope: precision solar full-disk photometry during solar cycles 23-25
Authors:
Ilaria Ermolli,
Fabrizio Giorgi,
Theodosios Chatzistergos
Abstract:
The Rome Precision Solar Photometric Telescope (Rome/PSPT) is a ground-based telescope engaged in precision solar photometry. It has a 27-year database of full-disk images of the photosphere and chromosphere beginning in 1996 and continuing to 2022. The solar images have been obtained daily, weather permitting, with approximately 2 arcsec/pixel scale in Ca II K line at 393.3 nm, G-band at 430.6 nm…
▽ More
The Rome Precision Solar Photometric Telescope (Rome/PSPT) is a ground-based telescope engaged in precision solar photometry. It has a 27-year database of full-disk images of the photosphere and chromosphere beginning in 1996 and continuing to 2022. The solar images have been obtained daily, weather permitting, with approximately 2 arcsec/pixel scale in Ca II K line at 393.3 nm, G-band at 430.6 nm, and continuum in the blue and red parts of the spectrum at 409.4 nm and 607.2 nm, respectively. Regular observations were also performed at the green continuum at 535.7 nm for a period of about 18 months. Since the first-light, Rome/PSPT operations have been directed at understanding the source of short-and long-term solar irradiance changes, spanning from one minute to several months, and from one year to a few solar cycles, respectively. However, Rome/PSPT data have also served to study a variety of other topics, including the photometric properties of solar disk features and of the supergranulation manifested by the chromospheric network. Moreover, they have been unique in allowing to connect series of historical and modern full-disk solar observations, especially the Ca II K line data. Here, we provide an overview of the Rome/PSPT telescope and of the solar monitoring carried out with it from its first light to the present, across solar cycles 23-25. We also briefly describe the main results achieved with Rome/PSPT data, and give an overview of new results being derived with the whole time series of observations covering the period 1996-2022.
△ Less
Submitted 18 November, 2022;
originally announced November 2022.
-
Full-disc Ca II K observations -- a window to past solar magnetism
Authors:
Theodosios Chatzistergos,
Natalie A. Krivova,
Ilaria Ermolli
Abstract:
First such observations were made in 1892 and since then various sites around the world have carried out regular observations, with Kodaikanal, Meudon, Mt Wilson, and Coimbra being some of the most prominent ones. By now, Ca II K observations from over 40 different sites allow an almost complete daily coverage of the last century. Ca II K images provide direct information on plage and network regi…
▽ More
First such observations were made in 1892 and since then various sites around the world have carried out regular observations, with Kodaikanal, Meudon, Mt Wilson, and Coimbra being some of the most prominent ones. By now, Ca II K observations from over 40 different sites allow an almost complete daily coverage of the last century. Ca II K images provide direct information on plage and network regions on the Sun and, through their connection to solar surface magnetic field, offer an excellent opportunity to study solar magnetism over more than a century. This makes them also extremely important, among others, for solar irradiance reconstructions and studies of the solar influence on Earth's climate. However, these data also suffer from numerous issues, which for a long time have hampered their analysis. Without properly addressing these issues, Ca II K data cannot be used to their full potential. Here, we first provide an overview of the currently known Ca II K data archives and sources of the inhomogeneities in the data, before discussing existing processing techniques, followed by a recap of the main results derived with such data so far.
△ Less
Submitted 24 October, 2022;
originally announced October 2022.
-
Scrutinizing the relationship between plage areas and sunspot areas and numbers
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Natalie A. Krivova,
Teresa Barata,
Sara Carvalho,
Jean-Marie Malherbe
Abstract:
Studies and reconstructions of past solar activity require data on sunspots as well as faculae/plage and network. Such data are also important for understanding the magnetic activity and variability of the Sun and Sun-like stars. The longest available direct faculae/plage datasets are white-light facular and Ca II K observations going back to 1874 and 1892, respectively. Prior to that time, the on…
▽ More
Studies and reconstructions of past solar activity require data on sunspots as well as faculae/plage and network. Such data are also important for understanding the magnetic activity and variability of the Sun and Sun-like stars. The longest available direct faculae/plage datasets are white-light facular and Ca II K observations going back to 1874 and 1892, respectively. Prior to that time, the only direct data available are for sunspots. We reassess the relationship between plage areas and sunspot records (areas and numbers) since 1892, to allow reconstructions of facular/plage areas which can be employed for studies going further back in time, i.e. over the period when solely sunspot observations are available. We use the plage areas derived from 38 consistently processed Ca II K archives as well as the plage area composite based on these archives. We find the relationship between plage and sunspot areas and numbers to be well represented by a power law function. We further find that the relationships depend on the bandwidth and the solar cycle strength. The reconstructions with a power law relationship are in good agreement with the original plage area series, whereas employment of a cycle-strength-dependent relationship improves the reconstructions only marginally. Performing the same analysis on other previously published plage area series, usually derived from a single archive with diverse processing techniques, returns different results when using different time series. This highlights the importance of applying a consistent processing to the various archives and demonstrates the uncertainties introduced by using previously published series. Our results have implications for past solar activity and irradiance reconstructions as well as for stellar activity studies, which sometimes assume a linear dependence between plage and sunspot areas.
△ Less
Submitted 15 September, 2022;
originally announced September 2022.
-
Reconstructing solar irradiance from historical Ca II K observations. I. Method and its validation
Authors:
Theodosios Chatzistergos,
Natalie A. Krivova,
Ilaria Ermolli,
Kok Leng Yeo,
Sudip Mandal,
Sami K. Solanki,
Greg Kopp,
Jean-Marie Malherbe
Abstract:
Knowledge of solar irradiance variability is critical to Earth's climate models and understanding the solar influence on Earth's climate. Direct solar irradiance measurements are only available since 1978. Reconstructions of past variability typically rely on sunspot data. These provide only indirect information on the facular and network regions, which are decisive contributors to irradiance vari…
▽ More
Knowledge of solar irradiance variability is critical to Earth's climate models and understanding the solar influence on Earth's climate. Direct solar irradiance measurements are only available since 1978. Reconstructions of past variability typically rely on sunspot data. These provide only indirect information on the facular and network regions, which are decisive contributors to irradiance variability on timescales of the solar cycle and longer. Our ultimate goal is to reconstruct past solar irradiance variations using historical full-disc Ca II K observations to describe the facular contribution independently of sunspot observations. Here, we develop the method and test it extensively by using modern CCD-based Ca II K observations and carry out initial tests on two photographic archives. We employ carefully reduced and calibrated Ca II K images from 13 datasets, such as those from the Meudon, Mt Wilson, and Rome observatories. We convert them to unsigned magnetograms and then use them as input to the adapted SATIRE model to reconstruct TSI variations over the period 1978-2019, for which direct irradiance measurements are available. The reconstructed TSI from the analysed Ca II K archives agrees well with direct TSI measurements and existing reconstructions. The model also returns good results on data taken with different bandpasses and images with low spatial resolution. Historical Ca II K archives suffer from numerous inconsistencies, but we show that these archives can still be used to reconstruct TSI with reasonable accuracy provided the observations are accurately processed. By using the unsigned magnetograms of the Sun reconstructed from high-quality Ca II K observations as input into the SATIRE model, we can reconstruct solar irradiance variations nearly as accurately as from directly recorded magnetograms.
△ Less
Submitted 13 September, 2021;
originally announced September 2021.
-
Modelling a multi-spacecraft coronal mass ejection encounter with EUHFORIA
Authors:
E. Asvestari,
J. Pomoell,
E. Kilpua,
S. Good,
T. Chatzistergos,
M. Temmer,
E. Palmerio,
S. Poedts,
J. Magdalenic
Abstract:
Coronal mass ejections (CMEs) are a manifestation of the Sun's eruptive nature. They can have a great impact on Earth, but also on human activity in space and on the ground. Therefore, modelling their evolution as they propagate through interplanetary space is essential. EUropean Heliospheric FORecasting Information Asset (EUHFORIA) is a data-driven, physics-based model, tracing the evolution of C…
▽ More
Coronal mass ejections (CMEs) are a manifestation of the Sun's eruptive nature. They can have a great impact on Earth, but also on human activity in space and on the ground. Therefore, modelling their evolution as they propagate through interplanetary space is essential. EUropean Heliospheric FORecasting Information Asset (EUHFORIA) is a data-driven, physics-based model, tracing the evolution of CMEs through background solar wind conditions. It employs a spheromak flux rope, which provides it with the advantage of reconstructing the internal magnetic field configuration of CMEs. This is something that is not included in the simpler cone CME model used so far for space weather forecasting. This work aims at assessing the spheromak CME model included in EUHFORIA. We employed the spheromak CME model to reconstruct a well observed CME and compare model output to in situ observations. We focus on an eruption from 6 January 2013 encountered by two radially aligned spacecraft, Venus Express and STEREO-A. We first analysed the observed properties of the source of this CME eruption and we extracted the CME properties as it lifted off from the Sun. Using this information, we set up EUHFORIA runs to model the event. The model predicts arrival times from half to a full day ahead of the in situ observed ones, but within errors established from similar studies. In the modelling domain, the CME appears to be propagating primarily southward, which is in accordance with white-light images of the CME eruption close to the Sun. In order to get the observed magnetic field topology, we aimed at selecting a spheromak rotation angle for which the axis of symmetry of the spheromak is perpendicular to the direction of the polarity inversion line (PIL). The modelled magnetic field profiles, their amplitude, arrival times, and sheath region length are all affected by the choice of radius of the modelled spheromak.
△ Less
Submitted 25 May, 2021;
originally announced May 2021.
-
Modeling stellar Ca II H & K emission variations. I. Effect of inclination on the S-index
Authors:
K. Sowmya,
A. I. Shapiro,
V. Witzke,
N. -E. Nèmec,
T. Chatzistergos,
K. L. Yeo,
N. A. Krivova,
S. K. Solanki
Abstract:
The emission in the near ultraviolet Ca II H & K lines is modulated by stellar magnetic activity. Although this emission, quantified via the S-index, has been serving as a prime proxy of stellar magnetic activity for several decades, many aspects of the complex relation between stellar magnetism and Ca II H & K emission are still unclear. The amount of measured Ca II H & K emission is suspected to…
▽ More
The emission in the near ultraviolet Ca II H & K lines is modulated by stellar magnetic activity. Although this emission, quantified via the S-index, has been serving as a prime proxy of stellar magnetic activity for several decades, many aspects of the complex relation between stellar magnetism and Ca II H & K emission are still unclear. The amount of measured Ca II H & K emission is suspected to be affected not only by the stellar intrinsic properties but also by the inclination angle of the stellar rotation axis. Until now such an inclination effect on S-index has remained largely unexplored. To fill this gap, we develop a physics-based model to calculate S-index, focusing on the Sun. Using the distributions of solar magnetic features derived from observations together with Ca II H & K spectra synthesized in non-local thermodynamic equilibrium, we validate our model by successfully reconstructing the observed variations of solar S-index over four activity cycles. Further, using the distribution of magnetic features over the visible solar disk obtained from surface flux transport simulations, we obtain S-index time series dating back to 1700 and investigate the effect of inclination on S-index variability, both on the magnetic activity cycle and the rotational timescales. We find that when going from an equatorial to a pole-on view, the amplitude of S-index variations decreases weakly on the activity cycle timescale and strongly on the rotational timescale (by about 22% and 81%, respectively, for a cycle of intermediate strength). The absolute value of S-index depends only weakly on the inclination. We provide analytical expressions that model such dependencies.
△ Less
Submitted 25 March, 2021;
originally announced March 2021.
-
Modelling solar irradiance from ground-based photometric observations
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Fabrizio Giorgi,
Natalie A. Krivova,
Cosmin Constantin Puiu
Abstract:
Total solar irradiance (TSI) has been monitored from space since 1978. The measurements show a prominent variability in phase with the solar cycle, as well as fluctuations on timescales shorter than a few days. However, the measurements were done by multiple and usually relatively short-lived missions making the possible long-term trend in the TSI highly uncertain. While the variability in the UV…
▽ More
Total solar irradiance (TSI) has been monitored from space since 1978. The measurements show a prominent variability in phase with the solar cycle, as well as fluctuations on timescales shorter than a few days. However, the measurements were done by multiple and usually relatively short-lived missions making the possible long-term trend in the TSI highly uncertain. While the variability in the UV irradiance is clearly in-phase with the solar cycle, the phase of the variability in the visible range has been debated. In this paper, we aim at getting an insight into the long-term trend of TSI since 1996 and the phase of the solar irradiance variations in the visible part of the spectrum. We use independent ground-based full-disc photometric observations in Ca~II~K and continuum from the Rome and San Fernando observatories to compute the TSI since 1996. We follow the empirical San Fernando approach based on the photometric sum index. We find a weak declining trend in the TSI of -7.8$^{+4.9}_{-0.8}\times10^{-3}$ Wm$^{-2}$y$^{-1}$ between the 1996 and 2008 activity minima, while between 2008 and 2019 the reconstructed TSI shows no trend to a marginally decreasing (but statistically insignificant) trend of -0.1$^{+0.25}_{-0.02}\times10^{-3}$ Wm$^{-2}$y$^{-1}$. The reference TSI series used for the reconstruction does not significantly affect the determined trend. The variation in the blue continuum (409.2 nm) is rather flat, while the variation in the red continuum (607.1 nm) is marginally in anti-phase, although this result is extremely sensitive to the accurate assessment of the quiet Sun level in the images. These results provide further insights into the long-term variation of the total solar irradiance. The amplitude of the variations in the visible is below the uncertainties of the processing, which prevents an assessment of the phase of the variations.
△ Less
Submitted 24 August, 2020;
originally announced August 2020.
-
Analysis of full disc Ca II K spectroheliograms III. Plage area composite series covering 1892-2019
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Natalie A. Krivova,
Sami K. Solanki,
Dipankar Banerjee,
Teresa Barata,
Marcel Belik,
Ricardo Gafeira,
Adriana Garcia,
Yoichiro Hanaoka,
Manjunath Hegde,
Jan Klimeš,
Viktor V. Korokhin,
Ana Lourenço,
Jean-Marie Malherbe,
Gennady P. Marchenko,
Nuno Peixinho,
Takashi Sakurai,
Andrey G. Tlatov
Abstract:
We derive the plage area evolution over the last 12 solar cycles employing data from all Ca II K archives available publicly in digital form known to us, including several as yet unexplored Ca II K archives. We analyse more than 290,000 full-disc Ca II K observations from 43 datasets spanning the period 1892-2019. All images were consistently processed with an automatic procedure that performs the…
▽ More
We derive the plage area evolution over the last 12 solar cycles employing data from all Ca II K archives available publicly in digital form known to us, including several as yet unexplored Ca II K archives. We analyse more than 290,000 full-disc Ca II K observations from 43 datasets spanning the period 1892-2019. All images were consistently processed with an automatic procedure that performs the photometric calibration (if needed) and the limb-darkening compensation. The processing also accounts for artefacts plaguing many of the images, including some very specific artefacts such as bright arcs found in Kyoto and Yerkes data. We have produced a plage area time-series from each analysed dataset. We found that the differences between the plage areas derived from individual archives are mainly due to the differences in the central wavelength and the bandpass used to acquire the data at the various sites. We have empirically cross-calibrated and combined the results obtained from each dataset to produce a composite series of plage areas. "Backbone" series are used to bridge all the series together. We have also shown that the selection of the backbone series has little effect on the final plage area composite. We have quantified the uncertainty of determining the plage areas with our processing due to shifts in the central wavelength and found it to be less than 0.01 in fraction of the solar disc for the average conditions found on historical data. We also found the variable seeing conditions during the observations to slightly increase the plage areas during activity maxima. We provide the so far most complete time series of plage areas based on corrected and calibrated historical and modern Ca II K images. Consistent plage areas are now available on 88% of all days from 1892 onwards and on 98% from 1907 onwards.
△ Less
Submitted 4 May, 2020;
originally announced May 2020.
-
Delving into the Historical Ca II K Archive from the Kodaikanal Observatory: the Potential of the Most Recent Digitised Series
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Sami K. Solanki,
Natalie A. Krivova,
Dipankar Banerjee,
Bibhuti K. Jha,
Subhamoy Chatterjee
Abstract:
Full-disc Ca II K photographic observations of the Sun carry direct information about the evolution of solar-plage regions for more than a century and are therefore a unique dataset for solar-activity studies. For a long time Ca II K observations were barely explored, but recent digitisations of multiple archives have allowed their extensive analysis. However, various studies have reported diverse…
▽ More
Full-disc Ca II K photographic observations of the Sun carry direct information about the evolution of solar-plage regions for more than a century and are therefore a unique dataset for solar-activity studies. For a long time Ca II K observations were barely explored, but recent digitisations of multiple archives have allowed their extensive analysis. However, various studies have reported diverse results partly due to the insufficient quality of the digitised data. Furthermore, inhomogeneities have been identified within the individual archives, which, at least partly, could be due to the digitisation. As a result, some of the archives, e.g. that from the Kodaikanal observatory, were re-digitised. The results obtained by different authors who analysed the data from the new digitisation of the Kodaikanal archive differ from each other as well as from those derived from the old digitisation. Since the data were processed and analysed using different techniques, it is not clear, however, whether the differences are due to the digitisation or the processing of the data. To understand the reasons for such discrepancies, we analyse here the data from the two most recent digitisations of this archive. We use the same techniques to consistently process the images from both archives and to derive the plage areas from them. Some issues have been identified in both digitisations, implying that they are intrinsic characteristics of the data. Moreover, errors in timing of the observations plague both digitisations. Overall, the most recent 16-bit digitisation offers an improvement over the earlier 8-bit one. It also includes considerably more data and should be preferred.
△ Less
Submitted 15 August, 2019;
originally announced August 2019.
-
Recovering the unsigned photospheric magnetic field from Ca II K observations
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Sami K. Solanki,
Natalie A. Krivova,
Fabrizio Giorgi,
Kok Leng Yeo
Abstract:
We reassess the relationship between the photospheric magnetic field strength and the Ca II K intensity for a variety of surface features as a function of the position on the disc and the solar activity level. This relationship can be used to recover the unsigned photospheric magnetic field from images recorded in the core of Ca II K line. We have analysed 131 pairs of high-quality, full-disc, nea…
▽ More
We reassess the relationship between the photospheric magnetic field strength and the Ca II K intensity for a variety of surface features as a function of the position on the disc and the solar activity level. This relationship can be used to recover the unsigned photospheric magnetic field from images recorded in the core of Ca II K line. We have analysed 131 pairs of high-quality, full-disc, near-co-temporal observations from SDO/HMI and Rome/PSPT spanning half a solar cycle. To analytically describe the observationally-determined relation, we considered three different functions: a power law with an offset, a logarithmic function, and a power law function of the logarithm of the magnetic flux density. We used the obtained relations to reconstruct maps of the line-of-sight component of the unsigned magnetic field (unsigned magnetograms) from Ca II K observations, which were then compared to the original magnetograms. We find that both power-law functions represent the data well, while the logarithmic function is good only for quiet periods. We see no significant variation over the solar cycle or over the disc in the derived fit parameters, independently of the function used. We find that errors in the independent variable, usually not accounted for, introduce attenuation bias. To address this, we binned the data with respect to the magnetic field strength and Ca II K contrast separately and derived the relation for the bisector of the two binned curves. The reconstructed unsigned magnetograms show good agreement with the original ones. RMS differences are less than 90 G. The results were unaffected by the stray-light correction of the SDO/HMI and Rome/PSPT data. Our results imply that Ca~II~K observations, accurately processed and calibrated, can be used to reconstruct unsigned magnetograms by using the relations derived in our study.
△ Less
Submitted 9 May, 2019;
originally announced May 2019.
-
Historical solar Ca II K observations at the Rome and Catania observatories
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Mariachiara Falco,
Fabrizio Giorgi,
Salvo L. Guglielmino,
Natalie A. Krivova,
Paolo Romano,
Sami K. Solanki
Abstract:
Here we present the little explored Ca II K archives from the Rome and the Catania observatories and analyse the digitised images from these archives to derive plage areas.
Here we present the little explored Ca II K archives from the Rome and the Catania observatories and analyse the digitised images from these archives to derive plage areas.
△ Less
Submitted 20 February, 2019;
originally announced February 2019.
-
Analysis of full disc Ca II K spectroheliograms. II. Towards an accurate assessment of long-term variations in plage areas
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Natalie A. Krivova,
Sami K. Solanki
Abstract:
Reconstructions of past irradiance variations require suitable data on solar activity. The longest direct proxy is the sunspot number, and it has been most widely employed for this purpose. These data, however, only provide information on the surface magnetic field emerging in sunspots, while a suitable proxy of the evolution of the bright magnetic features, specifically faculae/plage and network,…
▽ More
Reconstructions of past irradiance variations require suitable data on solar activity. The longest direct proxy is the sunspot number, and it has been most widely employed for this purpose. These data, however, only provide information on the surface magnetic field emerging in sunspots, while a suitable proxy of the evolution of the bright magnetic features, specifically faculae/plage and network, is missing. This information can potentially be extracted from the historical full-disc observations in the Ca II K line. We have analysed over 100,000 historical images from 8 digitised photographic archives of the Arcetri, Kodaikanal, McMath-Hulbert, Meudon, Mitaka, Mt Wilson, Schauinsland, and Wendelstein observatories, as well as one archive of modern observations from the Rome/PSPT. The analysed data cover the period 1893--2018. We first performed careful photometric calibration and compensation for the centre-to-limb variation, and then segmented the images to identify plage regions. This has been consistently applied to both historical and modern observations. The plage series derived from different archives are generally in good agreement with each other. However, there are also clear deviations that most likely hint at intrinsic differences in the data and their digitisation. We showed that accurate image processing significantly reduces errors in the plage area estimates. Accurate photometric calibration also allows precise plage identification on images from different archives without the need to arbitrarily adjust the segmentation parameters. Finally, by comparing the plage area series from the various records, we found the conversion laws between them. This allowed us to produce a preliminary composite of the plage areas obtained from all the datasets studied here. This is a first step towards an accurate assessment of the long-term variation of plage regions.
△ Less
Submitted 19 February, 2019;
originally announced February 2019.
-
The potential of Ca II K observations for solar activity and variability studies
Authors:
Ilaria Ermolli,
Theodosios Chatzistergos,
Natalie A. Krivova,
Sami K. Solanki
Abstract:
Several observatories around the globe started regular full-disc imaging of the solar atmosphere in the Ca II K line in the early decades of the 20th century. These observations are continued today at a few sites with either old spectroheliographs or modern telescopes equipped with narrow-band filters. The Ca II K time series are unique in representing long-term variations of the Sun's chromospher…
▽ More
Several observatories around the globe started regular full-disc imaging of the solar atmosphere in the Ca II K line in the early decades of the 20th century. These observations are continued today at a few sites with either old spectroheliographs or modern telescopes equipped with narrow-band filters. The Ca II K time series are unique in representing long-term variations of the Sun's chromospheric magnetic field. However, meaningful results from their analysis require accurate processing of the available data and robust merging of the information stored in different archives. This paper provides an overview of the historical and modern full-disc Ca II K observations, with focus on their quality and the main results obtained from their analysis over the last decade.
△ Less
Submitted 11 May, 2018;
originally announced May 2018.
-
Ca II K spectroheliograms for studies of long-term changes in solar irradiance
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Natalie A. Krivova,
Sami K. Solanki
Abstract:
We address the importance of historical full disc Ca II K spectroheliograms for solar activity and irradiance reconstruction studies. We review our work on processing such data to enable them to be used in irradiance reconstructions. We also present our preliminary estimates of the plage areas from five of the longest available historical Ca II K archives.
We address the importance of historical full disc Ca II K spectroheliograms for solar activity and irradiance reconstruction studies. We review our work on processing such data to enable them to be used in irradiance reconstructions. We also present our preliminary estimates of the plage areas from five of the longest available historical Ca II K archives.
△ Less
Submitted 10 May, 2018;
originally announced May 2018.
-
Analysis of full disc Ca II K spectroheliograms. I. Photometric calibration and CLV compensation
Authors:
Theodosios Chatzistergos,
Ilaria Ermolli,
Sami K. Solanki,
Natalie A. Krivova
Abstract:
Historical Ca II K spectroheliograms (SHG) are unique in representing long-term variations of the solar chromospheric magnetic field. They usually suffer from numerous problems and lack photometric calibration. Thus accurate processing of these data is required to get meaningful results from their analysis. In this paper we aim at developing an automatic processing and photometric calibration meth…
▽ More
Historical Ca II K spectroheliograms (SHG) are unique in representing long-term variations of the solar chromospheric magnetic field. They usually suffer from numerous problems and lack photometric calibration. Thus accurate processing of these data is required to get meaningful results from their analysis. In this paper we aim at developing an automatic processing and photometric calibration method that provides precise and consistent results when applied to historical SHG. The proposed method is based on the assumption that the centre-to-limb variation of the intensity in quiet Sun regions does not vary with time. We tested the accuracy of the proposed method on various sets of synthetic images that mimic problems encountered in historical observations. We also tested our approach on a large sample of images randomly extracted from seven different SHG archives. The tests carried out on the synthetic data show that the maximum relative errors of the method are generally <6.5%, while the average error is <1%, even if rather poor quality observations are considered. In the absence of strong artefacts the method returns images that differ from the ideal ones by <2% in any pixel. The method gives consistent values for both plage and network areas. We also show that our method returns consistent results for images from different SHG archives. Our tests show that the proposed method is more accurate than other methods presented in the literature. Our method can also be applied to process images from photographic archives of solar observations at other wavelengths than Ca II K.
△ Less
Submitted 13 October, 2017;
originally announced October 2017.
-
New reconstruction of the sunspot group number since 1739 using the direct calibration and "backbone" methods
Authors:
Theodosios Chatzistergos,
Ilya G. Usoskin,
Gennady A. Kovaltsov,
Natalie A. Krivova,
Sami K. Solanki
Abstract:
Group sunspot number (GSN) series constitute the longest instrumental astronomical database providing information on solar activity. It is a compilation of observations by many individual observers, and their inter-calibration has usually been performed using linear rescaling. There are multiple published series that show different long-term trends for solar activity. We aim at producing a GSN ser…
▽ More
Group sunspot number (GSN) series constitute the longest instrumental astronomical database providing information on solar activity. It is a compilation of observations by many individual observers, and their inter-calibration has usually been performed using linear rescaling. There are multiple published series that show different long-term trends for solar activity. We aim at producing a GSN series, with a non-linear non-parametric calibration. The only underlying assumptions are that the differences between the various series are due to different acuity thresholds of the observers, and that the threshold of each observer remains constant throughout the observing period. We use a daisy chain process with backbone (BB) observers and calibrate all overlapping observers to them. The calibration of each individual observer is performed with a probability distribution function (PDF) matrix, constructed considering all daily values for the overlapping period with the BB. The calibration of the BBs is done in a similar manner. Propagation of errors is modelled with Monte Carlo simulations. The final series extends back to 1739 and includes data from 314 observers. It suggests moderate activity during the 18th and 19th century, which is significantly lower than the high level of solar activity predicted by other recent reconstructions applying linear regressions. The new series provides a robust reconstruction, based on modern and non-parametric methods, of sunspot group numbers since 1739, and it confirms the existence of the Modern grand maximum of solar activity in the second half of the 20th century.
△ Less
Submitted 27 February, 2017; v1 submitted 20 February, 2017;
originally announced February 2017.
-
Dependence of the Sunspot-group Size on the Level of Solar Activity and its Influence on the Calibration of Solar Observers
Authors:
I. G. Usoskin,
G. A. Kovaltsov,
T. Chatzistergos
Abstract:
The distribution of the sunspot group size (area) and its dependence on the level of solar activity is studied. It is shown that the fraction of small groups is not constant but decreases with the level of solar activity so that high solar activity is largely defined by big groups. We study the possible influence of solar activity on the ability of a realistic observer to see and report the daily…
▽ More
The distribution of the sunspot group size (area) and its dependence on the level of solar activity is studied. It is shown that the fraction of small groups is not constant but decreases with the level of solar activity so that high solar activity is largely defined by big groups. We study the possible influence of solar activity on the ability of a realistic observer to see and report the daily number of sunspot groups. It is shown that the relation between the number of sunspot groups as seen by different observers with different observational acuity thresholds is strongly non-linear and cannot be approximated by the traditionally used linear scaling ($k-$factors). The observational acuity threshold [$A_{\rm th}$] is considered to quantify the quality of each observer, instead of the traditional relative $k-$factor. A nonlinear $c-$factor based on $A_{\rm th}$ is proposed, which can be used to correct each observer to the reference conditions. The method is tested on a pair of principal solar observers, Wolf and Wolfer, and it is shown that the traditional linear correction, with the constant $k-$factor of 1.66 to scale Wolf to Wolfer, leads to an overestimate of solar activity around solar maxima.
△ Less
Submitted 2 September, 2016;
originally announced September 2016.