This is our reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars" by T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gomez. We hope that both the comment and our reply lead to fruitful discussions which of the two presented scenarios is more likely.

How global faculae and network coverage relates to that of sunspots is relevant to the brightness variations of the Sun and Sun-like stars. We extend earlier studies that found the facular-to-sunspot-area ratio diminishes with total sunspot coverage. Chromospheric indices and the total magnetic flux enclosed in network and faculae, referred to here as 'facular indices', are modulated by the amount of facular and network present. We probed the relationship between various facular and sunspot indices through an empirical model that takes into account how active regions evolve. This model was incorporated into a total solar irradiance (TSI) model. The model presented here replicates most of the observed variability in the facular indices, and is better at doing so than earlier models. Contrary to recent studies, we found the relationship between the facular and sunspot indices to be stable over the past four decades. The model indicates that, like the facular-to-sunspot-area ratio, the ratio of the variation in chromospheric emission and total network and facular magnetic flux to sunspot area decreases with the latter. The TSI model indicates the ratio of the TSI excess from faculae and network to the deficit from sunspots also declines with sunspot area, with the consequence being that TSI rises with sunspot area more slowly than if the two quantities were linearly proportional to one another. The extrapolation of the TSI model to higher activity levels indicates that in the activity range where Sun-like stars are observed to switch from growing brighter with increasing activity to becoming dimmer instead, the activity-dependence of TSI exhibits a similar transition as sunspot darkening starts to rise more rapidly with activity than facular brightening. This bolsters the interpretation of this behavior of Sun-like stars as the transition from a faculae-dominated to a spot-dominated regime.

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.

Magnetic activity of the Sun and other stars causes their brightness to vary. We investigate how typical the Sun's variability is compared to other solar-like stars, i.e. those with near-solar effective temperatures and rotation periods. By combining four years of photometric observations from the Kepler space telescope with astrometric data from the Gaia spacecraft, we measure photometric variabilities of 369 solar-like stars. Most of the solar-like stars with well-determined rotation periods show higher variability than the Sun and are therefore considerably more active. These stars appear nearly identical to the Sun, except for their higher variability. Their existence raises the question of whether the Sun can also experience epochs of such high variability.

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.

Comparing solar and stellar brightness variations is hampered by the difference in spectral passbands used in observations as well as by the possible difference in the inclination of their rotation axes from the line of sight. We calculate the rotational variability of the Sun as it would be measured in passbands used for stellar observations. In particular, we consider the filter systems used by the CoRoT, $Kepler$, TESS, and $Gaia$ space missions. We also quantify the effect of the inclination of the rotation axis on the solar rotational variability. We employ the Spectral And Total Irradiance REconstructions (SATIRE) model to calculate solar brightness variations in different filter systems as observed from the ecliptic plane. We then combine the simulations of the surface distribution of the magnetic features at different inclinations using a surface flux transport model (SFTM) with the SATIRE calculations to compute the dependence of the variability on the inclination. For an ecliptic-bound observer, the amplitude of the solar rotational variability, as observed in the total solar irradiance (TSI) is 0.68 mmag (averaged over solar cycles 21-24). We obtained corresponding amplitudes in the $Kepler$ (0.74 mmag), CoRoT (0.73 mmag), TESS (0.62 mmag), $Gaia~ $ (0.74 mmag), $Gaia~ G_{RP}$ (0.62 mmag), and ), $Gaia~ G_{BP}$ (0.86 mmag) passbands. Decreasing the inclination of the rotation axis decreases the rotational variability. For a sample of randomly inclined stars, the variability is on average 15% lower in all filter systems considered in this work. This almost compensates for the difference in the amplitudes of the variability in TSI and $Kepler$ passbands, making the amplitudes derived from the TSI records an ideal representation of the solar rotational variability for comparison to $Kepler$ stars with unknown inclinations.

Magnetic features on the surfaces of cool stars cause variations of their brightness. Such variations have been extensively studied for the Sun. Recent planet-hunting space telescopes allowed measuring brightness variations in hundred thousands of other stars. The new data posed the question of how typical is the Sun as a variable star. Putting solar variability into the stellar context suffers, however, from the bias of solar observations being made from its near-equatorial plane, whereas stars are observed at all possible inclinations. We model solar brightness variations at timescales from days to years as they would be observed at different inclinations. In particular, we consider the effect of the inclination on the power spectrum of solar brightness variations. The variations are calculated in several passbands routinely used for stellar measurements. We employ the Surface Flux Transport Model (SFTM) to simulate the time-dependent spatial distribution of magnetic features on both near- and far-sides of the Sun. This distribution is then used to calculate solar brightness variations following the SATIRE (Spectral And Total Irradiance REconstruction) approach. We have quantified the effect of the inclination on solar brightness variability at timescales down to a day. Thus, our results allow making solar brightness records directly comparable to those obtained by the planet-hunting space telescopes. Furthermore, we decompose solar brightness variations into the components originating from the solar rotation and from the evolution of magnetic features.

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.

Young and active stars generally have regular, almost sinusoidal, patterns of variability attributed to their rotation, while the majority of older and less active stars, including the Sun, have more complex and non-regular light-curves which do not have clear rotational-modulation signals. Consequently, the rotation periods have been successfully determined only for a small fraction of the Sun-like stars observed by transit-based planet-hunting missions, such as CoRoT, Kepler, and TESS. This suggests that only a small fraction of such systems have been properly identified as solar-like analogs. We apply a new method for determining rotation periods of low-activity stars, like the Sun. The method is based on calculating the gradient of the power spectrum (GPS) of stellar brightness variations and identifying a tell-tale inflection point in the spectrum. The rotation frequency is then proportional to the frequency of that inflection point. In this paper test this GPS method against available photometric records of the Sun. We apply GPS, autocorrelation functions, Lomb-Scargle periodograms, and wavelet analyses to the total solar irradiance (TSI) time series obtained from the Total Irradiance Monitor (TIM) on the Solar Radiation and Climate Experiment (SORCE) and the Variability of solar IRradiance and Gravity Oscillations (VIRGO) experiment on the SOlar and Heliospheric Observatory (SoHO) missions. We analyse the performance of all methods at various levels of solar activity. We show that the GPS method returns accurate values of solar rotation independently of the level of solar activity. In particular, it performs well during periods of high solar activity, when TSI variability displays an irregular pattern and other methods fail. Our results suggest that the GPS method can successfully determine the rotational periods of stars with both regular and non-regular light-curves.

Accurate determination of stellar rotation periods is important for estimating stellar ages as well as for understanding stellar activity and evolution. While for about thirty thousand stars in the Kepler field rotation periods can be determined, there are over hundred thousand stars, especially with low photometric variability and irregular pattern of variations, for which rotational periods are unknown. Here, we investigate the effect of metallicity on the detectability of rotation periods. This is done by synthesising light curves of hypothetical stars, which are identical to our Sun, with the exception of the metallicity. These light curves are then used as an input to the period determination algorithms. We find that the success rate for recovering the rotation signal has a minimum close to the solar metallicity value. This can be explained by the compensation effect of facular and spot contributions. In addition, selecting solar-like stars with near-solar effective temperature, near solar photometric variability, and with metallicity between M/H = -0.35 and M/H = 0.35 from the Kepler sample, we analyse the fraction of stars for which rotational periods have been detected as a function of metallicity. In agreement with our theoretical estimate we found a local minimum for the detection fraction close to the solar metallicity. We further report rotation periods of 87~solar-like Kepler stars for the first time.

Context. Considerable effort has been put into using light curves observed by space telescopes such as CoRoT, Kepler and TESS for determining stellar rotation periods. While rotation periods of active stars can be reliably determined, the light curves of many older and less active stars (e.g. stars similar to the Sun) are quite irregular, which hampers determination of their periods. Aims. We examine the factors causing the irregularities in stellar brightness variations and develop a method for determining rotation periods of low activity stars with irregular light curves. Methods. We extend the Spectral And Total Irradiance Reconstruction (SATIRE) approach for modelling solar brightness variations to Sun-like stars. We calculate the power spectra of stellar brightness variations for various combinations of parameters defining the surface configuration and evolution of stellar magnetic features. Results. The short lifetime of spots in comparison to the stellar rotation period as well as the interplay between spot and facular contributions to brightness variations of stars with near solar activity cause irregularities in their light curves. The power spectra of such stars often lack a peak associated with the rotation period. Nevertheless, the rotation period can still be determined by measuring the period where the concavity of the power spectrum plotted in the log-log scale changes sign, i.e. by identifying the position of the inflection point. Conclusions. The inflection point of the (log-log) power spectrum is found to be a new diagnostic for stellar rotation periods that is shown to work even in cases where the power spectrum shows no peak at the rotation rate.

Context. Solar spectral irradiance (SSI) variability is one of the key inputs to models of the Earth's climate. Understanding solar irradiance fluctuations also helps to place the Sun among other stars in terms of their brightness variability patterns and to set detectability limits for terrestrial exo-planets. Aims. One of the most successful and widely used models of solar irradiance variability is SATIRE-S. It uses spectra of the magnetic features and surrounding quiet Sun computed with the ATLAS9 spectral synthesis code under the assumption of Local Thermodynamic Equilibrium (LTE). SATIRE-S has been at the forefront of solar variability modelling, but due to the limitations of the LTE approximation its output SSI has to be empirically corrected below 300 nm, which reduces the physical consistency of its results. This shortcoming is addressed in the present paper. Methods. We replace the ATLAS9 spectra of all atmospheric components in SATIRE-S with the spectra calculated using the non-LTE Spectral Synthesis Code (NESSY). We also use Fontenla et al. (1999) temperature and density stratification models of the solar atmosphere to compute the spectrum of the quiet Sun and faculae. Results. We compute non-LTE contrasts of spots and faculae and combine them with the SDO/HMI filling factors of the active regions to calculate the total and spectral solar irradiance variability during solar cycle 24. Conclusions. The non-LTE contrasts result in total and spectral solar irradiance in good agreement with the empirically corrected output of the LTE version. This suggests that empirical correction introduced into SATIRE-S output is well judged and that the corrected total and spectral solar irradiance obtained from the SATIRE-S model in LTE is fully consistent with the results of non-LTE computations.

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.

Contex. Stellar spectra synthesis is essential for the characterization of potential planetary hosts. In addition, comprehensive stellar variability calculations with fast radiative transfer are needed to disentangle planetary transits from stellar magnetically-driven variability. The planet-hunting space telescopes, such as CoRoT, Kepler, and TESS will bring vast quantities of data, rekindling the interest in fast calculations of the radiative transfer. Aims. We revisit the Opacity Distribution Functions (ODF) approach routinely applied to speedup stellar spectral synthesis. To achieve a considerable speedup relative to the current state-of-the-art, we further optimize the approach and search for the best ODF configuration. Furthermore, we generalize the ODF approach for fast calculations of flux in various filters often used in stellar observations. Methods. In a parameter-sweep-fashion, we generated ODF in the spectral range from UV to IR with different setups. The most accurate ODF configuration for each spectral interval was determined. For calculations of the radiative fluxes through filters we adapted the wavelength grid based on the transmission curve, whereafter the normal ODF procedure was performed. Results. Our optimum ODF configuration allows for a three fold speedup, compared to the previously used ODF configurations. The ODF generalization to calculate fluxes through filters results in a speedup of more than two orders of magnitude.

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.

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.

Changes in solar irradiance and in its spectral distribution are among the main natural drivers of the climate on Earth. However, irradiance measurements are only available for less than four decades, while assessment of solar influence on Earth requires much longer records. The aim of this work is to provide the most up-to-date physics-based reconstruction of the solar total and spectral irradiance (TSI/SSI) over the last nine millennia. The concentrations of the cosmogenic isotopes 14C and 10Be in natural archives have been converted to decadally averaged sunspot numbers through a chain of physics-based models. TSI and SSI are reconstructed with an updated SATIRE model. Reconstructions are carried out for each isotope record separately, as well as for their composite. We present the first ever SSI reconstruction over the last 9000 years from the individual 14C and 10Be records as well as from their newest composite. The reconstruction employs physics-based models to describe the involved processes at each step of the procedure. Irradiance reconstructions based on two different cosmogenic isotope records, those of 14C and 10Be, agree well with each other in their long-term trends despite their different geochemical paths in the atmosphere of Earth. Over the last 9000 years, the reconstructed secular variability in TSI is of the order of 0.11%, or 1.5 W/m2. After the Maunder minimum, the reconstruction from the cosmogenic isotopes is consistent with that from the direct sunspot number observation. Furthermore, over the nineteenth century, the agreement of irradiance reconstructions using isotope records with the reconstruction from the sunspot number by Chatzistergos et al. (2017) is better than that with the reconstruction from the WDC-SILSO series (Clette et al. 2014), with a lower chi-square-value.

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.

Context. Comparison studies of Sun-like stars with the Sun suggest an anomalously low photometric variability of the Sun compared to Sun-like stars with similar magnetic activity. Comprehensive understanding of stellar variability is needed, to find a physical reasoning for this observation. Aims. We investigate the effect of metallicity and effective temperature on the photometric brightness change of Sun-like stars seen at different inclinations. The considered range of fundamental stellar parameters is sufficiently small so the stars, investigated here, still count as Sun-like or even as solar twins. Methods. To model the brightness change of stars with solar magnetic activity, we extend a well established model of solar brightness variations, SATIRE (which stands for Spectral And Total Irradiance Reconstruction), which is based on solar spectra, to stars with different fundamental parameters. For that we calculate stellar spectra for different metallicities and effective temperature using the radiative transfer code ATLAS9. Results. We show that even a small change (e.g. within the observational error range) of metallicity or effective temperature significantly affects the photometric brightness change compared to the Sun. We find that for Sun-like stars, the amplitude of the brightness variations obtained for Str\"omgren (b + y)/2 reaches a local minimum for fundamental stellar parameters close to the solar metallicity and effective temperature. Moreover, our results show that the effect of inclination decreases for metallicity values greater than the solar metallicity. Overall, we find that an exact determination of fundamental stellar parameters is crucially important for understanding stellar brightness changes.

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.

Solar activity in the past millennia can only be reconstructed from cosmogenic radionuclide records in terrestrial archives. However, because of the diversity of the proxy archives, it is difficult to build a homogeneous reconstruction. Here we provide a new consistent multiproxy reconstruction of the solar activity over the last 9000 years, using available long-span datasets of 10Be and 14C in terrestrial archives. A new method, based on a Bayesian approach, was applied for the first time to solar activity reconstruction. A Monte Carlo search for the most probable value of the modulation potential was performed to match data from different datasets for a given time. We used six 10Be series from Greenland and Antarctica, and the global 14C production series. The 10Be series were resampled to match wiggles related to the grand minima in the 14C reference dataset. The GRIP and the EDML 10Be series diverge from each other during the second half of the Holocene, while the 14C series lies between them. A likely reason for this is the insufficiently precise beryllium transport and deposition model for Greenland. A slow 6-millennia variability with lows at ca. 5500 BC and 1500 AD of solar activity is found. Two components of solar activity can be statistically distinguished: the main 'normal' component and a component corresponding to grand minima. A possible existence of a component representing grand maxima is indicated, but it cannot be separated from the main component in a statistically significant manner. A new consistent reconstruction of solar activity over the last nine millennia is presented with the most probable values of decadal sunspot numbers and their realistic uncertainties. Independent components of solar activity corresponding to the main moderate activity and the grand-minimum state are identified; they may be related to different operation modes of the dynamo.

The solar brightness varies on timescales from minutes to decades. Determining the sources of such variations, often referred to as solar noise, is of importance for multiple reasons: a) it is the background that limits the detection of solar oscillations, b) variability in solar brightness is one of the drivers of the Earth's climate system, c) it is a prototype of stellar variability which is an important limiting factor for the detection of extra-solar planets. Here we show that recent progress in simulations and observations of the Sun makes it finally possible to pinpoint the source of the solar noise. We utilise high-cadence observations from the Solar Dynamic Observatory and the SATIRE model to calculate the magnetically-driven variations of solar brightness. The brightness variations caused by the constantly evolving cellular granulation pattern on the solar surface are computed with the MURAM code. We find that surface magnetic field and granulation can together precisely explain solar noise on timescales from minutes to decades, i.e. ranging over more than six orders of magnitude in the period. This accounts for all timescales that have so far been resolved or covered by irradiance measurements. We demonstrate that no other sources of variability are required to explain the data. Recent measurements of Sun-like stars by CoRoT and Kepler uncovered brightness variations similar to that of the Sun but with much wider variety of patterns. Our finding that solar brightness variations can be replicated in detail with just two well-known sources will greatly simplify future modelling of existing CoRoT and Kepler as well as anticipated TESS and PLATO data.

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.

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.

The variation in the radiative output of the Sun, described in terms of solar irradiance, is important to climatology. A common assumption is that solar irradiance variability is driven by its surface magnetism. Verifying this assumption has, however, been hampered by the fact that models of solar irradiance variability based on solar surface magnetism have to be calibrated to observed variability. Making use of realistic three-dimensional magnetohydrodynamic simulations of the solar atmosphere and state-of-the-art solar magnetograms from the Solar Dynamics Observatory, we present a model of total solar irradiance (TSI) that does not require any such calibration. In doing so, the modeled irradiance variability is entirely independent of the observational record. (The absolute level is calibrated to the TSI record from the Total Irradiance Monitor.) The model replicates 95% of the observed variability between April 2010 and July 2016, leaving little scope for alternative drivers of solar irradiance variability at least over the time scales examined (days to years).

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...

We present a new empirical model of total and spectral solar irradiance (TSI and SSI) variability entitled EMPirical Irradiance REconstruction (EMPIRE). As with existing empirical models, TSI and SSI variability is given by the linear combination of solar activity indices. In empirical models, UV SSI variability is usually determined by fitting the rotational variability in activity indices to that in measurements. Such models have to date relied on ordinary least squares regression, which ignores the uncertainty in the activity indices. In an advance from earlier efforts, the uncertainty in the activity indices is accounted for in EMPIRE by the application of an error-in-variables regression scheme, making the resultant UV SSI variability more robust. The result is consistent with observations and unprecedentedly, with that from other modelling approaches, resolving the long-standing controversy between existing empirical models and other types of models. We demonstrate that earlier empirical models, by neglecting the uncertainty in activity indices, underestimate UV SSI variability. The reconstruction of TSI and visible and IR SSI from EMPIRE is also shown to be consistent with observations. The EMPIRE reconstruction is of utility to climate studies as a more robust alternative to earlier empirical reconstructions.

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.

We present a reconstruction of the spectral solar irradiance since 1700 using the SATIRE-T2 (Spectral And Total Irradiance REconstructions for the Telescope era version 2) model. This model uses as input magnetograms simulated with a surface flux transport model fed with semi-synthetic records of emerging sunspot groups. We used statistical relationships between the properties of sunspot group emergence, such as the latitude, area, and tilt angle, and the sunspot cycle strength and phase to produce semi-synthetic sunspot group records starting in the year 1700. The semisynthetic records are fed into a surface flux transport model to obtain daily simulated magnetograms that map the distribution of the magnetic flux in active regions (sunspots and faculae) and their decay products on the solar surface. The magnetic flux emerging in ephemeral regions is accounted for separately based on the concept of extended cycles whose length and amplitude are linked to those of the sunspot cycles through the sunspot number. The magnetic flux in each surface component (sunspots, faculae and network, and ephemeral regions) was used to compute the spectral and total solar irradiance between the years 1700 and 2009. This reconstruction is aimed at timescales of months or longer although the model returns daily values. We found that SATIRE-T2, besides reproducing other relevant observations such as the total magnetic flux, reconstructs the total solar irradiance (TSI) on timescales of months or longer in good agreement with the PMOD composite of observations, as well as with the reconstruction starting in 1878 based on the RGO-SOON data. The model predicts an increase in the TSI of 1.2[+0.2, -0.3] Wm-2 between 1700 and the present. The spectral irradiance reconstruction is in good agreement with the UARS/SUSIM measurements as well as the Lyman-alpha composite.

Regular spaceborne measurements have revealed that solar brightness varies on multiple timescales, variations on timescales greater than a day being attributed to surface magnetic field. Independently, ground-based and spaceborne measurements suggest that Sun-like stars show a similar, but significantly broader pattern of photometric variability. To understand whether the broader pattern of stellar variations is consistent with the solar paradigm we assess relative contributions of faculae and spots to solar magnetically-driven brightness variability. We investigate how the solar brightness variability as well as its facular and spot contributions depend on the wavelength, timescale of variability, and position of the observer relative to the ecliptic plane. We perform calculations with the SATIRE model, which returns solar brightness with daily cadence from solar disc area coverages of various magnetic features. Moving the observer away from the ecliptic plane increases the amplitude of 11-year variability as it would be seen in Str\"omgren $(b+y)/2$ photometry, but decreases the amplitude of the rotational brightness variations as it would appear in Kepler and CoRoT passbands. The spot and facular contributions to the 11-year solar variability in the Str\"omgren $(b+y)/2$ photometry almost fully compensate each other so that the Sun appears anomalously quiet with respect to its stellar cohort. Such a compensation does not occur on the rotational timescale. The rotational solar brightness variability as it would appear in Kepler and CoRoT passband from the ecliptic plane is spot-dominated but the relative contribution of faculae increases for out-of-ecliptic viewing so that the apparent brightness variations are faculae-dominated for inclinations less than about $i=45^{\circ}$.

Reliable historical records of total solar irradiance (TSI) are needed for climate change attribution and research to assess the extent to which long-term variations in the Sun's radiant energy incident on the Earth may exacerbate (or mitigate) the more dominant warming in recent centuries due to increasing concentrations of greenhouse gases. We investigate potential impacts of the new Sunspot Index and Long-term Solar Observations (SILSO) sunspot-number time series on model reconstructions of TSI. In contemporary TSI records, variations on time scales longer than about a day are dominated by the opposing effects of sunspot darkening and facular brightening. These two surface magnetic features, retrieved either from direct observations or from solar activity proxies, are combined in TSI models to reproduce the current TSI observational record. Indices that manifest solar-surface magnetic activity, in particular the sunspot-number record, then enable the reconstruction of historical TSI. Revisions to the sunspot-number record therefore affect the magnitude and temporal structure of TSI variability on centennial time scales according to the model reconstruction methodologies. We estimate the effects of the new SILSO record on two widely used TSI reconstructions, namely the NRLTSI2 and the SATIRE models. We find that the SILSO record has little effect on either model after 1885 but leads to greater amplitude solar-cycle fluctuations in TSI reconstructions prior, suggesting many 18th and 19th century cycles could be similar in amplitude to those of the current Modern Maximum. TSI records based on the revised sunspot data do not suggest a significant change in Maunder Minimum TSI values, and comparing that era to the present we find only very small potential differences in estimated solar contributions to climate with this new sunspot record.

As a part of the long-term program at Kitt Peak National Observatory (KPNO), the Mn I 539.4 nm line has been observed for nearly three solar cycles using the McMath telescope and the 13.5 m spectrograph in double-pass mode. These full-disk spectrophotometric observations revealed an unusually strong change of this line's parameters over the solar cycle. Optical pumping by the Mg II k line was originally proposed to explain these variations. More recent studies have proposed that this is not required and that the magnetic variability might explain it. Magnetic variability is also the mechanism that drives the changes in total solar irradiance variations (TSI). With this work we investigate this proposition quantitatively by using using the model SATIRE-S. We applied exactly the same model atmospheres and value of the free parameter as were used in previous solar irradiance reconstructions to now model the variation in the Mn I 539.4 nm line profile and in neighboring Fe I lines. We compared the results of the theoretical model with KPNO observations. Our result confirms that optical pumping of the Mn I 539.4 nm line by Mg II k is not the main cause of its solar cycle change. It also provides independent confirmation of solar irradiance models which are based on the assumption that irradiance variations are caused by the evolution of the solar surface magnetic flux. The result obtained here also supports the spectral irradiance variations computed by these models.

Extending the knowledge about the properties of solar cycles into the past is essential for understanding the solar dynamo. This paper aims at estimating areas of sunspots observed by Schwabe in 1825-1867 and at calculating the tilt angles of sunspot groups. The sunspot sizes in Schwabe's drawings are not to scale and need to be converted into physical sunspot areas. We employed a statistical approach assuming that the area distribution of sunspots was the same in the 19th century as it was in the 20th century. Umbral areas for about 130,000 sunspots observed by Schwabe were obtained, as well as the tilt angles of sunspot groups assuming them to be bipolar. There is, of course, no polarity information in the observations. The annually averaged sunspot areas correlate reasonably with sunspot number. We derived an average tilt angle by attempting to exclude unipolar groups with a minimum separation of the two alleged polarities and an outlier rejection method which follows the evolution of each group and detects the moment it turns unipolar at its decay. As a result, the tilt angles, although displaying considerable scatter, place the leading polarity on average 5.85+-0.25 closer to the equator, in good agreement with tilt angles obtained from 20th-century data sets. Sources of uncertainties in the tilt angle determination are discussed and need to be addressed whenever different data sets are combined. The sunspot area and tilt angle data are provided online.

The solar brightness varies on timescales from minutes to decades. A clear identification of the physical processes behind such variations is needed for developing and improving physics-based models of solar brightness variability and reconstructing solar brightness in the past. This is, in turn, important for better understanding the solar-terrestrial and solar-stellar connections. We estimate the relative contributions of the continuum, molecular, and atomic lines to the solar brightness variations on different timescales. Our approach is based on the assumption that variability of the solar brightness on timescales greater than a day is driven by the evolution of the solar surface magnetic field. We calculated the solar brightness variations employing the solar disc area coverage of magnetic features deduced from the MDI/SOHO observations. The brightness contrasts of magnetic features relative to the quiet Sun were calculated with a non-LTE radiative transfer code as functions of disc position and wavelength. By consecutive elimination of molecular and atomic lines from the radiative transfer calculations, we assessed the role of these lines in producing solar brightness variability. We show that the variations in Fraunhofer lines define the amplitude of the solar brightness variability on timescales greater than a day and even the phase of the total solar irradiance variability over the 11-year cycle. We also demonstrate that molecular lines make substantial contribution to solar brightness variability on the 11-year activity cycle and centennial timescales. In particular, our model indicates that roughly a quarter of the total solar irradiance variability over the 11-year cycle originates in molecular lines. The maximum of the absolute spectral brightness variability on timescales greater than a day is associated with the CN violet system between 380 and 390 nm.

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.

Total solar irradiance and UV spectral solar irradiance have been monitored since 1978 through a succession of space missions. This is accompanied by the development of models aimed at replicating solar irradiance by relating the variability to solar magnetic activity. The NRLSSI and SATIRE-S models provide the most comprehensive reconstructions of total and spectral solar irradiance over the period of satellite observation currently available. There is persistent controversy between the various measurements and models in terms of the wavelength dependence of the variation over the solar cycle, with repercussions on our understanding of the influence of UV solar irradiance variability on the stratosphere. We review the measurement and modelling of UV solar irradiance variability over the period of satellite observation. The SATIRE-S reconstruction is consistent with spectral solar irradiance observations where they are reliable. It is also supported by an independent, empirical reconstruction of UV spectral solar irradiance based on UARS/SUSIM measurements from an earlier study. The weaker solar cycle variability produced by NRLSSI between 300 and 400 nm is not evident in any available record. We show that although the method employed to construct NRLSSI is principally sound, reconstructed solar cycle variability is detrimentally affected by the uncertainty in the SSI observations it draws upon in the derivation. Based on our findings, we recommend, when choosing between the two models, the use of SATIRE-S for climate studies.

Total and spectral solar irradiance are key parameters in the assessment of solar influence on changes in the Earth's climate. We present a reconstruction of daily solar irradiance obtained using the SATIRE-S model spanning 1974 to 2013 based on full-disc observations from the KPVT, SoHO/MDI and SDO/HMI. SATIRE-S ascribes variation in solar irradiance on timescales greater than a day to photospheric magnetism. The solar spectrum is reconstructed from the apparent surface coverage of bright magnetic features and sunspots in the daily data using the modelled intensity spectra of these magnetic structures. We cross-calibrated the various data sets, harmonizing the model input so as to yield a single consistent time series as the output. The model replicates 92% of the variability in the PMOD TSI composite including the secular decline between the 1996 and 2008 solar cycle minima. The model also reproduces most of the variability in observed Lyman-alpha irradiance and the Mg II index. The UV solar irradiance measurements from the UARS and SORCE missions are mutually consistent up to about 180 nm before they start to exhibit discrepant rotational and cyclical variability, indicative of unresolved instrumental effects. As a result, the agreement between model and measurement starts to deteriorate above this wavelength. As with earlier similar investigations, the reconstruction cannot reproduce the overall trends in SORCE/SIM SSI. We argue, from the lack of clear solar cycle modulation in the SIM record and the inconsistency between the total flux recorded by the instrument and TSI, that unaccounted instrumental trends are present. The reconstruction is consistent with observations from multiple sources, demonstrating its validity. It also provides further evidence that photospheric magnetism is the prime driver of variation in solar irradiance on timescales greater than a day.

We present a PSF for SDO/HMI and discuss the effects of its removal on the apparent properties of solar surface phenomena in HMI data. The PSF was retrieved from observations of Venus in transit by matching it to the convolution of a model of the venusian disc and solar background with a guess PSF. Observations recorded near in time to the transit of Venus were corrected for instrumental scattered light by the deconvolution with the PSF. Granulation contrast in restored HMI data is greatly enhanced relative to the original data and exhibit reasonable agreement with numerical simulations. Image restoration enhanced the apparent intensity and pixel averaged magnetic field strength of photospheric magnetic features significantly. For small-scale magnetic features, restoration enhanced intensity contrast in the continuum and core of the Fe I 617.3 nm line by a factor of 1.3, and the magnetogram signal by a factor of 1.7. For sunspots and pores, the enhancement varied strongly within and between features, being more acute for smaller features. Magnetic features are also rendered smaller, as signal smeared onto the surrounding quiet Sun is recovered. Image restoration increased the apparent amount of magnetic flux above the noise floor by a factor of about 1.2, most of the gain coming from the quiet Sun. Line-of-sight velocity due to granulation and supergranulation is enhanced by a factor of 1.4 to 2.1, depending on position on the solar disc. There are indications that the PSF varies across the FOV, with time and between the two CCDs. However, all these variations were found to be relatively small, such that the derived PSF can be applied to HMI data from both CCDs, over the period examined without introducing significant error. Image restoration affects the measured radiant, magnetic and dynamic properties of solar surface phenomena sufficiently to significantly impact interpretation.

We present a new model of total solar irradiance (TSI) based on magnetograms simulated with a surface flux transport model (SFTM) and the SATIRE (Spectral And Total Irradiance REconstructions) model. Our model provides daily maps of the distribution of the photospheric field and the TSI starting from 1878. We first calculate the magnetic flux on the solar surface emerging in active and ephemeral regions. The evolution of the magnetic flux in active regions is computed using a surface flux transport model fed with the observed record of sunspot group areas and positions. The magnetic flux in ephemeral regions is treated separately using the concept of overlapping cycles. To model the ephemeral region cycles, we assume that their length and amplitude are related to that of the sunspot cycles. We then use a version of the SATIRE model to compute the TSI. The area coverage and the distribution of different magnetic features as a function of time, which are required by SATIRE, are extracted from the simulated magnetograms and the modelled ephemeral region magnetic flux. Previously computed intensity spectra of the various types of magnetic features are employed. Our model reproduces the PMOD composite of TSI measurements starting from 1978 at daily and rotational timescales more accurately than the previous version of the SATIRE model computing TSI over this period of time. The simulated magnetograms provide a more realistic representation of the evolution of the magnetic field on the photosphere and also allow us to make use of information on the spatial distribution of the magnetic fields before the times when observed magnetograms were available. We find that the secular increase in TSI since 1878 is fairly stable to modifications of the treatment of the ephemeral region magnetic flux.

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.

The correlation between solar irradiance and the 11-year solar activity cycle is evident in the body of measurements made from space, which extend over the past four decades. Models relating variation in solar irradiance to photospheric magnetism have made significant progress in explaining most of the apparent trends in these observations. There are, however, persistent discrepancies between different measurements and models in terms of the absolute radiometry, secular variation and the spectral dependence of the solar cycle variability. We present an overview of solar irradiance measurements and models, and discuss the key challenges in reconciling the divergence between the two.

The Sun and stars with low magnetic activity levels, become photometrically brighter when their activity increases. Magnetically more active stars display the opposite behaviour and get fainter when their activity increases. We reproduce the observed photometric trends in stellar variations with a model that treats stars as hypothetical Suns with coverage by magnetic features different from that of the Sun. The presented model attributes the variability of stellar spectra to the imbalance between the contributions from different components of the solar atmosphere, such as dark starspots and bright faculae. A stellar spectrum is calculated from spectra of the individual components, by weighting them with corresponding disc area coverages. The latter are obtained by extrapolating the solar dependences of spot and facular disc area coverages on chromospheric activity to stars with different levels of mean chromospheric activity. We have found that the contribution by starspots to the variability increases faster with chromospheric activity than the facular contribution. This causes the transition from faculae-dominated variability and direct activity--brightness correlation to spot-dominated variability and inverse activity--brightness correlation with increasing chromospheric activity level. We have shown that the regime of the variability also depends on the angle between the stellar rotation axis and the line-of-sight and on the latitudinal distribution of active regions on the stellar surface. Our model can be used as a tool to extrapolate the observed photometric variability of the Sun to Sun-like stars at different activity levels, which makes possible the direct comparison between solar and stellar irradiance data.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

The brightness of the Sun varies on all time scales on which it has been observed, and there is increasing evidence that it has an influence on climate. The amplitudes of such variations depend on the wavelength and possibly on the time scale. Although many aspects of this variability are well established, the exact magnitude of secular variations (going beyond a solar cycle) and the spectral dependence of variations are under discussion. The main drivers of solar variability are thought to be magnetic features at the solar surface. The climate reponse can be, on a global scale, largely accounted for by simple energetic considerations, but understanding the regional climate effects is more difficult. Promising mechanisms for such a driving have been identified, including through the influence of UV irradiance on the stratosphere and dynamical coupling to the surface. Here we provide an overview of the current state of our knowledge, as well as of the main open questions.

The lack of long and reliable time series of solar spectral irradiance (SSI) measurements makes an accurate quantification of solar contributions to recent climate change difficult. Whereas earlier SSI observations and models provided a qualitatively consistent picture of the SSI variability, recent measurements by the SORCE satellite suggest a significantly stronger variability in the ultraviolet (UV) spectral range and changes in the visible and near-infrared (NIR) bands in anti-phase with the solar cycle. A number of recent chemistry-climate model (CCM) simulations have shown that this might have significant implications on the Earth's atmosphere. Motivated by these results, we summarize here our current knowledge of SSI variability and its impact on Earth's climate. We present a detailed overview of existing SSI measurements and provide thorough comparison of models available to date. SSI changes influence the Earth's atmosphere, both directly, through changes in shortwave (SW) heating and therefore, temperature and ozone distributions in the stratosphere, and indirectly, through dynamical feedbacks. We investigate these direct and indirect effects using several state-of-the art CCM simulations forced with measured and modeled SSI changes. A unique asset of this study is the use of a common comprehensive approach for an issue that is usually addressed separately by different communities. Omissis. Finally, we discuss the reliability of the available data and we propose additional coordinated work, first to build composite SSI datasets out of scattered observations and to refine current SSI models, and second, to run coordinated CCM experiments.

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.

This study aims at setting observational constraints on the continuum and line core intensity contrast of network and faculae, specifically, their relationship with magnetic field and disc position. Full-disc magnetograms and intensity images by the HMI instrument onboard SDO were employed. Bright magnetic features, representing network and faculae, were identified and the relationship between their intensity contrast at continuum and line core with magnetogram signal and heliocentric angle examined. Care was taken to minimize the inclusion of the magnetic canopy and straylight from sunspots and pores as network and faculae. In line with earlier studies, network features, on a per unit magnetic flux basis, appeared brighter than facular features. Intensity contrasts in the continuum and line core differ considerably, most notably, they exhibit opposite centre-to-limb variations. We found this difference in behaviour to likely be due to the different mechanisms of the formation of the two spectral components. From a simple model based on bivariate polynomial fits to the measured contrasts we confirmed spectral line changes to be a significant driver of facular contribution to variation in solar irradiance. The discrepancy between the continuum contrast reported here and in the literature was shown to arise mainly from differences in spatial resolution and treatment of magnetic signals adjacent to sunspots and pores. HMI is a source of accurate contrasts and low-noise magnetograms covering the full solar disc. For irradiance studies it is important to consider not just the contribution from the continuum but also from the spectral lines. In order not to underestimate long-term variations in solar irradiance, irradiance models should take the greater contrast per unit magnetic flux associated with magnetic features with low magnetic flux into account.

Context: The study of variations in total solar irradiance (TSI) is important for understanding how the Sun affects the Earth's climate. Aims: Full-disk continuum images and magnetograms are now available for three full solar cycles. We investigate how modelled TSI compares with direct observations by building a consistent modelled TSI dataset. The model, based only on changes in the photospheric magnetic flux can then be tested on rotational, cyclical and secular timescales. Methods: We use Kitt Peak and SoHO/MDI continuum images and magnetograms in the SATIRE-S model to reconstruct TSI over cycles 21-23. To maximise independence from TSI composites, SORCE/TIM TSI data are used to fix the one free parameter of the model. We compare and combine the separate data sources for the model to estimate an uncertainty on the reconstruction and prevent any additional free parameters entering the model. Results: The reconstruction supports the PMOD composite as being the best historical record of TSI observations, although on timescales of the solar rotation the IRMB composite provides somewhat better agreement. Further to this, the model is able to account for 92% of TSI variations from 1978 to 2009 in the PMOD composite and over 96% during cycle 23. The reconstruction also displays an inter-cycle, secular decline of 0.20 (+0.12 / -0.09) Wm-2 between cycle 23 minima, in agreement with the PMOD composite. Conclusions: SATIRE-S is able to recreate TSI observations on all timescales of a day and longer over 31 years from 1978. This is strong evidence that changes in photospheric magnetic flux alone are responsible for almost all solar irradiance variations over the last three solar cycles.

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.

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.

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.

The apparently unusual behaviour of the TSI during the most recent minimum of solar activity has been interpreted as evidence against solar surface magnetism as the main driver of the secular change in the TSI. We test claims that the evolution of the solar surface magnetic field does not reproduce the observed TSI in cycle 23. We use sensitive, 60-minute averaged MDI magnetograms and quasi-simultaneous continuum images as an input to our SATIRE-S model and calculate the TSI variation over cycle 23, sampled roughly twice-monthly. The computed TSI is then compared to the PMOD composite of TSI measurements and to the data from two individual instruments, SORCE/TIM and UARS/ACRIM II, that monitored the TSI during the declining phase of cycle 23 and over the previous minimum in 1996, respectively. Excellent agreement is found between the trends shown by the model and almost all sets of measurements. The only exception is the early, i.e. 1996 to 1998, PMOD data. Whereas the agreement between the model and the PMOD composite over the period 1999-2009 is almost perfect, the modelled TSI shows a steeper increase between 1996 and 1999 than implied by the PMOD composite. On the other hand, the steeper trend in the model agrees remarkably well with the ACRIM II data. A closer look at the VIRGO data, that make the basis of the PMOD composite after 1996, reveals that only one of the two VIRGO instruments, the PMO6V, shows the shallower trend present in the composite, whereas the DIARAD measurements indicate a steeper trend. We conclude that (1) the sensitivity changes of the PMO6V radiometers within VIRGO during the first two years have very likely not been correctly evaluated, and that (2) the TSI variations over cycle 23 and the change in the TSI levels between the minima in 1996 and 2008 are consistent with the solar surface magnetism mechanism.

The historical record of sunspot areas is a valuable and widely used proxy of solar activity and variability. The Royal Greenwich Observatory (RGO) regularly measured this and other parameters between 1874 and 1976. After that time records from a number of different observatories are available. These, however, show systematic differences and often have significants gaps. Our goal is to obtain a uniform and complete sunspot area time series by combining different data sets. A homogeneus composite of sunspot areas is essential for different applications in solar physics, among others for irradiance reconstructions. Data recorded simultaneously at different observatories are statistically compared in order to determine the intercalibration factors. Using these data we compile a complete and cross-calibrated time series. The Greenwich data set is used as a basis until 1976, the Russian data (a compilation of observations made at stations in the former USSR) between 1977 and 1985 and data compiled by the USAF network since 1986. Other data sets (Rome, Yunnan, Catania) are used to fill up the remaining gaps. Using the final sunspot areas record the Photometric Sunspot Index is calculated. We also show that the use of uncalibrated sunspot areas data sets can seriously affect the estimate of irradiance variations. Our analysis implies that there is no basis for the claim that UV irradiance variations have a much smaller influence on climate than total solar irradiance variations.

It is known that the tilt angles of active regions increase with their latitude (Joy's law). It has never been checked before, however, whether the average tilt angles change from one cycle to another. Flux transport models show the importance of tilt angles for the reversal and build up of magnetic flux at the poles which is, in turn, correlated with the strength of the next cycle. Here we analyse time series of tilt angle measurements and look for a possible relationship of the tilt angles with other solar cycle parameters, in order to glean information on the solar dynamo and to estimate their potential for predictions of solar activity. We employ tilt angle data from Mount Wilson and Kodaikanal observatories covering solar cycles 15 to 21. We analyse the latitudinal distribution of the tilt angles (Joy's law), their variation from cycle to cycle and their relationship to other solar cycle parameters, such as the strength, amplitude and length. The two main results are: 1. An anti-correlation between the mean normalized tilt angle of a given cycle and the strength (or amplitude) of that cycle, with a correlation coefficient of r=-0.95 and r=-0.93 for Mount Wilson and Kodaikanal data, respectively. 2. The product of the cycle averaged tilt angle and the strength of the same cycle displays a significant correlation with the strength of the next cycle (r=0.65 and r=0.70 for Mount Wilson and Kodaikanal data, respectively). An even better correlation is obtained between the source term of the poloidal flux in Babcock-Leighton-type dynamos (which contains the tilt angle) and the amplitude of the next cycle. The results of this study indicate that in combination with the cycle strength, the active region tilt angles play an important role in building up the polar fields at cycle minimum.

The variation of total solar irradiance (TSI) has been measured since 1978 and that of the spectral irradiance for an even shorter amount of time. Semi-empirical models are now available that reproduce over 80% of the measured irradiance variations. An extension of these models into the more distant past is needed in order to serve as input to climate simulations. Here we review our most recent efforts to model solar total and spectral irradiance on time scales from days to centuries and even longer. Solar spectral irradiance has been reconstructed since 1947. Reconstruction of solar total irradiance goes back to 1610 and suggests a value of about 1-1.5 Wm$^{-2}$ for the increase in the cycle-averaged TSI since the end of the Maunder minimum, which is significantly lower than previously assumed but agrees with other modern models. First steps have also been made towards reconstructions of solar total and spectral irradiance on time scales of millennia.

A gap in the total solar irradiance (TSI) measurements between ACRIM-1 and ACRIM-2 led to the ongoing debate on the presence or not of a secular trend between the minima preceding cycles 22 (in 1986) and 23 (1996). It was recently proposed to use the SATIRE model of solar irradiance variations to bridge this gap. When doing this, it is important to use the appropriate SATIRE-based reconstruction, which we do here, employing a reconstruction based on magnetograms. The accuracy of this model on months to years timescales is significantly higher than that of a model developed for long-term reconstructions used by the ACRIM team for such an analysis. The constructed `mixed' ACRIM - SATIRE composite shows no increase in the TSI from 1986 to 1996, in contrast to the ACRIM TSI composite.

Total solar irradiance reconstructed between 1978 and 2003 using solar surface magnetic field distributions is compared with three composites of total solar irradiance measurements. A good correspondence is found with the total solar irradiance composite from PMOD/WRC, with no bias between the three cycles. The agreement with the other composites (the ACRIM composite, mainly based on the Active Cavity Radiometer Irradiance Monitors I, II & III, and the IRMB composite from the Institut Royal Meteorologique Belgique) is significantly poorer. In particular, a secular increase in the irradiance exhibited by these composites is not present in the reconstructions. Hence any secular trend in total solar irradiance between 1978 and 2003 is not due to magnetic fields at the solar surface.

Variations of the solar UV irradiance are an important driver of chemical and physical processes in the Earth's upper atmosphere and may also influence global climate. Here we reconstruct solar UV irradiance in the range 115-400 nm over the period 1974-2007 by making use of the recently developed empirical extension of the SATIRE models employing SUSIM data. The evolution of the solar photospheric magnetic flux, which is a central input to the model, is described by the magnetograms and continuum images recorded at the Kitt Peak National Solar Observatory between 1974 and 2003 and by the MDI instrument on SoHO since 1996. The reconstruction extends the available observational record by 1.5 solar cycles. The reconstructed Ly-alpha irradiance agrees well with the composite time series by Woods et al (2000). The amplitude of the irradiance variations grows with decreasing wavelength and in the wavelength regions of special interest for studies of the Earth's climate (Ly-alpha and oxygen absorption continuum and bands between 130 and 350 nm) is one to two orders of magnitude stronger than in the visible or if integrated over all wavelengths (total solar irradiance).

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.

Various observatories around the globe started regular full-disk imaging of the solar atmosphere in the Ca II K line since the early decades of the 20th century. The archives made by these observations have the potential of providing far more detailed information on solar magnetism than just the sunspot number and area records to which most studies of solar activity and irradiance changes are restricted. We evaluate the image contents of three Ca II K spectroheliogram time-series, specifically those obtained by the digitization of the Arcetri, Kodaikanal, and Mt Wilson photographic archives. We describe the main problems afflicting these data and analyze their quality by expressing the image contents through several quantities. We compare the results obtained with those for similar present-day observations taken with the Meudon spectroheliograph and with the Rome-PSPT. We show that historic data suffer from stronger geometrical distortions and photometric uncertainties than similar present-day observations. The latter uncertainties mostly originate from the photographic calibration of the original data and from stray-light effects. We also show that the image contents of the three analyzed series vary in time. These variations are probably due to instrument changes and aging of the spectrographs used, as well as changes of the observing programs. Our results imply that the main challenge for the analysis of historic data is their accurate photometric calibration. This problem must be solved before they can provide reliable information about solar magnetism and activity over the last century. Moreover, inter-calibration of results obtained from independent time-series is required to reliably trace changes of solar properties with time from the analysis of such data.

We present a reconstruction of total solar irradiance since 1610 to the present based on variations of the surface distribution of the solar magnetic field. The latter is calculated from the historical record of the Group sunspot number using a simple but consistent physical model. Our model successfully reproduces three independent data sets: total solar irradiance measurements available since 1978, total photospheric magnetic flux from 1974 and the open magnetic flux since 1868 (as empirically reconstructed from the geomagnetic aa-index). The model predicts an increase in the total solar irradiance since the Maunder Minimum of about 1.3 \rm{Wm$^{-2}$}.

In this paper we investigate the dependence of umbral core brightness, as well as the mean umbral and penumbral brightness on the phase of the solar cycle and on the size of the sunspot. Albregtsen & Maltby (1978) reported an increase in umbral core brightness from the early to the late phase of solar cycle from the analysis of 13 sunspots which cover solar cycles 20 and 21. Here we revisit this topic by analysing continuum images of more than 160 sunspots observed by the MDI instrument on board the SOHO spacecraft for the period between 1998 March to 2004 March, i.e. a sizable part of solar cycle 23. The advantage of this data set is its homogeneity, with no seeing fluctuations. A careful stray light correction, which is validated using the Mercury transit of 7th May, 2003, is carried out before the umbral and penumbral intensities are determined. The influence of the Zeeman splitting of the nearby NiI spectral line on the measured 'continuum' intensity is also taken into account. We did not observe any significant variation in umbral core, mean umbral and mean penumbral intensities with solar cycle, which is in contrast to earlier findings for the umbral core intensity. We do find a strong and clear dependence of the umbral brightness on sunspot size, however. The penumbral brightness also displays a weak dependence. The brightness-radius relationship has numerous implications, some of which, such as those for the energy transport in umbrae, are pointed out.