Open Access
Issue
A&A
Volume 671, March 2023
Article Number A130
Number of page(s) 16
Section Catalogs and data
DOI https://doi.org/10.1051/0004-6361/202244811
Published online 17 March 2023

© The Authors 2023

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

Observations of solar limb darkening - or, more generally speaking, the center-to-limb variation (CLV) of continuum radiation and spectral lines - have been an important tool for probing the solar atmosphere since Schwarzschild (1906) used this phenomenon to show that the photosphere must be in radiative equilibrium (e.g., Hale & Adams 1907; Voigt 1950; Pagel 1955; Balthasar 1985; Alissandrakis et al. 2017). Perhaps one of the most used limb darkening studies is the one done by Neckel & Labs (1994), where the continuum limb darkening was measured for wavelengths between 330 and 1250 nm and presented in the form of fifth-degree polynomials. There are also several other such atlases that focus on various center-to-limb positions, for example the third solar spectrum (SS3)1 by Sten-flo (2015) and Ramelli et al. (2017); the CLV line studies by David (1961), White (1962), White & Suemoto (1968), Gur-tovenko et al. (1975), Balthasar (1988), Ding & Fang (1991), Faurobert et al. (2013), and Takeda & UeNo (2019), for example; or the atlas by Fathivavsari et al. (2014).

Center-to-limb variation measurements provide a good test for the accuracy of models of atmospheric structure and spectral line formation. For example, in Unsold (1955) the Milne-Eddington model is compared to observations of lines that are dominated by absorption and scattering processes, and similar checks were done more recently by Vernazza et al. (1976), Bj∅rgen et al. (2018), and Alsina Ballester et al. (2021), for example. Other studies have used such measurements to constrain the elemental abundance of spectral lines (e.g., Moore et al. 1966; Pereira et al. 2009; Bergemann et al. 2021). In these cases, the full-width at half-maximum (FWHM) and equivalent width are often reported, or synthesized spectra with different abundance values are compared to the observations. Additionally, studies of the Lyman-α and Mg II h&k lines have been performed by Bonnet et al. (1978) and Lemaire et al. (1981), for example, and more recently by Gunár et al. (2020, 2021), and Rachmeler et al. (2022) to create quiet-Sun reference profiles which could serve as an incident radiation boundary condition for the radiative transfer modeling of chromospheric and coronal structures. Additionally, CLV measurements of the polarization of these lines have been studied by Rachmeler et al. (2022), full disk integrated atlases were presented by Molaro & Monai (2012), and the formation height of several chromospheric lines was inferred by means of CLVs by Wittmann (1976).

Typically, CLV observations are recorded in one of three ways.

  1. One can place a slit or imager beyond the western limb of the Sun and turn off the telescope tracking. This causes the Sun to drift past the field of view (FOV) at a constant speed, and the resulting signal is then recorded (e.g., Wittmann 1976; Neckel & Labs 1994; Löfdahl 2016). This method requires many full disk scans in order to average out the small-scale structure across the solar diameter, as well as the absence of any larger structures that could contaminate the average. These kinds of measurements are typically done from west to east, but they can be achieved in other directions by modifying the tracking speed. This method is restricted to monochromatic measurements or spectrographs that give simultaneous data at several wavelengths.

  2. One can observe at several distinct disk positions (e.g., Ding & Fang 1991; Pierce & Slaughter 1982; Pereira et al. 2009; Lind et al. 2017). This method is easier to work with because the observer can choose the location without activity on the Sun for each μ (the cosine of the heliocentric angle). Additionally, the fact that tracking is active allows for each observation to stay pointed at one location for a longer time to form temporal averages. This comes at the cost of precision, as there is typically no overlap between consecutive pointings.

  3. One can create mosaics of monochromatic images (or sets of monochromatic images) where large amounts of pixels can be averaged, potentially excluding areas that contain active regions (Gunár et al. 2021). This method allows the user to choose between longer integration times per position at the expense of time or vice versa. This is the approach taken here.

In this study, we present a quiet-Sun CLV analysis of the Ca II H&K, Hα, Hβ, Mg I 5173 Å, and Ca II 8542 Å lines. These are strong, broad spectral lines with their cores having been formed in between the temperature minimum and the upper chromosphere. In addition to that, we also study the C I 5380 Å, Fe 16173 Å, Fe I 6301 Å, and 0 I 7772 Å lines which are mainly formed in the photosphere. We also observe the Fe I 5379 Å, Fe I 5172 Å, Ti II 5381 Å, and Ti I 5174 Å lines around the C I 5380 Å and Mg I 5173 Å lines, but we do not study them.

We carried out these observations with a series of semi-overlapping mosaics acquired along a solar radius. The main goal of this study is to present representative profiles of the lines mentioned above, as well as to test their robustness by comparing the spectral parameters extracted from these lines against the ones available in the literature. We expect that these profiles will be useful for abundance studies as well as reference profiles, and as approximate background profiles in the study of optically thin chromospheric structures as was done using a cloud model in Pietrow et al. (2022). Additionally, we aim to use our limb darkening curves to improve the intensity calibration of the (SSTRED; de la Cruz Rodríguez et al. 2015; Löfdahl et al. 2021) data-processing pipeline of the Swedish 1 m-Solar Telescope (SST Scharmer et al. 2003a,b).

2 Observations and data processing

Our data consist of ten mosaics, each spanning one solar radius, taken with the SST, using both the CRisp Imaging SpectroPolarimeter (CRISP, Scharmer et al. 2008) and the CHROMospheric Imaging Spectrometer (CHROMIS; Scharmer 2017) simultaneously. All mosaics were observed during two campaigns, one in June 2019 and another in June 2021, which henceforth are called Campaign I and Campaign II, respectively. Observations from Campaign I were taken from the disk center toward the solar north pole, while observations from Campaign II were taken from the solar south pole toward the disk center (see Fig. 1). In all cases, there was roughly a 30% area overlap between each consecutive pointing (see Fig. 2

To maximize the number of quiet-Sun pixels in each FOV, all observations were carried out during periods when there were no active regions along the mosaic path. The observations were carried out under variable seeing conditions. This is acceptable because we are primarily interested in taking large-scale averages along the mosaic and using the fine-scale structure solely for alignment purposes. The spectral lines from Campaign I were observed with one spectral scan per pointing during better seeing, while the lines from Campaign II were taken with three spectral scans per pointing in order to increase the likelihood of a good scan during poor seeing conditions. All the spectral lines were observed without polarimetry, except for Fe I 6301Å.

We list each of the observed lines in Tables 1 and 2. All wavelength positions are given relative to their line center. We note that not all lines are symmetrically sampled; in such a case, a value only has a + or − instead of a ± symbol in front of it.

The SSTRED pipeline (de la Cruz Rodríguez et al. 2015; Löfdahl et al. 2021) has been designed to process the data from the SST. It not only includes dark, flat field, and polarimetric calibrations, but it also performs image restoration, removing optical aberrations caused by turbulence in the atmosphere (and partially corrected for by the SST adaptive optics) using Multi-Object Multi-Frame Blind Deconvolution (MOMFBD; Löfdahl 2002; van Noort et al. 2005). The MOMFBD algorithm assumes that the point spread function (PSF) is constant within subfields of 5″ squared. The restored subfields are then mosaicked, resulting in a restored version of the entire FOV. MOMFBD methods work best when the contrast is high, the noise is low, and the exposure time is short. However, under poor seeing conditions (Fried parameter r0 of 5 cm or lower), the reconstruction can fail. This is characterized by the subfields becoming apparent as a checkered pattern of squares in the images. In this situation, the intensity may not be conserved, and thus the resulting spectra are expected to become affected and are unusable for the current purpose.

To overcome this problem under poor seeing conditions, we used a modified version of the SSTRED pipeline that performs all the data processing but does not apply MOMFBD. We believe that this can also be useful for time series with seeing dips as they allow for the temporal continuity of the data to be maintained, for data where spatial resolution is not relevant, and as a quick look tool. This feature has now been implemented into the SSTRED pipeline by a keyword that allows the user to set the amount of MOMFBD iterations to 0, thus only skipping the reconstruction while keeping the rest of the process.

The SSTRED pipeline makes use of the flat fields to produce a so-called cavity map that describes the variations of the central wavelength of the CRISP and CHROMIS instrumental profile over the FOV. In order to avoid unnecessary interpolation in the data, SSTRED does not apply the cavity-map corrections by default, but it leaves it to the end-user to apply it – which is typically done by correcting derived Doppler shifts. Since the large-scale variations gave a strong imprint on the derived CLV curves, even after the spatial averaging, we applied the cavity-map correction by interpolating the intensities for each pixel between different wavelength frames. In the case of spectroscopic inversions, the cavity map is often subtracted from the resulting velocities in order to correct for the line shifts without changing the profiles beforehand (e.g., Díaz Baso et al. 2019; Pietrow et al. 2020).

thumbnail Fig. 1

Screen capture of the primary image guider of the SST showing the HMI magnetogram with 25 pointings in a mosaic that spans from the solar south pole to the disk center, taken on 12 June 2021. The red marker represents the location at which the telescope was aimed at the moment that the image was taken.

thumbnail Fig. 2

Overview of mosaics taken during Campaign I, starting from the solar limb on the left and going to disk center on the right. The mosaics span roughly 1000 × 80″. Top: continuum near the blue wing of Fe I 6301 Å. Bottom: line center of Ca II 8542 Å. It can be seen that the Sun seems to extend past μ = 0 in the chromosphere.

Table 1

Summary of all data used in this paper.

3 Methods

In this section, we describe the methods used for the alignment of the mosaics, intensity, and position calibration of the data and the removal of p modes.

3.1 Mosaic alignment

The SST pointing and tracking system has features to simplify the making of mosaics. These features have allowed studies such as Scullion et al. (2014) where a large active region was observed in a 9 × 6 pointing mosaic that spanned a FOV of 280″ × 180″. Pastor Yabar et al. (2020) performed a partially overlapping mosaic with 17 pointings of the north pole region, and Vissers et al. (2021) presented a smaller 2 × 2 mosaic of a flaring region that spanned a FOV of roughly 110 × 110″. Our mosaics, which look at the quiet Sun, consist of 1 × 25 pointings and span a FOV of roughly 60 × 1000″.

Typically, mosaics of this type are aligned by using a cross-correlation algorithm (e.g., Cauzzi & Reardon 2012; Hammerschlag et al. 2013). However, these algorithms require several distinguishing features to operate, which for quiet-Sun observations is not the case, especially very close to the limb where the image contrast decreases. Therefore, the images were aligned by manually looking for matching patterns in the granulation and bright points in the overlapping areas between pointings, as well as with SDO/HMI continuum images (Lemen et al. 2012) which were spatially enhanced (Díaz Baso & Asensio Ramos 2018). Afterward, the alignment was refined by means of an ISPy mosaic alignment tool created for this project (Díaz Baso et al. 2021a). This tool allows the user to plot two consecutive mosaic pointings at an approximate location, which are moved, blinked, and further aligned.

A perfect match is not possible since our images are not destreched and taken under suboptimal seeing conditions. For this reason, we assume a conservative estimate for the alignment error by assuming that the pointing is accurate to within 16 pixels, which corresponds to roughly 1 and 0.5″ for CRISP and CHROMIS, respectively. Finally, we obtained our absolute solar coordinates with respect to the disk center from the alignment of our data with the enhanced SDO/HMI continuum observations.

Table 2

FPI tuning positions for each observed line.

3.2 Intensity calibration

The integration time is constant, which means that apart from the varying extinction, there is a linear relation between detector counts and solar intensity. We corrected our data for the atmospheric extinction variations, which are caused by the change in the Sun’s altitude during the acquisition of the mosaic, by regularly taking a disk center observation in between our observations. In Campaign I we took one disk center observation before starting the mosaic and one at the end. In Campaign II we added a third disk center observation halfway through the mosaic. A quadratic function was then fitted to these observations to obtain an approximation for the extinction at the time of the other pointings. This method has also been integrated into the SSTRED pipeline, while before it was only correcting for a single extinction value. After correcting for the extinction, we calibrated the intensity of all pointings with respect to the disk center pointing by simultaneously fitting several wing points of the averaged disk center profile to the solar atlas by Neckel & Labs (1984). We found that this works well and that no additional calibration between frames are necessary in order to assure a smooth transition between each overlapping frame. We do note that this method cannot be assumed to work with rapidly varying extinction, caused by clouds or calima for example.

3.3 Position calibration

For each mosaic, we used the Pythagorean theorem to calculate the distance of each pixel in the mosaic from the disk center, which can be used to collapse the mosaic into a 1D intensity array. We defined the limb in the same way as Neckel & Labs (1994), which is the inflection point of the sharp tail end of the intensity array in the continuum. This was done by taking the numerical derivative of a single-pixel column from the disk center to the limb, which allowed us to accurately pinpoint the inflection point, which we defined as μ = 0, to within 3 pixels (about 0.2″) for CRISP). Pixels beyond this point were discarded.

The continuum wavelengths used for each spectral line are given in Table 1. The broad lines do not have access to a continuum within their prefilter range, and therefore a compromise was needed. The largest wavelength difference between a line and the continuum wavelength used to define the limb position is for Hα, which uses the CHROMIS continuum filter at 4001 Å.

The wavelength difference between the continuum-defined wavelength and the chromospheric wavelength that it is applied to in itself is not expected to give a large error. Neckel & Labs (1994) reported that the position of the limb varies by 0.12″ between 3010 and 10990 Å, and up to 0.5″ for bad seeing at 400 nm (being smaller for longer wavelengths). In addition to that, the atmospheric differential refraction seemingly compresses the Sun along one axis. However, we estimate that the difference in apparent elevation of the solar limb in our observations that are caused by this effect should be less than 1″ in the extreme case of Ηα and 4001 Å. With this in mind, and considering seeing-induced smearing, we chose a limb distance of 3″ (corresponding to μ = 0.08) as a limit beyond which our results for the broad lines should be taken with caution.

The commonly cited CRISP and CHROMIS image scales of 0.058″ pixel–1 and 0.0375″ pixel–1, respectively, were derived from measurements by Norén (2013). These values are accurate enough when working inside a single FOV or a small mosaic, but not over the range of 25 pointings. Therefore, we redefine the angular scale for our 1D center-to-limb curves to be equal to the angular radius of the Sun as viewed from Earth for the date of observation. For campaigns I and II, this is 944 and 945″, respectively. Here we find image scales of 0.0584″ per pixel and 0.0384″ per pixel for the two instruments, respectively. Based on these scales, we can assign a μ value at each pixel in the mosaic. Finally, we binned down the data into 50 bins that each span Δμ = 0.02 in the μ scale.

thumbnail Fig. 3

Comparison of normalized average O I 7772 Å profiles before and after smoothing to remove effects from p modes. Left: un-smoothed limb darkening data of each wavelength point in the O I 7772 Å line on an arbitrary intensity scale. The curves are ordered by wavelength, with the top (red) curve representing the red continuum, the middle (green) curve representing the line core, and the bottom (blue) curve representing the blue continuum. The oscillatory behavior is most pronounced around the line core (green), between a μ of 0.3–0.8. Middle: original representative profiles, ranging from μ = 1 to μ = 0.02. The profile representing μ = 1 is shown in yellow at the top of the figure, with each consecutive profile beneath it in a different color. Right: same as the middle panel, but with most of the p-mode effects removed after limb darkening curves were smoothed.

3.4 Removing p-mode oscillations

After the examination of the limb darkening curve of each wavelength point, we found that our observations are affected by 5-min oscillations (p modes, Leighton 1960). The oscillations in intensity as a function of μ caused by this effect become especially pronounced in the steep parts of the line wings, where a small shift in wavelength becomes a relatively large change in intensity. Typically, p-mode oscillations are removed by averaging over time. Since we do not have a long enough time series for these mosaics, we instead chose to work with the limb darkening curves in the μ scale. This was done by applying a standard box filter, also known as a moving average, with a width of 5Δμ points to the central part (excluding the three outermost points on both sides) of each limb darkening curve (in μ scale) in order to preserve the edges of the limb darkening curves. The smoothed curves are illustrated in the left panel of Fig. 3, while the profiles before and after the smoothing are shown in the two remaining panels. We use both these smoothed μ-averaged profiles as well as the original profiles for all of the following steps in this work and include them in the supplemental material.

4 Results and discussion

We created a set of 50 averaged line profiles for each spectral line along the observed mosaic. Each profile represents a μ position between 0 and 1, spaced by the step of Δμ = 0.02 This means that while each of our profiles represents an equal distance in the μ scale, they span uneven distances in the angular scale, with the first bin spanning nearly 200″ and the last one less than 1″ (see Fig. 2). However, due to the spatial resolution of CRISP and CHROMIS, this means that even the smallest bins still contain an average of hundreds of pixels (see Fig. 2).

Figures 4 and 5 show the CLV curves of the selected wavelength points of the six broad lines. The curves are plotted in two ways: with calibrated atlas intensities in the left column and normalized to the intensity at μ = 1 in the right one. The color of the lines varies as shown in the colorbar above the images, with wavelength points on the blue side of the line core being drawn in blue, and red ones in red. Additionally, we have marked the line center with a black dashed curve and the line continuum as given by Neckel & Labs (1994) with a solid black line. We discuss these plots further in Sect. 4.3.

In order to compare our data to other observations, we calculated a set of parameters for each μ position of each spectral line. For photospheric lines, these are the equivalent width (Wμ, the shift of the line center relative to its position at the disk center, the FWHM, and the line depth. We define the equivalent width for a discretely sampled spectral line as follows: Wμ=i=0nλ(1I(λi,μ)/Ic)Δλi,$ {W_\mu } = \sum\limits_{i = 0}^{n\lambda } {\left( {1 - I{{\left( {{\lambda _i},\mu } \right)} \mathord{\left/ {\vphantom {{\left( {{\lambda _i},\mu } \right)} {{I_{\rm{c}}}}}} \right. \kern-\nulldelimiterspace} {{I_{\rm{c}}}}}} \right){\rm{\Delta }}{\lambda _i}} , $(1)

where Ι(λi,μ)/Ιc is the intensity at a given wavelength point normalized to the continuum intensity Ic, ηλ is the number of wavelength points, and Δλi; is the spacing between each consecutive point.

Our set of broad spectral lines cannot be fully captured within the prefilter range of the instruments. For these lines, we replaced the equivalent width and the FWHM with the “bound equivalent width” (WBμ) and the core width (CW), respectively. For the latter, we followed a convention introduced by Cauzzi et al. (2009). This modified FWHM-like parameter is taken halfway between the average intensity of both sides of the profile at a given separation parameter Αλl and the line core. In Cauzzi et al. (2009), the values for Ha and Ca II 8542 A were chosen to be 900 and 600 mA, respectively, based on a large-scale inspection of spectral profiles where the boundary between the chromospheric core and the photospheric wings was determined. We adopt these values and define our own for the remaining broad lines in Table 3. This allows us to study the behavior of the line cores without being affected by the broad photospheric wings.

For Ca II H & K, we needed yet another approach because of the highly irregular shape of the profiles. Instead of the equivalent width and FWHM, we show the relative shift of the K2 and H2 peaks compared to the line center and their separation as a function of μ, respectively. These peaks can be difficult to find in quiet-Sun profiles, especially the red peak, which tends to disappear in the averaged profiles (e.g., Druzhinin et al. 1987). However, we can still see a clear imprint of both peaks in the wavelength-dependent standard deviation of each averaged profile, which allows us to estimate their location within 50 mA. We refer the reader to Linsky & Avrett (1970), Linsky (2017), and Ayres (2019) for a comprehensive review about the Ca II H & K lines, although the latter two citations primarily focus on the K line so as to avoid effects caused by the Hϵ line.

All of these parameters are evaluated by fitting a univariate spline to each of the profiles. This is done to preserve any asymmetries in the lines that would be lost by fitting a Voigt or similar function. This method breaks down around μ = 0.2, when the continuum intensity starts to decrease rapidly and most lines start showing an emission-line-like behavior. This loss of smoothness in the line profile increases the degeneracy of possible spline fits, most of which have sharp peaks that greatly affect the inferred parameters. In addition, when the viewing angle is almost parallel to the solar surface, the averages contain a smaller number of pixels, so the variations observed in the spectral shapes are subsequently dominated by specific features (such as spicules) and to achieve a reliable average we would need many observations of these regions (Pazira et al. 2017). Therefore, we limit our study to μ = 0.2 and above in order to avoid introducing spurious results.

thumbnail Fig. 4

CLV curves of our chromospheric lines ordered by wavelength. In the left column, each CLV curve is shown in atlas calibrated intensities, in the right column each curve is normalized to one to better illustrate the differences in gradient of each wavelength point. The colors of the lines show their approximate wavelength position within the profile, with blue being the blue wing, green the line core, and red the red wing (see colorbar above the plots). The continuum intensity by Neckel & Labs (1994) is denoted by a solid black line, and the line core by a dashed black line. Figure 5 shows an enlarged version of the first two panels in the left column. We only plot every other wavelength point to preserve clarity in the figure. The black line at μ = 0.08 indicates the point to which we trust the accuracy of our data.

thumbnail Fig. 5

Enlarged version of the first two panels in the left column of Fig. 4.

Table 3

Separation parameters for broad lines.

4.1 Continuum measurements

In order to validate this method versus more classical approaches, we start with a comparison between the continuum limb darkening given by Neckel & Labs (1994) and our observations. For this purpose, we use the continuum points at 4001, 5382, 6172, 6301, and 7771 Å. In Fig. 6, we see that our data agree with the literature values for values of μ larger than 0.12, while below it our values go down steeper. This is due to the fact that Neckel & Labs (1994) excluded data taken within 7″ of the limb to avoid seeing related uncertainties. This corresponds to a μ of roughly 0.12 when using their solar radius value. This means that our data are in full agreement with Neckel & Labs (1994), but that their curves should not be used beyond the aforementioned limit. Neckel (2003, 2005) suggests that the extrapolated values from the Neckel & Labs (1994) data at μ = 0.0 give a meaningful “limb temperature” of 4546 K that is nearly constant over wavelength. This is not the case with our extrapolated curves, which give values between roughly 3000 and 4300 K for the five continuum points that we measured. However, it is important to note that the extrapolation of our fitted spline is very sensitive to the three points closest to the limb, and that changing the intensity of any of these three points by as little as 5% changes the resulting brightness temperature by 200 K.

The uncertainties on these curves are dominated by the 5-min oscillations, as the classical uncertainties become very small due to the large number of pixels over which each point is averaged. Therefore, we estimate an uncertainty of these values of 1%.

Besides the Neckel & Labs (1994) atlas, there are also several local thermodynamic equilibrium (LTE) codes that model limb darkening curves for a wide range of stars, including the Sun. These limb darkening curves can have a strong impact on stellar and exoplanet characterization (e.g., Mandel & Agol 2002; Charbonneau et al. 2002; Spake et al. 2018). In order to examine this assumption, we compare our data to three limb darkening approximations obtained by fitting the ATLAS stellar models2 using the LDCU python package3 (Espinoza & Jordán 2015) and assuming an effective temperature of 5778±58 K, a log(g) of 4.374±0.0005, and a turbulent velocity of 1.5±0.5 km s–1 (Gray 2008). The simplest of these approximations is the linear limb darkening law (Milne 1921): I(λ,μ)I(λ,μ=1)=1u(1μ),$ {{I\left( {\lambda ,\mu } \right)} \over {I\left( {\lambda ,\mu = 1} \right)}} = 1 - u\left( {1 - \mu } \right), $(2)

where I is the intensity as a function of λ and μ and u is the linear limb darkening coefficient. However, linear approximations tend to be too simplistic and cannot reproduce the behavior at the edge of the disk. The four parameter law (Claret 2000) is generally considered the most accurate and a good description of the stellar limb darkening and is as follows: I(λ,μ)I(λ,μ=1)=1i=14ci(1μi/2).$ {{I\left( {\lambda ,\mu } \right)} \over {I\left( {\lambda ,\mu = 1} \right)}} = 1 - \sum\limits_{i = 1}^4 {{c_i}\left( {1 - {\mu ^{{i \mathord{\left/ {\vphantom {i 2}} \right. \kern-\nulldelimiterspace} 2}}}} \right).} $(3)

Here ci represents the limb darkening coefficients. Unfortunately, the quality of the available transit data usually make a four-parameter limb darkening fit degenerate and this likely results in overfitting. For this reason, the squared law is the most widely used (Diaz-Cordoves & Gimenez 1992): I(λ,μ)I(λ,μ=1)=1c1(1μ)c2(1μ)2.$ {{I\left( {\lambda ,\mu } \right)} \over {I\left( {\lambda ,\mu = 1} \right)}} = 1 - {c_1}\left( {1 - \mu } \right) - {c_2}{\left( {1 - \mu } \right)^2}. $(4)

This law is generally considered capable of reproducing the shape of the CLV without being likely to result in overfitting.

Figure 7 shows the trends obtained from the ATLAS model overplotted on the data. It can be seen that while the CLV profiles calculated by the model (shown in red) are accurate within 4% of the disk-center intensity up to μ = 0.08 with the largest discrepancy occurring at 5832 Å towards the limb. It can also be seen that, as expected, the four parameter law most closely follows the calculated profile and the data, while the linear fit provides the least accurate representation of the intensity variations.

Analysis of the coefficients retrieved from the similarly widely used PHOENIX model (Husser et al. 2013) shows that the limb darkening is consistent with both the ATLAS model and the data within 4% of disk-center intensity out to μ = 0.08. Overall, it appears that the model parameters obtained from both the ATLAS and PHOENIX codes are in good agreement with solar data at most of the observed continuum wavelengths, although further investigation at longer wavelengths is required.

thumbnail Fig. 6

Limb darkening as a function of μ of five continuum points for our mosaic, plotted on top of the predicted limb darkening curve from Neckel & Labs (1994) (dashed red) and a univariate spline fitted to the data (red). The vertical black line shows the furthest measurement toward the limb done by Neckel & Labs (1994).

thumbnail Fig. 7

Limb darkening as a function of μ compared to the intensity variation generated for the Sun from the ATLAS stellar model (red) and the resulting linear (blue), quadratic (orange), and four parameter (green) limb darkening laws. The model and all three laws generally match the solar limb darkening out to μ = 0.1 at 4001-6301 Å. The linear law is the most different from the data at all wavelengths due to its simplicity.

4.2 Narrow line parameters

We present the narrow spectral lines in Figs. 9 to 12. Each figure consists of a set of six panels. The first four on the left side represent the equivalent width, line shift, FWHM, and line depth, respectively. The upper right panel shows a solar atlas profile that is convolved4 with the SST instrumental profile, as well as the average profile at μ = 1 and μ = 0.2. The last panel shows a 2D representation of the average profiles as a function of μ. All of these lines, except for O I 7772 Å fall within the spectral range of the IRSOL SS3 spectrum (Ramelli et al. 2017, Ramelli et al. 2019). In these cases, we have used this data alongside our own for comparison. However, no degradation was applied to the IRSOL SS3 spectra to compare the parameters under a different spectral resolution.

C I 5380 Å, shown in Fig. 9, is a relatively weak line that has been shown to virtually disappear at the far limb (Livingston et al. 1977). Our observations agree with this statement as we see that the line disappears almost completely around μ = 0.1. We find that our equivalent widths perfectly match those of the IRSOL SS3 atlas. However, this is not the case for the line shift. Here we see that, unlike the IRSOL data, our data exhibit a sharp increase around μ = 0.35. A similar discrepancy can be seen in the FWHM, where our data exhibit a peak around 0.35 Å. After inspecting the average profiles, we find that they take on an asymmetrical shape around μ = 0.3, but the reason for this is not clear as the mosaic does not seem to contain more activity in this region. A similar effect is discussed in Osipov & Vasilyeva (2019) where a series of iron lines are observed to have an increasing equivalent width closer to the disk center, but then they start to decrease when very close to the limb. Regardless, additional observations of this μ range are required. We also see the Fe I 5379 Å and Ti II 5381 Å lines in this plot.

Fe 16173 Å, shown in Fig. 10, is a strong iron line, commonly observed by SDO/HMI (Lemen et al. 2012). We once again find the equivalent widths of our data and that of the IRSOL SS3 atlas and Balthasar (1988) to be in agreement, while the line shift diverges slightly at lower μ. The differences in FWHM and line depth are caused by the different spectral resolutions of the compared data sets.

Fe I 6301 Å, shown in Fig. 11, is a commonly used photospheric line at the SST (e.g., Díaz Baso et al. 2021b) and Hinode/SP (Kosugi et al. 2007). Here we find that our observations match those of the IRSOL SS3 atlas for the equivalent width and line shift, but both data sets return equivalent width values that are lower than those reported in Balthasar (1988). This is likely due to the fact that our data do not include enough of a continuum. We again note the same divergence for the line shift as was seen in the previous two lines.

O I 7772 Å, shown in Fig. 12, is a strong oxygen line often used for oxygen abundance measurements (Bergemann et al. 2021). This line falls outside of the wavelength range of the IRSOL SS3 atlas. We compared our values to the equivalent width measurements presented in Pereira et al. (2009). We find our Wμ values to be systematically above those of Pereira et al. (2009), but we find a similar gradient in the line (see Fig. 8). The discrepancy in values is likely once again caused by a difference in the definition of the continuum and the bounds between which the equivalent width is calculated.

thumbnail Fig. 8

Enlarged version of the first panel of Fig. 12 showing the equivalent width of the O I 7772 Å line as a function of μ compared with measured values by Pereira et al. (2009). We plot the data (gray), the smoothed data (red), the measurements of Pereira et al. (2009) (blue), and our smoothed data shifted down to match the disk center observation of the blue points in light red. The latter was added to ease comparison in regards to the gradient between the two data sets.

4.3 Broad line parameters

The Ca II K line, shown in Fig. 13, has a more complicated shape compared with most of the other lines that we study. This is because the line source function decouples from temperature just after the temperature minimum, which results in two peaks on either side of the line core. The blue peak (K2V) is sensitive to a larger height range than the red peak (K2R) (e.g., Bj∅rgen et al. 2018). In the quiet Sun, we find that this often leads to the K2R peak disappearing, with about 25% of quiet-Sun Ca II K profiles being reversal-free (Rezaei et al. 2008). This effect is strengthened when averaging, as the already low intensity of the K2R peak gets smeared out. We can see this in the solar atlas (upper right panel in Fig. 13), where there is no clear K2R peak, while a K2V peak is visible. The same is true for our data at μ = 1. However, the location of both peaks can still be inferred from the wavelength-dependent standard deviation of each averaged profile5, which is largest at the location of the peaks. This location is plotted in the top left panel of the same figure. Additionally, we have overplotted observational measurements of the K2R peak location from Bj∅rgen et al. (2018) in the same panel for comparison. While our data are noisier than the literature values, they match the trend for almost the entire range apart from the last few points at the extreme limb, where the peaks seem to be moving closer together again. For the line depth, we have overplotted values from Ding & Fang (1991), which match very well with our data. In Fig. 4, we can see that apart from the noise, the entire profile has a similar behavior up until μ = 0.2, at which point the line wings dip down stronger. This is in line with the derived parameters discussed above, where most of the change in the profile seems to come from shifts in the K2 peaks.

As, shown in Fig. 14, behaves very similarly to Ca II K. We find a lower spread in the R and V peak shifts than we did for Ca II K, but a similar behavior with regard to both peaks moving apart with lower μ until the last few measurements. We find that the line core behaves in a similar way to Ca II K, and that the velocities seem to be slewing toward the blue. In Fig. 4 we also see a similar behavior to that of Ca II K, where most of the profile seems to undergo the same limb darkening effects until roughly μ = 0.2

Hβ, shown in Fig. 15, differs strongly from the IRSOL SS3 atlas values for the bound equivalent width. This is likely due to the chromium blend that can be found in the blue wing of the line, which has a differing depth due to the mismatch in spectral resolution between the two data sets. Nevertheless, our data once more are in strong agreement with the values provided by Ding & Fang (1991). In Fig. 4, we see a slightly larger spread in the limb darkening curves, but with largely the same shape. The line core seems to be dipping slightly faster than the wings.

Mg 15173 Å, shown in Fig. 16, is a line that forms around the temperature minimum. When comparing our data to the IRSOL SS3 atlas, we find that the bound equivalent width matches the IRSOL values closely. This is not the case for the line shift, where we see a curve toward the red at low μ values similar to the behavior of the photospheric lines. In Fig. 4, we see a much bigger spread in the limb darkening curves, with the far wings following the behavior of the continuum, and the core intensities being less affected. This data set also contains the Fe I 5172 Å and the Ti I 5174 Å lines.

Hα, shown in Fig. 17, has similar behavior to what we saw in the Hβ line, where our data strongly differ from the IRSOL data.

The difference in trend between our data and that of IRSOL is likely due to the telluric water line that is visible in the red wing of the IRSOL data, but it is completely washed out in ours. The line depth has a different gradient to that of the IRSOL data, but it matches the data from Ding & Fang (1991) reasonably well, except for the last point. In Fig. 4, we see that the limb darkening behavior of this line is quite different from Hβ. First, there are the wing points that dip strongly, but these can be explained by the fact that Hα samples further into the far wings. Second, the line core is more affected by the limb darkening than the near wings and thus acting oppositely to Hβ. Finally, there is a larger spread in the profiles’ limb darkening curves.

Ca II8542 Å, shown in Fig. 18, behaves similarly to the other chromospheric lines when it comes to the studied parameters. In Fig. 4, we see Ca II 8542 Å behaves in a similar manner to Hβ, with the line core being more affected by the limb darkening than the continuum. However, besides that, we see that the normalized line core closely follows the continuum curve until roughly μ = 0.3.

Finally, we would like to point to the fact that all of our narrow photospheric lines as well as Mg I 5173 Å show a positive (red) shift of their core with decreasing μ. This can be explained by the fact that the disk center has a blue-shifted velocity on average because of granular motions. This average blue shift decreases with μ as we start to look at the granules from the side. However, we find that all of the remaining broad chromospheric lines tend to shift toward the blue with decreasing μ. It is possible that this is because the chromospheric canopy tends to be horizontal and that these velocities enter our line of sight at higher μ angles. The implied preferential direction of the canopy could be explained by the projection effects of a radial structure. At lower μ the part pointing toward the disk center looks bigger than the part pointing away from it, which would result in a net blue shift in the line of sight. The fact that Mg I 5173 Å behaves more like a photospheric line is in line with this given that we do not see horizontal canopy structures in the line core. However, it is also possible that this is caused by an unforeseen bias in the averaging process, and further investigation into individual profile statistics is required.

thumbnail Fig. 9

Summary of the CLV data of the C I 5380 Aline. Left: four panels showing the equivalent width, line shift, FWHM, and line depth, respectively. Right: representative profiles of C I 5380 Å. In the upper panel, we plot a convolved profile from the solar atlas as well as our data at μ = 1 and μ = 0.2. The pair of black lines show the range over which the equivalent width has been calculated. The lower panel shows a 2D representation of all 50 profiles.

thumbnail Fig. 10

Same as Fig. 9, but for Fe I 6173 Å.

thumbnail Fig. 11

Same as Fig. 9, but for Fe I 6301 Å.

thumbnail Fig. 12

Same as Fig. 9, but for O 17772 Å. Additionally, we include the measured equivalent widths by Pereira et al. (2009) in the upper left panel.

thumbnail Fig. 13

Left: four panels showing the K2 peak location, line shift, K2 peak separation, and line depth, respectively. Right: representative profiles of Ca II K. In the upper panel, we plot a convolved profile from the solar atlas as well as our data at μ = 1 and μ = 0.2. The lower panel shows a 2D representation of all 50 profiles.

thumbnail Fig. 14

Same as Fig. 13, but for As.

thumbnail Fig. 15

Left: four panels showing the bound equivalent width, line shift, core width, and line depth, respectively. Right: representative profiles of Hβ. In the upper panel, we plot a convolved profile from the solar atlas as well as our data at μ = 1 and μ = 0.2. The lower panel shows a 2D representation of all 50 profiles.

thumbnail Fig. 16

Same as Fig. 15, but for Mg 15173 Å.

thumbnail Fig. 17

Same as Fig. 15, but for Hα.

thumbnail Fig. 18

Same as Fig. 15, but for Ca II 8542 Å.

5 Conclusions

In this work, we present the results derived from a series of high-resolution mosaics of the broad Ca II H&K Mg I 5173 Å, Hβ, Ha, and Ca II 8542 Å lines, as well as of the four narrow lines C I 5380 Å, Fe I 6173 Å, Fe I 6301Å, and O I 7772 Å. These mosaics, composed of 25 pointings, each spanning roughly 60″ × 60″, allowed us to create a series of 50 high-resolution reference profiles, spaced by 0.02 in the μ scale for each of the observed lines. Additionally, we created a second set of reference profiles that were smoothed using a box smoothing algorithm in order to remove imprints caused by the 5-min oscillations.

We tested the accuracy of the limb darkening curves obtained from these reference profiles against CLVs of atlas continuum intensities and synthetic continuum intensities, and found it to match in all cases (see Figs. 6 and 7). In fact, we were able to extend the continuum limb darkening curves of selected wavelengths past the μ = 0.12 cutoff that was imposed by Neckel & Labs (1994). We compared our continuum curves to stellar limb darkening models and found a generally good agreement above μ = 0.08, where our data become less reliable. Of the most commonly used limb darkening laws, the quadratic and four parameter laws showed the best match, while the linear law deviated the most due to its simplicity. These curves are available in the supplemental material.

We then investigated the narrow lines by measuring the equivalent width, line shift, FWHM, and line depth of each of these lines, as well as by comparing these parameters to existing literature where this was available. In general, we find that our data match existing literature values, and that discrepancies can be explained by differences in the continuum definition and integration range for the equivalent width.

The same was done for the regularly shaped broad lines, with a slight modification to the definition of the equivalent width and FWHM, following conventions started by Cauzzi et al. (2009). This allowed us to study the chromospheric behavior of these lines without it being drowned out by their large photospheric wings, but such a method is very sensitive to the spectral resolution of the studied line. For this reason we suggest the future comparisons using this parameter be done on lines with the same spectral resolution. We find that these chromospheric lines exhibit a small shift that tends toward the blue, a phenomenon that could potentially be caused by projection effects of the chromospheric canopy which tends to flow horizontally. We note that Mg I 5173 Å tends to behave much more like a photospheric line with a line shift toward the red. This could be due to the low formation height of the line, where it is still influenced by the granular velocities.

For the CaII H&K lines, we investigated the separation between the reversal peaks on both ends of the line core. These peaks are not always visible in our spectra but they were found by looking at the wavelength-dependent standard deviation. In both cases, we find that the peaks move away from one another with decreasing μ, and that our values match the observational values presented in Bj∅rgen et al. (2018). As with the other broad lines, we find the Ca II H & K lines tend to shift toward the blue for decreasing μ.

For future work of this kind, we would advise taking longer time series at each mosaic position or taking multiple time-independent mosaics of each line in order to average out the p modes and other local effects, since this is the main source of uncertainty in our data. This would come at the cost of temporal coherence between overlapping mosaic parts, which illustrates the trade-offs that one must consider while making these kinds of observations. A potential way around this would be to apply an integral field spectrometer that can simultaneously capture spatial and spectral information and to drift scan this instrument across the Sun, as was done with single wavelengths. However, the FOV of these kinds of instruments tends to be rather small, which will provide us with fewer pixels to average over.

We hope that these average profiles can be of use to the community for abundance studies and as an incident boundary condition for radiative transfer modelers. We also aim to apply these profiles to improve the intensity calibrations of the SSTRED pipeline.

Acknowledgements

We are grateful to Sepideh Kianfar and Carolina Robustini for participating in the SST observations and to Mats Löfdahl for the stimulating discussions on this topic. We also thank Horst Balthasar and the anonymous referee for their valuable suggestions during the peer-review process. This work was supported through the CHROMATIC project (2016.0019) of the Knut and Alice Wallenberg foundation and by the European Commission’s Horizon 2020 Program under grant agreements 824064 (ESCAPE – European Science Cluster of Astronomy & Particle Physics ESFRI Research Infrastructures) and 824135 (SOLARNET – Integrating High Resolution Solar Physics). This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (SUNMAG, grant agreement 759548). This research is supported by the Research Council of Norway, project number 325491 and through its Centers of Excellence scheme, project number 262622. Contributions by D.P. have been carried out within the framework of the NCCR PlanetS supported by the Swiss National Science Foundation under grants 51NF40_182901 and 51NF40_205606. D.P also acknowledges support of the Swiss National Science Foundation under grant number PCEFP2_194576. The Swedish 1-m Solar Telescope is operated on the island of La Palma by the Institute for Solar Physics of Stockholm University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. The Institute for Solar Physics was supported by a grant for research infrastructures of national importance from the Swedish Research Council (registration number 2017-00625). This research has made use of NASA’s Astrophysics Data System Bibliographic Services. We acknowledge the community effort devoted to the development of the following open-source packages that were used in this work: numpy (http://numpy.org), matplotlib (http://matplotlib.org), astropy (http://astropy.org). We extensively used the CRISPEX analysis tool (Vissers & Rouppe van der Voort 2012), the ISPy library (Díaz Baso et al. 2021a), Specutils (Earl et al. 2022), sTools (Kuckein et al. 2017), and SOAImage DS9 (Joye & Mandel 2003) for data reduction and visualization during the preparation of this article.

References

  1. Alissandrakis, C.E., Patsourakos, S., Nindos, A., & Bastian, T.S. 2017, A&A, 605, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Alsina Ballester, E., Belluzzi, L., & Trujillo Bueno, J. 2021, Phys. Rev. Lett., 127, 081101 [NASA ADS] [CrossRef] [Google Scholar]
  3. Ayres, T.R. 2019, in The Sun as a Guide to Stellar Physics, eds. O. Engvold, J.-C. Vial, & A. Skumanich, 27 [CrossRef] [Google Scholar]
  4. Balthasar, H. 1985, Sol. Phys., 99, 31 [NASA ADS] [CrossRef] [Google Scholar]
  5. Balthasar, H. 1988, A&AS, 72, 473 [NASA ADS] [Google Scholar]
  6. Bergemann, M., Hoppe, R., Semenova, E., et al. 2021, MNRAS, 508, 2236 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bj∅rgen, J.P., Sukhorukov, A.V., Leenaarts, J., et al. 2018, A&A, 611, A62 [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bonnet, R.M., Lemaire, P., Vial, J.C., et al. 1978, ApJ, 221, 1032 [NASA ADS] [CrossRef] [Google Scholar]
  9. Cauzzi, G., & Reardon, K. 2012, IAU Special Session, 6, E5.11 [NASA ADS] [Google Scholar]
  10. Cauzzi, G., Reardon, K., Rutten, R.J., Tritschler, A., & Uitenbroek, H. 2009, A&A, 503, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Charbonneau, D., Brown, T.M., Noyes, R.W., & Gilliland, R.L. 2002, ApJ, 568, 377 [NASA ADS] [CrossRef] [Google Scholar]
  12. Claret, A. 2000, A&A, 363, 1081 [NASA ADS] [Google Scholar]
  13. David, K.H. 1961, ZAp, 53, 37 [NASA ADS] [Google Scholar]
  14. de la Cruz Rodríguez, J., Löfdahl, M.G., Sütterlin, P., Hillberg, T., & Rouppe van der Voort, L. 2015, A&A, 573, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Díaz Baso, C.J., & Asensio Ramos, A. 2018, A&A, 614, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Díaz Baso, C.J., Martínez González, M.J., Asensio Ramos, A., & de la Cruz Rodríguez, J. 2019, A&A, 623, A178 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Díaz Baso, C.J., Vissers, G., Calvo, F., et al. 2021a, https://zenodo.org/record/5608441 [Google Scholar]
  18. Díaz Baso, C.J., de la Cruz Rodríguez, J., & Leenaarts, J. 2021b, A&A, 647, A188 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Diaz-Cordoves, J., & Gimenez, A. 1992, A&A, 259, 227 [Google Scholar]
  20. Ding, M., & Fang, C. 1991, Chinese J. Space Sci., 11, 241 [NASA ADS] [Google Scholar]
  21. Druzhinin, S.A., Pevtsov, A.A., & Teplitskaja, R.B. 1987, Astronomicheskij Tsirkulyar, 1512, 5 [NASA ADS] [Google Scholar]
  22. Earl, N., Tollerud, E., Jones, C., et al. 2022, https://zenodo.org/record/6207491 [Google Scholar]
  23. Espinoza, N., & Jordán, A. 2015, MNRAS, 450, 1879 [Google Scholar]
  24. Fathivavsari, H., Ajabshirizadeh, A., & Koutchmy, S. 2014, Ap&SS, 353, 347 [NASA ADS] [CrossRef] [Google Scholar]
  25. Faurobert, M., Ricort, G., & Aime, C. 2013, A&A, 554, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Gray, D.F. 2008, The Observation and Analysis of Stellar Photospheres (Cambridge, UK: Cambridge University Press) [Google Scholar]
  27. Gunár, S., Schwartz, P., Koza, J., & Heinzel, P. 2020, A&A, 644, A109 [EDP Sciences] [Google Scholar]
  28. Gunár, S., Koza, J., Schwartz, P., Heinzel, P., & Liu, W. 2021, ApJS, 255, 16 [CrossRef] [Google Scholar]
  29. Gurtovenko, E.A., Kostyk, R.I., Orlova, T.V., et al. 1975, Profiles of selected Fraunhofer lines for different positions center-limb on the solar disk (Kiev: Naukova dumka) [Google Scholar]
  30. Hale, G.E., & Adams, W.S. 1907, ApJ, 25, 300 [NASA ADS] [CrossRef] [Google Scholar]
  31. Hammerschlag, R.H., Sliepen, G., Bettonvil, F.C.M., et al. 2013, Opt. Eng., 52, 081603 [NASA ADS] [CrossRef] [Google Scholar]
  32. Husser, T.O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Joye, W.A., & Mandel, E. 2003, in Astronomical Society of the Pacific Conference Series, 295, Astronomical Data Analysis Software and Systems XII, eds. H.E. Payne, R.I. Jedrzejewski, & R.N. Hook, 489 [Google Scholar]
  34. Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3 [Google Scholar]
  35. Kuckein, C., Denker, C., Verma, M., et al. 2017, in Fine Structure and Dynamics of the Solar Atmosphere, 327, eds. S. Vargas Domínguez, A.G. Kosovichev, P. Antolin, & L. Harra, 20 [Google Scholar]
  36. Leighton, R.B. 1960, in Aerodynamic Phenomena in Stellar Atmospheres, 12, ed. R.N. Thomas, 321 [Google Scholar]
  37. Lemaire, P., Gouttebroze, P., Vial, J.C., & Artzner, G.E. 1981, A&A, 103, 160 [NASA ADS] [Google Scholar]
  38. Lemen, J.R., Title, A.M., Akin, D.J., et al. 2012, Sol. Phys., 275, 17 [Google Scholar]
  39. Lind, K., Amarsi, A.M., Asplund, M., et al. 2017, MNRAS, 468, 4311 [NASA ADS] [CrossRef] [Google Scholar]
  40. Linsky, J.L. 2017, ARA&A, 55, 159 [NASA ADS] [CrossRef] [Google Scholar]
  41. Linsky, J.L., & Avrett, E.H. 1970, PASP, 82, 169 [CrossRef] [Google Scholar]
  42. Livingston, W., Milkey, R., & Slaughter, C. 1977, ApJ, 211, 281 [NASA ADS] [CrossRef] [Google Scholar]
  43. Löfdahl, M.G. 2002, SPIE Conf. Ser., 4792, 146 [Google Scholar]
  44. Löfdahl, M.G. 2016, A&A, 585, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Löfdahl, M.G., Hillberg, T., de la Cruz Rodríguez, J., et al. 2021, A&A, 653, 653 [Google Scholar]
  46. Mandel, K., & Agol, E. 2002, ApJ, 580, L171 [Google Scholar]
  47. Milne, E.A. 1921, MNRAS, 81, 361 [NASA ADS] [CrossRef] [Google Scholar]
  48. Molaro, P., & Monai, S. 2012, A&A, 544, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Moore, C.E., Minnaert, M.G.J., & Houtgast, J. 1966, The solar spectrum 2935 A to 8770 A [Google Scholar]
  50. Neckel, H. 2003, Sol. Phys., 212, 239 [NASA ADS] [CrossRef] [Google Scholar]
  51. Neckel, H. 2005, Sol. Phys., 229, 13 [NASA ADS] [CrossRef] [Google Scholar]
  52. Neckel, H., & Labs, D. 1984, Solar Phys., 90, 90 [NASA ADS] [CrossRef] [Google Scholar]
  53. Neckel, H., & Labs, D. 1994, Sol. Phys., 153, 91 [CrossRef] [Google Scholar]
  54. Norén, A. 2013, ArXiv e-prints [arXiv:1311.5148] [Google Scholar]
  55. Osipov, S.N., & Vasilyeva, I.E. 2019, Kinematics Phys. Celest. Bodies, 35, 85 [CrossRef] [Google Scholar]
  56. Pagel, B.E.J. 1955, MNRAS, 115, 493 [NASA ADS] [CrossRef] [Google Scholar]
  57. Pastor Yabar, A., Martínez González, M.J., & Collados, M. 2020, A&A, 635, A210 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Pazira, H., Kiselman, D., & Leenaarts, J. 2017, A&A, 604, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Pereira, T.M.D., Asplund, M., & Kiselman, D. 2009, A&A, 508, 1403 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Pierce, A.K., & Slaughter, C. 1982, ApJS, 48, 73 [NASA ADS] [CrossRef] [Google Scholar]
  61. Pietrow, A.G.M., Kiselman, D., de la Cruz Rodríguez, J., et al. 2020, A&A, 644, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Pietrow, A.G.M., Druett, M.K., de la Cruz Rodriguez, J., Calvo, F., & Kiselman, D. 2022, A&A, 659, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Rachmeler, L.A., Bueno, J.T., McKenzie, D.E., et al. 2022, ApJ, 936, 67 [NASA ADS] [CrossRef] [Google Scholar]
  64. Ramelli, R., Setzer, M., Engelhard, M., et al. 2017, ArXiv e-prints [arXiv:1708.03284] [Google Scholar]
  65. Ramelli, R., Bianda, M., Setzer, M., et al. 2019, in Astronomical Society of the Pacific Conference Series, 526, Solar Polariation Workshop 8, eds. L. Belluzzi, R. Casini, M. Romoli, & J. Trujillo Bueno, 287 [Google Scholar]
  66. Rezaei, R., Bruls, J.H.M.J., Schmidt, W., et al. 2008, A&A, 484, 503 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Scharmer, G. 2017, in SOLARNET IV: The Physics of the Sun from the Interior to the Outer Atmosphere, 85 [Google Scholar]
  68. Scharmer, G.B., Bjelksjo, K., Korhonen, T.K., Lindberg, B., & Petterson, B. 2003a, SPIE Conf. Ser., 4853, 341 [NASA ADS] [Google Scholar]
  69. Scharmer, G.B., Dettori, P.M., Lofdahl, M.G., & Shand, M. 2003b, SPIE Conf. Ser., 4853, 370 [NASA ADS] [Google Scholar]
  70. Scharmer, G.B., Narayan, G., Hillberg, T., et al. 2008, ApJ, 689, L69 [NASA ADS] [CrossRef] [Google Scholar]
  71. Schwarzschild, K. 1906, Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse, 1906, 41 [Google Scholar]
  72. Scullion, E., Rouppe van der Voort, L., Wedemeyer, S., & Antolin, P. 2014, ApJ, 797, 36 [Google Scholar]
  73. Spake, J.J., Sing, D.K., Evans, T.M., et al. 2018, Nature, 557, 68 [NASA ADS] [CrossRef] [Google Scholar]
  74. Stenflo, J.O. 2015, A&A, 573, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Takeda, Y., & UeNo, S. 2019, Sol. Phys., 294, 63 [NASA ADS] [CrossRef] [Google Scholar]
  76. Unsold, A. 1955, Physik der Sternatmospharen, MIT besonderer Berucksichti- gung der Sonne (Berlin, Heidelberg: Springer) [Google Scholar]
  77. van Noort, M., Rouppe van der Voort, L., & Löfdahl, M.G. 2005, Sol. Phys., 228, 191 [NASA ADS] [CrossRef] [Google Scholar]
  78. Vernazza, J.E., Avrett, E.H., & Loeser, R. 1976, ApJS, 30, 1 [NASA ADS] [CrossRef] [Google Scholar]
  79. Vissers, G., & Rouppe van der Voort, L. 2012, ApJ, 750, 22 [Google Scholar]
  80. Vissers, G.J.M., Danilovic, S., Zhu, X., et al. 2021, A&A, 662, A88 [Google Scholar]
  81. Voigt, H.-H. 1950, ZAp, 27, 82 [NASA ADS] [Google Scholar]
  82. White, O.R. 1962, ApJS, 7, 333 [NASA ADS] [CrossRef] [Google Scholar]
  83. White, O.R., & Suemoto, Z. 1968, Sol. Phys., 3, 523 [NASA ADS] [CrossRef] [Google Scholar]
  84. Wittmann, A. 1976, A&A, 48, 121 [NASA ADS] [Google Scholar]

1

Which denotes the ratio between the spectrum at a given μ position and that at disk center.

4

The CRISP and CHROMIS instrumental profiles can be calculated with the “spec” module of ISPy (Díaz Baso et al. 2021a). However, for convenience, we have added precomputed instrumental profiles to the supplemental material.

5

The standard deviation of the intensity of the individual pixels inside each μ bin for each wavelength point in the spectral line.

All Tables

Table 1

Summary of all data used in this paper.

Table 2

FPI tuning positions for each observed line.

Table 3

Separation parameters for broad lines.

All Figures

thumbnail Fig. 1

Screen capture of the primary image guider of the SST showing the HMI magnetogram with 25 pointings in a mosaic that spans from the solar south pole to the disk center, taken on 12 June 2021. The red marker represents the location at which the telescope was aimed at the moment that the image was taken.

In the text
thumbnail Fig. 2

Overview of mosaics taken during Campaign I, starting from the solar limb on the left and going to disk center on the right. The mosaics span roughly 1000 × 80″. Top: continuum near the blue wing of Fe I 6301 Å. Bottom: line center of Ca II 8542 Å. It can be seen that the Sun seems to extend past μ = 0 in the chromosphere.

In the text
thumbnail Fig. 3

Comparison of normalized average O I 7772 Å profiles before and after smoothing to remove effects from p modes. Left: un-smoothed limb darkening data of each wavelength point in the O I 7772 Å line on an arbitrary intensity scale. The curves are ordered by wavelength, with the top (red) curve representing the red continuum, the middle (green) curve representing the line core, and the bottom (blue) curve representing the blue continuum. The oscillatory behavior is most pronounced around the line core (green), between a μ of 0.3–0.8. Middle: original representative profiles, ranging from μ = 1 to μ = 0.02. The profile representing μ = 1 is shown in yellow at the top of the figure, with each consecutive profile beneath it in a different color. Right: same as the middle panel, but with most of the p-mode effects removed after limb darkening curves were smoothed.

In the text
thumbnail Fig. 4

CLV curves of our chromospheric lines ordered by wavelength. In the left column, each CLV curve is shown in atlas calibrated intensities, in the right column each curve is normalized to one to better illustrate the differences in gradient of each wavelength point. The colors of the lines show their approximate wavelength position within the profile, with blue being the blue wing, green the line core, and red the red wing (see colorbar above the plots). The continuum intensity by Neckel & Labs (1994) is denoted by a solid black line, and the line core by a dashed black line. Figure 5 shows an enlarged version of the first two panels in the left column. We only plot every other wavelength point to preserve clarity in the figure. The black line at μ = 0.08 indicates the point to which we trust the accuracy of our data.

In the text
thumbnail Fig. 5

Enlarged version of the first two panels in the left column of Fig. 4.

In the text
thumbnail Fig. 6

Limb darkening as a function of μ of five continuum points for our mosaic, plotted on top of the predicted limb darkening curve from Neckel & Labs (1994) (dashed red) and a univariate spline fitted to the data (red). The vertical black line shows the furthest measurement toward the limb done by Neckel & Labs (1994).

In the text
thumbnail Fig. 7

Limb darkening as a function of μ compared to the intensity variation generated for the Sun from the ATLAS stellar model (red) and the resulting linear (blue), quadratic (orange), and four parameter (green) limb darkening laws. The model and all three laws generally match the solar limb darkening out to μ = 0.1 at 4001-6301 Å. The linear law is the most different from the data at all wavelengths due to its simplicity.

In the text
thumbnail Fig. 8

Enlarged version of the first panel of Fig. 12 showing the equivalent width of the O I 7772 Å line as a function of μ compared with measured values by Pereira et al. (2009). We plot the data (gray), the smoothed data (red), the measurements of Pereira et al. (2009) (blue), and our smoothed data shifted down to match the disk center observation of the blue points in light red. The latter was added to ease comparison in regards to the gradient between the two data sets.

In the text
thumbnail Fig. 9

Summary of the CLV data of the C I 5380 Aline. Left: four panels showing the equivalent width, line shift, FWHM, and line depth, respectively. Right: representative profiles of C I 5380 Å. In the upper panel, we plot a convolved profile from the solar atlas as well as our data at μ = 1 and μ = 0.2. The pair of black lines show the range over which the equivalent width has been calculated. The lower panel shows a 2D representation of all 50 profiles.

In the text
thumbnail Fig. 10

Same as Fig. 9, but for Fe I 6173 Å.

In the text
thumbnail Fig. 11

Same as Fig. 9, but for Fe I 6301 Å.

In the text
thumbnail Fig. 12

Same as Fig. 9, but for O 17772 Å. Additionally, we include the measured equivalent widths by Pereira et al. (2009) in the upper left panel.

In the text
thumbnail Fig. 13

Left: four panels showing the K2 peak location, line shift, K2 peak separation, and line depth, respectively. Right: representative profiles of Ca II K. In the upper panel, we plot a convolved profile from the solar atlas as well as our data at μ = 1 and μ = 0.2. The lower panel shows a 2D representation of all 50 profiles.

In the text
thumbnail Fig. 14

Same as Fig. 13, but for As.

In the text
thumbnail Fig. 15

Left: four panels showing the bound equivalent width, line shift, core width, and line depth, respectively. Right: representative profiles of Hβ. In the upper panel, we plot a convolved profile from the solar atlas as well as our data at μ = 1 and μ = 0.2. The lower panel shows a 2D representation of all 50 profiles.

In the text
thumbnail Fig. 16

Same as Fig. 15, but for Mg 15173 Å.

In the text
thumbnail Fig. 17

Same as Fig. 15, but for Hα.

In the text
thumbnail Fig. 18

Same as Fig. 15, but for Ca II 8542 Å.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.