Context. The study of magnetic activity in the Sun's polar regions is essential for understanding the solar cycle. However, measuring polar magnetic fields presents challenges due to projection effects and their intrinsically weak magnetic field strength. Faculae, bright regions on the visible solar surface associated with increased magnetic activity, offer a valuable proxy for measuring polar fields.

Aims. This research aims to analyze the magnetic activity of the Sun's polar regions through the use of polar faculae.

Methods. A neural network model (U-Net) was employed to detect polar faculae in images from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). The model was trained on synthetic data, eliminating the need for manual labeling, and was used to analyze 14 years of data from May 2010 to May 2024.

Results. The U-Net model demonstrates superior performance and efficiency over existing methods, enabling automated large-scale studies. We find that polar faculae numbers exhibit cyclical behavior with distinct minima and maxima, showing similar patterns between poles but with notable temporal delays (south pole: minimum early 2014, maximum late 2016; north pole: minimum late 2014, maximum mid-2019). Polar faculae magnetic fields remain consistent in magnitude (∼±75 G) across both poles and throughout the solar cycle. A strong linear correlation was found between the polar faculae count and the overall polar magnetic field strength. The spatio-temporal evolution reveals systematic migration of field polarity reversals from mid-latitudes toward the poles at rates of 3−8 m/s. During solar minimum, we observe a small relative increase in stronger-field faculae compared to solar maximum, suggesting either the coexistence of two magnetic distributions or subtle solar cycle dependence in faculae properties.

Context. Spatially resolved observations of the Sun and the astronomical sample size of stellar bodies are the respective key strengths of solar and stellar observations. However, the large difference in object brightness between the Sun and other stars has led to distinctly different instrumentation and methodologies between the two fields.Aims. We produced and analyzed synthetic full-disk spectra derived from 19 small area field-of-view optical observations of solar flares acquired by the Swedish 1-m Solar Telescope (SST) between 2011 and 2024. These were used to investigate what can and cannot be inferred about physical processes on the Sun from Sun-as-a-star observations.Methods. The recently released Numerical Empirical Sun-as-a-Star Integrator (NESSI) code provides synthetic full-disk integrated spectral line emission based on smaller field-of-view input while accounting for center-to-limb variations and differential rotation. We used this code to generate pseudo-Sun-as-a-star spectra from the SST observations.Results. We show that limited-area solar observations can be extrapolated to represent the full disk accurately in a manner close to what is achievable with Sun-as-a-star telescopes. Additionally, we identify nine spectral features, four of which are caused by instrumental effects. Most notably, we find a relation between the heliocentric angle of flares and the width of the excess emission left by them as well as a source of false positive coronal mass ejections-like signatures, and we defined an energy scaling law based on chromospheric line intensities that shows that the peak flare contrast roughly scales with the square root of the bolometric energy.Conclusions. The presented method of creating pseudo-Sun-as-a-star observations from limited field-of-view solar observations allows for the accurate comparison of solar flare spectra with their stellar counterparts while allowing for the detection of signals at otherwise unachievable noise levels.

Context. Sunspots survive on the solar surface for timescales ranging from days to months. This requires them to be in an equilibrium involving magnetic fields and hydrodynamic forces. Unfortunately, theoretical models of sunspot equilibrium are very simplified as they assume that spots are static and possess a self-similar and axially symmetric magnetic field. These assumptions neglect the role of small-scale variations of the magnetic field along the azimuthal direction produced by umbral dots, light bridges, penumbral filaments, and so forth.

Aims. We aim to study whether sunspot equilibrium is maintained once azimuthal fluctuations in the magnetic field, produced by the sunspot fine structure, are taken into account.

Methods. We apply the FIRTEZ Stokes inversion code to spectropolarimetric observations to infer the magnetic and thermodynamic parameters in two sunspots located at the disk center and observed with two different instruments: one observed from the ground with the 1.5-meter German GREGOR Telescope and another with the Japanese spacecraft Hinode. We compare our results with three-dimensional radiative magnetohydrodynamic simulations of a sunspot carried out with the MuRAM code.

Results. We infer clear variations in the gas pressure and density of the plasma directly related to fluctuations in the Lorentz force and associated with the filamentary structure in the penumbra. Similar results are obtained in the umbra despite its lack of an observed filamentary structure. Results from the two observed sunspots are in excellent qualitative and quantitative agreement with the numerical simulations.

Conclusions. Our results indicate that the magnetic topology of sunspots along the azimuthal direction is very close to magnetohydrostatic equilibrium, thereby helping to explain why sunspots are such long-lived structures capable of surviving on the solar surface for days or even full solar rotations.

Context. Inferences of the magnetic field in the solar atmosphere by means of spectropolarimetric inversions (i.e., Stokes inversion codes) yield magnetic fields that are non-solenoidal (∇ B ≠ 0). Because of this, results obtained by such methods are sometimes put into question. Aims. We aim to develop and implement a new technique that, in conjunction with Stokes inversion codes, can retrieve magnetic fields that are simultaneously consistent with observed polarization signals and with the null divergence condition. Methods. The method used in this work strictly imposes ∇ B=0 by determining the vertical component of the magnetic field (Bz) from the horizontal ones (BxBy). We implement this technique, which we refer to as solenoidal inversion, into the FIRTEZ Stokes inversion code and apply it to spectropolarimetric observations of a sunspot observed with the Hinode/SP instrument. Results. We show that the solenoidal inversion retrieves a vertical component of the magnetic field that is consistent in 80% of the analyzed three-dimensional (x, y, z) domain, with the vertical component of the magnetic field inferred from the non-solenoidal inversion. We demonstrate that the solenoidal inversion is capable of a better overall fitting to the observed Stokes vector than the non-solenoidal inversion. In fact, the solenoidal magnetic field fits Stokes V worse, but this is compensated by a better fit to Stokes I. We find a direct correlation between the worsening in the fit to the circular polarization profiles by the solenoidal inversion and the deviations in the inferred Bz with respect to the non-solenoidal inversion. Finally, we also show that the spatial distribution of the electric currents given by ∇x B does not change significantly after imposing the null divergence condition. Conclusions. In spite of being physically preferable, solenoidal magnetic fields are topologically very similar in 80% of the analyzed three-dimensional domain to the non-solenoidal fields obtained from spectropolarimetric inversions. These results support the idea that common Stokes inversion techniques fail to reproduce ∇ B=0 mainly as a consequence of the uncertainties in the determination of the individual components of the magnetic field. In the remaining 20% of the analyzed domain, where the Bz inferred by the solenoidal and non-solenoidal inversions disagree, it remains to be proven that the solenoidal inversion is to be preferred because even though the overall fit to the Stokes parameters improves, the fit to Stokes V worsens. It is in these regions where the application of the Stokes inversion constrained by the null divergence condition can yield new insights about the topology of the magnetic field in the solar photosphere.

Context. One-dimensional, semi-empirical models of the solar atmosphere are widely employed in numerous contexts within solar physics, ranging from the determination of element abundances and atomic parameters to studies of the solar irradiance and from Stokes inversions to coronal extrapolations. These models provide the physical parameters (i.e. temperature, gas pressure, etc.) in the solar atmosphere as a function of the continuum optical depth τc. The transformation to the geometrical z scale (i.e. vertical coordinate) is provided via vertical hydrostatic equilibrium. Aims. Our aim is to provide updated, one-dimensional, semi-empirical models of the solar atmosphere as a function of z, but employing the more general case of three-dimensional magneto-hydrostatic equilibrium (MHS) instead of vertical hydrostatic equilibrium (HE). Methods. We employed a recently developed Stokes inversion code that, along with non-local thermodynamic equilibrium effects, considers MHS instead of HE. This code is applied to spatially and temporally resolved spectropolarimetric observations of the quiet Sun obtained with the CRISP instrument attached to the Swedish Solar Telescope. Results. We provide average models for granules, intergranules, dark magnetic elements, and overall quiet-Sun as a function of both τc and z from the photosphere to the lower chromosphere. Conclusions. We demonstrate that, in these quiet-Sun models, the effect of considering MHS instead of HE is negligible. However, employing MHS increases the consistency of the inversion results before averaging. We surmise that in regions with stronger magnetic fields (i.e. pores, sunspots, network) the benefits of employing the magneto-hydrostatic approximation will be much more palpable.

The focus of this investigation is to quantify the conversion of magnetic to thermal energy initiated by a quiet Sun cancellation event and to explore the resulting dynamics from the interaction of the opposite-polarity magnetic features. We used imaging spectroscopy in the Hα line, along with spectropolarimetry in the Fe I 6173 Å and Ca II 8542 Å lines from the Swedish Solar Telescope (SST) to study a reconnection-related cancellation and the appearance of a quiet Sun Ellerman bomb (QSEB). We observed, for the first time, QSEB signature in both the wings and core of the Fe I 6173 Å line. We also found that, at times, the Fe I line-core intensity reaches higher values than the quiet Sun continuum intensity. From FIRTEZ-dz inversions of the Stokes profiles in Fe I and Ca II lines, we found enhanced temperature, with respect to the quiet Sun values, at the photospheric (log τ_c  = −1.5; ∼1000 K) and lower chromospheric heights (log τ_c  = −4.5; ∼360 K). From the calculation of total magnetic energy and thermal energy within these two layers, it was confirmed that the magnetic energy released during the flux cancellation can support heating in the aforesaid height range. Further, the temperature stratification maps enabled us to identify cumulative effects of successive reconnection on temperature pattern, including recurring temperature enhancements. Similarly, Doppler velocity stratification maps revealed impacts on plasma flow pattern, such as a sudden change in the flow direction. 

In this work, we study the accuracy that can be achieved when inferring the atmospheric information from realistic numerical magnetohydrodynamic simulations that reproduce the spatial resolution we will obtain with future observations made by the 4 m 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 x 10(-3) of I-c. 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 x 10(-4) of I-c. We, therefore, advise aiming for noise values on the order of or lower than 5 x 10(-4) of I-c 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.

Context. Electric currents play an important role in the energy balance of the plasma in the solar atmosphere. They are also indicative of non-potential magnetic fields and magnetic reconnection. Unfortunately, the direct measuring of electric currents has traditionally been riddled with inaccuracies.

Aims. We study how accurately we can infer electric currents under different scenarios.

Methods. We carry out increasingly complex inversions of the radiative transfer equation for polarized light applied to Stokes profiles synthesized from radiative three-dimensional magnetohydrodynamic (MHD) simulations. The inversion yields the magnetic field vector, B, from which the electric current density, j, is derived by applying Ampere’s law.

Results. We find that the retrieval of the electric current density is only slightly affected by photon noise or spectral resolution. However, the retrieval steadily improves as the Stokes inversion becomes increasingly elaborated. In the least complex case (a Milne-Eddington-like inversion applied to a single spectral region), it is possible to determine the individual components of the electric current density (j_x, j_y, j_z) with an accuracy of σ = 0.90 − 1.00 dex, whereas the modulus (∥j∥) can only be determined with σ = 0.75 dex. In the most complicated case (with multiple spectral regions, a large number of nodes, Tikhonov vertical regularization, and magnetohydrostatic equilibrium), these numbers improve to σ = 0.70 − 0.75 dex for the individual components and σ = 0.5 dex for the modulus. Moreover, in regions where the magnetic field is above 300 gauss, ∥j∥ can be inferred with an accuracy of σ = 0.3 dex. In general, the x and y components of the electric current density are retrieved slightly better than the z component. In addition, the modulus of the electric current density is the best retrieved parameter of all, and thus it can potentially be used to detect regions of enhanced Joule heating.

Conclusions. The fact that the accuracy does not worsen with decreasing spectral resolution or increasing photon noise, and instead increases as the Stokes inversion complexity grows, suggests that the main source of errors in the determination of electric currents is the lack of realism in the inversion model employed to determine variations in the magnetic field along the line of sight at scales smaller than the photon mean-free path, along with the intrinsic limitations of the model due to radiative transfer effects.

Ca ii K grains, i.e., intermittent, short-lived (about 1 minute), periodic (2-4 minutes), pointlike chromospheric brightenings, are considered to be the manifestations of acoustic waves propagating upward from the solar surface and developing into shocks in the chromosphere. After the simulations of Carlsson and Stein, we know that hot shocked gas moving upward interacting with the downflowing chromospheric gas (falling down after having been displaced upward by a previous shock) nicely reproduces the spectral features of the Ca ii K profiles observed in such grains, i.e., a narrowband emission-like feature at the blue side of the line core. However, these simulations are one-dimensional and cannot explain the location or the pointlike shape of the grains. Here, we report on the magnetic nature of these events. Furthermore, we report on similar events occurring at the largest flux concentrations, though they are longer-lived (up to 8 minutes) and exhibit the typical signature of steep velocity gradients traveling across the atmosphere. The spectral signatures of the studied events resemble their counterparts in sunspots, the umbral flashes. We then propose that magnetohydrodynamical waves are not only channeled through the magnetic field in sunspots, but they pervade the whole atmosphere. The propagation along magnetic fields can explain the pointlike appearance of the calcium grains observed in the quiet chromosphere.

We study the impact of the loss of axial symmetry around the optical axis on the polarimetric properties of a telescope with a segmented primary mirror when each segment is present in a different aging stage. The different oxidation stage of each segment as it is substituted in time leads to nonnegligible cross-talk terms. This effect is wavelength dependent, and it is mainly determined by the properties of the reflecting material. For an aluminum coating, the worst polarimetric behavior due to oxidation is found for the blue part of the visible. Contrarily, dust—as modeled in this work—does not significantly change the polarimetric behavior of the optical system. Depending on the telescope, there might be segment substitution sequences that strongly attenuate this instrumental polarization.