Search

Ca II K observations of the Sun have a great potential for probing the Sun's magnetism and activity, as well as for reconstructing solar irradiance. The Kodaikanal Solar Observatory (KoSO) in India, houses one of the most prominent Ca II K archives, spanning from 1904 to 2007, obtained under the same experimental conditions over a century, a feat very few other sites have achieved. However, the KoSO Ca II K archive suffers from several inconsistencies (e.g., missing/incorrect timestamps of observations and orientation of some images) which have limited the use of the archive. This study is a step towards bringing the KoSO archive to its full potential. We did this by developing an automatic method to orient the images more accurately than in previous studies. Furthermore, we included more data than in earlier studies (considering images that could not previously be analyzed by other techniques as well as 2845 newly digitized images), while also accounting for mistakes in the observational date/time. These images were accurately processed to identify plage regions along with their locations, enabling us to construct the butterfly diagram of plage areas from the entire KoSO Ca II K archive covering 1904-2007. Our butterfly diagram shows significantly fewer data gaps compared to earlier versions due to the larger set of data used in this study. Moreover, our butterfly diagram is consistent with Sp\"orer's law for sunspots, validating our automatic image orientation method. Additionally, we found that the mean latitude of plage areas calculated over the entire period is 20.5%+/-2.0 higher than that of sunspots, irrespective of the phase or the strength of the solar cycle. We also studied the North-South asymmetry showing that the northern hemisphere dominated plage areas during solar cycles 19 and 20, while the southern hemisphere dominated during solar cycles 21--23., Comment: 15 pages, 9 figures, Under Review in Solar Physics

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

Stars appear darker at their limbs than at their disk centers because at the limb we are viewing the higher and cooler layers of stellar photospheres. Limb darkening derived from state-of-the-art stellar atmosphere models systematically fails to reproduce recent transiting exoplanet light curves from the Kepler, TESS, and JWST telescopes -- stellar brightness obtained from measurements drops less steeply towards the limb than predicted by models. All previous models assumed atmosphere devoid of magnetic fields. Here we use our new stellar atmosphere models computed with the 3D radiative magneto-hydrodynamic code MURaM to show that small-scale concentration of magnetic fields on the stellar surface affect limb darkening at a level that allows us to explain the observations. Our findings provide a way forward to improve the determination of exoplanet radii and especially the transmission spectroscopy analysis for transiting planets, which relies on a very accurate description of stellar limb darkening from the visible through the infrared. Furthermore, our findings imply that limb darkening allows measuring the small-scale magnetic field on stars with transiting planets., Comment: 19 pages, 6 figures This is the submission version of manuscript to Nature Astronomy

Full disc observations of the Sun in the H$\alpha$ 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$\alpha$ 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$\alpha$ observations. We have consistently processed observations from 15 H$\alpha$ 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$\alpha$ 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$\alpha$ 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., Comment: 21 pages, 14 figures, accepted for publication in A&A

Magnetic features on the surface of stars, such as spots and faculae, cause stellar spectral variability on time-scales of days and longer. For stars other than the Sun, the spectral signatures of faculae are poorly understood, limiting our ability to account for stellar pollution in exoplanet transit observations. Here we present the first facular contrasts derived from magnetoconvection simulations for K0, M0 and M2 main-sequence stars and compare them to previous calculations for G2 main-sequence stars. We simulate photospheres and immediate subsurface layers of main-sequence spectral types between K0 and M2, with different injected vertical magnetic fields (0 G, 100 G, 300 G and 500 G) using MURaM, a 3D radiation-magnetohydrodynamics code. We show synthetic spectra and contrasts from the UV (300 nm) to the IR (10000 nm) calculated using the ATLAS9 radiative transfer code. The calculations are performed for nine viewing angles to characterise the facular radiation across the disc. The brightness contrasts of magnetic regions are found to change significantly across spectral type, wavelength and magnetic field strength, leading to the conclusion that accurate contrasts cannot be found by scaling solar values. This is due to features of different size, apparent structure and spectral brightness emerging in the presence of a given magnetic field for different spectral types., Comment: 17 pages, 12 figures, to be published in MNRAS

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 JASTP

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., Comment: 47 pages, 13 figures, accepted for publication in Frontiers in Astronomy and Space Sciences

Over the past years, thousands of stellar flares have been detected by harvesting data from large photometric surveys. These detections, however, do not account for potential sources of contamination such as background stars appearing in the same aperture as the primary target. We present a new method for identifying the true flare sources in large photometric surveys using data from the Kepler mission. Potential flares are identified in two steps: first, we search the light curves for at least two subsequent data points exceeding a 5{\sigma} threshold above the running mean. For these two cadences, we subtract the "quiet" stellar flux from the Kepler pixel data to obtain new images where the potential flare is the main light source. In the second step, we use a Bayesian approach to fit the point spread function of the instrument to determine the most likely location of the flux excess on the detector. We applied our method to 5862 main-sequence stars with near-solar effective temperatures. We found 2274 events exceeding the 5-sigma in at least two consecutive points in the light curves. Applying the second step reduced this number to 342 superflares. Of these, 283 flares happened on 178 target stars, 47 events are associated with fainter background stars, and in 10 cases, the flare location cannot be distinguished between the target and a background star. We also present cases where flares have been reported previously but our technique could not attribute them to the target star. We conclude that 1) identifying outliers in the light curves alone is insufficient to attribute them to stellar flares and 2) flares can only be uniquely attributed to a certain star when the instrument pixel-level data together with the point spread function are taken into account. As a consequence, previous flare statistics are likely contaminated by instrumental effects and unresolved astrophysical sources., Comment: accepted for publication in Astronomy and Astrophysics

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., Comment: 22 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

The Solar Ultraviolet Imaging Telescope (SUIT) is an instrument onboard the Aditya-L1 mission of ISRO that will measure and monitor the solar radiation emitted in the near-ultraviolet wavelength range (200-400 nm). SUIT will simultaneously map the photosphere and the chromosphere of the Sun using 11 filters sensitive to different wavelengths and covering different heights in the solar atmosphere and help us understand the processes involved in the transfer of mass and energy from one layer to the other. SUIT will also allow us to measure and monitor spatially resolved solar spectral irradiance that governs the chemistry of oxygen and ozone in the stratosphere of Earth's atmosphere. This is central to our understanding of the Sun climate relationship., Comment: 5 Figures, Published in Current Science

Study Analysis Group 21 (SAG21) of NASA's Exoplanet Exploration Program Analysis Group (ExoPAG) was organized to study the effect of stellar contamination on space-based transmission spectroscopy, a method for studying exoplanetary atmospheres by measuring the wavelength-dependent radius of a planet as it transits its star. Transmission spectroscopy relies on a precise understanding of the spectrum of the star being occulted. However, stars are not homogeneous, constant light sources but have temporally evolving photospheres and chromospheres with inhomogeneities like spots, faculae, plages, granules, and flares. This SAG brought together an interdisciplinary team of more than 100 scientists, with observers and theorists from the heliophysics, stellar astrophysics, planetary science, and exoplanetary atmosphere research communities, to study the current research needs that can be addressed in this context to make the most of transit studies from current NASA facilities like HST and JWST. The analysis produced 14 findings, which fall into three Science Themes encompassing (1) how the Sun is used as our best laboratory to calibrate our understanding of stellar heterogeneities ("The Sun as the Stellar Benchmark"), (2) how stars other than the Sun extend our knowledge of heterogeneities ("Surface Heterogeneities of Other Stars") and (3) how to incorporate information gathered for the Sun and other stars into transit studies ("Mapping Stellar Knowledge to Transit Studies"). In this invited review, we largely reproduce the final report of SAG21 as a contribution to the peer-reviewed literature., Comment: Invited review in press at RASTI. Based on the ExoPAG SAG21 report (arXiv:2201.09905v1) and refined via feedback from three reviewers. 75 pages, 30 figures, 5 tables

Astrometry has long been a promising technique for exoplanet detection. At the theoretical limits, astrometry would allow for the detection of smaller planets than previously seen by current exoplanet search methods, but stellar activity may make these theoretical limits unreachable. Astrometric jitter of a Sun-like star due to magnetic activity in its photosphere induces apparent variability in the photocenter of order $0.5\ \textrm{m}R_\odot$. This jitter creates a fundamental astrophysical noise floor preventing detection of lower-mass planets in a single spectral band. By injecting planet orbits into simulated solar data at five different passbands, we investigate mitigation of this fundamental astrometric noise using correlations across passbands. For a true solar analog and a planet at 1 au semimajor axis, the $6\sigma$ detection limit set by stellar activity for an ideal telescope at the best single passband is $0.01$ Earth masses. We found that pairs of passbands with highly correlated astrometric jitter due to stellar activity, but with less motion in the redder band, enable higher-precision measurements of the common signal from the planet. Using this method improves detectable planet masses at 1 au by up to a factor of $10$, corresponding to at best $0.005$ Earth masses for a Sun-like star with a perfect telescope. Given these results, we recommend that future astrometry missions consider proceeding with two or more passbands to reduce noise due to stellar activity., Comment: 11 pages, 9 figures. Published in AJ

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., Comment: 25 pages, 15 figures, accepted for publication in A&A

Variability observed in photometric lightcurves of late-type stars (on timescales longer than a day) is a dominant noise source in exoplanet surveys and results predominantly from surface manifestations of stellar magnetic activity, namely faculae and spots. The implementation of faculae in lightcurve models is an open problem, with scaling typically based on spectra equivalent to hot stellar atmospheres or assuming a solar-derived facular contrast. We modelled rotational (single period) lightcurves of active G2, K0, M0 and M2 stars, with Sun-like surface distributions and realistic limb-dependent contrasts for faculae and spots. The sensitivity of lightcurve variability to changes in model parameters such as stellar inclination, feature area coverage, spot temperature, facular region magnetic flux density and active band latitudes is explored. For our lightcurve modelling approach we used actress, a geometrically accurate model for stellar variability. actress generates 2-sphere maps representing stellar surfaces and populates them with user-prescribed spot and facular region distributions. From this, lightcurves can be calculated at any inclination. Quiet star limb darkening and limb-dependent facular contrasts were derived from MURaM 3D magnetoconvection simulations using ATLAS9. 1D stellar atmosphere models were used for the spot contrasts. We applied actress in Monte-Carlo simulations, calculating lightcurve variability amplitudes in the Kepler band. We found that, for a given spectral type and stellar inclination, spot temperature and spot area coverage have the largest effect on variability of all simulation parameters. For a spot coverage of 1%, the typical variability of a solar-type star is around 2 parts-per-thousand. The presence of faculae clearly affects the mean brightness and lightcurve shape, but has relatively little influence on the variability., Comment: 17 pages, 16 figures. Accepted for publication in MNRAS 22 April 2021

Size distribution of sunspots provides key information about the generation and emergence processes of the solar magnetic field. Previous studies on the size distribution have primarily focused on either the whole group or individual spot areas. In this paper, we investigate the organization of spot areas within sunspot groups. In particular, we analyze the ratio, $\rm{R}$, of the area of the biggest spot ($\rm{A_{big\_spot}}$) inside a group, to the total area of that group ($\rm{A_{group}}$). We use sunspot observations from Kislovodsk, Pulkovo and Debrecen observatories, together covering solar cycles 17 to 24. We find that at the time when the group area reaches its maximum, the single biggest spot in a group typically occupies about 60% of the group area. For half of all groups, $\rm R$ lies in the range between roughly 50% and 70%. We also find R to change with the group area, $\rm{A_{group}}$, such that $\rm{R}$ reaches a maximum of about 0.65 for groups with $\rm{A_{group}}\approx 200\mu$Hem and then remains at about 0.6 for lager groups. Our findings imply a scale invariant emergence pattern, providing an observational constraint on the emergence process. Furthermore, extrapolation of our results to larger sunspot groups may have a bearing on the giant unresolved starspot features found in Doppler images of highly active sun-like stars. Our results suggest that such giant features are composed of multiple spots, with the largest spot occupying roughly 55--75% of the total group area (i.e. of the area of the giant starspots seen in Doppler images)., Comment: Accepted for publication in Astronomy and Astrophysics

The advent of Gaia, capable of measuring stellar wobbles caused by orbiting planets, raised an interest to the astrometric detection of exoplanets. Another source of such wobbles (often also called jitter) is stellar magnetic activity. A quantitative assessment of the stellar astrometric jitter is important for a more reliable astrometric detection and characterisation of exoplanets. We calculate the displacement of the solar photocentre due to the magnetic activity for an almost 16-year period (February 2, 1999 - August 1, 2014). We also investigate how the displacement depends on the spectral passband chosen for observations, including the wavelength range to be covered by the upcoming Small-JASMINE mission of JAXA. This is done by extending the SATIRE-S model for solar irradiance variability to calculating the displacement of the solar photocentre caused by the magnetic features on the surface of the Sun. We found that the peak to peak amplitude of the solar photocentre displacement would reach 0.5 mas if the Sun were located 10 pc away from the observer and observed in the Gaia G filter. This is by far too small to be detected by the Gaia mission. However, the Sun is a relatively inactive star so that one can expect significantly larger signals for younger, and, consequently, more active stars. The model developed in this study can be combined with the simulations of emergence and surface transport of magnetic flux which have recently became available to model the astrometric jitter over the broad range of magnetic activities., Comment: Accepted for publication in ApJ

Kepler observations revealed that hundreds of stars with near-solar fundamental parameters and rotation periods have much stronger and more regular brightness variations than the Sun. Here we identify one possible reason for the peculiar behaviour of these stars. Inspired by solar nests of activity, we assume that the degree of inhomogeneity of active-region (AR) emergence on such stars is higher than on the Sun. To test our hypothesis, we model stellar light curves by injecting ARs consisting of spots and faculae on stellar surfaces at various rates and nesting patterns, using solar AR properties and differential rotation. We show that a moderate increase of the emergence frequency from the solar value combined with the increase of the degree of nesting can explain the full range of observed amplitudes of variability of Sun-like stars with nearly the solar rotation period. Furthermore, nesting in the form of active longitudes, in which ARs tend to emerge in the vicinity of two longitudes separated by $180^\circ$, leads to highly regular, almost sine-like variability patterns, rather similar to those observed in a number of solar-like stars., Comment: Typo in Equation 1 corrected in this version. Erratum accepted by ApJL

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., Comment: 17 pages, 8 figures, accepted for publication in Journal of Space Weather and Space Climate as part of the topical issue "Space climate: The past and future of solar activity"

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., Comment: 23 pages, 19 figures, accepted for publication in A&A

Long and consistent sunspot area records are important for understanding the long-term solar activity and variability. Multiple observatories around the globe have regularly recorded sunspot areas, but such individual records only cover restricted periods of time. Furthermore, there are also systematic differences between them, so that these records need to be cross-calibrated before they can be reliably used for further studies. We produce a cross-calibrated and homogeneous record of total daily sunspot areas, both projected and corrected, covering the period between 1874 and 2019. A catalogue of calibrated individual group areas is also generated for the same period. We have compared the data from nine archives: Royal Greenwich Observatory (RGO), Kislovodsk, Pulkovo, Debrecen, Kodaikanal, Solar Optical Observing Network (SOON), Rome, Catania, and Yunnan Observatories, covering the period between 1874 and 2019. Mutual comparisons of the individual records have been employed to produce homogeneous and inter-calibrated records of daily projected and corrected areas. As in earlier studies, the basis of the composite is formed by the data from RGO. After 1976, the only datasets used are those from Kislovodsk, Pulkovo and Debrecen observatories. This choice was made based on the temporal coverage and the quality of the data. In contrast to the SOON data used in previous area composites for the post-RGO period, the properties of the data from Kislovodsk and Pulkovo are very similar to those from the RGO series. They also directly overlap the RGO data in time, which makes their cross-calibration with RGO much more reliable. We have also computed and provide the daily Photometric Sunspot Index (PSI) widely used, e.g., in empirical reconstructions of solar irradiance., Comment: Accepted for publication in Astronomy and Astrophysics. An additional file (animation) is available at https://www.dropbox.com/s/armawopcxt8kmb9/Mandal_2020_butterfly_diagram.zip?dl=0

The variability of the spectral solar irradiance (SSI) over the course of the 11-year solar cycle is one of the manifestations of solar magnetic activity. There is a strong evidence that the SSI variability has an effect on the Earth's atmosphere. The faster rotation of the Sun in the past lead to a more vigorous action of solar dynamo and thus potentially to larger amplitude of the SSI variability on the timescale of the solar activity cycle. This could led to a stronger response of the Earth's atmosphere as well as other solar system planets' atmospheres to the solar activity cycle. We calculate the amplitude of the SSI and TSI variability over the course of the solar activity cycle as a function of solar age. We employ the relationship between the stellar magnetic activity and the age based on observations of solar twins. Using this relation we reconstruct solar magnetic activity and the corresponding solar disk area coverages by magnetic features (i.e. spots and faculae) over the last four billion years. These disk coverages are then used to calculate the amplitude of the solar-cycle SSI variability as a function of wavelength and solar age. Our calculations show that the young Sun was significantly more variable than the present Sun. The amplitude of the solar-cycle Total Solar Irradiance (TSI) variability of the 600 Myr old Sun was about 10 times larger than that of the present Sun. Furthermore, the variability of the young Sun was spot-dominated (the Sun being brighter at the activity minimum than in the maximum), i.e. the Sun was overall brighter at activity minima than at maxima. The amplitude of the TSI variability decreased with solar age until it reached a minimum value at 2.8 Gyr. After this point, the TSI variability is faculae-dominated (the Sun is brighter at the activity maximum) and its amplitude increases with age., Comment: 8 pages, 6 figures, accepted publication for Astronomy&Astrophysics

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 physics

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&A

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 Community

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&A

In G dwarfs, the surface distribution, coverage and lifetimes of starspots deviate from solar-like patterns as the rotation rate increases. We set up a numerical platform which includes the large-scale rotational and surface flow effects, aiming to simulate evolving surface patterns over an activity cycle for up to 8 times the solar rotation and flux emergence rates. At the base of the convection zone, we assume a solar projected butterfly diagram. We then follow the rotationally distorted trajectories of rising thin flux tubes to obtain latitudes and tilt angles. Using them as source distributions, we run a surface flux transport model with solar parameters. Our model predicts surface distributions of the signed radial fields and the starspots that qualitatively agree with observations., Comment: 3 pages, 3 figures; submitted to Proceedings of the 20th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun; ed. S. Wolk. This (newest) version has a reference update

The latitudinal distribution of starspots deviates from the solar pattern with increasing rotation rate. Numerical simulations of magnetic flux emergence and transport can help model the observed stellar activity patterns and the associated brightness variations. We set up a composite model for the processes of flux emergence and transport on Sun-like stars, to simulate stellar brightness variations for various levels of magnetic activity and rotation rates. Assuming that the distribution of magnetic flux at the base of the convection zone follows solar scaling relations, we calculate the emergence latitudes and tilt angles of bipolar regions at the surface for various rotation rates, using thin-flux-tube simulations. Taking these two quantities as input to a surface flux transport SFT model, we simulate the diffusive-advective evolution of the radial field at the stellar surface, including effects of active region nesting. As the rotation rate increases, (1) magnetic flux emerges at higher latitudes and an inactive gap opens around the equator, reaching a half-width of $20^\circ$ for $8\Omega_\odot$, (2) the tilt angles of freshly emerged bipolar regions show stronger variations with latitude. Polar spots can form at $8\Omega_\odot$ by accumulation of follower-polarity flux from decaying bipolar regions. From $4\Omega_\odot$ to $8\Omega_\odot$, the maximum spot coverage changes from 3 to 20%, respectively, compared to 0.4% for the solar model. Nesting of activity can lead to strongly non-axisymmetric spot distributions. On Sun-like stars rotating at $8\Omega_\odot$ ($P_{\rm rot}\simeq 3$ days), polar spots can form, owing to higher levels of flux emergence rate and tilt angles. Defining spots by a threshold field strength yields global spot coverages that are roughly consistent with stellar observations., Comment: 16 pages, 13 figures. Astron. & Astrophys. (in press); minor language corrections made

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 figures

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 figures

Observations of Sun-like stars over the last half-century have improved our understanding of how magnetic dynamos, like that responsible for the 11-year solar cycle, change with rotation, mass and age. Here we show for the first time how metallicity can affect a stellar dynamo. Using the most complete set of observations of a stellar cycle ever obtained for a Sun-like star, we show how the solar analog HD 173701 exhibits solar-like differential rotation and a 7.4-year activity cycle. While the duration of the cycle is comparable to that generated by the solar dynamo, the amplitude of the brightness variability is substantially stronger. The only significant difference between HD 173701 and the Sun is its metallicity, which is twice the solar value. Therefore, this provides a unique opportunity to study the effect of the higher metallicity on the dynamo acting in this star and to obtain a comprehensive understanding of the physical mechanisms responsible for the observed photometric variability. The observations can be explained by the higher metallicity of the star, which is predicted to foster a deeper outer convection zone and a higher facular contrast, resulting in stronger variability., Comment: Submitted to ApJ

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., Comment: 24 pages, 24 figures, accepted for publication in Astronomy and Astrophysics

Context. Stellar spectral variability on timescales of a day and longer, arising from magnetic surface features such as dark spots and bright faculae, is an important noise source when characterising extra-solar planets. Current 1D models of faculae do not capture the geometric properties and fail to reproduce observed solar facular contrasts. Magnetoconvection simulations provide facular contrasts accounting for geometry. Aims. We calculate facular contrast spectra from magnetoconvection models of the solar photosphere with a view to improve (a) future parameter determinations for planets with early G type host stars and (b) reconstructions of solar spectral variability. Methods. Regions of a solar twin (G2, log g=4.44) atmosphere with a range of initial average vertical magnetic fields (100 to 500~G) were simulated using a 3D radiation-magnetohydrodynamics code, MURaM, and synthetic intensity spectra were calculated from the ultraviolet (149.5~nm) to the far infrared (160000~nm) with the ATLAS9 radiative transfer code. Nine viewing angles were investigated to account for facular positions across most of the stellar disc. Results. Contrasts of the radiation from simulation boxes with different levels of magnetic flux relative to an atmosphere with no magnetic field are a complicated function of position, wavelength and magnetic field strength that is not reproduced by 1D facular models. Generally, contrasts increase towards the limb, but at UV wavelengths a saturation and decrease are observed close to the limb. Contrasts also increase strongly from the visible to the UV; there is a rich spectral dependence, with marked peaks in molecular bands and strong spectral lines. At disc centre, a complex relationship with magnetic field was found and areas of strong magnetic field can appear either dark or bright, depending on wavelength..., Comment: 17 pages, 11 figures

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., Comment: 19 pages, 12 figures, Astronomy and Astrophysics, accepted

Aims: Although the time of the Maunder minimum (1645--1715) is widely known as a period of extremely low solar activity, claims are still debated that solar activity during that period might still have been moderate, even higher than the current solar cycle #24. We have revisited all the existing pieces of evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods: We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. Results: The level of solar activity during the Maunder minimum is reassessed on the basis of all available data sets. Conclusions: We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely below that during the Dalton minimum around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected at a high confidence level., Comment: Accepted to Astron. Astrophys

We present a revised and extended total and spectral solar irradiance (SSI) reconstruction, which includes a wavelength-dependent uncertainty estimate, spanning the last three solar cycles using the SATIRE-S model. The SSI reconstruction covers wavelengths between 115 and 160,000 nm and all dates between August 1974 and October 2009. This represents the first full-wavelength SATIRE-S reconstruction to cover the last three solar cycles without data gaps and with an uncertainty estimate. SATIRE-S is compared with the NRLSSI model and SORCE/SOLSTICE ultraviolet (UV) observations. SATIRE-S displays similar cycle behaviour to NRLSSI for wavelengths below 242 nm and almost twice the variability between 242 and 310 nm. During the decline of last solar cycle, between 2003 and 2008, SSI from SORCE/SOLSTICE version 12 and 10 typically displays more than three times the variability of SATIRE-S between 200 and 300 nm. All three datasets are used to model changes in stratospheric ozone within a 2D atmospheric model for a decline from high solar activity to solar minimum. The different flux changes result in different modelled ozone trends. Using NRLSSI leads to a decline in mesospheric ozone, while SATIRE-S and SORCE/SOLSTICE result in an increase. Recent publications have highlighted increases in mesospheric ozone when considering version 10 SORCE/SOLSTICE irradiances. The recalibrated SORCE/SOLSTICE version 12 irradiances result in a much smaller mesospheric ozone response than when using version 10 and now similar in magnitude to SATIRE-S. This shows that current knowledge of variations in spectral irradiance is not sufficient to warrant robust conclusions concerning the impact of solar variability on the atmosphere and climate., Comment: 25 pages (18 pages in main article with 6 figures; 7 pages in supplementary materials with 6 figures) in draft mode using the American Meteorological Society package. Submitted to Journal of Atmospheric Sciences for publication

We analyze and model total solar irradiance variability on time scales from minutes to months, excluding variations due to p-mode oscillations, using a combination of convective and magnetic components. These include granulation, the magnetic network, faculae and sunspots. Analysis of VIRGO data shows that on periods of a day or longer solar variability depends on magnetic activity, but is nearly independent at shorter periods. We assume that only granulation affects the solar irradiance variability on time scales from minutes to hours. Granulation is described as a large sample of bright cells and dark lanes that evolve according to rules deduced from observations and radiation hydrodynamic simulations. Comparison of this model combined with a high time resolution magnetic-field based irradiance reconstruction, with solar data reveals a good correspondence except at periods of 10 to 30 hours. This suggests that the model is missing some power at these periods, which may be due to the absence of supergranulation or insufficient sensitivity of MDI magnetograms used for the reconstruction of the magnetic field-based irradiance reconstructions. Our model also shows that even for spatially unresolved data (such as those available for stars) the Fourier or wavelet transform of time series sampled at high cadence may allow properties of stellar granulation, in particular granule lifetimes to be determined.

Solar irradiance models that assume solar irradiance variations to be due to changes in the solar surface magnetic flux have been successfully used to reconstruct total solar irradiance on rotational as well as cyclical and secular time scales. Modelling spectral solar irradiance is not yet as advanced, and also suffers from a lack of comparison data, in particular on solar-cycle time scales. Here we compare solar irradiance in the 220 nm to 240 nm band as modelled with SATIRE-S and measured by different instruments on the UARS and SORCE satellites. We find good agreement between the model and measurements on rotational time scales. The long-term trends, however, show significant differences. Both SORCE instruments, in particular, show a much steeper gradient over the decaying part of cycle 23 than the modelled irradiance or that measured by UARS/SUSIM., Comment: 8 pages, 2 figures, conference proceedings to appear in Surveys in Geophysics

Aims: We investigate how well modeled solar irradiances agree with measurements from the SORCE satellite, both for total solar irradiance and broken down into spectral regions on timescales of several years. Methods: We use the SATIRE model and compare modeled total solar irradiance (TSI) with TSI measurements between 2003 and 2009. Spectral solar irradiance over 200-1630nm is compared with the SIM instrument on SORCE between 2004 and 2009 during a period of decline from moderate activity to the recent solar minimum in 10 nm bands and for three spectral regions of significant interest: the UV integrated over 200-300nm, the visible over 400-691nm and the IR between 972-1630 nm. Results: The model captures 97% of observed TSI variation. In the spectral comparison, rotational variability is well reproduced, especially between 400 and 1200 nm. The magnitude of change in the long-term trends is many times larger in SIM at almost all wavelengths while trends in SIM oppose SATIRE in the visible between 500 and 700nm and between 1000 and 1200nm. We discuss the remaining issues with both SIM data and the identified limits of the model, particularly with the way facular contributions are dealt with, the limit of flux identification in MDI magnetograms during solar minimum and the model atmospheres in the IR employed by SATIRE. It is unlikely that improvements in these areas will significantly enhance the agreement in the long-term trends. This disagreement implies that some mechanism other than surface magnetism is causing SSI variations, in particular between 2004 and 2006, if the SIM data are correct. Since SATIRE was able to reproduce UV irradiance between 1991 and 2002 from UARS, either the solar mechanism for SSI variation fundamentally changed around the peak of cycle 23, or there is an inconsistency between UARS and SORCE UV measurements. We favour the second explanation., Comment: 14 pages, 13 figures

Aims. We present a physically consistent reconstruction of the total solar irradiance for the Holocene. Methods. We extend the SATIRE models to estimate the evolution of the total (and partly spectral) solar irradiance over the Holocene. The basic assumption is that the variations of the solar irradiance are due to the evolution of the dark and bright magnetic features on the solar surface. The evolution of the decadally averaged magnetic flux is computed from decadal values of cosmogenic isotope concentrations recorded in natural archives employing a series of physics-based models connecting the processes from the modulation of the cosmic ray flux in the heliosphere to their record in natural archives. We then compute the total solar irradiance (TSI) as a linear combination of the jth and jth + 1 decadal values of the open magnetic flux. Results. Reconstructions of the TSI over the Holocene, each valid for a di_erent paleomagnetic time series, are presented. Our analysis suggests that major sources of uncertainty in the TSI in this model are the heritage of the uncertainty of the TSI since 1610 reconstructed from sunspot data and the uncertainty of the evolution of the Earth's magnetic dipole moment. The analysis of the distribution functions of the reconstructed irradiance for the last 3000 years indicates that the estimates based on the virtual axial dipole moment are significantly lower at earlier times than the reconstructions based on the virtual dipole moment. Conclusions. We present the first physics-based reconstruction of the total solar irradiance over the Holocene, which will be of interest for studies of climate change over the last 11500 years. The reconstruction indicates that the decadally averaged total solar irradiance ranges over approximately 1.5 W/m2 from grand maxima to grand minima.

Aims: We test the reliability of the observed and calculated spectral irradiance variations between 200 and 1600 nm over a time span of three solar rotations in 2004. Methods: We compare our model calculations to spectral irradiance observations taken with SORCE/SIM, SoHO/VIRGO and UARS/SUSIM. The calculations assume LTE and are based on the SATIRE (Spectral And Total Irradiance REconstruction) model. We analyse the variability as a function of wavelength and present time series in a number of selected wavelength regions covering the UV to the NIR. We also show the facular and spot contributions to the total calculated variability. Results: In most wavelength regions, the variability agrees well between all sets of observations and the model calculations. The model does particularly well between 400 and 1300 nm, but fails below 220 nm as well as for some of the strong NUV lines. Our calculations clearly show the shift from faculae-dominated variability in the NUV to spot-dominated variability above approximately 400 nm. We also discuss some of the remaining problems, such as the low sensitivity of SUSIM and SORCE for wavelengths between approximately 310 and 350 nm, where currently the model calculations still provide the best estimates of solar variability., Comment: 15 pages, 11 figures, accepted by A&A