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.

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.

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.

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.

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.

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.

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.

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.

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.

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.