Abstract
The enrichment of coronal loops and the slow solar wind with elements that have low First Ionization 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.
Introduction
The First Ionization Potential (FIP) effect is an enrichment of heavy elements with low-FIP such as Fe, Si, and Mg compared with photospheric abundances. It was initially measured in the solar wind and Solar Energetic Particles (SEPs) and later inferred from spectroscopic observations of the corona (see , ; ; ; , and references therein). The FIP bias, i.e., the ratio of coronal to photospheric abundances, is moreover mass independent. This means that processes below the transition region are strongly affecting the hydrostatic balance of the partially ionized chromosphere. Early on, explanations of the FIP effect involved diffusion, flows or some sort of turbulent mixing in the chromosphere to prevent any gravitational settling (see, e.g., ; Schwadron et al., 1999).
Schwadron et al. (1999) favored the hypothesis of turbulent wave heating in the chromosphere as a way to both prevent a mass-dependant fractionation and to obtain a low-FIP bias. Later on, proposed the ponderomotive acceleration, i.e., the time averaged gradient of electromagnetic fluctuations, as the origin of the FIP effect. The ponderomotive acceleration has this advantage that it may change sign and could explain the inverse FIP effect observed in low-mass stars (Wood and Linsky, 2010; ). Both processes rely on Alfvén waves propagating parallel and anti-parallel to the magnetic field, to trigger a turbulent cascade through non-linear interactions and heating. These wave populations naturally arise in coronal loops, where footprints motions excite the loop on both ends, but they are also expected in the open solar wind where reflection on large scale gradients (Velli et al., 1989; Zhou and Matthaeus, 1989) or compressible instabilities (Tenerani and Velli, 2013; Shoda et al., 2018b; ) create an sunward component from a purely anti-sunward wave packet.
The Ulysses spacecraft have shown clear composition differences between the fast and the slow wind components, the slow wind showing FIP biases close to the one observed in coronal loops (). After all, decades of observations (; Tu and Marsch, 1995; ) and modeling of fluctuations in the solar wind have now brought compelling evidence that turbulence is likely to be a fundamental ingredient of coronal heating and solar wind acceleration (see, e.g., Verdini and Velli, 2007; ; Shoda et al., 2018b; ). Perturbations, or waves are also essential to provide the additional acceleration giving birth to the fast wind ( km/s), through the ponderomotive force (; ; ; ). The low-FIP bias of the slow wind is often understood as the proof that the it originates in coronal loops, through exchange and magnetic reconfiguration (or reconnection) in the low corona (see, e.g., ). Yet, if turbulence is both responsible for the solar wind acceleration and the FIP effect, can we rule out a scenario where FIP fractionation occurs in purely open regions?
We investigate this question using a coupled modeling approach. First, we extract unidimensional profiles of open solar wind flux tubes and coronal loops using a multi-dimensional MHD code (). We then use the Shell-Atm code (; Verdini et al., 2009) to compute the propagation, cascade and dissipation of purely transverse perturbations of velocity and magnetic field along these profiles. We perform a parameter study where we vary the geometry, the initial perturbation amplitude and the initial injection scale, and discuss their effect on the turbulence properties in the chromosphere and transition region. We estimate the resulting FIP biases with an analytical fractionation model.
Background Wind and Transition Region Profiles
In this section, we describe the global MHD simulation used to extract the different wind profiles later input in our turbulence model. The code itself is a global MHD solver based on PLUTO () and our current implementation is described in details in and . The main purpose of the global simulation is to provide a realistic structure of the transition region, which plays an essential part in the present work. We rely on a single run, of a dipolar solar field of 5 G, which is typical of solar minimum configurations (see, e.g., ), and a uniform coronal heating prescribed with the following function:where , and . The value of corresponds to the global energy output of the solar wind.
The result of the simulation is shown in Figure 1. The left panel is a zoom on the transition region as a function of height and latitude. Contrarily with previous studies (; ), whose inner boundary was located at the base of the corona, we start the simulation at the photosphere with a temperature of 6000 K. We initialize the simulation from a one dimensional profile of the solar atmosphere obtained with the 1D version of the code described in and the boundary conditions are identical to . We use a spherical grid with a radial resolution that goes from to from the photosphere to 20 . There are 544 points in the radial and 512 in the latitudinal direction. The location of the transition region (TR) varies with the latitude and the temperature profiles. It is relatively constant in coronal holes at around , i.e., 2000 km above the photosphere. The minimum height corresponds to the highest coronal temperatures at the very edge of the helmet streamer. The height of the TR then increases within the core of the streamer up to , accordingly with lower coronal temperatures. These lower temperatures are likely due to increased density that result in increased cooling inside the loop (radiative losses are proportional to the squared density in the model, see ). This structure has already been described in the work of using a similar heating function.
FIGURE 1
The shock capturing ability of PLUTO, based on Riemann solvers, is essential to describe strong density and temperature gradients defining the transition region. have developed an additional module to propagate Alfvén waves in the corona. Yet, since this module is based on a Wentzel-Krimers-Brillouin approximation, it is not suited for a precise description of the transition region. To further study the role of turbulence in the composition of the solar wind, we thus rely on the Shell-Atm code, which solves the non-linear incompressible Alfvén wave equations on given profiles along magnetic fields lines, neglecting magneto-acoustic modes. We extracted four different solar wind profiles from our 2.5 D simulation, shown in Figure 1. In the following section, we describe the Shell-Atm simulations made from these profiles.
Shell Turbulence Model
Shell-Atm is a low dimension fluid turbulence code based on the shell technique (see ; ), that can model coronal loops (; ) and open solar wind solutions (Verdini et al., 2009). As fixed background profiles of flow speed, density and temperature, we use the four profiles extracted from our global MHD simulation. They are shown in Figure 2. Cases identified with CH correspond to open field regions in Figure 1. The CH1 profile is well in the coronal hole of the northern hemisphere in the fast wind region. CH2 is located close to the streamer and thus in the slow wind region of the simulation. LO cases are coronal loops, the LO1 case being at the open/close boundary with the minimum height of the transition region and LO2 being the smallest loop, well inside the helmet streamer.
FIGURE 2
Shell-Atm solves the reduced (incompressible) evolution of Alfvénic fluctuations, i.e. two coupled equations for the evolution of the parallel and anti-parallel wave populations. The transverse components of the fluctuations are discretized in the spectral space, starting from a scale , and then into 21 other shells at scales . We define the Elsässer perturbations at scale . The equations controlling the transport and dissipation of the fields are described in details in and Verdini et al. (2009). Non-linear interactions are operated in triads following the model of . Dissipation in Shell-Atm is obtained via explicit resistivity and viscosity. We set the magnetic Prandtl number to unity and the Lundquist number lies between in the chromosphere and in the corona. We write the heating per unit mass created by the cascade:Transverse motions’ amplitude is forced at the base of the domain in the chromosphere over three shells, starting at the scale . The forcing uses random fluctuation phases that are reset over a correlation time (see ). This is akin to setting a parallel frequency for the Alfvénic perturbations. is chosen around a few hundreds seconds to follow observations of transverse motions in the corona (see the review of ).
We chose to maintain the amplitude at the lower boundary, around 100 km above the photosphere, and on both boundaries for coronal loops configurations, using the following condition:This makes the inner boundary condition partially reflective, with a reflection coefficient being a function of time and of the balance of inward and outward wave populations. Imposing a fully reflecting inner boundary conditions leads indeed to an increase of the base velocity perturbations up to unrealistic values.
Next, we definethe energy flux of the perturbations at any given location in the domain. For both coronal holes and loops cases, we have the total energy injected per second and surface unit at a given timewhere the subscripts 0 and 1 denote the bottom and top of the computational domain respectively. In a statistical steady state, we expect to have the losses in the volume compensating for the input energy, i.e.,:whereis the integrated turbulent heating in the domain and W is the work of the force exerted by the perturbations on the solar wind flow:The component of the ponderomotive acceleration along field lines can be written (; ; ; ; ):where is the time averaged electric field due to the perturbations. Equation (10) also assumes that the ion cyclotron frequency is much larger than the wave frequency, which is easily verified in the inner heliosphere for typical Alfvén wave spectra and makes the force mass independent. In the remainder of the text, we will refer to Eq. 10 by ponderomotive force (per unit mass and unit volume), or ponderomotive acceleration indifferently.
As we consider one dimensional profile, we use the normalized area , which is unity at the base of all flux tubes, and we write Eq. 7 in .
Table 1 sums up all simulations made with the shell model for this work. For each profile geometry, we explore three sets of turbulence parameters referred as cases A, B, and C. In general, we chose the input parameters to create a turbulent heating close to what is observed in the (open) solar wind and as such used in the MHD described in Background Wind and Transition Region Profiles. We find that a transverse perturbation of km/s yields heating rate close to . — the minimum value to power the solar wind—in both coronal holes, which defines our reference set of turbulent inputs: cases A. Cases B use a higher forcing amplitude km/s, which logically results in a higher heating rate. Both cases A and B have the same injection scale km, which is the largest scale introduced in the system and corresponds to the size of supergranules. Cases C are set with a injection scale ten times lower, close to the size of granules. To get the heating in the coronal holes at the right order of magnitude, we had to decrease the input velocity perturbations to km/s. All cases have the same correlation time s.
TABLE 1
| Cases | (km/s) | (s) (s) | (km) | ( cgs) | ( cgs) |
|---|---|---|---|---|---|
| CH1 A | 7 | 600 | 35,000 | 0.84 | 0.49 |
| CH2 A | 7 | 600 | 35,000 | 1.8 | 1.0 |
| LO1 A | 7 | 600 | 35,000 | 20 | - |
| LO2 A | 7 | 600 | 35,000 | 10 | - |
| CH1 B | 10 | 600 | 35,000 | 3.1 | 1.6 |
| CH2 B | 10 | 600 | 35,000 | 4.6 | 2.9 |
| LO1 B | 10 | 600 | 35,000 | 26 | - |
| LO2 B | 10 | 600 | 35,000 | 21 | - |
| CH1 C | 3 | 600 | 3500 | 1.2 | 0.081 |
| CH2 C | 3 | 600 | 3500 | 1.2 | 0.092 |
| LO1 C | 3 | 600 | 3500 | 2.6 | - |
| LO2 C | 3 | 600 | 3500 | 2.3 | - |
Parameters and outputs for the Shell-Atm simulations.
In Figure 3, we show some properties of the solutions of the shell model runs A for km/s, s, and km. As the transverse perturbations are forced over three shells, the actual rms amplitude of the forcing is km/s. These parameters are close to the one used in Verdini et al. (2019) for a typical coronal hole solution. The results shown in Figure 3 and in Table 1 are obtained with a time average of the solution between 150,000 and 200,000 s, which represents around ten Alfvén crossing times for the coronal holes profiles and more for the coronal loops. All simulations have reached a pseudo-steady state during the period. The top panel shows the outward and inward rms fluctuations as a function of the curvilinear abscissa s along the profile. For coronal loops, we show in a logit scale, where L is the total length of the loop, to emphasize the behaviour in the chromosphere and the transition region at both boundaries.
FIGURE 3
In coronal holes simulations, one striking feature is the regime difference between a mostly balanced turbulence () below the transition region and an imbalanced turbulence beyond. This can be seen also in the bottom panel of Figure 3, where the inward and outward spectra are very comparable at (in blue), with a cascade covering four orders of magnitude and a spectral slope close to the usual Kolmogorov index . Higher in the corona, the outward wave dominates clearly while keeping a slope, while the inward component has a flatter spectrum with slope close to . In coronal loops, both populations are everywhere well represented and spectra seems perhaps closer to slopes. In the middle panel of Figure 3, we show the dissipation profile computed by the shell model. The plain black line is the true heating, while the gray line corresponds to the so-called phenomenological heating, a proxy often used in large scale fluid models (see, e.g. ; Verdini and Velli, 2007; ; Shoda et al., 2018a; ). It can be written:As said earlier, the total heating obtained in coronal holes solution is close to the input energy of the global MHD simulations, itself chosen to provide enough power in the solar wind (see Table 1). Nevertheless, the total heating increases as we go from a typical polar coronal hole, to a slow denser wind around the streamer and to coronal loops. This is expected, and increased heating in the low corona is one explanation for the denser slower wind, which originates somewhere close to the loops. Interestingly, we see that the phenomenological heating is in general an overestimation of the true heating in the shell model (as already noted in Verdini et al., 2019), except around the transition region. It does, however, seems to be a reasonable estimate in imbalanced turbulence region, i.e., in the open coronal wind regions.
The properties of cases B are very similar to the one of cases A, only with larger wave amplitudes and larger heating rates. For cases C, the wave amplitude is logically smaller, and the cascade extent is also shorter as we remain with a fixed value of the viscosity/resistivity . In all cases, the dissipation starts to be significant for. .
Ponderomotive Force and First Ionization Potential Fractionation
The main objective of this work is to compare the effect of the turbulence properties on the FIP fractionation for different background configurations. Among the key parameter of the current FIP models is the ponderomotive force, or the wave pressure exerted by the perturbations on the background flow. Figure 4 shows the ponderomotive force obtained in all Table 1 cases, close to the inner boundary, in the chromosphere and transition region. Let us first study the cases A of Table 1, shown in Figure 4 in the top panel. The profile of has in all cases a similar shape, with a slow decrease and a peak located around the transition region. The TR peak, which is responsible for most of the FIP fractionation in the model that follows, has an interesting ordering. The LO2 profile, has usually the highest averaged ponderomotive acceleration and a higher peak, but the maximum of the LO1 is usually of the order of the coronal holes configuration peaks. Moreover, it seems that the amplitude of the peak is a growing function of the height of the TR.
FIGURE 4
For cases B, we observe for the ponderomotive acceleration essentially a shift up of all the curves, conserving the hierarchy of the previous cases as a function of the geometry. For most cases A and B, the strength of the ponderomotive force is higher than the opposing gravity of the Sun, around , shown in black in Figure 4. The amplitude of the ponderomotive force is thus significant, especially at the transition region and Alfvén waves can have an influence on the coronal abundances. For cases C, however, the first striking feature is the net decrease of the ponderomotive acceleration below the Sun’s gravity pull. This means that although the turbulent heating is perfectly compatible with what is necessary to power the solar wind (compare especially case CH1 A with CH1 C), this set of parameter will likely not create a low-FIP bias through the ponderomotive acceleration. In these last cases, the turbulence injection scale is ten times smaller, i.e., the size of large granules. The initial amplitude is logically lower as the injection is made later in the cascade, which can explain the lower gradient of magnetic fluctuations.
We now compute the FIP fractionation created by the ponderomotive force in our models. We use for this the derivation of , that reads:Here, is the density of a given species j, the ionization fraction andwhere and are the collision frequencies of the species j with ions and neutrals respectively. We use the formulation of Schwadron et al. (1999), for the collision frequencies. The ionization fraction of each heavy ion is computed through an interpolation of Saha equilibria (; Saha, 1921) for K and the CHIANTI database and the ChiantiPy interface () for K. Finally, is the (quadratic) sum of the thermal speed and of a mass independent turbulent velocity , which represent a wave turbulence heating and is essential to avoid a mass dependant fractionation (see Schwadron et al., 1999).
Figure 5 shows the resulting FIP biases obtained with Eq. 12 and a constant km/s. We represent the minor ion density ratio between the low corona (7000 km above the photosphere) and the chromosphere (700 km above the photosphere), taking the Oxygen as a reference. The top and middle panel of Figure 5 shows a clear low-FIP bias in cases A and B. Elements Fe, Mg, and Si, show bias up to a factor 10 or more compared with Oxygen, in the red case, i.e., the smallest coronal loop. The profile of the FIP bias follows the hierarchy of the ponderomotive force amplitude shown in Figure 4. For CH1 and LO1 cases, the low-FIP bias is weak and roughly similar, while it is generally higher for the CH2 configuration. Finally, as intuited, cases C do not show any clear FIP biases. This is a direct result of the much weaker ponderomotive force obtained with these simulation parameters, and the heavy ion densities are the same in coronal holes and coronal loops.
FIGURE 5
In Figure 6, we repeat the same exercise, but choosing a turbulent velocityWe only sum the shell components for , where dissipation is significant. The average amplitude of is around 8 km/s for cases A, 12 km/s for cases B, and 5 km/s for cases C. The choice of is crucial to get FIP bias comparable with observations. It needs to be large enough to avoid the gravitational settling, and thus the mass dependence of the abundances in the corona. However, if is too large, it will kill the low-FIP bias. The low-FIP bias of cases A is thus slightly enhanced, because the average turbulent velocity given by Eq. 14 is lower than the value used in Figure 5. The FIP bias of case B is quenched and cases C start to show enrichment in low-FIP elements. Moreover, the hierarchy between the different magnetic configurations is much less visible in Figure 6. CH2 cases have for instance, properties very similar to LO2 cases. This shows that a very careful treatment of the turbulence is important to build FIP fractionation models.
FIGURE 6
Discussion
In this work, we have combined several models and tools to assess the efficiency of the Alfvénic turbulence in enriching abundances of low-FIP elements in the low corona. First, we used a global MHD model to get realistic profiles of the Alfvén speed and velocity gradients in both coronal holes and coronal loops, starting from the low chromosphere. They are indeed fundamental in the amplification, reflection and non-linear interaction of propagating Alfvén waves. Then, we took advantage of the Shell-Atm code to compute the interaction between counter-propagating Alfvén waves, the turbulent cascade and the dissipation. Our approach certainly lacks self-consistency, as, for instance, the heating obtained with the shell model is significantly larger in coronal loops than in coronal holes (at least for cases A and B), while the MHD simulation has a much more homogeneous heating. It is, however, the first time that both coronal loops and coronal holes are treated with a shell turbulence model. This model solves the fully non-linear incompressible Alfvén wave equations, including the cascading process and the frequency dependant reflections and interactions of Alfvén wave populations. This is an improvement in comparison with analytical non-WKB approach of previous works (see ; ; ). Moreover, the dissipation is treated physically at very high Lundquist numbers, which cannot be achieved in direct numerical simulations (see, e.g., ).
A few comments on the applicability of incompressible MHD to this problem are in order. Reduced MHD equations are derived assuming small gradients in the guide field (or even constant , see Strauss, 1976; Zank and Matthaeus, 1992; Zank and Matthaeus, 1993; ) as well as in the density, conditions that are not verified in the upper chromosphere and the transition region. The excluded nonlinear couplings in the parallel and perpendicular directions (essentially compressible interactions) are responsible for the coupling of Alfvén waves to slow and fast magneto-acoustic modes, leading to steepening, shock formation and dissipation of slow magneto-acoustic modes and fast modes, and the refraction of the latter downwards back to the chromosphere for large perpendicular wavenumbers. The parametric decay instability (PDI), involving compressible processes, is also excluded. The latter plays a role in the formation of the turbulent spectrum and the balance of inward and outward Alfvén wave population, but is effective most in the lower beta regions higher up in the corona (see Tenerani and Velli, 2013; Shoda et al., 2018b; ). As noted in Verdini et al. (2019), Shoda et al. (2019), density fluctuations created by the PDI increase the inward Alfvén waves amplitude and eventually increase the turbulent heating in the open wind regions.
Wave steepening, shocks and mode conversion, which do create a cascade in the direction parallel to the magnetic field, have typical spectra and are thus likely to be a secondary dissipation channel in comparison with the dissipation in the perpendicular plane, with spectra close to the Kolmogorov and Iroshnikov-Kraichnan phenomenology (see Figure 3). Alfvén wave coupling with magneto-acoustic modes may accelerate the perpendicular cascade, but because we run the simulation until a pseudo steady-state is reached, this is likely not a strong limitation of our approach. Numerous other modelers have made similar approximations to our own, as for example the studies by van Ballegooijen et al. (2011), , van Ballegooijen and Asgari-Targhi (2016), van Ballegooijen and Asgari-Targhi (2017), and . This study has no intention to be exhaustive, but we do examine a significant input parameter space in order to extract meaningful insights on the problem.
Our study shows that, assuming that turbulence is the dominant factor in the coronal heating and solar wind acceleration, a ponderomotive force can appear in the chromosphere and the transition region, and can be strong enough to create a low-FIP bias. This depends however on the turbulence parameters. Injecting energy at the scales of super granules provide the wave amplitude necessary for a low-FIP bias comparable with observations (see, e.g. ). The force is related to the a amplification of the waves in the corona and the strong gradient that appears at the transition region. Nevertheless, if the energy is injected at the scale of granules, the resulting ponderomotive acceleration seems too weak to explain the observations. Recent works have debated of the right injection scale parameters (van Ballegooijen and Asgari-Targhi, 2017; ), and the amplitude of the ponderomotive force could be a way of constraining solar wind turbulence models.
A second result of this study is that the low-FIP bias is not exclusive to coronal loops. Interestingly, we obtain significant low-FIP bias for open field configurations along the streamer (CH2). We also have usually little FIP bias for the LO1 loop, i.e., the loop at the very edge of the helmet streamer. This result might seem in opposition with previous works, such as the one of , . Their result is indeed the consequence of resonant Alfvén waves in the loops, likely triggered by reconnection, which amplify significantly the ponderomotive acceleration. Injecting transverse motions from the inner boundary, we do not observe such resonances, and further studies are needed to compare the relative importance of different triggers. We do show however that resonances are not needed to obtain a low-FIP bias in coronal loops and open slow wind regions. This study thus questions, anyhow, the necessity of interchange reconnection to explain the composition of the slow solar wind.
Finally, it is important to stress that our modeling of the solar chromosphere is very simple. We assume, following , and writing Eq. 12, that proton drag, collisions and turbulent mixing, are exactly compensating for the Sun’s gravity pull in the chromosphere. The actual balance of all these processes remains very hard to estimate. Shocks and compressible fluctuations are also believed to be significant in the chromosphere and even in the transition region (see , for a recent review) and are not accounted for in the shell simulations. Ongoing works are dedicated to build a compressible kinetic-fluid model, including heavy ions populations with various charge state and their interactions with protons and waves.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
VR is the main developer of the multi-D MHD model and performed all simulations. AV, ÉB, and MV are the main developers of the SHELL-ATM code. All authors have contributed to the data analysis, text editing and discussions.
Funding
This research was funded by the ERC SLOW SOURCE project (SLOW SOURCE–DLV-819189). MV was supported by the NASA HERMES Drive Science Center grant No. 80NSSC20K0604. We thank M. Laming for very useful discussions. Numerical computations were performed using GENCI grants numbers A0080410133 and A0070410293. This study has made use of the NASA Astrophysics Data System.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AlazrakiG.CouturierP. (1971). Solar wind accejeration caused by the gradient of alfven wave pressure. Astron. Astrophys.13, 380–389.
2
AntiochosS. K.LinkerJ. A.LionelloR.MikićZ.TitovV.ZurbuchenT. H. (2012). The structure and dynamics of the corona-heliosphere connection. Space Sci. Rev.172, 169–185. 10.1007/s11214-011-9795-7
3
BelcherJ. W. (1971). Alfvénic wave pressures and the solar wind. Astrophys. J.168, 509. 10.1086/151105
4
BelcherJ. W.DavisL.Jr. (1971). Large-amplitude Alfvén waves in the interplanetary medium, 2. J. Geophys. Res.76 (16), 3534–3563. 10.1029/ja076i016p03534
5
BochslerP.GeissJ.KunzS. (1986). Abundances of carbon, oxygen, and neon in the solar wind during the period from August 1978 to June 1982. Sol. Phys.103 (1), 177–201. 10.1007/bf00154867
6
BrunoR.CarboneV. (2013). The Solar wind as a turbulence laboratory. Living Rev. Sol. Phys.10, 2. 10.12942/lrsp-2013-2
7
BuchlinE.VelliM. (2007). Shell models of RMHD turbulence and the heating of solar coronal loops. Astrophys. J.662 (1), 701–714. 10.1086/512765
8
CarlssonM.De PontieuB.HansteenV. H. (2019). New view of the solar chromosphere. Annu. Rev. Astron. Astrophys.57 (1), 189–226. 10.1146/annurev-astro-081817-052044
9
ChandranB. D. G.HollwegJ. V. (2009). Alfvén wave reflection and turbulent heating in the solar wind from 1 solar radius to 1 Au: an analytical treatment. Astrophys. J.707 (2), 1659–1667. 10.1088/0004-637x/707/2/1659
10
ChandranB. D. G.PerezJ. C. (2019). Reflection-driven magnetohydrodynamic turbulence in the solar atmosphere and solar wind. J. Plasma Phys.85 (4), 905850409. 10.1017/s0022377819000540
11
DahlburgR. B.LamingJ. M.TaylorB. D.ObenschainK. (2016). Ponderomotive acceleration in coronal loops. Astrophys. J.831 (2), 160. 10.3847/0004-637x/831/2/160
12
DereK. P.LandiE.MasonH. E.Monsignori FossiB. C.YoungP. R. (1997). CHIANTI—an atomic database for emission lines. Astron. Astrophys. Suppl. Ser.125 (1), 149–173. 10.1051/aas:1997368
13
DeRosaM. L.BrunA. S.HoeksemaJ. T. (2012). Solar magnetic field reversals and the role of dynamo families. Astrophys. J.757 (1), 96. 10.1088/0004-637x/757/1/96
14
DmitrukP.MatthaeusW. H.MilanoL. J.OughtonS.ZankG. P.MullanD. J. (2002). Coronal heating distribution due to low‐frequency, wave‐driven turbulence. Astrophys. J.575 (1), 571–577. 10.1086/341188
15
FeldmanU. (1992). Elemental abundances in the upper solar atmosphere. Phys. Scripta46 (3), 202–220. 10.1088/0031-8949/46/3/002
16
GeissJ.GloecklerG.von SteigerR. (1995). Origin of the solar wind from composition Data. Space Sci. Rev.72 (1–2), 49–60. 10.1007/978-94-011-0167-7_8
17
GiulianiP.CarboneV. (1998). A note on shell models for MHD Turbulence. Europhys. Lett.43 (5), 527–532. 10.1209/epl/i1998-00386-y
18
GloecklerG.GeissJ. (1989). The abundances of elements and isotopes in the solar wind.AIP Conf. Proc. 183 (1), 49–71. 10.1063/1.37985
19
JacquesS. A. (1977). Momentum and energy transport by waves in the solar atmosphere and solar wind. Astrophys. J.215, 942. 10.1086/155430
20
LamingJ. M. (2004). A unified picture of the first ionization potential and inverse first ionization potential effects. Astrophys. J.614 (2), 1063–1072. 10.1086/423780
21
LamingJ. M. (2009). Non-WKB models of the first ionization potential effect: implications for solar coronal heating and the coronal helium and neon abundances. Astrophys. J.695 (2), 954–969. 10.1088/0004-637x/695/2/954
22
LamingJ. M. (2012). Non-WKB models of the first ionization potential effect: the role of slow mode waves. Astrophys. J.744 (2), 115. 10.1088/0004-637x/744/2/115
23
LamingJ. M. (2015). The FIP and Inverse FIP Effects in Solar and Stellar Coronae. Living Rev. Sol. Phys.12, 2. 10.1007/lrsp-2015-2
24
LeerE.HolzerT. E.FlåT. (1982). Acceleration of the solar wind. Space Sci. Rev.33, 161–200. 10.1007/BF00213253
25
LionelloR.LinkerJ. A.MikićZ. (2001). Including the transition region in models of the large‐scale solar corona. Astrophys. J.546 (1), 542–551. 10.1086/318254
26
LitwinC.RosnerR. (1998). Coronal scale-height enhancement by magnetohydrodynamic waves. Astrophys. J.506 (2), L143–L146. 10.1086/311656
27
MarschE.von SteigerR.BochslerP. (1995). Element fractionation by diffusion in the solar chromosphere. Astron. Astrophys.301, 261–276.
28
MeyerJ.-P. (1985a). The baseline composition of solar energetic particles. Astrophys. J. Suppl. Ser.57, 173–204.
29
MeyerJ.-P. (1985b). The baseline composition of solar energetic particles. Astrophys. J. Suppl. Ser.57, 151–171. 10.1086/191000
30
MignoneA.BodoG.MassagliaS.MatsakosT.TesileanuO.ZanniC.et al (2007). PLUTO: a numerical code for computational Astrophysics. Astrophys. J. Suppl. Ser.170 (1), 228–242. 10.1086/513316
31
NakariakovV. M.VerwichteE. (2005). Coronal Waves and Oscillations. Living Rev. Sol. Phys.2, 3. 10.12942/lrsp-2005-3
32
NigroG.MalaraF.CarboneV.VeltriP. (2004). Nanoflares and MHD turbulence in coronal loops: a hybrid shell model. Phys. Rev. Lett.92 (19), 194501. 10.1103/physrevlett.92.194501
33
OughtonS.MatthaeusW. H.DmitrukP. (2017). Reduced MHD in astrophysical applications: two-dimensional or three-dimensional?Astrophys. J.839 (1), 2. 10.3847/1538-4357/aa67e2
34
PerezJ. C.ChandranB. D. G. (2013). Direct numerical simulations of reflection-driven, reduced magnetohydrodynamic turbulence from the Sun to the Alfvén critical point. Astrophys. J.776 (2), 124. 10.1088/0004-637x/776/2/124
35
RévilleV.TeneraniA.VelliM. (2018). Parametric decay and the origin of the low-frequency alfvénic spectrum of the solar wind. Astrophys. J.866 (1), 38. 10.3847/1538-4357/aadb8f
36
RévilleV.VelliM.PanasencoO.TeneraniA.ShiC.BadmanS. T.et al (2020a). The role of Alfvén wave dynamics on the large-scale properties of the solar wind: comparing an MHD simulation with parker solar probe E1 Data. Astrophy. J. Suppl. Ser.246 (2), 24. 10.3847/1538-4365/ab4fef
37
RévilleV.VelliM.RouillardA. P.LavraudB.TeneraniA.ShiC.et al (2020b). Tearing instability and periodic density perturbations in the slow solar wind. Astrophy. J. Suppl. Ser.895 (1), L20. 10.3847/2041-8213/ab911d
38
SahaM. N. (1920). Ionisation in the solar chromosphere. Nature105 (2634), 232–233. 10.1038/105232b0
39
SahaM. N. (1921). On a physical theory of stellar spectra. Proc. Royal Soc. London Ser. A99, 135–153. 10.1098/rspa.1921.0029
40
SchwadronN. A.FiskL. A.ZurbuchenT. H. (1999). Elemental fractionation in the slow solar wind. Astrophy. J.521 (2), 859–867. 10.1086/307575
41
ShodaM.SuzukiT. K.Asgari-TarghiM.YokoyamaT. (2019). Three-dimensional simulation of the fast solar wind driven by compressible magnetohydrodynamic turbulence. Astrophy. J.880 (1), L2. 10.3847/2041-8213/ab2b45
42
ShodaM.YokoyamaT.SuzukiT. K. (2018a). A self-consistent model of the coronal heating and solar wind acceleration including compressible and incompressible heating processes. Astrophy. J.853 (2), 190. 10.3847/1538-4357/aaa3e1
43
ShodaM.YokoyamaT.SuzukiT. K. (2018b). Frequency-dependent Alfvén-wave propagation in the solar wind: onset and suppression of parametric decay instability. Astrophy. J.860 (1), 17. 10.3847/1538-4357/aac218
44
StraussH. R. (1976). Nonlinear, three-dimensional magnetohydrodynamics of noncircular tokamaks. Phys. Fluids19 (1), 134. 10.1063/1.861310
45
TeneraniA.VelliM. (2013). Parametric decay of radial Alfvén waves in the expanding accelerating solar wind. J. Geophys. Res. Space Phys.118 (12), 7507–7516. 10.1002/2013ja019293
46
TuC.-Y.MarschE. (1995). MHD structures, waves and turbulence in the solar wind: observations and theories. Space Sci. Rev.73, 1–210. 10.1007/BF00748891
47
van BallegooijenA. A.Asgari-TarghiM. (2017). Direct and inverse cascades in the acceleration region of the fast solar wind. Astrophys. J.835 (1), 10. 10.3847/1538-4357/835/1/10
48
van BallegooijenA. A.Asgari-TarghiM.CranmerS. R.DeLucaE. E. (2011). Heating of the solar chromosphere and corona by Alfvén wave turbulence. Astrophy. J.736 (1), 3. 10.1088/0004-637x/736/1/3
49
van BallegooijenA. A.Asgari-TarghiM. (2016). Heating and acceleration of the fast solar wind by Alfvén wave turbulence. Astrophy. J.821 (2), 106. 10.3847/0004-637x/821/2/106
50
VelliM.GrappinR.MangeneyA. (1989). Turbulent cascade of incompressible unidirectional Alfvén waves in the interplanetary medium. Phys. Rev. Lett.63 (17), 1807–1810. 10.1103/physrevlett.63.1807
51
VerdiniA.GrappinR.Montagud-CampsV. (2019). Turbulent Heating in the Accelerating Region Using a Multishell Model. Solar Phy.294, 65. 10.1007/s11207-019-1458-y
52
VerdiniA.VelliM. (2007). Alfvén waves and turbulence in the solar atmosphere and solar wind. Astrophy. J.662 (1), 669–676. 10.1086/510710
53
VerdiniA.VelliM.BuchlinE. (2009). Turbulence in the sub-alfvénic solar wind driven by reflection of low-frequency Alfvén waves. Astrophy. J.700 (1), L39–L42. 10.1088/0004-637x/700/1/l39
54
WoodB. E.LinskyJ. L. (2010). Resolving the ξ BOO binary withchandra, and revealing the spectral type dependence OF the coronal “FIP effect”. Astrophy. J.717 (2), 1279–1290. 10.1088/0004-637x/717/2/1279
55
ZankG. P.MatthaeusW. H. (1992). The equations of reduced magnetohydrodynamics. J. Plasma Phys.48 (1), 85–100. 10.1017/s002237780001638x
56
ZankG. P.MatthaeusW. H. (1993). Nearly incompressible fluids. II: magnetohydrodynamics, turbulence, and waves. Phys. Fluid. Fluid Dynam.5 (1), 257–273. 10.1063/1.858780
57
ZhouY.MatthaeusW. H. (1989). Non-WKB evolution of solar wind fluctuations: a turbulence modeling approach. Geophys. Res. Lett.16 (7), 755–758. 10.1029/gl016i007p00755
Summary
Keywords
turbulence, slow wind, transition region, chromosphere, Alfvén waves
Citation
Réville V, Rouillard AP, Velli M, Verdini A, Buchlin É, Lavarra M and Poirier N (2021) Investigating the Origin of the First Ionization Potential Effect With a Shell Turbulence Model. Front. Astron. Space Sci. 8:619463. doi: 10.3389/fspas.2021.619463
Received
20 October 2020
Accepted
04 January 2021
Published
08 February 2021
Volume
8 - 2021
Edited by
Luca Sorriso-Valvo, Institute for Space Physics, Sweden
Reviewed by
Mahboubeh Asgari-Targhi, Harvard University, United States
Giuseppina Nigro, Università Della Calabria, Italy
Updates
Copyright
© 2021 Réville, Rouillard, Velli, Verdini, Buchlin, Lavarra and Poirier.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Victor Réville, victor.reville@irap.omp.eu
This article was submitted to Space Physics, a section of the journal Frontiers in Astronomy and Space Sciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.