Analysis of historical solar observations and long-term changes in solar irradiance

The Sun is the main external driver of Earth's climate. Various mechanisms of the solar influence on climate have been proposed. The debate is ongoing, but variation in the radiative flux of the Sun is among the main candidates. Direct measurements of the solar irradiance exist for merely 40 years, which is a rather short period to derive conclusions about any possible long term changes in solar irradiance and their possible influence on climate. The main driver of the irradiance variations on time-scales of days to decades, and possibly longer, is believed to be the solar surface magnetism. Models have been developed that reconstruct the irradiance by using appropriate proxies of the magnetic activity of the Sun. Irradiance models require input data representing both dark and bright magnetic features emerging at the solar surface. The most widely ever used proxies are the sunspot areas (available since 1874), the sunspot number (available since 1700), and the sunspot group number (available since 1610). However, these records do not provide direct information on bright features. Their evolution has to be inferred from the sunspot data via certain assumptions whose justification is still very unclear. Therefore, there is a strong need for a more direct facular proxy. %in order to resolve the ambiguities of the models. Ca II K full-disc spectroheliograms are uniquely suited for that purpose. Observations in the Ca II K line started as early as in 1892 at various sites, providing a good temporal coverage of the whole 20th century. However, these data suffer from a variety of problems hindering their immediate applicability. The historical Ca II K observations are stored in photographic plates or films, which have a non-linear response to the incident radiation. Information on this relation has not been recorded for the majority of the historical observations. Furthermore, a plethora of artefacts have been introduced on these images during their various life stages. These artefacts need to be accurately accounted for in order to provide meaningful results from such data. We have developed a method to recover the relation for the response of the plates to the incident radiation by using information that is stored on the solar disc of the image. This method is based on the assumption that the darker parts of the quiet Sun regions remain unchanged in time. We have also developed a method of calculating the quiet Sun background of the images, which includes the centre-to-limb variation and takes into account various large-scale artefacts. We have shown that the accuracy of this method is greater than that of previously proposed techniques. We have also reassessed the relation between the magnetic field strength and the Ca II K contrast, by using a larger number of Ca II images than was done in earlier such studies. We find that this relation can be well described with a power law function, and the best fit parameters are unaffected by the activity level or the position on the disc. Hence this relation potentially allows a reconstruction of pseudo-magnetograms from the available Ca II K observations covering almost the whole 20th century, that can be used by irradiance models. The sunspot data (i.e. records of sunspot number), despite their extensive use, are not free of problems either. The process of cross-calibrating different records by individual observers has recently become a matter of debate. This debate brought to light that the majority of the methods used so far fail to take into consideration the non-linearity that arises due to different observing capabilities of the observers. We addressed the issue of the shape of the relation for the inter-calibration between different group sunspot number series. We have shown, with the aid of synthetic observations derived from the royal Greenwich observatory sunspot area records, that it is strongly non-linear, contrary to what is commonly assumed. We have developed a method to recalibrate the sunspot group number series with a non-linear non-parametric method based on daily statistics between the observers. Using Monte Carlo simulations we have accounted for the error propagation, which has not been done by any previous reconstruction. Our reconstruction of the group sunspot number favours moderate activity levels for the 18th and 19th century and supports the existence of the modern grand maximum of solar activity

Ca II K spectroheliograms for studies of long-term changes in solar irradiance

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.Comment: 5 pages, 3 figure

Analysis of full disc Ca II K spectroheliograms. II. Towards an accurate assessment of long-term variations in plage areas

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.Comment: 30 pages, 22 figures, accepted in A&

The potential of Ca II K observations for solar activity and variability studies

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.Comment: 6 pages, 2 figure

Long-term changes in solar activity and irradiance

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.Comment: 17 pages, 8 figures, accepted for publication in JAST

Recovering the unsigned photospheric magnetic field from Ca II K observations

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.Comment: 18 pages, 22 figures, accepted in A&

Rome Precision Solar Photometric Telescope: precision solar full-disk photometry during solar cycles 23–25

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 1 min to several months, and from 1 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

Full-disc Ca II K observations -- a window to past solar magnetism

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

Historical solar Ca II K observations at the Rome and Catania observatories

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.Comment: 5 pages, 3 figures, to be published in "Nuovo Cimento C" as proceeding of the Third Meeting of the Italian Solar and Heliospheric Communit

Delving into the Historical Ca II K Archive from the Kodaikanal Observatory: the Potential of the Most Recent Digitised Series

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.Comment: 28 pages, 17 figures, accepted for publication in solar physic