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

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.


Introduction
The variation in solar irradiance at timescales greater than a day is believed to be dominantly driven by photospheric magnetism (Solanki et al. 2013;Yeo et al. 2017a).Models developed to reproduce solar irradiance variability by relating it to magnetic activity on the solar surface provide the radiative forcing input required by climate simulations (Haigh 2007).Solar irradiance variability is modelled as the sum effect of the intensity deficit from sunspots and the excess from faculae and network, determined from observations of solar magnetism (Domingo et al. 2009;Yeo et al. 2014a).Most of the models aimed at reconstructing solar irradiance variability back to pre-industrial times, a period of particular interest to climate studies, rely on sunspot indices such as the total sunspot area, international sunspot number, and group sunspot number, as these are the only direct observations of solar magnetic features to go this far back in time (e.g.Lean 2000;Krivova et al. 2007Krivova et al. , 2010;;Dasi-Espuig et al. 2014, 2016;Coddington et al. 2016;Wu et al. 2018).Of course, inferring the effect of not just sunspots but also of faculae and network on solar irradiance from sunspot indices requires knowledge of how the amount of faculae and network present relates to sunspots.For example, the solar irradiance reconstruction by Dasi-Espuig et al. (2014, 2016) made use of the model by Cameron et al. (2010), which incorporates the empirical relationship between facular and sunspot area reported by Chapman et al. (1997), to calculate the amount of faculae and network present from the sunspot area and number.
The understanding of how faculae and network relate to sunspots is also relevant to that of the brightness variations of Sun-like stars.The synoptic programmes at the Fairborn, Lowell, and Mount Wilson observatories monitored the brightness and activity of a number of Sun-like stars as indicated by the Strömgren b and y (i.e.visible) photometry and Ca II H&K emission (as a proxy of activity).These observations revealed a dichotomy in the relationship between brightness and activity.While the brightness of less active, older stars rises with increasing activity, it diminishes for more active, younger stars (Lockwood et al. 2007;Hall et al. 2009;Shapiro et al. 2014;Radick et al. 2018).This switch in activity-dependence is interpreted as the transition from a faculae-dominated regime, where the intensity excess from faculae and network has a greater effect on brightness variations than the deficit from starspots, to a spot-dominated regime where the converse is true.The A&A 639, A139 (2020) threshold between the two regimes is estimated to be at log R HK (see definition in Noyes et al. 1984) of between −4.9 and −4.7 (Lockwood et al. 2007;Hall et al. 2009;Shapiro et al. 2014;Radick et al. 2018;Reinhold et al. 2019).The Sun has a mean log R HK of about −4.9 (see Lockwood et al. 2007, and Sect. 4.3) and appears to be faculae-dominated (Shapiro et al. 2016;Radick et al. 2018), suggesting that it lies not far below this threshold.The transition between the faculae-dominated and spot-dominated regimes can therefore be probed by looking at how solar faculae and network relate to sunspots and by extrapolating the apparent relationship to higher activity levels.Chapman et al. (1997) examined the relationship between total sunspot area and Ca II K plage area, taken as a proxy of facular area, over the declining phase of solar cycle 22.This makes use of the fact that chromospheric emission is strongly enhanced in plage and network features overlaying photospheric faculae and network (e.g.Schrijver et al. 1989;Harvey & White 1999;Loukitcheva et al. 2009;Kahil et al. 2017;Barczynski et al. 2018).Chapman et al. (1997) found that facular area conforms to a quadratic relationship with sunspot area.The coefficient of the second-order term is negative, such that the facular-to-sunspotarea ratio decreases as sunspot area increases.Investigations by Foukal (1993Foukal ( , 1996Foukal ( , 1998) ) and Shapiro et al. (2014), extending multiple solar cycles, returned similar results.However, while Chapman et al. (1997) made use of modern, relatively pristine Ca II K spectroheliograms, Foukal (1993Foukal ( , 1996Foukal ( , 1998) ) looked at facular areas based on historical Ca II K spectroheliograms, which suffer from calibration issues and defects (Ermolli et al. 2009), and white-light heliograms, where the intensity contrast of faculae is not only weak, but also diminishes towards the disc centre (Foukal et al. 2004).While Chapman et al. (1997) and Foukal (1993Foukal ( , 1996Foukal ( , 1998) made use of measured facular area, Shapiro et al. (2014) examined the total facular and network disc coverage from Ball et al. (2012), which was determined indirectly from full-disc magnetograms using an empirical relationship between the magnetogram signal and the facular filling factor (Fligge et al. 2000).A proper examination of the relationship between sunspot and facular area over multiple solar cycles is still lacking.The restoration and calibration of the various historical Ca II K spectroheliogram archives, which extend as far back as the beginning of the last century, will facilitate such investigations (Chatzistergos et al. 2018(Chatzistergos et al. , 2019a,b),b).
As noted in the previous paragraph, chromospheric emission is strongly enhanced in plage and network features overlaying photospheric faculae and network.It follows that chromospheric indices such as the 10.7 cm radio flux (F 10.7 , Tapping 2013), Ca II K 1 Å emission index (Bertello et al. 2016), Lyman α irradiance (Woods et al. 2000), and Mg II index (Heath & Schlesinger 1986;Snow et al. 2014) are modulated with plage and chromospheric network emission, and their relation to sunspot indices offers another avenue to probe how faculae and network relate to sunspots.In this article, we present such an effort, examining the relationship between the aforementioned chromospheric indices and the group sunspot number, international sunspot number and total sunspot area.At the same time, we also look at how the total magnetic flux enclosed in faculae and network as apparent in full-disc solar magnetograms (Yeo et al. 2014b), denoted F φ , relates to the various sunspot indices.The aim is to complement and extend the studies on the relationship between sunspot and facular area (i.e.Foukal 1993Foukal , 1996Foukal , 1998;;Chapman et al. 1997;Shapiro et al. 2014) by investigating how chromospheric emission and F φ relate to sunspots.The advantage is that while it is still a challenge to compare sunspot and facular area over multiple solar cycles, the F 10.7 goes back to 1947, and the other chromospheric indices and F φ to the 1970s (see Sect. 2), lending themselves to a multi-cycle comparison to sunspot indices.It is worth pointing out that there are sources of variability in chromospheric and coronal emission other than their enhancement over faculae and network.For example, the enhancement of solar 10.7 cm emission in compact sources associated with sunspots and in coronal loops (Tapping 1987).The F 10.7 and the various chromospheric indices are strongly, but not solely, modulated by faculae and network prevalence.
How faculae and network relate to sunspots is complex in that the amount of faculae and network present at a particular time is not indicated by prevailing sunspots alone.Active regions and their decay products dissipate slower than the sunspots they bear, with the result being that the magnetic flux of a given active region persists, manifesting as faculae and network, even after the embedded sunspots have decayed (van Driel-Gesztelyi & Green 2015).This means, at a given time, there can be faculae and network present that are not associated with the sunspotbearing active regions present at that moment, but with earlier active regions where the sunspots have already dissipated.Also, magnetic flux emerges on the solar surface, not just in active regions, but also in ephemeral regions (Harvey 1993(Harvey , 2000) ) and in the form of the internetwork magnetic field (Livingston & Harvey 1975;Borrero et al. 2017).Ephemeral regions and the internetwork magnetic field contribute to the magnetic network, but as they do not contain sunspots, their prevalence is not captured by monitoring sunspots.The internetwork magnetic field does not appear to vary over the solar cycle (Buehler et al. 2013;Lites et al. 2014), suggesting that its contribution to the magnetic network is invariant over cycle timescales.In contrast, the number of ephemeral regions varies along with the solar cycle (i.e. higher at cycle maxima and lower at minima, Harvey 1993Harvey , 2000)), indicating that their contribution to the magnetic network might similarly exhibit cyclic variability.Preminger & Walton (2005, 2006a,b, 2007), denoted here as PW, sought to reproduce solar irradiance, total photospheric magnetic flux, and various chromospheric and coronal indices, referred to here as the target indices, from total sunspot area by convolving it with the appropriate finite impulse response (FIR) filter.In other words, they modelled the relationship between total sunspot area and the target indices as a linear transformation.For a given target index, the FIR filter is given by the deconvolution of total sunspot area from the target index.As such, the FIR filter encapsulates the response of the target index to sunspots.The FIR filters indicate that the appearance of a sunspot would produce a response in the target indices over multiple rotation periods, and the time variation in this response is consistent with what is expected from the fact that active region magnetic flux persists, in the form of faculae and network, for some time after their sunspots have dissipated (cf.Sect.3.2).The application of the FIR filters to total sunspot area closely replicated the target indices, leading the authors to conclude that the relationship between total sunspot area and the various target indices is well represented by a linear transformation and does not change with time.A modification of the PW model was recently proposed by Dudok de Wit et al. (2018), which we discuss in Sect.3.3.
The more recent investigations by Svalgaard & Hudson (2010), Tapping & Valdés (2011), Livingston et al. (2012), and Tapping & Morgan (2017) reported that the relationship between the F 10.7 and the international sunspot number and total sunspot area appears to have been changing since solar cycle 23.In each of these studies, the authors modelled the F 10.7 implicitly assuming the level at a particular time is a function of A139, page 2 of 13 prevailing sunspots alone.For example, Tapping & Valdés (2011) and Tapping & Morgan (2017), hereinafter referred to collectively as TVM, described the F 10.7 as an exponentialpolynomial function of the sunspot indices, described here in Sect.3.1.Svalgaard & Hudson (2010) also examined the scatter plot of the F 10.7 and the international sunspot number (see Fig. 2 in their paper).This compares the F 10.7 at each time to the sunspot number at that time, which again implicitly assumes that the F 10.7 is a function of prevailing sunspots alone.We had noted that due to the way active regions evolve, the amount of faculae and network present at a given time is indicated not just by the sunspots present at that moment, but also by sunspots in the recent past.This would, of course, extend to the response of the F 10.7 to faculae and network.The analyses of Svalgaard & Hudson (2010), Livingston et al. (2012), and TVM, by treating the F10.7 as a function of prevailing sunspots alone, does not take this into account.Since solar 10.7 cm emission is enhanced not just over sunspots, but also over faculae and network, it is inconclusive if the findings of these authors point to secular variability in the relationship between the F 10.7 and sunspots.
Let us refer to chromospheric indices and F φ , which are both modulated by faculae and network prevalence, collectively as facular indices.We examined the relationship between sunspot and facular indices through an empirical model that extends the linear transformation approach put forward by PW.In the following, we describe the sunspot and facular indices considered (Sect.2) before presenting the model (Sect.3).The realism of any model of the relationship between sunspot and facular indices is, of course, indicated by how well the facular indices can be reconstructed from the sunspot indices with the model.We demonstrate the proposed model to be competent in this regard (Sect.4.1).Making use of the model reconstruction of the facular indices from the sunspot indices and an empirical model of solar irradiance variability based on it, we examine how chromospheric emission and F φ scale with sunspots (Sect.4.2), and what the apparent relationship implies for solar and stellar brightness variations (Sect.4.3).Finally, we provide a summary of the study in Sect. 5.

Data
In this study, we investigate how faculae and network relate to sunspots by examining the relationship between sunspot indices and what we term facular indices (defined in Sect.1), denoted S and F, respectively.Specifically, we compared the daily total sunspot area, S A , international sunspot number, S N and group sunspot number S G to the daily 10.7 cm radio flux, F 10.7 , Ca II K 1 Å emission index, F CaIIK , Lyman α irradiance, F Lα , Mg II index, F MgII and total magnetic flux enclosed in faculae and network, F φ (depicted in Fig. 1).We made use of the Penticton F 10.7 record (Tapping 2013) and the composite time series of S A by Balmaceda et al. (2009), F CaIIK by Bertello et al. (2016), F Lα by Machol et al. (2019) and F MgII provided by IUP (Institut für Umweltphysik, Universität Bremen).The F φ time series is taken from Yeo et al. (2014b), who isolated the faculae and network features in daily full-disc magnetograms dating back to 1974.While the Machol et al. (2019) F Lα composite goes back to 1947, the 1947 to 1977 segment (grey, Fig. 1h) is not provided by F Lα measurements but a model based on the 10.7 cm and 30 cm radio flux.For this reason, we exclude it from further consideration.
Both the S G and S N have been revised recently.Several competing revisions of the Hoyt & Schatten (1998)    From top to bottom: total sunspot area (S A ), international sunspot number (original; S N1 , revision; S N2 ), group sunspot number (original; S G1 , revision by Chatzistergos et al. 2017; S G2 ), 10.7 cm radio flux (F 10.7 ), Ca II K 1 Å emission index (F CaIIK ), Lyman α irradiance (F Lα ), Mg II index (F MgII ), and the total magnetic flux enclosed in faculae and network (F φ ).The S A , F 10.7 , F CaIIK , F Lα and F φ are in units of ppm of the solar hemisphere, solar flux units, Ångström, Wm −2 and Weber, respectively.The grey segment of the F Lα time series is excluded from the analysis.See Sect. 2 for details.
(2016), Cliver &Ling (2016), andChatzistergos et al. (2017).Of these, only the Chatzistergos et al. (2017) revision is suitable for the current study.We are interested in the daily S G , but the Svalgaard & Schatten (2016) revision is only available at monthly cadence due to the calibration method.The Usoskin et al. (2016) revision is moot here since it does not introduce any modifications to the Hoyt & Schatten (1998) time series after 1900, and therefore it compares similarly to the various facular indices as the original time series.As the Cliver & Ling (2016) revision only goes up to 1976, we cannot compare it to the F CaIIK , F Lα and F MgII as none of these go further back than 1976.To find any effect of the changes introduced to the S G on the analysis, we looked at both the Hoyt & Schatten (1998) and Chatzistergos et al. (2017) time series, distinguished here as S G1 and S G2 .For the same reason, we examined both the original and revised S N composites (Clette et al. 2014(Clette et al. , 2016;;Clette & Lefèvre 2016), denoted S N1 and S N2 .
To investigate the implications of the apparent relationship between sunspot and facular indices on solar and stellar brightness variations, we also make use of the PMOD total solar irradiance (TSI) composite (version 42_65_1709, A139, page 3 of 13 A&A 639, A139 (2020) Fröhlich 2000Fröhlich , 2006) ) and the Balmaceda et al. (2009) photometric sunspot index (PSI) composite.The PSI (Hudson et al. 1982;Fröhlich et al. 1994) indicates the proportional deficit in TSI due to sunspots.
We made use of the various data sets as available at the time of study, downloaded on 22 January 2020.The online sources are listed in the acknowledgements.

Models
We examined the relationship between the various sunspot and facular indices (Fig. 1) through an empirical model.The model, abbreviated to YSK, is an extension of the linear transformation approach proposed by PW (Preminger & Walton 2005, 2006a,b, 2007).In Sect.4.1, we examine how well we can replicate the facular indices from the sunspot indices with the YSK model and with the PW and TVM (Tapping & Valdés 2011;Tapping & Morgan 2017) approaches (serving as control).Before that, we first describe the TVM (Sect.3.1), PW (Sect.3.2), and YSK models (Sect.3.3).We denote the TVM model of facular index F as a function of sunspot index S as F * TVM (S ), and similarly that by PW and YSK as F * PW (S ) and F * YSK (S ).

TVM model
Tapping & Valdés (2011) and Tapping & Morgan (2017) examined the relationship between the F 10.7 and the S A and S N .
They smoothed the various time series and fit an exponentialpolynomial relationship of the form where f 1 to f 3 are fit parameters, with the condition that f 1 ≥ 0, and f 4 is fixed at 67 sfu.In the comparison to the S N , f 2 was also fixed at null.Here, we applied the TVM model (Eq.( 1)) to the S and F data sets (Sect.2) as they are (i.e.no smoothing) and without any constraints on f 1 to f 4 , apart from the f 1 ≥ 0 condition.

PW model
Preminger & Walton aimed to reproduce various chromospheric indices, including the four examined here, from the S A .The response of a given chromospheric index to S A is modelled as a linear transformation.Specifically, as the convolution of S A with the finite impulse response (FIR) filter, H PW,emp derived empirically by the deconvolution of S A from the chromospheric index.That is, where g 3 is a fit parameter.Preminger & Walton found H PW,emp to be consistent with what is expected from active region evolution.They described the form of H PW,emp with a model FIR filter, H PW,mod (Preminger & Walton 2007).As a function of time, t: where −0.25t ≤ t ≤ 14.25t and t is the synodic rotation period of the Sun, taken here to be 26.24days.As illustrated in Fig. 2a, this describes a sequence of lobes of diminishing amplitude, as modulated by the exponential envelope of time constant g 2 .The lobes are 0.5t wide and come in intervals of t , with the first lobe centred on t = 0.The sequence of lobes is truncated at t = 14.25t on the observation that there is no discernible signal in H PW,emp above this limit.Active regions emerge rapidly (days) and decay slowly (weeks to months).The response of F to active regions as they emerge is represented by the rising edge of the first lobe.As active regions decay, the effect on F not only diminishes with time but is also modulated by solar rotation, delineated by the sequence of lobes of diminishing amplitude.We note here that H PW,mod , defined such that it is fixed at unity at t = 0, describes the form, but not the amplitude of H PW,emp .
While the empirical FIR filters comply with what is expected from active region evolution, they do not contain any features that can be clearly attributed to ephemeral regions.This is the case even though it is known that the number of ephemeral regions varies roughly in phase with the sunspot cycle (Harvey 1993(Harvey , 2000)), which alludes to a connection between ephemeral regions and sunspots.We surmise that though both sunspots and ephemeral regions are manifestations of the same magnetic cycle, ephemeral region emergence is not coupled to sunspot emergence.As such, the PW model is essentially a model of the effect of active regions on F.
Preminger & Walton compared S A ⊗ H PW,mod , with g 2 fixed at certain arbitrary values, to the F 10.7 and Ap index (Foukal 1996).Otherwise, there was no attempt to model the relationship between chromospheric and sunspot indices with H PW,mod .Generalising the PW model (Eq.( 2)) to other sunspot indices and taking into account that H PW,mod describes the time-dependence of H PW,emp , we can rewrite Eq. ( 2) as where g 1 to g 3 are fit parameters.The g 1 term, missing from Eq. ( 3), scales the amplitude of the model FIR filter to the appropriate level.In our analysis, we applied this form of the PW model to the S and F data sets (Sect.2).

YSK model
The current model is an extension of the PW model (Eq.( 4)).It is given by A139, page 4 of 13 where and h 1 to h 4 are fit parameters.
The model FIR filter here, H YSK (Eq.( 6)) is identical to that in the PW model (Eq.( 3)), except the exponential envelope is given by exp (−|t|/h 3 ) instead of exp (−t/h 3 ).The effect on the model FIR filter is illustrated in Fig. 2b.While the envelope function adopted by PW skews the first lobe towards the negative time domain (black solid line), the proposed envelope function renders it symmetrical at about t = 0 (red dashed line).We introduced this modification on the observation that the empirical FIR filters derived by PW (see, for example, Fig. 1  In another departure from the PW model, the model FIR filter is applied to S h 2 instead of S .By applying the FIR filter to S , the PW model implicitly assumes that, active region evolution aside, F scales linearly with S , which is unlikely to be the case.It is known that chromospheric emission does not scale linearly with photospheric magnetic flux density (e.g.Schrijver et al. 1989;Harvey & White 1999;Loukitcheva et al. 2009;Kahil et al. 2017;Barczynski et al. 2018, see also Sect.4.2) and facular area is a quadratic function of sunspot area (Foukal 1993(Foukal , 1996(Foukal , 1998;;Chapman et al. 1997;Shapiro et al. 2014).This alludes to a nonlinear relationship between plage emission and facular magnetic flux, and between facular magnetic flux and S .We introduced the h 2 parameter to take this into account.
As noted in the introduction, Dudok de Wit et al. ( 2018) presented a modified version of the PW model.In their model, the variables are at 27-day (instead of daily) cadence so as to exclude solar rotation effects from the FIR filter.In addition, the g 3 term (Eq.( 2)) is allowed to vary with time.The authors found that with the coarser time resolution, most of the long-term (annual to decadal) variation in F is captured in g 3 (t) instead of the convolution of S A and the FIR filter.We did not adopt these modifications, as neither excluding solar rotation effects from the FIR filter nor splitting the variability in F into two separate terms is necessary for the purposes of the current study.

Model validation
We modelled the relationship between each facular index, F and each sunspot index, S .Taking each combination of F and S , we fitted the YSK model (Eqs.( 5) and ( 6)) and the TVM (Eq.( 1)) and PW models (Eq.( 4)), (serving as control).The fit parameters of the YSK model are listed in Table 1.
To validate the YSK model, we examined how well it reproduces the facular indices from the sunspot indices as compared to the two control models.For each model and each combination of F and S , we derive the following.We calculate the deviation from unity of the Pearson's correlation coefficient between F and the model reconstruction of F from S , F * (S ), denoted 1 − R 2 .This quantity indicates the variability in F that is not replicated in F * (S ).Taking the F-versus-F * (S ) scatter plot, we derive the ratio of the variance normal to and in the direction of the F = F * (S ) line, denoted as σ 2 ⊥ /σ 2 .The more F * (S ) replicates the variability and the scale of F, the lower the value of σ 2 ⊥ /σ 2 .To reveal how closely the long-term (annual to decadal) trend in F is reproduced in F * (S ), we took the three-year running mean of F and F * (S ), denoted as F 3Y and F * (S ) 3Y , and calculate the normalised residual, given by This is the difference between F 3Y and F * (S ) 3Y , normalised to the change in the former between the 2000 solar cycle maximum and 2008 minimum.Following convention, the epoch of solar cycle extrema is taken from the 13-month moving mean of the monthly S N2 .The normalisation expresses the residual as a proportion of solar cycle variability.The discrepancy between the measured and modelled long-term variability is also encapsulated in the root-mean-square of the normalised residual, abbreviated as RMS NR .We tabulate 1 − R 2 , σ 2 ⊥ /σ 2 and RMS NR in Table 2, and depict the normalised residue of the YSK model in Fig. 3.
For the YSK model, 1 − R 2 ranges from about 0.06 to 0.18 and RMS NR from 0.02 to 0.08 (Table 2), indicating that it reproduces about 82% to 94% of the variability in the various facular indices and their long-term trend to about 2% to 8% of solar cycle variability.In terms of 1 − R 2 , σ 2 ⊥ /σ 2 and RMS NR , the YSK and PW models replicate the facular indices better than the TVM model.The only exceptions are the F MgII & S N1 and F Lα & S N1 analyses, where only the YSK model registered a lower RMS NR than the TVM model.The strength of the YSK and PW models over the TVM model highlights how important it is, in such studies, to account for the fact that active region magnetic fluxes decay slower than sunspots (as similarly argued by Foukal 1998;Preminger & Walton 2007), and the suitability of the linear transformation approach proposed by PW for this purpose.
A139, page 5 of 13 A&A 639, A139 (2020) Table 2.For each combination of F and S , the agreement between F and the reconstruction from the TVM, PW and YSK models, as indicated by 1 − R 2 , σ 2 ⊥ /σ 2 and RMS NR , defined in Sect.4.1.The PW model registered a lower σ 2 ⊥ /σ 2 than the YSK model in the F CaIIK & S N2 and F φ & S G1 analyses, and a lower RMS NR for eight of the 25 combinations of F and S .Otherwise, the YSK model achieved a lower value of 1 − R 2 , σ 2 ⊥ /σ 2 and RMS NR than the PW model.Overall, the YSK model describes the relationship between the sunspot and facular indices better than both the TVM and PW models.
Secular variation in the relationship between F and S , if present, will imprint itself on the normalised residual.This will, however, be obscured by fluctuations in the normalised residual from data uncertainty and model limitations.To count any trend in the normalised residual as corresponding to secular variability in the relationship between F and S with confidence, it has to be apparent in multiple combinations of F and S , and significant compared to RMS NR .For the YSK model, in absolute terms, the normalised residual (Fig. 3) is below 10% (dashed lines) almost everywhere, meaning it is comparable to RMS NR (2% to 8%, Table 2).Looking at the F 10.7 & S A (red, Fig. 3b), F 10.7 & S N1 (green, Fig. 3b), and F Lα & S A analyses (red, Fig. 3f), the normalised residual rose by up to 20% between 2000 and 2005.However, this is not corroborated by how these two facular indices compare to the reconstruction from the other sunspot indices (Figs.3b and f), or by the F CaIIK (Fig. 3d), F MgII (Fig. 3h), and F φ analyses (Fig. 3j).In the case of the F CaIIK & S A analysis (red, Fig. 3d), the normalised residual rose to about 20% between 1980 and 1985, and again around 2005, but as before, this is not corroborated by any of the other combinations of F and S .Within the limits of the current analysis, there is no clear evidence of any secular variation in the relationship between F and S over the past four decades.We remind the reader that Svalgaard & Hudson (2010), Livingston et al. (2012), and TVM found the relationship between the F 10.7 and the S A and S N to have been changing since solar cycle 23.The analysis here, making use of a model demonstrated to be more physical than the TVM model and extended to include more facular and sunspot indices, indicates otherwise.We conclude that the contradictory results from these earlier studies might be an artefact of data uncertainty and model limitations.
The S A is a continuous quantity, at least to the resolution limit of the underlying sunspot area measurements, but the S N and S G , indicating the number of sunspots and sunspot groups, are discrete quantities.In the PW and YSK models, F at a given time is, in effect, given by the weighted sum of S over an extended period.The result is that the discrete nature of the S N and S G is not seen in the model reconstruction of F from them.In other words, this time-averaging suppresses the effect of the quantization noise in S N and S G on the model output.It is not straightforward to isolate and quantify the uncertainty introduced into the YSK model by S N and S G being discrete, but the impact on the current discussion is likely to be minimal.As is evident from Fig. 3, at least in terms of the three-year running mean, the modelling results from S N (green and blue) and S G (purple and cyan) are consistent with that from S A (red).
A139, page 6 of 13 Fig. 3. a: three-year running mean of F 10.7 (black) and the YSK model reconstruction of this facular index from the S A (red), S N1 (green), S N2 (blue), S G1 (purple), and S G2 (cyan).b: residual between the observed and reconstructed time series, normalised to the change in the former between the 2000 solar cycle maximum and 2008 minimum (Eq.( 7)).The dashed lines mark the 10% bound.c-j: The corresponding plots for the F CaIIK , F Lα , F MgII and F φ .

How faculae and network relate to sunspots
As noted in the introduction, the amount of faculae and network present on the solar disc at a given time is not indicated by prevailing sunspots alone due to how active regions evolve and the contribution by ephemeral regions and the internetwork magnetic field.In the YSK model, F * YSK (S ), active region evolution is taken into account by the convolution of the sunspot indices, S with the model FIR filter, H YSK (Eqs.( 5) and ( 6)).So whether in measurements or in the YSK model, a particular level of S does not map to a unique value of F. Nonetheless, we can still gain insight into how faculae and network relate to sunspots by looking at the overall trend in F with S .To this end, we compared the amplitude of the solar cycle in F * YSK (S ) and S .We opted to compare F * YSK (S ) instead of the measured F to S because of the following considerations.Since the S time series go further back in time than the F time series (Fig. 1), and F * YSK (S ) evidently extend as far as S , we can compare F * YSK (S ) and S over longer periods than when comparing F and S .More critically, how the cycle amplitude in F and S compare can be affected by the uncertainty in the decadal trend in the various time series.This uncertainty is irrelevant when comparing F * YSK (S ) and S .We recognise that the comparison of F * YSK (S ) and S is only valid as far as the YSK model is physical, but the robustness of the model, demonstrated in Sect.4.1, renders confidence in this approach.The S A is a more direct measure of sunspot prevalence than the S N and S G , which give sunspots and sunspot groups of different areas the same weighting.With this in mind, the focus here is on how F * YSK (S A ) and S A compare.We use ∆F as an abbreviation of the deviation in F from the 2008 solar cycle minimum level.Here, we define the cycle amplitude of F and S , which we denote as A (F) and A (S ), as the value of the three-year running mean of ∆F and S at cycle maxima.In Fig. 4, we chart A F * YSK (S A ) /A (S A ) against A (S A ) (red), revealing the trend in the F-to-S A cycle amplitude ratio with S A cycle amplitude as indicated by the YSK model.To compare the results from the various facular indices, we normalised the F-to-S A cycle amplitude ratio from each index to the cycle 22 value.The uncertainty in the F-to-S A cycle amplitude ratio, marked in the figure, is propagated from the uncertainty in the long-term trend in F * YSK (S ), RMS NR (Table 2).For the various facular indices, the F-to-S A cycle amplitude ratio decreases with increasing S A cycle amplitude (Fig. 4).The decline is steepest for the F CaIIK (Fig. 4b), followed by the F Lα and F MgII (where the trend with S A cycle amplitude is closely similar, Figs.4c,d), then the F φ (Fig. 4e), and finally the F 10.7 (Fig. 4a).For the F 10.7 , the trend is weak in relation to the uncertainty.
As a check, we repeated the above analysis with F * YSK (S N1 ) and S N1 , that is, we examined A F * YSK (S N1 ) /A (S N1 ) as a function of A (S N1 ).To allow a direct comparison to the F * YSK (S A ) and S A analysis, A (S N1 ) is calculated after rescaling this sunspot index to the scale of the S A using the quadratic polynomial fit to the S A -versus-S N1 scatter plot.This is repeated for F * YSK (S N2 ) and S N2 , F * YSK (S G1 ) and S G1 , and F * YSK (S G2 ) and S G2 .The results are drawn in grey in Fig. 4. The A F * YSK (S ) /A (S )versus-A (S ) profiles from the various sunspot indices lie largely within 1σ of one another, indicating that they are mutually consistent, affirming what we noted with the F * YSK (S A ) and S A analysis (red).Notably, for the F 10.7 (Fig. 4a), while the various profiles are within error of one another, they indicate conflicting trends with S A cycle amplitude.In other words, for this particular facular index, the underlying trend is too weak to be established by the current analysis.
The observation here that the F 10.7 departs from the other facular indices in terms how it compares to the S A is at least partly due to the following.We noted in the introduction that while the various chromospheric indices are strongly modulated by faculae and network prevalence due to the enhancement of A139, page 7 of 13 A&A 639, A139 (2020) chromopheric emission over these photospheric magnetic features, there are other sources of variability.Solar 10.7 cm emission is enhanced in compact sources that are associated with sunspots (Tapping 1987).This is not the case for the Ca II H&K, Lyman α, and Mg II h&k lines, and while sunspots can still be darker or brighter in these lines, depending on height in the solar atmosphere, the effect is much weaker.Clearly, this would have contributed to the divergence between the F 10.7 and the other chromospheric indices and F φ noted here.
Excluding the F 10.7 , the F-to-S A cycle amplitude ratio decreases more steeply with increasing S A cycle amplitude for the chromospheric indices (Figs.4b-d) than for the F φ (Fig. 4e).This means the chromospheric index-to-F φ cycle amplitude ratio also diminishes with rising S A cycle amplitude.We attribute this to how chromospheric emission relates to the photospheric magnetic field.Various studies have noted that at chromospheric passbands, the relationship between the intensity excess of chromospheric features and the underlying photospheric magnetic flux density can be described by a power law with an exponent that is below unity (e.g.Schrijver et al. 1989;Harvey & White 1999;Loukitcheva et al. 2009;Kahil et al. 2017;Barczynski et al. 2018;Chatzistergos et al. 2019b).In fact, part of the purpose of the h 2 term in the YSK model (Eq.( 5)) is to capture this relationship (cf.Sect.3.3).The power-law exponent lying below unity means the ratio of the intensity excess and photospheric magnetic flux density of chromospheric features declines with rising photospheric magnetic flux density.The F CaIIK -to-F φ ratio decreasing with increasing S A cycle amplitude is the extension of this behaviour to the disc-integrated Ca II K emission and magnetic flux.The same argument applies to the F Lα and F MgII .Of course, the more the power-law exponent deviates from unity, the more pronounced this effect is.
The power-law exponent has been reported for the Ca II K (Schrijver et al. 1989;Harvey & White 1999;Loukitcheva et al. 2009;Chatzistergos et al. 2019b) and Mg II k lines (Barczynski et al. 2018).However, since the passband of the Mg II k filtergrams used in the Barczynski et al. (2018) study differs from the spectral sampling of the Mg II h&k doublet in the derivation of the F MgII , the reported exponents are of limited relevance to the F MgII .The situation is similar for the F CaIIK .To the best of our knowledge, the power-law exponent at the Lyman α line has not been reported in the literature.The study by Barczynski et al. (2018), which examined six passbands formed at various heights in the solar atmosphere, noted the following.Going from the upper photosphere to the transition region, the power-law exponent decreases with height, up to the temperature minimum, before it starts to increase with height instead.The Mg II h&k doublet is formed higher in the chromosphere than the Ca II K line (Leenaarts et al. 2013a,b), and the Lyman α line is formed even higher, at the boundary to the transition region (Vernazza et al. 1981).So while the power-law exponent at each line as it would apply to the corresponding chromospheric index is not known, the Barczynski et al. (2018) study and the formation height of the various lines do suggest that it should be lower for the Ca II K line than for the Lyman α and Mg II h&k lines.This is consistent with the F-to-S A cycle amplitude ratio dropping more steeply with rising S A cycle amplitude for the F CaIIK (Fig. 4b) than for the F Lα (Fig. 4c) and F MgII (Fig. 4d).
As stated in the introduction, various studies have found the facular-to-sunspot-area ratio to decrease with increasing sunspot area (Foukal 1993(Foukal , 1996(Foukal , 1998;;Chapman et al. 1997;Shapiro et al. 2014).In this section, we see that like facular area, the variation in chromospheric emission and the total magnetic flux enclosed in faculae and network scales with sunspots in such a way that the ratio to sunspot area also decreases with increasing sunspot area (Figs.4b-e).

Implication for solar and stellar brightness variations
To examine the implications of the apparent relationship between facular and sunspot indices on solar and stellar brightness variations, we make use of the empirical model of this relationship, F * YSK (S ) to reconstruct the variation in TSI.For this purpose, we adopt what is termed the proxy approach (e.g.Lean & Foukal 1988;Chapman et al. 1996Chapman et al. , 2013;;Lean et al. 1997;Coddington et al. 2016;Yeo et al. 2017b).We refer to the TSI excess from faculae and network as facular brightening, and the deficit from sunspots as sunspot darkening, abbreviated as FB and SD, respectively.In the proxy approach, TSI variability is given by the multiple linear regression of a chromospheric index and a sunspot index, acting as proxies of FB and SD, to measured TSI.Here, we employ the reconstruction of F from S A , F * YSK (S A ) as the FB proxy and the PSI as the SD proxy.The PSI indicates the proportional deficit in TSI due to sunspots (Hudson et al. 1982;Fröhlich et al. 1994), and is calculated from the same sunspot area and position measurements as the S A .The TSI model is given by A139, page 8 of 13 Fig. 5. a: PMOD TSI composite (black) and the model based on the multiple linear regression of the reconstruction of the F CaIIK from S A , F * CaIIK,YSK (S A ) and the PSI to this TSI time series (red, Eq. ( 8)).b: difference between measurement and model (grey), and the corresponding three-year running mean (blue).The models based on F * Lα,YSK (S A ) and F * MgII,YSK (S A ) compare similarly to the PMOD TSI composite (Table 3), not shown to avoid repetition.
where k 1 to k 3 are fit parameters.What we did here is incorporate the YSK model into an existing TSI modelling technique.The objective being to examine what the observations we made earlier, based on the YSK model, about the relationship between F and S A (Sect.4.2) might imply for the effect of faculae, network, and sunspots on solar and stellar brightness variations.Using F * YSK (S A ) as the FB proxy also allowed us to reconstruct TSI over a longer period than if we used the significantly shorter measured F time series.
We derived three TSI models taking the reconstruction of the F CaIIK , F Lα and F MgII from S A as the FB proxy in turn.To minimise any bias from data uncertainty, the three models were optimised to the PMOD TSI composite in such a way that the k 2 term, the coefficient of the PSI in Eq. ( 8), is common.In Fig. 5a, we depict the model based on the F CaIIK reconstruction (red) along the PMOD composite (black).The fit parameters and the agreement between modelled and measured TSI, in terms of 1−R 2 , σ 2 ⊥ /σ 2 and RMS NR (defined in Sect.4.1), are summarised in Table 3.Looking at 1 − R 2 and RMS NR , the TSI models replicate about 69% of the variability in the PMOD composite and the long-term trend to about 11% to 12% of solar cycle variability (see also, Fig. 5b).The three models replicate the PMOD composite reasonably well, and almost equally so.For comparison, the TSI models by Dasi-Espuig et al. (2016) and Pelt et al. (2017), also based on sunspot area and position measurements, reproduced about 76% and 69% of the variability in the PMOD composite, respectively.
We did not model TSI with the F 10.7 for the following reasons.The proxy approach requires a quantity that represents the effect of faculae and network on TSI.The enhancement of solar 10.7 cm emission over sunspots makes the F 10.7 less suited for this purpose than the F CaIIK , F Lα and F MgII , where the intensity deficit or excess over sunspots provides a much smaller contribution to total emission (cf.Sect.4.2).While the main contribution to TSI is from the photosphere, the Ca II H&K and Mg II h&k doublets are formed in the chromosphere, the Lyman α line at the boundary to the transition region, and solar 10.7 cm emis-sion in the corona.The physical processes underlying how photospheric magnetism modulates the temperature structure and brightness of the solar atmosphere changes as we go up in height, such that the various indices would vary differently with faculae and network as compared to TSI.For example, in the lower photosphere, magnetic flux tubes are heated through their side walls by radiative heating, and in the upper photosphere and chromosphere through mechanical and Ohmic dissipations.The consequence of this is that the intensity contrast of faculae and network exhibits rather different disc centre-to-limb variation in the two atmospheric regimes (Yeo et al. 2013;Yeo & Krivova 2019).The corona, in terms of the structure of the magnetic field and heating processes, differs markedly from the photosphere and chromopshere.We expect the response of the F10.7 to faculae and network to depart from that of TSI more than the other chromospheric indices.Indeed, tests indicate the TSI reconstruction based on the F 10.7 to be poorer at reproducing the PMOD composite than the models based on the F CaIIK , F Lα and F MgII .
The proxy approach has its shortcomings, such as the limitations of using chromospheric indices to represent the effects of faculae and network on solar irradiance, as we just discussed (see also Yeo et al. 2014a;Yeo & Krivova 2019).While more physical and sophisticated techniques to model TSI exist, such as the 3D magnetohydrodynamics simulation-based model presented by Yeo et al. (2017a), there is no obvious way to incorporate the YSK model into such TSI modelling approaches.It is straightforward, however, to incorporate the YSK model into the proxy approach.And as noted earlier, the resulting TSI reconstruction reproduces observed TSI variability reasonably closely, and just as well as existing models based on similar data.This is sufficient for our objective, which is to examine what the observations we made with the YSK model in Sect.4.2 might imply for solar and stellar brightness variations through a TSI model that incorporates the YSK model.
For the same reason that a particular level of S does not correspond to a unique value of F, discussed in Sect.4.2, S does not uniquely map to TSI.Just as we did in that section to elucidate the overall trend in F with S , we also examined how FB, SD, the ratio of the two, and TSI vary with S A by looking at the solar cycle amplitude.
As with ∆F, ∆TSI represents the deviation in TSI from the 2008 solar cycle minimum level.We define the cycle amplitude of FB, SD, and TSI, denoted by A (FB), A (SD) and A (TSI), as the value of the three-year running mean of k 1 ∆F * YSK (S A ), k 2 PSI and ∆TSI at cycle maxima.In Fig. 6, we chart A (FB), |A (SD)|, A (FB) / |A (SD)| , and A (TSI) as a function of A (S A ).The mean of the values from the three TSI models is depicted (dots).The uncertainty is propagated from the RMS NR between F * YSK (S A ) and F (Table 2), and the formal regression error of k 1 and k 2 .
We found A (FB) to scale quadratically with A (S A ) (Fig. 6a) and A (SD) to scale linearly with A (S A ) (Fig. 6b) such that the ratio, A (FB) / |A (SD)| , decreases with A (S A ) (Fig. 6c).Consequently, A(∆TSI) increases with A (S A ) at a diminishing rate (Fig. 6d).This means, similar to what was noted for the facularto-sunspot-area ratio (Foukal 1993(Foukal , 1996(Foukal , 1998;;Chapman et al. 1997;Shapiro et al. 2014) and the ratio of the variation in facular indices to sunspot area (Sect.4.2), the ratio of facular and network brightening to sunspot darkening decreases as sunspot area increases.As a result, TSI rises more slowly with sunspot area than if the two quantities were linearly related to one another.This is why even though solar cycle 23 is markedly weaker than cycle 22, TSI rose to comparable levels over both cycles (de Toma et al. 2001).Comparing solar cycle 23 to cycle 22 (marked in Fig. 6d), the A (S A ) ratio is 0.70 and the A (TSI) ratio is 0.89.A139, page 9 of 13 A&A 639, A139 (2020) Table 3.The fit parameters of the empirical TSI models (k 1 to k 3 , Eq. ( 8)) derived taking the reconstruction of the F CaIIK , F Lα and F MgII from S A , denoted F * CaIIK,YSK (S A ), F * Lα,YSK (S A ) and F * MgII,YSK (S A ), as the FB proxy (i.e. the proxy of the TSI excess from faculae and network).
FB proxy Notes.The various models are optimised to the PMOD TSI composite in such a way that k 2 is identical.The agreement between model and measurement in terms of 1 − R 2 , σ 2 ⊥ /σ 2 and RMS NR , defined in Sect.4.1, is also tabulated.Due to the non-linear relationship between TSI and sunspot area, while cycle 23 is, in terms of S A , 30% weaker than cycle 22, the TSI cycle amplitude is only 11% weaker.We remind the reader that for Sun-like stars, going above a certain level of activity, they switch from growing brighter with rising activity to becoming dimmer instead, which is interpreted as the transition from a faculae-dominated to a spot-dominated regime (Lockwood et al. 2007;Hall et al. 2009;Shapiro et al. 2014;Radick et al. 2018).The Sun appears to be just below this transition, opening up the possibility of studying this phenomena by looking at how solar faculae and network relate to sunspots, and extrapolating the apparent relationship to higher activity levels.We attempt exactly that here by projecting the TSI model, which is based on the empirical model of the relationship between facular and sunspot indices, to higher activity levels.
We had examined A (FB), |A (SD)|, A (FB) / |A (SD)| , and A (TSI) as a function of activity as indicated by A (S A ) (Fig. 6).In the stellar studies, activity is characterised by the log R HK .Accordingly, we converted F * CaIIK,YSK (S A ) to the S -index scale (not to be confused with sunspot indices, abbreviated as S in this article) using the conversion relationship reported by Egeland et al. (2017), and then the result to log R HK following the procedure of Noyes et al. (1984).In this computation, we assumed a solar (B − V) colour index of 0.653 (Ramírez et al. 2012).In Fig. 7, we plot A (FB) (red dots), |A (SD)| (blue dots) and A (TSI) (black dots) again, this time as a function of log R HK .The log R HK level corresponding to each A (FB), A139, page 10 of 13 Fig. 8. a: log R HK over a simulated solar cycle.The dashed and dotted lines, as in Fig. 7, mark log R HK at T1 and T2.The scenario where the activity peak is below T1 (red), at T1 (green), between T1 and T2 (blue), at T2 (purple) and above T2 (cyan) are depicted.b: TSI corresponding to each scenario, similarly colour-coded, as given by the empirical relationship between A (TSI) and log R HK (black curve, Fig. 7).The difference to the 2008 minimum TSI level (solid black line) is drawn.
|A (SD)|, and A (TSI) value is given by the value of the threeyear running mean of log R HK at that solar cycle maximum.The data points sit just above log R HK of −4.9.The switch in the activity-dependence of the brightness of Sun-like stars is estimated to be between log R HK of −4.9 and −4.7 (Lockwood et al. 2007;Hall et al. 2009;Shapiro et al. 2014;Radick et al. 2018;Reinhold et al. 2019).To probe what happens to facular brightening, sunspot darkening, and TSI here, we fitted a straight line to A (FB) (red line) and a quadratic polynomial to |A (SD)| (blue curve), and extrapolated both fits to log R HK of −4.65.Then, we calculated the TSI level corresponding to the A (FB) line and the A (SD) curve, represented by the black curve.
While facular brightening scales linearly with log R HK (red line, Fig. 7), sunspot darkening (blue curve) rises increasingly rapidly with activity.At log R HK of about −4.80, sunspot darkening rises more rapidly with activity than facular brightening, with the consequence that TSI (black curve) goes from increasing with activity to decreasing instead (dashed lines).To aid this discussion, let us abbreviate this transition in TSI activitydependence as T1.Eventually, the negative correlation between TSI and activity results in TSI dropping below the 2008 minimum level at log R HK of about −4.68 (dotted lines), denoted as T2 (recall, A (TSI) is defined such that null corresponds to the 2008 minimum level).What happens to TSI at T1 and T2 is further illustrated in the following.Supposing log R HK varies over a solar cycle as pictured in Fig. 8a, using the A (TSI) extrapolation (black curve, Fig. 7), we computed the corresponding TSI (Fig. 8b) in the scenario that activity peaks below T1 (red), at T1 (green), between T1 and T2 (blue), at T2 (purple), and above T2 (cyan).When a given cycle is sufficiently strong, such that activity peaks above T1, TSI starts to dip around the cycle maximum, such that instead of varying along with the activity cycle, TSI exhibits a double-peaked form.It is worth emphasising that we based this analysis on the three-year running mean of the various quantities so as to elucidate the overall trend with activity.This discussion does not apply to the variability at shorter timescales.
The observation here that T1 occurs around the range of activity where the activity-dependence of the brightness of Sunlike stars exhibits a similar switch bolsters the interpretation of what we see in Sun-like stars as the transition from a faculaedominated to a spot-dominated regime.Since the TSI model is based on the empirical model of the relationship between fac-ular and sunspot indices, this result also implies that the latter accurately captures how solar faculae and network relate to sunspots.Converted to log R HK , the F CaIIK composite, which extends four solar cycles, has a mean and standard deviation of −4.90 and 0.03, respectively.This means T1, at log R HK of about −4.80, is just over three standard deviations away from the mean observed solar level.This underlines how delicately balanced facular brightening and sunspot darkening are on the Sun (cf.Shapiro et al. 2016).
Of course, any extrapolation of empirical relationships needs to be interpreted with care.In this case, the TSI model and the underlying model of the relationship between facular and sunspot indices closely replicate a range of observations (see Sect. 4.1, Fig. 5, and Table 3), and the turning point of the A (TSI) extrapolation (i.e.T1: dashed lines, Fig. 7) is close to the data points, where the extrapolation is the most reliable.As is evident from the figure, the 95% confidence interval of the A (TSI) extrapolation (shaded region) remains relatively small around T1.This gives us confidence that the extrapolation of the TSI model is, at least up to T1, somewhat reasonable.
The stellar studies are based on the Strömgren b and y photometry of Sun-like stars, not the bolometric photometry.However, as there are no direct observations of the Sun in Strömgren b and y, deriving an empirical model of solar Strömgren b and y as we did for TSI is currently not possible.This is a limitation of the current model.Studies have reported linear relationships for converting TSI to solar Strömgren b and y (see Radick et al. 2018 and references therein), but these are very approximate.Additionaly, since these relationships are linear, repeating the above analysis on the results of having applied them to the PMOD TSI composite would have no qualitative effect on the results.It would be possible to model Strömgren photometry with a proper model of spectral solar irradiance (e.g.Shapiro et al. 2016), but is outside the scope of the current study.This does not detract however, from the fact that the TSI model indicates the bolometric photometry of the Sun to behave in an analogous manner to the Strömgren b and y photometry of Sunlike stars.

Summary
How faculae and network relate to sunspots is of interest for the implications for solar and stellar brightness variations.In this study, we probed this relationship by looking at how chromospheric indices and the total magnetic flux enclosed in faculae and network, termed faculae indices, compare to sunspot indices.This makes use of the fact that chromospheric emission is strongly enhanced in plage and network features overlaying photospheric faculae and network, such that chromospheric indices are, to a greater or lesser degree, modulated with the amount of faculae and network present.Specifically, we compared the 10.7 cm radio flux, Ca II K 1 Å emission index, Lyman α irradiance, Mg II index and total facular and network magnetic flux to the total sunspot area, international sunspot number and group sunspot number.
We presented an empirical model of the relationship between facular and sunspot indices (Sect.3.3).The model is a modification of the Preminger & Walton (2005, 2006a,b, 2007) model, which describes the variation in chromospheric indices from active region evolution.The relationship between a given chromospheric index and sunspot index is described as the convolution of the latter with a suitable finite impulse response (FIR) filter.The modifications introduced here were aimed at more accurately reflecting what Preminger & Walton (2006a) noted of A139, page 11 of 13 A&A 639, A139 (2020) the FIR filters they had derived empirically from observations, and also to take into account the likely non-linear relationship between plage emission and facular magnetic flux, and between the latter and sunspot prevalence.Taking the current model, we reconstructed the facular indices from the sunspot indices.The model not only replicated most of the observed variability (up to 94%), but it is also better at doing so than the Preminger & Walton (2005, 2006a,b, 2007) model, as well as that of Tapping & Valdés (2011) and Tapping & Morgan (2017) (Sect. 4.1).Tapping & Valdés (2011) and Tapping & Morgan (2017), along with Svalgaard & Hudson (2010) and Livingston et al. (2012), found the relationship between the 10.7 cm radio flux and the total sunspot area and international sunspot number to have changed since solar cycle 23.The cited studies, in their analyses, had implicitly assumed that the 10.7 cm radio flux at a given time is a function of the sunspot area/number at that time alone (see Sects. 1 and 3.1).The magnetic flux in active regions persists, manifest as faculae and network, for some time after the sunspots they bear have decayed.As a consequence, the amount of faculae and network present at a given time, and therefore the response of the F10.7 to these magnetic structures, is indicated not just by prevailing sunspots, but also sunspots that had emerged in the recent past.Contrary to these studies, we found no clear evidence of any secular variation in the relationship between the facular and sunspot indices examined over the past four decades (Sect.4.1).The present analysis made use of a model that takes what we just noted about active region evolution into account and demonstrated to be more physical than the Tapping & Valdés (2011) and Tapping & Morgan (2017) model, and is extended to more facular and sunspot indices.Taking this into consideration, the conflicting results from the earlier studies is likely an artefact of data uncertainty and limitations in their analyses and models.
Various studies have noted that the facular-to-sunspot-area ratio diminishes with sunspot area (Foukal 1993(Foukal , 1996(Foukal , 1998;;Chapman et al. 1997;Shapiro et al. 2014).Comparing the reconstruction of the facular indices from the sunspot indices to the latter, we found the ratio of the variation in chromospheric emission and total facular and network magnetic flux to sunspot area to exhibit the same behaviour, decreasing with increasing sunspot area (Sect.4.2).
Making use of the fact that chromospheric indices are reasonable proxies of the effect of faculae and network on solar irradiance, we examined the implications for solar and stellar brightness variations by means of an empirical model of total solar irradiance (TSI).The TSI model, which incorporates our model of the relationship between chromospheric indices and total sunspot area, indicates that the ratio of the TSI excess from faculae and network and the deficit from sunspots also decreases with increasing sunspot area.The consequence is that TSI rises with sunspot area more slowly than if the two quantities are linearly related to one another (Sect.4.3).TSI rising with activity at a diminishing rate explains why TSI rose to comparable levels over solar cycles 22 and 23, even though cycle 23 is significantly weaker than cycle 22.The current study extended and improved upon what was noted earlier of the facular-to-sunspot-area ratio through an examination of independent data sets over multiple cycles (as noted in the introduction, a proper study of the relationship between facular and sunspot area over multiple solar cycles is still lacking).
We extrapolated the trend in facular and network brightening, sunspot darkening, and TSI with activity to higher activity levels, over the range where Sun-like stars are observed to switch from growing brighter with rising activity to becoming dimmer instead.The projection indicates that in this activity range, sunspot darkening will gradually rise faster with activity than faculae and network brightening, such that the activitydependence of TSI exhibits a similar switch.This bolsters the interpretation of the dichotomy in the relationship between the brightness and activity of Sun-like stars as the transition from a faculae-dominated to a spot-dominated regime.This result also suggests that our model of the relationship between facular and sunspot indices, which underlies the TSI model, accurately captures how solar faculae and network relate to sunspots. h

Fig. 1 .
Fig.1.Sunspot and facular indices examined in this study.From top to bottom: total sunspot area (S A ), international sunspot number (original; S N1 , revision; S N2 ), group sunspot number (original; S G1 , revision byChatzistergos et al. 2017; S G2 ), 10.7 cm radio flux (F 10.7 ), Ca II K 1 Å emission index (F CaIIK ), Lyman α irradiance (F Lα ), Mg II index (F MgII ), and the total magnetic flux enclosed in faculae and network (F φ ).The S A , F 10.7 , F CaIIK , F Lα and F φ are in units of ppm of the solar hemisphere, solar flux units, Ångström, Wm −2 and Weber, respectively.The grey segment of the F Lα time series is excluded from the analysis.See Sect. 2 for details.

Fig. 2 .
Fig. 2. The model FIR filter in the PW model, H PW,mod (black solid lines) and in the current model, H YSK (red dashed lines) for time constants of (a) 60 days and (b) 10 days.The black and red dotted lines follow the exponential envelope of the respective model FIR filters.Right panel: we limit the plot range to −10 ≤ t ≤ 10 days to highlight the difference in the two model FIR filters from the different envelope functions assumed.See Sect. 3 for details.
in Preminger & Walton 2006a) do not indicate any clear skewness in the first lobe at about t = 0.

Fig. 4 .
Fig. 4. Red: for each facular index, the F-to-S A solar cycle amplitude ratio, relative to cycle 22, as a function of S A cycle amplitude.The dots represent the individual cycle values, determined using the model reconstruction of F from S A .The error bars and curves denote the corresponding 1σ uncertainty and quadratic polynomial fit.Grey: the F-to-S cycle amplitude ratio-versus-S cycle amplitude profiles for the S N1 , S N2 , S G1 and S G2 .The various sunspot indices are rescaled to the S A scale to allow a direct comparison.Only the values from cycles 12 to 22, where the five sunspot index time series overlap, are depicted.The horizontal line marks the cycle 22 level (i.e.unity).See Sect.4.2 for details.

Fig. 6 .
Fig.6.From top to bottom: as a function of S A cycle amplitude, the cycle amplitude of (a) the effect of faculae and network (red), and (b) of sunspots on TSI (blue), (c) the ratio of the two (green), and (d) TSI (black).That is, A (FB), |A (SD)|, A (FB) / |A (SD)|, and A (TSI)-versus-A (S A ).The plot points and error bars represent the mean of the values from the three TSI models and the associated 1σ uncertainty (see Sect. 4.3).The latter is omitted in the case of A (FB) and A (SD), where it is so minute as to be obscured by the plot points.The values from the individual TSI models, lying within 1σ of the mean, are not drawn to avoid cluttering.The S A and PSI time series on which the TSI models are based, and therefore the plot points, cover solar cycles 12 to 24.The A (TSI) plot points corresponding to solar cycles 22 and 23 are labelled.The lines correspond to the linear or quadratic polynomial fit.

Fig. 7 .
Fig. 7. Cycle amplitude of FB (red plot points), SD (blue) and TSI (black), taken from Fig.6, as a function of log R HK .The corresponding 1σ uncertainty, not drawn, is generally much smaller than the plot symbols.The red line follows the linear fit to A (FB) and the blue line the quadratic to |A (SD)|, while the black line indicates the corresponding A (TSI) level.The shaded region encloses the 95% confidence interval of the A (TSI) curve.The dashed lines mark the turning point of the A (TSI) curve and the dotted lines where it goes below zero, denoted T1 and T2, respectively.See Sect.4.3 for the physical interpretation.The boxed area is blown up in the bottom panel.

Table 1 .
For each combination of F and S , the fit parameters of the YSK model (h 1 to h 4 , Eqs. (5) and (6)).