Observations of solar wind turbulence indicate the existence of multi-scale pressure-balanced structures (PBSs) in the solar wind. In this work, we conduct a numerical simulation to investigate multi-scale PBSs and in particular their formation in compressive MHD turbulence. By the use of a higher order Godunov code Athena,a driven compressible turbulence with an imposed uniform guide field is simulated. The simulation results show that both the magnetic pressure and the thermal pressure exhibit a turbulent spectrum with a Kolmogorov-like power law, and that in many regions of the simulation domain they are anti-correlated. The computed wavelet cross-coherence spectrum of the magnetic pressure and the thermal pressure, as well as their space series, indicate the existence of multi-scale PBSs, with the small PBSs being embedded in the large ones. These multi-scale PBSs are likely to be related with the highly oblique-propagating slow-mode waves, as the traced multi-scale PBS is found to be traveling in a certain direction at a speed consistent with that predicted theoretically for a slow-mode wave propagating in the same direction.

In the solar atmosphere, the jets are ubiquitous and found to be at various spatia-temporal scales. They are significant to understand energy and mass transport in the solar atmosphere. Recently, the high-speed transition region jets are reported from the observation. Here we conduct a numerical simulation to investigate the mechanism in their formation. Driven by the supergranular convection motion, the magnetic reconnection between the magnetic loop and the background open flux occurring in the transition region is simulated with a two-dimensional magnetohydrodynamics model. The simulation results show that not only a fast hot jet, much resemble the found transition region jets, but also a adjacent slow cool jet, mostly like classical spicules, is launched. The force analysis shows that the fast hot jet is continually driven by the Lorentz force around the reconnection region, while the slow cool jet is induced by an initial kick through the Lorentz force associated with the emerging magnetic flux. Also, the features of the driven jets change with the amount of the emerging magnetic flux, giving the varieties of both jets. These results will inspire our understanding of the formation of the prevalence of both the fast hot jet and slow cool jet from the solar transition region and chromosphere.

In the solar atmosphere, jets are prevalent and they are significant for the mass and energy transport. Here we conduct numerical simulations to investigate the mass and energy contributions of the recently observed high-speed jets to the solar wind. With a one-dimensional hydrodynamic solar wind model, the time-dependent pulses are imposed at the bottom to simulate the jets. The simulation results show that without other energy source, the injected plasmas are accelerated effectively to be a transonic wind with a substantial mass flux. The rapid acceleration occurs close to the Sun, and the resulting asymptotic speed, number density at 0.3 AU, as well as mass flux normalized to 1 AU are compatible with in situ observations. As a result of the high speed, the imposed pulses generate a train of shocks traveling upward. By tracing the motions of the injected plasma, it is found that these shocks heat and accelerate the injected plasmas successively step by step to push them upward and eventually allow them to escape. The parametric studies show that increasing the speed of the imposed pulses or their temperature gives a considerably faster, and hotter solar wind, while increasing their number density or decreasing their recurring period only bring a denser solar wind. These studies provide a possibility that the ubiquitous high-speed jets are a substantial mass and energy contributions to the solar wind.

Quasi-periodic disturbances of emission-line parameters are frequently observed in the corona. These disturbances propagate upward along the magnetic field with speeds $\sim100~\rm{km~s}^{-1}$. This phenomenon has been interpreted as evidence of the propagation of slow magnetosonic waves or argued to be signature of the intermittent outflows superposed on the background plasmas. Here we aim to present a new "wave + flow" model to interpret these observations. In our scenario, the oscillatory motion is a slow mode wave, and the flow is associated with a beam created by the wave-particle interaction owing to Landau resonance. With the help of a Vlasov model, we simulate the propagation of the slow mode wave and the generation of the beam flow. We find that weak periodic beam flows can be generated owing to Landau resonance in the solar corona, and the phase with strongest blueward asymmetry is ahead of that with strongest blueshift by about 1/4 period. We also find that the slow wave damps to the level of 1/e after the transit time of two wave periods, owing to Landau damping and Coulomb collisions in our simulation. This damping time scale is similar to that resulting from thermal-conduction in the magnetohydrodynamics regime. The beam flow is weakened/attenuated with increasing wave period and decreasing wave amplitude since Coulomb collision becomes more and more dominant over the wave action. We suggest that this "wave + flow" kinetic model provides an alternative explanation for the observed quasi-periodic propagating perturbations in various parameters in the solar corona.

Intensive studies have been conducted to understand the anisotropy of solar wind turbulence. However, the anisotropy of Els\"asser variables ($\textbf{Z}^\pm$) in 2D wave-vector space has yet to be investigated. Here we first verify the transformation based on the projection-slice theorem between the power spectral density PSD$_{2D}(k_\parallel,k_\perp )$ and the spatial correlation function CF$_{2D} (r_\parallel,r_\perp )$. Based on the application of the transformation to the magnetic field and the particle measurements from the WIND spacecraft, we investigate the spectral anisotropy of Els\"asser variables ($\textbf{Z}^\pm$), and the distribution of residual energy E$_{R}$, Alfv\'en ratio R$_{A}$ and Els\"asser ratio R$_{E}$ in the $(k_\parallel,k_\perp)$ space. The spectra PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{B}$, $\textbf{V}$, and $\textbf{Z}_{major}$ (the larger of $\textbf{Z}^\pm$) show a similar pattern that PSD$_{2D}(k_\parallel,k_\perp )$ is mainly distributed along a ridge inclined toward the $k_\perp$ axis. This is probably the signature of the oblique Alfv\'enic fluctuations propagating outwardly. Unlike those of $\textbf{B}$, $\textbf{V}$, and $\textbf{Z}_{major}$, the spectrum PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{Z}_{minor}$ is distributed mainly along the $k_\perp$ axis. Close to the $k_\perp$ axis, $\left| {E}_{R}\right|$ becomes larger while R$_{A}$ becomes smaller, suggesting that the dominance of magnetic energy over kinetic energy becomes more significant at small $k_\parallel$. R$_{E}$ is larger at small $k_\parallel$, implying that PSD$_{2D}(k_\parallel,k_\perp )$ of $\textbf{Z}_{minor}$ is more concentrated along the $k_\perp$ direction as compared to that of $\textbf{Z}_{major}$. The residual energy condensate at small $k_\parallel$ is consistent with simulation results in which E$_{R}$ is spontaneously generated by Alfv\'en wave interaction.

Magnetohydronamic turbulence is believed to play a crucial role in heating the laboratorial, space, and astrophysical plasmas. However, the precise connection between the turbulent fluctuations and the particle kinetics has not yet been established. Here we present clear evidence of plasma turbulence heating based on diagnosed wave features and proton velocity distributions from solar wind measurements by the Wind spacecraft. For the first time, we can report the simultaneous observation of counter-propagating magnetohydrodynamic waves in the solar wind turbulence. Different from the traditional paradigm with counter-propagating Alfv\'en waves, anti-sunward Alfv\'en waves (AWs) are encountered by sunward slow magnetosonic waves (SMWs) in this new type of solar wind compressible turbulence. The counter-propagating AWs and SWs correspond respectively to the dominant and sub-dominant populations of the imbalanced Els\"asser variables. Nonlinear interactions between the AWs and SMWs are inferred from the non-orthogonality between the possible oscillation direction of one wave and the possible propagation direction of the other. The associated protons are revealed to exhibit bi-directional asymmetric beams in their velocity distributions: sunward beams appearing in short and narrow patterns and anti-sunward broad extended tails. It is suggested that multiple types of wave-particle interactions, i.e., cyclotron and Landau resonances with AWs and SMWs at kinetic scales, are taking place to jointly heat the protons perpendicularly and parallel.

Transient brightenings in the transition region of the Sun have been studied for decades and are usually related to magnetic reconnection. Recently, absorption features due to chromospheric lines have been identified in transition region emission lines raising the question of the thermal stratification during such reconnection events. We analyse data from the Interface Region Imaging Spectrograph (IRIS) in an emerging active region. Here the spectral profiles show clear self-absorption features in the transition region lines of Si\,{\sc{iv}}. While some indications existed that opacity effects might play some role in strong transition region lines, self-absorption has not been observed before. We show why previous instruments could not observe such self-absorption features, and discuss some implications of this observation for the corresponding structure of reconnection events in the atmosphere. Based on this we speculate that a range of phenomena, such as explosive events, blinkers or Ellerman bombs, are just different aspects of the same reconnection event occurring at different heights in the atmosphere.

Measurements of solar wind turbulence reveal the ubiquity of discontinuities. In this study, we investigate how the discontinuities, especially rotational discontinuities (RDs), are formed in magnetohydrodynamic (MHD) turbulence. In a simulation of the decaying compressive three-dimensional (3-D) MHD turbulence with an imposed uniform background magnetic field, we detect RDs with sharp field rotations and little variations of magnetic field intensity as well as mass density. At the same time, in the de Hoffman-Teller (HT) frame, the plasma velocity is nearly in agreement with the Alfv\'{e}n speed, and is field-aligned on both sides of the discontinuity. We take one of the identified RDs to analyze in details its 3-D structure and temporal evolution. By checking the magnetic field and plasma parameters, we find that the identified RD evolves from the steepening of the Alfv\'{e}n wave with moderate amplitude, and that steepening is caused by the nonuniformity of the Alfv\'{e}n speed in the ambient turbulence.

In solar wind, magnetohydrodynamic (MHD) discontinuities are ubiquitous and often found to be at the origin of turbulence intermittency. They may also play a key role in the turbulence dissipation and heating of the solar wind. The tangential (TD) and rotational (RD) discontinuities are the two most important types of discontinuities. Recently, the connection between turbulence intermittency and proton thermodynamics has been being investigated observationally. Here we present numerical results from three-dimensional MHD simulation with pressure anisotropy and define new methods to identify and to distinguish TDs and RDs. Three statistical results obtained about the relative occurrence rates and heating effects are highlighted: (1) RDs tend to take up the majority of the discontinuities along with time; (2) the thermal states embedding TDs tend to be associated with extreme plasma parameters or instabilities, while RDs do not; (3) TDs have a higher average T as well as perpendicular temperature $T_\perp$. The simulation shows that TDs and RDs evolve and contribute to solar wind heating differently. These results will inspire our understanding of the mechanisms that generate discontinuities and cause plasma heating.

Superhalo electrons appear to be continuously present in the interplanetary medium, even at very quiet times, with a power-law spectrum at energies above $\sim$2 keV. Here we numerically investigate the generation of superhalo electrons by magnetic reconnection in the solar wind source region, using the MHD and test particle simulations for both single X-line reconnection and multiple X-line reconnection. We find that the direct current electric field, produced in the magnetic reconnection region, can accelerate electrons from an initial thermal energy of T $\sim10^5$ K up to hundreds of keV. After acceleration, some of the accelerated electrons, together with the nascent solar wind flow driven by the reconnection, propagate upwards along the newly-opened magnetic field lines into the interplanetary space, while the rest move downwards into the lower atmosphere. Similar to the observed superhalo electrons at 1 AU, the flux of the upward-traveling accelerated electrons versus energy displays a power-law distribution at $\sim$ 2 $-$ 100 keV, $f(E) \sim E^{-\delta}$, with a $\delta$ of $\sim$ 1.5 $-$ 2.4. For single (multiple) X-line reconnection, the spectrum becomes harder (softer) as the anomalous resistivity parameter $\alpha$ (uniform resistivity $\eta$) increases. These modeling results suggest that the acceleration in the solar wind source region may contribute to superhalo electrons.

The Sun is the only star that we can spatially resolve and it can be regarded as a fundamental plasma laboratory of astrophysics. The solar transition region (TR), the layer between the solar chromosphere and corona, plays an important role in solar wind origin and coronal heating. Recent high-resolution observations made by SOHO, TRACE, and Hinode indicate that the TR is highly nonuniform and magnetically structured. Through a combination of spectroscopic observations and magnetic field extrapolations, the TR magnetic structures and plasma properties have been found to be different in coronal holes and in the quiet Sun. In active regions, the TR density and temperature structures also differ in sunspots and the surrounding plage regions. Although the TR is believed to be a dynamic layer, quasi-steady flows lasting from several hours to several days are often present in the quiet Sun, coronal holes, and active regions, indicating some kind of plasma circulation/convection in the TR and corona. The emission of hydrogen Lyman lines, which originates from the lower TR, has also been intensively investigated in the recent past. Observations show clearly that the flows and dynamics in the middle and upper TR can greatly modify the Lyman line profiles.

We present observations of the power spectral anisotropy in wave-vector space of solar wind turbulence, and study how it evolves in interplanetary space with increasing heliocentric distance. For this purpose we use magnetic field measurements made by the Helios-2 spacecraft at three positions between 0.29 and 0.9 AU. To derive the power spectral density (PSD) in $(k_\parallel, k_\bot)$-space based on single-satellite measurements is a challenging task not yet accomplished previously. Here we derive the spectrum $\rm{PSD}_{\rm{2D}}$($\rm{k}_\parallel$, $\rm{k}_\bot$) from the spatial correlation function $\rm{CF}_{\rm{2D}}(r_\parallel, r_\bot)$ by a transformation according to the projection-slice theorem. We find the so constructed PSDs to be distributed in k-space mainly along a ridge that is more inclined toward the $\rm{k}_\bot$ than $\rm{k}_\parallel$ axis, a new result which probably indicates preferential cascading of turbulent energy along the $\rm{k}_\bot$ direction. Furthermore, this ridge of the distribution is found to gradually get closer to the $\rm{k}_\bot$ axis, as the outer scale length of the turbulence becomes larger while the solar wind flows further away from the Sun. In the vicinity of the $\rm{k}_\parallel$ axis, there appears a minor spectral component that probably corresponds to quasi-parallel Alfv\'enic fluctuations. Their relative contribution to the total spectral density tends to decrease with radial distance. These findings suggest that solar wind turbulence undergoes an anisotropic cascade transporting most of its magnetic energy towards larger $\rm{k}_\bot$, and that the anisotropy in the inertial range is radially developing further at scales that are relatively far from the ever increasing outer scale.

The origin of the solar wind is one of the most important unresolved prob- lems in space and solar physics. We report here the first spectroscopic signatures of the nascent fast solar wind on the basis of observations made by the EUV Imaging Spectrometer (EIS) on Hinode in a polar coronal hole, in which patches of blueshift are clearly present on dopplergrams of coronal emission lines with a formation temperature of lg(T/K)=6.0. The corresponding upflow is associated with open field lines in the coronal hole, and seems to start in the solar transition region and becomes more prominent with increasing temperature. This temperature-dependent plasma outflow is interpreted as evidence of the nascent fast solar wind in the polar coronal hole. The patches with significant upflows are still isolated in the upper transition region but merge in the corona, in agreement with the scenario of solar wind outflow being guided by expanding magnetic funnels.

We investigate the physical meaning of the prominent blue shifts of Ne VIII, which is observed to be associated with quiet-Sun network junctions (boundary intersections), through data analyses combining force-free-field extrapolations with EUV spectroscopic observations. For a middle-latitude region, we reconstruct the magnetic funnel structure in a sub-region showing faint emission in EIT-Fe 195. This funnel appears to consist of several smaller funnels that originate from network lanes, expand with height and finally merge into a single wide open-field region. However, the large blue shifts of Ne VIII are generally not associated with open fields, but seem to be associated with the legs of closed magnetic loops. Moreover, in most cases significant upflows are found in both of the funnel-shaped loop legs. These quasi-steady upflows are regarded as signatures of mass supply to the coronal loops rather than the solar wind. Our observational result also reveals that in many cases the upflows in the upper transition region (TR) and the downflows in the middle TR are not fully cospatial. Based on these new observational results, we suggest different TR structures in coronal holes and in the quiet Sun.

We study the emission heights of the coronal bright points (BPs) above the photosphere in the bipolar magnetic loops that are apparently associated with them. As BPs are seen in projection against the disk their true emission heights are unknown. The correlation of the BP locations on the Fe XII radiance map from EIT with the magnetic field features (in particular neutral lines) was investigated in detail. The coronal magnetic field was determined by an extrapolation of the photospheric field to different altitudes above the disk. It was found that most BPs sit on or near a photospheric neutral line, but that the emission occurs at a height of about 5 Mm. Some BPs, while being seen in projection, still seem to coincide with neutral lines, although their emission takes place at heights of more than 10 Mm. Such coincidences almost disappear for emissions above 20 Mm. We also projected the upper segments of the 3-D magnetic field lines above different heights, respectively, on to the x-y plane. The shape of each BP was compared with the respective field-line segment nearby. This comparison suggests that most coronal BPs are actually located on the top of their associated magnetic loops. Finally, we calculated for each selected BP region the correlation coefficient between the Fe XII intensity enhancement and the horizontal component of the extrapolated magnetic field vector at the same x-y position in planes of different heights, respectively. We found that for almost all the BP regions we studied the correlation coefficient, with increasing height, increases to a maximal value and then decreases again. The height corresponding to this maximum was defined as the correlation height, which for most bright points was found to range below 20 Mm.

The radiance and Doppler-shift distributions across the solar network provide observational constraints of two-dimensional modeling of transition-region emission and flows in coronal funnels. Two different methods, dispersion plots and average-profile studies, were applied to investigate these distributions. In the dispersion plots, we divided the entire scanned region into a bright and a dark part according to an image of Fe xii; we plotted intensities and Doppler shifts in each bin as determined according to a filtered intensity of Si ii. We also studied the difference in height variations of the magnetic field as extrapolated from the MDI magnetogram, in and outside network. For the average-profile study, we selected 74 individual cases and derived the average profiles of intensities and Doppler shifts across the network. The dispersion plots reveal that the intensities of Si ii and C iv increase from network boundary to network center in both parts. However, the intensity of Ne viii shows different trends, namely increasing in the bright part and decreasing in the dark part. In both parts, the Doppler shift of C iv increases steadily from internetwork to network center. The average-profile study reveals that the intensities of the three lines all decline from the network center to internetwork region. The binned intensities of Si ii and Ne viii have a good correlation. We also find that the large blue shift of Ne viii does not coincide with large red shift of C iv. Our results suggest that the network structure is still prominent at the layer where Ne viii is formed in the quiet Sun, and that the magnetic structures expand more strongly in the dark part than in the bright part of this quiet Sun region.

Aims. The physical implication of large blue shift of Ne viii in the quiet Sun region is investigated in this paper. Methods. We compare the significant Ne viii blue shifts, which are visible as large blue patches on the Doppler-shift map of a middlelatitude quiet-Sun region observed by SUMER, with the coronal magnetic-field structures as reconstructed from a simultaneous photospheric magnetogram by means of a force-free-field extrapolation. Results. We show for the first time that coronal funnels also exist in the quiet Sun. The region studied contains several small funnels that originate from network lanes, expand with height and finally merge into a single wide open-field region. However, the large blue shifts of the Ne viii line are not generally associated with funnels. A comparison between the projections of coronal loops onto the solar x-y-plane and the Ne viii dopplergram indicates that there are some loops that reveal large Ne viii blue shifts in both legs, and some loops with upflow in one and downflow in the other leg. Conclusions. Our results suggest that strong plasma outflow, which can be traced by large Ne viii blue shift, is not necessarily associated with the solar wind originating in coronal funnels but appears to be a signature of mass supply to coronal loops. Under the assumption that the measured Doppler shift of the Ne viii line represents the real outflow velocity of the neon ions being markers of the proton flow, we estimate the mass supply rate to coronal loops to be about 10\^{34} s\^{-1}.

We study the height variations of the sizes of chromospheric and transition-region features in a small coronal hole and the adjacent quiet Sun, considering images of the intensity, Doppler shift, and non-thermal motion of ultraviolet emission lines as measured by SUMER, together with the magnetic field as obtained by extrapolation from photospheric magnetograms. In order to estimate the characteristic sizes of the different features present in the chromosphere and transition region, we have calculated the autocorrelation function for the images as well as the corresponding extrapolated magnetic field at different heights. The HWHM of the autocorrelation function is considered to be the characteristic size of the feature shown in the corresponding image. Our results indicate that, in both the coronal hole and quiet Sun, the HWHM of the intensity image is larger than that of the images of Doppler-shift and non-thermal width at any given altitude. The HWHM of the intensity image is smaller in the chromosphere than in the TR, where the sizes of intensity features of lines at different temperatures are almost the same. But in the upper part of the transition region, the intensity size increases more strongly with temperature in the coronal hole than in the quiet Sun. We also studied the height variations of the HWHM of the magnetic field magnitude B and its component |Bz|, and found they are equal to each other at a certain height below 40 Mm in the coronal hole. The height variations of the HWHM of |Bz/B| seem to be consistent with the temperature variations of the intensity size. Our results suggest that coronal loops are much lower, and magnetic structures expand through the upper TR and lower corona much more strongly with height in the coronal hole than in the quiet Sun.