Quasi-biennial oscillations in the cross-correlation of properties of macrospicules

Jets, whatever small (e.g. spicules) or large (e.g. macrospicules) their size, may play a key role in momentum and energy transport from photosphere to chromosphere and at least to the low corona. Here, we investigate the properties of abundant, large-scale dynamic jets observable in the solar atmosphere: the macrospicules (MS). These jets are observationally more distinct phenomena than their little, and perhaps more ubiquitous, cousins, the spicules. Investigation of long-term variation of the properties of macrospicules may help to a better understanding of their underlying physics of generation and role in coronal heating. Taking advantage of the high temporal and spatial resolution of the Solar Dynamics Observatory, a new dataset, with several hundreds of macrospicules, was constructed encompassing a period of observations over six years. Here, we analyse the measured properties and relations between these properties of macrospicules as function of time during the observed time interval. We found that cross-correlations of several of these macrospicule properties display a strong oscillatory pattern. Next, wavelet analysis is used to provide more detailed information about the temporal behaviour of the various properties of MS. For coronal hole macrospicules, a significant peak is found at around 2-year period. This peak also exists partially or is shifted to longer period, in the case of quiet Sun macrospicules. These observed findings may be rooted in the underlying mechanism generating the solar magnetic field, i.e. the global solar dynamo.     

Analysis of a long-duration AR throughout five solar rotations: Magnetic properties and ejective events

Coronal mass ejections (CMEs), which are among the most magnificent solar eruptions, are a major driver of space weather and can thus affect diverse human technologies. Different processes have been proposed to explain the initiation and release of CMEs from solar active regions (ARs), without reaching consensus on which is the predominant scenario, and thus rendering impossible to accurately predict when a CME is going to erupt from a given AR. To investigate AR magnetic properties that favor CMEs production, we employ multi-spacecraft data to analyze a long duration AR (NOAA 11089, 11100, 11106, 11112 and 11121) throughout its complete lifetime, spanning five Carrington rotations from July to November 2010. We use data from the Solar Dynamics Observatory to study the evolution of the AR magnetic properties during the five near-side passages, and a proxy to follow the magnetic flux changes when no magnetograms are available, i.e. during far-side transits. The ejectivity is studied by characterizing the angular widths, speeds and masses of 108 CMEs that we associated to the AR, when examining a 124-day period. Such an ejectivity tracking was possible thanks to the multi-viewpoint images provided by the Solar-Terrestrial Relations Observatory and Solar and Heliospheric Observatory in a quasi-quadrature configuration. We also inspected the X-ray flares registered by the GOES satellite and found 162 to be associated to the AR under study. Given the substantial number of ejections studied, we use a statistical approach instead of a single-event analysis. We found three well defined periods of very high CMEs activity and two periods with no mass ejections that are preceded or accompanied by characteristic changes in the AR magnetic flux, free magnetic energy and/or presence of electric currents. Our large sample of CMEs and long term study of a single AR, provide further evidence relating AR magnetic activity to CME and Flare production.     

A statistical study of CME-associated flare during the solar cycle 24

We investigate on the relationship between flares and coronal mass ejections (CMEs) in which a flare started before and after the CME events which differ in their physical properties, indicating potentially different initiation mechanisms. The physical properties of two types flare-correlated CME remain an interesting and important question in space weather. We study the relationship between flares and CMEs using a different approach requiring both temporal and spatial constraints during the period from December 1, 2008 to April 30, 2017 in which the CMEs data were acquired by SOHO/LASCO (Solar and Heliospheric Observatory/Large Angle Spectrometric Coronagraph) over the solar cycle 24. The soft X-ray flare flux data, such as flare class, location, onset time and integrated flux, are collected from Geostationary Environmental satellite (GOES) and XRT Flare catalogs. We selected 307 CMEs-flares pairs applying simultaneously temporal and spatial constraints in all events for the distinguish between two associated CME-flare types. We study the correlated properties of coincident flares and CMEs during this period, specifically separating the sample into two types: flares that precede a CME and flares that follow a CME. We found an opposite correlation relationship between the acceleration and velocity of CMEs in the After- and Before-CMEs events. We found a log-log relation between the width and mass of CMEs in the two associated types. The CMEs and flares properties show that there were significant differences in all physical parameters such as (mass, angular width, kinetic energy, speed and acceleration) between two flare-associated CME types.     

Topological and statistical properties of nonlinear force-free fields

We use our semi-analytic solution of the nonlinear force-free field equation to construct three-dimensional magnetic fields that are applicable to the solar corona and study their statistical properties for estimating the degree of braiding exhibited by these fields. We present a new formula for calculating the winding number and compare it with the formula for the crossing number. The comparison is shown for a toy model of two helices and for realistic cases of nonlinear force-free fields; conceptually the formulae are nearly the same but the resulting distributions calculated for a given topology can be different. We also calculate linkages, which are useful topological quantities that are independent measures of the contribution of magnetic braiding to the total free energy and relative helicity of the field. Finally, we derive new analytical bounds for the free energy and relative helicity for the field configurations in terms of the linking number. These bounds will be of utility in estimating the braided energy available for nano-flares or for eruptions.     

A MHD-turbulence model for solar corona

The disposition of energy in the solar corona has always been a problem of great interest. It remains an open question how the low temperature photosphere supports the occurence of solar extreme phenomena. In this work, a turbulent heating mechanism for the solar corona through the framework of reduced magnetohydrodynamics (RMHD) is proposed. Two-dimensional incompressible long time simulations of the average energy disposition have been carried out with the aim to reveal the characteristics of the long time statistical behavior of a two-dimensional cross-section of a coronal loop and the importance of the photospheric time scales in the understanding of the underlying mechanisms. It was found that for a slow, shear type photospheric driving the magnetic field in the loop self-organizes at large scales via an inverse MHD cascade. The system undergoes three distinct evolutionary phases. The initial forcing conditions are quickly “forgotten” giving way to an inverse cascade accompanied with and ending up to electric current dissipation. Scaling laws are being proposed in order to quantify the nonlinearity of the system response which seems to become more impulsive for decreasing resistivity. It is also shown that few, if any, qualitative changes in the above results occur by increasing spatial resolution.     

Solar physics with the Square Kilometre Array

The Square Kilometre Array (SKA) will be the largest radio telescope ever built, aiming to provide collecting area larger than 1?km². The SKA will have two independent instruments, SKA-LOW comprising of dipoles organized as aperture arrays in Australia and SKA-MID comprising of dishes in South Africa. Currently the phase-1 of SKA, referred to as SKA1, is in its late design stage and construction is expected to start in 2020. Both SKA1-LOW (frequency range of 50–350?MHz) and SKA1-MID Bands 1, 2, and 5 (frequency ranges of 350–1050, 950–1760, and 4600–15,300?MHz, respectively) are important for solar observations. In this paper we present SKA’s unique capabilities in terms of spatial, spectral, and temporal resolution, as well as sensitivity and show that they have the potential to provide major new insights in solar physics topics of capital importance including (i) the structure and evolution of the solar corona, (ii) coronal heating, (iii) solar flare dynamics including particle acceleration and transport, (iv) the dynamics and structure of coronal mass ejections, and (v) the solar aspects of space weather. Observations of the Sun jointly with the new generation of ground-based and space-borne instruments promise unprecedented discoveries.     

Driving major solar flares and eruptions: A review

This review focuses on the processes that energize and trigger M- and X-class solar flares and associated flux-rope destabilizations. Numerical modeling of specific solar regions is hampered by uncertain coronal-field reconstructions and by poorly understood magnetic reconnection; these limitations result in uncertain estimates of field topology, energy, and helicity. The primary advances in understanding field destabilizations therefore come from the combination of generic numerical experiments with interpretation of sets of observations. These suggest a critical role for the emergence of twisted flux ropes into pre-existing strong field for many, if not all, of the active regions that produce M- or X-class flares. The flux and internal twist of the emerging ropes appear to play as important a role in determining whether an eruption will develop predominantly as flare, confined eruption, or CME, as do the properties of the embedding field. Based on reviewed literature, I outline a scenario for major flares and eruptions that combines flux-rope emergence, mass draining, near-surface reconnection, and the interaction with the surrounding field. Whether deterministic forecasting is in principle possible remains to be seen: to date no reliable such forecasts can be made. Large-sample studies based on long-duration, comprehensive observations of active regions from their emergence through their flaring phase are needed to help us better understand these complex phenomena.     

Variations of the solar acoustic high-degree mode frequencies over solar cycle 23

Using full-disk observations obtained with the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory (SOHO) spacecraft, we present variations of the solar acoustic mode frequencies caused by the solar activity cycle. High-degree (100 < ? < 900) solar acoustic modes were analyzed using global helioseismology analysis techniques over most of solar cycle 23. We followed the methodology described in details in [Korzennik, S.G., Rabello-Soares, M.C., Schou, J. On the determination of Michelson Doppler Imager high-degree mode frequencies. ApJ 602, 481–515, 2004] to infer unbiased estimates of high-degree mode parameters ([see also Rabello-Soares, M.C., Korzennik, S.G., Schou, J. High-degree mode frequencies: changes with solar cycle. ESA SP-624, 2006]). We have removed most of the known instrumental and observational effects that affect specifically high-degree modes. We show that the high-degree changes are in good agreement with the medium-degree results, except for years when the instrument was highly defocused. We analyzed and discuss the effect of defocusing on high-degree estimation. Our results for high-degree modes confirm that the frequency shift scaled by the relative mode inertia is a function of frequency and it is independent of degree.     

The dynamic chromosphere: Pushing the boundaries of observations and models

The interface between the bright solar surface and the million-degree corona continues to hold the key to many unsolved problems in solar physics. Advances in instrumentation now allow us to observe the dynamic structures of the solar chromosphere down to less than $0_{.}^{″}1$ with cadences of just a few seconds and in multiple polarisation states. Such observational progress has been matched by the ever-increasing sophistication of numerical models, which have become necessary to interpret the complex observations. With an emphasis on the quiet Sun, I will review recent progress in the observation and modelling of the chromosphere. Models have come a long way from 1D static atmospheres, but their predictions still fail to reproduce several key observed features. Nevertheless, they have given us invaluable insight into the physical processes that energise the atmosphere. With more physics being added to models, the gap between predictions and observations is narrowing. With the next generation of solar observatories just around the corner, the big question is: will they close the gap?     

The damping of small-amplitude oscillations in quiescent prominences

The presence of small-amplitude oscillations in prominences is well-known from long time ago. These oscillations, whose exciters are still unknown, seem to be of local nature and are interpreted in terms of magnetohydrodynamic (MHD) waves. During last years, observational evidence about the damping of these oscillations has grown and several mechanisms able to damp these oscillations have been the subject of intense theoretical modelling. Among them, the most efficient seem to be radiative cooling and ion-neutral collisions. Radiative cooling is able to damp slow MHD waves efficiently, while ion-neutral collisions, in partially ionised plasmas like those of solar prominences, can also damp fast MHD waves. In this paper, we plan to summarize our current knowledge about the time and spatial damping of small-amplitude oscillations in prominences.     

Magnetic connectivity and solar energetic proton event intensity profiles at deka-MeV energy

We present an analysis of the time-intensity profiles of 25 solar energetic proton events at 18.2 MeV, modelled by fitting an analytical function form (a modified Weibull function) to the observed intensities. Additionally relying on previous work that characterized the magnetic connectivity between the event-related solar flare and the observer in these events with three angular parameters, we investigate the fit function parameters, the connectivity parameters, and the iron-to-carbon ratio of the events for dependencies and correlations. We find that the fit parameter controlling the basic shape of the profile (parameter a) is not clearly dependent on the connectivity parameters or the Fe/C ratio, suggesting that the profile shapes of neither well and weakly connected nor generally “impulsive” and “gradual” events differ systematically during the early stages of the event at 1 AU. In contrast, the time scaling of the fit function (parameter b) is at least moderately correlated with both the magnetic connectivity parameters and the Fe/C ratio, in that well-connected and iron-rich events are typically shorter in relative duration than weakly connected and nominal-abundance events; intensity rise times display a similar correlation with the connectivity parameters. We interpret the former result as following from the combined effect of various transport processes acting on the particles in interplanetary space, while the latter is essentially consistent with established knowledge regarding the observed dependence of the time-intensity profile shapes of solar energetic particle events on their magnetic connectivity and heavy ion abundances. The desirability of modelling the particle transport effects in detail and extending the analysis to cover higher energies is indicated.     

Observational study of magnetic chirality of solar active regions

We present the evolution of magnetic field and relationship with the magnetic (current) helicity in solar active regions from a series of photospheric vector magnetograms obtained at Huairou Solar Observing Station near Beijing, and also longitudinal magnetograms by MDI of SOHO, white light and 171 Å images by TRACE and soft X-ray images by Yohkoh.The conclusions in the analysis of the formation process of complex and delta magnetic configuration in some super active regions are the following: (1) The magnetic shear and gradient provide the non-potentiality of the magnetic field of active regions reflecting the existence of electric current. (2) Some of large-scale delta active regions could be due to the emergence of highly sheared non-potential magnetic flux bundles from the subatmosphere with amount of magnetic helicity, in addition to the emergence of twisted magnetic ropes. (3) We also present some results on the study of the magnetic (current) helicity in solar active regions.     

Model of a fluxtube with a twisted magnetic field in the stratified solar atmosphere

We build a single vertical straight magnetic fluxtube spanning the solar photosphere and the transition region which does not expand with height. We assume that the fluxtube containing twisted magnetic fields is in magnetohydrostatic equilibrium within a realistic stratified atmosphere subject to solar gravity. Incorporating specific forms of current density and gas pressure in the Grad–Shafranov equation, we solve the magnetic flux function, and find it to be separable with a Coulomb wave function in radial direction while the vertical part of the solution decreases exponentially. We employ improved fluxtube boundary conditions and take a realistic ambient external pressure for the photosphere to transition region, to derive a family of solutions for reasonable values of the fluxtube radius and magnetic field strength at the base of the axis that are the free parameters in our model. We find that our model estimates are consistent with the magnetic field strength and the radii of Magnetic bright points (MBPs) as estimated from observations. We also derive thermodynamic quantities inside the fluxtube.     

“Oblique” Bernstein modes in solar preflare plasma: Generation of second harmonics

Process of second harmonics generation due to development of corresponding instability has been investigated for pure electron weakly oblique Bernstein mode. This mode was supposed to be modified by taking into account the influence of pair Coulomb collisions and weak large-scale electric field in flare loop. Investigated area was located near the loop foot-point in the “lower–middle” chromosphere of active region. It has been shown, that for the Fontenla–Avrett–Loeser model of solar atmosphere the investigated process of second harmonics generation starts at the extremely low threshold values of subdreicer electric field, well before the beginning of “preheating” phase of flare process.     

From solar nanoflares to stellar giant flares: Scaling laws and non-implications for coronal heating

In this study we explore physical scaling laws applied to solar nanoflares, microflares, and large flares, as well as to stellar giant flares. Solar flare phenomena exhibit a fractal volume scaling, V(L)  L^1.9, with L being the flare loop length scale, which explains the observed correlation between the total emission measure EM_p and flare peak temperature T_p in both solar and stellar flares. However, the detected stellar flares have higher emission measures EM_p than solar flares at the same flare peak temperature T_p, which can be explained by a higher electron density that is caused by shorter heating scale height ratios s_H/L ≈ 0.04–0.1. Using these scaling laws we calculate the total radiated flare energies E_X and thermal flare energies E_T and find that the total counts C are a good proxy for both parameters. Comparing the energies of solar and stellar flares we find that even the smallest observed stellar flares exceed the largest solar flares, and thus their observed frequency distributions are hypothetically affected by an upper cutoff caused by the maximum active region size limit. The powerlaw slopes fitted near the upper cutoff can then not reliably be extrapolated to the microflare regime to evaluate their contribution to coronal heating.     

Power-law spectra of energetic electrons in solar flares from the maximum entropy and dimensional considerations

The maximum entropy formalism and dimensional analysis are used to derive a power-law spectrum of accelerated electrons in impulsive solar flares, where the particles can contain a significant fraction of the total flare energy. Entropy considerations are used to derive a power-law spectrum for a particle distribution characterised by its order of magnitude of energy. The derivation extends an earlier one-dimensional argument to the case of an isotropic three-dimensional particle distribution. Dimensional arguments employ the idea that the spectrum should reflect a balance between the processes of energy input into the corona and energy dissipation in solar flares. The governing parameters are suggested on theoretical grounds and shown to be consistent with solar flare observations. The flare electron flux, differential in the non-relativistic electron kinetic energy E, is predicted to scale as $E^{-3}$. This scaling is in agreement with RHESSI measurements of the hard X-ray flux that is generated by deka-keV electrons, accelerated in intense solar flares.     

Asymmetric precipitation in a coronal loop as explanation of a singular observed spectrum

Almost 10 years of solar submillimeter observations have shown new aspects of solar activity, such as the presence of rapid solar spikes associated with the launch of coronal mass ejections and an increasing submillimeter spectral component in flares. We analyse the singular microwave–submillimeter spectrum of an M class solar flare on 20 December, 2002. Flux density observations measured by Sun patrol telescopes and the Solar Submillimeter Telescope are used to build the radio spectrum, which is fitted using Ramaty’s code. At submillimeter frequencies the spectrum shows a component different from the microwave classical burst. The fitting is achieved proposing two homogeneous sources of emission. This theoretical fitting is in agreement with differential precipitation through a magnetically asymmetric loop or set of loops. From a coronal magnetic field model we infer an asymmetric magnetic structure at the flare location. The model proposed to quantify the differential precipitation rates due to the asymmetry results in a total precipitation ratio Q₂/Q₁≈10⁴–10⁵, where Q₁(Q₂) represents the total precipitation in the loop foot with the high (low) magnetic field intensity. This ratio agrees with the electron total number ratio of the two sources proposed to fit the radio spectrum.     

Multi-scale reconnections in a complex CME

A series of three flares of GOES class M, M and C, and a CME were observed on 20 January 2004 occurring in close succession in NOAA 10540. Types II, III, and N radio bursts were associated. We use the combined observations from TRACE, EIT, Hα images from Kwasan Observatory, MDI magnetograms, GOES, and radio observations from Culgoora and Wind/ WAVES to understand the complex development of this event. We reach three main conclusions. First, we link the first two impulsive flares to tether-cutting reconnections and the launch of the CME. This complex observation shows that impulsive quadrupolar flares can be eruptive. Second, we relate the last of the flares, an LDE, to the relaxation phase following forced reconnections between the erupting flux rope and neighbouring magnetic field lines, when reconnection reverses and restores some of the pre-eruption magnetic connectivities. Finally, we show that reconnection with the magnetic structure of a previous CME launched about 8 h earlier injects electrons into open field lines having a local dip and apex (located at about six solar radii height). This is observed as an N-burst at decametre radio wavelengths. The dipped shape of these field lines is due to large-scale magnetic reconnection between expanding magnetic loops and open field lines of a neighbouring streamer. This particular situation explains why this is the first N-burst ever observed at long radio wavelengths.     

How far are we from a ‘Standard Model’ of the solar dynamo?

Over the last few years, dynamo theorists seem to be converging on a basic scenario as to how the solar dynamo operates. The strong toroidal component of the magnetic field is produced in the tachocline, from where it rises due to magnetic buoyancy to produce active regions at the solar surface. The decay of tilted bipolar active regions at the surface gives rise to the poloidal component, which is first advected poleward by the meridional circulation and then taken below the surface to the tachocline where it can be stretched to produce the toroidal component. The mathematical formulation of this basic model, however, involves the specification of some parameters which are still uncertain. We review these remaining uncertainties which have resulted in disagreements amongst various research groups and have made it impossible to still arrive at something that can be called a standard model of the solar dynamo.     

Statistical study of radio drifting pulsation structures with associated CMEs and other observations

Solar radio burst, especially the fine structures (FSs) and the drifting pulsation structures (DPSs), may be used as an important diagnostics tool to draw the evolution map of the flare loop in the initial phase of solar flares. In this work, 52 radio events were found accompanying with DPSs. They were all observed with the Solar Radio Spectrometers (0.625–7.6 GHz) of China during 1998–2004. Combining the radio observations with LASCO-C2, Goes-8 SXR, H_α, EUV and Trace observations, we analyzed all these events and obtained some statistic conclusions: First, 88% DPSs take place at the initial phase of the radio burst, and their rich spectrum characteristics are helpful to understand the events further. Second, 83% DPSs are associated with CMEs or ejection events, and all the events are accompanied by Goes SXR flare. Third, for CMEs and DPSs, which take the first step, there is no significant predominance of either of them. The relationship between the DPSs and CMEs is still not clear in this study because of the lack of spatial resolution in the centimeter–decimeter band. However, the EIT or Trace ejection happened during the onset/end time of DPSs. They are signatures of the initial phase of CMEs. Two events will be illustrated to explain this.