R. A. M. Van der Linden

AUTOMATED LASCO CME CATALOG FOR SOLAR CYCLE 23: ARE CMES SCALE INVARIANT?

In this paper, we present the first automatically constructed LASCO coronal mass ejection (CME) catalog, a result of the application of the Computer Aided CME Tracking software (CACTus) on the LASCO archive during the interval 1997 September-2007 January. We have studied the CME characteristics and have compared them with similar results obtained by manual detection (CDAW CME catalog). On average, CACTus detects less than two events per day during solar minimum, up to eight events during maximum, nearly half of them being narrow (<20?). Assuming a correction factor, we find that the CACTus CME rate is surprisingly consistent with CME rates found during the past 30 years. The CACTus statistics show that small-scale outflow is ubiquitously observed in the outer corona. The majority of CACTus-only events are narrow transients related to previous CME activity or to intensity variations in the slow solar wind, reflecting its turbulent nature. A significant fraction (about 15%) of CACTus-only events were identified as independent events, thus not related to other CME activity. The CACTus CME width distribution is essentially scale invariant in angular span over a range of scales from 20? to 120? while previous catalogs present a broad maximum around 30?. The possibility that the size of coronal mass outflows follow a power-law distribution could indicate that no typical CME size exists, i.e., that the narrow transients are not different from the larger well defined CMEs.

Heating the corona by nanoflares: simulations of energy release triggered by a kink instability

Context. The heating of solar coronal plasma to millions of degrees is likely to be due to the superposition of many small energy-releasing events, known as nanoflares. Nanoflares dissipate magnetic energy through magnetic reconnection. Aims. A model has been recently proposed in which nanoflare-like heating naturally arises, with a sequence of dissipation events of various magnitudes. It is proposed that heating is triggered by the onset of ideal instability, with energy release occurring in the nonlinear phase due to fast magnetic reconnection. The aim is to use numerical simulations to investigate this heating process. Methods. Three-dimensional magnetohydrodynamic numerical simulations of energy release are presented for a cylindrical coronal loop model. Initial equilibrium magnetic-field profiles are chosen to be linearly unstable, with a two-layer parameterisation of the current profile. The results are compared with calculations of linear instability, with line-tying, which are extended to account for a potential field layer surrounding the loop. The energy release is also compared with predictions that the field relaxes to a state of minimum magnetic energy with conserved magnetic helicity (a constant a force-free field). Results. The loop initially develops a helical kink, whose structure and growth rate are generally in accordance with linear stability theory, and subsequently a current sheet forms. During this phase, there is a burst of kinetic energy while the magnetic energy decays. A new relaxed equilibrium is established that corresponds quite closely to a constant a field. The fraction of stored magnetic energy released depends strongly on the initial current profile, which agrees with the predictions of relaxation theory. Conclusions. Energy dissipation events in a coronal loop are triggered by the onset of ideal kink instability. Magnetic energy is dissipated, leading to large or small heating events according to the initial current profile.

Coronal heating by magnetic reconnection in loops with zero net current

Context. The paper is concerned with heating of the solar corona by nanoflares: a superposition of small transient events in which stored magnetic energy is dissipated by magnetic reconnection. It is proposed that heating occurs in the nonlinear phase of an ideal kink instability, where magnetic reconnection leads to relaxation to a state of minimum magnetic energy. Aims. The aim is to investigate the nonlinear aspects of magnetic relaxation on a current loop with zero net axial current. The dynamical processes leading to the establishment of a relaxed state are explored. The efficiency of loop heating is investigated. Methods. A 3D magnetohydrodynamic numerical code is used to simulate the evolution of coronal loops which are initially in ideally unstable equilibrium. The initial states have zero net current. The results are interpreted by comparison both with linear stability analysis and with helicity-conserving relaxation theory. Results. The disturbance due to the unstable mode is strongly radially confined when the loop carries zero net current. Strong current sheets are still formed in the nonlinear phase with dissipation of magnetic energy by fast reconnection. The nonlinear development consists first of reconnection in a large scale current sheet, which forms near the quasi-resonant surface of the equilibrium field. Subsequently, the current sheet extends and then fragments, leading to multiple reconnections and effective relaxation to a constant α field. Conclusions. Magnetic reconnection is triggered in the nonlinear phase of kink instability in loops with zero net current. Initially, reconnection occurs in a single current sheet, which then fragments into multiple reconnection sites, allowing almost full relaxation to the minimum energy state. The loop is heated to high temperatures throughout its volume.

A nanoflare distribution generated by repeated relaxations triggered by kink instability

Context. It is thought likely that vast numbers of nanoflares are responsible for the corona having a temperature of millions of degrees. Current observational technologies lack the resolving power to confirm the nanoflare hypothesis. An alternative approach is to construct a magnetohydrodynamic coronal loop model that has the ability to predict nanoflare energy distributions. Aims. This paper presents the initial results generated by a coronal loop model that flares whenever it becomes unstable to an ideal MHD kink mode. A feature of the model is that it predicts heating events with a range of sizes, depending on where the instability threshold for linear kink modes is encountered. The aims are to calculate the distribution of event energies and to investigate whether kink instability can be predicted from a single parameter. Methods. The loop is represented as a straight line-tied cylinder. The twisting caused by random photospheric motions is captured by two parameters, representing the ratio of current density to field strength for specific regions of the loop. Instability onset is mapped as a closed boundary in the 2D parameter space. Dissipation of the loops magnetic energy begins during the nonlinear stage of the instability, which develops as a consequence of current sheet reconnection. After flaring, the loop evolves to the state of lowest energy where, in accordance with relaxation theory, the ratio of current to field is constant throughout the loop and helicity is conserved. Results. There exists substantial variation in the radial magnetic twist profiles for the loop states along the instability threshold. These results suggest that instability cannot be predicted by any simple twist-derived property reaching a critical value. The model is applied such that the loop undergoes repeated episodes of instability followed by energy-releasing relaxation. Hence, an energy distribution of the nanoflares produced is collated. This paper also presents the calculated relaxation states and energy releases for all instability threshold points. Conclusions. The final energy distribution features two nanoflare populations that follow different power laws. The power law index for the higher energy population is more than sufficient for coronal heating.

A two-dimensional magnetohydrodynamic stability model for helicity-injected devices with open flux

Models of the ideal magnetohydrodynamic (MHD) stability of spheromaks and spherical tokamaks are presented, including the effects of current on the open flux which plays a key role in helicity-injected current drive. The stability of spheromak equilibria with both open and closed flux and realistic current profiles representative of helicity-injected state is investigated, where a kink instability in the open flux is shown to dominate a tilt mode in the closed flux as the open flux current density is increased. A previous one-dimensional model is extended to more realistic two-dimensional equilibria which properly incorporate a region of closed magnetic flux as well as open flux penetrating the boundaries at electrodes. A new stability code SCOTS has been developed and benchmarked which can determine the growth rates of ideal MHD modes in this geometry. The coordinate system for this code has been developed such that it extends smoothly across the separatrix between closed and open flux, thus not imposing...

The Flare-Energy Distributions Generated by Kink-Unstable Ensembles of Zero-Net-Current Coronal Loops

It has been proposed that the million-degree temperature of the corona is due to the combined effect of barely detectable energy releases, called nanoflares, that occur throughout the solar atmosphere. Unfortunately, the nanoflare density and brightness implied by this hypothesis means that conclusive verification is beyond present observational abilities. Nevertheless, we investigate the plausibility of the nanoflare hypothesis by constructing a magnetohydrodynamic (MHD) model that can derive the energy of a nanoflare from the nature of an ideal kink instability. The set of energy-releasing instabilities is captured by an instability threshold for linear kink modes. Each point on the threshold is associated with a unique energy release; thus we can predict a distribution of nanoflare energies. When the linear instability threshold is crossed, the instability enters a nonlinear phase as it is driven by current sheet reconnection. As the ensuing flare erupts and declines, the field transitions to a lower energy state, which is modelled by relaxation theory; i.e., helicity is conserved and the ratio of current to field becomes invariant within the loop. We apply the model so that all the loops within an ensemble achieve instability followed by energy-releasing relaxation. The result is a nanoflare energy distribution. Furthermore, we produce different distributions by varying the loop aspect ratio, the nature of the path to instability taken by each loop and also the level of radial expansion that may accompany loop relaxation. The heating rate obtained is just sufficient for coronal heating. In addition, we also show that kink instability cannot be associated with a critical magnetic twist value for every point along the instability threshold.

The solar influences data analysis centre

Abstract Since 1981, the Royal Observatory of Belgium has operated the Sunspot Index Data Centre, the World Data Centre for the Sunspot Index. Recently, the Space Weather Forecast Centre of Paris–Meudon was transferred and added to the activities of the SIDC. Moreover, a complete archive of all images of the SOHO instrument EIT has become available at the SIDC. Given all these extensions, the new style SIDC has become a ‘Solar Influences Data Centre’ that analyses solar activity and provides services on three different time scales: 1. Fast warnings and real time monitoring . As the Regional Warning Centre (RWC) for Western Europe of the International Space Environment Service (ISES), the SIDC collects and redistributes solar, geomagnetic, and ionospheric data in Western Europe. Short-term predictions (3 days) and alerts are produced on a daily basis. 2. Forecasts and middle term analysis . The SIDC takes care of the calculation of a sunspot index, called the International Sunspot Number. We compute and broadcast the daily, monthly, yearly international sunspot numbers, with middle range predictions (up to 12 months). 3. Post-event analysis and long-term solar cycle analysis . Since the launch of SOHO, EIT offers a global view of the EUV corona over the whole rising phase of the solar activity cycle. Such a long-duration data series is unprecedented and allows the study of the evolution over the solar cycle of objects classes such as active regions, coronal holes, coronal mass ejections or flares.

Recent Developments of NEMO: Detection of EUV Wave Characteristics

Recent developments in space instrumentation for solar observations and increased telemetry have necessitated the creation of advanced pattern recognition tools for different classes of solar events. The Extreme Ultraviolet Imaging Telescope (EIT) onboard the SOHO spacecraft has uncovered a new class of eruptive events on the solar disk, which are often identified as signatures of the initiation of coronal mass ejections (CMEs). The development of an automatic detection tool of these signatures is an important task. The Novel EIT Wave Machine Observing (NEMO) code (http://sidc.be/nemo) is an operational tool that automatically detects EUV waves using a sequence of EUV images. NEMO applies techniques based on the general statistical properties of the underlying physical mechanisms of eruptive events. Originally, the technique was applied to images taken with the EIT telescope. In this work, the most recent updates of the NEMO code are presented. These updates include calculations of the area of the dimming region, a novel clustering technique for the extraction of dimming regions, and new criteria to identify eruptive dimmings based on their complex characteristics. The efficiency of NEMO has been significantly increased and now permits the extraction of dimming regions observed near the solar limb and also the detection of small-scale events. Furthermore, the catalogs of solar eruptive events based on the updated NEMO may include a larger number of physical parameters associated with the dimming regions.

Helicity-injected current drive and open flux instabilities in spherical tokamaks

The toroidal current driven by relaxation processes in a cylindrical spherical tokamak (ST) geometry with coaxial injected flux is estimated by use of the linear ideal stability boundary of equilibria with a high current on the open driven flux and a lower current on the closed flux. Instabilities with toroidal mode number n = 1 have been shown to play a vital role in the helicity injection current drive, being closely associated with the relaxation process which distributes current from a directly driven open flux to a closed flux. Previous results for spheromaks (1D and 2D equilibria) and STs (1D equilibria) have predicted stabilization, for a given open flux current, if the closed flux plasma current is sufficiently large, suggesting that the current drive mechanism is self-limiting. New results presented here for 2D ST equilibria are consistent with the 1D results, but new features appear in the stability maps as the axial length and toroidal field (TF) strength are varied in the equilibria. These include changes in the shape of stability boundaries and the estimated driven current, changes in the mode structure due to equilibrium changes and resonance effects which extend stability boundaries into the stable region. As the minimum and maximum of the safety factor q profile cross integer rational values, the resonant mode is destabilized, causing regions of enhanced instability in the current profile parameter space. The results show the effects on the stability of varying the geometric length ratio R/L and provide driven current estimates with varying imposed TF strengths; these results have implications both for existing STs and for the design of future devices.