Showing 1–50 of 82 results for author: Krivova, N A

Butterfly Diagram and Other Properties of Plage Areas from Kodaikanal Ca II K Photographs Covering 1904-2007

Authors: Bibhuti Kumar Jha, Theodosios Chatzistergos, Dipankar Banerjee, Ilaria Ermolli, Natalie A. Krivova, Sami K. Solanki, Aditya Priyadarshi

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

Submitted 23 August, 2024; originally announced August 2024.

Comments: 15 pages, 9 figures, Under Review in Solar Physics

The variation in the response of solar full-disc magnetographs

Authors: K. L. Yeo, S. K. Solanki, N. A. Krivova

Abstract: We compared magnetograms from the KPVT/SPM, SoHO/MDI, SOLIS/VSM, and SDO/HMI with the aim of probing the effect on measured solar magnetism of the variation in instrument response with time, magnetogram signal level, and position on the solar disc. Taking near-simultaneous observations from the various instruments, we examined the surface coverage by magnetic activity and the effect of cross-calib… ▽ More We compared magnetograms from the KPVT/SPM, SoHO/MDI, SOLIS/VSM, and SDO/HMI with the aim of probing the effect on measured solar magnetism of the variation in instrument response with time, magnetogram signal level, and position on the solar disc. Taking near-simultaneous observations from the various instruments, we examined the surface coverage by magnetic activity and the effect of cross-calibrating the various instruments under different assumptions. By comparing the surface coverage by magnetic activity in the observations from the various instruments, we traced the effect of the time variation in instrument response on the longitudinal magnetogram signal and disc-integrated unsigned magnetic flux. This yielded evidence of acute changes in the response of MDI and VSM with certain events such as the SoHO vacation in 1998 and the upgrade of the VSM CCD camera in 2009. Excluding these changes, the effect of instrument instability on the magnetogram signal and disc-integrated magnetic flux appears to be rather benign, with an associated uncertainty of less than 2%. We determined the magnetogram signal ratio between each instrument pairing as a function of magnetogram signal level and distance from disc centre and with it cross-calibrated the various instruments. We compared the result with that from repeating the cross-calibration with the overall magnetogram signal ratio. This allowed us to estimate the uncertainty in the magnetogram signal associated with the variation in instrument response with magnetogram signal level and distance from disc centre to be about 8% to 14%. The corresponding uncertainty in the disc-integrated magnetic flux is about 7% to 23 %. The results here will be useful to the interpretation of SPM, MDI, VSM, and HMI magnetograms. As examples, we applied our findings to selected results from earlier studies based on such data. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Stellar Surface Magnetic Fields Impact Limb Darkening

Authors: Nadiia M. Kostogryz, Alexander I. Shapiro, Veronika Witzke, Robert H. Cameron, Laurent Gizon, Natalie A. Krivova, Hans-G. Ludwig, Pierre F. L. Maxted, Sara Seager, Sami K. Solanki, Jeff Valenti

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

Submitted 29 February, 2024; originally announced March 2024.

Comments: 19 pages, 6 figures This is the submission version of manuscript to Nature Astronomy

Reconstruction of total solar irradiance variability as simultaneously apparent from Solar Orbiter and Solar Dynamics Observatory

Authors: K. L. Yeo, N. A. Krivova, S. K. Solanki, J. Hirzberger, D. Orozco Suárez, K. Albert, N. Albelo Jorge, T. Appourchaux, A. Alvarez-Herrero, J. Blanco Rodríguez, A. Gandorfer, P. Gutierrez-Marques, F. Kahil, M. Kolleck, J. C. del Toro Iniesta, R. Volkmer, J. Woch, B. Fiethe, I. Pérez-Grande, E. Sanchis Kilders, M. Balaguer Jiménez, L. R. Bellot Rubio, D. Calchetti, M. Carmona, A. Feller , et al. (20 additional authors not shown)

Abstract: Solar irradiance variability has been monitored almost exclusively from the Earth's perspective. {We present a method to combine the unprecedented observations of the photospheric magnetic field and continuum intensity from outside the Sun-Earth line, which is being recorded by the Polarimetric and Helioseismic Imager on board the Solar Orbiter mission (SO/PHI), with solar observations recorded fr… ▽ More Solar irradiance variability has been monitored almost exclusively from the Earth's perspective. {We present a method to combine the unprecedented observations of the photospheric magnetic field and continuum intensity from outside the Sun-Earth line, which is being recorded by the Polarimetric and Helioseismic Imager on board the Solar Orbiter mission (SO/PHI), with solar observations recorded from the Earth's perspective to examine the solar irradiance variability from both perspectives simultaneously.} Taking SO/PHI magnetograms and continuum intensity images from the cruise phase of the Solar Orbiter mission and concurrent observations from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) as input into the SATIRE-S model, we successfully reconstructed the total solar irradiance variability as apparent from both perspectives. In later stages of the SO mission, the orbital plane will tilt in such a way as to bring the spacecraft away from the ecliptic to heliographic latitudes of up to $33^{\circ}$. The current study sets the template for the reconstruction of solar irradiance variability as seen from outside the ecliptic from data that SO/PHI is expected to collect from such positions. {Such a reconstruction will be beneficial to factoring inclination into how the brightness variations of the Sun compare to those of other cool stars, whose rotation axes are randomly inclined. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Journal ref: A&A 679, A25 (2023)

Analysis of full-disc H alpha observations: Carrington maps and filament properties in 1909-2022

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Dipankar Banerjee, Teresa Barata, Ioannis Chouinavas, Mariachiara Falco, Ricardo Gafeira, Fabrizio Giorgi, Yoichiro Hanaoka, Natalie A. Krivova, Viktor V. Korokhin, Ana Lourenço, Gennady P. Marchenko, Jean-Marie Malherbe, Nuno Peixinho, Paolo Romano, Takashi Sakurai

Abstract: Full disc observations of the Sun in the H$αあるふぁ$ line provide information about the solar chromosphere and in particular about the filaments, which are dark and elongated features that lie along magnetic field polarity inversion lines. This makes them important for studies of solar magnetism. Since full disc H$αあるふぁ$ observations have been performed at various sites since 1800s, with regular photographic… ▽ More Full disc observations of the Sun in the H$αあるふぁ$ line provide information about the solar chromosphere and in particular about the filaments, which are dark and elongated features that lie along magnetic field polarity inversion lines. This makes them important for studies of solar magnetism. Since full disc H$αあるふぁ$ observations have been performed at various sites since 1800s, with regular photographic data having started in the beginning of the 20th century, they are an invaluable source of information on past solar magnetism. In this work we aimed at deriving accurate information about filaments from historical and modern full disc H$αあるふぁ$ observations. We have consistently processed observations from 15 H$αあるふぁ$ archives spanning 1909-2022. Our data processing includes photometric calibration of the data stored on photographic plates. We have constructed also Carrington maps from the calibrated H$αあるふぁ$ images. We find that filament areas are affected by the bandwidth of the observation. Thus, cross-calibration of the filament areas derived from different archives is needed. We have produced a composite of filament areas from individual archives by scaling all of them to the Meudon series. Our composite butterfly diagram shows very distinctly the common features of filament evolution, that is the poleward migration as well as a decrease in the mean latitude of filaments as the cycle progresses. We also find that during activity maxima, filaments on average cover about 1% of the solar surface. We see only a weak change in the amplitude of cycles in filament areas, in contrast to sunspot and plage areas. Analysis of H$αあるふぁ$ data for archives with contemporaneous Ca II K observations allowed us to identify and verify archive inconsistencies, which will also have implications for reconstructions of past solar magnetism and irradiance from Ca II K data. △ Less

Submitted 18 October, 2023; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 21 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 680, A15 (2023)

Modeling Stellar CaII H & K Emission Variations: Spot Contribution to the S-index

Authors: K. Sowmya, A. I. Shapiro, L. H. M. Rouppe van der Voort, N. A. Krivova, S. K. Solanki

Abstract: The S-index is a measure of emission in the CaII H & K lines and is a widely used proxy of stellar magnetic activity. It has been assumed until now that the S-index is mainly affected by bright plage regions in the chromosphere. In particular, the effect of starspots on the S-index has been neglected. In this study we revisit this assumption. For this we analyze high-resolution observations of sun… ▽ More The S-index is a measure of emission in the CaII H & K lines and is a widely used proxy of stellar magnetic activity. It has been assumed until now that the S-index is mainly affected by bright plage regions in the chromosphere. In particular, the effect of starspots on the S-index has been neglected. In this study we revisit this assumption. For this we analyze high-resolution observations of sunspots recorded in the CaII H spectral line at the Swedish 1-m Solar Telescope and determine the contrast of spots with respect to the quiet surroundings. We find that the CaII H line core averaged over whole sunspots (including superpenumbrae) is brighter than in the quiet surroundings and that the spot contrast in the line core is comparable to the facular contrast. This allows us to get a first estimate of the influence of spots on the S-index. We show that spots increase the S-index. While this increase is quite small for the Sun, it becomes significantly larger for more active stars. Further, we show that the inclusion of the contribution of spots to the S-index strongly affects the relationship between the S-index and stellar disk area coverages by spots and faculae, and present the new relations. △ Less

Submitted 23 September, 2023; v1 submitted 7 September, 2023; originally announced September 2023.

Comments: 11 pages, 5 figures, accepted for publication in ApJ Letters

Intensity contrast of solar network and faculae close to the solar limb, observed from two vantage points

Authors: K. Albert, N. A. Krivova, J. Hirzberger, S. K. Solanki, A. Moreno Vacas, D. Orozco Suárez, N. Albelo Jorge, T. Appourchaux, A. Alvarez-Herrero, J. Blanco Rodríguez, A. Gandorfer, P. Gutierrez-Marques, F. Kahil, M. Kolleck, R. Volkmer, J. C. del Toro Iniesta, J. Woch, B. Fiethe, I. Pérez-Grande, E. Sanchis Kilders, M. Balaguer Jiménez, L. R. Bellot Rubio, D. Calchetti, M. Carmona, A. Feller , et al. (21 additional authors not shown)

Abstract: The brightness of faculae and network depends on the angle at which they are observed and the magnetic flux density. Close to the limb, assessment of this relationship has until now been hindered by the increasingly lower signal in magnetograms. This preliminary study aims at highlighting the potential of using simultaneous observations from different vantage points to better determine the propert… ▽ More The brightness of faculae and network depends on the angle at which they are observed and the magnetic flux density. Close to the limb, assessment of this relationship has until now been hindered by the increasingly lower signal in magnetograms. This preliminary study aims at highlighting the potential of using simultaneous observations from different vantage points to better determine the properties of faculae close to the limb. We use data from the Solar Orbiter/Polarimetric and Helioseismic Imager (SO/PHI), and the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI), recorded at $\sim60^\circ$ angular separation of their lines of sight at the Sun. We use continuum intensity observed close to the limb by SO/PHI and complement it with the co-observed $B_{\rm LOS}$ from SDO/HMI, originating closer to disc centre (as seen by SDO/HMI), thus avoiding the degradation of the magnetic field signal near the limb. We derived the dependence of facular brightness in the continuum on disc position and magnetic flux density from the combined observations of SO/PHI and SDO/HMI. Compared with a single point of view, we were able to obtain contrast values reaching closer to the limb and to lower field strengths. We find the general dependence of the limb distance at which the contrast is maximum on the flux density to be at large in line with single viewpoint observations, in that the higher the flux density is, the closer the turning point lies to the limb. There is a tendency, however, for the maximum to be reached closer to the limb when determined from two vantage points. We note that due to the preliminary nature of this study, these results must be taken with caution. Our analysis shows that studies involving two viewpoints can significantly improve the detection of faculae near the solar limb and the determination of their brightness contrast relative to the quiet Sun. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Journal ref: A&A 678, A163 (2023)

Spectral variability of photospheric radiation due to faculae II: Facular contrasts for cool main-sequence stars

Authors: Charlotte M. Norris, Yvonne C. Unruh, Veronika Witzke, Sami K. Solanki, Natalie A. Krivova, Alexander I. Shapiro, Robert Cameron, Benjamin Beeck

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

Submitted 7 June, 2023; originally announced June 2023.

Comments: 17 pages, 12 figures, to be published in MNRAS

Long-term changes in solar activity and irradiance

Authors: Theodosios Chatzistergos, Natalie A. Krivova, Kok Leng Yeo

Abstract: The Sun is the main energy source to Earth, and understanding its variability is of direct relevance to climate studies. Measurements of total solar irradiance exist since 1978, but this is too short compared to climate-relevant time scales. Coming from a number of different instruments, these measurements require a cross-calibration, which is not straightforward, and thus several composite record… ▽ More The Sun is the main energy source to Earth, and understanding its variability is of direct relevance to climate studies. Measurements of total solar irradiance exist since 1978, but this is too short compared to climate-relevant time scales. Coming from a number of different instruments, these measurements require a cross-calibration, which is not straightforward, and thus several composite records have been created. All of them suggest a marginally decreasing trend since 1996. Most composites also feature a weak decrease over the entire period of observations, which is also seen in observations of the solar surface magnetic field and is further supported by Ca II K data. Some inconsistencies, however, remain and overall the magnitude and even the presence of the long-term trend remain uncertain. Different models have been developed, which are used to understand the irradiance variability over the satellite period and to extend the records of solar irradiance back in time. Differing in their methodologies, all models require proxies of solar magnetic activity as input. The most widely used proxies are sunspot records and cosmogenic isotope data on centennial and millennial time scale, respectively. None of this, however, offers a sufficiently good, independent description of the long-term evolution of faculae and network responsible for solar brightening. This leads to uncertainty in the amplitude of the long-term changes in solar irradiance. Here we review recent efforts to improve irradiance reconstructions on time scales longer than the solar cycle and to reduce the existing uncertainty in the magnitude of the long-term variability. In particular, we highlight the potential of using 3D magnetohydrodynamical simulations of the solar atmosphere as input to more physical irradiance models and of historical full-disc Ca II K observations encrypting direct facular information back to 1892. △ Less

Submitted 16 October, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 17 pages, 8 figures, accepted for publication in JASTP

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

Authors: Theodosios Chatzistergos, Natalie A. Krivova, Ilaria Ermolli

Abstract: First such observations were made in 1892 and since then various sites around the world have carried out regular observations, with Kodaikanal, Meudon, Mt Wilson, and Coimbra being some of the most prominent ones. By now, Ca II K observations from over 40 different sites allow an almost complete daily coverage of the last century. Ca II K images provide direct information on plage and network regi… ▽ More First such observations were made in 1892 and since then various sites around the world have carried out regular observations, with Kodaikanal, Meudon, Mt Wilson, and Coimbra being some of the most prominent ones. By now, Ca II K observations from over 40 different sites allow an almost complete daily coverage of the last century. Ca II K images provide direct information on plage and network regions on the Sun and, through their connection to solar surface magnetic field, offer an excellent opportunity to study solar magnetism over more than a century. This makes them also extremely important, among others, for solar irradiance reconstructions and studies of the solar influence on Earth's climate. However, these data also suffer from numerous issues, which for a long time have hampered their analysis. Without properly addressing these issues, Ca II K data cannot be used to their full potential. Here, we first provide an overview of the currently known Ca II K data archives and sources of the inhomogeneities in the data, before discussing existing processing techniques, followed by a recap of the main results derived with such data so far. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: 47 pages, 13 figures, accepted for publication in Frontiers in Astronomy and Space Sciences

Superflares on solar-like stars: A new method for identifying the true flare sources in photometric surveys

Authors: Valeriy Vasilyev, Timo Reinhold, Alexander I. Shapiro, Natalie A. Krivova, Ilya Usoskin, Benjamin T. Montet, Sami K. Solanki, Laurent Gizon

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

Submitted 28 September, 2022; originally announced September 2022.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 668, A167 (2022)

Scrutinizing the relationship between plage areas and sunspot areas and numbers

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Natalie A. Krivova, Teresa Barata, Sara Carvalho, Jean-Marie Malherbe

Abstract: Studies and reconstructions of past solar activity require data on sunspots as well as faculae/plage and network. Such data are also important for understanding the magnetic activity and variability of the Sun and Sun-like stars. The longest available direct faculae/plage datasets are white-light facular and Ca II K observations going back to 1874 and 1892, respectively. Prior to that time, the on… ▽ More Studies and reconstructions of past solar activity require data on sunspots as well as faculae/plage and network. Such data are also important for understanding the magnetic activity and variability of the Sun and Sun-like stars. The longest available direct faculae/plage datasets are white-light facular and Ca II K observations going back to 1874 and 1892, respectively. Prior to that time, the only direct data available are for sunspots. We reassess the relationship between plage areas and sunspot records (areas and numbers) since 1892, to allow reconstructions of facular/plage areas which can be employed for studies going further back in time, i.e. over the period when solely sunspot observations are available. We use the plage areas derived from 38 consistently processed Ca II K archives as well as the plage area composite based on these archives. We find the relationship between plage and sunspot areas and numbers to be well represented by a power law function. We further find that the relationships depend on the bandwidth and the solar cycle strength. The reconstructions with a power law relationship are in good agreement with the original plage area series, whereas employment of a cycle-strength-dependent relationship improves the reconstructions only marginally. Performing the same analysis on other previously published plage area series, usually derived from a single archive with diverse processing techniques, returns different results when using different time series. This highlights the importance of applying a consistent processing to the various archives and demonstrates the uncertainties introduced by using previously published series. Our results have implications for past solar activity and irradiance reconstructions as well as for stellar activity studies, which sometimes assume a linear dependence between plage and sunspot areas. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: 22 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 667, A167 (2022)

Faculae cancel out on the surfaces of active Suns

Authors: N. -E. Nèmec, A. I. Shapiro, E. Işık, K. Sowmya, S. K. Solanki, N. A. Krivova, R. H. Cameron, L. Gizon

Abstract: Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they… ▽ More Surfaces of the Sun and other cool stars are filled with magnetic fields, which are either seen as dark compact spots or more diffuse bright structures like faculae. Both hamper detection and characterisation of exoplanets, affecting stellar brightness and spectra, as well as transmission spectra. However, the expected facular and spot signals in stellar data are quite different, for instance they have distinct temporal and spectral profiles. Consequently, corrections of stellar data for magnetic activity can greatly benefit from the insight on whether the stellar signal is dominated by spots or faculae. Here, we utilise a surface flux transport model (SFTM) to show that more effective cancellation of diffuse magnetic flux associated with faculae leads to spot area coverages increasing faster with stellar magnetic activity than that by faculae. Our calculations explain the observed dependence between solar spot and facular area coverages and allow its extension to stars more active than the Sun. This extension enables anticipating the properties of stellar signal and its more reliable mitigation, leading to a more accurate characterisation of exoplanets and their atmospheres. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 10 pages, 6 Figures, accepted for publication in ApJL

Predictions of Astrometric Jitter for Sun-like Stars. III. Fast Rotators

Authors: K. Sowmya, N. -E. Nèmec, A. I. Shapiro, E. Işık, N. A. Krivova, S. K. Solanki

Abstract: A breakthrough in exoplanet detections is foreseen with the unprecedented astrometric measurement capabilities offered by instrumentation aboard Gaia space observatory. Besides, astrometric discoveries of exoplanets are expected from the planned space mission, Small-JASMINE. In this setting, the present series of papers focuses on estimating the effect of magnetic activity of G2V-type host stars o… ▽ More A breakthrough in exoplanet detections is foreseen with the unprecedented astrometric measurement capabilities offered by instrumentation aboard Gaia space observatory. Besides, astrometric discoveries of exoplanets are expected from the planned space mission, Small-JASMINE. In this setting, the present series of papers focuses on estimating the effect of magnetic activity of G2V-type host stars on the astrometric signal. This effect interferes with the astrometric detections of Earth-mass planets. While the first two papers considered stars rotating at the solar rotation rate, this paper focuses on stars having solar effective temperature and metallicity but rotating faster than the Sun, and consequently more active. By simulating the distribution of active regions on such stars using the Flux Emergence And Transport model, we show that the contribution of magnetic activity to the astrometric measurements becomes increasingly significant with increasing rotation rates. We further show that the jitter for the most variable periodic Kepler stars is high enough to be detected by Gaia. Furthermore, due to a decrease in the facula-to-spot area ratio for more active stars, the magnetic jitter is found to be spot-dominated for rapid rotators. Our simulations of the astrometric jitter has the potential to aid the interpretation of data from Gaia and upcoming space astrometry missions. △ Less

Submitted 28 June, 2022; v1 submitted 15 June, 2022; originally announced June 2022.

Comments: 21 pages, 20 figures, accepted for publication in ApJ

The relationship between bipolar magnetic regions and their sunspots

Authors: K. L. Yeo, S. K. Solanki, N. A. Krivova, J. Jiang

Abstract: The relationship between bipolar magnetic regions (BMRs) and their sunspots is an important property of the solar magnetic field, but it is not well constrained. One consequence is that it is a challenge for surface flux transport models (SFTMs) based on sunspot observations to determine the details of BMR emergence, which they require as input, from such data. We aimed to establish the relationsh… ▽ More The relationship between bipolar magnetic regions (BMRs) and their sunspots is an important property of the solar magnetic field, but it is not well constrained. One consequence is that it is a challenge for surface flux transport models (SFTMs) based on sunspot observations to determine the details of BMR emergence, which they require as input, from such data. We aimed to establish the relationship between the amount of magnetic flux in newly emerged BMRs and the area of the enclosed sunspots. Earlier attempts to constrain BMR magnetic flux were hindered by the fact that there is no proper database of the magnetic and physical properties of newly emerged BMRs currently available. We made use of the empirical model of the relationship between the disc-integrated facular and network magnetic flux and the total surface coverage by sunspots reported in a recent study. The structure of the model is such that it enabled us to establish, from these disc-integrated quantities, an empirical relationship between the magnetic flux and sunspot area of individual newly emerged BMRs, circumventing the lack of any proper BMR database. Applying the constraint on BMR magnetic flux derived here to an established SFTM retained its ability to replicate various independent datasets and the correlation between the model output polar field at the end of each cycle and the observed strength of the following cycle. The SFTM output indicates that facular and network magnetic flux rises with increasing sunspot magnetic flux at a slowing rate such that it appears to gradually saturate, analogous to earlier studies. The activity dependence of the ratio of facular and network flux to sunspot flux is consistent with the findings of recent studies: although the Sun is faculae-dominated, it is only marginally so as facular and network brightening and sunspot darkening appear to be closely balanced. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: In press

Journal ref: A&A (2021)

Reconstructing solar irradiance from historical Ca II K observations. I. Method and its validation

Authors: Theodosios Chatzistergos, Natalie A. Krivova, Ilaria Ermolli, Kok Leng Yeo, Sudip Mandal, Sami K. Solanki, Greg Kopp, Jean-Marie Malherbe

Abstract: Knowledge of solar irradiance variability is critical to Earth's climate models and understanding the solar influence on Earth's climate. Direct solar irradiance measurements are only available since 1978. Reconstructions of past variability typically rely on sunspot data. These provide only indirect information on the facular and network regions, which are decisive contributors to irradiance vari… ▽ More Knowledge of solar irradiance variability is critical to Earth's climate models and understanding the solar influence on Earth's climate. Direct solar irradiance measurements are only available since 1978. Reconstructions of past variability typically rely on sunspot data. These provide only indirect information on the facular and network regions, which are decisive contributors to irradiance variability on timescales of the solar cycle and longer. Our ultimate goal is to reconstruct past solar irradiance variations using historical full-disc Ca II K observations to describe the facular contribution independently of sunspot observations. Here, we develop the method and test it extensively by using modern CCD-based Ca II K observations and carry out initial tests on two photographic archives. We employ carefully reduced and calibrated Ca II K images from 13 datasets, such as those from the Meudon, Mt Wilson, and Rome observatories. We convert them to unsigned magnetograms and then use them as input to the adapted SATIRE model to reconstruct TSI variations over the period 1978-2019, for which direct irradiance measurements are available. The reconstructed TSI from the analysed Ca II K archives agrees well with direct TSI measurements and existing reconstructions. The model also returns good results on data taken with different bandpasses and images with low spatial resolution. Historical Ca II K archives suffer from numerous inconsistencies, but we show that these archives can still be used to reconstruct TSI with reasonable accuracy provided the observations are accurately processed. By using the unsigned magnetograms of the Sun reconstructed from high-quality Ca II K observations as input into the SATIRE model, we can reconstruct solar irradiance variations nearly as accurately as from directly recorded magnetograms. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 25 pages, 15 figures, accepted for publication in A&A

Journal ref: A&A 656, A104 (2021)

Predictions of Astrometric Jitter for Sun-like Stars. II. Dependence on Inclination, Metallicity, and Active-Region Nesting

Authors: K. Sowmya, N. -E. Nèmec, A. I. Shapiro, E. Işık, V. Witzke, A. Mints, N. A. Krivova, S. K. Solanki

Abstract: Ultra-precise astrometry from the Gaia mission is expected to lead to astrometric detections of more than 20,000 exoplanets in our Galaxy. One of the factors that could hamper such detections is the astrometric jitter caused by the magnetic activity of the planet host stars. In our previous study, we modeled astrometric jitter for the Sun observed equator-on. In this work, we generalize our model… ▽ More Ultra-precise astrometry from the Gaia mission is expected to lead to astrometric detections of more than 20,000 exoplanets in our Galaxy. One of the factors that could hamper such detections is the astrometric jitter caused by the magnetic activity of the planet host stars. In our previous study, we modeled astrometric jitter for the Sun observed equator-on. In this work, we generalize our model and calculate the photocenter jitter as it would be measured by the Gaia and Small-JASMINE missions for stars with solar rotation rate and effective temperature, but with various values of the inclination angle of the stellar rotation axis. In addition, we consider the effect of metallicity and of nesting of active regions (i.e. the tendency of active regions to emerge in the vicinity of each other). We find that, while the jitter of stars observed equator-on does not have any long-term trends and can be easily filtered out, the photocenters of stars observed out of their equatorial planes experience systematic shifts over the course of the activity cycle. Such trends allow the jitter to be detected with continuous measurements, in which case it can interfere with planet detectability. An increase in the metallicity is found to increase the jitter caused by stellar activity. Active-region nesting can further enhance the peak-to-peak amplitude of the photocenter jitter to a level that could be detected by Gaia. △ Less

Submitted 3 July, 2021; originally announced July 2021.

Comments: 24 pages, 20 figures in the main body and 4 in the appendix, accepted for publication in ApJ

Journal ref: 2021, ApJ, 919, 94

On the size distribution of spots within sunspot groups

Authors: Sudip Mandal, Natalie A. Krivova, Robert Cameron, Sami K. Solanki

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

Submitted 8 April, 2021; originally announced April 2021.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 652, A9 (2021)

Modelling the evolution of the Sun's open and total magnetic flux

Authors: N. A. Krivova, S. K. Solanki, B. Hofer, C. -J. Wu, I. G. Usoskin, R. Cameron

Abstract: Solar activity in all its varied manifestations is driven by the magnetic field. Particularly important for many purposes are two global quantities, the Sun's total and open magnetic flux, which can be computed from sunspot number records using models. Such sunspot-driven models, however, do not take into account the presence of magnetic flux during grand minima, such as the Maunder minimum. Here… ▽ More Solar activity in all its varied manifestations is driven by the magnetic field. Particularly important for many purposes are two global quantities, the Sun's total and open magnetic flux, which can be computed from sunspot number records using models. Such sunspot-driven models, however, do not take into account the presence of magnetic flux during grand minima, such as the Maunder minimum. Here we present a major update of a widely used simple model, which now takes into account the observation that the distribution of all magnetic features on the Sun follows a single power law. The exponent of the power law changes over the solar cycle. This allows for the emergence of small-scale magnetic flux even when no sunspots are present for multiple decades and leads to non-zero total and open magnetic flux also in the deepest grand minima, such as the Maunder minimum, thus overcoming a major shortcoming of the earlier models. The results of the updated model compare well with the available observations and reconstructions of the solar total and open magnetic flux. This opens up the possibility of improved reconstructions of sunspot number from time series of cosmogenic isotope production rate. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Journal ref: A&A 650, A70 (2021)

Modeling stellar Ca II H & K emission variations. I. Effect of inclination on the S-index

Authors: K. Sowmya, A. I. Shapiro, V. Witzke, N. -E. Nèmec, T. Chatzistergos, K. L. Yeo, N. A. Krivova, S. K. Solanki

Abstract: The emission in the near ultraviolet Ca II H & K lines is modulated by stellar magnetic activity. Although this emission, quantified via the S-index, has been serving as a prime proxy of stellar magnetic activity for several decades, many aspects of the complex relation between stellar magnetism and Ca II H & K emission are still unclear. The amount of measured Ca II H & K emission is suspected to… ▽ More The emission in the near ultraviolet Ca II H & K lines is modulated by stellar magnetic activity. Although this emission, quantified via the S-index, has been serving as a prime proxy of stellar magnetic activity for several decades, many aspects of the complex relation between stellar magnetism and Ca II H & K emission are still unclear. The amount of measured Ca II H & K emission is suspected to be affected not only by the stellar intrinsic properties but also by the inclination angle of the stellar rotation axis. Until now such an inclination effect on S-index has remained largely unexplored. To fill this gap, we develop a physics-based model to calculate S-index, focusing on the Sun. Using the distributions of solar magnetic features derived from observations together with Ca II H & K spectra synthesized in non-local thermodynamic equilibrium, we validate our model by successfully reconstructing the observed variations of solar S-index over four activity cycles. Further, using the distribution of magnetic features over the visible solar disk obtained from surface flux transport simulations, we obtain S-index time series dating back to 1700 and investigate the effect of inclination on S-index variability, both on the magnetic activity cycle and the rotational timescales. We find that when going from an equatorial to a pole-on view, the amplitude of S-index variations decreases weakly on the activity cycle timescale and strongly on the rotational timescale (by about 22% and 81%, respectively, for a cycle of intermediate strength). The absolute value of S-index depends only weakly on the inclination. We provide analytical expressions that model such dependencies. △ Less

Submitted 25 March, 2021; originally announced March 2021.

Comments: 23 pages, 13 figures in the main body and 2 in the appendix, accepted for publication in ApJ

Intensity contrast of solar network and faculae -- II. Implications for solar irradiance modelling

Authors: K. L. Yeo, N. A. Krivova

Abstract: We aim to gain insight into the effect of network and faculae on solar irradiance from their apparent intensity. Taking full-disc observations from the Solar Dynamics Observatory, we examined the intensity contrast of network and faculae in the continuum and core of the Fe I 6173 Å line and 1700 Å, including the variation with magnetic flux density, distance from disc centre, nearby magnetic field… ▽ More We aim to gain insight into the effect of network and faculae on solar irradiance from their apparent intensity. Taking full-disc observations from the Solar Dynamics Observatory, we examined the intensity contrast of network and faculae in the continuum and core of the Fe I 6173 Å line and 1700 Å, including the variation with magnetic flux density, distance from disc centre, nearby magnetic fields, and time. The brightness of network and faculae is believed to be suppressed by nearby magnetic fields from its effect on convection. The difference in intensity contrast between the quiet-Sun network and active region faculae, noted by various studies, arises because active regions are more magnetically crowded and is not due to any fundamental physical differences between network and faculae. These results highlight that solar irradiance models need to include the effect of nearby magnetic fields on network and faculae brightness. We found evidence that suggests that departures from local thermal equilibrium (LTE) might have limited effect on intensity contrast. This could explain why solar irradiance models that are based on the intensity contrast of solar surface magnetic features calculated assuming LTE reproduce the observed spectral variability even where the LTE assumption breaks down. Certain models of solar irradiance employ chromospheric indices as direct indications of the effect of network and faculae on solar irradiance. Based on past studies of the Ca II K line and on the intensity contrast measurements derived here, we show that the fluctuations in chromospheric emission from network and faculae are a reasonable estimate of the emission fluctuations in the middle photosphere, but not of those in the lower photosphere. The data set, which extends from 2010 to 2018, indicates that intensity contrast was stable to about 3% in this period. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Journal ref: A&A 624 A135 (2019)

The Dimmest State of the Sun

Authors: K. L. Yeo, S. K. Solanki, N. A. Krivova, M. Rempel, L. S. Anusha, A. I. Shapiro, R. V. Tagirov, V. Witzke

Abstract: How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm… ▽ More How the solar electromagnetic energy entering the Earth's atmosphere varied since pre-industrial times is an important consideration in the climate change debate. Detrimental to this debate, estimates of the change in total solar irradiance (TSI) since the Maunder minimum, an extended period of weak solar activity preceding the industrial revolution, differ markedly, ranging from a drop of 0.75 Wm-2 to a rise of 6.3 Wm-2. Consequently, the exact contribution by solar forcing to the rise in global temperatures over the past centuries remains inconclusive. Adopting a novel approach based on state-of-the-art solar imagery and numerical simulations, we establish the TSI level of the Sun when it is in its least-active state to be 2.0 +/- 0.7 Wm-2 below the 2019 level. This means TSI could not have risen since the Maunder minimum by more than this amount, thus restricting the possible role of solar forcing in global warming. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: Supporting information: https://agupubs.onlinelibrary.wiley.com/action/downloadSupplement?doi=10.1029%2F2020GL090243&file=grl61316-sup-0001-Text_SI-S01.pdf

Journal ref: Geophysical Research Letters, 47, e2020GL090243 (2020)

Predictions of astrometric jitter for Sun-like stars. I. The model and its application to the Sun as seen from the ecliptic

Authors: Alexander I. Shapiro, Sami K. Solanki, Natalie A. Krivova

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

Submitted 22 December, 2020; originally announced December 2020.

Comments: Accepted for publication in ApJ

Amplification of brightness variability by active-region nesting in solar-like stars

Authors: Emre Işık, Alexander I. Shapiro, Sami K. Solanki, Natalie A. Krivova

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

Submitted 11 December, 2020; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: Typo in Equation 1 corrected in this version. Erratum accepted by ApJL

Journal ref: ApJL 901, L12 (2020)

Modelling solar irradiance from ground-based photometric observations

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Fabrizio Giorgi, Natalie A. Krivova, Cosmin Constantin Puiu

Abstract: Total solar irradiance (TSI) has been monitored from space since 1978. The measurements show a prominent variability in phase with the solar cycle, as well as fluctuations on timescales shorter than a few days. However, the measurements were done by multiple and usually relatively short-lived missions making the possible long-term trend in the TSI highly uncertain. While the variability in the UV… ▽ More Total solar irradiance (TSI) has been monitored from space since 1978. The measurements show a prominent variability in phase with the solar cycle, as well as fluctuations on timescales shorter than a few days. However, the measurements were done by multiple and usually relatively short-lived missions making the possible long-term trend in the TSI highly uncertain. While the variability in the UV irradiance is clearly in-phase with the solar cycle, the phase of the variability in the visible range has been debated. In this paper, we aim at getting an insight into the long-term trend of TSI since 1996 and the phase of the solar irradiance variations in the visible part of the spectrum. We use independent ground-based full-disc photometric observations in Ca~II~K and continuum from the Rome and San Fernando observatories to compute the TSI since 1996. We follow the empirical San Fernando approach based on the photometric sum index. We find a weak declining trend in the TSI of -7.8$^{+4.9}_{-0.8}\times10^{-3}$ Wm$^{-2}$y$^{-1}$ between the 1996 and 2008 activity minima, while between 2008 and 2019 the reconstructed TSI shows no trend to a marginally decreasing (but statistically insignificant) trend of -0.1$^{+0.25}_{-0.02}\times10^{-3}$ Wm$^{-2}$y$^{-1}$. The reference TSI series used for the reconstruction does not significantly affect the determined trend. The variation in the blue continuum (409.2 nm) is rather flat, while the variation in the red continuum (607.1 nm) is marginally in anti-phase, although this result is extremely sensitive to the accurate assessment of the quiet Sun level in the images. These results provide further insights into the long-term variation of the total solar irradiance. The amplitude of the variations in the visible is below the uncertainties of the processing, which prevents an assessment of the phase of the variations. △ Less

Submitted 24 August, 2020; originally announced August 2020.

Comments: 17 pages, 8 figures, accepted for publication in Journal of Space Weather and Space Climate as part of the topical issue "Space climate: The past and future of solar activity"

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"

Authors: T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gómez

Abstract: 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. 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. △ Less

Submitted 9 July, 2020; originally announced July 2020.

How faculae and network relate to sunspots, and the implications for solar and stellar brightness variations

Authors: K. L. Yeo, S. K. Solanki, N. A. Krivova

Abstract: 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… ▽ More 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. △ Less

Submitted 5 October, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Analysis of full disc Ca II K spectroheliograms III. Plage area composite series covering 1892-2019

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Natalie A. Krivova, Sami K. Solanki, Dipankar Banerjee, Teresa Barata, Marcel Belik, Ricardo Gafeira, Adriana Garcia, Yoichiro Hanaoka, Manjunath Hegde, Jan Klimeš, Viktor V. Korokhin, Ana Lourenço, Jean-Marie Malherbe, Gennady P. Marchenko, Nuno Peixinho, Takashi Sakurai, Andrey G. Tlatov

Abstract: We derive the plage area evolution over the last 12 solar cycles employing data from all Ca II K archives available publicly in digital form known to us, including several as yet unexplored Ca II K archives. We analyse more than 290,000 full-disc Ca II K observations from 43 datasets spanning the period 1892-2019. All images were consistently processed with an automatic procedure that performs the… ▽ More We derive the plage area evolution over the last 12 solar cycles employing data from all Ca II K archives available publicly in digital form known to us, including several as yet unexplored Ca II K archives. We analyse more than 290,000 full-disc Ca II K observations from 43 datasets spanning the period 1892-2019. All images were consistently processed with an automatic procedure that performs the photometric calibration (if needed) and the limb-darkening compensation. The processing also accounts for artefacts plaguing many of the images, including some very specific artefacts such as bright arcs found in Kyoto and Yerkes data. We have produced a plage area time-series from each analysed dataset. We found that the differences between the plage areas derived from individual archives are mainly due to the differences in the central wavelength and the bandpass used to acquire the data at the various sites. We have empirically cross-calibrated and combined the results obtained from each dataset to produce a composite series of plage areas. "Backbone" series are used to bridge all the series together. We have also shown that the selection of the backbone series has little effect on the final plage area composite. We have quantified the uncertainty of determining the plage areas with our processing due to shifts in the central wavelength and found it to be less than 0.01 in fraction of the solar disc for the average conditions found on historical data. We also found the variable seeing conditions during the observations to slightly increase the plage areas during activity maxima. We provide the so far most complete time series of plage areas based on corrected and calibrated historical and modern Ca II K images. Consistent plage areas are now available on 88% of all days from 1892 onwards and on 98% from 1907 onwards. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Comments: 23 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 639, A88 (2020)

The Sun is less active than other solar-like stars

Authors: Timo Reinhold, Alexander I. Shapiro, Sami K. Solanki, Benjamin T. Montet, Nathalie A. Krivova, Robert H. Cameron, Eliana M. Amazo-Gomez

Abstract: 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 variabi… ▽ More 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. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Comments: Accepted for publication in Science. 3 (main) + 10 (supplementary) figures

Journal ref: Science 01 May 2020: Vol. 368, Issue 6490, pp. 518-521

Sunspot area catalogue revisited: Daily cross-calibrated areas since 1874

Authors: Sudip Mandal, Natalie A. Krivova, Sami K. Solanki, Nimesh Sinha, Dipankar Banerjee

Abstract: 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 r… ▽ More 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. △ Less

Submitted 30 April, 2020; originally announced April 2020.

Comments: Accepted for publication in Astronomy and Astrophysics. An additional file (animation) is available at https://www.dropbox.com/s/armawopcxt8kmb9/Mandal_2020_butterfly_diagram.zip?dl=0

Journal ref: A&A 640, A78 (2020)

Connecting measurements of solar and stellar brightness variations

Authors: N. -E. Nèmec, E. Işık, A. I. Shapiro, S. K. Solanki, N. A. Krivova, Y. Unruh

Abstract: 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… ▽ More 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. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 9 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 638, A56 (2020)

Power spectra of solar brightness variations at different inclinations

Authors: N. -E. Nèmec, A. I. Shapiro, N. A. Krivova, S. K. Solanki, R. V. Tagirov, R. H. Cameron, S. Dreizler

Abstract: 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, h… ▽ More 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. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: 11 pages, 9 figures plus 3 additional figures in the appendix, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 636, A43 (2020)

Solar-cycle irradiance variations over the last four billion years

Authors: Anna V. Shapiro, Alexander I. Shapiro, Laurent Gizon, Natalie A. Krivova, Sami K. Solanki

Abstract: 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 variabilit… ▽ More 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. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: 8 pages, 6 figures, accepted publication for Astronomy&Astrophysics

Rotation periods from the inflection point in the power spectrum of stellar brightness variations: II. The Sun

Authors: E. M. Amazo-Gómez, A. I. Shapiro, S. K. Solanki, N. A. Krivova, G. Kopp, T. Reinhold, M. Oshagh, A. Reiners

Abstract: 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-li… ▽ More 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. △ Less

Submitted 9 February, 2020; originally announced February 2020.

Journal ref: A&A 636, A69 (2020)

Effect of metallicity on the detectability of rotational periods in solar-like stars

Authors: V. Witzke, T. Reinhold, A. I. Shapiro, N. A. Krivova, S. K. Solanki

Abstract: 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… ▽ More 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. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 9 pages, 7 figures, accepted in A&A

Journal ref: A&A 634, L9 (2020)

Inflection point in the power spectrum of stellar brightness variations. I. The model

Authors: A. I. Shapiro, E. M. Amazo-Gómez, N. A. Krivova, S. K. Solanki

Abstract: 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 ex… ▽ More 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. △ Less

Submitted 18 October, 2019; originally announced October 2019.

Comments: 18 pages, 20 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 633, A32 (2020)

Readdressing the UV solar variability with SATIRE-S: non-LTE effects

Authors: R. V. Tagirov, A. I. Shapiro, N. A. Krivova, Y. C. Unruh, K. L. Yeo, S. K. Solanki

Abstract: 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 variabili… ▽ More 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. △ Less

Submitted 25 September, 2019; originally announced September 2019.

Comments: 8 pages, 6 figures

Journal ref: A&A 631, A178 (2019)

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

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Sami K. Solanki, Natalie A. Krivova, Dipankar Banerjee, Bibhuti K. Jha, Subhamoy Chatterjee

Abstract: Full-disc Ca II K photographic observations of the Sun carry direct information about the evolution of solar-plage regions for more than a century and are therefore a unique dataset for solar-activity studies. For a long time Ca II K observations were barely explored, but recent digitisations of multiple archives have allowed their extensive analysis. However, various studies have reported diverse… ▽ More Full-disc Ca II K photographic observations of the Sun carry direct information about the evolution of solar-plage regions for more than a century and are therefore a unique dataset for solar-activity studies. For a long time Ca II K observations were barely explored, but recent digitisations of multiple archives have allowed their extensive analysis. However, various studies have reported diverse results partly due to the insufficient quality of the digitised data. Furthermore, inhomogeneities have been identified within the individual archives, which, at least partly, could be due to the digitisation. As a result, some of the archives, e.g. that from the Kodaikanal observatory, were re-digitised. The results obtained by different authors who analysed the data from the new digitisation of the Kodaikanal archive differ from each other as well as from those derived from the old digitisation. Since the data were processed and analysed using different techniques, it is not clear, however, whether the differences are due to the digitisation or the processing of the data. To understand the reasons for such discrepancies, we analyse here the data from the two most recent digitisations of this archive. We use the same techniques to consistently process the images from both archives and to derive the plage areas from them. Some issues have been identified in both digitisations, implying that they are intrinsic characteristics of the data. Moreover, errors in timing of the observations plague both digitisations. Overall, the most recent 16-bit digitisation offers an improvement over the earlier 8-bit one. It also includes considerably more data and should be preferred. △ Less

Submitted 15 August, 2019; originally announced August 2019.

Comments: 28 pages, 17 figures, accepted for publication in solar physics

Opacity distribution functions for stellar spectra synthesis

Authors: M. Cernetic, A. I. Shapiro, V. Witzke, N. A. Krivova, S. K. Solanki, R. V. Tagirov

Abstract: 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 th… ▽ More 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. △ Less

Submitted 24 June, 2019; v1 submitted 7 June, 2019; originally announced June 2019.

Comments: 8 pages, 11 figures

Journal ref: A&A 627, A157 (2019)

Recovering the unsigned photospheric magnetic field from Ca II K observations

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Sami K. Solanki, Natalie A. Krivova, Fabrizio Giorgi, Kok Leng Yeo

Abstract: We reassess the relationship between the photospheric magnetic field strength and the Ca II K intensity for a variety of surface features as a function of the position on the disc and the solar activity level. This relationship can be used to recover the unsigned photospheric magnetic field from images recorded in the core of Ca II K line. We have analysed 131 pairs of high-quality, full-disc, nea… ▽ More We reassess the relationship between the photospheric magnetic field strength and the Ca II K intensity for a variety of surface features as a function of the position on the disc and the solar activity level. This relationship can be used to recover the unsigned photospheric magnetic field from images recorded in the core of Ca II K line. We have analysed 131 pairs of high-quality, full-disc, near-co-temporal observations from SDO/HMI and Rome/PSPT spanning half a solar cycle. To analytically describe the observationally-determined relation, we considered three different functions: a power law with an offset, a logarithmic function, and a power law function of the logarithm of the magnetic flux density. We used the obtained relations to reconstruct maps of the line-of-sight component of the unsigned magnetic field (unsigned magnetograms) from Ca II K observations, which were then compared to the original magnetograms. We find that both power-law functions represent the data well, while the logarithmic function is good only for quiet periods. We see no significant variation over the solar cycle or over the disc in the derived fit parameters, independently of the function used. We find that errors in the independent variable, usually not accounted for, introduce attenuation bias. To address this, we binned the data with respect to the magnetic field strength and Ca II K contrast separately and derived the relation for the bisector of the two binned curves. The reconstructed unsigned magnetograms show good agreement with the original ones. RMS differences are less than 90 G. The results were unaffected by the stray-light correction of the SDO/HMI and Rome/PSPT data. Our results imply that Ca~II~K observations, accurately processed and calibrated, can be used to reconstruct unsigned magnetograms by using the relations derived in our study. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: 18 pages, 22 figures, accepted in A&A

Journal ref: A&A 626, A114 (2019)

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

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Mariachiara Falco, Fabrizio Giorgi, Salvo L. Guglielmino, Natalie A. Krivova, Paolo Romano, Sami K. Solanki

Abstract: Here we present the little explored Ca II K archives from the Rome and the Catania observatories and analyse the digitised images from these archives to derive plage areas. Here we present the little explored Ca II K archives from the Rome and the Catania observatories and analyse the digitised images from these archives to derive plage areas. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: 5 pages, 3 figures, to be published in "Nuovo Cimento C" as proceeding of the Third Meeting of the Italian Solar and Heliospheric Community

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

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Natalie A. Krivova, Sami K. Solanki

Abstract: Reconstructions of past irradiance variations require suitable data on solar activity. The longest direct proxy is the sunspot number, and it has been most widely employed for this purpose. These data, however, only provide information on the surface magnetic field emerging in sunspots, while a suitable proxy of the evolution of the bright magnetic features, specifically faculae/plage and network,… ▽ More Reconstructions of past irradiance variations require suitable data on solar activity. The longest direct proxy is the sunspot number, and it has been most widely employed for this purpose. These data, however, only provide information on the surface magnetic field emerging in sunspots, while a suitable proxy of the evolution of the bright magnetic features, specifically faculae/plage and network, is missing. This information can potentially be extracted from the historical full-disc observations in the Ca II K line. We have analysed over 100,000 historical images from 8 digitised photographic archives of the Arcetri, Kodaikanal, McMath-Hulbert, Meudon, Mitaka, Mt Wilson, Schauinsland, and Wendelstein observatories, as well as one archive of modern observations from the Rome/PSPT. The analysed data cover the period 1893--2018. We first performed careful photometric calibration and compensation for the centre-to-limb variation, and then segmented the images to identify plage regions. This has been consistently applied to both historical and modern observations. The plage series derived from different archives are generally in good agreement with each other. However, there are also clear deviations that most likely hint at intrinsic differences in the data and their digitisation. We showed that accurate image processing significantly reduces errors in the plage area estimates. Accurate photometric calibration also allows precise plage identification on images from different archives without the need to arbitrarily adjust the segmentation parameters. Finally, by comparing the plage area series from the various records, we found the conversion laws between them. This allowed us to produce a preliminary composite of the plage areas obtained from all the datasets studied here. This is a first step towards an accurate assessment of the long-term variation of plage regions. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: 30 pages, 22 figures, accepted in A&A

Journal ref: A&A 625, A69 (2019)

Activity variation driven by flux emergence and transport on Sun-like stars

Authors: Emre Işık, Sami K. Solanki, Natalie A. Krivova, Alexander I. Shapiro

Abstract: 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… ▽ More 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. △ Less

Submitted 21 December, 2018; originally announced December 2018.

Comments: 3 pages, 3 figures; submitted to Proceedings of the 20th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun; ed. S. Wolk. This (newest) version has a reference update

Solar total and spectral irradiance reconstruction over the last 9000 years

Authors: C. -J. Wu, N. A. Krivova, S. K. Solanki, I. G. Usoskin

Abstract: 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 irradian… ▽ More 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. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Forward modelling of brightness variations in Sun-like stars I. Emergence and surface transport of magnetic flux

Authors: Emre Işık, Sami K. Solanki, Natalie A. Krivova, Alexander I. Shapiro

Abstract: 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 variatio… ▽ More 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Ωおーむ_\odot$, (2) the tilt angles of freshly emerged bipolar regions show stronger variations with latitude. Polar spots can form at $8Ωおーむ_\odot$ by accumulation of follower-polarity flux from decaying bipolar regions. From $4Ωおーむ_\odot$ to $8Ωおーむ_\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Ωおーむ_\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. △ Less

Submitted 23 October, 2018; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: 16 pages, 13 figures. Astron. & Astrophys. (in press); minor language corrections made

Journal ref: A&A 620, A177 (2018)

From Solar to Stellar Brightness Variations: The Effect of Metallicity

Authors: V. Witzke, A. I. Shapiro, S. K. Solanki, N. A. Krivova, W. Schmutz

Abstract: 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 chan… ▽ More 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ömgren (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. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: 12 pages, 12 figures, accepted in A&A

Journal ref: A&A 619, A146 (2018)

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

Authors: Ilaria Ermolli, Theodosios Chatzistergos, Natalie A. Krivova, Sami K. Solanki

Abstract: Several observatories around the globe started regular full-disc imaging of the solar atmosphere in the Ca II K line in the early decades of the 20th century. These observations are continued today at a few sites with either old spectroheliographs or modern telescopes equipped with narrow-band filters. The Ca II K time series are unique in representing long-term variations of the Sun's chromospher… ▽ More Several observatories around the globe started regular full-disc imaging of the solar atmosphere in the Ca II K line in the early decades of the 20th century. These observations are continued today at a few sites with either old spectroheliographs or modern telescopes equipped with narrow-band filters. The Ca II K time series are unique in representing long-term variations of the Sun's chromospheric magnetic field. However, meaningful results from their analysis require accurate processing of the available data and robust merging of the information stored in different archives. This paper provides an overview of the historical and modern full-disc Ca II K observations, with focus on their quality and the main results obtained from their analysis over the last decade. △ Less

Submitted 11 May, 2018; originally announced May 2018.

Comments: 6 pages, 2 figures

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

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Natalie A. Krivova, Sami K. Solanki

Abstract: We address the importance of historical full disc Ca II K spectroheliograms for solar activity and irradiance reconstruction studies. We review our work on processing such data to enable them to be used in irradiance reconstructions. We also present our preliminary estimates of the plage areas from five of the longest available historical Ca II K archives. We address the importance of historical full disc Ca II K spectroheliograms for solar activity and irradiance reconstruction studies. We review our work on processing such data to enable them to be used in irradiance reconstructions. We also present our preliminary estimates of the plage areas from five of the longest available historical Ca II K archives. △ Less

Submitted 10 May, 2018; originally announced May 2018.

Comments: 5 pages, 3 figures

The nature of solar brightness variations

Authors: A. I. Shapiro, S. K. Solanki, N. A. Krivova, R. H. Cameron, K. L. Yeo, W. K. Schmutz

Abstract: 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… ▽ More 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. △ Less

Submitted 26 November, 2017; v1 submitted 11 November, 2017; originally announced November 2017.

Comments: This is the submitted version of the paper published in Nature Astronomy

Analysis of full disc Ca II K spectroheliograms. I. Photometric calibration and CLV compensation

Authors: Theodosios Chatzistergos, Ilaria Ermolli, Sami K. Solanki, Natalie A. Krivova

Abstract: Historical Ca II K spectroheliograms (SHG) are unique in representing long-term variations of the solar chromospheric magnetic field. They usually suffer from numerous problems and lack photometric calibration. Thus accurate processing of these data is required to get meaningful results from their analysis. In this paper we aim at developing an automatic processing and photometric calibration meth… ▽ More Historical Ca II K spectroheliograms (SHG) are unique in representing long-term variations of the solar chromospheric magnetic field. They usually suffer from numerous problems and lack photometric calibration. Thus accurate processing of these data is required to get meaningful results from their analysis. In this paper we aim at developing an automatic processing and photometric calibration method that provides precise and consistent results when applied to historical SHG. The proposed method is based on the assumption that the centre-to-limb variation of the intensity in quiet Sun regions does not vary with time. We tested the accuracy of the proposed method on various sets of synthetic images that mimic problems encountered in historical observations. We also tested our approach on a large sample of images randomly extracted from seven different SHG archives. The tests carried out on the synthetic data show that the maximum relative errors of the method are generally <6.5%, while the average error is <1%, even if rather poor quality observations are considered. In the absence of strong artefacts the method returns images that differ from the ideal ones by <2% in any pixel. The method gives consistent values for both plage and network areas. We also show that our method returns consistent results for images from different SHG archives. Our tests show that the proposed method is more accurate than other methods presented in the literature. Our method can also be applied to process images from photographic archives of solar observations at other wavelengths than Ca II K. △ Less

Submitted 13 October, 2017; originally announced October 2017.

Comments: 24 pages, 24 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 609, A92 (2018)