Accurately constraining velocity information from spectral imaging observations using machine learning techniques.

Determining accurate plasma Doppler (line-of-sight) velocities from spectroscopic measurements is a challenging endeavour, especially when weak chromospheric absorption lines are often rapidly evolving and, hence, contain multiple spectral components in their constituent line profiles. Here, we present a novel method that employs machine learning techniques to identify the underlying components present within observed spectral lines, before subsequently constraining the constituent profiles through single or multiple Voigt fits. Our method allows active and quiescent components present in spectra to be identified and isolated for subsequent study. Lastly, we employ a Ca ɪɪ 8542 Å spectral imaging dataset as a proof-of-concept study to benchmark the suitability of our code for extracting two-component atmospheric profiles that are commonly present in sunspot chromospheres. Minimization tests are employed to validate the reliability of the results, achieving median reduced χ2-values equal to 1.03 between the observed and synthesized umbral line profiles. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

Characterization of shock wave signatures at millimetre wavelengths from Bifrost simulations.

Observations at millimetre wavelengths provide a valuable tool to study the small-scale dynamics in the solar chromosphere. We evaluate the physical conditions of the atmosphere in the presence of a propagating shock wave and link that to the observable signatures in mm-wavelength radiation, providing valuable insights into the underlying physics of mm-wavelength observations. A realistic numerical simulation from the three-dimensional radiative magnetohydrodynamic code Bifrost is used to interpret changes in the atmosphere caused by shock wave propagation. High-cadence (1 s) time series of brightness temperature (Tb) maps are calculated with the Advanced Radiative Transfer code at the wavelengths 1.309 mm and 1.204 mm, which represents opposite sides of spectral band 6 of the Atacama Large Millimeter/submillimeter Array (ALMA). An example of shock wave propagation is presented. The brightness temperatures show a strong shock wave signature with large variation in formation height between approximately 0.7 and 1.4 Mm. The results demonstrate that millimetre brightness temperatures efficiently track upwardly propagating shock waves in the middle chromosphere. In addition, we show that the gradient of the brightness temperature between wavelengths within ALMA band 6 can potentially be used as a diagnostics tool in understanding the small-scale dynamics at the sampled layers. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

Magnetoacoustic wave energy dissipation in the atmosphere of solar pores.

The suitability of solar pores as magnetic wave guides has been a key topic of discussion in recent years. Here, we present observational evidence of propagating magnetohydrodynamic wave activity in a group of five photospheric solar pores. Employing data obtained by the Facility Infrared Spectropolarimeter at the Dunn Solar Telescope, oscillations with periods of the order of 5 min were detected at varying atmospheric heights by examining Si ɪ 10827 Å line bisector velocities. Spectropolarimetric inversions, coupled with the spatially resolved root mean square bisector velocities, allowed the wave energy fluxes to be estimated as a function of atmospheric height for each pore. We find propagating magnetoacoustic sausage mode waves with energy fluxes on the order of 30 kW m-2 at an atmospheric height of 100 km, dropping to approximately 2 kW m-2 at an atmospheric height of around 500 km. The cross-sectional structuring of the energy fluxes reveals the presence of both body- and surface-mode sausage waves. Examination of the energy flux decay with atmospheric height provides an estimate of the damping length, found to have an average value across all five pores of Ld ≈ 268 km, similar to the photospheric density scale height. We find the damping lengths are longer for body mode waves, suggesting that surface mode sausage oscillations are able to more readily dissipate their embedded wave energies. This work verifies the suitability of solar pores to act as efficient conduits when guiding magnetoacoustic wave energy upwards into the outer solar atmosphere. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

High-resolution wave dynamics in the lower solar atmosphere.

The magnetic and convective nature of the Sun's photosphere provides a unique platform from which generated waves can be modelled, observed and interpreted across a wide breadth of spatial and temporal scales. As oscillations are generated in-situ or emerge through the photospheric layers, the interplay between the rapidly evolving densities, temperatures and magnetic field strengths provides dynamic evolution of the embedded wave modes as they propagate into the tenuous solar chromosphere. A focused science team was assembled to discuss the current challenges faced in wave studies in the lower solar atmosphere, including those related to spectropolarimetry and radiative transfer in the optically thick regions. Following the Theo Murphy international scientific meeting held at Chicheley Hall during February 2020, the scientific team worked collaboratively to produce 15 independent publications for the current Special Issue, which are introduced here. Implications from the current research efforts are discussed in terms of upcoming next-generation observing and high-performance computing facilities. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

The Role of Twist in Kinked Flux Rope Emergence and Delta-spot Formation.

It has been observationally well established that the magnetic configurations most favorable for producing energetic flaring events reside in δ-spots, a class of sunspots defined as having opposite-polarity umbrae sharing a common penumbra. They are frequently characterized by extreme compactness, strong rotation, and anti-Hale orientation. Numerous studies have shown that nearly all of the largest solar flares originate in δ-spots, making the understanding of these structures a fundamental step in predicting space weather. Despite their important influence on the space environment, surprisingly little is understood about the origin and behavior of δ-spots. In this paper, we perform a systematic study of the behavior of emerging flux ropes to test a theoretical model for the formation of δ-spots: the kink instability of emerging flux ropes. We simulated the emergence of highly twisted, kink-unstable flux ropes from the convection zone into the corona, and we compared their photospheric properties to those of emerged weakly twisted, kink-stable flux ropes. We show that the photospheric manifestations of the emergence of highly twisted flux ropes closely match the observed properties of δ-spots, and we discuss the resulting implications for observations. Our results strongly support and extend previous theoretical work that suggested that the kink instability of emerging flux ropes is a promising candidate to explain δ-spot formation, as it reproduces their key characteristics very well.

Structure of the solar photosphere studied from the radiation hydrodynamics code ANTARES.

The ANTARES radiation hydrodynamics code is capable of simulating the solar granulation in detail unequaled by direct observation. We introduce a state-of-the-art numerical tool to the solar physics community and demonstrate its applicability to model the solar granulation. The code is based on the weighted essentially non-oscillatory finite volume method and by its implementation of local mesh refinement is also capable of simulating turbulent fluids. While the ANTARES code already provides promising insights into small-scale dynamical processes occurring in the quiet-Sun photosphere, it will soon be capable of modeling the latter in the scope of radiation magnetohydrodynamics. In this first preliminary study we focus on the vertical photospheric stratification by examining a 3-D model photosphere with an evolution time much larger than the dynamical timescales of the solar granulation and of particular large horizontal extent corresponding to [Formula: see text] on the solar surface to smooth out horizontal spatial inhomogeneities separately for up- and downflows. The highly resolved Cartesian grid thereby covers [Formula: see text] of the upper convection zone and the adjacent photosphere. Correlation analysis, both local and two-point, provides a suitable means to probe the photospheric structure and thereby to identify several layers of characteristic dynamics: The thermal convection zone is found to reach some ten kilometers above the solar surface, while convectively overshooting gas penetrates even higher into the low photosphere. An [Formula: see text] wide transition layer separates the convective from the oscillatory layers in the higher photosphere.

SEARCH FOR A SIGNATURE OF TWIST-REMOVAL IN THE MAGNETIC FIELD OF SUNSPOTS IN RELATION WITH MAJOR FLARES.

We investigate the restructuring of the magnetic field in sunspots associated with two flares: the X6.5 flare on 6 December 2006 and the X2.2 flare on 15 February 2011. The observed changes were evaluated with respect to the so-called twist-removal model, in which helicity (twist) is removed from the corona as the result of an eruption. Since no vector magnetograms were available for the X6.5 flare, we applied the azimuthal symmetry approach to line-of-sight magnetograms to reconstruct the pseudo-vector magnetic field and investigate the changes in average twist and inclination of magnetic field in the sunspot around the time of the flare. For the X2.2 flare, results from the full vector magnetograms were compared with the pseudo-vector field data. For both flares, the data show changes consistent with the twist-removal scenario. We also evaluate the validity of the azimuthal symmetry approach on simple isolated round sunspots. In general, the derivations based on the azimuthal symmetry approach agree with true-vector field data though we find that even for symmetric sunspots the distribution of the magnetic field may deviate from an axially symmetric distribution.

HARD X-RAY EMISSION DURING FLARES AND PHOTOSPHERIC FIELD CHANGES.

We study the correlation between abrupt permanent changes of magnetic field during X-class flares observed by the Global Oscillation Network Group and Helioseismic and Magnetic Imager instruments, and the hard X-ray (HXR) emission observed by RHESSI, to relate the photospheric field changes to the coronal restructuring and investigate the origin of the field changes. We find that spatially the early RHESSI emission corresponds well to locations of the strong field changes. The field changes occur predominantly in the regions of strong magnetic field near the polarity inversion line (PIL). The later RHESSI emission does not correspond to significant field changes as the flare footpoints are moving away from the PIL. Most of the field changes start before or around the start time of the detectable HXR signal, and they end at about the same time or later than the detectable HXR flare emission. Some of the field changes propagate with speed close to that of the HXR footpoint at a later phase of the flare. The propagation of the field changes often takes place after the strongest peak in the HXR signal when the footpoints start moving away from the PIL, i.e., the field changes follow the same trajectory as the HXR footpoint, but at an earlier time. Thus, the field changes and HXR emission are spatio-temporally related but not co-spatial nor simultaneous. We also find that in the strongest X-class flares the amplitudes of the field changes peak a few minutes earlier than the peak of the HXR signal. We briefly discuss this observed time delay in terms of the formation of current sheets during eruptions.