Issue |
A&A
Volume 675, July 2023
|
|
---|---|---|
Article Number | A93 | |
Number of page(s) | 14 | |
Section | The Sun and the Heliosphere | |
DOI | https://doi.org/10.1051/0004-6361/202345890 | |
Published online | 04 July 2023 |
A study of the capabilities for inferring atmospheric information from high-spatial-resolution simulations
1
Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
e-mail: carlos.quintero@iac.es
2
Departamento de Astrofísica, Univ. de La Laguna, La Laguna, Tenerife, 38205
Spain
3
Instituto de Astrofísica e Ciências do Espaço, Departamento de Física, Universidade de Coimbra, Rua do Observatório s/n, 3040-004 Coimbra, Portugal
4
High Altitude Observatory, NCAR, PO Box 3000 Boulder, CO, 80307
USA
5
Astrophysics Research Centre, School of Mathematics and Physics, Queen’s University Belfast, Belfast, BT7 1NN Northern Ireland, UK
6
Institute for Solar Physics, Department of Astronomy, Stockholm University, AlbaNova University Centre, 10691 Stockholm, Sweden
7
National Solar Observatory, University of Colorado Boulder, 3665 Discovery Drive, Boulder, CO, 80303
USA
8
Instituto de Astrofísica de Andalucía (CSIC), Apdo. de Correos 3004, 18080 Granada, Spain
Received:
11
January
2023
Accepted:
19
May
2023
In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magneto-hydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4m class telescopes DKIST and EST. We first study multiple inversion configurations using the SIR code and the Fe I transitions at 630 nm until we obtain minor differences between the input and the inferred atmosphere in a wide range of heights. Also, we examine how the inversion accuracy depends on the noise level of the Stokes profiles. The results indicate that when the majority of the inverted pixels come from strongly magnetised areas, there are almost no restrictions in terms of the noise, obtaining good results for noise amplitudes up to 1 × 10−3 of Ic. At the same time, the situation is different for observations where the dominant magnetic structures are weak, and noise restraints are more demanding. Moreover, we find that the accuracy of the fits is almost the same as that obtained without noise when the noise levels are on the order of 1 × 10−4of Ic. We, therefore, advise aiming for noise values on the order of or lower than 5 × 10−4 of Ic if observers seek reliable interpretations of the results for the magnetic field vector reliably. We expect those noise levels to be achievable by next-generation 4m class telescopes thanks to an optimised polarisation calibration and the large collecting area of the primary mirror.
Key words: Sun: magnetic fields / magnetohydrodynamics (MHD) / techniques: polarimetric / radiative transfer / techniques: high angular resolution
© The Authors 2023
Open 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
Understanding the physics of solar phenomena implies the inference of information about their thermal and magnetic properties. In the case of the magnetic field, there are multiple ways to infer that information, from basic techniques such as deriving the field strength from the separation between the Stokes I Zeeman components (e.g., Landi Degl’Innocenti & Landolfi 2004), comparing the derivative of the Stokes I profile with Stokes V with the so-called weak-field approximation (see de la Cruz Rodríguez & van Noort 2017 or Centeno 2018 for a recent example), to more elaborate techniques such as the so-called inversion codes (see del Toro Iniesta & Ruiz Cobo 2016, for an introduction). Regarding the estimation of the magnetic properties of quiet Sun regions using inversion techniques, there were multiple works in the early 2000s. Those works used observations made with ground-based telescopes, for instance, with the Tenerife Infrared Spectropolarimeter (Collados et al. 2007) installed at the Vacuum Tower Telescope (von der Lühe 1998) and the Advanced Stokes Polarimeter (Elmore et al. 1992) installed in the Dunn Solar Telescope. The observations usually recorded the visible Fe I line pair at 630 nm or the infrared Fe I lines at 1.56 microns. The results of their analysis were published in multiple publications (see, among others, Lites 2002; Khomenko et al. 2003; Lites & Socas-Navarro 2004; Sánchez Almeida 2004; Domínguez Cerdeña et al. 2006).
After that period, the solar community saw an exponential increase in the number of works dedicated to elucidating the magnetic properties of the quiet Sun, mainly because the Hinode mission (Kosugi et al. 2007) was launched in 2006, in a period of minimum solar activity. The spacecraft offered high-stability observations (Shimizu et al. 2008) over long periods thanks to the excellent performance of the Solar Optical Telescope (Tsuneta et al. 2008). This stability allowed multiple types of quiet Sun observations to be made with short and long exposure times, over small and large fields of view, and at the disk centre and the limb and so on using the spectropolarimeter (Lites et al. 2013) installed on the optical telescope. Later, those data were complemented with higher spatial resolution observations taken with the Imaging Magnetograph eXperiment (Martínez Pillet et al. 2011) on board the Sunrise balloon (Solanki et al. 2010), which was launched for the first time in 2009. There is a long list of publications from quiet Sun observations made using the Hinode satellite and the Sunrise balloon, and so with the best intention of not leaving any critical work behind, we cite the recent review of Bellot Rubio & Orozco Suárez (2019), where those works are gathered and thoroughly explained. However, we cite a few particular publications that we consider relevant to the present article. One of the early works from the Hinode era was that of Lites et al. (2008), where the authors provided statistics regarding the properties of the magnetic field in the quiet Sun. This was followed by a plethora of related publications, with some works building on the dynamics of the quiet Sun magnetic fields (for instance, Centeno et al. 2007; Martínez González & Bellot Rubio 2009; Danilovic et al. 2010).
In parallel, we also saw a rise in the number of works focusing on the possible limitations of studying weak polarisation signals that are dominated by noise using synthetic data generated from semi-empirical atmospheres as well as from realistic 3D magneto-hydrodynamic (MHD) simulations (Borrero & Kobel 2011; de la Cruz Rodríguez et al. 2012; Danilovic et al. 2016; Beck et al. 2017; Milić et al. 2019; Campbell et al. 2021, among others,). Those works tested the ability to retrieve the original atmosphere used as input for the synthesis of the Stokes profiles under different conditions of spatial and spectral degradation, or noise. We believe these works are an excellent reference for observers that aim to, for example, prepare a proposal for observations, and we want to expand their work further. We aim to reproduce observations that will soon be made by next-generation 4 m class telescopes such as DKIST (Rimmele et al. 2020) and EST (Quintero Noda et al. 2022). To this end, we use state-of-the-art high-spatial-resolution numerical simulations with a pixel element with the size expected at the diffraction limit of those telescopes. Also, we aim to set the basis for highly accurate inversions where the user of inversion codes seeks to infer the stratification of the atmospheric parameters at different layers, obtaining precise fits of complex Stokes profiles through the use of gradients along the line of sight (LOS) on temperature, velocity, and the magnetic field vector. Thus, we aim to provide an inversion strategy from which future observers at the aforementioned next-generation telescopes can start from. Moreover, we are curious to see how the inversion technique behaves in such complex numerical scenarios when different noise levels are included and try to set a reference noise level that observers would need to aim for depending on their research interests.
2. Data and methodology
We use two different 3D snapshots in the present study. The first was computed with the MANCHA code (Khomenko & Collados 2006; Felipe et al. 2010; González-Morales et al. 2018) with the configuration described in Khomenko et al. (2017) and González-Morales et al. (2018). The realistic convection simulation snapshot used here is a simulation run initiated as purely hydrodynamic convection in a box spanning from approximately 1 Mm below the photosphere to 600 km above it. The pixel size in the horizontal and vertical directions is 10 and 7 km, respectively, while the snapshot has a squared field of view (FOV) of 15.88×15.88 arcsec2 (a total of 1152×1152 pix2). Once the convection has reached a stationary state, an initial magnetic field was seeded by ‘switching on’ the Biermann battery term in the induction equation (Khomenko et al. 2017). Its continuous action, together with the dynamo amplification, provided the spatially averaged magnetisation of the model of about ⟨|B|⟩ ≈ 102 G at the surface where the continuum optical depth at 500 nm is unity after a few hours of solar time. The snapshot we use in this work corresponds to the saturated dynamo regime.
The second snapshot corresponds to the ‘Sunspot’ run from Rempel (2012) computed with the MURaM code (Vögler et al. 2005). The pixel size in the horizontal and vertical domains is 12 and 8 km, respectively. We focused only on the bottom left corner of the original field of view. This small subset contains 1000×1000 pix2 and has a total size of 16.55×16.55 arcsec2.
The original atmospheres, computed on a regular vertical geometrical grid, were transformed to a 3D cube where the vertical domain corresponds to the continuum optical depth at 500 nm (log τ, from now on). For that purpose, we used the routine modelador.x included in the Stokes Inversion based on Response functions (SIR, Ruiz Cobo & del Toro Iniesta 1992) code. Thus, all the plots in this work correspond to the atmospheric parameters at a constant optical depth (while the geometrical height may vary from pixel to pixel). Figure 1 shows the spatial distribution of the temperature for the two snapshots at a height where log τ = 0. We also highlight a smaller region with a coloured square in each panel which corresponds to the spatial domain we use for optimising the inversion configuration (orange) we use as a reference setup in the following sections or a region we study in detail in this work (blue).
Fig. 1. Spatial distribution of the temperature where the continuum optical depth at 500 nm is unity. The top panel shows the results from a quiet-Sun run with 10 km spatial resolution computed with the MANCHA code. The bottom panel corresponds to a cut of the original Sunspot run with a 12 km spatial resolution generated with the MURaM code. The regions highlighted with coloured boxes designate the areas of interest on which we focus in this work. |
We mention above that the aim of this analysis is to understand our capabilities for inferring atmospheric information from observations made by future-generation 4m-class telescopes, such as DKIST and EST. Both telescopes will have access to the 630 nm region in the visible, which is traditionally used by many telescopes. Therefore, in this first work, we focus only on this spectral region, including the Fe I transitions at 6301.5 and 6302.5 Å while additional spectral lines will be studied in the future. Furthermore, we use the atomic parameters presented in Quintero Noda et al. (2021), and we do not add any spatial degradation as the pixel size at 630 nm for a 4m telescope is comparable to that of the simulations used in this work:
The synthesis of the Stokes profiles is done with SIR with a spectral sampling of 10 mÅ/pix over a spectral range that goes from −500 to 1500 mÅ from the rest wavelength of the line core of Fe I 6301.5 Å. The conditions mentioned above are easily achievable with the Visible Spectro-Polarimeter (de Wijn et al. 2022) installed in DKIST, for example. We do not add any spectral degradation resulting from the finite spectral resolution of an instrument; for example, we do not convolve the spectrum by a point spread function of a given width. Furthermore, we do not consider additional instrumental effects like those produced by long integration times, limited slit width, or limited spectral sampling in the case of filtergraph instruments. However, a noise with different amplitude values is added to the Stokes parameters later on. The atomic abundance used during the process corresponds to that presented in Asplund et al. (2009). In addition, although some authors have pointed out the necessity of including non-local thermodynamic equilibrium effects when modelling these spectral lines (e.g., Smitha et al. 2020), we are going to make the synthesis and the inversion using the same forward module. Therefore, we do not expect to have significant issues in this regard. All the inversions are done using the Python wrapper presented in Gafeira et al. (2021). We start the inversion from a modified version of the Harvard-Smithsonian reference atmosphere (HSRA, Gingerich et al. 1971), which is designed to represent the inversion of actual observations as closely as possible; that is, the guess atmosphere we start the inversion from is independent of the simulation atmospheres used during the synthesis. The changes in the HSRA atmosphere include a constant LOS velocity and magnetic field vector so that SIR can apply multiplicative perturbations on those atmospheric parameters. The temperature-electron pressure relation is the original one published in the cited publication.
2.1. General properties of the simulation runs
The selected snapshots differ in many aspects, and so we believe it is beneficial to cover scenarios that could resemble strongly magnetised quiet Sun areas and weakly magnetised regions. We illustrate those differences in Fig. 2, displaying the horizontal (B ⋅ sin(γ), being γ the inclination) and vertical (B ⋅ cos(γ)) components of the magnetic field at log τ = 0. The quiet Sun run from the MANCHA code (left column) is characterised by weaker magnetic fields that populate most of the FOV with an intricate pattern. There are also areas with almost zero field strength that can represent a challenge for the inversion process. In the case of the MURaM Sunspot run (rightmost column), there are strongly magnetised areas with field strengths of over 2000 G. We believe they are an extension of the nearby sunspot. However, at the same time, those strong fields are surrounded by weak magnetic fields. Thus, the selected region resembles a plage region or small pore observation. We assume that fitting the Stokes profiles from those areas void of a magnetic field will also represent a challenge for the inversion process.
Fig. 2. Spatial distribution of the horizontal (top) and longitudinal (bottom) component of the magnetic field where log τ = 0. The left panel corresponds to a quiet Sun run with 10 km spatial resolution generated with the MANCHA code, while the right panel shows a cut of the original Sunspot run with 12 km spatial resolution computed with the MURaM code. |
We complement the previous figure, adding the spatially averaged atmospheric parameters at every optical depth in Fig. 3. We computed those values for the complete FOV of the MANCHA (top) and MURaM (bottom) runs and for the small regions highlighted with coloured squares in Fig. 1. In the case of the MANCHA run, both areas seem to be relatively similar with weak magnetic fields (B < 100 G) and parallel to the solar surface (γ ∼ 90 degrees). Concerning the MURaM run, the averaged atmospheric parameters for the selected area and those for the entire FOV differ more significantly, with the atmospheric parameters of the enclosed region being more vertical with respect to the solar surface (γ ∼ 45 degrees). The enclosed region also shows a stronger magnetic field (B > 1000 G). Still, in all the cases, we can see large deviations at each optical depth, indicating that the high spatial resolution of these runs produces significant small-scale variations (from one pixel to the neighbouring one), particularly in the magnetic field vector.
Fig. 3. Evolution of the averaged atmospheric parameters with height. Top and bottom panels correspond to the MANCHA and MURaM runs presented in Fig. 1. Colours designate the averaged results over the small regions highlighted in the same figure. Bars display the standard deviation of the atmospheric parameters over the corresponding areas. |
2.2. Inversion configuration I: Nodes and inversion cycles
One of the goals of this work is to estimate the loss of accuracy we suffer when inverting spectropolarimetric data with various noise levels. It is, therefore, paramount to find an inversion configuration that provides good results where the inferred atmosphere is as close as possible to the input one. For that purpose, as the two snapshots are complex, multiple inversion tests are required. However, each snapshot has more than 1000×1000 pix2, and so it is not feasible to perform all those inversion tests with the complete FOV within reasonable time limits. In that sense, we perform all the inversion tests described in this section, only inverting a small patch of the MANCHA run with a size of 100×100 pix2 (see the orange square in Fig. 1). In addition, although several authors have confirmed over the years (among others, Borrero et al. 2013; Quintero Noda et al. 2014; Campbell et al. 2021) that using different initial atmospheres helps to improve the accuracy of the solution and reduce the chances of the code falling into a local χ2 minimum, in this work, we run every test with one atmosphere for simplicity and to reduce computational time. Therefore, we assume that the predictions we get here may be slightly improved using different guess atmospheres, which we plan to explore in more detail in future publications.
We start the study for the selection of nodes, that is, the heights where the guess atmosphere is perturbed and where the response functions (RFs) are evaluated using the values provided in Table 1. We follow a straightforward process to make the exercise as simple as possible and provide a reference guide for future observers interested in performing node-based inversions. We start from a simple configuration running one inversion cycle –see C1– changing the temperature at the top and bottom boundaries while the LOS velocity and magnetic field vector are constant with height. However, we do not aim to use these results as we know that the simulated atmospheres have physical parameters that strongly vary with height. Those variations with optical depth produce various asymmetries of the Stokes profiles (among others, Balasubramaniam et al. 1997; Martínez Pillet et al. 1997; Viticchié & Sánchez Almeida 2011; Sainz Dalda et al. 2012). For this reason, we need more than 1 node on the LOS velocity and the magnetic field vector stratifications to fit the Stokes parameters with the highest accuracy. We therefore increase the number of nodes used for the inversion as we add more inversion cycles as described in Table 1. In each new iteration, we compare the input and the inferred atmosphere to determine whether or not we are improving the results by adding more nodes. We then compute the average differences, in absolute value, over the FOV used for these tests (the orange region in Fig. 1) at each optical depth (the step is every log τ = 0.1) to estimate the general accuracy of the fits. The results are shown in Fig. 4.
Fig. 4. Average differences between the input and the inferred atmosphere at different optical depths. From left to right, we show the differences in temperature, LOS velocity, magnetic field strength, and inclination, respectively. Colours correspond to the various inversion configurations indicated in Table 1. |
Number of nodes used for the inversion of each atmospheric parameter.
Starting with the results for the temperature, we can see that the first configuration, with constant perturbations along the LOS for the velocity and magnetic field, provides relatively good results at log τ = 0 and around log τ = −2 and much larger deviations between and outside those layers. We believe this is due to a lack of freedom, as we can see improvements in subsequent inversion cycles until we reach five nodes in C3 and C4. Interestingly, there is an improvement in C4, although the nodes for temperature are the same as in the previous inversion cycle. We believe it is because the rest of the atmospheric parameters have more freedom, and the fits are consequently better. In the case of the LOS velocity, we also see a continuous improvement from C1 to C4. For the field strength, we again see that the differences become smaller as we move from C1 to C4, with an average difference of less than 30 G for the latter configuration. Finally, for the magnetic field inclination, we see improvements again as we move from C1 to C4, with average differences of lower than 15 degrees in the latter case.
After examining the previous results, we chose to continue our analysis using configuration C4 as it is better than the others in all the atmospheric parameters. However, we encourage the reader to examine the results in detail, particularly the fits of the Stokes profiles, when trying to define the inversion configuration. On the one hand, we need to add gradients to fit asymmetric profiles. On the other hand, deciding whether a linear gradient (2 nodes) is sufficient or more complex parabolic (3 nodes) or polynomial shapes (more than 3 nodes) are required is not entirely straightforward without an in-depth analysis of the results. Still, we recommend that users begin experimenting with the number of nodes presented here, as we find this to be a good baseline.
Finally, we want to mention that we do not pay particular attention to the magnetic field azimuth in this work. This is mainly because the expected magnetic fields are weak in strength and are strongly local (maybe occupying a few pixels in size), and so it would be challenging to solve the 180 degree ambiguity in order to interpret the inferred azimuth values reliably (e.g., Metcalf et al. 2006). Instead, could check the input atmosphere and, based on that, add or subtract 180 degrees to match the correct orientation of the magnetic field. Also, we could compare the cosine of the azimuth (that will range from 0 to 1), removing the effect of the 180 ambiguity. However, we prefer to leave the azimuth comparison out of this work, and we will address this subject –and the different options to reliably infer it in this complex numerical scenario– in future publications. Still, we want to emphasise that the atmospheric parameters are fitted using the node values presented in Table 1 to get good fits for the Stokes Q and U profiles, and hence the most accurate atmosphere for all the physical quantities.
2.3. Inversion configuration II: Profiles and atmospheres
We mentioned in the previous section that examining the fits of a few reference profiles is advisable in order to gain confidence when choosing the inversion configuration that works best for our goals. In that sense, we show in Fig. 5 the fit for three reference pixels where the code fits the profiles relatively poorly (top), relatively well (middle), and with high accuracy. Furthermore, we show the differences between the input (black) and the inferred Stokes parameter (colours) at the bottom of each panel. We again see significant differences for C1 (blue) and C2 (orange) and much lower deviations for C3 (grey) and C4 (red). In general, the latter two configurations behave very well. Still, in a few cases, C4 produces smaller differences, and in the rest of the cases, the results are similar, indicating that C4 provides the most accurate results for our scenario. Additionally, we can investigate further, comparing the input atmosphere that produced those profiles with that obtained by the code when using each inversion configuration.
Fig. 5. Stokes profiles corresponding to the three reference pixels (rows). From left to right, Stokes (I, Q, U, V), respectively. Black corresponds to the input profiles, while colours designate the Stokes parameters corresponding to the different node configurations presented in Table 1. Each panel includes the differences between the input and the inferred profiles below the horizontal solid black line. From top to bottom, each reference profile corresponds to an intergranular lane, edge of a granule, and top of a granule, respectively. We show the exact location in Fig. 11, where we provide a zoom onto the region highlighted in the top panel of Fig. 1. |
In Fig. 6, we present the comparison for temperature, LOS velocity, and magnetic field strength and inclination for the three reference pixels shown in Fig. 5. Starting with the temperature, we have relatively smooth stratifications that are well recovered by all the configurations. In this case, we cannot see significant differences between C3 and C4. This is probably because the number of nodes is the same in both cases. Something similar happens for the LOS velocity, with both configurations showing approximately the same behaviour and more accurate results than C1 and C2. Finally, we see slightly better accuracy for C4 in the magnetic field strength and inclination, although there are no significant differences. Therefore, based on these results, combined with the previous analysis, we continue using C4. However, we would not object to another user preferring to use C3 simply because the number of nodes is lower, considering that the fit of the Stokes profiles (the only measure of accuracy when working with observations) is slightly poorer than with C4.
Fig. 6. Atmospheres corresponding to the pixels presented in Fig. 5 (different rows). We show the temperature, LOS velocity, magnetic field strength, and inclination from left to right. Black corresponds to the input profiles, while colours designate the Stokes parameters corresponding to the different node configurations presented in Table 1. |
Finally, before examining the inversion configuration, we need to comment on the general accuracy of the fits and the atmospheres we are inferring. It is clear that in some cases (Fig. 6), the code struggles to reproduce the complex atmosphere from the 3D simulation. This lack of accuracy is something we expected in the case of having atmospheres that change over short height scales. The code uses several nodes (up to five in our case) that are equidistant in optical depth (see Fig. 7 as a reference), and so it is not possible to produce such complex stratifications.
Fig. 7. Reference location of nodes in SIR. The cases corresponding to 2, 3, and 5 nodes (see Table 1) have nodes located as shown by triangles, squares and diamonds, respectively. As a reference, we include the temperature of the HSRA atmosphere. |
At the same time, it is not possible from a radiative transfer point of view to have sensitivity to such small-scale variations with our observation configuration. Spectral lines are sensitive to a given range of optical depth layers that it usually not as narrow as the changes shown by the simulation. In this case, we only have two spectral lines that are sensitive within a similar range of heights (see, e.g., Quintero Noda et al. 2021). Therefore, it is impossible to reach the degree of complexity shown by the simulation, although the code generally does relatively well.
2.4. Inversion configuration III: Iteration cycles
At this point, it may not be clear as to whether or not we need to run multiple inversion cycles, increasing the number of nodes each time, for example, as we move from C1 to C4 in Table 1, or if we could run one inversion cycle using the configuration C4. Therefore, we compare both scenarios in Fig. 8.
Fig. 8. Average differences between the input and the inferred atmosphere at different optical depths. From left to right, we show the differences for temperature, LOS velocity, magnetic field strength and inclination, respectively. Line styles correspond to the case where we reach configuration C4 in Table 1 after running the previous cycles (solid) or to cases where we run a single inversion cycle using the nodes described in column C4 of Table 1 starting from the modified HSRA atmosphere (dashed). |
The results for the case where we run four inversion cycles, gradually increasing the number of nodes in each cycle, are presented in solid black. At the same time, the dashed line shows the result with one inversion cycle using the nodes specified by C4 and starting from the modified HSRA atmosphere. We start in both cases from a modified HSRA atmosphere where LOS velocity and magnetic field vector are constant with height. In the first example, the code changes the temperature at the top and bottom boundary in C1 and applies a constant perturbation in the rest of the atmospheric parameters. In the next cycle, the code improves the previous atmosphere, changing its temperature at the bottom, the top and the middle of the atmosphere, and at the bottom and the top for the rest of the parameters (adding gradients along the LOS for the velocity and the magnetic field vector). From there, in each inversion cycle, the code starts from a perturbed atmosphere, which delivered reasonable fits, and changes it again based on the RF of the previous atmosphere. Each iteration gradually adds more freedom in incremental steps. In the second study, that is, starting with C4 directly and running one inversion cycle, the code probably has too much freedom when starting from an atmosphere (modified HSRA) that is far from the target solution; this means the inversion iterations that the code can make will therefore allow a solution to be achieved that, in general, is poorer than in the case of gradually increasing the freedom of the atmospheric parameters. Therefore, we encourage the users of node-based inversion codes to gradually increase the number of nodes in a manner similar to that presented in Table 1.
2.5. Inversion configuration IV: Weights
We evaluate whether or not we can improve the previous results using a different configuration in this section. SIR is an inversion code that synthesises the Stokes profiles from different atmospheric models until the synthetic spectrum matches the observed one (see, for instance, Ruiz Cobo et al. 2022). The process implies solving the radiative transfer equation multiple times while perturbing the atmospheric parameters until an accurate solution is achieved, that is, one where the difference between synthetic and observed profiles is as small as possible for the entire analysed spectral range. That accuracy is defined by the following χ2 merit function:
where index S runs over the four Stokes parameters, i covers the wavelength samples up to nλ, and Nf stands for the number of degrees of freedom, that is, the difference between the number of observables (four Stokes profiles ×nλ) and the number of free parameters used in the inversion (the sum of all the values included in each column in Table 1). The weights ws(λi) are ‘best-fit normalisation’ factors used to add more emphasis to a specific Stokes profile, for example, Stokes V, or a specific spectral line, when fitting multiple spectral lines simultaneously, promoting better fits of those individual features. In this section, knowing that the target is to fit weak polarisation signals in a quiet Sun scenario, we study the correlation between the input and the inferred atmosphere when we add more weight to the Stokes Q and U profiles.
We compare the four different weight configurations presented in Table 2. In each case, we run the inversion using configuration C4, increasing the weight given to the Stokes Q and U parameters. The results are presented in Fig. 9. We can see that in the case of temperature, the differences are similar except for W4 (where the weights are much higher for the linear polarisation parameters), which shows greater differences in general. This higher difference means that the weight is so significant for Stokes Q and U that the code does not fit Stokes I with reasonable accuracy. Hence, the greater deviations in the temperature. We do not recommend reaching this point, even if we are primarily interested in the magnetic field vector, because, in the end, the Stokes Q, U, and V profiles are derivatives of the Stokes I parameter. In other words, a poor fit for Stokes I will also mean poor results for the rest of the Stokes parameters. We can see this trend for the rest of the atmospheric parameters, where the grey line produces a lower correlation for the magnetic field strength and LOS velocity. It is comparable to the other solutions for the inclination, which means that the code accurately fits Stokes Q and U. Still, we are losing accuracy in the rest of the atmospheric parameters. It is therefore essential to bear in mind that weights should be chosen carefully, in this case, emphasising Stokes Q and U because the magnetic field is weak and inclined (see Fig. 2) but without noticeable degradation of the fit of Stokes I because the error will extend to the inferred temperature, velocity, and field strength. Therefore, the most balanced solutions are obtained using W2 and W3. Of the two, W2 (red) is slightly better for the magnetic field strength, but W2 and W3 are similar for the inclination and the temperature. We therefore plan to use W2 combined with C4, although we recommend users experiment with weights until the inversions provide good fits on average for the pixels of interest.
Fig. 9. Average differences between the input and the inferred atmosphere at different optical depths. From left to right, the differences for temperature, LOS velocity, magnetic field strength, and inclination, respectively. Colours correspond to the inversion configurations indicated in Table 2. |
Weight options used for the inversion of each atmospheric parameter.
2.6. Inversion configuration V: Microturbulence
The microturbulence is a correction added when computing the Doppler broadening; in SIR, it is defined in the Doppler width, ΔλD in the following way:
where λ is the spectral line rest wavelength, ξ is the microturbulence, T the temperature, M the mass of the atom involved in the transition, and c and k are the speed of light and the Boltzmann constant, respectively.
There is no complete consensus as to why the microturbulent velocity is used when fitting observations (see, for instance, Sect. 3 in Milić et al. 2019). However, we still add it as a free parameter in inversions to account for the complicated vertical stratification in the LOS velocity, which cannot be fully resolved by node-based inversion codes, mainly in the present scenario where the numerical simulations were computed with 7 and 8 km vertical resolution, respectively. Therefore, in this section, we evaluate the accuracy of the fits when using or not microturbulence as a free parameter. We reiterate here that we did not incorporate any microturbulence correction during the synthesis process (the original simulation runs do not include one). However, similarly to what was shown in Milić et al. (2019) and Ruiz Cobo et al. (2022), we expect the results to be more accurate when using it as a free parameter.
Figure 10 shows the average differences between the input and the inferred atmosphere when inverting the microturbulence with a constant value with height (solid) or not (dashed). Results for temperature are more accurate in the whole domain, indicating that the fits of Stokes I are also more accurate. This improved accuracy also translates to better fits in general for the rest of the Stokes profiles and smaller differences for all the atmospheric parameters. Therefore, as previous works have found, the inversion code provides better results when using microturbulence as a free parameter. We believe the main reasons for this behaviour are either that the degrees of freedom in the temperature stratification are insufficient in number (which is unlikely for this case) or that there is a strong coupling or degeneracy between temperature and microturbulence. In any case, and given our commitment to exploring these results in future works, we continue using configuration C4 with one node for the microturbulence in this publication.
Fig. 10. Average differences between the input and the inferred atmosphere at different optical depths. From left to right, the differences for temperature, LOS velocity, magnetic field strength, and inclination, respectively. Line styles correspond to the configuration C4 in Table 1 and W2 in Fig. 9 when setting the microturbulence as a free parameter (solid) or not (dashed). |
2.7. Inversion configuration VI: Spatial distribution of the atmospheric parameters
Above, we focused on the average differences between the input and the inferred atmospheres to define the inversion configuration. In this section, we examine the accuracy of the results in more detail, studying the spatial distribution of the atmospheric parameters at different optical depth layers. We use three atmospheric layers at log τ = [0, −1, −2] as a reference, because they cover the typical range of maximum sensitivity for the spectral lines at 630 nm (e.g., Quintero Noda et al. 2021). Finally, we compare the input and the inferred atmospheres using the node configuration C4 and the weights W2 in Fig. 11.
Fig. 11. Comparison between the input atmosphere (odd rows) and the one inferred with SIR using configuration C4 and the weights W2 (even rows). From left to right, we show the temperature, LOS velocity, magnetic field strength, and inclination. From top to bottom, we display the spatial distribution of the atmospheric parameters at three reference layers at log τ500 = [0, −1, −2]. The FOV corresponds to the spatial domain within the square region (orange) in the top panel of Fig. 1. Coloured squares in the top leftmost panel designate the locations of the pixels studied in Figs. 5 and 6. Black, blue, and red correspond to the top, middle, and bottom rows in the cited figures, respectively. |
Starting with the temperature, we can see that the results at log τ = 0 are very accurate, with almost no apparent differences between the input and the inferred atmosphere. We see a similar behaviour for higher layers, that is, log τ = −1, although some pixels deviate from the input atmosphere. This deviation continues increasing at log τ = −2, where some areas appear with much lower temperatures in the inferred atmosphere (see around [0.4,1.2] arcsec).
In the case of the LOS velocity, we see that the inversion results closely resemble that of the input atmosphere in all layers, even at log τ = −2, where only minor differences can be found.
Magnetic field strength results at log τ = 0 are reasonably accurate with many of the small-scale features recovered in the inversion, although there are more significant differences at around [0.6, 0.9] arcsec. This behaviour is also found at log τ = −1, with the inversion reproducing the complex and thin structures found in the selected FOV but also providing stronger magnetic fields in areas where the input parameter was weak (see e.g., [0.4, 1.2] arcsec). Higher up, at log τ = −2, the results resemble the input atmosphere, but the general accuracy has dropped substantially. Finally, the magnetic field inclination is well recovered at log τ = 0, with the inferred atmosphere following the intricate patterns where the magnetic field changes rapidly, even from one polarity to the opposite one, on spatial scales of a few pixels. The same happens for higher layers, with a general landscape that resembles the complexity of the input atmosphere, although that substantially drops in accuracy at log τ = −2.
2.8. Inversion configuration VII: Accuracy of the inferred atmosphere
In order to investigate the accuracy of the fits in more detail, in Fig. 12 we show a histogram containing the differences for all the pixels in the selected FOV between the input and the inferred atmosphere at the same optical depth layers examined in the previous figure. Again, we are interested in the accuracy for inferring all the atmospheric parameters, and so we show, from left to right, the results for temperature, LOS velocity, magnetic field strength, and inclination.
Fig. 12. Histogram containing the differences between the inferred and the input atmosphere at selected optical depths for the region enclosed by the orange square in Fig. 1. The differences in temperature, LOS velocity, magnetic field strength, and inclination are shown from left to right. Each row corresponds, from top to bottom, to reference layers at log τ500 = [0, −1, −2]. The bin size is 10 K, 0.1 km s−1, 10 G, and 2 degrees, respectively. |
Starting with temperature, we obtain results that peak close to zero but show that the inversion process tends to obtain slightly hotter atmospheres than the input ones. However, this trend is less pronounced at log τ = −1, where the peak is much closer to zero. Furthermore, this trend is reverted at log τ = −2 where the code provides slightly cooler atmospheres. In any case, the width of the histogram is relatively narrow, having a Gaussian-like shape with wings that drop substantially before the ±100 K mark. As mentioned in the previous section, LOS velocity results are very close to the input atmosphere at most layers. This can be seen with a narrow Gaussian-like function for the differences that peaks close to zero in the three selected atmospheric layers. In the case of the field strength, differences also peak at zero, although we see an extended right tail at log τ = 0, indicating that multiple pixels show larger amplitude for the field strength than the original ones. However, this tail is less prominent in the upper layers, indicating that the inferred magnetic field strength is close to the original one in general. Finally, for the inclination, we again see that the differences display a Gaussian-like profile that peaks at zero, although it has extended wings, mainly at log τ = −2. In those cases, the differences could be larger than 30 degrees, which can be considered as low accuracy for the magnetic field inclination. Still, in general, mainly at lower layers, we have a prominent peak at zero, pointing out that the fits are generally accurate.
3. Results for noisy spectra
The strategy for the following sections is similar for each selected FOV and simulation run. First, we add a noise degradation, which is different for every wavelength and Stokes parameter. Next, we modify the amplitude of the noise using the following values: [0, 1 × 10−4, 5 × 10−4, 1 × 10−3, 5 × 10−3, 1 × 10−2] of Ic. We then invert the Stokes profiles using the configuration C4 and the weights W2 (see Tables 1 and 2) for each noise amplitude and compute the average differences over the whole spatial domain at each optical depth. We start with the highlighted area (see the orange box in the leftmost column of Fig. 2) and the complete FOV of the Mancha run, and then move to the selected area (see the blue box in the rightmost column of Fig. 2), and the complete FOV of the MURaM run. The results for the inversion of the weakly magnetised areas in the Mancha run, that is, the same region used in Sect. 2, are displayed in the first row of Fig. 13. Each coloured line represents a different noise level, from no noise in black to 1 × 10−2 of Ic in yellow.
Fig. 13. Average differences between the input and the inferred atmosphere at different optical depths when adding noise to the synthetic profiles. From left to right, we show the differences for temperature, LOS velocity, magnetic field strength, and inclination, respectively. Colours correspond to the noise values indicated in the leftmost column. Rows are associated with different reference FOVs, with the top, second, third, and the bottom row corresponding to the region inside the orange box, the complete Mancha FOV, the region inside the blue box, and the complete FOV of the MURaM run presented in Fig. 1, respectively. |
Starting with the temperature, we can see that the differences are similar to low noise values and begin to increase as we increase the noise, being noticeably higher at 5 × 10−3 of Ic and with larger noise values. Interestingly, the same pattern is found for the rest of the atmospheric parameters with similar and small differences for noise levels up to 1 × 10−3 of Ic. Only the inclination shows more variations for different noise values, where only a noise level of 1 × 10−4 of Ic (red) is very close to the case without noise (black). The same pattern is found when inverting the Stokes profiles over the entire FOV presented in the top row of Fig. 1, with similar differences (second row in Fig. 13) for most of the atmospheric parameters when the noise levels are below 1 × 10−3 of Ic. Only the inclination, once again, is more sensitive to the noise values and shows noticeable deviations when the noise is larger than ×10−4 of Ic.
The results for the inversion of the strongly magnetised area in the MURaM run are displayed in the third row of Fig. 13. In this case, the results are similar to the previous case, that is, differences are mainly noticeable when the noise is larger 1 × 10−3 of Ic. However, in this scenario, the same pattern can be seen for the inclination with more significant deviations only when reaching 5 × 10−3 of Ic and 1 × 10−2 of Ic. This applies to the case where we only choose a strongly magnetised area (see third row) like the one highlighted in blue in the rightmost column of Fig. 2, and to the case where we use the entire FOV. This points out that when the average field strength is relatively large on average (see the bottom row in Fig. 3), the accuracy of the results – even for the inclination – is still high when the noise levels are on the order of 1 × 10−3 of Ic. Finally, it is worth mentioning that the black (noiseless) and red (1 × 10−4 of Ic) lines are almost identical in all the examined scenarios, indicating that the maximum accuracy the code can provide for the present configuration (e.g., node selection, spectral lines, etc.) is reached with noise values of that order.
4. Summary
In this work, we study the accuracy with which we can infer atmospheric information from realistic numerical MHD simulations that reproduce the spatial resolution we can achieve with future observations using the 4m class telescopes DKIST and EST. We used two intrinsically different simulation runs from different numerical codes (MANCHA and MURaM), simulating different scenarios: a weakly magnetised area and a sunspot. The aim is to increase the applicability of the analysis. We studied multiple inversion configurations until we obtained a good correlation between the input and the inferred atmosphere in a wide range of heights. We hope the thorough explanations provided in this step of the process will help future users of inversion codes define an inversion configuration that matches their needs and gets the most from their observations. Also, to continue providing a reference for future observing proposals, we examine how the inversion accuracy depends on the noise level of the Stokes profiles. The results indicate that when the majority of the inverted pixels come from strongly magnetised areas, we can obtain accurate results even when the amplitude is up to 1 × 10−3 of Ic. However, the situation is more critical for observations where the dominant magnetic structures are weak and inclined with respect to the solar surface, finding that inaccurate results will start to appear for the inferred magnetic field inclination when signals are on the order of or greater than 5 × 10−4 of Ic. At the same time, we find that the accuracy of the fit is almost identical to the noiseless case when the noise levels are on the order of 1 × 10−4 of Ic. Hence, we only recommend values lower than 5 × 10−4 of Ic for high-spatial-resolution quiet Sun observations – such as the one represented here – if the users of inversion codes aim to reliably interpret the results for the magnetic field vector. We understand that this noise floor seems almost impossible with current observatories unless the spatial information is lost due to the necessary spatial binning to achieve those noise values within reasonable integration times. However, we know, at the same time, that next-generation 4 m class telescopes are designed to achieve those levels (although perhaps not at the diffraction limit) with high spatial resolution thanks to the large area of the primary mirror, which can work as a photon collecting bucket at the cost of reducing the spatial resolution.
We also want to emphasise that although there is almost no difference between the results of the inversion without noise and a noise level on the order of 1 × 10−4 of Ic, on many occasions, both inversions results deviate from the original solution. We argue that this is probably a limitation of node-based inversion codes such as SIR and the spectral line selection, which fail to reproduce the stark small-scale variations in height present in the simulated atmospheres. Therefore, in the case of the former, we may, in the future, explore different inversion alternatives (for instance, Asensio Ramos & Díaz Baso 2019), while in the case of the latter, we plan to add complementary spectral lines with different sensitivity to the atmospheric parameters, which could help to improve the accuracy of the results. Also, the community has been pushing hard to deliver evermore complex and rich numerical simulations of different solar phenomena, such as those presented in Danilovic (2017), and Cheung et al. (2019), and so we plan to use them as laboratory scenarios as well.
Finally, the authors look forward to sharing any data products used in this work, either the atmospheres transformed to optical depth scale or the output spectra at 630 nm, with anyone interested in using them. Please contact the corresponding author with any requests.
Acknowledgments
This work has also been supported by the Spanish Ministry of Economy and Competitiveness through the project ESP-2016-77548-C5-1-R, RTI2018-096886-B-C53 and PID2021-125325OB-C51. This work was supported by Fundação para a Ciência e a Tecnologia (FCT) through the research grants UIDB/04434/2020 and UIDP/04434/2020. This material is based upon work supported by the National Center for Atmospheric Research, which is a major facility sponsored by the National Science Foundation under Cooperative Agreement No. 1852977. 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). D. Orozco Suárez also acknowledges financial support through the Ramón y Cajal fellowships.
References
- Asensio Ramos, A., & Díaz Baso, C. J. 2019, A&A, 626, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Balasubramaniam, K. S., Keil, S. L., & Tomczyk, S. 1997, ApJ, 482, 1065 [NASA ADS] [CrossRef] [Google Scholar]
- Beck, C., Fabbian, D., Rezaei, R., & Puschmann, K. G. 2017, ApJ, 842, 37 [NASA ADS] [CrossRef] [Google Scholar]
- Bellot Rubio, L., & Orozco Suárez, D. 2019, Liv. Rev. Sol. Phys., 16, 1 [Google Scholar]
- Borrero, J. M., & Kobel, P. 2011, A&A, 527, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Borrero, J. M., Martínez Pillet, V., Schmidt, W., et al. 2013, ApJ, 768, 69 [Google Scholar]
- Campbell, R. J., Shelyag, S., Quintero Noda, C., et al. 2021, A&A, 654, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Centeno, R. 2018, ApJ, 866, 89 [Google Scholar]
- Centeno, R., Socas-Navarro, H., Lites, B., et al. 2007, ApJ, 666, L137 [CrossRef] [Google Scholar]
- Cheung, M. C. M., Rempel, M., Chintzoglou, G., et al. 2019, Nat. Astron., 3, 160 [NASA ADS] [CrossRef] [Google Scholar]
- Collados, M., Lagg, A., Díaz Garcí A, J. J., et al. 2007, in The Physics of Chromospheric Plasmas, eds. P. Heinzel, I. Dorotovič, & R. J. Rutten, ASP Conf. Ser., 368, 611 [Google Scholar]
- Danilovic, S. 2017, A&A, 601, A122 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Danilovic, S., Beeck, B., Pietarila, A., et al. 2010, ApJ, 723, L149 [Google Scholar]
- Danilovic, S., van Noort, M., & Rempel, M. 2016, A&A, 593, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- de la Cruz Rodríguez, J., Socas-Navarro, H., Carlsson, M., & Leenaarts, J. 2012, A&A, 543, A34 [Google Scholar]
- de la Cruz Rodríguez, J., & van Noort, M. 2017, Space. Sci. Rev., 210, 109 [CrossRef] [Google Scholar]
- de Wijn, A. G., Casini, R., Carlile, A., et al. 2022, Sol. Phys., 297, 22 [NASA ADS] [CrossRef] [Google Scholar]
- del Toro Iniesta, J. C., & Ruiz Cobo, B. 2016, Liv. Rev. Sol. Phys., 13, 4 [Google Scholar]
- Domínguez Cerdeña, I., Sánchez Almeida, J., & Kneer, F. 2006, ApJ, 646, 1421 [CrossRef] [Google Scholar]
- Elmore, D. F., Lites, B. W., Tomczyk, S., et al. 1992, in Polarization Analysis and Measurement, eds. D. H. Goldstein, & R. A. Chipman, SPIE Conf. Ser., 1746, 22 [NASA ADS] [CrossRef] [Google Scholar]
- Felipe, T., Khomenko, E., & Collados, M. 2010, ApJ, 719, 357 [Google Scholar]
- Gafeira, R., Orozco Suárez, D., Milić, I., et al. 2021, A&A, 651, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gingerich, O., Noyes, R. W., Kalkofen, W., & Cuny, Y. 1971, Sol. Phys., 18, 347 [Google Scholar]
- González-Morales, P. A., Khomenko, E., Downes, T. P., & de Vicente, A. 2018, A&A, 615, A67 [Google Scholar]
- Khomenko, E., & Collados, M. 2006, ApJ, 653, 739 [Google Scholar]
- Khomenko, E., Vitas, N., Collados, M., & de Vicente, A. 2017, A&A, 604, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Khomenko, E. V., Collados, M., Solanki, S. K., Lagg, A., & Trujillo Bueno, J. 2003, A&A, 408, 1115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kosugi, T., Matsuzaki, K., Sakao, T., et al. 2007, Sol. Phys., 243, 3 [Google Scholar]
- Landi Degl’Innocenti, E., & Landolfi, M. 2004, Polarization in Spectral Lines, 307 (Dordrecht: Kluwer Academic Publishers) [Google Scholar]
- Lites, B. W. 2002, ApJ, 573, 431 [NASA ADS] [CrossRef] [Google Scholar]
- Lites, B. W., & Socas-Navarro, H. 2004, ApJ, 613, 600 [NASA ADS] [CrossRef] [Google Scholar]
- Lites, B. W., Kubo, M., Socas-Navarro, H., et al. 2008, ApJ, 672, 1237 [Google Scholar]
- Lites, B. W., Akin, D. L., Card, G., et al. 2013, Sol. Phys., 283, 579 [NASA ADS] [CrossRef] [Google Scholar]
- Martínez González, M. J., & Bellot Rubio, L. R. 2009, ApJ, 700, 1391 [CrossRef] [Google Scholar]
- Martínez Pillet, V., Lites, B. W., & Skumanich, A. 1997, ApJ, 474, 810 [CrossRef] [Google Scholar]
- Martínez Pillet, V., del Toro Iniesta, J. C., Álvarez-Herrero, A., et al. 2011, Sol. Phys., 268, 57 [CrossRef] [Google Scholar]
- Metcalf, T. R., Leka, K. D., Barnes, G., et al. 2006, Sol. Phys., 237, 267 [Google Scholar]
- Milić, I., Smitha, H. N., & Lagg, A. 2019, A&A, 630, A133 [Google Scholar]
- Quintero Noda, C., Borrero, J. M., Orozco Suárez, D., & Ruiz Cobo, B. 2014, A&A, 569, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Quintero Noda, C., Barklem, P. S., Gafeira, R., et al. 2021, A&A, 652, A161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Quintero Noda, C., Schlichenmaier, R., Bellot Rubio, L. R., et al. 2022, A&A, 666, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rempel, M. 2012, ApJ, 750, 62 [Google Scholar]
- Rimmele, T. R., Warner, M., Keil, S. L., et al. 2020, Sol. Phys., 295, 172 [Google Scholar]
- Ruiz Cobo, B., & del Toro Iniesta, J. C. 1992, ApJ, 398, 375 [Google Scholar]
- Ruiz Cobo, B., Quintero Noda, C., Gafeira, R., et al. 2022, A&A, 660, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sainz Dalda, A., Martínez-Sykora, J., Bellot Rubio, L., & Title, A. 2012, ApJ, 748, 38 [NASA ADS] [CrossRef] [Google Scholar]
- Sánchez Almeida, J. 2004, in The Solar-B Mission and the Forefront of Solar Physics, eds. T. Sakurai, & T. Sekii, ASP Conf. Ser., 325, 115 [Google Scholar]
- Shimizu, T., Nagata, S., Tsuneta, S., et al. 2008, Sol. Phys., 249, 221 [Google Scholar]
- Smitha, H. N., Holzreuter, R., van Noort, M., & Solanki, S. K. 2020, A&A, 633, A157 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Solanki, S. K., Barthol, P., Danilovic, S., et al. 2010, ApJ, 723, L127 [NASA ADS] [CrossRef] [Google Scholar]
- Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, Sol. Phys., 249, 167 [Google Scholar]
- Viticchié, B., & Sánchez Almeida, J. 2011, A&A, 530, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vögler, A., Shelyag, S., Schüssler, M., et al. 2005, A&A, 429, 335 [Google Scholar]
- von der Lühe, O. 1998, New Astron. Rev., 42, 493 [Google Scholar]
All Tables
All Figures
Fig. 1. Spatial distribution of the temperature where the continuum optical depth at 500 nm is unity. The top panel shows the results from a quiet-Sun run with 10 km spatial resolution computed with the MANCHA code. The bottom panel corresponds to a cut of the original Sunspot run with a 12 km spatial resolution generated with the MURaM code. The regions highlighted with coloured boxes designate the areas of interest on which we focus in this work. |
|
In the text |
Fig. 2. Spatial distribution of the horizontal (top) and longitudinal (bottom) component of the magnetic field where log τ = 0. The left panel corresponds to a quiet Sun run with 10 km spatial resolution generated with the MANCHA code, while the right panel shows a cut of the original Sunspot run with 12 km spatial resolution computed with the MURaM code. |
|
In the text |
Fig. 3. Evolution of the averaged atmospheric parameters with height. Top and bottom panels correspond to the MANCHA and MURaM runs presented in Fig. 1. Colours designate the averaged results over the small regions highlighted in the same figure. Bars display the standard deviation of the atmospheric parameters over the corresponding areas. |
|
In the text |
Fig. 4. Average differences between the input and the inferred atmosphere at different optical depths. From left to right, we show the differences in temperature, LOS velocity, magnetic field strength, and inclination, respectively. Colours correspond to the various inversion configurations indicated in Table 1. |
|
In the text |
Fig. 5. Stokes profiles corresponding to the three reference pixels (rows). From left to right, Stokes (I, Q, U, V), respectively. Black corresponds to the input profiles, while colours designate the Stokes parameters corresponding to the different node configurations presented in Table 1. Each panel includes the differences between the input and the inferred profiles below the horizontal solid black line. From top to bottom, each reference profile corresponds to an intergranular lane, edge of a granule, and top of a granule, respectively. We show the exact location in Fig. 11, where we provide a zoom onto the region highlighted in the top panel of Fig. 1. |
|
In the text |
Fig. 6. Atmospheres corresponding to the pixels presented in Fig. 5 (different rows). We show the temperature, LOS velocity, magnetic field strength, and inclination from left to right. Black corresponds to the input profiles, while colours designate the Stokes parameters corresponding to the different node configurations presented in Table 1. |
|
In the text |
Fig. 7. Reference location of nodes in SIR. The cases corresponding to 2, 3, and 5 nodes (see Table 1) have nodes located as shown by triangles, squares and diamonds, respectively. As a reference, we include the temperature of the HSRA atmosphere. |
|
In the text |
Fig. 8. Average differences between the input and the inferred atmosphere at different optical depths. From left to right, we show the differences for temperature, LOS velocity, magnetic field strength and inclination, respectively. Line styles correspond to the case where we reach configuration C4 in Table 1 after running the previous cycles (solid) or to cases where we run a single inversion cycle using the nodes described in column C4 of Table 1 starting from the modified HSRA atmosphere (dashed). |
|
In the text |
Fig. 9. Average differences between the input and the inferred atmosphere at different optical depths. From left to right, the differences for temperature, LOS velocity, magnetic field strength, and inclination, respectively. Colours correspond to the inversion configurations indicated in Table 2. |
|
In the text |
Fig. 10. Average differences between the input and the inferred atmosphere at different optical depths. From left to right, the differences for temperature, LOS velocity, magnetic field strength, and inclination, respectively. Line styles correspond to the configuration C4 in Table 1 and W2 in Fig. 9 when setting the microturbulence as a free parameter (solid) or not (dashed). |
|
In the text |
Fig. 11. Comparison between the input atmosphere (odd rows) and the one inferred with SIR using configuration C4 and the weights W2 (even rows). From left to right, we show the temperature, LOS velocity, magnetic field strength, and inclination. From top to bottom, we display the spatial distribution of the atmospheric parameters at three reference layers at log τ500 = [0, −1, −2]. The FOV corresponds to the spatial domain within the square region (orange) in the top panel of Fig. 1. Coloured squares in the top leftmost panel designate the locations of the pixels studied in Figs. 5 and 6. Black, blue, and red correspond to the top, middle, and bottom rows in the cited figures, respectively. |
|
In the text |
Fig. 12. Histogram containing the differences between the inferred and the input atmosphere at selected optical depths for the region enclosed by the orange square in Fig. 1. The differences in temperature, LOS velocity, magnetic field strength, and inclination are shown from left to right. Each row corresponds, from top to bottom, to reference layers at log τ500 = [0, −1, −2]. The bin size is 10 K, 0.1 km s−1, 10 G, and 2 degrees, respectively. |
|
In the text |
Fig. 13. Average differences between the input and the inferred atmosphere at different optical depths when adding noise to the synthetic profiles. From left to right, we show the differences for temperature, LOS velocity, magnetic field strength, and inclination, respectively. Colours correspond to the noise values indicated in the leftmost column. Rows are associated with different reference FOVs, with the top, second, third, and the bottom row corresponding to the region inside the orange box, the complete Mancha FOV, the region inside the blue box, and the complete FOV of the MURaM run presented in Fig. 1, respectively. |
|
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.