The tearing mode instability is one important mechanism that may explain the triggering of fast magnetic reconnection in astrophysical plasmas such as the solar corona. In this paper, the linear stability analysis of the tearing mode is carried out for a current sheet in the presence of a guide field, and including the Hall effect. We show that the presence of a strong guide field does not modify the most unstable mode in the two-dimensional wave vector space orthogonal to the current gradient direction, which remains the fastest growing parallel mode. With the Hall effect, the inclusion of a guide field turns the non-dispersive propagation along the guide field direction to a dispersive one. The oblique modes have a wave-like structure along the normal direction of the current sheet and a strong guide field suppresses this structure while making the eigen-functions asymmetric.

In many astrophysical environments the plasma is only partially ionized, and therefore the interaction of charged and neutral particles may alter both the triggering of reconnection and its subsequent dynamical evolution. We derive the tearing mode maximum growth rate for partially ionized plasmas in the cases of weak and strong coupling between the plasma and the neutrals. In addition, critical scalings for current sheet aspect ratios are presented in terms of Lundquist number and ion-neutral collision frequencies. In the decoupled regime the standard tearing mode is recovered with a small correction depending on the ion-neutral collision frequency; in the intermediate regime collisions with neutrals are shown to stabilize current sheets, resulting in larger critical aspect ratios for ideal tearing to occur. Nonetheless, the additional electron-neutral collisions, hidden in the definition of the Lundquist number, can shrink the critical aspect ratios below the fully ionized case. In the coupled regime, the growth rate depends on the density ratio between ions and neutrals through the collision frequency between these two species. These provide critical aspect ratios for which the tearing mode instability transitions from slow to ideal, that depend on the neutral-ion density ratio.

The slow solar wind is typically characterized as having low Alfv\'enicity. However, Parker Solar Probe (PSP) observed predominately Alfv\'enic slow solar wind during several of its initial encounters. From its first encounter observations, about 55.3\% of the slow solar wind inside 0.25 au is highly Alfv\'enic ($|\sigma_C| > 0.7$) at current solar minimum, which is much higher than the fraction of quiet-Sun-associated highly Alfv\'enic slow wind observed at solar maximum at 1 au. Intervals of slow solar wind with different Alfv\'enicities seem to show similar plasma characteristics and temperature anisotropy distributions. Some low Alfv\'enicity slow wind intervals even show high temperature anisotropies, because the slow wind may experience perpendicular heating as fast wind does when close to the Sun. This signature is confirmed by Wind spacecraft measurements as we track PSP observations to 1 au. Further, with nearly 15 years of Wind measurements, we find that the distributions of plasma characteristics, temperature anisotropy and helium abundance ratio ($N_\alpha/N_p$) are similar in slow winds with different Alfv\'enicities, but the distributions are different from those in the fast solar wind. Highly Alfv\'enic slow solar wind contains both helium-rich ($N_\alpha/N_p\sim0.045$) and helium-poor ($N_\alpha/N_p\sim0.015$) populations, implying it may originate from multiple source regions. These results suggest that highly Alfv\'enic slow solar wind shares similar temperature anisotropy and helium abundance properties with regular slow solar winds, and they thus should have multiple origins.

In contrast with the fast solar wind, that originates in coronal holes, the source of the slow solar wind is still debated. Often intermittent and enriched with low FIP elements -- akin to what is observed in closed coronal loops -- the slow wind could form in bursty events nearby helmet streamers. Slow winds also exhibit density perturbations which have been shown to be periodic and could be associated with flux ropes ejected from the tip of helmet streamers, as shown recently by the WISPR white light imager onboard Parker Solar Probe (PSP). In this work, we propose that the main mechanism controlling the release of flux ropes is a flow-modified tearing mode at the heliospheric current sheet (HCS). We use MHD simulations of the solar wind and corona to reproduce realistic configurations and outflows surrounding the HCS. We find that this process is able to explain long ($\sim 10-20$h) and short ($\sim 1-2$h) timescales of density structures observed in the slow solar wind. This study also sheds new light on the structure, topology and composition of the slow solar wind, and could be, in the near future, compared with white light and in situ PSP observations.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

One of the discoveries made by Parker Solar Probe during first encounters with the Sun is the ubiquitous presence of relatively small-scale structures standing out as sudden deflections of the magnetic field. They were called switchbacks as some of them show up the full reversal of the radial component of the magnetic field and then return to regular conditions. Analyzing the magnetic field and plasma perturbations associated with switchbacks we identify three types of structures with slightly different characteristics: 1. Alfvenic structures, where the variations of the magnetic field components take place while the magnitude of the field remains constant; 2. Compressional, the field magnitude varies together with changes of the components; 3. Structures manifesting full reversal of the magnetic field (extremal class of Alfvenic structures). Processing of structures boundaries and plasma bulk velocity perturbations lead to the conclusion that they represent localized magnetic field tubes with enhanced parallel plasma velocity and ion beta moving together with respect to surrounding plasma. The magnetic field deflections before and after the switchbacks reveal the existence of total axial current. The electric currents are concentrated on the relatively narrow boundary layers on the surface of the tubes and determine the magnetic field perturbations inside the tube. These currents are closed on the structure surface, and typically have comparable azimuthal and the axial components. The surface of the structure may also accommodate an electromagnetic wave, that assists to particles in carrying currents. We suggest that the two types of structures we analyzed here may represent the local manifestations of the tube deformations corresponding to a saturated stage of the Firehose instability development.

Magnetic reconnection in current sheets converts magnetic energy into particle energy. The process may play an important role in the acceleration and heating of the solar wind close to the Sun. Observations from Parker Solar Probe provide a new opportunity to study this problem, as it measures the solar wind at unprecedented close distances to the Sun. During the 1st orbit, PSP encountered a large number of current sheets in the solar wind through perihelion at 35.7 solar radii. We performed a comprehensive survey of these current sheets and found evidence for 21 reconnection exhausts. These exhausts were observed in heliospheric current sheets, coronal mass ejections, and regular solar wind. However, we find that the majority of current sheets encountered around perihelion, where the magnetic field was strongest and plasma beta was lowest, were Alfv\'enic structures associated with bursty radial jets and these current sheets did not appear to be undergoing local reconnection. We examined conditions around current sheets to address why some current sheets reconnected, while others did not. A key difference appears to be the degree of plasma velocity shear across the current sheets: The median velocity shear for the 21 reconnection exhausts was 24% of the Alfv\'en velocity shear, whereas the median shear across 43 Alfv\'enic current sheets examined was 71% of the Alfv\'en velocity shear. This finding could suggest that large, albeit sub-Alfv\'enic, velocity shears suppress reconnection. An alternative interpretation is that the Alfv\'enic current sheets are isolated rotational discontinuities which do not undergo local reconnection.

Observations at 1 au have confirmed that enhancements in measured energetic particle fluxes are statistically associated with "rough" magnetic fields, i.e., fields having atypically large spatial derivatives or increments, as measured by the Partial Variance of Increments (PVI) method. One way to interpret this observation is as an association of the energetic particles with trapping or channeling within magnetic flux tubes, possibly near their boundaries. However, it remains unclear whether this association is a transport or local effect; i.e., the particles might have been energized at a distant location, perhaps by shocks or reconnection, or they might experience local energization or re-acceleration. The Parker Solar Probe (PSP), even in its first two orbits, offers a unique opportunity to study this statistical correlation closer to the corona. As a first step, we analyze the separate correlation properties of the energetic particles measured by the \isois instruments during the first solar encounter. The distribution of time intervals between a specific type of event, i.e., the waiting time, can indicate the nature of the underlying process. We find that the \isois observations show a power-law distribution of waiting times, indicating a correlated (non-Poisson) distribution. Analysis of low-energy \isois data suggests that the results are consistent with the 1 au studies, although we find hints of some unexpected behavior. A more complete understanding of these statistical distributions will provide valuable insights into the origin and propagation of solar energetic particles, a picture that should become clear with future PSP orbits.

Direct evidence of an inertial-range turbulent energy cascade has been provided by spacecraft observations in heliospheric plasmas. In the solar wind, the average value of the derived heating rate near 1 au is $\sim 10^{3}\, \mathrm{J\,kg^{-1}\,s^{-1}}$, an amount sufficient to account for observed departures from adiabatic expansion. Parker Solar Probe (PSP), even during its first solar encounter, offers the first opportunity to compute, in a similar fashion, a fluid-scale energy decay rate, much closer to the solar corona than any prior in-situ observations. Using the Politano-Pouquet third-order law and the von K\'arm\'an decay law, we estimate the fluid-range energy transfer rate in the inner heliosphere, at heliocentric distance $R$ ranging from $54\,R_{\odot}$ (0.25 au) to $36\,R_{\odot}$ (0.17 au). The energy transfer rate obtained near the first perihelion is about 100 times higher than the average value at 1 au. This dramatic increase in the heating rate is unprecedented in previous solar wind observations, including those from Helios, and the values are close to those obtained in the shocked plasma inside the terrestrial magnetosheath.

The solar wind shows periods of highly Alfv\'enic activity, where velocity fluctuations and magnetic fluctuations are aligned or anti-aligned with each other. It is generally agreed that solar wind plasma velocity and magnetic field fluctuations observed by Parker Solar Probe (PSP) during the first encounter are mostly highly Alfv\'enic. However, quantitative measures of Alfv\'enicity are needed to understand how the characterization of these fluctuations compares with standard measures from prior missions in the inner and outer heliosphere, in fast wind and slow wind, and at high and low latitudes. To investigate this issue, we employ several measures to quantify the extent of Alfv\'enicity -- the Alfv\'en ratio $r_A$, {normalized} cross helicity $\sigma_c$, {normalized} residual energy $\sigma_r$, and the cosine of angle between velocity and magnetic fluctuations $\cos\theta_{vb}$. We show that despite the overall impression that the Alfv\'enicity is large in the solar wind sampled by PSP during the first encounter, during some intervals the cross helicity starts decreasing at very large scales. These length-scales (often $> 1000 d_i$) are well inside inertial range, and therefore, the suppression of cross helicity at these scales cannot be attributed to kinetic physics. This drop at large scales could potentially be explained by large-scale shears present in the inner heliosphere sampled by PSP. In some cases, despite the cross helicity being constant down to the noise floor, the residual energy decreases with scale in the inertial range. These results suggest that it is important to consider all these measures to quantify Alfv\'enicity.

During Parker Solar Probe's first orbit, the solar wind plasma has been observed in situ closer than ever before, the perihelion on November 6th 2018 revealing a flow that is constantly permeated by large amplitude Alfv\'enic fluctuations. These include radial magnetic field reversals, or switchbacks, that seem to be a persistent feature of the young solar wind. The measurements also reveal a very strong, unexpected, azimuthal velocity component. In this work, we numerically model the solar corona during this first encounter, solving the MHD equations and accounting for Alfv\'en wave transport and dissipation. We find that the large scale plasma parameters are well reproduced, allowing the computation of the solar wind sources at Probe with confidence. We try to understand the dynamical nature of the solar wind to explain both the amplitude of the observed radial magnetic field and of the azimuthal velocities.

We carry out two-dimensional magnetohydrodynamic (MHD) simulations of an ensemble of Alfv\'enic fluctuations propagating in a structured, expanding solar wind including the presence of fast and slow solar wind streams. Using an appropriate expanding box model, the simulations incorporate the effects of fast-slow stream shear and compression and rarefaction self-consistently. We investigate the radial and longitudinal evolution of the cross-helicity, the total and residual energies and the power spectra of outward and inward Alfv\'enic fluctuations. The stream interaction is found to strongly affect the radial evolution of Alfv\'enic turbulence. The total energy in the Alfv\'en waves is depleted within the velocity shear regions, accompanied by the decrease of the normalized cross-helicity. The presence of stream-compression facilitates this process. Residual energy fluctuates around zero due to the correlation and de-correlation between the inward/outward waves but no net growth or decrease of the residual energy is observed. The radial power spectra of the inward/outward Alfv\'en waves show significant longitudinal variations. Kolmogorov-like spectra are developed only inside the fast and slow streams and when both the compression and shear are present. On the other hand, the spectra along the longitudinal direction show clear Kolmogorov-like inertial ranges in all cases.

During the Parker Solar Probe's (PSP) first perihelion pass, the spacecraft reached within a heliocentric distance of \(\sim 37~R_\odot\) and observed numerous magnetic and flow structures characterized by sharp gradients. To better understand these intermittent structures in the young solar wind, an important property to examine is their degree of correlation in time and space. To this end, we use the well-tested Partial Variance of Increments (PVI) technique to identify intermittent events in FIELDS and SWEAP observations of magnetic and proton-velocity fields (respectively) during PSP's first solar encounter, when the spacecraft was within 0.25 au from the Sun. We then examine distributions of waiting times between events with varying separation and PVI thresholds. We find power-law distributions for waiting times shorter than a characteristic scale comparable to the correlation time, suggesting a high degree of correlation that may originate in a clustering process. Waiting times longer than this characteristic time are better described by an exponential, suggesting a random memory-less Poisson process at play. These findings are consistent with near-Earth observations of solar wind turbulence. The present study complements the one by Dudok de Wit et al. (2020, present volume), which focuses on waiting times between observed "switchbacks" in the radial magnetic field.

Parker Solar Probe (PSP) observations during its first encounter at 35.7 $R_\odot$ have shown the presence of magnetic field lines which are strongly perturbed to the point that they produce local inversions of the radial magnetic field, known as switchbacks. Their counterparts in the solar wind velocity field are local enhancements in the radial speed, or jets, displaying (in all components) the velocity-magnetic field correlation typical of large amplitude Alfv\'en waves propagating away from the Sun. Switchbacks and radial jets have previously been observed over a wide range of heliocentric distances by Helios, WIND and Ulysses, although they were prevalent in significantly faster streams than seen at PSP. Here we study via numerical MHD simulations the evolution of such large amplitude Alfv\'enic fluctuations by including, in agreement with observations, both a radial magnetic field inversion and an initially constant total magnetic pressure. Despite the extremely large excursion of magnetic and velocity fields, switchbacks are seen to persist for up to hundreds of Alfv\'en crossing times before eventually decaying due to the parametric decay instability. Our results suggest that such switchback/jet configurations might indeed originate in the lower corona and survive out to PSP distances, provided the background solar wind is sufficiently calm, in the sense of not being pervaded by strong density fluctuations or other gradients, such as stream or magnetic field shears, that might destabilize or destroy them over shorter timescales.

One of the most striking observations made by Parker Solar Probe during its first solar encounter is the omnipresence of rapid polarity reversals in a magnetic field that is otherwise mostly radial. These so-called switchbacks strongly affect the dynamics of the magnetic field. We concentrate here on their macroscopic properties. First, we find that these structures are self-similar, and have neither a characteristic magnitude, nor a characteristic duration. Their waiting time statistics shows evidence for aggregation. The associated long memory resides in their occurrence rate, and is not inherent to the background fluctuations. Interestingly, the spectral properties of inertial range turbulence differ inside and outside of switchback structures; in the latter the $1/f$ range extends to higher frequencies. These results suggest that outside of these structures we are in the presence of lower amplitude fluctuations with a shorter turbulent inertial range. We conjecture that these correspond to a pristine solar wind.

We compare magnetic field measurements taken by the FIELDS instrument on Parker Solar Probe (PSP) during its first solar encounter to predictions obtained by Potential Field Source Surface (PFSS) modeling. Ballistic propagation is used to connect the spacecraft to the source surface. Despite the simplicity of the model, our results show striking agreement with PSPs first observations of the heliospheric magnetic field from 0.5 AU (107.5 Rs) down to 0.16 AU (35.7 Rs). Further, we show the robustness of the agreement is improved both by allowing the photospheric input to the model to vary in time, and by advecting the field from PSP down to the PFSS model domain using in situ PSP/SWEAP measurements of the solar wind speed instead of assuming it to be constant with longitude and latitude. We also explore the source surface height parameter (RSS) to the PFSS model finding that an extraordinarily low source surface height (1.3-1.5Rs) predicts observed small scale polarity inversions which are otherwise washed out with regular modeling parameters. Finally, we extract field line traces from these models. By overlaying these on EUV images we observe magnetic connectivity to various equatorial and mid-latitude coronal holes indicating plausible magnetic footpoints and offering context for future discussions of sources of the solar wind measured by PSP.

We use 2.5D Magnetohydrodynamic simulations to investigate the spectral signatures of the nonlinear disruption of a tearing unstable current sheet via the generation of multiple secondary current sheets and magnetic islands. During the nonlinear phase of tearing mode evolution, there develops a regime in which the magnetic energy density shows a spectrum with a power-law close to $B(k)^2\sim k^{-0.8}$. Such an energy spectrum is found in correspondence of the neutral line, within the diffusion region of the primary current sheet, where energy is conveyed towards smaller scales via a ``recursive'' process of fast tearing-type instabilities. Far from the neutral line we find that magnetic energy spectra evolve towards slopes compatible with the ``standard'' Kolmogorov spectrum. Starting from a self-similar description of the nonlinear stage at the neutral line, we provide a model that predicts a reconnecting magnetic field energy spectrum scaling as $k^{-4/5}$, in good agreement with numerical results. An extension of the predicted power-law to generic current sheet profiles is also given and possible implications for turbulence phenomenology are discussed. These results provide a step forward to understand the ``recursive'' generation of magnetic islands (plasmoids), which has been proposed as a possible explanation for the energy release during flares, but which, more in general, can have an impact on the subsequent turbulent evolution of unstable sheets that naturally form in the high-Lundquist number and collisionless plasmas found in most of the astrophysical environments.

Metis is the first solar coronagraph designed for a space mission capable of performing simultaneous imaging of the off-limb solar corona in both visible and UV light. The observations obtained with Metis aboard the Solar Orbiter ESA-NASA observatory will enable us to diagnose, with unprecedented temporal coverage and spatial resolution, the structures and dynamics of the full corona from 1.7 $R_\odot$ to about 9 $R_\odot$. Due to the uniqueness of the Solar Orbiter mission profile, Metis will be able to observe the solar corona from a close vantage point (down to 0.28 AU), achieving out-of-ecliptic views with the increase of the orbit inclination over time. Moreover, observations near perihelion, during the phase of lower rotational velocity of the solar surface relative to the spacecraft, will allow longer-term studies of the coronal features. Thanks to a novel occultation design and a combination of a UV interference coating of the mirrors and a spectral bandpass filter, Metis images the solar corona simultaneously in the visible light band, between 580 and 640 nm, and in the UV H I Lyman-{\alpha} line at 121.6 nm. The coronal images in both the UV Lyman-{\alpha} and polarised visible light are obtained at high spatial resolution with a spatial scale down to about 2000 km and 15000 km at perihelion, in the cases of the visible and UV light, respectively. A temporal resolution down to 1 second can be achieved when observing coronal fluctuations in visible light. The Metis measurements will allow for complete characterisation of the main physical parameters and dynamics of the electron and neutral hydrogen/proton plasma components of the corona in the region where the solar wind undergoes acceleration and where the onset and initial propagation of coronal mass ejections take place, thus significantly improving our understanding of the region connecting the Sun to the heliosphere.

Magnetohydrodynamic (MHD) theory and simulations have shown that reconnection is triggered via a fast "ideal" tearing instability in current sheets whose inverse aspect ratio decreases to $a/L\sim S^{-1/3}$, with $S$ is the Lundquist number defined by the half-length $L$ of the current sheet (of thickness $2a$). Ideal tearing, in 2D sheets, triggers a hierarchical collapse via stretching of X-points and recursive instability. At each step, the local Lundquist number decreases, until the subsequent sheet thickness either approaches kinetic scales or the Lundquist number becomes sufficiently small. Here we carry out a series of Hall-MHD simulations to show how the Hall effect modifies recursive reconnection once the ion inertial scale is approached. We show that as the ion inertial length becomes of the order of the inner, singular layer thickness at some step of the recursive collapse, reconnection transits from the plasmoid-dominant regime to an intermediate plasmoid+Hall regime and then to the Hall-dominant regime. The structure around the X-point, the reconnection rate, the dissipation property and the power spectra are also modified significantly by the Hall effect.

We investigate the dynamics of a closed corona cartesian reduced magnetohydrodynamic (MHD) model where photospheric vortices twist the coronal magnetic field lines. We consider two corotating or counter-rotating vortices localized at the center of the photospheric plate, and additionally more corotating vortices that fill the plate entirely. Our investigation is specifically devoted to study the fully nonlinear stage, after the linear stage during which the vortices create laminar and smoothly twisting flux tubes. Our main goal is to understand the dynamics of photospheric vortices twisting the field lines of a coronal magnetic field permeated by finite amplitude broadband fluctuations. We find that depending on the arrangement and handedness of the photospheric vortices an inverse cascade storing a significant amount of magnetic energy may occur or not. In the first case a reservoir of magnetic energy available to large events such as destabilization of a pre-CME configuration develops, while in the second case the outcome is a turbulent heated corona. Although our geometry is simplified our simulations are shown to have relevant implications for coronal dynamics and CME initiation.

Magnetic reconnection may be the fundamental process allowing energy stored in magnetic fields to be released abruptly, solar flares and coronal mass ejection (CME) being archetypal natural plasma examples. Magnetic reconnection is much too slow a process to be efficient on the large scales, but accelerates once small enough scales are formed in the system. For this reason, the fractal reconnection scenario was introduced (Shibata and Tanuma 2001) to explain explosive events in the solar atmosphere: it was based on the recursive triggering and collapse via tearing instability of a current sheet originally thinned during the rise of a filament in the solar corona. Here we compare the different fractal reconnection scenarios that have been proposed, and derive generalized scaling relations for the recursive triggering of fast, `ideal' - i.e. Lundquist number independent - tearing in collapsing current sheet configurations with arbitrary current profile shapes. An important result is that the Sweet-Parker scaling with Lundquist number, if interpreted as the aspect ratio of the singular layer in an ideally unstable sheet, is universal and does not depend on the details of the current profile in the sheet. Such a scaling however must not be interpreted in terms of stationary reconnection, rather it defines a step in the accelerating sequence of events of the ideal tearing mediated fractal cascade. We calculate scalings for the expected number of plasmoids for such generic profiles and realistic Lundquist numbers.

This paper briefly reviews a number of fundamental measurements that need to be made in order to characterize turbulence in space plasmas such as the solar wind. It has long been known that many of these quantities require simultaneous multipoint measurements to attain a proper characterization that would reveal the fundamental physics of plasma turbulence. The solar wind is an ideal plasma for such an investigation, and it now appears to be technologically feasible to carry out such an investigation, following the pioneering Cluster and MMS missions. Quantities that need to be measured using multipoint measurements include the two-point, two-time second correlation function of velocity, magnetic field and density, and higher order statistical objects such as third and fourth order structure functions. Some details of these requirements are given here, with a eye towards achieving closure on fundamental questions regarding the cascade rate, spectral anisotropy, characteristic coherent structures, intermittency, and dissipation mechanisms that describe plasma turbuelence, as well as its variability with plasma parameters in the solar wind. The motivation for this discussion is the current planning for a proposed Helioswarm mission that would be designed to make these measurements,leading to breakthrough understanding of the physics of space and astrophysical turbulence.

Ion fractional charge states, measured in situ in the heliosphere, depend on the properties of the plasma in the inner corona. As the ions travel outward in the solar wind and the electron density drops, the charge states remain essentially unaltered or "frozen in". Thus they can provide a powerful constraint on heating models of the corona and acceleration of the solar wind. We have implemented non-equilibrium ionization calculations into a 1D wave-turbulence-driven (WTD) hydrodynamic solar wind model and compared modeled charge states with the Ulysses 1994-5 in situ measurements. We have found that modeled charge state ratios of $C^{6+}/C^{5+}$ and $O^{7+}/O^{6+}$, among others, were too low compared with Ulysses measurements. However, a heuristic reduction of the plasma flow speed has been able to bring the modeled results in line with observations, though other ideas have been proposed to address this discrepancy. We discuss implications of our results and the prospect of including ion charge state calculations into our 3D MHD model of the inner heliosphere.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

We explore new opportunities for solar physics that could be realized by future missions providing sustained observations from vantage points away from the Sun-Earth line. These include observations from the far side of the Sun, at high latitudes including over the solar poles, or from near-quadrature angles relative to the Earth (e.g., the Sun-Earth L4 and L5 Lagrangian points). Such observations fill known holes in our scientific understanding of the three-dimensional, time-evolving Sun and heliosphere, and have the potential to open new frontiers through discoveries enabled by novel viewpoints.

Our knowledge of the Sun, its atmosphere, long term activity and space weather potential is severely limited by the lack of good observations of the polar and far-side regions. Observations from a polar vantage point would revolutionize our understanding of the mechanism of solar activity cycles, polar magnetic field reversals, the internal structure and dynamics of the Sun and its atmosphere. Only with extended (many day) observations of the polar regions can the polar flows be determined down to the tachocline where the dynamo is thought to originate. Rapid short period polar orbits, using in situ and remote sensing instrumentation, distributed over a small number of spacecraft, will provide continuous 360o coverage of the solar surface and atmosphere in both longitude and latitude for years on end. This unprecedented full coverage will enable breakthrough studies of the physical connection between the solar interior, the solar atmosphere, the solar wind, solar energetic particles and the inner heliosphere at large. A potential implementation, the Solar Polar Diamond Explorer (SPDEx) built upon the Solar Polar Imager mission design, involves up to four small spacecraft in a 0.48-AU orbit with an inclination of 75o. The orbit is achieved using solar sails or ion engines, both technologies already demonstrated in space.

The heating of coronal loops is investigated to understand the observational consequences in terms of the thermodynamics and radiative losses from the Sun as well as the magnetized coronae of stars with an outer convective envelope. The dynamics of the Parker coronal heating model are studied for different ratios of the photospheric forcing velocity timescale $t_p$ to the Alfv\'en crossing time along a loop $t_A$. It is shown that for $t_p/t_A \gtrsim$ 10--24 the heating rate and maximum temperature are largest and approximately independent of $t_p/t_A$, leading to a strong emission in X-rays and EUV. On the opposite decreasing $t_p/t_A$ to smaller values leads to lower heating rates and plasma temperatures, and consequently fading high-energy radiative emission once $t_p/t_A \lesssim$ 1--3. The average volumetric loop heating rate is shown to scale as $\ell_p u_p B_0^2/4\pi L^2$, where $\ell_p$ and $u_p$ are respectively the convective granule length-scale and velocity, $B_0$ is the intensity of the strong magnetic field threading the loop, and $L$ the loop length. These findings support a recent hypothesis explaining ultracool dwarf observations of stars with similar magnetic field strength but radically different topologies displaying different radiative emission.

Magnetohydrodynamic simulations have shown that a non-unique critical Lundquist number $S_c$ exists, hovering around $S_c \sim 10^4$, above which threshold Sweet-Parker type stationary reconnecting configurations become unstable to a fast tearing mode dominated by plasmoid generation. It is known that the flow along the sheet plays a stabilizing role, though a satisfactory explanation of the non-universality and variable critical Lundquist numbers observed is still lacking. Here we discuss this question using 2D linear MHD simulations and linear stability analyses of Sweet-Parker type current sheets in the presence of background stationary inflows and outflows at low Lundquist numbers ($S\le 10^4$). Simulations show that the inhomogeneous outflow stabilizes the current sheet by stretching the growing magnetic islands and at the same time evacuating the magnetic islands out of the current sheet. This limits the time during which fluctuations which begin at any given wave-length can remain unstable, rendering the instability non-exponential. We find that the linear theory based on the expanding-wavelength assumption works well for $S$ larger than $\sim 1000$. However we also find that the inflow and location of the initial perturbation also affect the stability threshold.

In this paper we study the scaling relations for the triggering of the fast, or "ideal", tearing instability starting from equilibrium configurations relevant to astrophysical as well as laboratory plasmas that differ from the simple Harris current sheet configuration. We present the linear tearing instability analysis for equilibrium magnetic fields which a) go to zero at the boundary of the domain and b) contain a double current sheet system (the latter previously studied as a cartesian proxy for the m=1 kink mode in cylindrical plasmas). More generally, we discuss the critical aspect ratio scalings at which the growth rates become independent of the Lundquist number $S$, in terms of the dependence of the $\Delta'$ parameter on the wavenumber $k$ of unstable modes. The scaling $\Delta'(k)$ with $k$ at small $k$ is found to categorize different equilibria broadly: the critical aspect ratios may be even smaller than $L/a \sim S^{\alpha}$ with $\alpha=1/3$ originally found for the Harris current sheet, but there exists a general lower bound $\alpha\ge1/4$.

The magnetic topology and field line random walk properties of a nanoflare-heated and magnetically confined corona are investigated in the reduced magnetohydrodynamic regime. Field lines originating from current sheets form coherent structures, called Current Sheet Connected (CSC) regions, extended around them. CSC field line random walk is strongly anisotropic, with preferential diffusion along the current sheets' in-plane length. CSC field line random walk properties remain similar to those of the entire ensemble but exhibit enhanced mean square displacements and separations due to the stronger magnetic field intensities in CSC regions. The implications for particle acceleration and heat transport in the solar corona and wind, and for solar moss formation are discussed.

Alfv\'enic fluctuations in the solar wind display many properties reflecting an ongoing nonlinear cascade, e.g. a well-defined spectrum in frequency, together with some characteristics more commonly associated with the linear propagation of waves from the Sun, such as the variation of fluctuation amplitude with distance, dominated by solar wind expansion effects. Therefore both nonlinearities and expansion must be included simultaneously in any successful model of solar wind turbulence evolution. Because of the disparate spatial scales involved, direct numerical simulations of turbulence in the solar wind represent an arduous task, especially if one wants to go beyond the incompressible approximation. Indeed, most simulations neglect solar wind expansion effects entirely. Here we develop a numerical model to simulate turbulent fluctuations from the outer corona to 1 AU and beyond, including the sub-Alfv\'enic corona. The accelerating expanding box (AEB) extends the validity of previous expanding box models by taking into account both the acceleration of the solar wind and the inhomogeneity of background density and magnetic field. Our method incorporates a background accelerating wind within a magnetic field that naturally follows the Parker spiral evolution using a two-scale analysis in which the macroscopic spatial effect coupling fluctuations with background gradients becomes a time-dependent coupling term in a homogeneous box. In this paper we describe the AEB model in detail and discuss its main properties, illustrating its validity by studying Alfv\'en wave propagation across the Alfv\'en critical point.

One of the main questions in magnetic reconnection is the origin of triggering behavior with on/off properties that accounts, once it is activated, for the fast magnetic energy conversion to kinetic and thermal energies at the heart of explosive events in astrophysical and laboratory plasmas. Over the past decade progress has been made on the initiation of fast reconnection via the plasmoid instability and what has been called "ideal" tearing, which sets in once current sheets thin to a critical inverse aspect ratio $(a/L)_c$: as shown by Pucci and Velli (2014), at $(a/L)_c \sim S^{-1/3}$ the time scale for the instability to develop becomes of the order of the Alfv\'en time and independent of the Lundquist number (here defined in terms of current sheet length $L$). However, given the large values of $S$ in natural plasmas, this transition might occur for thicknesses of the inner resistive singular layer which are comparable to the ion inertial length $d_i$. When this occurs, Hall currents produce a three-dimensional quadrupole structure of magnetic field, and the dispersive waves introduced by the Hall effect accelerate the instability. Here we present a linear study showing how the "ideal" tearing mode critical aspect ratio is modified when Hall effects are taken into account, including more general scaling laws of the growth rates in terms of sheet inverse aspect ratio: the critical inverse aspect ratio is amended to $a/L \simeq (di/L)^ {0.29} (1/S)^{0.19}$, at which point the instability growth rate becomes Alfv\'enic and does not depend on either of the (small) parameters $d_i/L, 1/S$. We discuss the implications of this generalized triggering aspect ratio for recently developed phase diagrams of magnetic reconnection.

To simulate the energy balance of coronal plasmas on macroscopic scales, we often require the specification of the coronal heating mechanism in some functional form. To go beyond empirical formulations and to build a more physically motivated heating function, we investigate the wave-turbulence-driven (WTD) phenomenology for the heating of closed coronal loops. Our implementation is designed to capture the large-scale propagation, reflection, and dissipation of wave turbulence along a loop. The parameter space of this model is explored by solving the coupled WTD and hydrodynamic evolution in 1D for an idealized loop. The relevance to a range of solar conditions is also established by computing solutions for over one hundred loops extracted from a realistic 3D coronal field. Due to the implicit dependence of the WTD heating model on loop geometry and plasma properties along the loop and at the footpoints, we find that this model can significantly reduce the number of free parameters when compared to traditional empirical heating models, and still robustly describe a broad range of quiet-sun and active region conditions. The importance of the self-reflection term in producing relatively short heating scale heights and thermal nonequilibrium cycles is also discussed.

Magnetic reconnection is thought to be the dynamical mechanism underlying many explosive phenomena observed both in space and in the laboratory, though the question of how fast magnetic reconnection is triggered in such high Lundquist ($S$) number plasmas has remained elusive. It has been well established that reconnection can develop over timescales faster than those predicted traditionally once kinetic scales are reached. It has also been shown that, within the framework of resistive Magnetohydrodynamics (MHD), fast reconnection is achieved for thin enough sheets via the onset of the so-called plasmoid instability. The latter was discovered in studies specifically devoted to the Sweet-Parker current sheet, either as an initial condition or an apparent transient state developing in nonlinear studies. On the other hand, a fast tearing instability can grow on an ideal, i.e., $S$-independent, timescale (dubbed "ideal" tearing) within current sheets whose aspect ratio scales with the macroscopic Lundquist number as $L/a\sim S^{1/3}$ -- much smaller than the Sweet-Parker one -- suggesting a new way to approach to the initiation of fast reconnection in collapsing current configurations. Here we present an overview of what we have called "ideal" tearing in resistive MHD, and discuss how the same reasoning can be extended to other plasma models commonly used that include electron inertia and kinetic effects. We then discuss a scenario for the onset of "ideal" fast reconnection via collapsing current sheets and describe a quantitative model for the interpretation of the nonlinear evolution of "ideally" unstable sheets in two dimensions.

Classical MHD reconnection theories, both the stationary Sweet-Parker model and the tearing instability, are known to provide rates which are too slow to explain the observations. However, a recent analysis has shown that there exists a critical threshold on current sheet's thickness, namely a/L~S^(-1/3), beyond which the tearing modes evolve on fast macroscopic Alfvenic timescales, provided the Lunquist number S is high enough, as invariably found in solar and astrophysical plasmas. Therefore, the classical Sweet-Parker scenario, for which the diffusive region scales as a/L~S^(-1/2) and thus can be up to ~100 times thinner than the critical value, is likely to be never realized in nature, as the current sheet itself disrupts in the elongation process. We present here two-dimensional, compressible, resistive MHD simulations, with S ranging from 10^5 to 10^7, that fully confirm the linear analysis. Moreover, we show that a secondary plasmoid instability always occurs when the same critical scaling is reached on the local, smaller scale, leading to a cascading explosive process, reminiscent of the flaring activity.

The evolution of a coronal loop is studied by means of numerical simulations of the fully compressible three-dimensional magnetohydrodynamic equations using the HYPERION code. The footpoints of the loop magnetic field are advected by random motions. As a consequence the magnetic field in the loop is energized and develops turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets: energy is deposited at small scales where heating occurs. Dissipation is non-uniformly distributed so that only a fraction of the coronal mass and volume gets heated at any time. Temperature and density are highly structured at scales which, in the solar corona, remain observationally unresolved: the plasma of our simulated loop is multi-thermal, where highly dynamical hotter and cooler plasma strands are scattered throughout the loop at sub-observational scales. Numerical simulations of coronal loops of 50000 km length and axial magnetic field intensities ranging from 0.01 to 0.04 Tesla are presented. To connect these simulations to observations we use the computed number densities and temperatures to synthesize the intensities expected in emission lines typically observed with the Extreme ultraviolet Imaging Spectrometer (EIS) on Hinode. These intensities are used to compute differential emission measure distributions using the Monte Carlo Markov Chain code, which are very similar to those derived from observations of solar active regions. We conclude that coronal heating is found to be strongly intermittent in space and time, with only small portions of the coronal loop being heated: in fact, at any given time, most of the corona is cooling down.

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

We study the linear and nonlinear evolution of the tearing instability on thin current sheets by means of two-dimensional numerical simulations, within the framework of compressible, resistive magnetohydrodynamics. In particular we analyze the behavior of current sheets whose inverse aspect ratio scales with the Lundquist number $S$ as $S^{-1/3}$. This scaling has been recently recognized to yield the threshold separating fast, ideal reconnection, with an evolution and growth which are independent of $S$ provided this is high enough, as it should be natural having the ideal case as a limit for $S\to\infty$. Our simulations confirm that the tearing instability growth rate can be as fast as $\gamma\approx 0.6\,{\tau_A}^{-1}$, where $\tau_A$ is the ideal Alfv\'enic time set by the macroscopic scales, for our least diffusive case with $S=10^7$. The expected instability dispersion relation and eigenmodes are also retrieved in the linear regime, for the values of $S$ explored here. Moreover, in the nonlinear stage of the simulations we observe secondary events obeying the same critical scaling with $S$, here calculated on the \emph{local}, much smaller lengths, leading to increasingly faster reconnection. These findings strongly support the idea that in a fully dynamic regime, as soon as current sheets develop, thin and reach this critical threshold in their aspect ratio, the tearing mode is able to trigger plasmoid formation and reconnection on the local (ideal) Alfv\'enic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas.

We investigate properties of the plasma fluid motion in the large amplitude low frequency fluctuations of highly Alfv\'enic fast solar wind. We show that protons locally conserve total kinetic energy when observed from an effective frame of reference comoving with the fluctuations. For typical properties of the fast wind, this frame can be reasonably identified by alpha particles, which, owing to their drift with respect to protons at about the Alfv\'en speed along the magnetic field, do not partake in the fluid low frequency fluctuations. Using their velocity to transform proton velocity into the frame of Alfv\'enic turbulence, we demonstrate that the resulting plasma motion is characterized by a constant absolute value of the velocity, zero electric fields, and aligned velocity and magnetic field vectors as expected for unidirectional Alfv\'enic fluctuations in equilibrium. We propose that this constraint, via the correlation between velocity and magnetic field in Alfv\'enic turbulence, is at the origin of the observed constancy of the magnetic field: while the constant velocity corresponding to constant energy can be only observed in the frame of the fluctuations, the correspondingly constant total magnetic field, invariant for Galilean transformations, remains the observational signature, in the spacecraft frame, of the constant total energy in the Alfv\'en turbulence frame.

Although it is widely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process, the detailed mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still not fully understood. Cranmer et al. (2007) developed a sophisticated, 1D, time-steady model of the solar wind with turbulence dissipation. By varying the coronal magnetic field, they obtain, for a single choice of wave properties, a realistic range of slow and fast wind conditions with a sharp latitudinal transition between the two streams. Using a 1D, time-dependent model of the solar wind of Lionello et al. (2014), which incorporates turbulent dissipation of Alfv\'en waves to provide heating and acceleration of the plasma, we have explored a similar configuration, obtaining qualitatively equivalent results. However, our calculations suggest that the rapid transition between slow and fast wind suggested by this 1D model may be disrupted in multidimensional MHD simulations by the requirement of transverse force balance.

The long-term evolution of large-amplitude Alfven waves propagating in the solar wind is investigated by performing two-dimensional MHD simulations within the expanding box model. The linear and nonlinear phases of the parametric decay instability are studied for both circularly polarized waves in parallel propagation and for arc-polarized waves in oblique propagation. The non-monochromatic case is also considered. In the oblique case, the direct excitation of daughter modes transverse to the local background field is found for the first time in an expanding environment, and this transverse cascade seems to be favored for monochromatic mother waves. The expansion effect reduces the instability growth rate, and it can even suppress its onset for the lowest frequency modes considered here, possibly explaining the persistence of these outgoing waves in the solar wind.

Although the mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still being actively investigated, it is largely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process. \citet{2010ApJ...708L.116V} presented a model for heating and accelerating the solar wind based on the turbulent dissipation of Alfv\'en waves. We first use a time-dependent model of the solar wind to reproduce one of \citeauthor{2010ApJ...708L.116V}'s solutions; then we extend its application to the case when the energy equation includes thermal conduction and radiation losses, and the upper chromosphere is part of the computational domain. Using this model, we explore parameter space and describe the characteristics of a fast-solar-wind solution. We discuss how this formulation may be applied to a 3D MHD model of the coron a and solar wind \citep{2009ApJ...690..902L}.

Recent high-resolution observations from the Solar Dynamics Observatory (SDO) have reawakened interest in the old and fascinating phenomenon of solar tornado-like prominences. This class of prominences was first introduced by Pettit (1932), who studied them over many years. Observations of tornado prominences similar to the ones seen by SDO had already been documented by Secchi (1877) in his famous "Le Soleil". High resolution and high cadence multiwavelength data obtained by SDO reveal that the tornado-like appearance of these prominences is mainly an illusion due to projection effects. We discuss two different cases where prominences on the limb might appear to have a tornado-like behavior. One case of apparent vortical motions in prominence spines and barbs arises from the (mostly) 2D counterstreaming plasma motion along the prominence spine and barbs together with oscillations along individual threads. The other case of apparent rotational motion is observed in prominence cavities and results from the 3D plasma motion along the writhed magnetic fields inside and along the prominence cavity as seen projected on the limb. Thus, the "tornado" impression results either from counterstreaming and oscillations or from the projection on the plane of the sky of plasma motion along magnetic field lines, rather than from a true vortical motion around an (apparent) vertical or horizontal axis. We discuss the link between tornado-like prominences, filament barbs, and photospheric vortices at their base.

We present simulations modeling closed regions of the solar corona threaded by a strong magnetic field where localized photospheric vortical motions twist the coronal field lines. The linear and nonlinear dynamics are investigated in the reduced magnetohydrodynamic regime in Cartesian geometry. Initially the magnetic field lines get twisted and the system becomes unstable to the internal kink mode, confirming and extending previous results. As typical in this kind of investigations, where initial conditions implement smooth fields and flux-tubes, we have neglected fluctuations and the fields are laminar until the instability sets in. But previous investigations indicate that fluctuations, excited by photospheric motions and coronal dynamics, are naturally present at all scales in the coronal fields. Thus, in order to understand the effect of a photospheric vortex on a more realistic corona, we continue the simulations after kink instability sets in, when turbulent fluctuations have already developed in the corona. In the nonlinear stage the system never returns to the simple initial state with ordered twisted field lines, and kink instability does not occur again. Nevertheless field lines get twisted, but in a disordered way, and energy accumulates at large scales through an inverse cascade. This energy can subsequently be released in micro-flares or larger flares, when interaction with neighboring structures occurs or via other mechanisms. The impact on coronal dynamics and CMEs initiation is discussed.

The dynamical response of the solar corona to surface and sub-surface perturbations depends on the chromospheric stratification, and specifically on how efficiently these layers reflect or transmit incoming Alfv\'en waves. While it would be desirable to include the chromospheric layers in the numerical simulations used to study such phenomena, that is most often not feasible. We defined and tested a simple approximation allowing the study of coronal phenomena while taking into account a parametrised chromospheric reflectivity. We addressed the problems of the transmission of the surface rotation to the corona and that of the generation of polar plumes by Alfv\'en waves (Pinto et al., 2010, 2011). We found that a high (yet partial) effective chromospheric reflectivity is required to properly describe the angular momentum balance in the corona and the way the surface differential rotation is transmitted upwards. Alfv\'en wave-driven polar plumes maintain their properties for a wide range of values for the reflectivity, but they become bursty (and eventually disrupt) when the limit of total reflection is attained.

We investigate the role of kinetic instabilities driven by a proton anisotropy on the onset of magnetic reconnection by means of 2-D hybrid simulations. The collisionless tearing of a current sheet is studied in the presence of a proton temperature anisotropy in the surrounding plasma. Our results confirm that anisotropic protons within the current sheet region can significantly enhance/stabilize the tearing instability of the current. Moreover, fluctuations associated to linear instabilities excited by large proton temperature anisotropies can significantly influence the stability of the plasma and perturb the current sheets, triggering the tearing instability. We find that such a complex coupling leads to a faster tearing evolution in a regime with larger perpendicular temperature when an ion-cyclotron instability is generated by the anisotropic proton distribution functions. On the contrary, in the presence of the opposite anisotropy, fire hose fluctuations excited by the unstable background protons with larger parallel temperature are not able to efficiently destabilize the current sheets, which remain stable for a long time after fire hose saturation. We discuss possible influences of this novel coupling on the solar wind and heliospheric plasma dynamics.

We demonstrate that major asymmetries in erupting filaments and CMEs, namely major twists and non-radial motions are typically related to the larger-scale ambient environment around eruptive events. Our analysis of prominence eruptions observed by the STEREO, SDO and SOHO spacecraft shows that prominence spines retain, during the initial phases, the thin ribbon-like topology they had prior to the eruption. This topology allows bending, rolling, and twisting during the early phase of the eruption, but not before. The combined ascent and initial bending of the filament ribbon is non-radial in the same general direction as for the enveloping CME. However, the non-radial motion of the filament is greater than that of the CME. In considering the global magnetic environment around CMEs, as approximated by the Potential Field Source Surface (PFSS) model, we find that the non-radial propagation of both erupting filaments and associated CMEs is correlated with the presence of nearby coronal holes, which deflect the erupting plasma and embedded fields. In addition, CME and filament motions respectively are guided towards weaker field regions, namely null points existing at different heights in the overlying configuration. Due to the presence of the coronal hole, the large-scale forces acting on the CME may be asymmetric. We find that the CME propagates usually non-radially in the direction of least resistance, which is always away from the coronal hole. We demonstrate these results using both low and high latitude examples.

We discuss observations of pseudostreamers and their 3D magnetic configuration as reconstructed with potential field source surface (PFSS) models to study their contribution to the solar wind. To understand the outflow from pseudostreamers the 3D expansion factor must be correctly estimated. Pseudostreamers may contain filament channels at their base in which case the open field lines diverge more strongly and the corresponding greater expansion factors lead to slower wind outflow, compared with pseudostreamers in which filament channels are absent. In the neighborhood of pseudostreamers the expansion factor does not increase monotonically with distance from the sun, and doesn't simply depend on the height of the pseudostreamer null point but on the entire magnetic field configuration.

Magnetic reconnection at the interface between coronal holes and loops, so-called interchange reconnection, can release the hotter, denser plasma from magnetically confined regions into the heliosphere, contributing to the formation of the highly variable slow solar wind. The interchange process is often thought to develop at the apex of streamers or pseudo-streamers, near Y and X-type neutral points, but slow streams with loop composition have been recently observed along fanlike open field lines adjacent to closed regions, far from the apex. However, coronal heating models, with magnetic field lines shuffled by convective motions, show that reconnection can occur continuously in unipolar magnetic field regions with no neutral points: photospheric motions induce a magnetohydrodynamic turbulent cascade in the coronal field that creates the necessary small scales, where a sheared magnetic field component orthogonal to the strong axial field is created locally and can reconnect. We propose that a similar mechanism operates near and around boundaries between open and closed regions inducing a continual stochastic rearrangement of connectivity. We examine a reduced magnetohydrodynamic model of a simplified interface region between open and closed corona threaded by a strong unipolar magnetic field. This boundary is not stationary, becomes fractal, and field lines change connectivity continuously, becoming alternatively open and closed. This model suggests that slow wind may originate everywhere along loop-coronal hole boundary regions, and can account naturally and simply for outflows at and adjacent to such boundaries and for the observed diffusion of slow wind around the heliospheric current sheet.

Context. Photospheric motions shuffle the footpoints of the strong axial magnetic field that threads coronal loops giving rise to turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets where energy is deposited at small-scales and the heating occurs. Previous studies show that current sheets thickness is orders of magnitude smaller than current state of the art observational resolution (~700 km). Aim. In order to understand coronal heating and interpret correctly observations it is crucial to study the thermodynamics of such a system where energy is deposited at unresolved small-scales. Methods. Fully compressible three-dimensional magnetohydrodynamic simulations are carried out to understand the thermodynamics of coronal heating in the magnetically confined solar corona. Results. We show that temperature is highly structured at scales below observational resolution and nonhomogeneously distributed so that only a fraction of the coronal mass and volume gets heated at each time. Conclusions. This is a multi-thermal system where hotter and cooler plasma strands are found one next to the other also at sub-resolution scales and exhibit a temporal dynamics.

Context: Supersonic nonthermal motions in molecular clouds are often interpreted as long-lived magnetohydrodynamic (MHD) waves. The propagation and amplitude of these waves is affected by local physical characteristics, most importantly the gas density and the ionization fraction. Aims: We study the propagation, reflection and dissipation of Alfv\'en waves in molecular clouds deriving the behavior of observable quantities such as the amplitudes of velocity fluctuations and the rate of energy dissipation. Methods: We formulated the problem in terms of Els\"asser variables for transverse MHD waves propagating in a one-dimensional inhomogeneous medium, including the dissipation due to collisions between ions and neutrals and to a nonlinear turbulent cascade treated in a phenomenological way. We considered both steady-state and time-dependent situations and solved the equations of the problem numerically with an iterative method and a Lax-Wendroff scheme, respectively. Results: Alfv\'en waves incident on overdense regions with density profiles typical of cloud cores embedded in a diffuse gas suffer enhanced reflection in the regions of the steepest density gradient, and strong dissipation in the core's interior. These effects are especially significant when the wavelength is intermediate between the critical wavelength for propagation and the typical scale of the density gradient. For larger wave amplitudes and/or steeper input spectra, the effects of the perpendicular turbulent cascade result in a stronger energy dissipation in the regions immediately surrounding the dense core. Conclusions: The results may help to interpret the sharp decrease of line width observed in the environments of low-mass cloud cores in several molecular transitions.

We present a mechanism for the formation of the low frequency 1/f magnetic spectrum based on numerical solutions of a shell reduced-MHD model of the turbulent dynamics inside the sub-Alfv\'enic solar wind. We assign reasonably realistic profiles to the wind speed and the density along the radial direction, and a radial magnetic field. Alfv\'en waves of short periodicity (600 s) are injected at the base of the chromosphere, penetrate into the corona and are partially reflected, thus triggering a turbulent cascade. The cascade is strong for the reflected wave while it is weak for the outward propagating waves. Reflection at the transition region recycles the strong turbulent spectrum into the outward weak spectrum, which is advected beyond the Alfv\'enic critical point without substantial evolution. There, the magnetic field has a perpendicular power-law spectrum with slope close to the Kolmogorov -5/3. The parallel spectrum is inherited from the frequency spectrum of large (perpendicular) eddies. The shape is a double power-law with slopes of -1 and -2 at low and high frequencies respectively, the position of the break depending on the injected spectrum. We suggest that the double power-law spectrum measured by Helios at 0.3 AU, where the average magnetic field is not aligned with the radial (contrary to our assumptions) results from the combination of such different spectral slopes. At low frequency the parallel spectrum dominates with its characteristic 1/f shape, while at higher frequencies its steep spectral slope (-2) is masked by the more energetic perpendicular spectrum (slope -5/3).

Coronal loops act as resonant cavities for low frequency fluctuations that are transmitted from the deeper layers of the solar atmosphere and are amplified in the corona, triggering nonlinear interactions. However trapping is not perfect, some energy leaks down to the chromosphere, thus limiting the turbulence development and the associated heating. We consider the combined effects of turbulence and leakage in determining the energy level and associated heating rate in models of coronal loops which include the chromosphere and transition region. We use a piece-wise constant model for the Alfven speed and a Reduced MHD - Shell model to describe the interplay between turbulent dynamics in the direction perpendicular to the mean field and propagation along the field. Turbulence is sustained by incoming fluctuations which are equivalent, in the line-tied case, to forcing by the photospheric shear flows. While varying the turbulence strength, we compare systematically the average coronal energy level (E) and dissipation rate (D) in three models with increasing complexity: the classical closed model, the semi-open corona model, and the corona-chromosphere (or 3-layer) model, the latter two models allowing energy leakage. We find that: (i) Leakage always plays a role (even for strong turbulence), E and D are systematically lower than in the line-tied model. (ii) E is close to the resonant prediction, i.e., assuming effective turbulent correlation time longer than the Alfven coronal crossing time (Ta). (iii) D is close to the value given by the ratio of photospheric energy divided by Ta (iv) The coronal spectra exibits an inertial range with 5/3 spectral slope, and a large scale peak of trapped resonant modes that inhibit nonlinear couplings. (v) In the realistic 3-layer model, the two-component spectrum leads to a damping time equal to the Kolmogorov time reduced by a factor u_rms/Va_corona

Sympathetic eruptions on the Sun have been observed for several decades, but the mechanisms by which one eruption can trigger another one remain poorly understood. We present a 3D MHD simulation that suggests two possible magnetic trigger mechanisms for sympathetic eruptions. We consider a configuration that contains two coronal flux ropes located within a pseudo-streamer and one rope located next to it. A sequence of eruptions is initiated by triggering the eruption of the flux rope next to the streamer. The expansion of the rope leads to two consecutive reconnection events, each of which triggers the eruption of a flux rope by removing a sufficient amount of overlying flux. The simulation qualitatively reproduces important aspects of the global sympathetic event on 2010 August 1 and provides a scenario for so-called twin filament eruptions. The suggested mechanisms are applicable also for sympathetic eruptions occurring in other magnetic configurations.

The Parker model for coronal heating is investigated through a high resolution simulation. An inertial range is resolved where fluctuating magnetic energy E_M (k_perp) \propto k_\perp^{-2.7} exceeds kinetic energy E_K (k_\perp) \propto k_\perp^{-0.6}. Increments scale as \delta b_\ell \simeq \ell^{-0.85} and \delta u_\ell \simeq \ell^{+0.2} with velocity increasing at small scales, indicating that magnetic reconnection plays a prime role in this turbulent system. We show that spectral energy transport is akin to standard magnetohydrodynamic (MHD) turbulence even for a system of reconnecting current sheets sustained by the boundary. In this new MHD turbulent cascade, kinetic energy flows are negligible while cross-field flows are enhanced, and through a series of "reflections" between the two fields, cascade more than half of the total spectral energy flow.

We present a series of numerical simulations aimed at understanding the nature and origin of turbulence in coronal loops in the framework of the Parker model for coronal heating. A coronal loop is studied via reduced magnetohydrodynamics simulations in Cartesian geometry. A uniform and strong magnetic field threads the volume between the two photospheric planes, where a velocity field in the form of a 1D shear flow pattern is present. Initially the magnetic field which developes in the coronal loop is a simple map of the photospheric velocity field. This initial configuration is unstable to a multiple tearing instability which develops islands with X and O points in the plane orthogonal to the axial field. Once the nonlinear stage sets in the system evolution is characterized by a regime of MHD turbulence dominated by magnetic energy. A well developed power law in energy spectra is observed and the magnetic field never returns to the simple initial state mapping the photospheric flow. The formation of X and O points in the planes orthogonal to the axial field allows the continued and repeated formation and dissipation of small scale current sheets where the plasma is heated. We conclude that the observed turbulent dynamics are not induced by the complexity of the pattern that the magnetic field-lines footpoints follow but they rather stem from the inherent nonlinear nature of the system.

The plasma density enhancements recently observed by the Large-Angle Spectrometric Coronagraph (LASCO) instrument onboard the Solar and Heliospheric Observatory (SOHO) spacecraft have sparked considerable interest. In our previous theoretical study of the formation and initial motion of these density enhancements it is found that beyond the helmet cusp of a coronal streamer the magnetized wake configuration is resistively unstable, that a traveling magnetic island develops at the center of the streamer, and that density enhancements occur within the magnetic islands. As the massive magnetic island travels outward, both its speed and width increase. The island passively traces the acceleration of the inner part of the wake. In the present paper a few spherical geometry effects are included, taking into account either the radial divergence of the magnetic field lines and the average expansion suffered by a parcel of plasma propagating outward, using the Expanding Box Model (EBM), and the diamagnetic force due to the overall magnetic field radial gradients, the so-called melon-seed force. It is found that the values of the acceleration and density contrasts can be in good agreement with LASCO observations, provided the spherical divergence of the magnetic lines starts beyond a critical distance from the Sun and the initial stage of the formation and acceleration of the plasmoid is due to the cartesian evolution of MHD instabilities. This result provides a constraint on the topology of the magnetic field in the coronal streamer.

To understand the nonlinear dynamics of the Parker scenario for coronal heating, long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry are carried out. A loop is modeled as a box extended along the direction of the strong magnetic field $B_0$ in which the system is embedded. At the top and bottom plates, which represent the photosphere, velocity fields mimicking photospheric motions are imposed. We show that the nonlinear dynamics is described by different regimes of MHD anisotropic turbulence, with spectra characterized by intertial range power laws whose indexes range from Kolmogorov-like values ($\sim 5/3$) up to $\sim 3$. We briefly describe the bearing for coronal heating rates.

The Parker or field line tangling model of coronal heating is investigated through long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry within the framework of reduced magnetohydrodynamics (RMHD). Slow photospheric motions induce a Poynting flux which saturates by driving an anisotropic turbulent cascade dominated by magnetic energy and characterized by current sheets elongated along the axial magnetic field. Increasing the value of the axial magnetic field different regimes of MHD turbulence develop with a bearing on coronal heating rates. In physical space magnetic field lines at the scale of convection cells appear only slightly bended in agreement with observations of large loops of current (E)UV and X-ray imagers.

In a recent work Grappin et al. [1] have shown that low- frequency movements can be transmitted from one footpoint to the other along a magnetic loop, thus mimicking a friction effect of the corona on the photosphere, and invalidating the line-tying approximation. We consider here successively the effect of high frequencies and turbulent damping on the process. We use a very simple atmospheric model which allows to study analytically the laminar case, and to study the turbulent case both using simple phenomenological arguments and a more sophisticated turbulence model [2]. We find that, except when turbulent damping is such that all turbulence is damped during loop traversal, coupling still occurs between distant footpoints, and moreover the coronal field induced by photospheric movements saturates at finite values.

MHD turbulence has long been proposed as a mechanism for the heating of coronal loops in the framework of the Parker scenario for coronal heating. So far most of the studies have focused on its dynamical properties without considering its thermodynamical and radiative features, because of the very demanding computational requirements. In this paper we extend this previous research to the compressible regime, including an energy equation, by using HYPERION, a new parallelized, viscoresistive, three-dimensional compressible MHD code. HYPERION employs a Fourier collocation -- finite difference spatial discretization, and uses a third-order Runge-Kutta temporal discretization. We show that the implementation of a thermal conduction parallel to the DC magnetic field induces a radiative emission concentrated at the boundaries, with properties similar to the chromosphere--transition region--corona system.

A model is presented for generation of fast solar wind in coronal holes, relying on heating that is dominated by turbulent dissipation of MHD fluctuations transported upwards in the solar atmosphere. Scale-separated transport equations include large-scale fields, transverse Alfvenic fluctuations, and a small compressive dissipation due to parallel shears near the transition region. The model accounts for proton temperature, density, wind speed, and fluctuation amplitude as observed in remote sensing and in situ satellite data.

The problem of three-dimensional combined magnetic and velocity shear driven instabilities of a compressible magnetized jet modeled with a plane neutral/current double vortex sheet in the framework of the resistive magnetohydrodynamics is addressed. The resulting dynamics given by the stream+current sheet interaction is analyzed and the effects of a variable geometry of the basic fields are considered. Depending on the basic asymptotic magnetic field configuration, a selection rule of the linear instability modes can be obtained. Hence, the system follows a two-stage path developing either through a fully three-dimensional dynamics with a rapid evolution of kink modes leading to a final turbulent state, or rather through a driving two-dimensional instability pattern that develops on parallel planes on which a reconnection+coalescence process takes place.

We study the formation and evolution of a turbulent spectrum of Alfv\'en waves driven by reflection off the solar wind density gradients, starting from the coronal base up to 17 solar radii, well beyond the Alfv\'enic critical point. The background solar wind is assigned and 2D shell models are used to describe nonlinear interactions. We find that the turbulent spectra are influenced by the nature of reflected waves. Close to the base, these give rise to a flatter and steeper spectrum for the outgoing and reflected waves respectively. At higher heliocentric distance both spectra evolve toward an asymptotic Kolmogorov spectrum. The turbulent dissipation is found to account for at least half of the heating required to sustain the background imposed solar wind and its shape is found to be determined by the reflection-determined turbulent heating below 1.5 solar radii. Therefore reflection and reflection-driven turbulence are shown to play a key role in the accelerationof the fast solar wind and origin of the turbulent spectrum found at 0.3 AU in the heliosphere.

The POLAR Investigation of the Sun (POLARIS) mission uses a combination of a gravity assist and solar sail propulsion to place a spacecraft in a 0.48 AU circular orbit around the Sun with an inclination of 75 degrees with respect to solar equator. This challenging orbit is made possible by the challenging development of solar sail propulsion. This first extended view of the high-latitude regions of the Sun will enable crucial observations not possible from the ecliptic viewpoint or from Solar Orbiter. While Solar Orbiter would give the first glimpse of the high latitude magnetic field and flows to probe the solar dynamo, it does not have sufficient viewing of the polar regions to achieve POLARIS' primary objective : determining the relation between the magnetism and dynamics of the Sun's polar regions and the solar cycle.

The Parker or field line tangling model of coronal heating is studied comprehensively via long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry within the framework of reduced magnetohydrodynamics (RMHD). Slow photospheric motions induce a Poynting flux which saturates by driving an anisotropic turbulent cascade dominated by magnetic energy. In physical space this corresponds to a magnetic topology where magnetic field lines are barely entangled, nevertheless current sheets (corresponding to the original tangential discontinuities hypothesized by Parker) are continuously formed and dissipated. Current sheets are the result of the nonlinear cascade that transfers energy from the scale of convective motions ($\sim 1,000 km$) down to the dissipative scales, where it is finally converted to heat and/or particle acceleration. Current sheets constitute the dissipative structure of the system, and the associated magnetic reconnection gives rise to impulsive ``bursty'' heating events at the small scales. This picture is consistent with the slender loops observed by state-of-the-art (E)UV and X-ray imagers which, although apparently quiescent, shine bright in these wavelengths with little evidence of entangled features. The different regimes of weak and strong MHD turbulence that develop, and their influence on coronal heating scalings, are shown to depend on the loop parameters, and this dependence is quantitatively characterized: weak turbulence regimes and steeper spectra occur in {\it stronger loop fields} and lead to {\it larger heating rates} than in weak field regions.

Long-time high-resolution simulations of the dynamics of a coronal loop in cartesian geometry are carried out, within the framework of reduced magnetohydrodynamics (RMHD), to understand coronal heating driven by motion of field lines anchored in the photosphere. We unambiguously identify MHD anisotropic turbulence as the physical mechanism responsible for the transport of energy from the large scales, where energy is injected by photospheric motions, to the small scales, where it is dissipated. As the loop parameters vary different regimes of turbulence develop: strong turbulence is found for weak axial magnetic fields and long loops, leading to Kolmogorov-like spectra in the perpendicular direction, while weaker and weaker regimes (steeper spectral slopes of total energy) are found for strong axial magnetic fields and short loops. As a consequence we predict that the scaling of the heating rate with axial magnetic field intensity $B_0$, which depends on the spectral index of total energy for given loop parameters, must vary from $B_0^{3/2}$ for weak fields to $B_0^{2}$ for strong fields at a given aspect ratio. The predicted heating rate is within the lower range of observed active region and quiet Sun coronal energy losses.

Context: The location of coronal heating in magnetic loops has been the subject of a long-lasting controversy: does it occur mostly at the loop footpoints, at the top, is it random, or is the average profile uniform? Aims: We try to address this question in model loops with MHD turbulence and a profile of density and/or magnetic field along the loop. Methods: We use the ShellAtm MHD turbulent heating model described in Buchlin & Velli (2006), with a static mass density stratification obtained by the HydRad model (Bradshaw & Mason 2003). This assumes the absence of any flow or heat conduction subsequent to the dynamic heating. Results: The average profile of heating is quasi-uniform, unless there is an expansion of the flux tube (non-uniform axial magnetic field) or the variation of the kinetic and magnetic diffusion coefficients with temperature is taken into account: in the first case the heating is enhanced at footpoints, whereas in the second case it is enhanced where the dominant diffusion coefficient is enhanced. Conclusions: These simulations shed light on the consequences on heating profiles of the complex interactions between physical effects involved in a non-uniform turbulent coronal loop.

We solve the problem of propagation and dissipation of Alfvenic turbulence in a model solar atmosphere consisting of a static photosphere and chromosphere, transition region, and open corona and solar wind, using a phenomenological model for the turbulent dissipation based on wave reflection. We show that most of the dissipation for a given wave-frequency spectrum occurs in the lower corona, and the overall rms amplitude of the fluctuations evolves in a way consistent with observations. The frequency spectrum, for a Kolmogorov-like slope, is not found to change dramatically from the photosphere to the solar wind, however it does preserve signatures of transmission throughout the lower atmospheric layers, namely oscillations in the spectrum at high frequencies reminiscent of the resonances found in the linear case. These may disappear once more realistic couplings for the non-linear terms are introduced, or if time-dependent variability of the lower atmospheric layer is introduced.

A simplified non-linear numerical model for the development of incompressible magnetohydrodynamics (MHD) in the presence of a strong magnetic field B0 and stratification, nicknamed Shell-Atm, is presented. In planes orthogonal to the mean field, the non-linear incompressible dynamics is replaced by 2D shell-models for the complex variables u and b, allowing one to reach large Reynolds numbers while at the same time carrying out sufficiently long time integrations to obtain a good statistics at moderate computational cost. The shell-models of different planes are coupled by Alfven waves propagating along B0. The model may be applied to open or closed magnetic field configurations where the axial field dominates and the plasma pressure is low; here we apply it to the specific case of a magnetic loop of the solar corona heated via turbulence driven by photospheric motions, and we use statistics for its analysis. The Alfven waves interact non-linearly and form turbulent spectra in the directions perpendicular and, via propagation, also parallel to the mean field. A heating function is obtained, and is shown to be intermittent; the average heating is consistent with values required for sustaining a hot corona, and is proportional to the aspect ratio of the loop to the power -1.5; characteristic properties of heating events are distributed as power-laws. Cross-correlations show a delay of dissipation compared to energy content.

We analyze low and high--latitude fast solar wind data from the Ulysses spacecraft from 1992 to 1994 using a a systematic method to analyse the anisotropic content of the magnetic field fluctuations. We investigate all available frequencies, 1-10^{-6} Hz, for both high and low--latitudes datasets and are able to quantify the relative importance of the anisotropic versus the isotropic fluctuations. We analyse, up to sixth order, longitudinal, transverse and mixed magnetic field correlations. Our results show that strongly intermittent and anisotropic events are present in the solar wind plasma at high frequencies/small scales, indicating the absence of a complete recovery of isotropy. Anisotropic scaling properties are compatible for high and low--latitude data, suggesting a universal behaviour in spite of the different rate of evolution of the fast solar wind streams in the two environments.

The propagation and the evolution of Alfvenic pulses in the solar coronal arcades is investigated by means of MHD numerical simulations. Significant transverse oscillations in coronal loops, triggered by nearby flare events, are often measured in EUV lines and are generally interpreted as standing kink modes. However, the damping times of these oscillations are typically very short (from one to a few periods) and the physical mechanism responsible for the decay is still a matter of debate. Moreover, the majority of the observed cases actually appears to be better modeled by propagating, rather than standing, modes. Here we perform 2.5-D compressible MHD simulations of impulsively generated Alfven waves propagating in a potential magnetic arcade (assumed as a simplified 2-D loop model), taking into account the stratification of the solar atmosphere with height from the photosphere to the corona. The results show a strong spreading of the initially localized pulses along the loop, due to the variations in the Alfven velocity with height, and correspondingly an efficient damping of the amplitude of the oscillations. We believe that simple explanations based on the effects of wave propagation in highly inhomogeneous media may apply to the majority of the reported cases, and that variations of the background density and Alfven speed along the loop should be considered as key ingredients in future models.

A convenient and widely used way to study the turbulent plasma in the solar corona is to do statistics of properties of events (or structures), associated with flares, that can be found in observations or in numerical simulations. Numerous papers have followed such a methodology, using different definitions of an event, but the reasons behind the choice of a particular definition (and not another one) is very rarely discussed. We give here a comprehensive set of possible event definitions starting from a one-dimensional data set such as a time-series of energy dissipation. Each definition is then applied to a time-series of energy dissipation issued from simulations of a shell-model of magnetohydrodynamic turbulence as defined in Giuliani and Carbone (1998), or from a new model of coupled shell-models designed to represent a magnetic loop in the solar corona. We obtain distributions of the peak dissipation power, total energy, duration and waiting-time associated to each definition. These distributions are then investigated and compared, and the influence of the definition of an event on statistics is discussed. In particular, power-law distributions are more likely to appear when using a threshold. The sensitivity of the distributions to the definition of an event seems also to be weaker for events found in a highly intermittent time series. Some implications on statistical results obtained from observations are discussed.

To model the interaction between the solar wind and the interstellar wind, magnetic fields must be included. Recently Opher et al. 2003 found that, by including the solar magnetic field in a 3D high resolution simulation using the University of Michigan BATS-R-US code, a jet-sheet structure forms beyond the solar wind Termination Shock. Here we present an even higher resolution three-dimensional case where the jet extends for $150AU$ beyond the Termination Shock. We discuss the formation of the jet due to a de Laval nozzle effect and it's su bsequent large period oscillation due to magnetohydrodynamic instabilities. To verify the source of the instability, we also perform a simplified two dimensional-geometry magnetohydrodynamic calculation of a plane fluid jet embedded in a neutral sheet with the profiles taken from our 3D simulation. We find remarkable agreement with the full three-dimensional evolution. We compare both simulations and the temporal evolution of the jet showing that the sinuous mode is the dominant mode that develops into a velocity-shear-instability with a growth rate of $5 \times 10^{-9} sec^{-1}=0.027 years^{-1}$. As a result, the outer edge of the heliosphere presents remarkable dynamics, such as turbulent flows caused by the motion of the jet. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how this outer boundary behaves. Already, however, we can say that the magnetic field effects are a major player in this region changing our previous notion of how the solar system ends.

There is currently a controversy as to whether Voyager 1 has already crossed the Termination Shock, the first boundary of the Heliosphere. The region between the Termination Shock and the Heliopause, the Helisheath, is one of the most unknown regions theoretically. In the Heliosheath magnetic effects are crucial, as the solar magnetic field is compressed at the Termination Shock by the slowing flow. Recently, our simulations showed that the Heliosheath presents remarkable dynamics, with turbulent flows and the presence of a jet flow at the current sheet that is unstable due to magnetohydrodynamic instabilities \cite{opher,opher1}. In this paper we review these recent results, and present an additional simulation with constant neutral atom background. In this case the jet is still present but with reduced intensity. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how the Heliosheath behaves. Already we can say that this region presents remarkable dynamics, with turbulent flows, indicating that the Heliosheath might be very different from what we previously thought.

A 3D model intermediate between cellular automata (CA) models and the reduced magnetohydrodynamic (RMHD) equations is presented to simulate solar impulsive events generated along a coronal magnetic loop. The model consists of a set of planes distributed along a magnetic loop between which the information propagates through Alfven waves. Statistical properties in terms of power-laws for energies and durations of dissipative events are obtained, and their agreement with X-ray and UV flares observations is discussed. The existence of observational biases is also discussed.