-
Thermal and kinetic coronal rain diagnostics with MgII h & k lines
Authors:
M. Kriginsky,
R. Oliver
Abstract:
We performed the temperature diagnostics using inversions of data from nine spectroscopic observations obtained with the IRIS spectrograph in the MgII h & k lines. The sensitivity to the temperature of the emission peaks of these lines was exploited to determine the temperature of the coronal rain plasma using inversions of the spectroscopic profiles. Additional relationships between different spe…
▽ More
We performed the temperature diagnostics using inversions of data from nine spectroscopic observations obtained with the IRIS spectrograph in the MgII h & k lines. The sensitivity to the temperature of the emission peaks of these lines was exploited to determine the temperature of the coronal rain plasma using inversions of the spectroscopic profiles. Additional relationships between different spectral features of these lines, derived from the use of 3D radiative transfer line synthesis applied to simulations, were employed in order to derive the line-of-sight (LoS) velocities in different parts of the coronal rain plasma. For the first time, spectroscopic inversions of coronal rain were successfully performed. Temperatures derived from the inversions yield coronal rain clump temperatures at the formation height of the emission peaks of the MgII h & k lines in the range between 5000 and 7000 K. This narrow range of values remains consistent among all the different observations used in this work. We obtained LoS velocities of up to 40 km/s, which are consistent with the motion of the plasma being mostly constrained to the plane of the sky, as the coronal rain was mostly detected shortly after its formation and the observations took place in the disc. Furthermore, velocity diagnostics led to the detection of larger velocities at higher layers of the coronal rain plasma in some cases. This increased velocity seems to indicate that at some point (at least) during the fall of coronal rain clumps towards the chromosphere, the material in the upper part of the coronal rain plasma is falling with greater velocity than the material below it. The conditions of the temperature and density of the coronal rain plasma where the Mg II h line forms appear to be slightly different that those of the Mg II k line, with the former found at slightly colder and denser parts of the plasma.
△ Less
Submitted 26 January, 2024;
originally announced January 2024.
-
Temperature diagnostics of chromospheric fibrils
Authors:
M. Kriginsky,
R. Oliver,
D. Kuridze
Abstract:
Context. Chromospheric fibrils are thin and elongated structures that connect nearby photospheric magnetic field concentrations of opposite polarities. Aims. We assess the possibilities and drawbacks related to the use of current instrumentation and inversion techniques to infer the thermodynamic structure of chromospheric fibrils. Methods. We employed spectroscopic observations obtained in the Ca…
▽ More
Context. Chromospheric fibrils are thin and elongated structures that connect nearby photospheric magnetic field concentrations of opposite polarities. Aims. We assess the possibilities and drawbacks related to the use of current instrumentation and inversion techniques to infer the thermodynamic structure of chromospheric fibrils. Methods. We employed spectroscopic observations obtained in the Ca ii 854.2 nm line with the CRISP instrument at the Swedish 1-m Solar Telescope and in coordination with observations in the ultraviolet Mg ii h & k lines taken with the IRIS satellite. We studied the temperature sensitivity of these chromospheric lines to properly invert their spectral profiles with the Stockholm inversion Code and determine the temperature, line-of-sight velocity, and microturbulent velocity of manually traced chromospheric fibrils present in the field of view. Results. Fibril-like structures show a very particular dependence of their temperature as a function of the position along their length. Their temperatures at the detected footpoints are, on average, 300 K higher than the temperature at the midpoint. The temperature variation appears to be almost symmetrical in shape, with partially traced fibrils showing a similar trend for the temperature variation. Additionally, the response of the Ca ii 854.2 nm line core to variations of the temperature for the inverted models of the atmosphere in fibril areas seems to be insufficient to properly resolve the aforementioned temperature structure. Only the addition of more temperature sensitive lines such as the Mg ii h & k lines would make it possible to properly infer the thermodynamic properties of chromospheric fibrils. Comparisons between the results obtained here and in previous studies focused on bright Ca ii K fibrils yield great similarities between these structures in terms of their temperature.
△ Less
Submitted 16 February, 2023;
originally announced February 2023.
-
Large ion-neutral drift velocities and plasma heating in partially ionized coronal rain blobs
Authors:
David Martínez-Gómez,
Ramón Oliver,
Elena Khomenko,
Manuel Collados
Abstract:
In this paper we present a numerical study of the dynamics of partially ionized coronal rain blobs. We use a two-fluid model to perform a high-resolution 2D simulation that takes into account the collisional interaction between the charged and neutral particles contained in the plasma. We follow the evolution of a cold plasma condensation as it falls through an isothermal vertically stratified atm…
▽ More
In this paper we present a numerical study of the dynamics of partially ionized coronal rain blobs. We use a two-fluid model to perform a high-resolution 2D simulation that takes into account the collisional interaction between the charged and neutral particles contained in the plasma. We follow the evolution of a cold plasma condensation as it falls through an isothermal vertically stratified atmosphere that represents the much hotter and lighter solar corona. We study the consequences of the different degrees of collisional coupling that are present in the system. On the one hand, we find that at the dense core of the blob there is a very strong coupling and the charged and neutral components of the plasma behave as a single fluid, with negligible drift velocities (of a few cm s^-1). On the other hand, at the edges of the blob the coupling is much weaker and larger drift velocities (of the order of 1 km s^-1) appear. In addition, frictional heating causes large increases of temperature at the transition layers between the blob and the corona. For the first time we show that such large drift velocities and temperature enhancements can develop as a consequence of ion-neutral decoupling associated to coronal rain dynamics. This can lead to enhanced emission coming from the plasma at the coronal rain-corona boundary, which possesses transition region temperature.
△ Less
Submitted 10 November, 2022;
originally announced November 2022.
-
Dark off-limb gap: manifestation of temperature minimum and dynamic nature of the chromosphere
Authors:
D. Kuridze,
P. Heinzel,
J. Koza,
R. Oliver
Abstract:
We study off-limb emission of the lower solar atmosphere using high-resolution imaging spectroscopy in the H$β$ and Ca II 8542 Å lines obtained with the CHROMospheric Imaging Spectrometer (CHROMIS) and the CRisp Imaging SpectroPolarimeter (CRISP) on the Swedish 1-m Solar Telescope. The H$β$ line wing images show the dark intensity gap between the photospheric limb and chromosphere which is absent…
▽ More
We study off-limb emission of the lower solar atmosphere using high-resolution imaging spectroscopy in the H$β$ and Ca II 8542 Å lines obtained with the CHROMospheric Imaging Spectrometer (CHROMIS) and the CRisp Imaging SpectroPolarimeter (CRISP) on the Swedish 1-m Solar Telescope. The H$β$ line wing images show the dark intensity gap between the photospheric limb and chromosphere which is absent in the Ca II images. We calculate synthetic spectra of the off-limb emissions with the RH code in the one-dimension spherical geometry and find good agreement with the observations. The analysis of synthetic line profiles shows that the gap in the H$β$ line wing images maps the temperature minimum region between the photosphere and chromosphere due to the well known opacity and emissivity gap of Balmer lines in this layer. However, observed gap is detected farther from the line core in the outer line wing positions than in the synthetic profiles. We found that an increased microturbulence in the model chromosphere is needed to reproduce the dark gap in the outer line wing, suggesting that observed H$β$ gap is the manifestation of the temperature minimum and the dynamic nature of the solar chromosphere. The temperature minimum produces a small enhancement in synthetic Ca II line-wing intensities. Observed off-limb Ca II line-wing emissions show similar enhancement below temperature minimum layer near the edge of the photospheric limb.
△ Less
Submitted 30 August, 2022;
originally announced August 2022.
-
The cyclic behaviour in the N-S asymmetry of sunspots and solar plages for the period 1910 to 1937 using data from Ebro catalogues
Authors:
V. de Paula,
J. J. Curto,
R. Oliver
Abstract:
The heliophysics catalogues published by the Ebro Observatory during 1910--1937 have been converted into a digital format in order to provide the data for computational processing. This has allowed us to study in detail the North-South (N-S) asymmetry of solar activity in that period, focusing on two different structures located at two different layers of the solar atmosphere: sunspots (Photospher…
▽ More
The heliophysics catalogues published by the Ebro Observatory during 1910--1937 have been converted into a digital format in order to provide the data for computational processing. This has allowed us to study in detail the North-South (N-S) asymmetry of solar activity in that period, focusing on two different structures located at two different layers of the solar atmosphere: sunspots (Photosphere) and solar plages (Chromosphere). The examination of the absolute and normalised N-S asymmetry indices in terms of their monthly sum of occurrences and areas has made possible to find out a cyclic behaviour in the solar activity, in which the preferred hemisphere changes systematically with a global period of 7.9 $\pm$ 0.2 yr. In order to verify and quantify accurately this periodicity and study its prevalence in time, we employed the RGO-USAF/NOAA sunspot data series during 1874--2016. Then, we examined each absolute asymmetry index time series through different techniques as the power spectrum analysis, the Complete Ensemble Empirical Mode Decomposition With Adaptive Noise algorithm, or the Morlet wavelet transform. The combined results reveal a cyclic behaviour at different time scales, consisting in two quite stable periodicities of 1.47 $\pm$ 0.02 yr and 3.83 $\pm$ 0.06 yr, which coexist with another three discontinuous components with more marked time-varying periods with means of 5.4 $\pm$ 0.2 yr, 9.0 $\pm$ 0.2 yr, and 12.7 $\pm$ 0.3 yr. Moreover, during 1910--1937, only two dominant signals with averaged periods of 4.10 $\pm$ 0.04 yr and 7.57 $\pm$ 0.03 yr can be clearly observed. Finally, in both signals, periods are slightly longer for plages in comparison with sunspots.
△ Less
Submitted 17 February, 2022;
originally announced February 2022.
-
Construction of coronal hole and active region magnetohydrostatic solutions in two dimensions: Force and energy balance
Authors:
J. Terradas,
R. Soler,
R. Oliver,
P. Antolin,
I. Arregui,
M. Luna,
I. Piantschitsch,
E. Soubrié,
J. L.,
Ballester
Abstract:
Coronal holes and active regions are typical magnetic structures found in the solar atmosphere. We propose several magnetohydrostatic equilibrium solutions that are representative of these structures in two dimensions. Our models include the effect of a finite plasma-$β$ and gravity, but the distinctive feature is that we incorporate a thermal structure with properties similar to those reported by…
▽ More
Coronal holes and active regions are typical magnetic structures found in the solar atmosphere. We propose several magnetohydrostatic equilibrium solutions that are representative of these structures in two dimensions. Our models include the effect of a finite plasma-$β$ and gravity, but the distinctive feature is that we incorporate a thermal structure with properties similar to those reported by observations. We developed a semi-analytical method to compute the equilibrium configuration. Using this method, we obtain cold and under-dense plasma structures in open magnetic fields representing coronal holes, while in closed magnetic configurations, we achieve the characteristic hot and over-dense plasma arrangements of active regions. Although coronal holes and active regions seem to be antagonistic structures, we find that they can be described using a common thermal structure that depends on the flux function. In addition to the force balance, the energy balance is included in the constructed models using an a posteriori approach. From the two-dimensional computation of thermal conduction and radiative losses in our models, we infer the required heating function to achieve energy equilibrium. We find that the temperature dependence on height is an important parameter that may prevent the system from accomplishing thermal balance at certain spatial locations. The implications of these results are discussed in detail.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
One dimensional prominence threads: I. Equilibrium models
Authors:
J. Terradas,
M. Luna,
R. Soler,
R. Oliver,
M. Carbonell,
J. L Ballester
Abstract:
Threads are the building blocks of solar prominences and very often show longitudinal oscillatory motions that are strongly attenuated with time. The damping mechanism responsible for the reported oscillations is not fully understood yet. To understand the oscillations and damping of prominence threads it is mandatory to investigate first the nature of the equilibrium solutions that arise under st…
▽ More
Threads are the building blocks of solar prominences and very often show longitudinal oscillatory motions that are strongly attenuated with time. The damping mechanism responsible for the reported oscillations is not fully understood yet. To understand the oscillations and damping of prominence threads it is mandatory to investigate first the nature of the equilibrium solutions that arise under static conditions and under the presence of radiative losses, thermal conduction and background heating. This provides the basis to calculate the eigenmodes of the thread models. The nonlinear ordinary differential equations for hydrostatic and thermal equilibrium under the presence of gravity are solved using standard numerical techniques and simple analytical expressions are derived under certain approximations. The solutions to the equations represent a prominence thread, i.e., a dense and cold plasma region of a certain length that connects with the corona through a prominence corona transition region (PCTR). The solutions can also match with a chromospheric-like layer if a spatially dependent heating function localised around the footpoints is considered. We have obtained static solutions representing prominence threads and have investigated in detail the dependence of these solutions on the different parameters of the model. Among other results, we have shown that multiple condensations along a magnetic field line are possible, and that the effect of partial ionisation in the model can significantly modify the thermal balance in the thread and therefore their length. This last parameter is also shown to be comparable to that reported in the observations when the radiative losses are reduced for typical thread temperatures.
△ Less
Submitted 11 June, 2021;
originally announced June 2021.
-
Magnetic field inference in active region coronal loops using coronal rain clumps
Authors:
M. Kriginsky,
R. Oliver,
P. Antolin,
D. Kuridze,
N. Freij
Abstract:
Aims. We aim to infer information about the magnetic field in the low solar corona from coronal rain clumps using high-resolution spectropolarimetric observations in the Ca ii 8542 A line obtained with the Swedish 1-m Solar Telescope. Methods. The weak-field approximation (WFA) provides a simple tool to obtain the line-of-sight component of the magnetic field from spectropolarimetric observations.…
▽ More
Aims. We aim to infer information about the magnetic field in the low solar corona from coronal rain clumps using high-resolution spectropolarimetric observations in the Ca ii 8542 A line obtained with the Swedish 1-m Solar Telescope. Methods. The weak-field approximation (WFA) provides a simple tool to obtain the line-of-sight component of the magnetic field from spectropolarimetric observations. We adapted a method developed in a previous paper in order to assess the different conditions that must be satisfied in order to properly use the WFA for the data at hand. We also made use of velocity measurements in order to estimate the plane-of-the-sky magnetic field component, so that the magnetic field vector could be inferred. Results. We have inferred the magnetic field vector from a data set totalling 100 spectral scans in the Ca ii 8542 A line, containing an off-limb view of the lower portion of catastrophically cooled coronal loops in an active region. Our results, albeit limited by the cadence and signal-to-noise ratio of the data, suggest that magnetic field strengths of hundreds of Gauss, even reaching up to 1000 G, are omnipresent at coronal heights below 9 Mm from the visible limb. Our results are also compatible with the presence of larger magnetic field values such as those reported by previous works. However, for large magnetic fields, the Doppler width from coronal rain is not that much larger than the Zeeman width, thwarting the application of the WFA. Furthermore, we have determined the temperature, T , and microturbulent velocity, $ξ$, of coronal rain clumps and off-limb spicules present in the same data set, and we have found that the former ones have narrower T and $ξ$ distributions, their average temperature is similar, and coronal rain has microturbulent velocities smaller than those of spicules
△ Less
Submitted 7 April, 2021;
originally announced April 2021.
-
Gravitational instability of solar prominence threads I. Curved magnetic fields without dips
Authors:
A. Adrover-González,
J. Terradas,
R. Oliver,
M. Carbonell
Abstract:
Prominence threads are dense and cold structures lying on curved magnetic fields that can be suspended in the solar atmosphere against gravity. The gravitational stability of threads, in the absence of non-ideal effects, is comprehensively investigated in the present work by means of an elementary but effective model. Based on purely hydrodynamic equations in one spatial dimension and applying lin…
▽ More
Prominence threads are dense and cold structures lying on curved magnetic fields that can be suspended in the solar atmosphere against gravity. The gravitational stability of threads, in the absence of non-ideal effects, is comprehensively investigated in the present work by means of an elementary but effective model. Based on purely hydrodynamic equations in one spatial dimension and applying line-tying conditions at the footpoints of the magnetic field lines, we derive analytical expressions for the different feasible equilibria and the corresponding frequencies of oscillation. We find that the system allows for stable and unstable equilibrium solutions subject to the initial position of the thread, its density contrast and length, and the total length of the magnetic field lines. The transition between the two types of solutions is produced at specific bifurcation points that have been determined analytically in some particular cases. When the thread is initially at the top of the concave magnetic field, that is at the apex, we find a supercritical pitchfork bifurcation, while for a shifted initial thread position with respect to this point the symmetry is broken and the system is characterised by an S-shaped bifurcation. The plain results presented in this paper shed new light on the behaviour of threads in curved magnetic fields under the presence of gravity and help to interpret more complex numerical magnetohydrodynamics (MHD) simulations about similar structures.
△ Less
Submitted 22 March, 2021;
originally announced March 2021.
-
Semi-empirical models of spicule from inversion of Ca II 8542 Å line
Authors:
D. Kuridze,
H. Socas-Navarro,
J. Koza,
R. Oliver
Abstract:
We study a solar spicule observed off-limb using high-resolution imaging spectroscopy in the Ca II 8542 Å line obtained with the CRisp Imaging SpectroPolarimeter (CRISP) on the Swedish 1-m Solar Telescope. Using a new version of the non-LTE code NICOLE specifically developed for this problem we invert the spicule single- and double-peak line profiles. This new version considers off-limb geometry a…
▽ More
We study a solar spicule observed off-limb using high-resolution imaging spectroscopy in the Ca II 8542 Å line obtained with the CRisp Imaging SpectroPolarimeter (CRISP) on the Swedish 1-m Solar Telescope. Using a new version of the non-LTE code NICOLE specifically developed for this problem we invert the spicule single- and double-peak line profiles. This new version considers off-limb geometry and computes atomic populations by solving the 1D radiative transfer assuming a vertical stratification. The inversion proceeds by fitting the observed spectral profiles at 14 different heights with synthetic profiles computed in the model by solving the radiative transfer problem along its length. Motivated by the appearance of double-peak Ca II 8542 Å spicule profiles, which exhibit two distinct emission features well separated in wavelength, we adopt a double-component scenario. We start from the ansatz that the spicule parameters are practically constant along the spicule axis for each component, except for a density drop. Our results support this ansatz by attaining very good fits to the entire set of 14$\times$4 profiles (14 heights and 4 times). We show that the double-component model with uniform temperature of 9560 K, exponential decrease of density with a height scale of $1000-2000$ km, and the counter-oriented line-of-sight velocities of components reproduce the double-peak line profiles at all spicule segments well. Analyses of the numerical response function reveals the necessity of the inversions of spectra at multiple height positions to obtain height-dependent, degeneracy-free reliable model with a limited number of free parameters.
△ Less
Submitted 7 December, 2020;
originally announced December 2020.
-
Ubiquitous hundred-Gauss magnetic fields in solar spicules
Authors:
M. Kriginsky,
R. Oliver,
N. Freij,
D. Kuridze,
A. Asensio Ramos,
P. Antolin
Abstract:
Aims. We use high-resolution spectropolarimetric observations in the Ca ii 8542 A line obtained with the SST to study the magnetic field in solar spicules. Methods. The equations that result from the application of the Weak Field Approximation (WFA) to the radiative transfer equations are used to infer the LOS component of the magnetic field (BLOS). Two restrictive conditions are imposed on the St…
▽ More
Aims. We use high-resolution spectropolarimetric observations in the Ca ii 8542 A line obtained with the SST to study the magnetic field in solar spicules. Methods. The equations that result from the application of the Weak Field Approximation (WFA) to the radiative transfer equations are used to infer the LOS component of the magnetic field (BLOS). Two restrictive conditions are imposed on the Stokes I and V profiles at each pixel before they can be used in a Bayesian inversion to compute its BLOS. Results. The LOS magnetic field component has been inferred in six data sets totalling 448 spectral scans in the Ca ii 8542 A line and containing both active region and quiet Sun areas, with values of hundreds of G being abundantly inferred. There seems to be no difference, from the statistical point of view, between the magnetic field strength of spicules in the quiet Sun or near an active region. On the other hand, the BLOS distributions present smaller values on the disk than off-limb, a fact that can be explained by the effect of superposition on the chromosphere of on-disk structures. We find that spicules in the vicinity of a sunspot have a magnetic field polarity (i.e. north or south) equal to that of the sunspot. This paper also contains an analysis of the effect of off-limb overlapping structures on the observed Stokes I and V parameters and the BLOS obtained from the WFA. It is found that this value is equal to or smaller than the largest LOS magnetic field components of the two structures. In addition, using random BLOS, Doppler velocities and line intensities of these two structures leads in ~ 50% of the cases to Stokes I and V parameters unsuitable to be used with the WFA. Conclusions. Our results present a scarcity of LOS magnetic field components smaller than some 50 G, which must not be taken as evidence against the existence of these magnetic field strengths in spicules.
△ Less
Submitted 1 June, 2020;
originally announced June 2020.
-
Spectral Characteristics and Formation Height of Off-Limb Flare Ribbons
Authors:
D. Kuridze,
M. Mathioudakis,
P. Heinzel,
J. Koza,
H. Morgan,
R. Oliver,
A. F. Kowalski,
J. C. Allred
Abstract:
Flare ribbons are bright manifestations of flare energy dissipation in the lower solar atmosphere. For the first time, we report on high-resolution imaging spectroscopy observations of flare ribbons situated off-limb in the H$β$ and Ca II 8542 Å lines and make a detailed comparison with radiative hydrodynamic simulations. Observations of the X8.2-class solar flare SOL2017-09-10T16:06 UT obtained w…
▽ More
Flare ribbons are bright manifestations of flare energy dissipation in the lower solar atmosphere. For the first time, we report on high-resolution imaging spectroscopy observations of flare ribbons situated off-limb in the H$β$ and Ca II 8542 Å lines and make a detailed comparison with radiative hydrodynamic simulations. Observations of the X8.2-class solar flare SOL2017-09-10T16:06 UT obtained with the Swedish Solar Telescope reveal bright horizontal emission layers in H$β$ line wing images located near the footpoints of the flare loops. The apparent separation between the ribbon observed in the H$β$ wing and the nominal photospheric limb is about 300 - 500 km. The Ca II 8542 Å line wing images show much fainter ribbon emissions located right on the edge of the limb, without clear separation from the limb. RADYN models are used to investigate synthetic spectral line profiles for the flaring atmosphere, and good agreement is found with the observations. The simulations show that, towards the limb, where the line of sight is substantially oblique with respect to the vertical direction, the flaring atmosphere model reproduces the high contrast of the off-limb H$β$ ribbons and their significant elevation above the photosphere. The ribbons in the Ca II 8542 Å line wing images are located deeper in the lower solar atmosphere with a lower contrast. A comparison of the height deposition of electron beam energy and the intensity contribution function shows that the H$β$ line wing intensities can be an useful tracer of flare energy deposition in the lower solar atmosphere
△ Less
Submitted 21 May, 2020;
originally announced May 2020.
-
Two-dimensional simulations of coronal rain dynamics. I. Model with vertical magnetic field and an unbounded atmosphere
Authors:
David Martínez-Gómez,
Ramón Oliver,
Elena Khomenko,
Manuel Collados
Abstract:
Aims. We aim to improve the understanding of the physical mechanisms behind the slower than free-fall motion and the two-stage evolution (an initial phase of acceleration followed by an almost constant velocity phase) detected in coronal rain events.
Methods. Using the Mancha3D code, we solve the 2D ideal MHD equations. We represent the solar corona as an isothermal vertically stratified atmosph…
▽ More
Aims. We aim to improve the understanding of the physical mechanisms behind the slower than free-fall motion and the two-stage evolution (an initial phase of acceleration followed by an almost constant velocity phase) detected in coronal rain events.
Methods. Using the Mancha3D code, we solve the 2D ideal MHD equations. We represent the solar corona as an isothermal vertically stratified atmosphere with a uniform vertical magnetic field and the plasma condensation as a density enhancement described by a 2D Gaussian profile. We analyse the temporal evolution of the descending plasma and study its dependence on parameters such as density and magnetic field strength.
Results. We confirm previous findings that the pressure gradient is the main force that opposes the action of gravity and slows down the blob descent and that larger densities require larger pressure gradients to reach the constant speed phase. We find that the shape of a condensation with a horizontal variation of density is distorted as it falls, due to the denser parts of the blob falling faster than the lighter ones. This is explained by the fact that the duration of the initial acceleration phase, and therefore the maximum falling speed attained by the plasma, increases with the ratio of blob to coronal density. We also find that the magnetic field plays a fundamental role in the evolution of the descending condensations. A strong enough magnetic field (greater than 10 G in our simulations) forces each plasma element to follow the path given by a particular field line, which allows to describe the evolution of each vertical slice of the blob in terms of 1D dynamics, without influence of the adjacent slices. In addition, under the typical conditions of the coronal rain events, the magnetic field prevents the development of the Kelvin-Helmholtz instability.
△ Less
Submitted 15 November, 2019;
originally announced November 2019.
-
Fundamental Transverse Vibrations of the Active Region Solar Corona
Authors:
Manuel Luna,
Ramon Oliver,
Patrick Antolin,
Inigo Arregui
Abstract:
Some high-resolution observations have revealed that the active-region solar corona is filled with myriads of thin strands even in apparently uniform regions with no resolved loops. This fine structure can host collective oscillations involving a large portion of the corona due to the coupling of the motions of the neighbouring strands. We study these vibrations and the possible observational effe…
▽ More
Some high-resolution observations have revealed that the active-region solar corona is filled with myriads of thin strands even in apparently uniform regions with no resolved loops. This fine structure can host collective oscillations involving a large portion of the corona due to the coupling of the motions of the neighbouring strands. We study these vibrations and the possible observational effects. Here we theoretically investigate the collective oscillations inherent to the fine structure of the corona. We have called them fundamental vibrations because they cannot exist in a uniform medium. We use the T-matrix technique to find the normal modes of random arrangements of parallel strands. We consider an increasing number of tubes to understand the vibrations of a huge number of tubes of a large portion of the corona. We additionally generate synthetic time-distance Doppler and line broadening diagrams of the vibrations of a coronal region to compare with observations. We have found that the fundamental vibrations are in the form of clusters of tubes where not all the tubes participate in the collective mode. The periods are distributed over a wide band of values. The width of the band increases with the number of strands but rapidly reaches an approximately constant value. The frequency band associated with the fine structure of the corona depends on the minimum separation between strands. We have found that the coupling between the strands is of large extent. The synthetic Dopplergrams and line-broadening maps show signatures of collective vibrations, not present in the case of purely random individual kink vibrations. We conclude that the fundamental vibrations of the corona can contribute to the energy budget of the corona and they may have an observational signature.
△ Less
Submitted 11 July, 2019;
originally announced July 2019.
-
Magneto-Rossby waves in the solar tachocline and the annual variations in solar activity
Authors:
Tamar Gachechiladze,
Teimuraz V. Zaqarashvili,
Eka Gurgenashvili,
Giorgi Ramishvili,
Marc Carbonell,
Ramon Oliver,
Jose Luis Balleste
Abstract:
Annual oscillations have been detected in many indices of solar activity during many cycles. Recent multi spacecraft observations of coronal bright points revealed slow retrograde toroidal phase drift (with the speed of 3 m/s of 1 yr oscillations, which naturally suggested their connection with Rossby-type waves in the interior. We have studied from a theoretical point of view the dynamics of glob…
▽ More
Annual oscillations have been detected in many indices of solar activity during many cycles. Recent multi spacecraft observations of coronal bright points revealed slow retrograde toroidal phase drift (with the speed of 3 m/s of 1 yr oscillations, which naturally suggested their connection with Rossby-type waves in the interior. We have studied from a theoretical point of view the dynamics of global magneto-Kelvin and magneto-Rossby waves in the solar tachocline with toroidal magnetic field. Using spherical coordinates, the dispersion relations of the waves and latitudinal structure of solutions were obtained analytically. We have also obtained the spectrum of unstable magneto-Rossby wave harmonics in the presence of the latitudinal differential rotation. Estimated periods and phase speeds show that the magneto-Rossby waves rather than the Kelvin waves match with the observations of 1 yr oscillations. On the other hand, Morlet wavelet analysis of Greenwich Royal Observatory sunspot areas for the solar cycle 23 has revealed multiple periodicities with periods of 450-460 days, 370-380 days, 310-320 days, 240-270 days, and 150-175 days in hemispheric and full disk data. Comparison of theoretical results with the observations allow us to conclude that the global magneto-Kelvin waves in the upper overshoot tachocline may be responsible for the periodicity of 450-460 days (1.3 yrs), while the remaining periods can be connected with different harmonics of global fast magneto-Rossby waves.
△ Less
Submitted 29 April, 2019;
originally announced April 2019.
-
Mapping the magnetic field of flare coronal loops
Authors:
D. Kuridze,
M. Mathioudakis,
H. Morgan,
R. Oliver,
L. Kleint,
T. V. Zaqarashvili,
A. Reid,
J. Koza,
M. G. Löfdahl,
T. Hillberg,
V. Kukhianidze,
A. Hanslmeier
Abstract:
Here we report on the unique observation of flaring coronal loops at the solar limb using high resolution imaging spectropolarimetry from the Swedish 1-meter Solar Telescope. The vantage position, orientation and nature of the chromospheric material that filled the flare loops allowed us to determine their magnetic field with unprecedented accuracy using the weak-field approximation method. Our an…
▽ More
Here we report on the unique observation of flaring coronal loops at the solar limb using high resolution imaging spectropolarimetry from the Swedish 1-meter Solar Telescope. The vantage position, orientation and nature of the chromospheric material that filled the flare loops allowed us to determine their magnetic field with unprecedented accuracy using the weak-field approximation method. Our analysis reveals coronal magnetic field strengths as high as 350 Gauss at heights up to 25 Mm above the solar limb. These measurements are substantially higher than a number of previous estimates and may have considerable implications for our current understanding of the extended solar atmosphere.
△ Less
Submitted 20 February, 2019;
originally announced February 2019.
-
Normal modes of transverse coronal loop oscillations from numerical simulations: I. Method and test case
Authors:
S. Rial,
I. Arregui,
R. Oliver,
J. Terradas
Abstract:
The purpose of this work is to develop a procedure to obtain the normal modes of a coronal loop from time-dependent numerical simulations with the aim of better understanding observed transverse loop oscillations. To achieve this goal, in this paper we present a new method and test its performance with a problem for which the normal modes can be computed analytically. In a follow-up paper, the app…
▽ More
The purpose of this work is to develop a procedure to obtain the normal modes of a coronal loop from time-dependent numerical simulations with the aim of better understanding observed transverse loop oscillations. To achieve this goal, in this paper we present a new method and test its performance with a problem for which the normal modes can be computed analytically. In a follow-up paper, the application to the simulations of \citet{rial2013} is tackled. The method proceeds iteratively and at each step consists of (i) a time-dependent numerical simulation followed by (ii) the Complex Empirical Orthogonal Function (CEOF) analysis of the simulation results. The CEOF analysis provides an approximation to the normal mode eigenfunctions that can be used to set up the initial conditions for the numerical simulation of the following iteration, in which an improved normal mode approximation is obtained. The iterative process is stopped once the global difference between successive approximate eigenfunctions is below a prescribed threshold. The equilibrium used in this paper contains material discontinuities that result in one eigenfunction with a jump across these discontinuities and two eigenfunctions whose normal derivatives are discontinuous there. After 6 iterations, the approximation to the frequency and eigenfunctions are accurate to $\lesssim 0.7$\% except for the eigenfunction with discontinuities, which displays a much larger error at these positions.
△ Less
Submitted 13 June, 2019; v1 submitted 1 February, 2019;
originally announced February 2019.
-
Energy transport and heating by torsional Alfvén waves propagating from the photosphere to the corona in the quiet Sun
Authors:
Roberto Soler,
Jaume Terradas,
Ramon Oliver,
Jose Luis Ballester
Abstract:
In the solar atmosphere, Alfvén waves are believed to play an important role in the transfer of energy from the photosphere to the corona and solar wind, and in the heating of the chromosphere. We perform numerical computations to investigate energy transport and dissipation associated with torsional Alfvén waves propagating in magnetic flux tubes that expand from the photosphere to the corona in…
▽ More
In the solar atmosphere, Alfvén waves are believed to play an important role in the transfer of energy from the photosphere to the corona and solar wind, and in the heating of the chromosphere. We perform numerical computations to investigate energy transport and dissipation associated with torsional Alfvén waves propagating in magnetic flux tubes that expand from the photosphere to the corona in quiet-Sun conditions. We place a broadband driver at the photosphere that injects a wave energy flux of $10^7$ erg cm$^{-2}$ s$^{-1}$ and consider Ohm's magnetic diffusion and ion-neutral collisions as dissipation mechanisms. We find that only a small fraction of the driven flux, $\sim 10^5$ erg cm$^{-2}$ s$^{-1}$, is able to reach coronal heights, but it may be sufficient to partly compensate the total coronal energy loss. The frequency of maximal transmittance is $\sim 5$ mHz for a photospheric field strength of 1 kG and is shifted to smaller/larger frequencies for weaker/stronger fields. Lower frequencies are reflected at the transition region, while higher frequencies are dissipated producing enough heat to balance chromospheric radiative losses. Heating in the low and middle chromosphere is due to Ohmic dissipation, while ion-neutral friction dominates in the high chromosphere. Ohmic diffusion is enhanced by phase mixing because of the expansion of the magnetic field. This effect has the important consequence of increasing the chromospheric dissipation and, therefore, reducing the energy flux that reaches the corona. We provide empirical fits of the transmission coefficient that could be used as input for coronal models.
△ Less
Submitted 4 December, 2018;
originally announced December 2018.
-
The Coronal Monsoon: Thermal Nonequilibrium Revealed by Periodic Coronal Rain
Authors:
F. Auchère,
C. froment,
E. Soubrié,
P. Antolin,
R. Oliver,
G. Pelouze
Abstract:
We report on the discovery of periodic coronal rain in an off-limb sequence of {\it Solar Dynamics Observatory}/Atmospheric Imaging Assembly images. The showers are co-spatial and in phase with periodic (6.6~hr) intensity pulsations of coronal loops of the sort described by Auchere et al. (2014) and Froment et al. (2015, 2017). These new observations make possible a unified description of both phe…
▽ More
We report on the discovery of periodic coronal rain in an off-limb sequence of {\it Solar Dynamics Observatory}/Atmospheric Imaging Assembly images. The showers are co-spatial and in phase with periodic (6.6~hr) intensity pulsations of coronal loops of the sort described by Auchere et al. (2014) and Froment et al. (2015, 2017). These new observations make possible a unified description of both phenomena. Coronal rain and periodic intensity pulsations of loops are two manifestations of the same physical process: evaporation / condensation cycles resulting from a state of thermal nonequilibrium (TNE). The fluctuations around coronal temperatures produce the intensity pulsations of loops, and rain falls along their legs if thermal runaway cools the periodic condensations down and below transition-region (TR) temperatures. This scenario is in line with the predictions of numerical models of quasi-steadily and footpoint heated loops. The presence of coronal rain -- albeit non-periodic -- in several other structures within the studied field of view implies that this type of heating is at play on a large scale.
△ Less
Submitted 6 February, 2018;
originally announced February 2018.
-
North-south asymmetry in Rieger-type periodicity during solar cycles 19-23
Authors:
Eka Gurgenashvili,
Teimuraz V. Zaqarashvili,
Vasil Kukhianidze,
Ramon Oliver,
Jose Luis Ballester,
Mausumi Dikpati,
Scott W. McIntosh
Abstract:
Rieger-type periodicity has been detected in different activity indices over many solar cycles. It was recently shown that the periodicity correlates with solar activity having a shorter period during stronger cycles. Solar activity level is generally asymmetric between northern and southern hemispheres, which could suggest the presence of a similar behavior in the Rieger-type periodicity. We anal…
▽ More
Rieger-type periodicity has been detected in different activity indices over many solar cycles. It was recently shown that the periodicity correlates with solar activity having a shorter period during stronger cycles. Solar activity level is generally asymmetric between northern and southern hemispheres, which could suggest the presence of a similar behavior in the Rieger-type periodicity. We analyse the sunspot area/number and the total magnetic flux data for northern and southern hemispheres during solar cycles 19-23 which had remarkable north-south asymmetry. Using wavelet analysis of sunspot area and number during the north-dominated cycles (19-20) we obtained the periodicity of 160-165 days in the stronger northern hemisphere and 180-190 days in the weaker southern hemisphere. On the other hand, south-dominated cycles (21-23) display the periodicity of 155-160 days in the stronger southern hemisphere and 175-188 days in the weaker northern hemisphere. Therefore, the Rieger-type periodicity has the north-south asymmetry in sunspot area/number data during solar cycles with strong hemispheric asymmetry. We suggest that the periodicity is caused by magnetic Rossby waves in the internal dynamo layer. Using the dispersion relation of magnetic Rossby waves and observed Rieger periodicity we estimated the magnetic field strength in the layer as 45-49 kG in more active hemispheres (north during the cycles 19-20 and south during the cycles 21-23) and 33-40 kG in weaker hemispheres. The estimated difference in the hemispheric field strength is around 10 kG, which provides a challenge for dynamo models. Total magnetic flux data during the cycle 20-23 reveals no clear north-south asymmetry which needs to be explained in the future.
△ Less
Submitted 26 July, 2017;
originally announced July 2017.
-
IRIS Burst Spectra Co-Spatial To A Quiet-Sun Ellerman-Like Brightening
Authors:
C. J. Nelson,
N. Freij,
A. Reid,
R. Oliver,
M. Mathioudakis,
R. Erdélyi
Abstract:
Ellerman bombs (EBs) have been widely studied over the past two decades; however, only recently have counterparts of these events been observed in the quiet-Sun. The aim of this article is to further understand small-scale quiet-Sun Ellerman-like brightenings (QSEBs) through research into their spectral signatures, including investigating whether the hot signatures associated with some EBs are als…
▽ More
Ellerman bombs (EBs) have been widely studied over the past two decades; however, only recently have counterparts of these events been observed in the quiet-Sun. The aim of this article is to further understand small-scale quiet-Sun Ellerman-like brightenings (QSEBs) through research into their spectral signatures, including investigating whether the hot signatures associated with some EBs are also visible co-spatial to any QSEBs. We combine H$α$ and Ca II $8542$ Å line scans at the solar limb with spectral and imaging data sampled by the Interface Region Imaging Spectrograph (IRIS). Twenty one QSEBs were identified with average lifetimes, lengths, and widths measured to be around $120$ s, $0.63$", and $0.35$", respectively. Three of these QSEBs displayed clear repetitive flaring through their lifetimes, comparable to the behaviour of EBs in Active Regions (ARs). Two QSEBs in this sample occurred co-spatial with increased emission in SDO/AIA $1600$ Å and IRIS slit-jaw imager $1400$ Å data, however, these intensity increases were smaller than reported co-spatial to EBs. One QSEB was also sampled by the IRIS slit during its lifetime, displaying increases in intensity in the Si IV $1393$ Å and Si IV $1403$ Å cores as well as the C II and Mg II line wings, analogous to IRIS bursts (IBs). Using RADYN simulations, we are unable to reproduce the observed QSEB H$α$ and Ca II $8542$ Å line profiles leaving the question of the temperature stratification of QSEBs open. Our results imply that some QSEBs could be heated to Transition Region temperatures, suggesting that IB profiles should be observed throughout the quiet-Sun.
△ Less
Submitted 17 July, 2017;
originally announced July 2017.
-
Dissipative instability in a partially ionised prominence slab
Authors:
I. Ballai,
B. Pinter,
R. Oliver,
M. Alexandrou
Abstract:
We investigate the nature of dissipative instability appearing in a prominence planar thread filled with partially ionised plasma in the incompressible limit. The importance of partial ionisation is investigated in terms of the ionisation factor and wavelength of waves propagating in the slab. To highlight the role of partial ionisation, we have constructed models describing various situations we…
▽ More
We investigate the nature of dissipative instability appearing in a prominence planar thread filled with partially ionised plasma in the incompressible limit. The importance of partial ionisation is investigated in terms of the ionisation factor and wavelength of waves propagating in the slab. To highlight the role of partial ionisation, we have constructed models describing various situations we can meet in solar prominence fine structure. Matching the solutions for the transversal component of the velocity and total pressure at the interfaces between the prominence slab and surrounding plasmas, we derived a dispersion relation whose imaginary part describes the evolution of the instability. Results are obtained in the limit of weak dissipation. We have investigated the appearance of instabilities in prominence dark plumes using single and two-fluid approximations. We show that dissipative instabilities appear for flow speeds that are less than the Kelvin-Helmholtz instability threshold. The onset of instability is determined by the equilibrium flow strength, the ionisation factor of the plasma, the wavelength of waves and the ion-neutral collisional rate. For a given wavelength and for ionisation degrees closer to a neutral gas, the propagating waves become unstable for a narrow band of flow speeds, meaning that neutrals have a stabilising effect. Our results show that the partially ionised plasma describing prominence dark plumes becomes unstable only in a two-fluid (charged particles-neutrals) model, that is for periods that are smaller than the ion-neutral collision time. The present study improves our understanding of stability of solar prominences and the role of partial ionisation in destabilising the plasma. We show the necessity of two-fluid approximation when discussing the nature of instabilities: waves in a single fluid approximation show a great deal of stability.
△ Less
Submitted 21 March, 2017;
originally announced March 2017.
-
The effects of magnetic-field geometry on longitudinal oscillations of solar prominences: Cross-sectional area variation for thin tubes
Authors:
M. Luna,
A. J. Diaz,
R. Oliver,
J. Terradas,
J. Karpen
Abstract:
Solar prominences are subject to both field-aligned (longitudinal) and transverse oscillatory motions, as evidenced by an increasing number of observations. Large-amplitude longitudinal motions provide valuable information on the geometry of the filament-channel magnetic structure that supports the cool prominence plasma against gravity. Our pendulum model, in which the restoring force is the grav…
▽ More
Solar prominences are subject to both field-aligned (longitudinal) and transverse oscillatory motions, as evidenced by an increasing number of observations. Large-amplitude longitudinal motions provide valuable information on the geometry of the filament-channel magnetic structure that supports the cool prominence plasma against gravity. Our pendulum model, in which the restoring force is the gravity projected along the dipped field lines of the magnetic structure, best explains these oscillations. However, several factors can influence the longitudinal oscillations, potentially invalidating the pendulum model. The aim of this work is to study the influence of large-scale variations in the magnetic field strength along the field lines, i.e., variations of the cross-sectional area along the flux tubes supporting prominence threads. We studied the normal modes of several flux tube configurations, using linear perturbation analysis, to assess the influence of different geometrical parameters on the oscillation properties. We found that the influence of the symmetric and asymmetric expansion factors on longitudinal oscillations is small.}{We conclude that the longitudinal oscillations are not significantly influenced by variations of the cross-section of the flux tubes, validating the pendulum model in this context.
△ Less
Submitted 11 July, 2016;
originally announced July 2016.
-
The role of Alfvén wave heating in solar prominences
Authors:
Roberto Soler,
Jaume Terradas,
Ramon Oliver,
Jose Luis Ballester
Abstract:
Observations have shown that magnetohydrodynamic waves over a large frequency range are ubiquitous in solar prominences. The waves are probably driven by photospheric motions and may transport energy up to prominences suspended in the corona. Dissipation of wave energy can lead to heating of the cool prominence plasma, so contributing to the local energy balance within the prominence. Here we disc…
▽ More
Observations have shown that magnetohydrodynamic waves over a large frequency range are ubiquitous in solar prominences. The waves are probably driven by photospheric motions and may transport energy up to prominences suspended in the corona. Dissipation of wave energy can lead to heating of the cool prominence plasma, so contributing to the local energy balance within the prominence. Here we discuss the role of Alfvén wave dissipation as a heating mechanism for the prominence plasma. We consider a slab-like quiescent prominence model with a transverse magnetic field embedded in the solar corona. The prominence medium is modelled as a partially ionized plasma composed of a charged ion-electron single fluid and two separate neutral fluids corresponding to neutral hydrogen and neutral helium. Friction between the three fluids acts as a dissipative mechanism for the waves. The heating caused by externally-driven Alfvén waves incident on the prominence slab is analytically explored. We find that the dense prominence slab acts as a resonant cavity for the waves. The fraction of incident wave energy that is channelled into the slab strongly depends upon the wave period, $P$. Using typical prominence conditions, we obtain that wave energy trapping and associated heating are negligible when $P \gtrsim 100$ s, so that it is unlikely that those waves have a relevant influence on prominence energetics. When $1$ s $\lesssim P \lesssim 100$ s the energy absorption into the slab shows several sharp and narrow peaks, that can reach up to 100%, when the incident wave frequency matches a cavity resonance of the slab. Wave heating is enhanced at those resonant frequencies. Conversely, when $P \lesssim 1$ s cavity resonances are absent, but the waves are heavily damped by the strong dissipation. We estimate that wave heating may compensate for about 10% of radiative losses of the prominence plasma.
△ Less
Submitted 23 May, 2016;
originally announced May 2016.
-
Rieger-type periodicity during solar cycles 14-24: estimation of dynamo magnetic field strength in the solar interior
Authors:
Eka Gurgenashvili,
Teimuraz Zaqarashvili,
Vasil Kukhianidze,
Ramon Oliver,
Jose Luis Ballester,
Giorgi Ramishvili,
Bidzina Shergelashvili,
Arnold Hanslmeier,
Stefaan Poedts
Abstract:
Solar activity undergoes a variation over time scales of several months known as Rieger-type periodicity, which usually occurs near maxima of sunspot cycles. An early analysis showed that the periodicity appears only in some cycles, and is absent in other cycles. But the appearance/absence during different cycles has not been explained. We performed a wavelet analysis of sunspot data from the Gree…
▽ More
Solar activity undergoes a variation over time scales of several months known as Rieger-type periodicity, which usually occurs near maxima of sunspot cycles. An early analysis showed that the periodicity appears only in some cycles, and is absent in other cycles. But the appearance/absence during different cycles has not been explained. We performed a wavelet analysis of sunspot data from the Greenwich Royal Observatory and the Royal Observatory of Belgium during cycles 14-24. We found that the Rieger-type periods occur in all cycles, but they are cycle-dependent: shorter periods occur during stronger cycles. Our analysis revealed a periodicity of 185-195 days during the weak cycles 14-15 and 24, and a periodicity of 155-165 days during the stronger cycles 16-23. We derived the dispersion relation of the spherical harmonics of the magnetic Rossby waves in the presence of differential rotation and a toroidal magnetic field in the dynamo layer near the base of the convection zone. This showed that the harmonic of fast Rossby waves with m=1 and n=4, where m (n) indicate the toroidal (poloidal) wavenumbers, respectively, perfectly fit with the observed periodicity. The variation of the toroidal field strength from weaker to stronger cycles may lead to the different periods found in those cycles, which explains the observed enigmatic feature of the Rieger-type periodicity. Finally, we used the observed periodicity to estimate the dynamo field strength during cycles 14-24. Our estimations suggest a field strength of 40 kG for the stronger cycles, and 20 kG for the weaker cycles.
△ Less
Submitted 13 May, 2016;
originally announced May 2016.
-
ALMA Observations of the Sun in Cycle 4 and Beyond
Authors:
S. Wedemeyer,
B. Fleck,
M. Battaglia,
N. Labrosse,
G. Fleishman,
H. Hudson,
P. Antolin,
C. Alissandrakis,
T. Ayres,
J. Ballester,
T. Bastian,
J. Black,
A. Benz,
R. Brajsa,
M. Carlsson,
J. Costa,
B. DePontieu,
G. Doyle,
G. Gimenez de Castro,
S. Gunár,
G. Harper,
S. Jafarzadeh,
M. Loukitcheva,
V. Nakariakov,
R. Oliver
, et al. (12 additional authors not shown)
Abstract:
This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interestin…
▽ More
This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interesting observations could be made already with the still limited solar observing modes foreseen for Cycle 4 and that ALMA has the potential to make important contributions to answering long-standing scientific questions in solar physics. With the proposal deadline for ALMA Cycle 4 in April 2016 and the Commissioning and Science Verification campaign in December 2015 in sight, several of the SSALMON Expert Teams composed strategic documents in which they outlined potential solar observations that could be feasible given the anticipated technical capabilities in Cycle 4. These documents have been combined and supplemented with an analysis, resulting in recommendations for solar observing with ALMA in Cycle 4. In addition, the detailed science cases also demonstrate the scientific priorities of the solar physics community and which capabilities are wanted for the next observing cycles. The work on this White Paper effort was coordinated in close cooperation with the two international solar ALMA development studies led by T. Bastian (NRAO, USA) and R. Brajsa, (ESO). This document will be further updated until the beginning of Cycle 4 in October 2016. In particular, we plan to adjust the technical capabilities of the solar observing modes once finally decided and to further demonstrate the feasibility and scientific potential of the included science cases by means of numerical simulations of the solar atmosphere and corresponding simulated ALMA observations.
△ Less
Submitted 11 April, 2016; v1 submitted 4 January, 2016;
originally announced January 2016.
-
Solar prominences embedded in flux ropes: morphological features and dynamics from 3D MHD simulations
Authors:
J. Terradas,
R. Soler,
M. Luna,
R. Oliver,
J. L. Ballester,
A. N. Wright
Abstract:
The temporal evolution of a solar prominence inserted in a three-dimensional magnetic flux rope is investigated numerically. Using the model of Titov Demoulin (1999) under the regime of weak twist, the cold and dense prominence counteracts gravity by modifying the initially force-free magnetic configuration. In some cases a quasi-stationary situation is achieved after the relaxation phase, charact…
▽ More
The temporal evolution of a solar prominence inserted in a three-dimensional magnetic flux rope is investigated numerically. Using the model of Titov Demoulin (1999) under the regime of weak twist, the cold and dense prominence counteracts gravity by modifying the initially force-free magnetic configuration. In some cases a quasi-stationary situation is achieved after the relaxation phase, characterized by the excitation of standing vertical oscillations. These oscillations show a strong attenuation with time produced by the mechanism of continuum damping due to the inhomogeneous transition between the prominence and solar corona. The characteristic period of the vertical oscillations does not depend strongly on the twist of the flux rope. Nonlinearity is the responsible for triggering the Kelvin-Helmholtz instability associated to the vertical oscillations and that eventually produces horizontal structures. Contrary to other configurations in which the longitudinal axis of the prominence is permeated by a perpendicular magnetic field, like in unsheared arcades, the orientation of the prominence along the flux rope axis prevents the development of Rayleigh-Taylor instabilities and therefore the appearance of vertical structuring along this axis.
△ Less
Submitted 22 December, 2015;
originally announced December 2015.
-
Dynamics of coronal rain and descending plasma blobs in solar prominences: II. partially ionized case
Authors:
R. Oliver,
R. Soler,
J. Terradas,
T. V. Zaqarashvili
Abstract:
Coronal rain clumps and prominence knots are dense condensations with chromospheric to transition region temperatures that fall down in the much hotter corona. Their typical speeds are in the range 30--150~km~s$^{-1}$ and of the order of 10--30~km~s$^{-1}$, respectively, i.e., they are considerably smaller than free fall velocities. These cold blobs contain a mixture of ionized and neutral materia…
▽ More
Coronal rain clumps and prominence knots are dense condensations with chromospheric to transition region temperatures that fall down in the much hotter corona. Their typical speeds are in the range 30--150~km~s$^{-1}$ and of the order of 10--30~km~s$^{-1}$, respectively, i.e., they are considerably smaller than free fall velocities. These cold blobs contain a mixture of ionized and neutral material that must be dynamically coupled in order to fall together, as observed. We investigate this coupling by means of hydrodynamic simulations in which the coupling arises from the friction between ions and neutrals. The numerical simulations presented here are an extension of those of \citet{oliver2014} to the partially ionized case. We find that, although the relative drift speed between the two species is smaller than 1~m~s$^{-1}$ at the blob center, it is sufficient to produce the forces required to strongly couple charged particles and neutrals. The ionization degree has no discernible effect on the main results of our previous work for a fully ionized plasma: the condensation has an initial acceleration phase followed by a period with roughly constant velocity and, in addition, the maximum descending speed is clearly correlated with the ratio of initial blob to environment density.
△ Less
Submitted 28 October, 2015;
originally announced October 2015.
-
Long-term variation in the Sun's activity caused by magnetic Rossby waves in the tachocline
Authors:
T. V. Zaqarashvili,
R. Oliver,
A. Hanslmeier,
M. Carbonell,
J. L. Ballester,
T. Gachechiladze,
I. G. Usoskin
Abstract:
Long-term records of sunspot number and concentrations of cosmogenic radionuclides (10Be and 14C) on the Earth reveal the variation of the Sun's magnetic activity over hundreds and thousands of years. We identify several clear periods in sunspot, 10Be, and 14C data as 1000, 500, 350, 200 and 100 years. We found that the periods of the first five spherical harmonics of the slow magnetic Rossby mode…
▽ More
Long-term records of sunspot number and concentrations of cosmogenic radionuclides (10Be and 14C) on the Earth reveal the variation of the Sun's magnetic activity over hundreds and thousands of years. We identify several clear periods in sunspot, 10Be, and 14C data as 1000, 500, 350, 200 and 100 years. We found that the periods of the first five spherical harmonics of the slow magnetic Rossby mode in the presence of a steady toroidal magnetic field of 1200-1300 G in the lower tachocline are in perfect agreement with the time scales of observed variations. The steady toroidal magnetic field can be generated in the lower tachocline either due to the steady dynamo magnetic field for low magnetic diffusivity or due to the action of the latitudinal differential rotation on the weak poloidal primordial magnetic field, which penetrates from the radiative interior. The slow magnetic Rossby waves lead to variations of the steady toroidal magnetic field in the lower tachocline, which modulate the dynamo magnetic field and consequently the solar cycle strength. This result constitutes a key point for long-term prediction of the cycle strength. According to our model, the next deep minimum in solar activity is expected during the first half of this century.
△ Less
Submitted 11 May, 2015;
originally announced May 2015.
-
Formation and evolution of coronal rain observed by SDO/AIA on February 22, 2012
Authors:
Z. Vashalomidze,
V. Kukhianidze,
T. V. Zaqarashvili,
R. Oliver,
B. Shergelashvili,
G. Ramishvili,
S. Poedts,
P. De Causmaecker
Abstract:
The formation and dynamics of coronal rain are currently not fully understood. Coronal rain is the fall of cool and dense blobs formed by thermal instability in the solar corona towards the solar surface with acceleration smaller than gravitational free fall. We aim to study the observational evidence of the formation of coronal rain and to trace the detailed dynamics of individual blobs. We used…
▽ More
The formation and dynamics of coronal rain are currently not fully understood. Coronal rain is the fall of cool and dense blobs formed by thermal instability in the solar corona towards the solar surface with acceleration smaller than gravitational free fall. We aim to study the observational evidence of the formation of coronal rain and to trace the detailed dynamics of individual blobs. We used time series of the 171 Å\, and 304 Å\, spectral lines obtained by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) above active region AR 11420 on February 22, 2012. Observations show that a coronal loop disappeared in the 171 Å channel and appeared in the 304 Å line$\text{}\text{}$ more than one hour later, which indicates a rapid cooling of the coronal loop from 1 MK to 0.05 MK. An energy estimation shows that the radiation is higher than the heat input, which indicates so-called catastrophic cooling. The cooling was accompanied by the formation of coronal rain in the form of falling cold plasma. We studied two different sequences of falling blobs. The first sequence includes three different blobs. The mean velocities of the blobs were estimated to be 50 km s$^{-1}$, 60 km s$^{-1}$ and 40 km s$^{-1}$. A polynomial fit shows the different values of the acceleration for different blobs, which are lower than free-fall in the solar corona. The first and second blob move along the same path, but with and without acceleration, respectively. We performed simple numerical simulations for two consecutive blobs, which show that the second blob moves in a medium that is modified by the passage of the first blob. Therefore, the second blob has a relatively high speed and no acceleration, as is shown by observations. The second sequence includes two different blobs with mean velocities of 100 km s$^{-1}$ and 90 km s$^{-1}$, respectively.
△ Less
Submitted 14 April, 2015;
originally announced April 2015.
-
Dissipative instability in partially ionised prominence plasmas
Authors:
Istvan Ballai,
Ramon Oliver,
Marios Alexandrou
Abstract:
We investigate the nature of dissipative instability at the boundary (seen here as tangential discontinuity) between the viscous corona and the partially ionised prominence plasma in the incompressible limit. The importance of the partial ionisation is investigated in terms of the ionisation fraction. Matching the solutions for the transversal component of the velocity and total pressure at the in…
▽ More
We investigate the nature of dissipative instability at the boundary (seen here as tangential discontinuity) between the viscous corona and the partially ionised prominence plasma in the incompressible limit. The importance of the partial ionisation is investigated in terms of the ionisation fraction. Matching the solutions for the transversal component of the velocity and total pressure at the interface between the prominence and coronal plasmas, we derive a dispersion relation whose imaginary part describes the evolution of the instability. Results are obtained in the limit of weak dissipation. Using simple analytical methods, we show that dissipative instabilities appear for flow speeds that are lower than the Kelvin-Helmholtz instability threshold. While viscosity tends to destabilise the plasma, the effect of partial ionisation (through the Cowling resistivity) will act towards stabilising the interface. For ionisation degrees closer to a neutral gas the interface will be unstable for larger values of equilibrium flow. The same principle is assumed when studying the appearance of instability at the interface between prominences and dark plumes. The unstable mode appearing in this case has a very small growth rate and dissipative instability cannot explain the appearance of flows in plumes. The present study improves our understanding of the complexity of dynamical processes at the interface of solar prominences and solar corona, and the role partial ionisation can have on the stability of the plasma. Our results clearly show that the problem of partial ionisation introduces new aspects of plasma stability with consequences on the evolution of solar prominences.
△ Less
Submitted 26 March, 2015;
originally announced March 2015.
-
On the support of neutrals against gravity in solar prominences
Authors:
J. Terradas,
R. Soler,
R. Oliver,
J. L. Ballester
Abstract:
Cool and dense prominences found in the solar atmosphere are known to be partially ionized because of their relative low temperature. In this Letter, we address the long-standing problem of how the neutral component of the plasma in prominences is supported against gravity. Using the multiple fluid approach we solve the time-dependent equations in two dimensions considering the frictional coupling…
▽ More
Cool and dense prominences found in the solar atmosphere are known to be partially ionized because of their relative low temperature. In this Letter, we address the long-standing problem of how the neutral component of the plasma in prominences is supported against gravity. Using the multiple fluid approach we solve the time-dependent equations in two dimensions considering the frictional coupling between the neutral and ionized components of the magnetized plasma representative of a solar prominence embedded in a hot coronal environment. We demonstrate that given an initial density enhancement in the two fluids, representing the body of the prominence, the system is able to relax in the vicinity of magnetic dips to a stationary state in which both neutrals and ionized species are dynamically suspended above the photosphere. Two different coupling processes are considered in this study, collisions between ions and neutrals and charge exchange interactions. We find that for realistic conditions ions are essentially static while neutrals have a very small downflow velocity. The coupling between ions and neutrals is so strong at the prominence body that the behavior is similar to that of a single fluid with an effective density equal to the sum of the ion and neutral species. We also find that the charge exchange mechanism is about three times more efficient sustaining neutrals than elastic scattering of ions with neutrals.
△ Less
Submitted 18 March, 2015;
originally announced March 2015.
-
Propagation and dispersion of sausage wave trains in magnetic flux tubes
Authors:
R. Oliver,
M. S. Ruderman,
J. Terradas
Abstract:
A localized perturbation of a magnetic flux tube produces a pair of wave trains that propagate in opposite directions along the tube. These wave packets disperse as they propagate, where the extent of dispersion depends on the physical properties of the magnetic structure, on the length of the initial excitation, and on its nature (e.g., transverse or axisymmetric). In Oliver et al. (2014) we cons…
▽ More
A localized perturbation of a magnetic flux tube produces a pair of wave trains that propagate in opposite directions along the tube. These wave packets disperse as they propagate, where the extent of dispersion depends on the physical properties of the magnetic structure, on the length of the initial excitation, and on its nature (e.g., transverse or axisymmetric). In Oliver et al. (2014) we considered a transverse initial perturbation, whereas the temporal evolution of an axisymmetric one is examined here. In both papers we use a method based on Fourier integrals to solve the initial value problem. Previous studies on wave propagation in magnetic wave guides have emphasized that the wave train dispersion is influenced by the particular dependence of the group velocity on the longitudinal wavenumber. Here we also find that long initial perturbations result in low amplitude wave packets and that large values of the magnetic tube to environment density ratio yield longer wave trains. To test the detectability of propagating transverse or axisymmetric wave packets in magnetic tubes of the solar atmosphere (e.g., coronal loops, spicules, or prominence threads) a forward modelling of the perturbations must be carried out. This is left for a future work.
△ Less
Submitted 4 February, 2015;
originally announced February 2015.
-
Morphology and dynamics of solar prominences from 3D MHD simulations
Authors:
J. Terradas,
R. Soler,
M. Luna,
R. Oliver,
J. L. Ballester
Abstract:
In this paper we present a numerical study of the time evolution of solar prominences embedded in sheared magnetic arcades. The prominence is represented by a density enhancement in a background stratified atmosphere and is connected to the photosphere through the magnetic field. By solving the ideal magnetohydrodynamic (MHD) equations in three dimensions we study the dynamics for a range of param…
▽ More
In this paper we present a numerical study of the time evolution of solar prominences embedded in sheared magnetic arcades. The prominence is represented by a density enhancement in a background stratified atmosphere and is connected to the photosphere through the magnetic field. By solving the ideal magnetohydrodynamic (MHD) equations in three dimensions we study the dynamics for a range of parameters representative of real prominences. Depending on the parameters considered, we find prominences that are suspended above the photosphere, i.e., detached prominences, but also configurations resembling curtain or hedgerow prominences whose material continuously connects to the photosphere. The plasma$-β$ is an important parameter that determines the shape of the structure. In many cases magnetic Rayleigh-Taylor (MRT) instabilities and oscillatory phenomena develop. Fingers and plumes are generated, affecting the whole prominence body and producing vertical structures in an essentially horizontal magnetic field. However, magnetic shear is able to reduce or even to suppress this instability.
△ Less
Submitted 23 December, 2014;
originally announced December 2014.
-
Propagation and dispersion of transverse wave trains in magnetic flux tubes
Authors:
R. Oliver,
M. S. Ruderman,
J. Terradas
Abstract:
The dispersion of small amplitude, impulsively excited wave trains propagating along a magnetic flux tube is investigated. The initial disturbance is a localized transverse displacement of the tube that excites a fast kink wave packet. The spatial and temporal evolution of the perturbed variables (density, plasma displacement, velocity, ...) is given by an analytical expression containing an integ…
▽ More
The dispersion of small amplitude, impulsively excited wave trains propagating along a magnetic flux tube is investigated. The initial disturbance is a localized transverse displacement of the tube that excites a fast kink wave packet. The spatial and temporal evolution of the perturbed variables (density, plasma displacement, velocity, ...) is given by an analytical expression containing an integral that is computed numerically. We find that the dispersion of fast kink wave trains is more important for shorter initial disturbances (i.e. more concentrated in the longitudinal direction) and for larger density ratios (i.e. for larger contrasts of the tube density with respect to the environment density). This type of excitation generates a wave train whose signature at a fixed position along a coronal loop is a short event (duration ~ 20 s) in which the velocity and density oscillate very rapidly with typical periods of the order of a few seconds. The oscillatory period is not constant but gradually declines during the course of this event. Peak values of the velocity are of the order of 10 km/s and are accompanied by maximum density variations of the order of 10-15% the unperturbed loop density.
△ Less
Submitted 18 February, 2014;
originally announced February 2014.
-
Synthetic hydrogen spectra of prominence oscillations
Authors:
Petr Heinzel,
Maciej Zapiór,
Ramon Oliver,
Jose Luis Ballester
Abstract:
Context. Prominence oscillations have been mostly detected using Doppler velocity, although there are also claimed detections by means of periodic variations in half-width or line intensity. However, scarce observational evidence exists about simultaneous detection of oscillations in several spectral indicators. Aims. Our main aim here is to explore the relationship between spectral indicators, su…
▽ More
Context. Prominence oscillations have been mostly detected using Doppler velocity, although there are also claimed detections by means of periodic variations in half-width or line intensity. However, scarce observational evidence exists about simultaneous detection of oscillations in several spectral indicators. Aims. Our main aim here is to explore the relationship between spectral indicators, such as Doppler shift, line intensity, and line half-width, and the linear perturbations excited in a simple prominence model. Methods. Our equilibrium background model consists of a bounded, homogeneous slab, which is permeated by a transverse magnetic field, having prominence-like physical properties. Assuming linear perturbations, the dispersion relation for fast and slow modes has been derived, as well as the perturbations for the different physical quantities. These perturbations have been used as the input variables in a one-dimensional radiative transfer code, which calculates the full spectral profile of the hydrogen H-alpha and H-beta lines. Results. We have found that different oscillatory modes produce spectral indicator variations in different magnitudes. Detectable variations in the Doppler velocity were found for the fundamental slow mode only. Substantial variations in the H-beta line intensity were found for specific modes. Other modes lead to lower and even undetectable parameter variations. Conclusions. To perform prominence seismology, analysis of the H-alpha and H-beta spectral line parameters could be a good tool to detect and identify oscillatory modes.
△ Less
Submitted 9 January, 2014;
originally announced January 2014.
-
The behavior of transverse waves in nonuniform solar flux tubes. II. Implications for coronal loop seismology
Authors:
Roberto Soler,
Marcel Goossens,
Jaume Terradas,
Ramon Oliver
Abstract:
Seismology of coronal loops using observations of damped transverse oscillations in combination with results from theoretical models is a tool to indirectly infer physical parameters in the solar atmospheric plasma. Existing seismology schemes based on approximations to the period and damping time of kink oscillations are often used beyond their theoretical range of applicability. These approximat…
▽ More
Seismology of coronal loops using observations of damped transverse oscillations in combination with results from theoretical models is a tool to indirectly infer physical parameters in the solar atmospheric plasma. Existing seismology schemes based on approximations to the period and damping time of kink oscillations are often used beyond their theoretical range of applicability. These approximations assume that the variation of density across the loop is confined to a nonuniform layer much thinner than the radius of the loop, but the results of the inversion problem often do not satisfy this preliminary hypothesis. Here, we determine the accuracy of the analytic approximations to the period and damping time, and its impact on seismology estimates, when largely nonuniform loops are considered. We find that the accuracy of the approximations when used beyond their range of applicability is strongly affected by the form of the density profile across the loop, that is observationally unknown and so must be arbitrarily imposed as part of the theoretical model. The error associated with the analytic approximations can be larger than 50% even for relatively thin nonuniform layers. This error directly affects the accuracy of approximate seismology estimates compared to actual numerical inversions. In addition, assuming different density profiles can produce noncoincident intervals of the seismic variables in inversions of the same event. The ignorance about the true shape of density variation across the loop is an important source of error that may dispute the reliability of parameters seismically inferred assuming an ad hoc density profile.
△ Less
Submitted 18 December, 2013;
originally announced December 2013.
-
Dynamics of coronal rain and descending plasma blobs in solar prominences: I. Fully ionised case
Authors:
R. Oliver,
R. Soler,
J. Terradas,
T. V. Zaqarashvili,
M. L. Khodachenko
Abstract:
Observations of active regions and limb prominences often show cold, dense blobs descending with an acceleration smaller than that of free fall. The dynamics of these condensations falling in the solar corona is investigated in this paper using a simple fully ionised plasma model. We find that the presence of a heavy condensation gives rise to a dynamical rearrangement of the coronal pressure that…
▽ More
Observations of active regions and limb prominences often show cold, dense blobs descending with an acceleration smaller than that of free fall. The dynamics of these condensations falling in the solar corona is investigated in this paper using a simple fully ionised plasma model. We find that the presence of a heavy condensation gives rise to a dynamical rearrangement of the coronal pressure that results in the formation of a large pressure gradient that opposes gravity. Eventually this pressure gradient becomes so large that the blob acceleration vanishes or even points upwards. Then, the blob descent is characterised by an initial acceleration phase followed by an essentially constant velocity phase. These two stages can be identified in published time-distance diagrams of coronal rain events. Both the duration of the first stage and the velocity attained by the blob increase for larger values of the ratio of blob to coronal density, for larger blob mass, and for smaller coronal temperature. Dense blobs are characterised by a detectable density growth (up to 60% in our calculations) and by a steepening of the density in their lower part, that could lead to the formation of a shock. They also emit sound waves that could be detected as small intensity changes with periods of the order of 100 s and lasting between a few and about ten periods. Finally, the curvature of the falling path is only relevant when a very dense blob falls along inclined magnetic field lines.
△ Less
Submitted 4 November, 2013;
originally announced November 2013.
-
MHD waves in two-dimensional prominences embedded in coronal arcades
Authors:
J. Terradas,
R. Soler,
A. J. Díaz,
R. Oliver,
J. L. Ballester
Abstract:
Solar prominence models used so far in the analysis of MHD waves in such structures are quite elementary. In this work, we calculate numerically magnetohydrostatic models in two-dimensional configurations under the presence of gravity. Our interest is in models that connect the magnetic field to the photosphere and include an overlying arcade. The method used here is based on a relaxation process…
▽ More
Solar prominence models used so far in the analysis of MHD waves in such structures are quite elementary. In this work, we calculate numerically magnetohydrostatic models in two-dimensional configurations under the presence of gravity. Our interest is in models that connect the magnetic field to the photosphere and include an overlying arcade. The method used here is based on a relaxation process and requires solving the time-dependent nonlinear ideal MHD equations. Once a prominence model is obtained, we investigate the properties of MHD waves superimposed on the structure. We concentrate on motions purely two-dimensional neglecting propagation in the ignorable direction. We demonstrate how by using different numerical tools we can determine the period of oscillation of stable waves. We find that vertical oscillations, linked to fast MHD waves, are always stable and have periods in the 4-10 min range. Longitudinal oscillations, related to slow magnetoacoustic-gravity waves, have longer periods in the range of 28-40 min. These longitudinal oscillations are strongly influenced by the gravity force and become unstable for short magnetic arcades.
△ Less
Submitted 19 September, 2013;
originally announced September 2013.
-
The behavior of transverse waves in nonuniform solar flux tubes. I. Comparison of ideal and resistive results
Authors:
Roberto Soler,
Marcel Goossens,
Jaume Terradas,
Ramon Oliver
Abstract:
Magnetohydrodynamic (MHD) waves are ubiquitously observed in the solar atmosphere. Kink waves are a type of transverse MHD waves in magnetic flux tubes that are damped due to resonant absorption. The theoretical study of kink MHD waves in solar flux tubes is usually based on the simplification that the transverse variation of density is confined to a nonuniform layer much thinner than the radius o…
▽ More
Magnetohydrodynamic (MHD) waves are ubiquitously observed in the solar atmosphere. Kink waves are a type of transverse MHD waves in magnetic flux tubes that are damped due to resonant absorption. The theoretical study of kink MHD waves in solar flux tubes is usually based on the simplification that the transverse variation of density is confined to a nonuniform layer much thinner than the radius of the tube, i.e., the so-called thin boundary approximation. Here, we develop a general analytic method to compute the dispersion relation and the eigenfunctions of ideal MHD waves in pressureless flux tubes with transversely nonuniform layers of arbitrary thickness. Results for kink waves are produced and are compared with fully numerical resistive MHD eigenvalue computations in the limit of small resistivity. We find that the frequency and resonant damping rate are the same in both ideal and resistive cases. The actual results for thick nonuniform layers deviate from the behavior predicted in the thin boundary approximation and strongly depend on the shape of the nonuniform layer. The eigenfunctions in ideal MHD are very different from those in resistive MHD. The ideal eigenfunctions display a global character regardless of the thickness of the nonuniform layer, while the resistive eigenfunctions are localized around the resonance and are indistinguishable from those of ordinary resistive Alfvén modes. Consequently, the spatial distribution of wave energy in the ideal and resistive cases is dramatically different. This poses a fundamental theoretical problem with clear observational consequences.
△ Less
Submitted 13 September, 2013;
originally announced September 2013.
-
Numerical Simulations of Magnetoacoustic-Gravity Waves in the Solar Atmosphere
Authors:
K. Murawski,
A. K. Srivastava,
J. A. McLaughlin,
R. Oliver
Abstract:
We investigate the excitation of magnetoacoustic-gravity waves generated from localized pulses in the gas pressure as well as in vertical component of velocity. These pulses are initially launched at the top of the solar photosphere that is permeated by a weak magnetic field. We investigate three different configurations of the background magnetic field lines: horizontal, vertical and oblique to t…
▽ More
We investigate the excitation of magnetoacoustic-gravity waves generated from localized pulses in the gas pressure as well as in vertical component of velocity. These pulses are initially launched at the top of the solar photosphere that is permeated by a weak magnetic field. We investigate three different configurations of the background magnetic field lines: horizontal, vertical and oblique to the gravitational force. We numerically model magnetoacoustic-gravity waves by implementing a realistic (VAL-C) model of solar temperature. We solve two-dimensional ideal magnetohydrodynamic equations numerically with the use of the FLASH code to simulate the dynamics of the lower solar atmosphere. The initial pulses result in shocks at higher altitudes. Our numerical simulations reveal that a small-amplitude initial pulse can produce magnetoacoustic-gravity waves, which are later reflected from the transition region due to the large temperature gradient. The atmospheric cavities in the lower solar atmosphere are found to be the ideal places that may act as a resonator for various oscillations, including their trapping and leakage into the higher atmosphere. Our numerical simulations successfully model the excitation of such wave modes, their reflection and trapping, as well as the associated plasma dynamics.
△ Less
Submitted 29 August, 2012;
originally announced August 2012.
-
The role of Rayleigh-Taylor instabilities in filament threads
Authors:
J. Terradas,
R. Oliver,
J. L. Ballester
Abstract:
Many solar filaments and prominences show short-lived horizontal threads lying parallel to the photosphere. In this work the possible link between Rayleigh-Taylor instabilities and thread lifetimes is investigated. This is done by calculating the eigenmodes of a thread modelled as a Cartesian slab under the presence of gravity. An analytical dispersion relation is derived using the incompressible…
▽ More
Many solar filaments and prominences show short-lived horizontal threads lying parallel to the photosphere. In this work the possible link between Rayleigh-Taylor instabilities and thread lifetimes is investigated. This is done by calculating the eigenmodes of a thread modelled as a Cartesian slab under the presence of gravity. An analytical dispersion relation is derived using the incompressible assumption for the magnetohydrodynamic (MHD) perturbations. The system allows a mode that is always stable, independently of the value of the Alfvén speed in the thread. The character of this mode varies from being localised at the upper interface of the slab when the magnetic field is weak, to having a global nature and resembling the transverse kink mode when the magnetic field is strong. On the contrary, the slab model permits another mode that is unstable and localised at the lower interface when the magnetic field is weak. The growth rates of this mode can be very short, of the order of minutes for typical thread conditions. This Rayleigh-Taylor unstable mode becomes stable when the magnetic field is increased, and in the limit of strong magnetic field it is essentially a sausage magnetic mode. The gravity force might have a strong effect on the modes of oscillation of threads, depending on the value of the Alfvén speed. In the case of threads in quiescent filaments, where the Alfvén speed is presumably low, very short lifetimes are expected according to the slab model. In active region prominences, the stabilising effect of the magnetic tension might be enough to suppress the Rayleigh-Taylor instability for a wide range of wavelengths.
△ Less
Submitted 20 March, 2012;
originally announced March 2012.
-
Prominence seismology
Authors:
I. Arregui,
J. L. Ballester,
R. Oliver,
R. Soler,
J. Terradas
Abstract:
Given the difficulty in directly determining prominence physical parameters from observations, prominence seismology stands as an alternative method to probe the nature of these structures. We show recent examples of the application of magnetohydrodynamic (MHD) seismology techniques to infer physical parameters in prominence plasmas. They are based on the application of inversion techniques using…
▽ More
Given the difficulty in directly determining prominence physical parameters from observations, prominence seismology stands as an alternative method to probe the nature of these structures. We show recent examples of the application of magnetohydrodynamic (MHD) seismology techniques to infer physical parameters in prominence plasmas. They are based on the application of inversion techniques using observed periods, damping times, and plasma flow speeds of prominence thread oscillations. The contribution of Hinode to the subject has been of central importance. We show an example based on data obtained with Hinode's Solar Optical Telescope. Observations show an active region limb prominence, composed by a myriad of thin horizontal threads that flow following a path parallel to the photosphere and display synchronous vertical oscillations. The coexistence of waves and flows can be firmly established. By making use of an interpretation based on transverse MHD kink oscillations, a seismological analysis of this event is performed. It is shown that the combination of high quality Hinode observations and proper theoretical models allows flows and waves to become two useful characteristics for our understanding of the nature of solar prominences.
△ Less
Submitted 22 January, 2012;
originally announced January 2012.
-
Wave Leakage and Resonant Absorption in a Loop Embedded in a Coronal Arcade
Authors:
S. Rial,
I. Arregui,
J. Terradas,
R. Oliver,
J. L. Ballester
Abstract:
We investigate the temporal evolution of impulsively generated perturbations in a potential coronal arcade with an embedded loop. As the initial configuration we consider a coronal loop, represented by a density enhancement, which is unbounded in the ignorable direction of the arcade. The linearized time-dependent magnetohydrodynamic equations have been numerically solved in field-aligned coordina…
▽ More
We investigate the temporal evolution of impulsively generated perturbations in a potential coronal arcade with an embedded loop. As the initial configuration we consider a coronal loop, represented by a density enhancement, which is unbounded in the ignorable direction of the arcade. The linearized time-dependent magnetohydrodynamic equations have been numerically solved in field-aligned coordinates and the time evolution of the initial perturbations has been studied in the zero-beta approximation. For propagation constrained to the plane of the arcade, the considered initial perturbations do not excite trapped modes of the system. This weakness of the model is overcome by the inclusion of wave propagation in the ignorable direction. The inclusion of perpendicular propagation produces two main results. First, damping by wave leakage is less efficient because the loop is able to act as a wave trap of vertical oscillations. Second, the consideration of an inhomogeneous corona enables the resonant damping of vertical oscillations and the energy transfer from the interior of the loop to the external coronal medium.
△ Less
Submitted 19 January, 2012;
originally announced January 2012.
-
Rossby waves and polar spots in rapidly rotating stars: Implications for stellar wind evolution
Authors:
T. V. Zaqarashvili,
R. Oliver,
J. L. Ballester,
M. Carbonell,
M. L. Khodachenko,
H. Lammer,
M. Leitzinger,
P. Odert
Abstract:
Rapidly rotating stars show short-period oscillations in magnetic activity and polar appearance of starspots. The aim of this paper is to study large-scale shallow water waves in the tachoclines of rapidly rotating stars and their connection to the periodicity and the formation of starspots at high latitudes. Shallow-water magnetohydrodynamic equations were used to study the dynamics of large-scal…
▽ More
Rapidly rotating stars show short-period oscillations in magnetic activity and polar appearance of starspots. The aim of this paper is to study large-scale shallow water waves in the tachoclines of rapidly rotating stars and their connection to the periodicity and the formation of starspots at high latitudes. Shallow-water magnetohydrodynamic equations were used to study the dynamics of large-scale waves at the rapidly rotating stellar tachoclines in the presence of toroidal magnetic field. Dispersion relations and latitudinal distribution of wave modes were derived. We found that low-frequency magnetic Rossby waves tend to be located at poles, but high-frequency magnetic Poincare waves are concentrated near the equator in rapidly rotating stars. These results have important implications for the evolution of the stellar wind in young Sun-like stars. Unstable magnetic Rossby waves may lead to the local enhancement of magnetic flux at high latitudes of tachoclines in rapidly rotating stars. The enhanced magnetic flux may rise upwards owing to the magnetic buoyancy in the form of tubes and appear as starspots at polar regions. Magnetic Rossby waves may also cause observed short-term periodicity in the stellar magnetic activity.
△ Less
Submitted 18 July, 2011;
originally announced July 2011.
-
Quasi-biennial oscillations in the solar tachocline caused by magnetic Rossby wave instabilities
Authors:
T. V. Zaqarashvili,
M. carbonell,
R. Oliver,
J. L. Ballester
Abstract:
Quasi-biennial oscillations (QBO) are frequently observed in the solar activity indices. However, no clear physical mechanism for the observed variations has been suggested so far. Here we study the stability of magnetic Rossby waves in the solar tachocline using the shallow water magnetohydrodynamic approximation. Our analysis shows that the combination of typical differential rotation and a toro…
▽ More
Quasi-biennial oscillations (QBO) are frequently observed in the solar activity indices. However, no clear physical mechanism for the observed variations has been suggested so far. Here we study the stability of magnetic Rossby waves in the solar tachocline using the shallow water magnetohydrodynamic approximation. Our analysis shows that the combination of typical differential rotation and a toroidal magnetic field with a strength > 10^5 G triggers the instability of the m=1 magnetic Rossby wave harmonic with a period of 2 years. This harmonic is antisymmetric with respect to the equator and its period (and growth rate) depends on the differential rotation parameters and the magnetic field strength. The oscillations may cause a periodic magnetic flux emergence at the solar surface and consequently may lead to the observed QBO in the solar activity features. The period of QBO may change throughout the cycle, and from cycle to cycle, due to variations of the mean magnetic field and differential rotation in the tachocline.
△ Less
Submitted 5 November, 2010;
originally announced November 2010.
-
Spatial Damping of Propagating Kink Waves in Prominence Threads
Authors:
R. Soler,
R. Oliver,
J. L. Ballester
Abstract:
Transverse oscillations and propagating waves are frequently observed in threads of solar prominences/filaments and have been interpreted as kink magnetohydrodynamic (MHD) modes. We investigate the spatial damping of propagating kink MHD waves in transversely nonuniform and partially ionized prominence threads. Resonant absorption and ion-neutral collisions (Cowling's diffusion) are the damping…
▽ More
Transverse oscillations and propagating waves are frequently observed in threads of solar prominences/filaments and have been interpreted as kink magnetohydrodynamic (MHD) modes. We investigate the spatial damping of propagating kink MHD waves in transversely nonuniform and partially ionized prominence threads. Resonant absorption and ion-neutral collisions (Cowling's diffusion) are the damping mechanisms taken into account. The dispersion relation of resonant kink waves in a partially ionized magnetic flux tube is numerically solved by considering prominence conditions. Analytical expressions of the wavelength and damping length as functions of the kink mode frequency are obtained in the Thin Tube and Thin Boundary approximations. For typically reported periods of thread oscillations, resonant absorption is an efficient mechanism for the kink mode spatial damping, while ion-neutral collisions have a minor role. Cowling's diffusion dominates both the propagation and damping for periods much shorter than those observed. Resonant absorption may explain the observed spatial damping of kink waves in prominence threads. The transverse inhomogeneity length scale of the threads can be estimated by comparing the observed wavelengths and damping lengths with the theoretically predicted values. However, the ignorance of the form of the density profile in the transversely nonuniform layer introduces inaccuracies in the determination of the inhomogeneity length scale.
△ Less
Submitted 24 September, 2010;
originally announced September 2010.
-
Seismology of Standing Kink Oscillations of Solar Prominence Fine Structures
Authors:
Roberto Soler,
Inigo Arregui,
Ramon Oliver,
Jose Luis Ballester
Abstract:
We investigate standing kink magnetohydrodynamic (MHD) oscillations in a prominence fine structure modeled as a straight and cylindrical magnetic tube only partially filled with the prominence material, and with its ends fixed at two rigid walls representing the solar photosphere. The prominence plasma is partially ionized and a transverse inhomogeneous transitional layer is included between the p…
▽ More
We investigate standing kink magnetohydrodynamic (MHD) oscillations in a prominence fine structure modeled as a straight and cylindrical magnetic tube only partially filled with the prominence material, and with its ends fixed at two rigid walls representing the solar photosphere. The prominence plasma is partially ionized and a transverse inhomogeneous transitional layer is included between the prominence thread and the coronal medium. Thus, ion-neutral collisions and resonant absorption are the considered damping mechanisms. Approximate analytical expressions of the period, the damping time, and their ratio are derived for the fundamental mode in the thin tube and thin boundary approximations. We find that the dominant damping mechanism is resonant absorption, which provides damping ratios in agreement with the observations, whereas ion-neutral collisions are irrelevant for the damping. The values of the damping ratio are independent of both the prominence thread length and its position within the magnetic tube, and coincide with the values for a tube fully filled with the prominence plasma. The implications of our results in the context of the MHD seismology technique are discussed, pointing out that the reported short-period (2 - 10 min) and short-wavelength (700 - 8,000 km) thread oscillations may not be consistent with a standing mode interpretation and could be related to propagating waves. Finally, we show that the inversion of some prominence physical parameters, e.g., Alfvén speed, magnetic field strength, transverse inhomogeneity length-scale, etc., is possible using observationally determined values of the period and damping time of the oscillations along with the analytical approximations of these quantities.
△ Less
Submitted 30 August, 2010; v1 submitted 12 July, 2010;
originally announced July 2010.
-
Three-Dimensional Propagation of Magnetohydrodynamic Waves in Solar Coronal Arcades
Authors:
S. Rial,
I. Arregui,
J. Terradas,
R. Oliver,
J. L. Ballester
Abstract:
We numerically investigate the excitation and temporal evolution of oscillations in a two-dimensional coronal arcade by including the three-dimensional propagation of perturbations. The time evolution of impulsively generated perturbations is studied by solving the linear, ideal magnetohydrodynamic (MHD) equations in the zero-beta approximation. As we neglect gas pressure the slow mode is absent…
▽ More
We numerically investigate the excitation and temporal evolution of oscillations in a two-dimensional coronal arcade by including the three-dimensional propagation of perturbations. The time evolution of impulsively generated perturbations is studied by solving the linear, ideal magnetohydrodynamic (MHD) equations in the zero-beta approximation. As we neglect gas pressure the slow mode is absent and therefore only coupled MHD fast and Alfven modes remain. Two types of numerical experiments are performed. First, the resonant wave energy transfer between a fast normal mode of the system and local Alfven waves is analyzed. It is seen how, because of resonant coupling, the fast wave with global character transfers its energy to Alfvenic oscillations localized around a particular magnetic surface within the arcade, thus producing the damping of the initial fast MHD mode. Second, the time evolution of a localized impulsive excitation, trying to mimic a nearby coronal disturbance, is considered. In this case, the generated fast wavefront leaves its energy on several magnetic surfaces within the arcade. The system is therefore able to trap energy in the form of Alfvenic oscillations, even in the absence of a density enhancement such as that of a coronal loop. These local oscillations are subsequently phase-mixed to smaller spatial scales. The amount of wave energy trapped by the system via wave energy conversion strongly depends on the wavelength of perturbations in the perpendicular direction, but is almost independent from the ratio of the magnetic to density scale heights.
△ Less
Submitted 12 February, 2010; v1 submitted 2 February, 2010;
originally announced February 2010.
-
The spatial damping of magnetohydrodynamic waves in a flowing partially ionised prominence plasma
Authors:
M. Carbonell,
P. Forteza,
R. Oliver,
J. L. Ballester
Abstract:
Solar prominences are partially ionised plasmas displaying flows and oscillations. These oscillations show time and spatial damping and, commonly, have been explained in terms of magnetohydrodynamic (MHD) waves. We study the spatial damping of linear non-adiabatic MHD waves in a flowing partially ionised plasma, having prominence-like physical properties. We consider single fluid equations for a…
▽ More
Solar prominences are partially ionised plasmas displaying flows and oscillations. These oscillations show time and spatial damping and, commonly, have been explained in terms of magnetohydrodynamic (MHD) waves. We study the spatial damping of linear non-adiabatic MHD waves in a flowing partially ionised plasma, having prominence-like physical properties. We consider single fluid equations for a partially ionised hydrogen plasma including in the energy equation optically thin radiation, thermal conduction by electrons and neutrals, and heating. Keeping the frequency real and fixed, we have solved the obtained dispersion relations for the complex wavenumber, k, and have analysed the behaviour of the damping length, wavelength and the ratio of the damping length to the wavelength, versus period, for Alfven, fast, slow and thermal waves.
△ Less
Submitted 27 January, 2010;
originally announced January 2010.