Source author record

Marco Velli

Marco Velli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2026arXiv

Properties of Magnetic Switchbacks in the Near-Sun Solar Wind

Magnetic switchbacks are fluctuations in the solar wind in which the interplanetary magnetic field sharply deflects away from its background direction so as to create folds in magnetic field lines while remaining of roughly constant magnitude. The magnetic field and velocity fluctuations are extremely well correlated in a way corresponding to Alfvénic fluctuations propagating away from the Sun. For a background field which is nearly radial this causes an outwardly propagating jet to form. Switchbacks and their characteristic velocity jets have recently been observed to be nearly ubiquitous by Parker Solar Probe with in situ measurements in the inner heliosphere within 0.3 AU. Their prevalence, substantial energy content, and potentially fundamental role in the dynamics of the outer corona and solar wind motivate the significant research efforts into their understanding. Here we review the in situ measurements of these structures (primarily by Parker Solar Probe). We discuss how they are identified and measured, and present an overview of the primary observational properties of these structures, both in terms of individual switchbacks and their collective arrangement into ``patches''. We identify both properties for which there is a strong consensus and those that have limited or qualified support and require further investigation. We identify and collate several open questions and recommendations for future studies.

preprint2023arXiv

Formation of Magnetic Switchbacks Observed by Parker Solar Probe

Magnetic switchbacks are rapid high amplitude reversals of the radial magnetic field in the solar wind that do not involve a heliospheric current sheet crossing. First seen sporadically in the seventies in Mariner and Helios data, switchbacks were later observed by the Ulysses spacecraft beyond 1 au and have been recently identified as a typical component of solar wind fluctuations in the inner heliosphere by the Parker Solar Probe spacecraft. Here we provide a simple yet predictive theory for the formation of these magnetic reversals: the switchbacks are produced by the shear of circularly polarized Alfvén waves by a transversely varying radial wave propagation velocity. We provide an analytic expression for the magnetic field variation, establish the necessary and sufficient conditions and show that the mechanism works in a realistic solar wind scenario.

preprint2023arXiv

Solaris: A Focused Solar Polar Discovery-class Mission to achieve the Highest Priority Heliophysics Science Now

Solaris is a transformative Solar Polar Discovery-class mission concept to address crucial outstanding questions that can only be answered from a polar vantage. Solaris will image the Sun's poles from ~75 degree latitude, providing new insight into the workings of the solar dynamo and the solar cycle, which are at the foundation of our understanding of space weather and space climate. Solaris will also provide enabling observations for improved space weather research, modeling and prediction, revealing a unique, new view of the corona, coronal dynamics and CME eruptions from above.

preprint2022arXiv

Conservation of Total Wave Action in the Expanding Solar Wind

The conservation of wave action in moving plasmas has been well-known for over half a century. However, wave action is not conserved when multiple wave modes propagate and coexist close to degeneration condition (Sound speed equals Alfvén speed, i.e. plasma $β\sim 1$). Here we show that the violation of conservation is due to wave mode conversion, and that the total wave action summed over interacting modes is still conserved. Though the result is general, we focus on MHD waves and identify three distinctive mode conversion mechanisms, i.e. degeneracy, linear mode conversion, and resonance, and provide an intuitive physical picture for the mode conversion processes. We use 1D MHD simulations with the Expanding Box Model to simulate the nonlinear evolution of monochromatic MHD waves in the expanding solar wind. Simulation results validate the theory; total wave action therefore remains an interesting diagnostic for studies of waves and turbulence in the solar wind.

preprint2022arXiv

Constraining Global Coronal Models with Multiple Independent Observables

Global coronal models seek to produce an accurate physical representation of the Sun's atmosphere which can be used, for example, to drive space weather models. Assessing their accuracy is a complex task and there are multiple observational pathways to provide constraints and tune model parameters. Here, we combine several such independent constraints, defining a model-agnostic framework for standardized comparison. We require models to predict the distribution of coronal holes at the photosphere, and neutral line topology at the model outer boundary. We compare these predictions to extreme ultraviolet (EUV) observations of coronal hole locations, white-light Carrington maps of the streamer belt and the magnetic sector structure measured \textit{in situ} by Parker Solar Probe and 1AU spacecraft. We study these metrics for Potential Field Source Surface (PFSS) models as a function of source surface height and magnetogram choice, as well as comparing to the more physical Wang-Sheeley-Arge (WSA) and the Magnetohydrodynamics Algorithm outside a Sphere (MAS) models. We find that simultaneous optimization of PFSS models to all three metrics is not currently possible, implying a trade-off between the quality of representation of coronal holes and streamer belt topology. WSA and MAS results show the additional physics they include addresses this by flattening the streamer belt while maintaining coronal hole sizes, with MAS also improving coronal hole representation relative to WSA. We conclude that this framework is highly useful for inter- and intra-model comparisons. Integral to the framework is the standardization of observables required of each model, evaluating different model aspects.

preprint2022arXiv

Electron-driven instabilities in the solar wind

The electrons are an essential particle species in the solar wind. They often exhibit non-equilibrium features in their velocity distribution function. These include temperature anisotropies, tails (kurtosis), and reflectional asymmetries (skewness), which contribute a significant heat flux to the solar wind. If these non-equilibrium features are sufficiently strong, they drive kinetic micro-instabilities. We develop a semi-graphical framework based on the equations of quasi-linear theory to describe electron-driven instabilities in the solar wind. We apply our framework to resonant instabilities driven by temperature anisotropies. These include the electron whistler anisotropy instability and the propagating electron firehose instability. We then describe resonant instabilities driven by reflectional asymmetries in the electron distribution function. These include the electron/ion-acoustic, kinetic Alfvén heat-flux, Langmuir, electron-beam, electron/ion-cyclotron, electron/electron-acoustic, whistler heat-flux, oblique fast-magnetosonic/whistler, lower-hybrid fan, and electron-deficit whistler instability. We briefly comment on non-resonant instabilities driven by electron temperature anisotropies such as the mirror-mode and the non-propagating firehose instability. We conclude our review with a list of open research topics in the field of electron-driven instabilities in the solar wind.

preprint2022arXiv

Influence of the heliospheric current sheet on the evolution of solar wind turbulence

The effects of the heliospheric current sheet (HCS) on the evolution of Alfvénic turbulence in the solar wind are studied using MHD simulations incorporating the expanding-box-model (EBM). The simulations show that near the HCS, the Alfvénicity of the turbulence decreases as manifested by lower normalized cross helicity and larger excess of magnetic energy. The numerical results are supported by a superposed-epoch analysis using OMNI data, which shows that the normalized cross helicity decreases inside the plasma sheet surrounding HCS, and the excess of magnetic energy is significantly enhanced at the center of HCS. Our simulation results indicate that the decrease of Alfvénicity around the HCS is due to the weakening of radial magnetic field and the effects of the transverse gradient in the background magnetic field. The magnetic energy excess in the turbulence may be a result of the loss of Alfvénic correlation between velocity and magnetic field and the faster decay of transverse kinetic energy with respect to magnetic energy in a spherically expanding solar wind.

preprint2022arXiv

Magnetic field intermittency in the solar wind: PSP and SolO observations ranging from the Alfven region out to 1 AU

$PSP$ and $SolO$ data are utilized to investigate magnetic field intermittency in the solar wind (SW). Small-scale intermittency $(20-100d_{i})$ is observed to radially strengthen when methods relying on higher-order moments are considered ($SF_q$, $SDK$), but no clear trend is observed at larger scales. However, lower-order moment-based methods (e.g., PVI) are deemed more appropriate for examining the evolution of the bulk of Coherent Structures (CSs), $PVI \ge 3$. Using PVI, we observe a scale-dependent evolution in the fraction of the dataset occupied by CSs, $f_{PVI \ge 3}$. Specifically, regardless of the SW speed, a subtle increase is found in $f_{PVI\ge3}$ for $\ell =20 d_i$, in contrast to a more pronounced radial increase in CSs observed at larger scales. Intermittency is investigated in relation to plasma parameters. Though, slower SW speed intervals exhibit higher $f_{PVI \geq 6}$ and higher kurtosis maxima, no statistical differences are observed for $f_{PVI \geq 3}$. Highly Alfvénic intervals, display lower levels of intermittency. The anisotropy with respect to the angle between the magnetic field and SW flow, $Θ_{VB}$ is investigated. Intermittency is weaker at $Θ_{VB} \approx 0^{\circ}$ and is strengthened at larger angles. Considering the evolution at a constant alignment angle, a weakening of intermittency is observed with increasing advection time of the SW. Our results indicate that the strengthening of intermittency in the inner heliosphere is driven by the increase in comparatively highly intermittent perpendicular intervals sampled by the probes with increasing distance, an effect related directly to the evolution of the Parker spiral.

preprint2022arXiv

Patches of magnetic switchbacks and their origins

Parker Solar Probe (PSP) has shown that the solar wind in the inner heliosphere is characterized by the quasi omni-presence of magnetic switchbacks ("switchback" hereinafter), local backward-bends of magnetic field lines. Switchbacks also tend to come in patches, with a large-scale modulation that appears to have a spatial scale size comparable to supergranulation on the Sun. Here we inspect data from the first ten encounters of PSP focusing on different time intervals when clear switchback patches were observed by PSP. We show that the switchbacks modulation, on a timescale of several hours, seems to be independent of whether PSP is near perihelion, when it rapidly traverses large swaths of longitude remaining at the same heliocentric distance, or near the radial-scan part of its orbit, when PSP hovers over the same longitude on the Sun while rapidly moving radially inwards or outwards. This implies that switchback patches must also have an intrinsically temporal modulation most probably originating at the Sun. Between two consecutive patches, the magnetic field is usually very quiescent with weak fluctuations. We compare various parameters between the quiescent intervals and the switchback intervals. The results show that the quiescent intervals are typically less Alfvénic than switchback intervals, and the magnetic power spectrum is usually shallower in quiescent intervals. We propose that the temporal modulation of switchback patches may be related to the "breathing" of emerging flux that appears in images as the formation of "bubbles" below prominences in the Hinode/SOT observations.

preprint2022arXiv

Preferential Proton over Electron heating from coherent structures during the first perihelion of Parker Solar Probe

The solar wind undergoes significant heating as it propagates away from the Sun; the exact mechanisms responsible for this heating remain unclear. Using data from the first perihelion of the Parker Solar Probe mission, we examine the properties of proton and electron heating occurring within magnetic coherent structures identified by means of the Partial Variance of Increments (PVI) method. Statistically, regions of space with strong gradients in the magnetic field, $PVI \geq 1$, are associated with strongly enhanced proton but only slightly elevated electron temperatures. Our analysis indicates a heating mechanism in the nascent solar wind environment facilitated by a nonlinear turbulent cascade that preferentially heats protons over electrons.

preprint2022arXiv

Statistical analysis of intermittency and its association with proton heating in the near Sun environment

We use data from the first six encounters of Parker Solar Probe and employ the Partial Variance of Increments ($PVI$) method to study the statistical properties of coherent structures in the inner heliosphere with the aim of exploring physical connections between magnetic field intermittency and observable consequences such as plasma heating and turbulence dissipation. Our results support proton heating localized in the vicinity of, and strongly correlated with, magnetic structures characterized by $PVI \geq 1$. We show that on average, such events constitute $\approx 19\%$ of the dataset, though variations may occur depending on the plasma parameters. We show that the waiting time distribution ($WT$) of identified events is consistent across all six encounters following a power-law scaling at lower $WTs$. This result indicates that coherent structures are not evenly distributed in the solar wind but rather tend to be tightly correlated and form clusters. We observe that the strongest magnetic discontinuities, $PVI \geq 6$, usually associated with reconnection exhausts, are sites where magnetic energy is locally dissipated in proton heating and are associated with the most abrupt changes in proton temperature. However, due to the scarcity of such events, their relative contribution to energy dissipation is minor. Taking clustering effects into consideration, we show that smaller scale, more frequent structures with PVI between, $1\lesssim PVI \lesssim 6$, play the major role in magnetic energy dissipation. The number density of such events is strongly associated with the global solar wind temperature, with denser intervals being associated with higher $T_{p}$.

preprint2022arXiv

The role of Alfvén wave dynamics on the large scale properties of the solar wind: comparing an MHD simulation with PSP E1 data

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énic 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én 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.

preprint2022arXiv

Thin current sheet formation: comparison between Earth's magnetotail and coronal streamers

Magnetic field line reconnection is a universal plasma process responsible for the magnetic field topology change and magnetic field energy dissipation into charged particle heating and acceleration. In many systems, the conditions leading to the magnetic reconnection are determined by the pre-reconnection configuration of a thin layer with intense currents -- otherwise known as the thin current sheet. In this study we investigate two such systems: Earth's magnetotail and helmet streamers in the solar corona. The pre-reconnection current sheet evolution has been intensely studied in the magnetotail, where in-situ spacecraft observations are available; but helmet streamer current sheets studies are fewer, due to lack of in-situ observations -- they are mostly investigated with numerical simulations and information that can be surmised from remote sensing instrumentation. Both systems exhibit qualitatively the same behavior, despite their largely different Mach numbers, much higher at the solar corona than at the magnetotail. Comparison of spacecraft data (from the magnetotail) with numerical simulations (for helmet streamers) shows that the pre-reconnection current sheet thinning, for both cases, is primarily controlled by plasma pressure gradients. Scaling laws of the current density, magnetic field, and pressure gradients are the same for both systems. We discuss how magnetotail observations and kinetic simulations can be utilized to improve our understanding and modeling of the helmet streamer current sheets.

preprint2021arXiv

Comparative study of electric currents and energetic particle fluxes in a solar flare and Earth magnetospheric substorm

Magnetic field-line reconnection is a universal plasma process responsible for the conversion of magnetic field energy to the plasma heating and charged particle acceleration. Solar flares and Earth's magnetospheric substorms are two most investigated dynamical systems where magnetic reconnection is believed to be responsible for global magnetic field reconfiguration and energization of plasma populations. Such a reconfiguration includes formation of a long-living current systems connecting the primary energy release region and cold dense conductive plasma of photosphere/ionosphere. In both flares and substorms the evolution of this current system correlates with formation and dynamics of energetic particle fluxes. Our study is focused on this similarity between flares and substorms. Using a wide range of datasets available for flare and substorm investigations, we compare qualitatively dynamics of currents and energetic particle fluxes for one flare and one substorm. We showed that there is a clear correlation between energetic particle bursts (associated with energy release due to magnetic reconnection) and magnetic field reconfiguration/formation of current system. We then discuss how datasets of in-situ measurements in the magnetospheric substorm can help in interpretation of datasets gathered for the solar flare.

preprint2021arXiv

Deciphering Solar Magnetic Activity: 140 Years Of The `Extended Solar Cycle' -- Mapping the Hale Cycle

We investigate the occurrence of the "extended solar cycle" (ESC) as it occurs in a host observational data spanning 140 years. Investigating coronal, chromospheric, photospheric and interior diagnostics we develop a consistent picture of solar activity migration linked to the 22-year Hale (magnetic) cycle using superposed epoch analysis (SEA) using previously identified Hale cycle termination events as the key time for the SEA. Our analysis shows that the ESC and Hale cycle, as highlighted by the terminator-keyed SEA, is strongly recurrent throughout the entire observational record studied, some 140 years. Applying the same SEA method to the sunspot record confirms that Maunder's butterfly pattern is a subset of the underlying Hale cycle, strongly suggesting that the production of sunspots is not the fundamental feature of the Hale cycle, but the ESC is. The ESC (and Hale cycle) pattern highlights the importance of 55\degree\ latitude in the evolution, and possible production, of solar magnetism.

preprint2021arXiv

Investigating the origin of the FIP effect with a shell turbulence model

The enrichment of coronal loops and the slow solar wind with elements that have low First Ionisation Potential, known as the FIP effect, has often been interpreted as the tracer of a common origin. A current explanation for this FIP fractionation rests on the influence of ponderomotive forces and turbulent mixing acting at the top of the chromosphere. The implied wave transport and turbulence mechanisms are also key to wave-driven coronal heating and solar wind acceleration models. This work makes use of a shell turbulence model run on open and closed magnetic field lines of the solar corona to investigate with a unified approach the influence of magnetic topology, turbulence amplitude and dissipation on the FIP fractionation. We try in particular to assess whether there is a clear distinction between the FIP effect on closed and open field regions.

preprint2020arXiv

Collisionless heat flux regulation via electron firehose instability in presence of a core and suprathermal population in the expanding solar wind

The evolution of the electron heat flux in the solar wind is regulated by the interplay between several effects: solar wind expansion, that can potentially drive velocity-space instabilties, turbulence and wave-particle interactions, and, possibly, collisions. Here we address the respective role played by the solar wind expansion and the electron firehose instability, developing in the presence of multiple electron populations, in regulating the heat flux. We carry out fully kinetic, Expanding Box Model simulations and separately analyze the enthalpy, bulk and velocity distribution function skewness contributions for each of the electron species. We observe that the key factor determining electron energy flux evolution is the reduction of the drift velocity of the electron populations in the rest frame of the solar wind. In our simulations, redistribution of the electron thermal energy from the parallel to the perpendicular direction after the onset of the electron firehose instability is observed. However, this process seems to impact energy flux evolution only minimally. Hence, reduction of the electron species drift velocity in the solar wind frame appears to directly correlate with efficiency for heat flux instabilities

preprint2020arXiv

Proton Temperature Anisotropy Variations in Inner Heliosphere Estimated with First Parker Solar Probe Observations

We report proton temperature anisotropy variations in the inner heliosphere with Parker Solar Probe (PSP) observations. Using a linear fitting method, we derive proton temperature anisotropy with temperatures measured by the Solar Probe Cup (SPC) from the SWEAP instrument suite and magnetic field observations from the FIELDS instrument suite. The observed radial dependence of temperature variations in the fast solar wind implies stronger perpendicular heating and parallel cooling than previous results from Helios measurements made at larger radial distances. The anti-correlation between proton temperature anisotropy and parallel plasma beta is retained in fast solar wind. However, the temperature anisotropies of the slow solar wind seem to be well constrained by the mirror and parallel firehose instabilities. The perpendicular heating of the slow solar wind inside 0.24 AU may contribute to its same trend up against mirror instability thresholds as fast solar wind. These results suggest that we may see stronger anisotropy heating than expected in inner heliosphere.

preprint2020arXiv

Tearing instability and periodic density perturbations in the slow solar wind

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.

preprint2019arXiv

Magnetic field kinks and folds in the solar wind

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én 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énic 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én 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.

preprint2019arXiv

Propagation of Alfvén waves in the expanding solar wind with the fast-slow stream interaction

We carry out two-dimensional magnetohydrodynamic (MHD) simulations of an ensemble of Alfvénic 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énic fluctuations. The stream interaction is found to strongly affect the radial evolution of Alfvénic turbulence. The total energy in the Alfvén 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én 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.

preprint2016arXiv

"Ideally" unstable current sheets and the triggering of fast magnetic reconnection

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.

preprint2016arXiv

Closed-Field Coronal Heating Driven by Wave Turbulence

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.

preprint2015arXiv

"Ideal" tearing and the transition to fast reconnection in the weakly collisional MHD and EMHD regimes

This paper discusses the transition to fast growth of the tearing instability in thin current sheets in the collisionless limit where electron inertia drives the reconnection process. It has been previously suggested that in resistive MHD there is a natural maximum aspect ratio (ratio of sheet length and breadth to thickness) which may be reached for current sheets with a macroscopic length L, the limit being provided by the fact that the tearing mode growth time becomes of the same order as the Alfvèn time calculated on the macroscopic scale (Pucci and Velli (2014)). For current sheets with a smaller aspect ratio than critical the normalized growth rate tends to zero with increasing Lundquist number S, while for current sheets with an aspect ratio greater than critical the growth rate diverges with S. Here we carry out a similar analysis but with electron inertia as the term violating magnetic flux conservation: previously found scalings of critical current sheet aspect ratios with the Lundquist number are generalized to include the dependence on the ratio $(d_e/L)^2$ where de is the electron skin depth, and it is shown that there are limiting scalings which, as in the resistive case, result in reconnecting modes growing on ideal time scales. Finite Larmor Radius effects are then included and the rescaling argument at the basis of "ideal" reconnection is proposed to explain secondary fast reconnection regimes naturally appearing in numerical simulations of current sheet evolution.

preprint2015arXiv

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

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.

preprint2015arXiv

Magnetic reconnection: Self-similar current sheet collapse triggered by "ideal" tearing

We study, by means of MHD simulations, the onset and evolution of fast reconnection via the "ideal" tearing mode within a collapsing current sheet at high Lundquist numbers ($S\gg10^4$). We first confirm that as the collapse proceeds, fast reconnection is triggered well before a Sweet-Parker type configuration can form: during the linear stage plasmoids rapidly grow in a few Alfvén times when the predicted "ideal" tearing threshold $S^{-1/3}$ is approached from above; after the linear phase of the initial instability, X-points collapse and reform nonlinearly. We show that these give rise to a hierarchy of tearing events repeating faster and faster on current sheets at ever smaller scales, corresponding to the triggering of "ideal" tearing at the renormalized Lundquist number. In resistive MHD this process should end with the formation of sub-critical ($S \leq10^4$) Sweet Parker sheets at microscopic scales. We present a simple model describing the nonlinear recursive evolution which explains the timescale of the disruption of the initial sheet.

preprint2015arXiv

Resistive magnetohydrodynamic simulations of the ideal tearing mode

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 $γ\approx 0.6\,{τ_A}^{-1}$, where $τ_A$ is the ideal Alfvénic 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énic timescales, as required to explain the explosive flaring activity often observed in solar and astrophysical plasmas.

preprint2015arXiv

The tearing mode instability of thin current sheets: the transition to fast reconnection in the presence of viscosity

This paper studies the growth rate of reconnection instabilities in thin current sheets in the presence of both resistivity and viscosity. In a previous paper, Pucci and Velli (2014), it was argued that at sufficiently high Lundquist number S it is impossible to form current sheets with aspect ratios L/a which scale as $L/a\sim S^α$ with $α> 1/3$ because the growth rate of the tearing mode would then diverge in the ideal limit $S\rightarrow\infty$. Here we extend their analysis to include the effects of viscosity, (always present in numerical simulations along with resistivity) and which may play a role in the solar corona and other astrophysical environments. A finite Prandtl number allows current sheets to reach larger aspect ratios before becoming rapidly unstable in pile-up type regimes. Scalings with Lundquist and Prandtl numbers are discussed as well as the transition to kinetic reconnection

preprint2014arXiv

Application of a Solar Wind Model Driven by Turbulence Dissipation to a 2D Magnetic Field Configuration

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én 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.

preprint2014arXiv

Validating a Time-Dependent Wave-Turbulence-Driven Model of 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én 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}.

preprint2013arXiv

Apparent Solar Tornado - Like Prominences

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.

preprint2012arXiv

Coronal Pseudostreamers: Source of Fast or Slow Solar Wind?

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.

preprint2012arXiv

Origins of Rolling, Twisting and Non-Radial Propagation of Eruptive Solar Events

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.

preprint2011arXiv

Coronal heating in coupled photosphere-chromosphere-coronal systems: turbulence and leakage

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

preprint2009arXiv

Three-dimensional evolution of magnetic and velocity shear driven instabilities in a compressible magnetized jet

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.