Issue |
A&A
Volume 621, January 2019
|
|
---|---|---|
Article Number | A78 | |
Number of page(s) | 16 | |
Section | The Sun | |
DOI | https://doi.org/10.1051/0004-6361/201833722 | |
Published online | 11 January 2019 |
Intensity contrast of solar plage as a function of magnetic flux at high spatial resolution
1
Max-Planck-Institut für Sonnensystemforschung, Justus-von-Liebig-weg 3, 37077 Göttingen, Germany
e-mail: kahil@mps.mpg.de
2
School of Space Research, Kyung Hee University, Yongin, Gyeonggi 446-701, Republic of Korea
Received:
26
June
2018
Accepted:
10
November
2018
Magnetic elements have an intensity contrast that depends on the type of region they are located in (for example quiet Sun, or active region plage). Observed values also depend on the spatial resolution of the data. Here we investigate the contrast-magnetic field dependence in active region plage observed near disk center with SUNRISE during its second flight in 2013. The wavelengths under study range from the visible at 525 nm to the near ultraviolet (NUV) at 300 nm and 397 nm. We use quasi-simultaneous spectropolarimetric and photometric data from the Imaging Magnetograph eXperiment (IMaX) and the Sunrise Filter Imager (SuFI), respectively. We find that in all wavelength bands, the contrast exhibits a qualitatively similar dependence on the line-of-sight magnetic field, BLOS, as found in the quiet Sun, with the exception of the continuum at 525 nm. There, the contrast of plage magnetic elements peaks for intermediate values of BLOS and decreases at higher field strengths. By comparison, the contrast of magnetic elements in the quiet Sun saturates at its maximum value at large BLOS. We find that the explanation of the turnover in contrast in terms of the effect of finite spatial resolution of the data is incorrect with the evidence provided by the high-spatial resolution SUNRISE data, as the plage magnetic elements are larger than the quiet Sun magnetic elements and are well-resolved. The turnover comes from the fact that the core pixels of these larger magnetic elements are darker than the quiet Sun. We find that plages reach lower contrast than the quiet Sun at disk center at wavelength bands formed deep in the photosphere, such as the visible continuum and the 300 nm band. This difference decreases with formation height and disappears in the Ca II H core, in agreement with empirical models of magnetic element atmospheres.
Key words: Sun: activity / Sun: photosphere / Sun: chromosphere / Sun: magnetic fields / techniques: photometric / techniques: polarimetric
© ESO 2019
1. Introduction
In addition to quiet-Sun (QS) network, small-scale magnetic elements on the solar surface manifest themselves in the form of plage regions or faculae. Faculae are seen near the solar limb in white light as bright structures often surrounding sunspots and pores. The same structures seen in chromospheric emission are called plages. Here, we have used the latter nomenclature to describe such magnetic flux concentrations seen at wavelengths sampling heights ranging from the low photosphere to the lower chromosphere. The network is seen at the edges of supergranules, while plages are part of active regions (ARs). Network and plage both harbor kG magnetic fields (Stenflo 1973).
At the photospheric level, magnetic elements are located within the intergranular lanes and have been successfully described by flux tubes (Spruit 1976; Solanki 1993). The magnetic field is swept to the dark intergranular lanes by granular motions and concentrated there by the magnetic flux expulsion mechanism (Parker 1963). The intensification of the magnetic field within these elements causes the internal gas pressure at a given height to drop to maintain total horizontal pressure equilibrium with the non-magnetic surroundings. This leads to an opacity depression, so that one can see into deeper layers of the photosphere within such elements. Their brightness is dependent on their size, the extent of the opacity depression (i.e., effectively the magnetic field strength), which together determine the effective heating from the surrounding hot walls, which correspond to the sides of the surrounding granules.
Small-scale magnetic elements are also linked to the heating of the chromopshere and corona. For the chromosphere, this is indicated by the correlation found between the brightness in the cores of Ca II H and K lines and the magnetic field (Skumanich et al. 1975; Schrijver et al. 1989; Ortiz & Rast 2005; Rezaei et al. 2007; Loukitcheva et al. 2009; Kahil et al. 2017), while the importance of small magnetic features for coronal heating has been pointed out by Zhang et al. (1998), Ishikawa & Tsuneta (2009), Zhou et al. (2010), and Chitta et al. (2017). Together, the network and faculae are responsible for the brightening of the Sun with higher magnetic activity on time scales longer than solar rotation (Krivova et al. 2003; Yeo et al. 2017).
The reported brightness of magnetic elements with respect to the quiet surrounding (i.e., their contrast) in the photosphere, along with their morphological structure depends on the wavelength band in which these structures are seen, their position on the solar disk, the selection method used to identify them, and the spatial resolution at which they are observed. Consequently, different observations sometimes gave seemingly contradictory results regarding their photometric and magnetic properties.
Below we give an overview of some of the studies that have looked into the contrast–magnetic field relationship, concentrating on the (near) disk center observations.
At a spatial resolution of 4″, Ortiz et al. (2002) used simultaneous full-disk magnetograms and continuum images (at 676.8 nm) acquired by the Michelson Doppler Imager (MDI) on-board the Solar and Heliospheric Observatory (SoHO). At disk center, the contrast of active-region faculae (identified by their higher magnetic field values) is negative while it is positive for smaller magnetic elements (the network). They attributed the contrast difference to the size of the features, whereby the radiative heating from the surrounding walls is more important and effective when the cross-section of a flux tube is smaller. This study was redone at higher resolution by Yeo et al. (2013) using HMI data. They obtained essentially the same qualitative results, with the differences that the contrast values at weak magnetogram signals in their identified network are negative.
Frazier (1971) used simultaneous observations of photospheric magnetic field inferred from the Fe I 5233.0 Å line and continuum filtergrams at Fe I 5250.2 Å in addition to line-core data in Fe I 5250.2 Å and Ca II K of two active regions at disk center. Whereas continuum contrast increases up until 200 G followed by a decrease at higher fields, the line-core and Ca II K contrasts keep on increasing with increasing magnetogram signal. For the core of another photospheric line, Fe I 6173 Å, Yeo et al. (2013) obtained that the line contrast remains positive at all studied field strengths, but it does decrease somewhat at higher field strengths. For Ca II H and K this qualitative behavior was confirmed by, for example Skumanich et al. (1975), Schrijver et al. (1989), Loukitcheva et al. (2009), and Kahil et al. (2017), although the exact dependence varied.
Title et al. (1992), Topka et al. (1992), and Lawrence et al. (1993) used AR data from the 50 cm Swedish Solar Vacuum Telescope (SVST), while Lawrence et al. (1993) also analyzed the quiet-Sun network. At disk center, and at a spatial resolution of 0.3″ their measured continuum brightness at visible wavelengths (676.8, 525, 557.6, and 630.2 nm) is negative (less than the mean QS intensity) for all magnetogram signals. The line core brightness at 676.8 nm increases until 600 G, then drops monotonically at higher fields in active regions. They reported higher (more positive) contrast in the quiet Sun (Lawrence et al. 1993).
Kobel et al. (2011) looked into the brightness relationship with the magnetogram signal of both, QS and AR plages. They used spectropolarimetric data at 630.2 nm from the Solar Optical Telecope on-board Hinode (spatial resolution of 0.3″). At disk center, the contrast of both regions shows the same behavior when plotted against the longitudinal field strength: contrast peaking at 700 G, and turning over at higher fields. This partly contradicts the results of Title et al. (1992); Topka et al. (1992) and Lawrence et al. (1993). Kobel et al. (2011) also obtained higher contrast in the QS compared to AR plages, which confirmed the results of Title et al. (1992) and Lawrence et al. (1993).
At a spatial resolution of 0.1″, Berger et al. (2007) studied the relationship between the G-band contrast of solar magnetic elements and the magnetic field at different heliocentric angles using the Swedich Solar Telescope (SST; Scharmer et al. 2003). For faculae close to disk center (μ = 0.97), the G-band contrast is positive for fields up until 700 G and becomes negative for higher fields, with sunspots and pores excluded from their study. Also using the SST, Narayan & Scharmer (2010) analysed the dynamical properties of solar magnetic elements close to disk center, using spectropolarimetric data collected by the CRisp Imaging SPectropolarimter (CRISP) in the Fe I line at 630.25 nm, with the magnetic field obtained via Milne-Eddington inversions. After masking the dark pores in their images, their scatterplot of the continuum contrast against the longitudinal component of the magnetic field shows an initial decrease of the contrast for fields less than 200 G, then an increase until 600 G followed by a drop below the QS reference for higher fields.
So far the observational results show a strong disagreement with the predictions of radiative magnetohydrodynamic simulations (MHD) at their original resolution. According to such simulations, in agreement with the flux tube model, the continuum contrast of small-scale magnetic elements should increase monotonically with the magnetic field strength (Vögler et al. 2005; their Fig. 13). In all observational studies the continuum contrast (in visible wavelengths) of quiet-Sun network and faculae decreases at higher fields, even when excluding pores in the case of faculae. This decrease was attributed to the limited spatial resolution of observations (even if they are carried out using space-borne instruments) compared to the higher resolution of MHD simulations (Röhrbein et al. 2011; Danilovic et al. 2013). At the spatial resolution of 0.15″ achieved by the Imaging Magnetograph eXperiment (IMaX) onboard SUNRISE, Kahil et al. (2017) confirmed this hypothesis by demonstrating that the continuum contrast at 525 nm of a quiet-Sun region saturated at the strongest fields. Inspecting this relationship in the larger plage features at the spatial resolution of IMaX is one of the primary aims of the current work.
A qualitative and quantitative study has been carried out by Kahil et al. (2017; hereafter Paper I) for the pixel-by-pixel relationship of contrast to magnetic field in quiet-Sun network at disk center. We used observations recorded by the balloon-borne observatory SUNRISE during its first flight in 2009 (SUNRISE I), in the UV and visible wavelength ranges. Here, we have looked into the same relationship, but for plage regions observed by SUNRISE during its second flight in 2013 (SUNRISE II). The data and analysis techniques are described in Sect. 2. In Sect. 3 we show our results and compare them to similar studies and to those obtained in Paper I for quiet-Sun data. In Sect. 4 we discuss our findings and come up with conclusions and insights for future investigations.
2. Observations and data reduction
The time series used in this study is recorded with the 1-m telescope on the balloon-borne observatory SUNRISE (Solanki et al. 2010; Barthol et al. 2011; Berkefeld et al. 2011) during its second flight on 2013 June 12 (Solanki et al. 2017). We use simultaneous spectropolarimetric and imaging data collected by the two scientific instruments onboard, an imaging spectropolarimeter (IMaX; Martínez Pillet et al. 2011) and a UV filter imager (SuFI; Gandorfer et al. 2011). The observations were made close to disk center (cosine of heliocentric angle μ = 0.93), targeting the active region (AR) NOAA 11768 (see Sect. 2.2 for more details on the AR field of view).
2.1. IMaX data
The spectropolarimetric data were acquired by the Imaging Magnetograph eXperiment (IMaX; Martínez Pillet et al. 2011) between 23:39 and 23:55 UT in their V8-4 mode. This mode consists of measuring the full Stokes vector (I, Q, U, V) at eight wavelength positions around the center of the photospheric Fe I 5250.2 Å line, with four accumulations at each wavelength. The spectral positions sampled by IMaX were located at Δλ = −120, −80, −40, 0, +40, +80, +120 mÅ from the line center of rest wavelength λ0 = 5250.2 Å with an additional one in the red continuum (at Δλ = +227 mÅ). With an exposure time of 250 ms for each of the four polarization measurements at each of the eight spectral positions, the total cadence achieved for this mode was 36.5 s. The 17 min time series consists of 28 sets of observations, that is, images at 8 wavelengths in all four Stokes profiles, each image covering 51″ × 51″ on the solar surface with a platescale of pixel−1.
During the flight, images were stabilized, by the use of a tip-tilt mirror controlled by the Correlating-Wave Front Sensor (CWS; Berkefeld et al. 2011). Since there is a time lag in the acquisition of the polarization states between two successive wavelength scans, they were interpolated with respect to time to compensate for the solar evolution time during the acquisition of each observation cycle. In addition, images were corrected for instrumental polarization, and for low-order wavefront aberrations to improve the spatial resolution (Martínez Pillet et al. 2011). At this stage, data are called phase-diversity (PD) reconstructed (level 2.2), with a noise level of 7 × 10−3Ic in Stokes V (Ic is the Stokes I continuum intensity) and a spatial resolution of 0.15″ − 0.18″.
Finally, a very simple approach was taken to correct for the instrumental straylight by subtracting 25% of the spatial mean Stokes I profiles from the individual profiles at each of the 8 sampled wavelengths. The Stokes Q, U and V data were not corrected for straylight since it is considered to be non-polarized (see Riethmüller et al. 2017 for more details on the stray-light correction). After stray-light correction, data are called level 2.3 data.
2.2. Inversions
After applying the above corrections, IMaX data were inverted to retrieve the physical parameters needed for our study. The Stokes inversion code SPINOR (Frutiger et al. 2000), which uses the STOPRO routines for the radiative transfer (Solanki 1987) was employed. The main retrieved parameters of relevance for this study are the magnetic field strength (B), the line-of-sight velocity (VLOS), and field inclination (γ). A simple one-component atmospheric model was considered, with three optical depth nodes at logτ = −2.5, −0.9, 0 for the temperature and a height-independent magnetic field vector, line-of-sight velocity and micro-turbulence. More details on the inversion strategy can be found in Solanki et al. (2017) and Kahil et al. (2017). The line-of-sight velocity maps obtained from the inversions were corrected at each pixel for the wavelength blueshift caused by the Etalon used for the spectral analysis, and a constant velocity (0.6 km s−1) was removed so that the spatially averaged velocity across the FOV is zero.
Figure 1 shows an example IMaX continuum map (at Δλ = +227 mÅ from λ0, left panel) and the corresponding magnetogram (longitudinal magnetic field, BLOS map) returned by the inversions. Figure 1 displays a photospheric zoo of different structures with various photometric and magnetic properties: a large pore, which has some penumbral structure attached to it, smaller pores with kG magnetic fields (black boxes), a flux emergence (FE) region (enclosed in the dashed green box), a small QS internetwork region (dashed red box) and plage (partly enclosed by the dashed yellow ellipse). The latter is located mainly in the vicinity of the pores and shows elongated structures or “ribbons” in the intergranular lanes as described in Berger et al. (2004). Such structures (wherever present in the FOV; i.e., not just in the yellow ellipse) will be analyzed in Sect. 3.4.
![]() |
Fig. 1. Left panel: continuum contrast at 525 nm. Right panel: its co-spatial and co-temporal longitudinal magnetic field map retrieved from the inversions. The black boxes (solid lines) contain small pores characterized by contrast below unity and kG magnetic fields. The dashed green box encloses an area with emerging flux. The dashed yellow ellipse outlines a region of plage composed of magnetic elements embedded in intergranular lanes. The dashed red box contains a quiet-Sun internetwork region with weak fields (average of 20 G) and mean contrast of unity (by definition). |
2.3. SuFI
Simultaneously with IMaX, the Sunrise Filter Imager (SuFI; Gandorfer et al. 2011) acquired high resolution images sampling the low-mid photosphere (300 nm, bandwidth of 4.4 nm) and the lower-chromosphere (397 nm; core of Ca II H, bandwidth of 0.11 nm) of the same region but with a narrower FOV of 15″ × 38″. The SuFI FOV covered a part of the leading polarity flux and flux-emergence regions, in addition to small pores and plage regions (see yellow and red boxes overlaid on IMaX FOV in Figs. 2c and d). The integration time for both 300 nm and 397 nm images is 500 ms, with a cadence of 7 s, and the corresponding plate scale is pixel−1. The spatial correction for wave-front aberrations was done using the Multi-Frame Blind Deconvolution technique (MFBD; van Noort et al. 2005). The spatial resolution achieved at 397 nm after reconstruction is approximately 70 km. The reconstructed data are then corrected for instrumental straylight using the solar limb profiles recorded during the SUNRISE I flight. A stray-light modulation transfer function (MTF), with which images are deconvolved, is obtained upon comparing the observed limb profiles to those in the literature at each of the observed wavelength bands (Feller et al., in prep.).
![]() |
Fig. 2. SuFI images coaligned to the IMaX FOV shown in Fig. 1. Panel a: SuFI 300 nm. Panel b: SuFI Ca II H. Panel c: IMaX continuum intensity (normalized to the mean quiet-Sun intensity) to which SuFI 300 nm (overlaid in red) is aligned. Panel d: IMaX line-core intensity derived from the inversions (see Sect. 2.5) and normalized to the local continuum (I+227), and to which calcium images (overlaid in yellow) are aligned. The yellow and red boxes are slightly misaligned due to the differential offsets in the SuFI wavelegnth channels. The gray scale is set to cover two times the rms range of each image. |
2.4. Alignment
Since we are comparing images recorded by SuFI with magnetograms obtained by IMaX, a precise alignment procedure has to be followed. We have applied here the same alignment procedure described in Paper I for the quiet-Sun data, mainly resampling SuFI images to the pixel size of IMaX images, and using a cross-correlation technique to compute the shifts in both x and y directions. The SuFI images at 300 nm were aligned with IMaX Stokes I continuum images, since both of these spectral regions form in the low photosphere and show a normal granulation pattern (although bright points are more prominent at 300 nm, see Fig. 2 and Riethmüller et al. 2010). The Ca II H data were aligned with IMaX Stokes I line-core images. Both wavelength bands share similar spatial structures, such as the reversed granulation pattern (dark granules, bright intergranular lanes), a result of forming higher in the atmosphere (see Figs. 2b and d). In the calcium images also other features besides reversed granulation are used (as the latter is best seen outside active regions). These include brightenings associated with strong magnetic fields, although these are somewhat diffuser in Ca II H than in the line core of IMaX due to the expansion of the magnetic field, which forms a low-lying magnetic canopy in plage (Jafarzadeh et al. 2017a). We note that since the IMaX line is weak, the IMaX line-core image shows a reversed granulation pattern only if normalized by the local continuum (see Fig. 2d).
2.5. Contrast
At each pixel, the Stokes I continuum intensity is determined as the highest intensity of the synthetic best-fit profile retrieved from the inversion of IMaX data. We will consider these values instead of the ones observed at Δλ = +227 mÅ, which can deviate from the true continuum if the neighboring line (Fe I 5250.6 Å) influences the 227 mÅ intensity (due to a wavelength shift or broadening of the line).
We have also referred to the inverted profiles to compute the IMaX line-core intensity as follows: we use the VLOS values returned by the inversions to derive the wavelength shift from the rest wavelength of the Fe I line (λ0 = 5250.2 Å) due to the Doppler effect. We added to it the shifts caused by the etalon blueshift and the convective blueshift that were already subtracted from the inverted VLOS maps, this gives the total shift of the line with respect to λ0. The line-core intensity is the Stokes I value at the total computed wavelength shift. This approach was adopted since at pixels with higher magnetic fields, where Zeeman splitting is important, a simple Gaussian fit to Stokes I profiles was found to underestimate the line-core intensity. Moreover, the profiles fitted to IMaX data points by the inversions are far superior to Gaussian fits, since the latter do not take into account the various effects that contribute to the shape of the spectral line (temperature, density, magnetic field, etc.).
The relative intensity at each wavelength band was computed at every pixel by normalizing to the mean quiet-Sun intensity of a very quiet region (of weak magnetic fields). The contrast is the relative intensity minus unity, but we will use this terminology for the computed relative intensity for a consistent comparison with the literature. Hence the mean quiet Sun level will be attributed to a contrast of unity, brighter (darker) pixels than the mean QS will have a contrast above (below) unity. For the IMaX continuum and line core, the QS region is outlined by the dashed red rectangle in both panels of Fig. 1.
Due to the narrower field of view (FOV) of SuFI, finding a quiet-Sun region in these data is more challenging, particularly as the QS region has to contain a large number of granules. Locating this region in the calcium images is even harder since one should look for dark regions with low fields, but due to the long fibrils found almost everywhere in the images (also protruding into neighboring quiet regions), most of the low field regions are rather bright in calcium emission, affecting the computed mean quiet-Sun value.
Luckily, quiet-Sun images are available as flatfields, which are sufficient for the purpose of obtaining the average quiet-Sun intensity. Obviously flatfields were recorded before and after the observations, but not during them. Fortunately, the photon flux at 397 nm does not depend on elevation angle (and hence time) so that the average quiet-Sun intensity does not vary with time and can be obtained straightforwardly from the flatfields. However, the 300 nm wavelength exhibits a cyclic variation due to absorption by the Earth’s atmosphere. We have used this day-to-night cycle variation to derive the mean quiet-Sun intensities at the times of our data acquisition. In Appendix A we provided more details on the computation of the mean quiet-Sun intensities in SuFI data.
2.6. Masking
We are mainly interested in studying the relationship between contrast and magnetic field of magnetic elements in plages. These features are distributed everywhere in the FOV around the large and small pores. Hence, pores were masked by using an intensity threshold in the IMaX continuum maps. We have computed the mean and standard deviation (σ) of each continuum image. After applying a low-pass filter consisting of the running mean over 33 × 33 pixels (to smooth the image), we discarded pixels with intensities in the smoothed image lower than the mean minus σ value. The smoothing greatly reduces the granular contrast. In particular intergranular lanes are much less dark, so that they lie above the threshold of mean minus σ in the smoothed image. However, pores remain darker and the smoothing also leads to the exclusion of their immediate surroundings which are darkened by the pores after smoothing. To ignore pixels around the pores, which do not belong to magnetic elements (see Sect. 3.4 for a further discussion on the masking technique), the surroundings of pores in addition to the flux-emergence area (black and green boxes in Fig. 1) were masked manually. The flux-emergence region can be identified more easily when examining movies of Stokes I and V images at the inner wavelength points, that is, at Δλ = −40, 0, +40 mÅ from line center.
All unmasked pixels of the 28 images composing the IMaX time series were included in the pixel-by-pixel scatterplots analyzed in the next sections. For studying the contrast at 525 nm vs. BLOS, the pixels in the whole IMaX FOV of 51″ × 51″ are taken, while in the UV wavelengths we consider the pixels in the smaller FOV of 15″ × 38″ with the corresponding cropped BLOS maps from the inverted IMaX data. The masks computed from IMaX continuum maps were also used to mask the magnetic features in the corresponding IMaX magnetograms and SuFI images at 300 nm and 397 nm.
3. Results
3.1. Scatterplots of IMaX continuum and line-core constrasts vs. BLOS
After masking the pores, we plot in Figs. 3 and 4 pixel-by-pixel scatterplots of the contrast in the IMaX continuum and line-core of Fe I 5250.2 Å versus the line-of-sight component of the magnetic field, BLOS. To show the trend of data points in the scatterplots, the contrast values were averaged into bins, each containing 500 data points. The averages for each bin are overplotted in red in both figures.
![]() |
Fig. 3. Scatterplot of the continuum contrast of Fe I 5250.4 Å vs. the LOS component of the magnetic field, BLOS, for the AR plage. Pores and the FE region are excluded. The horizontal dashed red line is the mean QS continuum intensity level (the mean contrast of the pixels in the dashed red box in Fig. 1). The red curve joins the average contrast values in bins of 500 data points each. |
![]() |
Fig. 4. Same as Fig. 3 but for the line core of Fe I 5250.4 Å. The blue curve is the logarithmic fit to the averaged line-core contrast values starting at 120 G. |
The large scatter around BLOS = 0 in Fig. 3 is due to the distribution of the magnetic field in the granulation (Schnerr & Spruit 2011). At low magnetic field values (BLOS < 150 G), the contrast in the visible continuum decreases with increasing BLOS since those fields are being dragged to the intergranular lanes by flux expulsion (Parker 1963). In an image, these pixels (not shown here) are located in the interiors of bright granules and dark intergranular lanes with weak magnetic fields. For fields above 200 G, which correspond to magnetic features located in the intergranular lanes and undergoing convective instability, the contrast increases with BLOS to reach contrast values above unity at 600 G. The resulting shape of the curve at low fields is referred to as the fishhook (Schnerr & Spruit 2011).
Above 600 G, the average continuum contrast is > 1 and continues to increase until it reaches a maximum at Bmax = 855 G, then decreases again, dropping below unity at around 1200 G. The pixels belonging to this range (around the peak and after the turnover) are mainly located in the intergranular lanes, as is typical of small-scale magnetic elements (Solanki 1993). The numerical values of the longitudinal magnetic field obtained above (at the minimum of the fishhook and at the highest contrast, Bmax) were computed by modeling the scattered data points using one of the non-parametric regression techniques described in the Appendix of Paper I. We prefer the non-parametric regression over a simple parametric model, such as the polynomial fit used in Kobel et al. (2011), to find the contrast and BLOS at the peak of their scatterplots, for two reasons. Firstly, the polynomial fit is highly dependent on the lower and upper limits of BLOS within which the averaged contrast values are fitted. Secondly, it is also dependent on the degree of the fitting polynomial.
The IMaX continuum contrast in plage exhibits the same qualitative behavior as the observations of Frazier (1971), Narayan & Scharmer (2010), Kobel et al. (2011), Kostik & Khomenko (2012), and Buehler et al. (2018) and the degraded simulations of Danilovic et al. (2013): a peak in the contrast at intermediate field strengths followed by a turnover at higher fields. The main difference is that the magnetic field value at which the contrast reaches its maximum (850 G) is higher than the values of 700 G and 650 G reported by Kobel et al. (2011) at their spatial resolution of 0.3″ and by Narayan & Scharmer (2010) at 0.15″, respectively. These values support the correlation between the position of the peak and the corresponding spatial resolution of the observations, as pointed out by Yeo et al. (2013).
The scatterplot of the line-core contrast versus BLOS in Fig. 4 shows a weaker fishhook shape compared to the visible continuum. The monotonic increase of the line-core contrast with BLOS was not seen in similar studies where line-core data were available (Title et al. 1992; Yeo et al. 2013). Their scatterplots of the line-core contrast drop at high magnetogram signals, even after masking out the pores. Following Paper I we have modeled this relationship with a logarithmic function of the form I(B)=β + α.logB, which is overplotted in blue in Fig. 4. We could not compare our fit results to any other study simply because the non-linear monotonic shape of the contrast in a similar wavelength band vs. BLOS is not obtained by lower spatial resolution studies. Nevertheless, we have performed them to serve as a reference to which future studies in the same wavelength band, or studies based on radiative MHD simulations, can be compared. The IMaX continuum contrast could not be modeled with a logarithmic function given the non-monotonic increasing shape. This behavior will be analyzed in Sect. 3.4.
3.2. Scatterplots of UV contrast vs. BLOS
Scatterplots of the contrast at 300 and 397 nm against BLOS are shown in Figs. 5 and 6, respectively. Contrast values are averaged following the method described in Sect. 3.1.
![]() |
Fig. 5. Scatterplot of the SuFI intensity contrast at 300 nm vs. BLOS. The red curve joins the averaged contrast values inside each bin, the blue curve is the logarithmic fit to the red curve starting from 100 G. |
![]() |
Fig. 6. Same as Fig. 5 but for the line core of Ca II H. The blue curve is the logarithmic fit to the red curve starting from 50 G. |
At 300 nm, the fishhook shape is still visible due to the relatively low formation height of this wavelength band of around 50 km, as computed by Jafarzadeh et al. (2017a). After an initial decrease with BLOS, the contrast at 300 nm reaches its minimum at around 100 G. Beyond that point, the contrast increases with field strength until it saturates at higher BLOS values.
The contribution function of the Ca II H bandpass, computed from an atmospheric model corresponding to AR plages and showing a mean formation height of 550 km (Jafarzadeh et al. 2017b), implies that this line has a contribution from the photosphere, but also from the low-to-mid chromosphere. The fishhook shape in the scatterplot of Fig. 6 is gone, due to the absence of granulation at the larger formation height of this line compared to the 300 nm bandpass and the 525 nm continuum. The features are on average brighter than the QS even at very weak fields, due to the presence of bright fibrils that are prominent even in weak magnetic field regions. Weak field regions that are dark can be seen mainly in the upper middle and lower left regions of the calcium images, while fibrils are found nearly everywhere (see Fig. 2). The contrast then increases with BLOS until it saturates at higher fields. Pixels in this region of the scatterplot belong to plages as identified in the IMaX continuum images.
In the past, a power-law function has been used to describe the relation between Ca II H brightness and magnetic flux density for network and IN features both, within limited fields of view (Rezaei et al. 2007; Loukitcheva et al. 2009) and for magnetic features in full-disk observations (Ortiz & Rast 2005; Chatzistergos 2017). Here, however, we have followed Paper I, where we found that a logarithmic function provides a superior representation to the data at the high spatial resolution reached by SUNRISE. The fit parameters are shown in Table 1.
Parameters of logarithmic fits to the contrast in the core of Fe I line at 5250.2 Å (starting from 120 G) and for 300 nm and 397 nm (starting from 100 G) vs. BLOS.
3.3. Comparison between wavelengths
We plot in Fig. 7 the binned contrast curves for the wavelengths in the UV at 300 and 397 nm, and in the visible continuum at 525.04 nm plus line core of Fe I at 525.02 nm derived from the inversions as explained in Sect. 2.5. We plot for comparison the curve corresponding to the line-core contrasts derived from the Gaussian fits to individual IMaX Stokes I profiles. Histograms of the measured contrasts at the studied wavelengths are displayed in Fig. 8. These histograms are for the common FOV (13″ × 34″) of IMaX and SuFI data sets for a one to one correspondence of the contrasts. All the pixels of the time series except for those belonging to pores are included in these histograms. They are plotted for different ranges of BLOS, these ranges correspond to magnetic field values where the scatterplot of the contrast at 525.04 nm vs. BLOS shows different behavior: the weak field region (10 G < BLOS < 200 G), the range where features have an average contrast above unity (600 G < BLOS < 1000 G), and the range where features start to be darker than the mean quiet-Sun intensity (BLOS > 1200 G).
![]() |
Fig. 7. Intensity contrasts of plage vs. BLOS in the wavelengths considered in the current study (see the legend in the figure). Plotted are the contrast values (plotted in red in Figs. 3–6). For comparison, also plotted is the binned line-core contrast at 525.02 nm derived from Gaussian fits to individual Stokes I profiles. |
![]() |
Fig. 8. Contrast histograms of pixels belonging to the common FOV of SuFI and IMaX after masking out pores. Colors correspond to different wavelengths, and panels correspond to different ranges of BLOS values (see text for a detailed explanation). The histograms are normalized such that the integral over each is equal to one. |
For weak fields (Fig. 8a), 525 nm continuum (green) and 300 nm (violet) have averaged contrasts < 1. Both bands exhibit lower contrasts than the rest of the wavelengths at all field strengths, with the contrast at 300 nm becoming larger than 525 nm with increasing BLOS (Figs. 8b and c).
The higher contrast measured in the UV (at 300 and 397 nm) compared to the visible continuum at 525.04 nm is explained by the stronger response of the Planck function to the temperature variations at shorter wavelengths plus the greater height of formation in the case of the Ca II H line core. For the Fe I line core, the reason for the larger contrast, is this line’s great temperature sensitivity and its greater formation height. For all values of BLOS, Ca II H (red) exhibits the highest contrasts.
Figures 7 and 8 also show the difference in the computed line-core intensity from the inversions (blue) and the Gaussian fits (orange), which is more notable at higher fields (Fig. 8c). As explained in Sect. 2.5, this difference is due to the underestimation of the line-core intensity at higher fields if derived from the Gaussian fits. The curves obtained by the two methods agree rather well, except at the very highest field strengths, suggesting that the Zeeman splitting plays a minor role at the spectral resolution of IMaX.
3.4. Why does the continuum contrast display a turnover?
In this subsection we will try to understand why the averaged continuum contrast of plage features at 525 nm peaks at intermediate field strengths and decreases at higher fields when observed by IMaX near disk center.
3.4.1. Treatment of magnetic signals around pores
The derived contrast as a function of BLOS near disk center is believed to vary according to the treatment of pixels around pores and sunspots (Kobel et al. 2011; Yeo et al. 2013). Therefore, care has to be taken when removing pores since plages are often located close to their boundaries, and inclusion of pixels which might belong to pores could affect the contrast–BLOS relationship in the plage areas. We have tested this by extending our masks around the pores, described in Sect. 2.6, to exclude more of their surroundings. We then again produced the scatterplot to see if the shape described above persists or not. The shape of the scatterplot (not shown here) and in particular the position of the contrast peak were hardly affected.
3.4.2. Masking technique
The criterion we used here to exclude pores from our data is based solely on the IMaX continuum intensity (see Sect. 2.6). We have also tried a mask based on a contrast-BLOS threshold, such as that employed by Kobel et al. (2011). The threshold used for creating their masks (contrast < 0.85 and BLOS > 900 G) aimed to remove the inner dark and magnetic parts of pores. In addition, they extended their masks around pores to exclude the pixels darkened by telescope diffraction. However, in our data this threshold corresponds to micropores that are part of the plage, and also to the dark magnetic pixels found in the ribbon-like features which result from their possible inclination with respect to the solar normal (see Sect. 3.4.3). Thus, adopting the same contrast-magnetic field threshold in our data will remove the pixels contributing to the downturn in the bottom right part of our scatterplot in Fig. 3 in addition to a few brighter pixels if those masks were extended, since micropores in our data tend to have brighter edges (see Figs. 10a and b). This points to a dependence of the photospheric continuum contrast–BLOS relationship on the criterion used to distinguish between magnetic elements and pores. This also implies that a careful examination of the magnetic and photometric distribution of small-scale magnetic elements has to be done before adopting any masking criterion, which is only possible if features are well resolved, as is the case with IMaX (see Sects. 3.4.3 and 3.6).
3.4.3. Spatial resolution
In lower spatial resolution studies (Title et al. 1992; Topka et al. 1992; Lawrence et al. 1993) using the SVST, the continuum contrast monotonically decreases with magnetogram signal (see Introduction). The inability to resolve magnetic regions, or the fact that they are partly blended with the non-magnetic surroundings, primarily with the dark intergranular lanes, were considered the reason for the negative contrasts at higher field strengths. This was confirmed by Danilovic et al. (2013) using MHD snapshots of a plage region. After degrading the simulations to the spatial resolution of Hinode/SP, their scatterplots turn from a monotonic relationship at the original MuRAM resolution to a peak and a turnover, in agreement with the findings of Kobel et al. (2011).
Danilovic et al. (2013) showed that smearing isolated bright magnetic structures embedded in a dark, nearly field-free environment shifts the data points with high contrast and BLOS at the original resolution to a region of the scatterplot with lower contrasts and intermediate field strengths. The accumulation of data points there gives rise to a peak upon averaging the contrast values. Bigger and less bright magnetic structures are only slightly affected by the smearing and remain dark with higher field strengths, causing the turnover in the right side of the scatterplot. This interpretation of the effect of finite spatial resolution is valid for regions with low average magnetic field such as the quiet-Sun network, where the magnetic features have little internal structure. Thus, it works well for the data described in Paper I. However, The plage features in our data are larger than those found in the quiet Sun, so that they should be affected less by spatial smearing, as test have confirmed. Instead the magnetic features in the plage show more complex photometric and magnetic distributions than the simple picture of magnetic bright points surrounded by dark field-free lanes.
To isolate these structures, we located the pixels with BLOS > 1200 G and binned contrast < 1 as seen in Fig. 3. We call them the “turnover” pixels. We also located the pixels centered at ± 200 G from Bmax (i.e., 600 G < BLOS < 1000 G) with an averaged contrast > 1, and call them the “peak” pixels. Figures 9a and b show the contours enclosing both families of pixels. We find that together they form magnetic structures of different sizes and appearences embedded in intergranular lanes and spread mainly around the large and small pores, with the turnover pixels (yellow contours) typically being surrounded by the peak pixels (red contours). Small structures have the shape of circular bright points, while slightly larger ones resemble the so called “striation” or “ribbon” in the G-band observations of Berger et al. (2004). Figures. 9c, d, and e are blowups of three features, picked to illustrate three distinctive types: a micropore (Fig. 9c), bright points (Fig. 9d), and a ribbon-like feature (Fig. 9e). The red and yellow contours now enclose the bright (contrast > 1) and dark (contrast < 1) parts of the peak pixels, while the green and orange contours enclose the bright and dark parts of the turnover pixels. Clearly, the smallest features, the bright points in Fig. 9d, contain only bright pixels, while the larger features display an increasing fraction of dark pixels as their area increases. This behavior is found to be typical.
![]() |
Fig. 9. Panels a and b: IMaX data for the same time as in Fig. 1, now with the absolute value of BLOS plotted in the right panel. The red contours enclose the “peak” pixels with > 1 averaged contrast in the continuum of 525 nm and 600 G < BLOS < 1000 G, while the yellow contours enclose the “turnover” pixels with averaged contrasts < 1 and BLOS larger than 1200 G. Panels c, d, e: blowups of three magnetic features of different sizes. Red and green contours enclose the bright (contrast > 1) parts of both families of pixels, while orange and yellow contours enclose the dark (contrast < 1) parts. |
In Fig. 10 we plot horizontal cuts through four features belonging to the three types of magnetic elements shown in Fig. 9. Displayed are the continuum intensity at 525 nm and the LOS magnetic field. The cut is perpendicular to the limb, with the nearest limb being to the right of the figure, and disk center to the left.
![]() |
Fig. 10. Profiles of BLOS (green curve) and continuum intensity contrast at 525 nm (blue curve) along cuts through four different magnetic structures in an example IMaX image (see text for details). The red dashed curves are the Gaussian fits to the BLOS profiles. The FWHM of these Gaussian fits are given at the top of each frame. The horizontal dashed blue line is where the contrast is equal to unity. The limb is to the right of the figure, while disk center is to the left |
In all four features the LOS component of the magnetic field peaks in their centers and decreases toward their edges. The continuum contrast at 525 nm, however, shows three distinctive patterns according to the sizes of the features (given above each frame): for large features (horizontal extent of 500 km–600 km) it is > 1 at the edges and < 1 at the center where BLOS peaks (Figs. 10a and b). For small bright points (100 km–200 km), it is > 1 in their centers and the brightness peak coincides with the position where BLOS peaks (Fig. 10c). In the elongated features with intermediate sizes (200 km–400 km), the continuum contrast is > 1 in the limb side of the feature, and < 1 in the disk center side, with a neutral contrast in the center (Fig. 10d). The latter pattern is seen in most of the features in the IMaX data. The feature size is determined from a Gaussian fit (dashed red curves in Fig. 10) to the BLOS profile.
The profiles of the continuum intensity in Figs. 10a and b resemble the symmetric double-humped profiles at disk center in both G-band observations and numerical simulations (Berger et al. 2004; Steiner 2005). While the profile in Fig. 10d resembles the asymmetric pattern of contrast in magnetic elements near the limb (Hirzberger & Wiehr 2005; Steiner 2005). Given that our observations were not carried out exactly at disk center (μ = 0.93), it could be that the shape of the intensity profiles seen across the mid-sized features are due to an inclined line of sight (the hot wall effect). But it could also be a result of inclined magnetic fields with respect to the solar normal, as proposed in Keller et al. (2004).
Variations of the intensity contrast and magnetic field distributions across IMaX features with different sizes explain the scatter we see in the plot of IMaX continuum contrast vs. BLOS in the range 600 G–2000 G, which corresponds to magnetic elements belonging to plage. We conclude that the scatter around the mean QS level is real and is not a consequence of poor spatial resolution or noise in the data. This point is strengthened later when inspecting the brightness-magnetic field relationship in quiet-Sun areas located within the FOV of the AR observations (Sect. 3.6).
3.5. Comparison with 2009 quiet-Sun data
In this section we present a qualitative and quantitative comparison of the contrast–magnetic field relationship near disk center in AR plage with that of quiet Sun, as obtained in Paper I. This comparison is obviously restricted to the wavelengths that were recorded during both flights (IMaX continuum and line core, 300 nm continuum, and Ca II H line core at 397 nm). The reason we didn’t compare with the quiet regions extracted from the active-regions scans is that the quiet Sun is located mainly outside the SuFI FOV. In the next section we will analyze the IMaX properties of quiet-Sun features from 2013 scans, and compare them to the 2009 results to test for consistency.
We show in Fig. 11 the averaged contrasts vs. BLOS obtained from the quiet-Sun observations in 2009 (updated Fig. 9 of Paper I). Some of these curves differ from those plotted in Paper I, due to an error in the SuFI straylight removal in that paper.
![]() |
Fig. 11. All binned contrast vs. BLOS curves from Figs. 3, 4, 7 and 8 of Paper I plotted together. Also plotted is the contrast of the 5250 Å line core obtained by averaging the intensities at the wavelength positions +40 and −40 mÅ from line core. The curves are identified by their color in the legend in the upper left part of the figure. The black dashed line marks the mean quiet-Sun intensity level, i.e., a contrast of unity. |
In the following, we describe the main differences to the plot in Paper I. Firstly, the QS Ca II H (black) reaches higher contrast than the line core of IMaX (green). This is basically due to the fact that the Ca II H contrast (SuFI observations) increased relative to Paper I, while the IMaX line-core contrast remained unchanged. Secondly, the contrast at 313 nm (pink) is still lower than 388 nm (yellow) but the difference is less pronounced than obtained in Paper I. Lastly, The contrast at 300 nm (light blue) is now higher than 388 nm (yellow) for BLOS above 600 G.
In Fig. 12, we overplot the averaged contrasts for the quiet Sun (red curves) observed with SUNRISE I and analyzed in Paper I (updated and consistent with the curves shown in Fig. 11) and for the AR plage (blue curves) analyzed in the current work, of the IMaX continuum (panel a) and line core (panel b) at 525 nm, the continuum of 300 nm (panel c), and the line core of Ca II H at 397 nm (panel d) against the photospheric magnetic field. The black dashed line corresponds to a contrast of 1.
![]() |
Fig. 12. Averaged contrast of IMaX continuum (panel a), IMaX line core (panel b), SuFI 300 nm (panel c), SuFI 397 nm (panel d) against BLOS for the quiet Sun observed by SUNRISE in 2009 (blue curves) and for plage (red curves) in 2013. |
In all 4 panels, the plage curves lie generally higher than the QS ones at low BLOS, but are lower at larger BLOS, with the latter effect being marginal for Ca II H. Very close to BLOS = 0, however, the QS curves are higher again for SuFI 300 nm and for IMaX continuum. For SuFI 300 nm and IMaX continuum the larger variation of the contrast curve in QS at small BLOS is likely due to the brighter granules and darker intergranular lanes than in the abnormal granulation present in the plage region (Narayan & Scharmer 2010). This “fishhook shape” (Schnerr & Spruit 2011) is extended to larger BLOS in plage due to the filling of intergranular lanes (and to some extent the sides of granules) with weak field. The greater brightness of the IMaX continuum in the QS is in accordance with the findings of Title et al. (1992), Lawrence et al. (1993) and Kobel et al. (2011). Whereas the transition to higher QS contrast in magnetic features occurs at 500 G for the IMaX continuum, it occurs already at around 100 G for the IMaX line core and the SuFI 300 nm channel. In addition, plage magnetic elements at 300 nm reach a contrast > 1 at around 250 G, while QS magnetic elements make this transition already at 150 G. The brightness in the Ca II H line core in both, quiet-Sun network and plage is similar for BLOS above 1000 G. The slightly higher QS contrast above approximately 1400 G is based on too few QS points to be truly significant.
In the literature, the higher contrasts at larger BLOS values reached by the visible continuum in the QS than in AR plages was attributed to the higher efficiency of convective energy transport in the quiet Sun. Going from QS to AR, i.e., with increasing magnetic flux averaged over the region, convection is suppressed in the surroundings of the stronger and larger magnetic features (ARs) leading to less efficient heating from the convective walls, and therefore less brightening with respect to their counterparts in the quiet Sun (Vögler 2005; Morinaga et al. 2008; Kobel et al. 2012; Criscuoli 2013; Riethmüller et al. 2014). The same reasoning may also explain the significantly lower contrast in plage at larger BLOS seen at 300 nm, given the similar formation height to the IMaX continuum.
In the cores of spectral lines (Figs. 12b and d), this difference is less significant compared to features formed in lower layers. This result is in accordance with the temperature models of network and plage flux tubes derived from the FTS observations by Solanki (1996), cf. Solanki (1987) and Solanki & Brigljevic (1992). In these works, flux tubes in the network were found to be hotter than in plage, especially in lower regions of the photosphere.
3.6. Comparison between 2009 and 2013 quiet-Sun data
As mentioned in Sect. 2.2, we can identify some quiet-Sun areas in IMaX scans (see for example red box in Fig. 1). We also delimited a couple of quiet-Sun areas in the upper right side of the FOV. We produced contrast–BLOS scatterplots of the pixels embedded in these boxes to test if we can reproduce the relationship obtained from observations of the solar surface when the Sun was mostly quiet, i.e SUNRISE observations from 2009. For brevity, we refer to QS data from 2009 as QS-2009 and QS boxes from scans recorded in 2013 as QS-2013. Since the quiet-Sun boxes in 2013 data are outside the SuFI FOV, we restrict this comparison to the IMaX continuum and line-core contrasts vs. BLOS.
In Fig. 13 we show the scatterplots of the binned IMaX continuum (panel a) and line-core (panel b) contrasts, in blue for QS-2009 and in red for QS-2013, vs. BLOS. The black dashed horizontal line is the mean QS intensity level. Due to the smaller number of data points in QS-2013, the corresponding curves of averaged contrasts look noisier. Hence, we use non-parametric smoothing on the scattered data points of QS-2013 to make the trend easy for the eye to follow. The smoothed curves are overplotted in yellow in Fig. 13.
![]() |
Fig. 13. IMaX continuum (panel a) and line-core (panel b) contrasts vs. BLOS. The red curves are the binned contrasts of quiet-Sun areas extracted from the 2013 ARs scans, and the blue curves are taken from Paper I for quiet-Sun observations in 2009. The yellow curves are the result of non-parametric smoothing applied to the scattered data points of quiet-Sun pixels of 2013. The black dashed line is the quiet-Sun level where contrast is equal to unity. |
According to Fig. 13, contrast–BLOS relationships in both, IMaX continuum and line core agree qualitatively in the quiet-Sun regions from both flights. In Paper I we came to the conclusion that magnetic fields in the quiet-Sun network at disk center are resolved by IMaX since the contrast in the visible continuum saturated at higher BLOS. The fact that the quiet-Sun features embedded in AR scans exhibit the same behavior implies that such fields are also resolved in SUNRISE II observations of active regions. Not surprisingly, the behavior of the averaged contrast variation with BLOS of magnetic elements in the quiet Sun is independent of solar magnetic activity level (Ortiz et al. 2006), so that any difference between the two flights would reflect instrumental differences.
Quantitatively, in panel a of Fig. 13, the minumum contrast of the fishhook in QS-2013 occurs at around 150 G, a value higer than the 80 G obtained for QS-2009. In addition, the features in QS-2013 start to become brighter than the mean QS at 650 G, i.e., at fields higher than the 500 G for QS-2009.
The averaged contrast of magnetic features in the QS-2009 data is higher than in the QS-2013 in both wavelength bands. Considering the basic picture of the “hot wall” and that observations in 2013 were done further from disk center (μ = 0.93) than in 2009 (μ = 0.98), we expect that magnetic elements are brighter (in the visible continuum) with increasing heliocentric angle, but the opposite behavior is seen in Fig. 13a. However, the difference between the two data sets is small. The line core of IMaX (Fig. 13b) displays the same behavior, which in this case is consistent with the findings of Yeo et al. (2013). The differences in the contrast could be interpeted as the effect of the higher average magnetic flux density in AR scans of 2013 data (50 G) compared to the quiet Sun observed by SUNRISE in 2009 (30 G). The convection is more hampered during the period of higher activity (2013), which leads to less effective heating inside the magnetic elements from their surroundings, and therefore less radiative energy emitted by these elements.
4. Discussion and conclusions
We have analyzed the properties of solar AR plage near disk center, as observed by SUNRISE during its second flight on June 2013. In particular, we have qualitatively and quantitatively described how the brightness in the visible and UV of the different components of magnetic features depends on the longitudinal field strength, BLOS. The latter is computed among other physical parameters from SPINOR inversions of IMaX Stokes profiles corrected for wavefront aberrations and straylight.
At the spatial resolution of 0.15″ achieved by IMaX, the scatterplot of the continuum contrast at 525 nm vs. BLOS (Fig. 3) shows that the contrast peaks at around 850 G beyond which it decreases with increasing BLOS, a behavior that agrees qualitatively with the findings of Berger et al. (2007), Narayan & Scharmer (2010) and Kobel et al. (2011). The findings of the latter citation were interpreted to be an effect of the finite spatial resolution of the data (Danilovic et al. 2013).
After locating the pixels in the BLOS range of 600 G-2000 G, we found that these pixels belong to the same features, which are well resolved by SUNRISE (see Fig. 9). The BLOS is always high in their centers and decreases toward their edges, while the contrast displays a more complex behavior that can be different from one feature to another, depending on its size. In most of the features, however, the contrast is > 1 in the limb direction, < 1 toward disk center, and neutral near their centers (Fig. 10d), whereas for larger size features, the centers are dark and brighten toward their edges (Figs. 10a, b). These complex profiles are attributed to the fact that plage features are large and they are located slightly off-disk-center in the analyzed data (Steiner 2005).
Consequently, one and the same plage magnetic element provides pixels in different parts of the contrast vs. BLOS plot, including bright and dark intermediate field pixels, and dark strong-field pixels. As a result, the turnover of brightness with BLOS is not an effect of image smearing induced by the limited spatial resolution but rather of the internal morphology of individual plage magnetic features. This conclusion is strengthened by the saturation of the continuun contrast at higher BLOS in quiet-Sun magnetic features (Paper I and Sect. 3.6), which are smaller than the features in the plage, so that their contrasts are more likely to suffer from insufficient spatial resolution. The network features also show much less internal structure.
We think that the simulations used so far to interpret the downturn of the contrast in observations as an artifact of the poor spatial resolution (Röhrbein et al. 2011; Danilovic et al. 2013), lack the larger magnetic features found in strong plage that dominate their continuum contrasts. By introducing higher average magnetic flux (≈400 G) and by extending the vertical depth, we expect the boxes to be deep enough to allow convection to produce larger features. In addition, the current simulations for plage regions are analyzed for lines-of-sight parallel to the solar normal (i.e., μ = 1). Therefore, we propose that for such simulations to reproduce our data, bigger and deeper simulation boxes are needed, in addition to performing spectral line synthesis with such simulations for lines-of-sight that are inclined to the local normal (μ = 0.93).
In the UV at 397 nm (core of the chromospheric Ca II H line) and in the core of the IMaX line, the averaged contrasts of plage features increase with increasing BLOS. Given the larger formation heights of these wavelengths compared to the IMaX continuum, this implies that the brightness at larger atmospheric heights is independent of the feature size, while lower in the photosphere, the size of the magnetic features plays a decisive role in determining their contrast (Solanki et al. 2001).
Comparison of the contrast–BLOS relationship in plage with quiet-Sun observations presented in Paper I confirms the findings of earlier studies: at the photospheric level and for larger field strengths, the contrast in the quiet Sun network is higher than active region plages (Title et al. 1992; Morinaga et al. 2008; Kobel et al. 2011) due to the less efficient convective energy transport in the latter (Morinaga et al. 2008; Kobel et al. 2012; Criscuoli 2013).
In the cores of spectral lines, this difference is found to be smaller. At the atmospheric layers sampled by these wavelengths (upper photosphere and lower chromosphere), radiative heating from convection becomes less efficient in determining the contrast of these features, and other processes dominate in transporting and dissipating energy (for example oscillations and waves).
Acknowledgments
The German contribution to SUNRISE and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The Spanish contribution was funded by the Ministerio de Economía y Competitividad under Projects ESP2013-47349-C6 and ESP2014-56169-C6, partially using European FEDER funds. The HAO contribution was partly funded through NASA grant number NNX13AE95G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea.
References
- Barthol, P., Gandorfer, A., Solanki, S. K., et al. 2011, Sol. Phys., 268, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Berger, T. E., van der Rouppe Voort, L. H. M., & Löfdahl, M. G. 2004, A&A, 428, 613 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Berger, T. E., Rouppe van der Voort, L., & Löfdahl, M. 2007, ApJ, 661, 1272 [NASA ADS] [CrossRef] [Google Scholar]
- Berkefeld, T., Schmidt, W., Soltau, D., et al. 2011, Sol. Phys., 268, 103 [NASA ADS] [CrossRef] [Google Scholar]
- Buehler, D., Lagg, A., van Noort, M., & Solanki, S. K. 2018, A&A, submitted [Google Scholar]
- Chatzistergos, T. 2017, Analysis of Historical Solar Observations and Long-term Changes in Solar Irradiance (Uni-edition) (Göttingen: Georg-August-Universität) [Google Scholar]
- Chitta, L. P., Peter, H., Solanki, S. K., et al. 2017, ApJS, 229, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Criscuoli, S. 2013, ApJ, 778, 27 [NASA ADS] [CrossRef] [Google Scholar]
- Danilovic, S., Röhrbein, D., Cameron, R. H., & Schüssler, M. 2013, A&A, 550, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Frazier, E. N. 1971, Sol. Phys., 21, 42 [NASA ADS] [CrossRef] [Google Scholar]
- Frutiger, C., Solanki, S. K., Fligge, M., & Bruls, J. H. M. J. 2000, A&A, 358, 1109 [NASA ADS] [Google Scholar]
- Gandorfer, A., Grauf, B., Barthol, P., et al. 2011, Sol. Phys., 268, 35 [NASA ADS] [CrossRef] [Google Scholar]
- Hirzberger, J., & Wiehr, E. 2005, A&A, 438, 1059 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ishikawa, R., & Tsuneta, S. 2009, A&A, 495, 607 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jafarzadeh, S., Solanki, S. K., Stangalini, M., et al. 2017a, ApJS, 229, 10 [NASA ADS] [CrossRef] [Google Scholar]
- Jafarzadeh, S., Rutten, R. J., Solanki, S. K., et al. 2017b, ApJS, 229, 11 [NASA ADS] [CrossRef] [Google Scholar]
- Kahil, F., Riethmüller, T. L., & Solanki, S. K. 2017, ApJS, 229, 12 [NASA ADS] [CrossRef] [Google Scholar]
- Keller, C. U., Schüssler, M., Vögler, A., & Zakharov, V. 2004, ApJ, 607, L59 [NASA ADS] [CrossRef] [Google Scholar]
- Kobel, P., Solanki, S. K., & Borrero, J. M. 2011, A&A, 531, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kobel, P., Solanki, S. K., & Borrero, J. M. 2012, A&A, 542, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kostik, R., & Khomenko, E. V. 2012, A&A, 545, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Krivova, N., Solanki, S. K., Fligge, M., & Unruh, C. Y. 2003, A&A, 399, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lawrence, J. K., Topka, K. P., & Jones, H. P. 1993, J. Geophys. Res., 98, 18 [Google Scholar]
- Loukitcheva, M., Solanki, S. K., & White, S. M. 2009, A&A, 497, 273 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martínez Pillet, V., del Toro Iniesta, J. C., Álvarez-Herrero, A., et al. 2011, Sol. Phys., 268, 57 [NASA ADS] [CrossRef] [Google Scholar]
- Morinaga, S., Sakurai, T., Ichimoto, K., et al. 2008, A&A, 481, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Narayan, G., & Scharmer, G. B. 2010, A&A, 524, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ortiz, A., & Rast, M. 2005, Mem. Soc. Astron. It., 76, 1018 [Google Scholar]
- Ortiz, A., Solanki, S. K., Domingo, V., Fligge, M., & Sanahuja, B. 2002, A&A, 388, 1036 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ortiz, A., Domingo, V., & Sanahuja, B. 2006, A&A, 452, 311 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Parker, E. N. 1963, ApJ, 138, 552 [NASA ADS] [CrossRef] [Google Scholar]
- Pickering, K. A. 2002, DIO, 12, 3 [NASA ADS] [Google Scholar]
- Rezaei, R., Schlichenmaier, R., Beck, C. A. R., Bruls, J. H. M. J., & Schmidt, W. 2007, A&A, 466, 1131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Riethmüller, T. L., Solanki, S. K., Pillet, V. M., et al. 2010, ApJ, 723, L169 [NASA ADS] [CrossRef] [Google Scholar]
- Riethmüller, T. L., Solanki, S. K., Berdyugina, S. V., et al. 2014, A&A, 568, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Riethmüller, T. L., Solanki, S. K., Barthol, P., et al. 2017, ApJS, 229, 16 [NASA ADS] [CrossRef] [Google Scholar]
- Röhrbein, D., Cameron, R., & Schüssler, M. 2011, A&A, 532, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Scharmer, G. B., Bjelksjo, K., Korhonen, T. K., Lindberg, B., & Petterson, B. 2003, Proc. SPIE, 4853, 341 [Google Scholar]
- Schnerr, R. S., & Spruit, H. C. 2011, A&A, 532, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schrijver, C. J., Cote, J., Zwaan, C., & Saar, S. H. 1989, ApJ, 337, 964 [Google Scholar]
- Skumanich, A., Smythe, C., & Frazier, E. N. 1975, ApJ, 200, 747 [NASA ADS] [CrossRef] [Google Scholar]
- Solanki, S. K. 1987, PhD Thesis, No. 8309, ETH Zürich, Switzerland [Google Scholar]
- Solanki, S. K. 1993, Space Sci. Rev., 63, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Solanki, S. K. 1996, IAU Symp., 176, 201 [NASA ADS] [Google Scholar]
- Solanki, S. K. 2001, in Magnetic Fields Across the Hertzsprung-Russell Diagram, eds. G. Mathys, S. K. Solanki, & D. T. Wickramasinghe, ASP Conf. Ser., 248, 45 [NASA ADS] [Google Scholar]
- Solanki, S. K., & Brigljevic, V. 1992, A&A, 262, L29 [NASA ADS] [Google Scholar]
- Solanki, S. K., Barthol, P., Danilovic, S., et al. 2010, ApJ, 723, L127 [NASA ADS] [CrossRef] [Google Scholar]
- Solanki, S. K., Riethmüller, T. L., Barthol, P., et al. 2017, ApJS, 229, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Spruit, H. C. 1976, Sol. Phys., 50, 269 [NASA ADS] [CrossRef] [Google Scholar]
- Steiner, O. 2005, ESA SP, 600, 1 [Google Scholar]
- Stenflo, J. O. 1973, Sol. Phys., 32, 41 [NASA ADS] [CrossRef] [Google Scholar]
- Title, A. M., Topka, K. P., Tarbell, T. D., et al. 1992, ApJ, 393, 782 [NASA ADS] [CrossRef] [Google Scholar]
- Topka, K. P., Tarbell, T. D., & Title, A. M. 1992, ApJ, 396, 351 [NASA ADS] [CrossRef] [Google Scholar]
- van Noort, M., Rouppe van der Voort, L., & Löfdahl, M. G. 2005, Sol. Phys., 228, 191 [NASA ADS] [CrossRef] [Google Scholar]
- Vögler, A. 2005, Mem. Soc. Astron. It., 76, 842 [NASA ADS] [Google Scholar]
- Vögler, A., Shelyag, S., Schüssler, M., et al. 2005, A&A, 429, 335 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Yeo, K. L., Solanki, S. K., & Krivova, N. A. 2013, A&A, 550, A95 [Google Scholar]
- Yeo, K. L., Solanki, S. K., Norris, C. M., et al. 2017, Phys. Rev. Lett., 119, 091102 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, J., Lin, G., Wang, J., Wang, H., & Zirin, H. 1998, Sol. Phys., 178, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Zhou, G. P., Wang, J. X., & Jin, C. L. 2010, Sol. Phys., 267, 63 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: The quiet Sun’s mean intensity in SuFI data
We use the dark corrected (level 0.1) data taken over nearly 4 hours (from June 12, 22:15:37 UT to June 13, 01:58:00 UT) to visualize the day-to-night cycles of the mean quiet-Sun intensity at 300 nm and 397 nm. Quiet Sun images are acquired for more than an hour at the beginning of the series, just prior to the 1 h spent observing the active region, and in three intervals after the end of the AR observations.
We plot for both wavelengths the averaged level 0.1 QS intensity (I) versus time (t), with the time at which the first image of the time series was recorded being referred to as t = 0. The day-night cycle is clearly visible in the left panel of Fig. A.1.
For fitting the day-to-night cycles we use the Beer-Lambert law, which gives the amount of flux (of original value I0) absorbed by a medium with a Rayleigh airmass of m:
τλ is the optical depth of the terrestrial atmosphere at a given wavelength λ. The Rayleigh airmass depends on the Sun’s elevation angle, ϕ (Pickering 2002), and for a spherically-symmetric atmosphere it is given by:
For each wavelength, I0 and τλ are determined from a fit to the I vs. m data based on Eq. (A.1). They are shown in Table A1. Plots of I versus m and the corresponding fits are shown in Figs. A.2 and A.3 for 300 nm and 397 nm, respectively. As expected, the day-to-night cycle variation is negligible at 397 nm since the average quiet-Sun intensity does not vary with time (or elevation angle). Therefore, the average of the flatfields is the average QS intensity.
At 300 nm, we use the best-fit parameters to evaluate I at the airmass values of our AR data. The I values determined in this way are the true mean quiet-Sun intensities. We normalize each image in our level 3.1 data with the corresponding evaluated mean after restoring the original flux. The latter is simply the product of every level 3.1 pixel value with the averaged flatfield to which our data were normalized (assuming that the straylight-correction and phase-diversity reconstruction do not affect the mean value of the image).
Best-fit parameters at each wavelength obtained by using Eq. (A.1) and applied to I vs. m scatterplots.
![]() |
Fig. A.1. Part of the day-to-night cycle of the photon flux at 300 nm and 397 nm for quiet-Sun images taken at disk center. The dashed blue lines delimit the 1 h observation of the active region, which was partly analyzed in this work. |
![]() |
Fig. A.2. Red data points: the 300 nm quiet-Sun photon flux vs. m, the air mass factor during the 4 h observing period. The best-fit blue curve was calculated according to Eq. (A.1) with the parameters given in Table A.1. |
All Tables
Parameters of logarithmic fits to the contrast in the core of Fe I line at 5250.2 Å (starting from 120 G) and for 300 nm and 397 nm (starting from 100 G) vs. BLOS.
Best-fit parameters at each wavelength obtained by using Eq. (A.1) and applied to I vs. m scatterplots.
All Figures
![]() |
Fig. 1. Left panel: continuum contrast at 525 nm. Right panel: its co-spatial and co-temporal longitudinal magnetic field map retrieved from the inversions. The black boxes (solid lines) contain small pores characterized by contrast below unity and kG magnetic fields. The dashed green box encloses an area with emerging flux. The dashed yellow ellipse outlines a region of plage composed of magnetic elements embedded in intergranular lanes. The dashed red box contains a quiet-Sun internetwork region with weak fields (average of 20 G) and mean contrast of unity (by definition). |
In the text |
![]() |
Fig. 2. SuFI images coaligned to the IMaX FOV shown in Fig. 1. Panel a: SuFI 300 nm. Panel b: SuFI Ca II H. Panel c: IMaX continuum intensity (normalized to the mean quiet-Sun intensity) to which SuFI 300 nm (overlaid in red) is aligned. Panel d: IMaX line-core intensity derived from the inversions (see Sect. 2.5) and normalized to the local continuum (I+227), and to which calcium images (overlaid in yellow) are aligned. The yellow and red boxes are slightly misaligned due to the differential offsets in the SuFI wavelegnth channels. The gray scale is set to cover two times the rms range of each image. |
In the text |
![]() |
Fig. 3. Scatterplot of the continuum contrast of Fe I 5250.4 Å vs. the LOS component of the magnetic field, BLOS, for the AR plage. Pores and the FE region are excluded. The horizontal dashed red line is the mean QS continuum intensity level (the mean contrast of the pixels in the dashed red box in Fig. 1). The red curve joins the average contrast values in bins of 500 data points each. |
In the text |
![]() |
Fig. 4. Same as Fig. 3 but for the line core of Fe I 5250.4 Å. The blue curve is the logarithmic fit to the averaged line-core contrast values starting at 120 G. |
In the text |
![]() |
Fig. 5. Scatterplot of the SuFI intensity contrast at 300 nm vs. BLOS. The red curve joins the averaged contrast values inside each bin, the blue curve is the logarithmic fit to the red curve starting from 100 G. |
In the text |
![]() |
Fig. 6. Same as Fig. 5 but for the line core of Ca II H. The blue curve is the logarithmic fit to the red curve starting from 50 G. |
In the text |
![]() |
Fig. 7. Intensity contrasts of plage vs. BLOS in the wavelengths considered in the current study (see the legend in the figure). Plotted are the contrast values (plotted in red in Figs. 3–6). For comparison, also plotted is the binned line-core contrast at 525.02 nm derived from Gaussian fits to individual Stokes I profiles. |
In the text |
![]() |
Fig. 8. Contrast histograms of pixels belonging to the common FOV of SuFI and IMaX after masking out pores. Colors correspond to different wavelengths, and panels correspond to different ranges of BLOS values (see text for a detailed explanation). The histograms are normalized such that the integral over each is equal to one. |
In the text |
![]() |
Fig. 9. Panels a and b: IMaX data for the same time as in Fig. 1, now with the absolute value of BLOS plotted in the right panel. The red contours enclose the “peak” pixels with > 1 averaged contrast in the continuum of 525 nm and 600 G < BLOS < 1000 G, while the yellow contours enclose the “turnover” pixels with averaged contrasts < 1 and BLOS larger than 1200 G. Panels c, d, e: blowups of three magnetic features of different sizes. Red and green contours enclose the bright (contrast > 1) parts of both families of pixels, while orange and yellow contours enclose the dark (contrast < 1) parts. |
In the text |
![]() |
Fig. 10. Profiles of BLOS (green curve) and continuum intensity contrast at 525 nm (blue curve) along cuts through four different magnetic structures in an example IMaX image (see text for details). The red dashed curves are the Gaussian fits to the BLOS profiles. The FWHM of these Gaussian fits are given at the top of each frame. The horizontal dashed blue line is where the contrast is equal to unity. The limb is to the right of the figure, while disk center is to the left |
In the text |
![]() |
Fig. 11. All binned contrast vs. BLOS curves from Figs. 3, 4, 7 and 8 of Paper I plotted together. Also plotted is the contrast of the 5250 Å line core obtained by averaging the intensities at the wavelength positions +40 and −40 mÅ from line core. The curves are identified by their color in the legend in the upper left part of the figure. The black dashed line marks the mean quiet-Sun intensity level, i.e., a contrast of unity. |
In the text |
![]() |
Fig. 12. Averaged contrast of IMaX continuum (panel a), IMaX line core (panel b), SuFI 300 nm (panel c), SuFI 397 nm (panel d) against BLOS for the quiet Sun observed by SUNRISE in 2009 (blue curves) and for plage (red curves) in 2013. |
In the text |
![]() |
Fig. 13. IMaX continuum (panel a) and line-core (panel b) contrasts vs. BLOS. The red curves are the binned contrasts of quiet-Sun areas extracted from the 2013 ARs scans, and the blue curves are taken from Paper I for quiet-Sun observations in 2009. The yellow curves are the result of non-parametric smoothing applied to the scattered data points of quiet-Sun pixels of 2013. The black dashed line is the quiet-Sun level where contrast is equal to unity. |
In the text |
![]() |
Fig. A.1. Part of the day-to-night cycle of the photon flux at 300 nm and 397 nm for quiet-Sun images taken at disk center. The dashed blue lines delimit the 1 h observation of the active region, which was partly analyzed in this work. |
In the text |
![]() |
Fig. A.2. Red data points: the 300 nm quiet-Sun photon flux vs. m, the air mass factor during the 4 h observing period. The best-fit blue curve was calculated according to Eq. (A.1) with the parameters given in Table A.1. |
In the text |
![]() |
Fig. A.3. Same as Fig. A.2, but for the 397 nm data. |
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.