Sergiy Shelyag

NUMERICAL MODELING OF FOOTPOINT-DRIVEN MAGNETO-ACOUSTIC WAVE PROPAGATION IN A LOCALIZED SOLAR FLUX TUBE

In this paper, we present and discuss results of two-dimensional simulations of linear and nonlinear magneto-acoustic wave propagation through an open magnetic flux tube embedded in the solar atmosphere expanding from the photosphere through to the transition region and into the low corona. Our aim is to model and analyze the response of such a magnetic structure to vertical and horizontal periodic motions originating in the photosphere. To carry out the simulations, we employed our MHD code SAC (Sheffield Advanced Code). A combination of the VALIIIC and McWhirter solar atmospheres and coronal density profiles were used as the background equilibrium model in the simulations. Vertical and horizontal harmonic sources, located at the footpoint region of the open magnetic flux tube, are incorporated in the calculations, to excite oscillations in the domain of interest. To perform the analysis we have constructed a series of time-distance diagrams of the vertical and perpendicular components of the velocity with respect to the magnetic field lines at each height of the computational domain. These time-distance diagrams are subject to spatio-temporal Fourier transforms allowing us to build ω-k dispersion diagrams for all of the simulated regions in the solar atmosphere. This approach makes it possible to compute the phase speeds of waves propagating throughout the various regions of the solar atmosphere model. We demonstrate the transformation of linear slow and fast magneto-acoustic wave modes into nonlinear ones, i.e., shock waves, and also show that magneto-acoustic waves with a range of frequencies efficiently leak through the transition region into the solar corona. It is found that the waves interact with the transition region and excite horizontally propagating surface waves along the transition region for both types of drivers. Finally, we estimate the phase speed of the oscillations in the solar corona and compare it with the phase speed derived from observations.

Acoustic wave propagation in the solar sub-photosphere with localised magnetic field concentration: effect of magnetic tension

Aims - We analyse numerically the propagation and dispersion of acoustic waves in the solar-like sub-photosphere with localised non-uniform magnetic field concentrations, mimicking sunspots with various representative magnetic field configurations. n nMethods - Numerical simulations of wave propagation through the solar sub-photosphere with a localised magnetic field concentration are carried out using SAC, which solves the MHD equations for gravitationally stratified plasma. The initial equilibrium density and pressure stratifications are derived from a standard solar model. Acoustic waves are generated by a source located at the height corresponding approximately to the visible surface of the Sun. By means of local helioseismology we analyse the response of vertical velocity at the level corresponding to the visible solar surface to changes induced by magnetic field in the interior. n nResults - The results of numerical simulations of acoustic wave propagation and dispersion in the solar sub-photosphere with localised magnetic field concentrations of various types are presented. Time-distance diagrams of the vertical velocity perturbation at the level corresponding to the visible solar surface show that the magnetic field perturbs and scatters acoustic waves and absorbs the acoustic power of the wave packet. For the weakly magnetised case, the effect of magnetic field is mainly thermodynamic, since the magnetic field changes the temperature stratification. However, we observe the signature of slow magnetoacoustic mode, propagating downwards, for the strong magnetic field cases.

Alfvén Waves in Simulations of Solar Photospheric Vortices

Using advanced numerical magneto-hydrodynamic simulations of the magnetized solar photosphere, including non-gray radiative transport and a non-ideal equation of state, we analyze plasma motions in photospheric magnetic vortices. We demonstrate that apparent vortex-like motions in photospheric magnetic field concentrations do not exhibit tornado-like behavior or a bath-tub effect. While at each time instance the velocity field lines in the upper layers of the solar photosphere show swirls, the test particles moving with the time-dependent velocity field do not demonstrate such structures. Instead, they move in a wave-like fashion with rapidly changing and oscillating velocity field, determined mainly by magnetic tension in the magnetized intergranular downflows. Using time-distance diagrams, we identify horizontal motions in the magnetic flux tubes as torsional Alfven perturbations propagating along the nearly vertical magnetic field lines with local Alfven speed.

ELLERMAN BOMBS—EVIDENCE FOR MAGNETIC RECONNECTION IN THE LOWER SOLAR ATMOSPHERE

The presence of photospheric magnetic reconnection has long been thought to give rise to short and impulsive events, such as Ellerman bombs (EBs) and Type II spicules. In this article, we combine high-resolution, high-cadence observations from the Interferometric BIdimensional Spectrometer and Rapid Oscillations in the Solar Atmosphere instruments at the Dunn Solar Telescope, National Solar Observatory, New Mexico, with co-aligned Solar Dynamics Observatory Atmospheric Imaging Assembly and Hinode Solar Optical Telescope (SOT) data to observe small-scale events situated within an active region. These data are then compared with state-of-the-art numerical simulations of the lower atmosphere made using the MURaM code. It is found that brightenings, in both the observations and the simulations, of the wings of the Hα line profile, interpreted as EBs, are often spatially correlated with increases in the intensity of the Fe I λ6302.5 line core. Bipolar regions inferred from Hinode/SOT magnetic field data show evidence of flux cancellation associated, co-spatially, with these EBs, suggesting that magnetic reconnection could be a driver of these high-energy events. Through the analysis of similar events in the simulated lower atmosphere, we are able to infer that line profiles analogous to the observations occur co-spatially with regions of strong opposite-polarity magnetic flux. These observed events and their simulated counterparts are interpreted as evidence of photospheric magnetic reconnection at scales observable using current observational instrumentation.

Magnetohydrodynamic code for gravitationally-stratified media

Aims. We describe a newly-developed magnetohydrodynamic (MHD) code with the capacity to simulate the interaction of any arbitrary perturbation (i.e., not necessarily limited to the linearised limit) with a magnetohydrostatic equilibrium background. Methods. By rearranging the terms in the system of MHD equations and explicitly taking into account the magnetohydrostatic equilibrium condition, we define the equations governing the perturbations that describe the deviations from the background state of plasma for the density, internal energy and magnetic field. We found it was advantageous to use this modified form of the MHD equations for numerical simulations of physical processes taking place in a stable gravitationally-stratified plasma. The governing equations are implemented in a novel way in the code. Sub-grid diffusion and resistivity are applied to ensure numerical stability of the computed solution of the MHD equations. We apply a fourth-order central difference scheme to calculate the spatial derivatives, and implement an arbitrary Runge-Kutta scheme to advance the solution in time. Results. We have built the proposed method, suitable for strongly-stratified magnetised plasma, on the base of the well-documented Versatile Advection Code (VAC) and performed a number of one- and multi-dimensional hydrodynamic and MHD tests to demonstrate the feasibility and robustness of the code for applications to astrophysical plasmas.

Forward modeling of acoustic wave propagation in the quiet solar subphotosphere

The results of numerical simulations of acoustic wave propagation and dispersion in the nonmagnetic solar subphotosphere are presented. Initial equilibrium density and pressure stratifications are taken from a standard solar model but modified to suppress convective instabilities in fully compressible two-dimensional ideal hydrodynamical modeling. Acoustic waves are generated by sources located below the height corresponding to the visible solar surface. The dynamic response of the solar interior to two acoustic source types, namely a harmonic source and one representing downward-propagating photospheric plumes, is studied. A large number of randomly distributed localized cooling sources with random amplitudes is also introduced. The methods used to analyze the simulation data are similar to ones used in observational studies in local helioseismology. Time-distance diagrams of the pressure and vertical velocity perturbations at the level corresponding to the solar surface show the appearance of wave packets propagating with different speeds, which are reflected at different depths beneath the subphotosphere. The (ω, kh) power spectra, derived from the vertical velocity data, show the existence of g-, f-, and p-modes; p-mode ridges are identifiable up to high radial orders of n ≈ 11; g-modes appear in the simulations, unlike in the real Sun, where they cannot propagate in the convectively unstable solar subphotosphere. Cross-correlation analysis of vertical velocity perturbations shows a good correspondence with the observed time-distance helioseismic data for quiet Sun. Thus, the ability of the implemented approach of forward modeling to investigate propagation of acoustic, internal, and surface gravity waves in a realistic solar interior model is shown.

Alfvén Wave heating of the solar chromosphere : 1.5D models

Physical processes that may lead to solar chromospheric heating are analyzed using high-resolution 1.5D non-ideal MHD modeling. We demonstrate that it is possible to heat the chromospheric plasma by direct resistive dissipation of high-frequency Alfven waves through Pedersen resistivity. However, this is unlikely to be sufficient to balance radiative and conductive losses unless unrealistic field strengths or photospheric velocities are used. The precise heating profile is determined by the input driving spectrum, since in 1.5D there is no possibility of Alfven wave turbulence. The inclusion of the Hall term does not affect the heating rates. If plasma compressibility is taken into account, shocks are produced through the ponderomotive coupling of Alfven waves to slow modes and shock heating dominates the resistive dissipation. In 1.5D shock coalescence amplifies the effects of shocks, and for compressible simulations with realistic driver spectra, the heating rate exceeds that required to match radiative and conductive losses. Thus, while the heating rates for these 1.5D simulations are an overestimate, they do show that ponderomotive coupling of Alfven waves to sound waves is more important in chromospheric heating than Pedersen dissipation through ion–neutral collisions.

The Dynamics of Rapid Redshifted and Blueshifted Excursions in the Solar Hα Line

We analyze high temporal and spatial resolution time-series of spectral scans of the Hα line obtained with the CRisp Imaging SpectroPolarimeter instrument mounted on the Swedish Solar Telescope. The data reveal highly dynamic, dark, short-lived structures known as Rapid Redshifted and Blueshifted Excursions (RREs, RBEs) that are on-disk absorption features observed in the red and blue wings of spectral lines formed in the chromosphere. We study the dynamics of RREs and RBEs by tracking their evolution in space and time, measuring the speed of the apparent motion, line of sight (LOS) Doppler velocity, and transverse velocity of individual structures. A statistical study of their measured properties shows that RREs and RBEs have similar occurrence rates, lifetimes, lengths, and widths. They also display non-periodic, nonlinear transverse motions perpendicular to their axes at speeds of 4–31 km s−1. Furthermore, both types of structures either appear as high speed jets and blobs that are directed outwardly from a magnetic bright point with speeds of 50–150 km s−1, or emerge within a few seconds. A study of the different velocity components suggests that the transverse motions along the LOS of the chromospheric flux tubes are responsible for the formation and appearance of these redshifted/blueshifted structures. The short lifetime and fast disappearance of the RREs/RBEs suggests that, similar to type II spicules, they are rapidly heated to transition region or even coronal temperatures. We speculate that the Kelvin–Helmholtz instability triggered by observed transverse motions of these structures may be a viable mechanism for their heating.

STELLAR SURFACE MAGNETO-CONVECTION AS A SOURCE OF ASTROPHYSICAL NOISE. I. MULTI-COMPONENT PARAMETERIZATION OF ABSORPTION LINE PROFILES

We outline our techniques to characterize photospheric granulation as an astrophysical noise source. A four-component parameterization of granulation is developed that can be used to reconstruct stellar line asymmetries and radial velocity shifts due to photospheric convective motions. The four components are made up of absorption line profiles calculated for granules, magnetic intergranular lanes, non-magnetic intergranular lanes, and magnetic bright points at disk center. These components are constructed by averaging Fe I 6302 A magnetically sensitive absorption line profiles output from detailed radiative transport calculations of the solar photosphere. Each of the four categories adopted is based on magnetic field and continuum intensity limits determined from examining three-dimensional magnetohydrodynamic simulations with an average magnetic flux of 200 G. Using these four-component line profiles we accurately reconstruct granulation profiles, produced from modeling 12 × 12 Mm2 areas on the solar surface, to within ~ ±20 cm s–1 on a ~100 m s–1 granulation signal. We have also successfully reconstructed granulation profiles from a 50 G simulation using the parameterized line profiles from the 200 G average magnetic field simulation. This test demonstrates applicability of the characterization to a range of magnetic stellar activity levels.

Ellerman bombs with jets: cause and effect

Ellerman Bombs (EBs) are thought to arise as a result of photospheric magnetic reconnection. We use data from the Swedish 1 m Solar Telescope to study EB events on the solar disk and at the limb. Both data sets show that EBs are connected to the foot points of forming chromospheric jets. The limb observations show that a bright structure in the Hα blue wing connects to the EB initially fueling it, leading to the ejection of material upwards. The material moves along a loop structure where a newly formed jet is subsequently observed in the red wing of Hα. In the disk data set, an EB initiates a jet which propagates away from the apparent reconnection site within the EB flame. The EB then splits into two, with associated brightenings in the inter-granular lanes. Micro-jets are then observed, extending to 500 km with a lifetime of a few minutes. Observed velocities of the micro-jets are approximately 5–10 km s−1, while their chromospheric counterparts range from 50 to 80 km s−1. MURaM simulations of quiet Sun reconnection show that micro-jets with properties similar to those of the observations follow the line of reconnection in the photosphere, with associated Hα brightening at the location of increased temperature.