High-resolution Wave Dynamics Research Articles

Influence of ambipolar and Hall effects on vorticity in three-dimensional simulations of magneto-convection: Vorticity with ambipolar and Hall effect

This paper presents the results of the analysis of three-dimensional simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation runs are compared: One including both ambipolar diffusion and the Hall effect; one including only ambipolar diffusion and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 h of simulated solar time after the local dynamo reaches the saturation regime. We analyse the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect. 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'.

Proper orthogonal and dynamic mode decomposition of sunspot data.

High-resolution solar observations show the complex structure of the magnetohydrodynamic (MHD) wave motion. We apply the techniques of proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) to identify the dominant MHD wave modes in a sunspot using the intensity time series. The POD technique was used to find modes that are spatially orthogonal, whereas the DMD technique identifies temporal orthogonality. Here, we show that the combined POD and DMD approaches can successfully identify both sausage and kink modes in a sunspot umbra with an approximately circular cross-sectional shape. 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'.

On the influence of magnetic topology on the propagation of internal gravity waves in the solar atmosphere.

The solar surface is a continuous source of internal gravity waves (IGWs). IGWs are believed to supply the bulk of the wave energy for the lower solar atmosphere, but their existence and role for the energy balance of the upper layers is still unclear, largely due to the lack of knowledge about the influence of the Sun's magnetic fields on their propagation. In this work, we look at naturally excited IGWs in realistic models of the solar atmosphere and study the effect of different magnetic field topographies on their propagation. We carry out radiation-magnetohydrodynamic simulations of a magnetic field free and two magnetic models-one with an initial, homogeneous, vertical field of 100 G magnetic flux density and one with an initial horizontal field of 100 G flux density. The propagation properties of IGWs are studied by examining the phase-difference and coherence spectra in the kh - ω diagnostic diagram. We find that IGWs in the upper solar atmosphere show upward propagation in the model with predominantly horizontal field similar to the model without magnetic field. In contrast to that the model with predominantly vertical fields show downward propagation. This crucial difference in the propagation direction is also revealed in the difference in energy transported by waves for heights below 0.8 Mm. Higher up, the propagation properties show a peculiar behaviour, which require further study. Our analysis suggests that IGWs may play a significant role in the heating of the chromospheric layers of the internetwork region where horizontal fields are thought to be prevalent. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

Acoustic-gravity wave propagation characteristics in three-dimensional radiation hydrodynamic simulations of the solar atmosphere.

There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of acoustic-gravity waves at various heights. We find considerable differences between the various models. The height dependence of wave power, in particular of high-frequency waves, varies by up to two orders of magnitude between the models, and the phase difference spectra of several models show unexpected features, including ±180° phase jumps. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

Transverse motions in sunspot super-penumbral fibrils.

Sunspots have played a key role in aiding our understanding of magnetohydrodynamic (MHD) wave phenomena in the Sun's atmosphere, and it is well known they demonstrate a number of wave phenomena associated with slow MHD modes. Recent studies have shown that transverse wave modes are present throughout the majority of the chromosphere. Using high-resolution Ca II 8542 Å observations from the Swedish Solar Telescope, we provide the first demonstration that the chromospheric super-penumbral fibrils, which span out from the sunspot, also show ubiquitous transverse motions. We interpret these motions as transverse waves, in particular the MHD kink mode. We compile the statistical properties of over 2000 transverse motions to find distributions for periods and amplitudes, finding they are broadly consistent with previous observations of chromospheric transverse waves in quiet Sun fibrils. The very presence of the waves in super-penumbral fibrils raises important questions about how they are generated, and could have implications for our understanding of how MHD wave energy is transferred through the atmosphere of a sunspot. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

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'.

Oscillations observed in umbra, plage, quiet-Sun and the polarity inversion line of active region 11158 using Helioseismic Magnetic Imager/Solar Dynamics Observatory data.

Using data from the Helioseismic Magnetic Imager, we report on the amplitudes and phase relations of oscillations in quiet-Sun, plage, umbra and the polarity inversion line (PIL) of an active region NOAA#11158. We employ Fourier, wavelet and cross-correlation spectra analysis. Waves with 5 min periods are observed in umbra, PIL and plage with common phase values of ϕ(v, I) = π/2, ϕ(v, Blos) = −(π/2). In addition, ϕ(I, Blos) = π in plage are observed. These phase values are consistent with slow standing or fast standing surface sausage wave modes. The line width variations, and their phase relations with intensity and magnetic oscillations, show different values within the plage and PIL regions, which may offer a way to further differentiate wave mode mechanics. Significant Doppler velocity oscillations are present along the PIL, meaning that plasma motion is perpendicular to the magnetic field lines, a signature of Alvènic waves. A time–distance diagram along a section of the PIL shows Eastward propagating Doppler oscillations converting into magnetic oscillations; the propagation speeds range between 2 and 6 km s−1. Lastly, a 3 min wave is observed in select regions of the umbra in the magnetogram data.This article is part of the Theo Murphy meeting issue ‘High-resolution wave dynamics in the lower solar atmosphere’.

The dynamics of 3-min wavefronts and their relation to sunspot magnetic fields.

We present a study of wave processes occurring in solar active region NOAA 11131 on 10 December 2010, captured by the Solar Dynamics Observatory in the 1600 Å, 304 Å and 171 Å channels. For spectral analysis, we employed pixelized wavelet filtering together with a developed digital technique based on empirical mode decomposition. We studied the ∼3-minute wave dynamics to obtain relationships with the magnetic structuring of the underlying sunspot. We found that during development of wave trains the motion path occurred along a preferential direction, and that the broadband wavefronts can be represented as a set of separate narrowband oscillation sources. These sources become visible as the waves pass through the umbral inhomogeneities caused by the differing magnetic field inclination angles. We found the spatial and frequency fragmentation of wavefronts, and deduced that the combination of narrowband spherical and linear parts of the wavefronts provide the observed spirality. Maps of the magnetic field inclination angles confirm this assumption. We detect the activation of umbral structures as the increasing of oscillations in the sources along the front ridge. Their temporal dynamics are associated with the occurrence of umbral flashes. Spatial localization of the sources is stable over time and depends on the oscillation period. We propose that these sources are the result of wave paths along the loops extending outwards from the magnetic bundles of the umbra. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

An overall view of temperature oscillations in the solar chromosphere with ALMA.

By direct measurements of the gas temperature, the Atacama Large Millimeter/submillimeter Array (ALMA) has yielded a new diagnostic tool to study the solar chromosphere. Here, we present an overview of the brightness-temperature fluctuations from several high-quality and high-temporal-resolution (i.e. 1 and 2 s cadence) time series of images obtained during the first 2 years of solar observations with ALMA, in Band 3 and Band 6, centred at around 3 mm (100 GHz) and 1.25 mm (239 GHz), respectively. The various datasets represent solar regions with different levels of magnetic flux. We perform fast Fourier and Lomb-Scargle transforms to measure both the spatial structuring of dominant frequencies and the average global frequency distributions of the oscillations (i.e. averaged over the entire field of view). We find that the observed frequencies significantly vary from one dataset to another, which is discussed in terms of the solar regions captured by the observations (i.e. linked to their underlying magnetic topology). While the presence of enhanced power within the frequency range 3-5 mHz is found for the most magnetically quiescent datasets, lower frequencies dominate when there is significant influence from strong underlying magnetic field concentrations (present inside and/or in the immediate vicinity of the observed field of view). We discuss here a number of reasons which could possibly contribute to the power suppression at around 5.5 mHz in the ALMA observations. However, it remains unclear how other chromospheric diagnostics (with an exception of Hα line-core intensity) are unaffected by similar effects, i.e. they show very pronounced 3-min oscillations dominating the dynamics of the chromosphere, whereas only a very small fraction of all the pixels in the 10 ALMA datasets analysed here show peak power near 5.5 mHz. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

Influence of ambipolar and Hall effects on vorticity in three-dimensional simulations of magneto-convection.

This paper presents the results of the analysis of three-dimensional simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation runs are compared: one including both ambipolar diffusion and the Hall effect; one including only ambipolar diffusion and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 h of simulated solar time after the local dynamo reaches the saturation regime. We analyse the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

A new method for detecting solar atmospheric gravity waves.

Internal gravity waves have been observed in the Earth's atmosphere and oceans, on Mars and Jupiter, and in the Sun's atmosphere. Despite ample evidence for the existence of propagating gravity waves in the Sun's atmosphere, we still do not have a full understanding of their characteristics and overall role for the dynamics and energetics of the solar atmosphere. Here, we present a new approach to study the propagation of gravity waves in the solar atmosphere. It is based on calculating the three-dimensional cross-correlation function between the vertical velocities measured at different heights. We apply this new method to a time series of co-spatial and co-temporal Doppler images obtained by SOHO/MDI and Hinode/SOT as well as to simulations of upward propagating gravity wave-packets. We show some preliminary results and outline future developments. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

Spectropolarimetric fluctuations in a sunspot chromosphere.

The instrumental advances made in this new era of 4 m class solar telescopes with unmatched spectropolarimetric accuracy and sensitivity will enable the study of chromospheric magnetic fields and their dynamics with unprecedented detail. In this regard, spectropolarimetric diagnostics can provide invaluable insight into magneto-hydrodynamic (MHD) wave processes. MHD waves and, in particular, Alfvénic fluctuations associated with particular wave modes were recently recognized as important mechanisms not only for the heating of the outer layers of the Sun's atmosphere and the acceleration of the solar wind, but also for the elemental abundance anomaly observed in the corona of the Sun and other Sun-like stars (also known as first ionization potential) effect. Here, we take advantage of state-of-the-art and unique spectropolarimetric Interferometric BIdimensional Spectrometer observations to investigate the relation between intensity and circular polarization (CP) fluctuations in a sunspot chromosphere. Our results show a clear link between the intensity and CP fluctuations in a patch which corresponds to a narrow range of magnetic field inclinations. This suggests the presence of Alfvénic perturbations in the sunspot. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

High-frequency oscillations in small chromospheric bright features observed with Atacama Large Millimetre/Submillimetre Array.

We report detection of oscillations in brightness temperature, size and horizontal velocity of three small bright features in the chromosphere of a plage/enhanced-network region. The observations, which were taken with high temporal resolution (i.e. 2 s cadence) with the Atacama large millimetre/ submillimetre array (ALMA) in Band 3 (centred at 3 mm; 100 GHz), exhibit three small-scale features with oscillatory behaviour with different, but overlapping, distributions of period on the order of, on average, 90 ± 22 s, 110 ± 12 s and 66 ± 23 s, respectively. We find anti-correlations between perturbations in brightness, temperature and size of the three features, which suggest the presence of fast sausage-mode waves in these small structures. In addition, the detection of transverse oscillations (although with a larger uncertainty) may also suggest the presence of Alfvénic oscillations which are likely representative of kink waves. This work demonstrates the diagnostic potential of high-cadence observations with ALMA for detecting high-frequency magnetohydrodynamic waves in the solar chromosphere. Such waves can potentially channel a vast amount of energy into the outer atmosphere of the Sun. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.

On the effect of oscillatory phenomena on Stokes inversion results.

Stokes inversion codes are crucial in returning properties of the solar atmosphere, such as temperature and magnetic field strength. However, the success of such algorithms to return reliable values can be hindered by the presence of oscillatory phenomena within magnetic wave guides. Returning accurate parameters is crucial to both magnetohydrodynamics (MHD) studies and solar physics in general. Here, we employ a simulation featuring propagating MHD waves within a flux tube with a known driver and atmospheric parameters. We invert the Stokes profiles for the 6301 Å and 6302 Å line pair emergent from the simulations using the well-known Stokes Inversions from Response functions code to see if the atmospheric parameters can be returned for typical spatial resolutions at ground-based observatories. The inversions return synthetic spectra comparable to the original input spectra, even with asymmetries introduced in the spectra from wave propagation in the atmosphere. The output models from the inversions match closely to the simulations in temperature, line-of-sight magnetic field and line-of-sight velocity within typical formation heights of the inverted lines. Deviations from the simulations are seen away from these height regions. The inversions results are less accurate during passage of the waves within the line formation region. The original wave period could be recovered from the atmosphere output by the inversions, with empirical mode decomposition performing better than the wavelet approach in this task. This article is part of the Theo Murphy meeting issue 'High-resolution wave dynamics in the lower solar atmosphere'.